《三维含时波动方程在球坐标下的求解.doc》由会员分享，可在线阅读，更多相关《三维含时波动方程在球坐标下的求解.doc（19页珍藏版）》请在金锄头文库上搜索。
1、毕 业 论 文 题 目： 三维含时波动方程在球坐标下的求解 学 号： 06110901020 姓 名： 朱启龙 教 学 院： 物理科学与技术学院 专业班级： 09物理本科（1）班 指导教师： 高慧昀 完成时间： 2013年 5 月 8 日 毕节学院教务处制三维含时波动方程在球坐标下的求解作者姓名：朱启龙 专业班级:09物本（1）班 学号：06110901020 指导老师： 高 慧 昀摘要：本文利用行波法和分离变量法两种方法研究了三维含时波动方程在球坐标系下的初值问题并进行求解。关键词：波动方程 ; 分离变量法 ; 行波法 ; 球坐标系The solution of the wave equat
2、ion in spherical coordinates in three-dimensional timeCandidate:ZhuQiLong Major: phsicesStudent No:06110901020 Advisor: GaoHuiYunAbstract：This paper line wave method of separation of variables method are two ways to study the three-dimensional wave equation and solve the initial value problem in sph
3、erical coordinates.Key words: wave equation; Method of separation of variables ; Traveling wave method ; Spherical coordinates目录引言11 行波法求解三维含时波动方程21.1球面平均法降维21.2 行波法求解32 分离变量法求解三维含时波动方程62.1 分离t,r, ,并求T62.2球函数方程的求解(,)82.3贝赛尔方程2.1.7的求解(r)113 小结13参考文献14致谢15IV引言三维含时波动方程描述了波（一种重要的描述自然界中的不同种类的波动现象的偏微分方程）在
4、均匀各项同性弹性介质中的传播。三维含时波动方程是随时间及空间的变化而变化的。对其初值问题的求解，如按我们熟悉常微分方程求解的方法是很难求解的。因此，我们熟悉的常微分方程的求解只适用于少数的某些定解问题，运用行波法（达朗贝尔法）就能解决这个难题；但是，行波法对三维情况也很难求解，故而引人球面平均法将三维情况降至一维以便于运用行波法求解。除此之外，分离变量法也适用于大量的各种各样定解问题，其基本思想是把偏微分方程分解为几个常微分方程，其中有的常微分方程带有附加条件而构成本征值问题1-3。本文利用行波法和分离变量法两种方法研究了三维含时波动方程在球坐标系下的初值问题并进行求解。1行波法求解三维含时波
5、动方程由于行波法对三维含时波动方程很难直接求解，故本章希望采用某种办法将三维问题转化为一维问题.为此引入球面平均法降维4.如三维含时波动方程： （0.1） （0.2） (0.3)在球坐标系下的算符的表达式为 (0.4)1.1球面平均法降维用表示以为中心，r为半径的球面，表示函数在球面上的平均值，则单位球面上的立体角元，则是只含有r,t的二元函数，下面均用代替，且时，既，将（0.1）式、(0.2)式、(0.3)式在球面上取平均值得: (1.1.1) (1.1.2) (1.1.3)由于为r和t的函数，所以在球坐标系下，由(0.4)式可得 (1.1.4 ) 将(1.14)式代入(1.1.1)式有既
6、(1.1.5)令，将其代入(1.1.2)、（1.1.3）、（1.1.5）式后有 (1.1.6)式（1.1.6）就是用球面平均法降维后所得的形式.1.2 行波法求解到目前为止，（1.1.6）式已完全属于一维r的含时齐次波动方程，此时利用行波法求解是非常方便的，其特征线方程为 解之得 令 (1.2.1)则代入(1.1.6)式中的有 （1.2.2）上式对求积分得继续对求积分得将z代入（1.1.6）式后得 （1.2.3）又将（1.2.3）式代入（1.1.6）式中的，有 （1.2.4） （1.2.5）将(1.2.5)变形后得代入（1.2.4）得当r=0时，由（1.1.6）式有则当，既时当，既时代入（1.
7、2.3）式得因此（1.1.6）式的解为又因为，则代入初始条件得 则可以写成 (1.2.6)也可以将其写为（1.2.7）将（1.2.6）和（1.2.7）式叫做三维含时波动方程的泊松公式，从这两式可以看出解在处，于时刻t的值由以为球心，为半径的球面上的初始扰动所确定.故称球面为点解值的决定区域，这是因为在球面上的初始扰动，经过时间t刚好都传播到了点处，换一种说法，在一点出的初始扰动，经过时间t刚好传播在球面上，这表示波以速度a传播4.2分离变量法求解三维含时波动方程分离变量法是将偏微分方程分解为几个常微分方程（有些关于空间的常微分方程带有附件条件构成本征值问题）之积的形式。分离变量法适用于大量的各
8、种各样的定解问题，本章主要针对n维空间的齐次问题（齐次方程和齐次边界问题）5-7。例如：三维空间的齐次问题2.1 分离t,r, ,并求T对波动方程（0.1）式分离时间变量t和空间变量r,令 （2.1.1）将其代入（1.1.1）式得变形得 两边分别是空间的函数和时间t的函数，则是不会相等的，若要相等，只有两边都等于相同的常数C.令C=-,且为实数，由后面的解答可以知道是本征值或者是两个本征值相加的值.于是有 于是得到两个微分方程 （2.1.2） （2.1.3）而方程的解为 （2.1.4）且方程（2.1.2）式称为亥姆霍兹方程.将（0.4）式代入（2.1.2）式有 （2.1.5）继续分离变量，令，
9、代入(2.1.5)式，并在等式两边同时乘以并移项得 左边是r的函数，右边是与的函数，则左右要相等，只有左边和右边同等于同一个常数D，令，则有于是得两个方程 （2.1.6） （2.1.7）由（2.1.4）式可知，直接分离时间变量和空间变量便可得到的函数表达式.2.2球函数方程的求解(,)方程（2.1.6）变形得球函数方程，既 （2.2.1）分离变量、，令,将其代入球函数方程（2.2.1）式有两边同时乘以并移项有两边分别是和的函数，则也不可能左右相等，只有左右都等于一个相同的常数E则又得两个微分方程 （2.2.2） （2.2.3）常微分方程（2.2.2）满足自然的周期条件，一个确定的极角可以加减2
10、的整数倍，既，则方程（2.2.2）和上述条件构成本征值问题6-8，本征值是 （2.2.4）将（2.2.4）式代入（2.2.2）式有解之得 （2.2.5）将（2.2.4）式代入（2.2.3）式有 （2.2.6）令，相当于，则有于是（2.2.6）可变为既 （2.2.7）当m=0时，上式变为 （2.2.8）（2.2.7）和（2.2.8）分别叫作l阶连带勒让得方程和l阶勒让得方程对l阶勒让得方程（2.2.7）运用级数解法在=0的领域上求解l阶勒让得方程（2.2.8），先将其变形得 （2.2.9）可见它的系数，都是有限值，二者在=0一定是解析的。所以，=0点一定是方程的常点。由常点领域上解的定理，且其解
11、的泰勒形式为则将上面的这些式子全部代入2.16,合并同幂项后各项的系数如下表.常数项X项项.项.令同幂项的系数和为零，找出系数的递推公式既 (2.2.10)如此，l阶勒让得方程的解可以写为 （2.2.11）将代入后得 （2.2.12）对l阶连带勒让得方程（2.2.7）的求解运用级数解法由于比较复杂，故作变换则有 其中、分别表示f对x求二阶导、一阶导，将以上三式代入（2.2.7）可把（2.2.7）变为 （2.2.13）（2.2.13）运用级数法来解会简单得多，这里不用此方法，只需要将勒让得方程（2.2.7）对x求导m次，可得其中，分别表示对x求次、次导，于是 （2.2.14）可见该式既为（2.2.13）式，所以（2.2.13）的解就是勒让得方程的解（这里用表示，以便进行区别）的m阶导数，既 （2.2.15）2.3
