傅里葉級數展開


閱讀本篇內容之前可先閱讀博客:三角函數定義歐拉公式

拉格朗日等數學家發現某些周期函數可以由三角函數的和來表示,比如下圖中,黑色的斜線是周期為 $2\pi$ 的函數,而紅色的曲線是三角函數

之和,可以看出兩者確實近似:

             

另一位數學家傅里葉猜測任意周期函數都可以寫成三角函數之和。

首先先證明一個結論:任意一個函數都可拆分成一個偶函數和一個奇函數的和。

證明

   易知

$$f(x) = \frac{f(x) - f(-x)}{2} + \frac{f(x) + f(-x)}{2} = g(x) + h(x)$$

   因為

$$g(-x) = \frac{f(-x) - f(x)}{2} = -g(x) \\
h(-x) = \frac{f(-x) + f(x)}{2} = h(x)$$

   所以 $g(x)$ 是奇函數,$h(x)$ 是偶函數,即 $f(x)$ 可以展開為奇函數和偶函數的和。

證畢

傅里葉就試圖將周期為 $T$ 的函數 $f(x)$ 展開為 $\sin x$ 和 $\cos x$ 函數和的形式。

那怎么保證組合出來的函數周期依然為 $T$ 呢?

我們知道函數 $f = \sin \omega t$ 的周期為 $T^{'} = \frac{2\pi}{\omega}$,要使得原函數能夠被三角函數表示,那么三角函數的粒度必然要小於原函數,

即三角函數的最小周期 $T^{'}$ 必須小於原函數 $f(x)$ 的最小周期 $T$,即

$$\frac{2\pi}{\omega} \leq T \\
\Rightarrow \frac{2\pi n}{\omega} = T, \; n > 0$$

這里的 $n$ 必須是一個整數,對於周期函數,它的分量的角頻率也必定是原函數角頻率的整數倍,這一點不做證明。

現在我們有了各種三角函數,如 $\sin x,\sin 2x,\sin 3x,\sin 4x$,對於最初舉的周期為 $2\pi$ 的函數,可以用這些三角函數慢慢逼近。

                

綜上,構造出來的三角函數展開式形式為

$$f(x) = \sum_{n = 1}^{\infty}\left ( a_{n}\cos \frac{2\pi n}{T}x + b_{n}\sin \frac{2\pi n}{T}x \right ) + C$$

其中 $\frac{2\pi}{T}$ 是原始信號的角頻率,因為 $n > 1$,可見分量的角頻率必然不小於原始信號的頻率,是原始信號頻率的整數倍。所以這里的 $n$ 

不要認為是索引序號,而是分量的角頻率與原始信號角頻率的倍數關系,如果沒有某個倍數對應的分量,那么這一項就是 $0$。

那怎么確定參數 $C,a_{n},b_{n}$ 呢?先來了解一下三角函數的正交性。

三角函數系列 $1,\cos x, \sin x, \cos 2x, \sin 2x,..., \cos nx, \sin nx,...$,在區間 $\left [ -\pi, +\pi \right ]$ 上是正交函數系,即

$$\int_{-\pi}^{\pi}\cos nxdx = 0, (n = 1,2,3,...)  \;\;\; \int_{-\pi}^{\pi}\sin nxdx = 0, (n = 1,2,3,...)   \\
\int_{-\pi}^{\pi}\cos mx \sin nx dx = 0, (m,n = 1,2,3,...)  \\
\int_{-\pi}^{\pi}\cos mx \cos nx dx = 0, (m,n = 1,2,3,..., m \neq n) \;\;\; \int_{-\pi}^{\pi}\sin mx \sin nx dx = 0, (m,n = 1,2,3,..., m \neq n) $$

比如,我們要求 $\cos \frac{2\pi n}{T}x$ 這一項的系數,兩邊都乘上它,然后兩邊在 $0-T$ 積分,注意這里的 $T$ 是原始信號的周期,但因為 $n > 1$,所以分解的

分量信號的周期都會比 $T$ 小而且是整數倍,那么依然可以用正交性,就會發現只留下自身那一項,基於此可求得

                                                     

求 $C$ 的話,兩邊直接在 $0-T$ 上積分即可,可求得

$$C = \frac{a_{0}}{2}$$

在求系數的公式中,$T,f(x)$ 都是已知的,所以只要給定一個周期信號,我們就能通過傅里葉級數分解為一個一個的正弦或余弦分量,求系數很容易,

代公式算即可,難的是為什么一個周期信號一定可以分解為三角函數求和的形式,傅里葉級數展開的形式是怎么來的?這里就不做證明了。

根據歐拉公式,我們來推導一下傅里葉級數的指數形式:

$$f(x) = \frac{a_{0}}{2} + \sum_{n = 1}^{\infty}\left ( a_{n}\cos n\omega x + b_{n}\sin n\omega x \right ) \\
= \frac{a_{0}}{2} + \sum_{n = 1}^{\infty}\left [ a_{n}\frac{e^{in\omega x} + e^{-in\omega x}}{2} + b_{n}\frac{e^{in\omega x} - e^{-in\omega x}}{2i} \right ] \\
= \frac{a_{0}}{2} + \sum_{n = 1}^{\infty}\left [ \frac{a_{n} - ib_{n}}{2}e^{in\omega x} + \frac{a_{n} + ib_{n}}{2}e^{-in\omega x} \right ]$$

其中,$\omega$ 是原始信號的角頻率,$n\omega$ 是分量的角頻率,也可寫成 $\omega_{n}$,因為

$$\frac{a_{n} - ib_{n}}{2} = \frac{1}{T}\int_{0}^{T}f(x)\cos n\omega xdx - i\frac{1}{T}\int_{0}^{T}f(x)\sin n\omega xdx \\
= \frac{1}{T}\int_{0}^{T}f(x)\left (\cos n\omega x - i\sin n\omega x \right )dx \\
= \frac{1}{T}\int_{0}^{T}f(x)e^{-in\omega x}dx, \; n = 1,2,3,\cdots $$

$$\frac{a_{n} + ib_{n}}{2} = \frac{1}{T}\int_{0}^{T}f(x)\cos n\omega xdx + i\frac{1}{T}\int_{0}^{T}f(x)\sin n\omega xdx \\
= \frac{1}{T}\int_{0}^{T}f(x)\left (\cos n\omega x + i\sin n\omega x \right )dx \\
= \frac{1}{T}\int_{0}^{T}f(x)e^{in\omega x}dx, \; n = 1,2,3,\cdots$$

$$\frac{a_{0}}{2} = \frac{1}{T}\int_{0}^{T}f(x)dx$$

驚訝得發現

$$\frac{a_{n} - ib_{n}}{2} \; \xrightarrow[]{let \; n = -n} \; \frac{a_{n} + ib_{n}}{2} \\
\frac{a_{n} - ib_{n}}{2} \; \xrightarrow[]{let \; n = 0} \;  \frac{a_{0}}{2}$$

若令

$$c_{n} = \frac{1}{T}\int_{0}^{T}f(x)e^{-in\omega x}dx \\
= \frac{1}{T}\int_{0}^{T}f(x)e^{-i\omega_{n} x}dx, \; n =0, \pm 1, \pm 2, \pm 3, \cdots$$

其中 $\omega_{n}$ 是分量的角頻率,則傅里葉級數展開可以寫出如下形式:

$$f(x) = \sum_{n = -\infty}^{+ \infty} c_{n}e^{-in\omega x} = \sum_{n = -\infty}^{+ \infty} c_{n}e^{-i\omega_{n} x}, \;\; n \in Z$$

指數形式的傅里葉級數展開比三角形式的簡潔,前者只需要求一個參數 $c_{n}$,而后者需要求三個參數 $a_{n},b_{n},a_{0}$。

由歐拉公式我們知道,$e^{j\theta}$ 表示復平面上的一個旋轉了 $\theta$ 弧度的向量,那么指數形式的傅里葉形式就可以理解為:無數旋轉的疊加。

或者說是:圓周運動的組合。

               

因為歐拉公式的存在,正弦函數或者余弦函數以及它們的組合就可以看成虛數 $e^{j\theta}$ 的實部或者虛部,也就提供了另一個角度來觀察函數 $f(x)$

的展開分量,原來的級數展開是時域上的各個信號相加,到復平面之后,它的分量就變成了一個個向量 $e^{j\omega t}$,在以各自的角頻率旋轉。比如上面

右圖中所示就是函數 $g(x)$ 分解為兩個正弦分類后在復平面上的情況:

$$g(x) = \sin x + \sin 2x    \\
\Rightarrow e^{ix} = \cos x + i \sin x \;\; e^{i2x} = \cos 2x + i \sin 2x  \\
\Rightarrow e^{ix} + e^{i2x} = (\cos x + \cos 2x) + i(\sin x + \sin 2x)$$

那么 $g(x)$ 的時域圖像就是向量 $e^{ix}$ 和 $e^{i2x}$ 在復平面上合成后的虛軸變化情況。

頻譜一般分為幅度譜和相位譜兩種,那怎么畫幅度譜呢?幅度譜是哪兩個量之間的關系呢?幅度譜分為兩種

1)單邊幅度譜

   只考慮頻率為正數部分得到的頻譜,由三角形式的傅里葉級數得到,根據三角函數的歸一公式是可以只保留正弦函數或余弦函數的,即

$$f(x) = \frac{a_{0}}{2} + \sum_{n = 1}^{\infty}\left ( a_{n}\cos n\omega x + b_{n}\sin n\omega x \right ) \\
= c_{0} + \sum_{n = 1}^{\infty}c_{n}\cos\left ( n\omega x + \varphi_{n} \right)$$

   其中

$$\varphi_{n} = -\arctan \frac{b_{n}}{a_{n}} \\
c_{n} = \sqrt{a_{n}^{2} + b_{n}^{2}}, n = 1,2,3,\cdots \\
c_{0} = \frac{a_{0}}{2}$$

   單邊幅度譜就是 $c_{n},n = 1,2,3,...$ 關於 $n\omega$ 的離散的圖像。

2)雙邊幅度譜

   由指數形式的傅里葉級數推導過程可知:

$$c_{n} = \frac{1}{T}\int_{0}^{T}f(x)e^{-in\omega x}dx, \; n =0, \pm 1, \pm 2, \pm 3, \cdots$$

   因為 $c_{n}$ 是一個復數,可以可以寫成模和輻角的形式,即

$$c_{n} = \left | c_{n} \right |e^{i\varphi_{n}}, \; n =0, \pm 1, \pm 2, \pm 3, \cdots$$

   又因為

$$c_{n} = \frac{a_{n} - ib_{n}}{2}, \; n =0, \pm 1, \pm 2, \pm 3, \cdots$$

   所以

$$\left |c_{n}  \right | = \frac{\sqrt{a_{n}^{2} + b_{n}^{2}}}{2} \\
\varphi_{n} = -\arctan \frac{b_{n}}{a_{n}}$$

   雙邊幅度譜就是 $c_{n},n = 0, \pm1,\pm2,\pm3,...$ 關於 $n\omega$ 的離散的圖像。

   可以看出,除了 $n = 0$,雙邊幅度譜的高度是單邊幅度譜的一半,在 $n = 0$ 兩者高度的絕對值相等。

 


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM