第九讲 差分方程 Matlab语言程序的设计 .pptVIP

三、线性常系数差分方程组 例3：汽车租赁公司的运营 作图观察数量变化趋势 k=0:10; plot(k,x) ，grid gtext(x1(k)), gtext(x2(k)), gtext(x3(k)) exno97 三、线性常系数差分方程组 例3：汽车租赁公司的运营 可以看到时间充分长以后3个城市汽车数量趋于180，300，120 可以考察这个结果与初始条件是否有关 若最开始600辆汽车都在A市，可以看到变化时间充分长以后，各城市汽车数量趋于稳定，与初始值无关 东北电力大学 理学院 徐 屹 Northeast Dianli College Matlab语言程序设计 第九讲 差分方程 一阶线性常系数差分方程 高阶线性常系数差分方程 2 线性常系数差分方程组 3 3 非线性差分方程 4 4 1 差分方程是在离散时段上描述现实世界中变化过程的数学模型 例1 某种货币1年期存款的年利率是r ，现存入M元，问年后的本金与利息之和是多少？ Xk+1=(1+r)xk , k = 0 , 1 , 2 ····· 以k=0时x0=M代入，递推n次可得n年后本息为 第九讲 差分方程 污水处理厂每天可将处理池的污水浓度降低一个固定比例q，问多长时间才能将污水浓度降低一半？ 记第k天的污水浓度为ck,则第k+1天的污水浓度为 ck+1=(1-q)ck,k=0,1,2,···· 从k=0开始递推n次得 以cn=c0/2代入即求解。 第九讲 差分方程 一、一阶线性常系数差分方程 濒危物种的自然演变和人工孵化 问题 Florida沙丘鹤属于濒危物种，它在较好自然环境下，年均增长率仅为1.94%，而在中等和较差环境下年均增长率分别为 -3.24% 和 -3.82%，如果在某自然保护区内开始有100只鹤，建立描述其数量变化规律的模型，并作 数值计算。 例1：濒危物种的自然演化和人工孵化 一、一阶线性常系数差分方程 例1：濒危物种的自然演化和人工孵化 在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]) exno91 一、一阶线性常系数差分方程 例1：濒危物种的自然演化和人工孵化 利用plot 绘图观察数量变化趋势 可以用不同线型和颜色绘图 r g b c m y k w 红 绿 兰 兰绿 洋红 黄 黑 白色 ： + o * . X s d 表示不同的线型 一、一阶线性常系数差分方程 例1：濒危物种的自然演化和人工孵化 plot(k,y1,k,y2,k,y3) 在同一坐标系下画图 subplot(2,2,1),plot(k,y1,r:) subplot(2,2,2),plot(k,y2,k-) subplot(2,2,3),plot(k,y3,g*) subplot(2,2,4),plot(k,y1,r:,k,y2,k-,k,y3,g*) 用gtext(‘r=0.0194’),gtext(‘r=-0.0324’),gtext(‘r=-0.0382’)在图上做标记。 exno92 一、一阶线性常系数差分方程 例1：濒危物种的自然演化和人工孵化 人工孵化是挽救濒危物种的措施之一，如果每年孵化5只鹤放入保护区，观察在中等自然条件下沙丘鹤的数量如何变化 Xk+1=aXk+5 ,a=1+r 如果我们想考察每年孵化多少只比较合适，可以令 Xk+1=aXk+b ,a=1+r 一、一阶线性常系数差分方程 例1：濒危物种的自然演化和人工孵化 function x=fhsqh(n,r,b) a=1+r; x=100; for k=1:n x(k+1)=a*x(k)+b; end fhsqh 一、一阶线性常系数差分方程 例1：濒危物种的自然演化和人工孵化 k=(0:20); %一个行向量 y1=fhsqh(20,-0.0324,5); % 也是一个行向量 round([ k’,y1’])% 对k,y1四舍五入，但是不改变变量的值 plot(k,y1) %k,y1是行向量列向量都可以 也可以观察200年的发展趋势，以及在较差条件下的发展趋势，也可以考察每年孵化数量变化的影响。 exno93,exno94 一、一阶线性常系数差分方程 例1：濒危物种的自然演化和人工孵化 一阶线性常系数差分方程的解、平衡点及其稳定性 自然环境下,b=0 人工孵化条件下 令xk=xk+1=x得