高手来帮帮忙啊 俺谢谢了啊 只要得个图形

dxylww

用MATLAB仿真一个系统,调用一个函数ste(),调用格式为[sigma,tp,ts]=ste(y,t)运行正确可以得到系统的阶越响应曲线, 但我的MATLAB7怎么都得不出来,谢谢哪位好心的朋友帮我一下 ,告诉我这到底怎么调用哪个函数,能把你得到的图形回贴发过来最好 ,真心感谢 !!!!!!~~~~~~~~~~~~~加我QQ350480320指导我也行 ,总之就是谢你了!!!!!!~~

程序如下:
function [sigma,tp,ts]=ste(y,t)
[mp,tf]=max(y);
cs=length(t);
yss=y(cs);
sigma=(mp-yss)/yss
tp=t(tf)
i=cs+1;n=0,
while n==0,
          i=i-1;
             if i==1,
                    n=1;
             elseif y(i)>1.05*yss,
                    n=1;
             end;
end;
t1=t(i);cs=length(t);j=cs+1;
n=0;
while n==0,
          j=j-1,
             if j==1,
                    n=1;
             elseif y(i)>0.95*yss,
                    n=1;
             end;
end;
t2=t(j);
if t2<tp
       if t1>t2
              ts=t1
       end
elseif t2>tp,
      if t2<t1,
             ts=t2
      else
             ts=t1
      end
end



global y t
s2=tf([1 3.841],[1 73.12]);
s1=tf(20700,[1 20 75 0]);
sope=s1*s2;sys=feedback(sope,1)
step(sys)
[y,t]=step(sys);
[sigma,tp,ts]=ste(y,t)