偏微分方程数值解

偏微分方程数值解_图文_图文.ppt1偏微分方程数值解5/9/2024参考数目GeorgeJ.Haltiner,RogerTerryWilliams,NumericalPredictionandDynamicMeteorology(2ndEdition),theUnitedStatesofAmerica,1979.2.CurtisF.GeraldandPatrickO.,AppliedNumericalAnalysis,PersonEducation,Inc.,2004.3.EugeniaKalnay,AtmosphericModeling,DataAssimilationandPredictability,thepressSyndicateoftheUniversityofCambridge,2003.4.AriehIserles,AFirstCourseintheNumericalAnalysisofDifferentialEquations,CambridgeUniversityPress,1996.5.李荣华，冯国忱.微分方程数值解.北京：人民教育出版社，1980.6.徐长发，李红.实用偏微分方程数值解法.华中科技大学出版社，2003.7.沈桐立，田永祥等.数值天气预报.北京：气象出版社，2007.2偏微分方程数值解5/9/2024数值天气预报—PDE数值解挪威气象学家V.Bjerknes（1904）提出数值预报的思想：通过求解一组方程的初值问题可以预报将来某个时刻的天气—思想；L.F.Richardson(1922)：开创了利用数值积分进行预报天气的先例，由于一些原因（如，计算稳定性问题“Courant,1928”）并没有取得预期的效果—尝试；Charney,Fjortoft,andVonNeumann(1950),借助于Princeton大学的的计算机(ENIAC)，利用一个简单的正压涡度方程（C.G.Rossby,1940）对500mb的天气形式作了24小时预报---成功；3偏微分方程数值解5/9/2024TheElectronicNumericalIntegratorandComputer(ENIAC).4偏微分方程数值解5/9/2024常微分方程的数值解大气科学中常微分方程和偏微分方程的关系1.大气行星边界层(近地面具有湍流运动特性的大气薄层,1—1.5km),埃克曼（V.W.Ekman）（瑞典）螺线的导出；2.1963年，美国气象学家Lorenz在研究热对流的不稳定问题时，使用高截断的谱方法，由Boussinesq流体的闭合方程组得到了一个完全确定的三阶常微分方程组，即著名的Lorenz系统。5偏微分方程数值解5/9/2024Lorenz系统dx/dt=a(y-x)dy/dt=x(b-z)-ydz/dt=xy-cz其中，a=10,(Prandtlnumber);b=28(Rayleighnumber);c=8/3;(x,y,z)_0=(0.01;0.01;1e-10)6偏微分方程数值解5/9/20247偏微分方程数值解5/9/20248偏微分方程数值解5/9/20249偏微分方程数值解5/9/202410偏微分方程数值解5/9/2024Franceshini将Navier-Stokes方程截断为五维的截谱模型如下：11偏微分方程数值解5/9/2024欧拉法—折线法常微分方程能直接进行积分的是少数，而多数是借助于计算机来求常微分方程的近似解；有限差分法是常微分方程中数值解法中通常有效的方法；建立差分算法的两个基本的步骤：1.建立差分格式，包括：a.对解的存在域剖分；b.采用不同的算法可得到不同的逼近误差—截断误差（相容性）；c.数值解对真解的精度—整体截断误差（收敛性）；d.数值解收敛于真解的速度；e.差分算法—舍人误差（稳定性）.12偏微分方程数值解5/9/20242.差分格式求解将积分方程通过差分方程转化为代数方程求解，一般常用递推算法。

在常微分方程差分法中最简单的方法是Euler方法，尽管在计算中不会使用，但从中可领悟到建立差分格式的技术路线，下面将对其作详细介绍:13偏微分方程数值解5/9/2024差分方法的基本思想“就是以差商代替微商”考虑如下两个Taylor公式：（1）（2）从（1）得到：14偏微分方程数值解5/9/2024从（2）得到：从（1）-（2）得到：从（1）+（2）得到：15偏微分方程数值解5/9/2024对经典的初值问题满足Lipschitz条件保证了方程组的初值问题有唯一解。16偏微分方程数值解5/9/2024一、算法构造：0tuT1.在求解域上等距离分割：2.在有：微分方程的精确解差分方程的精确解17偏微分方程数值解5/9/20243.应用时采用如下递推方式计算：4.例题对初值问题用Euler法求解，用即，18偏微分方程数值解5/9/20245.Euler法的几何意义0t在递推的每一步，设定过点作的切线，该切线的方程为：即：19偏微分方程数值解5/9/2024二、误差分析构造算法后，这一算法在实际中是否可行呢？也就是说是否使计算机仿真而不失真，这还需要进一步分析。1.局部截断误差--相容性为了分析分析数值方法的精确度，常常在成立的假定下，估计误差这种误差称为“局部截断误差”，如图。局部截断误差是以点的精确解为出发值，用数值方法推进到下一个点而产生的误差。20偏微分方程数值解5/9/2024整体截断误差是以点的初始值为出发值，用数值方法推进i+1步到点，所得的近似值与精确值的偏差：2.整体截断误差—收敛性称为整体截断误差。21偏微分方程数值解5/9/2024特例，若不计初始误差，即则即3.舍入误差—稳定性假设一个计算机仅表示4个数字（小数点后面），那么计算22偏微分方程数值解5/9/2024我们的要求是：最初产生的小误差在以后的计算中虽然会传递下去，但不会无限制的扩大，这就是稳定性所描述的问题。下面引进稳定性的概念：设由初值得到精确解，由初值得到精确解，若存在常数和充分小的步长使得则称数值方法是稳定的。tu023偏微分方程数值解5/9/2024计算例题其解析解为：x=00.20000.40000.60000.80001.0000y=1.00001.20001.37331.53151.68111.826924偏微分方程数值解5/9/202425偏微分方程数值解5/9/2024三、改进的Euler法将微分方程在区间上积分，得到用梯形法计算积分的近似值，有于是这是一个隐式格式，一般需要用迭代法来求，而用显式的Euler法提供初值。26偏微分方程数值解5/9/2024为了简化计算的过程，在此基础上进一步变为如下算法：此式称为“改进的Euler法。接下来讨论其几何意义预估校正其局部截断误差为这个问题将在下节讨论。27偏微分方程数值解5/9/2024tu028偏微分方程数值解5/9/2024Euler法、改进的Euler法和解析解的比较29偏微分方程数值解5/9/2024四、（龙格-库塔）Runge-Kutta方法简单的Euler法是建立在Taylor级数的一项展开；改进的Euler法是以两项Taylor级数为基础建立的，如：如果我们截取Taylor级数的更多项会得到什么样的求解方法呢？两个德国数学家（C.Runge&M.kutta）以这种思想为基础建立了求解微分方程的龙格-库塔方法。它是常微分方程数值解法中使用最为广泛的方法之一。30偏微分方程数值解5/9/2024一般地，一个K阶的Runge-Kutta方法可用下面的公式表示：其中，是待定的加权系数，是待定的系数。Euler法就是的R-K法。其系数的确定如下：将展开成的幂级数，并与微分方程的精确解在点的Taylor展开式相比较，使两者的前项相同，这样确定的R-K法，其局部截断误差为，根据所得关于待定系数的方程组，求出它们的值后代入公式，就成为一个阶R-K方法。31偏微分方程数值解5/9/2024例题以二阶R-K法为例说明上述过程把代入中，有32偏微分方程数值解5/9/2024经比较得到取为自由参数：从而得到不同的但都是二阶的R-K方法，对应的有中点法、Heun（亨）法以及改进的Euler法。33偏微分方程数值解5/9/2024基于相同的过程，通过比较五次Taylor多项式，得到更加复杂的结果，给出了包含13个未知数的11个方程。得到多组系数，其中常用的是以下四阶R-K法：改进的Euler法、R-K法以及解析解的比较：34偏微分方程数值解5/9/202435偏微分方程数值解5/9/2024五、线性多步（LinearMultistepMethod）法1.预备知识：插值多项式插值是离散函数逼近的重要方法，利用它可通过函数在有限个点处的取值状况，估算出函数在其他点处的近似值。

从几何上理解：对一维而言，已知平面上n＋1个不同点，要寻找一条n次多项式曲线通过这些点。插值多项式一般常见的是拉格朗日插值多项式。2.气象应用不均匀站点上的气象要素数据均匀网格点上的数据插值3.拉格朗日插值多项式拉格朗日插值多项式逼近可能是求插值节点不均匀的插值多项式的最简单的方法。实验观察结果或原始测量数据的分布通常是非均匀的。例如，四个点可以确定一个三次多项式，其拉格朗日形式为：36偏微分方程数值解5/9/20244.Adams-Bashforth（阿达姆斯—贝雪福斯）公式首先，用以下四个点对进行三次Langrage插值：则于是，有容易算出，例如，我们可以算得*2*137偏微分方程数值解5/9/2024将（*2）代入（*1）得到Adams-Bashforth公式：基于同样的计算过程，可以得到另外一个计算公式：这称为Adams-Moulto公式。预估校正38偏微分方程数值解5/9/2024偏微分方程数值解主讲：王曰朋39偏微分方程数值解5/9/2024一、区域的离散1.2.3.40偏微分方程数值解5/9/2024则函数可表示为：二、1.（一维）一、二阶导数的有限差分近似表达式41偏微分方程数值解5/9/20242.（二维）一、二阶偏导数的有限差分近似42偏微分方程数值解5/9/20243.抛物型方程初条：精确解为（以热传导或磁扩散方程为例）初值问题不论初始分布如何集中，它总在瞬间影响于无穷远，虽该影响随距离按指数衰减，然而它是以无限速度传播。此乃抛物型方程解的特征。43偏微分方程数值解5/9/2024三、热传导方程（抛物方程）1.热传导方程的介绍2.离散化（1）向前差分格式：44偏微分方程数值解5/9/2024计算：这是一个显式格式（四点格式）每一层各个节点上的值是通过一个方程组求解得到的。这可以从下面的计算过程看出来。45偏微分方程数值解5/9/2024系数矩阵为46偏微分方程数值解5/9/2024计算实例：47偏微分方程数值解5/9/202448偏微分方程数值解5/9/20242.向后差分格式当知道第n层上的时，要确定第n+1层上各点值必须通过求解一个线性代数方程组。49偏微分方程数值解5/9/2024其矩阵表达式如下：50偏微分方程数值解5/9/2024这是一个古典四点向后差分格式。计算实例51偏微分方程数值解5/9/202452偏微分方程数值解5/9/20243.Crank-Nicolson格式，亦称六点对称格式53偏微分方程数值解5/9/202454偏微分方程数值解5/9/202455偏微分方程数值解5/9/20244.Richardson格式这是一个五点三层差分显式格式56偏微分方程数值解5/9/2024讨论：假若由于的作用，导致差分方程的近似解设为：于是，我们可得到差分格式的误差方程如下：xtRichardson格式是不稳定的。57偏微分方程数值解5/9/20245.稳定性判别Von-Neumann稳定性在判断有限差分近似的稳定性方法中，以Von-Neumann方法使用较为广泛，它仅适用于线性常系数的有限差分近似。其过程如下：首先，要研究的差分方程可写为：如，其次，对进行变量分离：58偏微分方程数值解5/9/2024最后将代入所考察的有限差分方程。定义为放大系数59偏微分方程数值解5/9/2024下面说明，在什么条件下能使对所有的成立。从上式，我们看出，60偏微分方程数值解5/9/2024６.交替显隐式格式（１）显式预测隐式校正格式在n+1/2层上用古典显式格式计算出过度值，再在n+1层上用古典隐格式校正预测值，即：61偏微分方程数值解5/9/2024(2).跳点格式首先将网格点（j,n）按j+n等于偶数或奇数分成两组，分别称为偶数网点和奇数网点。从到的计算过程中，先在偶数网格点上用古典显式格式计算，再在奇数网点上用古典隐格式计算，即：62偏微分方程数值解5/9/2024三双曲型方程（a）一阶常系数线性双曲型方程（b）二阶常系数线性双曲型方程（波动方程）

其中a为常数主要对象为这些方程的定解条件，可仅有初始条件，也可以有初始条件和边界条件。其中a为常数

同椭圆型方程与抛物型方程相比，双曲型方程差分格式的性质与定解问题解析解的性质有更密切的关系。63偏微分方程数值解5/9/20241一阶线性双曲型方程（1）初值问题考虑由于u(x,t)沿x-t平面上方向为dx/dt=a的直线x

at=C（C为常数）的变化率为0，即故沿x-t平面上任一条斜率为1/a的直线x

at=C，u(x,t)为常数。平行直线族x

at=C就是方程(3.1)的特征线。(3.2)(3.1)64偏微分方程数值解5/9/2024利用特征线，可以求出初值问题(3.1)、(3.2)的解：由于u(x,t)在点处的值依赖与

(x)在点的值，故初始线t=0上的点称为解u(x,t)在点的依赖区域。与抛物型方程求解类似，对x-t平面进行矩形网格剖分，x方向的步长为h，t方向的步长为

，网点简记为(j,k)。65偏微分方程数值解5/9/2024（1）偏心格式和中心差分格式对方程(3.1)，利用差商代替导数的方法，可得前两个格式的局部截断误差阶为，分别称为左、右偏心格式。第三个格式的截断误差阶为，它称为中心差分格式。其中即66偏微分方程数值解5/9/2024从差分格式依赖区域和微分方程依赖区域的关系，可以得到差分格式收敛的必要条件：

差分格式的依赖区域包含微分方程的依赖区域（也称为CFL条件）。对于左偏心格式，CFL条件为：a>0，且。对于右偏心格式，CFL条件为：a<0，且。所以，当a>0（或a<0）时，左（或右）偏心格式才有实用值；中心差分格式无实用价值。对于中心格式，CFL条件为：。可以证明：左、右偏心格式的稳定条件与其CFL条件相同，但中心差分格式是恒不稳定的，67偏微分方程数值解5/9/2024对方程(3.1)，可利用特征线构造格式。设系数a>0，上网点A(j

1,k)，B(j,k)，C(j+1,k)处的解值已经算出，从点P(j,k+1)作特征线，它与线段AB交于点D。（2）Lax格式PCBAD由u(p)=u(D)，有这样，得到Lax格式：当，Lax格式稳定，截断误差阶为。68偏微分方程数值解5/9/2024（3）Lax—Wendroff格式对方程(3.1)，利用特征线作二次插值，即可得到Lax—Wendroff格式：当，Lax—Wendroff格式稳定，它的截断误差阶为。69偏微分方程数值解5/9/2024应该注意：边值条件的给法与其它两类方程不同。如果a>0，方程特征线向右倾，只能在x变化区域的左边界上给出边界条件：

（2）初边值问题如果a<0，方程特征线向左倾，只能在x变化区域的右边界上给出边界条件，即使x的变化区域为0

x

d

，也只能给出边界条件：a<0XOY0

x

da>0XOY0

x<70偏微分方程数值解5/9/2024设常数a>0，考虑下面模型问题：前面建立的几个显格式，都适用于这个问题。下面建立隐格式。71偏微分方程数值解5/9/2024连同初始条件与边界条件：（1）最简隐格式该格式的局部截断误差阶为。令，格式可改写为该格式可在0<x,t<内所有网点上显示地计算解之近似值。72偏微分方程数值解5/9/2024然后用中心差商逼近这些导数值，则可得到Wendroff格式：在点处，用（2）Wendroff格式PCHADGFBE73偏微分方程数值解5/9/2024连同初始条件与边界条件：该格式的局部截断误差阶为，且无条件稳定。令，格式可改写为该格式可在0<x,t<内所有网点上显示地计算解之近似值。74偏微分方程数值解5/9/20242二阶线性双曲型方程（波动方程）

考察对x-t平面进行矩形网格剖分，x方向的步长为h，t方向的步长为

，网点简记为(j,k)。用二阶中心差商代替(3.3)中的二阶导数，则得到网点(j,k)处的差分方程：

(3.3)(3.4)其中。或

75偏微分方程数值解5/9/2024该格式稳定的充分条件为。初始条件离散：初始条件离散：由消去，得上述差分格式与初始条件的截断误差阶均为。取为76偏微分方程数值解5/9/2024上述方法也可用于求解初边值问题：3交替方向隐格式77偏微分方程数值解5/9/2024任何模拟方法，都必须在最佳计算速度和数值精度之间寻找平衡点。要在各种可能的求解方法中找到一种统一地适用于计算材料学领域（或其它领域）的理想方法，一般是不现实的。

由于实际问题的具体特征、复杂性以及算法自身的适用范围决定了应用中必须选择、设计适合于自己特定问题的算法，因而掌握数值方法的思想至关重要。

科学计算（数值模拟）已经被公认为与理论分析、实验分析并列的科学研究三大基本手段之一，但三者之间相辅相成。78偏微分方程数值解5/9/2024谢谢大家！79偏微分方程数值解5/9/202480偏微分方程数值解5/9/2024第三部分偏微分方程的有限元方法一边值问题的变分原理1引论模型：在条件下求使得泛函达到最大的函数。（1）等周问题在长度一定的所有平面封闭曲线中，求所围面积为最大的曲线。81偏微分方程数值解5/9/2024定义：当求泛函在一个函数集合K中的极小（或极大）问题，则该问题称为变分问题。变分问题与微分方程的定解问题有一定的联系。（2）初等变分原理①一元二次函数的变分原理考察J(x)的极值情况。变分原理：设求，使与求解方程Lx=f等价。82偏微分方程数值解5/9/2024对称正定②多元二次函数的变分原理求J(x)取极小值的驻点，其中设设则J(x)可表示为：83偏微分方程数值解5/9/2024变分原理：设矩阵A对称正定，则下列两个命题等价：求，使(a)(b)是方程的解上述两个例子表明：其中

求二次函数的极小值问题和求线性代数方程（组）的解是等价的。84偏微分方程数值解5/9/2024（1）弦平衡的平衡原理与极小位能原理2两点边值问题的变分原理考察一根长为l的弦，两端固定在点A(0,0)和B(l,0)。当没有外力作用时，它的位置沿水平方向与X轴重合。设有强度为f(x)的外荷载垂直向下作用在弦上，于是弦发生形变。假定荷载很小，因而发生的形变也很小。用u(x)表示在荷载f(x)的作用下弦的平衡位置。

85偏微分方程数值解5/9/2024求弦的平衡位置归结为求解两点边值问题：设弦处于某一位置u=u(x)，可得到其总位能为

极小位能原理：其中T是弦的张力。平衡原理弦的平衡位置(记为)将在满足边值条件u(0)=0，u(l)=0的一切可能位置中，使位能取极小值。弦的平衡位置是下列变分问题的解86偏微分方程数值解5/9/2024在数学上，要将某个微分方程的定解问题转化为一个变分问题求解，必须针对已给的定解问题构造一个相应的泛函，并证明定解问题的解与泛函极值问题的解等价。有限元方法正是利用这种等价性（边值问题与变分问题的等价性），先将微分方程定解问题转化为变分问题（或变分方程）的求解问题，然后再设法近似求解变分问题（或变分方程）。87偏微分方程数值解5/9/2024（2）两点边值问题的变分原理①构造泛函考察二阶常微分方程边值问题：引入泛函算子则88偏微分方程数值解5/9/2024②变分问题与前述二阶常微分方程边值问题相应的变分问题是其中求，使89偏微分方程数值解5/9/2024③变分原理（变分问题与边值问题的等价性）设，是边值问题的解，则使J(u)达到极小值；反之，若使J(u)

达到极小值，则是边值问题的解。其中是强制边界条件，是自然边界条件，区别这两类边界条件在用有限元方法求解边值问题时很重要。90偏微分方程数值解5/9/2024（3）虚功原理对两点边值问题：其中虚功原理，且满足变分方程：设，以v乘方程两端，沿[a,b]积分，并利用，得变分方程对任意在力学里，表示虚功设，则

是边值问题解的充要条件是：91偏微分方程数值解5/9/2024

对于复杂的边界条件，边值问题的求解一般是困难的。若将微分方程化为相应的变分问题或变分方程，则只需处理强加边界条件，无需处理自然边界条件（自然边界条件已包含于变分问题中泛函的构造或已包含于给出的变分方程之中）。这一特点对研究微分方程离散化方法及其数值解带来了极大的方便。92偏微分方程数值解5/9/20243二阶椭圆边值问题的变分原理（1）极小位能原理模型方程其中G是平面有界区域。①构造泛函引入泛函算子则93偏微分方程数值解5/9/2024②变分问题与前述二阶椭圆边值问题相应的变分问题是求，使其中94偏微分方程数值解5/9/2024③变分原理（变分问题与边值问题的等价性）对第一边值问题，无论齐次或非齐次边界条件，泛函是一样的，只是边界条件要作为强加边值条件加在所取的函数类上。设，是二阶椭圆边值问题的解，则使J(u)达到极小值；反之，若使J(u)

达到极小值，则是二阶椭圆边值问题的解。其中对第二、三类边值问题，无论齐次或非齐次边界条件，二次泛函形式相对于第一边值问题有所改变，但函数类的选取与边界条件无关。95偏微分方程数值解5/9/2024（2）虚功原理问题其中设，以v乘方程两端后在G上积分，并利用Green公式，得变分方程96偏微分方程数值解5/9/2024虚功原理在力学里，表示虚功设是边值问题的解，则对任意，满足变分方程。反之，若，且对任意满足变分方程，则为边值问题的解。与极小位能原理类似，第一类边界条件为强加边界条件，第二、三类边界条件为自然边界条件。虚功原理比极小位能原理应用更广。97偏微分方程数值解5/9/2024目的：求解相应的变分问题或相应的变分方程。

Ritz方法是近似求解变分问题(即二次泛函极小值)的算法。Galerkin方法是近似求解变分方程的算法，这两种算法统称为Ritz-Galerkin方法。Ritz-Galerkin方法的基本思想以下用V表示等Sobolev空间，L表示微分算子，

(u,v)为由L及边值条件决定的双线性泛函。4Ritz-Galerkin方法用有限维空间的函数代替变分问题(或变分方程)中无限维空间的函数，从而在有限维函数空间中求变分问题(或变分方程)的近似解，并要求当有限维空间的维数不断增加时，有限维近似解逼近原变分问题(或变分方程)的解。98偏微分方程数值解5/9/2024由极小位能原理得出的变分问题为：Ritz方法：求变分问题的近似解。（1）Ritz方法求，使其中，设是V的n维子空间，是的一组基底(称为基函数)。中任一元素可表示为即选择适当的，使取极小值。求，使Ritz方法：99偏微分方程数值解5/9/2024展开令则满足解出代入，则得100偏微分方程数值解5/9/2024Ritz方法步骤为：根据最小位能原理构造相应于微分方程或物理问题的变分问题；取作为的一组基底，即用近似代替无穷维空间V；根据二次函数取极值的必要条件，得到中所满足的方程组：求解关于的线性代数方程组。101偏微分方程数值解5/9/2024由虚功原理得出的变分方程为：Galerkin方法：求变分方程的近似解。（2）Galerkin方法设是V的n维子空间，是的一组基底(称为基函数)。中任一元素可表示为即选择适当的，使取极小值。Galerkin方法：求，使对，满足102偏微分方程数值解5/9/2024由的任意性，取作为v，则得将代入变分方程，则解出代入，则得103偏微分方程数值解5/9/2024Galerkin步骤为：根据虚功原理构造相应于微分方程或物理问题的变分方程；取作为的一组基底，即用近似代替无穷维空间V；求解关于的线性代数方程组。取作为v，将代入变分方程，得到满足的方程组：104偏微分方程数值解5/9/2024有限元法广泛应用的原因Ritz-Galerkin方法应用的困难①基函数选取必须满足强加边界条件，因此选取困难；②计算量、存储量巨大；③方程组求解病态严重。充分发挥了变分形式和Ritz-Galerkin方法的优点；②摆脱了传统的基函数取法；③各种问题的结构程序格式统一。105偏微分方程数值解5/9/2024有限元方法基于变分原理，又具有差分方法的一些特点，并且适于较复杂的区域和不同粗细的网格。二椭圆型方程的有限元方法

差分法解偏微分方程，解得的结果就是准确解u在节点上的近似值；Ritz-Galerkin方法得到近似的解析解，但对一般区域，却往往难以实现。有限元方法与传统Ritz-Galerkin方法的差别在于有限维函数空间的构造方法。Ritz-Galerkin方法选用的基函数在整个定解区域上整体光滑，有限元则取分段或分片连续且局部非零的基函数。106偏微分方程数值解5/9/2024考虑两点边值问题：1一维问题的线性元将区间[a,b]分割为n个子区间。第i个单元记为，其长度。（1）试探函数与试探函数空间设则称为试探函数空间，称为试探函数。107偏微分方程数值解5/9/2024（2）用单元形状函数表示试探函数设在节点上试探函数在节点上的一组值为最简单的试探函数空间由分段线性函数组成。在第i个单元上的线性插值函数为即当时，的（线性）插值公式称为（线性）单元形状函数。108偏微分方程数值解5/9/2024把每个单元形状函数合并起来，就得到整个区间[a,b]上都有定义的函数：109偏微分方程数值解5/9/2024为使分段插值标准化，通常用仿射变换显然把变到，令则变为或110偏微分方程数值解5/9/2024定义基函数系（3）用节点基函数表示试探函数111偏微分方程数值解5/9/2024线性无关，它们可组成试探函数空间的基，常称为节点基函数。几何形状如图ab任一试探函数可表示为用这类插值型基函数，可以构造出适合各种边界条件的试探函数。112偏微分方程数值解5/9/2024若借助前述放射变换节点基函数可用变量

表示为113偏微分方程数值解5/9/2024①直接形成有限元方程(a)把表达式代入泛函；（4）从Ritz方法出发形成有限元方程(b)将泛函表达式中积分区间[a,b]变到[0,1]；(c)由达到极小值的条件得到含的有限元方程这儿(d)解出有限元方程的数值解，就得到使二次泛函取极小的近似函数(有限元解)114偏微分方程数值解5/9/2024有限元方程可用矩阵表示为其中称为总刚矩阵。115偏微分方程数值解5/9/2024工程中形成有限元方程时，通常先在每个单元上形成单元矩阵（称为单元刚度矩阵），然后由单元刚度矩阵形成总刚度矩阵（称为总体合成）。②用单元刚度分析形成有限元方程(a)把按单元组织，则在第i个单元上，令其中称为单元刚度矩阵。各元素可计算得到。116偏微分方程数值解5/9/2024再把扩展成n

n矩阵，使其第i

1行、第i行和第i

1列、第i列交叉位置的元素就是单元刚度矩阵的四个元素，其余全为零（只是第一行，第一列元素非零）。即记则其中称为总刚矩阵。117偏微分方程数值解5/9/2024(b)由达到极小值的条件(c)解出有限元方程的数值解，就得到使二次泛函取极小的近似函数(有限元解)得到有限元方程。118偏微分方程数值解5/9/2024（5）从Galerkin方法出发形成有限元方程把表达式代入变分方程对前面的两点边值问题，变分方程变为其中与Ritz方法相比，Galerkin方法形成的有限元方程其系数矩阵就是总刚矩阵。该方程即为Galerkin法形成的有限元方程。由Galerkin方法推导有限元方程更加方便直接，且适用面广。119偏微分方程数值解5/9/2024

若希望在每个单元上提高逼近的精确度，则可通过提高插值多项式次数来实现，在单元上可构造一、二、三及高次插值多项式，其方法有两种：2一维问题的高次元

整个问题计算的全过程除分析单元插值外，均与前面框架类似。①Lagrange型：在单元内部增加一些插值节点。②Hermite型：在节点引进一阶、二阶乃至更高阶导数。120偏微分方程数值解5/9/2024①线性元(Lagrange型)要求：在每一个单元上是一次多项式，在单元节点处连续。插值条件：在单元的两个端点取指定值。②二次元(Lagrange型)要求：在每一个单元上是二次多项式，在单元节点处连续。插值条件：在单元的两个端点及单元中点取指定值。③三次元（Hermite型）要求：在每一个单元上是三次多项式，在单元节点处连续。插值条件：在两个端点取指定的函数值和一阶导数值。121偏微分方程数值解5/9/2024采用高次元，有限元方程形成的方法和线性元类似，但工作量增加。一是计算积分的复杂性增加，二是矩阵的带宽增加。高次元的主要优点是收敛阶高，且提高了函数逼近的光滑性。122偏微分方程数值解5/9/2024假定区域G可以分割成有限个矩形的和，且每个小矩形（单元）的边和坐标轴平行。3二维问题的矩形元通过仿射变换采用矩形剖分后，任一个矩形总可变成单位正方形如果在上造出单元形状函数，就可得到试探函数。而上的形状函数可通过先在上造出形状函数，再通过仿射变化而得到。123偏微分方程数值解5/9/2024在上构造形状函数，也采用Lagrange型和Hermite型插值。Lagrange型：根据若干插值节点处的函数值决定插值函数。Hermite型：根据若干插值节点处的函数值、一阶偏导数乃至更高阶偏导数决定插值函数。124偏微分方程数值解5/9/2024（1）Lagrange型公式①双一次插值插值条件：给定顶点上的函数值求：双线性函数满足设令由为双线性函数，可求得125偏微分方程数值解5/9/2024令则通过仿射变换消去

、

，就得到上的形状函数。把这些函数按单元叠加，即对所有单元求和，就得到G上的试探函数。实际计算时，并不消去中间变量

、

，因为计算刚度矩阵元素（定积分）用

、

作自变量更为方便。126偏微分方程数值解5/9/2024插值条件：给定I

I上九个插值节点(0,0)、(1/2,0)、(1,0)、(0,1/2)、(1/2,1/2)、(1,1/2)、(0,1)、(1/2,1)、(1,1)的函数值。求：双二次函数满足②双二次插值127偏微分方程数值解5/9/2024故通过仿射变换消去

、

，就得到上的形状函数。令由为二次函数，可求得设128偏微分方程数值解5/9/2024插值条件：给定I

I上十六个插值节点(见图)。求：双三次函数满足设③双三次插值129偏微分方程数值解5/9/2024故令由为三次函数，可求得130偏微分方程数值解5/9/2024可以在四个顶点分别给定函数值、两个一阶偏导数的值和二阶混合偏导数的值(共十六条件)，确定一个双三次多项式的十六个系数。（2）Hermite型公式Lagrange型公式中不出现导数，这样的试探函数只属于。为了得到属于的试探函数，需要Hermite型插值公式。双三次多项式含有十六项：简单且常用的是不完全的双三次多项式插值。它去掉双三次多项式中的项。131偏微分方程数值解5/9/2024插值条件：给定I

I上四个插值节点。求：不完全双三次函数满足四个顶点处的函数值等于在该点的函数值；四个顶点处的值等于在该点的值；四个顶点处的值等于在该点的值。根据仿射变换则可将原插值问题转化为I

I上的插值问题。132偏微分方程数值解5/9/2024满足四个顶点处的函数值等于在该点的函数值；四个顶点处的值等于在该点的值乘以

x；四个顶点处的值等于在该点的值乘以

y。插值条件：给定I

I上四个插值节点(0,0)、(1,0)、(0,1)、(1,1)。求：不完全双三次函数类似于Lagrange型公式的构造，可以求得上的形状函数。133偏微分方程数值解5/9/2024

在三角形元的有限元方法中，先将定解区域G化分为若干个小三角形（称作单元）。然后在每个单元上构造插值型函数，并用分片函数（但整体连续的函数）代替变分问题或变分方程中所需求解的函数。4二维问题的三角形元用有限元求解二维椭圆边值问题时，应用最广的是三角形元。134偏微分方程数值解5/9/2024（1）三角剖分将定解区域化分成若干个小三角形单元时应注意：③为了保证有限元解的精确度和收敛性，并避免其离散后代数方程组系数矩阵的病态性，网格剖分中疏密的过渡不要太陡。错误为了保证有限元解有