get_ecg.m

function [ecg] = get_ecg(time,V_thresh,V_stim)
    
%%%%%Generates the ECG waveform

    x=time;
    V=V_stim;% example input  voltage
    threshold=V_thresh; % example threshold
    if(V<threshold)
          li=30/72;  
    
          a_pwav=0.25;
          d_pwav=0.09;
          t_pwav=0.16;  

          a_qwav=0.025;
          d_qwav=0.066;
          t_qwav=0.166;

          a_qrswav=1.6;
          d_qrswav=0.11;

          a_swav=0.25;
          d_swav=0.066;
          t_swav=0.09;

          a_twav=0.35;
          d_twav=0.142;
          t_twav=0.2;

          a_uwav=0.035;
          d_uwav=0.0476;
          t_uwav=0.433;
          pwav=p_wav(x,a_pwav,d_pwav,t_pwav,li);

     %qwav output
     qwav=q_wav(x,a_qwav,d_qwav,t_qwav,li);

     %qrswav output
     qrswav=qrs_wav(x,a_qrswav,d_qrswav,li);
     %swav output
     swav=s_wav(x,a_swav,d_swav,t_swav,li);

     %twav output
     twav=t_wav(x,a_twav,d_twav,t_twav,li);

     %uwav output
     uwav=u_wav(x,a_uwav,d_uwav,t_uwav,li);
     %ecg output
     ecg=pwav+qrswav+twav+swav+qwav+uwav;
    else

     ecg=0.89892*ones(1,length(x));



    end

end