第一章-数学物理中的偏微分方程.ppt

1、数学物理方程EquationsofMathematicalPhysics,2,数学物理方程指从物理学或其他各门自然科学、技术科学中的某些物理问题导出的偏微分方程(有时也包括积分方程、微分积分方程等)。它们反映了有关的未知变量关于时间的导数和与空间变量的导数之间的制约关系。连续介质力学、电磁学、量子力学等方面的基本方程都属于数学物理方程的范围。教学目的通过本课程的教学使学生获得有关偏微分方程的一些基本概念、基本方法，掌握三类典型方程定解问题的解法，进一步扩大学生的数学知识面，为后继课程提供必要的数学基础。,3,参考书目,数学物理方程,王明新,清华大学出版社。,数学物理方程，姜礼尚，高教出版社。,

2、工程技术中的偏微分方程，潘祖梁，浙江大学出版社。,1.1偏微分方程的一些基本概念,4,5,一.偏微分方程（partialdifferentialequation)（PDE）的基本概念,自变量,未知函数,偏微分方程的一般形式,6,PDE的阶：,PDE的解,古典解,广义解,概念,是指这样一个函数，它满足方程，并且在所考虑的区域内有m阶连续偏导数。,线性PDE,非线性PDE,半线性PDE,拟线性PDE,完全非线性PDE,自由项在偏微分方程中，不含有未知函数及其偏导数的项称为自由项,7,线性PDE：,PDE中对所含未知函数及其各阶导数的全体都是线性的。例如：,常系数线性PDE:,不然称为变系数的,齐次

3、线性PDE:,不然称为非齐次的,线性PDE的主部:具有最高阶数偏导数组成的部分,主部,8,PDE中对最高阶导数是线性的。例如：,半线性PDE：,完全非线性PDE：,PDE中对最高阶导数不是线性的。,拟线性PDE：,拟线性PDE中，最高阶导数的系数仅为自变量的函数。例如：,非线性PDE,9,举例（未知函数为二元函数）,1.,2.,变换,解为：,解为：,10,举例（未知函数为二元函数）,4.,3.,解为：,变换,解为：,11,5.,不易找出其通解，但还是可以找出一些特解,任意解析函数的实部和虚部均满足方程。,也是解,6.,特解都不易找到,KDV方程,举例（未知函数为二元函数）,12,7.,拟线性P

4、DE,8.,拟线性PDE,9.,半线性PDE,10.,半线性PDE,11.,完全非线性PDE,1.2三个典型的方程,13,14,举例(多元函数),拉普拉斯(Laplace)方程,热传导方程,波动方程,物理模型与定解问题的导出,15,16,弦振动方程的导出,17,一长为L的柔软均匀细弦，拉紧后，当它受到与平衡位置垂直的外力作用时，开始作微小横振动。假设这运动发生在同一平面内，求弦上各点位移随时间变化规律。,弦上各点作往返运动的主要原因在于弦的张力作用，弦在运动过程中各点的位移、加速度和张力都在不断变化，但它们遵循物理的运动规律。由此可以建立弦上各点的位移函数所满足的微分方程。,物理背景：波的传播

5、和弹性体振动。,弦振动方程的导出首先，考察弦横振动这个物理问题：给定一根两端固定的拉紧的均匀柔软的弦线，设其长度为l，它在外力作用下在平衡位置附近作微小的横振动，求弦上各点的运动规律。把实际问题提炼为数学模型时必须做一定的理想化假设，以便抓住问题的最本质特征。,基本假设：1.弦的质量是均匀的，弦的截面直径与长度相比可以忽略。弦可以视为一条曲线，线密度为常数。（细弦）2.弦在某一个平面内作微小横振动。弦的位置始终在一直线段附近，弦上各点在同一平面内垂直于该直线的方向上作微小振动。（微幅）3.弦是柔软的，它在形变时不抵抗弯曲。弦上各质点的张力方向与弦的切线方向一致，而弦的伸长变形与张力的关系服从虎

6、克定律。（横振动）基本规律：牛顿第二定律（冲量定律）,弦线上任意一点在t时刻沿y轴上的位移,研究对象：,在右图所示的坐标系，用u(x,t)表示弦上各点在时刻t沿垂直于x方向的位移。在这条弦上任意取一弦段(x,x+x)，它的弧长为：,由假设3，弦线张力T(x)总是沿着弦在x处的切线方向由于弦只在垂直x轴的方向进行横振动，因此可以把弦线的张力T(x)在x轴的方向的分量看成常数。对于图中选取的弦段而言，张力在x轴的垂直方向上的合力为：,在时间段(t,t+t)内该合力产生的冲量为：,另一方面，在时间段(t,t+t)内弦段(x,x+x)的动量变化为：,于是由冲量定理：,从而有：,进一步由t，x的任意性，

7、有,假定有垂直于x轴方向的外力存在，并设其线密度为F(x,t)，则弦段(x,x+x)上的外力为：,它在时间段(t,t+t)内的冲量为：,类似地，三维波动方程可以表示为：,于是有：,简化假设：,(2)振幅极小，张力与水平方向的夹角很小。,(1)弦是柔软的，弦上的任意一点的张力沿弦的切线方向。,牛顿运动定律：,横向：,纵向：,其中：,其中：,其中：,一维波动方程,令：,-非齐次方程,自由项,-齐次方程,忽略重力作用：,非均匀弦的强迫横振动方程,一维波动方程不仅可以描述弦的振动，还可以描述：弹性杆的纵向振动管道中气体小扰动的传播等等因此，一个方程反应的不止是一个物理现象，而是一类问题。,28,2+1

8、维波动方程或膜振动方程,一块均匀的拉紧的薄膜，离开静止水平位置作垂直于水平位置的微小振动，其运动规律满足,29,3+1维波动方程或声波方程,n+1维波动方程,1.4定解条件和定解问题,30,列出微分方程的目的是要从微分方程中求得具体问题的解或者研究解的性质。前面我们看到，弦振动方程描述的是弦作微小横振动时的位移函数u(x,t)所应满足的一般性规律。仅仅利用它并不能完全确定一条弦的具体运动状况。这是因为弦的运动还与其初始状态以及边界所处的状况有关系，因此对于具体的弦振动问题而言，还需要结合实际问题附加某些特定条件。例如:在前面的推导中，弦的两端被固定在x=0和x=l两点，即u(0,t)=0，u(

9、l,t)=0，这两个等式称为边界条件。此外，设弦在初始时刻t=0时的位置和速度为,这两个等式称为初始条件。边界条件和初始条件总称为定解条件。把微分方程和定解条件结合起来，就得到了与实际问题相对应的定解问题。对于弦振动方程而言，与上述定解条件结合后，其定解问题可以描述为：,定解条件,要在区域,上（见右上图）求上述定解问题的解，就是,要求这样的连续函数u(x,t)，它在区域00中满足波动方程(2.1)；在x轴上的区间0,l上满足初始条件(2.2)；并在边界x=0和x=l上满足边界条件(2.3)和(2.4)。,一般称形如(2.3)和(2.4)的边界条件为第一类边界条件，也叫狄利克雷（Dirichle

10、t）边界条件。,定解条件,波动方程的初始条件,1、初始条件描述系统的初始状态,系统各点的初位移系统各点的初速度,定解条件,(2)自由端：x=a端既不固定，又不受位移方向力的作用。,2、边界条件描述系统在边界上的状况,波动方程的三类边界条件,(1)固定端：对于两端固定的弦的横振动，其为：,或：,(3)弹性支承端：在x=a端受到弹性系数为k的弹簧的支承。,或,诺依曼（Neumann）边界条件,狄利克雷（Dirichlet）边界条件,同一类物理现象中，各个具体问题又各有其特殊性。边界条件和初始条件反映了具体问题的特殊环境和历史，即个性。,初始条件：够用来说明某一具体物理现象初始状态的条件。,边界条件

11、：能够用来说明某一具体物理现象边界上的约束情况的条件。,其他条件：能够用来说明某一具体物理现象情况的条件。,定解条件,定解问题,定解问题适定性概念,(1)初始问题：只有初始条件，没有边界条件的定解问题；(2)边值问题：没有初始条件，只有边界条件的定解问题；(3)混合问题：既有初始条件，也有边界条件的定解问题。,把某种物理现象满足的偏微分方程和其相应的定解条件结合在一起，就构成了一个定解问题。,定解问题的检验,解的存在性：定解问题是否有解；解的唯一性：是否只有一解；解的稳定性：定解条件有微小变动时，解是否有相应的微小变动。,37,经典的定解问题举例,维波动方程(弦振动方程)的初值问题,38,经典

12、的定解问题举例,热传导方程的初值问题,39,经典的定解问题举例,二维调和方程的边值问题,第一边值问题（Dirichlet）,第二边值问题（Neumann）,第三边值问题（Robin）,40,经典的定解问题举例,热传导方程的初、边值问题,41,何为适定性？,存在性唯一性连续依赖性（稳定性）,适定性,若PDE在附加条件及求解域的一定要求下，它的解在已知度量的某函数类中存在、唯一而且关于附加条件为稳定的，就称定解问题在相应的函数类中为适定的。,稳定性：只要定解条件的偏差足够小，相应的定解问题解的偏差也将非常小,除了研究定解问题的适定性外，数理方程中还经常研究的问题包括：解的正则性（光滑性）、解的渐近

13、性（包括衰减性）和定解问题的求解方法（精确解、渐近解、数值解）等。,定解问题适定性概念,43,定解问题的适定性,定解问题,PDE,定解条件,初值条件initialcondition,边值条件boundarycondition,初、边值条件,初值问题、边值问题、混合问题,热传导方程,44,45,热传导分析：设杆长方向为x轴，考虑杆上从x到x+dx的一段(代表)，其质量为dm=dx，热容量为cdm。设杆中的热流沿x轴正向，强度为q(x,t)，温度分布为u(x,t)，则问题：一根长为L的均匀导热细杆，侧面绝热，内部无热源。其热传导系数为k，比热为c，线密度为。求杆内温度变化的规律。,由能量守恒定律c

14、dmdu=dQ=q(x,t)-q(x+dx,t)dt=-qx(x,t)dxdt于是有cut=-qx由热传导定律q(x,t)=-kux(x,t)代入前面的式子，得到cut=kuxxut=a2uxx,46,推广：情况：内部有热源(或侧面不绝热)分析：设热源强度(单位时间在单位长度中产生的热量)为F(x,t)，代表段的吸热为Fdxdt方程：cut=kuxx+Fut=a2uxx+f，f=F/(c),稳定场方程,47,48,产生：在演化问题中，有时会到达一个不随时间变化的稳定状态，对应的方程称为稳定场方程。形式：在对应的演化方程中取消时间变量t，对t的导数为零。分类：无外界作用情况拉普拉斯方程：u=ut

15、t+uyy+uzz=0有外界作用情况泊松方程：u=utt+uyy+uzz=f(x,y,z)典型应用静电场方程：u=-/稳定温度分布：u=-F/k,数学物理方程的分类,在数学物理方程的建立过程中，我们主要讨论了三种类型的偏微分方程：波动方程；热传导方程；稳定场方程这三类方程描写了不同物理现象及其过程,49,50,二阶线性PDE方程的分类,两个自变量，齐次,主部,目的：,通过自变量的非奇异变换来简化方程的主部，从而据此分类。,非奇异,（1）,51,复合求导,52,系数之间的关系,（2）,（1）,（3）,53,其他系数之间的关系,（3*）,54,考虑,如若能找到两个相互独立的解,那么就作变换,从而有

16、,（4）,55,假设,是方程,的特解，则关系式,是常微分方程,（4）,（5）,的一般积分。反之亦然。,引理,由此可知，要求方程（4）的解，只须求出常微分方程（5）的一般积分。,56,定义,称常微分方程（5）为PDE（1）的特征方程。,称（5）的积分曲线为PDE（1）的特征曲线。,（6）,57,记,定义,方程(1)在点M处是,双曲型：,椭圆型：,抛物型：,若在点M处，有,若在点M处，有,若在点M处，有,58,双曲型PDE,右端为两相异的实函数,它们的一般积分为,由此令,，方程（1）可改写为,双曲型方程的第一标准型,双曲型方程的第二标准型,59,抛物型PDE,由此得到一般积分为,由此令,，其中,与

17、,独立(线性无关)的任意函数。,60,由于,由此推出,61,因此，方程（1）可改写为,抛物型方程的标准型,而,62,椭圆型PDE,右端为两相异的复数,由此推出两族复数积分曲线为,其中,63,由此令,从而方程（1）可改写为,，满足方程（4）,椭圆型方程的标准型,64,例1,抛物型方程,令,65,例2,双曲型方程,66,例3,Tricomi方程,椭圆型,双曲型,抛物型,67,叠加原理,弦振动方程的达朗贝尔解法,叠加原理,从本节开始我们讨论弦振动方程的各类定解问题。在此之前，先介绍叠加原理,在物理学研究中经常出现这样的现象：几种不同原因的综合所产生的效果等于这些不同原因单独（假设其他原因不存在）产生

18、的效果的累加。这就是叠加原理。,典型例子：力和加速度的关系，万有引力场的可叠加性复杂的声音各种单音的叠加电磁场中的叠加原理,则对于任意的常数C1、C2，函数,因此，弦振动方程满足叠加原理,线性方程都满足叠加原理,线性方程解（线性系统）具有叠加特性,几种不同的原因的综合所产生的效果等于这些不同原因单独产生的效果的累加。(物理上),弦振动方程的达朗贝尔解法,先从最简单的情形入手，即首先考察边界的影响可以忽略不计的情况（如果所考察的物体（弦线）长度很长，而我们所关注的又只是在较短时间内且距离边界较远的一段范围中的运动情况，那么边界条件的影响就可以忽略，并不妨把所考察的物体的长度视为无限长）。这样的情

19、况下，定解问题归结为如下形式：,这个定解问题中，定解条件只有初始条件，故通常称为初值问题(也称柯西（Cauchy）问题）。相应地，前一节中的定解问题(1.1)(1.4)由于既有初始条件，又有边界条件，故称为初边值问题或混合问题。,方程(1.5)中的自由项f(x,t)是由于外力作用产生的，因此方程(1.5)中f(x,t)恒为零的情况对应于自由振动；f(x,t)不为零的情况对应于强迫振动。,下面，我们求解上述初值问题。首先注意到微分方程及定解条件都是线性的。对于这种定解问题，同样存在叠加原理，即若u1(x,t)和u2(x,t)分别是下述初值问题,和,的解，那么u=u1(x,t)+u2(x,t)就一

20、定是原初值问题(1.5)、(1.6)的解(证明作为课后习题)。这样求解初值问题(1.5)、(1.6)就转化为分别求解齐次方程带非齐次边界条件的初值问题(I)和非齐次方程带齐次初始条件的初值问题(II),单独初始振动状态对振动过程的影响。,单独考虑外力因素对振动过程的影响。,首先，我们考察代表自由振动情况的初值问题(I)，它可以通过自变量变换的方法求解。引如新自变量：=x-at，=x+at。利用复合函数求导的法则，有,类似地，,从而，方程(1.7)就化为，这个方程可以直接求解。把它关于积分一次，再关于积分一次，就可以得到它的通解为u(,)=F()+G()，其中，F和G是任意两个可微分的单变量函数。代回原来的自变量，方程(1.7)的通解表示为u(x,t)=F(x-at)+G(x+at)。,利用这个通解表达式，就可以利用初始条件(1.8)来决定函数F和G，进而求出初值问题(I)的解。把上述通解表达式代入初始条件(1.8)，得到：,(1.12)式是一个简单的常微分方程，求解它得到,由(1.11)和(1.13)式联立求解可以得出函数F和G,把它们代入方程(1.7)的通解表达式就得到了初值问题(I)的解,这个公式(1.14)称为达朗贝尔公式。从以上推导过程可以看出：如果初值问题(I)有解，则解一定可以根据初始条件由达朗贝尔公式表达出来，因此