somebody plz help me with this code. I got an error after run this code is subscript indices must either be real positive integers or logical.

for i=1:length(left);
   
    [R_A(i) R_t(i)]=max(sigL(left(i):raight(i)));
    R_t(i)=R_t(i)-1+left(i) %add offset
   
    [Q_A(i) Q_t(i)]=min(sigL(left(i):R_t(i)));
    Q_t(i)=Q_t(i)-1+left(i)
  
    [S_A(i) S_t(i)]=min(sigL(left(i):raight(i)));
    S_t(i)=S_t(i)-1+left(i)
    
    [P_A(i) P_t(i)]=max(sigL(left(i):Q_t(i)));
    P_t(i)=P_t(i)-1+left(i)
    
    [T_A(i) T_t(i)]=max(sigL(S_t(i):raight(i)));
    T_t(i)=T_t(i)-1+left(i)+47;
end

[per isakson]

To make it easier for us to help, provide the data we need to run your code.

[Guillaume]

MANISHA KUMARI's comment originally posted as an Answer moved here:

sig=load('ecg_60hz_200.dat');
N=length(sig);
fs=200;
t=[0:N-1]/fs;
figure(1);subplot(4,2,1);plot(sig)
title('Original Signal')
%%
     %           Low Pass Filter
b=1/32*[1 0 0 0 0 0 -2 0 0 0 0 0 1];
a=[1 -2 1];
sigL=filter(b,a,sig);
subplot(4,2,3);plot(sigL)
title('Low Pass Filter')
subplot(4,2,4);zplane(b,a)
%%
     %           High Pass Filter
b=[-1/32 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 -1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1/32];
a=[1 -1];
sigH=filter(b,a,sigL);
subplot(4,2,5);plot(sigH)
title('High Pass Filter')
subplot(4,2,6);zplane(b,a)
%%
     %          Derivative Base Filter
b=[1/4 1/8 0 -1/8 -1/4];
a=[1];
sigD=filter(b,a,sigH);
subplot(4,2,7);plot(sigD)
title('Derivative Base Filter')
subplot(4,2,8);zplane(b,a)
%%
     %      be tavane 2 miresanim
sigD2=sigD.^2;
%%
     %      normalization
signorm=sigD2/max(abs(sigD2));
%%
  
h=ones(1,31)/31;
sigAV=conv(signorm,h);
sigAV=sigAV(15+[1:N]);
sigAV=sigAV/max(abs(sigAV));
figure(2);plot(sigAV)
title('Moving Average filter')
%%
treshold=mean(sigAV);
P_G= (sigAV>0.01);
figure(3);plot(P_G)
title('treshold Signal')
figure;plot(sigL)
%%
difsig=diff(P_G);
left=find(difsig==1);
raight=find(difsig==-1);
%%
     %      run cancel delay
     %      6 sample delay because of LowPass filtering
     %      16 sample delay because of HighPass filtering
left=left-(6+16);
raight=raight-(6+16);
%%
    % P-QRS-t
for i=1:length(left);
   
    [R_A(i) R_t(i)]=max(sigL(left(i):raight(i)));
    R_t(i)=R_t(i)-1+left(i) %add offset
   
    [Q_A(i) Q_t(i)]=min(sigL(left(i):R_t(i)));
    Q_t(i)=Q_t(i)-1+left(i)
  
    [S_A(i) S_t(i)]=min(sigL(left(i):raight(i)));
    S_t(i)=S_t(i)-1+left(i)
    
    [P_A(i) P_t(i)]=max(sigL(left(i):Q_t(i)));
    P_t(i)=P_t(i)-1+left(i)
    
    [T_A(i) T_t(i)]=max(sigL(S_t(i):raight(i)));
    T_t(i)=T_t(i)-1+left(i)+47
end
%%
figure;plot(t,sigL,t(Q_t),Q_A,'*g',t(S_t),S_A,'^k',t(R_t),R_A,'ob',t(P_t),P_A,'+b',t(T_t),T_A,'+r');
for i=1:((length(P_t))-1)
    
    HRV=P_t(i+1)-P_t(i)
end

[Guillaume]

per isakson's reply moved here as well:

You should use Comments to add to your question. From the outside it now looks as if the questions has been answered.

Cannot run your code without some values for the variable sig.

[MANISHA KUMARI]

new.txt is the ecg raw signal value.

In code there will be an editing of sig=load('new.txt'); instead of sig=load('ecg_60hz_200.dat');

Answer 1

[Guillaume]

%your code:
left=find(difsig==1);
raight=find(difsig==-1);

So left and the mispelled raight are a vector of indices of difsig. They're positive integers that can be any value from 1 to numel(difsig),

%[...]
left=left-(6+16);
raight=raight-(6+16);

From these positive integers you then subtract 10, so they're now integers from -9 to numel(difsig)-9,

%[...]
    [R_A(i) R_t(i)]=max(sigL(left(i):raight(i)));

which you're now using as indices. Clearly, there's a problem here since if left or the mispelled raight are less than 1, they can't be used to index.

Unfortunately, the comments don't really explain what you're trying to do, so it's difficult to tell you how to fix it other thanmaking sure that the indices you generate can't be less than 1.

Answer 2

[Image Analyst]

See the FAQ