1、 数学建模练习 第 7 次练习 姓名：程涛 胡铭 李冲 题目：最优捕鱼策略问题 最优捕鱼策略问题 摘 要 本文主要解决在在保护自然环境、实现可持续捕捞的前提下，针对鳀鱼提出 合理的捕捞方案以实现最大产量或最佳效益， 同时为渔业部门制定相关规定提出 参考依据。 针对问题一，在实现可持续捕捞的前提下，得到最高的年收获量。将鳀鱼分 为4个年龄段，分别建立微分方程，得到各龄组鳀鱼数量与时间的关系式，再以 每年开始时渔场中各年龄组鱼群条数不变为约束条件， 以总捕捞量最大为目标函 数，建立非线性规划模型，运用 MATLAB 软件求得当 3 龄鱼、4 龄鱼的捕捞强度 分别为 7.51、 17.88时， 捕捞
2、3龄鱼、 4龄鱼对应的条数为 、 条，此时捕捞总量最大为： 。并且得到年初各年龄组鱼群的条数分 别为 、5.36 、 、 条。 10 2.17 10 7 6.95 10 11 3.87 10 g 10 10 2.41 10 11 1.19 10 10 7 7.26 10 针对问题二，要求 5 年后鱼群的生产能力不能受到太大破坏，使得 5 年总 收获量最高。 以1龄鱼第5年末的数量与第1年初数量的相对误差小于10%作为 鱼群生产能力没有受到太大破坏的指标。 利用第一问中所得到的各年龄鱼数量随 时间的关系，以鱼群生产能力未受到太大破坏的指标为约束条件，5年总捕捞量 最大为目标函数，建立非线性规划模
3、型，运用LINGO软件求得当3龄鱼、4龄鱼 的捕捞强度分别为7.38、 17.58 时， 每年的捕捞量约为： 、 、 、 、 。此时 5 年捕捞总量最大约为： ，并且得到每年初各年龄组鱼群数量见表三。 11 2.39 10 11 2.2 g 10 11 4.05 10 g 12 1.6 10 g 11 3.87 10 g 11 3.9 10 g 关键词：最优捕鱼；微分方程；非线性优化模型；LINGO软件；MATLAB软件1 一 问题重述 1.1 问题背景 为了保护人类赖以生存的自然环境，可再生资源（如渔业，林业资源）的开 发必须适度。一种合理、简化的策略是，在实现可持续收获的前提下，追求最大
4、产量或最佳效益。 考虑对鳀鱼的最优捕捞策略： 假设这种鱼分 4 个年龄组，称为 1 龄鱼，2 龄鱼，3 龄鱼，4 龄鱼。各年龄 组每条鱼的平均重量分别为 5.07，11.55，17.86，22.99（克） ，各年龄组的自 然死亡率为0.8（/年） ，这种鱼为季节性集中产卵繁殖，平均每条4龄鱼的产卵 量为 （个） ，3 龄鱼的产卵量为这个数的一半，2龄和1龄鱼不产卵， 产卵和孵化期为每年的最后4个月，卵孵化并成活为1龄鱼，成活率（1龄鱼条 数与产卵量 之比）为 5 1.109 10 n () 11 11 1.22 10 1.22 10 n + 。 1.2 问题描述 渔业管理部门规定，每年只允许在
5、产卵孵化期前的8个月进行捕捞作业。如 果每年投入的捕捞能力（如渔船数，下网次数等）固定不变，这时单位时间捕捞 量将与各年龄组鱼群条数成正比， 比例系数不妨称捕捞强度系数。 通常使用13mm 网眼的拉网， 这种网只能捕捞3龄鱼和4龄鱼， 其两个捕捞强度系数之比为0.42 1。渔业上称这种方式为固定努力量捕捞。 1）建立数学模型分析如何实现可持续捕捞（即每年开始捕捞时渔场中各年 龄组鱼群条数不变） ，并且在此前提下得到最高的年收获量（捕捞总重量） 。 2）某渔业公司承包这种鱼的捕捞业务5年，合同要求5年后鱼群的生产能 力不能受到太大破坏。 已知承包时各年龄组鱼群的数量分别为： 122， 29.7，
6、 10.1， 3.29（ 条） ，如果仍用固定努力量的捕捞方式，该公司应采取怎样的策略才 能使总收获量最高。 9 10 二 问题分析 本为需要解决的根本问题在于：如何合理的制定捕捞方案既能达到最佳收 益，又能不太大破坏渔场生态系统，实现可持续收获。 对于问题一，在实现可持续捕捞的前提下，得到最大的年收获量。题中所给 的自然死亡率为0.8（/年）理解为瞬时死亡率。对于1龄鱼、2龄鱼由于不存在 捕捞，建立鱼群数量与时间的微分方程，得到1龄鱼、2龄鱼数量随时间的关系。 对于3龄鱼、4龄鱼，每年前8个月不仅要考虑自然死亡率，还要考虑捕捞量； 后四个月只需考虑死亡率，建立鱼群数量随时间的变化关系，得到3
7、龄鱼、4龄 鱼数量随时间的关系。假定产卵在9月初集中进行，由于只有3龄鱼、4龄鱼进 行产卵且每年只允许在产卵孵化期前的8个月进行捕捞作业。根据3龄鱼、4龄 鱼数量随时间的关系，求得每年产卵前3龄鱼、4龄鱼的数量，这样便可得到该2年的产卵量，结合成活率便得到年末1龄鱼数量。年末后低龄鱼变为高龄鱼，根 据年龄组鱼群条数不变得到约束条件，建立非线性优化模型，运用 MATLAB 便可 求得3龄鱼和4龄鱼的捕捞强度、年最大捕捞总量。 对于问题二，要求 5 年后鱼群的生产能力不能受到太大破坏，使得 5 年总 收获量最高。假设每年的捕捞强度系数不变。考虑到5年后鳀鱼总数量主要取决 于第1年初1龄鱼的数量，
8、将1龄鱼第5年末的数量与第1年初数量的变化率较 小作为鱼群生产能力没有受到太大破坏的指标， 利用第一问中所得到的各年龄鱼 数量随时间的关系，建立非线性规划模型，运用LINGO便可以求得5年最大总捕 捞量、捕捞强度。 三 模型假设 1）该鱼场的鳀鱼生态系统是独立的生态系统，不受其他生态系统的影响。 2）鳀鱼各年龄组的自然死亡率已经包含了鳀鱼群体与其他生物之间的竞争。 3）各年龄组的鱼经过一年后即进入高一年龄组，4龄鱼经过一年后仍为4龄鱼。 4）不考虑水质污染、气候异常等自然灾害对鳀鱼生长的影响。 5）3龄鱼、4龄鱼在9月初集中产卵。 6) 问题二中，假设五年3、4龄鱼的捕捞强度系数相同。 四 符
9、号说明 t ：表示时刻。 T ：表示年份。(= 1 , 6 ) T () i xt：表示t 时刻i龄鱼的数量。 ()=1,2,3,4 i i x ：表示初始时i龄鱼的数量。 ()=1,2,3,4 i i w :表示i龄鱼个体平均重量。 ( ) / 位：克 单条 r：表示各年龄鱼的自然死亡率。( 0 . 8 r = ) 4 n ：表示4龄个体平均产卵量。 ()5 4 = .109 10 n 1 个 n ：表示年产卵总量。 i f ：表示每年i龄鱼的捕捞数量。 ( ) 位： 单条 f ：表示年捕捞总重量。 ( ) 位：克 单 i k ：表示 龄鱼捕捞强度系数。 ()4 =3,4 i e：表示卵孵化
10、并成活为1龄鱼的成活率。 3五 模型建立与求解 5.1问题一模型建立与求解 5.1.1模型建立 由题已知 、 。并且对 1、2 龄鱼不能进行捕捞，其生 长只受自然死亡率的的影响。而题目中所给的各年龄组的自然死亡率为 0.8（/ 年） ，将死亡率理解为瞬时死亡率。由此求得1、2龄鱼的生长微分方程： 3 =17.86g w 4 =22.99g w (t) =( )(1 , i i dx rx t i dt =2 ) 2 ) 4 k并且，当 。 0( ) =( = 1 , ii tx t x i = 时 可得： （1） () ( 1 ,2) i rt ii xt xe i - = 对于 3、4 龄鱼
11、在前 8 个月不仅受到自然死亡率的影响，还要受到捕捞强度 系数的影响，并且捕捞强度系数的关系为： 。而当处于后 4 个月时进 入卵孵化期则只受自然死亡率的影响。 3 =0.42 k 捕捞期3、4龄鱼数量满足的微分方程： (t) =( )()( 3 , 4 ) i ii dx rkxti dt += 2 0, 3 t 并且，当 0( ) =( = 3 , ii tx t x i =时 4 ) 可得： () () ( 3,4) i rkt ii xt xe i -+ = 2 0, 3 t （2） 解得， 2 = 3 t 时 2 () 3 2 () = 3 i rk ii xx e -+由于每年的捕
12、捞只在 1 月到 8 月进行，且只能捕捞到 3、4 龄鱼，设在微小 的时间段dt 内，捕捞量为 ，则： i df =( ) iii df k x t dt 则 8 个月总的捕捞量为： 2 3 0 =( ii i fk x t ) d t 2( + ) 3 =1e (3 , + i kr ii i kx i kr - - 4 ) = （3） 产卵孵化期3、4龄鱼数量满足的微分方程： () =( ) (3 , 4 i ii dx t kx t i dt = ) 2 ,1 3 t 并且，当 2 (+) 3 22 () =e (= 3 , 4 ) 33 i rk ii txxi - = 时 解得： 4
13、 2 3 () e ( 3,4) i kr t ii xt x i - = 2 ,1 3 t （4） 其中总产卵量为： 3 22 () ( 3 43 44 0.5 rk rk nn x e n x e -+ -+ =+ 4 ) 3（5） 又卵孵化并成活为 1 龄鱼的成活率为： 11 11 1.22 10 1.22 10 n e = + （） （6） 则年末1 龄鱼的数量为： 1 (1) xn e = 建立非线性规划模型： 目标函数： 使得在该年内捕捞总重量最大，即： 33 44 max + fwfwf = 约束条件： 1） 第一年末， 由卵孵化并成活为 1 龄鱼的数量等于第一年初 1 龄鱼的数
14、量， 即： 1 = nx e 2）第一年末，i龄鱼的数量等于第一年初高一龄的数量，即： +1 (1)= ( =1,2) ii xxi 3）第一年末，3、4 龄鱼的数量之和等于第一年初4 龄鱼的数量，即： 34 (1)+ (1)= xx 4 x故模型可以总结为如下优化模型： 33 44 max + wf wf34 1 0.8 12 0.8 23 22 33 34 34= e= e= e+ e= =0.42 kk rr st nx xx xx xxx kk e - - - - 45.1.2 模型求解与结果分析 代入相关数据，由式(1) ，运用 MATLAB 求解得到：3 龄鱼、4 龄鱼的 捕捞强度
15、系数分别为k 7.51、 17.88，此时捕捞总重量最大为： ；对应的各年龄段鱼群数量见表一： (6) 3 = 4 k = 11 3.87 10 (g) f=各年龄段鱼群数量 表一 鱼年龄段 1 2 3 4 数量 11 1.19 10 10 5.36 10 10 2.41 10 7 7.26 10 3 龄鱼、4 龄鱼捕捞数量如下表二： 53、4龄鱼捕捞数量 表二 鱼年龄段 3 4 数量（条） 10 2.17 10 7 6.95 10 同时运用MATLAB 画出捕捞强度系数与总捕捞量的关系图,即图一: 图一 捕捞强度系数与总捕捞量 根据图一可以得出：在一定的范围内总捕获量随着捕捞强度的增大而增大
16、； 当超过某一最大之后便随着捕捞强度的增大而减小， 而那最高点即为所求的最高 的年收获量。这与实际吻合的比较好，当过度捕捞时会破坏鱼群生态的平衡。可 见该模型具有一定的合理性，能够为实际生产活动提供一定的参考。 5.2问题二模型建立与求解 5.2.1模型建立 为方便模型的建立，假设不同年份以相同捕捞强度系数捕捞，并且 3 龄鱼与 4 龄鱼捕捞强度系数之比为0.42:1。 由于五年后鳀鱼总数量主要取决于第一年初 1 龄鱼的数量，将 1 龄鱼第 5 年 末的数量与第 1 年初数量的变化率小于10%作为鱼群生产能力没有受到太大破 坏的指标。 在问题一的基础上，可以求出时间 时各龄段组鱼的数量随时间
17、的函数如下： ,+ 1 tT T 6 3 3 4 4 () 11 () 22 (+) (-) 33 2 (-) 3 33 (+) (-) 44 2 (-) 3 44 ()= ()e , + 1 ()= ()e , + 1 2 ()= ()e + 3 2 () =() e +,+ 1 3 2 ()= ()e + 3 2 () =() e +,+ 1 3 rtT rtT krtT krtT krtT krtT xtxT t TT xtxT t TT xtxT t T , T xtxT t T T xtxT t T , T xtxT t T T - - - - - - （7） 则有，第T 年内总的产
18、卵量为： 34 22 () ( 3 43 44 ()0 . 5 () () rk rk nT nx Te nx Te -+ -+ =+ ) 3（8） 第T 年内卵孵化并成活为 1 龄鱼的成活率为： 11 11 1.22 10 () 1.22 10 ( ) T nT e = + （） （9） 并且可以得到从第T 年初到第 年初这一年的时间内，捕捞 3 龄鱼和 4 龄鱼数量分别为： +1 T3 4 2( + ) 33 3 3 3 2( + ) 44 3 4 4 () () = 1e + () () = 1e + kr kr kx T fT kr kx T fT kr - - - - )（10） 则
19、第T 年内总的捕捞量为： （11） 33 44 () = () + () fT wfT wfT 建立非线性规划模型： 目标函数： 五年总的捕捞量最大，即： 5 =1 max ( ) T fT 约束条件： 1）第 +1 T 年年初的 () +1 i 龄鱼等于第T 年末的i龄鱼，即： * +1 (+ 1 ) =()(= 1 , 2 ) ii xTxTi 其中， 表示i龄鱼第T 年年末的数量。 * ( i xT 2） 第 年年初1龄鱼的数量等于第T 年末卵孵化并成活为1龄鱼的数量， 即： +1 T 1 (+ 1 ) =()() xT nTT e 3）第 年年初 4 龄鱼的数量等于第T 年年末 3 龄
20、鱼与 4 龄鱼之和，即： +1 T* 434 (+ 1 ) =() +() xT xT xT74）以1 龄鱼第 5 年末的数量与第 1 年初数量的相对误差小于10 作为鱼群 生产能力没有受到太大破坏的指标，即： %11 1 (6) (1) 0.10 (1) xx x - 又由题已知： 11 10 10 9 1234 (1)=1.22 10 (1)=2.97 10 (1)=1.01 10 (1)=3.29 10 xxxx 故模型总结为如下优化模型： 5 =1 * +1 1 * 434 11 1 max ( )(+ 1 ) =()(= 1 , 2 ) (+ 1 ) =()() (+ 1 ) =()
21、 +() (6) (1) 0.10 (1) T ii fT st xTxTi xT nTT xT xT xT xx x e - 5.2.2模型求解与评价 代入相关数据，由式(7 ，已知 的初值，逐步迭代，求出 的表达式，再运用 LINGO 软件求得：3 龄鱼、4 龄鱼的捕捞强 度分别为 7.38、 k ， 此时 5 年捕捞总重量最大为： ； 每年3 龄鱼、4 龄鱼捕捞量见下表三： ) (11) 17.58 (1) i x ()(= 2 , 6 ) i xTT 3 k = 4 = 12 1.6 10 g f= 每年 3、4龄鱼捕捞量和年总捕捞量 表三 年份 第一年 第二年 第三年 第四年 第五年
22、 3龄鱼捕捞 （g） 11 1.66 10 11 2.19 10 11 4.05 10 11 3.85 10 11 3.90 10 4 龄鱼捕捞(g） 10 7.30 10 8 7.34 10 8 9.70 10 9 1.79 10 9 1.70 10 捕捞总量（g） 11 2.39 10 11 2.20 10 11 4.05 10 11 3.87 10 11 3.90 10 第六年初 1 龄鱼的数量为: ,求得： 11 1 (6) 1.19 10 x = 11 1 (6) (1) 0.02 0.1 (1) xx x - = ， 故可以认为对渔群的生产能力破坏处于可接受范围， 即该模型不仅可以
23、保证持续捕捞又可实 现 5 年捕捞量最大。8 六 模型评价 本文所建模型的优点主要体现在：针对问题一、问题二，在合理的假设前提 下，分别建立了非线性规划模型，通过约束条件求解出最优解。运用 MATLAB 对其结果用图像直观的表现出来，且与实际吻合的比较好，可见本文所建立的模 型具有一定的合理性与参考价值。 但所建模型也存在一些不足： 为方便模型的建立， 假设鳀鱼的产卵时刻集中， 而在实际上鳀鱼的产卵在一定时间内时连续的；在问题二中，为了简化模型，假 定五年 3、4 龄鱼捕捞强度不变，而实际上捕捞强度与每年各年龄组鱼群条数成 正比，是个变化的量；并且，在为鱼群生产能力不能受到太大破坏给出指标时， 有一定的主观性； 这些都可能会导致模型所求的结果和实际数据之间存在一些偏 差。 参考文献 1 姜启源，谢金星，叶俊，数学模型（第三版）M.北京：高等教育出版社，2003。 2 谢金星，薛毅，优化建模与 LINDO/LINGO 软件M北京：清华大学出版社，2005。 9
