課程設計報告題目：題目：MATLA及控制系統仿真課程設計學院電子信息工程學院學科門類電氣信息類專業自動化學號2012449107姓名陳文華指導教師姜萍2016年1月16日河北大學2012級自動化MATLAE及控制系統仿真課程設計目錄TOC\o"1-5"\h\z一引言 2..\o"CurrentDocument"1.1實驗目的 2.\o"CurrentDocument"1.2實驗內容與要求 2.\o"CurrentDocument"1.2.1實驗內容 2.\o"CurrentDocument"1.2.2實驗要求 2.\o"CurrentDocument"二倒立擺控制系統設計 3.\o"CurrentDocument"2.1倒立擺的簡介 3.\o"CurrentDocument"2.2倒立擺的數學模型 3.\o"CurrentDocument"2.2.1本設計中所用到的各變量的取值及其意義 3\o"CurrentDocument"2.2.2動力學模型 3.\o"CurrentDocument"2.3模型轉化 5.\o"CurrentDocument"三基于狀態反饋的倒立擺系統設計 6\o"CurrentDocument"3.1系統的開環仿真 6.\o"CurrentDocument"3.1.1開環仿真的系統Simulink結構 6\o"CurrentDocument"3.1.2開環系統的分析 7.\o"CurrentDocument"3.2輸出反饋設計方法 7.\o"CurrentDocument"3.2.1輸出反饋仿真 7.\o"CurrentDocument"3.2.2輸出反饋系統的分析 8\o"CurrentDocument"3.3狀態反饋設計 8.\o"CurrentDocument"3.3.1基于狀態反饋控制器的倒立擺設計過程 8\o"CurrentDocument"3.3.2狀態反饋仿真 9.\o"CurrentDocument"333狀態反饋分析 10\o"CurrentDocument"3.4全維狀態觀測器的倒立擺控制系統設計與仿真 10\o"CurrentDocument"3.4.1基于全維狀態觀測器的倒立擺系統設計步驟 10\o"CurrentDocument"3.4.2系統仿真 10\o"CurrentDocument"3.4.3基于狀態觀測器的狀態反饋曲線分析 11\o"CurrentDocument"四鍋爐過熱汽溫控制系統設計及仿真 12\o"CurrentDocument"4.1蒸汽溫度控制的任務 12\o"CurrentDocument"4.2影響蒸汽溫度的因素 12\o"CurrentDocument"4.3蒸汽溫度系統開環模型建立 12\o"CurrentDocument"4.3.1減溫水量對蒸汽溫度的影響 12\o"CurrentDocument"4.3.2動態特性 12\o"CurrentDocument"4.4蒸汽溫度控制系統設計 12\o"CurrentDocument"4.4.1開環系統動態特性仿真及分析 12\o"CurrentDocument"4.4.2開環特性曲線分析 1.3\o"CurrentDocument"4.5單回路控制系統 13\o"CurrentDocument"4.5.1單回路控制系統仿真及分析 13\o"CurrentDocument"4.5.2系統PID參數的整定 1.3\o"CurrentDocument"4.5.3單回路控制系統仿真曲線分析 15\o"CurrentDocument"4.6串級控制系統 15\o"CurrentDocument"4.6.1串級控制系統仿真 15\o"CurrentDocument"4.6.2系統PID參數的整定 1.6\o"CurrentDocument"4.6.3串級系統響應曲線分析 18五總結 18附錄 19引言1.1實驗目的（1） 加強學生對控制理論及控制系統的理解，熟練應用計算機仿真常用算法和工具，完成控制系統計算機輔助設計的訓練。（2） 提高學生對控制系統的綜合及設計技能，擴大學生的知識面，培養學生獨立分析問題及解決問題的能力，為以后從事實際控制系統的設計工作打下基礎。1.2實驗內容與要求1.2.1實驗內容（1） 基于觀測器的倒立擺控制系統設計及仿真（2） 鍋爐過熱汽溫控制系統設計及仿真1.2.2實驗要求（1） 系統分析及數學模型建立（2） 開環系統仿真及動態特性分析（3） 控制方案設計及閉環系統仿真實驗（4） 實驗結果分析二倒立擺控制系統2.1倒立擺的簡介倒立擺系統是一個復雜的、高度非線性的、不穩定的高階系統，是學習和研究現代控制理論最合適的實驗裝置。倒立擺的控制是控制理論應用的一個典型范例,一個穩定的倒立擺系統對于證實狀態空間理論的實用性是非常有用的。由于倒立擺本身是自不穩定的系統，實驗建模存在一定的困難。但是經過假設忽略掉一些次要的因素后，倒立擺系統就是一個典型的運動的剛體系統，可以在慣性坐標系內應用經典力學理論建立系統的動力學方程關系。在此，我們首先應用動力學方程建立一級倒立擺的非線性數學模型；采用小偏差線性化的方法在平衡點附近局部線性化得到線性化的數學模型；然后應用狀態空間分析方法，采用狀態反饋為倒立擺系統建立穩定的控制律；最后應用狀態觀測器實現倒立擺系統的穩定控制。使倒立擺保持2.2倒立擺的數學模型使倒立擺保持倒立擺示意圖如圖2-1所示，通過對小車施加一定的驅動力,定的位姿圖2-1倒立擺示意圖2.2.1本設計中所用到的各變量的取值及其意義小車質量Mm小球的質量；I:倒擺的桿長；g：重力加速度；9:表示倒擺偏離垂直方向的角度；u是小車受到的水平方向的驅動力；2.2.2動力學模型小球受力分析如圖2-2所示，其中（Xc,yj表示小球的重心坐標通過受力分析，由牛頓第二運動定律，系統的運動滿足下面的方程:x軸方向：2 2d dM—2x m2 xG=udt2 dt2小球的重心坐標滿足： Mxm（xlsinn）二udt2 dt2整理后得：（Mm）x-ml（sinv）『?ml（cos打-u小球的力矩平衡方程：(FxCOOSv)l-(FySinn)l=(mgsinv)lFx=m—2冷=m[x-1(sin二)^2I(cosp]dtd2 2 . ..Fy=m2yG=m[T(cosd)v-1(sinr)^]dtxG=xlsinvyG=lcosvmxcosv-ml(sin二cos^)二2ml(cos2二)亠ml(sin二cosd)^2ml(sin-mgsin整理可得：mxcos：mlv-mgsin^最后得到倒立擺系統的動力學方程：2(M+m)x—ml(sin日)日+ml(cos8)T=uimxcos日+ml°=mgsin日

顯然該系統為明顯的非線性系統。但是對小車施加驅動力的目的是要保持小球在垂直方向的姿態，因此，我們關注的是小球在垂直方向附近的動態行為變化，為此將系統在該參考位(9二0)附近進行線性化處理。2.3模型轉化微分方程一狀態方程由倒擺系統的動力學模型(Mm)x—mg(sin對，ml(cos小-u一mxcosvmlv-mgsinvZ1取如下狀態變量： 乙-v,Z2-V-N,Z3二X,乙=X=zZ1可得到倒擺系統的狀態方程：一 J2ucos-(M m)gsinz1 ml(cos召sinz1)z2dz=dtd

dt2mlcosw-dz=dtd

dt2mlcosw-(Mm)lX■XJZ4uml(sinzjz22—mgcosz1sinz12Mm-mcosz12.4狀態方程的線性化:采用Jacobian矩陣線性化模型,「02.4狀態方程的線性化:采用Jacobian矩陣線性化模型,「0(Mm)gdz

dtMl0最終得到系統的線性化狀態方程為:-0〕_1IMl0-M」假定系統的輸出為倒擺的角度和小車的x軸坐標，則系統的輸出方程為:假定系統的輸出為倒擺的角度和小車的飛1飛1一10001日1=CZ「0010一1x!

三基于狀態反饋的倒立擺系統設計3.1系統的開環仿真3.1.1開環仿真的系統Simulink結構開環仿真的系統Simulink結構圖如圖3-1所示運行后觀察小車位置響應曲線如圖3-2所示，小球角度響應曲線如圖3-3所示。

3.1.2開環系統的分析由圖3-2和圖3-3所示，小球的角度會隨著小車的位移的增大而增大， 并不能自動調整在平衡點附近來回擺動。可見開環系統并不能維持系統的穩定性。3.2輸出反饋設計方法3.2.1輸出反饋仿真輸出反饋結構Simulink結構圖如圖3-4所示運行后系統波形倒擺的角度的響應曲線如圖 3-5，小車的位置的響應曲線圖3-6所示。

3.2.2輸出反饋系統的分析通過反復的調整和研究增益k1、k2對于系統誤差的敏感性，最終能夠穩定系統。然而系統的動態性能遠不能讓人滿意，對于k仁-50,k2=-2，系統只是臨界穩定，它仍在新的參考點附近反復震蕩。3.3狀態反饋設計3.3.1基于狀態反饋控制器的倒立擺設計過程（1）系統能控性判別，應用可控性判別矩陣 CM=ctrb（A,B），再判斷該矩陣的秩rank（CM=4，由開環系統分析部分已經得知系統狀態完全能控。（2） 閉環系統的極點配置。根據系統的動態性能，確定閉環系統的期望幾點clp，clp=[-1.5+3.0j-1.5-3.0j-5-4] 。（3）確定反饋增益。應用MATLAB勺place函數Ks=place（A,B,clp），確定反饋增益Ks,Ks=[-432.6154-176.2944-89.5077-64.1472]系統設計。由狀態反饋方框圖可得系統的狀態空間表達式為x=AxBu=AxB(Nrv-Ksx)=(A-BKs)xBNr此時，系統矩陣為A-BKs，(其中Ks為反饋增益矩陣)，控制矩陣為BNr(其中Nr?1/C2(A—BKs)」B),因為對小車的控制要求靜態終值Iimx3(t)=1，所t—JpC以Iimc2(sl-ABKS)JBNr=1。此時的系統設計sys二ss(A-B*Ks,B*Nr,C,D)s03.3.2狀態反饋仿真狀態反饋結構Simulink結構圖如圖3-7所示小車位置和狀態變量的響應曲線如圖 3-8所示

333狀態反饋分析從響應曲線可以看出，小車開始沿 x軸正向移動，大約3s后靜止在x=1m處。并且此時所有的狀態變量都趨于0，x(t)趨于平衡點。3.4全維狀態觀測器的倒立擺控制系統設計與仿真3.4.1基于全維狀態觀測器的倒立擺系統設計步驟系統能觀性判別。應用客觀性判別矩陣 N=obsv(A,C)，判別該矩陣的秩rank(N)=4，所以系統狀態完全能觀。狀態觀測器閉環極點配置。適當選擇觀測器的極點，使觀測器的動態速度是系統的兩倍以上，所觀測的極點op=2*clp。指定極點的觀測器增益L。同樣應用place函數:G=place(A'',op)，TOC\o"1-5"\h\zG=G,G=1.0e+00.*[-2.882-9.84010.0240.2382] 。系統設計。A A廠x=Ax+bu=Ax+b(v-Ksx)=Ax-bKsx+bvA A A A A] x=(A—GC)x+Gy+bu=(A—GC)x+GCx+b(v—Ksx)=(A—GC—bKs)x十GCx+bv<y=cx其中A11=A,A12 bKs,A21=A-GC-bKs.B^ -bNr3.4.2系統仿真基于狀態觀測器的狀態反饋Simulink結構圖如圖3-9所示

仿真結果狀態曲線圖如圖3-10,圖3-11顯示了系統狀態與觀測器得到的估計狀態之間的誤差曲線343基于狀態觀測器的狀態反饋曲線分析從響應曲線可以看出，小車開始沿x軸正向移動，并且此時所有的狀態變量都趨于0,x(t)趨于平衡點。四鍋爐過熱汽溫控制系統設計及仿真4.1蒸汽溫度控制的任務鍋爐出口過熱蒸汽溫度是蒸汽的重要質量指標，是整個鍋爐汽水通道中溫度最高的，直接關系到設備的安全和系統的生產效率。過高，使金屬強度降低，影響設備安全；過低，使全廠熱效率顯著下降，每下降5oC使熱效率下降1%。鍋爐過熱蒸汽溫度控制的基本任務就是維持過熱器出口溫度在允許范圍內， 保護設備安全，并使生產過程經濟、高效的持續運行。4.2影響蒸汽溫度的因素（1）減溫水量QW（控制量）（2）蒸汽流量D（3）煙氣熱量QH4.3蒸汽溫度系統開環模型建立4.3.1減溫水量對蒸汽溫度的影響過熱器具有多分布參數的對象，可以把管內蒸汽和金屬管壁看作多個單容對象串聯組成的多容對象。當減溫水流量發生變化后，需要通過這些串聯單容對象,最終引起出口蒸汽溫度變化。減溫器距離出口越遠延遲就越大。4.3.2動態特性本實驗采用的動態特性的高階模型為負荷為100%動態特性為（1）導前區:W2(s)（1）導前區:W2(s)二1.58(14s1)2（2）惰性區W（s）二2.45(15.8s1)44.4蒸汽溫度控制系統設計4.4.1開環系統動態特性仿真及分析開環系統動態特性如圖4-1所示Dmr Mi7*??Ji***1,?aInlOuMInlOuMSleplSubs怦m3—圖4-1開環系統動態特性Simulink結構圖

4.4.2開環特性曲線分析由圖4-2可知，系統在250秒左右穩定在3.84.5單回路控制系統4.5.1單回路控制系統仿真及分析單回路控制系統仿真如圖4-3所示4.5.2系統PID參數的整定（1） 取Ti=%,Td=0。P較大（Kp較小）工況穩定時投入自動；（2） 逐漸減小P（或增大Kp）每改變一次都給系統施加一次定值階躍，觀察輸出曲線，直至出現等幅振蕩（四，五次即可），如圖4-4所示，記錄此時的Kp=0.61,Pm=1/Kp=1.64,測出振蕩周期Tm=150

r--——IFAnAnnuv■AHTJ—JTJ—tf 111圖4-4等幅震蕩曲線r--——IFAnAnnuv■AHTJ—JTJ—tf 111圖4-4等幅震蕩曲線ZHMnl—TH! LrPID參數整定經驗公式計算PTiPI2.2Pm=0.5720.85Tm=10.2根據整定的參數，進行PID參數設置如圖4-5所示11E(S)u(s)=p(1tis)E(S)

E(S)可見振蕩較厲害，響應曲線品質不夠理想，在此基礎上繼續調整，增大積分時間、減小比例系數(均為增強穩定性)并嘗試加上微分作用。參數整定如圖4-7所示輸出仿真結果，如圖4-8所示4.5.3單回路控制系統仿真曲線分析由圖4-8可見，控制效果大大改善，有效抑制了超調并增強穩定性，快速達到平衡。4.6串級控制系統4.6.1串級控制系統仿真串級控制系統仿真Simulink結構圖如圖4-9所示

462系統PID參數的整定取Ti=%,Td=O。P較大(Kp較小)工況穩定時投入自動；逐漸減小P(或增大Kp)每改變一次都給系統施加一次定值階躍，觀察輸出曲線，直至出現等幅振蕩(四，五次即可)，如圖 4-10所示，記錄此時的Kp=2.4，Pm=1/Kp=0.42,測出振蕩周期Tm=100E(S)PIDE(S)經驗公式計算PTiPI2.2Pm0.85Tm根據整定的參數，進行PID參數設置如圖4-11所示1TiS

得到仿真特性曲線如圖4-12所示可見振蕩較厲害，響應曲線品質不夠理想，在此基礎上繼續調整，增大積分時間、減小比例系數（均為增強穩定性）調整后的 PID參數如圖4-13所示

463串級系統的響應曲線分析主控制器的輸出即副控制器的給定，而副控制器的輸出直接送往控制閥。主控制器的給定值是由工藝規定的，是一個定制，因此，主環是一個定值控制系統;而副控制器的給定值是由主控制器的輸出提供的，它隨主控制器輸出變化而變化，因此，副環是一個隨動控制系統。串級