窗函数法设计FIR数字滤波器.pdf

数字信号处理实验报告

实验4窗函数法设计FIR数字滤波器

一、实验目的

1．了解常用的几种窗函数，能正确选择适当的窗函数进行滤波器设计；

2．掌握窗函数法设计数字低通滤波器。

二、实验原理

1．常用的窗函数：

矩形窗函数为boxcar和rectwin，调用格式：

w=boxcar（N）w=rectwin(N)

其中N是窗函数的长度，返回值w是一个N阶的向量。

三角窗函数为triang，调用格式：

w=triang(N)

汉宁窗函数为hann，调用格式：

w=hann（N）

海明窗函数为hamming，调用格式：

w=hamming（N）

2．各个窗函数的性能比较

窗函第一瓣相主瓣宽阻带最

数对于主瓣小衰减

衰减（dB）

矩形-1321

窗

三角-2525

窗

汉宁-3144

窗

海明-4153

窗

三、实验内容

题一：生成四种窗函数：矩形窗、三角窗、汉宁窗、海明窗，并观察其

频率响应。

题二：根据下列技术指标，设计一个FIR数字低通滤波器：wp=0.2π，

ws=0.4π，ap=0.25dB，as=50dB，选择一个适当的窗函数，确定单位冲激

响应，绘出所设计的滤波器的幅度响应。

四、上机程序及运行结果

题一：n=30;

%矩形窗及其频响

window1=rectwin(n);

[h1,w1]=freqz(window1,1);

subplot(4,2,1);

stem(window1);title(矩形窗);subplot(4,2,2);

plot(w1/pi,20*log(abs(h1))/abs(h1(1)));title(矩形窗频响);

%三角窗及其频响

window2=triang(n);

[h2,w2]=freqz(window2,1);

subplot(4,2,3);stem(window2);title(三角窗);

subplot(4,2,4);

plot(w2/pi,20*log(abs(h2))/abs(h2(1)));title(三角窗频响);

%汉宁窗及其频响

window3=hann(n);

[h3,w3]=freqz(window3,1);

subplot(4,2,5);stem(window3);title(汉宁窗);

subplot(4,2,6);

plot(w3/pi,20*log(abs(h3))/abs(h3(1)));title(汉宁窗频响);

%海明窗频响

window4=hamming(n);

[h4,w4]=freqz(window4,1);

subplot(4,2,7);stem(window4);title(海明窗);

subplot(4,2,8);

plot(w4/pi,20*log(abs(h4))/abs(h4(1)));title(海明窗频响);

题二：

%理想

低通滤