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

2 views (last 30 days)
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

Answers (0)

Categories

Find more on ECG / EKG in Help Center and File Exchange

Products

Community Treasure Hunt

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

Start Hunting!