本文主要从傅里叶级数 FS，连续时间傅里叶变换 CTFT，离散时间傅里叶变换 DTFT，以及离散傅里叶变换 DFT 之间的区别与联系进行了比较详细的讨论，主要注重于公式形式上的推导，略去了相关的图像示例，对于傅里叶级数等数学上的严格证明，以及各种变换的性质也不属于本文的重点。本文主要参考自以下教材。

《数字信号处理–基于计算机的方法（第四版）》Sanjit K. Mitra ,电子工业出版社。

由于本文公式所占用的字符比较多，无法在一篇博客中完整发布，所以将其分为两篇博客。本篇主要介绍了傅里叶级数 FS，连续时间傅里叶变换 CTFT，以及离散时间傅里叶变换 DTFT。对于离散傅里叶变换 DFT 的内容，以及相关的总结，请移步以下链接。

傅里叶级数FS，连续时间傅里叶变换CTFT，离散时间傅里叶变换DTFT，离散傅里叶变换DFT，推导与联系（二）
https://blog.csdn.net/qq_33552519/article/details/130260657

对于离散傅里叶变换所衍生出的离散余弦变换和离散正弦变换，在我的另一篇文章中有所介绍，欢迎阅读与提出相关修改意见。

DFT：https://blog.csdn.net/qq_33552519/article/details/124917126

DCT：https://blog.csdn.net/qq_33552519/article/details/124917473

DST：https://blog.csdn.net/qq_33552519/article/details/124962042

1 周期函数的傅里叶级数

根据高数的知识，设 f ( t ) f\left( t \right) f(t) 是一个以 T T T 为周期的函数，且在 [ − T / 2 , T / 2 ] [ - T/2,T/2] [−T/2,T/2] 上有界可积，我们称数串

a n = 2 T ∫ − T / 2 T / 2 f ( t ) cos ⁡ ( n ω 0 t ) d t , b n = 2 T ∫ − T / 2 T / 2 f ( t ) sin ⁡ ( n ω 0 t ) d t , w h e r e ω 0 = 2 π / T . (1.1) \begin{gathered} {a_n} = \frac{2}{T}\int\limits_{ - T/2}^{T/2} {f\left( t \right)\cos \left( {n{\omega _0}t} \right)dt} , \\ {b_n} = \frac{2}{T}\int\limits_{ - T/2}^{T/2} {f\left( t \right)\sin \left( {n{\omega _0}t} \right)dt} , \\ where\;{\omega _0} = 2\pi /T. \\ \end{gathered} \tag{1.1} an=T2−T/2∫T/2f(t)cos(nω0t)dt,bn=T2−T/2∫T/2f(t)sin(nω0t)dt,whereω0=2π/T.(1.1)

为函数 f ( t ) f\left( t \right) f(t) 的傅里叶系数。以傅里叶系数为系数，所作的三角级数

f ( t ) ∼ a 0 2 + ∑ n = 1 ∞ ( a n cos ⁡ ( n ω 0 t ) + b n sin ⁡ ( n ω 0 t ) ) , (1.2) f\left( t \right) \sim \frac{{{a_0}}}{2} + \sum\limits_{n = 1}^\infty {\left( {{a_n}\cos \left( {n{\omega _0}t} \right) + {b_n}\sin \left( {n{\omega _0}t} \right)} \right)} ,\tag{1.2} f(t)∼2a0+n=1∑∞(ancos(nω0t)+bnsin(nω0t)),(1.2)

称为函数 f ( t ) f\left( t \right) f(t) 的傅里叶级数 (Fourier Series, FS)，其中 cos ⁡ ( n ω 0 t ) \cos \left( {n{\omega _0}t} \right) cos(nω0t) 与 sin ⁡ ( n ω 0 t ) \sin \left( {n{\omega _0}t} \right) sin(nω0t) 称为基本三角函数系。考虑欧拉公式，有

cos ⁡ ( n ω 0 t ) = e j n ω 0 t + e − j n ω 0 t 2 , sin ⁡ ( n ω 0 t ) = e j n ω 0 t − e − j n ω 0 t 2 j . (1.3) \begin{gathered} \cos \left( {n{\omega _0}t} \right) = \frac{{{e^{jn{\omega _0}t}} + {e^{ - jn{\omega _0}t}}}}{2}, \\ \sin \left( {n{\omega _0}t} \right) = \frac{{{e^{jn{\omega _0}t}} - {e^{ - jn{\omega _0}t}}}}{{2j}}. \\ \end{gathered} \tag{1.3} cos(nω0t)=2ejnω0t+e−jnω0t,sin(nω0t)=2jejnω0t−e−jnω0t.(1.3)

且根据式 (1.1)，有

a − n = a n , b − n = − b n . (1.4) \begin{gathered} {a_{ - n}} = {a_n}, \\ {b_{ - n}} = - {b_n}. \\ \end{gathered} \tag{1.4} a−n=an,b−n=−bn.(1.4)

于是

f ( t ) ∼ a 0 2 + ∑ n = 1 ∞ [ a n cos ⁡ ( n ω 0 t ) + b n sin ⁡ ( n ω 0 t ) ] = a 0 2 + ∑ n = 1 ∞ ( a n − j b n 2 e j n ω 0 t + a n + j b n 2 e − j n ω 0 t ) = a 0 − j b 0 2 e j 0 ω 0 t + ∑ n = 1 ∞ ( a n − j b n 2 e j n ω 0 t + a − n − j b − n 2 e − j n ω 0 t ) = ∑ n = − ∞ ∞ a n − j b n 2 e j n ω 0 t = ∑ n = − ∞ ∞ F n e j n ω 0 t . (1.5) \begin{aligned} f\left( t \right) &\sim \frac{{{a_0}}}{2} + \sum\limits_{n = 1}^\infty {\left[ {{a_n}\cos \left( {n{\omega _0}t} \right) + {b_n}\sin \left( {n{\omega _0}t} \right)} \right]} \\ &= \frac{{{a_0}}}{2} + \sum\limits_{n = 1}^\infty {\left( {\frac{{{a_n} - j{b_n}}}{2}{e^{jn{\omega _0}t}} + \frac{{{a_n} + j{b_n}}}{2}{e^{ - jn{\omega _0}t}}} \right)} \\ &= \frac{{{a_0} - j{b_0}}}{2}{e^{j0{\omega _0}t}} + \sum\limits_{n = 1}^\infty {\left( {\frac{{{a_n} - j{b_n}}}{2}{e^{jn{\omega _0}t}} + \frac{{{a_{ - n}} - j{b_{ - n}}}}{2}{e^{ - jn{\omega _0}t}}} \right)} \\ &= \sum\limits_{n = - \infty }^\infty {\frac{{{a_n} - j{b_n}}}{2}{e^{jn{\omega _0}t}}} = \sum\limits_{n = - \infty }^\infty {{F_n}{e^{jn{\omega _0}t}}} . \\ \end{aligned} \tag{1.5} f(t)∼2a0+n=1∑∞[ancos(nω0t)+bnsin(nω0t)]=2a0+n=1∑∞(2an−jbnejnω0t+2an+jbne−jnω0t)=2a0−jb0ej0ω0t+n=1∑∞(2an−jbnejnω0t+2a−n−jb−ne−jnω0t)=n=−∞∑∞2an−jbnejnω0t=n=−∞∑∞Fnejnω0t.(1.5)

代入式 (1.1)，可得

F n = a n − j b n 2 = 1 T ∫ − T / 2 T / 2 f ( t ) [ cos ⁡ ( n ω 0 t ) − j sin ⁡ ( n ω 0 t ) ] d t = 1 T ∫ − T / 2 T / 2 f ( t ) ( e j n ω 0 t + e − j n ω 0 t 2 − j e j n ω 0 t − e − j n ω 0 t 2 j ) d t = 1 T ∫ − T / 2 T / 2 f ( t ) e − j n ω 0 t d t . (1.6) \begin{aligned} {F_n} &= \frac{{{a_n} - j{b_n}}}{2} \\ &= \frac{1}{T}\int\limits_{ - T/2}^{T/2} {f\left( t \right)\left[ {\cos \left( {n{\omega _0}t} \right) - j\sin \left( {n{\omega _0}t} \right)} \right]dt} \\ &= \frac{1}{T}\int\limits_{ - T/2}^{T/2} {f\left( t \right)\left( {\frac{{{e^{jn{\omega _0}t}} + {e^{ - jn{\omega _0}t}}}}{2} - j\frac{{{e^{jn{\omega _0}t}} - {e^{ - jn{\omega _0}t}}}}{{2j}}} \right)dt} \\ &= \frac{1}{T}\int\limits_{ - T/2}^{T/2} {f\left( t \right){e^{ - jn{\omega _0}t}}dt} . \\ \end{aligned} \tag{1.6} Fn=2an−jbn=T1−T/2∫T/2f(t)[cos(nω0t)−jsin(nω0t)]dt=T1−T/2∫T/2f(t)(2ejnω0t+e−jnω0t−j2jejnω0t−e−jnω0t)dt=T1−T/2∫T/2f(t)e−jnω0tdt.(1.6)

这样我们就得到了更加统一简洁的复数形式的傅里叶级数公式，在信号与系统中也被称为周期函数的傅里叶变换对。习惯上将求解傅里叶系数的过程称为正变换，而傅里叶级数本身则称为反变换，即

F [ n ] = F { f ( t ) } = 1 T ∫ − T / 2 T / 2 f ( t ) e − j n ω 0 t d t , f ( t ) = F − 1 { F [ n ] } = ∑ n = − ∞ ∞ F n e j n ω 0 t . (1.7) \begin{gathered} F\left[ n \right] = \mathcal{F}\left\{ {f\left( t \right)} \right\} = \frac{1}{T}\int\limits_{ - T/2}^{T/2} {f\left( t \right){e^{ - jn{\omega _0}t}}dt} , \\ f\left( t \right) = {\mathcal{F}^{ - 1}}\left\{ {F\left[ n \right]} \right\} = \sum\limits_{n = - \infty }^\infty {{F_n}{e^{jn{\omega _0}t}}} . \\ \end{gathered} \tag{1.7} F[n]=F{f(t)}=T1−T/2∫T/2f(t)e−jnω0tdt,f(t)=F−1{F[n]}=n=−∞∑∞Fnejnω0t.(1.7)

注意在高数中因为考虑数学的严谨性，傅里叶级数是用波浪线 ∼ { \sim } ∼ 符号来表示的，表示傅里叶级数不一定能收敛到 f ( t ) f\left( t \right) f(t)，但在信号与系统更多考量的是时域与频域能量相等，即帕塞瓦尔定理，且所涉及的函数的性质大多也比较好，为避免麻烦一般用等号来表示。在后面的讨论中，除了一些特殊的函数如狄拉克函数以及由狄拉克函数表示的级数以外，在忽略掉周期内有限个间断点后，不妨假设函数 f ( t ) f\left( t \right) f(t) 都是连续的。

尽管周期函数 f ( t ) f\left( t \right) f(t) 的自变量 t t t 本身是连续的，但所得的傅里叶系数 F [ n ] F\left[ n \right] F[n] 却只在离散的频率点 n ω 0 n{\omega _0} nω0 上有定义。这种性质就类似于电子能级跃迁的光谱，其是由若干离散频率的电磁波分量所组成的。因此，傅里叶系数 F [ n ] F\left[ n \right] F[n] 常被称为周期信号 f ( t ) f\left( t \right) f(t) 的频谱， F [ n ] F\left[ n \right] F[n] 所在的空间也就被称为频域。另外，因为 e j n ω 0 t {e^{jn{\omega _0}t}} ejnω0t 在复平面上对应着一个单位圆，而不是实数域中的无限数轴，这使得复数域中的信号分析与处理相比于实数域会更加的简便，这种性质在后来的拉普拉斯变换与 z 变换中尤为重要。

2 非周期函数的傅里叶变换

如果考虑非周期函数 f ( t ) f\left( t \right) f(t)，我们可认为其周期 T → ∞ {T \to \infty } T→∞，于是有基频 ω 0 → 0 {{\omega _0} \to 0} ω0→0。引入变量 ω \omega ω，且有 ω n = n ω 0 {{\omega _n} = n{\omega _0}} ωn=nω0，那么

Δ ω n = ω n − ω n − 1 = ω 0 → 0. \Delta {\omega _n} = {\omega _n} - {\omega _{n - 1}} = {\omega _0} \to 0. Δωn=ωn−ωn−1=ω0→0.

于是有 Δ ω n → d ω {{\Delta {\omega _n} \to d\omega }} Δωn→dω，此时离散的 ω n {\omega _n} ωn 成为连续的 ω {\omega} ω 而无需再区分。这时

lim ⁡ T → ∞ f ( t ) = lim ⁡ T → ∞ ∑ n = − ∞ ∞ F n Δ ω n Δ ω n e j ω n t = lim ⁡ T → ∞ ∑ n = − ∞ ∞ ( 1 T ⋅ Δ ω n ∫ − T / 2 T / 2 f ( τ ) e − j ω n τ d τ ) e j ω n t Δ ω n = lim ⁡ T → ∞ ∑ n = − ∞ ∞ ( 1 2 π ∫ − T / 2 T / 2 f ( τ ) e − j ω n τ d τ ) e j ω n t Δ ω n = 1 2 π ∫ − ∞ ∞ ( ∫ − ∞ ∞ f ( τ ) e − j ω τ d τ ) e j ω t d ω . (2.1) \begin{aligned} \mathop {\lim }\limits_{T \to \infty } f\left( t \right) &= \mathop {\lim }\limits_{T \to \infty } \sum\limits_{n = - \infty }^\infty {\frac{{{F_n}}}{{\Delta {\omega _n}}}\Delta {\omega _n}{e^{j{\omega _n}t}}} \\ &= \mathop {\lim }\limits_{T \to \infty } \sum\limits_{n = - \infty }^\infty {\left( {\frac{1}{{T \cdot \Delta {\omega _n}}}\int\limits_{ - T/2}^{T/2} {f\left( \tau \right){e^{ - j{\omega _n}\tau }}d\tau } } \right){e^{j{\omega _n}t}}\Delta {\omega _n}} \\ &= \mathop {\lim }\limits_{T \to \infty } \sum\limits_{n = - \infty }^\infty {\left( {\frac{1}{{2\pi }}\int\limits_{ - T/2}^{T/2} {f\left( \tau \right){e^{ - j{\omega _n}\tau }}d\tau } } \right){e^{j{\omega _n}t}}\Delta {\omega _n}} \\ &= \frac{1}{{2\pi }}\int\limits_{ - \infty }^\infty {\left( {\int\limits_{ - \infty }^\infty {f\left( \tau \right){e^{ - j\omega \tau }}d\tau } } \right){e^{j\omega t}}d\omega } . \\ \end{aligned} \tag{2.1} T→∞limf(t)=T→∞limn=−∞∑∞ΔωnFnΔωnejωnt=T→∞limn=−∞∑∞ T⋅Δωn1−T/2∫T/2f(τ)e−jωnτdτ ejωntΔωn=T→∞limn=−∞∑∞ 2π1−T/2∫T/2f(τ)e−jωnτdτ ejωntΔωn=2π1−∞∫∞ −∞∫∞f(τ)e−jωτdτ ejωtdω.(2.1)

于是可以定义，设非周期函数 f ( t ) {f\left( t \right)} f(t) 在 ( − ∞ , + ∞ ) {\left( { - \infty , + \infty } \right)} (−∞,+∞) 上绝对可积，则 f ( t ) f\left( t \right) f(t) 的傅里叶变换 (Fourier Transform, FT) 与反变换分别为

F ( j ω ) = F { f ( t ) } = ∫ − ∞ ∞ f ( t ) e − j ω t d t , f ( t ) = F − 1 { F ( j ω ) } = 1 2 π ∫ − ∞ ∞ F ( j ω ) e j ω t d ω . (2.2) \begin{gathered} F\left( {j\omega } \right) = \mathcal{F}\left\{ {f\left( t \right)} \right\} = \int\limits_{ - \infty }^\infty {f\left( t \right){e^{ - j\omega t}}dt} , \\ f\left( t \right) = {\mathcal{F}^{ - 1}}\left\{ {F\left( {j\omega } \right)} \right\} = \frac{1}{{2\pi }}\int\limits_{ - \infty }^\infty {F\left( {j\omega } \right){e^{j\omega t}}d\omega } . \\ \end{gathered} \tag{2.2} F(jω)=F{f(t)}=−∞∫∞f(t)e−jωtdt,f(t)=F−1{F(jω)}=2π1−∞∫∞F(jω)ejωtdω.(2.2)

为了与后续的离散时间傅里叶变换做区分，有时还会特别称作连续时间傅里叶变换 (Continuous Time FT, CTFT)。注意以上的推导是不严格的，这里着重于其形式而略去相关证明。另外 F ( j ω ) F\left( {j\omega } \right) F(jω) 也可直接写作 F ( ω ) F\left( \omega \right) F(ω)，加上虚数符号主要是为了表达傅里叶变换的结果为复数，同时也与后面的离散时间傅里叶变换做区分。

相比傅里叶级数所得离散系数 F [ n ] F\left[ n \right] F[n]，傅里叶变换 F ( j ω ) F\left( {j\omega } \right) F(jω) 定义在连续的频率自变量 ω \omega ω上，且 ω \omega ω与 n n n 的量纲也有所区别， ω \omega ω 是由 lim ⁡ ω 0 → 0 n ω 0 {\mathop {\lim }\limits_{{\omega _0} \to 0} n{\omega _0}} ω0→0limnω0 这个极限来定义的。考虑定义在有限区间 [ − T / 2 , T / 2 ] \left[ { - T/2,\;T/2} \right] [−T/2,T/2] 的函数 f ( t ) f\left( t \right) f(t)，对其以 T T T 为周期进行延拓获得周期函数 f ~ ( t ) \tilde f\left( t \right) f~(t)，那么可得 f ( t ) f\left( t \right) f(t) 的傅里叶变换谱 F ( j ω ) F\left( {j\omega } \right) F(jω) 与 f ~ ( t ) \tilde f\left( t \right) f~(t) 的傅里叶级数系数 F ~ [ n ] \tilde F\left[ n \right] F~[n] 分别为

F ( j ω ) = ∫ − T / 2 T / 2 f ( t ) e − j ω t d t , F ~ [ n ] = 1 T ∫ − T / 2 T / 2 f ~ ( t ) e − j n Ω T t d t = 1 T F ( j ω ) ∣ ω = n Ω T , w h e r e Ω T = 2 π / T . (2.3) \begin{gathered} F\left( {j\omega } \right) = \int\limits_{ - T/2}^{T/2} {f\left( t \right){e^{ - j\omega t}}dt} , \\ \tilde F\left[ n \right] = \frac{1}{T}\int\limits_{ - T/2}^{T/2} {\tilde f\left( t \right){e^{ - jn{\Omega _T}t}}dt} = \frac{1}{T}{\left. {F\left( {j\omega } \right)} \right|_{\omega = n{\Omega _T}}}, \\ where\;{\Omega _T} = 2\pi /T. \\ \end{gathered} \tag{2.3} F(jω)=−T/2∫T/2f(t)e−jωtdt,F~[n]=T1−T/2∫T/2f~(t)e−jnΩTtdt=T1F(jω)∣ω=nΩT,whereΩT=2π/T.(2.3)

也就是说，周期函数的傅里叶系数 F ~ [ n ] \tilde F\left[ n \right] F~[n] 是对该周期函数在一个周期内进行傅里叶变换所得 F ( j ω ) F\left( {j\omega } \right) F(jω) 以 Ω T {\Omega _T} ΩT 为间隔均匀采样以及以 T T T 进行幅度缩放的结果。这里涉及到了采样的概念，从而引出了下面我们所要说的离散时间傅里叶变换。我们也将看到，时域与频域的性质总是镜像对称的。

3 离散时间傅里叶变换

在讨论采样与离散时间傅里叶变换之前，我们需要先介绍狄拉克函数，其通常是基于极限定义的，即

δ ( t ) = lim ⁡ ε → 0 { 1 / ε , t ∈ ( − ε / 2 , ε / 2 ) 0 , o t h e r s (3.1) \delta \left( t \right) = \mathop {\lim }\limits_{\varepsilon \to 0} \left\{ {\begin{array}{c} {1/\varepsilon ,}&{t \in \left( { - \varepsilon /2,\varepsilon /2} \right)} \\ {0,}&{others} \end{array}} \right.\tag{3.1} δ(t)=ε→0lim{1/ε,0,t∈(−ε/2,ε/2)others(3.1)

注意狄拉克函数的定义不是唯一的，以上采用了比较简单的矩形函数定义，也可使用其他如三角形函数的定义，等等。其只在 t → 0 t \to 0 t→0 有无穷大的值，而在其他位置为 0，因此狄拉克函数也被称为冲激函数，并且有

∫ − ∞ ∞ δ ( t ) d t = lim ⁡ ε → 0 ∫ − ε / 2 ε / 2 1 ε d t = 1. (3.2) \int\limits_{ - \infty }^\infty {\delta \left( t \right)dt} = \mathop {\lim }\limits_{\varepsilon \to 0} \int\limits_{ - \varepsilon /2}^{\varepsilon /2} {\frac{1}{\varepsilon }dt} = 1.\tag{3.2} −∞∫∞δ(t)dt=ε→0lim−ε/2∫ε/2ε1dt=1.(3.2)

假设函数 f ( t ) f\left( t \right) f(t) 连续，因为 lim ⁡ t → t 0 f ( t ) = f ( t 0 ) \mathop {\lim }\limits_{t \to {t_0}} f\left( t \right) = f\left( {{t_0}} \right) t→t0limf(t)=f(t0)，容易推导得其以下性质

∫ − ∞ ∞ f ( t ) δ ( t − t 0 ) d t = lim ⁡ ε → 0 ∫ t 0 − ε / 2 t 0 + ε / 2 f ( t ) 1 ε d t = f ( t 0 ) lim ⁡ ε → 0 ∫ t 0 − ε / 2 t 0 + ε / 2 1 ε d t = f ( t 0 ) = ∫ − ∞ ∞ f ( t 0 ) δ ( t − t 0 ) d t . (3.3) \begin{aligned} \int\limits_{ - \infty }^\infty {f\left( t \right)\delta \left( {t - {t_0}} \right)dt} &= \mathop {\lim }\limits_{\varepsilon \to 0} \int\limits_{{t_0} - \varepsilon /2}^{{t_0} + \varepsilon /2} {f\left( t \right)\frac{1}{\varepsilon }dt} \\ &= f\left( {{t_0}} \right)\mathop {\lim }\limits_{\varepsilon \to 0} \int\limits_{{t_0} - \varepsilon /2}^{{t_0} + \varepsilon /2} {\frac{1}{\varepsilon }dt} = f\left( {{t_0}} \right) \\ &= \int\limits_{ - \infty }^\infty {f\left( {{t_0}} \right)\delta \left( {t - {t_0}} \right)dt} . \\ \end{aligned} \tag{3.3} −∞∫∞f(t)δ(t−t0)dt=ε→0limt0−ε/2∫t0+ε/2f(t)ε1dt=f(t0)ε→0limt0−ε/2∫t0+ε/2ε1dt=f(t0)=−∞∫∞f(t0)δ(t−t0)dt.(3.3)

这个性质也称为狄拉克函数的采样性质。除此以外，狄拉克函数的卷积性质也十分常用，即

f ( t ) ∗ δ ( t − t 0 ) = ∫ − ∞ ∞ f ( τ ) δ ( t − t 0 − τ ) d τ = lim ⁡ ε → 0 ∫ t − t 0 − ε / 2 t − t 0 + ε / 2 f ( τ ) 1 ε d τ = f ( t − t 0 ) lim ⁡ ε → 0 ∫ t − t 0 − ε / 2 t − t 0 + ε / 2 1 ε d τ = f ( t − t 0 ) . (3.4) \begin{aligned} f\left( t \right) * \delta \left( {t - {t_0}} \right) &= \int\limits_{ - \infty }^\infty {f\left( \tau \right)\delta \left( {t - {t_0} - \tau } \right)d\tau } \\ &= \mathop {\lim }\limits_{\varepsilon \to 0} \int\limits_{t - {t_0} - \varepsilon /2}^{t - {t_0} + \varepsilon /2} {f\left( \tau \right)\frac{1}{\varepsilon }d\tau } \\ &= f\left( {t - {t_0}} \right)\mathop {\lim }\limits_{\varepsilon \to 0} \int\limits_{t - {t_0} - \varepsilon /2}^{t - {t_0} + \varepsilon /2} {\frac{1}{\varepsilon }d\tau } \\ &= f\left( {t - {t_0}} \right). \\ \end{aligned} \tag{3.4} f(t)∗δ(t−t0)=−∞∫∞f(τ)δ(t−t0−τ)dτ=ε→0limt−t0−ε/2∫t−t0+ε/2f(τ)ε1dτ=f(t−t0)ε→0limt−t0−ε/2∫t−t0+ε/2ε1dτ=f(t−t0).(3.4)

也就是说，狄拉克函数可以使得时间轴 t t t 的原点根据冲激产生的位置进行移动，从而可以以卷积的形式来表达函数的时移。

狄拉克函数提供了从连续时间信号到离散时间信号的分析方法。现实采样通常涉及采样周期 T T T 与采样时间 ε \varepsilon ε 两个概念。采样周期为两个相邻离散采样点的时间距离，一般我们只讨论等间隔的采样。采样时间主要受采样技术的限制，现有技术通常无法在一瞬间测量到信号比较准确的值。例如测量电压时通常需要对电容进行充电，而充电需要一定时间，这个过程可以用积分来表示，即

f ˉ ( n T ) = 1 ε ∫ n T − ε / 2 n T + ε / 2 f ( t ) d t . (3.5) \bar f\left( {nT} \right) = \frac{1}{\varepsilon }\int\limits_{nT - \varepsilon /2}^{nT + \varepsilon /2} {f\left( t \right)dt} .\tag{3.5} fˉ(nT)=ε1nT−ε/2∫nT+ε/2f(t)dt.(3.5)

可以发现狄拉克函数相当于一个理想的瞬时采样过程，使用狄拉克函数可以忽略采样时间这个概念，从而在信号分析中免去很多不必要的运算。

那么，假设有一连续时间信号 f a ( t ) {f_a}\left( t \right) fa(t)，以 T T T 为采样周期，对 f a ( t ) {f_a}\left( t \right) fa(t) 在 t = n T t = nT t=nT 时均匀采样，可得离散序列 f [ n ] f\left[ n \right] f[n]，其中

f [ n ] = f a ( n T ) , − ∞ < n < ∞ (3.6) f\left[ n \right] = {f_a}\left( {nT} \right),\quad - \infty < n < \infty \tag{3.6} f[n]=fa(nT),−∞<n<∞(3.6)

采样的过程可表示为周期狄拉克函数的加权叠加，从而获得连续时间信号 f p ( t ) {f_p}\left( t \right) fp(t)，即

f p ( t ) = ∑ n = − ∞ ∞ f a ( n T ) δ ( t − n T ) . (3.7) {f_p}\left( t \right) = \sum\limits_{n = - \infty }^\infty {{f_a}\left( {nT} \right)\delta \left( {t - nT} \right)} .\tag{3.7} fp(t)=n=−∞∑∞fa(nT)δ(t−nT).(3.7)

当然， f p ( t ) {f_p}\left( t \right) fp(t) 也等价于 f a ( t ) {f_a}\left( t \right) fa(t) 与周期冲激串 p ( t ) p\left( t \right) p(t) 的乘积，即

f p ( t ) = f a ( t ) p ( t ) , w h e r e p ( t ) = ∑ n = − ∞ ∞ δ ( t − n T ) . (3.8) \begin{gathered} {f_p}\left( t \right) = {f_a}\left( t \right)p\left( t \right), \\ where\;p\left( t \right) = \sum\limits_{n = - \infty }^\infty {\delta \left( {t - nT} \right)} . \\ \end{gathered} \tag{3.8} fp(t)=fa(t)p(t),wherep(t)=n=−∞∑∞δ(t−nT).(3.8)

注意连续时间信号与连续信号的区别，连续时间只是说自变量 t t t 是连续的，并没有说信号本身是连续的。另外，由于 f p ( t ) {f_p}\left( t \right) fp(t) 是由狄拉克函数的加权叠加来表示的，其只有在区间积分下才具有意义，例如

∫ n T − τ n T + τ f p ( t ) d t = f a ( n T ) = f [ n ] , τ < T / 2. (3.9) \int\limits_{nT - \tau }^{nT + \tau } {{f_p}\left( t \right)dt} = {f_a}\left( {nT} \right) = f[n],\;\tau < T/2.\tag{3.9} nT−τ∫nT+τfp(t)dt=fa(nT)=f[n],τ<T/2.(3.9)

那么，可推导得 f p ( t ) {f_p}\left( t \right) fp(t) 的傅里叶变换为

F p ( j ω ) = ∫ − ∞ ∞ f p ( t ) e − j ω t d t = ∫ − ∞ ∞ ( ∑ n = − ∞ ∞ f a ( n T ) δ ( t − n T ) ) e − j ω t d t = ∑ n = − ∞ ∞ f a ( n T ) ∫ − ∞ ∞ δ ( t − n T ) e − j ω t d t = ∑ n = − ∞ ∞ f a ( n T ) ∫ − ∞ ∞ δ ( t − n T ) e − j ω ( t − n T ) − j ω n T d ( t − n T ) = ∑ n = − ∞ ∞ f a ( n T ) e − j ω n T . (3.10) \begin{aligned} {F_p}\left( {j\omega } \right) &= \int\limits_{ - \infty }^\infty {{f_p}\left( t \right){e^{ - j\omega t}}dt} \\ &= \int\limits_{ - \infty }^\infty {\left( {\sum\limits_{n = - \infty }^\infty {{f_a}\left( {nT} \right)\delta \left( {t - nT} \right)} } \right){e^{ - j\omega t}}dt} \\ &= \sum\limits_{n = - \infty }^\infty {{f_a}\left( {nT} \right)\int\limits_{ - \infty }^\infty {\delta \left( {t - nT} \right){e^{ - j\omega t}}dt} } \\ &= \sum\limits_{n = - \infty }^\infty {{f_a}\left( {nT} \right)\int\limits_{ - \infty }^\infty {\delta \left( {t - nT} \right){e^{ - j\omega \left( {t - nT} \right) - j\omega nT}}d\left( {t - nT} \right)} } \\ &= \sum\limits_{n = - \infty }^\infty {{f_a}\left( {nT} \right){e^{ - j\omega nT}}} . \\ \end{aligned} \tag{3.10} Fp(jω)=−∞∫∞fp(t)e−jωtdt=−∞∫∞(n=−∞∑∞fa(nT)δ(t−nT))e−jωtdt=n=−∞∑∞fa(nT)−∞∫∞δ(t−nT)e−jωtdt=n=−∞∑∞fa(nT)−∞∫∞δ(t−nT)e−jω(t−nT)−jωnTd(t−nT)=n=−∞∑∞fa(nT)e−jωnT.(3.10)

因为我们更关心信号本身，所以当连续时间信号经过采样获得离散时间信号后，可以忽略采样周期 T T T，从而简化后续的运算。为了区别于一般的连续时间信号的傅里叶变换，我们用 Ω \Omega Ω 代替式 (3.10) 中的 ω \omega ω，即

F p ( j Ω ) = ∑ n = − ∞ ∞ f a ( n T ) e − j Ω n T . (3.11) {F_p}\left( {j\Omega } \right) = \sum\limits_{n = - \infty }^\infty {{f_a}\left( {nT} \right){e^{ - j\Omega nT}}} .\tag{3.11} Fp(jΩ)=n=−∞∑∞fa(nT)e−jΩnT.(3.11)

那么可定义序列 f [ n ] f\left[ n \right] f[n] 的离散时间傅里叶变换 (Discrete Time FT, DTFT) 为

F ( e j ω ) = ∑ n = − ∞ ∞ f [ n ] e − j ω n , w h e r e ω = Ω T . (3.12) F\left( {{e^{j\omega }}} \right) = \sum\limits_{n = - \infty }^\infty {f\left[ n \right]{e^{ - j\omega n}}} ,\;where\;\omega = \Omega T.\tag{3.12} F(ejω)=n=−∞∑∞f[n]e−jωn,whereω=ΩT.(3.12)

并且有

ω = 2 π ⇔ Ω = ω T = 2 π T = Ω T . (3.13) \omega = 2\pi \Leftrightarrow \Omega = \frac{\omega }{T} = \frac{{2\pi }}{T} = {\Omega _T}.\tag{3.13} ω=2π⇔Ω=Tω=T2π=ΩT.(3.13)

这样我们就有了分析离散时间信号的方法，这对于只能处理离散数据的数字电路尤为重要。然而，尽管我们在式 (3.12) 中有意忽略了采样周期 T T T，但这并不代表不需要对其进行进一步的讨论。因为采样是从连续时间到离散时间的转换过程，直觉上这会导致一部分信息丢失，那么该如何去证明以及衡量这种信息失真是我们接下来要讨论的重点。

可以发现，不同于一般连续时间信号的傅里叶变换，由周期狄拉克函数加权叠加所定义的连续时间信号 f p ( t ) {f_p}\left( t \right) fp(t) 的傅里叶变换 F p ( j Ω ) {F_p}\left( {j\Omega } \right) Fp(jΩ)，以及序列 f [ n ] f\left[ n \right] f[n] 的离散时间傅里叶变换 F ( e j ω ) F\left( {{e^{j\omega }}} \right) F(ejω) 都是具有明确周期性的，即

F p ( j Ω ) = ∑ n = − ∞ ∞ f a ( n T ) e − j Ω n T = ∑ n = − ∞ ∞ f a ( n T ) e − j Ω n 2 π Ω T = ∑ n = − ∞ ∞ f a ( n T ) e − j ( Ω + k Ω T ) n 2 π Ω T = F p ( j ( Ω + k Ω T ) ) , F ( e j ω ) = ∑ n = − ∞ ∞ f [ n ] e − j ω n = ∑ n = − ∞ ∞ f [ n ] e − j ( ω + 2 k π ) n = F ( e j ( ω + 2 k π ) ) = F ( e j ( ω + k Ω T T ) ) , w h e r e ω = Ω T , Ω T = 2 π / T . (3.14) \begin{aligned} {F_p}\left( {j\Omega } \right) &= \sum\limits_{n = - \infty }^\infty {{f_a}\left( {nT} \right){e^{ - j\Omega nT}}} = \sum\limits_{n = - \infty }^\infty {{f_a}\left( {nT} \right){e^{ - j\Omega n\tfrac{{2\pi }}{{{\Omega _T}}}}}} \\ &= \sum\limits_{n = - \infty }^\infty {{f_a}\left( {nT} \right){e^{ - j\left( {\Omega + k{\Omega _T}} \right)n\tfrac{{2\pi }}{{{\Omega _T}}}}}} = {F_p}\left( {j\left( {\Omega + k{\Omega _T}} \right)} \right), \\ F\left( {{e^{j\omega }}} \right) &= \sum\limits_{n = - \infty }^\infty {f\left[ n \right]{e^{ - j\omega n}}} = \sum\limits_{n = - \infty }^\infty {f\left[ n \right]{e^{ - j\left( {\omega + 2k\pi } \right)n}}} \\ &= F\left( {{e^{j\left( {\omega + 2k\pi } \right)}}} \right) = F\left( {{e^{j\left( {\omega + k{\Omega _T}T} \right)}}} \right), \\ where\;\omega &= \Omega T,\;{\Omega _T} = 2\pi /T. \\ \end{aligned} \tag{3.14} Fp(jΩ)F(ejω)whereω=n=−∞∑∞fa(nT)e−jΩnT=n=−∞∑∞fa(nT)e−jΩnΩT2π=n=−∞∑∞fa(nT)e−j(Ω+kΩT)nΩT2π=Fp(j(Ω+kΩT)),=n=−∞∑∞f[n]e−jωn=n=−∞∑∞f[n]e−j(ω+2kπ)n=F(ej(ω+2kπ))=F(ej(ω+kΩTT)),=ΩT,ΩT=2π/T.(3.14)

即 F p ( j Ω ) {F_p}\left( {j\Omega } \right) Fp(jΩ) 是以采样频率 Ω T {\Omega _T} ΩT 为周期的函数， F ( e j ω ) F\left( {{e^{j\omega }}} \right) F(ejω) 则相应地以 2 π 2\pi 2π 为周期，这种周期性是由采样过程中所使用的周期冲激串 p ( t ) p\left( t \right) p(t) 所引入的。除此以外，尽管我们通过 F p ( j Ω ) {F_p}\left( {j\Omega } \right) Fp(jΩ) 与 F ( e j ω ) F\left( {{e^{j\omega }}} \right) F(ejω) 定义了离散时间傅里叶变换的形式，但并没有将其与一般连续时间信号的傅里叶变换 F a ( j Ω ) {F_a}\left( {j\Omega } \right) Fa(jΩ) 建立起联系。由于连续时间函数 f a ( t ) {f_a}\left( t \right) fa(t) 的傅里叶变换 F a ( j Ω ) {F_a}\left( {j\Omega } \right) Fa(jΩ) 本身并没有周期性，即 F a ( j Ω ) {F_a}\left( {j\Omega } \right) Fa(jΩ) 有可能在整个数轴 Ω ∈ ( − ∞ , + ∞ ) \Omega \in \left( { - \infty , + \infty } \right) Ω∈(−∞,+∞) 上非周期性地延伸，那么如果 F p ( j Ω ) {F_p}\left( {j\Omega } \right) Fp(jΩ) 是由 F a ( j Ω ) {F_a}\left( {j\Omega } \right) Fa(jΩ) 推导而来的， F p ( j Ω ) {F_p}\left( {j\Omega } \right) Fp(jΩ) 的周期性将不可避免地导致 F a ( j Ω ) {F_a}\left( {j\Omega } \right) Fa(jΩ) 的部分信息被破坏。

为了证明这种信息失真，我们需要引入泊松求和公式，即对于一非周期连续时间函数 ϕ ( t ) \phi \left( t \right) ϕ(t)，有

∑ n = − ∞ ∞ ϕ ( t + n T ) = 1 T ∑ k = − ∞ ∞ Φ ( j k Ω T ) e j k Ω T t , w h e r e Ω T = 2 π / T . (3.15) \begin{gathered} \sum\limits_{n = - \infty }^\infty {\phi \left( {t + nT} \right)} = \frac{1}{T}\sum\limits_{k = - \infty }^\infty {\Phi \left( {jk{\Omega _T}} \right){e^{jk{\Omega _T}t}}} , \\ where\;{\Omega _T} = 2\pi /T. \\ \end{gathered} \tag{3.15} n=−∞∑∞ϕ(t+nT)=T1k=−∞∑∞Φ(jkΩT)ejkΩTt,whereΩT=2π/T.(3.15)

其中 Φ ( j Ω ) \Phi \left( {j\Omega } \right) Φ(jΩ) 为 ϕ ( t ) \phi \left( t \right) ϕ(t) 的连续时间傅里叶变换， T T T 为任意非零长度。要证明这一点，需要综合地使用到我们前面所述的周期函数的傅里叶级数以及连续时间函数的傅里叶变换的知识，以及关于狄拉克函数的性质，即

∑ n = − ∞ ∞ ϕ ( t + n T ) = ∑ n = − ∞ ∞ ϕ ( t ) ∗ δ ( t + n T ) = ϕ ( t ) ∗ ∑ n = − ∞ ∞ δ ( t + n T ) = ϕ ( t ) ∗ [ ∑ k = − ∞ − ∞ ( 1 T ∫ − T / 2 T / 2 δ ( τ ) e − j k Ω T τ d τ ) e j k Ω T t ] = ϕ ( t ) ∗ 1 T ∑ k = − ∞ ∞ e j k Ω T t = 1 T ∑ k = − ∞ ∞ ϕ ( t ) ∗ e j k Ω T t = 1 T ∑ k = − ∞ ∞ ( ∫ − ∞ ∞ ϕ ( τ ) e j k Ω T ( t − τ ) d τ ) = 1 T ∑ k = − ∞ ∞ ( ∫ − ∞ ∞ ϕ ( τ ) e − j k Ω T τ d τ ) e j k Ω T t = 1 T ∑ k = − ∞ ∞ Φ ( j k Ω T ) e j k Ω T t . (3.16) \begin{aligned} \sum\limits_{n = - \infty }^\infty {\phi \left( {t + nT} \right)} &= \sum\limits_{n = - \infty }^\infty {\phi \left( t \right) * \delta \left( {t + nT} \right)} \\ &= \phi \left( t \right) * \sum\limits_{n = - \infty }^\infty {\delta \left( {t + nT} \right)} \\ &= \phi \left( t \right) * \left[ {\sum\limits_{k = - \infty }^{ - \infty } {\left( {\frac{1}{T}\int\limits_{ - T/2}^{T/2} {\delta \left( \tau \right){e^{ - jk{\Omega _T}\tau }}d\tau } } \right){e^{jk{\Omega _T}t}}} } \right] \\ &= \phi \left( t \right) * \frac{1}{T}\sum\limits_{k = - \infty }^\infty {{e^{jk{\Omega _T}t}}} \\ &= \frac{1}{T}\sum\limits_{k = - \infty }^\infty {\phi \left( t \right) * {e^{jk{\Omega _T}t}}} \\ &= \frac{1}{T}\sum\limits_{k = - \infty }^\infty {\left( {\int\limits_{ - \infty }^\infty {\phi \left( \tau \right){e^{jk{\Omega _T}\left( {t - \tau } \right)}}d\tau } } \right)} \\ &= \frac{1}{T}\sum\limits_{k = - \infty }^\infty {\left( {\int\limits_{ - \infty }^\infty {\phi \left( \tau \right){e^{ - jk{\Omega _T}\tau }}d\tau } } \right){e^{jk{\Omega _T}t}}} \\ &= \frac{1}{T}\sum\limits_{k = - \infty }^\infty {\Phi \left( {jk{\Omega _T}} \right){e^{jk{\Omega _T}t}}} . \\ \end{aligned} \tag{3.16} n=−∞∑∞ϕ(t+nT)=n=−∞∑∞ϕ(t)∗δ(t+nT)=ϕ(t)∗n=−∞∑∞δ(t+nT)=ϕ(t)∗ k=−∞∑−∞ T1−T/2∫T/2δ(τ)e−jkΩTτdτ ejkΩTt =ϕ(t)∗T1k=−∞∑∞ejkΩTt=T1k=−∞∑∞ϕ(t)∗ejkΩTt=T1k=−∞∑∞ −∞∫∞ϕ(τ)ejkΩT(t−τ)dτ =T1k=−∞∑∞ −∞∫∞ϕ(τ)e−jkΩTτdτ ejkΩTt=T1k=−∞∑∞Φ(jkΩT)ejkΩTt.(3.16)

当 t = 0 t = 0 t=0 时，有

∑ n = − ∞ ∞ ϕ ( n T ) = 1 T ∑ k = − ∞ ∞ Φ ( j k Ω T ) . (3.17) \sum\limits_{n = - \infty }^\infty {\phi \left( {nT} \right)} = \frac{1}{T}\sum\limits_{k = - \infty }^\infty {\Phi \left( {jk{\Omega _T}} \right)} .\tag{3.17} n=−∞∑∞ϕ(nT)=T1k=−∞∑∞Φ(jkΩT).(3.17)

根据 CTFT 频移性质，即

1 2 π ∫ − ∞ ∞ Φ ( Ω + Ψ ) e j Ω t d Ω = 1 2 π ∫ − ∞ ∞ Φ ( Ω + Ψ ) e j ( Ω + Ψ ) t − j Ψ t d ( Ω + Ψ ) = e − j Ψ t 1 2 π ∫ − ∞ ∞ Φ ( Ω ) e j Ω t d Ω = ϕ ( t ) e − j Ψ t . (3.18) \begin{aligned} &\frac{1}{{2\pi }}\int\limits_{ - \infty }^\infty {\Phi \left( {\Omega + \Psi } \right){e^{j\Omega t}}d\Omega } \\ &= \frac{1}{{2\pi }}\int\limits_{ - \infty }^\infty {\Phi \left( {\Omega + \Psi } \right){e^{j\left( {\Omega + \Psi } \right)t - j\Psi t}}d\left( {\Omega + \Psi } \right)} \\ &= {e^{ - j\Psi t}}\frac{1}{{2\pi }}\int\limits_{ - \infty }^\infty {\Phi \left( \Omega \right){e^{j\Omega t}}d\Omega } \\ &= \phi \left( t \right){e^{ - j\Psi t}}. \\ \end{aligned} \tag{3.18} 2π1−∞∫∞Φ(Ω+Ψ)ejΩtdΩ=2π1−∞∫∞Φ(Ω+Ψ)ej(Ω+Ψ)t−jΨtd(Ω+Ψ)=e−jΨt2π1−∞∫∞Φ(Ω)ejΩtdΩ=ϕ(t)e−jΨt.(3.18)

于是有

∑ n = − ∞ ∞ ϕ ( n T ) e − j Ψ n T = 1 T ∑ k = − ∞ ∞ Φ ( j ( k Ω T + Ψ ) ) . (3.19) \sum\limits_{n = - \infty }^\infty {\phi \left( {nT} \right){e^{ - j\Psi nT}}} = \frac{1}{T}\sum\limits_{k = - \infty }^\infty {\Phi \left( {j\left( {k{\Omega _T} + \Psi } \right)} \right)} .\tag{3.19} n=−∞∑∞ϕ(nT)e−jΨnT=T1k=−∞∑∞Φ(j(kΩT+Ψ)).(3.19)

将 f a ( t ) {f_a}\left( t \right) fa(t) 替代 ϕ ( t ) \phi \left( t \right) ϕ(t)， k Ω T k{\Omega _T} kΩT 替代 Ψ \Psi Ψ，即有

∑ n = − ∞ ∞ f a ( n T ) e − j Ω n T = F p ( j Ω ) = 1 T ∑ k = − ∞ ∞ F a ( j ( Ω + k Ω T ) ) . (3.20) \sum\limits_{n = - \infty }^\infty {{f_a}\left( {nT} \right){e^{ - j\Omega nT}}} = {F_p}\left( {j\Omega } \right) = \frac{1}{T}\sum\limits_{k = - \infty }^\infty {{F_a}\left( {j\left( {\Omega + k{\Omega _T}} \right)} \right)} .\tag{3.20} n=−∞∑∞fa(nT)e−jΩnT=Fp(jΩ)=T1k=−∞∑∞Fa(j(Ω+kΩT)).(3.20)

因此可知，函数 f p ( t ) {f_p}\left( t \right) fp(t) 的连续时间傅里叶变换是 f a ( t ) {f_a}\left( t \right) fa(t) 的连续时间傅里叶变换经过以采样频率 Ω T {\Omega _T} ΩT 为周期进行平移以及幅度缩放后相加的结果。那么可预知的是，如果 F a ( j Ω ) {F_a}\left( {j\Omega } \right) Fa(jΩ) 在 ∣ Ω ∣ > Ω T / 2 \left| \Omega \right| > {\Omega _T}/2 ∣Ω∣>ΩT/2 时有非零值，就极有可能出现与 F a ( j ( Ω + k Ω T ) ) {F_a}\left( {j\left( {\Omega + k{\Omega _T}} \right)} \right) Fa(j(Ω+kΩT)) 重叠的情况，而重叠部分的频谱将无法再通过滤波进行分离，即我们无法从 F p ( j Ω ) {F_p}\left( {j\Omega } \right) Fp(jΩ) 完整地恢复 F a ( j Ω ) {F_a}\left( {j\Omega } \right) Fa(jΩ)，从而造成信息丢失。但反过来说，如果 F a ( j Ω ) {F_a}\left( {j\Omega } \right) Fa(jΩ) 只在 ∣ Ω ∣ < Ω T / 2 \left| \Omega \right| < {\Omega _T}/2 ∣Ω∣<ΩT/2 内有非零值，就不会与 F a ( j ( Ω + k Ω T ) ) {F_a}\left( {j\left( {\Omega + k{\Omega _T}} \right)} \right) Fa(j(Ω+kΩT)) 存在重叠的部分，在只知道 F p ( j Ω ) {F_p}\left( {j\Omega } \right) Fp(jΩ) 的情况下，我们可以通过理想低通滤波器取 F p ( j Ω ) {F_p}\left( {j\Omega } \right) Fp(jΩ) 在 Ω ∈ ( Ω T / 2 , Ω T / 2 ) \Omega \in \left( {{\Omega _T}/2,{\Omega _T}/2} \right) Ω∈(ΩT/2,ΩT/2) 的部分完整地恢复 F a ( j Ω ) {F_a}\left( {j\Omega } \right) Fa(jΩ)，由于傅里叶变换是可逆的，这意味着我们通过离散序列 f [ n ] f\left[ n \right] f[n] 无失真地恢复连续时间信号 f a ( t ) {f_a}\left( t \right) fa(t)。因为序列 f [ n ] f\left[ n \right] f[n] 的离散时间傅里叶变换 F ( e j ω ) F\left( {{e^{j\omega }}} \right) F(ejω) 只是 F p ( j Ω ) {F_p}\left( {j\Omega } \right) Fp(jΩ) 将 Ω T {\Omega _T} ΩT 缩放到 2 π 2\pi 2π 的结果，其性质只需参照以上推论即可。

因此可以发现，虽然直觉上对连续时间信号 f a ( t ) {f_a}\left( t \right) fa(t) 进行采样获得离散序列 f [ n ] f\left[ n \right] f[n] 的过程会存在信息失真，但这种情况并不是绝对的，这取决于连续时间信号的频谱有效范围以及所使用采样频率。设当 ∣ Ω ∣ > Ω m \left| \Omega \right| > {\Omega _m} ∣Ω∣>Ωm 时，有 F a ( j Ω ) = 0 {F_a}\left( {j\Omega } \right) = 0 Fa(jΩ)=0，那么采样频率 Ω T {\Omega _T} ΩT 只要满足

Ω m ⩽ Ω T / 2 ⇔ Ω T ⩾ 2 Ω m . (3.21) {\Omega _m} \leqslant {\Omega _T}/2 \Leftrightarrow {\Omega _T} \geqslant 2{\Omega _m}. \tag{3.21} Ωm⩽ΩT/2⇔ΩT⩾2Ωm.(3.21)

即可实现从序列 f [ n ] f\left[ n \right] f[n] 到 f a ( t ) {f_a}\left( t \right) fa(t) 的无失真恢复。这个结论也就是带限信号的抽样定理，式 (3.21) 也被称为奈奎斯特 (Nyquist) 条件， 2 Ω m 2{\Omega _m} 2Ωm 称为奈奎斯特率，它确定了从抽样形式中完全恢复 f a ( t ) {f_a}\left( t \right) fa(t) 的最小抽样频率。但要注意的是， F a ( j Ω ) {F_a}\left( {j\Omega } \right) Fa(jΩ) 在 ∣ Ω ∣ = Ω m \left| \Omega \right| = {\Omega _m} ∣Ω∣=Ωm 处不能包含冲激函数（也就是狄拉克函数），因为冲激函数是由极限形式定义的，这意味着 Ω → Ω m + {\Omega \to \Omega _m^ + } Ω→Ωm+ 时 F a ( j Ω ) {F_a}\left( {j\Omega } \right) Fa(jΩ) 并不为 0，此时采样频率要大于 2 Ω m 2{\Omega _m} 2Ωm 才能恢复出 f a ( t ) {f_a}\left( t \right) fa(t)。例如以 Ω T = 2 Ω m {\Omega _T} = 2{\Omega _m} ΩT=2Ωm 对 f a ( t ) = sin ⁡ ( Ω m t ) {f_a}\left( t \right) = \sin \left( {{\Omega _m}t} \right) fa(t)=sin(Ωmt) 进行采样时，有

f [ n ] = f a ( 2 π n Ω T ) = sin ⁡ ( π n ) = 0. f\left[ n \right] = {f_a}\left( {\frac{{2\pi n}}{{{\Omega _T}}}} \right) = \sin \left( {\pi n} \right) = 0. f[n]=fa(ΩT2πn)=sin(πn)=0.

这明显是不可能从 f [ n ] f\left[ n \right] f[n] 恢复 f a ( t ) {f_a}\left( t \right) fa(t) 的。

总之，我们分析了连续时间信号 f a ( t ) {f_a}\left( t \right) fa(t) 通过均匀采样获得离散序列 f [ n ] f\left[ n \right] f[n] 的过程以及连续时间傅里叶变换与离散时间傅里叶变换的内在联系。在这之后，我们主要注重于离散序列 f [ n ] f\left[ n \right] f[n] 本身的分析，而不再考虑采样频率所带来的问题。前面已经提到，对于离散序列 f [ n ] f\left[ n \right] f[n]，假设其绝对可和，定义其离散时间傅里叶变换为

F ( e j ω ) = ∑ n = − ∞ ∞ f [ n ] e − j ω n . (3.22) F\left( {{e^{j\omega }}} \right) = \sum\limits_{n = - \infty }^\infty {f\left[ n \right]{e^{ - j\omega n}}} .\tag{3.22} F(ejω)=n=−∞∑∞f[n]e−jωn.(3.22)

为了推导其反变换形式，因为 F ( e j ω ) F\left( {{e^{j\omega }}} \right) F(ejω) 是以 2 π 2\pi 2π 为周期的，不妨仿照连续时间傅里叶反变换形式，而只在 ( − π , π ) \left( { - \pi ,\pi } \right) (−π,π) 上进行积分，可得

1 2 π ∫ − π π F ( e j ω ) e j ω n d ω = 1 2 π ∫ − π π ( ∑ m = − ∞ ∞ f [ m ] e − j ω m ) e j ω n d ω = ∑ m = − ∞ ∞ f [ m ] ( 1 2 π ∫ − π π e − j ω ( n − m ) d ω ) . (3.23) \begin{aligned} &\frac{1}{{2\pi }}\int\limits_{ - \pi }^\pi {F\left( {{e^{j\omega }}} \right){e^{j\omega n}}d\omega } \\ = &\frac{1}{{2\pi }}\int\limits_{ - \pi }^\pi {\left( {\sum\limits_{m = - \infty }^\infty {f\left[ m \right]{e^{ - j\omega m}}} } \right){e^{j\omega n}}d\omega } \\ = &\sum\limits_{m = - \infty }^\infty {f\left[ m \right]\left( {\frac{1}{{2\pi }}\int\limits_{ - \pi }^\pi {{e^{ - j\omega \left( {n - m} \right)}}d\omega } } \right)} . \\ \end{aligned} \tag{3.23} ==2π1−π∫πF(ejω)ejωndω2π1−π∫π(m=−∞∑∞f[m]e−jωm)ejωndωm=−∞∑∞f[m] 2π1−π∫πe−jω(n−m)dω .(3.23)

当 n = m n = m n=m 时，明显有

1 2 π ∫ − π π e − j ω 0 d ω = 1. (3.24) \frac{1}{{2\pi }}\int\limits_{ - \pi }^\pi {{e^{ - j\omega 0}}d\omega } = 1.\tag{3.24} 2π1−π∫πe−jω0dω=1.(3.24)

当 n ≠ m n \ne m n=m 时，有

1 2 π ∫ − π π e − j ω ( n − m ) d ω = 1 2 π ∫ − π π cos ⁡ ( ω ( n − m ) ) − j sin ⁡ ( ω ( n − m ) ) d ω = 1 2 π ( n − m ) [ sin ⁡ ( ω ( n − m ) ) + j cos ⁡ ( ω ( n − m ) ) ] ∣ − π π = 0. (3.25) \begin{aligned} &\frac{1}{{2\pi }}\int\limits_{ - \pi }^\pi {{e^{ - j\omega \left( {n - m} \right)}}d\omega } \\ = &\frac{1}{{2\pi }}\int\limits_{ - \pi }^\pi {\cos \left( {\omega \left( {n - m} \right)} \right) - j\sin \left( {\omega \left( {n - m} \right)} \right)d\omega } \\ = &\frac{1}{{2\pi \left( {n - m} \right)}}\left. {\left[ {\sin \left( {\omega \left( {n - m} \right)} \right) + j\cos \left( {\omega \left( {n - m} \right)} \right)} \right]} \right|_{ - \pi }^\pi \\ = &0. \\ \end{aligned} \tag{3.25} ===2π1−π∫πe−jω(n−m)dω2π1−π∫πcos(ω(n−m))−jsin(ω(n−m))dω2π(n−m)1[sin(ω(n−m))+jcos(ω(n−m))]∣−ππ0.(3.25)

于是有

f [ n ] = 1 2 π ∫ − π π F ( e j ω ) e j ω n d ω . ((3.26) f\left[ n \right] = \frac{1}{{2\pi }}\int\limits_{ - \pi }^\pi {F\left( {{e^{j\omega }}} \right){e^{j\omega n}}d\omega } .\tag{(3.26} f[n]=2π1−π∫πF(ejω)ejωndω.((3.26)

从而可以总结得序列 f [ n ] f\left[ n \right] f[n] 的离散时间傅里叶变换对为

F ( e j ω ) = F { f [ n ] } = ∑ n = − ∞ ∞ f [ n ] e − j ω n . f [ n ] = F − 1 { F ( e j ω ) } = 1 2 π ∫ − π π F ( e j ω ) e j ω n d ω . (3.27) \begin{gathered} F\left( {{e^{j\omega }}} \right) = \mathcal{F}\left\{ {f\left[ n \right]} \right\} = \sum\limits_{n = - \infty }^\infty {f\left[ n \right]{e^{ - j\omega n}}} . \\ f\left[ n \right] = {\mathcal{F}^{ - 1}}\left\{ {F\left( {{e^{j\omega }}} \right)} \right\} = \frac{1}{{2\pi }}\int\limits_{ - \pi }^\pi {F\left( {{e^{j\omega }}} \right){e^{j\omega n}}d\omega } . \\ \end{gathered} \tag{3.27} F(ejω)=F{f[n]}=n=−∞∑∞f[n]e−jωn.f[n]=F−1{F(ejω)}=2π1−π∫πF(ejω)ejωndω.(3.27)

类似于连续时间下定义的狄拉克函数，我们也可定义离散时间下的冲激序列

δ [ n ] = { 1 , n = 0 , 0 , n ≠ 0. (3.28) \delta \left[ n \right] = \left\{ {\begin{array}{c} {1,}&{n = 0,} \\ {0,}&{n \ne 0.} \end{array}} \right.\tag{3.28} δ[n]={1,0,n=0,n=0.(3.28)

于是有

F { δ [ n − m ] } = ∑ n = − ∞ ∞ δ [ n − m ] e − j ω n = e − j ω m ∑ n = − ∞ ∞ δ [ n − m ] e − j ω ( n − m ) = e − j ω m . (3.29) \begin{gathered} \mathcal{F}\left\{ {\delta \left[ {n - m} \right]} \right\} = \sum\limits_{n = - \infty }^\infty {\delta \left[ {n - m} \right]{e^{ - j\omega n}}} \\ = {e^{ - j\omega m}}\sum\limits_{n = - \infty }^\infty {\delta \left[ {n - m} \right]{e^{ - j\omega \left( {n - m} \right)}}} = {e^{ - j\omega m}}. \\ \end{gathered} \tag{3.29} F{δ[n−m]}=n=−∞∑∞δ[n−m]e−jωn=e−jωmn=−∞∑∞δ[n−m]e−jω(n−m)=e−jωm.(3.29)

刚好对应了式 (3.23) 中的

1 2 π ∫ − π π F { δ [ n − m ] } e j ω n d ω = 1 2 π ∫ − π π e − j ω m e j ω n d ω = 1 2 π ∫ − π π e − j ω ( n − m ) d ω = δ [ n − m ] . (3.30) \begin{aligned} &\frac{1}{{2\pi }}\int\limits_{ - \pi }^\pi {\mathcal{F}\left\{ {\delta \left[ {n - m} \right]} \right\}{e^{j\omega n}}d\omega } \\ = &\frac{1}{{2\pi }}\int\limits_{ - \pi }^\pi {{e^{ - j\omega m}}{e^{j\omega n}}d\omega } \\ = &\frac{1}{{2\pi }}\int\limits_{ - \pi }^\pi {{e^{ - j\omega \left( {n - m} \right)}}d\omega } = \delta \left[ {n - m} \right]. \\ \end{aligned} \tag{3.30} ==2π1−π∫πF{δ[n−m]}ejωndω2π1−π∫πe−jωmejωndω2π1−π∫πe−jω(n−m)dω=δ[n−m].(3.30)

因此，离散冲激序列 δ [ n ] \delta \left[ n \right] δ[n] 在离散时间信号的分析中同样非常重要，其性质可以在计算序列的离散时间傅里叶变换时提供不少便利。

傅里叶级数FS，连续时间傅里叶变换CTFT，离散时间傅里叶变换DTFT，离散傅里叶变换DFT，推导与联系（一）相关推荐

连续时间信号，离散时间信号，模拟信号，数字信号的区别。
最近在准备考研复试,查询了很久连续时间信号.离散时间信号.模拟信号和数字信号的区别,发现网上各种说法不一,其中不乏一些混淆了概念的回答,非常容易误导大家,在这里以熊庆旭老师编写的<信号与系统&g ...
离散时间傅里叶变换，Z变换，离散傅里叶变换性质对比
本篇讲述数字信号处理中离散时间傅里叶变换,z变换,离散傅里叶变换性质的比较.公式是在word中使用mathtype一个一个公式打的,表格无法直接复制到web编辑器中,因此只能上传截图.整理不易,尊重版 ...
第一单元：连续时间信号与离散时间信号
前言 <信号与系统>是一门很难的课,也是许多学校考研要考的专业课,由于每周只有两节课,所以每次上完都要及时的去复习,这里参考的教材是奥本海姆著作,刘海棠译,北京:电子工业出版社,2013年 ...
傅里叶变换一维快速傅里叶变换（快速的一维离散傅里叶变换、分治法）
一.介绍 1.一维离散傅里叶变换DFT. DFT:(Discrete Fourier Transform)离散傅里叶变换是傅里叶变换在时域和频域上都呈离散的形式,将信号的时域采样变换为其DTFT的频域 ...
fft的c语言和matlab对比_傅里叶级数（FS）、傅里叶变换（FT）快速傅里叶变换（FFT）及量子傅里叶变换（QFT）之间推导关系...
1 引言傅里叶级数 (Fourier Series, FS) 是<高等数学>中遇到的一个重要的级数,它可以将任意一个满足狄利克雷条件的函数为一系列三角级数的和.最早由法国数学家傅里叶在研 ...
OPenCV：傅里叶变换、时域和频域、频谱和相位谱、傅里叶级数、离散傅里叶变换(DFT)、频域滤波、高通和低通滤波器、带通和带阻滤波器
日萌社人工智能AI:Keras PyTorch MXNet TensorFlow PaddlePaddle 深度学习实战(不定时更新) 快速了解傅立叶变换(播放PPT即能动态地显示gif图) ...
离散傅里叶变换----解释的最透彻的|包括定义物理意义
1. 傅里叶变换的集中形式及应用傅立叶变换是以时间为自变量的信号和以频率为自变量的频谱函数之间的一种变换关系. 由于自变量时间和频率可以是连续的,也可以是离散的,因此可以组成几种不同的变换对. (1 ...
连续时间采样及采样定理——MATLAB
一.实验目的 1.掌握连续时间信号离散化的方法(即采样),并能利用Matlab编程加以仿真实现; 2.掌握连续时间信号的傅立叶变换和离散时变换的仿真实现方法: 3.学会利用傅里叶变换和离散时间信号的傅 ...
第五章连续时间信号的采样
第五章连续时间信号的采样错题笔记: 冲激响应:对于连续时间 δ ( t ) \delta(t) δ(t),表示为 h ( t ) h(t) h(t):单位脉冲响应:对于离散序列 δ [ n ] \ ...

傅里叶级数FS，连续时间傅里叶变换CTFT，离散时间傅里叶变换DTFT，离散傅里叶变换DFT，推导与联系（一）

1 周期函数的傅里叶级数

2 非周期函数的傅里叶变换

3 离散时间傅里叶变换

傅里叶级数FS，连续时间傅里叶变换CTFT，离散时间傅里叶变换DTFT，离散傅里叶变换DFT，推导与联系（一）相关推荐

最新文章

热门文章