1 Introduction

This paper relates to a method for detecting the peak value of an electrocardiogram waveform, comprising the following steps: preprocessing the electrocardiogram waveform; calculating the cosine value of the sampling point of the electrocardiogram waveform to obtain a triangular cosine sequence; detecting the peak position according to the triangular cosine sequence; The position is corrected and confirmed, and points that are not peaks are excluded. This article also involves a peak detection system for ECG waveforms. This article can quickly and accurately detect the peak position of the ECG waveform, especially the QRS wave, and effectively reduce the workload of doctors.​

2 part code

%% 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

                       %from original signal                  

plot(rr4,rr5,'r*') %plot  R-R(n)(X-Axis) vs R-R(n-1)(Y-Axis)

 title('POINCARE PLOT'), xlabel('RR(n+1)') ,ylabel('RR(n)')

 %% Task 4-b

ki=length(rr3) ;

ahr=mean(hr);       % mean heart rate;

disp(['mean hrv = ' num2str(ahr)]); 

% disp is used to display the value(s);

SDNN = std(rr3); 

% SDNN, standard deviation for RR interval used in statical analysis;

disp(['SDNN = ' num2str(SDNN)]);

sq = diff(rr3).^2;

rms = sqrt(mean(sq)); % RMSSD,

disp(['RMSSD = ' num2str(rms)]);  

% RMS difference for RR interval used in statical analysis;

NN50 = sum(abs(diff(rr3))>.05); 

% NN50 no. of pairs of RR that is more than 50, used in statical analysis;

disp(['NN50 = ' num2str(NN50)]);

%% Try 4/3

%figure,

%comet(ecg_smooth(1:1000))

3 Simulation results

4 References

[1] Zhang Gaodeng, Dong Jun. Electrocardiogram waveform peak detection method and system:, CN104939821A[P].

Introduction to the blogger: Good at Matlab simulation in various fields such as intelligent optimization algorithm, neural network prediction, signal processing, cellular automata, image processing, path planning, drone, etc. You can communicate with Matlab code in private.

Some theories cite online literature, if there is any infringement, contact the blogger to delete it.