数学建模实验_稳定性模型讲课.docVIP

实验08 稳定性模型（4学时） （第7章 稳定性模型） 1.（验证）捕鱼业的持续收获 ——产量模型p215~219 产量模型： 其中， x(t)为t时刻渔场中的鱼量。 r是固有增长率。 N是环境容许的最大鱼量。 E是捕捞强度，即单位时间捕捞率。 要求： 运行下面的m文件，并把相应结果填空，即填入“_________”。%7.1 捕鱼业的持续收获——产量模型 %文件名：p.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)=