《捕鱼业持续收获效益问题》由会员分享，可在线阅读，更多相关《捕鱼业持续收获效益问题（14页珍藏版）》请在金锄头文库上搜索。
1、捕鱼业持续收获的效益问题天津师范大学数学科学学朱保军摘要为了保护人类赖以生存的自然环境，可再生资源(渔业)的开发必须适度。而在社会经济生活中，我们要使商业活动在一段时期内达到最大收益，因此我们要合理的开发资源，这时，我们不仅要考虑商业活动的当前经济效益，还要考虑生态效益及由此产生的对整体经济效益的影响。本文就是对渔业这类可再生资源的开发问题进行研究，利用相关的数学软件matlab进行求解。针对问题一，建立产量模型，在该模型中建立捕捞情况下鱼场鱼量方程x(t).一一一.一和捕捞万程E(t),对万程乂(。和(讨论其平衡点x0和稳定性F(X。)如下：Xo=N(1-),F(X0)=Er0r针对问题二，
2、对产量模型进行延伸可求出最大的持续产量情况下的稳定平衡点X0，最大持续产量hm和保持渔场鱼量稳定在X0的捕捞率E。再建立效益模型，可求出利润R(E),最大捕强度Er,最大利润下的渔场稳定鱼量xR及单位时间的持续产量hR。最后，建立捕捞过度模型，可求出盲目捕捞下的渔场稳定鱼量为Xs。稳定平衡点持续产量捕捞强度产量模型NX0=2rNhm=4*rE=一2效益模型NcXR-+22prNc2、hR_/1224PNJEr=2f,1-捕捞过度cxs=PEs=rA、1-PNJ关键字：产量模型效益模型捕捞过度matlab、问题提出为了保护人类赖以生存的自然环境，可再生资源(渔业)的开发必须适度。一种合理、简化的
3、策略是，在实现可持续收获的前提下，追求最大产量或最佳效益。设渔场鱼量方程为：dx=x(tx=rx(1-)-ExdtN.其中x(t)表示t时刻鱼塘中的鱼数，r表示鱼的自然增长率，N为鱼塘中最大可容纳的鱼量，E为捕捞率。.问题一，建立关于E(t)的方程，求x(t),E(t)的平衡点并讨论其稳定性。问题二，将所得结果与捕捞的效益模型与过度模型进行比较。综上所诉，本文要建立在捕捞情况下渔场鱼量遵从的方程，分析鱼量稳定的条件，并且在稳定的前提下讨论如何控制捕捞使持续产量或经济效益达到最大。最后研究所谓捕捞过度问题。二、模型假设1 .这种鱼在一年内的任何时间都会发生自然死亡，即死亡是一个连续的过程;2 .
4、捕捞也是一个连续的过程，不是在某一时刻突然发生；3 .鱼群的死亡率已考虑种群的相互竞争及环境等因素；三、字符说明x(t):时刻t时,渔场中的鱼量r:鱼的固有增长率N:环境容许的最大鱼量f(x)：单位时间的增长量h(x):单位时间的捕捞量F(x):F(x)=f(x)-h(x)X0、Xi:为方程F(x)的解%:稳止平衡点hm:最大持续产量E*:保持渔场鱼量稳定在*的捕捞率p:鱼的销售单价C：单位捕捞率的费用T：单位时间的收人S:单位时间的支出R(E):单位时间的利润Er：最大捕强度冷：最大利润下的渔场稳定鱼量hR:单位时间的持续产量Es：R(E)=0的解xs：盲目捕捞情况下的渔场稳定鱼量四、问题分
5、析问题一分析：主要是列出x(t)、E(t)的方程，可知是一阶微分方程的平衡点及稳定性分析，当x(t)=0时，X=Xg为方程x(t)的稳定点，若F(x)0,则x0对于万程是不稳定的。同理对于E(tW勺方程。问题二分析：列出产量模型主方程x(t),分析出最大的持续产量情况下的稳定平衡点*,最大持续产量hm和保持渔场鱼量稳定在*的捕捞率E*；再列出效益模型R(E),分析出利润R(E),最大捕强度Er,最大利润下的渔场稳定鱼量xR及单位时间的持续产量hR；最后，列出捕捞过度模型Es，分析出盲目捕捞下的渔场稳定鱼量为xs及其cp2c时的情况。sNNN五、模型建立问题一:产量模型：在无捕捞条件下x(t)的
6、增长服从logistic规律即x(t)=f(x)=rx(1-9(1)单位时间的捕捞量与渔场鱼量x(t)成正比，比例常数E表示单位时间捕捞率，即E(t尸cx(t)。h(x尸Ex(2)根据以上并记Fx=fX-hX得到捕捞情况下渔场鱼量满足的方程x(t)=F(x)=rx(1N)Ex(3)我们不需要解方程(3)以得到x(t)的动态变化，只希望知道渔场的稳定鱼量和保持稳定的条件，并由此来确定最大持续产量。为此可以直接求方程(3)的平衡点并分析其稳定性。令xFx=rx(1-一)Ex=0N得到两个平衡点xo=N(1-5),%=0(4)r不难算出Fx=E-r,Fxi=r-E所以若Er(5)有F(M)0,故几点
7、稳定，为点不稳定；Er,则结果正好相反。又因为E(t尸cx(t),故E(t)同x(t)的平衡点及稳定性。问题二：产量模型的延伸：进一步讨论渔场鱼量稳定在飞的前提下，如何控制捕捞强度E使持续产量最大的问题。用图解法非常简单地得到结果。作出抛物线y=f(x)和直线y=h(x)=Ex,如上图。注意到y=f(x)在原点的切线为y=rx,所以在条件(5)下y=Ex必与y=f(x)有交点P,P的横坐标就是平衡点x0P点的纵坐标h为稳定条件下单位时间的持续产量，由图立刻知道，当y=Ex与y=f(x)在抛物线顶点p相交时可获得最大的持续产量，此时的稳定平衡点为(6)且单位时间的最大持续产量为而由(4)式不难算
8、出保持渔场鱼量稳定在*的捕捞率为(8)E.=r综上所述，产量模型的结论是将捕捞率控制在固有增长率r的一半，更简单一些。可以说使渔场鱼量保持在最大鱼量N的一半时，能够获得最大的持续产里。效益模型：从经济角度看不应追求产量最大，而应考虑效益最佳T=ph(x)=pEx,S=cE(9)(11)(10)Er1-(12)单位时间利润问R=T-S=pEx-cE在稳定条件x=%下，以(4)式代入(10)式，得RE=T(E)-S(E)=pNE(1-)-cEr用微分法容易求出使利润达到最大的捕捞强度为hR为(13)hR=rxR(1-XRrN2c-22pN(14)将Er代入(4)式，可得最大利润下的渔场稳定量xr及
9、单位时间的持续产量NcxR=T-22p将(12)(14)式与产量模型中的(6)(8)式相比较可以看出，在最大效益原则下捕捞率和持续产量有所减少，而渔场应保持的稳定鱼量有所增加，并且减少或增加的部分随着捕捞成本c增长而变大，随着销售价格p的增长而变小，这结果是符合实际情况的。捕捞过度：上面的效益模型是以计划捕捞为基础的。如果渔场向众多盲目的经营者开放，那么即使只有微薄的利润，经营者也会蜂拥而去，这种情况称为盲目捕捞。这种捕捞方式将导致捕捞过度，下面讨论这个模型。(11)式给出了利润与捕捞强度的关系R(E),令R(E)=0的解为Es,可得Es=r1-I(15)PN；当E0,盲目的经营者们会加大捕捞
10、强度；若EaEs时，利润R(E)0)是(16)cPN即售价大于成本。并且由15)式可知，成本越低，售价越高，则Es越大。将(15)代入(4)式，得到盲目捕捞下的渔场稳定鱼量为(17)cXs=P入完全由成本-价格比决定，随着价格的上升和成本的下降，力将迅速减少，出现捕捞过度。比较（12）和（15）式知，Es=2Er,即盲目捕捞强度比最大效益下捕捞强度大一倍。cc*从（15）式和上图还可以得到，当彳52小时，（Er）Es2：时，EsE*,如图中Es2,称生态学捕捞过度。N综上所诉得稳定平衡点持续产量捕捞强度产量模型NX0一万,rNhm-，4e*=-2效益模型NcXR=c丁22p7f2,rNcK*1
11、Er=2/71-hR12.24PN)捕捞过度cXs=p1Es=rA、1-PNJ六、模型改进为了研究渔业的产量、效益及捕捞过度的问道，首先在对鱼的自然增长和捕捞情况的合理假设下，建立渔场鱼量的基本方程（3）,并利用平衡点稳定性分析确定了保持渔场鱼量稳定的条件。产量、效益和捕捞过度的3个模型在稳定的前提下步步深入，数学推导过程十分简单，却得到了在定性关系上与实际情况完全符合的结果。但如果改变对鱼的自然增长和人工捕捞的假设，模型及结果将随之变化七、参考文献11赵静，但琦，数学建模与数学实验（第3版）高等教育出版社2姜启源、谢金星、叶俊，数学模型（第三版）高等教育出版社，2003【3】刘来福，最优捕鱼
12、策略问题答案评述，数学的实践与认识，1997年1月八、附录图一clc;clear;ezplot(y=-xA2+4*x,0,4,0,5)text(3.5,0.4,y=f(x)text(4,-0.3,N)holdonx1=0:0.01:2y1=4*x1plot(x1,y1,r:)text(1.2,4.5,y=rx)holdonx2=0:0,01:3y2=2*x2plot(x2,y2,b)text(2.1,4.1,P*)holdonx3=0:0.01:4y3=x3plot(x3,y3,b)text(3.15,3,P)text(3.3,3.1,y=h(x)=Ex)holdonx4=3y4=0:0.05
13、:3plot(x4,y4,r:)text(3.05,-0.2,x0)holdonx5=2y5=0:0.05:4plot(x5,y5,r:)text(2.1,-0.3,x0*=N/2)holdony6=3plot(x6,y6,r:)text(-0.2,3,h)holdonx7=0:0.05:2y7=4plot(x7,y7,r:)text(-0.3,4,hm)clc;clear;ezplot(y=-xA2+4*x,0,4,0,5)text(4,-0.3,r)text(3.5,1,T(E)holdonx1=0:0.01:2y1=4*x1plot(x1,y1,r:)text(1.2,4.5,pNE)holdonx2=0:0.01:3y2=3*x2plot(x2,y2,b)text(1.6,4.3,S(E)=cE)holdonx3=0:0.01:4y3=x3plot(x3,y3,b)text(3.1,3,S(E)=cE)holdonx4=3y4=0:0.05:3plot(x4,y4,r:)text(3.05,-0.2,Es2)hold
