自动控制原理MATLAB仿真实验三控制系统的频域.docxVIP

实验三 MATLAB及仿真实验(控制系统的频域分析) 一 实验目的 利用计算机作出开环系统的波特图 观察记录控制系统的开环频率特性 控制系统的开环频率特性分析 二实验内容 1.用Matlab作Bode图.要求：画出对应 Bode图,并加标题. (1)G(s)二 2 25 — s2 +4s +25 num=[25]; den=[1 4 25]; G=tf( num,de n); margin(G);% 增益和相位裕度 title( Bode 图); grid on; BodeE： -10 -20 -30 〔mpj 否口 nlhlJdHS巨 -90 -135 ^180 Frequency (rad/s) (畀邑器?匸CL (2)G(s)9(s2 0.2s 1) (2) G(s) s(s2 1.2s 9) num=[9 1.8 9];den=co nv([1 0],[1 1.2 9]); G=tf( num,de n); margi n( G); BodeEititle( Bode 图); grid on; BodeEi 20 10 0 -10 20 -30 Frequency (rad/s)90 Frequency (rad/s) 45 O9IS呼Ha-45 O9IS呼Ha -45 2.用Matlab作Nyquist 图.要求画对应 Nyquist图，并加网格和标题 1 G (s)=— s +0.8s+1 n um=[1]; den=[1 0.8 1]; G=tf( num,de n); ny quist(G); title( 奈奎斯特图’)； grid on; M爰 XJBEEndui- .5o 2 de 4dB 10 20 dB I 0 0.5 Real Axis 3?典型二阶系统Gy）： 2〔 2，试绘制 取不同值时的Bode图。取 n =6』：工[0.1: 0.1:1.0]。 for e=[0.1:0.1:1] num=[36]; den=[1 12*e 36]; G=tf( num,de n); margi n( G); hold on; end 20Bode DiagramGm = Inf dB (at Inf rtud/s), Pm = -1B0 deg (d 0 rad/s)o^18010 20 Bode Diagram Gm = Inf dB (at Inf rtud/s), Pm = -1B0 deg (d 0 rad/s) o ^180 10_1 10 Frequency (隔松) 10J 35 某开环传函为： G(s)二 50 ，试绘制系统的 Nyquist曲线，并判断闭环系统稳 (s+5)(s-2) 定性，最后求出闭环系统的单位脉冲响应。 n um=[50]; den=co nv([1 5],[1 -2]); G=tf( num,de n); figure(1) ny quist(G); title( 奈奎斯特图’)； figure(2) impulse(G); title( 单位脉冲响应); 1_5?菱 XJBEmaiui-奈豪商祷閉-3.5 -3 -2.$ -2 -1.5 -1 -g 0Real AxisQpn=dlu鱼更加冲有应10 1S 萸S-4 3 2 1_5 ?菱 XJBEmaiui- 奈豪商祷閉 -3.5 -3 -2.$ -2 -1.5 -1 -g 0 Real Axis Qpn=dlu 鱼更加冲有应 10 1S 萸 S- 4 3 2 Time 泪econd却 由图可知曲线有（-1 , j0 ）左侧的极点，为一次副穿越，R=-2 ;且P=1不等 于R;所以系统不稳定 1 1 5. T2s2 2 Ts 1 T =0.1 1, 0.5, 0.1 作波特曲线图 for e=[0.1 0.5 1 2] n um=[1]; den=[0.01 0.2*e 1]; G=tf( num,de n); margi n( G); hold on; end 6.Bodfr DiagramGm = Inf dB (at Inf racWs). Pm = -1B0 deg (? 0 rndA)-90Frequency (rad/s)10J：31.6s 0.01s 1 0.1s 1, 6. Bodfr Diagram Gm = Inf dB (at Inf racWs). Pm = -1B0 deg (? 0 rndA) -90 Frequency (rad/s) 10J： 31.6 s 0.01s 1 0.1s 1 ,要求: (a)作波特图 num=[31.6]; den仁 co nv([1 0],[0.01 1]); den=co nv(de n1,[0.1 1]); G=tf( nu m,de n); m