信号与系统复习笔记——采样与通讯系统
采样定理
冲激串采样函数可表示为:
p ( t ) = ∑ n = − ∞ + ∞ δ ( t − n T ) p(t) = \sum_{n=-\infty}^{+\infty} \delta(t - nT) p(t)=n=−∞∑+∞δ(t−nT)
周期 T T T 称为采样周期,而 ω s = 1 T \omega_s = \frac{1}{T} ωs=T1 为采样频率。
对信号 x ( t ) x(t) x(t) 进行冲激串采样的结果为:
x p ( t ) = x ( t ) p ( t ) = ∑ n = − ∞ + ∞ x ( n T ) δ ( t − n T ) x_p(t) = x(t)p(t) = \sum_{n=-\infty}^{+\infty} x(nT) \delta(t - nT) xp(t)=x(t)p(t)=n=−∞∑+∞x(nT)δ(t−nT)
x p ( t ) x_p(t) xp(t) 本身仍是一个冲激串,幅度为 x ( t ) x(t) x(t) 在采样点上的值。
并且有:
X p ( j ω ) = 1 2 π [ X ( j ω ) ∗ P ( j ω ) ] X_p(j\omega) = \frac{1}{2\pi}[X(j\omega) \ast P(j\omega)] Xp(jω)=2π1[X(jω)∗P(jω)]
其中:
P ( j ω ) = 2 π T ∑ k = − ∞ + ∞ δ ( ω − k ω s ) P(j\omega) = \frac{2 \pi}{T} \sum_{k = -\infty}^{+\infty}\delta(\omega - k\omega_s) P(jω)=T2πk=−∞∑+∞δ(ω−kωs)
已知信号和单位冲激函数做卷积等于该信号的位移:
X ( j ω ) ∗ δ ( ω − ω 0 ) = X ( j ( ω − ω 0 ) ) X(j\omega) \ast \delta(\omega - \omega_0) = X(j(\omega - \omega_0)) X(jω)∗δ(ω−ω0)=X(j(ω−ω0))
于是有:
X p ( j ω ) = 1 T ∑ k = − ∞ + ∞ X ( j ( ω − k ω s ) ) X_p(j\omega) = \frac{1}{T} \sum_{k = -\infty}^{+\infty} X(j(\omega - k\omega_s)) Xp(jω)=T1k=−∞∑+∞X(j(ω−kωs))
即 X p ( j ω ) X_p(j\omega) Xp(jω) 频谱是 X ( j ω ) X(j\omega) X(jω) 频谱以周期 T T T 出现并叠加得到的结果,若频谱无重叠,那么则可以通过低通滤波器得到原信号频谱,这称为 采样定理 :
设 x ( t ) x(t) x(t) 是一个带限信号,在 ∣ ω ∣ > ω M |\omega| > \omega_M ∣ω∣>ωM 的时候 X ( j ω ) = 0 X(j\omega) = 0 X(jω)=0 。如果 ω s > 2 ω M \omega_s > 2\omega_M ωs>2ωM ,那么 x ( t ) x(t) x(t) 就可以被频率为 ω s \omega_s ωs 的周期采样串所采样并且唯一的确定。
其中2倍最大频率 2 ω M 2\omega_M 2ωM 称为 奈奎斯特频率 。
零阶保持采样
实际上传输窄和幅度较大的脉冲是相当困难的,我们可以采用零阶保持的方法来进行采样,即在采样的瞬间,直到下一次采样点来临,都是一直保持这一样本值:
这相当于,冲激采样之后,在串联一个零阶保持器:
零阶保持器具有一个从0开始的矩形窗口的单位脉冲响应,并且保持时间为采样频率 T T T 。
利用插值重建信号
插值,就是用一连续信号对一组样本值的拟合。常见的插值方法有:零阶保持插值、线性插值,多项式插值等。
对于冲激串采样的结果 x p ( t ) x_p(t) xp(t) 来说,通入截止频率为 ω M \omega_M ωM 的低通滤波器即可还原成源信号,这在时域上相当于做卷积:
x r ( t ) = x p ( t ) ∗ h ( t ) = ∫ − ∞ + ∞ h ( t − τ ) ∑ n = − ∞ + ∞ x ( n T ) δ ( τ − n T ) d τ = ∑ n = − ∞ + ∞ x ( n T ) ∫ − ∞ + ∞ h ( t − τ ) δ ( τ − n T ) d τ = ∑ n = − ∞ + ∞ x ( n T ) h ( t − n T ) x_r(t) = x_p(t) \ast h(t) = \int_{-\infty}^{+\infty}h(t - \tau) \sum_{n=-\infty}^{+\infty} x(nT) \delta(\tau - nT) d\tau = \sum_{n=-\infty}^{+\infty} x(nT) \int_{-\infty}^{+\infty}h(t - \tau) \delta(\tau - nT) d\tau = \sum_{n=-\infty}^{+\infty} x(nT) h(t - nT) xr(t)=xp(t)∗h(t)=∫−∞+∞h(t−τ)n=−∞∑+∞x(nT)δ(τ−nT)dτ=n=−∞∑+∞x(nT)∫−∞+∞h(t−τ)δ(τ−nT)dτ=n=−∞∑+∞x(nT)h(t−nT)
上式体现了低通滤波器的冲激响应曲线的拟合插值过程:
对一个理想的低通滤波器来说,其:
h ( t ) = T sin ( ω c t ) π t h(t) =T \frac{\sin(\omega_ct)}{\pi t} h(t)=Tπtsin(ωct)
这种插值方式通常称为 带限插值 。
理想滤波器物理上是不存在的,因此我们通常使用近似的非理想滤波器,会对信号还原造成一定损失,例如零阶保持器,在零阶保持采样中,可以将零阶保持器视为一个低通滤波器,其频率响应和理想滤波器相比如下:
若零阶保持器效果不好,可采用高阶保持器,这里的阶数指,插值后的函数曲线最高在该阶数的导数连续,例如零阶保持器表示在零阶导数都不连续,而线性插值(一阶保持器)表示在零阶导数连续,一阶导数不连续。
易知线性插值的冲激响应是如图的三角波形:
其频率响应和理想滤波器相比如下:
其频率响应为:
H ( j ω ) = 1 T [ sin ( ω T / 2 ) ω / 2 ] 2 H(j\omega) = \frac{1}{T} [\frac{\sin(\omega T/2)}{\omega/2}]^2 H(jω)=T1[ω/2sin(ωT/2)]2
欠采样
当 ω s > 2 ω M \omega_s > 2 \omega_M ωs>2ωM 的时候,采样信号的频谱是由原信号频谱重复组成,当 ω s > 2 ω M \omega_s > 2 \omega_M ωs>2ωM 的时候,采样信号的频谱不再是重复原信号频谱,而是发生 混叠现象 ,因此利用低通滤波器也不能把信号还原出来。
显然,对于采样点上的信号值通常能够准确的还原出来,即:
x r ( n T ) = x ( n T ) x_r(nT) = x(nT) xr(nT)=x(nT)
我们使用 x ( t ) = cos ω 0 t x(t) = \cos \omega_0 t x(t)=cosω0t 来说明欠采样的效果,当采用不同的采样频率的时候:
x r ( t ) = cos ω 0 t = x ( t ) ω 0 = ω s 6 x r ( t ) = cos ω 0 t = x ( t ) ω 0 = 2 ω s 6 x r ( t ) = cos [ ( ω s − ω 0 ) t ] ≠ x ( t ) ω 0 = 4 ω s 6 x r ( t ) = cos [ ( ω s − ω 0 ) t ] ≠ x ( t ) ω 0 = 5 ω s 6 \begin{matrix} x_r(t) = \cos \omega_0 t = x(t) & \omega_0 = \frac{\omega_s}{6} \\ x_r(t) = \cos \omega_0 t = x(t) & \omega_0 = \frac{2\omega_s}{6} \\ x_r(t) = \cos[(\omega_s - \omega_0)t] \neq x(t) & \omega_0 = \frac{4\omega_s}{6} \\ x_r(t) = \cos[(\omega_s - \omega_0)t] \neq x(t) & \omega_0 = \frac{5\omega_s}{6} \\ \end{matrix} xr(t)=cosω0t=x(t)xr(t)=cosω0t=x(t)xr(t)=cos[(ωs−ω0)t]=x(t)xr(t)=cos[(ωs−ω0)t]=x(t)ω0=6ωsω0=62ωsω0=64ωsω0=65ωs
当混叠发生的时候,原始频率 ω 0 \omega_0 ω0 就被混叠成一个较低的频率 ω s − ω 0 \omega_s - \omega_0 ωs−ω0 。当 ω 0 = ω s \omega_0 = \omega_s ω0=ωs 的时候,被重建的信号就是一个常数。
对于带相位的 x ( t ) = cos ( ω 0 t + ϕ ) x(t) = \cos (\omega_0 t + \phi) x(t)=cos(ω0t+ϕ) 来说,傅里叶变换的形状是相同的,只不过在 ω 0 \omega_0 ω0 的位置的冲激有个幅度因子 π e j ϕ \pi e^{j\phi} πejϕ 在 − ω 0 -\omega_0 −ω0 的时候冲激有个幅度因子 π e − j ϕ \pi e^{-j\phi} πe−jϕ ,当发生混叠的时候,此时:
x r ( t ) = cos [ ( ω s − ω 0 ) t − ϕ ] x_r(t) = \cos[(\omega_s - \omega_0)t - \phi] xr(t)=cos[(ωs−ω0)t−ϕ]
这个在相位上有个颠倒,称为 相位倒置 。
最后需要注意一点是,采样定理一定是采样频率大于信号最大频率的2倍,等于是不够的,考虑:
x ( t ) = cos ( ω s 2 t + ϕ ) x(t) = \cos(\frac{\omega_s}{2}t + \phi) x(t)=cos(2ωst+ϕ)
若通过频率 ω s \omega_s ωs 采样,得到的输出为:
x r ( t ) = cos ϕ cos ( ω s 2 t ) x_r(t) = \cos\phi\cos(\frac{\omega_s}{2}t) xr(t)=cosϕcos(2ωst)
可见只有当 2 π 2\pi 2π 倍相位偏移的时候,才能正确还原。