matlab根轨迹探究

G(S) = K/(S+10)

figure(1);
z = [];
p = [-10];
k = 1;
[num,den] = zp2tf(z,p,k);
subplot(1,3,1);
rlocus(num,den);
title('K = 1,G(S) = 1/S+10');
k1 = 10;
subplot(1,3,2);
[num1,den1] = zp2tf(z,p,k1);
rlocus(num1,den1);
title('K = 10,G(S) = 10/S+10');
k2=100;
[num2,den2] = zp2tf(z,p,k2);
subplot(1,3,3);
rlocus(num2,den2);
title('K = 100,G(S) = 100/S+10');//by 正茂.dog.jpg

G(S) = 10/(S+2)(S+6)

figure(2);
z = [];
p = [-2,-6];
k = 1;
[num,den] = zp2tf(z,p,k);
subplot(1,3,1);
rlocus(num,den);
title('K = 1,G(S) = 1/(S+2)(S+6)');
k1 = 10;
subplot(1,3,2);
[num1,den1] = zp2tf(z,p,k1);
rlocus(num1,den1);
title('K = 10,G(S) = 10/(S+2)(S+6)');
k2=100;
[num2,den2] = zp2tf(z,p,k2);
subplot(1,3,3);
rlocus(num2,den2);
title('K = 100,G(S) = 100/(S+2)(S+6)');//by 正茂.dog.jpg

G(S) =K(S+2)/(S+1)(S+3)^2

z = [-2];
p = [-1,-3,-3];
k = 1;
[num,den] = zp2tf(z,p,k);
subplot(1,3,1);
rlocus(num,den);
title('K = 1,G(S) = K*(S+2)/(S+1)(S+3)(S+3)');

k1 = 10;
subplot(1,3,2);
[num1,den1] = zp2tf(z,p,k1);
rlocus(num1,den1);
title('K = 10');

k2=100;
[num2,den2] = zp2tf(z,p,k2);
subplot(1,3,3);
rlocus(num2,den2);
title('K = 100');//by 正茂.dog.jpg

G(S )= K(S+3)/S(S+2)(S^2+10S+50)

figure(4);
z = [-3];
p = [0,-2,roots([1,10,50])'];
k = 1;
[num,den] = zp2tf(z,p,k);
subplot(1,3,1);
rlocus(num,den);
title('K = 1,G(S) = K(S+3)/S(S+2)(S^2+10S+50)');
k1 = 10;
[num1,den1] = zp2tf(z,p,k1);
subplot(1,3,2);
rlocus(num1,den1);
title('K = 10');
k2=100;
[num2,den2] = zp2tf(z,p,k2);
subplot(1,3,3);
rlocus(num2,den2);
title('K = 100');//BY ZEMAN

G(S)=K(S+1.5)(S^ 2 +4S+5)/S(S+2.5)(S^2+S+2.5)

figure(5);
z = [-1.5,roots([1 4 5])']';
p = [0,-2.5,roots([1 1 2.5])'];
k = 1;
[num,den] = zp2tf(z,p,k);
subplot(1,3,1);
rlocus(num,den);
title('K = 1,G(S) = K(S+1.5)(S^2+4S+5)/S(S+2.5)(S^2+S+2.5)');
k1 = 10;
subplot(1,3,2);
[num1,den1] = zp2tf(z,p,k1);
rlocus(num1,den1);
title('K = 10');
k2=100;
[num2,den2] = zp2tf(z,p,k2);
subplot(1,3,3);
rlocus(num2,den2);
title('K = 100');//BY ZEMAN