用matlab指数传递函数nyquist,用MATLAB制Nyquist图.ppt

用MATLAB制Nyquist图

幅值裕度和相位裕度 例三十二： 系统开环传递函数为 绘制k = 10，26，50时的带单位圆的奈奎斯特图， 单位圆的绘制是通过绘绘tjw的实部与虚部的轨迹而获得。 幅值裕度和相位裕度 n = 10; %取k = 10时的值 d = conv([1 2],[1 2 5]); w = [ 0:0.01: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 = 0.5时的奈氏值，结果不显示 n2 = 50; %取k = 50 d2 = d; %分母保留 [ a2,b2 ] = nyquist(n2,d2,w); plot(r,i,a,b,a1,b1,a2,b2),grid %绘出：r，i；a，b；a1，b1；a2，b2；的对应图形并加上栅格。 幅值裕度和相位裕度 执行以上程序后可在上图上得k = 10，26，50并加有单位圆的奈奎斯特图。一般来说由于此列的关系显示的图形不是一个正规的圆。 从图形上我们可以看出开环增益对闭环系统稳定性的影响: 当K值变化时，幅频特性成比例变化，而相频特性不受影响。 取k=26时，曲线恰好通过(-1,0j)点，这是临界稳定状态； 当k<26时，幅相曲线将从(-1,0j)点的右方穿过负实轴，不再包围(-1,0j)点，这时闭环系统是稳定的； 当k>26时，开环频率特性随着从变化到时，顺时针方向围绕(-1,0j)点一圈，即N=-1，可求得闭环系统在右半平面的极点数为： ，所以闭环系统不稳定。 * * 用MATLAB绘制Nyquist图 MATLAB命令 nyquist命令可以求得连续系统的奈奎斯特曲线。 命令格式： [ re，im，w ] = nyquist(num，den) [ re，im，w ] = nyquist(num，den，w) 当带有输出变量时，可得到相应的一组数据，不带输出变量时，则绘出奈奎斯特曲线。也可用指定频率向量w指定所要绘制的曲线范围。 * Y.Xing * * Y.Xing * (a) z=1,p=2 n=[0 0 1 1] m=[1 2 0 0] sys=tf(n,m) nyquist(sys) * Y.Xing * (b) z=2,p=1 * Y.Xing *