新·MATLAB实验.ppt

figure(2); subplot(2,1,1); plot(w/pi,db); axis([0,1,-50,10]); title(Designed Lowpass Filter Magnitude Response in dB) xlabel(frequency in pi units); ylabel(|H[exp(jw)]|); grid subplot(2,1,2); plot(w/pi,pha,r);grid title(Designed Lowpass Filter Phase Response in radians); xlabel(frequency in pi units);ylabel(arg(H[exp(jw)]); grid figure(3); subplot(2,1,1); plot(n/fs,x); title(Souce signal:x=sin(2*pi*100*t)+sin(2*pi*600*t) ); xlabel(time in s); ylabel(x); grid subplot(2,1,2); plot(n/fs,y); title(Filtered signal:y); xlabel(time in s); ylabel(y); grid 一、实验目的 1、熟悉窗函数法、频率取样法设计IIR数字滤波器的方法。 2、观察对实际正弦组合信号的滤波作用。 二、实验内容及要求 1、分别编制采用矩形窗、三角窗、汉宁窗、海明窗、布拉克曼窗、凯泽窗设计线性相位低通FIR数字滤波器的程序。 上机实验四：FIR低通数字滤波器的设计 要求的低通FIR数字滤波器指标如下： 通带截止频率 ，阻带起始频率 ，阻带最小衰减 As=50dB，求出滤波器的单位取样响应，幅频和相频响应，绘出它们的图形，并比较各种窗函数和窗口长度对滤波性能的影响。 2、编制采用频率取样法设计一个20阶的线性相位低通FIR数字滤波器的程序。 要求通带截止频率 ，阻带起始频率 。求出滤波器的单位取样响应，幅频和相频响应，绘出它们的图形，并比较阶数对滤波性能的影响。过渡带抽样点增加 (a)一点：T=0.5； (b)二点：T=0.66,T=0.33；比较对滤波性能的影响。 3、通过上述滤波器对以下正弦组合信号进行滤波，抽样频率为2KHz，绘出原始信号和滤波后的信号曲线，并比较滤波效果。 x(t)=sin(2*pi*100*t)+sin(2*pi*800*t) % Windows_Method.m % Linear Phase Lowpass Filter Design - Hamming window % Digital Filter Specifications: % wp = 0.2*pi; ws = 0.4*pi; wc=0.3*pi, (Rp=0.25dB), As=50dB % x(t)=sin(2*pi*100*t)+sin(2*pi*800*t),fs=2000; wp = 0.2*pi; ws = 0.4*pi; As=50; tr_width = ws - wp % M = ceil(1.8*pi/tr_width) + 1 % boxcar;As=21dB % M = ceil(6.1*pi/tr_width) + 1 % triang;As=25dB % M = ceil(6.2*pi/tr_width) + 1 % hanning;As=44dB M = ceil(6.6*pi/tr_width) + 1 % hamming;As=53dB % M = ceil(11*pi/tr_width) + 1 % blackman;As=74dB % M = ceil((As-7.95)/(14.36*tr_width/(2*pi))+1) + 1 % kaiser; % beta = 0.1102*(As-8.7) % kaiser;As50dB % beta = 0.584*(As-21)^0.4+0.07886*(As-21) % kaiser;21As50dB n=[0:1:M-1]; wc = (ws+wp)/2 % Ideal LowPass filter computation alpha = (M-1)/2; m = n - alpha + eps; hd = sin(wc*m) .