JOVANA
Library Glossary Getting Started Three Levels Fields How it works Mission
Join the mission
All guides

首批应用:振动的弦与受热的细棒

傅里叶级数在这里不再是代数练习,而开始求解真正的物理:一根被拨动而歌唱的弦,一根渐渐冷却的暖棒,二者都被同一个「把谐波加起来」的想法撬开。这就是从傅里叶级数通向偏微分方程的那座具体桥梁。

两个日常物件,各藏着一个方程

把一根吉他弦在两端固定,拉到一边,再松手——它便弹成一片歌唱的颤动。另外,给一根金属棒的中段加热,然后走开——暖意四散,热点冷却,给足时间,整根棒就漂向一个平坦的温度。这看起来是两个截然不同的故事:一个在鸣响,一个在归于沉寂。然而它们受几乎相同的数学支配,而傅里叶级数正是把二者一并撬开的那一件工具。本篇就是这一阶梯那个抽象断言——「任何信号都是谐波之和」——终于在真正的物理上兑现价值的地方。

每个物件都藏着一个偏微分方程——它联系的不是普通导数,而是某个含两个变量(这里是位置 x 与时间 t)的未知函数的好几个偏导数振动弦服从波动方程,写作 d^2u/dt^2 = c^2 d^2u/dx^2(常缩写为 u_tt = c^2 u_xx),其中 u(x,t) 是弦在位置 x、时刻 t 处的横向位移,c 是由张力与质量定下的波速。受热细棒服从热方程,du/dt = k d^2u/dx^2(u_t = k u_xx),此处 u(x,t) 是温度,k 是热扩散率。两个物理世界,两个偏微分方程,差别只在一个微妙处:一阶时间导数对上二阶时间导数。

另有两个事实把每个问题钉死为唯一的答案。其一,边界条件:两端都被夹住。弦在 x = 0 与 x = L 处被系住,故那里不能动,u(0,t) = u(L,t) = 0;棒的两端被维持在零温,方程形式相同。其二,初始条件:你在 t = 0 时刻起步的形状——弦被拨动后的轮廓,或棒的初始温度分布。偏微分方程加上边界条件再加上初始状态,构成一个边值问题,而这整套恰好只有一个解。我们全部的活儿,就是把它找出来。

分离变量法:让傅里叶变得必不可少的引擎

核心动作是分离变量法:赌答案能分解成「只依赖空间的部分」乘以「只依赖时间的部分」,即 u(x,t) = X(x) T(t)。这是一个猜测,而非定律——它无法表示两个变量的每一个函数——但对这些干净的、线性的、带着如此简单边界的问题,它恰恰是对的那个猜测。把它代入偏微分方程,一桩小小的奇迹发生了:所有含 x 的东西堆到一边,所有含 t 的东西堆到另一边。

  1. 把 u = X(x) T(t) 代入偏微分方程。对棒而言:X T' = k X'' T。两边同除以 k X T,得到 T'/(k T) = X''/X。
  2. 论证分离常数。左边只依赖 t,右边只依赖 x,可它们对一切 x 与 t 都相等——所以二者必同等于某一个常数,记作 -lambda。一个含两个变量的偏微分方程,就此裂成两个常微分方程。
  3. 在夹紧端条件 X(0) = X(L) = 0 下求解空间方程 X'' + lambda X = 0。只有特殊值 lambda = (n pi / L)^2 才给出非零解,而那个解就是 X_n(x) = sin(n pi x / L)。是边界亲手挑出了这些谐波。
  4. 对每个 n 求解相应的时间方程,再把所有积木 u_n(x,t) = X_n(x) T_n(t) 叠加成一个无穷和,并用初始形状钉住各系数——而这恰恰就是一个傅里叶级数。

请细看第 3 步刚发生了什么,因为这就是整条阶梯的秘密所在。我们并没有「选择」使用正弦——是夹紧的边界把正弦强加给了我们。一个在 x = 0 与 x = L 处都为零的函数,必由 sin(n pi x / L) 搭成;余弦不肯,因为 cos(0) 不为零。边界条件替我们挑出了那个正交族,正是你早先学过其欧拉-傅里叶公式的那套正弦基。这就是傅里叶级数与这些偏微分方程密不可分的缘由:被夹紧物件的物理,字面意义上挑选了谐波,而傅里叶的配方是求出每种各需多少的唯一办法。

弦:鸣响着的驻波

现在把弦做完。它的空间部分如前推得 X_n(x) = sin(n pi x / L)。它的时间部分来自波动方程的两个时间导数,给出 T_n'' + (c n pi / L)^2 T_n = 0——一个简谐振子的方程。于是 T_n(t) 是一个在时间上以频率 omega_n = c n pi / L 摇摆的正弦余弦组合。每个积木 u_n = sin(n pi x / L) [A_n cos(omega_n t) + B_n sin(omega_n t)] 都是一个驻波:一个在空间上固定的正弦拱形,其高度在时间里上下脉动,从不行进,只在原地呼吸。

完整的运动把它们全部叠加:u(x,t) = sum over n of sin(n pi x / L)[A_n cos(omega_n t) + B_n sin(omega_n t)]。最慢的模式 n = 1 定下你听到的音高——弦的基频 omega_1 = c pi / L——而更快的模式 n = 2, 3, …… 是它的泛音,正是这些谐波让小提琴和长笛在同一个音上有了不同的音色。要定出各个量 A_n,令 t = 0 并要求形状匹配被拨动的轮廓 f(x):于是 f(x) = sum A_n sin(n pi x / L),这正是初始形状的傅里叶正弦级数。这些 A_n 就是它的傅里叶系数;B_n 以同样方式来自初始速度。

VIBRATING STRING -- length L, ends fixed, plucked into shape f(x), released from rest

  PDE:        u_tt = c^2 u_xx
  Boundary:   u(0,t) = u(L,t) = 0          (tied at both ends)
  Initial:    u(x,0) = f(x),  u_t(x,0) = 0  (plucked, then let go)

  Separate u = X(x) T(t)  ->   X'' + lambda X = 0 ,  T'' + c^2 lambda T = 0
  Boundaries force         X_n = sin(n pi x / L),  lambda_n = (n pi / L)^2
  Time part (released from rest) ->  T_n = A_n cos(omega_n t),  omega_n = c n pi / L

  Solution (superpose the standing waves):
      u(x,t) = sum_{n>=1} A_n sin(n pi x / L) cos(c n pi t / L)

  Coefficients = Fourier sine series of the pluck:
      A_n = (2/L) integral_0^L f(x) sin(n pi x / L) dx
一头到尾解完的弦。边界挑出正弦模式,每个模式以自己的频率 omega_n 鸣响,而被拨动形状的傅里叶正弦系数 A_n 说明每个谐波起始时多响。从静止释放消去了 B_n sin(omega_n t) 那些项,时间上只剩下余弦。

细棒:同样的模式,却是消逝而非鸣响

现在轮到棒,留意它在何处与弦分道扬镳。空间部分一模一样——X_n(x) = sin(n pi x / L),正是同样的正弦模式,因为夹紧的边界相同。唯一变的是时间方程。热方程只有一个时间导数,而非两个,所以分离给出 T_n' = -k (n pi / L)^2 T_n——一个一阶方程,其解是一个衰减的指数,T_n(t) = exp(-k (n pi / L)^2 t)。一个导数对上两个导数,就是音乐与沉寂之间的全部差别:弦的模式永远振荡,棒的模式只是消融而去。

叠加起来便得到 u(x,t) = sum over n of B_n sin(n pi x / L) exp(-k (n pi / L)^2 t)。读一读这个指数:衰减率是 k(n pi / L)^2,它随 n 的平方增长。所以高次谐波——那些尖锐、扭动、含精细细节的模式——死得最快,而缓慢的 n = 1 鼓包逗留最久。这正是热为何把东西抹平:每一处带棱带角的特征都是高谐波的,而高谐波几乎瞬间蒸发,早在棒变得通体冰冷之前,就已留下一根光滑而平淡的棒。要求出 B_n,再次令 t = 0:初始温度 f(x) 必等于 sum B_n sin(n pi x / L),所以 B_n 又一次是起始形状的傅里叶正弦系数。

它为何管用、在哪里失效、又通向何方

让我们能够把无穷多个驻波相加、却仍然求解偏微分方程的那条逻辑,就是叠加原理:因为两个方程都是线性的,所以任何解之和仍是一个解。我们孤立地求解单个谐波——这很容易,因为单一模式无非是一个振子或一个衰减指数——然后相信把它们全部叠加就给出真实的运动,以初始形状的傅里叶系数作为混合权重。「先解一个模式再叠加」,正是你在线性常微分方程那里见过的分而治之,如今被抬升到了空间与时间的函数上。傅里叶级数恰恰是负责「分」的那台机器。

退后一步,好好品味你如今能做的事。物理学两个伟大的方程——波动与热——都向同一套有纪律的流程屈服:分离变量,让边界把正弦模式递到你手上,求解每个模式那个小小的时间方程,再让初始形状的傅里叶系数为求和加权。弦之所以鸣响,是因为它的模式振荡;棒之所以抹平,是因为它的模式衰减;而区分二者的那一个数,就是时间究竟出现一次还是两次。从这里起,道路朝两个方向攀升。当几何不再是一段平直的区间——一面振动的鼓膜、一根圆柱中的热——正弦模式便让位给别的正交族,比如你见过的、基于特殊函数的那些展开,它们全都受同一个「投影到一组基上」的思想支配。而当物体是无限长而非被夹紧时,离散的谐波模糊成一片连续区,傅里叶级数便化为傅里叶变换。你方才亲眼看着一个抽象的级数,成为了物理世界的语言。