波动方程或称波方程

波动方程波动方程或称波方程波方程 英语 wavewave equationequation 是一种重要的偏微分方程 主要 描述自然界中的各种的波动现象 包括横波和纵波 例如声波 光波 无线电 波和水波 波动方程抽象自声学 物理光学 电磁学 电动力学 流体力学等 领域 历史上许多科学家 如达朗贝尔 欧拉 丹尼尔 伯努利和拉格朗日等在研究 乐器等物体中的弦振动问题时 都对波动方程理论作出过重要贡献 波动方程是双曲形偏微分方程的最典型代表 其最简形式可表示为 关于位置 x 和时间t 的标量函数u 代表各点偏离平衡位置的距离 满足 这里 c c 通常是一个固定常数 代表波的传播速率 在常压 20 C 的空气中c 为 343 米 秒 参见音速 在弦振动问题中 c 依不同弦的密度大小和轴向张 力不同可能相差非常大 而在半环螺旋弹簧 一种玩具 英文商标为 Slinky 上 波速可以慢到 1 米 秒 在针对实际问题的波动方程中 一般都将波速表示成可随波的频率变化的量 这种处理对应真实物理世界中的色散现象 此时 c 应该用波的相速度代替 实际问题中对标准波动方程的另一修正是考虑波速随振幅的变化 修正后的方 程变成下面的非线性波动方程非线性波动方程 另需注意的是物体中的波可能是叠加在其他运动 譬如介质的平动 以气流中 传播的声波为例 上的 这种情况下 标量u 的表达式将包含一个马赫因子 对沿流动方向传播的波为正 对反射波为负 三维波动方程三维波动方程描述了波在均匀各向同性弹性体中的传播 绝大多数固体都是弹 性体 所以波动方程对地球内部的地震波和用于检测固体材料中缺陷的超声波 的传播能给出满意的描述 在只考虑线性行为时 三维波动方程的形式比前面 更为复杂 它必须同时考虑固体中的纵波和横波 式中 和 被称为弹性体的拉梅常数 也叫 拉梅模量 英文 Lam constants 或 Lam moduli 是描述各向同性固体弹性性质的参数 表示密度 是源函数 即外界施加的激振力 表示位移 注意在上述方程中 激振力和位移都是矢量 所以该方程也被称为矢量形式的矢量形式的 波动方程波动方程 其他形式的波动方程还能在量子力学和广义相对论理论中用到 目录目录 隐藏 1 标量形式的一维波动方程 o1 1 波动方程的推导 o1 2 初值问题的解 2 标量形式的三维波动方程 o2 1 球面波 2 1 1 时间箭头的讨论 o2 2 广义初值问题的解 3 标量形式的二维波动方程 4 边值问题 o4 1 一维情形 o4 2 多维情形 5 注释 6 参考文献 7 参看 8 外部链接 标量形式的一维波动方程标量形式的一维波动方程 编辑编辑 波动方程的推导波动方程的推导 编辑编辑 一维波动方程可用如下的方式推导 一列质量为m的小质点 相邻质点间用长 度h的弹簧连接 弹簧的弹性系数 又称 倔强系数 为k 其中u x 表示位于x的质点偏离平衡位置的距离 施加在位于x h 处的质点 m 上的力为 其中代表根据牛顿第二定律计算的质点惯性力 代表根据胡克 定律计算的弹簧作用力 所以根据分析力学中的达朗贝尔原理 位于x h 处质 点的运动方程为 式中已注明u x 是时间t 的显函数 若N 个质点间隔均匀地固定在长度L N h 的弹簧链上 总质量M N m 链 的总体劲度系数为K k N 我们可以将上面的方程写为 取极限 N h 就得到这个系统的波动方程 在这个例子中 波速 初值问题的解初值问题的解 编辑编辑 一维标量形式波动方程的一般解是由达朗贝尔给出的 原方程可以写成如下的 算子作用形式 从上面的形式可以看出 若F 和G 为任意函数 那么它们以下形式的组合 必然满足原方程 上面两项分别对应两列行波 行 与在 行动 中同音 F 表示经过该点 x 点 的右行波 G 表示经过该点的左行波 为完全确定F 和 G 的最终形式还需考虑如下初始条件 经带入运算 就得到了波动方程著名的达朗贝尔行波解达朗贝尔行波解 又称达朗贝尔公式 在经典的意义下 如果并且则 但是 行波函数F和G 也可以是广义函数 比如狄拉克 函数 在这种情况下 行 波解应被视作左行或右行的一个脉冲 基本波动方程是一个线性微分方程 也就是说同时受到两列波作用的点的振幅 就是两列波振幅的相加 这意味着可以通过把一列波分解成它的许求解中很有 效 标量形式的三维波动方程标量形式的三维波动方程 编辑编辑 三维波动方程初值问题的解可以通过求解球面波波动方程得到 求解结果可用 于推导二维情况的解 球面波球面波 编辑编辑 球面波方程的形式不随空间坐标系统的转动而变化 所以可以将它写成仅与距 源点距离r 相关的函数 方程的三维形式为 将方程变形为 此时 因变量ru 满足一维波动方程 于是可以利用达朗贝尔行波法将解写成 其中F 和G 为任意函数 可以理解为以速度c 从中心向外传播的波和从外面 向中心传播的波 这类从点源传出的波强度随距点源距离r 衰减 并且属于无无 后效波后效波 可以清晰地搭载信号 这种波仅在奇数维空间中存在 原因将在下一 小节中详细解释 幸运的是 我们生活的空间是三维的 所以我们可以清晰 地通过声波和电磁波 都属于球面波 来互相交流 时间箭头的讨论时间箭头的讨论 编辑编辑 上面方程的解里面 分成了两部分 一部分表示向外传播的波 一部分则是向 内 很明显 只要将 t 换成 t 就可以在这两部分之间转换 这体现了原始方 程对于时间是对称的 任意的一个解在时间轴上倒过来看仍然是一个解 然而 我们所观察到的实际的波 都是属于向外传播的 除非精心地加以调整 我们无法在自然界观察到向内的波 尽管它们也是波动方程的合法的解 关于这个现象 引起了不少讨论 有人认为 实际上它们即使存在 也无法加 以观察 想想如果四周的光向一个物体集中 则因为没有光到达我们的眼睛 我们不可能看见这个物体或者发现这个现象 见参考文献 2 广义初值问题的解广义初值问题的解 编辑编辑 波动方程中u 是线性函数 并且不随时间和空间坐标的平移而改变 所以我们 可以通过平移与叠加球面波获得方程各种类型的解 令 为任意 具有三个自变量的函数 球面波形F 为狄拉克 函数 数学语言是 F 是一 个在全空间积分等于 1 且非零区间收缩至原点的连续函数的弱极限 设 位一族球面波的源点 r 为距源点的径向距离 即 可定义 称为三维波动方程的影响函数影响函数 其意义为 点在t 0 时刻受到短促脉 冲 函数作用后向空间中传出的波的影响 系数分母 4 c 是为方便后续处理 而加上的 若u 是这一族波函数的加权叠加 且权函数为 则 从 函数的定义可知 u 还能写成 式中 和 是单位球面S 上点的坐标 d 为S 上的面积微元 该结 果的意义为 u t x y z 是以 x y z 为圆心 ct 为半径的球面上 的平 均值的t 倍 从上式易得 平均值是关于t 的偶函数 所以若 那么 以上得出的便是波动方程初值问题的解 从中可以看出 任意点P 在t 时刻受 到的波扰动只来自以P 为圆心 ct 为半径的球面上 而这个球的内部点在这 一时刻对P 点的状态完全没有影响 因为它们的影响之前就已经传过P 点了 换一个角度分析 假设三维空间中任意点P 在t 0 时刻受到一个脉冲扰动 那么由此发出的球面波在传过空间中的任意其它点Q 后 便再也不会对Q 的运动状态产生影响 这就是在物理学中也非常著名的惠更斯原理惠更斯原理 Huygens principle 也称为无后效现象无后效现象 表示传过的球面波不会留下任何后续效应 下面我们便可以解释上一小节中留下的问题了 事实上 前面所得到的球面波 解仅在奇数维空间中存在 偶数维空间中波动方程的解是弥散的 也就是说波 阵面掠过区域仍然会受其影响 以下面的二维波动方程 极坐标形式 注意和 上一小节三维形式的差别 为例 可以从三维形式的解通过降维法得到二维波动方程的影响函数 其中 设点M x y 到点 距离为d 那么从影响函数中可以看出 当t d c 即初始扰动已传过M 点后 M 仍在受到它的影响 二维球面波 柱面波 的这 一性质决定了它不能作为传递信号的工具 因为这种波 事实上包括所有偶数 维空间中的球面波 经过的点受到的是交织在一起的各个不同时刻的扰动 标量形式的二维波动方程标量形式的二维波动方程 编辑编辑 二维波动方程的直角坐标形式为 如前所述 我们可以从三维波动方程的解中将u 视为与其中一个自变量无关 降维法 来得到二维形式的解 将初始条件改写为 则三维形式的解就变成 其中 和 是单位球面上点的头两个坐标分量 d 是球面上的面积微元 此积分可变换为在 x y 为中心 ct 为半径的圆域D 上的积分 从这个结果也能得到上一小节最后的结论 二维波动方程解的一个例子是紧绷的鼓面的运动 边值问题边值问题 编辑编辑 一维情形一维情形 编辑编辑 一根自身绷紧 两端分别固定于x 0 和x L 的弹性弦在t 0 时刻 0 x L 上运动满足波动方程 在边界点处 可以要求u 满足各种边界条件 通常遇到 的边界条件都可归纳成下列形式 其中a b 非负 若要弦的两端固定不动 对应上面式子中a b 趋于无穷大 求解偏微分方程的分离变量法要求寻找以下形式的解 将上述假设形式代入原方程中可以得到 为使边值问题有非平凡解非平凡解 本征值 须满足 这是固有值问题的斯图姆 刘维尔理论的一个特例 若a b 为正数 则对应的 所有本征值均为正数 方程的解为三角函数 使u 和ut 满足平方可积条件的 解可以通过适当选取u 和ut 三角级数展开来求得 多维情形多维情形 编辑编辑 一维初始值 边值理论可以拓展至任意维空间中 考虑m 维空间 坐标简写为 x 中的域D B为D的边界 当 0 t 时 位于D 内的点x 满足波动方程 在D 的边界上 解u 须满足 其中n 是B 上指向域外的法向矢量 a 是定义在B 上的非负函数 要求u 在 B 上始终为 0 的边界条件相当于令a 趋于无穷 初始条件为 其中f 和g 是定义在D 内的函数 这个问题可以通过将f 和g 展开成域D 内 拉普拉斯算子满足边界条件的本征函数系的叠加来求解 这是分离变量法的一 般步骤 也就是求解在域D 内满足 在边界B 上满足 的本征函数系v 在二维情形下 上述本征函数系可以理解成绷紧地张在边界B 上的鼓面的自由 振动模态 若B 是一个圆 则这些本征函数是关于极角自变量 的三角函数 与关于极轴自变量r 的整阶贝塞尔函数的乘积 更详细的说明参见英文版条目 亥姆霍兹方程 在三维形式下 若边界是空间中的球面 那么本征函数是关于球坐标下两个极 角自变量的球面调和函数 乘以关于径向自变量 的半奇数阶贝塞尔函数 注释注释 编辑编辑 参考文献参考文献 编辑编辑 1 严镇军编 数学物理方程 第二版 中国科学技术大学出版社 合肥 2002 第 210 页 第 224 页 ISBN 7 312 00799 6 O 177 2 英 胡 普赖斯著 肖巍译 时间之矢与阿基米德之点 物理学 时间的新方向 上海科学技术出版社 上海 2001 ISBN 7 5323 5737 6 3 M F Atiyah R Bott L Garding Lacunas for hyperbolic differential operators with constant coefficients I Acta Math