连续时间信号的抽样及其重建

  现研究一连续信号进行抽样转换为数字信号,经数字信号处理器(DSP)或计算机处理后,再进行重建的过程,具体过程如下:

连续时间信号的抽样及其重建_第1张图片

其中采样/保持电路和A/D转换电路可以看做是一个理想抽样的过程,而D/A转换和平滑录播可以看做是一个理想内插的过程。

  假设理想抽样信号为
∑ n = − ∞ ∞ δ ( t − n T s ) \sum_{n=-\infty}^{\infty}\delta(t-nT_s) n=δ(tnTs)
其中 T s T_s Ts为抽样的周期。那么模拟信号 x a ( t ) x_a(t) xa(t)经理想抽样后得到的抽样信号 x ^ a ( t ) \hat{x}_a(t) x^a(t)
x ^ a ( t ) = x a ( t ) ⋅ ∑ n = − ∞ ∞ δ ( t − n T s ) \hat{x}_a(t)=x_a(t) \cdot \sum_{n=-\infty}^{\infty}\delta(t-nT_s) x^a(t)=xa(t)n=δ(tnTs)
设信号 x a ( t ) x_a(t) xa(t)的傅里叶变换为 X a ( j Ω ) X_a(j\Omega) Xa(jΩ),并且其最高频率为 Ω m \Omega_m Ωm,现研究抽样信号 x ^ a ( t ) \hat{x}_a(t) x^a(t)的傅里叶变换。
X ^ a ( j Ω ) = F [ x ^ a ( t ) ] = F [ x a ( t ) ⋅ ∑ n = − ∞ ∞ δ ( t − n T s ) ] = 1 2 π X ( j Ω ) ∗ 2 π T s ∑ n = − ∞ ∞ δ ( Ω − n Ω s ) = 1 T s ∑ n = − ∞ ∞ X ( j ( Ω − n Ω s ) ) \begin{aligned} \hat{X}_a(j\Omega)=F[\hat{x}_a(t)]&=F[x_a(t) \cdot \sum_{n=-\infty}^{\infty}\delta(t-nT_s)] \\ &=\frac{1}{2\pi}X(j\Omega)*\frac{2\pi}{T_s}\sum_{n=-\infty}^{\infty}\delta(\Omega - n\Omega_s) \\ &=\frac{1}{T_s}\sum_{n=-\infty}^{\infty}X(j(\Omega - n\Omega_s)) \end{aligned} X^a(jΩ)=F[x^a(t)]=F[xa(t)n=δ(tnTs)]=2π1X(jΩ)Ts2πn=δ(ΩnΩs)=Ts1n=X(j(ΩnΩs))
其中 Ω s = 2 π T s \Omega_s=\frac{2\pi}{T_s} Ωs=Ts2π

  从上式中就可以看出抽样信号的频谱是原信号频谱的周期延拓。

连续时间信号的抽样及其重建_第2张图片

要保证频谱在周期延拓时不发生混叠,那么就要求
Ω s − Ω m ≥ Ω m ⇒ Ω s ≥ 2 Ω m \Omega_s-\Omega_m\geq\Omega_m \Rightarrow\Omega_s\geq2\Omega_m ΩsΩmΩmΩs2Ωm
  把 Ω s = 2 Ω m \Omega_s=2\Omega_m Ωs=2Ωm称为奈奎斯特采样频率,这是频谱不发生混叠允许的最小采样频率,此时可以通过一低通滤波器将信号恢复出来,若频谱发生了混叠,则很难将信号重建出来。

  考虑信号的重建,由频谱图可知,通过一低通滤波器即可将信号完全的恢复出来,假设以频率 Ω s > 2 Ω m \Omega_s>2\Omega_m Ωs>2Ωm进行抽样,考虑这么一个低通滤波器:

连续时间信号的抽样及其重建_第3张图片

其傅里叶反变换为
h L P ( t ) = s i n c ( t T s ) h_{LP}(t)=sinc(\frac{t}{T_s}) hLP(t)=sinc(Tst)
其中 s i n c ( t ) = s i n ( π t ) π t sinc(t)=\frac{sin(\pi t)}{\pi t} sinc(t)=πtsin(πt)

  由频谱关系知
X ( j Ω ) = X ^ a ( j Ω ) ⋅ H L P ( j Ω ) X(j\Omega)=\hat{X}_a(j\Omega) \cdot H_{LP}(j\Omega) X(jΩ)=X^a(jΩ)HLP(jΩ)
所以
x a ( t ) = x ^ a ( t ) ∗ h L P ( t ) = x a ( t ) ⋅ ∑ n = − ∞ ∞ δ ( t − n T s ) ∗ h L P ( t ) = ∑ n = − ∞ ∞ x a ( n T s ) δ ( t − n T s ) ∗ s i n c ( t T s ) = ∑ n = − ∞ ∞ x a ( n T s ) s i n c ( 1 T s ( t − n T s ) ) \begin{aligned} x_a(t)&=\hat{x}_a(t) * h_{LP}(t) \\ &=x_a(t) \cdot \sum_{n=-\infty}^{\infty}\delta(t-nT_s) * h_{LP}(t) \\ &=\sum_{n=-\infty}^{\infty}x_a(nT_s)\delta(t-nT_s) *sinc(\frac{t}{T_s}) \\ &=\sum_{n=-\infty}^{\infty}x_a(nT_s)sinc(\frac{1}{T_s}(t-nT_s)) \end{aligned} xa(t)=x^a(t)hLP(t)=xa(t)n=δ(tnTs)hLP(t)=n=xa(nTs)δ(tnTs)sinc(Tst)=n=xa(nTs)sinc(Ts1(tnTs))
这就是信号的重建,这个过程也被称为理想内插过程。

你可能感兴趣的:(数字信号处理,数字信号处理)