[数学]差分方程模型.ppt

一、差分的概念及性质 二、差分方程的基本概念 三、线性差分方程解的结构 1、一阶常系数齐次线性差分方程的求解 2、一阶常系数非齐次线性差分方程的求解 差分方程是在离散时段上描述现实世界中变化过程的数学模型 例1、 某种货币1年期存款的年利率是r ，现存入M元，问n年后的本金与利息之和是多少？ Xk+1=(1+r)Xk , k = 0 , 1 , 2 ····· 以k=0时x0=M代入，递推n次可得n年后本息为 一阶线性常系数差分方程 模型建立 记第k年沙丘鹤的数量为Xk，年均增长率为r，则第k+1年鹤的数量为 Xk+1=(1+r)Xk k=0,1,2······ 已知X0=100, 在较好，中等和较差的自然环境下 r=0.0194, -0.0324,和-0.0382 , 利用Matlab编程，递推20年后观察沙丘鹤的数量变化情况 Matlab实现 首先建立一个关于变量n ，r的函数 function x=sqh(n,r) a=1+r; x=100; for k=1:n x(k+1)=a*x(k); end 在command窗口里调用sqh函数 利用plot 绘图观察数量变化趋势 可以用不同线型和颜色绘图 r g b c m y k w 分别表示 红 绿 兰 兰绿 洋红 黄 黑 白色 ： + o * . X s d 表示不同的线型 plot(k,y1,k,y2,k,y3) 在同一坐标系下画图 人工孵化是挽救濒危物种的措施之一，如果每年孵化5只鹤放入保护区，观察在中等自然条件下沙丘鹤的数量如何变化 Xk+1=aXk +5，a=1+r 如果想考察每年孵化多少只比较合适，可以令 Xk+1=aXk +b，a=1+r Matlab实现 function x=fhsqh(n,r,b) a=1+r; x=100; for k=1:n x(k+1)=a*x(k)+b; end k=(0:20); %一个行向量 y1=(20,-0.0324,5); 也是一个行向量 round( [ k ’, y 1 ’] ) 对k,y1四舍五入， 但是 不改变变量的值 plot(k, y1) k y1 是行向量列向量都可以 也可以观察200年的发展趋势，以及在较差条件下的发展趋势，也可以考察每年孵化数量变化的影响。 一阶线性常系数差分方程的解、平衡点及其稳定性 自然环境下,b=0 人工孵化条件下 令xk=xk+1=x得 差分方程的平衡点 k→∞时，xk→x，称平衡点是稳定的 高阶线性常系数差分方程 如果第k+1时段变量Xk+1不仅取决于第k时段变量Xk，而且与以前时段变量有关，就要用高阶差分方程来描述 一年生植物的繁殖 一年生植物春季发芽，夏天开花，秋季产种，没有腐烂，风干，被人为掠取的那些种子可以活过冬天，其中一部分能在第2年春季发芽，然后开花，产种，其中的另一部分虽未能发芽，但如又能活过一个冬天，则其中一部分可在第三年春季发芽，然后开花，产种，如此继续，一年生植物只能活1年，近似认为，种子最多可以活过两个冬天，试建立数学模型研究这种植物数量变化的规律，及它能一直繁殖下去的条件。 模型及其求解 记一棵植物秋季产种的平均数为c，种子能活过一个冬天的(1岁种子)比例为b，活过一个冬天没有发芽又活过一个冬天的（2岁种子）比例仍为b，1岁种子发芽率a1，2岁种子发芽率a2。 设c,a1,a2固定，b是变量，考察能一直繁殖的条件 记第k年植物数量为Xk，显然Xk与Xk-1,Xk-2有关，由Xk-1决定的部分是 a1bcXk-1,由Xk-2决定的部分是 a2b(1-a1)bcXk-2 Xk= a1bcXk-1 + a2b(1-a1)bcXk-2 实际上，就是Xk= pXk-1 + qXk-2 我们需要知道X0,a1,a2,c, 考察b不同时，种子繁殖的情况。在这里假设 X0=100,a1=0.5,a2=0.25,c=10,b=0.18~0.20 用matlab计算 function x=fz(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*X(1); for k=3:n X(k)=p*X(k-1)+q*X(k-2); end K=(0:20)’; Y1=fz(100,21,0.18); Y2=fz(100,21,