1764 字
9 分钟
泛函分析 08|正交展开、Fourier 与 Parseval

数学总导航 · 本课程目录 · ← 上一章 · 下一章 →

为什么调均衡器能压掉录音里的嗡嗡声?#

一段录音可能同时包含人声、低频嗡声和高频嘶声。直接在时间波形上看,它们叠在一起,难以分开。频率分析相当于问:这段声音与每种规则振动有多像?

先看易算的模型:f(t)=cos⁡(2π⋅2t)+0.3cos⁡(2π⋅5t)f(t)=\cos(2\pi\cdot2t)+0.3\cos(2\pi\cdot5t),观察区间为一秒。它由振幅为 1 的 2 Hz 波和振幅为 0.3 的 5 Hz 波组成。2 Hz、5 Hz 在这里只为便于手算,不是可听人声模型。

若分别与这些余弦相乘再积分,交叉项相互抵消,同频项留下。于是你能读出两种成分。真实均衡器的处理还涉及滤波设计与时域影响,但“先识别频率方向,再调整对应系数”正是核心入口。

把两种频率真的分开算#

对一秒内的 f(t)=cos⁡(4πt)+0.3cos⁡(10πt)f(t)=\cos(4\pi t)+0.3\cos(10\pi t),先与 cos⁡(4πt)\cos(4\pi t) 相乘再积分。用积化和差公式,交叉项为

∫01cos⁡(4πt)cos⁡(10πt)dt=12∫01[cos⁡(6πt)+cos⁡(14πt)]dt=0.\int_0^1\cos(4\pi t)\cos(10\pi t)dt =\frac12\int_0^1[\cos(6\pi t)+\cos(14\pi t)]dt=0.

同频项 ∫01cos⁡2(4πt)dt=1/2\int_0^1\cos^2(4\pi t)dt=1/2。所以读出 2 Hz 的系数时,需要再除以该方向的范数平方 1/21/2,得到 1。对 5 Hz 同理得到 0.3。

这就是正交投影公式,而不是仅靠图上找峰。观察窗口若不恰好包含整周期,交叉项一般不再恰好为零,实际分析就会出现频谱泄漏,需要考虑窗口选择。

为什么绕圆能检测频率#

把实信号乘以 e−2πiξte^{-2\pi i\xi t},是在对每个时刻施加一个反向旋转。如果信号中含有同样频率的复指数 e2πiξte^{2\pi i\xi t},两者相乘后恰好不再旋转,累积时会留下稳定贡献。

频率不匹配的成分仍在旋转,正负或不同方向的贡献会部分抵消。3Blue1Brown 的绕圆图像展示的正是这个机制;这里保留积分与归一化,避免把“重心平均”与未归一化 Fourier 积分混在一起。

重构为什么还需要完备性#

正交只能保证每次加入的新方向不重复,不能保证你已经拥有所有方向。如果只保留偶数频率,某些奇数频率信号的所有读数都为零,但信号本身并非零。

完备正交系统要求没有被整个系统漏掉的非零方向。满足这一点后,全部系数才足以在 L2L^2 意义下恢复原信号。Parseval 等式是在核对能量账:所有坐标能量相加,恰好等于原信号能量。

在跳跃边缘,误差小与看起来顺眼可以不同#

逼近一段突然开关的方波时,有限 Fourier 和仍是平滑波形,不可能在所有点瞬间变成跳跃。增加频率后,过渡区变窄,整体平方误差可以趋零,但边缘附近会有振铃。

因此重构质量要问清楚标准。如果只关心总能量误差,L2L^2 很合适;若产品不允许任何过冲,就得加入更强要求。数学定理的收敛方式直接影响你能从“逼近成功”推出什么实际结论。

带你做:把过滤掉的频率误差算清楚#

对一秒信号 f=cos⁡(4πt)+0.3cos⁡(10πt)f=\cos(4\pi t)+0.3\cos(10\pi t),如果只保留 2 Hz 成分,误差是 0.3cos⁡(10πt)0.3\cos(10\pi t),平方误差积分为 0.09/2=0.0450.09/2=0.045。

原信号总平方能量为 1/2+0.045=0.5451/2+0.045=0.545,因为交叉项为零。因此保留主频留下约 91.7%91.7\% 能量。这个百分比由明确的积分算出,不是看波形“好像差不多”。

最大时刻误差却可以达到 0.3。整体能量损失小与峰值偏差小是不同指标;如果被滤掉的是你真正需要的微弱报警音,能量小也不能说明它不重要。

再把观察窗口缩短,原先的正交积分未必仍为零。实际处理前应检查窗口、采样率和频率定义,而不是把一秒整周期的理想计算无条件用于任何录音片段。

回头整理:定义、工具与练习#

从调色盘理解频率#

声音里混着低音与高音。要问某个频率含量,不是看某一时刻有多高,而是让信号与该频率的波形对齐,累积它们的相似程度。这个累积就是内积。

3Blue1Brown 的傅里叶变换课程 用“把波形绕圆并观察重心”的图像解释频率检测;它对应复指数乘法后的积分。本章把这一直觉连接到 Hilbert 投影,并分别处理周期级数与实线上的变换。

正交坐标可无限延伸#

设 (en)(e_n) 为标准正交列,有限部分的最佳近似是

PNx=∑n=1N⟨x,en⟩en.P_Nx=\sum_{n=1}^N\langle x,e_n\rangle e_n.

由勾股定理,

∥x−PNx∥2=∥x∥2−∑n=1N∣⟨x,en⟩∣2≥0.\|x-P_Nx\|^2=\|x\|^2-\sum_{n=1}^N|\langle x,e_n\rangle|^2\ge0.

得到 Bessel 不等式。系数平方可和使投影序列为 Cauchy 列;Hilbert 完备性保证极限存在。

若这些方向的闭线性张成是整个 HH,称它们为完备正交系统,此时余量只能为零,得到

x=∑n⟨x,en⟩en,∥x∥2=∑n∣⟨x,en⟩∣2.x=\sum_n\langle x,e_n\rangle e_n,\qquad \|x\|^2=\sum_n|\langle x,e_n\rangle|^2.

第二式为 Parseval 等式。可分 Hilbert 空间存在可数标准正交基;非可分情形需要更一般的指标集。

周期 Fourier 级数#

在复 L2[0,2π]L^2[0,2\pi] 上取归一化内积 ⟨f,g⟩=(2π)−1∫fg‾\langle f,g\rangle=(2\pi)^{-1}\int f\overline g。则 ek(t)=eikte_k(t)=e^{ikt}(k∈Zk\in\mathbb Z)为完备正交系统。完备性是一个实质定理,可由三角多项式的稠密性证明。

f^(k)=12π∫02πf(t)e−iktdt,∥f−∑k=−NNf^(k)eikt∥2→0.\widehat f(k)=\frac1{2\pi}\int_0^{2\pi}f(t)e^{-ikt}dt, \quad \left\|f-\sum_{k=-N}^N\widehat f(k)e^{ikt}\right\|_2\to0.

对 f(t)=2+3cos⁡tf(t)=2+3\cos t,非零系数是 f^(0)=2\widehat f(0)=2、f^(±1)=3/2\widehat f(\pm1)=3/2,归一化能量为 4+9/2=17/24+9/2=17/2。

平方误差趋零没有保证什么#

L2L^2 收敛保证整体平方误差趋零,不自动保证每一点收敛,更不自动保证一致收敛。对分段光滑周期函数,可用额外的 Fourier 定理说明在跳跃点收敛到左右极限均值。跳跃附近还会出现 Gibbs 现象。

这些逐点定理有各自条件,不能仅从 Parseval 等式推导。

实线上的 Fourier 变换#

对 f∈L1(R)f\in L^1(\mathbb R),可定义 f^(ξ)=∫f(t)e−2πiξtdt\widehat f(\xi)=\int f(t)e^{-2\pi i\xi t}dt。实线上的 e2πiξte^{2\pi i\xi t} 不属于 L2(R)L^2(\mathbb R),所以不能直接把它们称为该空间内的普通正交基。

Plancherel 定理把变换从 L1∩L2L^1\cap L^2 延拓为 L2L^2 上的酉映射,并保持能量。一般 L2L^2 函数的变换由范数极限定义,原来的逐点积分未必绝对收敛。这里引用该定理,不把周期级数证明误用到非周期情形。

练习与解答#

  1. 标准正交列是否一定完备?
  2. f=sin⁡tf=\sin t 的归一化平方范数是多少?
  3. 最佳 NN 维正交近似误差如何计算?

解答:1. 不一定,只取 ℓ2\ell^2 的偶数编号单位向量就不完备。2. 1/21/2。3. ∥x∥2−∑n=1N∣⟨x,en⟩∣2\|x\|^2-\sum_{n=1}^N|\langle x,e_n\rangle|^2;只要这些方向正交,此式成立,不必先假设整个系统完备。


数学总导航 · 本课程目录 · ← 上一章 · 下一章 →

泛函分析 08|正交展开、Fourier 与 Parseval
https://f1yingwhite.github.io/posts/数学/泛函分析/08-正交展开、Fourier-与-Parseval/
作者
F1yingWhite
发布于
2026-11-03
许可协议
CC BY-NC-SA 4.0