数学建模实验案_稳定性模型.doc

实验08 稳定性模型（4学时） （第7章 稳定性模型） 1.（验证）捕鱼业的持续收获 ——产量模型p215~219 产量模型： 其中， x(t)为t时刻渔场中的鱼量。 r是固有增长率。 N是环境容许的最大鱼量。 E是捕捞强度，即单位时间捕捞率。 要求： 运行下面的m文件，并把相应结果填空，即填入“_________”。 %7.1 捕鱼业的持续收获——产量模型 %文件名：p215_217.m clear; clc; %无捕捞条件下单位时间的增长量：f(x)=rx(1-x/N) %捕捞条件下单位时间的捕捞量：h(x)=Ex %F(x)=f(x)-h(x)=rx(1-x/N)-Ex %捕捞情况下渔场鱼量满足的方程：x(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= ，见P216(4)式 x0=x(2); x1=x(1);%符号变量x的结构类型成为2×1sym %求F(x)的微分F(x) syms x; %定义符号变量x的结构类型为1×1sym dF=diff(Fx,x); %求导 dF=simple(dF) %简化符号表达式 %得 F(x)= %求F(x0)并简化 dFx0=subs(dF,x,x0); %将x=x0代入符号表达式dF dFx0=simple(dFx0) %得 F (x0)= %求F (x1) dFx1=subs(dF,x,x1) %得 F (x1)= %若 Er，有F(x0)0，F(x1)0，故x0点稳定，x1点不稳定（根据平衡点稳定性的准则）； %若 Er，则结果正好相反。 %在渔场鱼量稳定在x0的前提下（Er），求E使持续产量h(x0)达到最大hm。 %通过分析（见教材p216图1），只需求x0*使f(x)达到最大，且hm=f(x0*)。 syms r x N fx=r*x*(1-x/N); %fx=dF df=diff(fx,x); x0=solve(df,x) %得 x0*= ，见P217(6)式 hm=subs(fx,x,x0) %得 hm= ，见P217(7)式 %又由 x0*=N(1-E/r)，可得 E*= ，见P217(8)式 %产量模型的结论是： %将捕捞率控制在固有增长率的一半（E=r/2）时，能够获得最大的持续产量。 符号简化函数simple，变量替换函数sub的用法见提示。 ★ 给出填空后的M文件（见[ HYPERLINK \l _[215]7.1 捕鱼业的持续收获 215~217]）： %7.1 捕鱼业的持续收获——产量模型 %文件名：p215_217.m clear; clc; %无捕捞条件下单位时间的增长量：f(x)=rx(1-x/N) %捕捞条件下单位时间的捕捞量：h(x)=Ex %F(x)=f(x)-h(x)=rx(1-x/N)-Ex %捕捞情况下渔场鱼量满足的方程：x(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= -N*(-r+E)/r , x1= 0 ，见P216(4)式 x0=x(2); x1=x(1);%符号变量x的结构类型成为2×1sym %求F(x)的微分F(x) syms x; %定义符号变量x的结构类型为1×1sym dF=diff(Fx,x); %求导 dF=simple(dF) %简化符号表达式 %得 F(x)= r-2*r*x/N-E %求F(x0)并简化 dFx0=subs(dF,x,x0); %将x=x0代入符号表达式dF dFx0=simple(dFx0) %得 F (x0)= -r+E %求F (x1) dFx1=subs(dF,x,x1) %得 F (x1)= r-E %若 Er，有F(x0)0，F(x1)0，故x0点稳定，x1点不稳定（根据平衡点稳定性的准则）； %若 Er，则结果正好相反。 %在渔场鱼量稳定在x0的前提下（Er），求E使持续产量h(x0)达到最大hm。 %通过分析（见教材p216图1），只需求x0*使f(x)达到最大，且hm=f(x0*)。 syms r x N fx=r*x*(1-x/N); %fx=dF df=diff(fx,x); x0=solve(df,x)