1. 离散时间傅里叶变换的导出

针对离散时间非周期序列,为了建立它的傅里叶变换表示,我们将采用与连续情况下完全类似的步骤进行。

考虑某一序列 x [ n ] x[n] x[n],它具有有限持续期;也就是说,对于某个整数 N 1 N_1 N1 N 2 N_2 N2,在 N 1 N N 2 -N_1\leqslant N\leqslant N_2 N1NN2 以外, x [ n ] = 0 x[n]=0 x[n]=0。下图给出了这种类型的一个信号。

由这个非周期信号可以构成一个周期序列 <mover accent="true"> x ~ </mover> [ n ] \tilde x[n] x~[n],使 x [ n ] x[n] x[n] 就是 <mover accent="true"> x ~ </mover> [ n ] \tilde x[n] x~[n] 的一个周期。随着 N N N 的增大, x [ n ] x[n] x[n] 就在一个更长的时间间隔内与 <mover accent="true"> x ~ </mover> [ n ] \tilde x[n] x~[n] 相一致。而当 N N\to \infty N,对任意有限时间值 n n n 而言,有 <mover accent="true"> x ~ </mover> [ n ] = x [ n ] \tilde x[n]=x[n] x~[n]=x[n]

现在我们来考虑一下 <mover accent="true"> x ~ </mover> [ n ] \tilde x[n] x~[n] 的傅里叶级数表示式

<mlabeledtr> <mtext> (1) </mtext> <mover accent="true"> x ~ </mover> [ n ] = <munder> k = ( N ) </munder> a k e j k ( 2 π / N ) n </mlabeledtr> \tag{1}\tilde x[n] = \sum_{k=(N)}a_ke^{jk{(2\pi/N)}n} x~[n]=k=(N)akejk(2π/N)n(1)

<mlabeledtr> <mtext> (2) </mtext> a k = 1 N <munder> n = ( N ) </munder> <mover accent="true"> x ~ </mover> [ n ] e j k ( 2 π / N ) n </mlabeledtr> \tag{2}a_k = \frac{1}{N} \sum_{n=(N)} \tilde x[n]e^{-jk{(2\pi/N)}n} ak=N1n=(N)x~[n]ejk(2π/N)n(2)

因为在 N 1 N N 2 -N_1 \leqslant N \leqslant N_2 N1NN2 区间的一个周期上 <mover accent="true"> x ~ </mover> [ n ] = x [ n ] \tilde x[n]=x[n] x~[n]=x[n],因此我们将上式的求和区间就选在这个周期上

<mlabeledtr> <mtext> (3) </mtext> a k = 1 N <munderover> n = N 1 N 2 </munderover> x [ n ] e j k ( 2 π / N ) n = 1 N <munderover> n = + </munderover> x [ n ] e j k ( 2 π / N ) n </mlabeledtr> \tag{3}a_k = \frac{1}{N} \sum_{n=-N_1}^{N_2} x[n]e^{-jk{(2\pi/N)}n} = \frac{1}{N} \sum_{n=-\infty}^{+\infty} x[n]e^{-jk{(2\pi/N)}n} ak=N1n=N1N2x[n]ejk(2π/N)n=N1n=+x[n]ejk(2π/N)n(3)

现定义函数

<mlabeledtr> <mtext> (4) </mtext> X ( e j ω ) = <munderover> n = + </munderover> x [ n ] e j ω n </mlabeledtr> \tag{4}X(e^{j\omega})=\sum_{n=-\infty}^{+\infty}x[n]e^{-j\omega n} X(ejω)=n=+x[n]ejωn(4)

可见这些系数 a k a_k ak 正比于 X ( e j ω ) X(e^{j\omega}) X(ejω) 的各样本值,即

<mlabeledtr> <mtext> (5) </mtext> a k = 1 N X ( e j k ω 0 ) </mlabeledtr> \tag{5}a_k = \frac{1}{N}X(e^{jk\omega_0}) ak=N1X(ejkω0)(5)

式中, ω 0 = 2 π / N \omega_0=2\pi/N ω0=2π/N 用来记作在频域中的样本间隔。将(1) 和 (5)结合在一起, <mover accent="true"> x ~ </mover> [ n ] \tilde x[n] x~[n] 就可以表示为

<mlabeledtr> <mtext> (6) </mtext> <mover accent="true"> x ~ </mover> [ n ] = <munder> k = ( N ) </munder> 1 N X ( e j k ω 0 ) e j k ω 0 n = 1 2 π <munder> k = ( N ) </munder> X ( e j k ω 0 ) e j k ω 0 n ω 0 </mlabeledtr> \tag{6}\tilde x[n] = \sum_{k=(N)} \frac{1}{N}X(e^{jk\omega_0})e^{jk\omega_0n} = \frac{1}{2\pi}\sum_{k=(N)} X(e^{jk\omega_0})e^{jk\omega_0n}\omega_0 x~[n]=k=(N)N1X(ejkω0)ejkω0n=2π1k=(N)X(ejkω0)ejkω0nω0(6)

随着 N N\to \infty N <mover accent="true"> x ~ </mover> [ n ] \tilde x[n] x~[n] 趋近于 x [ n ] x[n] x[n],式(6)的极限就变成 x [ n ] x[n] x[n] 的表达式。再者,当 N N\to \infty N 时,有 ω 0 0 \omega_0\to 0 ω00,式(6)的右边就过渡为一个积分。

右边的每一项都可以看作是高度为 X ( e j k ω 0 ) e j k ω 0 n X(e^{jk\omega_0})e^{jk\omega_0n} X(ejkω0)ejkω0n 宽度为 ω 0 \omega_0 ω0 的矩形的面积。而且,因为这个求和是在 N N N ω 0 = 2 π / N \omega_0=2\pi/N ω0=2π/N 的间隔内完成的,所以总的积分区间总是有一个 2 π 2\pi 2π 的宽度。式(6)和式(4)就分别变成

<mlabeledtr> <mtext> (7) </mtext> <menclose notation="box"> <mstyle displaystyle="false" scriptlevel="0"> <mstyle displaystyle="true" scriptlevel="0"> x [ n ] = 1 2 π 2 π X ( e j ω ) e j ω n d ω </mstyle> </mstyle> </menclose> </mlabeledtr> \tag{7}\boxed{ x[n]=\frac{1}{2\pi}\int_{2\pi} X(e^{j\omega})e^{j\omega n}d\omega} x[n]=2π12πX(ejω)ejωndω(7)

<mlabeledtr> <mtext> (8) </mtext> <menclose notation="box"> <mstyle displaystyle="false" scriptlevel="0"> <mstyle displaystyle="true" scriptlevel="0"> X ( e j ω ) = <munderover> n = + </munderover> x [ n ] e j ω n </mstyle> </mstyle> </menclose> </mlabeledtr> \tag{8}\boxed{X(e^{j\omega})=\sum_{n=-\infty}^{+\infty}x[n]e^{-j\omega n}} X(ejω)=n=+x[n]ejωn(8)

(7)式和 (8)式被称为离散时间傅里叶变换对。函数 X ( j ω ) X(j\omega) X(jω) 称为 X ( t ) X(t) X(t)离散时间傅里叶变换,也通常被称为频谱

  • 例 1

  • 例 2

2. 周期信号的傅里叶变换

考虑如下信号
<mlabeledtr> <mtext> (9) </mtext> x [ n ] = e j ω 0 n </mlabeledtr> \tag{9} x[n] = e^{j\omega_0 n} x[n]=ejω0n(9)

其傅里叶变换是如下的冲激串

<mlabeledtr> <mtext> (10) </mtext> X ( e j ω ) = <munderover> l = + </munderover> 2 π δ ( ω ω 0 2 π l ) </mlabeledtr> \tag{10}X(e^{j\omega}) = \sum_{l=-\infty}^{+\infty}2\pi\delta(\omega-\omega_0-2\pi l) X(ejω)=l=+2πδ(ωω02πl)(10)

为了验证该式,必须求出其对应的反变换

<mlabeledtr> <mtext> (11) </mtext> 1 2 π 2 π X ( e j ω ) e j ω n d ω = 1 2 π 2 π <munderover> l = + </munderover> 2 π δ ( ω ω 0 2 π l ) e j ω n d ω </mlabeledtr> \tag{11} \frac{1}{2\pi}\int_{2\pi} X(e^{j\omega}) e^{j\omega n}d\omega = \frac{1}{2\pi}\int_{2\pi} \sum_{l=-\infty}^{+\infty}2\pi\delta(\omega-\omega_0-2\pi l) e^{j\omega n}d\omega 2π12πX(ejω)ejωndω=2π12πl=+2πδ(ωω02πl)ejωndω(11)

注意,在任意一个长度为 2 π 2\pi 2π 的积分区间内,在上式的和中真正包括的只有一个冲激,因此,如果所选的积分区间包含在 ω 0 + 2 π r \omega_0+2\pi r ω0+2πr 处的冲激,那么

<mlabeledtr> <mtext> (12) </mtext> 1 2 π 2 π X ( e j ω ) e j ω n d ω = e j ( ω 0 + 2 π r ) n = e j ω 0 n </mlabeledtr> \tag{12} \frac{1}{2\pi}\int_{2\pi} X(e^{j\omega}) e^{j\omega n}d\omega = e^{j(\omega_0+2\pi r) n} = e^{j\omega_0 n} 2π12πX(ejω)ejωndω=ej(ω0+2πr)n=ejω0n(12)

现在考虑一周期序列 x [ n ] x[n] x[n],周期为 N N N,其傅里叶级数为

<mlabeledtr> <mtext> (13) </mtext> x [ n ] = <munder> k = ( N ) </munder> a k e j k ( 2 π / N ) n </mlabeledtr> \tag{13} x[n] = \sum_{k=(N)} a_k e^{jk(2\pi/N)n} x[n]=k=(N)akejk(2π/N)n(13)

这时,傅里叶变换就是

<mlabeledtr> <mtext> (14) </mtext> X ( e j ω ) = <munderover> k = + </munderover> 2 π a k δ ( ω 2 π k N ) = <munderover> l = + </munderover> <munder> k = ( N ) </munder> 2 π a k δ ( ω k ω 0 2 π l ) </mlabeledtr> \tag{14} X(e^{j\omega}) = \sum_{k=-\infty}^{+\infty} {2\pi} a_k \delta (\omega-\frac{2\pi k}{N}) =\sum_{l=-\infty}^{+\infty} \sum_{k=(N)}{2\pi} a_k \delta (\omega - k\omega_0 - 2\pi l) X(ejω)=k=+2πakδ(ωN2πk)=l=+k=(N)2πakδ(ωkω02πl)(14)

这样,一个周期信号的傅里叶变换就能直接从它的傅里叶级数系数得到。

3. 离散时间傅里叶变换性质

为了方便,我们将 x [ n ] x[n] x[n] X ( e j ω ) X(e^{j\omega}) X(ejω) 这一对傅里叶变换用下列符号表示

x [ n ] <mover> <mstyle displaystyle="true" scriptlevel="0"> F </mstyle> </mover> X ( e j ω ) x[n] \overset{{\displaystyle {\mathcal {F}}}}{\leftrightarrow} X(e^{j\omega}) x[n]FX(ejω)

3.1. 离散时间傅里叶变换的周期性

<mlabeledtr> <mtext> (15) </mtext> <menclose notation="box"> <mstyle displaystyle="false" scriptlevel="0"> <mstyle displaystyle="true" scriptlevel="0"> X ( e j ( ω + 2 π ) ) = X ( e j ω ) </mstyle> </mstyle> </menclose> </mlabeledtr> \tag{15} \boxed{ X(e^{j(\omega+2\pi)}) = X(e^{j\omega})} X(ej(ω+2π))=X(ejω)(15)

3.2. 线性

x 1 [ n ] <mover> <mstyle displaystyle="true" scriptlevel="0"> F </mstyle> </mover> X 1 ( e j ω ) x_1[n] \overset{{\displaystyle {\mathcal {F}}}}{\leftrightarrow} X_1(e^{j\omega}) x1[n]FX1(ejω)

x 2 [ n ] <mover> <mstyle displaystyle="true" scriptlevel="0"> F </mstyle> </mover> X 2 ( e j ω ) x_2[n] \overset{{\displaystyle {\mathcal {F}}}}{\leftrightarrow} X_2(e^{j\omega}) x2[n]FX2(ejω)

<mlabeledtr> <mtext> (16) </mtext> <menclose notation="box"> <mstyle displaystyle="false" scriptlevel="0"> <mstyle displaystyle="true" scriptlevel="0"> a x 1 [ n ] + b x 2 [ n ] <mover> <mstyle displaystyle="true" scriptlevel="0"> F </mstyle> </mover> a X 1 ( e j ω ) + b X 2 ( e j ω ) </mstyle> </mstyle> </menclose> </mlabeledtr> \tag{16} \boxed{ ax_1[n]+bx_2[n] \overset{{\displaystyle {\mathcal {F}}}}{\leftrightarrow} aX_1(e^{j\omega})+bX_2(e^{j\omega})} ax1[n]+bx2[n]FaX1(ejω)+bX2(ejω)(16)

3.3. 时移与频移性质

x [ n ] <mover> <mstyle displaystyle="true" scriptlevel="0"> F </mstyle> </mover> X ( e j ω ) x[n] \overset{{\displaystyle {\mathcal {F}}}}{\leftrightarrow} X(e^{j\omega}) x[n]FX(ejω)

<mlabeledtr> <mtext> (17) </mtext> <menclose notation="box"> <mstyle displaystyle="false" scriptlevel="0"> <mstyle displaystyle="true" scriptlevel="0"> x [ n n 0 ] <mover> <mstyle displaystyle="true" scriptlevel="0"> F </mstyle> </mover> e j ω n 0 X ( e j ω ) </mstyle> </mstyle> </menclose> </mlabeledtr> \tag{17} \boxed{ x[n-n_0] \overset{{\displaystyle {\mathcal {F}}}}{\leftrightarrow} e^{-j\omega n_0}X(e^{j\omega})} x[nn0]Fejωn0X(ejω)(17)

<mlabeledtr> <mtext> (18) </mtext> <menclose notation="box"> <mstyle displaystyle="false" scriptlevel="0"> <mstyle displaystyle="true" scriptlevel="0"> e j ω 0 n x [ n ] <mover> <mstyle displaystyle="true" scriptlevel="0"> F </mstyle> </mover> X ( e j ( ω ω 0 ) ) </mstyle> </mstyle> </menclose> </mlabeledtr> \tag{18} \boxed{ e^{j\omega_0 n}x[n] \overset{{\displaystyle {\mathcal {F}}}}{\leftrightarrow} X(e^{j(\omega-\omega_0)})} ejω0nx[n]FX(ej(ωω0))(18)

3.4. 共轭及共轭对称性

x [ n ] <mover> <mstyle displaystyle="true" scriptlevel="0"> F </mstyle> </mover> X ( e j ω ) x[n] \overset{{\displaystyle {\mathcal {F}}}}{\leftrightarrow} X(e^{j\omega}) x[n]FX(ejω)

<mlabeledtr> <mtext> (19) </mtext> <menclose notation="box"> <mstyle displaystyle="false" scriptlevel="0"> <mstyle displaystyle="true" scriptlevel="0"> x [ n ] <mover> <mstyle displaystyle="true" scriptlevel="0"> F </mstyle> </mover> X ( e j ω ) </mstyle> </mstyle> </menclose> </mlabeledtr> \tag{19} \boxed{ x^*[n] \overset{{\displaystyle {\mathcal {F}}}}{\leftrightarrow} X^*(e^{-j\omega})} x[n]FX(ejω)(19)

共轭性质就能证明,若 x ( t ) x(t) x(t) 为实函数,那么 X ( j ω ) X(j\omega) X(jω) 就具有共轭对称性,即

<mlabeledtr> <mtext> (20) </mtext> <menclose notation="box"> <mstyle displaystyle="false" scriptlevel="0"> <mstyle displaystyle="true" scriptlevel="0"> X ( e j ω ) = X ( e j ω ) [ x [ n ] ] </mstyle> </mstyle> </menclose> </mlabeledtr> \tag{20} \boxed{ X(e^{j\omega}) = X^*(e^{-j\omega}) \qquad [x[n] 为实]} X(ejω)=X(ejω)[x[n]](20)

这就是说,离散傅里叶变换的实部是频率的偶函数,而虚部则是频率的奇函数

3.5. 差分与累加

<mlabeledtr> <mtext> (21) </mtext> <menclose notation="box"> <mstyle displaystyle="false" scriptlevel="0"> <mstyle displaystyle="true" scriptlevel="0"> x [ n ] x [ n 1 ] <mover> <mstyle displaystyle="true" scriptlevel="0"> F </mstyle> </mover> ( 1 e j ω ) X ( e j ω ) </mstyle> </mstyle> </menclose> </mlabeledtr> \tag{21} \boxed{ x[n]-x[n-1] \overset{{\displaystyle {\mathcal {F}}}}{\leftrightarrow} (1-e^{-j\omega}) X(e^{j\omega})} x[n]x[n1]F(1ejω)X(ejω)(21)

<mlabeledtr> <mtext> (22) </mtext> <menclose notation="box"> <mstyle displaystyle="false" scriptlevel="0"> <mstyle displaystyle="true" scriptlevel="0"> <munderover> m = n </munderover> x [ m ] <mover> <mstyle displaystyle="true" scriptlevel="0"> F </mstyle> </mover> 1 1 e j ω X ( e j ω ) + π X ( e j 0 ) <munderover> k = + </munderover> δ ( ω 2 π k ) </mstyle> </mstyle> </menclose> </mlabeledtr> \tag{22} \boxed{ \sum_{m=-\infty}^{n}x[m] \overset{{\displaystyle {\mathcal {F}}}}{\leftrightarrow} \frac{1}{1-e^{-j\omega}} X(e^{j\omega})+\pi X(e^{j0}) \sum_{k=-\infty}^{+\infty}\delta(\omega-2\pi k)} m=nx[m]F1ejω1X(ejω)+πX(ej0)k=+δ(ω2πk)(22)

3.6. 时间反转

<mlabeledtr> <mtext> (23) </mtext> <menclose notation="box"> <mstyle displaystyle="false" scriptlevel="0"> <mstyle displaystyle="true" scriptlevel="0"> x [ n ] <mover> <mstyle displaystyle="true" scriptlevel="0"> F </mstyle> </mover> X ( e j ω ) </mstyle> </mstyle> </menclose> </mlabeledtr> \tag{23} \boxed{ x[-n] \overset{{\displaystyle {\mathcal {F}}}}{\leftrightarrow} X(e^{-j\omega})} x[n]FX(ejω)(23)

3.7. 时域扩展

若令 是一个正整数,并且定义

<mlabeledtr> <mtext> (24) </mtext> x ( k ) [ n ] = { <mstyle displaystyle="false" scriptlevel="0"> x [ n / k ] </mstyle> <mstyle displaystyle="false" scriptlevel="0"> <mtext> 当  </mtext> n <mtext>   </mtext> <mtext>   </mtext> k <mtext>   </mtext> </mstyle> <mstyle displaystyle="false" scriptlevel="0"> 0 , </mstyle> <mstyle displaystyle="false" scriptlevel="0"> <mtext> 当  </mtext> n <mtext>   </mtext> <mtext>   </mtext> k <mtext>   </mtext> </mstyle> </mlabeledtr> \tag{24} x_{(k)}[n] = \begin{cases} x[n/k] &amp;\text 当\space n \space为\space k\space的整数倍 \\ 0, &amp;\text 当\space n \space不为\space k\space的整数倍 \end{cases} x(k)[n]={x[n/k]0, n  k  n  k (24)

<mlabeledtr> <mtext> (25) </mtext> <menclose notation="box"> <mstyle displaystyle="false" scriptlevel="0"> <mstyle displaystyle="true" scriptlevel="0"> x ( k ) [ n ] <mover> <mstyle displaystyle="true" scriptlevel="0"> F </mstyle> </mover> X ( e j k ω ) </mstyle> </mstyle> </menclose> </mlabeledtr> \tag{25} \boxed{ x_{(k)}[n] \overset{{\displaystyle {\mathcal {F}}}}{\leftrightarrow} X(e^{jk\omega})} x(k)[n]FX(ejkω)(25)

3.8. 频域微分

<mlabeledtr> <mtext> (26) </mtext> <menclose notation="box"> <mstyle displaystyle="false" scriptlevel="0"> <mstyle displaystyle="true" scriptlevel="0"> n x [ n ] <mover> <mstyle displaystyle="true" scriptlevel="0"> F </mstyle> </mover> j d X ( e j ω ) d ω </mstyle> </mstyle> </menclose> </mlabeledtr> \tag{26} \boxed{nx[n] \overset{{\displaystyle {\mathcal {F}}}}{\leftrightarrow} j\frac{dX(e^{j\omega})}{d\omega} } nx[n]FjdωdX(ejω)(26)

3.9. 帕斯瓦尔定理

<mlabeledtr> <mtext> (27) </mtext> <menclose notation="box"> <mstyle displaystyle="false" scriptlevel="0"> <mstyle displaystyle="true" scriptlevel="0"> <munderover> + </munderover> <mtext>   </mtext> x [ n ] <mtext>   </mtext> 2 = 1 2 π 2 π X ( e j ω ) 2 d ω </mstyle> </mstyle> </menclose> </mlabeledtr> \tag{27} \boxed{\sum_{-\infty}^{+\infty}|\space x[n] \space |^2 =\frac{1}{2\pi}\int_{2\pi}|X(e^{j\omega})|^2d\omega } + x[n] 2=2π12πX(ejω)2dω(27)

3.10. 卷积性质

<mlabeledtr> <mtext> (28) </mtext> <menclose notation="box"> <mstyle displaystyle="false" scriptlevel="0"> <mstyle displaystyle="true" scriptlevel="0"> y [ n ] = h [ n ] x [ n ] <mover> <mstyle displaystyle="true" scriptlevel="0"> F </mstyle> </mover> Y ( e j ω ) = H ( e j ω ) X ( e j ω ) </mstyle> </mstyle> </menclose> </mlabeledtr> \tag{28} \boxed{y[n]=h[n]*x[n] \overset{{\displaystyle {\mathcal {F}}}}{\leftrightarrow} Y(e^{j\omega})=H(e^{j\omega})X(e^{j\omega})} y[n]=h[n]x[n]FY(ejω)=H(ejω)X(ejω)(28)

两个信号在时域内的卷积就等于它们傅里叶变换的乘积。

3.11. 相乘性质

<mlabeledtr> <mtext> (29) </mtext> <menclose notation="box"> <mstyle displaystyle="false" scriptlevel="0"> <mstyle displaystyle="true" scriptlevel="0"> y [ n ] = x 1 [ n ] x 2 [ n ] <mover> <mstyle displaystyle="true" scriptlevel="0"> F </mstyle> </mover> Y ( e j ω ) = 1 2 π 2 π X 1 ( e j θ ) X 2 ( e j ( ω θ ) ) d θ </mstyle> </mstyle> </menclose> </mlabeledtr> \tag{29} \boxed{y[n]=x_1[n]x_2[n] \overset{{\displaystyle {\mathcal {F}}}}{\leftrightarrow} Y(e^{j\omega})=\frac{1}{2\pi}\int_{2\pi}X_1(e^{j\theta})X_2(e^{j(\omega-\theta)})d\theta} y[n]=x1[n]x2[n]FY(ejω)=2π12πX1(ejθ)X2(ej(ωθ))dθ(29)
两个信号在时域内的相乘就对应于频域内的周期卷积

4. 傅里叶变换性质和基本傅里叶变化列表

5. 离散时间傅里叶变换和连续时间傅里叶级数之间的对偶型

获取更多精彩,请关注「seniusen」!