Signal segmentation error: how to select a specific number of points before and after the ECG signal?

Question

Ioana Cretu 2021-2-7

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/738892-signal-segmentation-error-how-to-select-a-specific-number-of-points-before-and-after-the-ecg-signal

评论： Mathieu NOE 2021-2-23

I want to segment an ECG signal in heartbeats by taking a number of sample points before and after each R-peak. I did the following loop:

before=100;
after=140;
for i=1:length(qrs_i_raw)
    beat(i,:)=ECG(qrs_i_raw(i+1)-before:qrs_i_raw(i+1)+after);
    
end
%qrs_i_raw has the index location of the peak

it gives me this error

Index exceeds the number of array elements (650000).

I understand why I get this error, but how can I modify the loop so for the first and last beat it will take just the existing points before the first beat and the remaining points after the last beat?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Star Strider 2021-2-9

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/738892-signal-segmentation-error-how-to-select-a-specific-number-of-points-before-and-after-the-ecg-signal#answer_618872

See How do I find and plot the average of action potentials from a trace? for an appropriate approach. Substitute your EKG vector and time vector for the signals in that code, and adjust the ‘ofst’ value to be the one you want (typically about 450 ms for an EKG, so that value would depend on the sampling frequency of your EKG and the heart rate, so multiply the sampling frequency in Hz by 0.45 to get the number of samples to use for ‘ofst’). It may be necessary to adjust that for the heart rate, so typically that interval should be slightly less than ½ of the previous — or mean — R-R interval.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 2

Mathieu NOE 2021-2-9

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/738892-signal-segmentation-error-how-to-select-a-specific-number-of-points-before-and-after-the-ecg-signal#answer_618787

在 MATLAB Online 中打开

hello Ioana

I modified a bit your code , hope it works

before=100;
after=140;
nn = length(qrs_i_raw);
% first and last buffer - specific cases
beat(1,:)=ECG(1:qrs_i_raw(1)+after);
beat(nn,:)=ECG(qrs_i_raw(nn)-before:end);
% general cases
for i=2:nn-1
    beat(i,:)=ECG(qrs_i_raw(i)-before:qrs_i_raw(i)+after); % pls note index i is same on both sides of the equation
    
end
%qrs_i_raw has the index location of the peak