数值分析23常微分方程初值问题.pptVIP

* 常微分方程初值问题 求数值解的Euler方法 修正的Euler方法 局部截断误差与p阶精度 MATLAB求初值问题命令 《数值分析》 23 ? ? ? ? ? 例1.罗伦茨模型 求解区域: 0≤ t ≤ 80 初始条件：x(0)=0, y(0)= 0, z(0) = 0.0001 取 ? =8/3, ? =28, ? =10 2/16 例2.Lotka-Volterra捕食者与被捕食者模型 取 a =1, b=0.01, c =1, d =0.02 求解区域: 0≤ t ≤ 15 初始条件: x(0)=20, y(0)=20 3/16 ---- x(t) ---- y(t) 例3.人口模型(Logistic ): P’ =rP ( 1 – P/K ) 取: r = 0.02, K = 500 初始条件: P(0)=76 例4. 传染病模型: P’ =rP ( L – P ) L = 25000, r = 0.00001 初始条件: P(0)=2 4/16 常微分方程初值问题: 例 数值方法：0 = x0 x1 x2 ······ xN ≤2 求未知函数 y=y(x) x∈(x0, ∞) 求: y(x1), y(x2), y(x3), ······· y(xN) 的近似值 y1, y2, y3, ······, yN