偏微分方程数值解复习题(2011硕士)

合集下载

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

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

偏微分方程数值解试题(06B )参考答案与评分标准信息与计算科学专业一(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 的最小值点. (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 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分,变分方程3分,极小函数及其变分问题4分,三(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 (1)建立该边值问题的五点差分格式(五点棱形格式又称正五点格式),推导截断误差的阶。

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

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

一.填空(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)建立该边值问题的五点差分格式(五点棱形格式又称正五点格式),推导截断误差的阶。

数学偏微分方程的数值解法

数学偏微分方程的数值解法

数学偏微分方程的数值解法当然可以。

这里是根据标题“数学偏微分方程的数值解法”出的20道试题,包括选择题和填空题,每道题目都有详细的序号介绍:1. 选择题:偏微分方程的数值解法主要适用于哪类方程?A. 常微分方程B. 偏微分方程C. 代数方程D. 差分方程2.填空题:数值解法中常用的一种基础方法是______________。

3.选择题:有限差分法是一种适用于哪种类型的偏微分方程的数值解法?A. 椭圆型B. 抛物型C. 双曲型D. 超越型4.填空题:偏微分方程的数值解法通常要求将空间区域离散化为___ ___________。

5. 选择题:哪种方法适合处理偏微分方程的初始值问题?A. 有限元法B. 有限差分法C. 傅里叶变换法D. 辛普森法则6.填空题:数值解法中,常用的稳定性分析方法包括_____________ _。

7.选择题:对于偏微分方程的边值问题,常用的数值方法是_______ _______。

A. 有限体积法B. 辛普森法则C. 椭圆积分法D. 有限元法8.填空题:描述一种常见的数值解法的收敛性条件______________。

9.选择题:哪种方法在处理时间依赖性偏微分方程时特别有效?A. 隐式方法B. 显式方法C. 中心差分法D. 前进差分法10.填空题:数值解法中的矩阵求解通常利用______________方法。

11.选择题:在有限元法中,通常要对空间区域进行如何划分?A. 正交分解B. 三角剖分C. 曲面划分D. 直角分割12.填空题:有限差分法的精度通常与______________相关。

13.选择题:哪种方法适合处理非线性偏微分方程的数值解?A. 变分法B. 有限元法C. 辛普森法则D. 显式方法14.填空题:对于稳定性的要求,常用的数值方法需要满足_________ _____条件。

15.选择题:哪种方法可以有效地处理多维偏微分方程的数值解?A. 傅里叶变换法B. 辛普森法则C. 多重网格法D. 变分法16.填空题:在求解偏微分方程数值解时,通常需要考虑___________ ___问题。

偏微分方程数值习题解答

偏微分方程数值习题解答

偏微分⽅程数值习题解答李微分⽅程数值解习题解答 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 ==∑=具体计算使⽤标准坐标ξ.。

偏微分方程数值分析练习题

偏微分方程数值分析练习题

一、判断题(判断下列题目是否正确,如果正确请打“√”,错误请打“×”)1、Possion 方程属于椭圆型方程。

( √ )2、所有的对角矩阵、Hermite 矩阵和酉矩阵都是正规矩阵。

( √ )3、一阶双曲型方程并不是每个边界都给定条件,因此差分方程所需的边界条件往往要比微分方程的边界条件要多。

( )4、设),(t x u 是定解问题充分光滑的解,对于其差分格式,如果当0,→h τ时,其截断误差趋向于0,则称差分与定解问题是稳定的。

( x )5、Von Neumann 条件只是稳定性的必要条件。

( x )6、 如果当时间步长τ和空间步长h 趋于0时,差分格式的截断误差趋于0,就说这个差分格式是与微分方程相容的。

( √ )7、 如果当时间步长τ和空间步长h 无限缩小时,差分格式的解能无限逼近到微分方程的解,就说这个差分格式是收敛的。

( √ )8、 差分格式的相容性和收敛性是两个完全不同的概念,但是相容的差分格式必定也是收敛的。

( x )9、 差分格式的稳定性与差分格式本身有关,而与网格比的大小无关。

( x )10、 给定一个适定的线性初值问题,如果逼近它的差分格式是和它相容的,那么差分格式的收敛性是差分格式的稳定性的充分且必要条件。

( √ )11、 Von Neumann 条件只是稳定性的必要条件。

( x )12、 如果我们得到差分格式的截断误差为)()/(23h O h O +τ,那么差分格式在时间上是三阶精度的,在空间上是二阶精度的。

( x )13、 在有限差分数值解法中,无法得到一个无条件相容并且无条件稳定的显示差分格式。

( √ )二、填空题1.如果一个定解问题的解的存在性、唯一性和稳定性成立,称定解问题是 适定 的。

2.如果n n C ⨯∈u ,满足 U H U=I ,则u 称为酉矩阵。

3.Lax 等价定理是: 给一个适定的线性初值问题以及与其相容的差分格式,则差分格式的稳定性是差分格式收敛性的充分必要条件4.迎风格式的基本思想是:双曲线方程中关于空间偏导数用在特征线方向一侧的单边差商来代替。

偏微分方程考试题

偏微分方程考试题

数学物理方程及数值解 复习提要一、偏微分方程的建立 CH1 典型方程和定解条件 【内容提要】1. 方程的建立(步骤:确定物理量;微元法建立等式;化简得方程)主要方法:微元法; 泛定方程:(1) 波动方程(双曲型):弦振动方程:222222(,)(,)(),()u x t u x t F a a txρ∂∂==∂∂张力单位长度弦质量 传输线方程:222222222221,00i a LCi a a t x t x νν∂∂∂∂-=-=∂∂∂=∂;, 电磁场方程:22222211,,H E H E t t εμεμ∂∂=∇=∇∂∂22222222221(),με标量函数形式:∂∂∂∂=++∂∂∂=∂u u u z a u a t x y (2) 热传导方程/扩散方程(抛物型):ρ,其中22u Fa u f f t c ∂=∇+=∂ 导热杆(无热源)222u u a t x ∂∂=∂∂, 导热片(无热源)22222()u u u a t x y ∂∂∂=+∂∂∂ (3) 稳恒方程(椭圆型):Poisson 方程:,2u f ∇= Laplace 方程:,20u ∇=2.定解条件:初始条件及边界条件边界条件(1)第一类边界条件(Dirichlet 条件): 1(,)(,)D u M t f M t ∂=(2) 第二类边界条件(Neumann 条件):2Duf n ∂∂=∂ (3) 第三类边界条件(Robin 条件): 3()Duu f n σ∂∂+=∂ 3.定解问题的提法:⎧⎪⎧⎨⎨⎪⎩⎩偏微分方程(泛定方程)定解问题初始条件定解条件边界条件()Cauchy ⎧⎨⎩泛定方程(1)初始问题初始条件 ⎧⎨⎩泛定方程(2)边界问题(第一,二,三)边界条件⎧⎪⎨⎪⎩泛定方程(3)混合问题初始条件边界条件4.线性偏微分方程的基本性质(1).线性迭加原理212,11,,,,,,,:nnij i ij i n i j i i j iL a b c a b c f x x x x x x ==∂∂=++∂∂∂∑∑其中是算子的函数111(1,2)(),nnni i ii ii i i i i i L u f in L c u c L u c f=====⇒==∑∑∑命题:21110(1,2),,()0,nnii i i i i i i i i i k j u Lu i c u c L c u x x ∞===∂==⇒=∂∂∑∑∑一致敛命收题:(2.) 齐次化原理(冲量原理)Duhamel 原理:设(,,)x t ωτ是方程22222,,(,)(,)0,(,),a x t t x x x f x x t ωτωτωττω⎧∂∂=-∞<<+∞>⎪∂∂⎪⎨∂⎪==-∞<<+∞⎪∂⎩的解,⇒0(,,)d ,()t x t u x t ωττ=⎰是方程22222(,),,0(,0)(,0)0,0,u u a f x t x t tx u x u x x t ⎧∂∂=+-∞<<+∞>⎪∂∂⎪⎨∂⎪==-∞<<+∞⎪∂⎩的解。

(完整word版)偏微分方程数值解法答案

(完整word版)偏微分方程数值解法答案

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.peEllip5 用五点差分格式解拉普拉斯方程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:nyu0(j,1) = EllIni2Uxl(minx,miny+(j-1)*hy);u0(j,nx) = EllIni2Uxr(maxx,miny+(j-1)*hy);endfor j=1:nxu0(1,j) = EllIni2Uyl(minx+(j-1)*hx,miny);u0(ny,j) = EllIni2Uyr(minx+(j-1)*hx,maxy);endA = -4*eye((nx-2)*(ny-2),(nx-2)*(ny-2));b = zeros((nx-2)*(ny-2),1);for i=1:(nx-2)*(ny-2)if mod(i,nx-2) == 1if i==1A(1,2) = 1;A(1,nx-1) = 1;b(1) = - u0(1,2) - u0(2,1);elseif i == (ny-3)*(nx-2)+1A(i,i+1) = 1;A(i,i-nx+2) = 1;b(i) = - u0(ny-1,1) - u0(ny,2);elseA(i,i+1) = 1;A(i,i-nx+2) = 1;A(i,i+nx-2) = 1;b(i) = - u0(floor(i/(nx-2))+2,1);endendelseif mod(i,nx-2) == 0if i == nx-2A(i,i-1) = 1;A(i,i+nx-2) = 1;b(i) = - u0(1,nx-1) - u0(2,nx);elseif i == (ny-2)*(nx-2)A(i,i-1) = 1;A(i,i-nx+2) = 1;b(i) = - u0(ny-1,nx) - u0(ny,nx-1);elseA(i,i-1) = 1;A(i,i-nx+2) = 1;A(i,i+nx-2) = 1;b(i) = - u0(floor(i/(nx-2))+1,nx);endendelseif i>1 && i< nx-2A(i,i-1) = 1;A(i,i+nx-2) = 1;A(i,i+1) = 1;b(i) = - u0(1,i+1);elseif i > (ny-3)*(nx-2) && i < (ny-2)*(nx-2)A(i,i-1) = 1;A(i,i-nx+2) = 1;A(i,i+1) = 1;b(i) = - u0(ny,mod(i,(nx-2))+1);elseA(i,i-1) = 1;A(i,i+1) = 1;A(i,i+nx-2) = 1;A(i,i-nx+2) = 1;endendendendendul = A\b;for i=1:(ny-2)for j=1:(nx-2)u(i,j) = ul((i-1)*(nx-2)+j);endendformat short;2.peEllip5m 用工字型差分格式解拉普拉斯方程function u = peEllip5m(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:nyu0(j,1) = EllIni2Uxl(minx,miny+(j-1)*hy);u0(j,nx) = EllIni2Uxr(maxx,miny+(j-1)*hy);endfor j=1:nxu0(1,j) = EllIni2Uyl(minx+(j-1)*hx,miny);u0(ny,j) = EllIni2Uyr(minx+(j-1)*hx,maxy);endA = -4*eye((nx-2)*(ny-2),(nx-2)*(ny-2));b = zeros((nx-2)*(ny-2),1);for i=1:(nx-2)*(ny-2)if mod(i,nx-2) == 1if i==1A(1,nx) = 1;b(1) = - u0(1,1) - u0(3,1) - u0(1,3);elseif i == (ny-3)*(nx-2)+1A(i,i-nx+1) = 1;b(i) = - u0(ny,1) - u0(ny,3) - u0(ny-2,1);elseA(i,i-nx+3) = 1;A(i,i+nx-1) = 1;b(i) = - u0(floor(i/(nx-2))+1,1) - u0(floor(i/(nx-2))+3,1);endendelseif mod(i,nx-2) == 0if i == nx-2A(i,i+nx-1) = 1;b(i) = - u0(1,nx-2) - u0(1,nx) - u0(3,nx);elseif i == (ny-2)*(nx-2)A(i,i-nx+1) = 1;b(i) = - u0(ny,nx) - u0(ny,nx-2) - u0(ny-2,nx);elseA(i,i-nx+1) = 1;A(i,i+nx-3) = 1;b(i) = - u0(floor(i/(nx-2)),nx) - u0(floor(i/(nx-2))+2,nx);endendelseif i>1 && i< nx-2A(i,i+nx-1) = 1;A(i,i+nx-3) = 1;b(i) = - u0(1,i) - u0(1,i+2);elseif i > (ny-3)*(nx-2) && i < (ny-2)*(nx-2)A(i,i-nx+3) = 1;A(i,i-nx+1) = 1;b(i) = - u0(ny,mod(i,(nx-2))) -u0(ny,mod(i,(nx-2))+2);elseA(i,i-nx+3) = 1;A(i,i+nx-1) = 1;A(i,i+nx-3) = 1;A(i,i-nx+1) = 1;endendendendendul = A\b;for i=1:(ny-2)for j=1:(nx-2)u(i,j) = ul((i-1)*(nx-2)+j);endendformat short;3.peHypbYF 用迎风格式解对流方程function u = peYF(a,dt,n,minx,maxx,M)format long;h = (maxx-minx)/(n-1);if a>0for j=1:(n+M)u0(j) = IniU(minx+(j-M-1)*h);endelsefor j=1:(n+M)u0(j) = IniU(minx+(j-1)*h);endendu1 = u0;for k=1:Mif a>0for i=(k+1):n+Mu1(i) = -dt*a*(u0(i)-u0(i-1))/h+u0(i);endelsefor i=1:n+M-ku1(i) = -dt*a*(u0(i+1)-u0(i))/h+u0(i);endendu0 = u1;endif a>0u = u1((M+1):M+n);elseu = u1(1:n);endformat long;4.peHypbLax 用拉克斯-弗里德里希斯格式解对流方程function u = peHypbLax(a,dt,n,minx,maxx,M)format long;h = (maxx-minx)/(n-1);for j=1:(n+2*M)u0(j) = IniU(minx+(j-M-1)*h);endu1 = u0;for k=1:Mfor i=k+1:n+2*M-ku1(i) = -dt*a*(u0(i+1)-u0(i-1))/h/2+(u0(i+1)+u0(i-1))/2;endu0 = u1;endu = u1((M+1):(M+n));format short;4.peHypbLaxW 用拉克斯-温德洛夫格式解对流方程function u = peLaxW(a,dt,n,minx,maxx,M)format long;h = (maxx-minx)/(n-1);for j=1:(n+2*M)u0(j) = IniU(minx+(j-M-1)*h);endu1 = u0;for k=1:Mfor i=k+1:n+2*M-ku1(i) = dt*dt*a*a*(u0(i+1)-2*u0(i)+u0(i-1))/2/h/h - ...dt*a*(u0(i+1)-u0(i-1))/h/2+u0(i);endu0 = u1;endu = u1((M+1):(M+n));format short;6.peHypbBW 用比姆-沃明格式解对流方程function u = peBW(a,dt,n,minx,maxx,M)format long;h = (maxx-minx)/(n-1);for j=1:(n+2*M)u0(j) = IniU(minx+(j-2*M-1)*h);endu1 = u0;for k=1:Mfor i=2*k+1:n+2*Mu1(i) = u0(i)-dt*a*(u0(i)-u0(i-1))/h-a*dt*(1-a*dt/h)* ...(u0(i)-2*u0(i-1)+u0(i-2))/2/h;endu0 = u1;endu = u1((2*M+1):(2*M+n));format short;7.peHypbRich 用Richtmyer多步格式解对流方程function u = peRich(a,dt,n,minx,maxx,M)format long;h = (maxx-minx)/(n-1);for j=1:(n+4*M)u0(j) = IniU(minx+(j-2*M-1)*h);endfor k=1:Mfor i=2*k+1:n+4*M-2*ktmpU1 = -dt*a*(u0(i+2)-u0(i))/h/4+(u0(i+2)+u0(i))/2;tmpU2 = -dt*a*(u0(i)-u0(i-2))/h/4+(u0(i)+u0(i-2))/2;u1(i) = -dt*a*(tmpU1-tmpU2)/h/2+u0(i);endu0 = u1;endu = u1((2*M+1):(2*M+n));format short;8.peHypbMLW 用拉克斯-温德洛夫多步格式解对流方程function u = peMLW(a,dt,n,minx,maxx,M)format long;h = (maxx-minx)/(n-1);for j=1:(n+2*M)u0(j) = IniU(minx+(j-M-1)*h);endu1 = u0;for k=1:Mfor i=k+1:n+2*M-ktmpU1 = -dt*a*(u0(i+1)-u0(i))/h/2+(u0(i+1)+u0(i))/2;tmpU2 = -dt*a*(u0(i)-u0(i-1))/h/2+(u0(i)+u0(i-1))/2;u1(i) = -dt*a*(tmpU1-tmpU2)/h+u0(i);endu0 = u1;endu = u1((M+1):(M+n));format short;9.peHypbMC 用MacCormack多步格式解对流方程function u = peMC(a,dt,n,minx,maxx,M)format long;h = (maxx-minx)/(n-1);for j=1:(n+2*M)u0(j) = IniU(minx+(j-M-1)*h);endfor k=1:Mfor i=k+1:n+2*M-ktmpU1 = -dt*a*(u0(i+1)-u0(i))/h+u0(i);tmpU2 = -dt*a*(u0(i)-u0(i-1))/h+u0(i-1);u1(i) = -dt*a*(tmpU1-tmpU2)/h/2+(u0(i)+tmpU1)/2;endu0 = u1;endu = u1((M+1):(M+n));format short;10.peHypb2LF 用拉克斯-弗里德里希斯格式解二维对流方程的初值问题function u = pe2LF(a,b,dt,nx,minx,maxx,ny,miny,maxy,M)%啦-佛format long;hx = (maxx-minx)/(nx-1);hy = (maxy-miny)/(ny-1);for i=1:nx+2*Mfor j=1:(ny+2*M)u0(i,j) = Ini2U(minx+(i-M-1)*hx,miny+(j-M-1)*hy);endendu1 = u0;for k=1:Mfor i=k+1:nx+2*M-kfor j=k+1:ny+2*M-ku1(i,j) = (u0(i+1,j)+u0(i-1,j)+u0(i,j+1)+u0(i,j-1))/4 ...-a*dt*(u0(i+1,j)-u0(i-1,j))/2/hx ...-b*dt*(u0(i,j+1)-u0(i,j-1))/2/hy;endendu0 = u1;endu = u1((M+1):(M+nx),(M+1):(M+ny));format short;11.peHypb2FL 用拉克斯-弗里德里希斯格式解二维对流方程的初值问题function u = pe2FL(a,b,dt,nx,minx,maxx,ny,miny,maxy,M)format long;hx = (maxx-minx)/(nx-1);hy = (maxy-miny)/(ny-1);for i=1:nx+4*Mfor j=1:(ny+4*M)u0(i,j) = Ini2U(minx+(i-2*M-1)*hx,miny+(j-2*M-1)*hy);endendu1 = u0;for k=1:Mfor i=2*k+1:nx+4*M-2*kfor j=2*k-1:ny+4*M-2*k+2tmpU(i,j) = u0(i,j) - a*dt*(u0(i+1,j)-u0(i-1,j))/2/hx + ...(a*dt/hx)^2*(u0(i+1,j)-2*u0(i,j)+u0(i-1,j))/2;endendfor i=2*k+1:nx+4*M-2*kfor j=2*k+1:nx+4*M-2*ku1(i,j) = tmpU(i,j) - b*dt*(tmpU(i,j+1)-tmpU(i,j-1))/2/hy + ...(b*dt/hy)^2*(tmpU(i,j+1)-2*tmpU(i,j)+tmpU(i,j-1))/2;endendu0 = u1;endu = u1((2*M+1):(2*M+nx),(2*M+1):(2*M+ny));format short;12.peParabExp 用显式格式解扩散方程的初值问题function u = peParabExp(c,dt,n,minx,maxx,M)format long;h = (maxx-minx)/(n-1);for j=1:(n+2*M)u0(j) = PrIniU(minx+(j-M-1)*h);endu1 = u0;for k=1:Mfor i=k+1:n+2*M-ku1(i) = dt*c*(u0(i+1)-2*u0(i)+u0(i-1))/h/h+u0(i);endendu = u1((M+1):(M+n));format short;13.peParabTD 用跳点格式解扩散方程的初值问题function u = peParabTD(c,dt,n,minx,maxx,M)%跳点format long;h = (maxx-minx)/(n-1);for j=1:(n+2*M)u0(j) = PrIniU(minx+(j-M-1)*h);endu1 = u0;for k=1:Mfor i=k+1:n+2*M-kif mod(n+i,2) == 0u1(i) = u0(i) + c*dt*(u0(i+1) - 2*u0(i) + u0(i-1))/h/h;if i > 2u1(i-1) =(u0(i-1) + c*dt*(u1(i) + u1(i-2))/h/h)/(1 + 2*c*dt/h/h);endendendu0 = u1;endu = u1((M+1):(M+n));format short;14.peParabImp 用隐式格式解扩散方程的初边值问题function u = peParabImp(c,dt,n,minx,maxx,lbu,rbu,M)format long;h = (maxx-minx)/(n-1);u0(1) = lbu;u0(n) = rbu;for j=2:n-1u0(j) = PrIniU(minx+(j-1)*h);endu1 = u0;for k=1:MA = zeros(n-2,n-2);cb = - transpose(u0(2:(n-1)));cb(1) = cb(1) - dt*c*lbu/h/h;cb(n-2) = cb(n-2) - dt*c*rbu/h/h;A(1,1) = -2*dt*c/h/h -1;A(1,2) = dt*c/h/h ;for i=2:n-3A(i,i-1) = dt*c/h/h ;A(i,i) = - 2*dt*c/h/h -1 ;A(i,i+1) = dt*c/h/h ;endA(n-2,n-2) = -2*dt*c/h/h -1;A(n-2,n-3) = dt*c/h/h;u1(2:(n-1)) = A\cb;u0 = u1;endu = u1;format short;15.peParabKN 用克拉克-尼科尔森格式解扩散方程的初边值问题function u = peParabKN(c,dt,n,minx,maxx,lbu,rbu,M)format long;h = (maxx-minx)/(n-1);u0(1) = lbu;u0(n) = rbu;for j=2:n-1u0(j) = PrIniU(minx+(j-1)*h);endu1 = u0;for k=1:MA = zeros(n-2,n-2);cb = zeros(n-2,1);cb(1) = -u0(2) -(u0(3)-2*u0(2)+u0(1))*dt*c*lbu/h/h/2 - dt*c*lbu/h/h/2;cb(n-2) = -u0(n-1) -(u0(n)-2*u0(n-1)+u0(n-2))*dt*c*lbu/h/h/2 - dt*c*rbu/h/h/2;for i=2:n-3cb(i) = -u0(i+1) -(u0(i+2)-2*u0(i+1)+u0(i))*dt*c*lbu/h/h/2;endA(1,1) = -dt*c/h/h -1;A(1,2) = dt*c/h/h/2 ;for i=2:n-3A(i,i-1) = dt*c/h/h/2 ;A(i,i) = - dt*c/h/h -1 ;A(i,i+1) = dt*c/h/h/2 ;endA(n-2,n-2) = -dt*c/h/h -1;A(n-2,n-3) = dt*c/h/h/2;u1(2:(n-1)) = A\cb;u0 = u1;endu = u1;format short;16.peParabWegImp 用加权隐式格式解扩散方程的初边值问题function u = peParabWegImp(c,sita,dt,n,minx,maxx,lbu,rbu,M)format long;h = (maxx-minx)/(n-1);u0(1) = lbu;u0(n) = rbu;for j=2:n-1u0(j) = PrIniU(minx+(j-1)*h);endu1 = u0;for k=1:MA = zeros(n-2,n-2);cb = zeros(n-2,1);cb(1) = -u0(2) -(1 - sita)*(u0(3)-2*u0(2)+u0(1))*dt*c*lbu/h/h/2 ...- sita*dt*c*lbu/h/h;cb(n-2) = -u0(n-1) -(1 - sita)*(u0(n)-2*u0(n-1)+u0(n-2))*dt*c*lbu/h/h/2 ...- sita*dt*c*rbu/h/h;for i=2:n-3cb(i) = -u0(i+1) -(1 - sita)*(u0(i+2)-2*u0(i+1)+u0(i))*dt*c*lbu/h/h;endA(1,1) = -2*sita*dt*c/h/h -1;A(1,2) = sita*dt*c/h/h ;for i=2:n-3A(i,i-1) = sita*dt*c/h/h ;A(i,i) = - 2*sita*dt*c/h/h -1 ;A(i,i+1) = sita*dt*c/h/h ;endA(n-2,n-2) = - 2*sita*dt*c/h/h -1;A(n-2,n-3) = sita*dt*c/h/h;u1(2:(n-1)) = A\cb;u0 = u1;endu = u1;format short;17.peDKExp 用指数型格式解对流扩散方程的初值问题function u = peDKExp(a,b,dt,n,minx,maxx,M)format long;h = (maxx-minx)/(n-1);for j=1:(n+2*M)u0(j) = DKIniU(minx+(j-M-1)*h);endu1 = u0;coff = (exp(a*h/2/b)+exp(-a*h/2/b))/(exp(a*h/2/b)-exp(-a*h/2/b));coff = dt*coff*a*h/2;for k=1:Mfor i=k+1:n+2*M-ku1(i) = coff*(u0(i+1)-2*u0(i)+u0(i-1))+a*dt*(u0(i+1)-u0(i-1))/h/2+u0(i);endu0 = u1;endu = u1((M+1):(M+n));format short;18.peDKSam 用萨马尔斯基格式解对流扩散方程的初值问题function u = peDKSam(a,b,dt,n,minx,maxx,M)format long;h = (maxx-minx)/(n-1);for j=1:(n+2*M)u0(j) = DKIniU(minx+(j-M-1)*h);endu1 = u0;coff = dt*b/(1+a*h/b/2);for k=1:Mfor i=k+1:n+2*M-ku1(i) = coff*(u0(i+1)-2*u0(i)+u0(i-1))+a*dt*(u0(i)-u0(i-1))/h+u0(i);endu0 = u1;endu = u1((M+1):(M+n));format short;。

偏微分方程数值解复习题(2013硕士)

偏微分方程数值解复习题(2013硕士)

偏微分方程数值解期末复习(2012硕士)一、考题类型本次试卷共六道题目,题型及其所占比例分别为:填空题20%;计算题80% 二、按章节复习内容第一章知识点:Euler法、向前差商、向后差商、中心差商、局部截断误差、整体截断误差、相容性、收敛性、阶、稳定性、显格式、隐格式、线性多步法、第一特征多项式、第二特征多项式、稳定多项式、绝对稳定等;要求: 熟练一元函数的数值微分公式;会辨认差分格式, 计算线性多步法的局部截断误差和阶;第二章知识点:矩形网格、(正则,非正则)内点、边界点、偏向前(向后,中心)差商、五点差分格式、增设虚点法、积分插值法、线性椭圆型差分格式、极值原理、比较定理、五点差分格式的相容、收敛和稳定性等;要求: 熟练多元函数的数值微分公式;会建立椭圆型方程边值问题的差分格式;计算局部截断误差;了解极值原理讨论格式的收敛性和稳定性;第四章知识点:最简显格式、最简隐格式、CN格式、双层加权格式、Richardson 格式、网格比、传播因子法(分离变量法) 、传播因子、传播矩阵、谱半径、von Neumann条件、跳点格式、ADI格式等;要求: 会建立抛物型方程边值问题的经典差分格式;计算局部截断误差;会计算格式的传播因子或传播矩阵;会讨论格式的稳定性;第五章知识点:依赖区域、左偏心格式、右偏心格式、中心格式、LF格式、LW 格式、Wendroff格式、跳蛙格式、特征线、CFL条件等;要求: 建立双曲型方程边值问题的差分格式;计算局部截断误差;会计算格式的传播因子或传播矩阵;讨论格式的稳定性;第七章知识点:单元、线性元、线性基、(单元)刚度矩阵、(单元)荷载向量等;要求: 会用线性元(线性基)建立常微分方程边值问题的有限元格式三 练习题1、 试建立Euler 法(向后Euler 法或梯形法),并讨论格式的局部截断误差和阶。

P4,6,8+课件2、 已知一个线性二或三步法,试讨论格式的局部截断误差和阶。

偏微分方程数值解试卷

偏微分方程数值解试卷

一、(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 ∈∀ 成立。

最新偏微分方程数值解试题参考答案

最新偏微分方程数值解试题参考答案

偏微分方程数值解一(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. 什么是偏微分方程的数值解?- A. 用数值方法求解偏微分方程的方法。

- B. 用解析方法求解偏微分方程的方法。

- C. 用近似方法求解偏微分方程的方法。

- D. 用几何方法求解偏微分方程的方法。

2. 数值解法在解偏微分方程中的作用是什么?- 空格填写:______________。

3. 常见的偏微分方程数值解方法包括哪些?- A. 有限差分法。

- B. 蒙特卡罗法。

- C. 拉格朗日法。

- D. 傅里叶法。

4. 有限差分法(FDM)适用于解决什么类型的偏微分方程?- A. 抛物型。

- B. 椭圆型。

- C. 双曲型。

- D. 超越型。

5. 解偏微分方程的数值方法通常需要进行什么样的数值计算?- 空格填写:______________。

6. 数值解方法的稳定性和收敛性为解偏微分方程的数值解带来什么保证?- 空格填写:______________。

7. 在有限元方法(FEM)中,什么是“元”?- 空格填写:______________。

8. 数值解法中常用的时间步长(time step)是用来控制什么的?- A. 精度。

- B. 稳定性。

- C. 收敛性。

- D. 平滑性。

9. 请列举两种常见的偏微分方程数值解法。

10. 数值解法中常用的空间离散化方法有哪些?11. 解偏微分方程的数值方法中,误差分析如何进行?12. 请简要描述有限元方法(FEM)在解偏微分方程中的应用。

13. 有限差分法和有限元法在解决偏微分方程中有何异同?14. 数值解法如何处理非线性偏微分方程?15. 为什么说数值解法对大规模偏微分方程问题有优势?16. 数值解法的计算效率受哪些因素影响?17. 有限体积法在解偏微分方程中有何特点?18. 数值解法中如何选择合适的网格?19. 如何评估数值解法的稳定性?20. 解偏微分方程的数值方法中,收敛性是如何定义和验证的?希望这些题目能够满足你的需求!。

西北工业大学2012级研究生课程《偏微分方程数值解法》考试试题

西北工业大学2012级研究生课程《偏微分方程数值解法》考试试题

2u 2u 0 x 0.3 , 0 y 0.2 2 2 2x y y x u x 0 u x 0.3 2 y u u 1 , u x y 0 n y 0.2
用 h 0.1 作正方形网格剖分 。 (1)用五点菱形差分格式将微分方程在内点离散化; (2)用两种截断误差为 O( h 2 ) 的差分法将第三边界条件离散化; (3)整理后的差分方程组为
2 2u 2 u 四. (12 分)建立波动方程 2 a 的初值问题的显格式,推导截断误差. t x 2
五. (12 分)讨论求解 ut uxx +u yy cu, c const. 的古典显格式的稳定性。 六.对微分方程边值问题
du d ( p ( x ) ) q( x)u f ( x), x (a, b) dx dx du (a ) 1u (a) d1 , 1 0, 1 0 1 dx du (b) 2 dx 2u (b) d 2 , 2 0, 2 0
* 3.将有限元解 y h 用Fra bibliotek函数的形式表示出来:
八. (12 分)设有常微分方程边值问题
y y x 2 - 3, 0 x 1 y (0) 2, y (1) 1
1. 转化为相应的变分问题; 2. 将[ 0 , 1 ]二等分,采用线性元的有限元方法,导出有限元方程并求解。
U A U B U C U D 三. (12)给定初值问题
A C
B D
u u t x

ux,0 x
取时间步长 0.1 ,空间步长 h 0.2 。试合理选用一阶偏心差分格式(最简显格 式), 并以此格式求出解函数 u( x , t ) 在 x 0.2 , t 0.3 处的近似值。 1.所选用的差分格式是: 2.计算所求近似值:

偏微分方程数值解(试题)

偏微分方程数值解(试题)

偏微分⽅程数值解(试题)偏微分⽅程数值解试题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 +=,写出该差分格式的矩阵表⽰。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

偏微分方程数值解期末复习(2011硕士)
一、考题类型
本次试卷共六道题目,题型及其所占比例分别为:
填空题20%;计算题80%
二、按章节复习内容
第一章
知识点:Euler法、向前差商、向后差商、中心差商、局部截断误差、整体截断误差、相容性、收敛性、阶、稳定性、显格式、隐格式、线性多步法、第一特征多项式、第二特征多项式、稳定多项式、绝对稳定等;
要求: 会辨认差分格式, 判断线性多步法的误差和阶;
第二章
知识点:矩形网格、(正则,非正则)内点、边界点、偏向前(向后,中心)差商、五点差分格式、增设虚点法、积分插值法、线性椭圆型差分格式、极值原理、比较定理、五点差分格式的相容收敛和、稳定性等;
要求: 建立椭圆型方程边值问题的差分格式, 极值原理;
第四章
知识点:最简显格式、最简隐格式、CN格式、双层加权格式、Richardson 格式、网格比、传播因子法(分离变量法) 、传播因子、传播矩阵、谱半径、von Neumann条件、跳点格式、ADI格式、线性椭圆型差分格式、极值原理、比较定理、五点差分格式的相容收敛和稳定性等;
要求: 建立抛物型方程边值问题的差分格式, 计算局部截断误差;
第五章
知识点:左偏心格式、右偏心格式、中心格式、LF格式、LW格式、Wendroff 格式、跳蛙格式、特征线、CFL条件等;
要求: 建立双曲型方程边值问题的差分格式, 计算局部截断误差;
第七章
要求: 会用线性元(线性基)建立常微分方程边值问题的有限元格式
三 练习题
1、 已知显格式21131(
)22
n n n n u u h f f +++-=-
,试证明格式是相容的,并求它的阶。

P39+P41
2、用Taylor 展开原理构造一元函数一阶导数和二阶导数的数值微分公式。

提示:向前、向后和中心差商与一阶导数间关系,二阶中心差商与二阶导数之间的关系 课件
3、用数值微分方法或数值积分方法建立椭圆型方程
2
2
22
(,),
(,),u u f x y x y x
y
∂∂--
=∀∈Ω∂∂ :01,01
x y Ω≤≤≤≤ 内点差分格式。

P75+课件
4、构造椭圆型方程边值问题的差分格式. P101 (4)题
5、构建一维热传导方程2
20,(0)u u Lu a
a t x
∂∂=
-=>∂∂的数值差分格式(显隐格式等)。

参考P132-135相关知识点
6、设有逼近热传导方程2
2
(0)u u Lu a
f a const t
x
∂∂≡
-==>∂∂的带权双层格式
()()1
11111112
2(1)2k k
j
j
k k k k k k
j j j j j j u u a
u u u u u u h θθτ
++++-+-+-⎡⎤=
-++--+⎣⎦
其中[0,1]
θ
∈,试求其截断误差。

并证明当2
1212h
a θτ=-时,截断误差的阶最
高阶为
2
4
()
O h τ+。

P135+P165+课件
7、传播因子法证明抛物型方程2
2(0)u u Lu a
f a const t x
∂∂≡-==>∂∂的最简显隐和六
点CN 格式稳定性。

P156+课件 8、对一阶常系数双曲型方程的初边值问题
0,0,0,0,(,0)(),0,(0,)(),0,
u
u a t T x a t x u x x x u t t t T φψ∂∂⎧+=<≤<<∞>⎪∂∂⎪

=≤<∞⎪⎪=≤≤⎩
试建立左右偏心差分格式。

P185+课件 9、设有逼近双曲型方程
0u u a
t
x
∂∂+=∂∂的双层加权格式
()()1
1111(1)0k k
j
j
k k k k
j j j j u u a
u u u u h
θθτ
+++---⎡⎤+
-+--=⎣⎦, [0,1]θ∈ 试求其截断误差, 并说明当112
2ar
θ=
-
时截断误差为最高阶22()O h τ+. P194
10、对于两点边值问题(),(,)(),()d
du p qu f x a b dx dx
u a u b αβ⎧-
+=∈⎪⎨⎪==⎩
用等距结点线性元推导有
限元方程. 参考P267+P271相关知识点+课件。

相关文档
最新文档