Location: BG_LCC_diadic @ 0e12b0f14ccd / matlab_parameter_fitting / calc_APD_array.m

Author:
Shelley Fong <sfon036@UoA.auckland.ac.nz>
Date:
2024-11-13 11:18:14+13:00
Desc:
Updating cellml and sedml
Permanent Source URI:
https://staging.physiomeproject.org/workspace/891/rawfile/0e12b0f14ccd4923dc8b6f7fd13ba9f3b0ca1d20/matlab_parameter_fitting/calc_APD_array.m
Compare... Download Source 

function APD = calc_APD_array(t,V)

num_beats = floor(t(end));
APD = zeros(num_beats,1);

for i_beat = 1:num_beats
    Idx_beat = ((1+1000*(i_beat-1)):(1+1000*i_beat));
    t_beat = 1000*(t(Idx_beat)-i_beat+1)-300;
    V_beat = V(Idx_beat);
    V_peak = max(V_beat);
    V_end = V_beat(end);
    V_90 = 0.9*V_end + 0.1*V_peak;
    Idx_repolarised = find((V_beat <= V_90) & (t_beat > 20),1);
    Idx_cross = [Idx_repolarised-1 Idx_repolarised];
    APD(i_beat) = interp1(V_beat(Idx_cross),t_beat(Idx_cross),V_90);
end



end