资源描述：
《差分方程及微分方程数值解》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、东南大学《数学实验》报告实验内容：差分方程及微分方程数值解一实验目的熟悉迭代法及微分方程数值方法二预备知识（1）了解差分方程稳定性、周期分解、混沌等相关知识（2）了解欧拉方法、龙格-库特方法。三实验内容与要求（一）Volterra方程数值解方程其中a=1,b=0.1,c=0.5,d=0.04命令与结果在函数编辑器中输入：functiondxdt=euler(t,x)dxdt=[x(1)*(1-0.1*x(2))x(2)*(-0.5+0.02*x(1))];end四阶龙格-库塔公式：在命令窗口中输入：tspan=[015];x0=[2
2、5;2];[t,x]=ode45(@euler,tspan,x0);plot(t,x(:,1),'r-','LineWidth',0.5);holdon;plot(t,x(:,2),'g-','LineWidth',0.5);holdon;（1）相图四阶龙格-库塔公式：7axis([0150125])legend('x(1)','x(2)')gridontitle('TheNumericalSolutionOfCalculablewayoffourthrankRounge-kutt')plot(x(:,1),x(:,2))欧拉方法：
3、在命令窗口中输入：tspan=[015];x0=[25;2];[t,x]=ode23(@euler,tspan,x0);plot(t,x(:,1),'r-','LineWidth',0.5);holdon;plot(t,x(:,2),'g-','LineWidth',0.5);holdon;axis([0150125])legend('x(1)','x(2)')gridontitle('TheNumericalSolutionOfEulerEquation')>>plot(x(:,1),x(:,2))欧拉公式：（1）轨线图四/五阶龙
4、格-库塔公式：7欧拉公式：（二）差分阻滞增长模型7在t时刻单位时间内的人口数量的变化量仅仅与此时的人口数量x有关(等于右边的值),其中的r表示人口的固有增长率,N表示能容纳的最大人口数.经过单位时间，即一个繁殖周期的种群数量的增长量仅仅与前一个时期的种群数量有关。于是模型为即这是一个一阶非线性差分方程.对于给定的初值,我们可以从这个递推公式运用计算机很容易地计算出一些xk.这是在计算机出现以后的一个新的特点.但是我们更关心的是当时间趋于无穷时,即k趋于无穷时,xk的极限如何，即差分方程平衡点的稳定性问题。求得差分方程(2)的平衡点为
5、0和(b-1)/b.它们分别对应于差分方程(1)的平衡点0和N.我们将这个差分方程(2)在平衡点附近展开,有注意到b=1+r>1,平衡点0是不稳定的.略去高阶项得因此当
6、2-b
7、<1时,平衡点x*＝(b-1)/b是稳定的.
8、2-b
9、>1即b>3时,平衡点x*是不稳定的.（1）稳定情况（）稳定性分析b=1.6,在1到3范围之间稳定命令与结果（图）functionx=block(b,x0)b=1.6;ifnargin<2x0=0.5;稳定性分析b=2.2,在1到3范围之间稳定命令与结果（图）functionx=block(b,x0)b=
10、2.2;ifnargin<2x0=0.5;7endx=x0;n=1:10000;fori=2:10000x(i)=b*x(i-1)*(1-x(i-1));endplot(n,x,'.');end结果图：endx=x0;n=1:10000;fori=2:10000x(i)=b*x(i-1)*(1-x(i-1));endplot(n,x,'.');end结果图：（2）不稳定情况稳定性分析b=3.3在1到3范围之间不稳定命令与结果（图）稳定性分析b=3.63在1到3范围之间不稳定命令与结果（图）7（3）倍周期分解（1）求稳定平衡点可得两个
11、平衡点：0.47940.8236程序：functionx=block(b,x0)b=3.3;ifnargin<2x0=0.5;endx=x0;n=1:100;fori=2:100x(i)=b*x(i-1)*(1-x(i-1));end（2）数值结果ans=Columns1through160.50000.82500.47640.82320.48040.82370.47920.82360.47950.82360.47940.82360.47940.82360.47940.8236Columns17through320.47940.82
12、360.47940.82360.47940.82360.47940.8236（1）求稳定平衡点四个平衡点：0.50090.87500.38280.8269程序：functionx=block(b,x0)b=3.5;ifnargin<2
