心电图峰值检测附Matlab代码"/>
【心电信号】心电图峰值检测附Matlab代码
1 简介
本文涉及一种心电图波形峰值检测方法,包括如下步骤:对心电图波形进行预处理;计算所述心电图波形采样点的余弦值,得到三角余弦序列;根据所述三角余弦序列检测波峰位置;对上述波峰位置进行校正确认,排除不是峰值的点.本文还涉及一种心电图波形峰值检测系统.本文能够快速准确地检测心电图波形特别是QRS波的峰值位置,有效降低医生的工作负荷.
2 部分代码
%% SOLUTION 4:
clc;
clear all;
close all;
ecg=load('12-2-2016 17.53.13.txt');
f_s=250;
N=length(ecg);
t=[0:N-1]/f_s; %time period(total sample/Fs )
figure
plot(t,ecg,'r'); title('Raw ECG Data plotting ')
xlabel('time')
ylabel('amplitude')
w=50/(250/2);
bw=w;
[num,den]=iirnotch(w,bw); % notch filter implementation
ecg_notch=filter(num,den,ecg);
[e,f]=wavedec(ecg_notch,10,'db6');% Wavelet implementation
g=wrcoef('a',e,f,'db6',8);
ecg_wave=ecg_notch-g; % subtracting 10th level aproximation signal
更多推荐
【心电信号】心电图峰值检测附Matlab代码
发布评论