| G | 2 1 r(1 cosh)2 r 2 sin2 h
1 2r(r 1)(1 cosh) 1
从而获得原格式的稳定性条件 1 r 0
即 a 0 且 h a
， ②：关于时间、空间的一阶偏导数分别利用一阶向前 差商和一阶向后差商近似，即有
u
u( x j , tk1 ) u( x j , tk )
uk1 j
r
uk j1
(1
r)
ukj ,
u0j ( x j ),
jZ, k 0
也不难得到此格式的增长因子为
G 1 r r eih 1 r(1 cosh) ir sinh
| G | 2 1 r(1 cosh)2 r 2 sin2 h
1 2r(r 1)(1 cosh) 1
误差为 O( )
u
u( x j1 , tk ) u( x j , tk )
x ( x j ,tk )
h
误差为 O(h)
将上面的式子代入离散方程，可得
u( x j , tk1 ) u( x j , tk ) a u( x j1, tk ) u( x j , tk ) O( h)
h
j Z , k 0,
五、Lax-Wendroff 格式
u a u
t
x
2u t 2
t
a
u x
a
x
u t
a
x
a
u x
a2
2u x 2
再根据泰勒公式就有
u( x j , tk1 )
u( x j , tk )
u t
( x j ,tk )
2
2
2u t 2
( x j ,tk )
O( 3 )
u( x j , tk ) a
只要将其沿特征线投影到 x 轴上，得到投影点 (0 , 0) ，
则 u( x0 , t0 ) (0 ) .
考虑一维双曲型对流方程：
u( x, t
t)
a
u( x, x
t)
0,
x , t 0,
u( x, 0) ( x), x ,
1. 区域剖分(区域离散） 将原方程的上半平面求解区域分割成矩形一致网格。
u(
x
j
,
t0
)
(
x
j
),
j Z,
3. 处理方程 中的偏导数 对偏导数用不同的差商近似将建立不同的差分格式。 下面进行具体的讨论。
二、迎风格式（Upwind Scheme）
①：关于时间、空间的一阶偏导数都用向前差商近似，
u
u( x j , tk1 ) u( x j , tk )
t ( x j ,tk )
u
u( x j
, tk1 )
1 2
u( x j1 , tk ) u( x j1 , tk )
t ( x j ,tk )
关于空间的一阶偏导数仍用二阶中心差分，即
u
u( x j1 , tk ) u( x j1 , tk )
x ( x j ,tk )
2h
再将上述近似代入离散方程，可得
u(
x
③：前面讨论了关于时间和空间的一阶偏导数均用
一阶差商近似的情况，接下来容易想到可以对空间
的偏导数采用二阶中心差分来近似，从而有
u
u( x j , tk1 ) u( x j , tk )
t ( x j ,tk )
u
u( x j1 , tk ) u( x j1 , tk )
x ( x j ,tk )
u( x j , 0) ( x j ), j Z ,
将数值解 ukj 代替精确解 u( x j , tk ) 并忽略高阶小项， 则可以建立以下差分格式：
ukj
1
u
k j
a
uk j1
ukj
h
0,
j Z , k 0,
u0j
(xj
),
j Z,
可见上述格式的局部截断误差为 O( h)
设 ukj vkei xj ，相当于对数值解进行变量分离， 对数值格式稳定性的考察现在就转化为对振幅 是否 会放大进行讨论。
如果 vk1 G vk , 则 G 就称为增长因子，且 |G | 1
是数值格式稳定的充要条件，也称为 Von Neumann 条件。
现在研究上述格式的稳定性。
uk1 j
2
j Z , k 0,
,
tk1 )
a
u( x j1 ,
tk
) u( x j1 , tk 2h
)
O(
2
h2
),
u(
x
j
,
0)
(
x
j
),
j Z,
将数值解 ukj 代替精确解 u( x j , tk ) 并忽略高阶小项，
则可以建立以下蛙跳格式：
ukj
1
2
ukj
1
a
uk j1
uk j 1
2h
知，
当
dx a dt
时，就有
du dt
0
.
即存在一族特征线
x a t 其中 为任意常数，
使得在这样的特征线上有 du 0 ，也就是 u 值为常数。
dt
t
x at x a t 0
.
(x0 , t0 )
.
O
0
x
要获得在 x- t 平面上的任意一点 ( x0 , t0 ) 处的函数值，
u0j ( x j ), j Z ,
可以改写为
uk 1 j
1 (1 2
r )ukj1
1 (1 2
r
)
uk j1
,
u0j ( x j ),
j Z , k 0,
利用分解式可以得到其增长因子为
G 1 r eih 1 r eih cosh ir sinh
2
2
从而 | G | 2 cos2 h r 2 sin2 h 1 (r 2 1)sin2 h
(1
r )ukj
r
uk j1
,
u0j ( x j ),
r a , j Z, k 0,
h
易见， v k1ei x j (1 r )v kei x j r v ek i x j1
从而 G 1 r r eih 1 r(1 cosh) ir sinh
为使数值格式稳定，则增长因子 G 必须满足
对流方程的差分法
一、研究对象
1. 研究的对象—— 对流方程（一阶双曲型）.
u( x, t)
t
a
u( x, t) x
0,
x , t 0,
u( x, 0) ( x), x ,
易见，此方程有精确解 u( x, t ) ( x at ).
事实上，由
du u dx u dt x dt t
4.差分格式的求解
ukj
1
u
k j
a
uk j1
ukj
h
0,
j Z , k 0,
u0j
( x j ),
j Z,
uk1 j
(1
r )ukj
r
uk j1
,
u0j ( x j ),
r a , j Z , k 0,
h
—— 时间渐进显格式
5. 用谐波分析方法利用增长因子来讨论稳定性。
u x
( x j ,tk )
a2 2
2
2u x 2
( x j ,tk )
O( 3 )
上式中一阶、二阶偏导都用中心差分来近似，
u( x j , tk1 )
u( x j , tk ) a
ukj
1
1 2
( u kj 1
uk j1
)
a
uk j1
uk j1
2h
0,
j Z , k 0,
u0j ( x j ), j Z ,
易见，其局部截断误差为
O(
h2 )
O
h2
.
Lax-Friedrichs格式
ukj
1
1 2
( ukj 1
uk j1
)
a
uk j1
uk j 1
2h
0,
j Z , k 0,
0,
j Z , k 0,
u0j
(xj
),
j Z,
可见上述格式的局部截断误差为 O( 2 h2 )
上述格式还可简写为三层格式
uk1 j
r
(
uk j1
uk j1
)
uk 1 j
0,
u0j ( x j ),
j Z , k 0,
也不难得到此格式的增长因子为
G ir sinh 1 r 2 sin2 h
从而有稳定性条件要求 0 r 1
即 a 0 且 h a
综上，我们有以下迎风格式：
a < 0时 a > 0时
ukj
1
ukj
a
uk j1
u
k j
h
0,
j Z , k 0,
u0j
(xj
),
j Z,
且稳定性条件为
ukj
1
ukj
a
ukj
uk j 1
h
0,
j Z , k 0,
| G | 2 cos2 h r 2 sin2 h 1 (r 2 1)sin2 h
可见，当 | r | 1 时就有 | G | 1 ，从而数值格式稳定。 当根r 据 a局 部取截定断为误常差数时O(， Lahx2 )-FrOiedhri2chs是知一，阶格式。
h
下面介绍一个二阶格式，通过泰勒公式及原方程 变形而获得。