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

复指数形式、奇偶对称与半区间展开

把傅里叶级数重新打包成一个紧凑的复指数和,让函数的对称性替你省掉一半的活,再学会那个半区间技巧——它能把只定义在半个窗口上的函数,变成边值问题恰好需要的正弦或余弦级数。

一种紧凑形式:复指数级数

在上一篇里你构造了[[trigonometric-fourier-series|三角傅里叶级数]]:把一个周期函数写成 a_0/2 加上一串 a_n cos(n omega x) 与 b_n sin(n omega x) 之和,系数则由欧拉–傅里叶公式撬出来。它管用,但稍嫌臃肿——两族系数 a_n 与 b_n,还有两种要分清的项。有一种更利落的重新打包方式,物理学家和工程师几乎是本能地伸手去取,而它倚靠的是你在第一卷里已经见过的一个恒等式。

这个恒等式就是欧拉公式 e^{i theta} = cos(theta) + i sin(theta)。把它正反两用,你就能把每一个余弦和正弦都换成指数:cos(theta) = (e^{i theta} + e^{-i theta})/2,sin(theta) = (e^{i theta} - e^{-i theta})/(2i)。把这些代回三角级数,频率为 n 的那两项就坍缩成 +n 与 -n 处的指数。整个和重组成[[complex-exponential-fourier-series|复指数傅里叶级数]]:f(x) = 对所有整数 n(从负无穷到正无穷)求和 c_n e^{i n omega x}。只剩一族系数 c_n、一种项,而且指标现在也跑遍负整数。

系数公式相应地也很干净:c_n = (1/T) 乘以 f(x) e^{-i n omega x} dx 在一个周期上的积分,其中 T 是周期、omega = 2 pi / T。注意这是对每一个 n(正、负、零)都通用的同一个公式——c_0 不过是 f 在一个周期上的平均,即旧的 a_0/2。指数里的那个负号才是全部的秘密:e^{-i n omega x} 正是那个复指数,它凭正交性恰好挑出第 n 个分量,而把其余全部清零。我们下面就把这条正交性拆开来看。

负号为何奏效,以及 c_n 的含义

三角基之所以正交,是因为 cos 乘 sin、以及不匹配谐波的积分在一个周期上都为零。指数继承了同一事实的一个更干净的版本。e^{i n omega x} 乘 e^{-i m omega x} 在一个周期上的积分:当 n = m 时等于 T,当 n 不等于 m 时恰好为零——因为 e^{i(n-m) omega x} 是一个绕单位圆整数圈、平均为零的函数。这正是三角系正交性最省力的重述。于是用 e^{-i m omega x} 乘 f 再积分,就读出了那唯一的系数 c_m,并让其余每一项都噤声。

负指标在物理上是什么意思?没什么神秘的。第 n 项与第 -n 项总是成对出现,而当 f 是实信号时,它们互为复共轭:c_{-n} = c_n 的共轭。把这一对加回去,虚部相消,又长回频率为 n 的一个实的 cos 加 sin 振荡。回到旧系数的字典是:对正的 n,c_n = (a_n - i b_n)/2。所以负频率不过是实振荡那共轭一半的记账方式——为换来一个整洁的公式,这个代价值得付。

对称性白送你一半的工作量

在算出哪怕一个积分之前,先看看函数的对称性——它能当场把工作量减半。回想第一卷的词汇:若 f(-x) = f(x),函数是的(关于纵轴的镜像,像 cos 或 x^2);若 f(-x) = -f(x),则是的(绕原点的半圈旋转,像 sin 或 x^3)。一个关键事实——从定积分作为带符号面积来看很容易明白——是:奇函数在对称区间 -L 到 L 上的积分恰好为零(左半与右半相消),而偶函数的积分不过是右半积分的两倍。

现在把这一点喂进系数积分里。余弦是偶的、正弦是奇的,所以欧拉–傅里叶公式里那些乘积也继承了奇偶性。若 f 是偶函数,那么 f 乘 sin(n omega x) 是奇的,它在对称周期上的积分为零,于是每一个 b_n 都为零:偶函数是纯余弦级数。若 f 是奇函数,同样的论证消灭每一个 a_n,于是 f 是纯正弦级数。仅凭瞥一眼图像,你就能预先知道一半的系数为零——而那些积分你压根不用算。

半区间展开:造出你需要的对称性

在这里,对称性这个念头从一个省力技巧,摇身变为一个解题手段。在物理问题里,函数往往只在半个窗口上给出——比如一根杆的温度分布只对 0 <= x <= L 给定,对负的 x 则只字未提,因为杆根本不向那边延伸。此时还没有周期、没有左半、没有对称性。半区间展开这一招,就是*选择*如何把函数向左延拓,并且精确地这样选,以制造出让答案成为纯余弦或纯正弦的那种奇偶性。

两种选择,两种结果。把 f 反射成偶函数(关于 x = 0 镜像),再以周期 2L 重复:结果是偶的,正弦系数全部死去,于是得到半区间余弦级数,a_n = (2/L) 乘以 f(x) cos(n pi x / L) 从 0 到 L 的积分。换成把它反射成奇函数(绕原点旋转),再以周期 2L 重复:这下余弦死去,得到半区间正弦级数,b_n = (2/L) 乘以 f(x) sin(n pi x / L) 从 0 到 L 的积分。同一个在 [0, L] 上的 f,两条迥然不同的级数——在给你的那半个窗口上二者都正确。

Given f(x) on  [0, L]  only.

EVEN extension  ->  period 2L, mirror across x=0
   half-range COSINE series:
   f(x) = a0/2 + sum_{n>=1} a_n cos(n pi x / L)
   a_n = (2/L) * integral_0^L f(x) cos(n pi x / L) dx
   slope zero at the ends:  f'(0) = f'(L) = 0   (Neumann)

ODD extension   ->  period 2L, rotate about origin
   half-range SINE series:
   f(x) = sum_{n>=1} b_n sin(n pi x / L)
   b_n = (2/L) * integral_0^L f(x) sin(n pi x / L) dx
   value zero at the ends:   f(0) = f(L) = 0     (Dirichlet)
同一个在 [0, L] 上的 f,既可给出余弦级数(偶延拓,端点斜率为零),也可给出正弦级数(奇延拓,端点取值为零)——你挑那个与边界条件相符的。

为何边值问题替你定下选择

这个选择并非随意——物理替你挑定了它。当你用分离变量法求解杆的热方程或弦的波动方程时,方法交给你一族积木式的解,而两端的条件决定取哪一族。若两端被钉在零(温度被夹持为零,即狄利克雷条件 f(0) = f(L) = 0),就只有正弦合适,因为 sin(n pi x / L) 本来就在两端为零。若两端被绝热、没有热量流出(诺伊曼条件 f'(0) = f'(L) = 0),就只有余弦合适,因为 cos(n pi x / L) 在两端的斜率为零。

因此半区间展开恰恰是数据与积木之间的那座桥。要启动两端被冰镇的受热杆的冷却过程,你把初始温度分布展成半区间*正弦*级数,因为这样每个正弦模态都各自衰减、并在任何时刻都尊重被夹持的端点。倘若杆改为绝热,你就会把同一个分布展成半区间*余弦*级数。交到你手上的函数从未告诉你该用正弦还是余弦;是边界条件告诉你的。

一句诚实的告诫为这个循环收口。该级数只在 [0, L] 上等于你的 f;在外面,它等于你所发明的那个被反射、被重复的延拓,而非任何物理实在。而延拓的对称性可能引入一个角点:当 f(0) 不为 0 时,一个奇延拓会在原点造出一个跳跃,那里[[dirichlet-conditions|狄利克雷条件]]依然适用——级数收敛到跳跃的中点,并且由你上一篇见过的吉布斯现象,无论你保留多少项,都会在它附近过冲。半区间技巧确实强大,但它构造的是一个特定的周期世界,而你只对被问及的那一半负责。

一个具体画面:[0, L] 上的直线 f(x) = x

用一个尽可能朴素的分布把它落到实处:[0, L] 上的直线斜坡 f(x) = x。它的奇延拓是那条爬升、跌落、再爬升的锯齿波——一个在 L 的每个奇数倍处都有大小为 2L 的跳跃的函数,因为斜坡在 x = L 处收高,而下一份拷贝却起低。它的偶延拓则是那条爬升又下降的三角波,没有任何跳跃,只有尖锐的角点。同一个分布,两种截然不同的周期性格,全凭我们如何折叠那缺失的左半来选定。

  1. 正弦选择(奇,锯齿)。用一次分部积分算 b_n = (2/L) 乘以 x sin(n pi x / L) 从 0 到 L 的积分;结果是 b_n = (2L/(n pi)) 乘以 (-1)^{n+1}。系数只按 1/n 衰减——很慢——因为锯齿波带有真正的跳跃,而一个跳跃总会逼出一条 1/n 的尾巴。
  2. 余弦选择(偶,三角)。这里 a_0/2 = L/2(平均值),而 a_n = (2/L) 乘以 x cos(n pi x / L) 从 0 到 L 的积分,算得:奇数 n 时为 -4L/(n pi)^2,偶数 n 时为 0。这些系数按 1/n^2 衰减——快得多——因为三角波是连续的,只有角点;没有跳跃,也就没有那条慢吞吞的 1/n 尾巴。
  3. 读出寓意。两条级数都在 [0, L] 上重现 x,但余弦级数收敛得快得多、且没有吉布斯过冲,而正弦级数则爬得慢、并在端点附近过冲。系数衰减得越快,正是延拓越光滑的指纹——一个第一卷的念头重新现身:光滑度掌控着高频成分死去的快慢,这里是逐项地度量的。

这一对比就是整章的缩影。复指数形式是一个整洁的打包;奇偶捷径是工作量的免费减半;半区间展开则是为了在边值问题自己的语境里迎合它,而刻意制造出的对称性。而系数如何衰减,则悄悄回报出你所构造的延拓究竟有多粗糙。把这份直觉带向前去——当正弦余弦让位给后续各级上别的本征函数时,同样的正交性与对称性反射,会以推广的形式回归。