用MATLAB绘制Nyquist图.ppt

用MATLAB绘制Nyquist图
MATLAB命令
nyquist命令可以求得连续系统的奈奎斯特曲线。
命令格式:
[ re,im,w ] = nyquist(num,den)
[ re,im,w ] = nyquist(num,den,w)
当带有输出变量时,可得到相应的一组数据,不带输出变量时,则绘出奈奎斯特曲线。也可用指定频率向量w指定所要绘制的曲线范围。
11/13/2018

2
11/13/2018

3
(a) z=1,p=2
n=[0 0 1 1]
m=[1 2 0 0]
sys=tf(n,m)
nyquist(sys)
11/13/2018

4
(b) z=2,p=1
11/13/2018

5
幅值裕度和相位裕度
例三十二: 系统开环传递函数为

绘制k = 10,26,50时的带单位圆的奈奎斯特图,
单位圆的绘制是通过绘绘tjw的实部与虚部的轨迹而获得。
幅值裕度和相位裕度
n = 10; %取k = 10时的值
d = conv([1 2],[1 2 5]);
w = [ 0::10 ]; %确定频率范围
e = exp(j*w); %给出指数函数ejw
r = real(e); %求指数函数的实部,结果不显示
i = imag(e); %求函数ejw的虚部,结果不显示
[ a,b ] = nyquist(n,d,w); %求指定频率范围内的奈氏值,不显示结果
n1 = 26; %取k = 26
d1 = d; %保留原分母矢量
[ a1,b1 ] = nyquist(n1,d1,w); %求k = ,结果不显示
n2 = 50; %取k = 50
d2 = d; %分母保留
[ a2,b2 ] = nyquist(n2,d2,w);
plot(r,i,a,b,a