离散时间信号的傅里叶变换与频域分析
序列的傅里叶变换
-
定义
现有一个连续时间信号 x a ( t ) x_a(t) xa(t),采样后的信号为 x ^ a ( t ) \widehat{x}_a(t) x a(t),对应的频谱分别为 X a ( j Ω ) 、 X ^ a ( j Ω ) X_a(j\Omega)、\widehat{X}_a(j\Omega) Xa(jΩ)、X a(jΩ)
根据抽样定理, X ^ a ( j Ω ) \widehat{X}_a(j\Omega) X a(jΩ)为 X a ( j Ω ) X_a(j\Omega) Xa(jΩ)的以采样频率为周期的周期延拓:
X ^ a ( j Ω ) = 1 T ∑ k = − ∞ + ∞ X a ( j Ω − j k Ω s ) \widehat{X}_a(j\Omega)=\frac{1}{T}\sum_{k=-\infin}^{+\infin}X_a(j\Omega-jk\Omega_s) X a(jΩ)=T1k=−∞∑+∞Xa(jΩ−jkΩs)
现在对 x ^ a ( t ) \widehat{x}_a(t) x a(t)进行傅里叶变换:
X ^ a ( j Ω ) = ∑ − ∞ + ∞ x a ( n T ) e − j Ω n T \widehat{X}_a(j\Omega)=\sum_{-\infin}^{+\infin}x_a(nT)e^{-j\Omega nT} X a(jΩ)=−∞∑+∞xa(nT)e−jΩnT
进行归一化处理,令采样周期T=1s,采样频率f=1Hz, Ω s = 2 π r a d / s \Omega_s=2\pi rad/s Ωs=2πrad/s,则 ω = Ω f s \omega=\frac{\Omega}{f_s} ω=fsΩ,那么:
X ^ a ( e j ω ) = ∑ − ∞ + ∞ x a ( n ) e − j ω n \widehat{X}_a(e^{j\omega})=\sum_{-\infin}^{+\infin}x_a(n)e^{-j\omega n} X a(ejω)=−∞∑+∞xa(n)e−jωn
x ( n ) = 1 2 π ∫ − π π X ( e j ω ) e j ω n d ω x(n)=\frac{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 a ( n ) = x a ( n T ) x_a(n)=x_a(nT) xa(n)=xa(nT),上式即为序列的傅里叶正变换和反变换。 -
性质
- 线性:若
X
1
(
e
j
ω
)
=
F
T
[
x
1
(
n
)
]
,
X
2
(
e
j
ω
)
=
F
T
[
x
2
(
n
)
]
X_1(e^{j\omega})=FT[x_1(n)],X_2(e^{j\omega})=FT[x_2(n)]
X1(ejω)=FT[x1(n)],X2(ejω)=FT[x2(n)],则有
F T [ a x 1 ( n ) + b x 2 ( n ) ] = a X 1 ( e j ω ) + b X 2 ( e j ω ) FT[ax_1(n)+bx_2(n)]=aX_1(e^{j\omega})+bX_2(e^{j\omega}) FT[ax1(n)+bx2(n)]=aX1(ejω)+bX2(ejω)
- 线性:若
X
1
(
e
j
ω
)
=
F
T
[
x
1
(
n
)
]
,
X
2
(
e
j
ω
)
=
F
T
[
x
2
(
n
)
]
X_1(e^{j\omega})=FT[x_1(n)],X_2(e^{j\omega})=FT[x_2(n)]
X1(ejω)=FT[x1(n)],X2(ejω)=FT[x2(n)],则有
-
时移与频移:若 X ( e j ω ) = F T [ x ( n ) ] X(e^{j\omega})=FT[x(n)] X(ejω)=FT[x(n)],则有
F T [ x ( n − n 0 ) ] = e − j ω 0 X ( e j ω ) FT[x(n-n_0)]=e^{-j\omega_0}X(e^{j\omega}) FT[x(n−n0)]=e−jω0X(ejω)
F T [ e j ω 0 n x ( n ) ] = X [ e j ( ω − ω 0 ) ] FT[e^{j\omega_0n}x(n)]=X[e^{j(\omega-\omega_0)}] FT[ejω0nx(n)]=X[ej(ω−ω0)]- 周期性: X ( e j ω ) X(e^{j\omega)} X(ejω)是 w w w的周期函数,周期为 2 π 2\pi 2π
- 对称性:若
x
(
n
)
x(n)
x(n)满足
x
(
n
)
=
x
∗
(
−
n
)
x(n)=x^*(-n)
x(n)=x∗(−n)则称
x
(
n
)
x(n)
x(n)为共轭对称序列,一般用
x
e
(
n
)
x_e(n)
xe(n)表示;若
x
(
n
)
x(n)
x(n)满足
x
(
n
)
=
−
x
∗
(
−
n
)
x(n)=-x^*(-n)
x(n)=−x∗(−n)则称
x
(
n
)
x(n)
x(n)为共轭反对称序列,一般用
x
o
(
n
)
x_o(n)
xo(n)表示。
共轭对称序列实部为偶数,虚部为奇数;
共轭反对称序列实部是奇数,虚部为偶数;
任何序列都可以表示为一个共轭对称序列与一个共轭反对称序列之和:
x e ( n ) = 1 2 [ x ( n ) + x ∗ ( − n ) ] x_e(n)=\frac{1}{2}[x(n)+x^*(-n)] xe(n)=21[x(n)+x∗(−n)]
x o ( n ) = 1 2 [ x ( n ) − x ∗ ( − n ) ] x_o(n)=\frac{1}{2}[x(n)-x^*(-n)] xo(n)=21[x(n)−x∗(−n)]
x ( n ) = x e ( n ) + x o ( n ) x(n)=x_e(n)+x_o(n) x(n)=xe(n)+xo(n)
对上式两边进行傅里叶变换,在频域中也有类似的结论:
X e ( e j ω ) = 1 2 [ X ( e j ω ) + X ∗ ( e j ω ) ] = X R ( e j ω ) X_e(e^{j\omega})=\frac{1}{2}[X(e^{j\omega})+X^*(e^{j\omega})]=X_R(e^{j\omega}) Xe(ejω)=21[X(ejω)+X∗(ejω)]=XR(ejω)
X o ( e j ω ) = 1 2 [ X ( e j ω ) − X ∗ ( e j ω ) ] = X I ( e j ω ) X_o(e^{j\omega})=\frac{1}{2}[X(e^{j\omega})-X^*(e^{j\omega})]=X_I(e^{j\omega}) Xo(ejω)=21[X(ejω)−X∗(ejω)]=XI(ejω)
X ( e j ω ) = X e ( e j ω ) + X o ( e j ω ) X(e^{j\omega})=X_e(e^{j\omega})+X_o(e^{j\omega}) X(ejω)=Xe(ejω)+Xo(ejω)
其中 F T [ x ∗ ( − n ) ] = ∑ − ∞ + ∞ x ∗ ( − n ) e − j ω n = [ ∑ − ∞ + ∞ x ( − n ) e j ω n ] ∗ = X ∗ ( e j ω ) , X R ( e e j ω ) 、 X I ( e j ω ) FT[x^*(-n)]=\sum_{-\infin}^{+\infin}x^*(-n)e^{-j\omega n}=[\sum_{-\infin}^{+\infin}x(-n)e^{j\omega n}]^*=X^*(e^{j\omega}),X_R(e^{e^{j\omega}})、X_I(e^{j\omega}) FT[x∗(−n)]=∑−∞+∞x∗(−n)e−jωn=[∑−∞+∞x(−n)ejωn]∗=X∗(ejω),XR(eejω)、XI(ejω)分别为x(n)实部与虚部(包括j)的傅里叶变换。 - 时域卷积定理
若 F T [ x ( n ) ] = X ( e j ω ) , F T [ h ( n ) ] = H ( e j ω ) , y ( n ) = x ( n ) ∗ h ( n ) FT[x(n)]=X(e^{j\omega}),FT[h(n)]=H(e^{j\omega}),y(n)=x(n)*h(n) FT[x(n)]=X(ejω),FT[h(n)]=H(ejω),y(n)=x(n)∗h(n),则:
F T [ y ( n ) ] = X ( e j ω ) H ( e j ω ) FT[y(n)]=X(e^{j\omega})H(e^{j\omega}) FT[y(n)]=X(ejω)H(ejω) - 频域卷积定理
若 F T [ x ( n ) ] = X ( e j ω ) , F T [ h ( n ) ] = H ( e j ω ) , y ( n ) = x ( n ) h ( n ) FT[x(n)]=X(e^{j\omega}),FT[h(n)]=H(e^{j\omega}),y(n)=x(n)h(n) FT[x(n)]=X(ejω),FT[h(n)]=H(ejω),y(n)=x(n)h(n),则:
F T [ y ( n ) ] = 1 2 π X ( e j ω ) ∗ H ( e j ω ) FT[y(n)]=\frac{1}{2\pi}X(e^{j\omega})*H(e^{j\omega}) FT[y(n)]=2π1X(ejω)∗H(ejω) - 帕塞瓦尔定理
若 F T [ x ( n ) ] = X ( e j ω ) FT[x(n)]=X(e^{j\omega}) FT[x(n)]=X(ejω),则:
∑ − ∞ + ∞ ∣ x ( n ) ∣ 2 = 1 2 π ∫ − π + π ∣ X ( e j ω ) ∣ 2 d ω \sum_{-\infin}^{+\infin}|x(n)|^2=\frac{1}{2\pi}\int_{-\pi}^{+\pi}|X(e^{j\omega})|^2d\omega −∞∑+∞∣x(n)∣2=2π1∫−π+π∣X(ejω)∣2dω
序列的Z变换
序列 Z变换 收敛域 δ(n) 1 0≤|z|≤∞ δ(n-k) z-k 0≤|z|≤∞ u(n) 1/(1-z-1) |z|>1 RN(n) (1-z-N)/(1-z-1) |z|>0 nu(n) z/(z-1)2 |z|>1 anu(n) z/(z-a) |z|>|a| -anu(-n-1) z/(z-a) |z|<|a|