差分方程解法_第1页
差分方程解法_第2页
差分方程解法_第3页
差分方程解法_第4页
差分方程解法_第5页
已阅读5页,还剩5页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

差分方程解法第1页,课件共10页,创作于2023年2月一阶线性常系数差分方程濒危物种的自然演变和人工孵化问题Florida沙丘鹤属于濒危物种,它在较好自然环境下,年均增长率仅为1.94%,而在中等和较差环境下年均增长率分别为-3.24%和

-3.82%,如果在某自然保护区内开始有100只鹤,建立描述其数量变化规律的模型,并作数值计算。第2页,课件共10页,创作于2023年2月模型建立记第k年沙丘鹤的数量为xk,年均增长率为r,则第k+1年鹤的数量为

xk+1=(1+r)xkk=0,1,2······已知x0=100,在较好,中等和较差的自然环境下r=0.0194,-0.0324,和-0.0382我们利用Matlab编程,递推20年后观察沙丘鹤的数量变化情况第3页,课件共10页,创作于2023年2月Matlab实现首先建立一个关于变量n,r的函数functionx=sqh(n,r)a=1+r;x=100;fork=1:nx(k+1)=a*x(k);end第4页,课件共10页,创作于2023年2月在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'])第5页,课件共10页,创作于2023年2月利用plot绘图观察数量变化趋势可以用不同线型和颜色绘图rgbcmykw分别表示红绿兰兰绿洋红黄黑白色:+o*.Xsd表示不同的线型第6页,课件共10页,创作于2023年2月

plot(k,y1,k,y2,k,y3)在同一坐标系下画图plot(k,y2,':')>>plot(k,y2,'--')>>plot(k,y2,'r')>>plot(k,y2,'y')>>plot(k,y2,'y',k,y1,':')>>plot(k,y2,k,y1,':')>>plot(k,y2,'oy',k,y1,':')用gtext(‘r=0.0194’),gtext(‘r=-0.0324’),gtext(‘r=-0.0382’)在图上做标记。第7页,课件共10页,创作于2023年2月functionx=fhsqh(n,r,b)a=1+r;X=100;Fork=1:nX(k+1)=a*x(k)+b;end第8页,课件共10页,创作于2023年2月k=(0:20);%一个行向量y1=(20,-0.0324,5);也是一个行向量round([k’,y1’])对k,y1四舍五入,但是不改变变量的值

plot(k,y1)ky1是行向量列向量都可以也可以观察200年的发展趋势,以及在较差条件下的发展趋势,也可以考察每年孵化数量变化的影响。第9页,课件共10页,创作于2023年2月一阶线性常系数差分方程的解、平衡点及其稳定性

自然环境下,b=0人工孵化条件下

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论