Interpretation of a pseudocode for ECG signal compression using SQ segment compression

조회 수: 2 (최근 30일)
I've been trying to implement the algo mentioned in paper entitled as ECG signal compression with reconstruction via cubic spline approximation by Olga Malgina, Jana Milenkovi, Emil Plesnik, Matej Zajc, Andrej Košir and Jurij F. Tasi. http://ieeexplore.ieee.org/xpl/login.jsp?tp=&arnumber=5752122&url=http%3A%2F%2Fieeexplore.ieee.org%2Fiel5%2F5746576%2F5751960%2F05752122.pdf%3Farnumber%3D5752122
But couldn't interpret it properly... Its mentioned on section 3.2. Method 2 (SQ segment compression). Please follow the link to access the image showing algorithm https://docs.google.com/file/d/0Bxet-v6Iwb_xaXluQzc3eTh3R1k/edit?usp=sharing
here is the code I've tried to implement, I'm able to detect R peak but after that ain't able to do so
if true
% code
load('Sig100ch1.mat')
In_Sig = ans;
N = length(In_Sig);
% Phase space construction
x = In_Sig;
for k = 1:N-1
y(k) = x(k+1) - x(k);
end
% R Peak Detection
for k = 2:N-1
if ((x(k+1) - x(k))*(x(k) - x(k-1)) < 0) && x(k) > 0.65
R(k) = k;
end
end
R = R(R~=0);
% S-Peak detection
for k = 1:length(R) - 1
jj = (R(k)):(R(k)+(R(k+1)-R(k))/2);
Xj = x(jj);
Yj = y(jj);
for ii = 1:length(jj)
D(ii,k) = sqrt((x(k) - Xj(ii))^2 + (y(k) - Yj(ii))^2);
end
[~,ind] = max(D);
S(k) = jj(ind);
end
end
The problem is I'm not able to comprehend it after R peak detection.
I've submitted the Amplitude threshold method which is also used in SQ segment compression and its working fine. which can be accessed here in http://www.mathworks.in/matlabcentral/fileexchange/43451-ecg-signal-compression

답변 (0개)

카테고리

Help CenterFile Exchange에서 ECG / EKG에 대해 자세히 알아보기

제품

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by