1、偏微分方程数值解偏微分方程数值解(Numerical Solution of Partial Differential Equations)主讲：王曰朋主讲：王曰朋1精选课件参考数目1.George J.Haltiner,Roger Terry Williams,Numerical Prediction and Dynamic 2.Meteorology(2nd Edition),the United States of America,1979.2.Curtis F.Gerald and Patrick O.,Applied Numerical Analysis,Person Educati
2、on,Inc.,2004.3.Eugenia Kalnay,Atmospheric Modeling,Data Assimilation and Predictability,the press Syndicate of the University of Cambridge,2003.4.Arieh Iserles,A First Course in the Numerical Analysis of Differential Equations,Cambridge University Press,1996.5.李荣华，冯国忱.微分方程数值解.北京：人民教育出版社，1980.6.徐长发，李
3、红.实用偏微分方程数值解法.华中科技大学出版社，2003.7.沈桐立，田永祥等.数值天气预报.北京：气象出版社，2007.2精选课件数值天气预报PDE数值解1.挪威气象学家V.Bjerknes（1904）提出数值预报的思想：通过求解一组方程的初值问题可以预报将来某个时刻的天气思想；2.L.F.Richardson(1922)：开创了利用数值积分进行预报天气的先例，由于一些原因（如，计算稳定性问题“Courant,1928”）并没有取得预期的效果尝试；3.Charney,Fjortoft,and Von Neumann(1950),借助于Princeton大学的的计算机(ENIAC)，利用一个简
4、单的正压涡度方程（C.G.Rossby,1940）对500mb的天气形式作了24小时预报-成功；3精选课件The Electronic Numerical Integrator and Computer(ENIAC).4精选课件常微分方程的数值解大气科学中常微分方程和偏微分方程的关系 1.大气行星边界层(近地面具有湍流运动特性的大气薄层,11.5km),埃克曼（V.W.Ekman）（瑞典）螺线的导出；2.1963年，美国气象学家Lorenz在研究热对流的不稳定问题时，使用高截断的谱方法，由Boussinesq流体的闭合方程组得到了一个完全确定的三阶常微分方程组，即著名的Lorenz系统。5精选
5、课件Lorenz系统dx/dt=a(y-x)dy/dt=x(b-z)-y dz/dt=xy-c z 其中，a=10,(Prandtl number);b=28(Rayleigh number);c=8/3;(x,y,z)_0=(0.01;0.01;1e-10)6精选课件05101520253035404550-30-20-10010203040507精选课件-20020-30-20-100102030010203040508精选课件9精选课件10精选课件Franceshini 将Navier-Stokes方程截断为五维的截谱模型如下：112345221 33312441655142449357
6、Re53xxx xx xxxx xxxx xxxx xxxx x=-+=-+=-+=-=-11精选课件欧拉法折线法 常微分方程能直接进行积分的是少数，而多数是借助于计算机来求常微分方程的近似解；有限差分法是常微分方程中数值解法中通 常有效的方法；建立差分算法的两个基本的步骤：1.建立差分格式，包括：a.对解的存在域剖分；b.采用不同的算法可得到不同的逼近误差截断误差（相容性）；c.数值解对真解的精度整体截断误差（收敛性）；d.数值解收敛于真解的速度；e.差分算法舍人误差（稳定性）.12精选课件2.差分格式求解将积分方程通过差分方程转化为代数方程求解，一般常用递推算法。在常微分方程差分法中最简单
7、的方法是Euler方法，尽管在计算中不会使用，但从中可领悟到建立差分格式的技术路线，下面将对其作详细介绍:13精选课件差分方法的基本思想“就是以差商代替微商”23()1111()()()()()()()2!3!nnnu thu tu t hu t hut hut hO hn23()1111()()()()()()()2!3!nnnu thu tu t hu t hut hut hO hn考虑如下两个Taylor公式：（1）（2）从（1）得到：1()()()()iiiu tu tu tO hh14精选课件从（2）得到：1()()()()iiiu tu tu tO hh211()()()()2ii
8、iu tu tu tO hh从（1）-（2）得到：从（1）+（2）得到：2112()2()()()()iiiiu tu tu tu tO hh15精选课件对经典的初值问题0(,)(0)duf t udtuu(0,)tT1212(,)(,)f t uf t uL uu满足Lipschitz条件保证了方程组的初值问题有唯一解。16精选课件一、算法构造：0tuTnt0tit1iitth1.在求解域上等距离分割：2.在 有：1,iit t1()()(,)()iiu tu tf t uO hh1(,)iiiiuuf t uh1(,)iiiiuuhf t u微分方程的精确解差分方程的精确解17精选课件3.
9、应用时采用如下递推方式计算：0u1u2u3u0t1t2t4.例题对初值问题2(0)1yxyy(0,1)x5,N 0.2h 用Euler法求解，用即，001111223334445(,)1.000,1.200;(,)1.600,1.520;(,)2.320,1.984;(,)3.184,2.621;(,)4.221,3.465f xyyf x yyf xyyf xyyf xyy18精选课件5.Euler法的几何意义0t0tit1iitth在递推的每一步，设定()iiu tu1()iu t1iu2()O h()iu t(),iiu tu过点(,)iit u作 的切线，该切线的()u t方程为：11
10、()()iiiiiuuu ttt即：1(,)iiiiuuhf t u19精选课件二、误差分析构造算法后，这一算法在实际中是否可行呢？也就是说是否使计算机仿真而不失真，这还需要进一步分析。1.局部截断误差-相容性为了分析分析数值方法的精确度，常常在()iiuu t成立的假定下，估计误差111()iiieu tu这种误差称为“局部截断误差”，如图。21()()()()2iiihu tu thu tu2()(,()()2iiihu thf t u tu22111()()()2iiiheu tuuO h局部截断误差是以点it的精确解()iu t为出发值，用数值方法推进到下一个点1it而产生的误差。20
11、精选课件整体截断误差是以点0t的初始值0u为出发值，用数值方法推进i+1步到点1it，所得的近似值 与精确值 的偏差：2.整体截断误差收敛性1iu1()iu t111()iiiu tu称为整体截断误差。1(,)iiiiuuhf t u11()()(,()iiiiiu tu thf t u te21(,()(,)iiiiiih f t u tf t uDh21iiihLDh12210(1)1(1)(1)(1)iiihLDhhLhLhL2110(1)(1)1iiDhhLhLhL21精选课件0(1)TLTLDheeL特例，00若不计初始误差，即则1(1)TLiDheL1()iO h即3.舍入误差稳定
12、性假设一个计算机仅表示4个数字（小数点后面），那么10.3333,310.11119计算20.373410.1112011119110.33340.33333xxx20.33341190.6667133xxxx22精选课件我们的要求是：最初产生的小误差在以后的计算中虽然会传递下去，但不会无限制的扩大，这就是稳定性所描述的问题。下面引进稳定性的概念：设由初值0u0v得到精确解 ，nu由初值得到精确解 ，nv若存在常数C和充分小的步长0h使得00nnuvC uv则称数值方法是稳定的。tu00u0vntnunv23精选课件2(0)1dyxydxyy计算例题(0,1)x其解析解为：12yxx=0 0.
13、2000 0.4000 0.6000 0.8000 1.0000y=1.0000 1.2000 1.3733 1.5315 1.6811 1.826924精选课件00.10.20.30.40.50.60.70.80.9111.11.21.31.41.51.61.71.81.9xy25精选课件三、改进的Euler法将微分方程()(,()u tf t ut在区间1,iit t上积分，得到11()()(,()iitiitu tu tf t u t dt用梯形法计算积分的近似值，有1111(,()(,()(,()2iitiiiitf t u t dtf t u tf tu t于是1111()()(,(
14、)(,()2iiiiiiu tu th f t u tf tu t1111(,)(,)2iiiiiiuuh f t uf tu这是一个隐式格式，一般需要用迭代法来求 ，而用显式的Euler法提供初值。1iu26精选课件01(,)iiiiuuhf t u1 111/2(,)(,)kkiiiiiiuuhf t uf tu为了简化计算的过程，在此基础上进一步变为如下算法：111/2(,)(,)iiiiiiuuhf t uf tu1(,)iiiiuuhf t u此式称为“改进的Euler法。接下来讨论其几何意义预估校正其局部截断误差为3()O h这个问题将在下节讨论。27精选课件tu0it1it0(,
15、)iimf t u()u tiu1iu111(,)iimf tu122averagemmm1iu1()iu t28精选课件00.10.20.30.40.50.60.70.80.9111.11.21.31.41.51.61.71.81.9xyEuler法、改进的Euler法和解析解的比较2(0)1dyxydxyy(0,1)x12yx29精选课件四、（龙格-库塔）Runge-Kutta方法简单的Euler法是建立在Taylor级数的一项展开；改进的Euler法是以两项Taylor级数为基础建立的，如：23()()()()()2iiiihu thu thu tu tO h23()()()()()/2
16、()iiiiiu thu thu thu thu thO h()()()2iiiu thu tu th11(,)(,)2iiiiif tuf t uuh如果我们截取Taylor级数的更多项会得到什么样的求解方法呢？两个德国数学家（C.Runge&M.kutta）以这种思想为基础建立了求解微分方程的龙格-库塔方法。它是常微分方程数值解法中使用最为广泛的方法之一。30精选课件一般地，一个K阶的Runge-Kutta方法可用下面的公式表示：11111(,)(,)kiijjjiijjijijllluuw KKhf t uKhf tc h ua K2,3,jk其中，jw是待定的加权系数，2321,1,k
17、k kc cc aa是待定的系数。Euler法就是11,1kw的R-K法。其系数的确定如下：将1iu展开成h的幂级数，并与微分方程的精确解1()iu t在点it的Taylor展开式相比较，使两者的前1p项相同，这样确定的R-K法，其局部截断误差为1()pO h，根据所得关于待定系数的方程组，求出它们的值后代入公式，就成为一个p阶R-K方法。31精选课件例题以二阶R-K法为例说明上述过程11122112211(,)(,)iiiiiiuuw Kw KKhf t uKhf tc h ua K21()()2iiiihuuhu tu t2(,)(,)2iiiiihuhf t uf t u211(,)(,
18、)(,)(,)22iiitiiuiiiiuhf t uhf t uf t uf t u把12,K K代入1iu中，有112221(,),(,)iiiiiiiiuuwhf t uw hf tc h ua hf t u12221(,)(,)(,)(,)iiiiitiiuiiuwhf t uw h f t uc hf t ua hff t u21222221()(,)(,)(,)iiitiiuiiuh wwf t uh w c f t uw a ff t u32精选课件经比较得到122222111212www cw a1222212112www cac取 为自由参数：2c21 2,12 3c 121
19、 31 1(,)(0,1),(,),(,)4 42 2w w从而得到不同的但都是二阶的R-K方法，对应的有中点法、Heun（亨）法以及改进的Euler法。33精选课件基于相同的过程，通过比较五次Taylor多项式，得到更加复杂的结果，给出了包含13个未知数的11个方程。得到多组系数，其中常用的是以下四阶R-K法：1123412132431()6(,)11(,)2211(,)22(,)iiiiiiiiiiuuKKKKKhf t uKhf th uKKhf th uKKhf th uK改进的Euler法、R-K法以及解析解的比较：34精选课件00.10.20.30.40.50.60.70.80.9
20、111.11.21.31.41.51.61.71.8xy2(0)1dyxydxyy(0,1)x12yx35精选课件五、线性多步（Linear Multistep Method）法1.预备知识：插值多项式插值是离散函数逼近的重要方法，利用它可通过函数在有限个点处的取值状况，估算出函数在其他点处的近似值。从几何上理解：对一维而言，已知平面上n1个不同点，要寻找一条n次多项式曲线通过这些点。插值多项式一般常见的是拉格朗日插值多项式。2.气象应用不均匀站点上的气象要素数据均匀网格点上的数据插值 3.拉格朗日插值多项式拉格朗日插值多项式逼近可能是求插值节点不均匀的插值多项式的最简单的方法。实验观察结果或
21、原始测量数据的分布通常是非均匀的。例如，四个点可以确定一个三次多项式，其拉格朗日形式为：23413431212131421232412312434313234414243()()()()()()()()()()()()()()()()()()()()()()()()()xxxxxxxxxxxxP xffxxxxxxxxxxxxxxxxxxxxxxxxffxxxxxxxxxxxx36精选课件4.Adams-Bashforth（阿达姆斯贝雪福斯）公式首先，用以下四个点对(,()f t u t进行三次Langrage插值：111222333(,(,(),(,(,(),(,(,(),(,(,()nnn
22、nnnnnnnnntf t u ttf tu ttf tu ttf tu t则1233123123333132()()()(,)()(,()()()()()()()(,()()()()nnnnnnnnnnnnnnnnnnnnnnttttttf x ttf t u tttttttttttttf tu ttttttt于是，有1113()()(,()()()nnnntnnttntu tu tf t u t dtu tt dt容易算出，13123()55()59()37()9()24nnxnnnnxhx dxF xF xF xF x例如，我们可以算得11233123()()()1()(2)(3)()(
23、)()6nnnnxxhnnnnnnxxnnnnnnxxxxxxdxxxh xxh xxh dxxxxxxxh*2*137精选课件433011 55()(2)(3)5566424hhth th th dthhh将（*2）代入（*1）得到Adams-Bashforth公式：111223355(,)59(,)37(,)9(,)24nnnnnnnnnnhuuf x uf xuf xuf xu基于同样的计算过程，可以得到另外一个计算公式：11111229(,)19(,)5(,)(,)24nnnnnnnnnnhuuf xuf x uf xuf xu这称为Adams-Moulto公式。111223355(,
24、)59(,)37(,)9(,)24nnnnnnnnnnhuuf x uf xuf xuf xu11111229(,)19(,)5(,)(,)24nnnnnnnnnnhuuf xuf x uf xuf xu预估校正38精选课件偏微分方程数值解主讲：王曰朋39精选课件一、区域的离散1.2.3.40精选课件则函数可表示为：(,)(,)u x yu ih jk1,2,3i 1,2,3j 二、1.（一维）一、二阶导数的有限差分近似表达式41精选课件2.（二维）一、二阶偏导数的有限差分近似(,)(,)(,)()ijijijx yu xh yu x yuO hxh1(,)ijjjiix yuuuxh2(,)
25、(,)(,)()2ijijijx yu xh yu xh yuO hxh11(,)2ijjjiix yuuuxh211111111(,)1()22ijjjjjiiiix yuuuuux yhk 22()()O hO k42精选课件3.抛物型方程抛物型方程初条：精确解为（以热传导或磁扩散方程为例）（以热传导或磁扩散方程为例）22xutu)()0,(xfxu初值问题)0,(txd)(4)(exp41),(2ftxttxu不论初始分布如何集中，它总在瞬间影响于无穷远，虽该影响随距离按指数衰减，然而它是以无限速度传播。此乃抛物型方程解的特征。43精选课件三、热传导方程（抛物方程）1.热传导方程的介绍2
26、.222(0,)(,)0(,0)()uuatxutu L tu xf x离散化121122jjjjjiiiiiuuuuuakh0(0,)0juujk(,)0jNuu L t0(,0)()iiuu ihf ihf（1）向前差分格式：44精选课件计算：111(1 2)jjjjiiiiusus usu这是一个显式格式（四点格式）j1j i1i1i每一层各个节点上的值是通过一个方程组求解得到的。这可以从下面的计算过程看出来。22kash45精选课件1000121010002321100034321000112(1 2)(1 2)(1 2)(1 2)NNNNusus usuusus usuusus us
27、uusus usu1 21 21 21 2ssssssss系数矩阵为46精选课件22(,0)sin()(0,)(1,)01;00.1uutxu xxututxt 计算实例：47精选课件00.020.040.060.080.100.51-3-2-10123txu48精选课件2.向后差分格式111121122jjjjjiiiiiuuuuuakh11111(12)jjjjiiiisus usuu当知道第n层上的 时，要确定第n+1层上各点值 必须通过求解一个线性代数方程组。jiu1jiu11121011113212111432311111(12)(12)(12)(12)jjjjjjjjjjjjjjj
28、jNNNNsus usuusus usuusus usuusus usuu49精选课件111122111112121212jjjjjjNNjjNNuussuusssssuussuu其矩阵表达式如下：j1j 1ii1i50精选课件这是一个古典四点向后差分格式。计算实例22(,0)sin()(0,)(1,)01;00.1uutxu xxututxt 51精选课件00.020.040.060.080.100.5100.20.40.60.81txu52精选课件3.Crank-Nicolson格式，亦称六点对称格式1211111112222()2jjjjjjjjiiiiiiiiuuauuuuuukhh1
29、111111(12)(1 2)jjjjjjiiiiiisus ususus usu11121111/2/21/21/212jjjNjNussusssssussu53精选课件j1j 1i1i1211/2/21/21/2/21 2jjjNjNussusssssussu54精选课件00.020.040.060.080.100.5100.20.40.60.81txu55精选课件4.Richardson格式11211222jjjjjiiiiiuuuuuakh11112(2)jjjjjiiiiius uuuu1j 1i1ij1j i这是一个五点三层差分显式格式56精选课件讨论：假若由于 的作用，导致差分方
30、程的近似解设为：jiu 于是，我们可得到差分格式的误差方程如下：11112(2)jjjjjiiiiisxtRichardson格式是不稳定的。57精选课件5.稳定性判别 Von-Neumann 稳定性在判断有限差分近似的稳定性方法中，以Von-Neumann方法使用较为广泛，它仅适用于线性常系数的有限差分近似。其过程如下：首先，要研究的差分方程可写为：101nnmj mmj mm Nm Na ub u如，其次，对进行变量分离：jiu2expnnjpjppuVixl111(1 2)nnnnjjjjusus usu58精选课件expnnjjuVix最后将代入所考察的有限差分方程。111expexp
31、(1 2)expexpnnjjnnjjVi xsVi xs Vi xsVi x1exp(1 2)expnnnnVsVi hs VsVi h1(cossin)(1 2)(cossin)1 2(1 cos)nnVshihsVshihsh 定义为放大系数59精选课件21 4 sin2hs 下面说明，在什么条件下能使11nnVV对所有的成立。21 1 4 sin12hs 从上式，我们看出，21 1 4 sin2hs 212 sin2sh60精选课件12s.交替显隐式格式（）显式预测隐式校正格式在n+1/2层上用古典显式格式计算出过度值，再在n+1层上用12nju古典隐格式校正预测值，即：1/21121
32、1/21111122222nnnnnjjjjjnnnnnjjjjjuuuuuhuuuuuh61精选课件(2).跳点格式首先将网格点（j,n）按j+n等于偶数或奇数分成两组，分别称为偶数网点和奇数网点。从 到的计算过程中，先在偶数网格点上用古典显式格式计算，再在奇数网点上用古典隐格式计算，即：nt1nt111211111122121nnnnnjjjjjnnnnnjjjjjuuuuuhnjuuuuuhnj偶 数奇 数62精选课件222220uuatx0uuatx（a）一阶常系数线性双曲型方程（b）二阶常系数线性双曲型方程（波动方程）其中a为常数为为 这些方程的定解条件，可仅有初始条件，也可以有初始
33、条件和边界条件。其中a为常数 同椭圆型方程与抛物型方程相比，双曲型方程差分格式的性质与定解问题解析解的性质有更密切的关系。63精选课件考虑0dtdxatudtdxxutudtduxxuatu0 由于u(x,t)沿x-t平面上方向为dx/dt=a的直线 xat=C（C为常数）的变化率为0，即故沿x-t平面上任一条斜率为1/a的直线xat=C，u(x,t)为常数。平行直线族xat=C就是方程(3.1)的特征线。(3.2)(3.1)xxxu)()0,(64精选课件 利用特征线，可以求出初值问题(3.1)、(3.2)的解：)(),(atxtxu),(00tx 由于u(x,t)在点 处的值依赖与(x)在
34、点 的值，故初始线t=0上的点 称为解u(x,t)在点 的依赖区域。00atx 00atx),(00tx 与抛物型方程求解类似，对x-t平面进行矩形网格剖分，x方向的步长为h，t方向的步长为，网点 简记为(j,k)。),(kjtx65精选课件对方程(3.1)，利用差商代替导数的方法，可得011huuauukjkjkjkj011huuauukjkjkjkj02111huuauukjkjkjkj,2,1,0,2,1,0kj 前两个格式的局部截断误差阶为 ，分别称为左、右偏心格式。)(hO)(2hO 第三个格式的截断误差阶为 ，它称为中心差分格式。kjkjkjurarauu)1(11kjkjkjra
35、uurau11)1()(2111kjkjkjkjuuaruu其中hr/即66精选课件 从差分格式依赖区域和微分方程依赖区域的关系，可以得到差分格式收敛的必要条件：对于左偏心格式，CFL条件为：a0，且 。1ar1ar对于右偏心格式，CFL条件为：a0（或a0，上网点A(j1,k)，B(j,k)，C(j+1,k)处的解值已经算出，从点P(j,k+1)作特征线，它与线段AB交于点 D。ktt P1kttktt CBAD由u(p)=u(D)，有)(2)(2)(AuhahCuhahDu,2,1,0,2,1,0kj)1()1(21111kjkjkjuaruaru这样，得到Lax格式：当 ，Lax格式稳定
36、，截断误差阶为 。1ar)(22hhO68精选课件 对方程(3.1)，利用特征线作二次插值，即可得到LaxWendroff格式：当 ，LaxWendroff格式稳定，它的截断误差阶为 。1ar)(22hOkjxkjkjkjkjurauuaruu2221112)(269精选课件应该注意：边值条件的给法与其它两类方程不同。如果 a0，方程特征线向右倾，只能在 x变化区域的左边界上给出边界条件：)(),0(1ttu 如果 a0，方程特征线向左倾，只能在 x变化区域的右边界上给出边界条件，即使 x 的变化区域为0 x d，也只能给出边界条件：)(),(2ttdua0XOY0 x0，考虑下面模型问题：前
37、面建立的几个显格式，都适用于这个问题。TtttuxxxuxTtxuatu0)(),0(0)()0,(000下面建立隐格式。71精选课件连同初始条件与边界条件：,2,1,2,1,001111jkhuuauukjkjkjkjkjkjkjuaruararu111111该格式的局部截断误差阶为 。)(hO令 ，格式可改写为hr/该格式可在0 x,t 内所有网点上显示地计算解之近似值。,2,1,0,2,1,00kjuukkjj),(kj)1,(kj)1,1(kj72精选课件 然后用中心差商逼近这些导数值，则可得到Wendroff格式：在点 处，用)21,21(kjPP1kttktt CHADGFBEjh
38、x hjx)1(GEPtututu21FHPxuxuxu21,2,1,2,1,0022111111111jkhuuhuuauuuukjkjkjkjkjkjkjkj）（）（73精选课件连同初始条件与边界条件：11111()1kkkkjjjjaruuuuar 该格式的局部截断误差阶为 ，且无条件稳定。22()Oh令 ，格式可改写为/rh该格式可在0 x,t 内所有网点上显示地计算解之近似值。00,1,2,0,1,2,kkjjuujk74精选课件 222220uuaxatx 为正常数考察 对x-t平面进行矩形网格剖分，x方向的步长为h，t方向的步长为，网点 简记为(j,k)。),(kjtx 用二阶中
39、心差商代替(3.3)中的二阶导数，则得到网点(j,k)处的差分方程：(3.3)(3.4)(,0)()(,0)()uu xxxtxt 222221kktjxjauuh12222111()2(1)kkkkkjjjjjua ruua ruu,2,1,0,2,1,0;/kjhr其中 。或 1kk1k1jj1j75精选课件该格式稳定的充分条件为 。1ar初始条件)()0,(xxu离散：)()0,(txtu初始条件 离散：由0222022111)(21jxjtjjjuhauuu消去 ，得1jujjjjjararu)1(22211221)(22hO 上述差分格式与初始条件的截断误差阶均为 。jju0取为76
40、精选课件Ttxxuatu0,10022222上述方法也可用于求解初边值问题：10)()0,()()0,(xtxtuxxuTtttuttu0)(),1()(),0(77精选课件 78精选课件79精选课件80精选课件在条件下2122ssldsdsdydsdx求使得泛函达到最大的函数 。)(),(sysx2121),(ssdsdsdxydsdyxyxs 在长度一定的所有平面封闭曲线中，求所围面积为最大的曲线。81精选课件当求泛函在一个函数集合K中的极小（或极大）问题，则该问题称为变分问题。考察J(x)的极值情况。为实常数）fLLfxLxxJ,0(2)(2设求 ，使Rx 0)(min)(0 xJxJR
41、x与求解方程 Lx=f 等价。82精选课件对称正定nnnnnnaaaaaaaaaA212222111211求J(x)取极小值的驻点，其中 ninjniiijiijnxbxxaxxxJxJ1112121),()(设设 则J(x)可表示为：Tnxxxx),(21Tnbbbb),(21),(),(21)(xbxxAxJ83精选课件设矩阵A对称正定，则下列两个命题等价：求 ，使nRx 0)(min)(0 xJxJnRx(a)(b)是方程 的解0 xbxA：),(),(21)(xbxxAxJ其中 。84精选课件 考察一根长为 l 的弦，两端固定在点 A(0,0)和B(l,0)。当没有外力作用时，它的位置
42、沿水平方向与X轴重合。设有强度为f(x)的外荷载垂直向下作用在弦上，于是弦发生形变。假定荷载很小，因而发生的形变也很小。用 u(x)表示在荷载f(x)的作用下弦的平衡位置。u)0,0(A)0,(lBX85精选课件求弦的平衡位置归结为求解两点边值问题：设弦处于某一位置u=u(x)，可得到其总位能为 ldxufuTuJ02)2(21)(：0)(0)0(0)()(luulxxfxuT其中T是弦的张力。)(xuu 弦的平衡位置(记为 )将在满足边值条件 u(0)=0，u(l)=0 的一切可能位置中，使位能取极小值。弦的平衡位置 是下列变分问题的解)(xuu)(min)(20uJuJCu86精选课件 在
43、数学上，要将某个微分方程的定解问题转化为一个变分问题求解，必须针对已给的定解问题构造一个相应的泛函，并证明定解问题的解与泛函极值问题的解等价。87精选课件考察二阶常微分方程边值问题：0)(0)(),(buaubaxfqudxdupdxdLu),(),(21)(ufuuuJ),(),(21)(ufuLuuJbababafudxdxquudxdxdupdxd221badxfuquup)2(2122dxquvdxdvdxdupvuba),(引入泛函算子则88精选课件 与前述二阶常微分方程边值问题相应的变分问题是),(),(21)(ufuuuJ其中0)(,)()(,11aubaHxuxubaHE求 ，
44、使1EHu)(min)(1uJuJEHubadxfuquup)2(212289精选课件0)(,0)(),(buaubaxfqudxdupdxdLu设 ，是边值问题)(0ICf 2*Cu 的解，则 使 J(u)达到极小值；*u12EHCu 反之，若 使 J(u)达到极小值，则 是边值问题的解。*u),(),(21)(ufuuuJ其中 是强制边界条件，是自然边界条件，区别这两类边界条件在用有限元方法求解边值问题时很重要。0)(bu0)(au90精选课件对两点边值问题：0)(,0)(),(buaubaxfqudxdupdxdLu其中1*EHu，且满足变分方程：设 ，以v乘方程两端，沿a,b积分，并利
45、用 ，得变分方程1EHv0)(,0)(bvavdxquvdxdvdxdupvuba),(0),(),(vfvu对任意1EHv0),(),(*vfvu在力学里，表示虚功),(),(vfvu1EHv2*Cu 设 ，则 是边值问题解的充要条件是：*u91精选课件 。若将微分方程化为相应的变分问题或变分方程，则只需处理强加边界条件，无需处理自然边界条件（自然边界条件已包含于变分问题中泛函的构造或已包含于给出的变分方程之中）。这一特点对研究微分方程离散化方法及其数值解带来了极大的方便。92精选课件模型方程其中G是平面有界区域。0),(),(|2222uGyxyxfyuxuu),(),(21)(ufuuu
46、J),(),(21)(ufuuuJGdxdyyvyuxvxuvu)(),(引入泛函算子则93精选课件与前述二阶椭圆边值问题相应的变分问题是求 ，使)(10GHu)(min)()(10uJuJGHu其中0),(),(),(|),()(110yxuGHyxuyxuGH),(),(21)(ufuuuJ),(),(21ufuu94精选课件 对，泛函是一样的，只是边界条件要作为强加边值条件加在所取的函数类上。设 ，是二阶椭圆边值问题的解，则 使 J(u)达到极小值；)(0ICf)(2*GCu*u)()(102GHGCu 反之，若 使 J(u)达到极小值，则 是二阶椭圆边值问题的解。*u),(),(21)
47、(ufuuuJ其中 对，二次泛函形式相对于第一边值问题有所改变，但函数类的选取与边界条件无关。95精选课件问题000),(),(21|aaunuuGyxyxfu其中 设 ，以v乘方程两端后在G上积分，并利用Green公式，得变分方程)(1GHvE2),(auvdsdxdyyvyuxvxuvuG0),(),(vfvu0),(),(),(|),()(111yxuGHyxuyxuGHE)(1GHvE96精选课件在力学里，表示虚功),(),(vfvu)(2GCu 设 是边值问题的解，则对任意 ，满足变分方程。)(1GHvEu 反之，若 ，且对任意 满足变分方程，则 为边值问题的解。)()(12GHGC
48、uE)(1GHvEu 与极小位能原理类似，第一类边界条件为强加边界条件，第二、三类边界条件为自然边界条件。虚功原理比极小位能原理应用更广。97精选课件：求解相应的变分问题或相应的变分方程。这两种算法统称为Ritz-Galerkin方法。以下用V表示 等Sobolev空间，L表示微分算子，(u,v)为由L及边值条件决定的双线性泛函。110,HHHE 用有限维空间的函数代替变分问题(或变分方程)中无限维空间的函数，从而在有限维函数空间中求变分问题(或变分方程)的近似解，并要求当有限维空间的维数不断增加时，有限维近似解逼近原变分问题(或变分方程)的解。98精选课件由极小位能原理得出的变分问题为：ni
49、iincu1求 ，使Vu)(min)(uJuJVu其中 ，),(),(21)(ufuuuJ 设 是V 的n维子空间，是 的一组基底(称为基函数)。中任一元素 可表示为nVn,21nVnVnu即选择适当的 ，使 取极小值。nccc,21)(nuJ求 ，使nnVu)(min)(vJuJnVvn：99精选课件展开)(nuJ令nccc,21则 满足解出 代入 ，则得),1(nicinuniiincu1njfcjniiji,2,1),(),(1ninjnjjjjijicfcc111),(),(21njcuJjn,2,10)(),(),(21)(nnnnufuuuJ100精选课件 根据最小位能原理构造相应
50、于微分方程或物 理问题的变分问题；取 作为 的一组基底，即用 近似代替无穷维空间V；),1(nii,21nnspanVnV 根据二次函数取极值的必要条件，得到 nniiinVcu1ic中 所满足的方程组：求解关于 的线性代数方程组。icnjfcjniiji,2,1),(),(1101精选课件由虚功原理得出的变分方程为：niiincu1 设 是V 的n维子空间，是 的一组基底(称为基函数)。中任一元素 可表示为nVn,21nVnVnu：0),(),(vfvu0),(),(vfvun求 ，使对 ，满足nnVu nVvnu102精选课件nVv由 的任意性，取 作为v，则得n,21将 代入变分方程，则
