最优捕鱼策略

1、精选优质文档-倾情为你奉上最优捕鱼策略模型摘要本文涉及的问题是渔业资源可持续发展，即在我国一定渔场，在一段时间内，如何实现最大收获量的问题；同时保证渔场能实现稳定生产。我们的解题思路就是：以渔场生产过程中的两个相互制约的因素，分别是年捕捞量与再生产能力，从而确定最优发展策略:用微分方程分别描述各龄鱼1群数量虽时间变化的规律，并在此基础上确定总效益即总收获量为目标函数，以渔场的可持续捕获为约束条件，分别对长期生产和固定期生产建立规划模型。针对问题一：通过对4龄鱼在年末的不同状态(全部死亡；仍为4龄鱼)的考虑，从可可持续捕获条件出发，分别建立2个模型。最后求解在计算机上实现。针对问题二：确定一个整

2、体效益函数，综合考虑年捕捞能力和年再生产能力，用计算机数值解法进行搜索逐年确定各年的最优策略，从而得出五年的总最优策略。1. 先假设每年捕捞量强度相等，建立一个简单模型；2. 再假设每年捕捞强度不相等，建立一个复杂模型；3. 最后给出鱼群生产能力破还不大的含义（即鱼群减少率的上限）在它的约束之下再建立一个模型。关键词：微分方程；捕捞强度；再生产能力；规划模型一、 问题的重述为了保护人类赖以生存的自然环境，可再生资源（如渔业、林业资源）的开发必须适度以可持续发展。那么一种合理、简化的策略是在可持续收获的前提下，追求最大产量或最大效益。要求研究的问题是：在一段期间中，对某种鱼的最优捕捞策略。1.1

3、鱼的情况假设这种鱼分为4个年龄组：1，2，3，4龄鱼。各年龄组每条鱼的平均重量分别为5.07,11.55,17.86,22.99（克）；各年龄组的鱼自然死亡率均为0.8（1/年）；这种鱼为季节性集中产卵繁殖，平均每条4龄鱼的产量为1.109×105(个)，3龄鱼的产卵量为这个数的一半，2龄鱼和1龄鱼不产卵，产卵和孵化期为每年的最后4个月；卵孵化并成活为1龄鱼，成活率（1龄鱼条数与产卵总量n之比）为1.22×1011/（1.22×1011+n）。具体数据如下表：imi(g)r(1/年)ui(个/条)15.070.80211.550.80317.860.80.5545

4、×105422.990.80.5545×105其中，i表示i龄鱼，mi表示i龄鱼的重量，r表示i龄鱼的自然死亡率，ui表示平均每条i龄鱼的产卵量。又渔业管理部门规定，每年只允许在产卵孵化期前的8个月内进行捕捞作业。如果每年投入的捕捞能力（如渔船数、下网次数等）固定不变，这时单位时间捕捞量将于i成正比，比例系数称捕捞强度系数ki,通常使用13mm网眼的拉网，这种网只能捕捞3、4龄鱼，其两个捕捞强度系数之比为：k3:k4=0.42:1，k1=k2=0.渔业上称这种方式为固定努力量捕捞。1.2问题1) 建立数学模型分析如何实现可持续捕捞（即每年开始捕捞时渔场中各年龄组鱼条数不变）

5、，且在此前提下得到最高的年收获量（捕捞总量）。2) 某渔业公司承保这种鱼的捕捞业务5年，合同要求5年后鱼群的生产能力不能受太大的破坏。已知承保时各年龄组鱼群的数量分别为,如果122,29.7,10.1,3.29（*109条），如果仍用固定努力量的捕捞方式，该公司采取怎样的策略才能使总收获量最高。二、 模型假设1. 假设只考虑某种鱼群的繁殖和捕捞，期间不考虑鱼群的迁入和迁出，即假设大规模鱼群数量增长随时间是连续变化的。2. 根据模型已知条件，假设鱼群在一年的任何时间都会发生自然死亡，在一年的后四个月都会发生产卵。3. 4龄鱼的数量在第4年仍存活的数量占总的数量比例很小，这里可以假设其全部死亡，令

6、其退出完毕。4. 不考虑环境的影响，各年龄组的自然死亡率为0.8（1/年）。三、 符号说明a年平均固定死亡率，单位：1/年sj(t):t时刻j龄鱼的捕捞总重量，j=3或4t时间，t0,1H年总收获量，即捕捞总重量jj龄鱼，j=1,2,3,4i年数i0,5gjj龄鱼每条鱼的平均重量，单位：克si,j(t)T时刻j龄鱼第i年的捕捞重量，j=3,4k年平均捕捞率，单位：1/年ni第i年的产卵量n每年产卵量xi,j(t)j龄鱼在第i年时刻t的数量xj(t)j龄鱼在时刻t的数量Hi第i年总收获量，即捕捞总重量x0(j)j龄鱼在年初的数量x1(j)j龄鱼在年末的数量四、 问题的分析(一) 对于问题一的分析

7、1.对死亡率a的理解：我们定义平均死亡率a是单位时间鱼群死亡数量与现有鱼群数量的正比例系数，由假设可知，它是一个环境等其他因素无关的常数。由于鱼群的数量是连续变化的，而1、2龄鱼全年以上及3、4龄鱼在后4个月的数量只与死亡率有关，与其他因素无关，设鱼群量为x，则在时间t,t+t内，鱼群数量的减少=鱼群死亡数量，即xt+t-xt=-ax(t)t，xt+t-x(t)t =-ax(t).当t0时得：dx(t)dt=-ax(t), (1.1)2.对于捕捞强度系数k的理解：题目告诉我们，捕捞强度系数K一定，且只在捕捞期内（即每年的前8个月）捕捞3、4龄鱼，因此只会影响3、4龄鱼群数量，而不会影响其他的鱼

8、群数量。我们可以看3、4龄鱼鱼群的数量在捕捞期内不仅与k有关，而且还与死亡率a有关，类似第一点的分析，可以得到dx(t)dt=-axt-0.42kx(t) (1.2) dx(t)dt=-axt-kx(t) (1.3)其中x(t)表示3、4龄鱼的鱼群数量，t表示每年的前8个月，即t0,23.3.对于持续捕获的理解：随着时间的推移，各年龄组的鱼数量必将发生变化，但持续捕获要求每年开始捕捞时渔场中各年龄组鱼群条数不变，再根据鱼群的生长规律，我们可以得到关系式：上一年龄组鱼群年底的数量=下一年龄组鱼群年初的数量（1龄鱼除外），即x1j=x0(j+1) (1.4)4.对成活率m的应用又假设知，此种鱼在每

9、年的8月底一次产卵完毕，又已知3、4龄鱼每条产卵的个数，因此可将每年的产卵量n表示为：n=1.109×105×0.5x323+x423 (1.5)又已知成活率m=1.22×10111.22×1011+n (1.6)产卵量×成活率=1龄鱼每年年初的数量，即n×m=x01 （1.7）5.对最高收获量的描述根据第2点的分析，在t时刻的捕捞重量s(t) =3龄鱼捕捞重量s3(t)+4龄鱼捕捞重量s4(t).而s3(t)=0.42kx3(t)g3, s4(t)=kx4(t)g4，则st=0.42kx3(t)g3+ kx4(t)g4 （1.8）由

10、于捕捞被看成连续的作业，因此捕捞总收获量即年收获量可以用t时刻的捕捞重量s(t)关于t在捕捞期内的积分来刻画H=023stdt, (1.9)要求最高的年收获量即求H的最大值。6.对4龄鱼在年末情形的两个假设：（1）认为4龄鱼在年末与鱼群总数量相比十分微小，它们既不产卵又不会被捕捞，可以将它们忽略不计，令其退出系统。（2）近似的认为在年末未死亡的4龄鱼的各个特征（例如：重量、产卵个数等）均不发生改变，即仍回到4龄鱼组中。7.模型建立大纲（1）以第6点的第一个假设为基础，建立一个简单的模型I，其实实质上是联立以上分析的几个方程成为一个方程组，其中x04=x13.（2）以第6点的第二个假设为基础，改

11、进模型II，使得方程组中的一个方程x04=x13被方程x04=x13+x1(4)代替即可。（二）对于问题二的分析1.问题一的相似之处由于对各年龄组鱼群数量起到影响作用的各因素（例如：平均死亡率、成活率、捕捞期等）不变，因此，在每年各年龄组的鱼群数量变化情况与问题一类似。2.与问题一的不同之处，主要区别在于（1）问题一要求持续捕获，问题二只要求鱼量不受到太大的破坏，不限制各年龄组年初鱼群的数量，因此作为约束条件的方程组中各年龄组的鱼量肯定与年数有关，而不像问题一是常量。（2）问题一中的各变量呈周期变化，因此，只要考虑一个周期的变化情况即可，而问题二则不同，其各年的初始值在变化，因此，要考虑每一年

12、的捕获量，再讲5年求和，得到一个目标函数。综合以上两点，可以得到一个优化问题。3.根据优化问题我们又提出了三个模型模型I：简化使得每一年的捕捞强度系数相同，化为一元函数最优值的求解问题。模型II：考虑每一年的捕捞强度系数不同，得到一个多元函数最优值的求解问题。五、模型的建立问题一模型I：假设4龄鱼在年底推出系统和连续捕获前提 t0,1下如何得到高年收获量。由问题分析可知，可得以下优化问题：maxH=023s(t）dt,s(t)见问题分析， （2.1）dx1(t)dt=-ax1t t0,1dx2(t)dt=-ax2t t0,1dx3(t)dt=-ax3t-0.42k x3 t0,23 dx3(t

13、)dt=-ax3t t23,1dx4(t)dt=-ax4t-kx4 t0,23dx4(t)dt=-ax4t t23,1 解I得到x0(3)=1.22×1011e2a-1.22×10111.109×1050.5e-2×0.42k3+e-2-2×0.423k-23a (2.2)图：年度总捕获量H随捕捞强度K的变化曲线求解过程见附录。将（2.2）式代入目标函数（2.1）式中得到H关于k的医院函数，再利用一维搜索法求一元函数的最小值的方法上机求得H的最大值，k=17.4664，H=3.886e+011（吨），即k=17.4664时，H取最大。模型II：

14、进一步假设4龄鱼在年末的特征不变，仍当做4龄鱼，则在持续捕捞的情况下，求得最大捕获量。此模型类似模型I也可得优化问题，区别仅在式中，应改为x04=x13+x1(4)同理得x3=1.22×1011e2a-1.22×10111.109×1050.5e-23(0.42k-a)+e-2-2×0.423k-53a1-e-2k3a重复模型II的步骤解得k=17.4664，H=3.886e+011（吨）。问题二：根据题意，我们既要保证五年后鱼群的生产能力不能受到太大破坏，又要使公司总收获量最高。但是如果同时满足这两个条件，总的来说很难把握两者之间的度。因此，我们采用首

15、先使捕鱼收获量最大，其次再看其破坏程度是否影响鱼的生长，如果破坏程度能保证不影响鱼的生长，则使用该方法。即五年合同到期后鱼群尽可能接近可持续鱼的情况下使捕捞量达到最大。根据问题的分析可以得到以下优化方程： maxH=maxi=15(s3i+sei) dxi,1(t)dt=-axi,1t t0,1dxi,2(t)dt=-axi,2t t0,1dxi,3(t)dt=-axi,3t-0.42ki xi,3 t0,23 dxi,3(t)dt=-axi,3t t23,1dxi,4(t)dt=-axi,4t-ki xi,4 t0,23dxi,4(t)dt=-axi,4t t23,1 我们将此问题分成了两种

16、情况进行处理。情况一：假定五年中捕捞强度k一直保持不变，即ki=k，且第6年初的各年龄鱼数与平衡捕捞时的年初值x0(j)(j=1,2,3,4)的偏离不超过x0(j)的倍。Max f(k)=i=1517.86*0.420.42k+0.8x3(i-1)1-e-2(0.8+0.42k)3s.t x2i=x1(i-1)e-0.8x3i=x2(i-1)e-0.8x4i=x4i-1e-0.8+2k3+x3(i-1)e-0.8+0.84k3N3=1.5*1.109*1050.8x3(i-1)e-0.84k3e-1.63-e-0.8N4=3*1.109*1050.8x4(i-1)e-2k3e-1.63-e-0

17、.8x1i=1.22*1011(N3+N4)1.22*1011+(N3+N4)xi5-x0(i)x0(i)用Matlab进行求解（代码见附录），结果可得:K=17.1446 N1=0.9989 N2=0.9987 N3=0.9997 N4=0.9500此时，最高总收获量为1.6011e9吨。 情况二：假设每年的捕捞强度系数都不相同，即ki与i有关，且ki相互独立。第6年初的各年龄鱼数与平衡捕捞时的年初值x0(j)(j=1,2,3,4)的偏离不超过x0(j)的倍。Max f(k)=i=1517.86*0.420.42ki+0.8x3(i-1)1-e-2(0.8+0.42ki)3s.t x2i=x

18、1(i-1)e-0.8x3i=x2(i-1)e-0.8x4i=x4i-1e-0.8+2ki3+x3(i-1)e-0.8+0.84ki3N3=1.5*1.109*1050.8x3(i-1)e-0.84ki3e-1.63-e-0.8N4=3*1.109*1050.8x4(i-1)e-2ki3e-1.63-e-0.8x1i=1.22*1011(N3+N4)1.22*1011+(N3+N4)xi5-x0(i)x0(i)用Matlab进行求解（代码见附录），结果可得：k1=13.0018 k2=13.9406 k3=21.6540 k4=21.5971 k5=17.2248N1=0.9989 N2=0.

19、9500 N3=0.9500 N4=0.9500此时，最高总收获量为1.6290e9吨。六、模型的评价与改进6.1 模型的优点1）、每个模型的数值都可以通过Matlab计算解得，且解出来的值与实际情况相吻合，都可以用一般常识来解释，解的稳定性较好。2）、本题用到的一些数学方法及知识，比如说常微分方程求解、一元函数极值求解等都比较简单。3）、本模型通过了解鱼群数量关系，运用连续性随机变量的方法并建立微分方程解决问题，使鱼群死亡率在模型中更接近实际情况，实用性较强。4）、本模型的建立实现了可持续捕捞，在获得最高收获量时还能有效地保持生产平衡。5）、没有单一性，可以预测情况相差不大的其他物种之间的关

20、系和发展趋势。6.2 模型的缺点1）本题没有考虑社会资源的限制，这与实际情况不相符。2）没有考虑禁捕期内的3、4龄鱼的产卵过程的连续性，造成了一定的误差。3）忽略了一些其他因素对渔场的影响，得到的结果与实际有一定的差距。5.2 模型的改进1）该模型虽然提出了产卵与孵化过程应该是连续的过程，但在求解过程中实质上还是利用了离散的过程代替了连续过程。因此可以在此方向上进行改进，可以分成两种情况：1、每个时刻t，3、4龄鱼都可以固定的概率产卵；2、每个时期t，3、4龄鱼都以变概率产卵，这样又可以引入两个模型。2）对4龄鱼在每年年末的情况还可以再分一种情况，设其为与前面无关的第5类鱼，再根据有关的资料查

21、出产卵量与年龄的关系，得到一个数学表达式，对模型进行适当的改进。七、参考文献1孙祥,徐流美,吴清编，MATLAB7.0基础教程清华大学出版社 2005. 2吴振奎 王全文主编运筹学中国人民大学出版社 2005.3罗君，刘鹏周，鸣炜，最优捕鱼策略模型，数学的时间与认识，27卷：51-27,1997.4朱元鼎，王文浜，中国动物图谱（鱼类），科学出版社，北京，1973.5E.C.皮洛，数学生态学引论，科学出版社，北京，1978.附录微分方程求解过程：dx1(t)dt=-ax1t t0,1dx2(t)dt=-ax2t t0,1dx3(t)dt=-ax3t-0.42ki x3 t0,23 dx3(t)d

22、t=-ax3t t23,1dx4(t)dt=-ax4t-ki x4 t0,23dx4(t)dt=-ax4t t23,1 由（1）式两边对t积分x0xdx1(t)x1(t)=01-adt,得x1(1)=x0(1)e-a.同理对（2）式求积分，结果得x1(2)=x0(2)e-a.（3）式对t从0到积分23x0xdx3(t)x3(t)=023（-0.42k-a)dt,得x3(23)=x3(0)e-20.42k+at3对t0,23积分，得到x3t=x0(3)e-0.42k+at（4）式对t从23到1积分，得到x13=x0(3)e-(0.28k+a)同上分别对（5）、（6）式积分，得x4(23)=x0(

23、4)e-2k+at3 x4t=x0(4)e-k+at x14=x0(4)e-23(k+a)由模型假设，鱼产卵数 n=1.109*105*12x323+x423=1.109*105*x03e-20.42k+at3+x0(4)e-23(k+a)因为第二年年初1龄鱼的数量n*1.22*10111.22*1011+n.因为x0(3)=x0(1)e-2a,所以得到 x0(3)e2a=1.109*105*12x03e-0.28k+23a+x0(4）e-2(k+a)3)问题一代码（第一种情况）：%最优捕鱼策略ch431%文件名：ch431.mx=sym('x');E3=0.42*x;%捕捞率

24、d=1.22*1011;r=0.8;q=d*exp(-(3*r+2/3*E3)*(32529.55*exp(r)+65059.1*exp(-2/3*x)/(1-exp(-(r+2/3*x);N10=d*q/(d+q);N40=d*q/(d+q)*exp(-(3*r+2/3*E3)/(1-exp(-(r+2/3*x);a3=E3/(r+E3)*(1-exp(-2/3*(r+E3)*d*q/(d+q)*exp(-2*r);a4=x/(r+x)*(1-exp(-2/3*(r+x)*d*q/(d+q)*exp(-(3*r+2/3*E3)/(1-exp(-(r+2/3*x);M3=17.86;M4=22

25、.99;M=M3*a3+M4*a4;M1=-M;M10=char(M);M11=char(M1);fplot(M10,0,100)E4=fminbnd(M11,0,100);E3=0.42*E4;d=1.22*1011;r=0.8;q=d*exp(-(3*r+2/3*E3)*(32529.55*exp(r)+65059.1*exp(-2/3*E4)/(1-exp(-(r+2/3*E4);N10=d*q/(d+q);N40=d*q/(d+q)*exp(-(3*r+2/3*E3)/(1-exp(-(r+2/3*E4);a3=E3/(r+E3)*(1-exp(-2/3*(r+E3)*d*q/(d+q

26、)*exp(-2*r);a4=E4/(r+E4)*(1-exp(-2/3*(r+E4)*d*q/(d+q)*exp(-(3*r+2/3*E3)/(1-exp(-(r+2/3*E4);M3=17.86;M4=22.99;M=M3*a3+M4*a4;E3E4Max=M%第一种情况，四龄鱼过一年全挂了N1=N10%各年龄组鱼的数量N2=N1*exp(-r)N3=N2*exp(-r)N4=N40问题二代码（情况一、情况二）：function fishf()global x0x=zeros(6,4);r1=0;r2=100;r,v1=fminbnd(fish1,r1,r2);%解第一问ee=exp(-0

27、.8);x10=f1(r);x0(1)=x10;x0(2)=x10*ee;x0(3)=x10*ee*ee;x0(4)=x10/(1-exp(-0.8-2*r/3)*exp(-2.4-0.84*r/3);rw1=-v1x0pause;r2,v2=fmincon(fish2,10,0,100,fish2con);%解第二问，5年的r不变时r2w2=-v2x=f2(r2);x_x0=x(6,:)./x0pause;r0=10 10 10 10 10;vlb=0 0 0 0 0;vub=150 150 150 150 150;r3,v3=fmincon(fish3,r0,vlb,vub,fish3co

28、n);%解第二问，5年r不同时r3w3=-v3x=f3(r3);x_x0=x(6,:)./x0function f=fish1(r)c1=0.8/(3*1.109e5*(exp(-2.4-2.84*r/3)/(1-exp(-0.8-2*r/3);c2=0.5*exp(-1.6-0.84*r/3)+(exp(-2.4-2.84*r/3)/(1-exp(-0.8-2*r/3);c3=17.86*0.42/(0.42*r+0.8)*exp(-1.6)*(1-exp(-0.8-2*r/3)*(1-exp(-1.6-2*r/3);c4=22.99/(r+0.8)*(exp(-2.4-0.84*r/3)/

29、(1-exp(-0.8-2*r/3)*(1-exp(-1.6-2*r/3);f=-1.22e9*r*(1-c1/c2)*(c3+c4);returnfunction f=f1(r)c1=0.8/(3*1.109e5*(exp(-2.4-2.84*r/3)/(1-exp(-0.8-2*r/3);c2=0.5*exp(-1.6-0.84*r/3)+(exp(-2.4-2.84*r/3)/(1-exp(-0.8-2*r/3);f=1.22e11*(1-c1/c2);returnfunction x=f2(r)global x0ee=exp(-0.8);x(1,:)=122 29.7 10.1 3.2

30、9*1e9;for i=1:5,n=3*1.109e5/0.8*(exp(-1.6/3)-exp(-0.8);n3=n*x(i,3)*exp(-0.84*r/3);n4=n*x(i,4)*exp(-2*r/3);x(i+1,1)=1.22e11*(0.5*n3+n4)/(1.22e11+0.5*n3+n4);x(i+1,2)=x(i,1)*ee;x(i+1,3)=x(i,2)*ee;x(i+1,4)=x(i,4)*exp(-0.8-2*r/3)+x(i,3)*exp(-0.8-0.84*r/3);endreturnfunction g,eg=fish2con(r)global x0ee=exp

31、(-0.8);x(1,:)=122 29.7 10.1 3.29*1e9;for i=1:5,n=3*1.109e5/0.8*(exp(-1.6/3)-exp(-0.8);n3=n*x(i,3)*exp(-0.84*r/3);n4=n*x(i,4)*exp(-2*r/3);x(i+1,1)=1.22e11*(0.5*n3+n4)/(1.22e11+0.5*n3+n4);x(i+1,2)=x(i,1)*ee;x(i+1,3)=x(i,2)*ee;x(i+1,4)=x(i,4)*exp(-0.8-2*r/3)+x(i,3)*exp(-0.8-0.84*r/3);enda=0.05g=(1-a)*x

32、0-x(6,:)-(1+a)*x0+x(6,:);eg=;returnfunction g,eg=fish3con(r)global x0ee=exp(-0.8);x(1,:)=122 29.7 10.1 3.29.*1e9;for i=1:5,n=3*1.109e5/0.8*(exp(-1.6/3)-exp(-0.8);n3=n*x(i,3)*exp(-0.84*r(i)/3);n4=n*x(i,4)*exp(-2*r(i)/3);x(i+1,1)=1.22e11*(0.5*n3+n4)/(1.22e11+0.5*n3+n4);x(i+1,2)=x(i,1)*ee;x(i+1,3)=x(i,2)*ee;x(i+1,4)=x(i,4)*exp(-0.8-2*r(i)/3)+x(i,3)*exp(-0.8-0.84*r(i)/3);enda=0.05g=(1-a)*x0-x(6,:)-(1+a)*x0+x(6,:);eg=;returnfunction x=f3(r)global x0ee=exp(-0.8);x(1,:)=122 29.7 10.1 3.29.*1e9;for i=1:5,n=3*1.109e5/0.8*(exp(-1.6/3)-exp(-0.8);n3=n*x(i,3)*exp(-0.84*r(i)/3);n4=n*x(i,4)*exp(-2*r(i)