波动方程的行波解

波动方程的行波解

我们曾经讨论过波动方程 ∂ 2 u ∂ t 2 − a 2 ∂ 2 u ∂ x 2 = 0 \frac{\partial^{2} u}{\partial t^{2}}-a^{2} \frac{\partial^{2} u}{\partial x^{2}}=0 t22ua2x22u=0的通解,这里,把它改写成
u ( x , t ) = f ( x − a t ) + g ( x + a t ) u(x, t)=f(x-a t)+g(x+a t) u(x,t)=f(xat)+g(x+at)
其中 f f f g g g 是任意二阶可微函数。这个解式表明,波动方程的通解由两个波组成: f ( x − a t ) f(x-a t) f(xat) 代表沿 x x x 轴向右传播的波,当 t = 0 t=0 t=0 时,波形为 f ( x ) f(x) f(x),而后以恒定速率 a a a 向右传播,而保持波形不变; g ( x + a t ) g(x+a t) g(x+at) 则代表沿 x x x 轴向左传播的波,当 t = 0 t=0 t=0 时,波形为 g ( x ) g(x) g(x),而后也以同样的恒定速率 a a a 向左传播,保持波形不变。单独的 f ( x − a t ) f(x-a t) f(xat) g ( x + a t ) g(x+a t) g(x+at) 都是波动方程的解。它们独立传播,互不干扰,这正是因为波动方程是线性齐次方程,其解具有叠加性。

原则上说,函数 f f f g g g 应该由定解条件确定,但如果把问题简化为一维无界弦上的波的传播问题,那么 f f f g g g 便完全由初始条件决定。
{ u u = a 2 u x x ( − ∞ < x < + ∞ , t > 0 ) u ( x , 0 ) = φ ( x ) ( − ∞ < x < + ∞ ) u t ( x , 0 ) = ψ ( x ) ( − ∞ < x < + ∞ ) \left\{\begin{array}{lll} u_{u}=a^{2} u_{x x} & (-\infty<x<+\infty, t>0) \\ u(x, 0)=\varphi(x) & (-\infty<x<+\infty) \\ u_{t}(x, 0)=\psi(x) & (-\infty<x<+\infty) \end{array}\right. uu=a2uxxu(x,0)=φ(x)ut(x,0)=ψ(x)(<x<+,t>0)(<x<+)(<x<+)
其中 φ ( x ) , ψ ( x ) \varphi(x), \psi(x) φ(x),ψ(x) 为已知函数,“无限长”杆的自由振动,或电阻为零的 “无限长”传输线上电流、电压的变化,都可以提出相同的定解问题。我们已经得到了波动方程的通解 u ( x , t ) = f ( x − a t ) + g ( x + a t ) u(x, t)=f(x-a t)+g(x+a t) u(x,t)=f(xat)+g(x+at),接下来则需将初始条件带入求出函数 φ ( x ) \varphi(x) φ(x) ψ ( x ) \psi(x) ψ(x)
f ( x ) + g ( x ) = ϕ ( x ) a [ f ′ ( x ) − g ′ ( x ) ] = − ψ ( x ) \begin{gathered} f(x)+g(x)=\phi(x) \\ a\left[f^{\prime}(x)-g^{\prime}(x)\right]=-\psi(x) \end{gathered} f(x)+g(x)=ϕ(x)a[f(x)g(x)]=ψ(x)
a [ f ′ ( x ) − g ′ ( x ) ] = − ψ ( x ) a\left[f^{\prime}(x)-g^{\prime}(x)\right]=-\psi(x) a[f(x)g(x)]=ψ(x)积分,可以得到
f ( x ) − g ( x ) = − 1 a ∫ 0 x ψ ( ξ ) d ξ + C f(x)-g(x)=-\frac{1}{a} \int_{0}^{x} \psi(\xi) \mathrm{d} \xi+C f(x)g(x)=a10xψ(ξ)dξ+C
其中 C C C 是积分常数。将这个结果和 f ( x ) + g ( x ) = ϕ ( x ) f(x)+g(x)=\phi(x) f(x)+g(x)=ϕ(x)联立,即可求得
f ( x ) = 1 2 ϕ ( x ) − 1 2 a ∫ 0 x ψ ( ξ ) d ξ + C 2 , g ( x ) = 1 2 ϕ ( x ) + 1 2 a ∫ 0 x ψ ( ξ ) d ξ − C 2 . f(x)=\frac{1}{2} \phi(x)-\frac{1}{2 a} \int_{0}^{x} \psi(\xi) \mathrm{d} \xi+\frac{C}{2}, \quad g(x)=\frac{1}{2} \phi(x)+\frac{1}{2 a} \int_{0}^{x} \psi(\xi) \mathrm{d} \xi-\frac{C}{2} . f(x)=21ϕ(x)2a10xψ(ξ)dξ+2C,g(x)=21ϕ(x)+2a10xψ(ξ)dξ2C.
因此波动方程的解为
u ( x , t ) = f ( x − a t ) + g ( x + a t ) = 1 2 [ ϕ ( x − a t ) + ϕ ( x + a t ) ] + 1 2 a ∫ x − a t x + a t ψ ( ξ ) d ξ \begin{aligned} u(x, t) &=f(x-a t)+g(x+a t) \\ &=\frac{1}{2}[\phi(x-a t)+\phi(x+a t)]+\frac{1}{2 a} \int_{x-a t}^{x+a t} \psi(\xi) \mathrm{d} \xi \end{aligned} u(x,t)=f(xat)+g(x+at)=21[ϕ(xat)+ϕ(x+at)]+2a1xatx+atψ(ξ)dξ
这样,就求得了一维无界区间上波动方程定解问题的解。它称为一维波动方程定解问题的行波解,或达朗贝尔(d’Alembert) 解。这个解具有清楚的物理意义: 第一项表示由初位移 u ( x , t ) ∣ t = 0 = ϕ ( x ) \left.u(x, t)\right|_{t=0}=\phi(x) u(x,t)t=0=ϕ(x) 激发的行波, t > 0 t>0 t>0 以后分成相等的两部分,独立地向左、右传播,速率为 a a a;第二项表示由速度 ∂ u / ∂ t ∣ t = 0 = ψ ( x ) \partial u /\left.\partial t\right|_{t=0}=\psi(x) u/tt=0=ψ(x) 激发的行波[干涉项],在 t > 0 t>0 t>0 时刻,它左右对称地扩展到 [ x − a t , x + a t ] [x-a t, x+a t] [xat,x+at] 的范围,所以传播速率也是 a a a

d’Alembert 解法的思路容易理解,先求出通解,然后从中挑选特解。但是对一般偏微分方程而言,这是十分困难的,所以这种方法,并不是求解数理方程的常用方法,因此我们决定引入另一种解法–分离变量法。

下一章节:分离变量法

更多相关内容传送门:数学物理方法专栏阅读指南

posted @   Moster2469  阅读(550)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· winform 绘制太阳,地球,月球 运作规律
· 超详细:普通电脑也行Windows部署deepseek R1训练数据并当服务器共享给他人
· 上周热点回顾(3.3-3.9)
· AI 智能体引爆开源社区「GitHub 热点速览」
· 写一个简单的SQL生成工具
点击右上角即可分享
微信分享提示