样条函数与小波函数在偏微分方程数值解中的应用

张然  

【摘要】：物理学、力学和工程技术方面的许多问题都可归结为按初始条件和边界条件求解偏微分方程的初值-边值问题。这些问题中能采用解析法按照边值条件求解的仅限于极少数情况。所以，一般只能采用近似方法求解。近些年来，随着电子计算机的广泛应用，数值解法逐渐成为解初值-边值问题的一种非常有效的方法，在数值分析中占有重要的地位。 自1946年schoenberg提出样条函数的概念以来，样条函数方法得到了迅速的发展和广泛的应用，同时也成为了函数逼近理论中重要的研究对象。上世纪80年代，一批数学家、物理学家、地理学家为小波概念奠定了坚实的基础，形成了小波理论。小波函数在微分、积分方程数值解，数字信号(图象)处理等诸多领域均有应用。小波由于同时具有空间域与频率域的局部性，因此是描述、检取函数奇性的有效工具。样条与小波的思想对数值数学的发展有很大影响。 本文将进一步探讨样条函数与小波函数在偏微分方程数值解中的应用。主要做了以下几项工作： 1．在第二章中，我们应用一种边界型求积技术，推导出用于求解偏微分方程边值问题的一类边界元方案。作为数值解法，边界元法只需将区域的边界分割成边界单元，因此使所考虑问题的维数降低了一维。与对整个区域进行分割的有限元法和有限差分法相比，它具有输入数据少，计算时间短，特别适用于无限域问题和三维问题等优点。此外，边界元法只对边界离散，离散化误差仅来源于边界，区域内的有关物理量可由解析式的离散形式直接求得，因此提高了计算精度；求解时要改变内点的数量或位置也非常方便；对于只需求出边界值的问题，关于区域内的物理量可以不必进行计算，能提高计算效率。 设Ω为R~n中的有界闭区域，定义二阶微分算子L为 Lu=sum from i,j=1 to n(a_(ij)(X)(?)~2u/(?)x_i(?)x_j)+sum from i=1 to n(b_i(X)(?)u/(?)x_i+c(X)u,) (1) 、‘、、·‘rJ亏、 其中a:J(X)任万老(n):b:(X)，e(X)‘万二(‘2);且毋脾)(0全一)表示几上的函数集合， 集合中每个元素f(x)=f(x，，…，x。)都具有连续的偏导数D(il，…“)f，。三:、三。，k= 1.…，n· L的伴随算子M定义为 、艺闺 几 “。一艺 1，J=1 aZ(。a‘，(X)) 口x:口工， a(。b、(X)) axl +e(X)v. (2) 考虑下面的边值问题 { Lu(X)一夕(X)， X‘n， X‘a几一， X〔a几2， (3) 一一 一叭如一撇 二= UQ二 其中微分算子L由(l)式定义，并且口只、U口几2=a几 在互2.4，我们利用降维展开方法推导出一类用于处理边值问题(3)的边界元法.主 要结果如下: 定理1设几〔R”是。维有界闭区城，其边界曲面ag分片光滑，特别，当，‘=2 时，机2是一条有限长简单闭曲线.令。二试X)任CZ(几)，且满足L。=o，又令”二v(X) 为材。一占(X一XO)的解，其中X。为几中任意给定点.互为伴随的徽分耳子L与AI 的定义分别由(l)和(2)给出. (i)如果X。任几\口几，那么边值问题(3)在点Xo处的解为 一卜五。}喀一、}ds一五。…·客(鑫一(·)鬓)会ld·，(4) 其中。:=;:(X)二艺界，a!，(X)口。/ax，+艺界，vaa:，(X)/ax，一b:(X)v· (11)、于x0。。。，若;i。。n一亡，:(:)擎一。，。且11二。(。):一1一。，其中华 艺一U置口几‘一u口1‘ 为x:在曲面{X;{X}二时上的外法方向导数，则 C“”“。，一五。 。艺;: 瓮]d一五。}·睿(睿一爵)佘ld一、 其中口 一，二关.do，而d“”凡””应于X0”立体角·这里K· =几自口B。厂B‘= {X:!X一X01时是以X0为中心的小球. 下面考虑外边值问题，即当几的补卯=R”\几有界时，几上的边值问题.为得到问 题(3)的唯一解，我们还需要诚x)满足辐射条件 二(X)}}，。(X)1}二。(X)1}v。(X)}[二o(}X}一”)(}X卜co)， (6) 其中v(X)是Mv(X)=0的基本解. 类似于定理1，对外边值问题(3)加果诚X)满足辐射条件(6)，我们可以得到相应的 定理(见39页).之后，我们以平面调和方程、高维调和方程和Helmhol比方程为例给出 了定理l与定理2的应用. 在互2.5中，我们以求解二维Helmho一tz方程(即令(3)中L二△+儿2)的边值为例， 利用样条插值和边界型求积公式(例如[s]中给出的公式)把前面推导的边界积分方程离 散成代数方程，给出了边界元法的具体求解过程. 假定{xj}孔，为边界a。的一个给定剖分·下面，我们来求解边界上所有未知的 二;二。(凡)和勺=a。(xj)/on.由(5)式，我们知道对点凡〔a几，j=1，…，n 异·，一五。·(X)会d“一五。粉 口、(X) a牡 ds， (7) 其中马 的解. 一，二天;d“，“d”“K·的对应于凡的立体”，、为方程‘“十“’)· 二占(X一价) 为将(7)离散化，我们将二和q作插值展开 二(X)会艺。‘叻{‘’(X) 己之1 其中《‘’和《2，为Lagrange插值基函数. 如下代数方程 。(X)‘艺。功{，，(X)，(8) 己=l 将(s)代入(7)，我们就将积分方程(7)离散为 口，.声 二了一Uj+夕U更 Z7r碑‘.口“ 五。，，1)(X)豁d“一客。!瓜‘，2’‘X，、d“，，- 1，…，倪·(9) 为解该线性方程组，我们记 、飞

【关键词】：
【学位授予单位】：吉林大学
【学位级别】：博士
【学位授予年份】：2004
【分类号】：O241
【DOI】：CNKI:CDMD:1.2004.100282
【目录】：

全文下载： CAJ格式
不支持迅雷等加速下载工具，请取消加速工具后下载

CAJViewer7.0阅读器支持所有CNKI文件格式，AdobeReader仅支持PDF格式