clear,clc
m=0;
for i=1:200
w=0.5*i;
for j=1:80
k=0.01*j*w;
num=[k*w*w];
den=[1,0.4*w,w*w,k*w*w];
t=[0:0.01:100];
[c,x,t]=step(num,den,t);
ma=max(c);
tp=find(c==ma);
if (tp<100)
if (1<ma<1.05)
m=m+1;
ck(m)=k;
cw(m)=w;
end
end
end
end
plot(cw,ck,'.');
axis([0,100,0,40])
很多步骤我不太清楚具体是干嘛的，求帮助，不胜感激，
