%exp4_11.m
%已知系統(tǒng)的傳遞函數(shù)為：G(s)=K/(s^3+52s^2+100^s)
%求當(dāng)K分別取1300和5200時(shí)，系統(tǒng)的極坐標(biāo)頻率特性圖
clear
close all
clc
k1=1300;
k2=5200;
w=8:1:80;
num1=k1;
num2=k2;
den=[1 52 100 0];
figure(1)
subplot(211)
nyquist(num1,den,w);
subplot(212)
pzmap(num1,den);
figure(2)
subplot(211)
nyquist(num2,den,w);
subplot(212)
[rm,im]=nyquist(num2,den);
plot(rm,im)
xlabel('real')
ylabel('image')
title('w from 負(fù)無窮 to 零')
figure(3)
[numc,denc]=cloop(num2,den);
subplot(211)
step(numc,denc)
subplot(212)
[numc1,denc1]=cloop(num1,den);
step(numc1,denc1)