收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

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

张然  
【摘要】:物理学、力学和工程技术方面的许多问题都可归结为按初始条件和边界条件求解偏微分方程的初值-边值问题。这些问题中能采用解析法按照边值条件求解的仅限于极少数情况。所以,一般只能采用近似方法求解。近些年来,随着电子计算机的广泛应用,数值解法逐渐成为解初值-边值问题的一种非常有效的方法,在数值分析中占有重要的地位。 自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) 为解该线性方程组,我们记 、飞?


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 齐东旭;样条函数磨光法及其应用[J];吉林大学学报(理学版);1977年02期
2 张家驹;用样条函数解两类插值问题[J];计算数学;1984年01期
3 何天晓;朱功勤;;第三型三角剖分下具有最小矩形局部支集的样条函数[J];淮北煤炭师范学院学报(自然科学版);1984年Z1期
4 张其林,沈祖炎;非理想压杆稳定分析的样条函数法[J];应用力学学报;1994年02期
5 于巍;;样条函数的B-样条基和截断幂基表示之间的转换[J];高师理科学刊;2008年06期
6 张辉;陈建军;;利用样条函数导出数值积分公式[J];晋中学院学报;2010年03期
7 李岳生;样条函数的共轭插值(Ⅱ)——微分算子样条[J];中山大学学报(自然科学版);1981年01期
8 李大成;应用样条函数积分计算矿体储量的探讨[J];地质科技情报;1982年S1期
9 孙永生;一个广义样条函数类上的极值问题[J];科学通报;1982年11期
10 王省富;刘克轩;;几类插值有理样条函数[J];西北工业大学学报;1988年01期
11 陈发来;Schoenberg变差缩减算子的迭代极限[J];中国科学技术大学学报;1992年02期
12 毛瑞祥;刘中林;;斜板的双样条最小二乘配点法静力研究[J];长安大学学报(自然科学版);1993年03期
13 杨少伟,张乃苍;样条函数平曲线加宽的研究[J];重庆交通学院学报;1994年03期
14 丛爽,宋瑞祥;改进型B样条网络函数及其应用[J];中国科学技术大学学报;2000年06期
15 张伟林,牛忠荣;样条函数线法在壳体非线性分析中的应用[J];安徽建筑工业学院学报(自然科学版);2002年01期
16 何启志;指数样条函数在零息票债券定价中的应用[J];合肥工业大学学报(自然科学版);2005年05期
17 程正兴;一类新的光顺问题[J];上海第二工业大学学报;1986年01期
18 王小林;复插值样条函数及其收敛性[J];数学物理学报;1987年02期
19 段奇;李筛和;;一类带边界混合约束的拟半节点二元样条S_2~1(Δ_(mn)~(2))插值[J];山东大学学报(工学版);1989年02期
20 胡永谟;样条函数中的多项式Φ_n(λ,z)的一个应用[J];大学数学;1992年04期
中国重要会议论文全文数据库 前10条
1 赵建华;;局部坐标系下的代数与正弦混合样条[A];数学及其应用文集——中南模糊数学和系统分会第三届年会论文集(下卷)[C];1995年
2 冯华仲;莘仁俊;;样条函数在燃烧室出口温度数据处理中的应用[A];中国航空学会第七届动力年会论文摘要集[C];2010年
3 王文江;郑会永;;光顺样条滤波的数学原理及其应用[A];1996年中国控制会议论文集[C];1996年
4 武兰河;王道斌;;层合板壳稳定问题的样条函数解法[A];第七届全国结构工程学术会议论文集(第Ⅰ卷)[C];1998年
5 杨绿峰;赵艳林;李桂青;;闭口薄壁杆件约束扭转计算的样条里兹法[A];第八届全国结构工程学术会议论文集(第Ⅰ卷)[C];1999年
6 吴强;;GSM(1,1)模型的结构及数值解[A];数学及其应用文集——中南模糊数学和系统分会第三届年会论文集(上卷)[C];1995年
7 祁守宝;曹庆;;非圆曲线的线切割自动编程软件开发[A];第八届全国电加工学术年会论文集[C];1997年
8 王复明;刘文廷;许学标;;有限深弹性空间地基的半数值解[A];第四届全国结构工程学术会议论文集(下)[C];1995年
9 罗永峰;罗小华;陈扬骥;;一种改进的非线性旋转薄壳单元理论[A];庆贺刘锡良教授执教五十周年暨第一届全国现代结构工程学术报告会论文集[C];2001年
10 闫庐松;佟萌萌;张洋洋;季春霖;杨琳;姚广仁;;非线性自校正技术及其在涂层测厚仪设计中的应用[A];2004中国控制与决策学术年会论文集[C];2004年
中国博士学位论文全文数据库 前10条
1 邵新平;非线性微分方程的样条函数求解方法[D];浙江大学;2010年
2 张然;样条函数与小波函数在偏微分方程数值解中的应用[D];吉林大学;2004年
3 徐敏;样条函数与径向基函数逼近某些问题研究[D];大连理工大学;2010年
4 王彦博;数值微分及其应用[D];复旦大学;2005年
5 徐惠霞;B样条多重乘积理论与有理曲线曲面多项式逼近技术的研究[D];浙江大学;2008年
6 苏成;分域样条虚边界元法及其在高层建筑结构分析中的应用[D];华南理工大学;1997年
7 姜自武;样条函数与径向基函数的若干研究[D];大连理工大学;2010年
8 郭微光;机载超宽带合成孔径雷达运动补偿技术研究[D];中国人民解放军国防科学技术大学;2003年
9 常锦才;多元样条的力学意义与自适应数据分析的某些研究[D];大连理工大学;2008年
10 姜政毅;多元可加细函数及其应用[D];吉林大学;2009年
中国硕士学位论文全文数据库 前10条
1 廖肇源;样条函数在微分方程数值解中的若干应用[D];大连理工大学;2006年
2 李辉;三次样条插指函数的研究[D];北京交通大学;2008年
3 苏延辉;球面和环面上的样条函数[D];大连理工大学;2006年
4 孙增献;估计利率期限结构[D];中国科学技术大学;2009年
5 陈光;服装CAD[D];吉林大学;2005年
6 郭新伟;分布式捷联姿态基准中局域网数据传输及变形测量研究[D];哈尔滨工程大学;2005年
7 张琳;单锚腿系泊系统数值模拟[D];清华大学;2005年
8 张渔勇;桥梁结构列车多轴移动荷载识别研究[D];武汉理工大学;2007年
9 李志明;有理三次样条及空间闭曲线插值问题的研究[D];合肥工业大学;2008年
10 李崇民;线性微分方程边值问题数值求解的多步差分法[D];兰州大学;2008年
中国重要报纸全文数据库 前1条
1 ;李春祥:提出固有频率五种新模型[N];科技日报;2003年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978