信号的时域抽样与重建
离散时间信号大多数由连续时间信号进行抽样获的。在信号进行处理的过程中,要使有限带宽信号被抽样后能够不失真地还原出原模拟信号,抽样信号的周期Ts及抽样频率Fs的取值必须符合奈奎斯特定理。
如果Fs的取值大于两倍的信号最高频率fm,只要经过一个低通滤波器,抽样信号就能不失真地还原出原模拟信号。反之,则频谱将发生混叠,抽样信号将无法不失真地还原出原模拟信号。
实例
1、对连续信号进行采样
已知一个连续时间信号,其中:f0=1HZ,取最高有限带宽频率fm=5f0。分别显示原连续时间信号波形和3种情况下抽样信号的波形。(Fs=fm、Fs=2fm、Fs=3fm)
f0 = 1;
T = 1/f0;
fm = 5*f0;
t = 0:0.001:5;
x = sin(2*pi*f0*t) + (1/3)*sin(6*pi*f0*t);
subplot(411);plot(t,x);
for i=1:3;
Fs = i*fm;
Ts = 1/Fs;
n = 0:Ts:5; %离散信号的抽样周期
xs = sin(2*pi*f0*n) + (1/3)*sin(6*pi*f0*n);
subplot(4,1,1+i);stem(n,xs,'filled');
axis([min(n),max(n),1.1*min(xs),1.1*max(xs)]);
end
2、连续信号和抽样信号的频谱
根据理论可知,信号的频谱图可以直观地反映出抽样信号能否恢复原模拟信号波形。通过对上例中3种情况求振幅频谱,来进一步分析和证明时域抽样定理。
f0=1;T=1/f0;
t=-2:0.1:2;
N=length(t);
x=sin(2*pi*f0*t)+1/3*sin(6*pi*f0*t);
fm=5*f0;Tm=1/fm;
wm=2*pi*fm;
k=0:N-1;w1=k*wm/N;
X=x*exp(-j*t'*w1)*Tm;
subplot(4,1,1);plot(w1/(2*pi),abs(X))
axis([0,max(4*fm),1.1*min(abs(X)),1.1*max(abs(X))]);
for i=1:3;
if i<=2 c=0;else c=1; end
fs=(i+c)*fm;Ts=1/fs;
n=-2:Ts:2;
xs=sin(2*pi*f0*n)+1/3*sin(6*pi*f0*n);
N=length(n);
ws=2*pi*fs;
k=0:N-1;w=k*ws/N;
Xs=xs*exp(-j*n'*w)*Ts;
subplot(4,1,1+i);plot(w/(2*pi),abs(Xs));
axis([0,max(4*fm),1.1*min(abs(Xs)),1.1*max(abs(Xs))]);
end
3、利用时域卷积重建例题信号
时域卷积是用时域抽样信号xs(t)与理想滤波器系统的单位冲激响应h(t)进行卷积积分来求解。卷积积分的公式通过推导化简为内插公式
f0=1;T0=1/f0; %基波频率、周期
fm=5*f0;Tm=1/fm; %最高频率为基波的5倍
t=0:0.01:3*T0;
x=sin(2*pi*f0*t)+1/3*sin(6*pi*f0*t);
subplot(4,1,1);plot(t,x)
for i=1:3;
fs=i*fm;Ts=1/fs; %设置不同采样频率、周期
n=0:(3*T0)/Ts;
t1=0:Ts:3*T0;
xs=sin(2*pi*f0/fs*n)+1/3*sin(6*pi*f0/fs*n); %抽样信号
T_N=ones(length(n),1)*t1-n'*Ts*ones(1,length(t1));
x1=xs*sinc(2*pi*fs*T_N); %内插公式计算
subplot(4,1,1+i);plot(t1,x1)
axis([min(t1),max(t1),1.1*min(x1),1.1*max(x1)])
end