% hàm tìm thời gian quá độ
function[tset]=settling_time(g,stardard)
%This function compute settling time, input data is tranfer function, error
%stardard (0.02 or 0.05 or other) if stardard=null, I get: stardard=2%
%——————————————————-|%
%|(C)2005 Bui Trung Hieu |%
%|Website: http://www.khvt.com |%
%|Email: hieu@khvt.com |%
%——————————————————-|%

gtxl=freqresp(g,0);
if(nargin==1)
stardard=0.02;
end

ymin=gtxl*(1-stardard);
ymax=gtxl*(1+stardard);
[y,t]=step(g);
for k=length(y):-1:1
if ((y(k)>ymax)|y(k)f(x)
%——————————————————-|%
%|(C)2005 Bui Trung Hieu |%
%|Website: http://www.khvt.com |%
%|Email: hieu@khvt.com |%
%——————————————————-|%

exponent=g.iodelay;
[z,p,fac]=zpkdata(g,’v');
f=sym(’1′);
syms x
k=1;
while(k<=length(z))
f=f*(x-z(k));
k=k+1;
end
k=1;
while(k<=length(p))
f=f/(x-p(k));
k=k+1;
end
f=f*fac*exp(-exponent*x);
%%End function