球面平均法和泊松公式精编版.docxVIP

§7.2球面平均法和泊松公式

一本节主要思想

(1)对三维波动方程的初值问题，先假设已知空间某一点M3,y,z)的振动U(X,y,z,t)，然后以M点为发射子波的波源，根据球面波的对称性，可根据加权平均的思想来考察球面Sm的平均振动U(r,t)，从而将问题归结为两个自变量的一维波动方程，最后采用极限的思r

想，令r=0，即可得到M点的振动u(x,y,z,t).

(2)对二维波动方程的初值问题，采用将其上升到三维空间的思想，根据已有的三维波动方程的泊松公式，获得此问题的三维解，再采用降维法，最终获得二维解.

二三维波动方程的初值问题，泊松公式

1三维波动方程初值问题解的泊松公式

考察三维波动方程的初值问题：

u=a2Au t0,一8x,y,z^(1)

u\*(X,y,z),u\=w(X,y,z) 一8v,y,z布2)

kt=0 11=0

以SM表示以点M(x,y,z)为心，半径为r的球面，以血表示S的面元，则S的面元r 1 r

dS=r2dw.

[d①=sin0dOd甲

注：｛ ndS=r2d①

[dS=r2sin0d0d甲

下面用加权平均的思想(即球面平均法)求函数(x,y,z,t)在球面SM上的平均值u(r,t):u(r,t)=-^jju(，门,匚,t)dS(3)

4冗r2sm

r

(,n,SeSM,=X+ra,^=y+ra2,匚=z+ra3,(4)

a=sin0cos甲,a=sin0sin甲,a=cos0.(5)

注：(4)、(5)实际上是球的参数方程的表达式

(3)式也可写成

u(r,t)=-!-jju(x+ra,y+ra,z+ra,t)dw，

4兀sm 1 2 3

注：由于dS=r2sin0d0d^=r2d①不含dr,故可将上放入积分号中，与r2正好约掉.r2

由此可知u =u(x,y,z,t)，所以，为了求u可以先求u

r=0

下面就来讨论如何求u:

注意到r=i；(x一)2+(y-n)2+(z一S2，容易算得

8_du8,_dux一dxdr8xdrr

82ur2-(x-g)2du(x-g)2d2u = + .dx2 r3 dr r2dr2

d2ud2u

同理，可求出旬2和兼，将他们相加，得到

.-2dud2u1d2 -

Au= + = (ru).

rdrdr2 rdr2

注：

ddu-—(ru)=r——+udrdr

TOC\o1-5\h\zd2 -d,du_、_du d2u

dr2 (ru)=—(r一+u)=2一+rdr2 dr dr

dr2

1d2/-、2dud2u所以，—(ru)=—k+=

rdr2rdr dr2

将方程(1)两端取球面平均，即得

utta2d

u

tt

=a2Au= (ru)

rdr2

等式两边同乘以r得

TOC\o1-5\h\z\oCurrentDocumentd dz

—(ru)-a2^—(ru)=0 (6)

dt2 dr2

注:与时间t无关所以可将r放入算子瓦中这是关于岳的一维波动方程，其通解为

ru=f(r-at)+f(r+at) (7)

1 2

令r=0，得

0=f1(-at)+f2(at)

从而知f(a)=-f(-a)，故有

12

ru=f(r+at)-f(at-r) (8)

12

注：①a的取值范围应是(-8,0)

②我们所关心的是r-aK0的情况，因为此时表明振动已传至所考察的球面，若

—at0,表明振动还未传至所考察的球面，这样就无法考察心点的振动.

为求，将上式对r求导，得

r=0

8/ _du

—(ru)=u+r——

dr dr

=f(r+at)+f(at-r)

(9)

令r=0,即得

u(尤，y,z,t)=u

=2f(at).

r=0 2

(10)

所以，为求u，只需求?为此将⑻式再对t求导

d

—(ru)=a|f(r+at)-f(at-r)J.

d,-、 厂…

dr l2、

由(9)、(11)两式得

d 1dz_X

—(ru)+(ru=dr adt

在上式中令t=0，并注意到

-“rd-

2fr=〒他书一*侦)_dr adt _

t=