实验二：数字滤波器的设计

实验二：数字滤波器的设计

一、实验目的：

掌握IIR数字低通滤波器的设计方法。

二、实验原理：

 

 

 

 

 

% 清空命令窗口
设计代码：clc;
% 设定参数
Fs = 1000; % 采样频率
fp = 100; % 通带截止频率
fs = 300; % 阻带截止频率
rp = 3; % 通带最大衰减
rs = 20; % 阻带最小衰减
% 将频率转换为角频率
wp = 2 * pi * fp / Fs;
ws = 2 * pi * fs / Fs;
% 转换为模拟滤波器的通带和阻带边界频率
Ts = 1 / Fs;
wap = 2 / Ts * tan(wp / 2);
was = 2 / Ts * tan(ws / 2);
% 计算巴特沃斯滤波器阶数和截止频率
[n, wn] = buttord(wap, was, rp, rs, 's');
% 获取巴特沃斯滤波器的零点、极点和增益
[z, p, k] = buttap(n);
% 转换零点-极点-增益表示为分子-分母表示
[bp, ap] = zp2tf(z, p, k);
% 将低通巴特沃斯滤波器转换为所需的频率
[bs, as] = lp2lp(bp, ap, wn);
% 通过双线性变换将模拟滤波器转换为数字滤波器
[bz, az] = bilinear(bs, as, Fs);
% 计算滤波器的频率响应
[h, w] = freqz(bz, az, 1024, Fs);
% 绘制滤波器的振幅频率响应
figure;
subplot(4,1,1);
plot(w, abs(h));

title('低通滤波器');
xlabel('频率'); ylabel('振幅');
% 绘制滤波器的幅度响应的对数尺度图
subplot(4,1,2);
plot(w, 20 * log10(abs(h)));
title('Magnitude Frequency');
xlabel('频率'); ylabel('dB');

(2) 用所设计的滤波器对实际心电图信号采样序列进行仿真滤波处理，并分别打印出滤波前后的心电图信号波形图，改变参数Fs，观察总结滤波作用与效果。

心电图信号采样序列x(n)： 人体心电图信号在测量过程中往往受到工业高频干扰， 所以必须经过低通滤波处理后， 才能作为判断心脏功能的有用信息。 下面给出一实际心电图信号采样序列样本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}

% 原始心电图信号
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];
% 滤波
y = filter(bz, az, x);
% 画出原始心电图信号
subplot(4,1,3);
plot(x);
title('原始心电图信号');
% 画出滤波后心电图信号
subplot(4,1,4);
plot(y);
title('滤波后信号');

 

五、实验总结

1、实验结论

由上图可知，滤波器的通带在100hz时的衰减为1.19dB左右，小于技术要求的3dB；滤波器阻带的通带在300Hz时的衰减为20dB左右，大于技术要求的20dB。所以滤波器基本满足要求。

将滤波器应用到仿真的实际心电图信号中，得到的滤波后的心电图信号如上图中的第4个图，可见经过滤波器后比原始信号更加平滑。

2、心得体会

在此次实验中，我学习如何利用Matlab进行数字滤波处理，了解了如何选择滤波器的参数，包括采样频率、通带截止频率、阻带截止频率等，以及如何根据这些参数计算出滤波器的阶数和截止频率；学会了如何利用 MATLAB 中的函数来实现模拟滤波器到数字滤波器的转换，并且了解了双线性变换的原理和方法；通过计算滤波器的频率响应并绘制振幅频率响应图和幅度响应的对数尺度图，使我对滤波器在频率域上的性能有了更直观的了解；应用设计好的滤波器对心电图信号进行滤波处理，观察了滤波前后信号的变化，掌握了数字滤波在信号处理中的实际应用，进一步丰富和巩固了自己的专业技能。