1实验8 1捕鱼业的持续收获产量模型实验7 2种群的相互竞争1

1、河北大学数学模型实验实验报告班级专业15计科2班姓名张宇轩学号20151101006实验地点C1-229指导老师司建辉成绩实验项目实验8-1捕鱼业的持续收获一一产量模型实验7-2种群的相互竞争(1 )实验7-3种群的相互竞争(2 )、实验目的学会利用MATLABS行实验，学会使用符号表达式，熟练掌握用函数solve求解代数方程 组，复习在MATLAB画图方法，合理使用函数text、axis、grid等。二、实验要求f(x)=rx(1-x/N) h(x)=Ex实验8-1捕鱼业的持续收获 一一产量模型 运行下面的m文件，并把相应结果填空，即填入“ clear;clc;%无捕捞条件下单位时间的增长量

2、： %甫捞条件下单位时间的捕捞量：%F(x)=f(x)-h(x)=rx(1-x/N)-Exx(t)=F(x)%甫捞情况下渔场鱼量满足的方程： %满足F(x)=0的点x为方程的平衡点%求方程的平衡点 syms r x N E; % 定义符号变量Fx=r*x*(1-x/N)-E*x; %创建符号表达式x=solve(Fx,x) % 求解 F(x)=0 (求根)%得到两个平衡点，记为：% x0=, x1=x0=x (2);x1=x(1);%符号变量x的结构类型成为V2X 1sym %求 F(x)的微分 F(x)syms x; %定义符号变量x的结构类型为 dF=diff(Fx,x);dF=simpl

3、e(dF) %简化符号表达式F(x)=%求 F(x0)并简化 dFx0=subs(dF,x,x0); % 将 x=x0 代入符号表达式 dF dFx0=sim ple(dFx0) %寻 F (x0)=%求 F (x1) dFx1=subs(dF,x,x1)F (x1)=%若 Er，有F(x0)0，故xO点稳定，x1点不稳定(根据平衡点稳定 性的准则)；x0 的前提下( Er ，则结果正好相反。 %在渔场鱼量稳定在 %通过分析(见教材 syms r x N fx=r*x*(1-x/N);df=diff(fx,x); x0=solve(df,x) %得 x0*=hm=subs(fx,x,x0) %

4、得 hm=，可得 E*=%又由 x0*=N(1-E/r) %产量模型的结论是：%将捕捞率控制在固有增长率的一半(E=r/2 )时，能够获得最大的持续产量。simple 的格式：S 尝试多种不同的算法简化,以显示 S 表达式的长度最短的简化形式。 sub 的格式： 提示 符号简化函数 simple(S) 对符号表达式 变量替换函数Subs(S,OLD,NEW) 将符号表达式S中的OLD变量替换为NEW变量。实验7-2 种群的相互竞争( 1) 补充如下指出的程序段，然后运行该 m 文件，对照教材上的相应结果。 clear;clc;%甲乙两个种群满足的增长方程： % x1(t)=f(x1,x2)=r

5、1*x1*(1-x1/N1-k1*x2/N2) % x2(t)=g(x1,x2)=r2*x2*(1-k2*x1/N1-x2/N2)%求方程的平衡点，即解代数方程组% f(x1,x2)=0 g(x1,x2)=0 编写出该程序段。 提示 ( 1 ) 使用符号表达式；( 2) 用函数 solve 求解代数方程组；P的结构类型为,( 3)调整解(平衡点)的顺序放入P 中(见下面注释所示)，P 的第 1 列对应%得 4 个平衡点：% P(1)=P1% P(2)=P2% P(3)=P3% P(4)=P4x1 ，第 2 列对应 x2。0)N2)N1 ，O，N1*(-1+k1)/(-1+k2*k1) ， N2

6、*(-1+k2)/(-1+k2*k1),或同时大于 1 。O， O)%平衡点位于第一象限才有意义，故要求P3： k1,k2 同时小于%判断平衡点的稳定性(参考教材p224)fx1=diff(f,x1);fx2=diff(f,x2);gx1=diff(g,x1);gx2=diff(g,x2);A=fx1,fx2;gx1,gx2syms x1 x2;P=subs(-(fx1+gx2),x1,x2, P( :,1) ,P (:,2); p=simple(p);%简化符号表达式p q=subs(det(A),x1,x2, P(:,1), P( :,2); q=si mp le(q);P P q经测算可

7、得该表的第4列，即稳定条件。%寻到教材p225表1的前3列，实验7-3种群的相互竞争(2 ) 求微分方程组X2(t尸2兀2(1的数值解，分别画出教材P 227 有关数据参见教材P227中“计算与验证”。提示(1 )求微分方程组的数值解可参考教材P140(2 )在 figure(1) 中画图 2(a)，在 figure(2)中的图2(a)、(b、 (c)。的程序。中画图2(b)，在figure)中画图2(c)。在程序中，figure(图形编号)用于定位对应图形。(3 )使用text(x,y,标识文本)，坐标点(x,y)在“标识文本”的左边，调整(x,y) 值，使“标识文本”放在图中的适当位置。a

8、xis(xmin xmax ymin ymax)控制坐标的刻度范围。grid on 打开网格，grid off 关闭网格。hold on 把要画的图形保持在之前在同一 figure上所画的图形中(同一坐(4 )用(5 )用(6 )用 标系)。2(c)中的两“点线”直线，一条的两个端点为(0,1)和(1,0)，另一条的两(7 )图 个端点为(0,2)和(1.6,0)。三、实验内容产量模型1.实验8-1捕鱼业的持续收获 syms r x N E;Fx=r*x*(1-x/N)-E*x; x=solve(Fx,x)x0=x(2); x1=x(1); syms x; dF=diff(Fx,x); dF=

9、simple(dF) dFx0=subs(dF,x,x0); dFx0=simple(dFx0) dFx1=subs(dF,x,x1) syms r x N fx=r*x*(1-x/N); df=diff(fx,x); x0=solve(df,x) hm=subs(fx,x,x0) 2实验 7-2 种群的相互竞争( 1 ) syms r1 r2 x1 x2 N1 N2 k1 k2; F=r1*x1*(1-x1/N1-k1*x2/N2); G=r2*x2*(1-k2*x1/N1-x2/N2); x1,x2= solve(F,G) P=x1(2,4,1,3),x2(2,4,1,2); fx1=di

10、ff(F,x1); fx2=diff(F,x2); gx1=diff(G,x1); gx2=diff(G,x2); A=fx1,fx2;gx1,gx2; syms x1 x2; p=subs(-(fx1+gx2),x1,x2,P(:,1),P(:,2); p=simple(p); q=subs(det(A),x1,x2,P(:,1),P(:,2); q=simple(q); P p q3.实验 7-3 种群的相互竞争( 2) M文件中：function y=fun(t,x) r=1;d=0.5;a=0.1;b=0.02; y=(r-a*x(2).*x(1);(-d+b*x(1).*x(2);

11、ts=0:0.1:15; x0=25,2; t,x=ode45(fun,ts,x0); t,x plot(t,x),grid,gtext(x(t),gtext(y(t),pause, plot(x(:,1),x(:,2),grid t,x=ode45(fun,ts,x0);四、实验结果及其分析1 .实验8-1捕鱼业的持续收获x =0-N*（-叶E）/r产量模型dF = r-2*r*x/N-EdFxO =-r+EdFx1 =r-ExO =1/2*NE=r/2分析：将捕捞率控制在固有增长率的一半（2.实验7-2种群的相互竞争（1）x1 =E=r/2 ）时，能够获得最大的持续产量。0N10N1*(-

12、1+k1)/(-1+k2*k1)x2 =00N2N2*(-1+k2)/(-1+k2*k1)A =r1*(1-x1/N1-k1*x2/N2)-r1*x1/N1,-r1*x1*k1/N2-r2*x2*k2/N1, r2*(1-k2*x1/N1-x2/N2)-r2*x2/N2Jins 0,HL-r-r2.fl*-k+iil+E2*lr2+r2,-rl+r*hl-E-i+kJ*k2), -Cl-k3-hI也I怨即*rrZFC-】+1tLM笳 r-*c23.实验7-3种群的相互竞争（2 ）ans =025.00002.00000.100027.08182.00410.200029.33442.01700

13、.300031.76892.03940.400034.39612.07260.500037.22582.11780.600040.26732.17670.700043.50122.25340.800046.93602.35030.900050.60722.46831.000054.53012.61061.100058.69992.78191.200063.09172.98911.300067.66043.24111.400072.34093.54841.500077.04793.92381.600081.67594.38191.700086.09964.93911.800090.17325.6

14、1401.900093.73116.42682.000096.58737.40002.100098.53608.55772.200099.30559.92342.300098.614311.50852.400096.285113.30672.500092.247215.28822.600086.585317.39472.700079.534919.54272.800071.536421.62252.900063.084823.53003.000054.623625.18193.100046.544126.51633.200039.186027.49213.300032.793228.09783

15、.400027.336828.37663.500022.737528.37643.600018.913428.14263.700015.777127.71783.800013.235427.14263.900011.187326.45564.00009.527825.69114.10008.175824.87404.20007.068424.02454.30006.159223.15804.40005.418522.28534.50004.812921.41624.60004.312420.55824.70003.896719.71644.80003.550018.89484.90003.26

16、0118.09615.00003.018117.32195.10002.815416.57355.20002.646015.85155.30002.504715.15625.40002.387214.48755.50002.289813.84545.60002.209913.22955.70002.145512.63935.80002.095312.07405.90002.057811.53306.00002.031711.01546.10002.016210.52076.20002.010610.04796.30002.01439.59646.40002.02709.16536.50002.

17、04838.75396.60002.07818.36146.70002.11647.98706.80002.16327.63006.90002.21887.28977.00002.28336.96547.10002.35736.65657.20002.44106.36227.30002.53516.08217.40002.64015.81557.50002.75665.56187.60002.88555.32047.70003.02765.09107.80003.18374.87297.90003.35494.66568.00003.54244.46888.10003.74784.28198.

18、20003.97234.10468.30004.21743.93658.40004.48473.77738.50004.77613.62658.60005.09373.48388.70005.43983.34908.80005.81713.22178.90006.22833.10179.00006.67662.98889.10007.16532.88269.20007.69782.78319.30008.27812.69009.40008.91002.60329.50009.59802.52259.600010.34682.44799.700011.16202.37949.800012.049

19、42.31669.900013.01532.259610.000014.06652.208410.100015.21022.162910.200016.45412.123310.300017.80632.089810.400019.27552.062710.500020.87082.042210.600022.60162.028710.700024.47792.022810.800026.51022.024910.900028.70932.035511.000031.08442.055811.100033.64162.087511.200036.39822.130811.300039.36802.186911.400042.55982.257311.500045.97762.344711.600049.62052.452511.700053.48292.584911.800057.55392.746811.900061.81802.944112.000066.25493.183412.100070.84283.469212.200075.52223.812912.300080.19004.234812.400084.72134.754512.500088.96865.391812.600092.76226.165812