实验八(实验课)-论文指导设计-.pptxVIP

实验八

IIR和FIR滤波器过滤信号的实现及比较

一、实验目的

■

(1)

掌握数字滤波器的计算机仿真方法。

■

(2)

通过观察对实际心电图信号的滤波作用，

获得数字滤波的感性知识。

二、实验原理与方法

若h(n)是因果FIR滤波器的单位脉冲响应，其长度为N，当输入为

x(n)时，输出序列y(n)可表示为：

取h(n)为低通滤波器N=21截止频率ωc=0.2π，采用海明窗设计出h(n)

若H(z)是IIR滤波器，在通常内频率低于0.2π时，最大衰减小于1dB；在阻带内[0.3π,π]频率区间上，最小衰减大于15dB。用双线性变换法设计得到Butterworth滤波器系统函数H(z)为：

二、实验原理与方法

式中:

滤波器H(z)由三个二阶滤波器H1(z)、H2(z)和H3(z)级联组成，如下图所示：

二、实验原理与方法

可得：

当k=1时，

实验用MATLAB语言工具函数简介

利用IIR滤波器或FIR滤波器对输入信号x(n)进行滤波时，对于有限长脉冲响应滤波器（FIR）,若已知单位脉冲响应h(n),根据滤波器输出和输入之间的关系y(n)=x(n)*h(n)，滤波器的输出可利用函数conv()直接卷积得到（函数的用法见实验二）。

若已知滤波器的传递函数H(z)，在MATLAB的数字信号处理工具箱

中提供了函数filter(),该函数利用数字滤波器对输入信号x(n)进行滤波时其实现采用直接Ⅱ型结构，适用于IIR和FIR两种滤波器。

实验用MATLAB语言工具函数简

设滤介波器的系统传递函数：

函数调用格式：

y=filter(b,a,x)

[y,zf]=filter(b,a,x)

y=filter(b,a,x,zi)

其中格式[y,zf]=filter(b,a,x)除了得到矢量y外，

还可以得到x的最终状态矢量zf；

格式y=filter(b,a,x,zi)可以在zi中指定x的初始状态。

即：滤波器系数

，输入序列矢量为x。

三、实验内容及要求

心电图信号采样序列x(n)

人体心电图信号在测量过程中往往受到工业高频干扰，所以必须

经过低通滤波处理后，才能作为判断心脏功能有用信息。

下面给出一实际心电图信号采样序列样本x(n)，其中存在高频干扰

三、实验内容及要求

编写FIR滤波器仿真程序，计算其对心电图信号采样序列x(n)的响应序列y1(n)。

编写IIR滤波器仿真程序，计算其对心电图信号采样序列x(n)的响应序列y2(n)。

在通用计算机上运行仿真滤波器程序，在屏幕上打印出x(n)

y1(n)和y2(n)，并进行比较。

四、实验报告要求

(1)对比滤波前后的心电图信号波形，说明数字滤波器的滤波作用。

(2)比较FIR和IIR的滤波性能，运算量。

五、相关程序：

xn=[4-20-4-6...

-4-2-4-6-6...

-4-4-6-6-2...

61280-16...

-38-60-84-90-66...

-32-4-2-4-8...

12121066...

64000...

00-2-40...

00-2-20...

0-2-2-2-20]

生成序列x(n)

五、相关程序：

利用窗函数法设计FIR系统

hn=fir1(N,Wc/pi,hamming(N+1));

经过三个系统H(z)合为一个大系统为

H4=H1.*(H2);

H=H4.*(H3);

利用H(z)设计IIR系统

A=0.2318;

b1=[A,2*A,A];

a1=[1,-0.4164,0.6006];

y1=filter(b1,a1,xn);

五、相关程序：

xk=fft(xn,64);

yk=fft(yn,64);

[h,w]=freqz(hn,1,512);

求H(ejω)

FIR

IIR

求输出y(n)■

FIR

IIR

[H1,W]=freqz(b1,a1,512);

y=filter(b,a,x)

y=filter(hn,1,xn);

yn=filter(bz,az,xn);

求输入和输出的幅度曲线

可利用FFT

要求：

设计FIR滤波器和IIR滤波器分别处理心电图序列

FIR滤波器

N=21截止频率ωc=0.2π，采用海明窗设计出h(n)

IIR滤波器

要求：

画出：

FIR和IIR两种滤波器的系统幅度曲线图；