数学建模实验报告最优捕鱼策略

数学建模实验报告

最优捕鱼策略

一.实验目的：

1、了解与熟练掌握常系数线性差分方程的解法；

2、通过最优捕鱼策略建模案例，使用MATLAB软件认识与掌握差分方程模型在实际生活方面的重要作用。

二.实验内容：（最优捕鱼策略）

生态学表明，对可再生资源的开发策略应在事先可持续收获的前提下追求最大经济效益。考虑具有4个年龄鱼：1龄鱼，… ，4龄鱼的某种鱼。该鱼类在每年后4个月季节性集中产卵繁殖。而据规定，捕捞作业只允许在前每年投入的捕捞能力固定不变，单位时间捕捞量与个年龄鱼群条数的比例称为捕捞强度系数。使用只能捕捞3、4龄鱼的13mm网眼的拉网，其两个捕捞强度系数比为0.42:1.渔业上称这种方式为固定力量捕捞。

该鱼群本身有如下数据：

1． 各年龄组鱼的自然死亡率为0.8（1/年），5.07,11.55,17.86,22.99（单位：g）；

2． 1龄鱼和2龄鱼不产卵，产卵期间，平均每条4龄鱼产卵量为（个），3龄鱼为其一半；

3． 卵孵化的成活率为1.22ⅹ1011/（1.22ⅹ1011 + n）（n有如下问题需要解决：

1）分析如何实现可持续捕获（即每年开始捕捞时各年龄组鱼群不变）前提下得到最高收获量；

2）合同要求某渔业公司在5年合同期满后鱼群的生产能力不能受到太大的破坏，承包时各年龄组鱼群数量为122,29.7,10.1,3.29（ⅹ量的捕捞方式下，问该公司应采取怎样的捕捞策略，才能使总收获量最高。

三. 模型建立

假设a、鱼群总量的增加虽然是离散的，但对大规模鱼群而言，我们可以假设鱼群总量的变化随时间是连续的； b、龄鱼到来年分别长一岁成为wilyes11收集 博客(与学习无关)：http://www.77cn.com.cn/u/1810231802

8个月进行，平均质量分别为

1.109ⅹ105

； ，并在此9条），在固定努力i + 1龄鱼，i = 1，

其为产卵总量）10

数学建模实验报告

2，3； c、4龄鱼在年末留存的数量占全部数量的比例相对很小，可假设全部死 亡。d、连续捕获使各年龄组的鱼群数量呈周期性变化，周期为1年，可以只考虑鱼群数量在1年内的变化情况。（且可设xi（t）：在t时刻i龄鱼的条数，i = 1，2，3，4；n：每年的产卵量；k：4龄鱼捕捞强度系数；2ai0：每年初i龄鱼的数量，i = 1，2，3，4；）

进而可建立模型如下：

max（total（k））=17.86 0

2/3

0.42kx3(t)dt 22.99

2/3

kx4(t)dt

1.22 1011

11dx(t)

1 0.8x1(t) t∈[0，1]，x1（0）= n ×1.22 10 n

dtdx2(t)

0.8x2(t)dt t∈[0，1]，x2（0）= x1（1） dx3(t)

(0.8 0.42k)x3(t)

dt t∈[0，2/3]，x3（0）= x2（1） dx3(t)22

0.8x3(t)

s.t. dt t∈[2/3，1]，x3（3-）= x3（3+） dx4(t)

(0.8 k)x4(t)

dt t∈[0，2/3]，x4（0）= x3（1） dx4(t)22

0.8x4(t)

dt t∈[2/3，1]，x4（3-）= x4（3+） 22 n 1.109 105[0.5x3() x4( )]

33

四. 模型求解（含经调试后正确的源程序）

1. 先建立一个buyu.m的M文件： function y=buyu(x);

global a10 a20 a30 a40 total k; syms k a10;

x1=dsolve('Dx1=-0.8*x1','x1(0)=a10'); t=1;

a20=subs(x1);

x2=dsolve('Dx2=-0.8*x2','x2(0)=a20'); t=1;

a30=subs(x2);

x31=dsolve('Dx31=-(0.8+0.4*k)*x31','x31(0)=a30');

1

数学建模实验报告

t=2/3;

a31=subs(x31);

x32=dsolve('Dx32=-0.8*x32','x32(2/3)=a31'); t=1;

a40=subs(x32);

x41=dsolve('Dx41=-(0.8+k)*x41','x41(0)=a40'); t=2/3;

a41=subs(x41);

x42=dsolve('Dx42=-0.8*x42','x42(2/3)=a41'); t=2/3;

a31=subs(x31);

nn=1.109*10^5*(0.5*a31+a41);

Equ=a10-nn*1.22*10^11/(1.22*10^11+nn); S=solve(Equ,a10); a10=S(2,1); syms t; k=x;

t3=subs(subs(int(0.42*k*x31,t,0,2/3))); t4=subs(subs(int(k*x41,t,0,2/3))); total=17.86*t3+22.99*t4; y=subs((-1)*total)

2. 再建立一个buyu1.m的M文件： global a10 a20 a30 a40 total; [k,mtotal]=fminbnd('buyu',0,20); ezplot(total,0,25); xlabel(''); ylabel(''); title('');

format long; k

total=-mtotal; a10=eval(a10) a20=eval(a20) a30=eval(a30) a40=eval(a40) format short clear

果分析

． 鱼总量与时间图：

2

五． 结1

数学建模实验报告

2. 可以看出捕捞强度对收获量的影响：

实验输出数据：

y =

-3.6757e+011 y =

-3.9616e+011 y =

-4.0483e+011 y =

-4.0782e+011 y =

-4.0802e+011 y =

-4.0805e+011 y =

-4.0805e+011 y =

3

数学建模实验报告

-4.0805e+011 y =

-4.0805e+011 y =

-4.0805e+011 y =

-4.0805e+011 y = 0 y =

-4.0667e+011 k =

18.25976795085083 total =

4.080548655562244e+011 a10 =

1.195809275167686e+011 a20 =

5.373117428928620e+010 a30 =

2.414297288420686e+010 a40 =

8.330238542343275e+007

则k=18.25976795085083时，最高年收获量为（克），此时每年年初1，2，31.195809275167686×1011 5.373117428928620×1010 2.414297288420686×1010 8.330238542343275×107

4total=4.080548655562244

114

×10，年龄组鱼的数量分别为：

数学建模实验报告

六． 实验总结

本次实验的目的是了解差分方程（递推关系）的建立及求解，以及掌握用差分方程（递推关系）来求解现实问题的方法。实验中假设鱼群总量的变化随时间是连续的，从而利用微分方的知识建立最优捕鱼策略问题的优化模型。通过实验加深了对概念和方法的理解，了解了差分方程的程序解法。

学生签名：

年 月 日

七．教师评语及成绩

教师签名：

年 月 日

5