实验五、信号取样和取样定理.docVIP

实验五、信号取样与取样定理 1实验目的 1）学会利用MATLAB完成信号取样及对取样信号的频谱进行分析； 2）学会利用MATLAB改变取样间隔，观察取样后信号的频谱变化； 3）学会利用MATLAB对取样后的信号进行重建。 2实验原理及实例分析 （信号取样及取样定理的基本原理见教材的第4.8节） 采用数值计算法计算傅立叶变换的理论依据是： 若信号为时限信号，当时间间隔T取得足够小时，上式可演变为： 上式用MATLAB表示为： F=f * exp(j * t’ * w) * T 其中F为信号的傅里叶变换，w为频率，T为时间步长。 2.1 信号取样 例1：已知升余弦脉冲信号为，，令参数、，当时，编程实现该信号经过冲激脉冲取样后得到的取样信号及其频谱。 解：MATLAB程序如下： clc close all clear all Ts = 1; dt = 0.01; t1 = -4:dt:4; ft = 1/2 * (1+cos(t1)) .* (uCT(t1+pi) - uCT(t1-pi)); subplot(221);plot(t1,ft,Linewidth,2);grid; axis([-4 4 -0.1 1.1]);xlabel(t(sec));title(f(t)); w = -3*pi :0.01:3*pi; F1 = dt * ft * exp(-j * t1 * w); subplot(222);plot(w,abs(F1),Linewidth,2);grid; xlabel(\omega);title(|F(j\omega)|); axis([-10 10 -0.2 1.1*pi]); t2 = -4:Ts:4; fst = 1/2 * (1+cos(t2)) .* (uCT(t2+pi) - uCT(t2-pi)); subplot(223);plot(t1,ft,:);hold on stem(t2,fst);hold off; grid; axis([-4 4 -0.1 1.1]);xlabel(t(sec));title(f_s(t)); Fsw = Ts * fst * exp(-j * t2 * w); subplot(224);plot(w,abs(Fsw),Linewidth,2);grid; axis([-10 10 -0.2 1.1*pi]); xlabel(\omega);title(|F_s(\omega)|); 产生的图形如图1所示： 图1 例1程序产生的结果 2.2取样定理 例2：用例1来验证取样定理。 解：例1中信号的频谱大部分集中在之间，设其截止频率为，代入参数可得，因而Nyquist间隔为。在例1的程序中，可通过修改的值得到不同的结果。当时，程序运行结果如图2所示。当时，程序运行结果如图3所示，从图3可以看出产生了较为严重的频谱混叠现象。 图2 时的图形 图3 时的图形 2.3信号重建 例3：对例1中的升余弦脉冲信号，假设其截止频率，取样间隔，采用截止频率的低通滤波器对抽样信号滤波后重建信号，并计算重建信号与原升余弦脉冲信号的绝对误差。 解：MATLAB程序如下： clc close all clear all wm = 2; wc = 1.2 * wm; Ts = 1; N = 100; nTs = (-N:N) * Ts; fs = 1/2 * (1+cos(nTs)) .* (uCT(nTs+pi) - uCT(nTs-pi)); t = -4:0.01:4; ftr = zeros(1,length(t)); L = length(-N:N); for n = 1:L m = n - ((L-1) / 2 + 1); xt = fs(n) * sinc((wc/pi) * (t - m * Ts)); ftr = ftr + Ts * wc / pi * xt; end t1 = -4:0.01:4; f1 = 1/2 * (1+cos(t1)) .* (uCT(t1+pi) - uCT(t1-pi)); subplot(311);plot(t1,f1,:);hold on stem(nTs,fs);hold off; grid; axis([-4 4 -0.1 1.3]); xlabel(nTs);title(f(nTs)); subplot(312); plot(t,ftr,Linewidth,2);grid; axis([-4 4 -0.1 1.3]); xlabel(t);title(Reconstructed signal fr(t)); error = abs(ftr - f1); su