FIR滤波器设计实验：不同窗函数与频率抽样法对比.pdfVIP

实验三FIR滤波器的设计

一．实验目的：

研究FIR滤波器设计的有关问题

二．实验要求：

（1）当N=11、41、81、121时，用矩形窗设计一个低通FIR滤波器，其剪切频率w=0.3pi。

画出其幅度响应H（w），用哈明窗重复上述设计，比较两个窗。

（2）用凯泽尔窗，哈明窗，矩形窗设计具有指定指标的滤波器。

（3）用频率抽样法设计一个理想低通滤波器，通带截止频率为0.375pi，阻带截止频率为

0.5pi。分别为不加，加一个，加两个过渡点。

三．实验程序：

1．设计FIR滤波器：

functiondsp203()

closeall;

N=11;

wc=0.3;

B=fir1(N,wc,boxcar(N+1));

[Hz,wz]=freqz(B,1);

Hz=20*log10(abs(Hz));

subplot(2,2,1)

plot(wz/pi,Hz);

gridon

title(矩形窗函数法幅频特性N=11);

N=41;

wc=0.3;

B=fir1(N,wc,boxcar(N+1));

[Hz,wz]=freqz(B,1);

Hz=20*log10(abs(Hz));

subplot(2,2,2);

plot(wz/pi,Hz);

title(矩形窗函数法幅频特性N=41);

gridon

N=81;

wc=0.3;

B=fir1(N,wc,boxcar(N+1));

[Hz,wz]=freqz(B,1);

Hz=20*log10(abs(Hz));

subplot(2,2,3)

plot(wz/pi,Hz);

title(矩形窗函数法幅频特性N=81);

gridon

N=121;

wc=0.3;

B=fir1(N,wc,boxcar(N+1));

[Hz,wz]=freqz(B,1);

Hz=20*log10(abs(Hz));

subplot(2,2,4)

plot(wz/pi,Hz);

title(矩形窗函数法幅频特性N=121);

gridon

2．设计指定指标的滤波器：

functiondsp204()

closeall;

wc=0.45;

dw=0.1;

N=100;

F=[0,0.4,0.5,1];A=[1,1,0.0001,0];

B=fir2(N,F,A,rectwin(N+1));

[Hz,wz]=freqz(B,1);

Hz=20*log10(abs(Hz));

subplot(3,1,1)

plot(wz/pi,Hz);

title(窗函数法设计低通滤波器rectanglewindows);

gridon

N=100;a=10;

B=fir2(N,F,A,kaiser(N+1,a));

[Hz,wz]=freqz(B,1);

Hz=20*log10(abs(Hz));

subplot(3,1,2)

plot(wz/pi,Hz);

title(窗函数法设计低通滤波器kaiserwindows);

gridon

N=100;

B=fir2(N,F,A,hamming(N+1));

[Hz,wz]=freqz(B,1);

Hz=20*log10(abs(Hz));

subplot(3,1,3)

plot(wz/pi,Hz);

title(窗函数法设计低通滤波器hammingwindows);

gridon

F=[0,0.4,0.5,1];A=[0,0,1,1];

N=100;

figure;

B=fir2(N,F,A,rectwin(N+1));

[Hz,wz]=freqz(B,1);

Hz=20*log10(abs(Hz));

subplot(3,1,1)

plot(wz/pi,Hz);

title(窗函数法设计高通滤波器rectanglewindows);

gridon

N=100;a=10;

B=fir2(N,F,A,kaiser(N+1,a));

[Hz,wz]=freqz(B,1);

Hz=20*log10(abs(Hz));

subplot(3,1,2);

plot(wz/pi,Hz);

title(窗函数法设计高通滤波器kaiserwindows);

gridon