一类渔业经济模型的复杂动力学行为

Complex Dynamic Behavior of a Kind of Fishery Economic Model

摘要: 考虑到多个主体仅最大化个人利益的情形，本文建立了一个渔业经济模型。结合非线性动力学知识，利用单参数分岔图，吸引子，系统的间歇行为，以及吸引盆拓扑结构的变化分析了系统的动力学行为。较大的市场需求，可控的价格弹性和生物的内在增长率，较小的成本参数将会维持渔业资源的可持续发展，维护各主体的收益。

Abstract: Considering the situation that many agents only maximize their personal interests, this paper establishes a fishery economic model. Combined with the knowledge of nonlinear dynamics, the dynamic behavior of the system is analyzed by using one-parameter bifurcation diagram, attractors, intermittent behavior of the system, and the change of the topological structure of the basins. Larger market demand, controllable price elasticity and internal growth rate of organisms, and smaller cost parameters will maintain the sustainable development of fishery resources and maintain the return of each agent.

1. 引言

目前许多渔业资源正面临着衰退和不可持续的局面。对渔业经济学的研究，通常用动态模型来描述资源数量的变化 [1] - [6] 。Levhari和Mirman [1] 建立了关于渔获量的微分博弈的离散模型，获得了Cournot-Nash均衡的表达式。Bischi和Lamantia [3] 研究了一个古诺寡头垄断博弈的多主体渔业经济模型。顾恩国和陈宝香 [4] 研究了具有不同策略的多主体渔业经济模型。Raymond等 [5] 研究了渔业模型的动力学问题，不仅证明了共存平衡点的稳定性而且通过数值模拟证明了生物经济平衡。Zhou等人 [6] 基于梯度调整机制分析了在受其他参数影响下多个主体如何获取最大利润。

利用渔业经济模型预测和防止渔业资源的衰退和保持可持续发展是十分必要的，因此大多数学者引入了非线性系统的动态行为分析 [7] [8] [9] [10] [11] 。Din等 [7] 研究了一类离散奇异生物经济系统的稳定性与两种类型的分岔(即flip分叉和Neimark-Sacker分叉)。此外，吸引盆中还有其他复杂的动力学现象，如洞、吸引子共存、全局分叉等，一直是许多学者感兴趣的问题。比如Abraham和Laura等人 [8] 使用临界曲线方法直观地引入二维不可逆离散动力系统中的混沌和分叉。李慧 [9] 研究了鞍结分岔和两种间歇现象，并且利用不可逆映射理论和临界曲线等方法研究了模型的全局动力学。Liu等 [11] 利用临界曲线分析了吸引盆拓扑结构的变化以及吸引盆“缺口”的原因。此外，最大李亚普诺夫指数图，相图等，也是分析复杂动态行为的有力工具(见参考文献 [12] [13] [14] [15] )。

已有的文献大多采用线性需求函数，而线性水平一般不能很好地代表实际中较为复杂多变的情形，本文采用非线性需求函数使得模型更接近实际市场。并且将捕捞总量看成系统变量进行研究，结合非线性动力学知识，利用吸引子，系统的间歇行为，以及吸引盆拓扑结构的变化分析了系统的动力学行为。

2. 模型建立

假设有n个参与者和m个市场，其中 $n,m>1$ 。n个参与者收获鱼，每个参与者在m个市场上出售鱼。通过参考文献 [10] 的非线性需求函数，假定逆需求函数为：

$p_i=\frac{A_i-\left(x_{1i}+x_{2i}+\cdots +x_{ni}\right)}{B_i\left(x_{1i}+x_{2i}+\cdots +x_{ni}\right)}$ ， $i=1,2,\cdots ,m$ (1)

其中： $A_i$ 表示市场i的价格等于0时的最大瞬时需求； $B_i$ 表示需求相对于价格的敏感度的正参数； $x_{ki}$ 表示t时间内，参与者 $k\left(k=1,2,\cdots ,n\right)$ 收获并在市场i上出售的鱼的数量。

令 $X\left(t\right)$ 是在时间t时鱼群的总数量， $h_k\left(t\right)$ 是在时间t时参与者 $k\left(k=1,2,\cdots ,n\right)$ 收获(和出售)的鱼量。其中： $h_k\left(t\right)=x_{k1}\left(t\right)+x_{k2}\left(t\right)+\cdots +x_{km}\left(t\right)$ 。

假设参与者 $k\left(k=1,2,\cdots ,n\right)$ 的成本函数 [16] 为：

$C_k=c_k+\frac{{\gamma }_k{h}_k^2}{X}$ (2)

其中： $c_k$ 表示固定成本， ${\gamma }_k$ 表示技术参数且 ${\gamma }_k>0$ 。

假设 $s_i\left(t\right)=x_{1i}\left(t\right)+x_{2i}\left(t\right)+\cdots +x_{ni}\left(t\right)$ 为时间周期t内市场 $i\left(i=1,2,\cdots ,m\right)$ 供应(和销售)的鱼量；且

$H\left(t\right):=h_1\left(t\right)+h_2\left(t\right)+\cdots +h_n\left(t\right)=s_1\left(t\right)+s_2\left(t\right)+\cdots +s_m\left(t\right)$ (3)

假设渔民拥有完善的种群信息，在这种情况下，捕捞的鱼群的动态是由一维系统控制的：

$X\left(t+1\right)=X\left(t\right)\left(1+\alpha -\beta X\left(t\right)\right)-H\left(t\right)$ (4)

其中： $\alpha$ 表示内在增长率， $K=\frac{\alpha }{\beta }$ 表示环境承载力。

在完全不合作情形下，每个参与者 $k\left(k=1,2,\cdots ,n\right)$ 只是最大化他们自己的利益来决定他们的收获政策。则参与者 $k\left(k=1,2,\cdots ,n\right)$ 在时期t的利润函数为：

${\pi }_k\left(t\right)={\displaystyle \underset{i=1}{\overset{m}{\sum }}\frac{A_i-\left(x_{1i}+x_{2i}+\cdots +x_{ni}\right)}{B_i\left(x_{1i}+x_{2i}+\cdots +x_{ni}\right)}}{x}_{ki}-c_k-\frac{{\gamma }_k{h}_k{\left(t\right)}^2}{X}$ (5)

假设所有个体采取最佳策略反映，他们决定其捕捞策略使他们的利润最大化。每一时段t所有个体决定 $t+1$ 时段的捕捞量，必须形成其竞争对手的捕捞量的预期。则第 $k\left(k=1,2,\cdots ,n\right)$ 个个体在 $t+1$ 时段的捕捞量为：

$h_k\left(t+1\right)=\arg \max {\pi }_k\left(h_k\left(t\right),h_{-k}^e\left(t+1\right)\right)$ (6)

其中： $h_{-k}^e\left(t+1\right)$ 表示第 $k\left(k=1,2,\cdots ,n\right)$ 个个体关于其他个体的决策预期。

假设 $h_{-k}^e\left(t+1\right)=h_{-k}\left(t\right)$ ，则个体 $k\left(k=1,2,\cdots ,n\right)$ 的优化问题的解可以由以下方程(7)得到

$\frac{\partial {\pi }_k\left(t\right)}{\partial x_{ki}}=0$ (7)

则

$x_{ki}=\frac{A_i{s}_{i}X\left(t\right)}{\left(A_i-s_i\right)X\left(t\right)-2{\gamma }_k{h}_k\left(t\right)B_i{s}_i}$ (8)

对所有的i相加得到

$h_k={\displaystyle \underset{i=1}{\overset{m}{\sum }}{x}_{ki}}=\frac{{\displaystyle \underset{i=1}{\overset{m}{\sum }}{A}_i}H\left(t\right)X\left(t\right)}{\left({\displaystyle \underset{i=1}{\overset{m}{\sum }}{A}_i}-H\left(t\right)\right)X\left(t\right)-2{\gamma }_k{h}_k\left(t\right)H\left(t\right){\displaystyle \underset{i=1}{\overset{m}{\sum }}{B}_i}}$ (9)

令 ${\displaystyle \underset{i=1}{\overset{m}{\sum }}{A}_i}=a$ ， ${\displaystyle \underset{i=1}{\overset{m}{\sum }}{B}_i=b}$ ，则

$h_k=\frac{aHX\left(t\right)}{\left(a-H\right)X\left(t\right)-2{\gamma }_k{h}_k\left(t\right)Hb}$ (10)

对所有的 $k\left(k=1,2,\cdots ,n\right)$ 相加得到

$H=\frac{aX\left(t\right)\left(n-1\right)}{2bcH}-\frac{nX\left(t\right)}{2bc}$ (11)

其中： $c={\displaystyle \underset{k=1}{\overset{n}{\sum }}{\gamma }_k}$

结合方程(4)建立如下渔业资源储量与捕捞量之间相互作用的模型为：

$T:=\{\begin{array}{l}X\left(t+1\right)=X\left(t\right)\left(1+\alpha -\beta X\left(t\right)\right)-H\left(t\right)\\ H\left(t+1\right)=\frac{aX\left(t\right)\left(n-1\right)}{2bcH\left(t\right)}-\frac{nX\left(t\right)}{2bc}\end{array}$ (12)

为了使得模型有意义，对于 $t>0$ 限制在区域 $\left(X,H\right)\in R_{+}^2$ 研究。根据实际意义，所有时段总捕捞量小于总资源储量，即 $\left(X,H\right)\in \left\{\left(X,H\right)|0\le H\le X\right\}$ 。

3. 不动点的局部稳定性

系统T的平衡点满足

$X\left(t\right)\left(\alpha -\beta X\left(t\right)\right)=\frac{aX\left(t\right)\left(n-1\right)}{2bcH\left(t\right)}-\frac{nX\left(t\right)}{2bc}$ (13)

正不动点 $E\left(X^{*},H^{*}\right)$ 的坐标满足

$\alpha -\beta X^{*}\left(t\right)=\frac{\left(n-1\right)a}{2bc{H}^{*}\left(t\right)}-\frac{n}{2bc}$ (14)

化简(14)可以得到：

$2bc{\beta }^2{X}^{*}{}^3-\left(4bc\alpha \beta +n\beta \right)X^{*}{}^2+\left(2bc{\alpha }^2+n\alpha \right)X^{*}-a\left(n-1\right)=0$ (15)

记 $U=2bc{\beta }^2$ ， $V=-\left(4bc\alpha \beta +n\beta \right)$ ， $W=2bc{\alpha }^2+n\alpha$ ， $Q=-a\left(n-1\right)$ ， $X^{*}=Y^{*}-\frac{V}{3U}$ 。则(15)式简化为 $Y^{*}{}^3+p{Y}^{*}+q=0$ ，其中

$p=\frac{3UW-V^2}{3U^2},q=\frac{27U^{2}Q-9UVW+2V^3}{27U^3}$

记 $\Delta ={\left(\frac{q}{2}\right)}^2+{\left(\frac{p}{3}\right)}^3$ ，有以下的命题1。

命题1：

(1) 当 $\Delta >0$ 时，系统(13)有唯一的正平衡点 $\left(X^{*},H^{*}\right)$ 。

其中

$X^{*}=Y^{*}+\frac{2bc\alpha +n}{3bc\beta }$

$Y^{*}=\sqrt[3]{-\frac{q}{2}+\sqrt{\Delta }}+\sqrt[3]{-\frac{q}{2}-\sqrt{\Delta }}$

$H^{*}=\frac{\left(n-1\right)a}{2bc\left(\alpha -\beta X^{*}\right)+n}$

(2) 当 $\Delta =0$ 时，系统(13)有2个正平衡点 $\left(X_1^{*},H_1^{*}\right)$ 和 $\left(X_2^{*},H_2^{*}\right)$

其中：

$X_1^{*}=Y_1^{*}+\frac{2bc\alpha +n}{3bc\beta }$ ， $X_2^{*}=Y_2^{*}+\frac{2bc\alpha +n}{3bc\beta }$

$Y_1^{*}=2\sqrt[3]{-\frac{q}{2}}$ ， $Y_2^{*}=-\sqrt[3]{-\frac{q}{2}}$

$H_1^{*}=\frac{\left(n-1\right)a}{2bc\left(\alpha -\beta X_1^{*}\right)+n}$ ， $H_2^{*}=\frac{\left(n-1\right)a}{2bc\left(\alpha -\beta X_2^{*}\right)+n}$

(3) 当 $\Delta <0$ 时，系统(13)有3个正平衡点 $\left(X_1^{*},H_1^{*}\right)$ ， $\left(X_2^{*},H_2^{*}\right)$ 和 $\left(X_3^{*},H_3^{*}\right)$

其中

$X_1^{*}=Y_1^{*}+\frac{2bc\alpha +n}{3bc\beta }$ ， $X_2^{*}=Y_2^{*}+\frac{2bc\alpha +n}{3bc\beta }$ ， $X_3^{*}=Y_3^{*}+\frac{2bc\alpha +n}{3bc\beta }$

$Y_1^{*}=\sqrt[3]{-\frac{q}{2}+\sqrt{\Delta }}+\sqrt[3]{-\frac{q}{2}-\sqrt{\Delta }}$ ， $Y_2^{*}=\omega \sqrt[3]{-\frac{q}{2}+\sqrt{\Delta }}+{\omega }^2\sqrt[3]{-\frac{q}{2}-\sqrt{\Delta }}$ ， $Y_3^{*}={\omega }^2\sqrt[3]{-\frac{q}{2}+\sqrt{\Delta }}+\omega \sqrt[3]{-\frac{q}{2}-\sqrt{\Delta }}$

(这里 $\omega =\frac{-1+\sqrt{3}i}{2}$ )

$H_1^{*}=\frac{\left(n-1\right)a}{2bc\left(\alpha -\beta X_1^{*}\right)+n}$ ， $H_2^{*}=\frac{\left(n-1\right)a}{2bc\left(\alpha -\beta X_2^{*}\right)+n}$ ， $H_3^{*}=\frac{\left(n-1\right)a}{2bc\left(\alpha -\beta X_3^{*}\right)+n}$

不动点 $E\left(X^{*},H^{*}\right)$ 的局部稳定性依赖于系统(12)的Jacobi矩阵的特征值：

$J\left(E\right)=\left(\begin{array}{cc}1+\alpha -2\beta X^{*}& -1\\ \frac{a\left(n-1\right)}{2bc{H}^{*}}-\frac{n}{2bc}& -\frac{a{X}^{*}\left(n-1\right)}{2bc{\left(H^{*}\right)}^2}\end{array}\right)$ (16)

如果 $J\left(E\right)$ 的特征值的模小于1，则 $E\left(X^{*},H^{*}\right)$ 稳定。但是满足命题3.1条件的正平衡点 $E\left(X^{*},H^{*}\right)$ 存在时，很难解析地计算出Jacobi矩阵的特征值，所以严格分析稳定性是非常复杂的。因此有以下命题3.2：

命题2：如果满足条件：

$-1<\frac{a{H}^{*}-a{X}^{*}+H^{*}{}^{2}n-a{H}^{*}n-a\alpha X^{*}+an{X}^{*}+2a\beta X^{*}{}^2+a\alpha n{X}^{*}-2a\beta n{X}^{*}{}^2}{2bc{H}^{*}{}^2}<3$

则不动点 $E\left(X^{*},H^{*}\right)$ 是局部渐近稳定的。

证明：由于均衡点处的Jacobi矩阵的特征值很难获得具体的解析表达式，因此不能通过特征值在复平面内的具体位置来分析稳定性。然而可以通过数值分析不动点的稳定性，不动点的特征方程为：

$P\left(z\right)=z^2-\text{tr}J\left(E\right)z+\det J\left(E\right)$ (17)

其中

$\text{tr}J\left(E\right)=\alpha -2\beta X^{*}-\frac{a{X}^{*}\left(n-1\right)}{2bc{H}^{*}{}^2}+1$

$\det J\left(E\right)=-\frac{a{H}^{*}-a{X}^{*}+n{H}^{*}{}^2-an{H}^{*}-a\alpha X^{*}+an{X}^{*}+2a\beta X^{*}{}^2+a\alpha n{X}^{*}-2a\beta n{X}^{*}{}^2}{2bc{H}^{*}{}^2}$

不动点 $E\left(X^{*},H^{*}\right)$ 若满足Jury条件：

(i) $P\left(1\right)=1-\text{tr}J\left(E\right)+\det J\left(E\right)>0$ 。

(ii) $P\left(-1\right)=1+\text{tr}J\left(E\right)+\det J\left(E\right)>0$ 。

(iii) $1-\det J\left(E\right)>0$ 。

则不动点 $E\left(X^{*},H^{*}\right)$ 是局部稳定的，即

(i) $2\beta X^{*}-\alpha -\eta +\frac{a{X}^{*}\left(n-1\right)}{2bc{H}^{*}{}^2}>0$ 。

(ii) $\alpha -2\beta X^{*}-\eta -\frac{a{X}^{*}\left(n-1\right)}{2bc{H}^{*}{}^2}+2>0$ 。

(iii) $\eta +1>0$ 。

其中

$\eta =\frac{a{H}^{*}-a{X}^{*}+H^{*}{}^{2}n-a{H}^{*}n-a\alpha X^{*}+an{X}^{*}+2a\beta X^{*}{}^2+a\alpha n{X}^{*}-2a\beta n{X}^{*}{}^2}{2bc{H}^{*}{}^2}$

则有

(i') $\alpha -2\beta X^{*}<\frac{a{X}^{*}\left(n-1\right)}{2bc{H}^{*}{}^2}-\eta$ 。

(ii') $\alpha -2\beta X^{*}>\frac{a{X}^{*}\left(n-1\right)}{2bc{H}^{*}{}^2}+\eta -2$ 。

(iii') $\eta >-1$ 。

将(iii')代入(i')和(ii')，可以得到：

$\alpha -2\beta X^{*}>\frac{a{X}^{*}\left(n-1\right)}{2bc{H}^{*}{}^2}-1-2=\frac{a{X}^{*}\left(n-1\right)}{2bc{H}^{*}{}^2}-3$

所以有

$\frac{a{X}^{*}\left(n-1\right)}{2bc{H}^{*}{}^2}-3<\alpha -2\beta X^{*}<\frac{a{X}^{*}\left(n-1\right)}{2bc{H}^{*}{}^2}-\eta$

即当

$-1<\frac{a{H}^{*}-a{X}^{*}+H^{*}{}^{2}n-a{H}^{*}n-a\alpha X^{*}+an{X}^{*}+2a\beta X^{*}{}^2+a\alpha n{X}^{*}-2a\beta n{X}^{*}{}^2}{2bc{H}^{*}{}^2}<3$

不动点 $E\left(X^{*},H^{*}\right)$ 是局部渐近稳定的。

图1显示了b-c平面上的纳什均衡点的稳定区域。参数为 $\alpha =0.0305$ ， $\beta =0.662$ ， $a=0.3$ 。稳定区域的大小在很大程度上受瞬时市场需求和内在增长率的影响。如图1(b)和图1(c)所示，稳定区域的范围随着 $\alpha$ 的增大而变小，而随着a的增大而增大。

(a) $a=0.3,\alpha =0.0305$ (b) $a=0.8,\alpha =0.0305$ (c) $a=0.8,\alpha =0.8$

Figure 1. Stability regions of the Nash equilibrium point on the b-c plane

图1. b-c平面上的纳什均衡点的稳定域

4. 数值模拟

上一节使用Jury准则分析了定点的稳定性；本节将利用数值模拟研究参数对系统稳定性的影响。我们研究了两个代理在同一市场上运营的情况。假设固定成本为 $c_k=0$ ，并且两个代理都采用类似的方法捕获鱼类生物并在同一市场上销售。

系统关于参数b的一维分岔图如图2(a)所示，可以看出随着参数b的不断增加，系统先发生Neimark-Sacker分岔，然后发生逆倍周期分岔之后进入混沌。参数选取与图2(a)相同，但是将初值选取为(1.1, 0.93)，得到系统T关于参数b的单参数分岔图，如图2(b)。对比图2(a)和图2(b)，发现在跳跃处两个分岔图具有不同的状态，可以看出系统可能具有多稳态结构。当初始的渔业生物资源存储量和捕捞量较小时，系统在(1.79, 1.81)之间的复杂动力学行为是比较简单的周期行为，但是当初始值改变，初始的渔业生物量和捕捞量变大时，跳跃处的周期行为将变得复杂。利用吸引盆分析单参数分岔图发生跳跃的原因是因为吸引子在这期间共存。值得注意的是，系统在 $b=1.7913$ 时，单参数分岔有明显的跳跃，这些地方都可能有共存吸引子。这说明价格敏感度系数应该保持在一定的范围内才会使渔业经济达到稳定，否则会让渔业经济失去平衡，陷入不可逆的混沌状态。

为了研究单参数分岔图中 $b=1.7913$ 时的跳跃现象，选取与图2(a)相同的参数集，得到关于b的单参数分岔图跳跃处的吸引盆和吸引子，如图2(c)~图2(f)。当参数 $b=1.791$ 时，系统T有一个5-环，在图中用黑色表示，用浅蓝色标记为它的吸引盆，如图2(c)。调整价格敏感参数，当 $b=1.7913$ 时，系统中有两种吸引子共存，一种是用黑色标记5-环，另一种是用红色标记的30-环，其吸引盆用黄色标记，如图2(d)。继续增大价格敏感参数b，红色的30-环会经过flip分岔变成60-环，……，之后进入混沌，如图2(e)，当 $b=1.806$ 时变为周期为15的混沌吸引子，与此同时，5-环也从最初的不动点变为5-周期的不变环。随着

Figure 2. (a) When the initial conditions are set to $X=0.62$ and $H=0.61$ , the one-dimensional bifurcation diagram with respect to b; (b) When the initial conditions are set to $X=1.10$ and $H=0.93$ , the one-dimensional bifurcation diagram with respect to b; (c) $b=1.79100$ , 5-cycle and its basin; (d) $b=1.79130$ , 5-cycle coexists with 30-cycle and their basins; (e) $b=1.80600$ , 5-cycle coexists with period-15 chaotic attractor; (f) $b=1.80700$ , only period-5 invariant cycle remains in the system

图2. (a) 初始条件选择为(0.62, 0.61)时，关于b的单参数分岔图；(b) 初始条件选择为(1.1, 0.93)时，关于b的单参数分岔图；(c) $b=1.79100$ 时，5-环及其对应的吸引盆；(d) $b=1.79130$ ，5-环与30-环及其对应的吸引盆；(e) $b=1.80600$ ，5-周期不变环与15-周期混沌吸引子及其对应的吸引盆；(f) $b=1.80700$ ，5-周期不变环和它的吸引盆

价格敏感参数b不断增加红色的混沌吸引子与黑色的不变环也变得越来越大，各自对应的吸引盆范围也变大了。直到b增加到1.807时，由于15-周期的混沌吸引子与其吸引盆边界发生接触15-环及其黄色的吸引域都消失，最后系统只有5-周期的不变环，如图2(f)所示。价格敏感参数对系统稳定性有一定的影响，当价格敏感参数控制在适当的范围，才会让渔业经济系统保持在有序的状态，渔类生物资源存储量和捕捞量也会是可持续发展的。

如图3(a)所示，当成本参数c在不断变大时，系统会从最初的混沌状态经由flip分岔，再次进入混沌状态。当参数 $c\in \left(1.425,1.6\right)$ 和 $c\in \left(1.935,2.175\right)$ 时保持周期稳定，成本参数c不断增大，会让渔业经济系统在 $c=2.2$ 时，不可逆的进入一个无序的状态，并且会很难预测出此渔业经济系统未来的走向。同时，系统在 $c\approx 1.4257$ 会突然从混沌状态分别进入10-周期状态。接下来，将分析参数c在1.4257附近变化时系统的动力学行为。

如图3(b)所示， $c=1.4255$ 时，系统只存在一个混沌吸引子。图3(c)是图3(b)对应的时间序列图，这是PM-III型间歇混沌发生前的时间序列，表现为周期运动与混沌运动相互交替。 $c=1.4257$ 是PM-III型间歇混沌发生的临界值。当 $c=1.426$ ，如图3(d)是PM-III型间歇混沌发生后情形，系统经历短暂混沌运

Figure 3. (a) The one-dimensional bifurcation diagram with respect to bifurcation parameter c; (b) $c=1.42570$ , the attractor of the system is a chaotic attractor; (c) Time series diagram corresponding to (b); (d) $c=1.42600$ , the attractor of the system is a 10-cycle

图3. (a) 关于分岔参数c的单参数分岔图；(b) $c=1.42570$ ，系统只有一片混沌吸引子存在；(c) 图(b)所对应的时间序列图；(d) $c=1.42600$ ，系统只有一个吸引的10-环

动之后进入周期-10状态。间歇混沌现象意味着系统具有自我调节能力，系统在稳定与混沌之间不断循环，但随着时间的推移，系统失去了长期稳定性，这意味着这种自我调节能力无法使系统保持长期稳定发展。

5. 结语

基于非线性需求函数，本文建立了多主体渔业经济模型，其中多主体之间以自身利益最大化作为自己的收获决策。利用系统的雅可比矩阵和Jury判据讨论了不动点的存在性和局部稳定性。分析了系统的多稳态性，发现较大的市场需求，可控的价格弹性和生物的内在增长率，较小的成本参数将会维持渔业资源的可持续发展，各个主体的可观收益。并且发现由于吸引子的共存会让单参数分岔图有跳跃现象。此外，间歇混沌现象意味着系统具有自调节能力，系统不断在稳定以及混乱之间循环往复。通过分析吸引盆的形状和大小变化，将有效防止资源衰竭，维护主体利益。

NOTES

^*通讯作者。