实验三 带阻滤波抑制心电50HZ工频干扰

实验三带阻滤波抑制心电50HZ工频干扰

生物医学工程系罗融编一、实验目的：

1．学习带阻滤波器的设计方法。

2．学习用MATLAB语言编写带阻滤波器设计程序。

3．学习用带阻滤波器处理生理信号。

二、实验内容：

1．设计IIR带阻滤波器

中心频率f0=50Hz，通带边界频率f p=[48.04 52.04]，阻带边界频率f BS=[49.04

51.04]，通带最大波动 Amax=1dB ，阻带最小衰减Amin=20dB，抽样频率

fs=250Hz。要求绘出幅频特性曲线图与相频特性曲线图。（可用buttord，butter 函数设计IIR带阻滤波器，用freqz绘图）

2．设计FIR带阻滤波器

抽样频率fs=250Hz，阻带边界频率f BS=[49 51]，n=800，用窗函数法设计线性相位FIR带阻滤波器。要求绘出幅频特性曲线图与相频特性曲线图。（可用fir1函数设计FIR带阻滤波器，用freqz绘图）

3．用所设计的IIR带阻滤波器处理含有50Hz干扰的心电信号，要求显示输入输出的时域波形及它们的谱。（可用filter函数滤波，用fft函数计算频

谱，用plot，stem作图）。

4．用所设计的FIR带阻滤波器处理含有50Hz干扰的心电信号，要求显示输入输出的时域波形及它们的谱。（可用filter函数滤波，用fft函数计算频

谱，用plot，stem作图）。含有50HZ干扰的心电信号由数据文件

shiyansan.mat提供，用load shiyansan命令后，shiyansan数据文件中的变

量ecg50即在matlab工作空间中，可用plot（ecg50）语句观察该含有50HZ

干扰的心电信号。

三、报告要求：

报告格式要求同实验一。报告内容应包含实验名称，实验目的，实验内容，实验程序代码及结果，实验结果分析与讨论等，并对两种滤波器的特点进行比较分析。附录：

1）窗函数法：