《2022年第十章-偏微分方程数值解法》由会员分享,可在线阅读,更多相关《2022年第十章-偏微分方程数值解法(30页珍藏版)》请在金锄头文库上搜索。
1、精品资料欢迎下载第十章 偏微分方程数值解法偏微分方程问题,其求解十分困难。除少数特殊情况外,绝大多数情况均难以求出精确解。因此,近似解法就显得更为重要。本章仅介绍求解各类典型偏微分方程定解问题的差分方法。1 差分方法的基本概念1.1 几类偏微分方程的定解问题椭圆型方程:其最典型、最简单的形式是泊松(Poisson )方程),(2222yxfyuxuu特别地,当0),(yxf时,即为拉普拉斯( Laplace )方程,又称为调和方程02222yuxuuPoisson 方程的第一边值问题为),(),(),(),(),(2222yxyxuyxyxfyuxuyx其 中为 以为 边 界 的 有 界 区
2、域 ,为 分 段 光 滑 曲 线 ,称为定解区域,),(yxf,),(yx分别为,上的已知连续函数。第二类和第三类边界条件可统一表示为精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 1 页,共 30 页精品资料欢迎下载),(),(yxuuyxn其中n为边界的外法线方向。当0时为第二类边界条件,0时为第三类边界条件。抛物型方程:其最简单的形式为一维热传导方程220 (0)uuaatx方程可以有两种不同类型的定解问题:初值问题xxxuxtxuatu)()0 ,(, 0022初边值问题22120 0, 0( ,0)( )0(0, )( ),( , )(
3、)0uuatTxltxu xxxlutg tu l tgttT其中)(x,)(1tg,)(2tg为已知函数,且满足连接条件精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 2 页,共 30 页精品资料欢迎下载)0()(),0()0(21glg边界条件)(),(),(),0(21tgtlutgtu称为第一类边界条件。第二类和第三类边界条件为)()()()(22101tgutxutgutxulxxTt0其中0)(1t,0)(2t。当0)()(21tt时,为第二类边界条件,否则称为第三类边界条件。双曲型方程 : 最简单形式为一阶双曲型方程0xuatu物理中
4、常见的一维振动与波动问题可用二阶波动方程22222xuatu描述,它是双曲型方程的典型形式。方程的初值问题为精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 3 页,共 30 页精品资料欢迎下载xxtuxxuxtxuatut)()()0,(,0022222边界条件一般也有三类,最简单的初边值问题为2222201200, 0( ,0)( ),( )0(0, )( ),( , )( )0tuuatTxltxuu xxxxltutg tu l tg ttT1.2 差分方法的基本概念差分方法又称为有限差分方法或网格法,是求偏微分方程定解问题的数值解中应用最广
5、泛的方法之一。它的基本思想是:先对求解区域作网格剖分,将自变量的连续变化区域用有限离散点(网格点)集代替;将问题中出现的连续变量的函数用定义在网格点上离散变量的函数代替;通过用网格点上函数的差商代替导数,将含连续变量的偏微分方程定解问题化成只含有限个未知数的代数方程组(称为差分格式)。如果差分格式有解,且当网格无限变小时其解收敛于原微分方程定解问题的解,则差分格式的解就作为原问题的近似解(数值解)。因此,用差分方法求偏微分方程定解问题一般需要解决以下问题:(1)选取网格;(2)对微分方程及定解条件选择差分近似,列出差分格式;(3)求解差分格式;(4)讨论差分格式解对于微分方程解的收敛性及误差估
6、计。下面,用一个简单的例子来说明用差分方法求解偏微分方程问题的一般过程及差分方法的基本概念。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 4 页,共 30 页精品资料欢迎下载设有一阶双曲型方程初值问题。)()0 ,(, 00xxuxtxuatu(1) 选取网格:2h -h 0 h 2h 3h 首先对定解区域 0,),(txtxD作网格剖分,最简单常用一种网格是用两族分别平行于x轴与t轴的等距直线khxxk,(0, 1, 2,0,1,2,)jttjkj将D分 成 许多小矩形区域。这些直线称为网格线,其交点称为网格点,也称为节点,h和分别称作x方向和
7、t方向的步长。这种网格称为矩形网格。(2) 对微分方程及定解条件选择差分近似,列出差分格式:如果用向前差商表示一阶偏导数,即精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 5 页,共 30 页精品资料欢迎下载2211(,)12(,)(, )(, )(,)2(,)(, )(,)2kjkjkjkjkjxxtkjkjkjtxtu xtu x tuhuxh txhu x tu x tuuxtt其中1,021。方程0uuatx在节点),(jktx处可表示为htxutxuatxutxujkjkjkjk),(),(),(),(11),(2),(21222jkxj
8、ktthxuahtxu), 2, 1 , 0, 2, 1, 0(),(jktxRjk其中(,0)() (0,1,2,)kku xxk。由于当,h足够小时,在式中略去),(jktxR,就得到一个与方程相近似的差分方程1,1,0kjkkkjjjuuuuah精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 6 页,共 30 页精品资料欢迎下载此处,jku,可看作是问题的解在节点),(jktx处的近似值。同初值条件), 2, 1, 0()(0,kxukk结合,就得到求问题的数值解的差分格式。式)(),(2),(2),(1222hOthxuahtxutxRjk
9、xjktjk称为差分方程的截断误差。如果一个差分方程的截断误差为)(pqhOR, 则称差分方程对t是q阶精度,对x是p阶精度的。显然,截断误差的阶数越大,差分方程对微分方程的逼近越好。若网格步长趋于0 时,差分方程的截断误差也趋于0,则称差分方程与相应的微分方程是相容的。这是用差分方法求解偏微分方程问题的必要条件。如果当网格步长趋于0 时,差分格式的解收敛到相应微分方程定解问题的解,则称这种差分格式是收敛的。2 椭圆型方程第一边值问题的差分解法本节以 Poisson 方程为基本模型讨论第一边值问题的差分方法。2.1 差分格式的建立考虑 Poisson 方程的第一边值问题),(),(),(),(
10、),(2222yxyxuyxyxfyuxuyx精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 7 页,共 30 页精品资料欢迎下载取,h分别为x方向和y方向的步长,如图所示,以两族平行线khxxk,( ,0, 1, 2,)jyyjk j将定解区域剖分成矩形网格。节点的全体记为(,), ,kjkjRxyxkh yjkj为整数。定解区域内部的节点称为内点,记内点集R为h。边界与网格线的交点称为边界点,边界点全体记为h。与节点),(jkyx沿x方向或y方向只差一个步长的点),(1jkyx和),(1jkyx称为节点),(jkyx的相邻节点。如果一个内点的四
11、个相邻节点均属于,称为正则内点,正内点的全体记为) 1(,至少有一个相邻节点不属于的内点称为非正则内点,非正则内点的全体记为)2(。问题是要求出第一边值问题在全体内点上的数值解。为 简 便 , 记),(),(jkyxjk,( , )(,)kju k ju x y,),(,jkjkyxff。对正则精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 8 页,共 30 页精品资料欢迎下载内点) 1(),(jk,由二阶中心差商公式44422(42( , )2(4)222(4)22( , )(1, )( , )( , )(1, )12(1, )2 ( , )(1
12、, )(12( ,1)2 ( , )( ,1)(,12xk jkxkjyk ju kju k ju k ju kjuhhhuxhu kju k ju kjhuxhuu k ju k ju k juxyyPoisson方 程2222(,)uufxyxy在 点),(jk处可表示为22(1, ) 2 ( , )(1, )( ,1) 2 ( , )( ,u kju k ju kju k ju k ju kh其中(),(12),(12),(22)4(21)4(244hOyxuyhxuhjkRjkxjkx为其截断误差表示式,略去),(jkR,即得与方程相近似的差分方程jkjkjkjkjkjkjkfuuuh
13、uuu,21,1,2, 1, 122精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 9 页,共 30 页精品资料欢迎下载式中方程的个数等于正则内点的个数,而未知数jku,则除了包含正则内点处解u的近似值外,还包含一些非正则内点处u的近似值,因而方程个数少于未知数个数。在非正则内点处Poisson 方程的差分近似不能按上式给出,需要利用边界条件得到。边界条件的处理可以有各种方案,下面介绍较简单的两种。(1)直接转移用最接近非正则内点的边界点上的u值作为该点上u值的近似,这就是边界条件的直接转移。例如,点),(jkP为非正则内点,其最接近的边界点为Q点
14、,则有)2(,),()()(jkQQuujk上式可以看作是用零次插值得到非正则内点处u的近似值,容易求出,其截断误差为)(hO。将上式代入,方程个数即与未知数个数相等。(2)线性插值这种方案是通过用同一条网格线上与点P相邻的边界点与内点作线性插值得到非正则内点),(jkP处u值的近似。由点R与T的线性插值确定)(Pu的近似值jku,,得)()(,TdhdRdhhujk其中RPd,其截断误差为)(2hO。将其与方程相近似的差分方程联立,得到方程个数与未知数个数相等的方程组,求解此方程组可得 Poisson 方程第一边值问题的数值解。上面所给出的差分格式称为五点菱形格式,精选学习资料 - - -
15、- - - - - - 名师归纳总结 - - - - - - -第 10 页,共 30 页精品资料欢迎下载jkjkjkjkjkjkjkfuuuhuuu,21,1,2, 1, 122实际计算时经常取h,此时五点菱形格式可化为1,1,1,1,21(4)kjkjk jk jk jk juuuuufh简记为21hjkjkfu,其中jkjkjkjkjkjkuuuuuu,1,1, 1, 1,4。例 1 用五点菱形格式求解拉普拉斯(Laplace )方程第一边值问题)1lg(),(),(),(222222yxyxuyxyxfyuxu其中 1,0),(yxyx。取31h。(0,0 ) (1,0) (2,0)
16、(3,0) 解 网格中有四个内点,均为正则内点。由五点菱形格式,得方程组精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 11 页,共 30 页精品资料欢迎下载0)4(10)4(10)4(10)4(12,21 ,23,22, 12,322, 11 ,13, 12,02,221 ,20,22,21, 11 ,321, 10, 12, 11,01 ,22uuuuuhuuuuuhuuuuuhuuuuuh代 入 边 界 条 件1,02,00,10,21,32,33,13,21625lg,lg991013lg,lg992534lg,lg993740lg,lg9
17、9uuuuuuuu其解为2756919.01,1u,4603488.01,2u,3467842.02,1u,5080467.02,2u当h时,对1,1,1,1,21(4)kjkjk jk jk jk juuuuufh利用点),(jk,)1, 1(jk,)1, 1(jk构造的差分格式,称为五点矩形格式,简记为精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 12 页,共 30 页精品资料欢迎下载221hjkjkfu,其中jkjkjkjkjkjkuuuuuu,1, 11, 11, 11, 1,4, 其截断误差为244444244( , )( , )6()
18、12k jhuuuR k jO hOxxyy五点菱形格式与矩形格式的截断误差均为)(2hO,称它们具有二阶精度。如果用更多的点构造差分格式,其截断误差的阶数可以提高,如利用菱形格式及矩形格式所涉及的所有节点构造出的九点格式就是具有四阶精度的差分格式。3 抛物型方程的差分解法以一维热传导方程)0(022axuatu为基本模型讨论适用于抛物型方程定解问题的几种差分格式。3.1 差分格式的建立首先对xt平面进行网格剖分。分别取,h为x方向与t方向的 步 长 , 用 两 族 平 行 直 线(0, 1, 2,)kxxkh k,)2, 1 , 0( jjttj,将xt平面剖分成矩形网格,节点为(,) (0
19、, 1, 2,0,1,2)kjxtkj。为简便,记),(),(jktxjk,),(),(jktxujku,)(kkx,)(11jjtgg,精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 13 页,共 30 页精品资料欢迎下载221122( ),( ),( )jjjjjjggttt。(一)微分方程的差分近似在网格内点),(jk处,对tu分别采用向前、向后及中心差商公式(,)(,)2(,)( ,1)(,)( )( ,)(,1)( )( ,1)(,1)()2kjkjkjuu kju kjOtuu kju kjOtuu kju kjOt一维热传导方程)0(
20、022axuatu可分别表示为(), 1(),(2), 1(),() 1,(2hOhjkujkujkuajkujku(), 1(),(2), 1() 1,() 1,(2hOhjkujkujkuajkujku222( ,1)( ,1)(1, )2 ( , )(1, )(2u k ju k ju kju k ju kjaOhh由此得到一维热传导方程的不同差分近似精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 14 页,共 30 页精品资料欢迎下载022,1, 1,1,huuuauujkjkjkjkjk022,1, 11,huuuauujkjkjkjkj
21、k,1,11,1,2202k jk jkjk jkjuuuuuah上述差分方程所用到的节点各不相同。其截断误差分别为)(2hO,)(2hO和)(22hO。因此,它们都与一维热传导方程相容。如果将式22( ,1)( ,1)(1, )2 ( , )(1, )(2u k ju k ju kju k ju kjaOh中的jku,用)(211,1,jkjkuu代替,则可得到又一种差分近似022, 11,1, 11,1,huuuuauujkjkjkjkjkjk差分方程用到四个节点。由Taylor 公式容易得出)()(2121,1,Ouuujkjkjk故其的截断误差为2222)(hOhO。因而不是对任意的精
22、选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 15 页,共 30 页精品资料欢迎下载0, h,此差分方程都能逼近热传导方程)0(022axuatu,仅当( )o h时,才成立。综上可知,用不同的差商公式可以得到微分方程的不同的差分近似。构造差分格式的关键在于使其具有相容性、收敛性和稳定性。前面三个方程都具有相容性,而此方程则要在一定条件下才具有相容性。(二)初、边值条件的处理为用差分方法求解定解问题初值问题xxxuxtxuatu)()0 ,(, 0022初边值问题22120 0, 0(,0)()0(0,)( ),( , )( )0uuatTxltx
23、u xxxlutgtu l tgttT还需对定解条件进行离散化。对初始条件及第一类边界条件,可直接得到kkkxuu)0,(0,(0, 1,0,1, )kkn或精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 16 页,共 30 页精品资料欢迎下载jjjnjjjgtluugtuu2,1,0),(),0()1,1,0(mj其中Tmhln,对第二、三类边界条件)()()()(22101tgutxutgutxulxxTt0需用差分近似。下面介绍两种较简单的处理方法。(1)在左边界)0( x处用向前差商近似偏导数xu,在右边界)(lx处用向后差商近似xu,即)
24、(),1(),()(),0(),1(),(),0(hOhjnujnuxuhOhjujuxujnj),1 ,0(mj则得边界条件的差分近似为精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 17 页,共 30 页精品资料欢迎下载jjnjjnjnjjjjjguhuuguhuu2,2,1,1,01,0,1),1,0(mj其截断误差为)(hO。(2)用中心差商近似xu, 即2(0,)2( , )(1, )( 1, )()2(1, )(1, )()2jn juujujO hxhuu nju njO hxh),1 ,0(mj则得边界条件的差分近似为jjnjjnj
25、njjjjjguhuuguhuu2,2,1,11,01,1,12),1 ,0(mj其截断误差为)(2hO。误差的阶数提高了,但出现定解区域外的节点),1(j和), 1(jn,这就需要将解拓展到定解区域外。可以通过用内节点上的u值插值求出ju,1和jnu,1,也可以假定热传导方程在边界上也成立,将差分方程扩展到边界节点上,由此消去ju,1精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 18 页,共 30 页精品资料欢迎下载和jnu,1。(三)几种常用的差分格式以热传导方程的初边值问题221200, 0( , 0)()0(0,)( ),( , )( )
26、0uuatTxltxu xxxlutgtu l tgttT为例给出几种常用的差分格式。(1)古典显式格式令2har,则022,1,1,1,huuuauujkjkjkjkjk可改写成jkjkjkjkruurruu,1,11,)21(将其与初始条件及第一类边界条件kkkxuu)0,(0,(0, 1,0,1, )kkn或jjjnjjjgtluugtuu2,1,0),(),0()1,1,0(mj结合,我们得到求解此问题的一种差分格式精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 19 页,共 30 页精品资料欢迎下载), 1 , 0(,), 1 , 0(,
27、 1 , 0, 1, 2 , 1()21 (2,1,00, 1, 11,mjgugunkumjnkruurruujjnjjkkjkjkjkjk由于第 0 层)0(j上节点处的u值已知)(0,kku,由此即可算出u在第一层) 1( j上节点处的近似值1,ku。重复使用此式,可以逐层计算出所有的jku,,因此此差分格式称为古典显式格式。又因式中只出现相邻两个时间层的节点,故此式是二层显式格式。(2)古典隐式格式将式022,1, 11,huuuauujkjkjkjkjk整理并与初始条件及第一类边界条件式联立,得差分格式如下,1,1,1,11,1,00,1,2(2)(1,2, ,1,0,(0,1, ,
28、 ),(0,1, , )k jk jkjk jkjkkjjn jjuur uuuknjuknugugjm其中2har。虽然第 0 层上的u值仍为已知,但不能由上式直接计算以上各层节点上的值jku,,必须通过解下列线性方程组精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 20 页,共 30 页精品资料欢迎下载,11,1,11,11,1,11,(2)(1 2 )k jkjk jkjkjk jkjur uuurur uru1,11,2,12,2,12,1,11,12012001212jjjjnjnnjnjrruurgrrruuuuuurgrrrrr)1,
29、1 ,0(mj才能由jku,计算1, jku,故此差分格式称为古典隐式格式。此方程组是三对角方程组,且系数矩阵严格对角占优,故解存在唯一。(3)Richardson 格式Richardson 格式是将式,1,11,1,2202k jk jkjk jkjuuuuuah整理后与初始条件及第一类边界条件式联立。其计算公式为), 1 , 0(,), 1 , 0(, 1, 2, 1()2(22,1,00, 1, 11,1,mjgugunkunkuuuruujjnjjkkjkjkjkjkjk精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 21 页,共 30 页
30、精品资料欢迎下载这种差分格式中所涉及的节点出现在1, 1jjj三层上,故为三层显式格式。 Richardson格式是一种完全不稳定的差分格式,因此它在实际计算中是不能采用的。(4)杜福特 - 弗兰克尔( DoFort-Frankel)格式DoFort-Frankel格式也是三层显式格式,它是由式022, 11,1, 11,1,huuuuauujkjkjkjkjkjk与初始条件及第一类边界条件式结合得到的。具体形式如下:), 1 , 0(,), 1 , 0(, 2 , 1(2121)(2122,1, 00,1, 1, 11,mjgugunkunkurruurrujjnjjkkjkjkjkjk用这
31、种格式求解时,除了第0 层上的值0,ku由初值条件得到,必须先用二层格式求出第1 层上的值1 ,ku,然后再按上式逐层计算), 3, 2(,mjujk。(5)六点隐式格式对二阶中心差商公式222(,)(1,)2(,)(1,)kjuu kju kju kjxh精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 22 页,共 30 页精品资料欢迎下载如果用22xu在点) 1,(jk与点),(jk处的二阶中心差商的平均值近似22xu在)21,(jk处的值,即22,22, 11, 11,1, 1)21,(22huuuuuuxujkjkjkjkjkjk同时tu在
32、点)21,(jk处的值也用中心差商近似,即,1,1(,)2kjkjkjuuux这样又得到热传导方程的一种差分近似2,2(2,11,11,1,12,1,jkjkjkjkjkjkjkuuuuuhauu其截断误差为)(22hO,将上式与初始条件及第一类边界条件式联立并整理,得差分格式,1,1,1,11,11,00,1,2(2)(222(1,2,1,0,1,1)(0,1, ),(0,1,)k jk jkjk jkjkjk jkkkjjn jjrruuuuuuuuknjmuknugugjm精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 23 页,共 30 页
33、精品资料欢迎下载此格式涉及到六个节点,它又是隐式格式,故称为六点隐式格式。与古典隐式格式类似,用六点格式由第j层的值jku,计算第1j层的值1, jku时,需求解三对角方程组1,12,12,11,11/20/21/ 200/ 21/2/21jjnjnjrrurrruuurrrrr)1,1 ,0(mj此方程组的系数矩阵严格对角占优,故仍可用追赶法求解。例 2 用古典显式格式求初边值问题309), 3(,0),0(30)0,(30, 300222ttutuxxxuxtxutu的数值解,取5.0, 1h。解 这 里21,0.5aarh,2)(xx,精选学习资料 - - - - - - - - - 名
34、师归纳总结 - - - - - - -第 24 页,共 30 页精品资料欢迎下载9)(, 0)(21tgtg。由格式), 1 , 0(,), 1 , 0(, 1, 2, 1()21(2,1,00,1, 11,mjgugunkunkruurruujjnjjkkjkjkjkjk可得到)6, 1 , 0(9, 0)3,2, 1 ,0()5, 1 ,0, 2, 1()(5.0, 3,0220, 1, 11,juukkxujkuuujjkkjkjkjk将初值0,ku代入上式,即可算出2) 04(5 . 0)(5 . 00, 00,21 , 1uuu5) 19(5 . 0)(5 . 00, 10, 31
35、,2uuu将边界条件01, 0u,91 , 3u及上述结果代入又可求得9, 5 . 5, 5 . 2, 02, 32, 22, 12, 0uuuu如此逐层计算,得全部节点上的数值解为0,31,32,33,30,40,2.75,5.75,9,0,uuuuu精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 25 页,共 30 页精品资料欢迎下载4 双曲型方程的差分解法对二阶波动方程22222xuatu如果令1uvt,xuv2,则方程可化成一阶线性双曲型方程组xvtvxvatv12221记Tv),(21vv,则方程组可表成矩阵形式xAxatvvv0102矩
36、阵A有两个不同的特征值a,故存在非奇异矩阵P,使得aaPAP001作变换TwwP),(21vw,方程组可化为xtww方程组由二个独立的一阶双曲型方程联立而成。因此本节主要讨论一阶双曲型方程的差分解法。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 26 页,共 30 页精品资料欢迎下载4.1 几种简单的差分格式考虑一阶双曲型方程的初值问题xxxuxtxuatu)()0,(,00将xt平面剖分成矩形网格,取x方向步长为th,方向步长为,网格线为(0, 1, 2,)kxxkh k,)2, 1 , 0( jjttj。为简便,记),(),(jktxjk,)
37、(),(),(kkjkxtxujku。以不同的差商近似偏导数,可以得到方程的不同的差分近似0, 1,1,huuauujkjkjkjk0, 1,1,huuauujkjkjkjk02, 1, 1,1,huuauujkjkjkjk截断误差分别为)(hO,)(hO与)(2hO。结合离散化的初始条件,可以得到几种简单的差分格式精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 27 页,共 30 页精品资料欢迎下载), 2, 1 , 0, 2, 1, 0()(0, 1,1,jkuuuaruukkjkjkjkjk), 2 , 1 , 0, 2, 1, 0()(0,
38、 1,1,jkuuuaruukkjkjkjkjk,1,1,1,0()(0, 1, 2,0,1,2, )2k jk jkjkjkkaruuuukju其中rh。如果已知第j层节点上的值jku,,按上面三种格式就可求出第1j层上的值1,jku。因此,这三种格式都是显式格式。如果对tu采用向后差商,xu采用向前差商,则方程可化成0)(),(), 1() 1,(),(hOhjkujkuajkujku相应的差分格式为精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 28 页,共 30 页精品资料欢迎下载), 2 , 1 , 0, 2, 1, 0()(0,1,1,
39、 1,1,jkuuuaruukkjkjkjkjk此差分格式是一种隐式格式,必须通过解方程组才能由第j层节点上的值jku,求出第1j层节点上的值1,jku。例对初值问题xxxuxtxutu)()0,(, 00其中0002101)(xxxx用 差 分 格 式求 其数 值 解)4,3,2, 1(,jujk, 取21hr。 解 记),2,1,0(kkhxk,由初始条件,2,10021,2,11)(kkkxkk精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 29 页,共 30 页精品资料欢迎下载按差分格式2, 1 ,0, 2, 1,0()(0, 1,1,jkuuuaruukkjkjkjkjk计算公式为jkjkjkuuu, 1,1,2123计算结果略。如果用差分格式, 2 , 1 , 0, 2, 1, 0()(0, 1,1,jkuuuaruukkjkjkjkjk求解,计算公式为)(21,1,1,jkjkjkuuu计算结果见表略。与准确解11(,)200xtuxtxtx比较知,按前一个差分格式所求得的数值解不收敛到初值问题的解,而后一个差分格式的解收敛到原问题的解。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 30 页,共 30 页