自动控制大作业--西安电子科技大学--matlab实现，运行完美实现！(3)

二阶微分环节150Magnitude (dB)Phase (deg)100500-501801359045010-1100101102103Frequency (rad/sec) 振荡环节50Magnitude (dB)Phase (deg)0-50-100-1500-45-90-135-18010-1100101102103Frequency (rad/sec) 原总波特图0Magnitude (dB)Phase (deg)-50-100-1500-45-90-135-18010-1100101102103Frequency (rad/sec)

源程序：

num=[0,20,20,10];

den=conv([1,1,0],[1,10]); nyquist(num,den)

输出：

Nyquist Diagram2015105Imaginary Axis0-5-10-15-20-1-0.500.5Real Axis11.52

源程序：

num=[2000,2000];

den=conv([1 0.5 0],[1 14 400]); nichols(num,den) v = [-270 -90 -40 40]; axis(v) ngrid

输出：

Nichols Chart40 0 dB30 0.25 dB 0.5 dB20 1 dB 3 dB 6 dB0Open-Loop Gain (dB)10-10-20-30-40-270-225-180Open-Loop Phase (deg)-135-90

源程序：

num = [0 2000 2000]; den = conv([1 0.5 0],[1 14 400]); h=tf(num,den); bode(h);

num = [0 2000 2000]; den = conv([1 0.5 0],[1 14 400]); h=tf(num,den);

[gm,pm,wg,wc]=margin(h)

输出：

gm =

2.7493

pm =

73.3527

wg =

19.8244 wc =

5.3477 Bode Diagram10050Magnitude (dB)Phase (deg)0-50-100-150-90-135-180-225-27010-210-1100101102103Frequency (rad/sec)

源程序：

num=[0 0 0 1];

den=conv([.5 1 0],[1 1]); sysp=tf(num,den) sys=feedback(sysp,1) w=logspace(-1,1); bode(sys,w)