试验四IIR数字滤波器设计试验目的掌握用双线性变换法设计.doc

实验四 IIR数字滤波器设计 一、实验目的 (1)掌握用双线性变换法设计IIR数字低通和高通滤波器。 (2)设计低通滤波器对实际心电图信号进行滤波。 (3)设计低通滤波器对含有啸叫噪声的音乐信号进行消噪。 *(4)设计IIR数字低通和高通滤波器对某个DTMF(双音多频)信号进行频带分离。 二、实验环境 1.Windows98以上操作系统 2.安装MATLAB6.0以上版本 三、实验原理 1．选频型数字滤波器的种类有 低通 、 高通 、 带通 和 带阻 滤波器。 2. 从实现方法上，数字滤波器通常分为 IIR 和 FIR 滤波器。 3．IIR滤波器的设计目的是根据技术指标，找到 H(z)分子/分母系数b,a ；IIR滤波的MATLAB语句为y= filter(b,a,x) ; 四、实验内容 1. 人体心电图信号在测量过程中往往受到工业高频干扰， 必须经过低通滤波处理后才能作为判断心脏功能的有用信息。给出一实际心电图信号采样序列样本x(n)， 其中存在高频干扰。 试以x(n)作为输入序列， 滤除其中的干扰成分。 x(n)= {-4, -2, 0, -4, -6, -4, -2, -4, -6, -6, -4, -4, -6, -6, -2, 6, 12, 8, 0, -16, -38, -60, -84, -90, -66, -32, -4, -2, -4, 8, 12, 12, 10, 6, 6, 6, 4, 0, 0, 0, 0, 0, -2, -4, 0, 0, 0, -2, -2, 0, 0, -2, -2, -2, -2, 0} 。 低通滤波器 设计指标：ωp=0.2πrad,ωs=0.3πrad,ap=1dB,as=15dB。 已设计出H(z) (p300) 其中 A=0.09036 ；B1=1.2686,C1=-0.7051 ；B2=1.0106,C2=-0.3583；B3=0.9044,C3=-0.2155 *IIR滤波的Matlab语句：y=filter(b,a,x) b,a----Hk(z)分子/分母系数；x---输入信号x(n)；y---滤波结果y(n)。 **求频响特性的Matlab语句：[H,w]=freqz(b,a,N) b,a----Hk(z)分子/分母系数；N---频率点数；H—滤波器的频响特性H(w); w---数字频率轴w。 (1)用IIR低通滤波器(原理设计)对实际的心电信号进行滤波 (IIRECG.m) %程序：(IIRECG.m) clear clc x=[-4,-2,0,-4,-6,-4,-2,-4,-6,-6, ... -4, -4, -6, -6,-2,6,12,8,0,-16, ... -38,-60,-84,-90,-66,-32,-4, -2, -4, 8, ... 12, 12, 10, 6, 6, 6, 4, 0, 0, 0, ... 0, 0, -2, -4, 0, 0, 0, -2, -2, 0, ... 0, -2, -2, -2, -2, 0];%心电信号x(n)[含高频噪声] A=0.09036; %IIR低通滤波器系数 B1=1.2686;C1=-0.7051; B2=1.0106;C2=-0.3583; B3=0.9044;C3=-0.2155; b=[A 2*A A]; a1=[1 -B1 -C1]; N=128; [H1,w]=freqz(b,a1,N);%频响特性H1(w) y1=filter(b,a1,x); %y1(n)是x(n)经H1(z)滤波的结果 a2=[1 -B2 -C2]; [H2,w]=freqz(b,a2,N);%频响特性H2(w) y2=filter(b,a2,y1); %y2(n)是y1(n)经H2(z)滤波的结果 a3=[1 -B3 -C3]; [H3,w]=freqz(b,a3,N);%频响特性H3(w) H=H1.*H2.*H3; %总的频响特性H(w)=H1(w)H2(w)H3(w) mag=abs(H); %滤波器的幅频特性 db=20*log10((mag+eps)/max(mag));%幅频特性(dB) y3=filter(b,a3,y2); %y3(n)是y2(n)经H3(z)滤波的结果 X=fft(x,N); %对x(n)做N点fft,得其频