(高等数学) 偏微分方程.doc_第1页
(高等数学) 偏微分方程.doc_第2页
(高等数学) 偏微分方程.doc_第3页
(高等数学) 偏微分方程.doc_第4页
(高等数学) 偏微分方程.doc_第5页
已阅读5页,还剩36页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

第十四章 偏微分方程 物理、力学、工程技术和其他自然科学经常提出大量的偏微分方程问题.由于实践的需要和一些数学学科(如泛函分析,计算技术)的发展,促进了偏微分方程理论的发展,使它形成一门内容十分丰富的数学学科. 本章主要介绍一阶偏微分方程、线性方程组及二阶线性偏微分方程的理论.在二阶方程中,叙述了极值原理、能量积分及惟一性定理.阐明了一些解的性质和物理意义,介绍典型椭圆型、双曲型、抛物型方程的常用解法:分离变量法,基本解,格林方法,黎曼方法,势位方法及积分变换法.最后,扼要地介绍了有实用意义的数值解法:差分方法和变分方法.1 偏微分方程的一般概念与定解问题 偏微分方程及其阶数 一个包含未知函数的偏导数的等式称为偏微分方程.如果等式不止一个,就称为偏微分方程组.出现在方程或方程组中的最高阶偏导数的阶数称为方程或方程组的阶数. 方程的解与积分曲面 设函数u在区域D内具有方程中所出现的各阶的连续偏导数,如果将u代入方程后,能使它在区域D内成为恒等式,就称u为方程在区域D中的解,或称正规解. 在n+1维空间中是一曲面,称它为方程的积分曲面. 齐次线性偏微分方程与非齐次线性偏微分方程 对于未知函数和它的各阶偏导数都是线性的方程称为线性偏微分方程.如就是线性方程.在线性方程中,不含未知函数及其偏导数的项称为自由项,如上式的f(x,y).若自由项不为零,称方程为非齐次的.若自由项为零,则称方程为齐次的. 拟线性方程与半线性方程 如果一个方程,对于未知函数的最高阶偏导数是线性的,称它为拟线性方程.如就是拟线性方程,在拟线性方程中,由最高阶偏导数所组成的部分称为方程的主部.上面方程的主部为如果方程的主部的各项系数不含未知函数,就称它为半线性方程.如就是半线性方程. 非线性方程 不是线性也不是拟线性的方程称为非线性方程.如就是一阶非线性偏微分方程. 定解条件 给定一个方程,一般只能描写某种运动的一般规律,还不能确定具体的运动状态,所以把这个方程称为泛定方程.如果附加一些条件(如已知开始运动的情况或在边界上受到外界的约束)后,就能完全确定具体运动状态,称这样的条件为定解条件.表示开始情况的附加条件称为初始条件,表示在边界上受到约束的条件称为边界条件. 定解问题 给定了泛定方程(在区域D内)和相应的定解条件的数学物理问题称为定解问题.根据不同定解条件,定解问题分为三类. 1 初值问题 只有初始条件而没有边界条件的定解问题称为初值问题或柯西问题. 2 边值问题 只有边值条件而没有初始条件的定解问题称为边值问题. 3 混合问题 既有边界条件也有初始条件的定解问题称为混合问题(有时也称为边值问题). 定解问题的解 设函数u在区域D内满足泛定方程,当点从区域D内趋于给出初值的超平面或趋于给出边界条件的边界曲面时,定解条件中所要求的u及它的导数的极限处处存在而且满足相应的定解条件,就称u为定解问题的解. 解的稳定性 如果定解条件的微小变化只引起定解问题的解在整个定义域中的微小变化,也就是解对定解条件存在着连续依赖关系,那末称定解问题的解是稳定的.定解问题的适定性 如果定解问题的解存在与惟一并且关于定解条件是稳定的,就说定解问题的提法是适定的.2 一阶偏微分方程一、 柯西-柯娃列夫斯卡娅定理 一阶偏微分方程的通解 一阶偏微分方程的一般形式 在有些书中写作是或,其中如解出p1,可得:p1 = f (x1 , x2 , xn , u , p2 , pn ) 当方程的解包含某些“任意元素”(指函数),如果适当选取“任意元素”时,可得方程的任意解(某些“奇异解”除外),则称这样的解为通解. 在偏微分方程的研究中,重点在于确定方程在一些附加条件(即定解条件)下的解,而不在于求通解. 一阶方程的柯西问题称为柯西问题,式中为已知函数,对柯西问题有如下的存在惟一性定理. 柯西柯娃列夫斯卡娅定理 设 f ( x1 , x2 xn , u , p2 pn ) 在点 ( x10 , x20 xn0 , u0 , p20 pn0 ) 的某一邻域内解析,而在点( x20 xn0 ) 的某邻域内解析,则柯西问题在点 ( x10 xn0 ) 的某一邻域内存在着惟一的解析解. 这个定理应用的局限性较大,因它要求f及初始条件都是解析函数,一般的定解问题未必能满足这种条件. 对高阶方程也有类似定理.二、 一阶线性方程 1一阶齐次线性方程 特征方程特征曲线初积分(首次积分) 给定一阶齐次线性方程 (1)式中ai为连续可微函数,在所考虑的区域内的每一点不同时为零(下同).方程组 ( i = 1,2n ) 或 (2)称为一阶齐次线性偏微分方程的特征方程.如果曲线l: xi = xi (t) ( i=1,2n )满足特征方程(2),就称曲线l为一阶齐次线性方程的特征曲线. 如果函数 ( x1 , x2 xn )在特征曲线上等于常数,即 ( x1(t) , x2(t) xn(t) ) = c就称函数 ( x1, x2 xn )为特征方程(2)的初积分(首次积分). 齐次方程的通解 1o 连续可微函数u = ( x1, x2 xn ) 是齐次线性方程(1)的解的充分必要条件是: ( x1, x2 xn )是这个方程的特征方程的初积分. 2o 设i ( x1 , x2 xn ) ( i = 1,2 n) 是特征方程(2)在区域D上连续可微而且相互独立的初积分(因此在D内的每一点,矩阵的秩为n) ,则u = ( 1 ( x1 , x2 xn ) n-1 ( x1 , x2 xn ) )是一阶齐次线性方程(1)的通解,其中为n个变量的任意连续可微函数. 柯西问题 考虑方程的柯西问题式中 ( x2 xn )为已知的连续可微函数. 设i ( x1 , x2 xn ) ( i = 1,2 n) 为特征方程的任意n个相互独立的初积分,引入参变量 (),从方程组解出x2 xn 得则柯西问题的解为u = ( 2 ( 1 , 2 n-1 ) n ( 1 , 2 n-1 ) ) 2非齐次线性方程 它的求解方法与拟线性方程相同.三、 一阶拟线性方程 一阶拟线性方程为其中ai及R为x1 , x2 xn , u的连续可微函数且不同时为零. 一阶拟线性方程的求解和它的特征方程或为原拟线性方程的特征方程.如果曲线l: xi = xi (t) ( i=1,2n ) , u = u(t) 满足特征方程,则称它为拟线性方程的特征曲线. 设 i ( x1 xn,u ) ( i = 1,2 n) 为特征方程的n个相互独立的初积分,那末对于任何连续可微函数, ( 1 ( x1 xn , u) , 2 ( x1 xn , u) n ( x1 xn , u) ) = 0都是拟线性方程的隐式解. 柯西问题 考虑方程的柯西问题为已知的连续可微函数. 设 1 ( x1 , x2 xn , u) n ( x1 , x2 xn , u) 为特征方程的n个相互独立的初积分,引入参变量 , 从解出 x2 xn , u 则由给出柯西问题的隐式解.四、 一阶非线性方程 完全解通解奇异解 一阶非线性方程的一般形式为 若一阶偏微分方程的解包含任意n个独立的常数,则称这样的解为完全解(全积分). 若V ( x1, x2 xn , u , c1 , c2cn ) = 0为方程的完全解,从消去ci ,若得一个解,则称它为方程的奇异解(奇积分). 以两个独立变量为例说明完全解与通解、奇异解的关系,设方程有完全解V (x,y,z,a,b)=0 ( a,b为任意常数),则方程等价于从方程组消去a,b所得的方程. 利用常数变易法把a,b看作x, y的函数,将V (x,y,z,a,b)=0求关于x, y的偏导数,得那末与V=0联立可确定a,b.有三种情况: 1 ,将其与V(x,y,z,a,b)=0联立可确定不含任意常数的奇异解. 2 如,即回到完全解. 3 当时,必有,这时,如果不属于情形2 ,则a与b存在函数关系:b=(a),这里为任意可微函数,并从方程V(x,y,z,a,b)=0和消去a,b,可确定方程的通解. 定理 偏微分方程的任何解包含在完全解内或通解内或奇异解内. 特征方程特征带特征曲线初积分 在一阶非线性方程:中,设F对所有变量的二阶偏导数存在且连续,称或为非线性方程的特征方程.设特征方程的解为xi=xi(t), u=u(t), pi=pi(t) (i=1,2,n)称它为非线性方程的特征带.在x1,x2xn,u空间的曲线xi=xi(t), u=u(t) (i=1,2,n)称为非线性方程的特征曲线.如果函数在特征方程的任一解xi=xi(t) (i=1,2n), u=u(t), pi=pi(t) (i=1,2n)上等于常数,即那末函数称为特征方程的初积分. 求完全解的拉格朗日恰比方法 考虑两个变量的情况. 对于方程F(x,y,z,p,q)=0,选择使雅可比式的一个初积分G(x,y,z,p,q).解方程组 (a为任意常数)得p(x,y,z,a)及q(x,y,z,a).则方程dz=pdx+qdy的通解V(x,y,z,a,b)=0(b是积分dz=pdx+qdy出现的任意常数)就是方程F(x,y,z,p,q)=0的完全解. 例 求方程的完全解. 解 方程的特征方程为这里成立所以特征方程的一个初积分为z2p2 x2 . 解方程组 (a为任意常数)得 积分微分方程得完全解 (b为任意常数) 某些容易求完全解的方程 1 仅含p,q的方程F(p,q)=0 G=p是特征方程的一个初积分.从F(p,q)=0与p=a(a为任意常数)得q=(a),积分dz=adx+(a)dy得完全解z=ax+(a)y+b (b为任意常数) 2 不显含x,y的方程F(z,p,q)=0 特征方程为因此qdp-pdq=0,显然为一个初积分,由F(z,p,q)=0,q=pa(a为任意常数)解得p=(z,a).于是由dz=(z,a)dx+a(z,a)dy得 (b为任意常数)可确定完全解. 3 变量分离形式的方程特征方程为可取初积分Gi=fi(xi,pi) , (i=1,2n).从fi(xi,pi)=ai (i=1,2n)解出pi=i(xi,ai)得完全解式中ai,b为任意常数,且. 克莱罗方程 方程称为克莱罗方程,其完全解为对ci微分得 (i=1,2,n)与完全解的表达式联立消去ci即得奇异解. 例 求方程zxpyqpq=0的完全解和奇异解. 解 这是克莱罗方程,它的完全解是z=ax+by+ab 对a,b微分,得x=b,y=a,消去a,b得奇异解z=xy 发甫方程 方程P(x,y,z)dx+Q(x,y,z)dy+R(x,y,z)dz=0 (1)称为发甫方程,如果P,Q,R二次连续可微并满足适当条件,那末方程可积分.如果可积分成一关系式时,则称它为完全可积. 1 方程完全可积的充分必要条件 当且仅当P,Q,R满足条件 (2)时,存在一个积分因子(x,y,z),使dU1=(Pdx+Qdy+Rdz)从而方程的通解为U1(x,y,z)=c 特别,当时,存在一个函数U(x,y,z)满足从而 dU=Pdx+Qdy+Rdz 所以方程的通解为U(x,y,z)=c 所以完全可积的发甫方程的通解是一单参数的曲面族. 定理 设对于发甫方程(1)在某区域D上的完全可积条件(2)成立,则对D内任一点M(x,y,z)一定有方程的积分曲面通过,而且只有一个这样的积分曲面通过. 2 方程积分曲面的求法 设完全可积条件(2)成立.为了构造积分曲面,把z看成x,y的函数(设R(x,y,z)0),于是原方程化为由此得方程组发甫方程(1)与此方程组等价. 把方程(3)中的y看成参变量,积分后得一个含有常数的通解然后用未知函数代替常数,将代入方程(4),在完全可积的条件下,可得的一个常微分方程,其通解为c为任意常数,代回中即得发甫方程的积分曲面z=(x,y,(y,c) 由于发甫方程关于x,y,z的对称性,在上面的讨论中,也可把x或y看成未知函数,得到同样的结果. 例 求方程yzdx+2xzdy+xydz=0的积分曲面族. 解 容易验证完全可积条件成立,显然存在一个积分因子,用它乘原方程得积分后得积分曲面族xy2z=c 也可把方程化为等价的方程组把y看成参变量,积分得通解用未知函数代替,将代入方程得积分后有所以原方程的积分曲面族是xy2z=c五、 一阶线性微分方程组 一阶线性偏微分方程组的一般形式 两个自变量的一阶线性方程组的形式是或 (1)其中Aij,Bij,Cij,Fi,aij,bij,fi是(x,t)的充分光滑函数. 特征方程特征方向特征曲线称为方程组(1)的特征方程.在点(x,t)满足特征方程的方向称为该点的特征方向.如果一条曲线l,它上面的每一点的切线方向都和这点的特征方向一致,那末称曲线l为特征曲线. 狭义双曲型方程与椭圆型方程 如果区域D内的每一点都存在n个不同的实的特征方向,那末称方程组在D内为狭义双曲型的. 如果区域D内的每一点没有一个实的特征方向,那末称方程组在D内为椭圆型的. 狭义双曲型方程组的柯西问题 1 化方程组为标准形式对角型 因为det(aij-ij)=0有n个不同的实根1(x,t) n(x,t),不妨设那末常微分方程的积分曲线li (i=1,2,n)就是方程组(1)的特征曲线. 方程的非零解(k(1) k(n)称为对应于特征方向k的特征矢量. 作变换可将方程组化为标准形式对角型 所以狭义双曲型方程组可化为对角型,而一般的线性微分方程组(1)如在区域D内通过未知函数的实系数可逆线性变换可化为对角型的话,(此时不一定要求 i都不相同),就称这样的微分方程组在D内为双曲型的. 2 对角型方程组的柯西问题 考虑对角型方程组的柯西问题i(x)是a,b上的连续可微函数.设ij,i,i在区域D内连续可微,在D内可得相应的积分方程组式中为第i条特征曲线li上点(x,t)与点(xi,0)之间的一段,(xi,0)为li与x轴上a,b的交点.上式可以更确切地写为 (i=1,2n)式中xi=xi(x,t,t)为过点(x,t)的第i条特征曲线,利用逐次逼近法可解此积分方程.为此令序列vi(k) (k=0,1,2)一致收敛于积分方程的连续可微解vi(x,t) (i=1,2n),这个vi(x,t)也就是对角型方程组的柯西问题的解. 设在区域D内对角型方程组的柯西问题的解存在,那末解与初值有下面的关系: (i) 依赖区间:过D中任意点M(x,t)作特征曲线l1,ln,交x轴于B,A,称区间A,B为M点的依赖区间(图14.1(a),解在M点的值由区间A,B的初值确定而与A,B外的初值无关. (ii) 决定区域:过点A,B分别作特征曲线ln,l1,称ln,l1 与区间A,B围成的区域D1为区间A,B的决定区域(图14.1(b),在区域D1中解的值完全由A,B上的初值决定. (iii) 影响区域:过点A,B分别作特征曲线l1,ln,称l1,ln与A,B围成的区域D2为区间A,B的影响区域(图14.1(c).特别当区间A,B缩为一点A时,A点的影响区域为D3(图14.1(d).在区域D2中解的值受A,B上的初值影响,而在区域D2外的解的值则不受A,B上的初值影响.图14.1 线性双曲型方程组的边值问题 以下列线性方程组来说明: (1) 1 第一边值问题(广义柯西问题) 设在平面(x,t)上给定曲线段,它处处不与特征方向相切.过A,B分别引最左和最右的特征曲线l1及l2.要求函数u(x,t),v(x,t)在,l1及l2围成的闭区域上满足方程组,且在上取给定的函数值(图14.2(a). 2 第二边值问题(古沙问题) 设l1是过P点的第一族特征线,l2是第二族特征线,在l1的一段PA上给定v(x,t)的数值,在l2的一段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.3 图14.2 边值问题的近似解特征线法 以上定解问题,可用逐步逼近法求解,也可用特征线法求解的近似值.以第一边值问题为例说明. 在曲线AB上取n个分点A1,A2, An,并记A为A0,B为An+1,过A0按A0的第二特征方向作直线与过A1按A1的第一特征方向作直线相交于B0;过A1按A1第二特征方向作直线与过A2按A2的第一特征方向作直线相交于B1最后得到Bn(图14.3).用如下的近似公式来确定方程组(1)的解u(x,t),v(x,t)在Bi (i=0,1,2,n)的数值:于是在一个三角形网格的节点上得到u,v的数值.再经过适当的插值,当n相当大,Ai、Ai+1的距离相当小时,就得到所提问题的足够近似的解. 特殊形式的拟线性方程组可化约系统 一般的拟线性方程组的问题比较复杂,目前研究的结果不多,下面介绍一类特殊形式的拟线性方程组可化约系统.如果方程组中所有的系数只是u,v的函数,称它为可化约系统. 考虑满足条件的方程组的解u=u(x,t),v=v(x,t).x,t可以表示成u,v的函数,且原方程化为这是关于自变量u,v的线性方程组.这样就把求拟线性方程组满足的解,化为解线性方程组的问题.而此线性方程组满足条件 的解,在(x,t)平面上的象即为原来拟线性方程组的解.3 二阶偏微分方程一、 二阶偏微分方程的分类、标准形式与特征方程 考虑二阶偏微分方程 (1)式中aij(x)=aij(x1,x2,xn)为x1,x2,xn的已知函数. 特征方程特征方向特征曲面特征平面特征锥面代数方程称为二阶方程(1)的特征方程;这里a1,a2,an是某些参数,且有.如果点x=(x1,x2,xn)满足特征方程,即则过x的平面的法线方向l:(a1,a2,an)称为二阶方程的特征方向;如果一个(n)维曲面,其每点的法线方向都是特征方向,则称此曲面为特征曲面;过一点的(n)维平面,如其法线方向为特征方向,则称这个平面为特征平面,在一点由特征平面的包络组成的锥面称为特征锥面. n个自变量方程的分类与标准形式 在点P(x1,x2,xn),根据二次型 (ai为参量)的特征根的符号,可将方程分为四类: (i) 特征根同号,都不为零,称方程在点P为椭圆型. (ii) 特征根都不为零,有n个具有同一种符号 ,余下一个符号相反,称方程在点P为双曲型. (iii) 特征根都不为零,有个具有同一种符号(nm1),其余m个具有另一种符号,称方程在点P为超双曲型. (iv) 特征根至少有一个是零,称方程在点P为抛物型. 若在区域D内每一点方程为椭圆型,双曲型或抛物型,则分别称方程在区域D内是椭圆型、双曲型或抛物型. 在点P作自变量的线性变换可将方程化为标准形式: 椭圆型: 双曲型: 超双曲型: 抛物型:式中为不包含二阶导数的项. 两个自变量方程的分类与标准形式 方程的一般形式为 (2) a11,a12,a22为x,y的二次连续可微函数,不同时为零. 方程a11dy2a12dxdy+a22dx2=0称为方程(2)的特征方程.特征方程的积分曲线称为二阶方程(2)的特征曲线. 在某点P(x0,y0)的邻域D内,根据=a122a11a12的符号将方程分类: 当0时,方程为双曲型; 当=0时,方程为抛物型; 当0,存在两族实特征曲线,作变换,和方程化为标准形式或 (ii) 抛物型: 因=0,只存在一族实的特征曲线,取二次连续可微函数,使,作变换,方程化为标准形式(iii) 椭圆型:因0,对任意和任意的ai有 定理1 设u(x)为D内椭圆型方程的解,它在D内二次连续可微,在上连续,且不是常数,如f(x)0(或f(x)0),则u(x)不能在D的内点取非正最小值(或非负最大值). 如果过边界S上的任一点P都可作一球,使它在P点与S相切且完全包含在区域D内,则有 定理2 设u(x)为椭圆型方程在D内二次连续可微,在上连续可微的解,且不是常数,并设f(x)0(或f(x)0).若u(x)在边界S上某点M处取非正最小值(或非负最大值),只要外法向导数在点M存在,则 (或) 2 定解问题 (i)第一边值问题(狄利克莱问题) (S) (ii)第二边值问题(诺伊曼问题) (S)其中 N为S的外法线方向. (iii)第三边值问题(混合问题) (S)a(),b(),()在S上连续,N是S的外法线方向,a()0,b()0,且a2()+b2()0. 3 解的惟一性问题 设c(x)及b()不同时恒等于零,如果定解问题Lu=f,lu=的解存在,则是惟一的,设c(x)及b()都恒等于零,如果定解问题Lu=f,lu=的解存在,则除相差一个常数外,解是惟一的. 抛物型方程的极值原理与解的惟一性定理 设为柱体,在柱体内部 考虑抛物型方程式中aij(x,t),bi(x,t),c(x,t),f(x,t)在上连续,且正定. 1 强极值原理 设u(x,t)为抛物型方程Lu=f(x,t)在D(0,T)内连续可微在上连续的解.并设f(x)=0,若u(x,t)在D(0,T的某点(x0,t0)取非负的最大值,即则对任意满足下列条件的点P(x,t),都有u(x,t)=m:点P(x,t)满足tt0 ,且可用完全在D(0,T 内的连续曲线x=x(t)与点(x0,t0)相连. 如在的侧边界:S0,T上(S是D的边界)任一点P都可作一球,使它在P点与相切且完全在D(0,T)内,则有 定理 设u(x,t)在上连续,在D(0,T内满足抛物型方程Lu=f,且不是常数,设 f0,若u(x,t)在上某点M处取非正最小值,只要外法向导数在点M存在,则 2 柯西问题与混合问题 柯西问题的初值条件是 混合问题按下列的定解条件分别称为 (i)第一边值问题:,; (ii)线性边值问题:,其中N为的外法线方向为已知函数,a0,b0,a2+b20 3 解的惟一性定理 如果抛物型方程Lu=f的混合问题的解存在,那末它是惟一的.如果柯西问题存在有界的解,那末在有界函数类中,解是惟一的. 波动方程的能量积分与解的惟一性定理 1 波动方程的柯西问题与混合问题 设波动方程为 柯西问题的初值条件是 如果在有界区域Q:D(0,T中考虑波动方程,记的侧边界为,则混合问题的定解条件是(i) 第一边值问题(ii) 第二边值问题(iii) 第三边值问题式中N为 的外法线方向,(x),(x)为D上的已知函数,(,t)及(,t) (,0tT)为定义在上的已知函数,(,t)0. 2 解的惟一性定理 波动方程的混合问题与柯西问题的解如果存在必定惟一. 惟一性定理可用下面能量积分证明. 3 能量积分 积分* 是的简写,下同.称为波动方程的能量积分. 满足齐次波动方程及u|=0(或)的函数u(x,t)成立: 能量守恒原理 E(t)=E(0). 能量不等式 式中 满足齐次波动方程及的函数,在上面能量不等式E(t)中增加一项,上面关系仍成立. 对于柯西问题,在特征锥 (R为大于零的常数)中考虑齐次波动方程的解u,记特征锥与t=t0的截面为,关于能量积分成立下面的能量不等式式中 t是 t=常数的超平面与以为上底所作的柱体(母线平行于Ot轴)的交截面.三、 三种典型方程 1.波动方程 研究下面形式的波动方程式中f(x,y,z,t)为已知函数. 许多物体的运动规律可用波动方程来描述.如弦振动可用一维波动方程描述;膜的振动可用二维波动方程描述;声波和电磁波的振荡可用三维波动方程描述. 齐次方程柯西问题的解 设齐次波动方程的柯西问题满足下面初始条件:并设三次连续可微,二次连续可微,那末解u的表达式分别为 1 三维(克希霍夫公式)式中Sat表示球面:,dS表示球面的面积元素. 2 二维(泊松公式)式中Kat表示圆: . 3 一维(达兰贝尔公式)利用降维法可从高维的解推得低维的解. 非齐次方程柯西问题的解 非齐次波动方程柯西问题的解等于上面齐次方程柯西问题的解添加一项所谓推迟势 .1 三维式中积分区域是以(x,y,z)为中心,at为半径的球体, 2 二维式中 . 一维 解的物理意义 波动方程解的表达式具有明确的物理意义. 1 波的传播 以弦振动为例,在达兰贝尔公式中,形如(x-at)的解描写了弦振动以常速度a向右传播,称(x-at)为右传播波,(x+at)为左传播波,a为波速. 2 依赖区间 过点P(x,t)作两条特征线x=c1,x+at=c2交x轴于x1,x2,则区间x1,x2称为点P的依赖区间,由达兰贝尔公式可见解在P点的值只与x1,x2上初始条件有关,而与区间外(x),(x)的值无关. 3 决定区域 过x轴上两点x1,x2(x10)称为x1,x2的决定区域(图14.4(a),在区域中解的数值由x1,x2上的初始条件完全决定.任意改变初始条件在x1,x2外的数值,解在此区域中不会发生任何变化.图14.4 4影响区域 过x轴上两点分别作特征线x=x1,x=x2+at称区域 x1-atxx2 (t0)为x1,x2的影响区域(图14.4(b).在此区域中,解的数值受到x1,x2上初始条件的影响,而在此区域外,解的值不受x1,x2上的初始条件影响,当区域x1,x2缩为一点x0时,点x0的影响区域为x轴上区间(图14.4(c)x0xx0+at (t0) 对二维波动方程,点(x0,y0,t0)的依赖区域为t=0上的圆.(xx0)2+(yy0)2a2t02在t=0上圆(xx0)2+(yy0)2a2t02的决定区域是以(x0,y0,t0)为顶点的圆锥体区域(图14.5(a).(xx0)2+(yy0)2a2(tt0)2 (tt0)初始平面t=0上一点(x0,y0,0)的影响区域为圆锥体(图14.5(b).(xx0)2+(yy0)2a2t2 (t0) (1)初始平面t=0上某一区域的影响区域,就是由此区域上每一点所作的圆锥体(1)的包络面所围成的区域.图14.5 对三维波动方程,点(x0,y0,z0,t0)的依赖区域为t=0上的球面(xx0)2+(yy0)2+(zz0)2=a2t02初始平面t=0上的球体(xx0)2+(yy0)2+(zz0)2a2t02的决定区域是以它为底,以(x0,y0,z0,t0)为顶点的圆锥体区域(xx0)2+(yy0)2+(zz0)2a2(tt0)2 (tt0)在初始平面t=0上点(x0,y0,z0,0)的影响区域为锥面(xx0)2+(yy0)2+(zz0)2=a2t2 (t0) (2)初始平面上某一区域的影响区域就是它上面的每一点所作的锥面(2)的包络面围成的区域. 二维与三维波的传播存在着下述本质区别. 5惠更斯原理 对三维波动方程,点(x0,y0,z0,0)的影响区域为(xx0)2+(yy0)2+(zz0)2a2t2 (t0)若在某一有界区域有一个初始扰动,在时刻t受到此初始扰动的影响区域就是所有以点为中心,以at为半径的球面全体,当t足够大时,这种球面族有内外包络面,称外包络面为传播波的前阵面,内包络面为后阵面.前阵面以外的部分表示扰动尚未传到的区域,后阵面以内的部分是波已传过并恢复了原来状态的区域,前后阵面之间的区域就是受到扰动影响的部分,在三维,波的传播有清晰的前阵面与后阵面,称为惠更斯原理或称无后效现象. 6波的弥散 对二维波动方程,点(x0,y0)的影响区域为(xx0)2+(yy0)2a2t2 若在有界区域内有一个初始扰动,则波的传播只有前阵面而无后阵面,所以当的初始扰动传到某点后,扰动对此点的影响不会消失,不过随时间的增加而逐渐减弱.这种现象称为波的弥散,或说波具有后效现象. 2.热传导方程 热传导方程的一般形式为式中f(x,t)为连续有界函数. 热传导方程是描述热的传导过程,分子的扩散过程等物理规律的. 对于n维热传导方程的柯西问题的初值条件为式中为连续有界函数,方程的解的表达式为 3.拉普拉斯方程 研究重力场、静力场、磁场以及一些物理现象(如振动、热传导、扩散)的平衡或稳定过程,通常得到椭圆型方程,最典型的方程为拉普拉斯方程u=0及泊松方程u=式中为已知函数,为拉普拉斯算子, 圆或球的狄利克莱问题解的泊松积分 当区域为圆或球时,分别采用极坐标(r,)或球坐标(r,)较为方便. u=0的极坐标形式为 u=0的球坐标形式为 狄利克莱问题解的泊松积分为 1 区域是圆时,u=0, u|r=a=(),解为泊松积分式中()为已知连续函数,()=(+2). 2 区域是球时,u=0, u|r=a=(),解为泊松积分式中()为已知连续函数, 调和函数的性质 二维拉普拉斯方程的连续解称为调和函数,它具有以下重要性质: 1 设函数u(x,y)在以S为边界的有界区域D内调和,在上有连续一阶偏导数,则式中为外法向导数. 2 算术平均值定理 设函数u(x,y)在圆的内部调和,在闭圆上连续,则u(x,y)在圆心的值等于它在圆周上的值的算术平均值. 3 每一个调和函数u(x,y)对x,y无穷次可微. 4 哈拉克第一定理(一致收敛定理) 设uk(x,y), (k=1,2)在有界区域D内调和,在上连续,如果uk(x,y)在D的边界上一致收敛,则在D内也一致收敛,并且极限函数在D内调和. 5 哈拉克第二定理(单调性定理) 设调和函数列uk(x,y), (k=1,2,)在D的某一内点收敛,且对于任意k,uk+1(x,y)uk(x,y)则uk(x,y) 在D内处处收敛于某调和函数,同时在D的每一有界闭子区域上一致收敛. 6 刘维尔定理 如函数u(x,y)在全平面上调和且不是常数,则它不可能有上界和下界. 7 可去奇点定理 设u(x,y)在A点的一个邻域(除A点外)调和且有界,但在A点没有定义,则可定义函数u(x,y)在A点的值,使u在整个A点的邻域(包括A点)内是调和函数. 李雅普诺夫闭曲面与内、外边值问题 设S为En的有限闭曲面,如果满足下列条件,那末S称为李雅普诺夫闭曲面: (i)曲面到处有切面. (ii)存在常数d0,对曲面上每一点P,可作一个以P为中心、d为半径的球,使曲面在此球内的部分和任意一条与P点法线平行的直线相交不多于一点 (iii)曲面上任意二点P1及P2的法线的夹角(P1,P2)满足式中A,为正常数,01,是点P1与P2之间的距离. (iv)从空间任意一点P0看曲面的任一部分的立体角有界,即|k (k为常数)(从点P0看曲面S的立体角为式中表示矢量 ,NP表示S在点P的外法线矢量,dSP表示点P的面积元素.)设D为En的有界区域,其边界S为李雅普诺夫闭曲面.求在D内满足u=0而在S 上满足给定边界条件的解称为内边值问题;求在D外满足u=0而在S上满足给定边界条件的解称为外边值问题. 狄利克莱问题与诺伊曼问题的解 狄利克莱问题 u=0, 诺伊曼问题 u=0,式中MSS,为S上的已知连续函数,为外法向导数. 1 狄利克莱问题的解可表示为面积分.式中v(P)称为面密度,面积分u(M)称为双层位势,rPM为点M与变点P之间的距离,rPM为矢量 ,NP为S在点P的外法线矢量,v(M)满足第二类弗雷德霍姆积分方程(第十五章1): (i)内边值问题 (ii)外边值问题 2 诺伊曼问题的解可表示为面积分式中(P)称为面密度,面积分u(M)称为单层位势,(P)满足第二类弗雷德霍姆积分方程: (i)内边值问题 (ii)外边值问题 定理: 狄利克莱的内外边值问题及诺伊曼的外边值问题有惟一解,而诺伊曼的内边值问题解存在的充分必要条件是: 泊松方程 在区域D内,泊松方程u=(为已知连续函数)有特解: 三维: 体势位 二维: 对数势位式中rPM为点M与变点P之间的距离. 如果已知泊松方程的一个特解U(M),则=uU满足=0,从而泊松方程的边值问题可化为拉普拉斯方程相应的边值问题.四、 基本解与广义解 共轭微分算子与自共轭微分算子 算子称为二阶线性微分算子,式中aij,bi,c为x1,x2,xn的二次连续可微函数.由公式决定的算子L*称为L的共轭微分算子.如果L=L*,则称L为自共轭微分算子. 格林公式 1 算子L的格林公式是式中S为区域D的边界,N为S的外法线矢量,ei为xi的轴的矢量(0,0,0,0),cos(N,ei)表示矢量N与ei的夹角的余弦, 2 三维拉普拉斯算子的格林公式其中是外法向导数. 3 算子的格林公式式中L*为L的共轭微分算子,N为外法线矢量,i,j分别为x轴,y轴上的单位矢量. 基本解 1 方程Lu=f的基本解: 设M,M0为En中的点,满足方程的解U(M,M0)称为方程Lu=f的基本解,有时也称为方程Lu=0的基本解,式中(M-M0)称为n维狄拉克函数(-函数). 基本解U(M,M0)满足 (i) LU(M,M0)=0,当MM0, (ii) 对任意充分光滑的函数f(M),于是U(M,M0)满足Lu=f(M) . 所以有时也就把满足条件(i)、(ii)的函数U(M,M0)定义为方程Lu=f(M)的基本解. (a)u=0的基本解 二维: 三维: n维: 式中表示点M与M0之间的距离. (b)n维空间的多重调和方程mu=0的基本解 (c)热传导方程的基本解 (d)波动方程的基本解 一维: 二维: 三维: 2 柯西问题的基本解 (i)称满足的解为波动方程柯西问题的基本解,它的形式为 一维: 二维: 三维: (ii)称满足的解为热传导方程柯西问题的基本解,它的形式是 同样方法可以定义其他定解问题的基本解. 由定义可见,基本解表示由集中量(如点热源,点电荷等)所产生的解,下段介绍的格林函数,黎曼函数也具有这种特点,统称它们为点源函数,或影响函数. 广义解 在区域D中给定二阶线性方程式中f在D上连续. 1 设un(x)为D上充分光滑(如二阶连

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论