心电图(ECG)是临床最常用的检查之一, 是利用心电图机从体表记录心脏每一心动周期所产生的电活动变化图形的技术。根据获得的心电图测试数据文件”ECGrawdata.txt”,其中数据第一列为采样时间,第二列和第三列是对同一种心电信号进行不同的导联方式所得到的采样数据。试完成: (1)截取心电信号任意4个周期数据,画出序列x1(n)波形,在图中标出采样频率; (2)对x1(n)加入工频干扰(频率50Hz的正弦信号),画出加入噪声的心电信号序列x2(n)波形; (3)已知M点平滑滤波器为 ,当M分别为3、5、7时,对信号x2(n) 进行滤波,画出对应的输出y (n),分析M的大小对滤波器的影响; (4)对于三个不同的M值,分别画出滤波器对应的幅频特性曲线。
首先导入所需要的心电图数据
%% 导入ECGrawdata.txt文件中的数据
val = importdata('ECGrawdata.txt');
% 提取ECGrawdata.txt文件中的数据
data = val.data;
data1 = data(:,1);
data2 = data(:,2);
% 采样时间
t = 0:0.004:20.476;
figure(1);
subplot(2,1,1)
plot(t,data1,'.');
title('心电图信号的时