《有限差分法在微分方程中的应用》课程论文

课程论文

有限差分法在微分方程中的应用

本学期学习了《微分方程数值解》，本书中有限差分法给我留下的印象比较深刻，下边说说自己在方面的一点理解，请老师指正。

1.有限差分法的基本思想：

当系统的数学模型建立后，我们面对的主要问题就是微分积分方程的求解。基本思想是用离散的只含有限个未知量的差分方程组去近似地代替连续变量的微分方程和定解条件，并把差分方程组的解作为微分方程定解问题的近似解。将原方程及边界条件中的微分用差分来近似，对于方程中的积分用求和或及机械求积公式来近似代替，从而把原微分积分方程和边界条件转化成差分方程组。

2.有限差分法求解偏微分方程的步骤：

区域离散，即把所给偏微分方程的求解区域细分成由有限个格点组成的网格，这些离散点称作网格的节点；

近似替代，即采用有限差分公式替代每一个格点的导数。

逼近求解，换而言之，这一过程可以看作是用一个插值多项式及其微分来代替偏微分方程的解的过程。

从原则上说，这种方法仍然可以达到任意满意的计算精度。因为方程的连续数值解可以通过减小独立变量离散取值的间格，或者通过离散点上的函数值进行插值计算来近似得到。理论上，当网格步长趋近于零时，差分方程组的解应该收敛于精确解，但由于机器字节的限制，网格步长不可能也没有必要取得无限小，那么差分法的收敛性或者说算法的稳定性就显得至关重要。因此，在运用有限差分法时，除了要保证精度外，还必须要保证其收敛性。

3.构造差分法的几种形式：

主要草用的是泰勒级数展开的方法。其基本差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等。其中前两种形式为一阶计算精度，后一种为二阶计算精度。

4.有限差分法的应用:

4.1抛物线形的差分法中的一维常系数抛物线型方程 考虑最简单的以为常系数抛物线型方程

2

2()u u

Lu a f x t x

∂∂=-=∂∂ (,)x t ∈Ω 其中Ω是（x.t ）平面内的给定区域，可以是有节区域或无解区域；a>0是常数，L 是微分算

子。

根据定解条件的不同，可以将上述方程分为两类： 1.初值问题

在区域{(,)|,0}x t x t Ω=-∞<<+∞>上求解方程满足初始条件

(,0)(),u x x x φ=-∞<<+∞的解。

2，初边值问题（混合问题）

在区域{(,)|0,0}x t x l x T Ω=<<<≤内求方程满足初始条件

(,0)(),0u x x x l φ=≤≤和下列边界条件之一的解。

第一边届条件

1,(0,)u t μ=

2(,)u l t μ= 0t T ≤≤

第二边界条件

1,(0,)x u t β=

2(,)x u l t β= 0t T ≤≤

第三边界条件

101(())|()x x u t u r t α=-=

202(())|()x x u t u r t α=-= 0,1,2,0i i t T α≥=≤≤

用适当的差商代替方程中相应的偏导数，可得到以下几种最简差分格式： 古典显示格式：

1

11

(1)2

2k k k k k

j j

j j j k

k h j j u u u u u L u a

f h

τ

++---+≡

-=

古典隐式格式：

1111

11

(2)

112

2k k k k k j j

j j j k k h j

j u u u u u L u

a

f h

τ

+++++-++--+≡

-=

加权六点隐式格式：

11111111(3)122

22(1)(1)k k k k k k k k

j j

j j j j j j k k k

h j

j j u u u u u u u u L u a a f f h h θθθθτ

+++++-+-+⎡⎤--+-+≡

-+-=+-⎢⎥⎢⎥⎣⎦

4.2椭圆型方程边值问题的差分法

考虑如下两点边值问题

au bu cu f '''-++=

01(0),(1)u u u u == (0,1)x ∈Ω=

其中0u ，1u 为常数，系数a=a(x),b=b(x),c=c(x),f=f(x)为一致的充分光滑函数，且满足a(x)>0，c(x)>0.

首先将区间Ω离散化，我们采用剖分部分，取正整数M ，将区间M Ω等分，的M+1个节点：

0101M x x x =<<<=

其中1

(0,1,

),j x jh j M h M

===

。设U 为定义在节点(0,1,)j x jh j M ==上的网格函

数并用i U 近似()j u x 。下边可以得到两点边值问题的有限差分：

11

11

2

22j j j j j h j j

j

j j j U U U U U A U a b c U f h

h

-++--+--++=

0011,U u U u == 0,1,1

j M =- 上述方程组成为差分方程，它的解就是两点边值问题的差分解，当0,1,1j M =-时忧伤

出差分方程可得：

221111

()(2)()22

j j j j j j j j j j a hb U a h c U a hb U h f -+-+++--=

由此得出线性代数方程： AU g = 其中1,2,

1,1,2,1(),()T T M M U U U U g g g g --==，因此(1)(1)M M A R --∈