波动方程或称波方程.docx

波动方程或称波方程.docx

《波动方程或称波方程.docx》由会员分享，可在线阅读，更多相关《波动方程或称波方程.docx（9页珍藏版）》请在冰豆网上搜索。

波动方程或称波方程

波动方程或称波方程（英语:

waveequation）是一种重要的偏微分方程,主要描述自然界中的各种的波动现象,包括横波和纵波，例如吏波、光波、无线电波和生波。

波动方程抽象自声学、物理光学、电磁学、电动力学、流体力学等领域。

历史上许多科学家，如达朗贝尔、欧拉、丹尼尔•伯努利和拉格朗日等在研究兀蟲等物体中的弦振动问题时，都对波动方程理论作出过重要贡献。

波动方程是双曲形偏微分方程的最典型代表,其最简形式可表示为：

关于位置X和时间t的标量函数U（代表各点偏离平衡位置的距离）满足：

这里c通常是一个固定輕代表波的传播速率。

在常压、20°C的空气中C为343米/秒（参见音速）o在弦振动问题中，c依不同弦的密度大小和轴向张力不同可能相差非常大。

而在半环螺旋弹簧（一种玩具，英文商标为Slinky）上,波速可以慢到1米/秒。

在针对实际问题的波动方程中，一般都将波速表示成可随波的频圣变化的量，这种处理对应真实物理世界中的色散现象。

此时，c应该用波的相速度代替:

实际问题中对标准波动方程的另一修正是考虑波速随輕的变化，修正后的方程变成下面的非线性波动方程：

92U2n

丽=c（u）V-u

另需注意的是物体中的波可能是叠加在其他运动（譬如介质的平动，以气流中传播的声波为例）上的。

这种情况下，标量u的表达式将包含一个马赫因子（对沿流动方向传播的波为正，对反射波为负）。

三维波动方程描述了波在均匀各向同性弹性体中的传播。

绝大多数固体都是弹性体，所以波动方程对地球内部的地震波和用于检测固体材料中缺陷的超声波的传播能给出满意的描述。

在只考虑线性行为时，三维波动方程的形式比前面更为复杂，它必须同时考虑固体中的纵波和横波:

pii=f+（入+2;z）V（V•u）-pVx（Vxu）

式中:

・入和“被称为弹性体的拉鯉敎（也叫“拉梅模量”，英文LdmQconstants或LamQmoduli）,是描述各向同性固体弹性性质的参数;

・P表示密度;

・f是源函数（即外界施加的激振力）；

•U表示位移；

注意在上述方程中，激振力和位移都是矢量,所以该方程也被称为矢量形式的波动方程。

其他形式的波动方程还能在量子力学和广义相对论理论中用到。

［隐藏］

.1标量形式的一维波动方程

O1.1波动方程的推导

O1.2初值问题的解

•2标量形式的三维波动方程

O2.1球面波

-2.1・1时间箭头的讨论

o2.2广义初值问题的解

-3标量形式的二维波动方程

-4边值问题

o4・1一维悄形

O4.2多维情形

•5注释

•6参考文献

•1参看

•8外部链接

标量形式的一维波动方程［编辑］

波动方程的推导［编辑］

一维波动方程可用如下的方式推导：

一列质量为也的小质点，相邻质点间用长度力的弹簧连接。

弹簧的弹性系数（乂称“倔强系数”）为层

u（x）u（x+h）u（x+2h）

其中u（x）表示位于X的质点偏离平衡位置的距离。

施加在位于X马处的质点加上的力为：

Fffaoke=凡+2/1+Fx=k\u（x+2ft,i）-u（x+ft,t）]+k\ii（x,t）—u（x+h,

其中Fwe诚“代表根据牛顿笫二定律计算的质点惯性力,FHooke代表根据胡克定律计算的弹簧作用力。

所以根据分析力学中的达朗贝尔原理,位于卅力处质点的运动方程为：

式中已注明ud）是时间t的显函数。

若艸个质点间隔均匀地固定在长度Z二Nh的弹簧链上，总质量朋二Nm,链的总体劲度系数为K=k/N,我们可以将上面的方程写为：

d2u（x+h、t）KL2u（x+2h>t）—2u（x+h』）+u{x}t）毎="m稚

取极限沖Tg,I。

就得到这个系统的波动方程:

d2u（x}t）KL2d2u{x}t）

汛2

M加2

C=\1—

在这个例子中，

波速VMo

初值问题的解［编辑］

一维标量形式波动方程的一般解是山达朗贝尔给出的。

原方程可以写成如下的算子作用形式：

dd］\dd］n

dt-cd^\\dt+cdi\u=^

从上面的形式可以看出，若尸和G为任意函数，那么它们以下形式的组合

=F（x—ct）+G（x+ct）

必然满足原方程。

上面两项分别对应两列行波（"行"与在"行动"中同音）一一尸表示经过该点（x点）的右行波，G表示经过该点的左行波。

为完全确定尸和G的最终形式还需考虑如下初始条件：

“3,0）=f（x）

=g（x）

经带入运算，就得到了波动方程著名的达朗贝尔行波解，乂称达朗贝尔公式:

在经典的意义下，如果€C"并且g

（2）€,x）€竺但是,

行波函数尸和G也可以是广义函数,比如狄拉克6函数。

在这种悄况下，行波解应被视作左行或右行的一个业。

基本波动方程是一个线性微分方程,也就是说同时受到两列波作用的点的振幅就是两列波振幅的相加。

这意味着可以通过把一列波分解成它的许求解中很有效。

标量形式的三维波动方程［编辑］

三维波动方程初值问题的解可以通过求解球面波波动方程得到。

求解结果可用于推导二维情况的解。

球面波［编辑］

球面波方程的形式不随空间坐标系统的转动而变化，所以可以将它写成仅与距源点距离r相关的函数。

方程的三维形式为:

将方程变形为：

（™）tt-c?

（ru）rr=0;

此时，因变量ru满足一维波动方程，于是可以利用达朗贝尔行波法将解写成:

讥匸r）=~F（r-

其中尸和G为任意函数，可以理解为以速度c从中心向外传播的波和从外面向中心传播的波。

这类从点源传岀的波强度随距点源距离r衰减，并且属于无后效波，可以清晰地搭载信号。

这种波仅在餾维空间中存在（原因将在下一小节中详细解释）。

幸运的是，我们生活的空间是三维的，所以我们可以清晰地通过声波和电磁波（都属于球面波）来互相交流。

时间箭头的讨论［编辑］

上面方程的解里面，分成了两部分，一部分表示向外传播的波，一部分则是向内。

很明显，只要将t换成-1,就可以在这两部分之间转换。

这体现了原始方程对于时间是对称的，任意的一个解在时间轴上倒过来看仍然是一个解。

然而，我们所观察到的实际的波，都是属于向外传播的。

除非精心地加以调整，我们无法在自然界观察到向内的波，尽管它们也是波动方程的合法的解。

关于这个现象，引起了不少讨论。

有人认为，实际上它们即使存在，也无法加以观察。

想想如果四周的光向一个物体集中，则因为没有光到达我们的眼睛，我们不可能看见这个物体或者发现这个现象（见参考文献［2］）。

广义初值问题的解［编辑］

波动方程中“是线性函数，并且不随时间和空间坐标的平移而改变。

所以我们可以通过平移与叠加球面波获得方程各种类型的解。

令为任意具有三个自变量的函数，球面波形尸为狄拉克§函数（数学语言是:

尸是一个在全空间积分等于1且非零区间收缩至原点的连续函数的過极陛）。

设（Jn,G位一族球面波的源点，r为距源点的径向距离，即：

r2=（⑦一£）2+（妙一〃）2+（2—02-

可定义

称为三维波动方程的影响函数，其意义为（^,n,O点在t=o时刻受到短促脉冲§函数作用后向空间中传出的波的影响，系数超4nc是为方便后续处理而加上的。

若U是这一族波函数的加权叠加，且权函数为4）,则

从各函数的定义可知，"还能写成

"（佔g

人2）=命〃卩（金+cto,y+龙丁）血，

式中C1、0和的意义为：

的t倍：

Y是单位球面S上点的坐标，d3为S上的面积微元。

该结果x,y,z）是以（$y,z）为圆心，^为半径的球面上"的平均值

迩卫』z）=

从上式易得

班0卫』2）=0,如0,禺皿）=0（禺盼）・

平均值是关于t的偶函数,所以若

v{t,x.y}z）=缶（叫砌），

那么

v（0,j,y,z）=讽莎匕2）,仇（0卫崩“）=0.

以上得岀的便是波动方程初值问题的解。

从中可以看岀，任意点尸在广时刻受到的波扰动只来自以尸为圆心，cr为半径的球面上，而这个球的内部点在这一吐刻对尸点的状态完全没有影响（因为它们的影响之前就已经传过尸点了）。

换一个角度分析，假设三维空间中任意点P在t=0时刻受到一个脉冲扰动6,那么山此发出的球面波在传过空间中的任意其它点Q后，便再也不会对Q的运动状态产生影响，这就是在物理学中也非常著名的惠更斯原理（Huygens'principle），也称为无后效现象，表示传过的球面波不会留下任何后续效应。

下面我们便可以解释上一小节中留下的问题了。

事实上，前面所得到的球面波解仅在奇数维空间中存在。

偶数维空间中波动方程的解是弥散的,也就是说波阵面掠过区域仍然会受其影响。

以下面的二维波动方程（极坐标形式,注意和上一小节三维形式的差别）为例：

_c2{urr+—ur）=0

T

可以从三维形式的解通过隆继法得到二维波动方程的影响函数：

其中

F=丿（工一02+尸

设点Masy）到点（4,n）距离为丛那么从影响函数中可以看出，当t>d/c即初始扰动已传过於点后，於仍在受到它的影响。

二维球面波（柱面波）的这一性质决定了它不能作为传递信号的工具，因为这种波（事实上包括所有偶数维空间中的球面波）经过的点受到的是交织在一起的各个不同时刻的扰动。

标量形式的二维波动方程［编蚩］

二维波动方程的直角坐标形式为:

如前所述，我们可以从三维波动方程的解中将U视为与其中一个自变量无关（降维法）来得到二维形式的解。

将初始条件改写为

饥（0角0=o,弘（0卫為）=0匕讪,

则三维形式的解就变成

其中Cl和B是单位球面上点的头两个坐匡分量，d3是球面上的面积微元。

此积分可变换为在（x,y）为中心，ct为半径的圆域。

上的积分：

从这个结果也能得到上一小节最后的结论。

二维波动方程解的一个例子是紧绷的鼓面的运动。

边值问题［编辑］

一根自身绷紧，两端分别固定于沪0和沪£的弹性弦在乃>0时刻，0<-y

在边界点处，可以要求u满足各种边界条件。

通常遇到的边界条件都可归纳成下列形式：

一血仏0）+au（i,0）=0,

ux（t,L）+bu（t,L）=0,

其中/b非负。

若要弦的两端固定不动，对应上面式子中/b趋于无穷大。

求解偏微分方程的分离_变量迭要求寻找以下形式的解：

u（t,x）=T（t）v（x）.

将上述假设形式代入原方程中可以得到:

为使边值问题有非平凡解，本征值X须满足

v1-FAf=0,

—d（0）+av（0）=0,d（Z）+bv（L）=0.

这是固有值问题的斯图姆-刘维尔理论的一个特例。

若日、b为正数，则对应的所有本征值均为正数，方程的解为三角函数。

使口和仏满足平方可积条件的解可以通过适当选取U和4三角级数展开来求得。

多维情形［编辑］

一维初始值-边值理论可以拓展至任意维空间中。

考虑也维空间（坐标简写为X）中的域0,万为。

的边界。

当0〈方时，位于0内的点x满足波动方程。

在0的边界上，解u须满足

其中”是万上指向域外的法向矢量，a是定义在B上的非负函数。

要求u在万上始终为0的边界条件相当于令$趋于无穷。

初始条件为

饥（0应）=/（X）,ut=g（>）,

其中f和呂是定义在0内的函数。

这个问题可以通过将f和g展开成域0内拉普拉斯算子满足边界条件的本征函数系的輕加来求解（这是分离变量法的一般步骤）。

也就是求解在域。

内满足

V・Vv+Av=0,

在边界万上满足

的本征函数系Fo

在二维情形下，上述本征函数系可以理解成绷紧地张在边界万上的鼓面的自山振动模态。

若万是一个圆，则这些本征函数是关于极角自变量0的三角函数与关于极轴自变量r的整阶贝塞尔函数的乘积。

更详细的说明参见英文版条LI玄姆霍兹方程。

在三维形式下，若边界是空间中的球面，那么本征函数是关于球坐标下两个极角自变量的球面调和函数,乘以关于径向自变量Q的半奇数阶贝塞尔函数。

注释［编辑］

参考文献［编辑］

•[1]严镇军编，《数学物理方程》，第二版，中国科学技术大学出版社，合肥，2002,第210页'第224页，ISBN7-312-00799-6/0•177

•[2][英]胡•普赖斯著，肖巍译，《时间之矢与阿基米德之点一物理学时间的新方向》，上海科学技术出版社，上海，2001,ISBN7-5323-5737-6

•[3jM.F・Atiyah,R・Bott,L・Garding,LacunasforhyperbolicdifferentialoperatorswithconstantcoefficientsI、ActaMath・,124（1970）,109-189.

•[4]M.F.Atiyah,R・Bott,andL・Garding,LacunasforhyperbolicdifferentialoperatorswithconstantcoefficientsII,ActaMath・,131（1973）,145-206.

•[5]R・Courant,D・Hilbert,MethodsofMathematicalPhysics,volII.Interscience（Wiley）NewYork,1962.