get_yss.m

来自「其中提到遺傳學的程式碼與應用提供給次淚相向的研究者參考下載」· M 代码 · 共 15 行

function yss=Get_yss(t,y,FT);
%% yss=Get_yss(t,y,FT);
%% Get steady-state value of a waveform (y)
%% t-time, y-waveform, FT-Final simulation Time
%% PenChen Chou. 2004/1/15. 2004-4-16. 2005-4-16.

if nargin~=3, 
    error('Error! No enough input arguments (need 3)!'); 
end
N=length(y);II=find(t>0.9*FT);y=y(II);
if abs(std(y))<0.02 & abs(y(end-3)-y(end))<0.1
   yss=y(end);
else
   yss=inf;
end