[離散時間信號處理學習筆記] 6. 離散時間傅里葉變換


要理解這節課的內容需要先對傅里葉變換有一定程度的了解,這里主要分析的是離散時間傅里葉變換,這部分算是從傅里葉變換到離散傅里葉變換的過渡內容。推薦閱讀[傅里葉變換及其應用學習筆記] 課程概覽中離散傅里葉變換開頭的相關課程。

 

離散時間傅里葉變換

離散時間傅里葉變換(Discrete-Time Fourier Transform),簡稱DTFT,DTFT是從傅里葉變換(FT)中的來的。

 

從FT到DTFT

FT的分析對象是時域上的連續時間函數$x(t)$,DTFT的分析對象是對時域上的序列$x[n]$。兩者間有如下關系:

$x[n] = x(n), -\infty<n<\infty$

$x[n]$相當於$x(t)$在$n$上的取樣,不過$x[n]$終究是離散序列,為了使它跟傅里葉變換扯上關系,有必要把$x[n]$轉化成連續時間函數。而數學上可以用原函數與脈沖函數的乘積來表示取樣。

$\begin{align*}
x_s(t) &= x(t) \sum_{n=-\infty}^{\infty}\delta(t-n)\\
&=\sum_{n=-\infty}^{\infty}x(n)\delta(t-n)\\
&=\sum_{n=-\infty}^{\infty}x[n]\delta(t-n)
\end{align*}$

此時,上述取樣仍然為連續時間函數,對它進行傅里葉變換

$\begin{align*}
X_s(e^{j\omega}) = \mathcal{F}x_s(t) &=\sum_{n=-\infty}^{\infty}x[n]\mathcal{F}\delta(t-n)\\
&= \sum_{n=-\infty}^{\infty}x[n]e^{-j\omega n}
\end{align*}$

即得到

$\color{red}{X(e^{j\omega}) = \displaystyle{ \sum_{n=-\infty}^{\infty}x[n]e^{-j\omega n} }}$

這個式子就被稱為DTFT。其中變量$\omega$為頻率。那么相對地,通過$X(e^{j\omega})$來還原序列$x[n]$的式子就被稱為IDTFT(可以從傅里葉級數的式子進行推導得到)

$\color{red}{\displaystyle{x[n] = \frac{1}{2\pi}\int_{-\pi}^{\pi}X(e^{j\omega})e^{j\omega n}d\omega }}$

對於DTFT,套用推導傅里葉級數的思想:在時域上有一個序列$x[n]$,組成它冪級數$e^{j\omega n}$的頻率分布在$(-\pi,\pi)$之間,在頻域上呈現出一個函數$X(e^{j\omega})$。

 

一般來說,序列在進行DTFT后得到的是一個變量為$\omega$的復函數,這點和頻率響應一樣,它可以表示為實部與虛部的形式

$X(e^{j\omega}) = X_R(e^{j\omega})+jX_I(e^{j\omega})$

也可以表示為幅度與相位的形式

$X(e^{j\omega}) = |X(e^{j\omega})|e^{j\angle X(e^{j\omega})}$

 

DTFT與頻率響應

頻率響應有如下定義

$H(e^{j\omega}) = \displaystyle{\sum_{n=-\infty}^{\infty}h[n]e^{-j\omega n} }$

對比上述傅里葉變換,可以發現頻率響應就是單位脈沖響應的DTFT,那么通過對頻率響應進行IDTFT即可得到單位脈沖響應

$\displaystyle{h[n] = \int_{-\pi}^{\pi}H(e^{j\omega})e^{j\omega n}d\omega}$

 

DTFT的收斂性

收斂性問題我們在傅里葉變換課程能也有討論過,這里從序列這邊展開討論。DTFT不是對任何序列都適用的,顯然它對序列有一定的要求。

首先假設有一函數$X(e^{j\omega})$,其對應的序列為$x[n]$

  • 如果$\displaystyle{\sum_{n=-\infty}^{\infty}|x[n]|<\infty}$,則DTFT得到的函數會絕對收斂於$X(e^{j\omega})$,$\displaystyle{X(e^{j\omega}) = X_M(e^{j\omega})=\lim_{M\to\infty}\sum_{n=-M}^{M}x[n]e^{-j\omega n}}$
  • 如果$\displaystyle{\sum_{n=-\infty}^{\infty}|x[n]|^2<\infty}$,則DTFT得到的函數會均方收斂於$X(e^{j\omega})$,$\displaystyle{\lim_{M\to\infty}\int_{-\pi}^{\pi}|X(e^{j\omega})-X_M(e^{j\omega})|^2d\omega = 0}$
  • 對於不收斂的序列$x[n]$,有時也能得到其DTFT,如$u[n]$進行DTFT后可以得到$U(e^{j\omega}) = \frac{1}{1-e^{-j\omega}}+\displaystyle{\sum_{r=-\infty}^{\infty}\pi\delta(\omega+2\pi r)}$

 

 

關於收斂的意義,下面的均方收斂的例子能很好地展示其收斂過程。

有一低通濾波器的頻率響應如下

$H_{lp}(e^{j\omega})=\left\{\begin{matrix}
1, &|\omega|<\omega_c \\
0, & \omega_c<|\omega|\leqslant\pi
\end{matrix}\right.$

單位脈沖響應是頻率響應的逆傅里葉變換

$\begin{align*}
h_{lp}[n] &=\frac{1}{2\pi}\int_{-\omega_c}^{\omega_c}e^{j\omega n}d\omega \\
&= \left [\frac{1}{2\pi jn}e^{j\omega n}\right ]_{-\omega_c}^{\omega_c}\\
&= \frac{1}{2\pi jn}(e^{j\omega_c n}-e^{-j\omega_c n})\\
&= \frac{sin\omega_c n}{\pi n} \quad -\infty<n<\infty
\end{align*}$

可以看到當$n\to\infty$時,這個序列僅以$\frac{1}{n}$趨於0,因此$h_{lp}[n]$不是絕對可加的,也就是說

$\displaystyle{ \sum_{n=-\infty}^{\infty}\frac{sin\omega_c n}{\pi n}e^{-j\omega n} }$

並不絕對收斂於$H_{lp}(e^{j\omega})$。為了獲得直觀上的理解,先來考慮作為有限項和的

$\color{red}{H_{M}(e^{j\omega}) = \displaystyle{\sum_{n=-M}^{M}\frac{sin\omega_c n}{\pi n}e^{-j\omega n}}}$

當$M=1,3,7,19$時分別如下圖所示

 

1     3

 

7     19

可以發現隨着$M$增大,有限項和越來越趨近與原來的頻率響應。但同時也有一個問題,在$\omega_c$點附近的波形一直都是振盪的,並且幅度沒有變小,而是隨着$M$的增大越發集中於$\omega_c$。也就是說對$h_{lp}[n]$進行DTFT后得到的函數不是絕對收斂於$H_{lp}(e^{j\omega})$的,但它是均方收斂的。

 

 

DTFT的對稱性質

對稱性質主要用於簡化計算。

我們現實中所碰到的序列多是實數序列,不過在對信號進行傅里葉分析時,則不可避免地把序列$x[n]$的范圍擴展到復數域。對於復數序列$x[n]$,它可以被分為實數部分$x_R[n]$與虛數部分$x_I[n]$。

complex_series

cos_series     sin_series

 

共軛對稱與共軛反對稱

如果一個序列的實數域對稱,虛數域反(原點)對稱,則稱該序列為共軛對稱序列(conjugate-symmetric sequence),該序列用$x_e[n]$來表示,有$x_e[n] = x_e^{*}[-n]$

 

RealSymmetric     ImaginaryAntisymmetric

如果一個序列的實數域反(原點)對稱,虛數域對稱,則稱該序列為共軛反對稱序列(conjugate-antisymmetric sequence),該序列用$x_o[n]$來表示,有$x_o[n] = -x_o^{*}[n]$

 

RealAntisymmetric     ImaginarySymmetric

任何序列$x[n]$都能表示成一個共軛對稱序列和一個共軛反對稱序列之和

$x[n] = x_e[n]+x_o[n]$

其中

$\left\{\begin{matrix}
x_e[n] &=&\frac{1}{2}(x[n]+x^*[-n]) &=&x_e^*[-n] \\
x_o[n] &=&\frac{1}{2}(x[n]-x^*[-n]) &=&-x_o^*[-n]
\end{matrix}\right.$

利用共軛對稱的定義容易證明上述等式是成立的,利用共軛對稱的圖示也能很好地理解。

 

同理,復函數也可以有共軛對稱這一特征。DTFT得到的復函數$X(e^{j\omega})$也能分解為共軛對稱和共軛反對稱函數之和

$X(e^{j\omega}) = X_e(e^{j\omega})+X_o(e^{j\omega})$

其中

$\left\{\begin{matrix}
X_e(e^{j\omega}) &=&\frac{1}{2}[X(e^{j\omega})+X^*(e^{j\omega})] &=&X_e^*(e^{j\omega}) \\
X_o(e^{j\omega}) &=&\frac{1}{2}[X(e^{j\omega})-X^*(e^{j\omega})] &=&-X_o^*(e^{j\omega})
\end{matrix}\right.$

 

 

DTFT的對稱性質

DTFT有如下表的對稱性質

序列x[n] DTFT X(e)
1. x*[n] X*(e-jω)
2. x*[-n] X*(e)
3. Re{x[n]} Xe(e)  (X(e)的共軛對稱部分)
4. jIm{x[n]} Xo(e)  (X(e)的共軛反對稱部分)
5. xe[n]  (x[n]的共軛對稱部分) XR(e)=Re{X(e)}
6. xo[n]  (x[n]的共軛反對稱部分) jXI(e)=jIm{X(e)}
以下性質僅適用於x[n]為實序列
7. 任意實x[n] X(e)=X*(e-jω)  (共軛對稱)
8. 任意實x[n] XR(e)=XR(e-jω)  (實部為偶函數(對稱))
9. 任意實x[n] XI(e)=-XI(e-jω)  (虛部為奇函數(反對稱))
10. 任意實x[n] |X(e)|=|X(e-jω)|  (幅度為偶函數(對稱))
11. 任意實x[n] ∠X(e)=∠X(e-jω)  (相位為奇函數(反對稱))
12. xe[n]  (x[n]的偶部(對稱部分)) XR(e)
13. xo[n]  (x[n]的奇部(反對稱部分)) jXI(e)

 

上述性質中,只要證明了性質1,同理即可證明性質2,后面的性質可以通過上述共軛對稱之和的等式中容易證明。

性質1證明

對$x^{*}[n]$進行DTFT

$\begin{align*}
\sum_{n=-\infty}^{\infty}x^*[n]e^{-j\omega n}
&=\sum_{n=-\infty}^{\infty}x_R[n](cos\omega n-jsin\omega n)-j\sum_{n=-\infty}^{\infty}x_I[n](cos\omega n-jsin\omega n)\\
&=\sum_{n=-\infty}^{\infty}x_R[n]cos\omega n-j\sum_{n=-\infty}^{\infty}x_R[n]sin\omega n-j\sum_{n=-\infty}^{\infty}x_I[n]cos\omega n-\sum_{n=-\infty}^{\infty}x_I[n]sin\omega n\\
&=\left(\sum_{n=-\infty}^{\infty}x_R[n]cos\omega n-\sum_{n=-\infty}^{\infty}x_I[n]sin\omega n\right)-j\left( \sum_{n=-\infty}^{\infty}x_R[n]sin\omega n+\sum_{n=-\infty}^{\infty}x_I[n]cos\omega n\right)\\
&=\left(\sum_{n=-\infty}^{\infty}|x[n]|cos\phi cos\omega n-\sum_{n=-\infty}^{\infty}|x[n]|sin\phi sin\omega n\right )\\
&\quad-j\left(\sum_{n=-\infty}^{\infty}|x[n]|cos\phi sin\omega n+\sum_{n=-\infty}^{\infty}|x[n]|sin\phi cos\omega n\right)\quad letting\ \phi=\angle x[n]\\
&=\sum_{n=-\infty}^{\infty}|x[n]|cos(\phi+\omega n)-j\sum_{n=-\infty}^{\infty}|x[n]|sin(\phi+\omega n)\\
\end{align*}\\$

對$x[n]$進行DTFT

$\begin{align*}
\sum_{n=-\infty}^{\infty}x[n]e^{-j\omega n}
&= \sum_{n=-\infty}^{\infty}|x[n]|e^{j\angle x[n]}e^{-j\omega n} \\
&= \sum_{n=-\infty}^{\infty}|x[n]|e^{j(\phi-\omega n)}\quad letting\ \phi=\angle x[n]\\
&=\sum_{n=-\infty}^{\infty}|x[n]|cos(\phi-\omega n)+j\sum_{n=-\infty}^{\infty}|x[n]|sin(\phi-\omega n)\\
\end{align*}$

對比兩個結果即可發現$x^*[n]$進行DTFT后得到的是$X^*(e^{-j\omega})$

 

 

DTFT的相關定理

本書是把DTFT描述成序列的傅里葉變換,用的也是傅里葉變換的符號$\mathcal{F}$。

DTFT的相關定理與傅里葉變換相關定理相差無幾,如有必要也能用同樣的方法推導得出,下面只進行定理的羅列

序列x[n],y[n] DTFT X(e),Y(e)
1. $ax[n]+by[n]$ $aX(e^{j\omega})+bY(e^{j\omega})$
2. $x[n-n_d]$ ($n_d$為整數) $e^{-j\omega n_d}X(e^{j\omega})$
3. $e^{j\omega_0 n}x[n]$ $X(e^{j(\omega-\omega_0)})$
4. $x[-n]$ $X(e^{-j\omega})$
5. $nx[n]$ $j\frac{dX(e^{j\omega})}{d\omega}$
6. $x[n]*y[n]$ $X(e^{j\omega})Y(e^{j\omega})$
7. $x[n]y[n]$ $\frac{1}{2\pi}\int_{-\pi}^{\pi}X(e^{j\theta})Y(e^{j(\omega-\theta)})d\omega$
帕斯瓦爾定理:
8. $\displaystyle{\sum_{n=-\infty}^{\infty}|x[n]|^2=\frac{1}{2\pi}\int_{-\pi}^{\pi}|X(e^{j\omega})|^2d\omega }$
9. $\displaystyle{\sum_{n=-\infty}^{\infty}x[n]y^*[n]=\frac{1}{2\pi}\int_{-\pi}^{\pi}X(e^{j\omega})Y^*(e^{j\omega})d\omega }$

其中1為線性性質,2與3為移位定理,4為對偶性質,5為微分定理,6與7為卷積定理。

 

下面列出的是一些常用序列的DTFT

序列 DTFT
1. $\delta[n]$ $1$
2. $\delta[n-n_0]$ $e^{-j\omega n_0}$
3. $1\quad -\infty<n<\infty$ $\displaystyle{ \sum_{k=-\infty}^{\infty}2\pi\delta(\omega+2\pi k) }$
4. $a^nu[n]\quad (|a|<1)$ $\frac{1}{1-ae^{-j\omega}}$
5. $u[n]$ $\frac{1}{1-e^{-j\omega}}+\displaystyle{\sum_{k=-\infty}^{\infty}\pi\delta(\omega+2\pi k)}$
6. $(n+1)a^nu[n]\quad (|a|<1)$ $\frac{1}{(1-ae^{-j\omega})^2}$
7. $\frac{r^nsin\omega_p(n+1)}{sin\omega_p}u[n]\quad (|r|<1)$ $\frac{1}{1-2rcos\omega_p e^{-j\omega}+r^2e^{-j2\omega}}$
8. $\frac{sin\omega_c n}{\pi n}$ $X(e^{j\omega})=\left\{\begin{matrix}1, &|\omega|<\omega_c \\ 0, &\omega_c<|\omega|\leqslant\pi\end{matrix}\right.$
9. $x[n]=\left\{\begin{matrix}1, &0\leqslant n\leqslant M \\ 0, &else\end{matrix}\right.$ $\frac{sin[\omega(M+1)/2]}{sin(\omega/2)}e^{-j\omega M/2}$
10. $e^{j\omega_0 n}$ $\displaystyle{\sum_{k=-\infty}^{\infty}2\pi\delta(\omega-\omega_0+2\pi k)}$
11. $cos(\omega_0 n+\phi)$ $\displaystyle{ \sum_{k=-\infty}^{\infty}[\pi e^{j\phi}\delta(\omega-\omega_0+2\pi k)+\pi e^{-j\phi}\delta(\omega+\omega_0+2\pi k)] }$


免責聲明!

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



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