最优捕鱼策略实验报告
学号:104080298 姓名:宁亚会班级:10D
摘要
为了保护人类赖以生存的自然环境,可再生资源(如渔业、林业资源等)的开发必须适
度。而在社会经济生活中,我们要使商业活动在一段时期内达到最大收益,因此我们要合理的
开发资源,这时,我们不仅要考虑商业活动的当前经济效益,还要考虑生态效益及由此产生的
对整体经济效益的影响。本文就是对渔业这类可再生资源的开发问题进行研究,利用相关的
数学软件进行求解。
对于问题一,我们考虑渔场生产过程中的各年龄组鱼群数量的制约因素,将其分为两大类,第1,2龄鱼群为一类,该鱼群数量变化在一年内只受自然死亡率制约,写出鱼群数量满
足的微分方程;第3,4龄鱼群为一类,其数量变化在前8个月受捕捞强度和自然死亡率影响,后4个月只受自然死亡率的制约,分阶段写出写出鱼群数量满足的微分方程;根据微分方程,
求出在某时刻各鱼群的数量表达式(类似于人口增长模型)。因为捕捞是连续的,所以任意
一个时刻的捕捞量为捕捞强度乘以鱼群的数量,又捕捞只在前8个月进行,则年捕捞量为前
8个月各时刻鱼群数量的积分。最后建立年总捕捞量的函数与生产过程中满足的关系式,转化为非线性规划模型,利用lingo和matlab软件分别求解。
对于问题二,题中已给出各年龄组鱼群的初始值,我们利用问题一中所得到的迭代方程,
可迭代地求出第i年初各年龄组鱼群的数量;再根据问题一中的捕捞量表达式,可写出5年的捕捞总量表达式,以5年捕捞总量最大为前提,利用matlab软件求解出此时的捕捞强
度,然后再验证在此捕捞强度下会不会使5年后鱼群的生产能力有太大的破坏。
最后,我们得出以下结论:可持续捕获条件下,捕捞强度为17.36292时,达到最大捕
捞总质量3.887076 1011g ;5年后鱼群的生产能力不会有太大的破坏条件下,捕捞强度为
k 17.5,17.8,达到最大最大捕捞总质量1.6056 1012
一.问题重述
生态学表明,对可再生资源的开发策略应在事先可持续收获的前提下追求最大经济效益,考虑具有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龄鱼产卵量为1.109 105(个),3龄鱼为其一半;卵孵化的成活率为1.22 101;(1.22 1011 n)(n为产卵总量);
二、问题分析
对于问题一,要实现可持续捕获,即每年开始捕捞时渔场中各年龄组鱼群条数不变,因此我们要算出每年初各龄鱼组的数量。
对自然死亡率,捕捞强度系数和卵成活率的解释,即1 , 2龄鱼仅受自然死亡率的影响;
而3, 4龄鱼不仅受自然死亡率的影响,还受捕捞强度系数的影响;因为该种鱼的最高寿命为4,所以在后四月中4龄鱼都不存活;,而对于1龄鱼的数量,是3,4龄鱼在前年的后4 年产卵所存活下来的数量;对于捕捞量,题中规定只在1到8月才能捕捞,而且1, 2龄鱼
不被捕捞,所以主要来源于对3, 4龄鱼的捕捞。根据这些关系可列出一系列的方程,其中
捕捞量作为目标函数,其他的作为约束条件,建立一个非线性规划模型,再然后用lin go 软件和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龄鱼的生长只受自然死亡率的影响, 程
可得:X i(t) X o e",X o为每年年初
T年的i龄鱼在T+1年变为i+1龄鱼, 由此可知1, 2龄鱼的生长的微分方程满足方
dxi
dt
rxi,i 1,2,3,4
则有:
r x「(T 1) e X i (T)
而对于3, 4龄鱼的生长,在前八个月,他们的生长不仅受自然生长率的影响,还受捕捞强度系数的影响,而后四个月仅受自然生长率的影响。
我们以一年为一个时间单位,则这一时间单位可以分为两个阶段,见图( 1):
2/3 1
I :捕捞期
图(1)
因此,
①.前八个月3、4龄鱼生长的微分方程满足:
捞量为k/dt),则年捕捞量为:
约束条件:
X2ne
X32r ne
(
2
*0.42k 3r)/ (-k
X4n e 3(1 e 3
X
1
r)
)
r
dx i
dt
(r kJX j i 3,4
可得:X i(t) x°e (r k i)t, X0为每年年初i龄鱼的数量
由于每年的捕捞只在1到8月进行,并且只能捕到3, 4龄鱼,所以任意一个时刻的捕
2
]k i X i(t)dt
宀x°(1 e韓
k i r
r)
)
②.后四个月3、4龄鱼生长的微分方程满足方程
dxi
rxi, i
dt
1,234
可得:x i (t) x0e rt ,x0为每年年初
2
(r 0.42k)
其中产卵量n 0.5mx3e 3
-(r
3
m%e 3
k)
又孵化存活率
11
1.22*10
1.22*1011
(8)
所以年初1龄鱼的总量X1(T 1) n(T) (9)
根据以上分析,我们可以建立非线性规划模型:
2
17.86 Jk3x3(t)dt
目标函数:max
2
22.99 jk4X4(t)dt (10)
0.5mx3e
2
(0.42k r)
3
2
-(
k
r)
m%e 3
II :产卵期
4.2模型求解
可将目标函数和约束条件转化为: 目标函数为:
约束条件:
1.22*10^1 1.22*仍1 n
2
8 8
64
(才 0.28k p)
(0.28k 一)
e 3 3
n 1.22*10A11*(m* (0.5*e
3 2
) 1)
(_k 0.8)
1 e 3
下面利用 LINGO 软件求解,输入以下命令:
max =17.86*0.42*k/(0.8+0.42*k)*1.22*10A 11/(1.22*10A
11+ n)* n* @exp (-1.6)*
(1- @exp (-2/3*(0.8+0.42*k)))+22.99*k/(0.8+k)*1.22*10A11/(1.22*10A11+ n) *n* @exp(-0.28*k-2.4)/(1-
@exp (-2/3*k-0.8))*(1-
@exp(-2/3*(0.8+k)));
n=1.22*10人11*(1.109*10人5*(0.5* @exp (-0.28*k-6.4/3)+
@exp (-(0.28+2/3)*k-
8.8/3)/(1- @exp (-2/3*k-0.8)))-1); 直接运行,输出结果为:
Local optimal soluti on found at iterati on: 92
Objective value:
0.3887076E+12
Variable Value Reduced Cost
K 17.36292 -1.034723
N 0.6078067E+13
0.000000
Row Slack or Surplus Dual Price 1 0.3887076E+12 1.000000 2
0.000000
0.1258406E-02
k 17.36292, f 0.3887076E+12
下面再利用matlab 软件求解:
首先,建立 fun 1.m 、max1.m> picture.m 文件 picture1.m
k=li nspace(1,20,20);
n=1.22*10A11*(1.109*10A5*(0.5*exp(-0.28*k-6.4/3)+exp(-(0.28+2/3)*k-8.8/3)/(1-exp(-2/3 plot (k,n)
fun 1.m 文件:
max 17.86
0.42k
0.42k 0.8 ne '"(I
2(0.8 0.42 k)
e 3
) 22.99
0.8 k .3(0.8 k)
e 3 )
(0.28 k
ne
(Ik 0.8)
e 3 )(1
fun cti on n=fun 1(k) n=1.22*10八11*(1.109*10八5*(0.5*exp(-0.28*k-6.4/3)+exp(-(0.28+2/3)*k-8.8/3”(1-exp(-2/3
maxl.m 文件:
fun cti on y=max1( n,k) y=17.86*0.42*k/(0.8+0.42*k)*1.22*10A11/(1.22*10A11+ n)*
n*exp(-1.6)*(1-exp(-2/3*(0.8+( 99*k/(0.8+k)*1.22*10A11/(1.22*10A11+ n)* n*exp(-0.28*k-2.4)/(1-exp(-2/3*k-0.8))*(1-exp(-
我们对k 1,20 ,
运行picture.m 文件,画出n关于k的图像:
由图像,,我们看出k Z , n],这与事实是相符的。
然后,在命令窗口中输入以下命令:
程序1.m
for k=1:1:20
n=fu n1(k);
f=max1( n,k);
end
运行出结果我们观察到k 1,17 ,f Z ,k 17,20 ,f ] ,因此,改变k 的取值区间即步长,建立了以下文件,从而求得更精确的结果。
程序2.m
for k=17:0.1:18
n=fun1(k);
f=max1(n,k);
k
f
end
从运行结果看:
k =17.3000 或k =17.4000 ,取最大值f =3.8871e+011
最后,建立主程序3.m k 17.3,17.4
程序3.m
for k=17.3:0.01:17.4
n=fun1(k);
f=max1(n,k);
k
f
end
从运行结果看:
k 17.3,17.4 ,f =3.8871e+011
问题一所求结果为:k 17.3,17.4时,f取最大值3.8871e 011
4.3 结果分析
3龄鱼的捕捞强度为7.29/ 年;4龄鱼的捕捞强度为17.36/ 年:最优可持续捕捞量
3.8871* e11,可持续捕捞的鱼群大小(条数) :
1 龄 1.1960 1011
2 龄 5.3740 1010
10
7
3 龄 2.4147 10
4 龄 8.3967 10
分析结果发现,4龄鱼在年末存活的数量占全部数量的比例相对很小。
五.问题二的模型建立与求解
5.1模型建立
针对渔业公司的 始值的前提下,可迭代求出各龄鱼群第
n(t 1) %(t 1)
整个生存过程满足的关系式为
x 2(t 1)
X 3(t 1) X 4(t 1) k 3 0.42k 4
这样max 就变为关于k 的函数,易于求解。
用matlab 软件求解 1、利用 matlab 软件,建立fun.m 文件 fun.m 文件:
function y=fun(k3,k,m,l1,l2,l3)
5年捕捞计划,我们利用已得到的迭代方程在已知各个年龄组的鱼的初
i 年的鱼量的分布的函数。
(k 3 r)2/3 (k 3 r)2/3
m/2x 3(t)e
m&(t)e
1.22 1011 n(t 1) 1.22 1011 n(t 1) e /(t) 目标函数为
e 乂⑴
x 3(t)e (k 3
2/3 r)
他2/3 r)
x(t)e
max 5
(X 3(t))17.86
t 1
(k j r)2/3
5
(x 4(t) 22.99
t 1
k 4 k 4
r (1
e (k 4
r)2/3
(15)
5.2模型求解 将目标函数转化为: max
0.42k 0.42k r
-r 0.42 k
e 3
(10.1 29.7e
2r
122
e
1.22 1011 n 2 2r n e
1.22 10 n 2 *
1.22 1011 1.22 1011
n 3 2r nr e
)
k
2 r k
r
2
0.42k r £k
2
2r
0.42k 2r 3 0.42k 1 k
1
e
3
(3.29
10.1e
3
3.29e
3
29.7e
3
10.1e 3 3
k r
2r Ik
3r
2
0.42 k
3r
2 2 0.42 k k
3r 1 4 0.42k k
3r 2k
3%
3
122
e 3
29.7e
3
3
10.1e 3
3
3.29e
1.22 1011 n 2
3r - 0.42k
4r -
2
0.42 k k
4r - i
4
0.42 k k
3r - 0.42k 2k
3
122e
3
3
29.7e
3
3
10.1e
3
1.22 1011 n 2 e
2
(0.42k r)
其中 n(2) 10.2*m/2*e 3
2
-(k r)
3.29* m* e 3
(0.28k 5
r)
n(3) 29.7* m/2* e 3
(0.28 k
10.1* m* e
3.29* m* e
3.29e
3r 8k
3
))