Nyquist图 源代码

G=tf([1 5],[1 2 3 0]); rlocus(G) grid figure bode(G) grid figure nyquist(G) grid figure Gc=feedback(G,1) step(Gc) grid

% ex4_3 % 放大系数对频率响应的影响 figure('pos',[30 100 260 260],'color','w'); axes('pos',[0.15 0.2 0.7 0.7]); G1=tf([1],[5 1]); nyquist(…

% ex4_21 % 离散控制系统的频率响应图。 figure('pos',[50 50 750 250],'color','w'); axes('pos',[0.05 0.2 0.25 0.7]); G=tf(5,[1 6 11 6 0]); bo…

% ex4_5 % 几种类型系统的奈魁斯特图 figure('pos',[30 100 260 260],'color','w'); axes('pos',[0.15 0.2 0.7 0.7]); wn=1; w=[0,logspace(-2,5,…

% ex4_23 % 含时滞的一阶惯性环节的衰减频率响应 figure('pos',[30 100 600 260],'color','w'); axes('pos',[0.055 0.2 0.25 0.7]); G1=tf([5],[5 1]); …

% ex4_13 % 奈魁斯特稳定性判别示例 figure('pos',[30,30,300,250],'color','w'); axes('pos',[0.2,0.2,0.7,0.7]); T1=3; T2=2; T3=1; for k=[…

%fig6_40 figure('pos',[100,100,200,200],'color','w'); G=tf(1,1); set(G,'inputdelay',2); Gc=tf([1 1],1) nyquist(G*Gc,[0:.1:2]'…

% ex4_2 % 时滞对奈魁斯特图的影响。 figure('pos',[30 100 260 260],'color','w'); axes('pos',[0.15 0.2 0.7 0.7]); G1=tf([1],[5 1]); nyquist(…

for K=0.2:0.5:3 num=K; den=[5 3 1]; G=tf(num,den); nyquist(G); hold on; end grid

%exp4_11.m %已知系统的传递函数为：G(s)=K/(s^3+52s^2+100^s) %求当K分别取1300和5200时，系统的极坐标频率特性图 clear close all clc k1=1300; k2=5200; w=8:1:…