• 定义：对于一足够光滑函数 u ，若时间步长 t ，空间步长 x
趋近于0时，差分方程截断误差 Rnj 对于每一点 x j , tn 都趋近于
0，则该差分方程 Lunj 0 逼近微分方程 Lu 0，即差分方程与
微分方程是相容的。
• 差分方程相容性可以通过Taylor展开方法来证明。例如，扩散方
而是趋近于某值，或结论并不是对每个点 x j , t都n 成立，则差分方
程就不满足相容性条件，差分方程也就不逼近于微分方程。
③
相容性条件不仅要求差分方程截断误差R
n j
趋近于0，而且要求差分方
程定解条件截断误差rjn 也同时趋近于0。
④ 差分格式有两种不同形式的相容性，即无条件相容和有条件相容。
定义：在某一个时刻tn存在计算误差

n，若在
j
tn1时刻满足：

n1 j

k

n j
或

n j
k

0 j
条件，则差分方程是稳定的。
这里定义：
1

n j





n j
2
x
2

2

是某种定义的范数。
下面我们用几个简单的例子来说明差分方程稳定性概念。
（1）对流方程FTFS差分方程为：
• Lax定理：对于适定和线性的初值问题微分方程，若逼近它的差分方程和它
是相容的，则差分方程稳定性是差分方程收敛性的充分和必要条件。
• Lax定理可以形象地表示为：
e0j
O(x, t)
(d)
在t=0时，差分方程的初始条件应该是完全准确的，即：
u
0 j
(x j
),e0j

u0

u0j

0
即：
max j
en1 j
O(x, t)
(e)
即差分方程离散化误差和截断误差是相同数量级，因此，若
R
n j
→0，则：
lim t0
max
j
e
n1 j

0.8u
n j 1ຫໍສະໝຸດ ；当a=1，Δx=0.1，r=1.0，则有：
u
n1 j

u
n j1
；
当a=1，Δx=0.1，r=2.0，则有：
u
n1 j

u
u j

2u
n j1
。
图b中给出了上述不同条件下差分方程计算误差的图解。从图中 可以发现，当r=1.0时，差分方程解和微分方程解是一致的；当 r=0.8时，在差分方程解的两端有耗散现象，当r=2.0时，差分方 程解会出现振荡，并且在t=nΔt继续增加时，振荡也继续加剧， 直到计算完全失败。 数值分析表明，FTBS差分方程只有在r 1.0时计算才是稳定，当r ＞1.0时差分方程计算是不稳定。
程的FTCS差分格式为：
u n 1 j

u
n j
t

un j 1

2u
n j
x2

u
n j 1

0
•
把
u n1 j
作为t的函数，在 tn
邻域展开成Taylor级数，把
un j 1
和 u nj1作为x的函数，在 x j 邻域展开成Taylor级数：
u n 1 j

u
n j


n
j
x3


(x4
)
un j 1

u
n j


u x
n j
x

1 2

2u x2
n
j
x2

1 6

3u x3
n
j
x3

(x4
)
u u 将
u
n1 j
、
nj 1和
n j 1
代入FTCS格式中，即可得到：
u
un1 j

(1 r)unj

ru
, n
j 1
u0j
(xj )
(a)
设求解区域内任意一点 x p ，t p ，它的微分方程精确解为u，
差分方程解为
u
n j
，则离散化误差为
enj
u unj
，把差分方
程和微分方程相减可得离散化误差方程：
en1 j

(1
r)enj

renj1

O(x, t)
计算力学基础
第二章 有限差分方法
2.4 差分方程的相容性、收敛性和稳定性
一个微分方程采用不同的方法可以得到不同的差分方程。那么， 我们要问，对于这些不同的差分方程是否都同样有效，同样可靠，而 且能得到同样的计算结果呢？
答案是否定的。事实上，不同的差分方程和原方程有完全不同的 对应关系，它们具有各自不同的性质，因此，数值结果也完全不同。 在这些差分方程中有些差分方程是有效的、可靠的；些差分方程只有 在一定的条件下是有效的、可靠的；有些差分方程则是完全无效的、 不可靠的。所以，如何判断和分析差分方程有效性和可靠性就成为非 常必要和现实的问题了。
为离散化误差。
定义2：节点 xp , t p 为微分方程求解区域 内任意一点，当
当
x xp,t tp
时，差分方程数值解
u
n j
趋近于微分方程
精确解
u
，即
enj

u

u
n j

0
，则差分方程收敛于微分方程。
差分方程收敛性有两种证明方法，直接证明法和数值试验法。
一、直接证明法
对流方程 u a u 0 的FTBS差分格式为： t x
目完全淹没了，所求得差分方程数值解已经没有任何意义了，
因此，FTFS差分方程是不稳定的。
（2）对流方程FTBS差分格式的误差传播方程为：

n1 j


n j

r

n j

n j1

1
r

n j

r
n j1
(b)
当a＞0，a t x
1
时，

n1 j

max
j

n j
u t
n j
t

1 2

2u t 2
n j
t 2

1 6

3u t 3
n j
t 3

(t 4
)
un j 1

u
n j


u x
n
j
x

1 2

2u x2
n
j
x2

1 6

3u x3


0
(f)
x0
由此可知，FTBS格式在a＞0， a t 1 时，是收敛的。 x
二、数值试验法
数值试验法基本思想是用差分方程求出FTBS数值解，然后和 微分方程精确解进行比较，确定差分方程是否收敛。
直接证明法比较简单，但是只有很少几个差分方程可以采用直 接证明法来证明其收敛性，而数值试验法又非常麻烦，一般来 说，很难用数值试验结果严格证明差分方程是否收敛。总的说 来，不管是采用直接证明法，还是数值试验法，要证明差分方 程收敛性都是比较困难的。
t
1。但是，对于
x
不同的a，Δt，Δx，FTBS差分格式的稳定条件是不同的（见
图b）。

n1 j


n j
r

n j

n j1

1
r

n j

r
n j1
(b)
通过对（b）式的数值分析可知：
当a=1，Δx=0.1，r=0.8，则有：
u n 1 j

0.2u
n j
① 相容性是对求解区域内任意一点差分方程逼近于微分方程的
程度，相容性是有限差分算法（包括有限体积算法）首先必
须满足的有效性条件。
② 相容性要求对于求解区域内任意点 x j , tn ，在 t, x 同时趋近于0，
截断误差Rnj 趋近于0。如果 t, x 不是同时趋近于0或并不趋近于0，

t

2u x2
n j


1 2
2u t 2
t

1 6
3u t 3
t 2


(t
3
)


(x
2
)

n
j
•
当 t, x 0时，上等式右侧所有项都趋近0，差分方程趋近
于原微分方程，即FTCS差分方程和原方程是相容的。
•
关于差分方程相容性需要作以下说明：

n j
，而计算到n+100时
刻，（xj，tn+100）点的计算误差将发展到
r

1
100

n j

r100
n j1
，
假定只有在节点（xj，tn）上存在误差

n j
，其他各节点的计算
误差为零，则若取r=0.8，则

n100 j

1.8100

n j

3.37
1025

n j
。
由此可以看出，这个计算误差必定会将差分方程精确解原来面
关于差分方程收敛性需要作以下说明： (1) 差分方程收敛性表示差分方程数值解和微分方程精确解