本文主要从傅里叶级数 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\left(t\right)$ 是一个以 $T$ 为周期的函数，且在 $[-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 , l 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} , l \\\\ where\\;{\\omega _0} = 2\\pi /T. \\\\ \\end{gathered} \\tag{1.1} an​=T2​−T/2∫T/2​f(t)cos(nω0​t)dt,bn​=T2​−T/2∫T/2​f(t)sin(nω0​t)dt,lwhereω0​=2π/T.​(1.1)

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

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∑∞​(an​cos(nω0​t)+bn​sin(nω0​t)),(1.2)

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

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ω0​t)=2ejnω0​t+e−jnω0​t​,sin(nω0​t)=2jejnω0​t−e−jnω0​t​.​(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∑∞​[an​cos(nω0​t)+bn​sin(nω0​t)]=2a0​​+n=1∑∞​(2an​−jbn​​ejnω0​t+2an​+jbn​​e−jnω0​t)=2a0​−jb0​​ej0ω0​t+n=1∑∞​(2an​−jbn​​ejnω0​t+2a−n​−jb−n​​e−jnω0​t)=n=−∞∑∞​2an​−jbn​​ejnω0​t=n=−∞∑∞​Fn​ejnω0​t.​(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/2​f(t)[cos(nω0​t)−jsin(nω0​t)]dt=T1​−T/2∫T/2​f(t)(2ejnω0​t+e−jnω0​t​−j2jejnω0​t−e−jnω0​t​)dt=T1​−T/2∫T/2​f(t)e−jnω0​tdt.​(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/2​f(t)e−jnω0​tdt,f(t)=F−1{F[n]}=n=−∞∑∞​Fn​ejnω0​t.​(1.7)

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

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

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

如果考虑非周期函数 $f\left(t\right)$，我们可认为其周期 $T\to \infty$，于是有基频 ω 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→∞lim​f(t)​=T→∞lim​n=−∞∑∞​Δωn​Fn​​Δωn​ejωn​t=T→∞lim​n=−∞∑∞​ ​T⋅Δωn​1​−T/2∫T/2​f(τ)e−jωn​τdτ ​ejωn​tΔωn​=T→∞lim​n=−∞∑∞​ ​2π1​−T/2∫T/2​f(τ)e−jωn​τdτ ​ejωn​tΔωn​=2π1​−∞∫∞​ ​−∞∫∞​f(τ)e−jωτdτ ​ejωtdω.​(2.1)

于是可以定义，设非周期函数 $f\left(t\right)$ 在 $\left(-\infty ,+\infty \right)$ 上绝对可积，则 $f\left(t\right)$ 的傅里叶变换 (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\left(j\omega \right)$ 也可直接写作 $F\left(\omega \right)$，加上虚数符号主要是为了表达傅里叶变换的结果为复数，同时也与后面的离散时间傅里叶变换做区分。

相比傅里叶级数所得离散系数 $F\left[n\right]$，傅里叶变换 $F\left(j\omega \right)$ 定义在连续的频率自变量 $\omega$上，且 $\omega$与 $n$ 的量纲也有所区别，$\omega$ 是由 lim ⁡ ω 0 → 0 n ω 0 {\\mathop {\\lim }\\limits_{{\\omega _0} \\to 0} n{\\omega _0}} ω0​→0lim​nω0​ 这个极限来定义的。考虑定义在有限区间 $\left[-T/2,T/2\right]$ 的函数 $f\left(t\right)$，对其以 $T$ 为周期进行延拓获得周期函数 f ~ ( t ) \\tilde f\\left( t \\right) f~​(t)，那么可得 $f\left(t\right)$ 的傅里叶变换谱 $F\left(j\omega \right)$ 与 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/2​f(t)e−jωtdt,F~[n]=T1​−T/2∫T/2​f~​(t)e−jnΩT​tdt=T1​F(jω)∣ω=nΩT​​,whereΩT​=2π/T.​(2.3)

也就是说，周期函数的傅里叶系数 F ~ [ n ] \\tilde F\\left[ n \\right] F~[n] 是对该周期函数在一个周期内进行傅里叶变换所得 $F\left(j\omega \right)$ 以 Ω T {\\Omega _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\to 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​ε1​dt=1.(3.2)

假设函数 $f\left(t\right)$ 连续，因为 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→t0​lim​f(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​=ε→0lim​t0​−ε/2∫t0​+ε/2​f(t)ε1​dt=f(t0​)ε→0lim​t0​−ε/2∫t0​+ε/2​ε1​dt=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τ=ε→0lim​t−t0​−ε/2∫t−t0​+ε/2​f(τ)ε1​dτ=f(t−t0​)ε→0lim​t−t0​−ε/2∫t−t0​+ε/2​ε1​dτ=f(t−t0​).​(3.4)

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

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

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)=ε1​nT−ε/2∫nT+ε/2​f(t)dt.(3.5)

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

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

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\left(t\right)$ 的乘积，即

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$ 是连续的，并没有说信号本身是连续的。另外，由于 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$，从而简化后续的运算。为了区别于一般的连续时间信号的傅里叶变换，我们用 $\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\left[n\right]$ 的离散时间傅里叶变换 (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$，但这并不代表不需要对其进行进一步的讨论。因为采样是从连续时间到离散时间的转换过程，直觉上这会导致一部分信息丢失，那么该如何去证明以及衡量这种信息失真是我们接下来要讨论的重点。

可以发现，不同于一般连续时间信号的傅里叶变换，由周期狄拉克函数加权叠加所定义的连续时间信号 f p ( t ) {f_p}\\left( t \\right) fp​(t) 的傅里叶变换 F p ( j Ω ) {F_p}\\left( {j\\Omega } \\right) Fp​(jΩ)，以及序列 $f\left[n\right]$ 的离散时间傅里叶变换 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ΩT​2π​=n=−∞∑∞​fa​(nT)e−j(Ω+kΩT​)nΩT​2π​=Fp​(j(Ω+kΩT​)),=n=−∞∑∞​f[n]e−jωn=n=−∞∑∞​f[n]e−j(ω+2kπ)n=F(ej(ω+2kπ))=F(ej(ω+kΩT​T)),=ΩT,ΩT​=2π/T.​(3.14)