数学建模实验报告最优捕鱼策略
最优捕鱼策略问题
最优捕鱼策略问题摘要本文以最优捕鱼策略为主题，在logistic模型基础上建立了可持续发展捕鱼策略模型,并借助计算机Matlab，运用二分法近似求得了模型最优解。
在此基础上提出了灵敏度函数S，并由此判断死亡率w和捕捞强度E的变化对产量变化的影响。
最后根据实际生产需求，分析死亡率w对最大产量Qm的影响。
对于问题1，我们首先考虑不存在捕捞情况下的模型，再加入捕捞强度分析，最后根据问题1的条件（每年开始捕捞时渔场中各种年龄组鱼群条数不变）建立方程组，得到可持续发展捕鱼策略模型，解得方程组后在w=0.8时绘图得到最大产量Qm=3.8871*10^11。
对于问题2，我们引用了灵敏度函数S(ω,Q),起意义为ω变化率与Q变化率的比值，例如S=0.1，即表示当死亡率变化1%的时候，产量Q变化0.1%。
发现在问题1取得最大产量的情况下，死亡率每增加1%，最大产量减少1.743%。
并给出了不同死亡率w和产量下S的函数。
对于问题3，方法与问题2相似，灵敏度函数S(E,Q)在问题1的情况下，捕捞强度系数E每增加1%，产量Q减少0.0010%。
并给出了不同捕捞强度E和产量Q下S的函数。
对于问题4,我们取不同的死亡率w，得到不同的最大产量Q,利用MATLAB用函数拟合的方法得到了相似度很高的4阶拟合函数Qm（w）仿照问题2求解了灵敏度函数S(E,Qm)，发现了在问题1求得最大产量的时候，死亡率的波动对最大产量的影响是相对较大的。
现实生产中可表现为一段时间内大量鱼群的死亡对渔民的收获量会造成比较大的损失。
为此我们找到了影响较小的点，当把死亡率控制在0.957附近时，鱼群的突然大数目死亡短时间内对渔民造成的损失最小。
对此我们提出了一些策略。
关键词：可持续发展捕鱼策略模型，灵敏度分析，函数拟合，微分方程。
一、问题重述以鳀鱼为例，制定一种最优的捕鱼策略，要求实现可持续捕捞，并且在此前提下得到最高的年收获量，并进一步考虑自然死亡率和捕捞强度系数，提出相关建议。
最优捕鱼策略KYM)
存活率 pi ~同一时段的 xi+1与 xi之比
x* 1, p0 , p0 p1,L p0 p1L p99 T
（与pi 的定义 xi1(k 1) pi xi (k) 比较）
Leslie模型的应用：公园大象管理
南非的一家大型自然公园放养了大约11000头大象，管理部门希望 为大象创造一个健康的生存环境，将大象的总数控制在11000头左右。每 年，公园的管理人员都要统计当年大象的总数。过去20年里，公园每年都 要处理一些大象，以便保持大象总数维持在11000头左右，通常都是采用 捕杀或者迁移的方法来实现。统计表明，每年约处理600-800头大象。
02
,
p0 p1 p2
03
,L
,
p0 p1 p2 L
100
p99 )
1是A的一个特征值，相应的特征向量为
a0 a1 a2 L
1
0
0L
a99 a100
0
0
er=1,1,L ,1T
A
M 1LM
0
0 0
1 0
0L 1L
0 0
M M M O O
0
0
M
与特征值0相应的特r征向量为 X0=M e
(ii)
lim
t
X
(t
பைடு நூலகம்
)
/
0t
=
cX0
,
其中c是与X(0)有关的常数,即当t充分大时，
X (t) c0t X 0.
定理1的证明：
由非负矩阵的谱性质，L矩阵有最大的特征值是单重的且为正数。
设为 0 0, 则下面只需求出一个相应的特征向量 。
最优捕鱼策略 建模论文
最优捕鱼策略摘 要为了保护人类赖以生存的自然环境，实现资源的可持续发展。
可再生资源（如渔业、林业资源等）的开发必须适度。
因此本文针对可持续捕鱼提出的两个问题建立了两个优化模型。
模型1针对问题1，已知各年龄组鱼群之间的数量变化规律、自然死亡率、3,4龄鱼的捕捞强度系数之比、捕捞和产卵时间范围，要满足在实现可持续捕捞的前提下得到最高的年收获量。
即保证每年开始捕捞时渔场中各年龄组鱼群条数不变的情况下得到最高年收获量。
以3，4龄鱼的年产量为目标函数，各龄鱼在年初和年末的条数为约束条件，建立规划模型，利用和数学软件进行求解得到最高收获量为。
lingo matlab 113.88707610×模型2针对问题2，根据题意渔业公司承包这种鱼的捕捞，并且要求5年后鱼群的生产能力不能受到太大破坏，但又要使收获量最大。
首先，题中已给出各年龄组鱼群的初始值，我们利用模型1求出第6年初各年1龄鱼的数量；其次，根据问题1中的捕捞量表达式，可写出5年的捕捞总量表达式，以5年捕捞总量最大为前提，利用matlab 软件求解出此时的捕捞强度系数；再次，计算出第一年初与第六年1龄鱼的数量之比为，得到在此捕捞强度下不会使5年后鱼群的生产能力有太大的破坏；0.05423%最后，当捕捞强度系数为()17.5,17.8k ∈时，得到鱼的最高收获量为。
121.605610×关键词：自然死亡率 捕捞强度系数 ling 和数学软件 最高收获量 优化o matlab1 问题重述1.1 问题背景为保护人类赖以生存的自然环境，可再生资源（如渔业，林业等资源）的开发必须适度。
一种合理、简化的策略是，在实现可持续收获的前提下追求最大产量或最佳效益。
1.2 基本条件1 鱼分为4个年龄组：1龄鱼，2龄鱼，3龄鱼，4龄鱼。
各年龄组每条鱼的平均重量分别为5.07，11.55，17.86，22.99（克）；2 每年龄组鱼的自然死亡率均为0.8（1∕年）；3 这种鱼为季节性集中产卵繁殖，平均每条4龄鱼的产卵量为1.109×（个），3龄鱼的产卵量为51012×1.109×(个)，2龄鱼和1龄鱼不产卵，产卵和孵化期为每年的最后四个月；5104 卵孵化并成活为1龄鱼，成活率为（1龄鱼条数与产卵总量n 之比）；11111.2210/(1.2210)n ××+5 每年只允许在产卵孵化期的前8个月进行捕捞；6 每年投入的捕捞能力(如渔船数，下网次数等)固定不变，单位时间捕捞量与各年龄组鱼群条数成正比，比例系数为捕捞强度系数；7 使用13mm 网眼的拉网，以固定努力量的方式捕捞。
数学建模—最佳捕鱼方案
三、 符号说明
;当k 1 x ：表示 i 龄鱼第 j 年的年初（或年末）的鱼量（ k 0或1， 当k 0时， 表示年初 时表示年末。 i 1,2,3,4; j 1,2, ） 条 ； r ：表示各年龄组鱼群的死亡率： 0.8(1 年) ； ：表示 4 龄鱼的捕捞强度系数，则 3 龄鱼的捕捞强度系数为 0.42 ； n ：产卵总量 个 ； Z：捕鱼总重量 g ； xij t ：表示第 j 年 t 时刻 i 龄鱼的数量 条 ； j ：表示第 j 年的捕鱼总量；
4
年 收 获 总 量 （ g）
4.2 4.15 4.1 4.05 4 3.95 3.9 3.85
x 10
11
3.8 0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
死 亡 率 （ 1/ 年 ）
由上图可直观地看出：死亡率与年收获总量成正比例关系，即当死亡率增加时， 年收获总量则减少；反之，增加。由此可知，死亡率对年收获总量有显著的影响。 2．对模型中捕捞强度系数 的灵敏度分析 模型中其它因素不变， 只考虑 从 10 变到 19 时最大的年收获总量的变化情况， 分析 的变化对模型的影响（见下图）
年 收 获 总 量 （ g）
3.95 x 10
11
3.9
3.85
3.8
3.75
3.7
3.65
3.6
3.55
3.5 10
11பைடு நூலகம்
12
13
14
15
16
17
18
19
4龄 鱼 的 捕 捞 强 度 系 数
由上图可直观地看出：捕捞强度系数也是影响年收获总量的重要因素，年收获总量 随捕捞强度系数的增加而增加。只是增长速率逐渐减慢。 七、 模型评价与推广 模型的评价： 优点：1. 本文建立的模型与实际相联系，考虑到一些实际情况，从而使模型较贴近实 际；通用性.，推广性较强。 2．模型方便、直观，可以实现计算机模拟。 缺点： 1．模型虽然考虑到了很多因素，但为了建立模型，忽略了一些影响因素，具有 一定的局限性。 2．在建模过程中，简化了一些因素，得到了最优方案可能与实际有一定的出入。 模型的推广： 模型建立思想不但适合捕鱼方面，而且适合其它相关方面，只需稍加改动即可。
数学建模课程设计_最佳捕鱼方案
数学建模论文姓名: 文勇学号：201315020220论文标题：最佳捕鱼方案1.问题的提出一个水库，由个人承包，为了提高经济效益，保证优质鱼类有良好的生活环境，必须对水库的杂鱼做一次彻底清理，因此放水清库。
水库现有水位平均为15米，自然放水每天水位降低0.5米，经与当地协商，水库水位最低降至5米，这样预计需要二十天时间，水位可达到目标。
据估计水库内尚有草鱼25000余公斤，鲜活草鱼在当地市场上，若日供应量在500公斤以下，其价格为30元/公斤；日供应量在500—1000公斤，其价格降至25元/公斤，日供应量超过1000公斤时，价格降至20元/公斤以下，日供应量到1500公斤，已处于饱和，捕捞草鱼的成本水位于15米时，每公斤6元；当水位降至5米时，为3元/公斤。
同时随着水位的下降草鱼死亡和捕捞造成损失增加，至最低水位5米时损失率为10%。
承包人提出了这样一个问题：如何捕捞鲜活草鱼投放市场，效益最佳？2.问题分析通过简单的分析和思考，该问题可以归为一个数学规划问题。
条件（1）（2）是针对目前状况的约束，条件（3）是通过卖鱼可以获得的利润，条件（4）是对成本的约束。
在四个条件约束的情况下，我们可以建立模型。
由于对损失率的理解不同，我们进行了不同的假设，并在这些假设下建立了模型一和模型二、三。
模型一中，损失率是基于水库草鱼的总量，草鱼的损失是一些定值的累加。
而在模型二、三中，为了更接近现实生活中的情况及人们的认知观，我们对第n天草鱼的损失率的理解是基于第n-1天剩下的草鱼而言。
模型二将不考虑日供应量超过1500kg的情况，而模型三考虑。
模型三的建立采用多目标的规划方法进行求解。
3.条件假设1、日供应量不受外界条件的变化而变化，是一定的。
2、当天售出的草鱼数量等于当天捕捞的草鱼。
3、水位的变化除了每天的自然放水，不考虑蒸发等其他的情况。
4、假设在放水清库的过程中，随着水位的下降,捕捞成本成呈递减等差数列，而草鱼的损失成递增等差数列。
全国数模竞赛题选讲1-最优捕鱼策略(96A)
Jingsaitixuanjiang
Jingsaitixuanjiang
假设这种鱼分 4 个年龄组，称 1 龄鱼，…，4 龄
0.8 3
k
2 3
s 40 ) e
2 3
0.8
1.22 10
k
2 3
11
10
11
a (1 e
)( e
0.42 k
s 30 2 e
s 40 ) e
0.8
]
F 3 s 30 F 4 s 40
F3
1 .6 [ 1 .2 2 1 0
3 a (1
重复 ⑵ ，根据递推关系算出下一年的
s 12 , s 22 , s 32 , s 42 ;
⑷ 再重复 ⑵、⑶ 当计算到年初与年末的各龄鱼 的数量一致时，即鱼群稳定为止，根据
G P 3 m 3 P4 m 4 算出年捕获量；
⑸ 另定 k 值，重复 ⑴ ~ ⑷； ⑹ 根据年捕获量最大原则，最后确定最佳的 k 值。
s 21 s 1 s 10 e
3）3 龄鱼即上一年末 2 龄鱼
0.8
s 31 s 2 s 20 e
4）4 龄鱼即上一年末 3 龄鱼
0.8
s 41 s 3 s 30 e
0.8
e
0.24 k
2 3
Jingsaitixuanjiang
3 、鱼群持续变化的递推关系
捞方式，该公司应采取怎样的策略才能使总收获量
最优捕鱼策略(1)
第二步 得出最终模型 • 根据可持续捕捞的要求, 给出约束条件及其目标函数
最优捕鱼策略(1)
由于每年各龄鱼的演化规律相同，且捕捞模式相
同，综上可得：
第k年底i 龄鱼的数量Ni1(k)对第k年初i 龄鱼的数量Ni0(k) 的
递推关系
(4最优捕鱼策略(1)
由各龄鱼之间的年龄增长关系,并假定产卵在年底一次完成,利用关系 式(4)得
从而第k+1年初i 龄鱼的数量Ni0 (k+1)与第k年初i 龄鱼的数量Ni0 (k) 的递
最优捕鱼策略(1)
3rew
演讲完毕，谢谢听讲!
再见，see you again
2020/11/17
最优捕鱼策略(1)
最优捕鱼策略(1)
2020/11/17
最优捕鱼策略(1)
（1）建立数学模型分析如何实现可持续捕捞（即每年开始捕捞时渔场中
各年龄组鱼群条数不变)，并且在此前提下得到最高年收获量（捕捞总重 量）。 （2）某渔业公司承包这种鱼的捕捞业务5年，合同要求5年后鱼群的生产 能力不能受到太大破坏。
已 知 承 包 时 各 年 龄 组 鱼 群 数 量 分 别 为 ： 122 ， 29.7 ， 10.1 ， 3.29 （×109条）。如果仍用固定努力量的捕捞方式，该公司采用怎样的策略才 能使总收获量最高。
Qk —k年度鱼产卵总量
p —鱼卵的成活率
Mi—第i 龄鱼的平均重量(i=1，2，3，4) Ei —第i 龄鱼的捕捞强度系数 ai —对i 龄鱼的年捕捞量(i=3，4） W—年总收获量，即W=M3a3+M4a4 WW — 5年的总收获量为，即
数学建模——最优捕鱼模型
最优捕鱼模型一.问题的重述捕鱼业在当今社会中十分重要的行业，捕鱼量的大小决定着捕鱼的经济效益，其中捕鱼量与捕鱼时间有着密切关联. 所以如何利用数学模型了解捕鱼量与捕鱼时间之间的关系，是一个具有现实意义的问题.现假设在一个鱼塘中投放若干鱼苗，鱼苗尾数随着时间的增长而减少，且相对减少率为常数；每尾鱼的重量随着时间增长而增加，且由于喂养引起的每尾鱼重量增加率与鱼的表面积成正比，由于消耗引起的减少率与其重量本身成正比. 分析如下问题：问题一：建立尾数和时间的微分方程并求解；问题二：建立每尾鱼重量和时间的微分方程并求解；问题三：用控制网眼的方法不捕小鱼，从一定时刻开始捕捞，用尾数的相对减少率表示捕捞能力，分析开始捕鱼的最佳时刻，使得捕获量最大，并建立相关模型.二.问题分析1.针对问题一，根据相对减少率的数学定义，可以建立鱼尾数和时间的微分方程；2.针对问题二，将鱼体假设为球体，得出鱼的表面积与它重量的关系，使得鱼的重量完全成为一个关于时间的函数，进一步建立出鱼重量与时间的微分方程；3.针对问题三，将捕捞行为看作连续的过程，瞬时捕捞量与瞬时捕鱼尾数、每尾鱼瞬时重量呈正相关关系，瞬时捕鱼尾数与捕捞能力有关，每尾鱼瞬时重量可由对问题二的解答得出，总捕捞量即为瞬时捕捞量关于时间的积分.三.基本假设1．假设自然因素不会对鱼的尾数产生影响；2．假设在整个捕捞过程中鱼没有繁衍行为；3．假设每尾鱼都均衡生长；4．假设在捕捞过程中鱼的条数连续；5．假设鱼为球体.四．符号表示五．模型建立与求解模型一. 鱼苗尾数的相对减少率为常数r . 由相对减少率的定义得()()()t t t t n n rn t +∆-=-∆ 即()()()00lim lim t t t t t t n n rn t +∆∆→∆→-=-∆ 即()t dn rn dt=- 解得0rt n n e -=模型二. 假设鱼为球体，体积为V ，表面积为S ，半径为R ，重量为G ，初始重量为0G ，鱼的密度为ρ;且每尾鱼的重量随着时间增长而增加，其中由于喂养引起的每尾鱼重量增加率与鱼表面积成正比（比例系数为1k ），由于消耗引起的减少率与其重量本身成正比（比例系数为2k ）. 由343V R π=，2=4S R π，G V ρ=得2233S G ρ⎛⎫= ⎝⎭令23=b ρ⎛⎫ ⎝⎭又由于12=-dG k S k G dt,=0t ，0G G =所以231-11322+k t k b k b G e k k ⎡⎤⎫=⎢⎥⎪⎭⎣⎦模型三. 控制网眼不捕小鱼，鱼塘中瞬时鱼尾数用(t)n 表示，捕捞能力（E ）可以用尾数的相对减少率1dn n dt表示，从T 时刻开始捕捞，使得捕捞量W 能够最大.其中减少量包括自然减少量（即第一模型中的减少量）和捕捞量.此时，-(t)0(t)=-at n n e En-0-0(e )11=-=-=a e at at d n dn E n dt n dt所以，--00(t)==1+(1+)at aT T Tan e an W En dt dt e a a a ∞∞=⎰⎰ 则，在此模型下，捕捞时间越早，捕捞量越大.模型四. 建立在模型三的基础上，捕捞量的大小不仅取决于鱼尾数(t)n ，还取决于鱼的重量G .即(t)TW En Gdt ∞=⎰所以，231--0113(t)22=+1+at k t T T an e k b k b W En Gdt e dt a k k ∞∞⎡⎤⎫=⎢⎥⎪⎭⎣⎦⎰⎰ 可根据此函数求得最大捕捞量所对应的时刻T .感谢下载!欢迎您的下载，资料仅供参考。
数学建模案例——最佳捕鱼方案
最佳捕鱼方案摘要：本文解决的是一个最佳捕鱼方案设计的单目标线性规划问题，目的是制定每天的捕鱼策略，使得总收益最大。
根据题设条件，结合实际情况，我们设计了成本与损失率随天数的增加成反比变化的函数曲线（见图三所示），并导出总收益的表达式： 212121111i i i i i i i i W w p s q m =====⨯-⨯∑∑∑。
由于价格是关于供应量的分段函数（见图一所示），我们引入“0－1”变量法编写程序（程序见附录一），并用数学软件LINGO 求解，得到最大收益(W)为441291.4元，分21天捕捞完毕。
其中第1～16天，日捕捞量在1030～1070公斤之间，第17～21天的日捕捞量为1610～1670公斤之间（具体数值见正文）。
由结果分析，我们对模型提出了优化方向，例如人工放水来降低成本。
关键词：“0-1”整数规划，单目标线性规划，离散型分布。
一． 问题重述一个水库，由个人承包，为了提高经济效益，保证优质鱼类有良好的生活环境，必须对水库里的杂鱼做一次彻底清理，因此放水清库。
水库现有水位平均为15米，自然放水每天水位降低0.5米，经与当地协商水库水位最低降至5米，这样预计需要二十天时间，水位可达到目标。
据估计水库内尚有草鱼二万五千余公斤，鲜活草鱼在当地市场上，若日供应量在500公斤以下，其价格为30元/公斤；日供应量在500—1000公斤，其价格降至25元/公斤，日供应量超过1000公斤时，价格降至20元/公斤以下，日供应量到1500公斤处于饱和。
捕捞草鱼的成本水位于15米时，每公斤6元；当水位降至5米时，为3元/公斤。
同时随着水位的下降草鱼死亡和捕捞造成损失增加，至最低水位5米时损失率为10%。
承包人提出了这样一个问题：如何捕捞鲜活草鱼投放市场，效益最佳？二． 模型假设1．池塘中草鱼的生长处于稳定状态，不考虑种群繁殖以及其体重增减，即在捕捞过程中草鱼总量保持在25，000公斤不变。
2．第一天捕捞时水位为15m ，每天都在当天的初始水位捕捞草鱼，水库水位每天按自然放水0.5m 逐渐降低，20天后刚好达到最低要求水位5m 。
最优捕鱼策略_数学建模
精心整理西安邮电大学（理学院）数学建模报告摘要为了保护人类赖以生存的自然环境，可再生资源（如渔业、林业资源）的开发必须适度。
本文实际上就是为了解决渔业上最优捕鱼策略问题，即在可持续捕捞的前提下，追求捕捞量的最大化。
问题一采用条件极值列方程组的方法求解，即1龄鱼的数量由3龄鱼和4龄鱼的产卵孵化而来；2，3龄鱼的数量分别由上一年1龄鱼，2龄鱼生长而来；4龄鱼由上一年的3龄鱼和上一年末存活的4龄鱼组成。
最后得到：只要每年1-8月份3、4龄鱼捕捞总量小于、，就可以实现总捕捞量最大为；对结果分析得到捕捞的对象主要是3龄鱼，当3龄与4龄鱼的捕捞系数发生变化时，总的捕捞量变化不大。
???问题二给出年初各龄鱼的数量，要求在5年后鱼群的生产能力没有受到太大条），如果仍用固定努力量的捕捞方式，该公司采取怎样的策略才能使总收获量最高。
二、模型假设1、这种鱼分为四个年龄组：1龄鱼，2龄鱼，3龄鱼，4龄鱼；2、各年龄组每条鱼的平均重量分别为5.07克,11.55克,17.86克,22.99克；3、各年龄组鱼的自然死亡率均为0.8（1/年）；m……i龄鱼每条鱼的平均重量in……9月底该种鱼总共产卵数量*n……卵孵化成幼鱼进入1龄鱼阶段的数量k……对i龄鱼活鱼的捕捞强度系数i四、问题分析针对问题一：如何在满足可持续捕捞的前提下，实现每一年捕鱼的最大量（重量），文中给出各龄鱼在年底转化的具体情况：1龄鱼数量由3龄鱼和4龄鱼的产卵孵化而来；2，3龄鱼的数量分别由上一年龄段的鱼经自然死亡以及捕捞生长而来；4龄鱼是由上一年段3龄鱼经自然死亡以及捕捞后生长的和原有的4龄鱼组成的，并且规定只在每年的前八个月出船捕捞。
那么根据以上信息我们可以建立动态整型规划模型，即以每年的前八个月作为动态规划中的8种状态，在满足文中的可持续捕捞的约束条件下，先确定这前八个月中，每个月的捕捞量，最后求得这八个月总捕捞量的最大值；当然我们还可以建立微分方程模型，把每一龄鱼的数量变化看成是随时间连续变化的，将每一龄鱼的初始数量减去第八个月末的数量⎪⎩⎪⎨≤≤-=---129,1,1,1,,j c x x i j i j i i i j i j i 这个等式说明了该模型中我们把每一个月看做一个时间单位，鱼的数量随时间的变化是离散的，当每个月月初各龄鱼的数量固定时，该月要捕捞的总的活鱼数量也就固定了。
最优捕鱼策略试验报告
最优捕鱼策略实验报告学号:104080298 姓名:宁亚会班级:10D摘要为了保护人类赖以生存的自然环境，可再生资源(如渔业、林业资源等)的开发必须适度。
而在社会经济生活中,我们要使商业活动在一段时期内达到最大收益,因此我们要合理的开发资源，这时,我们不仅要考虑商业活动的当前经济效益,还要考虑生态效益及由此产生的对整体经济效益的影响。
本文就是对渔业这类可再生资源的开发问题进行研究，利用相关的数学软件进行求解。
对于问题一，我们考虑渔场生产过程中的各年龄组鱼群数量的制约因素，将其分为两大类，第1,2龄鱼群为一类，该鱼群数量变化在一年内只受自然死亡率制约，写出鱼群数量满足的微分方程;第3,4龄鱼群为一类，其数量变化在前8个月受捕捞强度和自然死亡率影响，后4个月只受自然死亡率的制约，分阶段写出写出鱼群数量满足的微分方程;根据微分方程，求出在某时刻各鱼群的数量表达式(类似于人口增长模型)。
因为捕捞是连续的，所以任意一个时刻的捕捞量为捕捞强度乘以鱼群的数量，又捕捞只在前8个月进行，则年捕捞量为前8个月各时刻鱼群数量的积分。
最后建立年总捕捞量的函数与生产过程中满足的关系式，转化为非线性规划模型，利用lingo和matlab软件分别求解。
对于问题二，题中已给出各年龄组鱼群的初始值，我们利用问题一中所得到的迭代方程，可迭代地求出第i年初各年龄组鱼群的数量;再根据问题一中的捕捞量表达式，可写出5年的捕捞总量表达式，以5年捕捞总量最大为前提，利用matlab软件求解出此时的捕捞强度，然后再验证在此捕捞强度下会不会使5年后鱼群的生产能力有太大的破坏。
最后，我们得出以下结论:可持续捕获条件下，捕捞强度为17.36292时，达到最大捕113.88707610，捞总质量g; 5年后鱼群的生产能力不会有太大的破坏条件下，捕捞强度为121.605610，k,17.5,17.8，达到最大最大捕捞总质量，，一(问题重述生态学表明，对可再生资源的开发策略应在事先可持续收获的前提下追求最大经济效益，考虑具有4个年龄组:1龄鱼，……，4龄鱼的某种鱼。
最优捕鱼问题
最优捕鱼策略优化模型摘要“最优捕鱼策略” 的数学模型通过鱼在单位时间内的死亡率来年调整捕鱼强度系数对现有的鱼进行捕捞并获取最大的产量。
由于鱼的生长具有周期性，每一种鱼的数量的改变对整个循环都有影响，因此必须综合考虑，以使每个种年龄段的鱼的数量不破坏的情况下的到最大产量，利用数学知识联系实际问题，作出相应的解答和处理。
问题一：根据已经掌握的人口模型，将鱼的死亡同人口增长联系起来，每种鱼的死亡也有相应的关系，从开始到一个循环的结束，死亡量由大到小，而死亡率保持不变。
通过对死亡率的分析讨论发现)()(t x k r dtdx+-= 经过不定积分可知tk r t e x x )()0()(+-=在此基础上对死亡和捕获量进行综合分析，从而避开了考虑具体的谁先谁后的问题。
通过使用了非线性等式的约束来实现可持续收获，采用了微分方程和非线性规划方法来解决该优化问题。
利用了MATLAB 软件工具求的每年年初的各年龄组鱼的量、最大捕捞量和捕捞强度系数。
得到了各年龄组鱼群的年初的量分别为111019599.1⨯，1110537395.0⨯，,102414672.011⨯7103959.8⨯（单位为条）。
最优的捕捞强度系数为四龄鱼的捕捞强度系数：()年/136279.174=k ，最大量为111088708.3max ⨯=（克）。
在第二问中，模型中通过对鱼群的循环周期考虑可知四年一个循环但模型中将5年作为一个周期来建立模型，这样可以得到最大捕捞量，综合题目一中的模型最终捕在保证破坏最少的情况下的最大产量，由于捞强度系数为未知量，在实现5年后鱼群的生产能力不受到太大破坏的前提下，通过最后一年的量与初始量相等建立模型并利用MATLAB 软件进行求解，求出最大捕捞量，收获的最大量。
求得的捕捞强度系数分别为18.217266（1/年），总收获量为1210604751.1⨯ 克，即160.4751万吨。
关键词：微分方程. 最大捕捞量. 捕捞强度系数. 死亡率. 非线性规划一.问题的提出（略）二.问题分析该问题是一个涉及到微分方程的优化问题，初步分析为非线性规划问题。
建模论文 最优化捕鱼措施
最佳捕鱼策略摘要渔业作为一种再生资源产业，在可持续发展的时代主题下，保证其持续稳产是形势所趋。
本文利用微分方程和非线性规划理论，探讨在可持续收获的条件下，如何通过调整捕捞强度系数，实现捕鱼量的最大化。
针对问题一，首先推导出鱼群产卵、自然死亡、年龄随时间变化等诸因素与各年龄组鱼群数量的数学表达式，结合可持续捕捞，形成一组约束条件，以年捕获量最大作为目标函数，建立非线性规划模型。
用Lingo 编程求解得到：当捕捞强度系数k 取17.36时，年捕获量最大，为3.88×1011克。
然后利用Matlab 画出了在保证可持续捕获的前提下，年度捕获量随捕捞强度系数k 变化的图象，并经过多次计算，验证了结果的准确性和稳定性。
针对问题二，在问题一模型的基础之上，修改约束条件。
首先采用每年的捕捞努力量固定，但各年彼此之间的捕捞努力量不尽相同的方式，然后采用每年的捕捞努力量都保持不变的方式，并将两个模型比较得出采用模型二收益更大。
鉴于此问是多元非线性规划问题，且数据较大，为了得到全局最优解，我们采用Matlab 进行求解，最终得到结果为：1k2k3k4k5kGG13.8815.8818.3633.095.52121.7210⨯得到最大的捕获量为1.72⨯1012克，从而制定出最佳捕鱼策略。
此外，在模型的推广中，改变模型一的假设，在认为4龄鱼一年后仍为4龄鱼的基础上，对问题一进行了改进，得出的结果虽相差甚微，但是思路更具逻辑性。
关键词：微分方程 多元非线性规划 马尔萨斯人口增长模型一、 问题重述为了保护人类赖以生存的自然环境，可再生资源（如渔业，林业资源）的开发必须适度。
一种合理、简化的策略是，在实现可持续收获的前提下，追求最大产量或最佳效益。
考虑对鳀鱼的最优捕捞策略，该种鱼的基本信息如表1所示；表1. 鳀鱼的基本信息1龄鱼 2龄鱼 3龄鱼 4龄鱼 平均重量 5.0711.5517.8622.99自然死亡率 0.8产卵量 00.5545×1051.109×105这种鱼为季节性集中产卵繁殖，产卵和孵化期为每年的最后4个月，卵孵化并成活为1龄鱼，成活率（1龄鱼条数与产卵量n 之比）为1.22ⅹ1011/（1.22ⅹ1011 + n ）.渔业管理部门规定，每年只允许在产卵孵化期前的8个月进行捕捞作业。
最佳捕鱼策略——数学建模论文
最佳捕鱼策略摘要为了实现鳀鱼持续的经济效益，可持续的捕捞方案必不可少。
本文建立了最优化模型，求出了在可持续条件下最大的鳀鱼年收获量以及自然死亡率和捕捞强度系数对模型的影响，并向渔业管理部门提出的鳀鱼资源利用的政策建议。
针对问题一，以一年为周期，年初各个年龄组鳀鱼的数量由上一年相关年龄组的数量决定，分别建立微分方程，得到各个年龄组鳀鱼数量与时间的关系式。
以可持续条件下各个年龄组鳀鱼数量相同为约束条件，以捕捞的3、4龄鱼最大数量为目标函数建立最优化模型。
采用Lingo17.0对模型进行求解，得到年初1龄鱼的数量为1110195994.1⨯条，年初2龄鱼的数量为1010373946.5⨯条，年初3龄鱼的数量为1010414670.2⨯条，年初4龄鱼的数量为710395523.8⨯条，年收获量最大值为1110887536.3⨯克。
针对问题二，由模型I 得出年收获量是自然死亡率和捕捞强度系数的关系。
将捕捞强度系数赋一固定值，用Matlab 软件得出了在4龄鱼的捕捞强度系数为5的情况下，年收获量和自然死亡率成反向关系。
针对问题三，由前述得到的年收获量与自然死亡率和捕捞强度系数的关系，运用Matlab2016求解得到当4龄鱼的捕捞强度系数(k)以0.01为步长，从0到20分布时对应的F(k)的数值，并以k 的取值为横坐标，对应的F(k)为纵坐标，绘制捕获量F(m)随捕捞强度系数变化的曲线图，得出年收获量与捕捞强度系数成正向关系。
最后，本文从提高捕捞技术、保护鳀鱼苗种和生存环境、开发产业链等四个方面对鳀鱼资源的综合利用提出了建议。
关键词：年收获量最优化模型1问题重述和分析本题是最优化问题，此问涉及的各个变量为：每条1龄鱼、2龄鱼、3龄鱼、4龄鱼的平均重量分别是 5.1g、11.6g、17.9g、23.0g，自然死亡率为0.8，各个年龄组鳀鱼产卵量情况，产卵孵化期为每年后4月，3龄鱼和4龄鱼捕捞强度系数比为0.42：1，卵的存活率等。
【实验】数学建模实验报告最优捕鱼策略
【关键字】实验最优捕鱼策略一.实验目的：1、了解与熟练掌握常系数线性差分方程的解法；2、通过最优捕鱼策略建模案例，使用MA TLAB软件认识与掌握差分方程模型在实际生活方面的重要作用。
二.实验内容：（最优捕鱼策略）生态学表明，对可再生资源的开发策略应在事先可持续收获的前提下追求最大经济效益。
考虑具有4个年龄鱼：1龄鱼，… ，4龄鱼的某种鱼。
该鱼类在每年后4个月季节性集中产卵繁殖。
而据规定，捕捞作业只允许在前8个月进行，每年投入的捕捞能力固定不变，单位时间捕捞量与个年龄鱼群条数的比率称为捕捞强度系数。
使用只能捕捞3、4龄鱼的网眼的拉网，其两个捕捞强度系数比为0.42:1.渔业上称这种方式为固定力量捕捞。
该鱼群本身有如下数据：1．各年龄组鱼的自然死亡率为0.8（1/年），其平均质量分别为5.07,11.55,17.86,22.99（单位：g）；2．1龄鱼和2龄鱼不产卵，产卵期间，平均每条4龄鱼产卵量为1.109ⅹ105（个），3龄鱼为其一半；3．卵孵化的成活率为1.22ⅹ1011/（1.22ⅹ1011 + n）（n为产卵总量）；有如下问题需要解决：1）分析如何实现可持续捕获（即每年开始捕捞时各年龄组鱼群不变），并在此前提下得到最高收获量；2）合同要求某渔业公司在5年合同期满后鱼群的生产能力不能受到太大的破坏，承包时各年龄组鱼群数量为122,29.7,10.1,3.29（ⅹ109条），在固定努力量的捕捞方式下，问该公司应采取怎样的捕捞策略，才能使总收获量最高。
三. 模型建立假设a、鱼群总量的增加虽然是离散的，但对大规模鱼群而言，我们可以假设鱼群总量的变化随时间是连续的；b、龄鱼到来年分别长一岁成为i + 1龄鱼，i = 1，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.86t∈[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（+）四. 模型求解（含经调试后正确的源程序）1.先建立一个buyu.m的M文件：function y=buyu(x);global 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');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;ktotal=-mtotal;a10=eval(a10)a20=eval(a20)a30=eval(a30)a40=eval(a40)format shortclear五．结果分析1．鱼总量与时间图：2.可以看出捕捞强度对收获量的影响：实验输出数据：y =-3.6757e+011y =-3.9616e+011y =-4.0483e+011y =-4.0782e+011y =-4.0802e+011y =-4.0805e+011y =-4.0805e+011y =-4.0805e+011y =-4.0805e+011y =-4.0805e+011y =-4.0805e+011y =y =-4.0667e+011k =18.25976795085083total =4.080548655562244e+011 a10 =1.195809275167686e+011a20 =5.373117428928620e+010a30 =2.414297288420686e+010a40 =8.330238542343275e+007则k=18.25976795085083时，最高年收获量为total=4.080548655562244×1011（克），此时每年年初1，2，3，4年龄组鱼的数量分别为：1.195809275167686×10115.373117428928620×10102.414297288420686×10108.330238542343275×107六．实验总结本次实验的目的是了解差分方程（递推关系）的建立及求解，以及掌握用差分方程（递推关系）来求解现实问题的方法。
最佳捕鱼策略
2015高教社杯全国大学生数学建模竞赛承诺书我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。
如有违反竞赛规则的行为，我们将受到严肃处理。
我们参赛选择的题号是（从A/B/C/D中选择一项填写）：我们的参赛报名号为（如果赛区设置报名号的话）：所属学校（请填写完整的全名）：参赛队员(打印并签名) ：1.2.3.指导教师或指导教师组负责人(打印并签名)：日期：年月日最佳捕鱼策略目录摘要 (2)1.问题的重述 (4)2.模型的假设与符号的说明 (5)2. 1 模型假设 (5)2. 2 符号说明 (5)3 模型建立和求解 (6)3.1问题1 (6)3. 2问题2 (11)4 模型评价及改进 (16)摘要针对本文要求的最大收益时的捕捞强度系数和各年龄鱼群相应各年数量，本文利用微分方程建模，得到各年龄鱼群相应各年数量与捕捞强度系数的关系如式（8），然后进行后续建模计算。
对于问题一利用平均死亡率和瞬时死亡率的关系我们得到瞬时死亡率值，本文根据其定义利用微分方程建立鱼群数量的变化模型。
对于3、4龄鱼，还应考虑前8个月的捕捞，由此建立分段函数模型。
第二年的1龄鱼是由上一年的3龄鱼、4龄鱼产卵所得，我们注意到，产卵鱼群应该分为产卵后活到下一年的鱼和产卵后还未活到下一年的鱼。
本文假设每时刻鱼死亡的概率均相等，则在产卵期的四个月里的某一时刻t，鱼死亡的概率服从几何概率分布，由此得到产卵后还未活到下一年的鱼的数量，最终得到3龄鱼和4龄鱼的产卵总数，便得到第二年1龄鱼的数量。
最后得到的捕捞的3、4龄鱼的数量和每条3、4龄鱼的重量列出捕捞总量即目标函数，结合1、2、3、4龄鱼的数量的等式关系作为约束条件，利用matlab计算进行优化计算得到捕捞量最大（为3.9512×10^11g）时的四龄鱼的捕捞强度系数(为16.05)和1、2、3、4龄鱼初始数量（分别为117.23×10^9、35.17×10^9、10.55×10^9、0.04×10^9条）。
最优捕鱼策略模型
求解 (1) (2) (3)式 ,得到 :
fi ( t) = e - kt fi ( 0) , i = 1, 2,
fi ( t) =
e-
( k + xi) t fi ( 0) ,
( 0 ≤ t≤
2 3
)
e- k(t-
2 3
)
fi (
2 3
),
(
2 3
≤ t≤1 )
可求得
fi (
2 3
)
= e-
2 3
nj:第 j年产卵总数量 ;
η j
=
1.1.2222××1100111η1 j为第
j年卵的存活率 。
2. 问题分析 。
问题一 :要求在可持续捕捞的基础上 ,实现年捕捞重量最大 。在鱼
的重量一定的情况下 ,捕捞量由 3龄鱼和 4 龄鱼捕捞的数量 (即题目中
的捕捞强度系数 )确定 ,另外还需要求出各年龄组的鱼的数量 。这里设
问题二 :这里同样需要求最大捕捞量 ,注意到这里仍然是固定努力
1996年大学生数学建模竞赛试题(最优捕鱼)市公开课获奖课件省名师示范课获奖课件
渔业管理部门要求，每年只允许在产卵卵化期前旳８个月内进行捕捞 作业。假如每年投入旳捕捞能力（如渔船数、下网次数等）固定不变，这 时单位时间捕捞量将与各年龄组鱼群条数成正比。百分比系数不妨称捕捞 强度系数。一般使用１３mm网眼旳拉网，这种网只能捕捞３龄鱼和４龄鱼， 其两个捕捞强度系数之比为0.42:1。渔业上称这种方式为固定努力量捕捞。
x1
(t
1)
1.22 1011 n(t 1) 1.22 1011 n(t 1)
x2 (t 1) er x1 (t) x3 (t 1) er x2 (t) x4 (t 1) x3 (t)e (k3 2 / 3r)
第六年1龄鱼数量占第一年1龄鱼数量旳百分比为：
q1 x1(6) 100% 1.1956 100% 98%
: 孵化存活率
三、模型建立
问题一模型
1．1，2龄鱼及3，4龄鱼后四个月旳生长只受自然死亡率旳影响，由此可知1，2龄
鱼旳生长旳微分方程满足方程（1）：
dxi dt
rxi , i
1, 2,3, 4
可得：x i
(t
)
x0ert
,
x0为每年年初i龄鱼的数量
T年旳i龄鱼在T+1年变为i+1龄鱼， 则有：xi1(T 1) er xi(T )
(0.28k 6.4 )
(m*(0.5*e 3
(
e
2 3
k
0.28k 8.8 3
( 2 k 0.8)
)
)
1)
1e 3
问题一：
model: max=17.86*0.42*k/(0.8+0.42*k)*1.22*10^11/(1.22*10^11+n)*n*@exp(-1.6)*(1-@exp(-
最优捕鱼策略-数学建模
西安邮电大学（理学院）数学建模报告最优捕鱼策略专业名称：信息与计算科学班级： 1302班学生姓名：张梦倩学号（8位）： 07131057指导教师：支晓斌摘要为了保护人类赖以生存的自然环境，可再生资源（如渔业、林业资源）的开发必须适度。
本文实际上就是为了解决渔业上最优捕鱼策略问题，即在可持续捕捞的前提下，追求捕捞量的最大化。
问题一采用条件极值列方程组的方法求解，即1龄鱼的数量由3龄鱼和4龄鱼的产卵孵化而来；2，3龄鱼的数量分别由上一年1龄鱼，2龄鱼生长而来；4龄鱼由上一年的3龄鱼和上一年末存活的4龄鱼组成。
最后得到：只要每年1-8月份3、4龄鱼捕捞总量小于、，就可以实现总捕捞量最大为；对结果分析得到捕捞的对象主要是3龄鱼，当3龄与4龄鱼的捕捞系数发生变化时，总的捕捞量变化不大。
问题二给出年初各龄鱼的数量，要求在5年后鱼群的生产能力没有受到太大的破坏的前提下，使5年的总收获量最大，即在5年内鱼群能够可持续繁殖和生长。
本题以5年的总捕获量为目标函数，以5年后各龄鱼的数量没有发生太大的变化为条件，建立承包期总产量模型。
最终得到的捕捞策略如表1-1。
只要各年龄鱼每年的捕捞数量小于表1-1中的数量，就可以实现5年后鱼群的生产能力没有发生太大的变化。
一、问题重述为了保护人类赖以生存的自然环境，可再生资源（如渔业、林业资源）的开发必须适度。
一种合理、简化的策略是，在实现可持续收获的前提下，追求最大产量或最佳效益。
考虑对某种鱼（鲳鱼）的最优捕捞策略：假设这种鱼分4个年龄组：称1龄鱼，……，4龄鱼。
各年龄组每条鱼的平均重量分别为5.07,11.55,17.86,22.99（克）；各年龄组鱼的自然死亡率均为0.8（1/年）；这种鱼为季节性集中产卵繁殖，平均每条4龄鱼的产卵量为1.109×105(个）；3龄鱼的产卵量为这个数的一半，2龄鱼和1龄鱼不产卵，产卵和孵化期为每年的最后4个月；卵孵化并成活为1龄鱼，成活率（1龄鱼条数与产卵总量n之比）为1.22×1011/(1.22×1011+n).渔业管理部门规定，每年只允许在产卵卵化期前的8个月内进行捕捞作业。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
最优捕鱼策略
一.实验目的：
1、了解与熟练掌握常系数线性差分方程的解法；
2、通过最优捕鱼策略建模案例，使用MATLAB软件认识与掌握差分方程模型在实际生活方面的重要作用。
二.实验内容：（最优捕鱼策略）
生态学表明，对可再生资源的开发策略应在事先可持续收获的前提下追求最大经济效益。
考虑具有4个年龄鱼：1龄鱼，…，4龄鱼的某种鱼。
该鱼类在每年后4个月季节性集中产卵繁殖。
而据规定，捕捞作业只允许在前8个月进行，每年投入的捕捞能力固定不变，单位时间捕捞量与个年龄鱼群条数的比例称为捕捞强度系数。
使用只能捕捞3、4龄鱼的13mm网眼的拉网，其两个捕捞强度系数比为0.42:1.渔业上称这种方式为固定力量捕捞。
该鱼群本身有如下数据：
1．各年龄组鱼的自然死亡率为0.8（1/年），其平均质量分别为
5.07,11.55,17.86,22.99（单位：g）；
2．1龄鱼和2龄鱼不产卵，产卵期间，平均每条4龄鱼产卵量为1.109ⅹ105（个），3龄鱼为其一半；
3．卵孵化的成活率为1.22ⅹ1011/（1.22ⅹ1011 + n）（n为产卵总量）；
有如下问题需要解决：
1）分析如何实现可持续捕获（即每年开始捕捞时各年龄组鱼群不变），并在此前提下得到最高收获量；
2）合同要求某渔业公司在5年合同期满后鱼群的生产能力不能受到太大的破坏，承包时各年龄组鱼群数量为122,29.7,10.1,3.29（ⅹ109条），在固定努力量的捕捞方式下，问该公司应采取怎样的捕捞策略，才能使总收获量最高。
三. 模型建立
假设a、鱼群总量的增加虽然是离散的，但对大规模鱼群而言，我们可以假设鱼群总量的变化随时间是连续的；b、龄鱼到来年分别长一岁成为i + 1龄鱼，i = 1，2，3；c、4龄鱼在年末留存的数量占全部数量的比例相对很小，可假设全部死亡。
d、连续捕获使各年龄组的鱼群数量呈周期性变化，周期为1年，可以只考
虑鱼群数量在1年内的变化情况。
（且可设x i （t ）：在t 时刻i 龄鱼的条数，i = 1，2，3，4；n ：每年的产卵量；k ：4龄鱼捕捞强度系数；2a i0：每年初i 龄鱼的数量，i = 1，2，3，4；）
进而可建立模型如下：
max （total （k ））=17.86⎰⎰+3
/203
/2043)(99.22)(42.0dt t kx dt t kx
)(8.0)(11t x dt
t dx -= t ∈[0，1]，x1（0）= n ×n +⨯⨯1111
1022.11022.1 )(8.0)(22t x dt t dx -= t ∈[0，1]，x2（0）= x1（1）
)()42.08.0()(33t x k dt t dx +-= t ∈[0，2/3]，x3（0）= x2（1）
s.t. )(8.0)(33t x dt t dx -= t ∈[2/3，1]，x3（32-）= x3（32+）
)()8.0()(44t x k dt t dx +-= t ∈[0，2/3]，x4（0）= x3（1）
)(8.0)(44t x dt t dx -= t ∈[2/3，1]，x4（32-）= x4（32+）
四. 模型求解（含经调试后正确的源程序）
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');
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
五．结果分析
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 =
-4.0805e+011
y =
-4.0805e+011
y =
-4.0805e+011
y =
-4.0805e+011
y =
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时，最高年收获量为total=4.080548655562244×1011（克），此时每年年初1，2，3，4年龄组鱼的数量分别为：
1.195809275167686×1011
5.373117428928620×1010
2.414297288420686×1010
8.330238542343275×107
六．实验总结
本次实验的目的是了解差分方程（递推关系）的建立及求解，以及掌握用差分方程（递推关系）来求解现实问题的方法。
实验中假设鱼群总量的变化随时间是连续的，从而利用微分方的知识建立最优捕鱼策略问题的优化模型。
通过实验加深了对概念和方法的理解，了解了差分方程的程序解法。
学生签名：
年月日
七．教师评语及成绩
教师签名：
年月日。