本帖最后由 ☞屌丝☞Yxm☆ 于 2013-3-26 04:37 编辑

下面一段程序是我网上找到的，但是不全，求各位大神帮忙补充下好么，我在做练习额！
急用！！！！
》》》红 {MOD}@处表示不全，需要补充的。《《《《
num=[k];
den=[T1*T2,T1+T2,1];
g0=tf(num,den);
t=(0:0.1:70);
y1=step(feedback(kp*g0,1),t);
n=length(t);
[y1max,ind]=max(y1);
y1ss=y1(n);
y1ssxw=k/T1*T2;
ess=y1ssxw-y1ss;
tp=t(ind);
mp=(y1max-y1ss)/y1ss*100;
for k=1:n
   if y1(k)<=0.9*y1ss&y1(k+1)>=0.9*y1ss
        tr=t(k);
break
end
   end
while y1(n)>0.98*y1ss&&y1(n)<1.02*y1ss;
n=n-1;
end
ts=t(n);
num1=[k*kp*td k*kp];
den1=[T1*T2 T1+T2 1];
y2=step(feedback(tf(num,den)));
[y2max,ind]=max(y2);
y2ss=y2(n);
mp2=100*(y2max-y2ss)/y2ss;
tp2=t(ind);
for k=1:n
if y2(k)<=0.9*y2ss&y2@@@@@@@；
tr2=t(k);
break
end
end
ess2=1-y2ss;
while y2(n)>0.98*y2ss&@@@@@@@;
n=n-1;
end
ts2=t(n);
num3=[kp*ti 1];
den3=[ti];
g1=tf(num3,den3);
y3=step(feedback(g0*g1@@@@@@@@@@);
[y3max,ind]=max(y3);
y3ss=y3(n)
@@@