【心电信号】心电图峰值检测附Matlab代码

心电图峰值检测附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