偏微分方程数值解法试题与答案

偏微分方程数值解试题(０6B ）参考答案与评分标准信息与计算科学专业一(10分）、设矩阵A 对称，定义)(),(),(21)(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),()()()(2000x 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 ,)(),(21)0()1()(00x J x Ax x x J >+==+ϕϕ,因此0x 是)(x J 的最小值点. (４分）评分标准：)(λϕ的展开式3分, 每问3分，推理逻辑性１分二（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 ∈,乘方程两端,积分应用分部积分得到 (３分))().(),(v f fvdx dx quv dxdv dx du p v u a b a ba ==+=⎰⎰,),(1b a H v E ∈∀ 即变分问题的Galerkin 形式. （3分)令⎰-+=-=b a dx fu qu dxdu p u f u u a u J ])([21),(),(21)(22，则变分问题的Ritz 形式为求),(1*b a H u E ∈,使)(m in )(1*u J u J EH u ∈= （4分) 评分标准:空间描述与积分步骤３分，变分方程3分，极小函数及其变分问题４分,三(20分)、对于边值问题⎪⎩⎪⎨⎧-====⨯=∈=∂∂+∂∂====x u u u u G y x y u x u y y x x 1||,0|,1|)1,0()1,0(),(,010102222 (１)建立该边值问题的五点差分格式(五点棱形格式又称正五点格式),推导截断误差的阶。

一．填空（1553=⨯分）1．若步长趋于零时，差分方程的截断误差0→lmR ，则差分方程的解lm U 趋近于微分方程的解lm u . 此结论_______（错或对）； 2．一阶Sobolev 空间{})(,,),()(21Ω∈''=ΩL f f f y x f H y x关于内积=1),(g f _____________________是Hilbert 空间；3．对非线性（变系数）差分格式，常用 _______系数法讨论差分格式的_______稳定性； 4．写出3x y =在区间]2,1[上的两个一阶广义导数：_________________________________， ________________________________________；5．隐式差分格式关于初值是无条件稳定的. 此结论_______（错或对）。

二．（13分）设有椭圆型方程边值问题用1.0=h 作正方形网格剖分 。

（1）用五点菱形差分格式将微分方程在内点离散化； （2）用截断误差为)(2h O 的差分法将第三边界条件离散化； （3）整理后的差分方程组为 三．（12）给定初值问题xut u ∂∂=∂∂ ， ()10,+=x x u 取时间步长1.0=τ，空间步长2.0=h 。

试合理选用一阶偏心差分格式（最简显格式）, 并以此格式求出解函数),(t x u 在2.0,2.0=-=t x 处的近似值。

1．所选用的差分格式是： 2．计算所求近似值：四．（12分）试讨论差分方程()ha h a r u u r u u k l k l k l k l ττ+-=-+=++++11,1111逼近微分方程0=∂∂+∂∂xu a t u 的截断误差阶R 。

思路一：将r 带入到原式，展开后可得格式是在点（l+1/2,k+1/2）展开的。

思路二：差分格式的用到的四个点刚好是矩形区域的四个顶点，可由此构造中心点的差分格式。

偏微分方程数值解一（10分）、设矩阵A 对称正定，定义)(),(),(21)(n R x x b x Ax x J ∈-=，证明下列两个问题等价：(1)求n R x ∈0使)(min )(0x J x J n Rx ∈=;(2)求下列方程组的解：b Ax = 解: 设n R x ∈0是)(x J 的最小值点,对于任意的n R x ∈,令),(2),()()()(2000x Ax x b Ax x J x x J λλλλϕ+-+=+=, (3分)因此0=λ是)(λϕ的极小值点,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分) 反之,若nR x ∈0满足bAx =0,则对于任意的x ,)(),(21)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 dxdu 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==∈=b u a u b a H u u H 为求解函数空间,检验函数空间.取),(10b a H v ∈,乘方程两端,积分应用分部积分得到 (3分))().(),(v f fvdx dx quv dxdv dx du p v u a b a ba ==+=⎰⎰,),(1b a H v ∈∀ 即变分问题的Galerkin 形式. (3分)令⎰-+=-=b a dx fu qu dxdup u f u u a u J ])([21),(),(21)(22,则变分问题的Ritz 形式为求),(1*b a H u ∈,使)(m in )(10*u J u J H u ∈= (4分) 评分标准:空间描述与积分步骤3分,变分方程3分,极小函数及其变分问题4分,三（20分）、对于边值问题⎪⎩⎪⎨⎧=⨯=∈-=∂∂+∂∂∂0|)1,0()1,0(),(,12222G u G y x yux u （1）建立该边值问题的五点差分格式（五点棱形格式又称正五点格式），推导截断误差的阶。

偏微分⽅程数值习题解答李微分⽅程数值解习题解答 1-1 如果0)0('=?，则称0x 是)(x J 的驻点（或稳定点）.矩阵A 对称（不必正定），求证0x 是)(x J 的驻点的充要条件是：0x 是⽅程组 b Ax =的解证明:由)(λ?的定义与内积的性线性性质,得),()),((21)()(0000x x b x x x x A x x J λλλλλ?+-++=+=),(2),()(200x Ax x b Ax x J λλ+-+=),(),()(0'x Ax x b Ax λλ?+-=必要性:由0)0('=?,得,对于任何n R x ∈,有0),(0=-x b Ax ,由线性代数结论知,b Ax b Ax ==-00,0充分性: 由b Ax =0,对于任何n R x ∈,0|),(),()0(00'=+-==λλ?x Ax x b Ax即0x 是)(x J 的驻点. §1-2补充: 证明)(x f 的不同的⼴义导数⼏乎处处相等.证明:设)(2I L f ∈,)(,221I L g g ∈为)(x f 的⼴义导数,由⼴义导数的定义可知,对于任意)()(0I C x ∞∈?,有-=ba ba dx x x f dx x x g )()()()('1?? ??-=ba ba dx x x f dx x x g )()()()('2?? 两式相减,得到)(0)()(021I C x g g ba ∞∈?=- 由变分基本引理,21g g -⼏乎处处为零,即21,g g ⼏乎处处相等.补充:证明),(v u a 的连续性条件(1.2.21) 证明: 设'|)(|,|)(|M x q M x p ≤≤,由Schwarz 不等式||||.||||||||.|||||)(||),(|'''''v u M v u M dx quv v pu v u a ba +≤+=?11*||||.||||2v u M ≤,其中},max{'*M M M =习题:1 设)('x f 为)(x f 的⼀阶⼴义导数,试⽤类似的⽅法定义)(x f 的k 阶导数,...2,1(=k ) 解:⼀阶⼴义导数的定义,主要是从经典导数经过分部积分得到的关系式来定义,因此可得到如下定义:对于)()(2I L x f ∈,若有)()(2I L x g ∈,使得对于任意的)(0I C ∞∈?,有 ?-=bak kba dx x x f dx x x g )()()1()()()(??则称)(x f 有k 阶⼴义导数,)(x g 称为)(x f 的k 阶⼴义导数,并记kk dxfd x g =)(注:⾼阶⼴义导数不是通过递推定义的,可能有⾼阶导数⽽没有低阶导数.2.利⽤)(2I L 的完全性证明))()((1I H I H m 是Hilbert 空间.证明:只证)(1I H 的完全性.设}{n f 为)(1I H 的基本列,即0||||||||||||0''01→-+-=-m n m n m n f f f f f f因此知}{},{'n n f f 都是)(2I L 中的基本列(按)(2I L 的范数).由)(2I L 的完全性,存在)(,2I L g f ∈,使0||||,0||||0'0→-→-g f f f n n ,以下证明0||||1→-f f n (关键证明dxdfg =)由Schwarz 不等式,有00||||.|||||)())()((|??f f x x f x f n ba n -≤-?00'''|||||||||)())()((|??f f dx x x g x f n ba n -≤-?对于任意的)()(0I C x ∞∈?,成⽴=∞a ba n n dx x x f dx x x f )()()()(lim ??=∞→ba b a nn dx x x g dx x x f )()()()(lim '??由?-=ba n ba ndx x x f dx x x f )()()()(''??取极限得到dx x x f dx x x g ba ba ??-=)()()()('??即')(f x g =,即)(1I H f ∈,且0||||||||||||0''01→-+-=-f f f f f f n n n故)(1I H 中的基本列是收敛的,)(1I H 是完全的. 3.证明⾮齐次两点边值问题证明:边界条件齐次化令)()(0a x x u -+=βα,则0u u w -=满⾜齐次边界条件.w 满⾜的⽅程为00Lu f Lu Lu Lw -=-=,即w 对应的边值问题为==-=0)(,0)('b w a w Lu f Lw (P) 由定理知,问题P 与下列变分问题等价求)(min )(,**12*1w J w J H C w Ew E ∈=∈其中),(),(21)(0*w Lu f w w a w J --=.⽽Cu u a u Lu u J u u Lu f u u u u a w J +-+=-----=),(),()(~),(),(21)(000000*⽽200)()(),(),(C b u b p u u a u Lu +-=-β从⽽**)()()(~)(C b u b p u Jw J +-=β则关于w 的变分问题P 等价于:求α=∈)(,12*a u H C u使得)(min )()(*1u J u J a u H u α=∈=其中)()(),(),(21)(b u b p u f u u a u J β--=4就边值问题（1.2.28）建⽴虚功原理解:令)(0a x u -+=βα,0u u w -=,则w 满⾜)(,0)('00==-=-=b w a w Lu f Lu Lu Lw等价于:1E H v ∈?0),(),(0=--v Lu f v Lw应⽤分部积分,+-=-=-b a b a b a dx dx dv dx dw p v dx dw p vdx dx du p dx d v dx dw p dx d |)()),((还原u ,)()(),(),(),(),(),(),(),(),(000b v b p v f v u a v u a v Lu v f v u a v Lu f v w a β--=-+-=--于是,边值问题等价于:求α=∈)(,1a u H u ,使得1E H v ∈?,成⽴0)()(),(),(=--b v b p v f v u a β注:形式上与⽤v 去乘⽅程两端,应⽤分部积分得到的相同. 5试建⽴与边值问题等价的变分问题.解:取解函数空间为)(20I H ,对于任意)(20I H v ∈⽤v 乘⽅程两端,应⽤分部积分,得到0),(),(44=-+=-v f u dx ud v f Lu⽽??-==b a b a b a dx dxdvdx u d v dx u d vdx dx u d v dx u d .|),(33334444 dx dxv d dx u d dx dx vd dx u d dx dv dx u d b a b a b a ??=+-=2222222222| 上式为),(][2222v f dx uv dx vd dx u d b a =+?定义dx uv dxvd dx u d v u a ba ][),(2222+=?,为双线性形式.变分问题为:求)(20I H u ∈,)(20I H v ∈?),(),(v f v u a =1-41.⽤Galerkin Ritz -⽅法求边值问题==<<=+-1)1(,0)0(102"u u x x u u 的第n 次近似)(x u n ,基函数n i x i x i ,...,2,1),sin()(==π?解:(1)边界条件齐次化:令x u =0,0u u w -=,则w 满⾜齐次边界条件,且)1(,0)0(20==-=-=w w x x Lu Lu Lw第n 次近似n w 取为∑==n i i i n c w 1,其中),...2,1(n i c i =满⾜的Galerkin Ritz -⽅程为n j x x c a j ni i j i ,...,2,1),(),(21=-=∑= ⼜xd jx ix ij dx x j x i dxx j x i ij dx a j i jij i ?-=+=+=ππππππππ)cos()cos(2)sin()sin()cos()cos()(),(1010210''-+πππjx ix sin sin 21由三⾓函数的正交性,得到≠=+=j i j i i a j i ,0,212),(22π??⽽]1)1[()(2)sin()1(),(3102--=-=-?jj j dx x j x x x x ππ? 于是得到+-=-=为偶数为奇数j j j j a x x c j j j j 0 )1()(8),(),(2232ππ最后得到∑+=-+---+=]21[1233])12(1[)12(])12sin[(8)(n k n k k x k x x u ππ 2.在题1中,⽤0)1(=u 代替右边值条件,)(x u n 是⽤Galerkin Ritz -⽅法求解相应问题的第n 次近似,证明)(x u n 按)1,0(2L 收敛到)(x u ,并估计误差.证明:n u 对应的级数绝对收敛,由}{sin x i π的完全性知极限就是解)(x u ,其误差估计为338nR n π≤3.就边值问题(1.2.28)和基函数),...,2,1()()(n i a x x i i =-=?,写出Galerkin Ritz -⽅程解:边界条件齐次化,取)(0a x u -+=βα,0u u w -=, w 对应的微分⽅程为)(,0)('00==-=-=b w a w Lu f Lu Lu Lw对应的变分⽅程为0),(),(0=--v Lu f v w a)]([)(000a x q dx dpqu dx du p dx d Lu -++-=+-=βαβ+-=-ba b a dx x pv b v b p v dxdp )()()(' 变分⽅程为dx v qu x pv b v b p v f v w a ba ?--+=])([)()(),(),(0'ββ取n i a x x i i ,...,2,1,)()(=-=?,则Galerkin -Ritz ⽅程为∑-++--+=-=ba i ba i i nj j jidxa x x q dx a x i x pb b p fc a )]()[()()()()(),(),(11βαβ?β??+=ba j i j i j i dx q p a ][),(''取1,0,1===f q p ,具体计算1=n , )(1),(11a b dx a ba -==221)(21)()()(21a b a b a b a b d -=---+-=ββ, )(211a b c -=,即解)(2101a x u u -+= 2=n :22111)()(2),(),(),(a b dx a x a a b a ba -=-=-=3222)(34)(4),(a b dx a x a ba -=-=3223222)(31)()()(31)(2)()(a b a b a b a b dxa x ab dx a x d ba b a -=---+-=---+-=??ββββ得到⽅程组为 --=----3221322)(31)(21c )(34)()(a b a b c a b a b a b a b特别取1,0==b a ,有= 31213411121c c求解得到1,21,6131122=-=-=c c c其解为202)(21)(a x a x u u ---+=C h2 椭圆与抛物型⽅程有限元法§1.1 ⽤线性元求下列边值问题的数值解: 10,2sin242"<<=+-x x y y ππ0)1(,0)0('==y y此题改为4/1,0)1()0(,1"====+-h y y y y解: 取2/1=h ,)2,1,0(==j jh x j ,21,y y 为未知数. Galerkin 形式的变分⽅程为),(),(v f v Lu =,其中+-=10210"4),(uvdx vdx u v Lu π,?=1)(2sin 2),(dx x xv v f π⼜dx v u dx v u v u vdx u =+-=-10''10''10'10"|因此dx uv v u v u a )4(),(12''?+=π在单元],[1i i i x x I -=中,应⽤仿射变换(局部坐标)hx x i 1--=ξ节点基函数为)3,2,1(,0,,,1)(111=≤≤-=≤≤-=-=--+i other x x x h x x x x x h x x x i i i i i i i ξξξξ?-+++=++=1022210222222'111)1(41]41[]4[),(1021ξξπξξπ?πd h d hh dxa x x x x取2/1=h ,则计算得124),(211π??+=a122)1(41[),(210221πξξξπ??+-=-+-=?d h h a-+++=10101)1)(2121(2sin )0(2sin [2),(ξξξπξξξπ?d d h h f ??-++=1010)1(4)1(sin 2sin ξξξπξξξπd d hξξξπ?d h f ?+=102)2121(2sin 2),(代数⽅程组为= ),(),(),(),(),(),(212122212111f f y y a a a a 代如求值.取4/1=h ,未知节点值为4321,,,u u u u ,⽅程为4,3,2,1),(),(41==∑=j f ua j i iji应⽤局部坐标ξ表⽰,-+++=10221022])1(41[)41(),(ξξπξξπ??d hh d h h a j j248]88[21022πξξπ+=+=?dξξξπ??d hh a j j ])1(41[),(1021?-+-=++964)1(164212πξξξπ+-=-+-=?d 964),(21π??+-=-j j a系数矩阵为}964,248,964{222πππ+-++-=diag A取1=f ,41)1(),(1010=-+=??ξξξξ?d h d h f j-+++=+10110)1)]((2sin[2)](2sin[2),(ξξξπξξξπd h x h d h x h f j j j -++++=1010)1)](4 41(2sin[21)]44(2sin[42ξξξπξξξπd j d j++?=+++++-+=100110|)]8)1([cos(821]8)1(sin[21]8)1(sin[]8)(sin[21ξππξξπξξξπξπj d j d j j+2.就⾮齐次第三边值条件22'11')()(,)()(βαβα=+=+b u b u a u a u导出有限元⽅程.解:设⽅程为f qu pu Lu =+-='')( 则由),()]()[()()]()[()(),(|),)((''1122'''''v pu a u a v a p b u b v b p v pu v pu v pu b a----=-=αβαβ变分形式为:),(1b a H v ∈?)()()()(),()()()()()()(),(),(1212''a v a p b v b p v f a v a u a p b v b u b p v qu v pu ββαα-+=-++)(),(0b u u a u u N ==记)()()()(),()()()()()()()(),(),(),(1212''a v a p b v b p v f v F a v a u a p b v b u b p v qu v pu v u A ββαα-+=-++=则上述变分形式可表⽰为)(),(v F v u A =设节点基函数为),...,2,1,0)((N j x j =? 则有限元⽅程为),...,1,0()(),(0N j F u A j Ni i j i ==∑=具体计算使⽤标准坐标ξ.。

偏微分方程习题及答案【篇一：偏微分方程数值解法期末考试题答案】题答案及评分标准学年学期：专业：班级：课程：教学大纲：使用教材：教材作者：出版社：数学与应用数学数学偏微分方程数值解法《偏微分方程数值解法》教学大纲（自编，2006）《偏微分方程数值解法》陆金甫、关治清华大学出版社一、判断题（每小题1分，共10分）1、（o）2、（o）3、（x）4、（x）5、（o）6、（o）7、（o）8、（x）9、（x） 10、（o）二、选择题（每小题2分，共10分） 11、（d） 12、（a） 13、（c） 14、（b）15、（c）三、填空题（每小题2分，共20分）?2?216、2?2??x1?x2?2?2 17、a=[4 5 9;23 5 17;11 23 1] 18、y=exp(-t/3)*sin(3*t) ?xn19、help 20、zeros(m,n)21、inva(a)*b或者a/b22、a=sym([cos(x-y),sin(x+y);exp(x-y),(x-1)^3])?(s)?1?(s)?c[??(s)]2?023、a[?2(s)]2?2b?224????v(?)ed? 25、i?xu(xj,tn?1)?u(xj,tn)?四、计算题：（每小题12分，共36分）?u?u?0（x?r,t?0）的有限差分方程（两层显示26、写成对流方程?a?t?x格式，用第n层计算第n+1层），并把有限差分方程改写为便于计算的迭代格式???/h为网格比。

解：在点(xj,tn)处，差分方程为?1un?unjj??anunj?1?ujh?0（j?0,?1,?2,，n?0,1,2,）（8分）便于计算的形式为?1nnn???/h （4分） un?u?a?(u?ujjj?1j)，?u?2u?a2的有限差分方程（中心差分格式，用第n层27、写出扩散方程?t?x计算第n+1层），并把有限差分方程改写为便于计算的迭代格式，???/h2为网格比。

《偏微分方程数值解》试卷参考答案与评分标准专业班级信息与计算科学开课系室考试日期 2006.4.14命题教师王子亭偏微分方程数值解试题(06A)参考答案与评分标准信息与计算科学专业一（10分）、设矩阵A 对称正定，定义)(),(),(21)(n R x x b x Ax x J ∈-=，证明下列两个问题等价：(1)求n R x ∈0使 )(min )(0x J x J nRx ∈=;(2)求下列方程组的解：b Ax =解: 设n R x ∈0是)(x J 的最小值点,对于任意的n R x ∈,令),(2),()()()(2000x Ax x b Ax x J x x J λλλλϕ+-+=+=, (3分)因此0=λ是)(λϕ的极小值点,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分) 反之,若nR x ∈0满足bAx =0,则对于任意的x ,)(),(21)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 dxdu 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)()(),,(|{110==∈=b u a u b a H u u H 为求解函数空间,检验函数空间.取),(10b a H v ∈,乘方程两端,积分应用分部积分得到 (3分))().(),(v f fvdx dx quv dxdv dx du pv u a b a ba ==+=⎰⎰,),(10b a H v ∈∀ 即变分问题的Galerkin 形式. (3分)令⎰-+=-=b a dx fu qu dxdup u f u u a u J ])([21),(),(21)(22,则变分问题的Ritz 形式为求),(10*b a H u ∈,使)(min )(1*u J u J H u ∈= (4分)评分标准:空间描述与积分步骤3分,变分方程3分,极小函数及其变分问题4分,三（20分）、对于边值问题⎪⎩⎪⎨⎧=⨯=∈-=∂∂+∂∂∂0|)1,0()1,0(),(,12222G u G y x yux u （1）建立该边值问题的五点差分格式（五点棱形格式又称正五点格式），推导截断误差的阶。

偏微分方程数值解法（带程序）例1 求解初边值问题22,(0,1),012,(0,]2(,0)12(1),[,1)2(0,)(1,)0,0u ux t t x x x u x x x u t u t t ⎧⎪⎪⎪⎧⎨⎪⎪⎪⎨⎪⎪⎪⎪⎩⎩∂∂=∈>∂∂∈=-∈==>要求采用树脂格式 111(2)n n n n nj j j j j u u u u u λ++-=+-+，2()tx λ∆=∆，完成下列计算： （1） 取0.1,0.1,x λ∆==分别计算0.01,0.02,0.1,t =时刻的数值解。

（2） 取0.1,0.5,x λ∆==分别计算0.01,0.02,0.1,t =时刻的数值解。

（3） 取0.1, 1.0,x λ∆==分别计算0.01,0.02,0.1,t =时刻的数值解。

并与解析解22()22181(,)sin()sin()2n t n u n x t n x e n ππππ∞-==∑进行比较。

解：程序function A=zhongxinchafen(x,y,la) U=zeros(length(x),length(y)); for i=1:size(x,2)if x(i)>0&x(i)<=0.5 U(i,1)=2*x(i); elseif x(i)>0.5&x(i)<1 U(i,1)=2*(1-x(i)); end endfor j=1:length(y)-1for i=1:length(x)-2U(i+1,j+1)=U(i+1,j)+la*(U(i+2,j)-2*U(i+1,j)+U(i,j)); end endA=U(:,size(U,2))function u=jiexijie1(x,t) for i=1:size(x,2) k=3;a1=(1/(1^2)*sin(1*pi/2)*sin(1*pi*x(i))*exp(-1^2*pi^2*t));a2=a1+(1/(2^2)*sin(2*pi/2)*sin(2*pi*x(i))*exp(-2^2*pi^2*t));while abs(a2-a1)>0.00001a1=a2;a2=a1+(1/(k^2)*sin(k*pi/2)*sin(k*pi*x(i))*exp(-k^2*pi^2*t));k=k+1;endu(i)=8/(pi^2)*a2;endclc; %第1题第1问clear;t1=0.01;t2=0.02;t3=0.1;x=[0:0.1:1];y1=[0:0.001:t1];y2=[0:0.001:t2];y3=[0:0.001:t3];la=0.1;subplot(131)A1=zhongxinchafen(x,y1,la);u1=jiexijie1(x,t1)line(x,A1,'color','r','linestyle',':','linewidth',1.5);hold online(x,u1,'color','b','linewidth',1);A2=zhongxinchafen(x,y2,la);u2=jiexijie1(x,t2)line(x,A2,'color','r','linestyle',':','linewidth',1.5);line(x,u2,'color','b','linewidth',1);A3=zhongxinchafen(x,y3,la);u3=jiexijie1(x,t3)line(x,A3,'color','r','linestyle',':','linewidth',1.5);line(x,u3,'color','b','linewidth',1); title('例1(1)');subplot(132);line(x,u1,'color','b','linewidth',1);line(x,u2,'color','b','linewidth',1);line(x,u3,'color','b','linewidth',1);title('解析解');subplot(133);line(x,A1,'color','r','linestyle',':','linewidth',1.5);line(x,A2,'color','r','linestyle',':','linewidth',1.5);line(x,A3,'color','r','linestyle',':','linewidth',1.5);title('数值解');clc; %第1题第2问clear;t1=0.01;t2=0.02;t3=0.1;x=[0:0.1:1];y1=[0:0.005:t1];y2=[0:0.005:t2];y3=[0:0.005:t3];la=0.5;subplot(131);A1=zhongxinchafen(x,y1,la);u1=jiexijie1(x,t1)line(x,A1,'color','r','linestyle',':','linewidth',1.5);hold online(x,u1,'color','b','linewidth',1);A2=zhongxinchafen(x,y2,la);u2=jiexijie1(x,t2)line(x,A2,'color','r','linestyle',':','linewidth',1.5);line(x,u2,'color','b','linewidth',1);A3=zhongxinchafen(x,y3,la);u3=jiexijie1(x,t3)line(x,A3,'color','r','linestyle',':','linewidth',1.5); line(x,u3,'color','b','linewidth',1);title('例1(2)'); subplot(132);line(x,u1,'color','b','linewidth',1); line(x,u2,'color','b','linewidth',1);line(x,u3,'color','b','linewidth',1);title('解析解'); subplot(133);line(x,A1,'color','r','linestyle',':','linewidth',1.5); line(x,A2,'color','r','linestyle',':','linewidth',1.5);line(x,A3,'color','r','linestyle',':','linewidth',1.5);title('数值解');clc; %第1题第3问 clear;t1=0.01;t2=0.02;t3=0.1;x=[0:0.1:1];y1=[0:0.01:t1];y2=[0:0.01:t2];y3=[0:0.01:t3];la=1.0; subplot(131);A1=zhongxinchafen(x,y1,la);u1=jiexijie1(x,t1)line(x,A1,'color','r','linestyle',':','linewidth',1.5);hold on line(x,u1,'color','b','linewidth',1);A2=zhongxinchafen(x,y2,la);u2=jiexijie1(x,t2) line(x,A2,'color','r','linestyle',':','linewidth',1.5); line(x,u2,'color','b','linewidth',1);A3=zhongxinchafen(x,y3,la);u3=jiexijie1(x,t3) line(x,A3,'color','r','linestyle',':','linewidth',1.5); line(x,u3,'color','b','linewidth',1);title('例1(3)'); subplot(132);line(x,u1,'color','b','linewidth',1); line(x,u2,'color','b','linewidth',1);line(x,u3,'color','b','linewidth',1);title('解析解'); subplot(133);line(x,A1,'color','r','linestyle',':','linewidth',1.5); line(x,A2,'color','r','linestyle',':','linewidth',1.5);line(x,A3,'color','r','linestyle',':','linewidth',1.5);title('数值解'); 运行结果：表1：取0.1,0.1,x λ∆==0.01,0.02,0.1,t =时刻的解析解与数值解表2：取0.1,0.5,x λ∆==0.01,0.02,0.1,t =时刻的解析解与数值解表3：取0.1, 1.0,x λ∆==0.01,0.02,0.1,t =时刻的解析解与数值解图1：取0.1,0.1,x λ∆==0.01,0.02,0.1,t =时刻的解析解与数值解图2：取0.1,0.5,x λ∆==0.01,0.02,0.1,t =时刻的解析解与数值解图3：取0.1, 1.0,x λ∆==0.01,0.02,0.1,t =时刻的解析解与数值解例2 用Crank-Nicolson 格式完成例1的所有任务。

偏微分方程考试题及答案一、单项选择题（每题2分，共10分）1. 偏微分方程的一般形式是什么？A. \( \frac{\partial u}{\partial x} + \frac{\partialu}{\partial y} = 0 \)B. \( \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} = 0 \)C. \( \frac{\partial u}{\partial x} + \frac{\partial^2u}{\partial y^2} = 0 \)D. \( \frac{\partial^2 u}{\partial x^2} + \frac{\partialu}{\partial y} = 0 \)答案：B2. 以下哪个方程不是线性偏微分方程？A. \( \frac{\partial^2 u}{\partial x^2} + 3u = 0 \)B. \( \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} = 1 \)C. \( \frac{\partial^2 u}{\partial x^2} + 2\frac{\partial u}{\partial y} = 0 \)D. \( u^2 \frac{\partial^2 u}{\partial x^2} + \frac{\partial u}{\partial y} = 0 \)答案：D3. 波动方程的解通常表示为两个函数的和，这两个函数分别是？A. \( f(x+ct) \) 和 \( g(x-ct) \)B. \( f(x+ct) \) 和 \( g(x+ct) \)C. \( f(x-ct) \) 和 \( g(x+ct) \)D. \( f(x+ct) \) 和 \( h(x-ct) \)答案：A4. 拉普拉斯方程的解是调和函数，以下哪个条件不是调和函数必须满足的？A. \( \Delta u = 0 \)B. \( \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2u}{\partial y^2} = 0 \)C. \( \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2u}{\partial y^2} + \frac{\partial^2 u}{\partial z^2} = 0 \) D. \( \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2u}{\partial y^2} = 1 \)答案：D5. 以下哪个条件不是偏微分方程解的存在性和唯一性定理所要求的？A. 初始条件B. 边界条件C. 系数的连续性D. 变量的离散性答案：D二、填空题（每题3分，共15分）1. 偏微分方程 \( \frac{\partial^2 u}{\partial x^2} -\frac{\partial^2 u}{\partial y^2} = 0 \) 是一个 ________ 方程。

1. 课本2p 有证明2. 课本812,p p 有说明3. 课本1520,p p 有说明4. Rit2法，设n u 是u 的n 维子空间，12,...n ϕϕϕ是n u 的一组基底，n u 中的任一元素n u 可表为1nn i i i u c ϕ==∑，则,1111()(,)(,)(,)(,)22j nnn 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 jJ u c ∂=∂，从而得到12,...n c c c 满足1(,)(,),1,2...niji j i a c f j n ϕϕϕ===∑，通过解线性方程组，求的i c ，代入1nn i i i u c ϕ==∑，从而得到近似解n u 的过程称为Rit2法简而言之，Rit2法：为得到偏微分方程的有穷维解，构造了一个近似解，1nn i ii u c ϕ==∑，利用,1111()(,)(,)(,)(,)22j nnn 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 法：为求得1nn i ii u c ϕ==∑形式的近似解，在系数i c 使n u 关于n V u ∈，满足(,)(,)n a u V f V =，对任意nV u ∈或（取,1j V j nϕ=≤≤）1(,)(,),1,2...nijij i a cf j n ϕϕϕ===∑的情况下确定i c ，从而得到近似解1nn i i i u c ϕ==∑的过程称Galerkin 法为 Rit2-Galerkin 法方程：1(,)(,)nijij i a cf ϕϕϕ==∑5. 有限元法：将偏微分方程转化为变分形式，选定单元的形状，对求解域作剖分，进而构造基函数或单元形状函数，形成有限元空间，将偏微分方程转化成了有限元方程，利用有效的有限元方程的解法，给出偏微分方程近似解的过程称为有限元法。

李微分方程数值解习题解答1-1如果（P（0）= 0 ,则称X。

是丿（X）的驻点（或稳走点）•矩阵A对称（不必正走）, 求证兀0是丿（兀）的驻点的充要条件是：兀0是方程组= b的解证明:由恥）的走义与内积的性线性性质，得必要性:由0 （0） = 0,得，对于任何X w R",有（Ax Q-b,x） = 0r由线性代数结论知，充分性：由Ax0 = b,对于任何'G R",即兀是丿0）的驻点.§1-2补充：证明广（劝的不同的广义导数几乎处处相等. 证明:设f e L2（Z）,碍e E⑴为/（x）的广义导数，由广义导数的定义可知，对于任意0（劝《：（/）,有两式相减，得到由变分基本引理胡-&几乎处处为零间g理几乎处处相等•补充证明血巧的连续性条件（证明：设I p(x) \<M.\q{x) \<M ,^Sch warz不等式I tz(w,v) 1=1 ^\puv + quv)dx \< M II u II .11 v II II w II .11 v II< 2M*II u II, .11 v II*其中W = max{ M,M }习题:1设fg为/g的一阶广义导数,试用类似的方法走义f⑴的k阶导数伙=1,2,...）解:一阶广义导数的定义，主要是从经典导数经过分部积分得到的关系式来走义，因此可得到如下走义:对于/⑴e厶S若有g（x） e皿）,使得对于任意的^eC；（Z）f有则称/（X）有k阶广义导数,g（x）称为/（X）的k阶广义导数，并记gd）卑dx注:高阶广义导数不是通过递推走义的，可能有高阶导数而没有低阶导数.2.利用必）的完全性证明是Hilbert空间• 证明:只证刃⑴的完全性•设｛/“｝为刃⑴的基本列间因此知｛九｝,｛/；｝都是厶（/）中的基本列（按厂（/）的范数）•由从/）的完全性，存在f,g € %）,使II f n -门0,11 f n - g ll°T 0,以下证明皿-门TO（关键证明防孚）ax由Schwarz不等式f有对于任意的0(兀)e C；(/)f成立由 f fn (x)(p(x)dx = -£ f n (x)(p (x)dx 取极限得到［g(x)0(x)Jx = -^f(x)(p(x)dx 即g(x) = f \即仔H'(D且故H (/)中的基本列是收敛的,H (/)是完全的.3•证明非齐次两点边值问题证明:边界条件齐次化令u0(x) = a + p(x — d)f则w = u—如满足齐次边界条件亠满足的方程为5 =Lu — Lu(} = f - Lu{} f即w对应的边值问题为f Lw = f-Lu0⑹[w(d) = 0, vv (/?) = 0由走理知，问题p与下列变分问题等价求哄cm；, J(w.) = rninJ\w) 其中丿(w) = |a(w,w)-(y- Lu0,w)・而而(Lu Q,u)-a(u Q9u) = -p(b)j3u(b) + C2从而J*(w) = 7(w)- p(b)陋b) + C则关于w的变分问题p等价于：求仏 e C2C\H l,u(a) = a使得其中丿(u) = - p(b)陋b)4就边值问题(解:令％ =& + J3(x-a)t w = u- u()t则"满足等价于:Vve/7；应用分部积分，还原"，于是，边值问题等价于:求"e H',u(a) = 使得VveH；,成立注:形式上与用i去乘方程两端，应用分部积分得到的相同.5试建立与边值问题 等价的变分问题.解:取解函数空间为比(D 对于任意「e 比(/)用v 乘方程两端,应用分部积分，得到rjz ,d 4u 、 z d'u - dhi lh c h u M 而加小【乔皿二乔以T 冇 上式为+ uv ^dx =(/' v)走义a(",y) = f ［丫 : : : +他］必, 为双线性形式. Ja dx- dx^变分问题为:求u e /7j(/)f Vv e &(/)1-41.•用Ritz - Galerkin 方法求边值问题(p t (x) = sin(iTZX ),/ = 1,2,…屮3 -bd u dv 1 2 •——dx ]il dx dx 次近似色(%) f 基解：(1)边界条件齐次化:令如=x,w = u-u0,则W 满足齐次边界条件，且第n次近似W”取为W” = £c,0 , M中 c.(z = l,2,...n)满足的叭 - Galerkin方程为由三角函数的正交性得到2而八®十（_ 25）心丙［―于是得到最后得到2•在题1中專（1）=0代替右边值条件川3是用Ritz - Galerkin方法求解相应问题的第斤次近似"证明色（兀）按厂（0,1）收敛至U叭x），并估计误差• 证明:他对应的级数绝对收敛，由｛血沁｝的完全性知极限就是解畑直误差估计为3•就边值问题（(p.(x) = (x-ay(i = 1,2,•・；)，写出Ritz - Galerkin方程解：边界条件齐次化，取竝=CC p(X _ 6Z)，W = U _ U() 9 W对应的微分方程为对应的变分方程为变分方程为d(w,y) = (/,*) + /3p(b)v(b)— £[/^?v(x) —qu^lx取0(x) = (x-ay9i= 1,2,…，心则Ritz ・GaleTkin方程为取P = 1,纟=0, / = 1，具体计算〃=1，d(0], 0]) = f \dx = (b _ a)£ =丄(b - a)2 + 0(b - a) - 0(b — a) = L (b - a)，，2 2C[ =q(b —%即解络=%0 + —(x — a)n = 2：d(%,0)= (b_a), Q(0,輕)=^2(x-a)dx = (b-a)2得到方程组为特别取a=oe=i z有求解得到* *=-詁=1其解为均=况（）+（X 一Cl）（X 一Q）~Ch2椭圆与抛物型方程有限元法§1.1用线性元求下列边值问题的数值解：此题改为一y + y = 1, y(0) = y(l) = 0/ = 1/4 解：取/i = l/2,x y. = = 0,1,2)』,儿为未知数.Galerkin形式的变分方程为(Lu.v) = (/», 其中龙2(Lu,v) = vdx + — ^uvdx,(/,v) = £2sin —xv{x)dx又一vdx = -uv Ijj ^uvdx =9因此a（u.v） = £（wv -^-—uv）dx在单元/. = [x,,兀]中，应用仿射变换侷部坐标）Eh节点基函数为取/7 = 1 / 2,则计算得G ®, %) = 4 +咅丄厶代数方程组为代如求值・取力= 1/4,未知节点值为坷上2，均，弘4方程为应用局部坐标f表示，系数矩阵为A =他{—4 + £ ,8 +壬,_4 + £} 取/=1, (/ ®)=町:磁+皿(1—§)込占2.就非齐次第三边值条件导出有限元方程.解:设方程为厶“ = -(pu) +qu=f则由变分形式为:V「记I IA(u,v) = (pu ) + (qu°) + a2p(b)u(b)v(b) - a l p(a)u(a)v(a) F(v) = (/» + p(b)/32v(b) — p(a)^v(a)则上述变分形式可表示为4仏巧= F(v)设节点基函数为03(丿=0,1,2, (7V)则有限元方程为具体计算使用标准坐标二。

偏微分方程习题及答案【篇一：偏微分方程数值解法期末考试题答案】题答案及评分标准学年学期：专业：班级：课程：教学大纲：使用教材：教材作者：出版社：数学与应用数学数学偏微分方程数值解法《偏微分方程数值解法》教学大纲（自编，2006）《偏微分方程数值解法》陆金甫、关治清华大学出版社一、判断题（每小题1分，共10分）1、（o）2、（o）3、（x）4、（x）5、（o）6、（o）7、（o）8、（x）9、（x） 10、（o）二、选择题（每小题2分，共10分） 11、（d） 12、（a） 13、（c） 14、（b）15、（c）三、填空题（每小题2分，共20分）?2?216、2?2??x1?x2?2?2 17、a=[4 5 9;23 5 17;11 23 1] 18、y=exp(-t/3)*sin(3*t) ?xn19、help 20、zeros(m,n)21、inva(a)*b或者a/b22、a=sym([cos(x-y),sin(x+y);exp(x-y),(x-1)^3])?(s)?1?(s)?c[??(s)]2?023、a[?2(s)]2?2b?224????v(?)ed? 25、i?xu(xj,tn?1)?u(xj,tn)?四、计算题：（每小题12分，共36分）?u?u?0（x?r,t?0）的有限差分方程（两层显示26、写成对流方程?a?t?x格式，用第n层计算第n+1层），并把有限差分方程改写为便于计算的迭代格式???/h为网格比。

解：在点(xj,tn)处，差分方程为?1un?unjj??anunj?1?ujh?0（j?0,?1,?2,，n?0,1,2,）（8分）便于计算的形式为?1nnn???/h （4分） un?u?a?(u?ujjj?1j)，?u?2u?a2的有限差分方程（中心差分格式，用第n层27、写出扩散方程?t?x计算第n+1层），并把有限差分方程改写为便于计算的迭代格式，???/h2为网格比。

《偏微分方程数值解法》试题（专业：凝聚态物理学号：60 姓名：鄢建军）1.考虑定解问题（1）用迎风格式（）求解1,0(,0)0,0t xu uxu xx+=⎧⎪≤⎧⎨=⎨⎪>⎩⎩。

利用迎风格式编写Fortran程序语言，运行结果如下：Fig 1.迎风格式求解结果（2）用Beam-Warming 格式（）求解。

利用Beam—Warming格式编写Fortran程序语言，运行结果如下：Fig 2. Beam—Warming 格式求解结果（3） 比较两种方法结果的异同。

将两种格式运行的结果绘制在一起，要求时间步长和空间步长在两种格式中都相同，运行结果如下图所示：Fig 3. 迎风格式和Beam-Warming 格式求解结果比较从两种格式的运行结果来看，都存在边缘的误差现象，相比而言，Beam-Warming 格式的运行结果差一些。

但是理论上分析，迎风格式的截断误差为()h οτ+，而Beam-Warming 格式的截断误差为22()h h οττ++。

稳定性上来分析，迎风格式的稳定性较好，要求1(/)a h λλτ≤=，Beam-Warming 格式的稳定性条件为2(/)a h λλτ≤=。

2.考虑定解问题2121110,04(,0)sin ,0(0,)(,)0u u a x l t t u x x x l l u t u l t π⎧∂∂-=<<⎪∂∂⎪⎪=<<⎨⎪⎪==⎪⎩实际计算时，取下列参数：a=1；1l=2.计算进行到合适的时刻为止。

要求：（1） 用加权隐式格式（）求解该问题，研究不同θ值对解的影响。

用加权隐式格式求解该问题，研究不同θ值对解的影响。

采用的差分格式为12212[(1)]0n nj jnn x j x j u u a u u hθδθδτ++---+=， 其截断误差为：2221(),21(),2o h o h τθτθ⎧+≠⎪⎪⎨⎪+=⎪⎩，稳定性条件为：112,(0),1221,(1).2a λθθθ⎧≤≤<⎪⎪-⎨⎪≤≤⎪⎩绝对稳定 当θ取不同值，得出的结果不同。

一、（10分)简叙偏微分方程数值解研究的内容。

建立一个偏微分方程数值格式，需要研究和讨论哪些问题？一个好的数值格式应该到达哪些要求？二、(10分）简叙用差分方法求解偏微分方程数值解的方法和步骤。

三、（10分）请解释收敛性、稳定性、相容逼近性三个概念。

它们之间有何区别和联系？四．对如下两点边值问题：⎪⎩⎪⎨⎧==<<=+-,0)1( )0(,10 ,222u u x f u dx ud 将区间[0,1］作剖分：N N x x x x <<<<=-1100 , （1） 在上述网格剖分下建立中心差分格式; （10分） （2） 写出局部截断误差；并问当网格剖分满足什么条件时，上述差分格式的误差阶为2； （10分）（3） 利用极值定理证明解的存在唯一性。

（10分)四、(1）中心差分格式：)()()(2))()()()((21111u R x f x u h x u x u h x u x u h h i i i ii i i i i i i +=+---+--+++ （2）局部截断误差：)(3)(2331h O dx u d hh u R iii i +⎥⎦⎤⎢⎣⎡--=+ 显然当剖分为等距时,误差阶为2 (3）用极值定理证存在唯一性：只需证明齐次（边值与右端恒零）问题只有平凡解.实际上,设iu 是齐次问题的解，则由极值定理，i u 既不能在区间内部取正的极大,也不能取负的极小，因此0=i u 。

七．已知两点边值问题：⎪⎩⎪⎨⎧==<<=-0)1( )0(10 ,22u u x f dx ud (A) （1） 写出问题（A ）的虚功原理（含证明)； （25分) （2） 写出建立线性有限元方程组h h h b U A =的主要步骤； （10分）（3） 对于线性有限元方程组h h h b U A =，设h A 产生扰动h A ∆,h b 产生扰动h b ∆，由此引起h U 产生扰动h u ∆，说明h u ∆与h A ∆，h b ∆之间的关系.(5分）七． （1）⎪⎩⎪⎨⎧==<<=-0)1( )0(10 ,22u u x f dx u d （A ）检验函数空间为：{}0)1(,0)0(),()()(1==∈=v v I H x v x v V ,]1,0[=I ，本题中检验函数空间与 试探函数空间相同.对V v ∈∀ 有：)(122---⎰dx v f dx ud ⇒0101=-⎰⎰dx v f dx dx vd dx u d dx fv dx dx vd dx u d v u a ⎰⎰==110v)(f, , ),( （2分)从而问题（A ）变分问题可描述为：求V U u =∈使得),(),(v f v u a = 对V v ∈∀ 成立。

二、改进的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 +=+ 校正:)].,(),([2111+++++=n n n n n n y x f y x f hy y(1.15)这个计算公式也可以表示为11(,),(,),1().2p n n nc n n p n p cy y hf x y y y hf x y y y y ++⎧=+⎪⎪=+⎪⎨⎪=+⎪⎪⎩例1 取步长0.1h =，分别用Euler 方法及改进的Euler 方法求解初值问题d (1),01,d (0) 1.yy xy x xy ⎧=-+≤≤⎪⎨⎪=⎩ 解 这个初值问题的准确解为()1(21)xy 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 ，得110.1[1(101)]0.9,10.1[0.9(10.10.9)]0.9019,1(0.90.9019)0.900952p c y y y ⎧=-⨯⨯+⨯=⎪⎪=-⨯⨯+⨯=⎨⎪⎪=+=⎩ 20.900950.1[0.90095(10.10.90095)]0.80274,0.900950.1[0.80274(10.20.80274)]0.80779,1(0.802740.80779)0.805262p c y y y ⎧=-⨯⨯+⨯=⎪⎪=-⨯⨯+⨯=⎨⎪⎪=+=⎩ 这样继续计算下去，其结果列于表9.1.从表9.1可以看出,Euler 方法的计算结果只有2位有效数字,而改进的Euler 方法确有3位有效数字,这表明改进的Euler 方法的精度比Euler 方法高.例2 试用Euler 方法、改进的Euler 方法及四阶经典R-K 方法在不同步长下计算初值问题⎪⎩⎪⎨⎧=≤≤+-=1)0(,10),1(d d y x xy y xy在0.2、0.4、0.8、1.0处的近似值，并比较它们的数值结果.解 对上述三种方法,每执行一步所需计算)1(),(xy y y x f +-=的次数分别为1、2、4。

偏微分方程数值解一（10分）、设矩阵A 对称正定，定义)(),(),(21)(n R x x b x Ax x J ∈-=，证明下列两个问题等价：(1)求n R x ∈0使)(min )(0x J x J n Rx ∈=;(2)求下列方程组的解：b Ax = 解: 设n R x ∈0是)(x J 的最小值点,对于任意的n R x ∈,令),(2),()()()(2000x Ax x b Ax x J x x J λλλλϕ+-+=+=, (3分)因此0=λ是)(λϕ的极小值点,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分) 反之,若nR x ∈0满足bAx =0,则对于任意的x ,)(),(21)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 dxdu 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==∈=b u a u b a H u u H 为求解函数空间,检验函数空间.取),(10b a H v ∈,乘方程两端,积分应用分部积分得到 (3分))().(),(v f fvdx dx quv dxdv dx du p v u a b a ba ==+=⎰⎰,),(1b a H v ∈∀ 即变分问题的Galerkin 形式. (3分)令⎰-+=-=b a dx fu qu dxdup u f u u a u J ])([21),(),(21)(22,则变分问题的Ritz 形式为求),(1*b a H u ∈,使)(m in )(10*u J u J H u ∈= (4分) 评分标准:空间描述与积分步骤3分,变分方程3分,极小函数及其变分问题4分,三（20分）、对于边值问题⎪⎩⎪⎨⎧=⨯=∈-=∂∂+∂∂∂0|)1,0()1,0(),(,12222G u G y x yux u （1）建立该边值问题的五点差分格式（五点棱形格式又称正五点格式），推导截断误差的阶。

偏微分⽅程数值解（试题）偏微分⽅程数值解试题1、考虑⼀维的抛物型⽅程：2200, [0,], 0t T (,), (,)(,0)()x x u ux t xu x t u u x t u u x x ππνπ?==??=∈≤≤??=== （1）导出时间离散是⼀阶向前Euler 格式，空间离散是⼆阶精度的差分格式；（2）讨论（1）中导出的格式的稳定性；（3）若时间离散为⼆阶精度的蛙跳格式，112n n n t t u u u t t+-=?-=空间离散是⼆阶精度的中⼼差分，问所导出的格式稳定吗？为什么？2、考虑Poission ⽅程2(,)1, (,)0, in AB and AD (,)0, in BC and CDu x y x y unu x y -?=∈Ω== 其中Ω是图1中的梯形。

使⽤差分⽅法来离散该⽅程。

由于梯形的对称性，可以考虑梯形的⼀半，如图2，图2 从物理空间到计算区域的⼏何变换图1 梯形为了求解本问题，采⽤如下⽅法：将Ω的⼀半投影到正⽅形区域?Ω，然后在?Ω上使⽤差分⽅法来离散该⽅程。

在计算区域?Ω上⽤N N ?个⽹格点，空间步长为1/(1)N ξη?=?=-。

（1）引⼊⼀个映射T 将原区域Ω（带有坐标,x y ）变换到单位正⽅形?Ω（带有坐标,ξη）。

同时导出在新区域上的⽅程和边界条件。

（2）在变换区域，使⽤泰勒展开导出各导数项在区域内部和边界点上的差分格式。

3、对线性对流⽅程0 constant >0u u a a t x+=，其⼀阶迎风有限体积法离散格式为 1?n j u +=?nj u a t x-（?n j u 1?n j u --）（1）写出0a <时的⼀阶迎风有限体积法的离散格式；（2）写出a 为任意符号的常数的⼀阶迎风有限体积法的守恒形式。

（3）使⽤0 u uu t x+=说明⼀阶迎风有限体积法不是熵保持的格式。

4、对⼀维Poission ⽅程, (0,1)(0)(1)0x xx u xe x u u ?-=∈?==? 将[]01，分成(1)n +等分，写出⽤中⼼差分离散上述⽅程的差分格式，并问：（1）该差分格式与原微分⽅程相容吗？为什么？（2）该差分格式稳定吗？为什么？（3）该差分格式是否收敛到原微分⽅程的解？为什么？（4）取(1)6n +=，写出该差分格式的矩阵表⽰。