Advertisement

数字信号处理——第三章、离散傅里叶变换(DFT)

阅读量:

DFT的实质是将有限长序列傅里叶变换的有限点离散采样 ,从而实现了频域离散化

3.1 离散傅里叶变换的定义及物理意义

3.1.1 DFT的定义

长度M的有限长序列x(n)的N点\textcolor{Red}{傅里叶变换} 为:
X(k)=DFT[x(n)]=\sum_{n=0}^{N-1}x(n)W_N^{kn} \,\,\,\,\,\,\,\,\,\,\,k=0,1,2,3,\ldots,N-1
\textcolor{Red}{离散傅里叶逆变换为:}
x(n)=IDFT[X(k)]=\frac{1}{N}\sum_{k=0}^{N-1}X(k)W_N^{-kn} \,\,\,\,\,\,\,\,\,\,\,n=0,1,2,3,\ldots,N-1\\ \\ W_N=e^{-j\frac{2\pi}{N}}

3.1.2 DFT与傅里叶变换和Z变换的关系

X(k)=X(z)|_{z=e^{j\frac{2\pi }{N}k}}\,\,\,\,\,\,\,\,\,k=0,1,2,3,\ldots,N-1\\ X(k)=X(e^{jw})|_{w=\frac{2\pi }{N}k}\,\,\,\,\,\,\,\,\,k=0,1,2,3,\ldots,N-1\\

3.1.3 DFT的隐含周期性

X(k)和x(n)隐含周期性,且周期均为N,对任意整数m,总有:
W_N^k=W_N^{k+mN}\,\,\,\,\,\,\,\,\,\,k,m为整数,N为自然数
对于任何周期为N的周期序列\hat{x}(n)都可以看作有限长序列x(n)的延拓,x(n)为\hat{x}(n)的一个周期
\hat{x}(n)=\sum_{m=-\infty}^{\infty}x(n+mN)\\ x(n)=\hat{x}(n)R_N(n)

3.2 离散傅里叶变换的基本性质

性质 序列 离散傅里叶变换 备注
线性性质 ax_1(n)+bx_2(n) aX_1(k)+bX_2(k)
序列循环移位性质 x(n) y(n)=x((n+m))_NR_N(n)
时域循环移位定理 \hat{x}(n+m) X(k)W_N^{-km}=X(k)e^{j\frac{2\pi}{N}mk}
频域循环移位定理 W_N^{nl}x(n) \hat{X}(k+l)
周期卷积和定理 \sum_{m=0}^{N-1}x_1(m)x_2{(n-m)} X_1(k)X_2(k)
频域卷积和定理 x_1{(n)}x_2(n) \frac{1}{N}\sum_{l=0}^{N-1}x_1(l)x_2{(k-l)}
复共轭序列 x^*(n) X^*(N-k)\,\,\,\,\,\,0<=k<=N-1 x^*(n)是x(n)的复共轭序列,长度为N

3.3 频率域采样

设序列的Z变换为:
X(z)=\sum_{n=-\infty}^{\infty}x(n)z^{-n}
且X(z)的收敛于包含单位圆。在单位圆 上对X(z)进行等间隔采集N点得到:
X(k)=X(z)|_{z=e^{j\frac{2\pi}{N}k}}=\sum_{n=0}^{N-1}x(n)W_N^{kn}\,\,\,\,\,\,\,\,\,\,\,\,(3.3.1)
将X(k)看做长度为N的有限长序列x_N(n)的DFT,即:x_N(n)=IDFT(X(k))

x_N(n)的周期延拓序列为\hat{x}(n),故X(k)是\hat{X}(k)的主序列
\hat{x}(n)=x_N((n))_N=IDFS[\hat{X}(k)]=\frac{1}{N}\sum_{k=0}^{N-1}X(k)W_N^{-kn}\\ 将3.3.1代入

x_N(n)=\hat{x}(n)R_N(n)=\sum_{i=-\infty}^{\infty}x(n+iN)R_N(n)

上式说明,X(z)在单位圆上的N点等间隔采样X(k)的N点IDFT是原序列x(n)以N为周期的周期性延拓序列的主值序列。如果序列的长度为M,则只有当频域\textcolor{Red}{采样点数N>=M} 时,才有:
x_N(n)=IDFT[X(k)]=x(n)
即可有频域采样X(k)恢复原序列,否则长生时域混叠现象。

3.4 DFT的应用

L点循环卷积:

y_c(n)=\sum_{m=0}^{L-1}h(m)x((n-m))_LR_L(n)

栅栏效应

在DFT中,采样点之间的频谱时看不到的,这就是栅栏效应。为了把原来被“栅栏”挡住的频谱分量检测出来,就必须提高频率的分辨率。

对于有限长序列,可以在原序列的尾部补零;对于无限长序列,可以最大截取长度及DFT变换区间长度,使得频域采样间隔变小。

截断效应

在实际遇到的序列可能时无限长的,用DFT对其进行谱分析时,需要通过加窗函数来对序列进行截短。加粗样式

全部评论 (0)

还没有任何评论哟~