资源描述：
《用窗口法设计FIR数字滤波器.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、信号、系统与信号处理实验Ⅱ实验报告实验名称：用窗口法设计FIR数字滤波器一、实验目的了解一个实际滤波器的设计过程，加深掌握用窗口法设计FIR数字滤波器的原理和窗函数对数字滤波器性能的影响二、实验内容与要求1．编写用改进的升余弦窗设计FIR线性相应低通数字滤波器的程序,已知=wc=0.5p，N=21。这里wcw为低通滤波器的截止频率，N为滤波器的长度,滤波器的阶数为N-1。2．调试运行程序,要求在屏幕上显示出单位脉冲响应h(n)的数值。画出其幅度响应
2、H(ejw)
3、及20log10
4、H(ejw)
5、的曲线。3.画出窗函数w(n)及其频谱
6、W(ejw)
7、和20log10
8、W(e
9、jw)
10、的曲线。三、实验程序与结果clearall;N=21;wc=0.5*pi;b=fir1(N-1,wc/pi,hamming(N))freqz(b,1,512);h(n)的数值即b：b=Columns1through110.00000.0036-0.0000-0.01220.00000.0343-0.0000-0.08580.00000.31060.4991Columns12through210.31060.0000-0.0858-0.00000.03430.0000-0.0122-0.00000.00360.0000clearall;N=21;wc=0.5*pi;
11、r=(N-1)/2;n=0:N-1;hdn=sin(wc*(n-r))/pi./(n-r);ifrem(N,2)~=0hdn(r+1)=wc/pi;endwn=hamming(N);h=hdn.*wn';H=fft(h,512);W=fft(wn,512);w=2*[0:511]/512;figure;subplot(3,1,1);plot(wn);subplot(3,1,2);plot(w,abs(W));ylabel('
12、W(eiw)
13、')subplot(3,1,3);plot(w,20*log10(abs(W)));ylabel('20lg
14、W(eiw)
15、');fi
16、gure;subplot(2,1,1);plot(w,abs(H));ylabel('
17、H(eiw)
18、')subplot(2,1,2);plot(w,20*log10(abs(H)));ylabel('20lg
19、H(eiw)
20、');一、仿真结果分析图一是FIR低通滤波器的幅频和相频特性图图二是滤波器设计过之后的幅度特性图三是窗函数的幅度特性FIR滤波器的设计过程在于运用海明窗口进行设计，按照定义h(n)=hd(n)*w(n)出发来求出w(n)的值。五、实验问题解答与体会思考题（1）如何根据20log10
21、H(e^jw)
22、曲线得到阻带最小衰减？最小阻带衰减As=20lg(
23、
24、H(e^j0)
25、/
26、H(e^jws)
27、)。所以在20log10
28、H(e^jw)
29、曲线中，w=0和w=1所对应的20lg
30、H(e^jw)
31、值之差就是最小阻带衰减。（2）如何根据20log10
32、W(eiw)
33、曲线得到旁瓣峰值幅度和主瓣宽度？20lg
34、W(e^jw)
35、曲线中第一个波谷对应w值的两倍就是主瓣宽度；旁瓣峰值幅度定义为窗函数的幅频函数
36、W(e^jw)
37、的最大旁瓣相对主瓣最大值的衰减。从20lg
38、W(e^jw)
39、曲线看即为w=0后开始的第一个波峰所对应的20lg
40、W(e^jw)
41、的值减去w=0的20lg
42、W(e^jw)
43、的值就是旁瓣峰值幅度
