- 1、本文档共10页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
最优捕鱼策略问题资料
课 时 授 课 计 划
实验序号: 12
一、课 题 实验八 最优捕鱼策略问题
二、课 型: 上机实验
三、目的要求:
能表述最优捕鱼策略问题的分析过程;
能表述模型的建立方法;
会建立最优捕鱼策略问题的优化模型;
会利用MATLAB或MATHMATICA求解最优捕鱼策略问题。
四、重点、难点:
重点:会利用MATLAB或MATHMATICA求解最优捕鱼策略问题
难点:会建立最优捕鱼策略问题的优化模型;
五、教学方法及手段:实验
六、参考资料:
1、 刘来福,曾文艺,数学模型与数学建模,北京师范大学出版社,1997第一版,2002第二版谭永基,俞文(鱼此),数学模型,复旦大学出版社,1997王树禾,数学模型基础,中国科学技术大学出版社,1996 李尚志等,数学建模竞赛教程,江苏教育出版社,1996 李大潜,中国大学生数学建模竞赛,高等教育出版社,1998
八、授课记录:
授课日期 2008. 班 次 信息051、2
九、授课效果分析:
教学效果良好
十、教学进程(教学内容、教学环节及时间分配等)
复习
导入课题
假设鱼群总量的变化随时间是连续的,从而利用微分方的知识建立最优捕鱼策略问题的优化模型
教学内容
为了保护人类赖以生存的自然环境,可再生资源(如渔业,林业资源)的开发必须适度。一种合理、简化的策略是,在实现可持续收获的前提下,追求最大产量或最佳效益。
考虑对鳀鱼的最优捕捞策略:假设这种鱼分4个年龄组,称为1龄鱼,2龄鱼,3龄鱼,4龄鱼。各年龄组每条鱼的平均重量分别为5.07,11.55,17.86,22.99(克),各年龄组的自然死亡率为0.8(/年),这种鱼为季节性集中产卵繁殖,平均每条4龄鱼的产卵量为1.109ⅹ1011(个),3龄鱼的产卵量为这个数的一半,2龄和1龄鱼不产卵,产卵和孵化期为每年的最后4个月,卵孵化并成活为1龄鱼,成活率(1龄鱼条数与产卵量n之比)为1.22ⅹ1011/(1.22ⅹ1011 + n)。
渔业管理部门规定,每年只允许在产卵孵化期前的8个月进行捕捞作业。如果每年投入的捕捞能力(如渔船数,下网次数等)固定不变,这时单位时间捕捞量将与各年龄组鱼群条数成正比,比例系数不妨称捕捞强度系数。通常使用13mm网眼的拉网,这种网只能捕捞3龄鱼和4龄鱼,其两个捕捞强度系数之比为0.42∶1。渔业上称这种方式为固定努力量捕捞。
建立数学模型分析如何实现可持续捕捞(即每年开始捕捞时渔场中各年龄组鱼群条数不变),并且在此前提下得到最高的年收获量(捕捞总重量)。
实验步骤
1、问题假设
鱼群总量的增加虽然是离散的,但对大规模鱼群而言,我们可以假设鱼群总量的变化随时间是连续的。
查阅有关鳀鱼的资料发现,鳀鱼一般在每年8月开始产卵,从而可以假设鱼群每年在8月底瞬间产卵完毕,卵在12月底全部孵化完毕。
龄鱼到来年分别长一岁成为i + 1龄鱼,i = 1,2,3。
4龄鱼在年末留存的数量占全部数量的比例相对很小,可假设全部死
亡。
连续捕获使各年龄组的鱼群数量呈周期性变化,周期为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
t∈[0,1],x1(0)= n ×
t∈[0,1],x2(0)= x1(1)
t∈[0,2/3],x3(0)= x2(1)
s.t. t∈[2/3,1],x3(-)= x3(+)
t∈[0,2/3],x4(0)= x3(1)
t∈[2/3,1],x4(-)= x4(+)
4、模型的求解
用MATLAB软件编程解微分方程组,先求得一元函数total(k)的表达式,画出total(k)函数的图形,
程序:
[buyu.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*x
您可能关注的文档
最近下载
- 胰岛素抵抗和代谢综合征课件.pptx VIP
- 胰岛素抵抗和代谢综合征.ppt VIP
- GB/T 17747.1-2011_天然气压缩因子的计算 第1部分:导论和指南.pdf
- 《基础护理学》第7章 休息与活动(含答案).docx
- 城市中心区综合性公园使用现状调查研究————以成都市人民公园为例.docx
- 产品档案管理制度及流程.pdf
- 中华民族一家亲,同心共筑中国梦.pptx VIP
- “社工+志愿者”联动模式的思考及对策研究--以惠州市河背社区志愿者项目为例.docx
- 国家开放大学,地域文化,人文武隆形考一 (3).pptx VIP
- (黑龙江省)新课标高中信息技术会考试题 学科整合 试题及答案..doc VIP
文档评论(0)