1 东南大学《数学实验》报告 实验内容：差分方程及微分方程数值解 一 实验目的 熟悉迭代法及微分方程数值方法 二 预备知识 （1）了解差分方程稳定性、周期分解、混沌等相关知识 （2）了解欧拉方法、龙格-库特方法。 三 实验内容与要求 （一）Volterra 方程数值解 方程0,,,,,dcbadxycxdtdybxyaxdtdx 其中a=1,b=0.1,c=0.5,d=0.04 .2)0(,25)0(yx 命令与结果 在函数编辑器中输入： function dxdt = euler( t,x ) dxdt=[ x(1)*(1-0.1*x(2)) x(2)*(-0.5+0.02*x(1))]; end 四阶龙格-库塔公式： 在命令窗口中输入： tspan=[0 15]; x 0=[25;2]; [t,x ]=ode45(@eu ler,tspan,x 0); plot(t,x (:,1),'r-','LineWidth',0.5); hold on; plot(t,x (:,2),'g-','LineWidth',0.5); hold on; （1） 相图 四阶龙格-库塔公式： 2 axis([0 15 0 125]) legend('x(1)','x(2)') grid on title('The Numerical Solution Of Calculable way of fourth rank Rounge-kutt') plot(x(:,1),x(:,2)) 欧拉方法： 在命令窗口中输入： tspan=[0 15]; x0=[25;2]; [t,x]=ode23(@euler,tspan,x0); plot(t,x(:,1),'r-','LineWidth',0.5); hold on; plot(t,x(:,2),'g-','LineWidth',0.5); hold on; axis([0 15 0 125]) legend('x(1)','x(2)') grid on title('The Numerical Solution Of Euler Equation') >> plot(x(:,1),x(:,2)) 欧拉公式： （2） 轨线图 四/五阶龙格-库塔公式： 3 欧拉公式： （二）差分阻滞增长模型 在t 时刻单位时间内的人口数量的变化量仅仅与此时的人口数量x 有关(等于右边的值),)1()(Nxrxtx 4 其中的r 表示人口的固有增长率,N 表示能容纳的最大人口数. 经过 单位时间，即一个繁殖周期的种群数量的增长量仅仅与前一个时期的种群数量有关。于是模型为 即 这是一个一阶非线性差分方程.对于给定的初值,我们可以从这个递推公式运用计算机很容易地计算出一些 xk.这是在计算机出现以后的一个新的特点.但是我们更关心的是当时间趋于无穷时,即k 趋于无穷时, xk 的极限如何，即差分方程平衡点的稳定性问题。 求得差分方程(2)的平衡点为 0 和(b-1)/b.它们分别对应于差分方程(1)的平衡点 0 和 N. 我们将这个差分方程(2)在平衡点附近展开,有 kkbxxx1:0处 注意到 b=1+r>1,平衡点 0 是不稳定的. *))(...
