函数作为向量
傅里叶级数不仅仅是计算工具——它是一个几何操作:函数在基上的正交投影。
函数空间
在 [0,1] 上平方可积函数的集合形成向量空间 L²[0,1]。加法和标量乘法按点进行。两个函数 f、g 的内积:
⟨f, g⟩ = ∫₀¹ f(t) · g*(t) dt
其中 g* 是 g 的复共轭。这满足所有内积公理。
傅里叶基的正交性
函数 φ_k(t) = e^{i2πkt} 形成 L²[0,1] 的正交归一基:
⟨φ_k, φ_m⟩ = ∫₀¹ e^{i2πkt} · e^{−i2πmt} dt = ∫₀¹ e^{i2π(k−m)t} dt = δ_{km}
(如果 k = m 则等于 1,否则等于 0——通过在完整周期上积分纯振荡。)
傅里叶系数作为内积
x(t) 的第 k 个傅里叶系数:
c_k = ⟨x, φ_k⟩ = ∫₀¹ x(t) · e^{−i2πkt} dt
这是 x 在基向量 φ_k 上的投影。系数衡量 φ_k 在 x 中的分量大小。
投影到子空间
傅里叶级数截断为 2N+1 项将 x 投影到由 {φ_{−N}、…、φ_N} 张成的子空间。截断级数是 x 在该有限维子空间上的正交投影。
根据贝塞尔不等式,投影最小化 L² 误差:
‖x − Σ_{k=−N}^{N} c_k φ_k‖² ≤ ‖x − Σ_{k=−N}^{N} a_k φ_k‖² 对于 a_k 的任何选择
傅里叶截断是来自该子空间的 L² 最佳近似。它最小化均方误差(差的 L² 范数的平方)。
矩形窗口 → Sinc 核
时域中的矩形窗口(仅保留 |k| ≤ N 的系数)对应系数索引中的 rect 函数的乘法。
一个域中的乘法对应于另一个域中的卷积。
矩形窗口的傅里叶变换(在离散系数空间中)是狄利克雷核——一个周期性的 sinc 类似函数:
D_N(f) = Σ_{k=−N}^{N} e^{i2πfk} = sin(π(2N+1)f) / sin(πf)
当我们截断傅里叶级数时,我们用 D_N(f) 卷积理想响应 H_ideal(f)。
为什么会出现吉布斯现象
狄利克雷核具有衰减缓慢的大旁瓣。在 H_ideal(f) 的阶跃不连续处附近,这些旁瓣产生振铃——它们在跳跃的一侧相干叠加,产生约 9% 的超调。
数学常数:∫₀^π sin(t)/t dt = Si(π) ≈ 1.8519。吉布斯超调高度 = (2/π)·Si(π) − 1 ≈ 0.0895 = 8.95%。这与 N 无关。
窗口几何
平滑窗口(汉明、汉宁、凯撒)的傅里叶变换具有更低的旁瓣。用具有较小旁瓣的核卷积 H_ideal(f) 会产生较少的振铃。权衡:较低的旁瓣总是伴随着较宽的主瓣,加宽过渡带。
吉布斯常数
吉布斯超调是一个定积分,而不是 N 的函数。
单位阶跃的 N 项部分傅里叶和的第一个最大值出现在距离不连续处约 f ≈ 1/(2N) 处。当 N → ∞ 时,这个最大值接近 1/(2)·(2/π)·Si(π) ≈ 1.0895。
超调:跳跃高度的 0.0895 或约 8.95%。
作为频域核的窗口
每个窗口函数都有一个傅里叶变换,描述用于平滑理想频率响应的核。
核的关键几何参数:
1. 主瓣宽度:确定过渡带宽度(更宽的主瓣 → 更宽的过渡)。
2. 峰值旁瓣水平:确定通带和阻带纹波(更低的旁瓣 → 更少的纹波)。
这两个参数不是独立的。对于给定的窗口长度 2N+1,降低旁瓣高度需要加宽主瓣——总是如此。
凯撒窗口给用户一个旋钮(α)来连续权衡旁瓣高度与主瓣宽度,而不是在固定窗口类型之间跳跃。
设计见解
过渡带宽度 ΔF ≈ 主瓣宽度 / N。纹波 δ ≈ 旁瓣水平。两个公式都是近似的;凯撒的方程使它们精确。