微分方程数值解 ■主要问题： ●如何对定义域进行离散？ ●如何对微分方程进行离散？ ●收敛性问题，即步长充分小时，所得到的数值解能否逼近问题的真解 ●误差估计 ·稳定性问题，即舍入误差在以后各步的计算中，是否会无限制扩大 ●计算效率 ·并行计算 6
微分方程数值解 主要问题： ⚫ 如何对定义域进行离散? ⚫ 如何对微分方程进行离散？ ⚫ 收敛性问题，即步长充分小时，所得到的数值解能否逼近问题的真解 ⚫ 误差估计 ⚫ 稳定性问题，即舍入误差在以后各步的计算中，是否会无限制扩大 ⚫ 计算效率 ⚫ 并行计算 ⚫ …… 6
Euler公式 ■ 对定义域[a,b]作等距剖分，即 :a h=5-a.1=0.1..m m ■向前差商公式 -x)=y'(x)+3y"5) h y'=fu,以 y(b) h y(a)=a WN x)-=fx,x》+2"5) ）=+》+行) W2 →y+1=y+hf(x,y,) =a ,,1x=b 7
Euler公式 对定义域 作等距剖分，即 向前差商公式 7 [ , ] a b : , , 0,1, , I i b a x a hi h i m m − = + = = 1 1 2 1 1 ( ) ( ) '( ) ''( ) 2 ( ) ( ) ( , ( )) ''( ) 2 ( ) ( ) ( , ( )) ''( ) 2 ( , ) i i i i i i i i i i i i i i i i i i y x y x h y x y h y x y x h f x y x y h h y x y x hf x y x y y y hf x y + + + + − = + − = + = + + = +
Euler公式 定义：在假设第步计算是精确的前提下，考虑载断误 差T1=y(x)-y1,称T+1为局部载断误差。若T=O(hP) ，则称方法是卫阶相容的，简称相容 ■向前差商公式的局部载断误差： T+1=y(x+1)-y+ h2 =x)+f(,x》+2"5)- h2 =x)+f(x(x》+2y"5)--f(x》 h2 =2) =0(h2) 8
Euler公式 定义：在假设第 步计算是精确的前提下，考虑截断误 差 ，称 为局部截断误差。若 ，则称方法是 阶相容的，简称相容 向前差商公式的局部截断误差： 8 i 1 1 1 ( ) T y x y i i i + + + = − Ti+1 1 1 1 2 1 2 2 2 ( ) ( ) ( , ( )) ''( ) 2 ( ) ( , ( )) ''( ) ( , )) 2 ''( ) 2 ( ) i i i i i i i i i i i i i i i i T y x y h y x hf x y x y y h y x hf x y x y y hf x y h y O h + + + + = − = + + − = + + − − = = 1 1 ( ) p T O h i + + = p
Euler公式 整体载断误差和收敛性：考虑局部载断误差的积果和 传播 e=y(x)- y(x)+hGx》+2")-片-(x,y ≤y(x,)-y+hf(x,y(x)》-f(x,y+T ≤le,+hLy(x)-y+|Tl s(1+hL)e,+T,T=max T, ≤(1+hL)(1+hL)e+T）+T =(1+hL)2le,+(1+hL)+1)T ≤(1+hL)2(1+hL)le,-2+T）+(1+hL)+1)7 =(1+hL)3le-2+(1+hL)2+(1+hL)+17 9
Euler公式 整体截断误差和收敛性：考虑局部截断误差的积累和 传播 9 ( ) 1 1 1 2 1 1 1 ( ) | ( ) ( , ( )) ''( ) ( , )) | 2 ( ) ( , ( )) ( , ) ( ) (1 ) , max (1 ) (1 ) (1 i i i i i i i i i i i i i i i i i i i i i i j j i e y x y h y x hf x y x y y hf x y y x y h f x y x f x y T e hL y x y T hL e T T T hL hL e T T + + + + + − = − = + + − − − + − + + − + + + = + + + + = + ( ) ( ) ( ) ( ) 2 1 2 2 3 2 2 ) (1 ) 1 (1 ) (1 ) (1 ) 1 (1 ) (1 ) (1 ) 1 i i i hL e hL T hL hL e T hL T hL e hL hL T − − − + + + + + + + + + = + + + + + +
Euler.公式 ■整体截断误差和收敛性： ≤… ≤(1+hL)e+(1+hL)'+…+(1+hL)+1)T =+y"e小+二+7 1-(1+hL) ≤1+hg,++h"7 hL ≤+"（+a) 是1阶方法 sew{+)}set(a*冠) 当a=00对，7=o)→e[+)>0,即Euler公式 是收敛的 10
Euler公式 整体截断误差和收敛性： 当 时， ，即Euler公式 是收敛的 10 ( ) 1 0 1 1 0 1 1 0 1 0 ( 1) ( ) 0 0 (1 ) (1 ) (1 ) 1 1 (1 ) (1 ) 1 (1 ) (1 ) (1 ) (1 ) i i i i i i i i hL b a L hL e hL hL T hL hL e T hL hL hL e T hL T hL e hL T T e e e e hL hL + + + + + + + − + + + + + + + − + = + + − + + + + + + + + 0 e h = → 0, 0 2 ( ) 0 ( ) 0 b a L T T O h e e hL − = + → 是1阶方法