汉宁窗实现fir理想低通滤波器.pdfVIP

数字信号处理

实验报告

实验一

一、实验目的

通过汉宁窗函数设计低通的FIR滤波器，使学生掌握FIR滤波器的工作原理和汉宁窗函数

对设计的影响。

二、实验代码

M=31;

alpha=(M-1)/2;

n0=0:M-1;

hd=(sin((pi/2)*(n0-alpha))./(pi*(n0-alpha)));

figure(1)

subplot(3,1,1);

stem(n0,hd)

title(理想脉冲响应);

%axis([-1M01.2]);

ylabel(hd(n));

%text(M+1,-1.2,n);

window=hann(M);

subplot(3,1,2);

stem(window);

title(汉宁窗);

h=hd.*(window);

subplot(3,1,3);

stem(n0,h);

[H,m]=freqz(window,1);

figure(2)

subplot(2,1,1);

Hf=abs(H);

Hx=angle(H);

plot(m,Hf);

title(离散幅频特性曲线)

subplot(2,1,2);

plot(m,Hx);

title(离散相频特性曲线)

（函数补充：matlab提供了专门用于求离散系统频响特性的幻术freqz（），调用freqz（）

格式如下：[H,w]=freqz(B,A,N).B和A分别为离散系统的系统函数分子分母多项式的系统向

量，N为正整数，返回量H则包含了离散系统频响在0——pi范围内N个频率等分点的值，向量w

则包含范围内N个频率等分点。调用中若N默认，默认值为512.）

三、实验结果

理想脉冲响应

0.5

)

n0

(

d

h

-0.5

051015202530

汉宁窗

1

0.5

0

05101520253035

0.5

0

-0.5

051015202530

离散幅频特性曲线

15

10

5

0

00.511.522.533.5

离散相频特性曲线

4

2

0

-2

-4

00.511.522.533.5

四、实验流程图