get_risetime.m

其中提到遺傳學的程式碼與應用提供給次淚相向的研究者參考下載

function riseTime=Get_riseTime(inp_cmd,yss,t,y,FT)
%% riseTime=Get_riseTime(inp_cmd,yss,t,y,FT);
%% Find Rise time of a waveform y
%% yss-Steady-state value of y, FT-Final Time
%% PenChen Chou. 2004/1/15. 2005-4-16.

if nargin~=5, 
    error('Error! No enough input arguments (need 5)!'); 
end

if abs(yss)/abs(inp_cmd)>3 | abs(yss)/abs(inp_cmd)<0.333
  riseTime=FT;
else
  if yss==eps
    riseTime=FT;
  else
    II=find(y>=0.1*yss);
    if isempty(II)
      riseTime=FT;
    else
      Ix=II(1); clear II
      II=find(y>=0.9*yss);
      if isempty(II)
         riseTime=FT;
      else
         riseTime=t(II(1))-t(Ix);
      end; %keyboard;pause
    end
  end
end