d’Alembert方程(或称为波动方程)是描述一维波动现象的基本偏微分方程,形式如下:
\[ \frac{\partial^2 u}{\partial t^2} = c^2 \frac{\partial^2 u}{\partial x^2} \]
其中:
- \( u(x, t) \) 是波的位移函数,表示位置 \( x \) 和时间 \( t \) 时的振幅。
- \( c \) 是波速(常数,由介质的性质决定)。
d’Alembert 方程(波动方程)通常表示为: \[ \Box \phi = \left( \frac{1}{c^2} \frac{\partial^2}{\partial t^2} - \nabla^2 \right) \phi = 0 \] 其中:
- \(\Box\) 是 d’Alembert 算子(波动算子),
- \(\nabla^2\) 是 拉普拉斯算子(空间二阶导数),
- \(c\) 是波速(如光速)。
d’Alembert 方程: \[ \Box \phi = 0 \] 可以写成: \[ \left( -\frac{1}{c^2} \frac{\partial^2}{\partial t^2} + \frac{\partial^2}{\partial x^2} + \frac{\partial^2}{\partial y^2} + \frac{\partial^2}{\partial z^2} \right) \phi = 0 \]
在 球坐标系 \((r, \theta, \phi)\) 中,拉普拉斯算子 \(\nabla^2\) 为: \[ \nabla^2 = \frac{1}{r^2} \frac{\partial}{\partial r} \left( r^2 \frac{\partial}{\partial r} \right) + \frac{1}{r^2 \sin \theta} \frac{\partial}{\partial \theta} \left( \sin \theta \frac{\partial}{\partial \theta} \right) + \frac{1}{r^2 \sin^2 \theta} \frac{\partial^2}{\partial \phi^2} \] 因此,d’Alembert 方程在球坐标下可以写成: \[ \left( \frac{1}{c^2} \frac{\partial^2}{\partial t^2} - \nabla^2 \right) \phi = 0 \] 但 矩阵形式 在球坐标中较复杂,通常使用张量分析或微分几何的方法。
解的形式(d’Alembert解)
方程的通解由d’Alembert公式给出:
\[ u(x, t) = f(x - ct) + g(x + ct) \]
这里:
- \( f(x - ct) \) 表示以速度 \( c \) 向右传播的行波。
- \( g(x + ct) \) 表示以速度 \( c \) 向左传播的行波。
- \( f \) 和 \( g \) 是任意二阶可微函数,由初始条件或边界条件确定。
初始条件问题
若给定初始位移和初始速度: \[ u(x, 0) = \phi(x), \quad \frac{\partial u}{\partial t}(x, 0) = \psi(x), \] 则解为: \[ u(x, t) = \frac{1}{2} \left[ \phi(x - ct) + \phi(x + ct) \right] + \frac{1}{2c} \int_{x - ct}^{x + ct} \psi(s) , ds. \]
物理意义
该方程适用于:
- 弦的横向振动(如吉他弦)。
- 声波在均匀介质中的传播。
- 电磁波在真空或均匀空间中的传播(此时 \( c \) 为光速)。
高维推广
在三维空间中,波动方程变为: \[ \frac{\partial^2 u}{\partial t^2} = c^2 \nabla^2 u, \] 其中 \( \nabla^2 \) 是拉普拉斯算子。
数学性质
- 满足线性叠加原理(解的线性组合仍是解)。
- 具有行波解和驻波解(如通过分离变量法得到)。
如果需要更具体的推导或应用场景,可以进一步说明!
1. 从球面波方程出发
假设有一个从原点发出的球对称波,其波前满足: \[ (ct)^2 = x^2 + y^2 + z^2 \] 这表示在时间 \( t \),波前到达的位置 \((x, y, z)\) 满足该方程。
2. 定义波函数 \( u(x, y, z, t) \)
假设波函数 \( u \) 在波前上是常数(例如 \( u = \text{const} \) 在 \((ct)^2 = x^2 + y^2 + z^2\) 上),那么我们可以尝试构造一个满足波动方程的函数。
3. 考虑一般行波解
在三维空间中,球对称波的一般解可以表示为: \[ u(r, t) = \frac{f(r - ct)}{r} + \frac{g(r + ct)}{r} \] 其中 \( r = \sqrt{x^2 + y^2 + z^2} \),\( f \) 和 \( g \) 分别代表向外和向内传播的波。
4. 验证是否满足波动方程
三维波动方程为: \[ \frac{\partial^2 u}{\partial t^2} = c^2 \nabla^2 u \] 其中 \( \nabla^2 \) 是拉普拉斯算子: \[ \nabla^2 u = \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} + \frac{\partial^2 u}{\partial z^2} \]
计算时间二阶导数
对 \( u(r, t) = \frac{f(r - ct)}{r} \)(仅考虑向外传播的波): \[ \frac{\partial u}{\partial t} = \frac{-c f’(r - ct)}{r}, \quad \frac{\partial^2 u}{\partial t^2} = \frac{c^2 f’’(r - ct)}{r} \]
计算空间二阶导数(拉普拉斯算子)
由于 \( u \) 是球对称的,拉普拉斯算子在球坐标下为: \[ \nabla^2 u = \frac{1}{r^2} \frac{\partial}{\partial r} \left( r^2 \frac{\partial u}{\partial r} \right) \] 计算 \( \frac{\partial u}{\partial r} \): \[ \frac{\partial u}{\partial r} = \frac{f’(r - ct)}{r} - \frac{f(r - ct)}{r^2} \] 再计算 \( \nabla^2 u \): \[ \nabla^2 u = \frac{1}{r^2} \frac{\partial}{\partial r} \left( r f’(r - ct) - f(r - ct) \right) = \frac{f’’(r - ct)}{r} \]
对比波动方程
\[ \frac{\partial^2 u}{\partial t^2} = c^2 \nabla^2 u \implies \frac{c^2 f’’(r - ct)}{r} = c^2 \left( \frac{f’’(r - ct)}{r} \right) \] 两边相等,验证了 \( u(r, t) = \frac{f(r - ct)}{r} \) 满足三维波动方程。
5. 从 \((ct)^2 = x^2 + y^2 + z^2\) 到波动方程
更直接的方法是利用 特征曲面(characteristic surfaces) 的概念。波动方程的特征曲面满足: \[ (ct)^2 = x^2 + y^2 + z^2 \] 这表明波前是一个光锥(light cone),其法向量满足: \[ c^2 dt^2 - dx^2 - dy^2 - dz^2 = 0 \] 这对应于 相对论中的时空间隔,并暗示了波动方程的形式: \[ \frac{\partial^2 u}{\partial t^2} = c^2 \left( \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} + \frac{\partial^2 u}{\partial z^2} \right) \]
总结
- 从 \((ct)^2 = x^2 + y^2 + z^2\) 可以推导出三维波动方程,因为它描述了波前的传播。
- 球对称解 \( u(r, t) = \frac{f(r - ct)}{r} \) 满足波动方程。
- 更一般地,波动方程的解可以表示为行波的叠加。
如果需要更详细的数学步骤(如拉普拉斯算子的计算),可以进一步展开!
- \( r = \sqrt{x^2 + y^2 + z^2} \) 是径向距离,
- \( c \) 是波速,
- \( f \) 和 \( g \) 是任意(二阶可微)函数,分别代表 向外传播 和 向内传播 的波。
1. 物理意义
(1) 向外传播的波 \( \frac{f(r - ct)}{r} \)
- 含义:表示波从原点(\( r=0 \))向外辐射,例如:
- 点光源发出的光波,
- 爆炸产生的声波。
- 为什么是 \( r - ct \):
- \( f(r - ct) \) 表示波以速度 \( c \) 向外传播,\( r - ct = \text{常数} \) 是波前的方程。
- 例如,\( t=0 \) 时波在 \( r = r_0 \),经过时间 \( t \),波到达 \( r = r_0 + ct \)。
- 为什么除以 \( r \):
- 在三维空间中,波的能量分布在球面上,面积 \( \propto r^2 \),因此振幅必须 \( \propto \frac{1}{r} \) 以保持能量守恒(否则总能量会无限增长)。
(2) 向内传播的波 \( \frac{g(r + ct)}{r} \)
- 含义:表示波从无穷远处 向原点汇聚,例如:
- 聚焦的声波或光波,
- 黑洞的引力波反射(理论模型)。
- 为什么是 \( r + ct \):
- \( g(r + ct) \) 表示波以速度 \( c \) 向内传播,\( r + ct = \text{常数} \) 是波前的方程。
- 例如,\( t=0 \) 时波在 \( r = r_0 \),经过时间 \( t \),波到达 \( r = r_0 - ct \)(向原点靠近)。
- 物理合理性:
- 实际物理问题中,向内传播的波通常需要特定边界条件(如完美反射)或源项支持,否则较少出现。
2. 为什么解是这种形式?
(1) 波动方程在球坐标下的化简
三维波动方程: \[ \frac{\partial^2 u}{\partial t^2} = c^2 \nabla^2 u \] 在球对称情况下(\( u \) 仅依赖 \( r \) 和 \( t \)),拉普拉斯算子为: \[ \nabla^2 u = \frac{1}{r^2} \frac{\partial}{\partial r} \left( r^2 \frac{\partial u}{\partial r} \right) \] 令 \( v(r, t) = r u(r, t) \),则方程化简为: \[ \frac{\partial^2 v}{\partial t^2} = c^2 \frac{\partial^2 v}{\partial r^2} \] 这是一维波动方程,其通解为: \[ v(r, t) = f(r - ct) + g(r + ct) \] 因此: \[ u(r, t) = \frac{f(r - ct)}{r} + \frac{g(r + ct)}{r} \]
(2) 边界条件决定 \( f \) 和 \( g \)
- 初始条件:给定 \( u(r, 0) \) 和 \( \frac{\partial u}{\partial t}(r, 0) \),可确定 \( f \) 和 \( g \)。
- 物理约束:
- 在 \( r=0 \) 处,\( u \) 应为有限值,因此 \( f(-ct) = -g(ct) \)(否则分母为零时发散)。
- 若问题无内向波(如自由辐射),则 \( g = 0 \)。
3. 直观类比:一维 vs 三维
- 一维波动方程:
\[
u(x, t) = f(x - ct) + g(x + ct)
\]
- \( f \) 和 \( g \) 直接代表右行和左行波,无衰减。
- 三维波动方程:
\[
u(r, t) = \frac{f(r - ct)}{r} + \frac{g(r + ct)}{r}
\]
- 振幅随 \( \frac{1}{r} \) 衰减,反映能量分布在球面上的扩散。
4. 经典例子
(1) 向外传播的脉冲波
设初始扰动仅在 \( r = r_0 \): \[ u(r, 0) = \frac{\delta(r - r_0)}{r}, \quad \frac{\partial u}{\partial t}(r, 0) = 0 \] 解为: \[ u(r, t) = \frac{\delta(r - r_0 - ct)}{r} \quad \text(仅外向波) \]
(2) 驻波(内向+外向叠加)
若 \( f(r - ct) = -g(r + ct) \),则: \[ u(r, t) = \frac{f(r - ct) - f(r + ct)}{r} \] 在 \( r=0 \) 处有限(如球对称谐振腔内的声波)。
总结
- \( f(r - ct) \):外向波(物理常见,如光辐射、声波发散)。
- \( g(r + ct) \):内向波(需特定条件支持,如反射或聚焦)。
- 因子 \( \frac{1}{r} \) 体现三维几何导致的振幅衰减。
如果需要具体问题的计算(如给定初始条件求 \( f \) 和 \( g \)),可以进一步讨论!
1. 相对论约定(4维时空)
在 狭义相对论 或 场论 中,通常采用 四维坐标系 \((x^0, x^1, x^2, x^3)\),其中:
- \(x^0 = ct\)(时间坐标,\(c\) 为光速,\(t\) 为时间),
- \(x^1 = x\),\(x^2 = y\),\(x^3 = z\)(空间坐标)。
此时偏导数的定义为: \[ \partial_\mu \phi = \frac{\partial \phi}{\partial x^\mu}, \quad \mu = 0, 1, 2, 3 \]
具体形式
-
时间导数(\(\partial_0 \phi\)): \[ \partial_0 \phi = \frac{\partial \phi}{\partial x^0} = \frac{1}{c} \frac{\partial \phi}{\partial t} \]
- 多出的因子 \(1/c\) 是因为 \(x^0 = ct\)。
-
空间导数(\(\partial_i \phi\),\(i = 1, 2, 3\)): \[ \partial_1 \phi = \frac{\partial \phi}{\partial x}, \quad \partial_2 \phi = \frac{\partial \phi}{\partial y}, \quad \partial_3 \phi = \frac{\partial \phi}{\partial z} \]
- 这是直角坐标系下的标准空间偏导。
2. 非相对论约定(时间与空间分离)
在经典物理或某些非相对论性场论中,通常直接区分时间和空间:
- 时间坐标:\(t\),
- 空间坐标:\(x, y, z\)。
此时:
-
时间导数: \[ \partial_0 \phi = \frac{\partial \phi}{\partial t} \]
- 注意:此处 \(\partial_0\) 直接对应 \(t\),无 \(c\) 因子。
-
空间导数: \[ \partial_i \phi = \frac{\partial \phi}{\partial x^i}, \quad i = 1, 2, 3 \]
- 与相对论约定相同。
3. 一般坐标系中的表示
在曲线坐标系(如球坐标、柱坐标)中,\(\partial_i \phi\) 需用协变导数或坐标变换表示。例如:
- 球坐标 \((r, \theta, \varphi)\): \[ \partial_1 \phi = \frac{\partial \phi}{\partial r}, \quad \partial_2 \phi = \frac{1}{r} \frac{\partial \phi}{\partial \theta}, \quad \partial_3 \phi = \frac{1}{r \sin\theta} \frac{\partial \phi}{\partial \varphi} \]
4. 示例:Klein-Gordon 方程
在相对论场论中,Klein-Gordon 方程用 \(\partial_\mu\) 表示为: \[ (\partial^\mu \partial_\mu + m^2) \phi = 0 \] 展开后: \[ \frac{1}{c^2} \frac{\partial^2 \phi}{\partial t^2} - \nabla^2 \phi + m^2 \phi = 0 \] 其中 \(\nabla^2 \phi = \partial_i \partial_i \phi = \frac{\partial^2 \phi}{\partial x^2} + \frac{\partial^2 \phi}{\partial y^2} + \frac{\partial^2 \phi}{\partial z^2}\)。
5. 总结
- \(\partial_0 \phi\):
- 相对论约定:\(\frac{1}{c} \frac{\partial \phi}{\partial t}\),
- 非相对论约定:\(\frac{\partial \phi}{\partial t}\)。
- \(\partial_i \phi\):
- 直角坐标:\(\frac{\partial \phi}{\partial x^i}\),
- 曲线坐标:需考虑度规(如球坐标中的 \(1/r\) 因子)。
如果需要具体问题中的计算(如某一场方程的展开),可以进一步说明!