Matlab大作业-几种应用领域.pdfVIP

Matlab大作业-几种应用领域

Matlab大作业

陈大Ben_

I.在养殖业,最优捕鱼策略应用：

最优捕鱼策略问题

一．问题

为了保护人类赖以生存的自然环境，可再生

资源（如渔业、林业资源）的开发必须适度。一

种合理、简化的策略是，在实现可持续收获的前

提下，追求最大产量或最佳效益。

考虑对某种鱼的最优捕捞策略：

假设这种鱼分4个年龄组，称1龄鱼，...，

1龄鱼。各年龄组每条鱼的平均重量分别为

5.07,11.55,17.86,22.99(克)，各年龄组鱼的自然死

亡率均0.8(1/年)，这种鱼为季节性集中产卵繁

殖，平均每条龄鱼的产卵量为1.109×1011

(个)，龄鱼的产卵量为这个数的一半，龄鱼和龄

鱼不产卵。产卵和孵化期为每年的最后4个月，

综上所述，原问题实质上是给出了各年龄组

鱼群之间数量的变化规律，并给出了它们的自然

死亡率及捕捞和产卵的时间分布，并固定3、4

龄鱼捕捞能力的比值，要求选择一定的捕捞能力

系数，使得各年龄组鱼的数量在各年开始的第一

天条数不变（第一问），5年后鱼群的生产能力

不会有太大的破坏（第二问），并在此条件下，

求到最大捕获量。

二．建模

1.这种鱼在一年内的任何时间都会发生自

然死亡，即死亡是一个连续的过程；

2.捕捞也是一个连续的过程，不是在某一时

刻突然发生；3.鱼群的死亡率已考虑种群的相

互竞争及环境等因素；

4.3、4龄鱼产卵集中在9月初期,到第二年

初完成孵化；

5.龄鱼到来年分别长一岁成i+1龄鱼，其中

上一年存活下来的4龄鱼仍是4龄鱼。

三．程序

[buyu.m]

functiony=buyu(x)

globala10a20a30a40totalk;

symska10;

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.42*k)*x31

,x31(0)=a30);

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,x4

1(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);

eq1=a10-nn*1.22*10^11/(1.22*10^11+

nn);

S=solve(eq1,a10);

a10=S(2);symst;

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;

k=x;

y=subs(-total);

[buyu1]

globala10a20a30a40total;

[k,mtotal]=fminbnd(buyu,16,18);

ezplot(total,0,25);

xlabel(