用窗口法设计FIR数字滤波器(杭电).pdfVIP

用窗口法设计FIR数字滤波器

（杭电）

信号、系统与信号处理实验H

实验报告

实验名称：用窗口法设计FIR数字滤波器

一、实验目的

了解一个实际滤波器设计过程，加深掌握

用窗口法设计FIR数字滤波器的原来和窗口函数

对数字滤波器性能的影响。

二、实验内容与要求

（）编写用改进的升余弦窗设计线性

1FIR

相应低通数字滤波器的程序，已知，，这

0.5N21

c

里为低通滤波器的截止频率，为滤波器的长度，

N

c

滤波器的阶数为。

N-1

（）调试运行程序，要求在屏幕上显示出

2

单位冲激响应的数值，画出其幅度响应及

hn|Hejw20lg

的曲线。

Hejw

（）画出窗口函数的幅度及的曲

3jwjw

nWe20lgWe

线。

思考：

（）如何根据曲线得到阻带最小衰

120lgHejw

减？

（）如何根据曲线得到旁瓣峰值幅

220lgWejw

度和主瓣宽度?

三、实验程序与结果

clc;

clearall;

N=21;

wc=0.5*pi;

r=(N-1)/2;

n=0:N-1;

hdn=sin(wc*(n-r))/pi./(n-r);ifrem(N,2)~=0;

hdn(r+1)=wc/pi;

end

wn=hamming(N);h=hdn.*wn;

H=fft(h,512);

w=2*[0:511]/512;

subplot(3,1,1);

stem(n,h);

title(h(n));

xlabel(n);

ylabel(幅度);subplot(3,1,2);plot(w,abs(H));

title(H(exp(jw)));xlabel(w/pi);

ylabel(幅度);subplot(3,1,3);

plot(w,20*log10(abs(H)));title(H(exp(jw)));

幅度()

xlabel(w/pi);ylabel(