2. 一维波动方程的数学基础

有限差分仿真的本质，是把连续的偏微分方程“翻译”成离散网格上可逐点计算的代数递推式。翻译是否忠实，取决于我们对原方程数学结构的理解：方程是齐次还是非齐次、边界是固定还是自由、初始条件是位移型还是速度型、解是行波还是驻波，这些都会直接改变差分格式的构造方式与稳定性条件。本节围绕一维波动方程，梳理四组关键数学事实，并说明它们各自如何落到差分格式的设计上。

齐次与非齐次：方程右端项决定递推式是否多一项

取弦振动为物理原型，设 u(x,t) 为横向位移，波速为 c，则一维齐次波动方程为

∂²u/∂t² = c² · ∂²u/∂x²        (齐次，右端为 0)

若弦上还作用有分布外力 f(x,t)（单位质量对应的加速度量纲），则得到非齐次方程

∂²u/∂t² = c² · ∂²u/∂x² + f(x,t)   (非齐次)

两者的差别在离散后非常直观。采用空间步长 Δx、时间步长 Δt，记 u_iⁿ ≈ u(iΔx, nΔt)，对时间二阶导用中心差商、空间二阶导用中心差商，得到显式格式：

u[i]^(n+1) = 2u[i]^n − u[i]^(n−1) + r²·(u[i+1]^n − 2u[i]^n + u[i−1]^n) + Δt²·f[i]^n
其中 r = c·Δt/Δx 为库朗数

齐次情形下最后一项为零，递推只依赖上一时刻与当前时刻的位移；非齐次情形下每个节点每步都要叠加 Δt²·f_iⁿ，相当于把外力作为“已知激励”注入。工程上这意味着：源项采样必须与时间层对齐，若 f 随时间剧烈变化，Δt 还需同时满足源项的分辨率要求，而不能只按波速选取。

初始条件与边界条件：定解信息的三种类型

波动方程对时间二阶，因此需要两个初始条件；对空间二阶，因此两端各需一个边界条件。常见组合如下：

• 初始位移型：u(x,0)=φ(x)，u_t(x,0)=ψ(x)。若 ψ=0，表示弦从静止释放。
• 固定端（Dirichlet）：u(0,t)=u(L,t)=0，边界节点值恒为零，无需递推。
• 自由端（Neumann）：u_x(0,t)=0，需用虚拟节点或单侧差商消去边界外的未知量。
• 弹性/阻抗边界：u_x 与 u_t 的线性组合为零，对应吸收或反射条件。

差分实现时有两个易错点。其一，速度初始条件不能直接代入位移递推式，必须用虚拟时间层 n=−1 消去：由 u_t(x,0)=ψ 得 u_i⁻¹=u_i¹−2Δt·ψ_i，再代入 n=0 的递推式，得到只含 u¹ 与 u⁰ 的启动公式。其二，Neumann 边界若用一阶单侧差商，边界处精度会降到一阶，可能污染整体解；常用做法是引入界外虚拟节点 u₋₁=u₁，保持二阶精度。

行波解与分离变量解：两种视角对应两类验证手段

齐次方程的通解可由达朗贝尔公式给出：

u(x,t) = ½[φ(x−ct) + φ(x+ct)] + (1/2c)·∫[x−ct, x+ct] ψ(s) ds

它表明任意扰动都分解为右行波 φ(x−ct) 与左行波 φ(x+ct) 的叠加，波形在传播中不变。这一视角直接支撑了数值格式的两个设计判断：显式格式的依赖域必须覆盖特征线 x±ct，这正是库朗条件 r≤1 的几何来源；而数值色散会使不同波数分量以不同速度传播，导致波形逐渐“散开”，因此用单脉冲传播测试最容易暴露格式的色散误差。

分离变量解则给出另一幅图景。令 u=X(x)T(t)，代入齐次方程得 T″/(c²T)=X″/X=−k²，于是

X(x) = sin(kx) 或 cos(kx)，T(t) = sin(ωt) 或 cos(ωt)，ω = ck

在两端固定条件下，k_n=nπ/L，解为各阶驻波模态的叠加，频率 ω_n=nπc/L。这一视角的实用价值在于：它给出了精确的模态频率与振型，可作为数值结果的解析基准。例如取 L=1、c=1，则基频为 π rad/s；若数值解的一阶模态频率明显偏离该值，说明空间分辨率不足或时间推进存在相位误差。

把数学结论落到格式参数上

综合上述结构，可以给出一套可操作的参数选取与验证流程。以长度 L=1 m、波速 c=100 m/s、两端固定、初始位移为半正弦 φ(x)=sin(πx)、初始速度为零的示例为例（以下数值为教学设定，非实测）：

1. 空间分辨率：若关心到第 5 阶模态，最高波数对应波长 2L/5=0.4 m，按每波长至少 10 个节点，取 Δx=0.04 m，即 N=25 个区间。
2. 时间步长：由库朗条件取 r=0.9，得 Δt=rΔx/c=0.9×0.04/100=3.6×10⁻⁴ s。若存在源项 f，还需检查其最快变化周期是否被 Δt 充分采样。
3. 启动步：因 ψ=0，虚拟层给出 u_i⁻¹=u_i¹，代入 n=0 递推式可解出 u_i¹=u_i⁰+½r²(u_i+1⁰−2u_i⁰+u_i−1⁰)，无需迭代。
4. 结果判读：运行若干周期后，对中心点时间序列做频谱分析，一阶峰应落在 c/(2L)=50 Hz 附近；若峰值偏移超过 1%～2%，优先减小 Δx 或改用更高阶空间差分，而不是单纯减小 Δt。

归结起来，齐次与非齐次决定递推式是否含源项，初边值条件决定启动步与边界闭合方式，行波解给出稳定性与色散的几何判据，分离变量解提供可对照的模态基准。四者共同构成差分格式设计的理论依据，也是后续编写程序、设置参数和排查异常时的判断出发点。