- 1、本文档共14页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
最优捕鱼策略(大).doc
最优捕鱼策略
为了保护人类赖以生存的自然环境,可再生资源(如渔业、林业资源)的开发必须适度。一种合理、简化的策略是,在实现可持续收获的前提下,追求最大产量或最佳效益。
考虑对某种鱼(鯷鱼)的最优捕捞策略:
假设这种鱼分4个年龄组,称1龄鱼,…,4龄鱼。各年龄组每条鱼的平均重量分别为5.07,11.55,17.86,22.99(克),各年龄组鱼自然死亡率均为0.8(1/年),这种鱼为季节性集中产卵繁殖,平均每条4龄鱼的产卵量为1.109(105(个),3龄鱼的产卵量为这个数的一半,2龄鱼和1龄鱼不产卵。产卵和孵化期为每年的最后3个月,卵孵化并成活为1龄鱼,成活率(1龄鱼条数与产卵总量n之比)为1.22(1011/(1.22(1011 ( n)。
渔业管理部门规定,每年只允许在产卵孵化期前的8个月内进行捕捞作业。如果每年投入的捕捞能力(如渔船数、下网次数等)固定不变,这时单位时间捕捞量将与各年龄组鱼群条数成正比,比例系数不妨称捕捞强度系数。通常使用13mm网眼的拉网,这种网只能捕捞3龄鱼和4龄鱼,其两个捕捞强度系数之比为0.42(1,渔业上称这种方式为固定努力量捕捞。
建立数学模型分析如何实现可持续捕获(即每年开始捕捞时渔场中各年龄组鱼群条数不变),并且在此前提下得到最高的年收获量(捕捞总重量)。
2)某渔业公司承包这种鱼的捕捞业务5年,合同要求5年后鱼群的生产能力不能受到太大破坏。已知承包时各年龄组的鱼群的数量分别为:122,29.7,10.1,32.9((109条),如果仍用固定努力量的捕捞方式,该公司应采取怎样的策略才能使总收获量最高.
基本假设:
1、? 鱼群生活在稳定的环境中,不考虑鱼群的迁入和迁出,也不考虑鱼群的空间分布;
2、? 1龄鱼、2龄鱼、3龄鱼、4龄鱼均可以在一年即一个周期的任意时间内死亡;
3、成活的i龄鱼(i=1,2,3)每经过一年即一个周期变为(i+1)龄鱼,而4龄鱼不变;
4、 假设相邻两个年龄组的鱼群在相邻两年之间的变化是连续的,即第T年底第i年龄组的鱼的条数等于第T+1年初第i+1年龄组的鱼的条数;
各年龄组鱼的平均重量和自然死亡率稳定,不考虑由于饲养技术、环境等因素引起变化;
6、 只考虑采用固定努力量捕捞方式下的捕捞策略。
符号说明:
符号 意义 说明 单位 r 对4龄鱼的捕捞强度 变量 1/年 t 时间 变量 年 t时刻i龄鱼的数量(i=1,2,3,4) 变量 条 n3 3龄鱼在第四季度的的产卵量 未知量 个 n4 4龄鱼在第四季度的的产卵量 未知量 个 未知量
模型建立:
一、各龄鱼的变化规律
解之得:
一年内的捕捞量f(r)为(单位:g):
对于平衡捕捞,有:
由可得:
由
将
对承包五年的情况,分两种情形处理。
第一种:假定五年中捕捞强度r?一直保持不变,且第6年初的各龄鱼数与平衡捕捞时的年初值的ε倍(如ε=0.95即95%):
第二种:假定五年中各年的捕捞强度不同,设ri为第i年的捕捞强度,且第6年初的各龄鱼数与平衡捕捞时的年初值的ε倍(如ε=0.95即95%):
附: MATLAB程序
%最优捕鱼策略.
function fishf()
global x0
x=zeros(6,4);
r1=0;r2=100;
[r,v1]=fminbnd(@fish1,r1,r2);%解第1问.
ee=exp(-0.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);
r
w1=-v1
x0
pause;
[r2,v2]=fmincon(@fish2,10,[],[],[],[],0,100,@fish2con);%解第2问,五年的r不变时.
r2
w2=-v2
x=f2(r2);
x_x0=x(6,:)./x0
pause;
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,@fish3con);;%解第2问,五年的r不相同时.
r3
w3=-v3
x=f3(r3);
x_x0=x(6,:)./x0
function f=fish1(r)
c1=0.8/(3*1.109e5*(exp(-1.6/3)-exp(-0.8)));
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
文档评论(0)