研究有限差分格式稳定性的其他方法-报告

2015年秋季学期研究生课程考核（读书报告、研究报告）考核科目:偏微分方程数值解法学生所在院（系）:理学院数学系学生所在学科:数学学生姓名:Hiter学号:1XS012000学生类别:考核结果阅卷人第1页共9页研究有限差分格式稳定性的其他方法摘要偏微分方程的求解一直是大家比较关心的一个问题，而有限差分格式则是求解偏微分方程时常用并且有效的一个方法。因此，研究有限差分格式的性质就显得尤为重要。在课上我们已经跟着老师学习了运用Fourier方法研究有限差分格式的稳定性，但是在很多研究有限差分格式稳定性的问题中仅仅会用Fourier方法是不够的，所以在本篇论文中，将会介绍其他三种常用的研究有限差分格式稳定性的方法，分别是：Hirt启示型方法、直接方法（或称矩阵方法）和能量不等式方法。关键字：偏微分方程；有限差分格式；稳定性AbstractThesolutionofpartialdifferentialequationshasbeenmoreconcernedwithaproblem,andthefinitedifferenceschemeisacommonandeffectivemethodforsolvingpartialdifferentialequations.Therefore,itisveryimportanttostudythecharacterofthefinitedifferencescheme.WehavefollowedtheteachertolearntheuseofFouriermethodoffinitedifferenceschemestability,butinalotofresearchonthestabilityoffinitedifferenceschemeisonlyusedFouriermethodisnotenough,sointhispaper,willintroducetheotherthreekindsofcommonlyusedinthestudyoffinitedifferenceschemestabilitymethod,respectivelyis:Hirtenlightenmentmethod,directmethod(ormatrixmethod)andenergyinequalitymethod.Keywords:partialdifferentialequation;finitedifferencescheme;stability1前言微分方程的定解问题就是在满足某些定解条件下求微分方程的解。在空间区域的边界上要满足的定解条件称为边值条件。如果问题与时间有关，在初始时刻所要满足的定解条件，称为初值条件。不含时间而只带边值条件的定解问题，称为边值问题。与时间有关而只带初值条件的定解问题，称为初值问题。同时带有两种定解条件的问题，称为初值边值混合问题。定解问题往往不具有解析解，或者其解析解不易计算。所以要采用可行的数值解法。有限差分方法就是一种数值解法，它的基本思想是先把问题的定义域进行网格剖分，然后在网格点上，按适当的数值微分公式把定解问题中的微商换成差商，从而把原问题离散化为差分格式，进而求出数值解。此外，还要研究差分格式的解的存在性和唯一性、解的求法、解法的数值稳定性、差分格式的解与原定解问题的真解的误差估计、差分格式的解当网格大小趋于零时是否趋于真解（即收敛性），等等。有限差分方法具有简单、灵活以及通用性强等特点，容易在计算机上实现。在课上我们已经跟着老师学习了运用Fourier方法研究有限差分格式的稳定性，但是在很多研究有限差分格式稳定性的问题中仅仅会用Fourier方法是不够的，所以在本篇论文中，将会介绍其他三种常用的研究有限差分格式稳定性的方法，分别是：Hirt启示型方法、直接方法和能量不等式方法。2Hirt启示性方法2.1方法概述Hirt启示性方法是一种近似分析方法。主要是把差分格式在某确定点上作泰勒级数近似第2页共9页展开，把高阶误差略去，只留下最低阶的误差项。如果差分格式是相容的，那么这样得到的新的微分方程（称之为第一微分近似或修正微分方程）与原来的微分方程相比只增加了一些含小参数的较高阶导数的附加项。Hirt方法就是利用第一微分近似的适应性来研究差分格式的稳定性。Hirt方法的判别准则是这样的：如果第一微分近似是适定的，那么原来微分方程的差分格式是稳定的，否则不稳定。其实所述的微分格式是原来微分方程问题的相容的差分格式，那么也可以看作第一微分近似问题的相容的差分格式。如果第一微分近似问题是不适定的，那么它的差分格式将不稳定[1]。2.2操作方法先给出几个方程0,,0,0tRxaxuatu（2.1）,2,1,0,,2,1,0,011njhuuauunjnjnjnj（2.2）011huuauunjnjnjnj（2.3）考虑对流方程（2.1）的差分格式（2.3），在点），（njtx进行Taylor技术展开，有)(][2][),(),(2221hOxuhxuhtuutxunjnjnjnj)(][2][),(),(2221Otuhtutuutxunjnjnjnj利用对流方程（2.1），有22222)(xuaxuattu因此，在点），（njtx上，有差分方程（2.3）可以得到)(2222222hOxuaahxuatu）（略去高阶误差项，得出第一微分方程近似22222xuaahxuatu）（要使上面的抛物型方程有意义，必须有0222aah而上面的不等号改为等号，则就化为原来的对流方程。在这两种情况下，相应的问题是适定的。即第一微分近似适定的条件是第3页共9页0222aah由此得出差分格式（2.3）的稳定性条件是1a，其中h。此结论与Fourier方法分析得到的结论是一致的。下面我们再来分析逼近对流方程（2.1）（仍设0a）的差分格式（2.2）的稳定性。模仿上面的推导可以得到它的第一微分近似是22222xuaahxuatu）（可以看出22xu的系数小于0，因此第一微分近似是不适定的，从而推出差分格式（2.2）是不稳定的。3直接方法关于抛物型方程初值问题的差分格式的稳定性问题，可以用直接方法（或称矩阵方法）来研究。下面用具体例子来说明这个方法的基本思想及使用方法。考虑常系数扩散方程的初值问题0,0),(),0(),0(),()0,(0),,0(,0,022ttlutulxxuxutlxaxuatu（3.1）采用显示差分格式来逼近，即0,01,2,1),(1,2,1,0,200002111nuuJjxuuJjnhuuuauuJnjjnjnjnjnjnj（3.2）其中lJh。先把差分格式（3.2）写成1,2,1,)21(111Jjuauauaunjnjnjnj（3.3）其中2h。可以把（3.3）写成向量形式，即]00[]][21212121[][00122111121211JnnJnJnnnJnJnnuuauuuuaaaaaaaaaauuuu（3.4）如果令第4页共9页TnJnnnuuuu),,,(121并考虑到00nJnuu，则（3.4）式可以写成nnAuu1（3.5）其中]21212121[aaaaaaaaaaA（3.6）从显示格式出发，得到方程组（3.5）式，也可以理解为较为一般的形式，即对于逼近初值问题（3.2）的其他二层格式也可以化为（3.5）式的形式。当然此时A不是（3.6）式所表示的形式。如果差分格式是二层隐式格式。则A为CB1这种形式。因此（3.5）式这种形式可理解为既包含二层显示格式又包含二层隐士格式的较为一般的形式。引入误差向量~nnnuuz，其中nu是差分方程（3.5）的精确值（理论值），~nu是差分方程（3.5）经数值求解得到的值（包括了舍入误差等）。显然，nz满足nnAzz1（3.7）从而推出0zAznn（3.8）差分格式（3.5）的稳定性就要求0,nKzn（3.9）其中为向量的2-范数。由于02zAznn因此（3.9）式成立的充分必要条件为MAn2（3.10）上述采用2-范数，当然也可以采用其他类型的范数。对于稳定性条件（3.10），可以仿Fourier方法中的推导，得到一些结论：（1）谱半径条件MA1)(（3.11）是差分格式稳定的一个必要条件，其中M为常数。（2）如果矩阵A是一个正规矩阵，则（3.11）式也是格式稳定的一个充分条件。下面讨论差分格式（3.5），（3.6）的稳定性。矩阵（3.6）是对称矩阵，所以只要使条件（3.11）第5页共9页成立即可。现在来计算A的特征值。令)1(J阶方阵]0110110110[S则A可以表示为SaIaA)21(其中I为)1(J阶单位矩阵。由此可知，关键是求出S的特征值和特征向量。设和TJ),,,(121分别为S的特征值和特征向量，wSw写成分量的形式有02,,1,0,0021JjjjwwJj（3.12）先求出jw，再求出S的特征值。由于S为对称矩阵，所以其特征值为实数。由Gerschgorin定理知，1Jkjkjkkss其中kjs为矩阵S的元素。由此得到2。（3.12）式的第一式为常系数线性差分方程。设其解具有如下形式：0,jjw将它代入（3.12）式的第一式，便得到关于的一元二次方程012此方程称为（3.12）式的第一式的特征方程。由于2，所以其解为2)2(12i其中1i。可以看到1)2(1)2(222第6页共9页取2)2(1sin,2cos，则ie。因此差分方程（3.12）的解可以表示为Jjeaeawijijj,,1,0,21由00w，得到021aa。再由0Jw，得到021iJiJeaea，从而有02）（iJiJeea由此可推0sin2Ja。02a，有1,2,1,JkkJ。所以得到Jk，可以得到Jkkcos2。注意到Jh1，则S的特征值为khkcos2。从而得到A的特征值为1,2,1,2sin41cos2212Jkkhakhaak当21a时，1)(A。因此显示格式的稳定性条件为21a。下面讨论隐式格式0)(20002111111nJnjjnjnjnjnjnjuuxuuhuuuauu的稳定性。可以把隐式格式写成向量形式nnuBu11其中TnJnnnuuuu),,,(121，SaIaB)21(。利用前面已经求得的S的特征值，可以得到B的特征值1,,2,1),cos1(21)(21)(JjkhaSaaBkk由此可知，1)(Bk，从而有1)(1Bk。注意的B为对称矩阵，所以1B也为对称矩阵，利用直接方法结论（2）知，扩散方程隐式格式是无条件稳定的。从上面的叙述看来，利用直接方法来分析抛物型方程的初值问题的差分格式并不困难。但在实际应用中却存在着一定的限制。上面讨论稳定性的两个例子中式依据了特殊矩阵S才求出了)1(J阶矩阵A、1B的特征值。一般说来，计算高阶矩阵的特征值是相当困难的，因此直接方法应用也就很困难了。第7页共9页4能量不等式方法4.1方法概述在讨论线性常系数差分格式的稳定性问题时，建立了判别差分格式的稳定性准则