【正文】
12 14 16 1 08642024681013579111315N y q u is t D ia g r a mR e a l A x isImaginary Axis 可以看出 隨著 k 的增大,該系統(tǒng)的 Nyquist 曲線向外擴展,但不影響穩(wěn)定性,如圖所示。 s=num2str(k)。 bode(num,den)。 for ks=:: den=[a 2*T*ks 1]。 結(jié)果及打印曲線: 6 5 4 3 2 1 0 1 . 51 0 . 500 . 511 . 50 . 10 . 30 . 50 . 70 . 90 . 10 . 30 . 50 . 70 . 9123456R o o t L o c u sR e a l A x is ( s e c o n d s1)Imaginary Axis (seconds1) 8 0 0 . 1 0 . 2 0 . 3 0 . 4 0 . 5 0 . 6 0 . 7 0 . 8 0 . 900 . 10 . 20 . 30 . 40 . 50 . 60 . 70 . 80 . 91S t e p R e s p o n s eT im e ( s e c o n d s )Amplitude Select a point in the graphics window selected_point = k = p = 【 】已知一振蕩 環(huán)節(jié)的傳遞函數(shù)為 G(s)=12122 ?? TssT ?, 求當 T=5, ζ =、 、 、…、 時的幅頻特性曲線和對數(shù)幅頻相頻特性曲線。 W1=feedback(W,1)。 sgrid(z,wn)。 rlocus(num,den)。 rlocus(num,den)。 end 結(jié)果及打印曲線: 0 2 4 6 8 10 1200 . 20 . 40 . 60 . 811 . 21 . 41 . 61 . 8S t e p R e s p o n s eT im e ( s e c o n d s )Amplitude 0 2 4 6 8 10 124321012345Im p u ls e R e s p o n s eT im e ( s e c o n d s )Amplitude 6 se4? 11002?s K 【 】 試繪制如圖所示的系統(tǒng)的根軌跡( K 由 0→ +∞變化)。 hold on。 num=w2。 end (2)用直接求根法判別系統(tǒng)穩(wěn)定性: d=[1 20 9 100]。(d20)amp。0 2 4 0。0 2 4]。1 3]。 G=L1/L2 % 傳遞函數(shù) simplify(G) H1=G*1。 解題過程及程序: syms t。 den=[3 2 4 8]。 MATLAB的 CONTROL SYSTEM工具箱和 SIMULINK仿真軟件,分析系統(tǒng)的性能。作為自動化專業(yè)的學生很有必要學會應用這一強大的工具,并掌握利用 MATLAB對控制理論內(nèi)容進行分析和研究的技能,以達到加深對課堂上所講內(nèi)容理解的目的。設(shè)置《 自動控制理論 A》課程設(shè)計的目的是使 MATLAB成為學生的基本技能,熟悉 MATLAB這一解決具體工程問題的標準軟件,能熟練地應用 MATLAB軟件解決控制理論中的復雜和工程實際問題,并給以后的模糊控制理論、最優(yōu)控制理論和多變量控制理論等奠定基礎(chǔ)。 MATLAB軟件設(shè)計控制系統(tǒng)以滿足具體的性能指標要求。 解題過程及程序: num=[3 4 6]。 G2=tf(num,den,Ts) 結(jié)果及打印曲線: Transfer function: 3 z^2 + 4 z + 6 3 z^3 + 2 z^2 + 4 z + 8 Sampling time (seconds): 4 3 【 】 某系統(tǒng)在零初始條件的單位階躍響應為 h(t)=(1 tt ee ?? ?2 )1(t),試求系統(tǒng)傳遞函數(shù)及零初始條下的單位脈沖響應。 L2=1/s。 a=[2 4。1 3 5。1 3 5 0。 if((d10)amp。 else WARNDLG(‘The system is unstable’,’Stability Analysis’)。 w2=wn*wn。 step(num,den)。 hold on。 den=conv(den1,[100 1])。 den=conv([1 1],[