差分方程建模专题讲座.pptVIP

Matlab实现首先建立一个关于变量n,r的函数functionx=sqh(n,r)a=1+r;x=100;fork=1:nx(k+1)=a*x(k);end第22页，共56页，星期日，2025年，2月5日在command窗口里调用sqh函数k=(0:20);y1=sqh(20,0.0194);y2=sqh(20,-0.0324);y3=sqh(20,-0.0382);round([k,y1,y2,y3])%一个行向量%也是一个行向量%对变量四舍五入，但是不改变变量的值第23页，共56页，星期日，2025年，2月5日利用plot绘图观察数量变化趋势可以用不同线型和颜色绘图rgbcmykw分别表示红绿兰兰绿洋红黄黑白色：+o*.Xsd表示不同的线型第24页，共56页，星期日，2025年，2月5日plot(k,y1,k,y2,k,y3)在同一坐标系下画图plot(k,y1,r,k,y2,y,k,y3,--)gtext(r=0.0194)gtext(r=-0.0324)gtext(r=-0.0382)在图上做标记。第25页，共56页，星期日，2025年，2月5日人工孵化是挽救濒危物种的措施之一，如果每年孵化5只鹤放入保护区，观察在中等自然条件下沙丘鹤的数量如何变化Xk+1=aXk+5,a=1+r如果我们想考察每年孵化多少只比较合适，可以令Xk+1=aXk+b,a=1+r第26页，共56页，星期日，2025年，2月5日观察200年的发展趋势，以及在较差条件下的发展趋势可以考察每年孵化数量变化的影响。第27页，共56页，星期日，2025年，2月5日functionx=fhsqh(n,r,b)a=1+r;x=100;Fork=1:nx(k+1)=a*x(k)+b;end第28页，共56页，星期日，2025年，2月5日k=(0:200);n=200;r=-0.0324b=5y=fhsqh(n,r,b)plot(k,y)第29页，共56页，星期日，2025年，2月5日一阶线性常系数差分方程的解、平衡点及其稳定性

自然环境下,b=0人工孵化条件下令xk=xk+1=x得差分方程的平衡点k→∞时，xk→x,称平衡点是稳定的第30页，共56页，星期日，2025年，2月5日高阶线性常系数差分方程如果第k+1时段变量xk+1不仅取决于第k时段变量xk，而且与以前时段变量有关，就要用高阶差分方程来描述第31页，共56页，星期日，2025年，2月5日一年生植物的繁殖一年生植物春季发芽，夏天开花，秋季产种，没有腐烂，风干，被人为掠取的那些种子可以活过冬天，其中一部分能在第2年春季发芽，然后开花，产种，其中的另一部分虽未能发芽，但如又能活过一个冬天，则其中一部分可在第三年春季发芽，然后开花，产种，如此继续，一年生植物只能活1年，而近似的认为，种子最多可以活过两个冬天，试建立数学模型研究这种植物数量变化的规律，及它能一直繁殖下去的条件。第32页，共56页，星期日，2025年，2月5日模型及其求解记一棵植物春季产种的平均数为c,种子能活过一个冬天的(1岁种子)比例为b,活过一个冬天没有发芽又活过一个冬天的（2岁种子）比例仍为b,1岁种子发芽率a1，2岁种子发芽率a2。设c,a1,a2固定，b是变量，考察能一直繁殖的条件记第k年植物数量为xk，显然xk与xk-1,xk-2有关，由xk-1决定的部分是xk-1Cba1,由xk-2决定的部分是xk-2cb(1-a1)ba2xk=a1bcxk-1+a2b(1-a1)bcxk-2第33页，共56页，星期日，2025年，2月5日实际上，就是xk=pxk-1+qxk-2,我们需要知道x0,a1,a2,c,考察b不同时，种子繁殖的情况。在这里假设x0=100,a1=0.5,a2=0.25,c=10,b=0.18---0.2这样可以用matlab计算了xk=a1bcxk-1+a2b(1-a1)bcxk-2第34页，共56页，星期日，2025年，2月5日functionx=zwfz(x0,n,b)c=10;a1=0.5;a2=0.25;p=a1*b*c;q=a2*b*(1-a1)*b*c;X(1)=x0;X(2)=p*(x1);fo