tìm thời gian quá độ của hệ thống điều bằng hàm trong matlab

[canhbg]

Chào các bác!
em có đoạn tìm thời gian quá độ của hệ thống điều khiển như sau, sử dụng hàm trong matlab
% tim thoi gian qua do
% de chay can dat ham dang text(.txt)
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 |%
|%
%-------------------------------------------------------|%

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)<ymin)
tset=t(k);
break;
end
if (k==1)
error('Loi nhap du lieu!')
tset=null;
end
end
clear y t ymin ymax

nhưng em khi em tạo M file rồi chạy vẫn chưa được
các bác, anh Bùi Hiếu giúp em với
1. LÀM THẾ NÀO để xác định được
1.1 thời gian quá độ (tqd t settling time)
1.2 thời gian hệ thống ổn định ( t rise time) (viết hàm như thế nào các bác nhỉ)
1.3 độ quá điều chỉnh (Mp)
1.4 sai số xác lập (Ess)
Rất mong nhận được hồi âm
Cảm ơn các bác nhiều