文档库 最新最全的文档下载
当前位置:文档库 › 偏微分方程理论学习

偏微分方程理论学习

偏微分方程理论学习
偏微分方程理论学习

偏微分方程理论学习

一. 偏微分方程发展简介

1. 常微分方程

十七世纪微积分创立之后,常微分方程理论立刻就发展起来,当时应用常微分方程,解决几何与理学中的新问题。结果是在天体理学中不仅能得到并解释早先已经知晓的那些事实,而且得到了性的发现(例如,海王星的发现就是在对微分方程分析的基础上作出的)。

2. 偏微分方程

偏微分方程的研究要晚得多,对物理学中出现的偏微分方程研究在十八世纪中叶导致了分析学的一个新的分支------数学物理方程的建立。

J.达朗贝尔(D’Alembert )(1717-1783)、L.欧拉(Euler )(1707-1783)、D.伯努利(Bernoulli )(1700-1782)、J.拉格朗日(Lagrange )(1736-1813)、P.拉普拉斯(Laplace )(1749-1827)、S.泊松(Poisson )(1781-1840)、J.傅里叶(Fourier )(1768-1830)等人的工作为这一学科分支奠定了基础。它们在考察具体的数学物理问题中,所提出的思想与方法,竟适用于众多类型的微分方程,成为十九世纪末偏微分方程一般理论发展的基础。

十九世纪,偏微分方程发展的序幕是由法国数学家傅里叶拉开的,他于1822年发表的《热的解析理论》是数学史上的经典文献之一。傅里叶研究的主要是吸热或放热物体内部任何点处的温度随空间和时间的变化规律。在对物体的物理性状作出一定的限制(如均匀、各向同性)后,他根据物理原理推导出了三维空间的热传导方程

其中k 是一个参数,其值依赖于物体的质料。傅里叶当时解决的是如下特殊的热传导问题:设所考虑的物体为两端保持在温度0度、表面绝热且无热流通过的柱轴。在此情形下求解上述热传导方程,因为柱轴只涉及一维空间,所以这个问题也就是求解偏微分方程

???

????<<=>==??=??,0),()0,(,0,0),(,0),0(T T 222l x x f x T t t l T t T x k x ,

其中后面两项分别是边界条件和初始条件。傅里叶为解这个方程用了分离变量法,他得到满足方程和边界条件的级数解为

为了满足初始条件,必须有

这就促使傅里叶不得不考虑任给一个函数,能否将它表示成三角级数的问题。傅里叶得出的结论是:每个函数都可以表示成

这样,每个n b 可由上式乘以,...)2,1(sin =n nx ,再从0到π积分而得到。他还指出这个程序可以应用于表达式

接着,他考虑了任何函数)(x f 在区间),(ππ-的表达式,利用对称区间上的任何函数可以表示成一个奇函数和一个偶函数之和这一事实,傅里叶可以将区间),(ππ-上的任何函数)(x f 表示为

其系数由

确定,这就是我们通常所称的傅里叶级数。

为了处理无穷区域上的热传导问题,傅里叶同时还导出了现在所谓的“傅里叶积分”:

需要指出的是,傅里叶从没有对“任意”函数可以展成傅里叶级数这一断言给出过任何完全的证明,它也没有说出一个函数可以展开为三角级数必须满足的条件。然而傅里叶本人对此充满信心,因为他的信念有几何上的根据。傅里叶的工作不仅发展了偏微分方程的理论,而且使函数概念得以改进,同时也标志着人们从解析函数或可展成泰勒级数的函数中解放出来。傅里叶的前辈都曾坚持一个函数必须是可用单个式子表示的,而傅里叶级数却可以表示那些在区间),0(π或),(ππ-的不同部分有不同解析式的函数,不论这些表示式相互是否连续地接合着。特别是,一个傅里叶级数是在一整段区间上表示一个函数的,而一个泰勒级数仅在函数的解析点附近表示该函数。

事实上,傅里叶的主要思想早在1807年他提交巴黎科学院的一篇关于热传导的论文中就出现了,但是这篇论文在拉格朗日等人评审后遭到拒绝。1811年,他又提交了经过修改的论文,以争取科学院为热传导问题所设立的高额奖金。这次他虽然获了奖,但仍因受到缺乏严格性的批评而未能将论文发表在当时科学院的《报告》里。1824年,傅里叶成为科学院的秘书,这回他终于能够把他1811年的论文原封不动地发表在《报告》里,而这已经是在他的名著《热的解析理论》出版两年以后的事情了。

十九世纪偏微分方程的另一个重要发展是围绕着位势方程来进行的,这方面的代表人物格林(G .. Green)是一位磨坊工出身、自学成才的英国数学家。位势方程也称拉普拉斯方程:

拉普拉斯曾采用球面调和函数法解这个方程,不过他得到一个错误的结论,认为这个方程当被吸引的点(x,y,z)位于物体内部时也成立。这个错误由泊松加以更正。泊松指出,如果点(x,y,z)在吸引体内

部,则满足方程πρ4V -=?,其中ρ是吸引体密度,它也是x,y,z 的一个函数。拉普拉斯和泊松的方法都只适用于特殊的几何体,格林则认识到函数V 的重要性,并赋予它“位势”(potential)的名称,与前人不同的是,格林发展了函数V 的一般理论。他求解位势方程的方法与用特殊函数的级数方法相反,称为奇异点方法。他在1828年私人印刷出版的小册子《关于数学分析应用于电磁学理论的一篇论文》中,建立了许多推动位势论的进一步发展极为关键的定理与概念,其中以格林公式

???????-??=?-?σd n

U V n V U dv U V V U )()( (n 为物体表面指向外部的法向,dv 是体积元,d σ是面积元)和作为一种带奇异性的特殊位势的格林函数概念影响最为深远。

格林是剑桥数学物理学派的开山祖师,他的工作培育了汤姆逊(W.Thomson)、斯托克斯(G.Stokes)、麦克斯韦(J.C.Maxwell)等强有力的后继者,他们是十九世纪典型的数学物理学家。他们的主要目标,是发展求解重要物理问题的一般数学方法,而他们手中的主要武器就是偏微分方程,以至于在十九世纪,偏微分方程几乎变成了数学物理的同义词。

剑桥数学物理学派的贡献使经历了一个多世纪沉寂后英国数学在十九世纪得以复兴,麦克斯韦1864年导出的电磁场方程 ,)(1rot t

E c H ??=ε ,)(1rot t

H c E ??-=μ ,)(ρε=E div

0)(=H div μ

是十九世纪数学物理最壮观的胜利,正是根据对这组方程的研究,麦克斯韦预言了电磁波的存在,不仅给科学和技术带来巨大的冲击,同时也是偏微分方程威名大振。爱因斯坦在一次纪念麦克斯韦的演讲中说:“偏微分方程进入理论物理学时是婢女,但逐渐变成了主妇,”他认为这是从十九世纪开始的,而剑桥数学物理学派尤其是麦克斯韦在这一转变中起了重要的作用。

除了麦克斯韦方程,十九世纪导出的著名偏微分方程组还有粘性流体运动的纳维(C.L.M.H. Navier)-斯托克斯和弹性介质的柯西方程等。所有这些方程都不存在普遍解法。不过,十九世纪的数学家们已经逐渐认识到在偏微分方程的情形,无论是单个方程还是方程组,通解实际上不如初始条件和边界条件已给出的特殊问题的解有用。因此他们在求解定结问题方面作了大量工作。

对18、19世纪建立起来类型众多的微分方程,数学家们求显式解的努力往往归于失败,这种情况促使他们转而证明解的存在性。最先考虑微分方程解的存在性问题的数学家是柯西。他指出:在求显式解无效的场合常常可以证明解的存在性。他在19世纪20年代对形如y)f(x,y'=的常微分方程给出了第一个存在性定理,这方面的工作被德国数学家李普希茨(R. Lipschitz)、法国数学家刘维尔(J.Liouville)和皮卡(C.E. Picard)等追随。柯西也是讨论偏微分方程解的存在性的第一人,他在1848年的一系列论文中论述了如何将任意阶数大于1的偏微分方程化为偏微分方程组,然后讨论了偏微分方程组解的存在性并提出了证明存在性的强函数方法。柯西的工作后来被俄国女数学家柯瓦列夫斯卡娅(C.B. Ковалевская)独立地发展为包括拟线性方程和高阶组在内非常一般的形式。有关偏微分方程解的存在唯一性定理在现代文献中就称为“柯西-柯瓦列夫斯卡娅定理”。

柯瓦列夫斯卡娅是历史上为数不多的杰出女数学家之一。她出生于莫斯科一个贵族家庭,17岁时就在彼得堡一位海军学校教师指导下掌握了微积分。然而当时俄国的大学拒收女生,为了求学深造,他只好出走德国,先在海德堡大学学习一年,后来慕名到柏林求见威尔斯特拉斯。初次见面,威尔斯特拉斯出了一堆难题考她,估计她多半做不出来,但一周以后,当柯瓦列夫斯卡娅如期带着

完满的答卷回来见他时,这位名重一时的数学家对她的数学才能不再怀疑。当时的柏林大学跟俄国的大学一样不收女生,威尔斯特拉斯决定为柯瓦列夫斯卡娅单独授课,每星期日下午一次,四年不曾中断。在这四年时间里,柯瓦列夫斯卡娅不仅学完了大学的全部数学课程,而且还写出了三篇重要论文,其中一篇就是前面提到的关于偏微分方程解存在性的研究。这些工作是那么出色,以至于哥廷根大学在没有经过考试和答辩的情况下破格授予她博士学位,使她成为历史上第一位女数学博士。

由于18世纪的大量开发,常微分方程的求解在19世纪反而局限于用分离变量法解偏微分方程时所得到的那些方程,并且多半使用级数解,这引导出一串特殊函数,如贝塞尔(Bessel)函数、高斯(Gauss)超几何函数等等。在十九世纪后半叶,对常微分方程研究的理论方面变得突出,并且在常微分方程解析理论和定性理论两个大的方向上开拓了常微分研究的新局面,其中重大发展都与庞加莱(H. Poincare)的名字联系着。

庞加莱从27岁起任巴黎大学教授,直到他去世。他是欧拉、柯西之后最多产的数学家,并且在研究领域的广泛方面很少有人能与他相比。每年他在巴黎大学讲授一门不同的科目,而在每一门科目中,他都留着他自己的创造印记。

庞加莱、克莱因和希尔伯特,是在19和20世纪数学交界线上高耸着的三个巨大身影。他们放射着19世纪数学的光辉,同时照耀着通往20世纪数学的道路。在19世纪末,数学发展呈现出一派生机勃勃的景象,这与18世纪形成了鲜明的对比。无论从内部需要还是外部应用看,数学家们似乎都有做不完的问题。1900年8月5日,庞加莱宣布巴黎国际数学家大会开幕,正是在这次会议期间,希尔伯特充满信心地走上讲台,以他著名的23个问题揭开了20世纪数学的序幕。

当研究在解决物理问题的过程中出现的具体微分方程时,往往会产生一些极具普遍性、起初并没有严格的数学根据而应用于范围广泛物理问题的方法。例如,傅里叶方法、里茨(Ritz)方法、伽辽金(Галёркин)方法、摄动理论方法等就是这一类方法。这些方法应用的有效性成为试图对它们进行严格论证的原因之一。这就导致新的数学理论、新的研究方向的建立(傅里叶积分理论、本证函数展开理论和广义函数论等等)。

二、偏微分方程理论的两个特点

1. 偏微分方程理论与应用、与物理问题的直接联系

偏微分方程理论产生于那些归结为考察某些具体偏微分方程的具体物理问题的研究,这些方程便得到数学物理方程的称谓。

数学在物理中应用的历史较长,18世纪是数学与经典力学相结合的黄金时期,19世纪数学应用的重点转移到电学与电磁学,并且由于剑桥学派的努力而形成了数学物理分支。进入20世纪以后,随着物理科学的发展,数学相继在应用于相对论、量子力学以及基本粒子理论等方面取得了一个又一个突破,极大地丰富了数学物理的内容,同时也反过来刺激了数学自身的进步。

在20世纪初狭义相对论和广义相对论的创立过程中,数学都建有奇功。1907年,德国数学家闵可夫斯基(Minkowski)提出了“闵可夫斯基空间”,即将时间与空间融合在一起的四维时空13

R,。闵可夫斯基几何为爱因斯坦狭义相对论提供了适合的数学模型。有了闵可夫斯基时空模型后,爱因斯坦又进一步研究引力场理论以建立广义相对论。1912年夏他已经概括出新的引力理论的基本物理原理,但为了实现广义相对论的目标,还必须寻求理论的数学结构,爱因斯坦为此花费了3年的时间,最后在数学家格罗斯曼(Grossmann)介绍下掌握了发展相对论引力学说所必需的数学工具-----以黎曼几何为基础的绝对微分学,亦即爱因斯坦后来所称的张量分析。在1915年11月25日发表的一篇论文中,爱因斯坦终于导出了广义协变的引力场方程

g就是黎曼度量张量。爱因斯坦指出:“由于这组方程,广义相对论作为一种逻辑结构终于大功告

成!”

根据爱因斯坦的理论,时空整体是不均匀的,只是在微小的区域内可以近似地看作均匀。在数学上,广义相对论的时空可以理解为一种黎曼空间,非均匀时空连续区可借助于现成的黎曼度量来描述。这样,广义相对的数学表述第一次揭示了非欧几何的现实意义,成为历史上数学应用最伟大的例子之一。

20世纪数学物理的另一项经典成果是量子力学数学基础的确立。20世纪初,普朗克(M. Planck)、爱因斯坦、玻尔(N. Bohr)等创立了量子力学,但是到1925年为止,还没有一种量子理论能以统一的结构来概括这一领域已经积累的知识,当时的量子力学可以说是本质上相互独立的、有时甚至相互矛盾部分的混合体。1925年有了重要进展,由海森堡(W. Heisenberg)建立的矩阵力学和由薛定谔发展的波动力学形成了两大量子理论,而进一步将这两大理论融合为统一的体系,便成为当时科学界的当务之急。恰恰在这时,数学又起了意想不到的但却是决定性的作用。1927年,希尔伯特和冯·诺依曼、诺德海姆(L. Nordheim)合作发表了论文《论量子力学基础》,开始了用积分方程等分析工具使量子力学统一化的努力。在随后两年中,冯·诺依曼又进一步利用他从希尔伯特关于积分方程的工作中提炼出来的抽象希尔伯特空间理论,去解决量子力学的特征值问题,并最终将希尔伯特的谱理论推广到量子力学中经常出现的无界算子情形,从而奠定了量子力学的严格数学基础。1932年,冯·诺依曼发表了总结性著作《量子力学的数学基础》,完成了量子力学的公理化。

抽象的数学成果最终成为其他科学新理论的仿佛是量身定做的工具,在20世纪下半叶又演出了精彩的一幕,这就是大范围微分几何在统一场论中的应用。广义相对论的发展,逐渐促使科学家们去寻求电磁场与引力场的统一表述,这方面第一个大胆的尝试是数学家外尔(H. Weyl)在1918年提出的规范场理论,外尔自己称之为“规范不变几何”。统一场论的探索后来又扩展到基本粒子间的强相互作用和弱相互作用。1954年,物理学家杨振宁和米尔斯(R. L. Mills)提出的“杨-米尔斯理论”,揭示了规范不变性可能是所有四种(电磁、引力、强、弱)相互作用的共性,开辟了用规范场论来统一自然界这4中相互作用的新途径。数学家们很快就注意到杨-米尔斯理论所需要的数学工具早已存在,物理规范势实际上就是微分几何中纤维丛上的联络,20世纪30、40年代以来已经得到深入的研究。不仅如此,人们还发现规范场的杨-米尔斯方程是一组在数学上有重要意义的非线性偏微分方程。1975年以来,对杨-米尔斯方程的研究取得了许多重要成果,展示了统一场论的诱人前景,同时也推动了数学自身的发展。

数学不仅在物理、化学等传统学科中有着广泛而重要的应用,数学在生物学中应用自20世纪初以来得到了很大发展。1926年,意大利数学家伏尔泰拉(V. V olterra)提出著名的伏尔泰拉方程。从此微分方程又成为建立各种生物模型的重要工具。用微分方程建立生物模型在20世纪50年代曾获得轰动性成果,这就是描述神经脉冲传导过程的数学模型霍奇金(Hodgkin)-哈斯利(Huxley)方程(1972)和描述视觉系统侧抑制作用的哈特莱茵(Hartline)-拉特利夫(Ratliff)方程(1958),它们都是复杂的非线性方程组,引起了数学家和生物学家的浓厚兴趣。这两项工作分别获得1963年和1967年度诺贝尔医学生理学奖。

与物理、化学和生物、甚至于经济领域现象有关的数学问题提出,导致现象的数学理想化,或者换句话说,导致建立描述所研究的各类现象基本规律的数学模型。对于一系列现象的数学模型的建立在于归结为以基本物理、经济规律为基础的方程,这些模型仅仅考虑到现象的本质特点而忽略一系列次要的特点。例如,动量守恒、能量守恒、质量守恒等就是这样的规律。用这种方法可以得到在电动力学、声学、弹性力学、流体动力学以及其他连续介质力学的分支所研究的物理现象的方程。用数学方法研究数学模型不仅可以得到物理现象的定量特征,以给定的精度计算实际过程,还可能洞察物理现象的本质,有时还可以预言新的效果。

3. 偏微分方程理论与其他数学分支如泛函分析、函数论、拓扑学、代数、复分析的紧密联系。

偏微分方程理论广泛应用数学这些领域中的基本概念、基础思想和基本方法,并且它本身也给这些学科分支的研究问题的范围与方向以影响。弦振动的研究就是这种相互影响的经典范例。弦振动是达朗贝尔于1747年建立的,它还得到了表达这个方程通解的公式。欧拉得出弦振动方程柯西(Cauchy )问题解的公式:这个公式今天称为达朗贝尔公式。D.伯努利断言:弦振动方程的任何解均可表示为三角级数。欧拉同达朗贝尔、D.伯努利关于弦振动方程解的性质的争论,对数学物理、分析学,特别是三角级数理论的发展具有重要意义。

J .傅里叶在1822年进一步研究了用三角级数表示函数的问题,这是与热传到问题有关的。随后在L.狄利克雷(Dirichlet )(1805-1859)的工作中最先指出了把函数展开成三角级数的充分条件。最先出现在数学物理问题中的把函数表示成三角级数的问题在很大程度上促成了现代的集论与函数论的建立。

偏微分方程在几何上的应用产生了微分几何,古典微分几何多是局部性即小范围的。黎曼几何在空间每一点附近建立局部的二次微分型式

∑=.d 2j i ij dx dx g s

19世纪末,意大利数学家里奇(C.G .Ricci )发展了黎曼关于微分型式不变量的研究,开创了所谓“绝对微分学”,即现在的张量分析,系统地研究里面度量在坐标变换之下的不变性。1917年,里奇的学生列维(T. Levi)-奇维塔(Civita)引进“列维-奇维塔平移”,将欧氏空间的平行概念推广到弯曲空间,是黎曼几何具有了明显的几何意义。后来外尔发现平行性与空间的度量性质无关,从而建立所谓仿射联络(1918),摆脱度量定义平移与曲率,从而建立更广泛的几何理论。1920年以后嘉当(E. Cartan)发展了一般的联络理论与活动标架法。嘉当联络是纤维丛概念的先声,但在20世纪30年代以前,黎曼几何的研究,包括嘉当的工作,主要是小范围的。1925年,霍普夫(Hopf)注意到黎曼空间的微分几何结构与拓扑结构的关系,微分几何开始经历从局部到整体的转移。

整体微分几何以研究微分几何(小范围)性质与大范围性质之间的联系为目标。由于纤维丛的概念反映了流形的固有的图片性质,它提供了从局部研究想整体研究过渡的合适机制。因此整体微分几何的研究与微分拓扑学便有不解之缘,纤维丛与示性类的引入,使整体微分几何的研究出现了突破,陈省身在这方面有奠基性的贡献。

微分几何本来就是分析在几何中的应用,整体微分几何则表现出与现代分析更深刻的联系,特别是非线性偏微分方程理论的运用,引出了整体微分几何的重大成果。典型的例子是丘成桐1976年解决了微分几何领域里著名的“卡拉比猜想”。这是给定里奇曲率求黎曼度量的问题,其中需要求解高难度的非线性偏微分方程。丘成桐还解决了一系列与非线性偏微分方程有关的其他几何问题,并证明了广义相对论中的正质量猜想等等。由于这些工作,1982年丘成桐荣获菲尔兹奖。

三、偏微分方程理论的内容

偏微分方程是数学的中心,不论是纯粹数学还是应用数学。它们通常发生于因变量作为以空间和时间为自变量的连续变化函数的数学模型中。它们最引人注目的特性是其普适性,这一特性使我们能够从流体和固体力学、电磁学、概率、金融到众多应用领域中找到偏微分方程理论中每个数学概念的来源。而且这种可应用性随着现代软件适用于这种方程的离散逼近格式的灵活性和威力的增加而日益增长。同样戏剧性的是所有这些应用领域中方程的产生方式能容易地成为非常重要和深刻的基本数学问题的研究动机,并且反过来从这些研究的成果中获益。

不管它是否作为一种物理现象的模型,偏微分方程的分析有许多目的。其中的一个主要目的是适定性。粗略地说,一个偏微分方程是适定的,若它有解(存在性)、解唯一(唯一性)、且对输入数据的微小改变的响应也是很小的改变(连续依赖性)。前两个准则是一个有意义的物理模型所要求的,第三个准则是实验观察的基础。考虑适定性时,还应记得对有实际意义的问题通畅不可能求得显示解,从而逼近格式,特别是数值解在应用中就具有特别的重要性。因此,适定性问题与偏微分方程科学计算的如下中心问题有密切联系:对一个问题给定一定精度的数据,数值解计算输出有多少精度?正因为这个问题对现代定量科学的重要性,适定性成偏微分方程理论的核心内容。

本课程作为偏微分方程的入门课程,主要研究椭圆、抛物和双曲这三类线性偏微分方程的适定性问题,它们分别以拉普拉斯方程、热传导方程和波动方程为代表。具体地,对于某些规则的求解区域试图求出满足特定线性偏微分方程和定解条件的具体解,这就解决了存在性问题;再利用方程本身所具有的特殊性质,将证明所求的解是唯一的,也就是解决了唯一问题;关于连续依赖性问题,需要在不同的函数空间中考虑,我们将在连续函数空间和平方可积函数空间中分别讨论解对于输入数据的连续依赖问题。

四、偏微分方程理论的学习

有鉴于上述内容,偏微分方程的学习以三类线性偏微分方程的适定性问题为主线。同时,考虑到偏微分方程理论的两个特点:一是与应用、与物理的紧密联系;二是与数学其他分支的联系。

对于每一类方程,首先需要对方程的导出有一定的了解。事实上,同一个方程可能有许多不同的来源,这一方面是偏微分方程理论具有广泛应用的原因之一。同时对于不同的来源进行类比研究可以更好地揭示物理过程的某些特性,因为某个具体物理特性在某个物理过程还没有被观察到或没有引起注意,而在另外某个物理过程已经被观察注意到了,如果这两个物理过程服从同一个偏微分方程,则在原来的物理过程中应该也具有这个特性。其次,在对数学模型进行研究之后,需要有意识地将数学解代回原来的物理意义中,去理解、解释物理现象。这一方面可以验证数学模型的有效性,另一方面可以更好地理解已知的物理现象,甚至于预言新的物理现象。

只有这样才能显示出偏微分方程理论是从物理等具体学科中发展出来,同时又服务于这些具体学科,逐步提高分析、解决实际问题的能力。

至于与数学其他学科的联系,首先,求解过程中将会用到许多微积分或数学分析的概念、思想和定理,解的表达形式也是有积分形式的或级数形式的,解空间的结构则用到许多线性代数的知识;其次,唯一性和连续依赖性则需要许多实变和泛函分析的理论和方法。事实上,实变函数和泛函分析的许多思想方法都是来源于偏微分方程理论研究的。

数学理想化

数学求解过程

物理解释

参考文献

1.李文林.数学史概论.北京:高等教育出版社,200

2.8

2.陈祖墀.偏微分方程.合肥:中国科学技术大学出版社,2004.8

3.N.Asmar.Partial Differential Equations with Fourier Series and

Boundary Value Problems.陈祖墀,宣本金译.北京:机械工业出版社,2006.10

偏微分方程数值解期末试题及标准答案

偏微分方程数值解试题(06B ) 参考答案与评分标准 信息与计算科学专业 一(10分)、设矩阵A 对称,定义)(),(),(2 1)(n R x x b x Ax x J ∈-=,)()(0x x J λλ?+=.若0)0('=?,则称称0x 是)(x J 的驻点(或稳定点).矩阵A 对称(不必正定),求证0x 是)(x J 的驻点的充要条件是:0x 是方程组 b Ax =的解 解: 设n R x ∈0是)(x J 的驻点,对于任意的n R x ∈,令 ),(2),()()()(2 000x Ax x b Ax x J x x J λλλλ?+-+=+=, (3分) 0)0('=?,即对于任意的n R x ∈,0),(0=-x b Ax ,特别取b Ax x -=0,则有0||||),(2000=-=--b Ax b Ax b Ax ,得到b Ax =0. (3分) 反之,若n R x ∈0满足b Ax =0,则对于任意的x ,)(),(2 1)0()1()(00x J x Ax x x J >+==+??,因此0x 是)(x J 的最小值点. (4分) 评分标准:)(λ?的展开式3分, 每问3分,推理逻辑性1分 二(10分)、 对于两点边值问题:?????==∈=+-=0 )(,0)(),()('b u a u b a x f qu dx du p dx d Lu 其中]),([,0]),,([,0)(min )(]),,([0min ],[1b a H f q b a C q p x p x p b a C p b a x ∈≥∈>=≥∈∈ 建立与上述两点边值问题等价的变分问题的两种形式:求泛函极小的Ritz 形式和Galerkin 形式的变分方程。 解: 设}0)(),,(|{11=∈=a u b a H u u H E 为求解函数空间,检验函数空间.取),(1b a H v E ∈,乘方程两端,积分应用分部积分得到 (3分) )().(),(v f fvdx dx quv dx dv dx du p v u a b a b a ==+=??,),(1 b a H v E ∈? 即变分问题的Galerkin 形式. (3分)

偏微分方程数值解实验报告

偏微分方程数值解实验报告

1、用有限元方法求下列边值问题的数值解:''()112x -y +y =2s i n ,0∈∈??∈(0,)?, 其中取1ν= 要求画出解曲面。迭代格式如下: 1221212111111111122142212n n n n n n j j j j j j n n n n n n j j j j j j V V V V V V h h V V V V V V h h τ++++++++++-+-??-()-()()-()??++?????? ??-+-+??=+??????

1、 %Ritz Galerkin方法求解方程 function u1=Ritz(x) %定义步长 h=1/100; x=0:h:1; n=1/h; a=zeros(n-1,1); b=zeros(n,1); c=zeros(n-1,1); d=zeros(n,1); %求解Ritz方法中内点系数矩阵 for i=1:1:n-1 b(i)=(1/h+h*pi*pi/12)*2; d(i)=h*pi*pi/2*sin(pi/2*(x(i)+h))/2+h*pi*pi/2*sin(pi/2*x(i+1))/2; end %右侧导数条件边界点的计算 b(n)=(1/h+h*pi*pi/12); d(n)=h*pi*pi/2*sin(pi/2*(x(i)+h))/2; for i=1:1:n-1 a(i)=-1/h+h*pi*pi/24; c(i)=-1/h+h*pi*pi/24; end %调用追赶法 u=yy(a,b,c,d) %得到数值解向量 u1=[0,u] %对分段区间做图 plot(x,u1) %得到解析解 y1=sin(pi/2*x); hold on plot(x,y1,'o') legend('数值解','解析解') function x=yy(a,b,c,d) n=length(b); q=zeros(n,1); p=zeros(n,1); q(1)=b(1); p(1)=d(1); for i=2:1:n

(完整版)偏微分方程的MATLAB解法

引言 偏微分方程定解问题有着广泛的应用背景。人们用偏微分方程来描述、解释或者预见各种自然现象,并用于科学和工程技术的各个领域fll。然而,对于广大应用工作者来说,从偏微分方程模型出发,使用有限元法或有限差分法求解都要耗费很大的工作量,才能得到数值解。现在,MATLAB PDEToolbox已实现对于空间二维问题高速、准确的求解过程。 偏微分方程 如果一个微分方程中出现的未知函数只含一个自变量,这个方程叫做常微分方程,也简称微分方程;如果一个微分方程中出现多元函数的偏导数,或者说如果未知函数和几个变量有关,而且方程中出现未知函数对几个变量的导数,那么这种微分方程就是偏微分方程。 常用的方法有变分法和有限差分法。变分法是把定解问题转化成变分问题,再求变分问题的近似解;有限差分法是把定解问题转化成代数方程,然后用计算机进行计算;还有一种更有意义的模拟法,它用另一个物理的问题实验研究来代替所研究某个物理问题的定解。虽然物理现象本质不同,但是抽象地表示在数学上是同一个定解问题,如研究某个不规则形状的物体里的稳定温度分布问题,由于求解比较困难,可作相应的静电场或稳恒电流场实验研究,测定场中各处的电势,从而也解决了所研究的稳定温度场中的温度分布问题。 随着物理科学所研究的现象在广度和深度两方面的扩展,偏微分方程的应用范围更广泛。从数学自身的角度看,偏微分方程的求解促使数学在函数论、变分法、级数展开、常微分方程、代数、微分几何等各方面进行发展。从这个角度说,偏微分方程变成了数学的中心。

一、MATLAB方法简介及应用 1.1 MATLAB简介 MATLAB是美国MathWorks公司出品的商业数学软件,用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境,主要包括MATLAB和Simulink两大部分。 1.2 Matlab主要功能 数值分析 数值和符号计算 工程与科学绘图 控制系统的设计与仿真 数字图像处理 数字信号处理 通讯系统设计与仿真 财务与金融工程 1.3 优势特点 1) 高效的数值计算及符号计算功能,能使用户从繁杂的数学运算分析中解脱出来; 2) 具有完备的图形处理功能,实现计算结果和编程的可视化; 3) 友好的用户界面及接近数学表达式的自然化语言,使学者易于学习和掌握; 4) 功能丰富的应用工具箱(如信号处理工具箱、通信工具箱等) ,

差分法求解偏微分方程MAAB

南京理工大学 课程考核论文 课程名称:高等数值分析 论文题目:有限差分法求解偏微分方程 姓名:罗晨 学号: 成绩: 有限差分法求解偏微分方程 一、主要内容 1.有限差分法求解偏微分方程,偏微分方程如一般形式的一维抛物线型方程:具体求解的偏微分方程如下: 2.推导五种差分格式、截断误差并分析其稳定性; 3.编写MATLAB程序实现五种差分格式对偏微分方程的求解及误差分析;

4.结论及完成本次实验报告的感想。 二、推导几种差分格式的过程: 有限差分法(finite-differencemethods )是一种数值方法通过有限个微分方程近似求导从而寻求微分方程的近似解。有限差分法的基本思想是把连续的定解区域用有限个离散点构成的网格来代替;把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似;把原方程和定解条件中的微商用差商来近似,积分用积分和来近似,于是原微分方程和定解条件就近似地代之以代数方程组,即有限差分方程组,解此方程组就可以得到原问题在离散点上的近似解。 推导差分方程的过程中需要用到的泰勒展开公式如下: ()2100000000()()()()()()()......()(()) 1!2!! n n n f x f x f x f x f x x x x x x x o x x n +'''=+-+-++-+-(2-1) 求解区域的网格划分步长参数如下: 11k k k k t t x x h τ ++-=?? -=?(2-2) 2.1古典显格式 2.1.1古典显格式的推导 由泰勒展开公式将(,)u x t 对时间展开得 2,(,)(,)( )()(())i i k i k k k u u x t u x t t t o t t t ?=+-+-?(2-3) 当1k t t +=时有 21,112,(,)(,)( )()(())(,)()() i k i k i k k k k k i k i k u u x t u x t t t o t t t u u x t o t ττ+++?=+-+-??=+?+?(2-4) 得到对时间的一阶偏导数 1,(,)(,)()=()i k i k i k u x t u x t u o t ττ+-?+?(2-5) 由泰勒展开公式将(,)u x t 对位置展开得 223,,21(,)(,)()()()()(())2!k i k i k i i k i i u u u x t u x t x x x x o x x x x ??=+-+-+-??(2-6) 当11i i x x x x +-==和时,代入式(2-6)得

Maab求解微分方程组及偏微分方程组

第四讲 Matlab 求解微分方程(组) 理论介绍:Matlab 求解微分方程(组)命令 求解实例:Matlab 求解微分方程(组)实例 实际应用问题通过数学建模所归纳得到的方程,绝大多数都是微分方程,真正能得到代数方程的机会很少.另一方面,能够求解的微分方程也是十分有限的,特别是高阶方程和偏微分方程(组).这就要求我们必须研究微分方程(组)的解法:解析解法和数值解法. 一.相关函数、命令及简介 1.在Matlab 中,用大写字母D 表示导数,Dy 表示y 关于自变量的一阶导数,D2y 表示y 关于自变量的二阶导数,依此类推.函数dsolve 用来解决常微分方程(组)的求解问题,调用格式为: X=dsolve(‘eqn1’,’eqn2’,…) 函数dsolve 用来解符号常微分方程、方程组,如果没有初始条件,则求出通解,如果有初始条件,则求出特解. 注意,系统缺省的自变量为t 2.函数dsolve 求解的是常微分方程的精确解法,也称为常微分方程的符号解.但是,有大量的常微分方程虽然从理论上讲,其解是存在的,但我们却无法求出其解析解,此时,我们需要寻求方程的数值解,在求常微分方程数值解方面,MATLAB 具有丰富的函数,我们将其统称为solver ,其一般格式为: [T,Y]=solver(odefun,tspan,y0) 说明:(1)solver 为命令ode45、ode23、ode113、ode15s 、ode23s 、ode23t 、ode23tb 、ode15i 之一. (2)odefun 是显示微分方程'(,)y f t y =在积分区间tspan 0[,]f t t =上从0t 到f t 用初始条件0y 求解. (3)如果要获得微分方程问题在其他指定时间点012,,,,f t t t t L 上的解,则令tspan 012[,,,]f t t t t =L (要求是单调的). (4)因为没有一种算法可以有效的解决所有的ODE 问题,为此,Matlab 提供了多种求解器solver ,对于不同的ODE 问题,采用不同的solver.

偏微分方程数值解实验报告

精品文档 偏微分方程数值解 上 机 实 验 报 告 (一)实验一 一、上机题目: 用线性元求解下列边值问题的数值解:

精品文档 ′′22?? ?? ??,0

精品文档 (二)实验二 四、上机题目: 求解 Helmholtz 方程的边值问题: u k 2u 1 ,于(0,1)*(0,1) u0,于1{ x0,0y1} U{0x1, y 1} 1{ x0,0y1} U{0x1, y1} u 0,于2{0x1, y 0} U { x1,0y1} n 其中 k=1,5,10,15,20 五、实验程序:

偏微分方程数值及matlab实验报告.docx

偏微分方程数值实验报告八 实验题目:利用有限差分法求解 u ( x) u(x) f (x), u( 1) 0, u(1) 0. 真解为 u( x) e x 2 (1 x 2 ) 实现算法:对于两点边值问题 d 2u f , x l , dx 2 (1) u(a),u(b) , 其中 l ( a, b) (a b), f 为 l [ a,b] 上的连续函数, , 为给定常数 . 其相应的有限差分法的算法如下: 1.对求解区域做网格剖分,得到计算网格 .在这里我们对区间 l 均匀剖分 n 段,每个剖分单元 b a 的剖分步长记为 h . n 2.对微分方程中的各阶导数进行差分离散,得到差分方程 .运用的离散方法有: 方法一 :用待定系数和泰勒展开进行离散 d 2u( x i ) i 1 u( x i 1) i u( x i ) i 1 u( x i 1) d( x i )2 方法二:利用差商逼近导数 d 2u( x i ) u( x i 1 ) 2u( x i ) u( x i 1 ) d( x i )2 h 2 将(2) 带入 (1)可以得到 u(x i 1) 2u(x i ) u(x i 1 ) ) R i (u) , h 2 f ( x i 其中 R i (u) 为无穷小量,这时我们丢弃 R i (u) ,则有在 x i 处满足的计算公式: u(x i 1) 2u( x i ) u( x i 1 ) 1,..., n 1 h 2 f ( x i ), i 3.根据边界条件,进行边界处理 .由 (1)可得 u 0 , u n (2) (3) (4) 称(3)(4)为逼近 (1) 的差分方程,并称相应的数值解向量 U n 1 为差分解, u i 为 u( x i ) 的近似值 . 4.最后求解线性代数方程组,得到数值解向量U n 1 .

偏微分方程数值解例题答案

二、改进的Euler 方法 梯形方法的迭代公式(1.10)比Euler 方法精度高,但其计算较复杂,在应用公式(1.10)进行计算时,每迭代一次,都要重新计算函数),(y x f 的值,且还要判断何时可以终止或转下一步计算.为了控制计算量和简化计算法,通常只迭代一次就转入下一步计算.具体地说,我们先用Euler 公式求得一个初步的近似值1+n y ,称之为预测值,然后用公式(1.10)作一次迭代得1+n y ,即将1+n y 校正一次.这样建立的预测-校正方法称为改进的Euler 方法: 预测: ),,(1n n n n y x hf y y +=+ 校正 : )].,(),([2 111+++++=n n n n n n y x f y x f h y y (1.15) 这个计算公式也可以表示为 11(,), (,), 1(). 2p n n n c n n p n p c y y hf x y y y hf x y y y y ++?=+??=+?? ?=+??? 例1 取步长0.1h =,分别用Euler 方法及改进的Euler 方法求解初值问题 d (1),01, d (0) 1. y y xy x x y ?=-+≤≤???=? 解 这个初值问题的准确解为()1(21)x y x e x =--. 根据题设知 ).1(),(xy y y x f +-= (1) Euler 方法的计算式为 )],1([1.01n n n n n y x y y y +?-=+ 由1)0(0==y y , 得 ,9.0)]101(1[1.011=?+??-=y ,8019.0)]9.01.01(9.0[1.09.02=?+??-=y 这样继续计算下去,其结果列于表9.1. (2) 改进的Euler 方法的计算式为 110.1[(1)],0.1[(1)], 1(), 2p n n n n c n p n p n p c y y y x y y y y x y y y y ++?=-?+?=-?+??? ?=+??? 由1)0(0==y y ,得

偏微分方程式之求解

第六章偏微分方程式之求解 在化工的领域中,有不少程序之动态是由以偏微分方程式(Partial differential equation;PDE)所描述的,例如热与质量在空间中的传递等。这些用以描述实际问题的PDE,除非具有某些特定的方程式型态及条件,否则甚难以手算的方式找出解析解。而在数值求解方面,最常被采用的方法为有限差分法(finite difference)何有限元素法(finite element)。然对于某些不熟悉数值分析及程序编写的化工人而言,欲充分了解以偏微分方程式所描述之系统动态是相当不容易的,更遑论进一步的设计与分析了。 值得庆幸的是,MATLAB 的环境中提供了一个求解PDE 问题的工具箱,让使用者得以利用简单的指令或图形接口工具输入欲解的PDE,并求解。使得PDE 之数值解在弹指之间完成,使用者不在为数值法所苦恼,轻松掌握偏微分方程式系统的动态,并可进一步进行后续之设计工作。 本章将以循渐进的方式,介绍PDE 工具箱及其用法,并以数个典型的化工范例进行示范,期能使初学者很快熟悉PDE工具箱,并使用它来设计与分析以偏微方方程式所描述的程序系统。 6.1 偏微分方程式之分类 偏微分方程式可根据其阶数(order),线性或非线性型态,以及边界条件进行分类。 6.1.1依阶数的分类 偏微分方程式是以偏微分项中之最高次偏微分来定义其阶数,例如: 一阶偏微分方程式: xy 二阶偏微分方程式: 三阶偏微分方程式: 6.1.2 依非线性程度分类

偏微分方程式亦可以其线性或非线性情况,区分为线性 (linear),似线性 (quasilinear),以及非线性三类。 例如,以下之二阶偏微分方程式 (Constantinides and Mostoufi,1999) 可依系数 ( )之情况,进行如下表之归类 类别 情况 线性 似线性 系数 ( )为定值,或仅为 (x,y)函数 系数 ( )为依变数 (dependent variable)u 或其比方程式中之偏微 分低阶之偏微分项的函数,如 ( ) (x,y,u, u x, u y) 非线性 系数 ( )中,具有与原方程式之偏微分同阶数之变数,如 () (x,y,u, 2u x 2 , 2 u y 2, 2u x y) 另外,对于线性二阶偏微分方程式,可进一步将其分类为椭圆型 (elliptic) ,拋 物线型 (parabolic),以及双曲线型 (hyperbolic) 。具体上来说,此类偏微分方程 式二阶线性之 一般式为 系数a,b,c,d,e 和 f 是定值或为 u 的函数。若 g=0,则上式为其次是偏微分方程 式。式子 ( )之分类及代表性例子,请见下表 (c ~ u) a ~ u f 2 热传导或扩散方程式 u 2 u xt a() 2 u y 2 22 uu b( ) c( ) 2 x y x 2 d() 0 22 uu b c 2 x y y 2 d u e u fu g 0 xy 方程式类别 判断式 椭圆型 b 2 4ac 0 拋物线型 b 2 4ac 0 代表性范例 Laplace 方程式, Poisson 方程式, 22 uu 22 xy 22 uu 22 xy f (x,y) x 2 t (c ~ u) a ~ u f

偏微分方程求解方法及其比较

偏微分方程求解方法及其比较 发表时间:2008-12-11T09:32:01.530Z 来源:《科海故事博览科教创新》2008年第10期供稿作者:曹海洋吕淑娟王淑芬 [导读] 近些年来,无限维动力系统得到了很大的发展.随着对它研究的深入和计算能力的迅速提高,使得与之相关的数值研究越来越被人们关注.谱方法作为一种数值求解偏微分方程的方法,它具有无穷阶收敛性.因此,谱方法也就引起人们更多的关注. 摘要:近些年来,无限维动力系统得到了很大的发展.随着对它研究的深入和计算能力的迅速提高,使得与之相关的数值研究越来越被人们关注.谱方法作为一种数值求解偏微分方程的方法,它具有无穷阶收敛性.因此,谱方法也就引起人们更多的关注. 关键词:谱方法;偏微分;收敛;逼近; 1偏微分方程及其谱方法的介绍 偏微分方程主要借助于未知函数及其导数来刻画客观世界的物理量的一般变化规律。理论上,对偏微分方程解法的研究已经有很长的历史了。最初的研究工作主要集中在物理,力学,几何学等方面的具体问题,其经典代表是波动方程,热传导方程和位势方程(调和方程)。通过对这些问题的研究,形成了至今仍然使用的有效方法,例如,分离变量法,fourier变换法等。早期的偏微分方程研究主要集中在理论上,而在实际操作中其研究方法和研究结果都难以得到广泛的应用。求解的主要方法为:有限差分法,有限元法,谱方法。 谱方法起源于Ritz-Galerkin方法,它是以正交多项式(三角多项式,切比雪夫多项式,勒让得多项式等)作为基函数的Galerkin方法、Tau 方法或配置法,它们分别称为谱方法、Tau方法或拟谱方法(配点法),通称为谱方法。谱方法是以正交函数或固有函数为近似函数的计算方法。从函数近似角度看.谱方法可分为Fourier方法.Chebyshev或Legendre方法。前者适用于周期性问题,后两者适用于非周期性问题。而这些方法的基础就是建立空间基函数。 下面介绍几种正交多项式各种节点的取值方法及权重。 1) Chebyshev-Gauss: 2) Chebyshev-Gauss-Radau: x0 =1, 3) Chebyshev-Gauss-Lobatto: x0 =1, xN =1, 4)Legendre-Gauss: xj 是的零点且 5) Legendre-Gauss-Radau: xj 是的N+1个零点且 6) Legendre-Gauss-Lobatto: x0=-1,xN=1其它N-1个点是的零点且 下面介绍谱方法中最重要的Jacobi正交多项式其迭代公式为: 其中: Jacobi正交多项式满足正交性: 而Chebyshev多项式是令时Jacobi多项式的特殊形式,另外Legendre多项式是令时Jacobi多项式的特殊形式。 2 几种典型的谱方法 谱方法是以正交函数或固有函数为近似函数的计算方法。谱近似可以分为函数近似和方程近似两种近似方式。从函数近似角度看.谱方法可分为Fourier方法.Chebyshev或Legendre方法。前者适用于周期性问题,后两者适用于非周期性问题。从方程近似角度看,谱方法可分为在物理空间离散求解的Collocation法、在谱空间进行离散求解的Galerkin法,以及先在物理空间离散求积,再变换到谱空间求解的Pseudo-spectral法。Collocation法适用于非线性问题.Galerkin法适用于线性问题,而Pseudo-spectral法适用于展开方程时的非线性项的处理。谱方法的特点是对光滑函数指数性逼近的谱精度;以较少的网格点得到较高的精度;无相位误差;适合多尺度的波动性问题;计算精度高于其他方法。快速傅立叶变化的提出大大促进了谱方法的发展,迄今已有各种的谱方法计算格式被提出.并被应用于天文学、电磁学、地理学等各种问题的计算。 下面介绍一下应用于各个区域的几种谱方法: 1)以Fourier谱方法为例介绍谱方法解方程的主要过程 以一阶波动方程为例: 其中u(x,t)为方程的解,L是包含u和u关于空间变量的导数的算子,除了方程以有初始条件和适当的边界条件。 故可设其中为试探空间的基函数,ak(t)为展开系数,对于傅立叶谱方法中的共轭有: 其中从而利用其正交性和周期性可以减少工作量,另外再结合边界条件就可以求出来。 2) Galerkin方法是谱方法中十分经典的解偏微分方程的方法,但还有其局限性,而利用Hermite谱方法中依赖时间的权函数对经典的Galerkin方法进行拓展后的新的方法能适用范围扩大了很多。它能很好的应用在微分方程最优控制问题有限元方法的分析中,并且如果能够灵活运用利用Chebyshev方法、Galerkin方法和配置方法,则会形成更强的计算方法。如将Tau方法的思想成功地应用于奇数阶微分方程Petrov-Galerkin谱方法。 3)在无界区域上谱方法和拟谱方法发展了以Hermite函数和Laguerre函数为基函数的正交逼近和插值理论,在这些结果的基础上发展了全空间和半空间上数理方程的谱方法和拟谱方法,从而形成一种新的能更好解决误解区域问题的方法,此种方法被很好的应用于统计物理、量子力学和流体力学中。 4) 我们利用非一致带权Sobolev空间中的Jacobi多项式正交逼近和Jacobi-Gauss型插值理论,提出以Jacobi多项式为基函数的Jacobi谱方法和拟谱方法用来解决一些奇异问题和计算某些特定的无界区域问题。 5)有限谱方法是基于有限点、有限项的局域谱方法。这种方法要求近似函数应具有等同隔网格和非周期性的性质。有限谱方法分为基于非

偏微分方程数值解法答案

1. 课本2p 有证明 2. 课本812,p p 有说明 3. 课本1520,p p 有说明 4. Rit2法,设n u 是u 的n 维子空间,12,...n ???是n u 的一组基底,n u 中的任一元素n u 可 表为1n n i i i u c ?==∑ ,则,11 11()(,)(,)(,)(,)22j n n n n n n i j i j j i j j J u a u u f u a c c c f ???=== -=-∑∑是12,...n c c c 的二次函数,(,)(,)i j j i a a ????=,令 () 0n j J u c ?=?,从而得到12,...n c c c 满足1 (,)(,),1,2...n i j i j i a c f j n ???===∑,通过解线性方程组,求的i c ,代入1 n n i i i u c ?==∑, 从而得到近似解n u 的过程称为Rit2法 简而言之,Rit2法:为得到偏微分方程的有穷维解,构造了一个近似解,1 n n i i i u c ?== ∑, 利用,11 11()(,)(,)(,)(,)22j n n n n n n i j i j j i j j J u a u u f u a c c c f ???===-=-∑∑确定i c ,求得近似解n u 的过程 Galerkin 法:为求得1 n n i i i u c ? == ∑形式的近似解,在系数i c 使n u 关于n V u ∈,满足(,)(,) n a u V f V =,对任 意 n V u ∈或(取 ,1j V j n ?=≤≤) 1 (,)(,),1,2...n i j i j i a c f j n ???===∑的情况下确定i c ,从而得到近似解1 n n i i i u c ?==∑的过程称 Galerkin 法为 Rit2-Galerkin 法方程: 1 (,)(,)n i j i j i a c f ???==∑ 5. 有限元法:将偏微分方程转化为变分形式,选定单元的形状,对求解域作剖分,进而构 造基函数或单元形状函数,形成有限元空间,将偏微分方程转化成了有限元方程,利用 有效的有限元方程的解法,给出偏微分方程近似解的过程称为有限元法。 6. 解:对求解区间进行网格剖分,节点01......i n a x x x x b =<<<<=得到相邻节点1,i i x x -

偏微分方程数值及matlab实验报告

偏微分方程数值实验报告八 实验题目:利用有限差分法求解 . 0)1(,0)1(),()()(==-=+''-u u x f x u x u 真解为 ) 1()(22 x e x u x -=-实现算法:对于两点边值问题 , )(,)(,,d 22βα==∈=-b u a u l x f dx u (1) 其中),(b a l =f b a ),(<为],[b a l =上的连续函数,βα,为给定常数. 其相应的有限差分法的算法如下: 1.对求解区域做网格剖分,得到计算网格.在这里我们对区间l 均匀剖分n 段,每个剖分单元的剖分步长记为n a b h -= .2.对微分方程中的各阶导数进行差分离散,得到差分方程.运用的离散方法有:方法一:用待定系数和泰勒展开进行离散 )()()()(d ) (d 11112 2++--++≈i i i i i i i i x u x u x u x x u ααα方法二:利用差商逼近导数 2 112 2) ()(2)()(d )(d h x u x u x u x x u i i i i i -++-≈(2) 将(2)带入(1)可以得到 )()() ()(2)(2 11u R x f h x u x u x u i i i i i +=+-- -+, 其中)(u R i 为无穷小量,这时我们丢弃)(u R i ,则有在i x 处满足的计算公式: 1,...,1)() ()(2)(2 11-==+-- -+n i x f h x u x u x u i i i i ,(3) 3.根据边界条件,进行边界处理.由(1)可得 β α==n u u ,0(4) 称(3)(4)为逼近(1)的差分方程,并称相应的数值解向量1-n U 为差分解,i u 为)(i x u 的近似值.4.最后求解线性代数方程组,得到数值解向量1 -n U .

《偏微分方程概述及运用matlab求解偏微分方程常见问题》要点

北京航空航天大学 偏微分方程概述及运用matlab求解微分方 程求解常见问题 姓名徐敏 学号57000211 班级380911班 2011年6月

偏微分方程概述及运用matlab求解偏微分 方程常见问题 徐敏 摘要偏微分方程简介,matlab偏微分方程工具箱应用简介,用这个工具箱解方程的过程是:确定待解的偏微分方程;确定边界条件;确定方程所在域的几何形状;划分有限元;解方程 关键词MATLAB 偏微分方程程序 如果一个微分方程中出现的未知函数只含有一个自变量,这个方程叫做常微分方程,也简称微分方程:如果一个微分方程中出现多元函数的偏导数,或者说如果未知函数和几个变量有关,而且方程中出现未知函数对几个变量的导数,那么这种微分方程就是偏微分方程。 一,偏微分方程概述 偏微分方程是反映有关的未知变量关于时间的导数和关于空间变量的导数之间制约关系的等式。许多领域中的数学模型都可以用偏微分方程来描述,很多重要的物理、力学等学科的基本方程本身就是偏微分方程。早在微积分理论刚形成后不久,人们就开始用偏微分方程来描述、解释或预见各种自然现象,并将所得到的研究方法和研究成果运用于各门科学和工程技术中,不断地取得了显著的成效,显示了偏微分方程对于人类认识自然界基本规律的重要性。逐渐地,以物

理、力学等各门科学中的实际问题为背景的偏微分方程的研究成为传统应用数学中的一个最主要的内容,它直接联系着众多自然现象和实际问题,不断地提出和产生出需要解决的新课题和新方法,不断地促进着许多相关数学分支(如泛函分析、微分几何、计算数学等)的发展,并从它们之中引进许多有力的解决问题的工具。偏微分方程已经成为当代数学中的一个重要的组成部分,是纯粹数学的许多分支和自然科学及工程技术等领域之间的一座重要的桥梁。 在国外,对偏微分方程的应用发展是相当重视的。很多大学和研究单位都有应用偏微分方程的研究集体,并得到国家工业、科学部门及军方、航空航天等方面的大力资助。比如在国际上有重大影响的美国的Courant研究所、法国的信息与自动化国立研究所等都集中了相当多的偏微分方程的研究人员,并把数学模型、数学方法、应用软件及实际应用融为一体,在解决实际课题、推动学科发展及加速培养人才等方面都起了很大的作用。 在我国,偏微分方程的研究起步较晚。但解放后,在党和国家的大力号召和积极支持下,我国偏微分方程的研究工作发展比较迅速,涌现出一批在这一领域中做出杰出工作的数学家,如谷超豪院士、李大潜院士等,并在一些研究方向上达到了国际先进水平。但总体来说,偏微分方程的研究队伍的组织和水平、研究工作的广度和深度与世界先进水平相比还有很大的差距。因此,我们必须继续努力,大力加强应用偏微分方程的研究,逐步缩小与世界先进水平的差距 二,偏微分方程的内容

偏微分方程实验1

《偏微分方程数值解》 课程实验报告(一)

1 实验题目 给定初值问题(注:共两次上机时间4学时) 1 02 )0(42'≤≤?? ?=--=x y x y y 其精确解为 12)(2+-=-x e x y x 。取h=0.1, 分别用显式Euler 法、隐式梯形法、二级二阶Runge-Kutta 、三级三阶Runge-Kutta 、四级四阶Runge-Kutta 计算数值解,并与精确解比较。 2 求解方法 1、每种方法的迭代公式: (1)显式Euler 法 p[i]=p[i-1]+h*(*f)(x[i-1],p[i-1]); (2)隐式梯形法 p[i]=p[i-1]+h*(*f)(x[i-1],p[i-1]);; (3)二级二阶Runge-Kutta K1=(*f)(x[i-1],p[i-1]); K2=(*f)(x[i-1]+h/2,p[i-1]+h*K1/2); p[i]=p[i-1]+h*K2; (4)三级三阶Runge-Kutta K1=(*f)(x[i-1],p[i-1]); K2=(*f)(x[i-1]+h/2,p[i-1]+h*K1/2); K3=(*f)(x[i-1]+h/2,p[i-1]+h*K2/2); K4=(*f)(x[i-1]+h,p[i-1]+h*K3); p[i]=p[i-1]+h*(K1+2*K2+2*K3+K4)/6; (5)四级四阶Runge-Kutta K1=(*f)(x[i-1],p[i-1]);

K2=(*f)(x[i-1]+h/2,p[i-1]+h*K1/2); K3=(*f)(x[i-1]+h/2,p[i-1]+h*K2/2); K4=(*f)(x[i-1]+h,p[i-1]+h*K3); p[i]=p[i-1]+h*(K1+2*K2+2*K3+K4)/6; 3 程序源代码 #include #include #define N 11 double fun(double x,double y) //微分方程{ return -2*y-4*x; } double p_fun(double x) //原函数 { return exp(-2*x)-2*x+1; } double* Exact(double x0, double h) //精确解{ int i; double *p=new double[N]; for(i=0;i

偏微分方程数值解法

“十二五”国家重点图书出版规划项目 信息与计算科学丛书 67 偏微分方程数值解法 陈艳萍鲁祖亮刘利斌编著

内 容 简 介 本书试图用较少的篇幅描述偏微分方程的几种数值方法. 主要内容包括:Sobolev空间初步, 椭圆边值问题的变分问题, 椭圆问题的有限差分方法, 抛物型方程的有限差分方法, 双曲型方程的有限差分方法, 椭圆型方程的有限元方法, 抛物及双曲方程的有限元方法, 椭圆型方程的混合有限元方法, 谱方法等. 本书内容丰富, 深入浅出, 尽可能地用简单的方法来描述一些理论结果, 并根据作者对有限差分、有限元、混合有限元、谱方法的理解和研究生教学要求, 全面、客观地评价各种数值计算方法,并列举一些数值计算的例子, 阐述许多新的学术观点. 本书可作为高等学校数学系高年级本科生和研究生的教材或参考书, 也可作为计算数学工作者和从事科学与工程计算的科研人员的参考书. 图书在版编目(CIP)数据 偏微分方程数值解法/陈艳萍, 鲁祖亮, 刘利斌编著. —北京:科学出版社, 2015.1 (信息与计算科学丛书67) ISBN 978-7-03-000000-0 Ⅰ. ①偏… Ⅱ. ①陈… ②鲁… ③刘… Ⅲ. ① Ⅳ.① 中国版本图书馆CIP数据核字(2014) 第000000号 责任编辑: 王丽平/责任校对: 彭涛 责任印制: 肖钦/封面设计: 陈敬 出版 北京东黄城根北街16号 邮政编码: 100717 https://www.wendangku.net/doc/8e1535139.html, 印刷 科学出版社发行 各地新华书店经销 * 2015年1月第一版开本: 720×1000 1/16 2015年1月第一次印刷印张: 14 字数: 280 000 定价: 88.00元 (如有印装质量问题, 我社负责调换)

几种常见的偏微分方程数值求解问题

一.椭圆型问题 1.1单位圆盘的泊松方程 泊松方程是最简单的椭圆型PDE问题。 该问题的公式为,边界上U=0。 该问题的精确解为 1使用命令行函数 首先必须创建MATLAB函数,使二维几何模型参数化。 M文件circle.m返回单位圆边界点的坐标。该文件内容为:nbs=4; if nargin==0, x=nbs; %边界线段个数 return end d=[ 00 0 0 1 1 1 1 1 1 1 1 00 0 0 ]; bs1=bs(:)’; if find(bs1<1 | bs1>nbs), error(‘Non existent boundary sement number’) end x=zeros(size(s)); y=zeros(size(s)); [m,n]=size(bs); if m==1 & n==1, bs=bs*ones(size(s)); %扩展bs elseif m~=size(s,1) | n~=size(s,1), error(‘bs must be scalar or of same size as s’); end if ~isempty(s), %边界线段1 ii=find(bs==1); x(ii)=1*cos((pi/2)*s(ii)-pi);

y(ii)=1*sin((pi/2)*s(ii)-pi); %边界线段2 ii=find(bs==2); x(ii)=1*cos((pi/2)*s(ii)-(pi/2)); y(ii)=1*sin((pi/2)*s(ii)- (pi/2)); %边界线段3 ii=find(bs==3); x(ii)=1*cos((pi/2)*s(ii)); y(ii)=1*sin((pi/2)*s(ii)); %边界线段4 ii=find(bs==4); x(ii)=1*cos((pi/2)*s(ii)-(3*pi/2); y(ii)=1*sin((pi/2)*s(ii)- (3*pi/2); end 然后用另一函数circleb1.m描述边界条件。 function[q,g,h,r]=circleb1(p,e,u,time) b1=[ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 48 48 48 48 48 48 48 48 49 49 49 49 48 48 48 48 ]; if any(size(u)) [q,g,h,r]=pdeexpd(p,e,u,time,b1); else [q,g,h,r]=pdeexpd(p,e,time,b1); end 现在可以用命令行进行工作: [p,e,t]=initmesh(‘circleg’,’Hmax’,1); error=[];err=1; while err>0.001, *p,e,t+=refinemesh(‘circleg’,p,e,t); u=assempde(‘circleb1’,p,e,t,1,0,1);

偏微分方程上机实验报告.doc

上机实验2:五点差分格式法 偏微分方程(Matlab )实验报告 ——五点差分格式法 一、 实验题目 设G 是形如下图的十字形域,由五个相等的单位正方形组成,用五点差分格式求下列边值问题的数值解: 22 2 21,u u G x y ??+=-???于u=0,于G 二、 实验原理 取定沿X 轴和Y 轴方向的步长1h 和2h ,() 12 22 1 2 h h h =+,作两族与坐 标轴平行的直线:x=i 1h ,y=j 2h ,,0,1,2,i j =±± 若(,i j x y )为正则内点,沿x,y 方向分别用二阶中心差商代替 xx yy u u 和则得 1,1,,1,1 2 212 22[ ]i j ij i j i j ij i j ij u u u u u u f h h +-+--+-+-+ = 特别取正方形网格:12h h h ==,则原差分方程可简化为 2 1,,11,,11()44 ij i j i j i j i j ij h u u u u u f --++-+++= 三、 实验程序 1)function uxy = EllIni2Uxl(x,y) format long ;

uxy = 0; 2)function uxy = EllIni2Uxr(x,y) format long; uxy = y*(2-y); 3)function uxy = EllIni2Uyl(x,y) format long; uxy = 0; 4)function uxy = EllIni2Uyr(x,y) format long; if x < 1 uxy = x; else uxy = 2 - x; end 5)function u = peEllip5(nx,minx,maxx,ny,miny,maxy) format long; hx = (maxx-minx)/(nx-1); hy = (maxy-miny)/(ny-1); u0 = zeros(nx,ny); for j=1:ny u0(j,1) = EllIni2Uxl(minx,miny+(j-1)*hy); u0(j,nx) = EllIni2Uxr(maxx,miny+(j-1)*hy); end for j=1:nx u0(1,j) = EllIni2Uyl(minx+(j-1)*hx,miny); u0(ny,j) = EllIni2Uyr(minx+(j-1)*hx,maxy); end A = -4*eye((nx-2)*(ny-2),(nx-2)*(ny-2)); b = ones((nx-2)*(ny-2),1).*(-1); for i=1:(nx-2)*(ny-2) if mod(i,nx-2) == 1 if i==1 A(1,2) = 1; A(1,nx-1) = 1; b(1) = - u0(1,2) - u0(2,1); else if i == (ny-3)*(nx-2)+1 A(i,i+1) = 1; A(i,i-nx+2) = 1;

相关文档