【正文】
系統(tǒng)的 bode圖: hold on margin(gn) 102101100101102 1 8 0 1 5 0 1 2 0 9 0Phase (deg)B o d e D i a g r a mG m = In f d B ( a t In f r a d / s e c ) , P m = 5 3 . 8 d e g ( a t 3 . 1 9 r a d / s e c )F r e q u e n c y ( r a d / s e c ) 6 0 4 0 2 00204060Magnitude (dB) 圖 22 校正后系統(tǒng)的 bode 圖 畫出其單位階躍響應圖: figure,step(feedback(g0,1)) hold on,step(feedback(gn,1),39。 45176。g0=tf(n,d) Transfer function: 4 s^2 + s margin(g0) 4 0 3 0 2 0 1 00102030405060Magnitude (dB)102101100101 1 8 0 1 3 5 9 0Phase (deg)B o d e D i a g r a mG m = In f d B ( a t In f r a d / s e c ) , P m = 1 4 . 2 d e g ( a t 1 . 9 7 r a d / s e c )F r e q u e n c y ( r a d / s e c ) 圖 21 原系統(tǒng)的 bode 圖 Pm=176。 wc=3)求所需的超前相位角 pha: + ans = 29 pha=ans pha = 29 4)計算α: alpha=(1+sin(pha*pi/180))/(1sin(pha*pi/180)) alpha = 5)求 wm: 10*log10(alpha) ans = 101100101102 1 8 0 1 3 5 9 0Phase (deg)B o d e D i a g r a mG m = In f d B ( a t In f r a d / s ) , P m = 2 5 . 2 d e g ( a t 4 . 2 5 r a d / s )F r e q u e n c y ( r a d / s ) 6 0 4 0 2 00204060S y s t e m : g 0F r e q u e n c y ( r a d / s ) : 5 . 6 6M a g n i t u d e ( d B ) : 4 . 6Magnitude (dB) 圖 12 尋求 wc的 bode圖 此時, wm=6rad/s,不滿足要求; 6)重新設計,利用ω c求α: 取ω m即新的ω c為 6rad/s,讀取此時的 L(ω ),由圖得 L(ω )=,則運行下列程 序,求α: