exp4_11.m

matlab与控制系统仿真的PPt多媒体教程和源码

%exp4_11.m
%已知系统的传递函数为：G(s)=K/(s^3+52s^2+100^s)
%求当K分别取1300和5200时，系统的极坐标频率特性图
clear
close all
clc
k1=1300;
k2=5200;
w=1:0.1:80;%指定角频率范围
num1=k1;
num2=k2;
den=[1 52 100 0];
figure(1)
nyquist(num1,den,w);%自动绘制
figure(2)
nyquist(num2,den,w);%对应K2
figure(3)
[rm,im]=nyquist(num2,den);%返回实部和虚部，没有指定角频率范围
plot(rm,im)
xlabel('real')
ylabel('image')
text(-18,1200,'w from 负无穷 to 零')