最优捕鱼策略实验报告
学号：104080298 姓名：宁亚会 班级：10D
摘要
为了保护人类赖以生存的自然环境，可再生资源（如渔业、林业资源等）的开发必须适度。而在社会经济生活中,我们要使商业活动在一段时期内达到最大收益,因此我们要合理的开发资源，这时,我们不仅要考虑商业活动的当前经济效益,还要考虑生态效益及由此产生的对整体经济效益的影响。本文就是对渔业这类可再生资源的开发问题进行研究，利用相关的数学软件进行求解。
对于问题一，我们考虑渔场生产过程中的各年龄组鱼群数量的制约因素，将其分为两大类，第1,2龄鱼群为一类，该鱼群数量变化在一年内只受自然死亡率制约，写出鱼群数量满足的微分方程；第3,4龄鱼群为一类，其数量变化在前8个月受捕捞强度和自然死亡率影响，后4个月只受自然死亡率的制约，分阶段写出写出鱼群数量满足的微分方程；根据微分方程，求出在某时刻各鱼群的数量表达式（类似于人口增长模型）。因为捕捞是连续的，所以任意一个时刻的捕捞量为捕捞强度乘以鱼群的数量，又捕捞只在前8个月进行，则年捕捞量为前8个月各时刻鱼群数量的积分。最后建立年总捕捞量的函数与生产过程中满足的关系式，转化为非线性规划模型，利用lingo 和matlab 软件分别求解。
对于问题二，题中已给出各年龄组鱼群的初始值，我们利用问题一中所得到的迭代方程，可迭代地求出第i 年初各年龄组鱼群的数量；再根据问题一中的捕捞量表达式，可写出5年的捕捞总量表达式，以5年捕捞总量最大为前提，利用matlab 软件求解出此时的捕捞强度，然后再验证在此捕捞强度下会不会使5年后鱼群的生产能力有太大的破坏。
最后，我们得出以下结论：可持续捕获条件下，捕捞强度为17.36292时，达到最大捕捞总质量11
3.88707610⨯g ； 5年后鱼群的生产能力不会有太大的破坏条件下，捕捞强度为
()17.5,17.8k ∈，达到最大最大捕捞总质量121.605610⨯
一．问题重述
生态学表明，对可再生资源的开发策略应在事先可持续收获的前提下追求最大经济效益，考虑具有4个年龄组：1龄鱼，……，4龄鱼的某种鱼。该鱼类在每年后4个月季节性集中产卵繁殖。而按规定，捕捞作业只允许在前8个月进行，每年投入的捕捞能力固定不变，单位时间捕捞量与各年龄组鱼群条数的比例称为捕捞强度系数。使用只能捕捞3，4龄鱼的13mm 网眼的拉网，其两个捕捞强度系数比为0.42：1。渔业上称这种方式为固定努力量捕捞。
该鱼群本身有如下数据：
各年龄组鱼的自然死亡率为0.8（1/年），其平均质量分别为5.07，11.55，17.86，22.99（单位：g ）；1，2龄鱼不产卵，平均每条4龄鱼产卵量为5
10109.1⨯（个），3龄鱼为其一半；卵孵化的成活率为)1022.1(10
22.11111
n +⨯⨯（n 为产卵总量）；
二、问题分析
对于问题一，要实现可持续捕获，即每年开始捕捞时渔场中各年龄组鱼群条数不变，因此我们要算出每年初各龄鱼组的数量。
对自然死亡率，捕捞强度系数和卵成活率的解释，即1，2龄鱼仅受自然死亡率的影响；而3，4龄鱼不仅受自然死亡率的影响，还受捕捞强度系数的影响；因为该种鱼的最高寿命为4，所以在后四月中4龄鱼都不存活；，而对于1龄鱼的数量，是3，4龄鱼在前年的后4年产卵所存活下来的数量；对于捕捞量，题中规定只在1到8月才能捕捞，而且1，2龄鱼不被捕捞，所以主要来源于对3，4龄鱼的捕捞。根据这些关系可列出一系列的方程，其中捕捞量作为目标函数，其他的作为约束条件，建立一个非线性规划模型，再然后用lingo 软件和matlab 软件进行求解。
对于问题二，合同要求5年后鱼群的生产能力不能受到太大破坏，又要使总收益最高，这就有可能发生满足了前者满足不了后者之类的情况。我们处理方法是先确定一个策略使其收益最高，再检验此捕鱼策略是否能保证5年后鱼群的生产能力不受到太大的破坏，若它让鱼群的生产能力受到了严重破坏，我们再求另外一种策略。但从理论分析可知，5年后将在鱼群尽可能接近可持续鱼群的情况下来使捕捞量达到最大。对于破坏大小，我们采用1龄鱼群数量变化率来衡量，即以第六年初1龄鱼群数量的变化量与承包时鱼群数量初值之比表示，因为2，3，4龄鱼群的数量在很大程度上受承包初1龄鱼影响，根据关系，可以知道5年后2，3，4龄鱼群的数量肯定会有较大变化。只要该比值小于5%，我们就认为鱼群的生产能力没有受到太大破坏。
题中已经给了我们各年龄组的初始值，而问题一中也已得出一组迭代方程，我们利用这些迭代方程，求出各年的鱼量分布；同样可以根据问题一中捕捞量的表达式求出5年的总捕捞量，以此来确定我们的最优捕捞策略。再然后我们通过验证来确定其5年后鱼群的生产能力有没有受到太大破坏。
综上所述，原问题实质上是给出了各年龄组鱼群之间数量的变化规律，并给出了它们的自然死亡率及捕捞和产卵的时间分布，并固定3、4龄鱼捕捞能力的比值，要求选择一定的捕捞能力系数，使得各年龄组鱼的数量在各年开始的第一天条数不变（第一问），5年后鱼群的生产能力不会有太大的破坏（第二问），并在此条件下，求到最大捕获量。
三．模型假设
由于问题本身存在很多不确定因素，为了使问题简化，作如下假设：
1．这种鱼在一年内的任何时间都会发生自然死亡，即死亡是一个连续的过程。 2．捕捞也是一个连续的过程，不是在某一时刻突然发生。 3．1、2龄鱼体形太小，不能被捕。
4．3、4龄鱼在一年中的后4个月的第一天集中一次产卵
5．i 龄鱼到来年分别长一岁成为i+1龄鱼，i=1，2，3，其中上一年存活下来的4龄鱼仍是4龄鱼
四．问题一的模型建立与求解
4.1模型建立
1，2龄鱼的生长只受自然死亡率的影响，由此可知1，2龄鱼的生长的微分方程满足方程
4,3,2,1,=-=i rxi dt
dxi
龄鱼的数量为每年年初可得：i ,)(00x e x t x rt i -= )1(
T 年的i 龄鱼在T+1年变为i+1龄鱼，
)()1(1T x e T x i r i -+=+则有： )2(
