常微分方程初值问题.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 5/16 连续问题: 离散问题: n = 1, 2, 3, ······ 6/16 求解常微分方程初值问题的Euler方法 7/16 例 取步长 h = 2/10, 2/20, 2/30, 2/40, 用欧拉法求解. 解: 记 f (x, y) = y－ x y2, xn=0.2n (n = 0, 1, 2,···, N ) 由Euler公式: yn+1 = yn + h( yn－ xn yn2 ) (n = 0, 1, ···,N) y0 = 1 N 10 20 30 40 h 0.2 0.1 0.0667 0.05 误差 0.1059 0.0521 0.0342 0.0256 8/16 n=input(input n:= ); f=inline(y-x.*y.^2); h=2/n;x=h:h:2; y0=1; y(1)=y0+h*f(0,y0); for k=1:n-1 xk=x(k);yk=y(k); y(k+1)=yk+h*f(xk,yk); end t=-1:h/2:5; yy=1./(t-1+2*exp(-t)); plot(0,1,ro,x,y,ro,t,yy) y1=1./(x-1+2*exp(-x)); error=max(abs(y-y1))) o : 数值解 ----: 准确解 9/16 解析解: y’ = f (x, y) 梯形公式: 左矩形公式 用数值积分方法离散化常微分方程 10/16 算法: k1 = f(xn , yn) , k2 = f( xn+1 , yn+ h k1) 由梯形公式推出的预-校方法: 11/16 n=input(input n:= ); f=inline(y-x.*y.^2); h=2/n;x=h:h:2; y0=1; k1=f(0,y0);k2=f(h,y0+h*k1); y(1)=y0+0.5*h*(k1+k2); for k=1:n-1 xk=x(k);yk=y(k); k1=f(xk,yk);k2=f(xk+h,yk+h*k1); y(k+1)=yk+0.5*h*(k1+k2); end t=-1:h/2:5; yy=1./(t-1+2*exp(-t)); plot(0,1,ro,x,y,ro,t,yy) y1=1./(x-1+2*exp(-x)); error=max(abs(y-y1)) input n:= 10 error = 0.0123 12/16 o : 数值解 ----: 准确解 *