数字信号处理实验报告--滤波器对数字信号的影响.docVIP

《数字信号处理》 实验报告 实 验 名 称： 滤波器对数字信号的影响 专 业 班 级 ： 学 生 姓 名 ： 学 号 ： 指 导 教 师 ： 实 验 时 间： 2011年6月15日 1 实验目的： 证明各式滤波器对包括不同频率信号的响应 2 运行环境： PC机，windows2000及其以上，matlab软件。 3 实验设计： Coeff_LF = 1; Coeff_MF = 0; Coeff_HF = 1; Coeff_NOISE = 0; F1 = 350; F2 = 600; F3 = 450; FS = 2000; EndTime = 1/FS * 1023; Time = 0:1/FS:EndTime; N = length(Time) 4 实验结果 1 经过滤后信号的频域形成一条直线 2带阻信号经过滤后时域和频域的波形图， 3 在时域信号BS滤波器的波形图 5 结论 当Rs的值越大，信号失真越小，越精确。说明滤波器的性能越好。 6 心得体会: 本次课程设计中，本次课程设计中，通过实验对各种滤波器的认识进步加强，提高了动手能力。并提高了对MATLAB的认识，进一步了解了其功能的强大和运用和范围的广泛。 7 实验代码： %实验设计程序如下 function FilterSignalDemo Coeff_LF = 1; % Coefficient of low frequency signal F1(350 Hz) Coeff_MF = 0; % Coefficient of middle frequency signal F3(450 Hz) Coeff_HF = 1; % Coefficient of High frequency signal F2(600 Hz) Coeff_NOISE = 0; % Coefficient of noise, such as 0(without noise), 1, sqrt(2), 3 and so on. F1 = 350; % Hz F2 = 600; % Hz F3 = 450; % Hz FS = 2000; % Hz sample frequency YBP = filter(b, a, SinSignal); figure subplot(2,1,1) plot(Time, YBP) xlim([0 0.05]) xlabel(Time) title(Time domain signals after BP filter) FFTSinSignalSum = 2 / N * fft(YBP) subplot(2,1,2) FnSinSum = abs(FFTSinSignalSum); plot(FW(2:length(FW))*FS/pi/2, FnSinSum(2:length(FW))) hold on plot((FW(2:length(FW))-2*pi)*FS/pi/2, FnSinSum(2:length(FW))) xlim([-FS/2 FS/2]) ylim([0 1]) xlabel(Frequency) %-------------------BandStop Filter--------------------------- figure Wp = [400 500]/1000; Ws = [420 480]/1000; Rp = 3; Rs = 30; [n,Wn] = cheb1ord(Wp,Ws,Rp,Rs) [b,a] = cheby1(n,Rp,Wn,stop); freqz(b,a,512,2000); title(Chebyshev Type I BandStop Filter) YBS = filter(b, a, SinSignal); figure subplot(2,1,1) plot(Time, YBS) xlim([0 0.05]) xlabel(Time) title(Time domain signals after BS filter) FFTSinSignalSum = 2 / N * fft(YBS) subplot(2,1,2) FnSinSum = abs(FFTSinSignalSum); plot(FW(2:length(FW))*FS/pi/2, FnSinSum(2:length