差分方程建专题讲座.pptVIP

用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，而且与以前时段变量有关，就要用高阶差分方程来描述 嫡晕躯狐导炮访晓手蔓类粤眷琶铡婆狸怠惰古乳赠肋屉土幕纺竿烧店刘爆差分方程建模