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

来源 :吉林大学 | 被引量 : 0次 | 上传用户:zyf20011027
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
物理学、力学和工程技术方面的许多问题都可归结为按初始条件和边界条件求解偏微分方程的初值-边值问题。这些问题中能采用解析法按照边值条件求解的仅限于极少数情况。所以,一般只能采用近似方法求解。近些年来,随着电子计算机的广泛应用,数值解法逐渐成为解初值-边值问题的一种非常有效的方法,在数值分析中占有重要的地位。 自1946年schoenberg提出样条函数的概念以来,样条函数方法得到了迅速的发展和广泛的应用,同时也成为了函数逼近理论中重要的研究对象。上世纪80年代,一批数学家、物理学家、地理学家为小波概念奠定了坚实的基础,形成了小波理论。小波函数在微分、积分方程数值解,数字信号(图象)处理等诸多领域均有应用。小波由于同时具有空间域与频率域的局部性,因此是描述、检取函数奇性的有效工具。样条与小波的思想对数值数学的发展有很大影响。 本文将进一步探讨样条函数与小波函数在偏微分方程数值解中的应用。主要做了以下几项工作: 1.在第二章中,我们应用一种边界型求积技术,推导出用于求解偏微分方程边值问题的一类边界元方案。作为数值解法,边界元法只需将区域的边界分割成边界单元,因此使所考虑问题的维数降低了一维。与对整个区域进行分割的有限元法和有限差分法相比,它具有输入数据少,计算时间短,特别适用于无限域问题和三维问题等优点。此外,边界元法只对边界离散,离散化误差仅来源于边界,区域内的有关物理量可由解析式的离散形式直接求得,因此提高了计算精度;求解时要改变内点的数量或位置也非常方便;对于只需求出边界值的问题,关于区域内的物理量可以不必进行计算,能提高计算效率。 设Ω为Rn中的有界闭区域,定义二阶微分算子L为 Lu=sum from i,j=1 to n(aij(X)(?)2u/(?)xi(?)xj)+sum from i=1 to n(bi(X)(?)u/(?)xi+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=1aZ(。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)为解该线性方程组,我们记、飞?
其他文献
为实现电力资源在电力市场下的有效配置和社会福利的最大化,主要从电价、机制和合约3个角度为电力市场下可中断负荷管理提供解决方法。根据可避免成本设计可中断电价,运用机
朗读是阅读学习常用方法之一,在小学语文阅读教学过程中占有举足轻重的地位。现实中,小学语文课堂中的朗读教学在重视程度、时间安排和评价考核方面均有不同程度的缺失,朗读
对GNSS整周模糊度估计的LAMBDA算法做了简要介绍,分析搜索空间大小的合理确定跟解算的可靠性以及模糊度搜索效率之间的关系。通过实验证明常被采用的两种搜索空间确定方法的
近些年来,随着国民经济快速的发展,国家的能源消费量也在高速增长,但以煤炭、石油为主的传统化石能源由于存在高污染和利用率低等问题,对环境造成破坏的同时,还造成了巨大的能源浪
对现阶段普通高校田径课开设现状及体育教育专业学生的专业选项情况进行调查分析,提出高校体育田径教学改革的对策与建议。
现今世界,国家与国家之间的联系愈来愈紧密,无论是经济还是文化方面,各国人们都开始互相学习,互相交流,学习外语已变成一个见怪不怪的现象了。而在学习一门语言过程中,语言文
汽车驾驶员驾驶行为差错,是汽车行驶过程中普遍存在的一种客观现象。统计表明,汽车驾驶行为差错导致的交通事故占事故总数的85%~92%;相关统计资料也表明,驾驶员因素导致的交通事故
子宫内膜炎(Endometritis)是奶牛常见疾病,是奶牛不孕不育的主要原因。奶牛患子宫内膜炎导致人工授精次数增加,配种和受孕间隔延长,妊娠率、受孕率降低等,使奶牛淘汰率提高,造成养牛
本论文将以城市形象的相关理论为指导,通过对在城市形象建设中取得良好效果的三个城市:大庆,大连和焦作等案例进行细致的分析和研究,归纳出城市形象评价的各个要素、构建出城
互联网技术革新及企业业务规模不断扩大,各企业在信息化建设方面的投入也随之增加。信息化优势凸显的同时也使企业生产系统数目不断攀升及功能愈演庞大复杂,然而传统的维护方