MATLAB绘制控制系统的极坐标图
极坐标绘图函数 nyquist图
函数用法
nyquist(sys)
nyquist(sys,w)
[re,im,w]= nyquist(sys)
第一种使用方法：给定开环系统的模型对象sys，自动选择w的绘制频率
第二种使用方法：给定开环系统的模型对象sys，手动选择w的绘制频率
第三种使用方法：返回极坐标图的参数向量，不做图
re为复变函数G(jw)的实部向量，im为虚部向量
绘制如下函数的极坐标图
G o p e n ( s ) = 10 s 2 + 2 s + 5 G_{open} (s)= \frac{10}{s^{2}+2s+5 } Gopen(s)=s2+2s+510
MATLAB代码如下：
num=[0 0 10]%分子，作多项式模型
den=[1 2 5]%分母
sys=tf(num,den)%构建系统对象
nyquist(sys)
输出：
num = 1×3
0 0 10
den = 1×3
1 2 5
sys =
10 s 2 + 2 s + 5 \frac{10}{s^{2}+2s+5 } s2+2s+510
Continuous-time transfer function.
图像：
绘制效果不理想的改进方法
1.使用axis()改变坐标显示范围
如：axis([-1,1.5,-2,2])
2.给定角频率指定范围
w=0:0.1:100;
nyquist(sys,w)