我們在前一章中,已經說明了一個 LTI 系統的頻率響應,也就是其單位脈衝響應的離散時間傅立葉轉換,可用來表示對於特定頻率弦波所產生的震幅和相角改變。此外,對於任意的離散時間訊號 x[n],其「離散時間傅立葉轉換」(簡稱 DTFT)及其反轉換可以表示如下:
由上述的第一個式子即可推導出第二個式子。在第二個式子中,由於 X(ejw) 的週期性是 2p,因此積分區間可以是任何長度為 2p 的區間。
X(ejw) = Ske-jwk x[k] x[n] = (2p)-1∫2pX(ejw)ejwn dw 由第二個式子,可以我們可以將 x[n] 拆解成無窮個基本函數 ejwn 的線性組合,而每一個基本函數的震幅則是由 X(ejw) 來控制。換句話說,X(ejw) 即代表 x[n] 在連續頻率 w 的分量,可以進一步說明如下。
由基本微積分可知,對於一個函數 f(x) 的積分,我們可以使用累加來逼近:
因此如果 x[n] 是實數,我們就可以將其拆解開來,表示如下:
∫02p f(x) dx = limN→∞Sk=0N-1 f(k•(2p/N))•(2p/N) 換句話說,x[n] 已經被拆解成 N 個餘弦函數的線性組合,而這些餘弦函數的角頻率是從 0 到 2p(N-1)/N,震幅則是 |X(ejw)|/N 在每個角頻率的值。(如果 x[n] 是複數,我們也可以進行此種拆解,只不過物理意義較不明顯。)
x[n] = Re{(2p)-1∫2pX(ejw)ejwn dw} = (2p)-1∫2p Re{X(ejw)ejwn} dw = (2p)-1∫2p |X(ejw)| cos(wn + q) dw, q=∠X(ejw) = N-1 Sk=0N-1 [|X(ejw)| cos(wn + q)]w=2pk/N dw, N→∞ 由上一節及本節的討論,我們就可以瞭解,DTFT 具有兩個非常重要的意義:
- 如果 h[n] 是一個 LTI 系統的脈衝響應(Impulse Response),那麼 H(ejw)
= Skh[k]e-jwk 就可以代表此系統對於角頻率等於 w 的訊號所造成的增益 |H(ejw)| 和相位 ∠H(ejw)。 - 如果訊號 x[n] 是一個任意訊號,那麼 X(ejw)
= Skx[k]e-jwk 就可以代表此訊號在不同角頻率 w 的分量大小 |H(ejw)| 和相位 ∠H(ejw)。 以下列出離散時間傅立葉轉換的一些常用到的重要性質。
如果 x[n] 是一個實數序列,那麼我們可以將對應的 DTFT X(ejw) 拆解成實部和虛部的組合,如下:
- 線性性質:
z[n] = a x[n] + b y[n] ←→ Z(ejw) = a X(ejw) + b Y(ejw) - 週期性:
X(ej(w + 2 k p)) = X(ejw) 換句話說,DTFT 的週期為 2p。- 延遲性質:
y[n] = x[n-k] ←→ Y(ejw) = ejwk X(ejw) - 旋積:
y[n] = x[n]*h[n] ←→ Z(ejw) = X(ejw)Y(ejw) 也就是說,在時域的旋積等於在頻域的乘積。其中
X(ejw) = Ske-jwk x[k] = Skx[k] cos(wk) - j Skx[k] sin(wk) = XR(ejw) + j XR(ejw) XR(ejw) = Skx[k] cos(wk) 這些函數滿足下列性質:
XI(ejw) = - Skx[k] sin(wk)
- 共軛對稱:X*(ejw) = X(e-jw)
- XR(ejw) 是偶函數
- XI(ejw) 是奇函數
- |X(ejw)| 是偶函數
- ∠X(ejw) 是奇函數
Audio Signal Processing and Recognition (音訊處理與辨識)