非线性偏微分方程-偏微分方程数值方法

非线性偏微分方程偏微分方程数值方法非线性偏微分方程偏微分方程数值方法非线性偏微分方程定义:各阶微分项有次数高于一的，该微分方程即为非线性微分方程(一)主要研究内容非线性偏微分方程是现代数学的一个重要分支，无论在理论中还是在实际应用中，非线性偏微分方程均被用来描述力学、控制过程、生态与经济系统、化工循环系统及流行病学等领域的问题。利用非线性偏微分方程描述上述问题充分考虑到空间、时间、时滞的影响，因而更能准确的反映实际。本方向主要研究非线性偏微分方程、H-半变分不等式、最优控制系统的微分方程理论及其在电力系统的应用。1.非线性偏微分方程的研究:我们主要研究偏微分方程解的存在唯一性(和多解性)及稳定性;偏微分方程的初值问题、初边值问题的整体解(包括周期解和概周期解)的存在性及渐近性;平衡解的存在性，尤其是当问题依赖于某些参数时平衡解的分叉结构，以及平衡解的稳定性问题;非线性方程的数值解。2.H-半变分不等式的研究:建立具有极大单调算子扰动的多值(S)型和伪单调型映象的广义度理论，广义不动点指标理论和具有非凸、不可微泛函的非线性发展型H-半变分不等式理论，由此来研究含间断项的非线性偏微分方程。3.最优控制系统的微分方程理论及其在电力系统的应用:主要研究与电力生产有关的控制系统的理论和应用。首先提出了对Banach空间中抽象非线性发展方程所描述的最优控制系统的研究。引进非光滑分析，研究最优控制系统的微分方程，利用变分不等式理论研究多值问题、数值计算等，所获理论成果应用于电力系统的许多最优控制问题(如:电力系统励磁调节器传递函数的辨识、牛顿最优潮流的数学模型等)。(二)研究方向的特色1.变分不等式理论与能量泛函的凸性密切相关，由于现代科学技术的需要，特别是研究自由边界和固体力学问题的需要，传统的方法往往都无法解决这类问题，人们对H-半变分不等式进行研究，研究涉及现代分析及应用、偏微分方程以及科学计算等众多领域中亟待解决和发展的重要课题。2.该研究是现代数学与电力生产的交叉学科研究课题，它对电力生产及管理有着十分重要的理论指导意义和实际应用价值，为控制系统设计、分析和计算都可提供一些重要的理论依据。在应用数学学科的这一研究领域中本课题属于国内外前沿性研究工作。(三)可取得的突破1.深入研究空间、时间、时滞对解的性质的影响，诸如静态解、周期解的存在性、解的存在性、渐近性等问题;寻求它们在含间断项的非线性偏微分方程方面的突破。2.寻求和发现新的处理非单调、非凸不可微能量泛函的方法(如建立Ishikawa迭代序列收敛准则)，建立发展型方程G-收敛准则，寻求可行的光滑方法将算子方程光滑化，创建新的先验估计方法。3.应用现代数学所获得的理论，研究最有控制系统的微分方程，为控制系统设计、分析和计算提供一些重要的理论依据和方法。1747年，法国的达朗贝尔等由弦振动的研究而开创偏微分方程论。1760~1761年，法国的拉格朗日系统地研究了变分法及其在力学上的应用。随机微分方程数值解在随机微分方程数值解这个领域，近几年来国内涉足它的人开始逐渐增多。它也是一门建立在随机分析与微分方程数值解之间的新兴学科。作为一个初学者，我想从它的框架简单谈一下自己的认识，以供讨论。从研究的问题本身来说它主要分为:1随机常微分方程数值方法2随机偏微分方程数值方法3随机延时微分方程数值方法4倒向随机微分方程数值方法仅这四个方面就已经涵盖目前非常重要的一些技术领域的应用。另外从数值方法上分，它可以分为:1强逼近问题2弱逼近问题还有更强的顺向逼近。国内最早涉足这个领域的是山大的彭实戈老师，已经在倒向随机微分方程理论及随机最优控制方面取得了惊人的突破。国外方面，在美国做随机常微分方程的很少(只有Hchurz，lamba几个)，做随机偏微分方如Allen,Cao等等)。在欧洲做随机常微分方程的很多(如Talay，程的较多(Higham，Milstein等)。另外澳洲也有专门研究随机常微分方程的(如Burrage)。随机微分方程(SDE)是a微分方程在哪些一个或更多期限是a随机过程因而造成是本身一个随机过程的解答。一般，SDEs合并空白噪声哪些能被重视作为衍生物苏格兰的植物学家RobertBrown的行动(或熏肉香肠过程);然而，值得一提的是，任意波动的其他类型是可能的，例如跳跃过程(参见[1]).内容1背景1.1术语1.2随机微积分1.3数值解2用途在物理2.1笔记关于"Langevin等式"3用途在可能性和财政数学4解答的存在和独特5参考6参见背景在SDEs的最早期的工作被完成描述苏格兰的植物学家RobertBrown的行动爱因斯坦's著名纸和同时由Smoluchowski。然而，其中一更加早期的工作与苏格兰的植物学家RobertBrown的行动有关相信Bachelier(1900)在他的论文'猜想理论'。这工作被跟随了Langevin.最新Ito和Stratonovich在更加坚实的数学立足处投入了SDEs。术语在物理学，SDEs通常被写当Langevin等式。这些有时缠扰不清称"Langevin等式"即使有许多可能的形式。这些包括包含一个确定部分和一另外任意的一个常微分方程空白噪声期限。第二个形式是福克战斗机Planck等式.福克战斗机Planck等式是描述时间演变的一个偏微分方程概率分布作用.第三个形式是在数学和财务最频繁使用(如下所示)的随机微分方程。这于Langevin形式是相似的，但它在有差别的形式通常被写。这个形式频繁地使用由数学家和在定量财务。SDEs进来二品种，对应于随机微积分的二个版本。随机微积分苏格兰的植物学家RobertBrown的行动或熏肉香肠过程数学上被发现是格外复杂的。熏肉香肠过程non-differentiable;因此，它要求微积分它自己的规则。使用随机微积分的二个版本，Ito随机微积分并且Stratonovich随机微积分.当你应该使用一或其他时，它是有些模棱两可的。方便地，你在解答可能再欣然转换ItoSDE成等效StratonovichSDE和后面成援助;然而，使用的你一定小心当的微积分SDE最初写下时。数值解随机微分方程的特别是数值解和随机偏微分方程相对地讲是一个年轻领域。几乎为常微分方程的解答使用的所有算法为SDEs非常不足将运作，有非常恶劣的数字汇合。用途在物理在物理，SDEs在Langevin形式典型地被写并且被称为"Langevin等式"。例如，一般被结合的套优先处理的SDEs在形式经常被写:那里是套未知数，fi并且gi是任意作用和ηm是，经常被称为的时间的任意作用"噪声命名"。这个形式通常是能用的，因为有变换的标准技术高次等式成数通过介绍新的未知数结合了优先处理的等式。如果gi是常数，系统被认为受叠加性噪声支配，否则它被认为受乘噪声支配。这个期限是有些引入歧途的，因为它来意味一般案件，即使看起来暗示有限的案件，:.叠加性噪声是简单的二个案件。正确解答可能使用平凡经常被发现微积分.特别是，平凡连锁法则微积分能使用。然而，在乘噪声情况下，Langevin等式不是明确定义的个体独自，并且必须指定它是否应该解释Langevin等式作为ItoSDE或StratonovichSDE。在物理，解答主要方法将发现概率分布作用作为时间功能使用等值福克战斗机Planck等式(FPE)。福克战斗机Planck等式是确定的偏微分方程.它告诉怎样概率分布作用及时相似地演变于怎样Schrdinger等式给量子波函数的时间演变或扩散等式给化工集中的时间演变。二者择一地数值解可以获得蒙特卡洛模仿。其他技术包括道路综合化那在比喻画在统计物理之间和量子力学(例如，福克战斗机Planck等式可以被变换成Schrdinger等式通过重新调节几可变物)或通过写下常微分方程为统计片刻概率分布作用。笔记关于"Langevin等式"""在"Langevin等式"是有些不合文法命名原则。每个单独物理模型有它自己的Langevin等式。或许，"Langevin等式"或"伴生的Langevin等式"更将好遵守共同的英国用法。用途在可能性和财政数学记法用于概率论例如(和在概率论的许多应用，财政数学)是轻微地不同的。这个记法做异乎寻常的自然时间的任意作用ηm在物理公式化更加明确。也是用于出版物的记法数字方法为解决随机微分方程。用严密的数学用语，ηm不能仅被选择作为一个通常作用，而是作为a广义函数.数学公式化比物理公式化对待这复杂化以较少二义性。一个典型的等式是形式那里B表示a熏肉香肠过程(标准苏格兰的植物学家RobertBrown的行动)。应该解释这个等式作为一个不拘形式的方式表达对应积分方程上面等式描绘行为连续的时间随机过程xt作为平凡的总和Lebesgue积分式并且Itō积分式.A启发式(但是非常随机微分方程的有用的)解释那在小规模间隔时间长度δ随机过程xt改变它的价值由是的数量通常分布与期望μ(xt,t)δ并且变化σ(xt,t)δ并且是过程的过去行为的独立。这如此是，因为熏肉香肠过程的增加是独立和通常分布。作用μ指漂泊系数，当时σ叫扩散率。随机过程xt叫a扩散过程和通常是aMarkov过程.SDE的正式解释被给根据什么构成解答对SDE。有解答对SDE，一种强的解答和一种微弱的解答的二个主要定义。两个要求过程的存在xt那解决SDE的积分方程版本。二句谎言之间的区别在部下的概率空间(ΩFPr)。一种微弱的解答包括a概率空间并且满足积分方程的过程，而一种强的解答是满足等式的过程和被定义在一个特定概率空间。一个重要例子是等式为几何学苏格兰的植物学家RobertBrown的行动哪些是等式为a的价格的动力学股票在黑Scholes定价财政数学的模型选择。也有更加一般的随机微分方程，系数μ并且σ取决于不仅过程的现值xt，而且在过程的早先价值和可能在其他过程的当前或早先价值也是。在那个案件解答过程，x不是Markov过程，并且它称Itō过程而不是扩散过程。当系数仅依靠礼物和通过价值x定义的等式称随机延迟微分方程。解答的存在和独特和以确定普通和偏微分方程，知道是重要的特定SDE是否有一种解答，并且是否它是独特的。下列是一个典型的存在和独特定理为Itō采取价值的SDEsn-尺寸欧几里德的空间Rn并且由驾驶m-尺寸苏格兰的植物学家RobertBrown的行动B;证明在ksendal(2003年，?5.2)也许被发现。让T0，和让是可测函数为哪些那里存在常数C并且D这样为所有t?[0,T]和所有x并且y?Rn的地方让Z是独立的一个随机变量σ-引起的代数Bs,s?0，和与有限二次矩:然后随机微分方程或初值问题xt=Z;有Pr-几乎肯定独特t-连续的解答(t,ω)|?xt(ω)这样x是适应对滤清FtZ引起Z并且Bs,s?t和参考adomian，乔治(1983)。随机系统数学在科学和工程学(169)。奥兰多，FL:学术出版社公司。adomian，乔治(1986)。非线性随机操作员等式.奥兰多，FL:学术出版社公司。adomian，乔治(1989)。在物理的非线性随机系统理论和应用数学和它的应用(46)。Dordrecht:Kluwer学术出版者小组。ksendal，BerntK。(2003).随机微分方程:介绍以应用.柏林:Springer。国际标准书号3-540-04758-1.Teugels，J。并且SundB。(eds。)(2004)。保险统计计算科学百科全书.Chichester:威里，523-527。C.W.Gardiner(2004)。随机方法手册:为物理、化学和自然科学.Springer，415。托马斯?Mikosch(1998)。基本的随机微积分:以财务视线内.新加坡:世界科学出版，212。国际标准书号981-02-3543-7.Bachelier，L.，(1900)。Théoriedelaspeculation(用法语)，PhD论文.NUMDAM:用英语在1971书'股市'Eds的任意字符。P.H.Cootner。高性能科学计算研究一、研究内容一般地，构成实际应用物理过程的各个不同阶段的物理模型，可分别由不同类型的时间相关或无关的偏微分方程在给定的物理区域上描述。如何针对不同偏微分方程的问题设计合适的网格和离散格式，如何设计可扩展的并行算法及其并行实现技术，在离散网格上给出方程的近似解，是我们研究的两个主要方面。本项目的研究以科学计算的共性问题为核心，包括具有最优复杂性的计算方法研究和能发挥计算机浮点计算峰值性能的实现技术研究，同时应用本项目科学计算的共性问题的研究成果，解决一批我国具有重大需求的科学计算问题。1.创新计算方法的基础理论研究计算数学是研究可在计算机上运行的数值算法的构造及其数学理论的学科。过去五十多年科学计算发展的历史表明:基础计算方法的重要突破如有限元方法、多重网格方法、快速傅里叶变换等都极大地改变了科学计算的面貌。我们将研究有限元新型算法包括多重网格与区域分解算法、均匀化多尺度算法、自适应高精度算法和各类方法的耦合，动力系统的保结构算法，守恒律高分辨率差分格式，各类快速算法包括非规则网格的快速傅里叶变换等，同时研究新的应用领域大规模高速集成电路中电磁信息计算中的计算方法。研究重点在并行自适应算法与理论，保结构计算方法的理论与应用，大规模高速集成电路中电磁信息计算。1.1并行自适应算法与理论这里自适应方法主要是指网格自适应方法，是一类渗透到了偏微分方程数值解、非线性逼近论、偏微分方程约束的最优工程设计、网格产生等科目研究的方法。现在网格自适应方法主要分为三种主要的类型，分别叫做h-方法、p-方法和r-方法。其中h-方法是对网格进行自适应的局部加密和稀疏化，p-方法是在网格的不同位置使用不同的基函数，r-方法是进行网格点的重新分布，又叫做移动网格方法。将h-方法和p-方法结合可以得到h-p方法，也可以将r-方法和p-方法结合得到r-p方法。网格自适应方法最根本的目标在于使用最少的计算资源来解决问题，从而可以在现有的硬件资源条件下扩大计算的规模和提高计算的精度。针对当前国际研究发展的趋势和本项目应用问题的需求，我们主要的研究内容集中在下面的二个方面:网格方法在偏微分方程数值解中的应用研究摘要:该文的主要目的是研究无网格方法,并将其应用于偏微分方程的数值解过程中.与传统的网格方法不同,无网格方法的核心是用"点云"离散求解区域,并基于当地点云离散结构,引入二次极小曲面逼近空间导数.该文先以代表定常不可压位势绕流的Laplace方程为例,研究了Laplace方程的无网格离散形式,并运用GMRES高效算法对其快速求解,数值模拟了典型的圆柱绕流;并通过不同点云尺度的数值模拟,显示出点云尺度对计算精度的影响.在此基础上,将该方法推广应用到解算Euler方程组.针对守恒型Euler方程组的无网格离散形式,借鉴非结构网格方法附加耗散模型,采用五步Runge-Kutta显式时间推进格式求解.并且基于点云离散结构,引入了当地时间步长、残值光顺等加速收敛技术,数值模拟了对称和非对称翼型绕流,获得较好的计算结果.该文还对基于点云结构的无网格计算软件的面向对象设计模式进行了研究,着重于提高软件的复用性和Matlab偏微分方程工具箱简介1.概述本文只给出该工具箱的函数列表，读者应先具备偏微分方程的基本知识，然后根据本文列出的函数查阅Matlab的help，便可掌握该工具箱的使用。2.偏微分方程算法函数列表adaptmesh生成自适应网络及偏微分方程的解assemb生成边界质量和刚度矩阵assema生成积分区域上质量和刚度矩阵assempde组成偏微分方程的刚度矩阵及右边hyperbolic求解双曲线型偏微分方程parabolic求解抛物线型偏微分方程pdeeig求解特征型偏微分方程pdenonlin求解非线性型微分方程poisolv利用矩阵格式快速求解泊松方程3.图形界面函数pdecirc画圆pdeellip画椭圆pdemdlcv转化为版本1.0式的*.m文件pdepoly画多边形pderect画矩形pdetool偏微分方程工具箱的图形用户界面4.几何处理函数csgchk检查几何矩阵的有效性csgdel删除接近边界的小区decsg将固定的几何区域分解为最小区域initmesh产生最初的三角形网络jigglemesh微调区域内的三角形网络poimesh在矩形区域上产生规则的网络refinemesh细化三角形网络wbound写一个边界描述文件wgeom写一个几何描述文件pdecont画轮廓图pdemesh画偏微分方程的三角形网络pdeplot画偏微分方程的三角形网络pdesurf画表面图命令5.通用函数pdetriq三角形单元的品性度量poiasma边界点对快速求解泊松方程的"贡献"矩阵poicalc规范化的矩阵格式的点索引poiindex规范化的矩阵格式的点索引sptarn求解一般的稀疏矩阵的特征值问题tri2grid由三角形格式转化为矩形格式《偏微分方程中多尺度问题的数值解法》偏微分方程数值方法理论及其应用、有限元方法、多重网格法与区域分解法"偏微分方程数值求解中的自适应网格方法研究"人工边界方法:无界区域上的偏微分方程数值解"有限元高精度理论及算法"、"具有奇异解的偏微分方程的数值解法"、"无界域上偏微分方程的数值解法"、"多尺度有限元方法及其快速算法"、"快速数值计算算法及软件"偏微分方程数值解法2所谓的偏微分方程(PDE)是指含两个以上自变量的微分方程。偏微分方程的求解一般说来太过复杂，所以现在还没有一个对所有偏微分进行求解的理论，所谓的求解偏微分方程也只是对某些人们比较熟悉的类型进行求解。对于一个形如A(x,y)Uxx+B(x,y)Uxy+C(x,y)Uyy=f(x,y,U,Ux,Uy)inΩ的偏微分方程其中Ω是给定的平面有界区域。如果B^2-4AC0椭圆型B^2-4AC=0抛物线型B^2-4AC0双曲线型如果ABC是常数，方程被称为拟线性方程。以上三类方程，人们有较成熟的解法。这三类方程也有物理意义，比如椭圆型方程常见于电磁场的分布，抛物线型方程常见于扩散，双曲线型常见于波动，后两者还常会带有对时间的求导项。这些方程，往往在一定的条件下才能有定解:Dirichlet条件，又称第一类边界条件，设定初值Neumann条件，又称第二类边界条件，设定边值条件很多情况下，两者都有，称为混合边界条件。我的课题中涉及到一个物质随着流动相在色谱柱里运动的方程，能够描述物质浓度波在柱内的运动和变形，因此会包括一阶时间项和二阶空间项，有个专有名词--对流扩散方程，是种抛物线型和双曲线型的混合型方程。偏微分方程数值解法差分方法有限元方法拟谱方法自适应格点方法小波分析方法解偏微分方程解决的方向:微分算子的计算或表达时间的差分离散边界的处理收敛性分析误差的估计稳定性分析微分算子的自适应计算时间和空间的自适应计算差分法从定解问题的微分或积分形式出发,用数值微商或数值积分公式导出相应的线性代数方程组.构造逼近微分方程定解问题的差分格式:直接差分化法,积分插值法以及有限体积法或广义差分法.差分解的存在唯一性,收敛性以及稳定性的研究.这些理论问题为对差分解作出先验估计.基于极值定理以及能量不等式作估计.有限元法从定解问题的变分形式出发,用Ritz-Galerkin方法导出相应的线性代数方程组.中文译名?偏微分方程的多尺度小波方法本书系《小波分析及其应用》第6卷，是一本论文集。小波分析是目前国际上公认的最新时-频分析工具，由于其具有自适应性和数学显微镜性质，而成为众多学科共同关注的焦点。从数学角度讲，小波分析对函数逼近、调和分析、统计学、微分和积分方程的数值解等均产生直接的影响。本书作为小波分析与一般偏微分方程(PDE-partialDifferentialEquation)技术的桥梁，将多尺度分解的概念引入到了PDE的数值求解，可有效的分析较复杂问题。书中内容分为6部分:(1)回顾了基于多层预调节及多网格技术的有限元法，多尺度空间分解框架，域内椭圆形问题的多尺度解法。(2)快速小波算法(压缩与自适应方面):D维二阶椭圆形PDE的自适应解的小波配置方法，求解非线性PDE的自适应小波分析，基于小波包最佳基的动态自适应概念在对流扩散PDE中的应用，求解椭圆算子方程中的非线性近似与自适应技术。(3)积分方程的小波求解，包括强椭圆边界积分方程的多尺度Galerkin法。(4)小波多尺度求解PDE的软件工具与数值实例。(5)多尺度分析在湍流中的应用。(6)偏微分算子的小波分析。本书收集的14篇论文代表了当前小波在偏微分方程应用中的最新进展，可供小波理论及应用、PDE等应用数学领域的科研人员学习参考。(力学系马坚伟)小波分析方法小波分析方法解偏微分方程思路:Galerkin方法为基础;半群方法为基础.基于偏微分方程或积分方程的信号处理,流体动力学的问题就能用此方程描述.这些问题解的特征为光滑的(smooth),非振荡的(non-oscillatory),shock.方法为:算子和解投影到小波基上.基函数的消失矩特性使得解和算子能够稀疏表达,因此就能给出快速,自适应算法.这些算法基于在光滑区域用较少的小波系数,在奇异区域得用较多的小波系数.解这类方程重要的一步为时间的离散.因为进化方程的扩散项,标准的显格式容许小的时间步长.另外,隐格式容许大的时间步长,但在每一步得解线性方程组,这就给应用带来了困难.B.Alpert,G.Beylkin,Tchamitchian(1990-2005)用的方法:Wavelet-Galerkinmethod,Taylor-Galerkinmethod,配点方法,非标准小波表示.JohnWeiss用小波Galerkin方法(Daubechies,1992,1993).用的是时间差分,空间离散.计算比较复杂,但精度好.小波Galerkin方法Galerkin配点方法:通过投影将连续算子离散化为矩阵形式,此方法的困难在于二重积分的数值计算;为解决这困难,研究者提出了函数基用小波基,此方法被称为小波Galerkin方法.在作数值逼近计算时,因为用了小波基,因此很多算子可用稀疏矩阵表示,那么小波Galerkin方法就为作快速数值计算提供了算法.总的来说,小波Galerkin方法在作逼近分析时比Adomian分解方法更可靠,在作数值逼近计算时比Galerkin方法速度更快.算法复杂性为另外,得分析稳定性;不同小波基础的误差估计;时间空间的自适应.Legendre多小波的非标准表示的优点:算子矩阵稀疏;子区间元素相同;维数低;可线性化非线性项.Legendre多小波不连续,微分算子的处理方法:通过尺度方程导出系数方程组,解此方程组可得到算子矩阵;用传统的弱导数通过积分计算算子矩阵.此小波处理边界有优势.边界的处理?构造多分辨分析,使得小波基满足边界条件.用插值小波,配点方法.变系数的处理?时间空间的自适应?应用小波分析求解微分方程研究作者:来源:信息与计算科学系责任编辑:xinxi课题主持人:孙涛项目组成员:孙涛、李震、武斌、赵燕项目研究时间:2010.5-2012.5项目研究内容:主要研究应用小波分析进行微分方程的求解特别是偏微分方程的数值求解。预期目标是研究应用小波理论进行微分方程求解的已有成果，分析比较各种方法在理论与应用上的优缺点，同时对其在适用范围、计算精度、计算复杂性、收敛性以及稳定性等方面进行对比，从而有针对性的对各种方法进行改进或完善;对将小波方法应用于偏微分方程数值求解的数学思想进行研究，形成基本的小波方法;对小波方法求解偏微分方程的小波基的特点进行分析，明确用于偏微分方程数值解法的小波基的数学特性，设计用小波方法求解偏微分方程的一般数学方法。研究成果形式:论文和研究报告。偏微分方程是需要常微分方程和随机微分(随机过程)两门课做基础的需同时具备边界条件和初始条件。只给边界条件，一般无法解。如题目无初始条件，可自定(设)一些初始条件。只有范围的结果,但不能求出精确的解.给了边界就能.稳定性分析是针对某一特定的差分算法来说的。而并不是对偏微分方程来说的。一般是用Fouier分析的办法来做。你可以看一下余德浩,汤华中编的科学出版社出版的"微分方程数值解法"里面216页有一些相关的东西。比较常用的差分算法有Lax_Wendroff格式以及MacCormack格式。另外，你如果想要解析解的话，估计可能要用特征线法。或者分离变量法看一下。微分方程数值解?NumericalSolutionsofDifferentialEquations课程编号:S080800XJ001课程属性:学科基础课学时/学分:40/2预修课程:高等数学(包括数学分析与线性代数)、数学物理方程、计算方法、程序设计。教学目的和要求:本课程为数学、物理、力学、化学及工程科学等专业硕士研究生的选修课。主要讲授常微分方程和偏微分方程差分方法的算法、稳定性和收敛性理论，内容包括常微分方程初值与边值问题的数值解法，抛物型、双曲型及椭圆型偏微分方程的差分方法等。通过本课程学习，希望学生掌握数值求解微分方程的一些基本方法，为进一步学习计算数学的专业课或在各自的专业工作中应用科学计算这一重要研究手段打下基础。内容提要:第一章常微分方程初、边值问题数值解法Euler方法;Runge-Kutta方法;线性多步方法;稳定性、收敛性和误差估计;常微分方程边值问题的数值方法。第二章抛物型方程的差分方法差分格式建立的基础;显式、隐式差分格式;差分格式的稳定性和收敛性;高维抛物型方程的差分方法;交替方向隐式差分方法。第三章双曲型方程的差分方法一维双曲型方程的特征线方法;一阶线性双曲型方程(组)的差分方法;双曲型守恒律方程及守恒型差分格式;二阶波动方程的差分方法。第四章椭圆型方程的差分方法Poisson方程第一边值问题的差分方法;Poisson方程的有限体积方法;差分方法的收敛性和误差估计;椭圆型差分方程的迭代解法;多重网格方法。教材余德浩、汤华中，《微分方程数值解法》，科学出版社，北京，2002。张文生，《科学计算中的偏微分方程有限差分法》，科学出版社，北京，2006。主要参考书:[1]J.W.Thomas.NumericalPartialDifferentialEquations:FiniteDifferenceMethods.Springer-VerlagNewYorkInc.1995.[2]胡健伟、汤怀民，《微分方程数值方法》，科学出版社，北京，1999。偏微分方程数值解的两类主要方法:差分方法和有限元方法二课程性质、目的与任务《偏微分方程数值解》是信息与计算科学专业的一门专业课，学生通过学习一些典型、通用的偏微数值方法，掌握用差分法，有限元法求解偏微分方法的基本理论，理解这些方法构造的基本思想，学会编制差分法和有限元法的计算程序，同时通过学习一些基本概念和基本理论(如稳定性、收敛性、误差估计培养一定的理论分析能力。等)三教学基本内容与基本要求教学基本内容包括:1.抛物型方程的有限差分方法2.双曲型方程的有限差分方法3.椭圆型方程的有限差分方法4.变分原理5.有限单元法6.有限元方法理论基础教学基本要求1.掌握差分法和有限元法的基本理论2.了解用差分法和有限元法计算偏微分方程的误差估计方法3.能独立编制差分法和有限元法的计算程序多尺度问题中的偏微分方程数值方法《微分方程数值解》第一章绪论一、学习目的通过本章的学习，了解偏微分方程中的三大类方程，以及偏微分的一些基本概念。计划8学时。二、课程内容第一节数学物理方程中的三大类方程(一)抛物型方程典型方程:热传导方程，由空间物体的热传导问题导出。利用物理中传热学的傅里叶实验定律。(二)双曲型方程典型方程:波动方程，由两端固定的细弦振动导出。利用胡克定律、牛顿第二定律等。(三)椭圆型方程典型方程:调和方程(Laplace方程)，由静电场的电位势或没有热源的热传导等导出。第二节数学物理方程中的基本概念何为线性的或非线性的，给出一个方程怎么判断它是哪类方程，定解问题的三种提法等。三、重点、难点提示和教学手段本章重点是三类方程的导出和偏微分方程中的基本概念。难点是导出过程的理论推导。四、思考与练习掌握、吸收所学知识。第二章常微分方程初值问题数值解法一、学习目的通过本章的学习，对常微分方程初值问题的几个典型方法了解、掌握，并能编写程序。计划8学时。二、课程内容2.1欧拉法(一)欧拉法的格式:用差商代替微商。(二)收敛性研究通过分析截断误差确定格式的收敛速度。(三)稳定性研究格式对初值误差的连续依赖性。2.2梯形法、隐式格式的迭代计算用梯形公式近似计算积分得到常微分方程的梯形公式，而且是一个隐式格式。估算梯形法的整体截断误差。2.3单步法、Runge-Kutta法用泰勒级数构造一般的单步法，几种不同的Runge-Kutta法，以及各自的优缺点。其中经典的四阶Runge-Kutta法尤为重要。2.4线性多步法用Lagrange插值近似小分割上的曲线，得到线性多步法。Adams外插、内插公式等。2.5误差的事后估计法、步长的自动选择何为误差的事后估计法，以及如何利用事后估计法得到的截断误差作为步长h自动选择的标准。2.6高阶常微分方程(组)的数值方法怎样把高阶微分方程转化为一阶的方程组，然后怎么对方程组利用前面所介绍的方法进行近似计算。三、重点、难点提示和教学手段本章重点是利用各种方法求方程的近似解。难点是方法的推导以及局部和整体截断误差的估计。四、思考与练习复习所学内容，计算课堂上没有推导的几种格式的截断误差，然后完成布置的作业。第三章抛物型方程的差分格式一、学习目的掌握有关差分格式以及稳定性的一些基本概念，会构造差分格式并可用两种方法分析差分格式的稳定性。了解差分格式稳定性的定义及其含义。计划14学时。二、课程内容3.1差分格式建立的基础对所考虑的方程的初边值问题进行网格剖分，建立差分格式。学习三种差商代替微商的方法。学会用算子形式表示差分格式。3.2显示差分格式一维常系数热传导方程的古典显式格式，以及系数依赖于x的一维热传导方程的显式格式。计算各自的截断误差。3.3隐式差分格式由向后差商得到古典隐式格式，推导常用的Crank-Nicolson隐式格式和加权的六点隐式格式，知道前两种是六点加权隐式格式的特殊形式。系数依赖于x,t的一维热传导方程的隐式格式。3.4解三对角形方程组的追赶法在求解隐式差分方程时形成一个线性代数方程组，它的系数矩阵是三对角形矩阵，因此要学会用追赶法求这类方程组，分为追和赶两步。3.5差分格式的稳定性和收敛性学习-图方法、矩阵法、Fourier级数法(VonNeumann方法)分析差分格式的稳定性。重点用矩阵法和VonNeumann方法分析前面学习的几种差分格式，并比较后两种方法的优劣。对于收敛性利用Lax等价性定理转化为对稳定性的研究。3.6非线性抛物型方程的差分解法举例包括Richtmyer线性方法和Less三层差分格式。对于Less三层差分格式要对第二层利用其他方法求出。3.7二维抛物型方程的差分格式初边值问题要进行三维方向的网格剖分，其中方法与一维的类似，也有显式和隐式之分，以及稳定性分析等，其中显式简单，但效果没有隐式好。3.8交替方向的隐式差分格式(ADI格式)为了提高精度和满足无条件稳定的差分格式，把每一时间层的计算分成几步进行，而使每步具有一维格式的特点，提出以下几种格式:Peaceman-Rachford格式、Douglas-Rachford格式、Mitchell-Fairweather格式等。三、重点、难点提示和教学手段本章重点是利用各种方法求方程的近似解，用矩阵法和VonNeumann方法进行稳定性分析。难点是格式的理论导出和稳定性分析。四、思考与练习复习所学内容，编写一定的程序，用两种方法进行稳定性分析，然后完成布置的作业。第四章椭圆型方程的差分格式一、学习目的掌握椭圆型方程的五点、九点差分格式，掌握极值原理，收敛性分析和误差估计。计划14个学时。二、课程内容4.1正方形区域中的Laplace方程Dirichlet边值问题的差分模拟对Dirichlet边值问题从x和y轴方向进行网格剖分得到Laplace方程的五点差分格式。然后转化为解一个线性矩阵。4.2Neumann边值问题的差分模拟由于边值问题通过告诉它的法向量在边值的值，这样关键就是如何把这个条件转化为边值上的解。利用中心差商代替微商把导数边值转为一般的边值条件。4.3混合边值问题区域的一部分是Dirichlet条件而另一部分是Neumann条件，那么对于Neumann条件利用类似上节的方法处理边值问题。4.4非矩形区域当区域不是规则的矩形时，我们对这种区域的邻接边界的内部结点需要特别的处理，它到边界的距离可以是非整数倍的分割。也可以得到Laplace方程的五点差分格式，它是前面五点格式的推广。4.5极坐标形式的差分格式有的时候所求区域是圆环、环形域或扇形域，采用极坐标形式更为方便，此时应该把一般的Poisson方程转化为极坐标的形式。在极坐标情况下会出现奇异点，故需要附加条件，对它需特别处理。4.6矩形区域上的Poisson方程的五点差分逼近的敛速分析利用极值原理分析五点差分格式的敛速估计。4.7一般二阶线性椭圆型方程差分逼近及其性质研究通过一些实例学习二阶线性椭圆型方程的差分格式。4.8椭圆型差分方程的迭代解法由于前面介绍的各种边值问题的差分格式最终都是解一个大型的线性方程组。那么怎样求这个大型的线性方程组?本节介绍三种迭代法(Jacobi迭代、Gauss-Seidel迭代、超松弛迭代)。通过比较Jacobi迭代和Gauss-Seidel迭代的敛速发现，Gauss-Seidel迭代是Jacobi迭代的两倍。虽然前两者都收敛，但是他们的速度还是比较慢，如果选择适当的松弛因子，利用超松弛方法可以大大提高敛速，故如何选择最佳松弛因子是关键。4.9多重网格法简介学习为何引入多重网格法，有哪些优点?包括二重网格法和多重网格法等。三、重点、难点提示和教学手段本章重点是:差分格式的建立，极值原理及数值解的收敛性分析。教学难点:边界条件的处理及非均匀部分差分格式的建立。四、思考与练习复习所学内容，编写一定的程序，然后完成布置的作业。第五章双曲型方程的差分格式一、学习目的掌握一阶拟线性双曲型方程(组)的特征线法，一阶双曲型方程(组)的差分方法，以及二阶线性双曲型方程的差分方法。计划12个学时。二、课程内容5.1一阶拟线性双曲型方程的特征线法对于一阶(拟)线性双曲型方程，通过一条特征曲线，把一个偏微分问题转化为常微分问题，然后再对此常微分方程进行近似求解，就可以求出原问题的近似解。给出这样的方程要知道怎么求它的特征曲线、特征方程以及特征关系。5.2一阶拟线性双曲型方程(组)的特征线法对于一阶(拟)线性双曲型方程组，首先求出它的正规形式，以及它的两个特征曲线、特征方程和特征关系。同上节一样也是把偏微分方程沿着特征方向转化为常微分方程组的形式。再用欧拉法对常微分方程组近似求解。5.3一阶双曲型方程的差分格式如果通过向前差商代替对t方向的微商，用中心差商代替对行x方向的微商,经过验证发现这是一个恒不稳定的格式，所以通过改进此格式得到Lax-Friedrichs格式。根据方程系数的不同对x方向的微商向前差商或向后差商就是Courant-Isaacson-Rees格式。如果对时间层进行中心差商代替就是跳蛙格式。还有Lax-Wendroff格式和隐式的Crank-Nicolson格式。5.4一阶双曲型方程组的差分格式类似一阶双曲型方程的差分格式，也有Lax-Friedrichs格式和Courant-Isaacson-Rees格式，以及Courant-Friedrichs-Lewy条件。5.5二阶双曲型方程的差分格式有显式格式和隐式格式之分，但是由于双曲型方程的初值问题比较复杂，因此对告诉初始时刻速度的初始条件需要像处理Neumann问题进行差商代替微商。三、重点、难点提示和教学手段本章重点是:一阶拟线性双曲型方程(组)的特征线法，一阶双曲型方程(组)的差分方法，以及二阶线性双曲型方程的差分方法。教学难点:一阶拟线性双曲型方程组的特征线法的导出。四、思考与练习复习所学内容，编写一定的程序，然后完成布置的作业。第六章非线性双曲型守恒律方程的差分格式一、学习目的掌握何为双曲型守恒律、弱解的定义，和几种典型的差分格式:守恒型差分格式，单调差分格式，TVD差分格式等。计划8个学时。二、课程内容6.1非线性双曲型守恒律简介、弱解的定义怎样判断一个方程组的对应Jacobi矩阵的特征值和特征向量决定此方程组是(严格)双曲型守恒律的?以及什么是弱解，为什么需要提出弱解的概念，有什么优点?6.2一阶拟线性双曲型方程(组)的特征线法把一阶的线性双曲型方程中的Lax-Friedrichs格式和Lax-Wendroff格式推广到双曲型守恒律方程，可以证明它们都是守恒型差分格式。6.3单调差分格式此前的守恒型差分格式虽然收敛到弱解，但是不能保证极限是唯一物理解，所以提出单调差分格式，这种格式若收敛，则收敛到唯一物理解。给出满足什么条件才是单调差分格式。另外这种格式只有一阶精度，对于高精度还在研究中。6.4TVD差分格式由于前面提到的单调差分格式精度不高，所以为了能够得到精度较高且能得到唯一物理解的差分格式，由A.Harten于1983年提出了总变差减少差分格式(TVD)。本节给出什么格式是TVD格式，以及保单调格式，并且证明前面讨论的差分格式在一定条件下都是TVD差分格式。6.5对一维方程组的推广把Lax-Wendroff格式和Lax-Friedrichs格式可以推广到一维方程组的情况。三、重点、难点提示和教学手段本章重点是:课本中对双曲型守恒律方程提出的各种差分方法。教学难点:各种差分格式的理论导出。四、思考与练习复习所学内容，然后完成布置的作业。第六章有限元方法简介一、学习目的对有限元方法有所了解，知道它提出的理论依据，以及它的优缺点。计划4个学时。二、课程内容7.1二阶常微分方程边值问题的有限元解法对于常微分边值问题，介绍如何利用有限元方法把边界条件当成最小化泛函的一部分，以及怎样在给定的集合上某一个函数类中找到一个使泛函I达到极小值的函数。7.2偏微分方程边值问题的有限元解法与常微分边值问题类似，也是归结为泛函求极小的一种解法，不过偏微分问题相对