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

1、有限差分法在微分方程中的应用本学期学习了微分方程数值解，本书中有限差分法给我留下的印象比较深刻，下边说说自己在方面的一点理解，请老师指正。1 .有限差分法的基本思想：当系统的数学模型建立后，我们面对的主要问题就是微分积分方程的求解。 基本思想是用离散的只含有限个未知量的差分方程组去近似地代替连续变量的 微分方程和定解条件，并把差分方程组的解作为微分方程定解问题的近似解。将 原方程及边界条件中的微分用差分来近似， 对于方程中的积分用求和或及机械求 积公式来近似代替，从而把原微分积分方程和边界条件转化成差分方程组。2 .有限差分法求解偏微分方程的步骤：区域离散，即把所给偏微分方程的求解区域细分成由

2、有限个格点组成的网格，这些离散点称作网格的节点；近似替代，即采用有限差分公式替代每一个格点的导数。逼近求解，换而言之，这一过程可以看作是用一个插值多项式及其微分来代 替偏微分方程的解的过程。从原则上说，这种方法仍然可以达到任意满意的计算精度。因为方程的连续数值 解可以通过减小独立变量离散取值的问格， 或者通过离散点上的函数值进行插值 计算来近似得到。理论上，当网格步长趋近于零时，差分方程组的解应该收敛于 精确解，但由于机器字节的限制，网格步长不可能也没有必要取得无限小， 那么 差分法的收敛性或者说算法的稳定性就显得至关重要。因此，在运用有限差分法时，除了要保证精度外，还必须要保证其收敛性。3

3、.构造差分法的几种形式：主要草用的是泰勒级数展开的方法。其基本差分表达式主要有三种形式：一 阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等。其中前两种形式 为一阶计算精度，后一种为二阶计算精度。4 .有限差分法的应用：4.1 抛物线形的差分法中的一维常系数抛物线型方程 考虑最简单的以为常系数抛物线型方程2Lu- au2 = f (x) (x ,t )FtFx其中Q是(x.t)平面内的给定区域，可以是有节区域或无解区域；a>0是常数，L是微分算子。根据定解条件的不同，可以将上述方程分为两类：1.初值问题在区域( =( x,t) | -co <x <+oc,t >0

4、上求解方程满足初始条件u(x,0) = 4(x),<x < +OC 的解。2,初边值问题(混合问题)在区域C=(x,t)|0<x<l,0 <x <T内求方程满足初始条件u(x,0) =4(x),0 Wx Ml和下列边界条件之一的解。第一边届条件u(0,t)=1,u(l,t) =20 MtMT第二边界条件Ux(0,t) ”,Ux(l,t)，20 <t <T第三边界条件(u -：i(t)ux) |xz0 = ri(t)(u -： 2(t)ux)|x£ = r2(t)= i -0j = 1, 2,。订用适当的差商代替方程中相应的偏导数，可得到

5、以下几种最简差分格式: 古典显示格式：k 1 k k _ uj - u jLh uj =akk kuj 1 -2u j uj 4 kh27古典隐式格式:k 1 k k 1 _ uj- u jLh uj -aTk 1 八 k 1 k 1uj 1 -2uj uj k 1h二 L加权六点隐式格式:k 1 k - klkl k 1k k k kLh3)uk 三 juL_ aeUj2u2 +Uj,+'5¥22 + Uj= f二+( fjt Ih2h2 i j4.2椭圆型方程边值问题的差分法 考虑如下两点边值问题-au . bu cu = fu(0) =u0U1) =u1x J =(0,

6、1)其中U0 , ui为常数，系数a=a(x),b=b(x),c=c(x),f=f(x)为一致的充分光滑函数，且满足 a(x)>0, c(x)>0.首先将区间 维离散化,我们采用剖分部分，取正整数M,将区间MC等分，的M+1个节点：0 = x0 ：二 x1 ：二 | H ：二 xM = 1 1.其中xj = jh (j =0,1J|M ), h = M。设U为定义在下点xj = jh( j = 0,1,1 |M )上的网格函数并用Ui近似u(xj)。下边可以得到两点边值问题的有限差分：I U i4 -2U j U j 1 U j 1 -U i4AhUjL-ajM也-yUj =fjh

7、2hU0=u0,U1=u1j=0,1(lM- 1上述方程组成为差分方程，它的解就是两点边值问题的差分解，当 j =0,1,| M -1时忧伤 出差分方程可得：1212,- 2hbj)Uj(2aj h&)Uj -2hbj)Uj1 =h2fj由此得出线性代数方程：AU = g其中 U =(U1,U2"Um J,g =§2,|® J ,因此 awr(M4)(M22a +h2h1Ya2 +2hb2) A = (aj)=，1,、 (a1 2hh)2a2 h2b1-(a3 2hb3)，1、-(a2 - hb2)2a3 h2b32aM _1 h2bMh2f1，1,.、(

8、a 2hb)u。h2f2 h fM _221h fM 1 ' (aM J - - hbM j)u1显然，当h充分小的时候矩阵 A为按行对角占优的矩阵，即1,I2M -且存在i0,使得：,0 >Z闻|j 04.3:双曲线方程的有限差分法我们考虑线性对流方程的初值问题a =0, (x,t) 'J.t二 xu(x,0) = ' (x),-二二 x ；二用差分法求解该微分方程的过程和用差分法求解抛物线型方程相类似。以下我们看看用差分法求解该方程的几种格式，不一一写出具体步骤，只列出结果。1。 Courant-Isaacson-Rees格式(迎风格式)其中 a：、1。a|

9、a)2k 1 kk kk kuj uj +a +uj uj= a_uj 卡uj =0 hh2.Lax-Friedrichs 格式。k 11 k kuj- 5 (uj1 uj) uk 1-uk 12 auj1 uj-=02h3。Lax-wendroff 格式。2ukk1 =u： 一23：.1 .u：) r2(u：.1.2ujk u：)4 .蛙跳(Leap-Frog)格式。k kk ku u- u uuj书 uj+ a uj4 uj_ 022h 一5 .Crank-Nicolson 格式。k 1 k r,k1 k 1 k k、u j=uj-(uj1-ujluj1- ujj.)4只是一些概括性的知识

10、点，下以上是本书中差分法在解决不同类型的初边值问题中的应用, 边我们仅用一个具体的实例来说明差分法的求解问题过程。考虑扩散方程的第一初边值问题：.2当 0 :二 x ：二 1,t 0_u二u二 2，.-t二 xu(x,0) =sinx,当0 _x _1时u(0,t) =u(1,t) =0,当t _0用分离变量法可得其解析式为：u(x,t) = e_2tsin二x,0 £xE1； t -0取j =10,h =0.1,毛=jh(j =0,1J|J), £为时间步长，r=：为网比，对于不同的r,h用加权六点隐式格式:k 1 k .(3) k _ uj - uj Lh uj k 1

11、 八 k 1 k 1-aBum-2) uj+a(16)h2kk ku j 1 -2u j u jh2-fjk1(1-) fjk计算上述问题的解u(x,t)在(0.5,0.5)处的近似值，计算结果如下表所示，上述问题的解析式在该点的值为 u (0.5,0.5) =0.00719188。用加权六点隐式格式计算解u (0.5,0.5)得近似值re =0古典显示格式0 =0.5六点对称格式0 =0.8加权六点隐式格式日=1古典隐式格式0.250.007046460.007486960.007758880.007943340.50.006616560.007481470.008030890.008409901.01.1056e+0.0060.007459540.008580110.009378182.0-2.058e+0030.007371960.009695640.011449145.00.001376620.006766860.013110770.0186116510.0-0.000000000.004733130.018664550.03295445 5.总结：