1、偏微分方程數(shù)值解的研究摘要：本文從偏微分方程問題的起源著手，介紹從18世紀起，達朗貝爾首先研究弦振動方 程。然后介紹邊界積分方程組的兩類常見的數(shù)值解方法，投影法和機械求積法。最后介紹兩 種非線性偏微分方程的展開解法，Exp-function方法和(G/G)展開法。關(guān)鍵詞：投影法：機械求積法；Exp-function方法；(G /G)展開法；0偏微分方程問題的起源微積分創(chuàng)立不久之后，從18世紀初開始，達朗貝爾(Jean Le Rond dAlembert, 1717 一 1783)首先研究了弦振動方程。傅立葉(J. Fourier, 1768 一 1830)第一個 發(fā)現(xiàn)并求解了熱傳導(dǎo)方程，其平衡

2、態(tài)則用拉普拉斯方程刻畫。這幾類方程是最經(jīng) 典的偏微分方程，由它們而產(chǎn)生大量的求解方法以及通解，但是，不存在求解偏 微分方程的一般方法。雖然19世紀初期柯西(A.L.Cauehy，1789 一 1851)、拉格 朗日(J. L. Lagrange，1736 一 1813)等人解決了一階偏微分方程的求解問題，但 這種化為一階常微分方程組求解的基本方法，對于二階偏微分方程并不適用。從數(shù)學(xué)研究來講，數(shù)學(xué)家們一直相信通解的存在，并且以此為出發(fā)點，先求 通解，最后確定常數(shù)和函數(shù)。1820年柯西首先證明了一階常微分方程初值問題 解的存在性和唯一性。而且很快從實數(shù)推廣到復(fù)數(shù)范圍了，對于偏微分方程他只 需要證明

3、一階偏微分方程組在復(fù)數(shù)范圍內(nèi)的解的存在性。正是在這里，他創(chuàng)造性 地發(fā)明了優(yōu)函數(shù)方法：首先應(yīng)用幕級數(shù)展開式給出形式解，然后與某個己知收斂 的幕級數(shù)相比較來證明形式解的收斂性，從而證明了解析解的存在性。30多年之后，魏爾斯托拉斯(K. WeierstraSS，1815 一 1897)的學(xué)生俄國數(shù)學(xué) 家科瓦列夫斯卡婭獨立地證明了偏微分方程組柯西問題解析解的存在唯一性，即 柯西一科瓦列夫斯卡婭定理。柯西一科瓦列夫斯卡婭定理是偏微分方程理論中第 一個普遍的存在定理。這樣普遍適用的、不依賴于方程類型的性質(zhì)的一直引人關(guān) 注，是偏微分方程理論的一個主要發(fā)展方向；但另一方面，18世紀以來對幾類 特殊方程的研究結(jié)

4、果表明，把微分方程分成不同的類型：雙曲型、拋物型和橢圓 型方程。這些是二階線性偏微分方程的三種基本類型，數(shù)學(xué)物理中常見的波動方 程、熱傳導(dǎo)方程及調(diào)和方程分別為它們的代表，其定解問題的提法、解的性質(zhì)以 及求解方法，多數(shù)可以推廣到創(chuàng)門所代表的三種類型上去。有些方程在區(qū)域的一部分可能是雙曲型的，而在另外某一部分可能是橢圓型 的，在分界部分或者退化為拋物型的或者是不確定的，這樣的方程稱為混合型。 1923年F.G特里科米(F. Tricomi)首先研究了這樣的方程，并得到了深刻的結(jié)果。關(guān)于各種各樣的定解問題的提法，既有其物理上的依據(jù)，在數(shù)學(xué)上又能按照 一定的準則被證明是合理的，這些準則中最常用的是解的

5、存在性、唯一性以及解 關(guān)于定解數(shù)據(jù)的連續(xù)依賴性。如果一個定解問題的解是存在、唯一且連續(xù)地依賴 于定解條件，則這個定解問題就稱為適定的，此概念是阿達瑪(Jacuqes Hdamaard，1565 一 1963)提出來的。雙曲型方程柯西問題的現(xiàn)代理論,是由阿達瑪對二階雙曲型方程柯西問題的 先驅(qū)工作開始的。他通過構(gòu)造在特征劈錐面上具有奇性的解-基本解來求解柯 西問題，并采用發(fā)散積分的有限部分的方法來克服所遇到的奇性困難。他的工作 經(jīng)過里斯及索伯烈夫等人的發(fā)展，對廣義函數(shù)論的建立是一個重要的推動，而阿 達瑪?shù)姆椒ㄔ趶V義函數(shù)論的框架中也得到了更清晰和完善的表達。黎曼(GF.B.Rimenan，1826

6、一 1866)依據(jù)所謂的狄里克雷原理斷言狄里克雷 問題有解，魏爾斯托拉斯指出了黎曼這個論斷的邏輯缺陷，并舉出了反例。多年 之后，希爾伯特(D.Hilbert，1862 一 1943)給出了狄里克雷原理的完整無缺的證明。 (狄里克雷問題解的存在性證明，有龐加萊一佩隆掃除法，施瓦茲交替法，差分 法等等。)這方面的研究極大地推動了泛函分析的發(fā)展，也使得變分法成為研究 偏微分方程的強有力的工具。拉普拉斯方程的狄里克雷問題解的存在性，波動方程的特征初值問題解的存 在性等具體方程的一般理論問題，在19世紀也同樣受到關(guān)注并分別推廣到更為 一般的一類方程。20世紀30年代彼得羅夫斯基(1.G.Petrovky

7、，1901 一 1973)的工作成為偏微 分方程理論發(fā)展的又一個里程碑，他對古典分型理論進行了推廣，并同時研究了 每一類方程的特點。1邊界積分方程數(shù)值解在積分方程的解析解很難求或無法求出時，借助數(shù)值方法求解它的近似解就 顯得非常重要。從已有的成果來看，方法的應(yīng)用先于理論分析，也就是說，不少 方法已被提出或應(yīng)用，但其數(shù)學(xué)基礎(chǔ)的研究則不甚深入。隨著應(yīng)用的愈高要求， 近年人們已把著眼點轉(zhuǎn)向方法的數(shù)學(xué)原理研究上。我國對奇異積分方程的數(shù)值 解法的研究起步較晚，但自82年全國積分方程會議強調(diào)開展這方面的研究后， 很快引起了一些學(xué)者的興趣和關(guān)注，特別是數(shù)值方法的理論分析方面，近年來我 國學(xué)者做了不少工作。這

8、里僅考慮方程的常見數(shù)值解法。1.1投影法常見的投影法通常有配置法，Galerkin有限元法和最小二乘法。設(shè)中.(j = 1,2,.,n)是化的基函數(shù)，近似解v ex可表示為這組基函數(shù)的線性組合 hh hv =n c 中(x)( 1)j=1投影法都要歸結(jié)為求解以c, = (j = 1,2,., n)為未知量的線性方程組。配置法：設(shè)e D是一組適當(dāng)選取的配置點(數(shù)目為n個)，令A(yù)v - f在配置點上恒等于零就得到n個線性方程工Ap (t )c = f (t),i = 1,2,.,n(2)j i j ij=1一般說來，系數(shù)矩陣A%(t)是不對稱的。配置點的選取要取決于積分算子的性 質(zhì)和邊界的幾何狀態(tài)

9、，它對于離散方程(2)的適定性有決定性的影響。Galerkin方法：在這種方法中，基函數(shù)系具有有限支集，系數(shù)七是通過下 面的Galerkin有限方程得到的工(Ap ,甲)c = (f，甲),i = 1,2,.,n(3)j i jij=1如果A是自伴的，則上面方程組的系數(shù)矩陣是對稱的。最小二乘法：這種方法要求| Av - f|取得最小值。如果方程的解存在唯12(. D)一，則最小二乘法將導(dǎo)致一個對稱，正定的線性方程組工(Ap.俱 M = (f, Api), i = 1,2,., nj=1以上方法都可以看作Galerkin-Petrov方法的特例。它們采用相同的試探函數(shù)，但 選擇不同的檢驗函數(shù)。所

10、得到的線性代數(shù)方程組的系數(shù)矩陣都是滿陣。1.2機械求積法選取一個適當(dāng)?shù)那蠓e公式來近似積分算子Av有A (s) = hLk(k,t )v(t ),t = jh,h = (b一a)/n, j = 1,2,n(4)nj j jj=1然后賦值就可以求出離散矩陣的每一個元素，離散方程為hk(s ,t )v(t ) = f (s ),s = ih,i = 1,2,n(5)i j j i ij=1機械求積法對連續(xù)核可以借助Anselone聚緊收斂理論得到它的收斂性證明，但 是對邊界積分方程由于核的奇性，則必須另辟蹊徑。2兩種非線性偏微分方程的展開解法2.1 Exp-function 方法描寫淺水波的KdV方

11、程，它的一般形式如下：u + uux +=0KdV方程應(yīng)用非常廣泛，它也是許多領(lǐng)域的孤波現(xiàn)象的模型：如等離子體 聲波，彈性桿中縱向色散波，低溫下非線性晶格的聲子波包的熱激發(fā)等。我們考慮如下的KdV方程：u + 6uu +uu = 0利用變換u = u（&）,& = kx + wt代入上述方程，則方程轉(zhuǎn)化為常微分方程：wu + 6kuu + k 3u= 0。現(xiàn)假設(shè)方程的解具有如下形式:d a exp（詣）f pb exp（用）我們平衡最高階線性項u”以及最高階非線性項uu:c exp(d + 7q)& + + c exp(-c-7p)&=q=一 七涼用)8c exp(2d + q)& + c e

12、xp(- 2c- p)& uu = 3寸4乙q b.exp( j& )3于是，我們有：d + 7q = 2d + 6q以及-c-7p = - 2c-6p從而我們得到c,d,p,q的關(guān)系為：p = c,q = d。取 p = c = 1, q = d = 1則上式為 u = aiexp(&) + aiexp(&) exp(&) + b + b exp(-&) TOC o 1-5 h z 0-1帶回原方程中我們得到exp(-2&) + c exp(-3&) = 03c exp(3&) + c exp(2&) + c exp(&) + c + c exp(-&) + c3210-1-2A = (ex

13、p(&) + b + b exp(-& )4其中0-1c = -wa + 6ka2b + wab k3a b k3a 6ka a301 01 01 001 0c =-6ka2-2wa -8k3a + 2wab - 12kaa+ 8k3ab + 2wab2 -2wa b -4k3a b + 6ka2b2 + 12ka2b1 -11 00 00 01 01 -1c = -18k3a bb + 6wa b b + 18ka2b b + wa b3 + wa b2 6ka2b + k3a b311 0 -11 0 -11 0 -11 00 00 01 0一 k3a b2 + 6ka ab2 18ka

14、a 一 wa b 一 23k3a b 一 5k3a b 一 5wa b +12 0 00 1 00-10-10-1-10-10ka a b 一 12ka a b1 0 -11 -1 0c = -4wa b2 +12ka2b2 + 32k3a b2 - 4k3a b2 -12ka2 - 4wa b-24ka a b + 24ka abb + 4wa b 2 + 4wa b?b + 4k 3a b2b0-1010-101 -11 0 -11 0 -1c = 6ka2b b + 18kaa b2 12ka a b + 5wab b + 12ka ab b wa b3 k3a b3 + wa b b2

15、 6ka a b2-10 011 0101 11 0111101 01 001 00 1 0 TOC o 1-5 h z + wa b2 + 5k3abb2 18ka2b 23k3a b 6wa b b + 18k3a b b kaab 12kaa b 0 110 110011011 011 01110c = 2wa b2b + 4k3a b2b 4k3a b2 b + 2wa b3 8k3a b2 6ka2b221011010 101111016ka2b2 2wa b2 12ka2b + 8k3a b3 + 2wa b2b + 12ka ab1 011111 10 1 0111c = 6ka ab + k 3a b3 k 3a b b 2 + wa b3 wa bb2 6ka 2 b b310 10 11 010 11 0110 1上述式子中假若系數(shù)為0,得到如下的代數(shù)方程：C3= 0，C2= 0，0 0時:C 2 11 ; 2a12au &) = a c1sinhM 立肉+c2coshn 甲肉 c cosh(1肉+ c sinh1var肉12 c