原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
实验 08 稳定性模型（ 4 学时）
（第7 章 稳定性模型）
1. （验证）捕鱼业的持续收获 ——产量模型 p215~219
产量模型：
x
x(t ) F ( x) rx 1 Ex
N
其中，
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 为方程的平衡点
% 求方程的平衡点
1
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)=
2
% 若 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 文件（见 [215~217] ）：
%7.1 捕鱼业的持续收获
文档评论（0）