matlab最优捕鱼策略.pptxVIP

matlab最优捕鱼策略会计学第1页/共21页 （1）建立数学模型分析如何实现可持续捕捞（即每年开始捕捞时渔场中各年龄组鱼群条数不变)，并且在此前提下得到最高年收获量（捕捞总重量）。（2）某渔业公司承包这种鱼的捕捞业务5年，合同要求5年后鱼群的生产能力不能受到太大破坏。 已知承包时各年龄组鱼群数量分别为：122，，，（×109条）。如果仍用固定努力量的捕捞方式，该公司采用怎样的策略才能使总收获量最高。第2页/共21页模型的假设（1）假设只考虑一种鱼的繁殖和捕捞，鱼群增长过程中不考虑鱼的迁入与迁出.（2）假设各年龄组的鱼在一年内的任何时间都会发生自然死亡，产卵可在后四个月内任何时间发生.（3）假设3、4龄鱼全部具有生殖能力，或者虽然雄鱼不产卵，但平均产卵量掩盖了这一差异.（4）假设各年龄组的鱼经过一年后，即进入高一级的年龄组，但4龄鱼经过一年后仍视为4龄鱼.（5）假设对鱼的捕捞用固定努力量捕捞方式，每年的捕捞强度系数保持不变，且捕捞只在前八个月进行.第3页/共21页符号说明Ni(t)—t 时刻i 龄鱼的数量 Ni0 (k) —第k 年初i 龄鱼的数量Ni1 (k)___第k年底i 龄鱼的数量 (i=1，2，3，4) r—鱼的自然死亡率c —4龄鱼的平均产卵量 (则c/2为3龄鱼的平均产卵量)Qk —k年度鱼产卵总量p —鱼卵的成活率Mi—第i 龄鱼的平均重量(i=1，2，3，4) Ei —第i 龄鱼的捕捞强度系数ai —对i 龄鱼的年捕捞量(i=3，4）W—年总收获量，即W=M3a3+M4a4WW — 5年的总收获量为，即第4页/共21页由已知条件，可得（E为捕捞努力量）模型的建立第一步 得出基本模型 给出第k年底i 龄鱼的数量Ni1(k)与第k年初i 龄鱼的数量Ni0(k)之间的递推关系 给出年度捕鱼量 给出第k+1年初i 龄鱼的数量Ni0(k+1)与第k年初i 龄鱼的数量Ni0(k+1)的递推关系第二步 得出最终模型 根据可持续捕捞的要求, 给出约束条件及其目标函数解得从而第5页/共21页第一步（时间以年为单位，考虑一年内各龄鱼数量的演化） 已知r为自然死亡率，其定义为单位时间内死亡的鱼的数量与鱼的总量之比。由于不捕捞1、2龄鱼，所以在[t，t+Δt]内，根据死亡率的定义，变形得(1)第6页/共21页对于3、4龄鱼由于捕捞在前8个月进行，因此在前8个月内，捕捞与死亡均影响鱼的变化，因而微分方程变形为(2)由(2)式解得从而对于3、4龄鱼由于后四个月无捕捞，只有自然死亡，所以在后四个月其数量演化的方程为(3)解得从而第7页/共21页 由于仅在前八个月捕捞,且仅捕捞3龄鱼和4领鱼，而且捕捞强度系数表示的是单位时间内捕捞量与各年龄组鱼群总量成正比的比例系数，所以对i 龄鱼的年捕捞量为从而一年内捕鱼总收获量为第8页/共21页 由于每年各龄鱼的演化规律相同，且捕捞模式相同，综上可得： 第k年底i 龄鱼的数量Ni1(k)对第k年初i 龄鱼的数量Ni0(k) 的递推关系(4) 第k年的年度捕鱼收获量(5)第9页/共21页由各龄鱼之间的年龄增长关系,并假定产卵在年底一次完成,利用关系式(4)得从而第k+1年初i 龄鱼的数量Ni0 (k+1)与第k年初i 龄鱼的数量Ni0 (k) 的递推关系为(6)(5)式是每年捕鱼的总收获量，式 (6)刻划了鱼群各年龄组每年的变化情况，它们一起构成了基本模型。第10页/共21页第二步 （1）为了实现可持续的最大捕捞（即每年开始捕捞时渔场中各年龄组鱼群条数不变），即要求的前提下获得最高年收获量。结合基本模型，即可得到年度产量最优模型：(7)其中约束条件第二个等号说明各组鱼群条数及产卵量均与k无关。第11页/共21页优化模型（7）中的约束条件与k无关，故可把k丢掉，并利用E34，把目标函数和约束条件同时化简得(7’)注意到四个约束条件中含五个变量,因此从约束方程组可用符号计算软件解出Ni0(i=1，2，3，4),它们都是E4的函数,从而目标函数就是E4的一元函数.问题最终归结为一元函数的极值问题.该模型也可完全通过数值迭代求解！即： E4从0开始，逐渐增加，逐个计算W，挑出使W最大的E4。第12页/共21页模型的求解供参考的MATLAB计算程序%最优捕鱼策略ch431%文件名：ch431.mx=sym(x); E3=0.42*x;d=1.22*10^11;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