波动方程的行波解
波动方程的行波解
我们曾经讨论过波动方程
∂
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
∂t2∂2u−a2∂x2∂2u=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(x−at)+g(x+at)
其中
f
f
f 和
g
g
g 是任意二阶可微函数。这个解式表明,波动方程的通解由两个波组成:
f
(
x
−
a
t
)
f(x-a t)
f(x−at) 代表沿
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(x−at) 和
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(x−at)+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)=−a1∫0xψ(ξ)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)−2a1∫0xψ(ξ)dξ+2C,g(x)=21ϕ(x)+2a1∫0xψ(ξ)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(x−at)+g(x+at)=21[ϕ(x−at)+ϕ(x+at)]+2a1∫x−atx+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/∂t∣t=0=ψ(x) 激发的行波[干涉项],在
t
>
0
t>0
t>0 时刻,它左右对称地扩展到
[
x
−
a
t
,
x
+
a
t
]
[x-a t, x+a t]
[x−at,x+at] 的范围,所以传播速率也是
a
a
a。
d’Alembert 解法的思路容易理解,先求出通解,然后从中挑选特解。但是对一般偏微分方程而言,这是十分困难的,所以这种方法,并不是求解数理方程的常用方法,因此我们决定引入另一种解法–分离变量法。
下一章节:分离变量法
更多相关内容传送门:数学物理方法专栏阅读指南
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· winform 绘制太阳,地球,月球 运作规律
· 超详细:普通电脑也行Windows部署deepseek R1训练数据并当服务器共享给他人
· 上周热点回顾(3.3-3.9)
· AI 智能体引爆开源社区「GitHub 热点速览」
· 写一个简单的SQL生成工具