周期函数可以展开成在完备正交信号空间中的无穷级数,所展开的无穷级数统称为傅里叶级数。设有周期函数 x ~ a ( t ) \widetilde{x}_a(t) x a(t),它的周期是T,角频率 Ω = 2 π F = 2 π T \Omega=2{\pi}F=\dfrac{2\pi}{T} Ω=2πF=T2π,它可分解为
x ~ a ( t ) = a 0 2 + ∑ n = 1 ∞ a n c o s ( n Ω t ) + ∑ n = 1 ∞ b n s i n ( n Ω t ) = A 0 2 + ∑ n = 1 ∞ A n c o s ( n Ω t + φ n ) \widetilde{x}_a(t)=\dfrac{a_0}{2}+\sum\limits_{n=1}^{\infty}a_ncos(n{\Omega}t)+\sum\limits_{n=1}^{\infty}b_nsin(n{\Omega}t)=\dfrac{A_0}{2}+\sum\limits_{n=1}^{\infty}{A_n}cos(n{\Omega}t+\varphi_n) x a(t)=2a0+n=1∑∞ancos(nΩt)+n=1∑∞bnsin(nΩt)=2A0+n=1∑∞Ancos(nΩt+φn)
傅里叶级数 a n = 2 T ∫ − T 2 T 2 x ~ a ( t ) c o s ( n Ω t ) d t , n = 0 , 1 , 2 , . . . a_n=\dfrac{2}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}\widetilde{x}_a(t)cos(n{\Omega}t)dt,{\qquad}n=0,1,2,... an=T2∫−2T2Tx a(t)cos(nΩt)dt,n=0,1,2,...
b n = 2 T ∫ − T 2 T 2 x ~ a ( t ) s i n ( n Ω t ) d t , n = 1 , 2 , . . . b_n=\dfrac{2}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}\widetilde{x}_a(t)sin(n{\Omega}t)dt,{\qquad}n=1,2,... bn=T2∫−2T2Tx a(t)sin(nΩt)dt,n=1,2,...
A n = a n 2 + b n 2 , φ n = − a r c t a n ( b n a n ) A_n=\sqrt{a_n^2+b_n^2},{\qquad}\varphi_n=-arctan(\dfrac{b_n}{a_n}) An=an2+bn2,φn=−arctan(anbn)
x ~ a ( t ) = ∑ n = − ∞ ∞ X n e j n Ω t \widetilde{x}_a(t)=\sum\limits _{n=-\infty}^{\infty}X_ne^{jn{\Omega}t} x a(t)=n=−∞∑∞XnejnΩt
傅里叶级数
X n = 1 T ∫ − T 2 T 2 x ~ a ( t ) e − j n Ω t d t , n = 0 , ± 1 , ± 2 , . . . X_n=\dfrac{1}{T}\int_{\frac{-T}{2}}^{\frac{T}{2}} \widetilde{x}_a(t)e^{-jn{\Omega}t}dt,{\qquad}n=0,\pm1,\pm2,... Xn=T1∫2−T2Tx a(t)e−jnΩtdt,n=0,±1,±2,...
如果周期性脉冲的重复周期足够长,使得后一个脉冲到来之前,前一个脉冲的作用实际上早已消息,这样的周期信号可以作为非周期信号处理。为了描述非周期信号的频谱特性,引入频谱密度的概念。当 T → ∞ T\to\infty T→∞时,频谱密度函数为
X a ( j ω ) = lim T → ∞ X n T = ∫ − ∞ ∞ x a ( t ) e − j ω t d t X_a(j\omega)=\lim_{T\to\infty}X_nT=\int_{-\infty}^{\infty} x_a(t)e^{-j{\omega}t}dt Xa(jω)=T→∞limXnT=∫−∞∞xa(t)e−jωtdt
傅里叶逆变换为 x a ( t ) = 1 2 π ∫ − ∞ ∞ X a ( j ω ) e j ω t d ω x_a(t)=\dfrac{1}{2\pi}\int _{-\infty}^{\infty}X_a(j\omega)e^{j{\omega}t}d\omega xa(t)=2π1∫−∞∞Xa(jω)ejωtdω
傅里叶变换存在的充分条件是在无线区间内 x a ( t ) x_a(t) xa(t)绝对可积,即 ∫ − ∞ ∞ ∣ x a ( t ) ∣ d t < ∞ \int_{-\infty}^{\infty}|x_a(t)|dt<\infty ∫−∞∞∣xa(t)∣dt<∞但它非必要条件。
现在考虑一个周期为T的周期函数 x ~ a ( t ) \widetilde{x}_a(t) x a(t),它可以展开成指数形式的傅里叶级数 x ~ a ( t ) = ∑ n = − ∞ ∞ X n e j n Ω t \widetilde{x}_a(t)=\sum\limits_{n=-\infty}^{\infty}X_{n}e^{jn{\Omega}t} x a(t)=n=−∞∑∞XnejnΩt
对该式的等号两端取傅立叶变换,得
X ( j ω ) = 2 π ∑ n = − ∞ ∞ X n δ ( ω − n Ω ) X(j\omega)=2{\pi}\sum\limits_{n=-\infty}^{\infty}X_n\delta({\omega-n{\Omega})} X(jω)=2πn=−∞∑∞Xnδ(ω−nΩ)
上式表明,周期信号的傅里叶变换由无穷多个冲奇函数组成,这些冲激函数位于信号的各谐波角频率 n Ω ( n = 0 , ± 1 , ± 2 , . . . ) n\Omega(n=0,\pm1,\pm2,...) nΩ(n=0,±1,±2,...)处,其强度为各相应幅度 F n F_n Fn的 2 π 2\pi 2π倍。
为了解决傅里叶变换的局限性,引入复频率进行复频域分析。因果信号 x a ( t ) ε ( t ) x_a(t)\varepsilon(t) xa(t)ε(t),记 s = α + j ω s=\alpha+j\omega s=α+jω,其拉普拉斯变换为
X a ( s ) = ∫ 0 − ∞ x ( t ) e − s t d t X_a(s)=\int_{0_-}^{\infty}x(t)e^{-st}dt Xa(s)=∫0−∞x(t)e−stdt
其逆变换为
x ( t ) = 1 2 π j ∫ α − j ∞ α + j ∞ X a ( s ) e a t d s x(t)=\dfrac{1}{2{\pi}j}\int_{\alpha-j\infty}^{\alpha+j\infty}X_a(s)e^{at}ds x(t)=2πj1∫α−j∞α+j∞Xa(s)eatds
为使拉普拉斯积分式绝对且一致收敛,因果函数 x a ( t ) ε ( t ) x_a(t)\varepsilon(t) xa(t)ε(t)必须满足:(1)在有限区间 a < t < b a<t<b a<t<b内(其中 0 ≤ a < b < ∞ 0{\leq}a<b<\infty 0≤a<b<∞)可积,(2)对于某个 α 0 \alpha_0 α0有 lim t → ∞ ∣ x a ( t ) ε ( t ) ∣ e − α t = 0 , α < α 0 \lim_{t\to\infty}|x_a(t)\varepsilon(t)|e^{-{\alpha}t}=0,\alpha<\alpha_0 t→∞lim∣xa(t)ε(t)∣e−αt=0,α<α0此时对于 R e [ s ] = α < α 0 Re[s]=\alpha<\alpha_0 Re[s]=α<α0,拉普拉斯积分式绝对且一致收敛。
我们往往可以通过周期采样一个模拟信号来得到序列,这样一来,序列中的第n个数的数值就等于模拟信号 x a ( t ) x_a(t) xa(t)在时刻nT处的值,即 x ( n ) = x a ( t ) , − ∞ < n < ∞ x(n)=x_a(t),{\quad}-\infty<n<\infty x(n)=xa(t),−∞<n<∞T称为采样周期,其倒数就是采样频率。
定义非周期序列 x ( n ) x(n) x(n)的离散时间傅里叶变换为
X ( e j ω ) = ∑ n = − ∞ ∞ x ( n ) e − j ω n X(e^{j\omega})=\sum\limits_{n=-\infty}^{\infty}x(n)e^{-j{\omega}n} X(ejω)=n=−∞∑∞x(n)e−jωn
其逆变换为 x ( n ) = 1 2 π ∫ − π π X ( e j ω ) e j ω n d ω x(n)=\dfrac{1}{2\pi}\int_{-\pi}^{\pi}X(e^{j\omega})e^{j{\omega}n}d\omega x(n)=2π1∫−ππX(ejω)ejωndω
它把序列 x ( n ) x(n) x(n)表示成频率在 2 π 2\pi 2π的区间范围内,由 X ( e j ω ) X(e^{j\omega}) X(ejω)确定每一个复正弦分量相对大小的无限小复正弦的叠加。
由于离散时间傅里叶变换不是对所有的序列都收敛,所以将其推广到z域。Z变换 X ( z ) X(z) X(z)定义为
X ( z ) = ∑ n = − ∞ ∞ x ( n ) z − n X(z)=\sum\limits_{n=-\infty}^{\infty}x(n)z^{-n} X(z)=n=−∞∑∞x(n)z−n
其逆变换为
x ( n ) = 1 2 π j ∮ c X ( z ) z n − 1 d z x(n)=\dfrac{1}{2{\pi}j}\oint_cX(z)z^{n-1}dz x(n)=2πj1∮cX(z)zn−1dz
由定义式可见当 z = e j w z=e^{jw} z=ejw时,Z变换就变为了离散时间傅里叶变换,或者说离散时间傅里叶变换就是在单位圆上进行Z变换。
将序列 x ( n ) x(n) x(n)进行周期延拓我们就可以到到周期序列 x ~ ( n ) \widetilde{x}(n) x (n),其周期为N,即有 x ~ ( n ) = x ~ ( n + l N ) , ( l 为 任 意 整 数 ) \widetilde{x}(n)=\widetilde{x}(n+lN),\quad(l为任意整数) x (n)=x (n+lN),(l为任意整数)
通过周期信号的傅里叶级数可以推导出周期序列的离散傅里叶级数为 X ~ ( k ) = ∑ n = 0 N − 1 x ~ ( n ) e − j 2 π N k n \widetilde{X}(k)=\sum\limits_{n=0}^{N-1}\widetilde{x}(n)e^{-j\frac{2\pi}{N}kn} X (k)=n=0∑N−1x (n)e−jN2πkn
逆变换为 x ~ ( n ) = 1 N ∑ k = 0 N − 1 X ~ ( k ) e j 2 π N k n \widetilde{x}(n)=\dfrac{1}{N}\sum\limits_{k=0}^{N-1}\widetilde{X}(k)e^{j\frac{2\pi}{N}kn} x (n)=N1k=0∑N−1X (k)ejN2πkn
设长度为N的有限长序列 x ( n ) x(n) x(n) 的区间为 [ 0 , N − 1 ] [0,N-1] [0,N−1],其余各处皆为0,即 x ( n ) = { x ( n ) , 0 ≤ n ≤ N − 1 0 , n为其余值 x(n)= \begin{cases} x(n),& \text{$0{\leq}n{\leq}N-1$}\\\\ 0,& \text{n为其余值} \end{cases} x(n)=⎩⎪⎨⎪⎧x(n),0,0≤n≤N−1n为其余值可以将其看作周期序列 x ~ ( n ) \widetilde{x}(n) x (n) 的一个周期,其离散傅里叶变换定义为
X ( k ) = ∑ n = 0 N − 1 x ( n ) e − j 2 π N k n 0 ≤ k ≤ N − 1 X(k)=\sum\limits_{n=0}^{N-1}x(n)e^{-j\frac{2\pi}{N}kn}{\qquad}0{\leq}k{\leq}N-1 X(k)=n=0∑N−1x(n)e−jN2πkn0≤k≤N−1
反变换为 x ( n ) = 1 N ∑ k = 0 N − 1 X ( k ) e j 2 π N k n 0 ≤ n ≤ N − 1 x(n)=\dfrac{1}{N}\sum\limits_{k=0}^{N-1}X(k)e^{j\frac{2\pi}{N}kn}{\qquad}0{\leq}n{\leq}N-1 x(n)=N1k=0∑N−1X(k)ejN2πkn0≤n≤N−1
可以看出DFT即在区间 [ 0 , 2 π ) [0,2\pi) [0,2π)内对DTFT进行N点等间隔采样,也就是在单位圆上对Z变换进行N点等间隔采样。
与周期信号的傅里叶变换类似,我们可以将周期序列的离散傅里叶变换看作是频域的脉冲幅值正比于序列DFS系数的一个脉冲串。 x ~ ( n ) \widetilde{x}(n) x (n) 离散傅里叶变换定义为脉冲串
X ~ ( e j ω ) = ∑ k = − ∞ ∞ 2 π N X ~ ( k ) δ ( ω − 2 π k N ) \widetilde{X}(e^{j\omega})=\sum\limits_{k=-\infty}^{\infty}\dfrac{2\pi}{N}\widetilde{X}(k)\delta(\omega-\dfrac{2{\pi}k}{N}) X (ejω)=k=−∞∑∞N2πX (k)δ(ω−N2πk)
可见 X ~ ( e j ω ) \widetilde{X}(e^{j\omega}) X (ejω)是周期的,周期为 2 π 2\pi 2π。这是因为 X ~ ( k ) \widetilde{X}(k) X (k)的周期为N,是 2 π N \frac{2\pi}{N} N2π的整数倍。
整理这篇文章的初衷是女朋友叫我帮她找这些变换的关系图,我在百度上找了好久,始终没有一个较为完善的版本,于是就萌发出了自己写一篇的想法。在整理的过程中,发现自己对很多变换的定义和关系掌握的其实很差,遇到了很多问题,也因此又翻阅了当时信号与系统的教材和奥本海姆的离散时间信号处理。
文章中我只简单的列举了各个变换的定义式,具体条件和性质大家可以自行翻阅教材。关于变换关系图,如果有错误,欢迎在留言区指正。