Fourier 级数 (Fourier Series)
Fourier 级数(傅里叶级数)是分析学中处理周期现象的核心工具。它的基本思想是:任何符合一定条件的周期函数,都可以表示为一系列正弦和余弦函数的叠加。这一理论不仅在纯数学(如偏微分方程、数论)中至关重要,也是现代信号处理、量子力学及通信工程的数学基石。
三角函数系 {1,cosx,sinx,cos2x,sin2x,…,cosnx,sinnx,…} 在区间 [−π,π] 上具有 正交性:
- ∫−ππcosnxsinmxdx=0 (对所有 n,m)
- ∫−ππcosnxcosmxdx=⎩⎨⎧0,π,2π,n=mn=m=0n=m=0
- ∫−ππsinnxsinmxdx={0,π,n=mn=m=0
设 f(x) 是以 2π 为周期的可积函数,其 Fourier 级数为:
f(x)∼2a0+∑n=1∞(ancosnx+bnsinnx)
其中系数定义为:
an=π1∫−ππf(x)cosnxdx,n=0,1,2,…
bn=π1∫−ππf(x)sinnxdx,n=1,2,…
Fourier 级数是否收敛于 f(x) 是一个深刻的问题。Dirichlet 给出了最实用的判定准则。
若 f(x) 是周期为 2π 的函数,且在一个周期内满足以下 Dirichlet 条件:
- 分段连续:仅有有限个第一类间断点。
- 分段单调:仅有有限个极值点。
则 f(x) 的 Fourier 级数 S(x) 在每一点 x 处均收敛,且:
S(x)=2f(x+0)+f(x−0)
- 在 连续点 x 处,S(x)=f(x);
- 在 间断点 x 处,S(x) 收敛于左右极限的平均值。
- 物理直觉:在突变点,级数倾向于取“中间值”,这体现了三角级数作为光滑函数逼近非光滑函数时的平滑特性。
- 一致收敛性:若 f(x) 全线连续且分段光滑,则其余 Fourier 级数在整个实轴上 一致收敛 于 f(x)。
对于定义在 [0,π] 上的非周期函数 f(x),我们可以通过 延拓 将其展开。
定义 F(x)={f(x),−f(−x),0≤x≤π−π≤x<0。
此时 an=0,级数仅包含正弦项:
f(x)∼∑n=1∞bnsinnx,bn=π2∫0πf(x)sinnxdx
定义 F(x)={f(x),f(−x),0≤x≤π−π≤x<0。
此时 bn=0,级数仅包含常数项和余弦项:
f(x)∼2a0+∑n=1∞ancosnx,an=π2∫0πf(x)cosnxdx
若函数定义在 [a,b] 上,通常先平移至 [0,L],再根据需要进行周期延拓。
- 普通周期延拓:直接令 f(x+L)=f(x),通常在端点处会引入间断。
- 反射延拓:如上述奇/偶延拓,通常能保持端点的连续性,从而提高收敛速度。
当 f(x) 存在跳变间断点时,Fourier 级数的部分和 SN(x) 在间断点附近会产生“振荡”和“过冲”。
设 f(x) 在 x0 处有一个跳变间断点,跳变高度为 h=∣f(x0+0)−f(x0−0)∣。其第 N 阶部分和 SN(x) 在 x0 附近的第一个极大值与函数值之差趋于:
limN→∞SN(xmax)=f(x0+0)+2h[π2∫0πtsintdt−1]
其中常量 G=π2Si(π)≈1.17897…。
- 过冲百分比:(G−1)/2≈0.08949⋯≈9%。
- 本质分析:Gibbs 现象揭示了用光滑正交基逼近非光滑信号时的能量分布特性。即使增加谐波项数,过冲幅度也不会消失,仅是过冲点向间断点无限靠拢。这体现了 L2 空间均方收敛并不保证点态一致收敛。
在图像处理中,Gibbs 现象表现为“振铃效应”(Ringing artifacts)。在 JPEG 压缩等频域算法中,必须通过窗函数 (Windowing) 或低通滤波来抑制这种由于频谱截断引起的伪影。
从 Hilbert 空间的视角看,Fourier 展开是正交基下的投影。
若 f(x) 在 [−π,π] 上平方可积,则其平均功率等于各谐波分量功率之和:
π1∫−ππ∣f(x)∣2dx=2a02+∑n=1∞(an2+bn2)
- 求级数和:它是计算如 ∑n21,∑n41 等数项级数的强力工具。
- 误差分析:用于评估用前 N 项逼近原函数时的均方误差。
设 f(x)=x,x∈[0,π]。
- 将其进行 偶延拓 展开为余弦级数。
- 利用结果求 ∑n=1∞(2n−1)21。
解析:
-
偶延拓:在 [−π,π] 上 F(x)=∣x∣。
a0=π2∫0πxdx=π。
an=π2∫0πxcosnxdx=π2[nxsinnx0π−∫0πnsinnxdx]=πn22(cosnπ−1)=πn22((−1)n−1)。
- n 为偶数时,an=0;
- n 为奇数时,an=−πn24。
展开式:x=2π−π4∑k=1∞(2k−1)2cos(2k−1)x,x∈[0,π]。
-
求和:令 x=0,由于偶延拓后在 0 处连续,f(0)=0。
0=2π−π4∑k=1∞(2k−1)21⟹∑k=1∞(2k−1)21=8π2。
已知 f(x)=x2,x∈[−π,π],其展开式为 f(x)=3π2+4∑n=1∞n2(−1)ncosnx。证明:
∑n=1∞n41=90π4
解析:
应用 Parseval 等式:π1∫−ππ(x2)2dx=2a02+∑n=1∞an2。
左边:π1∫−ππx4dx=5π2π5=52π4。
右边:a0=32π2,an=n24(−1)n。
52π4=21(32π2)2+∑n=1∞n416=92π4+16∑n=1∞n41。
整理得:16∑n41=4518π4−10π4=458π4⟹∑n41=90π4。
设 f(x)=x, x∈(−π,π) 且 f(x+2π)=f(x)。
- 求其 Fourier 展开式。
- 利用结果计算 ∑n=1∞nsinn。
解析:
- f(x) 是奇函数,故 an=0。
bn=π2∫0πxsinnxdx=π2[−nxcosnx+n2sinnx]0π=n2(−1)n+1。
展开式为:x∼2∑n=1∞n(−1)n+1sinnx,x∈(−π,π)。
- 令 x=1。由于 x=1 是 f(x) 的连续点,由 Dirichlet 定理:
1=2∑n=1∞n(−1)n+1sinn⟹∑n=1∞n(−1)n+1sinn=21。
注意:若要求 ∑nsinn,通常考虑 f(x)=2π−x 的展开。
设 f(x)={1,0,0<x<π−π<x<0。求其 Fourier 级数,并讨论 x=0 处的收敛值。
解析:
a0=π1∫0π1dx=1。
an=π1∫0πcosnxdx=0 (n≥1)。
bn=π1∫0πsinnxdx=πn1(1−cosnπ)=πn1−(−1)n。
- n 为偶数时,bn=0;
- n 为奇数时,bn=πn2。
展开式:f(x)∼21+π2∑k=1∞2k−1sin(2k−1)x。
收敛性:在 x=0 处,f(x) 有第一类间断点。
收敛值 S(0)=2f(0+0)+f(0−0)=21+0=21。
利用 f(x)=x(π−x) 在 [0,π] 上的正弦级数展开,计算 ∑n=1∞(2n−1)61。
解析:
- 奇延拓后 bn=π2∫0π(xπ−x2)sinnxdx=πn38 (n 为奇数),bn=0 (n 为偶数)。
- 应用 Parseval 等式:π2∫0π∣x(π−x)∣2dx=∑bn2。
左边:π2∫0π(x2π2−2x3π+x4)dx=π2[3π5−2π5+5π5]=15π4。
右边:∑k=1∞(π(2k−1)38)2=π264∑k=1∞(2k−1)61。
整理得:∑k=1∞(2k−1)61=960π6。
- 基础:将 f(x)=sin2x 在 [−π,π] 上展开为 Fourier 级数。
- 周期延拓:将 f(x)=ex 在 [0,1] 上分别进行奇延拓和偶延拓,写出其正弦级数和余弦级数的系数表达式。
- 收敛性分析:讨论函数 f(x)=∑n=1∞n1.5sinnx 的连续性与可微性。
- Parseval 实战:利用 f(x)=∣x∣,x∈[−π,π] 的展开式,通过 Parseval 等式求 ∑n=1∞(2n−1)41。
- 挑战 (Riemann-Lebesgue):证明若 f(x) 在 [a,b] 上 Riemann 可积,则 limn→∞∫abf(x)sinnxdx=0。并解释该引理在 Fourier 级数系数趋于零中的作用。
Fourier 级数不仅仅是一个数学公式,它代表了从时间域到频率域的视角转换。掌握它的关键在于深刻理解 Dirichlet 条件下的点收敛特性,以及 Parseval 等式所揭示的能量守恒本质。
练习库同步 (Analysis Exercise Sync)