差分方程及微分方程数值解

差分方程及微分方程数值解

部门： xxx

时间： xxx

整理范文，仅供参考，可下载自行编辑

东南大学《数学实验》报告实验内容：差分方程及微分方程数值解

一实验目的

熟悉迭代法及微分方程数值方法

二预备知识

<1）了解差分方程稳定性、周期分解、混沌等相关知识 <2）了解欧拉方法、龙格-库特方法。

三实验内容与要求

<一）Volterra方程数值解

Numerical

Calculable

rank

（2）轨线图

四/五阶龙格-库塔公式：

Numerical

Euler

欧拉公式：

（二）差分阻滞增长模型

在t时刻单位时间内的人口数量的变化量仅仅与此时的人口数量x有关(等于右边的值>,其中的r表示人口的固有增长率,N表示能容纳的最大人口数.b5E2RGbCAP

经过单位时间，即一个繁殖周期的种群数量的增长量仅仅与前一个时期的种群数量有关。于是模型为

即

这是一个一阶非线性差分方程.对于给定的初值,我们可以从这个递推公式运用计算机很容易地计算出一些xk.这是在计算机出现以后的一个新的特点.但是我们更关心的是当时间趋于无穷时,即k趋于无穷时, xk的极限如何，即差分方程平衡点的稳定性问题。p1EanqFDPw

求得差分方程(2>的平衡点为0和(b-1>/b.它们分别对应于差分方程(1>的平衡点0和N.

我们将这个差分方程(2>在平衡点附近展开,有

注意到b=1+r>1,平衡点0是不稳定的.

略去高阶项得

因此当|2-b|<1时,平衡点x*＝ (b-1>/b是稳定的. |2-b|>1即 b>3时,平衡点x*是不稳定的.DXDiTa9E3d

<1）稳定情况<）

<

function x=block(b,x0>

图）function x=block(b,x0>

结果

结果

图：

<2）不稳定情况

稳定性分析

b=3.3在1到3范围之间 不稳定

命令与结果<图）

稳定性分析

b=3.63在1到3范围之间 不稳定

命令与结果<图）

<3）倍周期分解

（1） 求稳定平衡点 可得两个平衡点：0.4794 0.8236 程序： <1）求稳定平衡点

四个平衡点：0.5009 0.8750 0.3828 0.8269 程序：

申明：

所有资料为本人收集整理，仅限个人学习使用，勿做商业用途。