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

例如，如认为第一季度的销售量大体按线性增长，可设销售量 由 x=[[1:5],ones(5,1)];y=[11 12 13 15 16];z=x\y 求得 ， . 预测第六年起第一季度的销售量为 由于数据少，用回归分析效果不一定好. 如认为销售量并非逐年等量增长而是按前一年或前几年同期销售量的一定比例增长的，则可建立相应的差分方程模型.仍以第一季度为例，为简单起见不再引入上标， 最小二乘法求一组总体吻合较好的数据 编写Matlab程序如下： y0=[11 12 13 15 16]; y=y0(3:5);x=[y0(2:4),y0(1:3),ones(3,1)]; z=x\y 求得 ， 虽然这一差分方程恰好使所有统计数据吻合，但这只是一个巧合 不难看出，如分别对每一季度建立一差分方程，则根据统计数据拟合出的系数可能会相差甚大，但对同一种商品，这种差异应当是微小的，故应根据统计数据建立一个共用于各个季度的差分方程.为此，将季度编号为 编写Matlab程序如下： y0=[11 16 25 12 12 18 26 14 13 20 27 15 15 24 30 15 16 25 32 17]; y=y0(9:20); x=[y0(5:16),y0(1:12),ones(12,1)]; z=x\y 求得 ， ， ， 还是较为可信的 全国赛题选讲 3 （1） 1996年 A题 最优捕鱼策略 （2） 2009年 A题 制动器试验台的控制方法 用Matlab求解差分方程问题 2 一、一阶线性常系数差分方程 二、高阶线性常系数差分方程 三、线性常系数差分方程组 一、一阶线性常系数差分方程 濒危物种的自然演变和人工孵化 问题 Florida沙丘鹤属于濒危物种，它在较好自然环境下，年均增长率仅为1.94%，而在中等和较差环境下年均增长率分别为 -3.24% 和 -3.82%，如果在某自然保护区内开始有100只鹤，建立描述其数量变化规律的模型，并作 数值计算。 模型建立 记第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函数 k=(0:20); y1=sqh(20,0.0194); y2=sqh(20,-0.0324); y3=sqh(20,-0.0382); round([k,y1,y2,y3]) %一个行向量 %也是一个行向量 %对变量四舍五入，但 是不改变变量的值 利用plot 绘图观察数量变化趋势 可以用不同线型和颜色绘图 r g b c m y k w 分别表示 红绿兰兰绿洋红黄黑白色 ： + o * . X s d 表示不同的线型 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) 在图上做标记。 人工孵化是挽救濒危物种的措施之一，如果每年孵化5只鹤放入保护区，观察在中等自然条件下沙丘鹤的数量如何变化 Xk+1=aXk +5 ,a=1+r 如果我们想考察每年孵化多少只比较合适，可以令 Xk+1=aXk +b ,a=1+r 观察200年的发展趋势，以及在较差条件下的发展趋势 可以考察每年孵化数量变化的影响。 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:200); n=200; r=-0.0324 b=5 y=fhsqh(n,r,b) plot(k,y) 一阶线性常系数差分方程的解、平衡点及其稳定性 自然环境下,b=0 人工孵化条件下 令xk=xk+1=x得 差分方程的平衡点 k→∞时，xk→x,称平衡点是稳定的 高阶线性常系数差分方程 如果第k+1时段变量xk+1不仅取决于第k时段变量xk，而且与以前时段变量有关，就要用高阶差分方程来描述 一年生植物的繁殖 一年生植物春季发芽，夏天开花，秋季产种，没有腐烂，风干，被人