分段周期函數用單一周期函數逼近的數學方法：傅里葉級數

問題背景

給定一個連續的分段周期函數，尋找一個單一的解析周期函數對其進行逼近。

核心理論：傅里葉級數

根據傅里葉理論，任何滿足狄利克雷條件的周期函數（包括分段連續函數）都可以展開為正弦/餘弦級數之和：

$$f(x) \sim \frac{a_0}{2} + \sum_{n=1}^{\infty} \left[ a_n \cos(n\omega x) + b_n \sin(n\omega x) \right]$$

系統系數計算

對於周期為 $T = 2\pi/\omega$ 的函數：

• $a_0 = \frac{\omega}{\pi} \int_{0}^{T} f(x) dx$
• $a_n = \frac{\omega}{\pi} \int_{0}^{T} f(x) \cos(n\omega x) dx$
• $b_n = \frac{\omega}{\pi} \int_{0}^{T} f(x) \sin(n\omega x) dx$

分段函數在各連續區間上分別積分再求和即可。

工程應用要點

1. 截斷階數：實際應用取前 $N$ 項，$N$ 越大逼近越精確，但計算量增加。
2. 吉布斯現象：在不連續點附近會有超調振盪，無法通過增加項數消除，需用濾波器平滑。
3. 頻域分析：傅里葉係數直接反映各諧波分量能量分布，可用於信號處理、壓縮。

替代逼近方法

學習資源建議

• 高等數學教材「傅里葉級數」章節。
• 《信號與系統》(Oppenheim) 第 3 章。
• MATLAB/Python 實驗：fft、scipy.signal.periodogram 驗證係數。