1、数学物理方程 与特殊函数,第1讲,第一章 一些典型方程和定解条件的推导,1.1 基本方程的建立,例1 弦的振动 设有一根均匀柔软的细弦, 平衡时沿直线拉紧, 而且除受不随时间而变的张力作用外, 不受外力影响. 下面研究弦作微小横向振动的规律. 所谓“横向“是指全部运动出现在一个平面上, 而且弦上的点沿垂直于x轴的方向运动. 所谓“微小“是指的振动的幅度及弦在任意位置处切线的倾角都很小, 以致它们的高于一次方的项都可略而不计.,x,O,u,a,T,a,T,M,M,N,N,x,x+dx,ds,或,或,最后得,(1.3)式称为一维波动方程.,如果在振动过程中, 弦上另外还受到一个与弦的振动方向平行的
2、外力, 且假定在时刻t弦上x点处的外力密度为F(x,t), 则按照同样的推动办法, 可得到弦的强迫振动方程为,其中 表示t时刻单位质量 的弦在x点处所受的外力密度,方程(1.3)与(1.3)的差别在于(1.3)的右端多了一个与未知函数u无关的项f(x,t), 这个项称为自由项. 包括有非零自由项的方程称为非齐次方程, 自由项恒等于零的方程称为齐次方程. (1.3)为齐次一维波动方程, (1.3)为非齐次一维波动方程.,例2 传输线方程 对于直流电或低频的交流电, 电路的基尔霍夫定律指出同一支路中电流相等, 但对于较高频率的电流(指频率还没有高到能显著地辐射电磁波的情况), 电路中导线的自感和电
3、容的效应不可忽略, 因而同一支路中电流未必相等.,今考虑一来一往的高频传输线, 它被当作具有分布参数的导体. 在具有分布参数的导体中, 电流通过的情况, 可以用电流强度i和电压v来描述, 此处i与v都是x,t的函数, 记作i(x,t)与v(x,t).,L每一回路单位的串联电感; C每一单位长度的分路电容.,(1.4)对t求偏导, 将(1.5)对x求偏导, 得,(1.4)对x求偏导, 将(1.5)对t求偏导, 得,最后得,这两个方程称为高频传输线方程. 若令,这两个方程与(1.3)完全相同. 由此可见, 同一个方程可以用来描述不同的物理现象,复习高等数学内容 以i,j,k分别表示沿x,y,z轴正
4、向的单位向量, 则任何一个三维空间的点或者向量a可表示为 a=xi+yj+zk, 或表示为a=(x,y,z) 则二向量a1=(x1,y1,z1)和a2=(x2,y2,z2)的数量积或点积为: a1a2=x1x2+y1y2+z1z2 它们的向量积或叉积为,用M来代表三维空间中的一点x,y,z, 则一个数量场函数就是一个三元函数, 用f(M)表示. 一个向量场函数F(M), 表示每给一点, 都有一个向量与之对应, F(M)可表示为 F(M)=P(M)i+Q(M)j+R(M)k 其中P(M),Q(M),R(M)都是数量函数. 一个数量场函数f(M)在某点x,y,z的沿l方向的方向导数为,其中a,b,
5、g为方向l的方向角,数量场函数f(M)的梯度是一个向量场函数, 记为grad f(M),函数在某点的梯度是这样一个向量, 它的方向与取得最大方向导数的方向一致, 而它的模为方向导数的最大值.,向量场函数F(M)=P(M)i+Q(M)j+R(M)k在点M处的散度是一个标量场函数, 记为div F,在这里可看作稳定流动的不可压缩流体在点M的源头强度在单位体积内所产生的流体的质量. 如果为负, 表示点M处流体在消失.,向量场函数F(M)=P(M)i+Q(M)j+R(M)k在点M处的旋度是一个向量场函数, 记为rot F,向量微分算子的定义为,它又称为(Nabla)算子或哈密尔顿(Hmilton)算子
6、. 运用向量微分算子, 设数量场函数为f(x,y,z), 向量场函数为F(x,y,z)=Pi+Qj+Rk,其中,称为拉普拉斯(Laplace)算子.,高斯公式和斯托克斯公式可分别写成,其中Fn为向量场函数在曲面元法线方向上的投影, 而Ft为向量场函数在曲线切线方向上的投影.,例3 电磁场方程 电磁场的特性可用电场强度E与磁场强度H以及电感应强度D与磁感应强度B来描述. 它们可用麦克斯韦(Maxwell)方程组描述为:,J传导电流面密度,r电荷的体密度.,e介电常数,m导磁率,s导电率.,而rotrot H=graddiv H-2H=2H,最后得到H所满足的方程为,同理, 若消去H即得E所满足的
7、方程,如s=0, 则,(1.15)与(1.16)称为三维波动方程. 用标量形式表示可写成,u是E(或H)的任意一个分量, 而,从(1.11)和(1.12)还可以推出静电场的电位所满足的微分方程. 因E=-grad u, 由 div D=e div E=r 可得 div grad u=-r/e, 或,此非齐次方程称为泊松(Poisson)方程. 如果静电场是无源的, 即r=0, 则(1.18)变成2u=0 (1.19) 这个方程称为拉普拉斯(Laplace)方程.,例4 热传导方程 在物体中任取一闭曲面S, 它所包围的区域记作V. 假设在时刻t区域V内点M(x,y,z)处的温度为u(x,y,z,
8、t), n为曲面元素DS的法向(从V内指向V外). 由传热学中傅里叶实验定律可知, 物体在无穷小时间段dt内, 流过一个无穷小面积dS的热量dQ与时间dt, 曲面面积dS, 以及物体温度u沿曲面dS的法线方向的方向导数三者成正比,热场,k为热传导系数.这里为常数。从时刻t1到t2, 通过曲面S流入区域V的全部热量为(负号表明热量是由高温向低温流动),上式的负号表示热流流向是温度梯度的相反方向。,因S为闭曲面, 式中的二重积分可利用高斯公式化为三重积分, 即,流入的热量使V内温度发生了变化, 在时间间隔t1,t2内区域V内各点温度从u(x,y,z,t1)变化到u(x,y,z,t2), 则在t1,
9、t2内温度升高所需要的热量为,c物体的比热, r物体的密度.,流入的热量应等于物体吸收的热量, 因此有,由于t1,t2及V是任取的, 上式成立的条件就是被积函数相等, 即,其中,方程(1.21)称为三维热传导方程.,若物体内有热源, 其强度为F(x,y,z,t), 则相应的热传导方程为,其中,作为特例, 如果所考虑的物体是一根细杆(或一块薄板), 或者即使不是细杆(或薄板)而其中的温度u只与x,t(或x,y,t)有关, 则(1.21)就变成一维或二维热传导方程,如果考虑稳恒温度场, 即在热传导方程中物体的温度趋于某种平衡状态, 这时温度u已与时间无关, 此时方程(1.21)就变成拉普拉斯方程(
10、1.19). 由此可见稳恒温度场内的温度u也满足 拉普拉斯方程. 在研究气体或液体的扩散过程时, 若扩散系数是常数, 则所得的扩散方程与热传导方程完全相同.,1.2 初始条件与边界条件,上面讨论的是如何将一个具体问题所具有的物理规律用数学式子表达出来. 除此以外, 我们还需要把这个问题所具有的特定条件也用数学形式表达出来, 这是因为任何一个具体的物理现象都是处在特定条件之下的. 例如弦振动问题, 上节所推导出来的方程是一切柔软均匀的弦作微小横向振动的共同规律, 在推导这个方程时没有考虑到弦在初始时刻的状态以及弦所受的约束情况.如果我们不是泛泛地研究弦的振动, 势必就要考虑到弦所具有的特定条件.
11、,因为任何一个振动物体在某时刻的振动状态总是和此时刻以前的状态有关, 从而就与初始时刻的状态有关. 另外, 弦的两端所受的约束也会影响弦的振动, 端点处的物理条件不同就会产生不同的影响, 因而弦的振动也不同. 所以对弦振动问题来说, 除了建立振动方程以外, 还需列出它所处的特定条件. 对热传导方程, 拉普拉斯方程也是如此.,提出的条件应该能够用来说明某一具体物理现象的初始状态或者边界上的约束情况. 用以说明初始状态的条件称为初始条件; 用以说明边界上的约束情况的条件称为边界条件.,下面具体说明初始条件和边界条件的表达形式. 先谈初始条件, 对于弦振动问题来说, 初始条件就是弦在开始时刻的位移及
12、速度, 若以j(x), y(x)分别表示初位移和初速度, 则初始条件可以表达为,而对热传导方程来说, 初始条件是指在开始时刻物体温度的分布情况, 若以j(M)表示t=0时物体内任一点M处的温度, 则热传导方程的初始条件就是 u(M,t)|t=0=j(M). (1.23) 泊松方程与拉普拉斯方程都是描述稳恒状态的, 与初始状态无关, 所以不提初始条件.,再谈边界条件. 还是先从弦振动问题说起, 从物理学得知, 弦在振动时, 其端点(以x=a表示这个端点)所受约束情况, 通常有以下三种类型: 第一, 固定端, 即弦在振动过程中这个端点始终保持不动. 对应于这种状态的边界条件为 u|x=a=0, (
13、1.24) 或 u(a,t)=0.,第二, 自由端, 即弦在这个端点不受位移方向的外力, 从而在这个端点弦在位移方向的张力应该为零. 由前面的推导过程可知, 此时对应的边界条件为,即,或,对热传导方程来说, 也有类似的情况. 以S表示某物体V的边界, 如果在导热过程中边界S的温度为已知的函数f(x,y,z,t), 则这时的边界条件为 u|S=f (1.27) 这里的f是定义在S上(一般依赖于t)的函数. 如果在导热过程中, 物体V与周围介质处于绝热状态, 则在S上的热量流速始终为零, 则在边界S上必满足,所以, 当物体和外界有热交换时, 相应的边界条件为,即,其中s=k1/k,综上所述可知,
14、不论对弦振动问题, 还是对热传导问题, 它们所对应的边界条件, 从数学角度看不外有三种类型: 一是在边界S上直接给出了未知函数u的数值,即 u|S=f1 (1.30) 这种形式的边界条件称为第一类边界条件. 二是在边界S上给出了未知函数u沿S的外法线方向的方向导数, 即,这种形式的边界条件称为第二类边界条件.,三是在边界S上给出了未知函数u及其沿S的外法向的方向导数某种线性组合的值, 即,这种形式的边界条件称为第三类边界条件. 需要注意的是(1.30),(1.31),(1.32)右端的fi(i=1,2,3)都是定义在边界S上(一般说来, 也依赖于t)的已知函数. 不论哪一类型的边界条件, 当它
15、的数学表达式中的自由项(即不依赖于u的项)恒为零时, 则这种边界条件称为是齐次的, 否则称为非齐次的.,1.3 定解问题的解法,前面推导了三种不同类型的偏微分方程并讨论了与它们相应的初始条件与边界条件的表达方式. 由于这些方程中出现的未知函数的偏导数的最高阶都是二阶, 而且它们对于未知函数及其各阶偏导数来说都是线性的, 所以这种方程称为二阶线性偏微分方程. 在工程技术上二阶偏微分方程遇到最多.,如果一个函数具有某偏微分方程中所需要的各阶连续偏导数, 并且代入该方程中能使它变成恒等式, 则此函数称为该方程的解(古典解). 由于每一个物理过程都处在特定的条件之下, 所以我们的任务是要求出偏微分方程
16、的适合某些特定条件的解. 初始条件和边界条件都称为定解条件, 把某个偏微分方程和相应的定解条件结合在一起, 就构成了一个定解问题.,只有初始条件, 没有边界条件的定解问题称为始值问题(或柯西(Cauchy)问题); 反之, 没有初始条件, 只有边界条件的定解问题称为边值问题. 即有初始条件也有边界条件的定解问题称为混合问题.,一个定解问题提得是否符合实际情况, 当然必须靠实践来证实. 然而从数学角度来看, 可以从三方面加以检验: 1) 解的存在性, 即看所归结出来的定解问题是否有解; 2) 解的惟一性, 即看是否只有一个解; 3) 解的稳定性, 即看当定解条件有微小变动时, 解是否相应地只有微
17、小的变动. 如果确实如此, 此解便称为稳定的, 否则所得的解就无实用价值.,如果一个定解问题存在惟一且稳定的解, 则此问题称为适定的. 在以后讨论中我们把着眼点放在讨论定解问题的解法上, 而很少讨论它的适定性, 这是因为讨论定解问题的适定性往往十分困难, 而本书所讨论的定解问题都是经典的, 它们的适定性都是经过证明了的.,一个含n个自变量的二阶线性偏微分方程的最一般的形式应该是,其中Aik, Bi, C, f都只是x1,x2,xn的已知函数, 与未知函数无关. 在两个自变量的情形, (1.33)式可写为,线性偏微分方程(1.33)满足叠加原理, 即若ui是方程 Lui=fi (i=1,2,) 的解, 而且级数,收敛, 并且能够逐项微分两次, 其中Ci(i=1,2,)为任意常数, 则u一定是方程,的解(假定这个方程右端的级数是收敛的).,特别是, 如果ui(i=1,2,)是二阶线性齐次方程 Lu=0 的解, 则只要,收敛, 并且可以逐项微分两次, u一定也是此方程的解. 这个结论的证明是非常容易的(留作练习, 但它却是下一章要讲的分离变量法的出发点.,作业 习题一 第17页 第1,3,5题,
