clear;
wn=0.25648*2*pi;
clpoles=[];param=[];
y0=[];
for x=0:0.02:50
num=[(2*x+1)*wn (2*x+1)*wn*wn wn^3];
den=[1 (2*x+1)*wn (2*x+1)*wn*wn wn^3];
clp=roots(den);
clpoles=[clpoles;clp'];param=[param;x];
end
disp([param,clpoles])
plot(clpoles,'-');
v=[-2 1 -2 2];
axis(v);axis('square')
sgrid([0.707] ,[])
%hold on
%ind=[1 2 3];%就是你需要的点的序号
%plot(clpoles(2,3),'*')
%text(clpoles(2,2),'o','color','g')
%if(-3<= real(clpoles)<=0)
% [y0]=clpoles;
%end
%hold on
%plot([y0],'r.','markersize',20)
%k=5;den=conv(conv([1 1],[1 3]),[1 12]);
%clpoles=[];param=[];
%for alpha=2:10
% num=[0,0,k,k*alpha];
% clpoly=num+den;
% clp=roots(clpoly);
% clpoles=[clpoles;clp'];param=[param;alpha];
%end
%disp([param,clpoles])
%plot(clpoles,'-')
%v=[-2 1 -2 2];
%axis(v);axis('square')
%sgrid([0.707] ,[])