偏微分课程课件4_双曲型方程的差分方法(I)
(高等数学)偏微分方程
第十四章 偏微分方程物理、力学、工程技术和其他自然科学经常提出大量的偏微分方程问题.由于实践的需要和一些数学学科（如泛函分析，计算技术）的发展，促进了偏微分方程理论的发展，使它形成一门内容十分丰富的数学学科.本章主要介绍一阶偏微分方程、线性方程组及二阶线性偏微分方程的理论.在二阶方程中，叙述了极值原理、能量积分及惟一性定理.阐明了一些解的性质和物理意义，介绍典型椭圆型、双曲型、抛物型方程的常用解法：分离变量法，基本解，格林方法,黎曼方法,势位方法及积分变换法.最后，扼要地介绍了有实用意义的数值解法：差分方法和变分方法.§1 偏微分方程的一般概念与定解问题[偏微分方程及其阶数] 一个包含未知函数的偏导数的等式称为偏微分方程.如果等式不止一个，就称为偏微分方程组.出现在方程或方程组中的最高阶偏导数的阶数称为方程或方程组的阶数.[方程的解与积分曲面] 设函数u 在区域D 内具有方程中所出现的各阶的连续偏导数，如果将u 代入方程后，能使它在区域D 内成为恒等式，就称u 为方程在区域D 中的解，或称正规解. ),,,(21n x x x u u = 在n +1维空间),,,,(21n x x x u 中是一曲面，称它为方程的积分曲面. [齐次线性偏微分方程与非齐次线性偏微分方程] 对于未知函数和它的各阶偏导数都是线性的方程称为线性偏微分方程.如()()()()y x f u y x c yuy x b x u y x a ,,,,=+∂∂+∂∂就是线性方程.在线性方程中，不含未知函数及其偏导数的项称为自由项，如上式的f (x,y ).若自由项不为零，称方程为非齐次的.若自由项为零，则称方程为齐次的.[拟线性方程与半线性方程] 如果一个方程，对于未知函数的最高阶偏导数是线性的，称它为拟线性方程.如()()()()()()0,,,,,,,,,,,,22222122211=+∂∂+∂∂+∂∂+∂∂∂+∂∂u y x c y uu y x b x u u y x a yu u y x a y x u u y x a x u u y x a就是拟线性方程，在拟线性方程中，由最高阶偏导数所组成的部分称为方程的主部.上面方程的主部为()()()22222122211,,,,,,yuu y x a y x u u y x a x u u y x a ∂∂+∂∂∂+∂∂如果方程的主部的各项系数不含未知函数，就称它为半线性方程.如()()()()0,,,,,,2222=∂∂+∂∂+∂∂+∂∂y yu y x d x y u y x c yu y x b x u y x a就是半线性方程.[非线性方程] 不是线性也不是拟线性的方程称为非线性方程.如1)()1(222=∂∂+∂∂+yux u u就是一阶非线性偏微分方程.[定解条件] 给定一个方程，一般只能描写某种运动的一般规律，还不能确定具体的运动状态，所以把这个方程称为泛定方程.如果附加一些条件（如已知开始运动的情况或在边界上受到外界的约束）后，就能完全确定具体运动状态，称这样的条件为定解条件.表示开始情况的附加条件称为初始条件，表示在边界上受到约束的条件称为边界条件.[定解问题] 给定了泛定方程（在区域D 内）和相应的定解条件的数学物理问题称为定解问题.根据不同定解条件，定解问题分为三类.1︒ 初值问题 只有初始条件而没有边界条件的定解问题称为初值问题或柯西问题. 2︒ 边值问题 只有边值条件而没有初始条件的定解问题称为边值问题.3︒ 混合问题 既有边界条件也有初始条件的定解问题称为混合问题（有时也称为边值问题）.[定解问题的解] 设函数u 在区域D 内满足泛定方程，当点从区域D 内趋于给出初值的超平面或趋于给出边界条件的边界曲面时，定解条件中所要求的u 及它的导数的极限处处存在而且满足相应的定解条件，就称u 为定解问题的解.[解的稳定性] 如果定解条件的微小变化只引起定解问题的解在整个定义域中的微小变化，也就是解对定解条件存在着连续依赖关系，那末称定解问题的解是稳定的.[定解问题的适定性] 如果定解问题的解存在与惟一并且关于定解条件是稳定的，就说定解问题的提法是适定的.§2 一阶偏微分方程一、 柯西-柯娃列夫斯卡娅定理[一阶偏微分方程的通解] 一阶偏微分方程的一般形式 是0),,,,,,,,(2121=∂∂∂∂∂∂nn x ux u x u u x x x F或()0,,,,,,,211=n n p p p u x x F ，其中()n i x up ii ,,2,1 =∂∂=如解出p 1，可得：p 1 = f (x 1 , x 2 ,…, x n , u , p 2 ,…, p n )当方程的解包含某些“任意元素”(指函数)，如果适当选取“任意元素”时，可得方程的任意解(某些“奇异解”除外)，则称这样的解为通解.在偏微分方程的研究中，重点在于确定方程在一些附加条件(即定解条件)下的解，而不在于求通解.[一阶方程的柯西问题]()()⎪⎩⎪⎨⎧==∂∂=n x x n n x x u p p u x x x f x u,,|,,,,,,,22211011 ϕ 称为柯西问题，式中),,(2n x x ϕ为已知函数，对柯西问题有如下的存在惟一性定理.[柯西－柯娃列夫斯卡娅定理] 设 f ( x 1 , x 2 ,, x n , u , p 2 ,, p n ) 在点 ( x 10 , x 20 ,, x n 0 , u 0 , p 20 ,, p n 0 ) 的某一邻域内解析，而),,(2n x x ϕ在点( x 20 ,, x n 0 ) 的某邻域内解析，则柯西问题在点 ( x 10 ,, x n 0 ) 的某一邻域内存在着惟一的解析解.这个定理应用的局限性较大，因它要求f 及初始条件都是解析函数，一般的定解问题未必能满足这种条件.对高阶方程也有类似定理.二、 一阶线性方程1． 一阶齐次线性方程[特征方程∙特征曲线∙初积分(首次积分)] 给定一阶齐次线性方程在有些书中写作0),,,,,,,,,(121=∂∂∂∂∂∂nn x u x u t u u x x x t F()()0,,,,,,211211=∂∂++∂∂nn n n x u x x x a x u x x x a (1) 式中a i 为连续可微函数，在所考虑的区域内的每一点不同时为零(下同).方程组()n i ix x x a tx ,,,d d 21 = ( i = 1,2,, n ) 或()()()n n n n n x x x a x x x x a x x x x a x ,,,d ,,,d ,,,d 2121222111 === (2)称为一阶齐次线性偏微分方程的特征方程.如果曲线l : x i = x i (t ) ( i =1,2,, n )满足特征方程(2)，就称曲线l 为一阶齐次线性方程的特征曲线.如果函数ψ ( x 1 , x 2 ,, x n )在特征曲线),,2,1()(n i t x x i i ==上等于常数，即ψ ( x 1(t ) , x 2(t ) ,, x n (t ) ) = c就称函数ψ ( x 1, x 2,, x n )为特征方程(2)的初积分(首次积分). [齐次方程的通解]1o 连续可微函数u = ψ ( x 1, x 2,, x n ) 是齐次线性方程(1)的解的充分必要条件是： ψ ( x 1, x 2,, x n )是这个方程的特征方程的初积分.2o 设ψi ( x 1 , x 2 ,, x n ) ( i = 1,2,, n 1-) 是特征方程(2)在区域D 上连续可微而且相互独立的初积分(因此在D 内的每一点，矩阵⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂---n n n n n n x x x x x x x x x 121112221212111ψψψψψψψψψ 的秩为n 1-) ，则u = ω ( ψ1 ( x 1 , x 2 ,, x n ) ,, ψn -1 ( x 1 , x 2 ,, x n ) )是一阶齐次线性方程(1)的通解，其中ω为n 1-个变量的任意连续可微函数. [柯西问题] 考虑方程的柯西问题()()⎪⎩⎪⎨⎧==∂∂==∑n x x ni i n i x x u x u x x x a ,,|0,,,2121011 ϕ 式中ϕ ( x2 ,, x n )为已知的连续可微函数.设ψi ( x 1 , x 2 ,, x n ) ( i = 1,2,, n 1-) 为特征方程的任意n 1-个相互独立的初积分，引入参变量 i ψ (1,,2,1-=n i )，从方程组()()()⎪⎪⎩⎪⎪⎨⎧===--120112201212011,,,,,,,,,n n n n n x x x x x x x x x ψψψψψψ 解出x 2 ,, x n 得()()⎪⎩⎪⎨⎧==--12112122,,,,,,n n nn x x ψψψωψψψω 则柯西问题的解为u = ϕ ( ω2 ( ψ1 , ψ2 ,, ψn -1 ) ,, ωn ( ψ1 , ψ2 ,, ψn -1 ) )2． 非齐次线性方程它的求解方法与拟线性方程相同.三、 一阶拟线性方程一阶拟线性方程为()()∑==∂∂ni n i n i u x x x R x uu x x x a 12121,,,,,,,, 其中a i 及R 为x 1 , x 2 ,, x n , u 的连续可微函数且不同时为零. [一阶拟线性方程的求解和它的特征方程]()()⎪⎩⎪⎨⎧===u x x x R t un i u x x x a t x n n i i,,,,d d ),,2,1(,,,,d d 2121 或()()()u x x R uu x x a x u x x a x n n n n n ,,,d ,,,d ,,,d 11111 ===为原拟线性方程的特征方程.如果曲线l : x i = x i (t ) ( i =1,2,, n ) , u = u (t ) 满足特征方程，则称它为拟线性方程的特征曲线.设 ψi ( x 1 ,, x n ,u ) ( i = 1,2,, n ) 为特征方程的n 个相互独立的初积分，那末对于任何连续可微函数ω，ω ( ψ1 ( x 1,, x n , u ) , ψ2 ( x 1,, x n , u ) ,, ψn ( x 1,, x n , u ) ) = 0都是拟线性方程的隐式解.[柯西问题] 考虑方程的柯西问题()()()⎪⎩⎪⎨⎧==∂∂==∑n x x ni n i ni x x u u x x x R x u u x x x a ,,|,,,,,,,,212121011 ϕ ϕ为已知的连续可微函数.设 ψ1 ( x 1 , x 2 ,, x n , u ) ,, ψn ( x 1 , x 2 ,, x n , u ) 为特征方程的n 个相互独立的初积分，引入参变量 n ψψψ,,,21 , 从()()()⎪⎪⎩⎪⎪⎨⎧===nn n n n u x x x u x x x u x x x ψψψψψψ,,,,,,,,,,,,2012201212011解出 x 2 ,, x n , u()()()⎪⎪⎩⎪⎪⎨⎧===n n n n n u x x ψψψωψψψωψψψω,,,,,,,,,21212122 则由()()()()()()()0,,,,,,,,,,,,,,,,,,,,,,2121221221121=-≡n n n n n n u x x x u x x x u x x x V ψψψωψψψωϕψψω给出柯西问题的隐式解.四、 一阶非线性方程[完全解·通解·奇异解] 一阶非线性方程的一般形式为()()n i x up p p p u x x x F ii n n ,,2,10,,,,,,,,2121 =∂∂== 若一阶偏微分方程的解包含任意n 个独立的常数，则称这样的解为完全解(全积分). 若V ( x 1, x 2 ,, x n , u , c 1 , c 2,, c n ) = 0为方程的完全解，从()n i c VV i,,2,10,0 ==∂∂= 消去c i ，若得一个解，则称它为方程的奇异解(奇积分).以两个独立变量为例说明完全解与通解、奇异解的关系，设方程()yzq x z p q p z y x F ∂∂=∂∂==,,0,,,,有完全解V (x ,y ,z ,a ,b )=0 ( a ,b 为任意常数),则方程等价于从方程组()⎪⎪⎩⎪⎪⎨⎧=∂∂+∂∂=∂∂+∂∂=0,00,,,,q z Vy V p z V x V b a z y x V 消去a ,b 所得的方程.利用常数变易法把a ,b 看作x , y 的函数，将V (x ,y ,z ,a ,b )=0求关于x , y 的偏导数，得00=∂∂⋅∂∂+∂∂⋅∂∂+∂∂+∂∂=∂∂⋅∂∂+∂∂⋅∂∂+∂∂+∂∂ybb V y a a V q z V y V xbb V x a a V p z V x V那末0,0=∂∂⋅∂∂+∂∂⋅∂∂=∂∂⋅∂∂+∂∂⋅∂∂yb b V y a a V x b b V x a a V 与V=0联立可确定a ,b .有三种情况：1︒ 0≡∂∂≡∂∂bVa V ，将其与V (x ,y ,z ,a ,b )=0联立可确定不含任意常数的奇异解. 2︒ 如0=∂∂=∂∂=∂∂=∂∂yb x b y a x a ，即回到完全解. 3︒ 当0/,0/≡∂∂≡∂∂b Va V 时，必有()()0,,=∂∂y x b a ，这时，如果不属于情形2︒ ，则a 与b 存在函数关系：b=ω(a )，这里ω为任意可微函数，并从方程V (x ,y ,z ,a ,b )=0和()∂∂∂∂ωV a Vba +'=0消去a ,b ，可确定方程的通解.定理 偏微分方程的任何解包含在完全解内或通解内或奇异解内. [特征方程·特征带·特征曲线·初积分] 在一阶非线性方程：()F x x x u p p p n n 12120,,,,,,,, =中，设F 对所有变量的二阶偏导数存在且连续，称()n i uFp x F t p p F p t u p Ft x i i i ni iii i ,,2,1)(d d d d ,1 =∂∂+∂∂-=∂∂=∂∂=∂∂∑=或u F p x F p u F p x F p p Fp up F x p F xp F x n nnni i i nn ∂∂+∂∂-==∂∂+∂∂-=∂∂=∂∂==∂∂=∂∂∑=d d d d d d 11112211为非线性方程的特征方程.设特征方程的解为x i =x i (t ), u=u (t ), p i =p i (t ) (i =1,2,…,n )称它为非线性方程的特征带.在x 1,x 2,, x n ,u 空间的曲线x i =x i (t ), u=u (t ) (i=1,2,…,n )称为非线性方程的特征曲线.如果函数()n n p p p u x x x G ,,,,,,,,2121 在特征方程的任一解x i =x i (t ) (i =1,2,, n ), u=u (t ), p i =p i (t ) (i =1,2,, n )上等于常数，即()()()()()()()()G x t x t x t u t p t p t p t C n n 1212,,,,,,,, =那末函数()n n p p p u x x x G ,,,,,,,,2121 称为特征方程的初积分.[求完全解的拉格朗日－恰比方法] 考虑两个变量的情况.对于方程F (x ,y ,z ,p ,q )=0，选择使雅可比式()()0,,≠∂∂q p G F 的一个初积分G (x ,y ,z ,p ,q ).解方程组()()F x y z p q G x y z p q a,,,,,,,,==⎧⎨⎪⎩⎪0(a 为任意常数) 得p (x ,y ,z ,a )及q (x ,y ,z ,a ).则方程d z=p d x+q d y的通解V (x ,y ,z ,a ,b )=0(b 是积分d z=p d x+q d y 出现的任意常数)就是方程F (x ,y ,z ,p ,q )=0的完全解.例 求方程()z p q x y 22222+=+的完全解.解 方程的特征方程为()()()qy x z y qp q p z x p q p z z q z y p z x 22222222222d 22d 2d 2d 2d +-=+-=+== 这里成立zpxx p z z p d d d =+ 所以特征方程的一个初积分为z 2p 2 －x 2 .解方程组 ()()z p q x y z p x a22222222+-+=-=⎧⎨⎪⎩⎪ (a 为任意常数) 得 p a x zq y az=+=-22, 积分微分方程dz a x zdx y azdy =++-22 得完全解z x x a y y a a x x a y y ab 22222=++-++++-+ln(b 为任意常数)[某些容易求完全解的方程] 1︒ 仅含p ,q 的方程F (p ,q )=0G =p 是特征方程的一个初积分.从F (p ,q )=0与p=a (a 为任意常数)得q=ψ(a )，积分d z=a d x+ψ(a )d y得完全解z=ax+ψ(a )y+b (b 为任意常数)2︒ 不显含x ,y 的方程F (z ,p ,q )=0 特征方程为zFqqz F p p q F q p F p z q F y p F x ∂∂-=∂∂-=∂∂+∂∂=∂∂=∂∂d d d d d 因此q d p-p d q =0，显然G qp=为一个初积分，由F (z ,p ,q )=0，q=pa (a 为任意常数)解得p=ψ(z ,a ).于是由d z=ψ(z ,a )d x+a ψ(z ,a )d y得()⎰++=b ay x a z z,d ψ (b 为任意常数)可确定完全解.3︒ 变量分离形式的方程()f x p i i i i n,=∑=10特征方程为n n n n i i iin n n x f p x f p p f p z p f x p f x ∂∂-==∂∂-=∂∂=∂∂==∂∂∑=d d d d d 1111111 可取初积分G i =f i (x i ,p i ) , (i =1,2,, n ).从f i (x i ,p i )=a i (i =1,2,, n )解出p i =ϕi (x i ,a i )得完全解()∑⎰=+=ni i i i i b x a x z 1d ,ϕ式中a i ,b 为任意常数，且a i i n=∑=10.[克莱罗方程] 方程()z p x f p p p i i n i n=+=∑121,,,称为克莱罗方程，其完全解为()z c x f c c c i i n i n=+=∑121,,,对c i 微分得x fc i i=-∂∂ (i =1,2,…,n ) 与完全解的表达式联立消去c i 即得奇异解.例 求方程z －xp －yq －pq =0的完全解和奇异解. 解 这是克莱罗方程，它的完全解是z=ax+by+ab对a,b 微分，得x=－b,y=－a ，消去a ,b 得奇异解z=－xy[发甫方程] 方程P (x,y,z )d x+Q (x,y,z )d y+R (x,y,z )d z=0 (1)称为发甫方程，如果P,Q,R 二次连续可微并满足适当条件，那末方程可积分.如果可积分成一关系式时，则称它为完全可积.1︒ 方程完全可积的充分必要条件 当且仅当P,Q,R 满足条件0)()()(=∂∂-∂∂+∂∂-∂∂+∂∂-∂∂yP x Q R x R z P Q z Q y R P (2) 时，存在一个积分因子μ(x,y,z )，使d U 1=μ(P d x+Q d y+R d z )从而方程的通解为U 1(x,y,z )=c特别，当0,0,0=∂∂-∂∂=∂∂-∂∂=∂∂-∂∂yP x Q x R z P z Q y R 时，存在一个函数U (x,y,z )满足 zU R y U Q x U P ∂∂=∂∂=∂∂=,,从而 d U=P d x+Q d y+R d z 所以方程的通解为U (x,y,z )=c所以完全可积的发甫方程的通解是一单参数的曲面族.定理 设对于发甫方程(1)在某区域D 上的完全可积条件(2)成立，则对D 内任一点M (x,y,z )一定有方程的积分曲面通过，而且只有一个这样的积分曲面通过. 2︒ 方程积分曲面的求法设完全可积条件(2)成立.为了构造积分曲面，把z 看成x,y 的函数(设R (x,y,z )≠0)，于是原方程化为y RQ x R P z d d d --=由此得方程组()()()()⎪⎪⎩⎪⎪⎨⎧≡-=∂∂≡-=∂∂4,,3,,11z y x Q R Q y z z y x P R P xz发甫方程(1)与此方程组等价.把方程(3)中的y 看成参变量，积分后得一个含有常数 c 的通解 ()cy x z ~;,ϕ= 然后用未知函数()~cy 代替常数 c ，将()()z x y c y =ϕ,;~代入方程(4)，在完全可积的条件下，可得()~cy 的一个常微分方程，其通解为 ()()~,cy y c =ψ c 为任意常数，代回()()z x y cy =ϕ,;~中即得发甫方程的积分曲面 z=ϕ(x,y,ψ(y,c ))由于发甫方程关于x,y,z 的对称性，在上面的讨论中，也可把x 或y 看成未知函数，得到同样的结果.例 求方程yz d x+2xz d y+xy d z=0的积分曲面族.解 容易验证完全可积条件成立，显然存在一个积分因子μ=1xyz，用它乘原方程得 0d d 2d =++zz y y x x 积分后得积分曲面族xy 2z=c也可把方程化为等价的方程组⎪⎪⎩⎪⎪⎨⎧-=∂∂-=∂∂y z yz x z xz 2 把y 看成参变量，积分xzx z -=∂∂得通解 zx c= 用未知函数()~cy 代替 c ，将()y c zx ~=代入方程y z y z 2-=∂∂得 ()()yy cy y c ~2d ~d -= 积分后有()~cy c y =2所以原方程的积分曲面族是xy 2z=c五、 一阶线性微分方程组[一阶线性偏微分方程组的一般形式] 两个自变量的一阶线性方程组的形式是()n i F u C x u B t u A i n j j ij n j n j jij j ij ,,2,10111 ==++∂∂+∂∂∑∑∑=== 或()n i f u b x u a t u i n j j ij n j j ij i,,2,1011 ==++∂∂+∂∂∑∑== (1) 其中A ij ,B ij ,C ij ,F i ,a ij ,b ij ,f i 是(x,t )的充分光滑函数. [特征方程·特征方向·特征曲线]⎩⎨⎧=≠==-j i j i t xa ij ij ij ,1,0,0)d d det(δδ称为方程组(1)的特征方程.在点(x,t )满足特征方程的方向txd d 称为该点的特征方向.如果一条曲线l ，它上面的每一点的切线方向都和这点的特征方向一致，那末称曲线l 为特征曲线. [狭义双曲型方程与椭圆型方程] 如果区域D 内的每一点都存在n 个不同的实的特征方向，那末称方程组在D 内为狭义双曲型的.如果区域D 内的每一点没有一个实的特征方向，那末称方程组在D 内为椭圆型的. [狭义双曲型方程组的柯西问题] 1︒ 化方程组为标准形式——对角型因为det(a ij -δij λ)=0有n 个不同的实根λ1(x,t ) ,, λn (x,t )，不妨设),(),(),(21t x t x t x n λλλ<<<那末常微分方程()()n i t x txi ,,2,1,d d ==λ 的积分曲线l i (i =1,2,…,n )就是方程组(1)的特征曲线. 方程()()aijk ij k i i n-==∑λδλ1的非零解(λk (1) ,, λk (n ))称为对应于特征方向λk 的特征矢量. 作变换()()n i u v nj jj i i ,,2,11==∑=λ可将方程组化为标准形式——对角型()()()()n i t x v t x a x v t x t v i nj j ij ii i ,,2,1,,,1=+=∂∂+∂∂∑=βλ 所以狭义双曲型方程组可化为对角型，而一般的线性微分方程组(1)如在区域D 内通过未知函数的实系数可逆线性变换可化为对角型的话，(此时不一定要求 λi 都不相同)，就称这样的微分方程组在D 内为双曲型的. 2︒ 对角型方程组的柯西问题 考虑对角型方程组的柯西问题()()()()()()n i x x v t x v t x a x v t x tv i inj i j ij i i i,,2,10,,,,1 =⎪⎩⎪⎨⎧=+=∂∂+∂∂∑=ϕβλ ϕi (x )是[a,b ]上的连续可微函数.设αij ,βi ,λi 在区域D 内连续可微，在D 内可得相应的积分方程组()()()n i tv x t x v il i n j j ij i i i ,,2,1d ,~1 =⎥⎦⎤⎢⎣⎡++=⎰∑=βαϕ 式中 l i 为第i 条特征曲线l i 上点(x,t )与点(x i ,0)之间的一段，(x i ,0)为l i与x 轴上[a,b ]的交点.上式可以更确切地写为()()[]()[]()[]()[]⎰∑⎭⎬⎫⎩⎨⎧+⋅+==t n j i i i j i ij i i i t x x t x x v t x x a t x x t x v 01d ,,,,,,,,,0,,,τττβττττϕ(i =1,2,, n )式中x i =x i (x ︒,t ︒,t )为过点(x ︒,t ︒)的第i 条特征曲线，利用逐次逼近法可解此积分方程.为此令()()()[]()()()()[]()[]()()[]()[]()()()()[]()[]()()[]()[]()n i t x x t x x v t x x a t x x t x v n i t x x t x x v t x x a t x x t x v n i t x x t x v i i tnj i k j i ij i i k ii i tnj i j i ij i i ii i i ,,2,1d ,,,,,,,,,0,,,,,2,1d ,,,,,,,,,0,,,,,2,10,,,}{}{01101010=+⋅+==+⋅+===⎰∑⎰∑=-=τττβττττϕτττβττττϕϕ序列{v i (k )} (k =0,1,2 ,)一致收敛于积分方程的连续可微解v i (x,t ) (i =1,2,, n )，这个v i (x,t )也就是对角型方程组的柯西问题的解.设在区域D 内对角型方程组的柯西问题的解存在，那末解与初值有下面的关系：(i) 依赖区间：过D 中任意点M (x,t )作特征曲线l 1,l n ，交x 轴于B,A ，称区间[A,B ]为M 点的依赖区间(图14.1(a ))，解在M 点的值由区间[A,B ]的初值确定而与[A,B ]外的初值无关. (ii) 决定区域：过点A,B 分别作特征曲线l n ,l 1,称l n ,l 1 与区间[A,B ]围成的区域D 1为区间[A,B ]的决定区域(图14.1(b ))，在区域D 1中解的值完全由[A,B ]上的初值决定.(iii) 影响区域：过点A,B 分别作特征曲线l 1,l n ，称l 1,l n 与[A,B ]围成的区域D 2为区间[A,B ]的影响区域(图14.1(c )).特别当区间[A,B ]缩为一点A 时，A 点的影响区域为D 3(图14.1(d )).在区域D 2中解的值受[A,B ]上的初值影响，而在区域D 2外的解的值则不受[A,B ]上的初值影响.图14.1[线性双曲型方程组的边值问题] 以下列线性方程组来说明：()⎪⎪⎩⎪⎪⎨⎧<++=∂∂+∂∂++=∂∂+∂∂2122221111λλλλc v b u a x v t v c v b u a xu t u (1) 1︒ 第一边值问题(广义柯西问题) 设在平面(x,t )上给定曲线段⋂AB ，它处处不与特征方向相切.过A,B 分别引最左和最右的特征曲线l 1及l 2.要求函数u (x,t ),v (x,t )在⋂AB ，l 1及l 2围成的闭区域D 上满足方程组，且在⋂AB 上取给定的函数值(图14.2(a )).2︒ 第二边值问题(古沙问题) 设l 1是过P 点的第一族特征线，l 2是第二族特征线，在l 1的一段PA 上给定v (x,t )的数值，在l 2的一段PB 上给定u (x,t )的数值，过A 点作第二族特征线，过B 点作第一族特征线相交于Q .求在闭区域PAQB 上方程组的解(图14.2(b )).3︒ 第三边值问题 设AB 为非特征曲线的曲线弧，AC 为一特征线弧，且在AB 与AC 之间不存在过A 点的另外特征曲线,过C 点作第二族特征线与过B 点的第一族特征线交于E 点，在AC 上给定v (x,t )的数值，在AB 上给定u (x,t )的数值，求ACEBA 所围成的闭区域D 上的方程组的解(图14.2(c)).图14.2[边值问题的近似解——特征线法] 以上定解问题，可用逐步逼近法求解，也可用特征线法求解的近似值.以第一边值问题为例说明.在曲线AB 上取n 个分点A 1,A 2,, A n ，并记A 为A 0，B 为A n +1，过A 0按A 0的第二特征方向作直线与过A 1按A 1的第一特征方向作直线相交于B 0；过A 1按A 1第二特征方向作直线与过A 2按A 2的第一特征方向作直线相交于B 1 ,最后得到B n (图14.3).用如下的近似公式来确定方程组(1)的解u (x,t ),v (x,t )在B i (i =0,1,2,…,n )的数值：()()()()()()(){}()[]()()()()()()(){}()[]u B u A B A a A u A b A v A c A A v B v A B A a A u A b A v A c A A i i i i i i i i i i i i i i i i i i i i -=++⨯+-=++⨯+⎧⎨⎪⎩⎪+++++++--11111111112122212121211λλ图14.3于是在一个三角形网格的节点上得到u,v 的数值.再经过适当的插值，当n 相当大，A i 、A i +1的距离相当小时，就得到所提问题的足够近似的解.[特殊形式的拟线性方程组——可化约系统] 一般的拟线性方程组的问题比较复杂，目前研究的结果不多，下面介绍一类特殊形式的拟线性方程组——可化约系统.如果方程组⎪⎪⎩⎪⎪⎨⎧=∂∂+∂∂+∂∂+∂∂=∂∂+∂∂+∂∂+∂∂0022221111x v D t v C x u B tu A xv D t v C x u B t uA 中所有的系数只是u,v 的函数，称它为可化约系统. 考虑满足条件()()0,,≠∂∂t x v u 的方程组的解u=u (x,t ),v=v (x,t ).x,t 可以表示成u,v 的函数，且()()()()()()()()v u t x u t x v v u t x u x t v v u t x v tx u v u t x v x t u ,,,,,,,,,,∂∂∂∂=∂∂∂∂∂∂-=∂∂∂∂∂∂-=∂∂∂∂∂∂=∂∂ 原方程化为⎪⎪⎩⎪⎪⎨⎧=∂∂+∂∂-∂∂-∂∂=∂∂+∂∂-∂∂-∂∂0022221111u t D u x C v t B vx A ut D u x C v t B v xA 这是关于自变量u,v 的线性方程组.这样就把求拟线性方程组满足()()0,,≠∂∂t x v u 的解，化为解线性方程组的问题.而此线性方程组满足条件()()0,,≠∂∂v u t x 的解，在(x,t )平面上的象即为原来拟线性方程组的解.§3 二阶偏微分方程一、 二阶偏微分方程的分类、标准形式与特征方程考虑二阶偏微分方程()0),,,,,,(111,2=∂∂∂∂+∂∂∂∑=nnnj i j i ij x u x u u x x F y x u x a (1) 式中a ij (x )=a ij (x 1,x 2,…,x n )为x 1,x 2,…,x n 的已知函数.[特征方程·特征方向·特征曲面·特征平面·特征锥面]代数方程()01,=∑=nj i jiijaa x a称为二阶方程(1)的特征方程；这里a 1,a 2,…,a n 是某些参数，且有012≠∑=ni i a .如果点x ︒=(x 1︒,x 2︒,…,x n ︒)满足特征方程，即()01,o =∑=nj i jiijaa x a则过x ︒的平面()01o=-∑=nk kk k x x a 的法线方向l :(a 1,a 2,…,a n )称为二阶方程的特征方向；如果一个(n 1-)维曲面，其每点的法线方向都是特征方向，则称此曲面为特征曲面；过一点的(n 1-)维平面，如其法线方向为特征方向，则称这个平面为特征平面，在一点由特征平面的包络组成的锥面称为特征锥面.[n 个自变量方程的分类与标准形式] 在点P (x 1︒,x 2︒,…,x n ︒)，根据二次型()∑=nj i jinijaa x x x a 1,o o 2o 1,,, (a i 为参量)的特征根的符号，可将方程分为四类：(i) 特征根同号，都不为零，称方程在点P 为椭圆型.(ii) 特征根都不为零，有n 1-个具有同一种符号 ，余下一个符号相反，称方程在点P 为双曲型.(iii) 特征根都不为零，有m n -个具有同一种符号(n >m >1)，其余m 个具有另一种符号，称方程在点P 为超双曲型.(iv) 特征根至少有一个是零，称方程在点P 为抛物型.若在区域D 内每一点方程为椭圆型，双曲型或抛物型，则分别称方程在区域D 内是椭圆型、双曲型或抛物型.在点P 作自变量的线性变换可将方程化为标准形式：椭圆型：∑==+∂∂ni ix u1220Φ双曲型：∑==+∂∂-∂∂n i ix ux u 22120Φ超双曲型：()10112222>>=+∂∂-∂∂∑∑=+=m n x ux u m i nm i ii Φ抛物型：()00122>=+∂∂∑-=m x umn i iΦ 式中Φ为不包含二阶导数的项.[两个自变量方程的分类与标准形式] 方程的一般形式为0,,,,222222122211=⎪⎪⎭⎫⎝⎛∂∂∂∂+∂∂+∂∂∂+∂∂y u x u u y x F y u a y x u a x u a (2) a 11,a 12,a 22为x ,y 的二次连续可微函数，不同时为零. 方程a 11d y 22-a 12d x d y +a 22d x 2=0称为方程(2)的特征方程.特征方程的积分曲线称为二阶方程(2)的特征曲线. 在某点P (x 0,y 0)的邻域D 内，根据Δ=a 122－a 11a 12的符号将方程分类： 当Δ>0时，方程为双曲型； 当Δ=0时，方程为抛物型； 当Δ<0时，方程为椭圆型.在点P 的邻域D 内作变量替换，可将方程化为标准形式：（i ） 双曲型：因Δ>0，存在两族实特征曲线11),(c y x =ϕ，22),(c y x =ϕ，作变换),(1y x ϕξ=，),(2y x ϕη=和,,ηηξ-=+=s t s 方程化为标准形式),,,,(2222tus u u t s t u s u ∂∂∂∂=∂∂-∂∂Φ或),,,,(12ηξηξΦηξ∂∂∂∂=∂∂∂uu u u （ii ） 抛物型: 因Δ=0，只存在一族实的特征曲线c y x =),(ϕ，取二次连续可微函数),(y x ψ，使0),(),(≠∂∂y x ψϕ，作变换),(y x ϕξ=，),(y x ψη=，方程化为标准形式),,,,(222ηξηξΦη∂∂∂∂=∂∂uu u u （iii ） 椭圆型：因Δ<0，不存在实特征曲线，设c y x i y x y x =+=),(),(),(21ϕϕϕ为11221121212d d a a a a a x y -+=的积分，y x ϕϕ,不同时为零，作变量替换),(1y x ϕξ=，),(2y x ϕη=,方程化为标准形式),,,,(32222ηξηξΦηξ∂∂∂∂=∂∂+∂∂uu u u u二、 极值原理·能量积分·定解问题的惟一性定理椭圆型方程、抛物型方程的极值原理及双曲型方程的能量守恒原理是相应方程的解所具有的最基本性质之一，在定解问题的研究中起着重要的作用. [椭圆型方程的极值原理与解的惟一性定理]1︒ 极值原理 设D 为n 维欧氏空间E n 的有界区域，S 是D 的边界，在D 内考虑椭圆型方程()()()()x x x x f u c x ub x x u a Lu ni i i n j i j i ij =+∂∂+∂∂∂≡∑∑==11,2式中a ij (x ),b i (x ),c (x ),f (x )在D 上连续，c (x )≤0且二次型()∑=nj i j i ij a a a 1,x 正定，即存在常数μ>0，对任意x D ∈和任意的a i 有()∑∑==≥ni i nj i jiija aa a 121,μx定理1 设u (x )为D 内椭圆型方程的解，它在D 内二次连续可微，在D 上连续，且不是常数，如f (x )≤0（或f (x )≥0），则u (x )不能在D 的内点取非正最小值（或非负最大值）. 如果过边界S 上的任一点P 都可作一球，使它在P 点与S 相切且完全包含在区域D 内，则有 定理2 设u (x )为椭圆型方程在D 内二次连续可微，在D 上连续可微的解，且不是常数，并设f (x )≤0（或f (x )≥0）.若u (x )在边界S 上某点M 处取非正最小值（或非负最大值），只要外法向导数错误！未定义书签。
第九章 偏微分方程差分方法
第9章 偏微分方程的差分方法含有偏导数的微分方程称为偏微分方程。
由于变量的增多和区域的复杂性，求偏微分方程的精确解一般是不可能的，经常采用数值方法求方程的近似解。
偏微分方程的数值方法种类较多，最常用的方法是差分方法。
差分方法具有格式简单，程序易于实现，计算量小等优点，特别适合于规则区域上偏微分方程的近似求解。
本章将以一些典型的偏微分方程为例，介绍差分方法的基本原理和具体实现方法。
9.1椭圆型方程边值问题的差分方法9.1.1 差分方程的建立最典型的椭圆型方程是Poisson （泊松）方程G y x y x f yux u u ∈=∂∂+∂∂-≡∆-),(),,()(2222 （9.1）G 是x ,y 平面上的有界区域，其边界Γ为分段光滑的闭曲线。
当f (x ,y )≡0时，方程（9.1）称为Laplace(拉普拉斯）方程。
椭圆型方程的定解条件主要有如下三种边界条件第一边值条件 ),(y x u α=Γ （9.2） 第二边值条件),(y x nuβ=∂∂Γ （9.3） 第三边值条件 ),()(y x ku nuγ=+∂∂Γ （9.4） 这里，n 表示Γ上单位外法向，α(x,y ),β(x,y ),γ(x,y )和k (x,y )都是已知的函数，k (x,y )≥0。
满足方程（9.1）和上述三种边值条件之一的光滑函数u (x ,y )称为椭圆型方程边值问题的解。
用差分方法求解偏微分方程，就是要求出精确解u (x ,y )在区域G 的一些离散节点（x i ，y i ）上的近似值u i ,j ≈(x i ，y i ）。
差分方法的基本思想是，对求解区域G 做网格剖分，将偏微分方程在网格节点上离散化，导出精确解在网格节点上近似值所满足的差分方程，最终通过求解差分方程，通常为一个线性方程组，得到精确解在离散节点上的近似值。
设G ={0<x <a , 0<y <b }为矩形区域，在x ,y 平面上用两组平行直线x =ih 1, i =0,1,…,N 1, h 1=a /N 1 y =jh 2, j =0,1,…,N 2, h 2=b /N 2将G 剖分为网格区域，见图9-1。
第十章 偏微分方程数值解3
2 kh 2 2 2 kh 2 kh 1 2ar sin 4a r sin 2 cos 2 2 2 kh 1 4ar(1 ar) sin 2
2
当a>0时，恒有 G ( , k )
2
1 ，格式（10.4.7）不稳定 ；
2
当a<0且ar 1时，G ( , k )
2
1 1 1 1 u n1 2u n 1 u n1 u n1 2u n 1 u n1 a j j j j j j 2 2 2 h h
1 1 a 2 2 u n1 2(1 a 2 2 )u n 1 a 2 2 u n1 j j j 1 1 4u n 2(1 a 2 2 )u n 1 a 2 2 (u n1 u n1 ) j j j j
要求：差分格式，定解条件，区域划分，边值条件的处理；
截断误差阶，稳定性结论
N为正整数。在t = 0上的结点称为边界结点，其余所有属于
c x , 0 t T
内的结点称为内部结点。
10.4.2 一阶双曲方程的差分法
a）迎风格式 ut (xj, tn)用向前差商代替，ux（xj, tn）用向前或向后 差商代替，得
n 1 j
u
或
u
n j
u
n 1 j
期末考试闭卷部分： 1. 第二章 非线性方程求根 ：二分法、迭代法、牛顿法和弦截法 要求：根的存在，公式，收敛性条件的判别 2. 第三章 解线性方程组的直接法：掌握Gauss消元法进行到底 的条件，矩阵三角分解定理的条件和结论，向量和
矩阵的范数，方程组的条件数与病态方程组的求解
3. 第四章 解线性方程组的迭代法：雅可比迭代法，高斯-赛德尔 迭代法；要求：求解公式，收敛条件。
2-双曲型方程的差分方法
其截断误差是
n 1 n 1 n n u u u u a j 1 j 1 j 1 j 1 0 2 2 h 2 h
T O( h )
2 2
其增长因子是
1 1 2 ia sin kh G 1 1 2 ia sin kh
2 2 2 1 1 a sin kh 4 G 1 2 1 2 2 1 4 a sin kh 2
),
a0 a0
1 n n n un u a ( u u j j j 1 j ),
也可写成统一形式
1 n n n n n n 1 1 un u a ( u u ) a ( u 2 u u j j j 1 j 1 j 1 j j 1 ) 2 2
u ( P) u (Q) u (C ) a u (C ) u ( B) 1 a (1 a ) u ( B) 2u (C ) u ( D) 2
对应差分格式即为Lax-Wendroff格式
2 2 a a n 1 n n n n n n uj uj u j 1 u j 1 u j 1 2u j u j 1 2 2
代入前面的表达式有
u
n 1 j
u
n j
a
u
n j 1
u
n j 1
2h
u u a x t j
n
2h 2
n n n 2 2 2 a2 u 2 u u O ( h h ) j j 1 j 1
得到二阶精度的显式格式，即Lax-Wendroff格式
隐式格式
u u
n j
n 1 j
二阶双曲方程显、隐差分法
一、研究对象
1. 研究的对象——二阶双曲型方程.
2 2 u( x , t ) 2 u( x , t ) a f ( x , t ), 0 x 1, 0 t T , 2 2 x t u u ( x , 0) ( x ), ( x , 0) ( x ), 0 x 1, t u(0, t ) ( t ), u(1, t ) ( t ), 0 t T ,
k u 将数值解 i 代替精确解 u( xi , tk ) 并忽略高阶小项， 则第四步，可以建立以下显差分格式：
k k k uik 1 2uik uik 1 2 ui 1 2ui ui 1 a f ( xi , t k ), 1 i m 1, 1 k n 1, 2 2 h 0 ui1 ui0 ( xi ), 0 i m , ui ( xi )， k k u0 ( t k )，um ( t k ), 1 k n.
从而得增长因子为
G 1 2r sin
2
h
2
4r sin
2
h
2
( r sin
2
h
2
1)
如果 r 1 ，则
G 1 2r sin
2
h
2
i 4r sin
2
h
2
(1 r sin
2
h
2
)
从而 | G | 1 ，满足Von Neumann 条件。 但此时由于 | G | 1 ，所以Von Neumann条件只 是差分格式稳定的必要条件而非充分条件。当 r <1
第十四章SECTION4偏微分方程的数值解法
§4 偏微分方程的数值解法一、 差分法差分法是常用的一种数值解法.它是在微分方程中用差商代替偏导数，得到相应的差分方程，通过解差分方程得到微分方程解的近似值. 1. 网格与差商在平面 (x ,y )上的一以S 为边界的有界区域D 上考虑定解问题.为了用差分法求解，分别作平行于x 轴和y 轴的直线族.⎩⎨⎧====jhy y ihx x i i (i ,j =0,±1,±2,…,±n ) 作成一个正方形网格，这里h 为事先指定的正数，称为步长；网格的交点称为节点，简记为(i ,j ).取一些与边界S 接近的网格节点，用它们连成折线S h ，S h 所围成的区域记作D h .称D h 内的节点为内节点，位于S h 上的节点称为边界节点（图14.7）.下面都在网格D h + S h 上考虑问题：寻求各个节点上解的近似值.在边界节点上取与它最接近的边界点上的边值作为解的近似值，而在内节点上，用以下的差商代替偏导数：()()[]()()[]()()()[]()()()[]()()()[]y x u h y x u y h x u h y x u hy x u h y x u y x u h y x u hy u y h x u y x u y h x u h x u y x u h y x u hyu y x u y h x u h x u ,),(,,1,,2,1,,2,1,,1,,122222222++-+-+≈∂∂∂-+-+≈∂∂-+-+≈∂∂-+≈∂∂-+≈∂∂注意， 1︒ 式中的差商()()[]y x u y h x u h ,,1-+称为向后差商，而()()[]y h x u y x u h,,1--称为向前差商，()()[]y h x u y h x u h,,21--+称为中心差商.也可用向前差商或中心差商代替一阶偏导数.2︒ x 轴与y 轴也可分别采用不同的步长h ,l ，即用直线族⎩⎨⎧====jhy y ihx x j i (i,j =0, ±1, ±2 , ) 作一个矩形网格.图14.72. 椭圆型方程的差分方法[五点格式] 考虑拉普拉斯方程的第一边值问题()()⎪⎪⎩⎪⎪⎨⎧=∈=∂∂+∂∂y x u D y x y ux u S ,,02222μ式中μ(x ,y )为定义在D 的边界S 上的已知函数.采用正方形网格，记u (x i ,y j )=u ij ，在节点(i ,j )上分别用差商 u u u h u u u h i j ij i j i j ij i j -+-+-+-+11211222,,,,,代替2222,yux u ∂∂∂∂，对应的差分方程为u u u h u u u hi j ij i j i j ij i j -+-+-++-+=112112220,,,, （1） 或u u u u u ij i j i j i j i j =+++-+-+141111,,,,即任一节点(i ,j )上u ij 的值等于周围相邻节点上解的值的算术平均，这种形式的差分方程称为五点格式，在边界节点上取()()()h j i ij S j i y x u ∈=,,**μ（2） 式中(x i *,y j *)是与节点(i ,j )最接近的S 上的点.于是得到了以所有内节点上的u ij 值为未知量的若干个线性代数方程，由于每一个节点都可列出一个方程，所以未知量的个数与方程的个数都等于节点的总数，于是，可用通常的方法（如高斯消去法）解此线性代数方程组，但当步长不很大时，用高斯消去法将会遇到很大困难，可用下面介绍的其他方法求解. 若h →0时，差分方程的解收敛于微分方程的解，则称差分方程为收敛的.在计算过程中，由于进行四则运算引起舍入误差，每一步计算的舍入误差都会影响以后的计算结果，如果这种影响所产生的计算偏差可以控制，而不至于随着计算次数的增加而无限增大，则称差分方程是稳定的.[迭代法解差分方程] 在五点格式的差分方程中，任意取一组初值{u ij }，只要求它们在边界节点(i ,j )上取以已知值μ(x i *,y j *)，然后用逐次逼近法（也称迭代法）解五点格式：()()()()()[]() ,2,1,0411,1,,1,11=+++=+-+-+n u u u u u n j i n j i n j i n j i n ij 逐次求出{u ij (n )}.当(i+1,j )，(i －1,j )，(i ,j －1)，(i ,j+1)中有一点是边界节点时，每次迭代时，都要在这一点上取最接近的边界点的值.当n →∞时，u ij (n )收敛于差分方程的解，因此n 充分大时，{u ij (n )}可作差分方程的近似解，迭代次数越多，近似解越接近差分方程的解.[用调节余数法求节点上解的近似值] 以差商代替Δu 时，用节点(i+1,j )，(i －1,j )，(i ,j+1)，(i ,j －1)上u 的近似值来表示u 在节点(i ,j )的值将产生的误差，称此误差为余数R ij ，即()()()()()ij j i j i j i j i j i R y x u h y x u h y x u y h x u y h x u =--+++-++,4,,,,设在(i ,j )上给u ij 以改变量δu ij ，从上式可见R ij 将减少4δu ij ，而其余含有u (x i ,y j )的差分方程中的余数将增加δu ij ，多次调整δu ij 的值就可将余数调整到许可的有效数字的范围内，这样可获得各节点上u (x ,y )的近似值.这种方法比较简单，特别在对称区域中计算更简捷.例 求Δu =0在内节点A ,B ,C ,D 上解的近似值.设在边界节点1,2,3,4上分别取值为1,2,3,4（图14.8） 解 记u (A )=u A ，点A ,B ,C ,D 的余数分别为－4u A + u B + u c +5=R A u A －4 u B + u D +7=R Bu A－4 u c + u D +3=R Cu B + u c －4u D +5=R D以边界节点的边值的算术平均值作为初次近似值，即u A (0)=u B (0)=u C (0)=u D (0)=2.5则相应的余数为：R A =0, R B =2, R C = －2, R D =0最大余数为±2.先用δu C =－0.5把R C 缩减为零，u C 相应地变为2，这时R A , R D 也同时缩减(－0.5)，新余数是R A =－0.5,R B =2,0=C R , R D =－0.5.类似地再变更δu B =0.5，从而 u B 变为3，则得新余数为0====D C B A R R R R .这样便可消去各节点的余数，于是u 在各节点的近似值为：u A =2.5, u B =3, u C =2, u D =2.5现将各次近似值及余数列表如下：次数调 整 值第n 次近似值及余数u A R A u B R B u C R C u D R D 0 1 2δu C = －0.5 δu B = 0.5 2.5 2.5 2.5 0 －0.5 0 2.5 2.5 3 2 2 0 2.5 2 2 －2 0 0 2.5 2.5 2.5 0 －0.5 0 结果近似值2.5322.5[解重调和方程的差分方法] 在矩形D (x 0≤x ≤x 0+a ,y 0≤y ≤y 0+a )中考虑重调和方程024*******=∂∂+∂∂∂+∂∂=yuy x u x u u ∆取步长h an=，引直线族图14.8⎩⎨⎧+=+=jh y y ihx x 00 (i , j = 0, 1, 2,, n ) 作成一个正方形网格.用差商代替偏导数()()()()()[]{()()()()[]()()()()[]}h y x u h y x u y h x u y h x u h y h x u h y h x u h y h x u h y h x u h y x u h y x u y h x u y h x u y x u 2,2,,2,2,,,,2,,,,8201,-+++-++---++-+-++++--+++-++=上式表明了以(x ,y )为中心时，u (x ,y )的函数值与周围各点函数值的关系，但对于邻近边界节点的点(x ,y )，如图14.9中的A ，就不能直接使用上式，此时将划分网格的直线族延伸，在延伸线上定出与边界距离为h 的点，称这些点为外邻边界节点，如图14.9以A 为中心时，点E ,C 为边界节点，点J ,K 为E ,C 的外邻边界节点，用下法补充定义外邻边界节点J 处函数的近似值u J ，便可应用上面的公式.1︒ 边界条件为()()()S P P x uP u SS ∈==21,μ∂∂μ 时，定义u J =u A －2μ2(E )h .2︒ 边界条件为()()()S P P xuP u SS ∈=∂∂=2221,μμ时，定义u J =2μ1(E )－u A －h 2μ2(E ). [其他与Δu 有关的网格] 1︒ 三角网格（图14.10(a )）取P 0(x ,y )为中心，它的周围6个邻近节点分别为：()()⎪⎪⎭⎫⎝⎛-+⎪⎪⎭⎫ ⎝⎛---⎪⎪⎭⎫⎝⎛+-⎪⎪⎭⎫⎝⎛+++h y h x P h y h x P y h x P h y h x P h y h x P y h x P 23,2,23,2,,23,223,2,,654321则 R u h u u u h i i +∆+∆=⎪⎭⎫⎝⎛-∑=226102161632式中u i =u (P i ), u 0=u (P 0)，R 表示余项. 2︒ 六角网格（图14.10(b )）取P 0(x ,y )为中心，它的三个邻近节点分别为图14.9()⎪⎪⎭⎫ ⎝⎛-+-⎪⎪⎭⎫ ⎝⎛++h y h x P y h x P h y h x P 23,2,,23,2321则 R u u u h i i +∆=⎪⎭⎫⎝⎛-∑=0312334.图14.103︒ 极坐标系中的网格（图14.10(c )）取P 0(r ,θ)为中心，它的四个邻近节点分别为()()()()l r P h r P l r P h r P ++--θθθθ,,,,,,4321而拉普拉斯方程01122222=∂∂+∂∂+∂∂=θ∆u r r u r ru u的相应的差分方程为()()()011221110222134222312=⎪⎭⎫ ⎝⎛+--++++u l r h u u rh u u l r u u h 3. 抛物型方程的差分方法 考虑热传导方程的边值问题()()()()()()⎪⎪⎪⎩⎪⎪⎪⎨⎧≥==<<=><<=∂∂-∂∂0,,,,00,0,0,0,021222t t t b u t t u bx x x u t b x x u a t u μμϕ 将[0,b ]分为n 等份，每段长为∆x bn=.引两族平行线（图14.11）图14.11x =x i =i ∆x (i =0,1,2,, n )y =y j =j ∆t (j =0,1,2,, ∆t 取值见后)作成一个长方形的网格，记u (x i ,t j )为u ij ，节点(x i ,t j )为(i ,j )，在节点(i ,j )上分别用(),2,1,1,,2,1Δ2,Δ2,1,11,=-=+---++j n i x u u u t u u ji ij j i ij j i 代替22,xut u ∂∂∂∂，于是边值问题化为差分方程()()()()⎪⎪⎪⎩⎪⎪⎪⎨⎧===-===-==+----++ ,2,1,0,Δ,Δ1,,2,1,Δ,2,1,0,1,,2,10Δ2Δ21002,1,121,j t j u t j u n i x i u j n i x u u u a tu u nj j i j i ij j i ijj i μμϕ 记()22x ta ∆∆=λ，差分方程可写成 () ,2,1,1,,2,121,1,11,=-=+-+=-++j n i u u u u ji ij j i j i λλλ （1）由此可按t 增加的方向逐排求解.在第0排上u i 0的值由初值ϕ(i ∆x )确定，j +1排u i ,j +1的值可由第j 排的三点(i +1,j ),(i ,j ),(i －1,j )上的值u i +1,j , u ij ,u i -1,j 确定，而u 0,j +1,u n ,j +1已由边界条件μ1((j +1)∆t )及μ2((j +1)∆t )给定，于是可逐排计算一切节点上的u ij 值.当ϕ(x ), μ1(x )和μ2(x )充分光滑，且λ≤12时，差分方程收敛而且稳定.所以利用差分方程（1）计算时，必须使λ≤12，即()22Δ21Δx at ≤.热传导方程还可用差分方程()0Δ2Δ21,11,1,121,=+---+-++++x u u u a t u u j i j i j i ij j i 代替，此时如已知前j 排u ij 的值，为求第j +1排的u i ,j +1 必须解包含n －1个未知量u u j n j 1111,,,,+-+ 的线性代数方程组，这种差分方程称为隐式格式的差分方程，前面所提的差分方程称为显式格式差分方程.隐式格式差分方程对任意的λ都是稳定的.4. 双曲型方程的差分方法 考虑弦振动方程的第一边值问题()()()()()()⎪⎪⎪⎩⎪⎪⎪⎨⎧≥==<<=∂∂=><<=∂∂-∂∂0,,,,00),()0,(,0,0,0,02122222t t t b u t t u b x x t x u x x u t b x x u a tu μμψϕ 用矩形网格，列出对应的差分方程：()()()()⎪⎪⎪⎩⎪⎪⎪⎨⎧===-=∆=∆-==-==+--+--+-+ ,2,1,0,Δ,Δ1,,2,1),(,Δ,2,1,1,,2,1,0Δ2)(Δ22100102,1,1221,1,j t j u t j u n i x i t u u x i u j n i x u u u a t u u u nj j i i i j i ij j i j i ij j i μμψϕ 记ω=a tx∆∆与上段一样，利用u u n 022,和在第0排及第1排的已知数值（初始条件）u i 0 , u i 1可计算u i 2，然后用已知的u i 1 , u i 2及u u n 033,可计算u i 3,类似地可确定一切节点上的u ij 值.当ϕ(x ),ψ(x ),μ1(x )和μ2(x )充分光滑，且ω≤1时，差分方程收敛且稳定，所以要取∆∆t ax ≤1.二、 变分方法1. 自共轭边值问题将§3定义的共轭微分算子的概念推广到一般方程.设D 是n E 中的有界区域，S 为其边界，在D 上考虑2k 阶线性微分方程()x f x x uaLu km mi i i ni m m i i n n n=∂∂≡∑∑==++201111 ∂ 的齐次边值问题()r j u l Sj ,,2,10==式中f (x )是D 内的已知函数，l j u 是线性微分算子. 将 ⎰DvLud Ω分部积分k 次得()⎰∑⎰⎪⎪⎭⎫ ⎝⎛+=Ω=S j j j D S v R u R v u vLu d ~,Λd k 1 式中Λ(u ,v )是一个D 上的积分，其被积函数包含u ,v 的k 阶导数；R j 和 R j 是定义在边界S 上的两个线性微分算子.再将Λ(u ,v )分部积分k 次得()()⎰∑⎰⎪⎪⎭⎫⎝⎛-Ω=Λ=S k j j j D S u R v R v uL v u d ~d ,1***式中L*是一个2k 阶的微分算子，称为L 的共轭微分算子.若L=L*，则称L 为自共轭微分算子.从上面可推出格林公式()()⎰∑⎰=-=Ω-Skj jjjjDS u R v R v R u R v uL vLu 1***d ~~d 如从l j u |S =l j v |S =0可推出在边界S 上()∑==-kj jjjju R v R v R u R 1**0~~ 则称l j u |S =0为自共轭边界条件.如果微分算子及边界条件都是自共轭的，则称相应的边值问题为自共轭边值问题，此时有()0d ][=Ω-⎰DuLv vLu每个边值问题对应于某希尔伯特空间H （例如L 2(D )，见第九章§7）中的一个算子A ，其定义域M A 是H 中一线性稠密集合，它由足够次连续可微且满足边界条件的函数组成，在M A 上，Au 的数值与Lu 的数值相同，从而求解边值问题化为解算子方程Au f = 的问题.设A 为定义在实的希尔伯特空间H 中的某线性稠密集合M A 上的线性算子.若对于M A 的任意非零元素,,v u 成立(Au ,v )=(u ,Av )则称A 为对称算子.若对任意非零元素u 成立()0,>u Au 则称A 为正算子.如成立更强的不等式(Au ,u )≥r ||u ||2 (r>0)则称A 为正定算子.此处(u ,v )表示希尔伯特空间的内积，||u ||2=(u ,u ). 2. 变分原理与广义解定理 设A 是正定算子，u 是方程Au =f 在M A 上的解的充分必要条件是： u 使泛函F (u )=(Au ,u )－2(f ,u )取极小值.上述将边值问题化为等价的求泛函极值问题的方法称为能量法.在算子的定义域不够大时，泛函F (u )的极值问题可能无解.不过对于正定算子，可以开拓集合M A ，使在开拓了的集合上，泛函的极值问题有解.为开拓M A ，在M A 上引进新的内积[u ,v ]=(Au ,v )，定义模||u ||2=[u ,u ]=(Au ,u )，在模||u ||的意义下，补充极限元素，得到一个新的完备希尔伯特空间H 0，在H 0上，泛函F (u )仍然有意义，而泛函的极值问题有解.但必须注意，此时使泛函F (u )取极小的元素u 0不一定属于M A ，因此它不一定在原来的意义下满足方程Au=f 及边界条件.称u 0为广义解. 3. 极小化序列与里兹方法在处理变分问题中，极小化序列起着重要的作用.考虑泛函F (u )=(Au ,u )－2(f ,u )以d 表示泛函的极小值.设在希尔伯特空间中存在一列元素{u n } (n =1,2 ,)，使()d u F n n =∞→lim则称{u n }为极小化序列.定理 若算子A 是正定的，则F (u )的每一个极小化序列既按H 空间的模也按H 0的模收敛于使泛函F (u )取极小的元素.这个定理不但指出利用极小化序列可求问题的解，而且提供一种近似解的求法，即把极小化序列中的每一个元素当作问题的近似解.设算子A 是正定的，构造极小化序列的里兹方法的主要步骤是：(1) 在线性集合M A 中选取H 0中完备的元素序列{ϕi } ， (i =1,2 ,) 并要求对任意的n ,ϕ1,ϕ2,…,ϕn 线性无关.称这样的元素为坐标元素.(2) 令u a n k k k n==∑ϕ1 ，其中a k 为待定系数.代入泛函F (u )，得自变量a 1,a 2,…,a n 的函数()()()∑∑==-=nj jjn k j kjkj n f a A a a u F 11,,2,ϕϕϕ(3) 为使函数F (u n )取极小，必须()()n j a u F jn ,,2,10 ==∂∂，从而求出a k (k =1,2,…,n ).序列{u n }即为极小化序列，u n 可作为问题的近似解. 4. 里兹方法在特征值问题上的应用 算子方程Au －λu =0的非零解λ称为算子A 的特征值，对应的非零解u 称为λ所对应的特征函数. 对线性算子A ，若存在常数K ，使对任何M A 的元素ϕ成立(A ϕ,ϕ)≥K ||ϕ||2则称A 为下有界算子，正定算子是下有界的（此时K =0）.记(A ϕ,ϕ)/||ϕ||2的下确界为d . 定理1 设A 为下有界对称算子，若存在不为零的元素ϕ0∈M A ，使()d A =200,ϕϕϕ则d 就是A 的最小特征值，ϕ0为对应的特征函数.于是求下有界对称算子的最小特征值问题化为变分问题，即在希尔伯特空间中求使泛函(A ϕ,ϕ)/||ϕ||2取极小的元素，或在||ϕ||=1的条件下求使泛函(A ϕ,ϕ)取极小的元素.定理2 设A 是下有界对称算子，λ1≤λ2≤…≤λn 是它的前n 个特征值，ϕ1,ϕ2,…,ϕn 是对应的标准正交特征函数，如果存在不为零的元素1+n ϕ，在附加条件(ϕ,ϕ)=1, (ϕ,ϕ1)=0, (ϕ,ϕ2)=0, …, (ϕ,ϕn )=0下使泛函(A ϕ,ϕ)取极小，则ϕn +1是算子A 的特征函数，对应的特征值()11,++=n n A ϕϕλ就是除λ1 ,，λn 外的最小的一个特征值.于是求第n +1个特征值就化为变分问题，即在附加条件(ϕ,ϕ)=1, (ϕ,ϕ1)=0, (ϕ,ϕ2)=0 ,, (ϕ,ϕn )=0 下求使泛函(A ϕ,ϕ)取极小的元素.为了利用里兹方法求特征值，在M A 中选取一列在H 0中完备的坐标元素序列{ϕi },(i =1,2 ,)， 令u a n k k k n==∑ϕ1，确定a k ，使在条件 (u n ,u n )=1下，(Au n ,u n )取极小，这个问题化为求n个变元a 1,a 2,…,a n 的函数()()∑==nm k m k k m n n a a A u Au 1,,,ϕϕ在条件()()∑===nm k m k m k n n a a u u 1,1,,ϕϕ下的极值问题，一般可用拉格朗日乘数法解（见第九章§3，t ），此时()()()()()()()()()()()()0,,,,,,,,,,,,11222121111111=------n n n n n n n n n n A A A A A A ϕϕλϕϕϕϕλϕϕϕϕλϕϕϕϕλϕϕϕϕλϕϕϕϕλϕϕ的最小的根即为特征值的近似值，如果将上式的根按大小排列，就依次得后面的特征值的近似值，但精确度较差. 对一般算子方程Au －λBu=0如果A 为下有界对称算子，B 为正定算子，则()()()()()()()()()()()()0,,,,,,,,,,,,11222121111111=------n n n n n n n n n n B A B A B A B A B A B A ϕϕλϕϕϕϕλϕϕϕϕλϕϕϕϕλϕϕϕϕλϕϕϕϕλϕϕ的根就是特征值的近似值. 5. 迦辽金方法用里兹方法解数学物理问题有很多限制，最主要的限制是要求算子正定，但很多问题不一定满足这个条件，迦辽金方法弥补了这个缺陷. 迦辽金方法的主要步骤是：(1) 在M A 中选取在空间H 中完备的元素序列{ϕi } (i =1,2 ,)，其中任意n 个元素线性无关，称{ϕi } (i =1,2,…)为坐标元素序列. (2) 把方程的近似解表示为u a n k k k n==∑ϕ1式中a k 是待定常数，把u n 代入方程Au=f 中的u ，两端与ϕj (j =1,2,…,n )求内积，得 a k 的n 个代数方程()()()n j f A a j n k j k k,,2,1,,1 ==∑=ϕϕϕ(3) 求出a k ，代回u n 的表达式，便得方程的近似解u n .在自共轭边值问题中，当算子是正定时，由迦辽金方法和里兹方法得到的关于a k 的代数方程组是相同的.。
偏微分-第三章
clear N=1000;dx=0.01; dt=0.001;c=dt/dx; x=linspace(0,3,300)'; u(1:300,1)=1; u(60:300,1)=0; h=plot(x,u(:,1),'linewidth',3); axis([0,3,-0.2,1.5]); set(h,'EraseMode','xor','MarkerSize',1) for k=2:N set(h,'XData',x,'YData',u(:,1)); drawnow; u(2:199,2)=u(2:199,1)+c/2*(u(1:198,1)u(3:200,1))+1/2*c*c*(u(3:200,1)-2*u(2:199,1)+u(1:198,1)); u(2:199,1)=u(2:199,2); end
偏微分课程课件5_双曲型方程的差分方法(II)
uvn j
uv t
n j
2
2
2uv t 2
n j
O(
3
)
uv n1 j
uv
n j
A
uv x
n j
2
2
A2
2uv n
x
2
j
O(
3)
用中心差商代替偏导数
uvn j
A
uvn uvn
j1
j1
2h
2
2
A2
2
x
uv n j
h2
O(
3
2h2
h2 )
舍去截断误差, 有LW差分格式.
1 2
a
nj （（u
n ）2
j 1
（u
n ）2）
j1
21
（unj 1）2
1 2
(unj 1 )2
（unj1）2 )
1 2
anj
( unj 1 )2
(
un j 1
)2
(4)用h乘上式两边并对 j 求和，记离散模
un
2
h
（unj）2h
||
un1
||h2 ||
un
||h2
1 2
a
n j
((unj 1
t
x
2u t 2
（ t
a(x)
u ） x
a(x)
x
(u ) t
= a(x) (a(x) u ) a(x) (a(x) u )
x
x
x x
25
代入Taylor展开式,于是有
u(
x
j
,
tn 1 )
u(x
j
,
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
2
2
2
2
1 4a(1 a)sin2 kh
2
a 0?
a 1,|G( ,k)| 1,Von Neumann 条件满足
条件稳定
7
a0
un1 j
u
n j
a
un j1
unj
0,
h
vn1 ((1 a ) aeikh )vn
| G(k, ) |2 (1 a a cos kh)2 a2 2 sin2 kh
第三章 双曲型方程定解问题 的有限差分法
3.1 一阶线性常系数双曲型方程 3.2 一阶线性常系数双曲型方程组 3.3 变系数双曲型方程及方程组
1
(一) 一阶线性常系数双曲型方程
u t
a
u x
0,
x R,
t0
u(x,0) u0(x), x R.
u定义在xt平面上的一个区域内.
下面考察方程解u在定义域内直线x=at+C上的变化规律.
a 0： 1 a 0
a 0:
jh an jh不收敛
微分方程解的依赖区域不属于差分方程解的依赖区域
21
Lax-Wendroff格式的C.F.L条件
x jn x j atn x jn
jh nh jh an jh nh
nh an nh |a| 1,其中 .
h
B:(x j1 ,tn ) ((j 1)h,n ) C:(x j ,tn ) (jh,n ) D:(x j1 ,tn ) ((j 1)h,n ) P:(x j ,tn1 ) (jh,(n 1) )
30
过P做特征线 x-at C（, C jh-a(n 1)）, 与BC相交于Q(jh-a ,n )
u(x j1,tn )
u(x j1,tn )
O(h2 ),
2u n
x
2
j
1 h2
u(x j1,tn )
2u(x j ,tn )
u(x j1,tn )
O(h2 )
16
un1 j
unj
a 2
h
un j1
un j1
a2 2
2
h2
(u nj 1
2unj
un j1
)
G( ,k) 1 a2 2( cos kh 1) ia sin kh | G( ,k)|2 1 4a2 2(1 a2 2 )sin4 kh
(a|a|)(unj
un j 1
)
1 2
(a|a|)(unj1
unj
)
,
un1 j
unj
1 2
a (unj1
un j 1
)
1 2
|a|
(unj 1
2unj
un j 1
).
11
2 Lax-Friedrichs格式
中心差分格式
un1 j
unj
a
un j1
un j -1
0
2h
1
1
T (x j , tn ) t u(x j , tn ) a h ox u(x j , tn )
24
（2）中心差分格式
un1 j
unj
a
un j1
un j 1
0，
2h
G( ,k ) 1 a (eikh eikh )
2
1 ai( sin kh)
|G( ,k)|2 1 a2 2sin2kh，
sin kh 0, |G( ,k)| 1
不稳定，C.F.L.条件仍为 | a| 1，
=1+4a(1+a )sin2 kh
2
a 0,| a | 1
1 a 0
| G( , k) | 1.
条件稳定
a 0? 8
证明：
课堂练习
a 0：
un1 j
unj
a
un j1
unj
0,
h
a
0:
un1 j
unj
a unj
un j 1
0,
h
绝对不稳定
9
a 0:
un1 j
unj
a
unj
P
AB C D
Q
u(P) u(Q)
1.由B、C线性插值求u(Q): a 0时的迎风格式 2.由B、D线性插值求u(Q): Lax-Fridriches格式 3.由B、C、D二次插值求u(Q): Lax-Wendroff格式 4.由A、B、C二次插值求u(Q): Beam-Warming格式
B，C点插值
2
2
cos kh ia sin kh,
| G( , k) |2 cos2 kh a2 2 sin1 则稳定
14
Lax-Friedrichs格式可以不考虑特征线走向， 但截断误差比迎风格式的截断误差大。
15
3. Lax-Wendroff格式（2阶精度）
1 2
(unj 1
un j 1
)
a
un j1
un j1
0
2h
T (x,t) O( h2 ) O( h2 )
取unj v en ikjh
v n1
1 2
(e ikh
eikh )
a
2
(e ikh
eikh )
vn
G( , k ) 1 (eikh eikh ) a (eikh eikh )
19
x jn x j atn x j
jh nh jh an jh
a 0 0 a 1
其中 .
h
a 0 jh an jh 不收敛
n
P
D
D' C
D'
20
右偏心格式C.F.L.条件
un1 j
unj
a
un j1
unj
0
h
x j x j atn x jn
jh jh an jh nh
(x ( xk1
xk )( x xk1 ) xk )( xk 1 xk 1 )
yk
(x ( xk
xk 1 )( x xk 1 ) xk 1 )( xk xk 1 )
yk 1
(x ( xk1
xk1 )( x xk xk 1 )( xk 1
) xk
)
.
29
P
B CD
设B,C,D处的值已知，下面来确定点P处的值u(P)
P点数值解依赖于DC内节点上的函数值---依赖区域 18
P点数值解依赖于DC内节点上的函数值---依赖区域
t
n
P
D
D' C x
jn
j
C.F.L.条件,
P 点微分方程解依赖区域 D' 应在
差分方程的依赖区域端点构成的区间DC内, 否则 unj 与u( xj , tn )
没有关系。即差分格式的依赖区域应该包含微分方程解的依赖区域。
a unj
un j 1
0,
h
un1 j
unj
a
un j1
unj
0,
h
a 0, a 0.
n1 n
j1 j j1
a0
n1
n j1 j
a0
j1
5
a0
un1 j
unj
a unj
un j 1
0
h
un1 j
unj
a
(u
n j
un j1
)
取unj v neikjh
G( ,k) 1 a(1 eikh ) 1 a(1 cos kh) ai( sin kh)
2
|G( ,k)| 1 |a| 1
T (x,t ) O( 2 h2 )
17
4. Courant-Friedrichs-Lewy条件（C.F.L.条件）
左偏心格式
un1 j
unj
a unj
un j1
0
h
t
P( j, k ) (xj ,tk )
( j,k)
D( xjk ,0)
C( xj ,0) x
22
C.F.L.条件 差分格式的依赖区域包含微分方程的依赖区域
C.F.L.条件是格式收敛的必要条件.
23
（1）Lax-Wendroff格式：
G(k, ) 1 a2 2( cos kh-1)-ia sin kh |G(k, )| 1 |a| 1
C.F.L.条件|a| 1也是稳定性条件
由Lax等价定理，C.F.L.条件是收敛充分条件。
un j1
0,
h
a 0:
un1 j
unj
a
un j1
unj
0
h
差分格式与微分方程的特征线走向一致，条件稳定。
n1 n
j1 j j1
a0
n1
n j1 j
a0
j1
10
a
0:
un1 j
unj
a
unj
un j 1
0,
h
a 0:
un1 j
unj
a
un j1
unj
0
h
迎风格式统一形式
un1 j
unj
1 2
C.F.L.条件下不收敛（非充分条件）
25
课堂练习
1. 试给出一阶双曲型方程左偏心格式、右偏 心格式、中心差分格式的C.F.L.条件。
右偏心格式C.F.L.条件
a 0： x j x j atn x jn 1 a 0
a 0:
jh an jh不收敛
n1
n j1 j j1
a0
左偏心格式C.F.L.条件
特征线与t=0交于点 ( x0 at0 ,0). u( x0 , t0 ) u( x0 at0 ,0) u0( x0 at0 ).