信号与系统matlab学习（五）

信号的时域抽样与重建

离散时间信号大多数由连续时间信号进行抽样获的。在信号进行处理的过程中，要使有限带宽信号被抽样后能够不失真地还原出原模拟信号，抽样信号的周期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