数学建模的课件.ppt

数学建模实验（八） 微分方程 例1 dx/dt=rx(1-x/k)，x(0)=0.01 1.符号解 S=dsolve(‘Dx=r*x*(1-x/k)’,’x(0)=0.01’) S =k/(1+100*exp(-r*t)*k-exp(-r*t)) %转变为函数值 r=0.3;k=8; s=subs(S) t=0:40; ss=subs(s,’t’,t); plot(t,ss),grid 常微分方程求解 2. 数值解 fun=@(t,x)0.3.*x.*(1-x/8);%定义赋值函数 ode45(fun,[0,10],0.01)%直接解方程作图 [tout,xout]=ode45(fun,[0,10],0.1) )%方程数值解。 输入参数说明： fun:，被积分的函数名,如果函数用文件 fun.m定义，则用@fun或‘fun’调用。 function y=fun(t, x) y=0.3*x*(1-x/8); [0, 40]=[t的初始值, t的终值] 。 0.01：x的初始值。 输出参数说明： tout: 离散的自变量值， xout: 离散的函数值。 例1 研究种群增长的Logistic微分方程模型dn/dt= r (1-n) n 的动力学行为。 数值试验（1）取r=0.8 对不同的初值N0=[0.01 0.2 0.5 0.8] 观察解的变化 。 建立函数文件: funlog.m function y=funlog(t,x) y=0.8*x*(1-x); 运行程序： xb=[0.01 0.2 0.5 0.8]; for i=1:4 x0=xb(i); [tt,xx]=ode45(funlog,[0,10],x0); plot(tt,xx),grid, hold on end 被积函数 积分区间 初值 数值实验（2） 对固定的初值N0=0.1，取不同的r=[0.1, 0.8, 1, 2, 5] 观察解的变化 。 建立函数文件: funlog.m function y=funlog(t,x) global r y=r*x*(1-x); 运行程序： rb=[0.1, 0.8, 1, 2, 5]; global r for i=1:5 r=rb(i); [tt,xx]=ode45(funlog,[0,10],0.1); plot(tt,xx),grid, hold on end 蓝鲸的内禀增长率每年估计为5%，估计蓝鲸的最大环境载为150,000条。磷虾是蓝鲸喜欢的一种食物。磷虾的最大饱和种群为 500吨/英亩。当缺少捕食者，环境不拥挤时，磷虾种群以每年25%的速率增长。磷虾500吨/英亩可以提高蓝鲸2%的年增长率， 同时150,000条蓝鲸将减少磷虾10%的年增长率。确定鲸鱼与磷虾是否可以长期共存。 (a) 组建一个蓝鲸和磷虾的动态模型。模拟两个种群随时间的变化。假设初始状态为蓝鲸 5,000条、磷虾750吨/英亩。 习题1 (b) 假设捕捞使得鲸鱼只剩下它的平衡态的5%。而磷虾保持平衡态的数量。描述一旦停止捕捞将发生什么情况。鲸鱼恢复需要多长时间？磷虾群体将发生什么变化？ (c)确定鲸鱼与磷虾是否可以长期共存。 (d)检验(b)的结论对鲸鱼剩余量5%的灵敏度，给出鲸鱼种群恢复时间对它所受伤害程度的依赖关系。 解（a）设蓝鲸x1千条，磷虾x2吨/英亩 模型： dx1/dt=0.05 x1 (1-x1/150)+0.02/500 x1 x2 dx2/dt=0.25x2 (1-x2/500)-0.1/150x1 x2 x1 (0)=5, x2(0)=750 Matlab指令 function y=fun(t,x) y=[0.05*x(1)*(1-x(1)/150)+0.02/500*x(1)*x(2); 0.25*x(2)*(1-x(2)/500)-0.1/150*x(1)*x(2)]; end [tt,xx]=ode45(fun,[0,150],[5 750]); plot(tt,xx),grid 解（b）令 0=0.05 x1 (1-x1/150)+0.02/500 x1 x2 0=0.25x2 (1-x2/500)-0.1/150x1 x2 得平衡态 （0，0） x1* =181.0345 x2 * = 258.6207 dx1/dt=0.05 x1 (1-x1/150)+0.02/500 x1 x2 dx2/dt=0.25x2 (1-x2/500)-0.1/150x1 x2 x1 (0)=0.05*181, x2(0)=258.62 鲸鱼群体恢复需要9