数学总导航 · 本课程目录 · ← 上一章 · 下一章 →
为什么调均衡器能压掉录音里的嗡嗡声?#
一段录音可能同时包含人声、低频嗡声和高频嘶声。直接在时间波形上看,它们叠在一起,难以分开。频率分析相当于问:这段声音与每种规则振动有多像?
先看易算的模型:f(t)=cos(2π⋅2t)+0.3cos(2π⋅5t),观察区间为一秒。它由振幅为 1 的 2 Hz 波和振幅为 0.3 的 5 Hz 波组成。2 Hz、5 Hz 在这里只为便于手算,不是可听人声模型。
若分别与这些余弦相乘再积分,交叉项相互抵消,同频项留下。于是你能读出两种成分。真实均衡器的处理还涉及滤波设计与时域影响,但“先识别频率方向,再调整对应系数”正是核心入口。
把两种频率真的分开算#
对一秒内的 f(t)=cos(4πt)+0.3cos(10πt),先与 cos(4πt) 相乘再积分。用积化和差公式,交叉项为
∫01cos(4πt)cos(10πt)dt=21∫01[cos(6πt)+cos(14πt)]dt=0.同频项 ∫01cos2(4πt)dt=1/2。所以读出 2 Hz 的系数时,需要再除以该方向的范数平方 1/2,得到 1。对 5 Hz 同理得到 0.3。
这就是正交投影公式,而不是仅靠图上找峰。观察窗口若不恰好包含整周期,交叉项一般不再恰好为零,实际分析就会出现频谱泄漏,需要考虑窗口选择。
为什么绕圆能检测频率#
把实信号乘以 e−2πiξt,是在对每个时刻施加一个反向旋转。如果信号中含有同样频率的复指数 e2πiξt,两者相乘后恰好不再旋转,累积时会留下稳定贡献。
频率不匹配的成分仍在旋转,正负或不同方向的贡献会部分抵消。3Blue1Brown 的绕圆图像展示的正是这个机制;这里保留积分与归一化,避免把“重心平均”与未归一化 Fourier 积分混在一起。
重构为什么还需要完备性#
正交只能保证每次加入的新方向不重复,不能保证你已经拥有所有方向。如果只保留偶数频率,某些奇数频率信号的所有读数都为零,但信号本身并非零。
完备正交系统要求没有被整个系统漏掉的非零方向。满足这一点后,全部系数才足以在 L2 意义下恢复原信号。Parseval 等式是在核对能量账:所有坐标能量相加,恰好等于原信号能量。
在跳跃边缘,误差小与看起来顺眼可以不同#
逼近一段突然开关的方波时,有限 Fourier 和仍是平滑波形,不可能在所有点瞬间变成跳跃。增加频率后,过渡区变窄,整体平方误差可以趋零,但边缘附近会有振铃。
因此重构质量要问清楚标准。如果只关心总能量误差,L2 很合适;若产品不允许任何过冲,就得加入更强要求。数学定理的收敛方式直接影响你能从“逼近成功”推出什么实际结论。
带你做:把过滤掉的频率误差算清楚#
对一秒信号 f=cos(4πt)+0.3cos(10πt),如果只保留 2 Hz 成分,误差是 0.3cos(10πt),平方误差积分为 0.09/2=0.045。
原信号总平方能量为 1/2+0.045=0.545,因为交叉项为零。因此保留主频留下约 91.7% 能量。这个百分比由明确的积分算出,不是看波形“好像差不多”。
最大时刻误差却可以达到 0.3。整体能量损失小与峰值偏差小是不同指标;如果被滤掉的是你真正需要的微弱报警音,能量小也不能说明它不重要。
再把观察窗口缩短,原先的正交积分未必仍为零。实际处理前应检查窗口、采样率和频率定义,而不是把一秒整周期的理想计算无条件用于任何录音片段。
回头整理:定义、工具与练习#
从调色盘理解频率#
声音里混着低音与高音。要问某个频率含量,不是看某一时刻有多高,而是让信号与该频率的波形对齐,累积它们的相似程度。这个累积就是内积。
3Blue1Brown 的傅里叶变换课程 用“把波形绕圆并观察重心”的图像解释频率检测;它对应复指数乘法后的积分。本章把这一直觉连接到 Hilbert 投影,并分别处理周期级数与实线上的变换。
正交坐标可无限延伸#
设 (en) 为标准正交列,有限部分的最佳近似是
PNx=n=1∑N⟨x,en⟩en.由勾股定理,
∥x−PNx∥2=∥x∥2−n=1∑N∣⟨x,en⟩∣2≥0.得到 Bessel 不等式。系数平方可和使投影序列为 Cauchy 列;Hilbert 完备性保证极限存在。
若这些方向的闭线性张成是整个 H,称它们为完备正交系统,此时余量只能为零,得到
x=n∑⟨x,en⟩en,∥x∥2=n∑∣⟨x,en⟩∣2.第二式为 Parseval 等式。可分 Hilbert 空间存在可数标准正交基;非可分情形需要更一般的指标集。
周期 Fourier 级数#
在复 L2[0,2π] 上取归一化内积 ⟨f,g⟩=(2π)−1∫fg。则 ek(t)=eikt(k∈Z)为完备正交系统。完备性是一个实质定理,可由三角多项式的稠密性证明。
f(k)=2π1∫02πf(t)e−iktdt,f−k=−N∑Nf(k)eikt2→0.对 f(t)=2+3cost,非零系数是 f(0)=2、f(±1)=3/2,归一化能量为 4+9/2=17/2。
平方误差趋零没有保证什么#
L2 收敛保证整体平方误差趋零,不自动保证每一点收敛,更不自动保证一致收敛。对分段光滑周期函数,可用额外的 Fourier 定理说明在跳跃点收敛到左右极限均值。跳跃附近还会出现 Gibbs 现象。
这些逐点定理有各自条件,不能仅从 Parseval 等式推导。
实线上的 Fourier 变换#
对 f∈L1(R),可定义 f(ξ)=∫f(t)e−2πiξtdt。实线上的 e2πiξt 不属于 L2(R),所以不能直接把它们称为该空间内的普通正交基。
Plancherel 定理把变换从 L1∩L2 延拓为 L2 上的酉映射,并保持能量。一般 L2 函数的变换由范数极限定义,原来的逐点积分未必绝对收敛。这里引用该定理,不把周期级数证明误用到非周期情形。
练习与解答#
- 标准正交列是否一定完备?
- f=sint 的归一化平方范数是多少?
- 最佳 N 维正交近似误差如何计算?
解答:1. 不一定,只取 ℓ2 的偶数编号单位向量就不完备。2. 1/2。3. ∥x∥2−∑n=1N∣⟨x,en⟩∣2;只要这些方向正交,此式成立,不必先假设整个系统完备。
数学总导航 · 本课程目录 · ← 上一章 · 下一章 →