[傅里葉變換及其應用學習筆記] 二. 周期性,三角函數表示復雜函數


    這份是本人的學習筆記,課程為網易公開課上的斯坦福大學公開課:傅里葉變換及其應用。

 

這節課目的

如何用像$sin$,$cos$這些簡單的函數來表示復雜周期函數。

 

信號周期化

並不是所有現象都是周期性的,而且即使是周期性的現象(時間周期性),最終都會終結。而$sin$,$cos$這些數學函數是無始無終的,那么我們該怎么做?

我們采用了一種叫信號周期化的方法:

設有如下信號(左)

image

我們可以把它無限復制,這樣就成了一個周期信號,然后研究我們感興趣的部分(單一周期內的信號)。

由於有了信號周期化這種做法,我們的傅里葉研究將相當廣泛。

 

設定周期

為了方便我們后面的學習,在此設定周期為1,后面的學習會遵循該設定,即

$f(t+1) = f(t)$

因此信號模型為$sin(2\pi t)$與$cos(2\pi t)$。

 

結論

首先引出結論,周期為1的信號,可以由$sin(2\pi t)$或$cos(2\pi t)$組成

 

一個周期,多個頻率

舉個例子

下圖分別為$sin(2\pi t)$,$sin(4\pi t)$,$sin(6 \pi t)$的圖形

sin2pit   sin4pit   sin6pit

$sin(2\pi t)$的周期是1,頻率是1。

$sin(4\pi t)$的周期是1/2,頻率是2,但是1也可以是它的周期。

$sin(6\pi t)$的周期是1/3,頻率是3,但是1也可以是它的周期。

把他們組合起來(相加)得到$sin(2\pi t)+sin(4 \pi t)+sin(6\pi t)$,圖形如下

sin(246pit)

這個復雜的圖形的周期還是1,它是由周期為1,頻率不同的sin函數組成的。

 

上面的例子只是不同頻率的組合,我們還可以改變他們的振幅,相位。這表明我們通過$sin$已經可以組成非常多的信號

$\displaystyle{\sum^n_{k=1}}A_k sin(2\pi kt+\varphi_k)$

注:k=1的項被稱為基波(fundamental wave),k>1的項被稱為諧波(harmonic)

 

公式推導

對sin進行分解

$sin(2\pi kt + \varphi_k)=sin(2\pi kt)cos\varphi_k+cos(2\pi kt)sin\varphi_k$

因此有

$\begin{align*}
&\quad \sum^n_{k=1}A_ksin(2\pi kt + \varphi_k)\\
&=\sum^n_{k=1}A_ksin(2\pi kt)cos\varphi_k+cos(2\pi kt)sin\varphi_k\\
&=\sum^n_{k=1}(a_kcos(2\pi kt)+b_ksin(2\pi kt))
\end{align*}$

$a_k,b_k$與相位$\varphi_k$和振幅$A_k$有關。

另外,我們還可以添加一個常量來表示其中不變的部分:

$\frac{a_0}{2}+\displaystyle{\sum^n_{k=1}}(a_kcos(2\pi kt)+b_ksin(2\pi kt))$

該常量$\frac{a_0}{2}$被稱為直流分量(DC component)。

 

復指數式

上面的式子還可以推導成復指數的方式

有如下歐拉公式:

$e^{2\pi ikt} = cos(2\pi kt)+isin(2\pi kt), i=\sqrt{-1}$

$cos(2\pi kt) = \frac{e^{2\pi ikt} + e^{-2\pi ikt}}{2}$

$sin(2\pi kt) = \frac{e^{2\pi ikt} - e^{-2\pi ikt}}{2i}$

 

通過歐拉公式對上述式子進行展開,得

$\begin{align*}
&\quad a_kcos(2\pi kt)+b_ksin(2\pi kt)\\
&= \frac{a_ke^{2\pi ikt}+a_ke^{-2\pi ikt}}{2}+\frac{b_ke^{2\pi ikt}-b_ke^{-2\pi ikt}}{2i}\\
&= \frac{a_ke^{2\pi ikt}+a_ke^{-2\pi ikt}}{2}+\frac{-b_kie^{2\pi ikt}+b_kie^{-2\pi ikt}}{2}\\
&= \frac{a_k-b_ki}{2}e^{2\pi ikt}+\frac{a_k+b_ki}{2}e^{-2\pi ikt}
\end{align*}$

分成$\frac{a_k-b_ki}{2}e^{2\pi ikt}$與$\frac{a_k+b_ki}{2}e^{-2\pi ikt}$兩部分。按照我們前面的推論,$k$作為調整頻率的系數,是一個正整數,現在如果我們把復指數上的符號移動到$k$上,$k$就稱為了覆蓋正負的整數,那么上面的式子就變成

$\begin{align*}
a_kcos(2\pi kt)+b_ksin(2\pi kt) 
&= \underbrace{\frac{a_k-b_ki}{2}e^{2\pi ikt}+\frac{a_k+b_ki}{2}e^{-2\pi ikt}}_{k>0}\\
&= \underbrace{\frac{a_k-b_ki}{2}e^{2\pi ikt}}_{k>0}+\underbrace{\frac{a_{-k}+b_{-k}i}{2}e^{2\pi ikt}}_{k<0}
\end{align*}$

把$\frac{a_k-b_ki}{2}$和$\frac{a_{-k}+b_{-k}i}{2}$取出來用$C_k$表示,則有,

$C_k=
\begin{cases}
&\frac{a_k-b_ki}{2} \text{ , } k>0 \\
&\frac{a_{-k}+b_{-k}i}{2} \text{ , } k<0
\end{cases}$

即$C_k$為復數且滿足以下條件,

$C_{-k}=\bar{C_k}$

 

有了上述條件,式子可以寫成

$\begin{align*}
&\quad \sum^n_{k=1}(a_kcos(2\pi kt)+b_ksin(2\pi kt))\\
&=\sum^n_{k=-n}C_ke^{2\pi ikt}
\end{align*}$

 

上述推導引出一個結論:對於一個真實的信號(值為實數),當它轉換為上述復數形式時,它的系數對稱存在,即有$k$必然會有$-k$,且$C_k$與$C_{-k}$共軛。反過來,如果系數滿足上述條件,那么此信號也是真實信號。

 

通用性

我們已經從sin的組合推導到了復指數之和的形式。那么說回來,這種三角函數的組合形式是否可以用到更大的范圍?它是否適用於一般周期函數?

下面,我們假設這個推斷是成立的,三角函數之和適用於一般周期函數,則有,

$f(t)=\displaystyle{\sum^n_{k=-n}}C_ke^{2\pi ikt}$

取出該多項式其中的一項$C_me^{2\pi imt},-n \leqslant m \leqslant n$,

$C_me^{2\pi imt} = f(t)-\displaystyle{\sum^n_{k\neq m}}C_k e^{2\pi ikt}$

等號兩邊同時乘以$e^{-2\pi imt}$,得

$\begin{align*}
& C_m = e^{-2\pi imt}f(t)-\sum^n_{k\neq m}C_k e^{-2\pi imt}e^{2\pi ikt}\\
&\quad \  = e^{-2\pi imt}f(t)-\sum^n_{k\neq m}C_k e^{2\pi i(k-m)t}
\end{align*}$

對等號兩邊同時積分

$\displaystyle{\int_{0}^{1}}C_mdt=C_m$

 

$\begin{align*}
&\quad \int_{0}^{1}(e^{-2\pi imt}f(t)-\sum^n_{k\neq m}C_k e^{2\pi i(k-m)t})dt \\
&= \int_{0}^{1}e^{-2\pi imt}f(t)dt - \sum^n_{k\neq m}C_k\int_{0}^{1} e^{2\pi i(k-m)t}dt \\
&= \int_{0}^{1}e^{-2\pi imt}f(t)dt - \sum^n_{k\neq m}C_k\left.\frac{1}{2\pi i(k-m)}e^{2\pi i(k-m)t}\right|^1_0 \\
&= \int_{0}^{1}e^{-2\pi imt}f(t)dt - \sum^n_{k\neq m}C_k\frac{1}{2\pi i(k-m)}(e^{2\pi (k-m)t}-e^0) \\
&= \int_{0}^{1}e^{-2\pi imt}f(t)dt - \sum^n_{k\neq m}C_k\frac{1}{2\pi i(k-m)}(cos2\pi(k-m)+isin2\pi(k-m) - 1) \quad spread \ with \ Euler \ Formular \\
&= \int_{0}^{1}e^{-2\pi imt}f(t)dt - \sum^n_{k\neq m}C_k\frac{1}{2\pi i(k-m)}(1+0-1) \quad k \ and \ m \ is \ interger \\
&= \int_{0}^{1}e^{-2\pi imt}f(t)dt
\end{align*}$

即,

$C_m = \displaystyle{\int_{0}^{1}}e^{-2\pi imt}f(t)dt$


免責聲明!

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



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