實用文檔研究生《現代控制理論及其應用》課程小論文一級倒立擺的建模與控制分析學院：機械工程學院班級：機研131姓名：尹潤豐學號：2013212020162014 年6月2日實用文檔目錄1.問題描述及狀態空間表達式建立............................-1-1.1問題描述.................................................................-1-1.2狀態空間表達式的建立.....................................................-1-1.2.1直線一級倒立擺的數學模型...........................................-1-1.2.2直線一級倒立擺系統的狀態方程......................................-5-2.應用MATLAB分析系統性能..................................-6-2.1直線一級倒立擺閉環系統穩定性分析.........................................-6-2.2系統可控性分析..........................................................-7-2.3系統可觀測性分析........................................................-8-3.應用matlab進行綜合設計.................................-8-3.1狀態反饋原理.............................................................-8-3.2全維狀態反饋觀測器和simulink仿真........................................-9-4.應用Matlab進行系統最優控制設計.........................-11-5.總結....................................................-13-實用文檔問題描述及狀態空間表達式建立1.1問題描述倒立擺是機器人技術、控制理論、計算機控制等多個領域、多種技術的有機結合，其被控系統本身又是一個絕對不穩定、高階次、多變量、強耦合的非線性系統，可以作為一個典型的控制對象對其進行研究。倒立擺系統作為控制理論研究中的一種比較理想的實驗手段，為自動控制理論的教學、實驗和科研構建一個良好的實驗平臺，以用來檢驗某種控制理論或方法的典型方案，促進了控制系統新理論、新思想的發展。下對于倒立擺系統，經過小心的假設忽略掉一些次要的因素后，它就是一個典型的運動的剛體系統，可以在慣性坐標系內應用經典力學理論建立系統的動力學方程。下面采用其中的牛頓—歐拉方法建立直線一級倒立擺系統的數學模型。1.2狀態空間表達式的建立 直線一級倒立擺的數學模型實用文檔圖1.1直線一級倒立擺系統1.1所示。本文中倒立擺系統描述中涉及的符號、物理意義及相關數值如表圖1.2是系統中小車的受力分析圖。其中，N和P為小車與擺桿相互作用力的水平和垂直方向的分量。實用文檔圖1.2系統中小車的受力分析圖圖1.3是系統中擺桿的受力分析圖。Fs是擺桿受到的水平方向的干擾力,Fh是擺桿受到的垂直方向的干擾力，合力是垂直方向夾角為α的干擾力Fg。圖1.3擺桿受力分析圖分析小車水平方向所受的合力，可以得到以下方程：MxFfxN11設擺桿受到與垂直方向夾角為α的干擾力Fg，可分解為水平方向、垂直方向的干擾力，所產生的力矩可以等效為在擺桿頂端的水平干擾力FS、垂直干擾力Fh產生的力矩。Fh12FSFgsinFgcos對擺桿水平方向的受力進行分析可以得到下面等式：2NFSmd2xlsin13dt實用文檔即：Nmxmlcosml2sinFfsin14對圖1.3擺桿垂直方向上的合力進行分析，可以得到下面方程：md2lcosPmgFh2l15dt即PmgFgcosmlsinml2cos1 6力矩平衡方程如下：FglsincosFglcossinPlsinNlcosI017代入P和N，得到方程：2Fglsincos2FglcossinIml2cos2mglsinml22sin2mlxcos018設，（φ是擺桿桿與垂直向上方向之間的夾角，單位是弧度），代入上式。假設φ<<1，則可進行近似處理：2cos1,sin,d0,cos21,sin2dt由于：I1ml23方程化為：2Fgsincos4mlmgmx193令：FfFgsincos，則19可化為：2Ff4mlmgmx1103110即是化簡后的直線一級倒立擺系統微分方程。帶入實際數據后，微分方程為：29.43x2Ff111m當忽略了F時，系統的微分方程如式（1-12）所示f29.43x112忽略干擾力后，直線一級倒立擺系統是單輸入二輸出的四階系統，考慮干擾力后，實用文檔直線一級倒立擺系統是二輸入二輸出的四階系統。其內部的4個狀態量分別是小車的位移x、小車的速度x、擺桿的角度θ、擺桿的角速度。系統輸出的觀測量為小車的位移x、擺桿的角度θ。其控制量為小車的加速度將微分方程（1-12）化為關于加速度輸入量和角度輸出量的傳遞函數：s3113Rss229.4 直線一級倒立擺系統的狀態方程實驗所使用的直線一級倒立擺系系統是加速度x作為系統的控制輸入，所以根據式（1-12）建立系統的狀態方程為：xxxx3g3x4l4l整理后得到系統狀態方程：x0100x000001xxx00010003g03g4l4lxyx1000x0x00100將實際參數代入得到一級倒立擺系統的狀態空間方程為：實用文檔x0100x0x0000x10001x00029.403xyx1000x00010x0010000000110000A001BC010D00000029.403應用MATLAB分析系統性能2.1直線一級倒立擺閉環系統穩定性分析構建如圖1.4所示閉環系統，則系統的閉環極點為（ -5.1381）、（5.1381）:圖1.4閉環系統結構圖由于有實部為正的極點，所以閉環系統不穩定，必須設計控制器使系統穩定。可以通過MATLABSimulink中對其進行仿真，判斷其穩定性。構建圖1.4所示系統的仿真程序e1，加入1m/s2的階躍信號實用文檔由上圖也能清楚的知道一級倒立擺系統是不穩定的。2.2 系統可控性分析系統的可控性可根據秩判據進行可控性判斷。 線性定常連續系統完全可控的充分必要條件是：rank(B AB An1B) n，其中 n為系統矩陣 A的階次，M (BAB An1B)為系統的可控性矩陣。matlab程序及運行結果如下：A=[0100;0000;0001;0029.40];B=[0;1;0;3];T=ctrb(A,B);rank(T)ans=4由于rank（Ic）=4，可見該系統是完全可控的。實用文檔2.3 系統可觀測性分析系統的可控性可根據秩判據進行可控性判斷。線性定常連續系統完全可控的充分必要條件是：CCArank N rank n或rank(CT ATCT(AT)2CT (AT)n1CT) nCAn1其中n為系數矩陣A的階次。matlab程序及運行結果如下：A=[0100;0000;0001;0029.40];C=[1000;0010];T0=obsv(A,C);rank(T0)ans=4由于rank（T0）=4，故該系統是可觀測的。應用matlab進行綜合設計3.1狀態反饋原理設ｎ維線性定常系統：x Ax Bu,y Cx其中x,u,y分別是n維、p維、q維向量；A、B、C分別是n*n維，n*p維，n*q維實數矩陣。狀態反饋系統的控制量u取為狀態x的線性函數：實用文檔u v Kx其中，v為p參考輸入向量，K為p*n維實反饋增益矩陣。加入狀態反饋后系統的結構圖如圖 3.1所示：圖3.1系統的全狀態反饋結構圖則系統狀態反饋的動態方程為：ABKxBv,yCx3.2全維狀態反饋觀測器和 simulink 仿真狀態反饋的的實現是利用狀態反饋使系統的閉環極點位于所希望的極點位置。而狀態反饋任意配置閉環極點的充分必要條件是被控系統可控。直線一級倒立擺系統是可控的。設系統期望極點為1234=2343i43i，則系統期望特征多項式為：a*ss1s2s3s4列寫狀態反饋系統的特征多項式：detSIABK令兩個特征多項式各項系數對應相等， 則可解出K陣。由matlab求出狀態反饋矩陣K，編程如下：A=[0100;0000;0001;0029.40];B=[0;1;0;3];K=acker(A,B,[-2-3-4+3i-4-3i])K=實用文檔-5.1020 -5.8844 35.1673 6.2948系統加入0.1m/s2的階躍輸入，在構成的狀態反饋調節器控制下，MATLAB中進行系統的階躍響應仿真，編程如下：A=[0100000000010029.40];B1=[0103];C=[1000010000100001];D1=[0000]';dt=0.005;ieof=801;fori=1:ieof;U(:,i)=[0.1];T(i)=i*dt;end;%%離散化op=[-2 % 期望極點-3-4+3i-4-3i];K=place(A,B1,op)Ak0=[(A-B1*K)];Bk0=[B1];Ck0=[C];Dk0=[D1];lqrop=eig(Ak0);x=[0000]';dt=0.005;%離散時[dA,dB]=c2d(Ak0,Bk0,dt);% 經離散化得到離散狀態方程Ak1=[(A-B1*K)];Bk1=[B1];Ck1=[C];Dk1=[D1];sys=ss(Ak1,Bk1,Ck1,Dk1);實用文檔[Y,X]=lsim(sys,U,T);plot(T,-Y),grid;legend('Cart','VCart','single','Vs');圖3.2 極點配置為[-2-3-4+3i-4-3i] 時的全狀態反饋仿真圖橫軸時間單位秒，從圖中可以看出，系統穩定。應用Matlab進行系統最優控制設計最優控制問題就是尋找一個控制系統的最優控制方案或最優控制規律，使系統能最優地達到預期的目標。對于線性連續系統，提出二次型目標函數：J1xT(tf)Sx(tf)1tfxTtQtxtuTtRtutdt22t0其中，R(t)正定，S及Q（t)半正定，且設它們為對稱矩陣， t0,tf固定。當tf趨近無窮時，在y(t) x(t)情況下，該問題即為無限時間輸出調節器問題。 此時穩態誤差項趨于零，在此題目中假設二次型最優控制性能指標為：J [xTQx uTRu]dt0實用文檔50000003000其中：Q050R=10000010Matlab編程如下：A=[0100;0000;0001;0029.40];B=[0;1;0;3];C=[1000;0100;0010;0001];D=[0000];Q=[500000;03000;00500;00010];R=1;[K,P,e]=lqr(A,B,Q,R)K=-22.3607-17.469770.104113.2462在simulink 下進行仿真模型的建立，如圖 4.1：圖4.1LQR仿真模型將K輸入后，進行仿真，結果如圖 4.2：實用文檔圖4.2LQR仿真結果由圖可見，在二次型最優控制下系統穩定性得到明顯改善。總結通過對一級倒立擺的分析