鍋爐過(guò)熱氣溫控制MATLAB及控制系統(tǒng)仿真

裝訂線 課程設(shè)計(jì)報(bào)告 題目：MATLAB及控制系統(tǒng)仿真課程設(shè)計(jì) 學(xué) 院 電子信息工程學(xué)院 學(xué)科門類 電氣信息類 專 業(yè) 自動(dòng)化 學(xué) 號(hào) 2012449107 姓 名 陳文華 指導(dǎo)教師 姜萍 2016年 1 月 16 日目 錄一 引言21.1 實(shí)驗(yàn)?zāi)康?1.2 實(shí)驗(yàn)內(nèi)容與要求21.2.1實(shí)驗(yàn)內(nèi)容21.2.2實(shí)驗(yàn)要求2二 倒立擺控制系統(tǒng)設(shè)計(jì)32.1倒立擺的簡(jiǎn)介32.2倒立擺的數(shù)學(xué)模型32.2.1本設(shè)計(jì)中所用到的各變量的取值及其意義32.2.2動(dòng)力學(xué)模型32.3模型轉(zhuǎn)化5三 基于狀態(tài)反饋的倒立擺系統(tǒng)設(shè)計(jì)63.1系統(tǒng)的開(kāi)環(huán)仿真63.1.1開(kāi)環(huán)仿真的系統(tǒng)Simulink結(jié)構(gòu)63.1.2開(kāi)環(huán)系統(tǒng)的分析73.2輸出反饋設(shè)計(jì)方法73.2.1輸出反饋仿真73.2.2輸出反饋系統(tǒng)的分析83.3狀態(tài)反饋設(shè)計(jì)83.3.1基于狀態(tài)反饋控制器的倒立擺設(shè)計(jì)過(guò)程83.3.2狀態(tài)反饋仿真93.3.3狀態(tài)反饋分析103.4全維狀態(tài)觀測(cè)器的倒立擺控制系統(tǒng)設(shè)計(jì)與仿真103.4.1基于全維狀態(tài)觀測(cè)器的倒立擺系統(tǒng)設(shè)計(jì)步驟103.4.2系統(tǒng)仿真10 3.4.3基于狀態(tài)觀測(cè)器的狀態(tài)反饋曲線分析11四 鍋爐過(guò)熱汽溫控制系統(tǒng)設(shè)計(jì)及仿真124.1蒸汽溫度控制的任務(wù)124.2影響蒸汽溫度的因素124.3蒸汽溫度系統(tǒng)開(kāi)環(huán)模型建立124.3.1減溫水量對(duì)蒸汽溫度的影響124.3.2動(dòng)態(tài)特性124.4蒸汽溫度控制系統(tǒng)設(shè)計(jì)124.4.1開(kāi)環(huán)系統(tǒng)動(dòng)態(tài)特性仿真及分析124.4.2開(kāi)環(huán)特性曲線分析134.5單回路控制系統(tǒng)134.5.1單回路控制系統(tǒng)仿真及分析134.5.2系統(tǒng)PID參數(shù)的整定134.5.3單回路控制系統(tǒng)仿真曲線分析154.6串級(jí)控制系統(tǒng)154.6.1串級(jí)控制系統(tǒng)仿真154.6.2系統(tǒng)PID參數(shù)的整定164.6.3串級(jí)系統(tǒng)響應(yīng)曲線分析18五 總結(jié)19附 錄20一 引言1.1 實(shí)驗(yàn)?zāi)康模?）加強(qiáng)學(xué)生對(duì)控制理論及控制系統(tǒng)的理解，熟練應(yīng)用計(jì)算機(jī)仿真常用算法和工具，完成控制系統(tǒng)計(jì)算機(jī)輔助設(shè)計(jì)的訓(xùn)練。（2）提高學(xué)生對(duì)控制系統(tǒng)的綜合及設(shè)計(jì)技能，擴(kuò)大學(xué)生的知識(shí)面，培養(yǎng)學(xué)生獨(dú)立分析問(wèn)題及解決問(wèn)題的能力，為以后從事實(shí)際控制系統(tǒng)的設(shè)計(jì)工作打下基礎(chǔ)。1.2 實(shí)驗(yàn)內(nèi)容與要求1.2.1實(shí)驗(yàn)內(nèi)容（1）基于觀測(cè)器的倒立擺控制系統(tǒng)設(shè)計(jì)及仿真（2）鍋爐過(guò)熱汽溫控制系統(tǒng)設(shè)計(jì)及仿真1.2.2實(shí)驗(yàn)要求（1）系統(tǒng)分析及數(shù)學(xué)模型建立（2）開(kāi)環(huán)系統(tǒng)仿真及動(dòng)態(tài)特性分析（3）控制方案設(shè)計(jì)及閉環(huán)系統(tǒng)仿真實(shí)驗(yàn)（4）實(shí)驗(yàn)結(jié)果分析二 倒立擺控制系統(tǒng)2.1倒立擺的簡(jiǎn)介倒立擺系統(tǒng)是一個(gè)復(fù)雜的、高度非線性的、不穩(wěn)定的高階系統(tǒng)，是學(xué)習(xí)和研究現(xiàn)代控制理論最合適的實(shí)驗(yàn)裝置。倒立擺的控制是控制理論應(yīng)用的一個(gè)典型范例,一個(gè)穩(wěn)定的倒立擺系統(tǒng)對(duì)于證實(shí)狀態(tài)空間理論的實(shí)用性是非常有用的。由于倒立擺本身是自不穩(wěn)定的系統(tǒng)，實(shí)驗(yàn)建模存在一定的困難。但是經(jīng)過(guò)假設(shè)忽略掉一些次要的因素后，倒立擺系統(tǒng)就是一個(gè)典型的運(yùn)動(dòng)的剛體系統(tǒng)，可以在慣性坐標(biāo)系內(nèi)應(yīng)用經(jīng)典力學(xué)理論建立系統(tǒng)的動(dòng)力學(xué)方程關(guān)系。在此，我們首先應(yīng)用動(dòng)力學(xué)方程建立一級(jí)倒立擺的非線性數(shù)學(xué)模型；采用小偏差線性化的方法在平衡點(diǎn)附近局部線性化得到線性化的數(shù)學(xué)模型；然后應(yīng)用狀態(tài)空間分析方法，采用狀態(tài)反饋為倒立擺系統(tǒng)建立穩(wěn)定的控制律；最后應(yīng)用狀態(tài)觀測(cè)器實(shí)現(xiàn)倒立擺系統(tǒng)的穩(wěn)定控制。2.2倒立擺的數(shù)學(xué)模型倒立擺示意圖如圖2-1所示，通過(guò)對(duì)小車施加一定的驅(qū)動(dòng)力，使倒立擺保持一定的位姿。圖2-1倒立擺示意圖2.2.1本設(shè)計(jì)中所用到的各變量的取值及其意義小車質(zhì)量M；m：小球的質(zhì)量；l：倒擺的桿長(zhǎng)；g：重力加速度；：表示倒擺偏離垂直方向的角度；u是小車受到的水平方向的驅(qū)動(dòng)力；2.2.2動(dòng)力學(xué)模型小球受力分析如圖2-2所示，其中表示小球的重心坐標(biāo)YFyFxGxl水平方向受到的合外力豎直方向受到的合力圖2-2小球受力分析示意圖通過(guò)受力分析，由牛頓第二運(yùn)動(dòng)定律，系統(tǒng)的運(yùn)動(dòng)滿足下面的方程：x軸方向：小球的重心坐標(biāo)滿足：整理后得：小球的力矩平衡方程： 整理可得：最后得到倒立擺系統(tǒng)的動(dòng)力學(xué)方程： 顯然該系統(tǒng)為明顯的非線性系統(tǒng)。但是對(duì)小車施加驅(qū)動(dòng)力的目的是要保持小球在垂直方向的姿態(tài)，因此，我們關(guān)注的是小球在垂直方向附近的動(dòng)態(tài)行為變化，為此將系統(tǒng)在該參考位(0)附近進(jìn)行線性化處理。 2.3模型轉(zhuǎn)化微分方程狀態(tài)方程由倒擺系統(tǒng)的動(dòng)力學(xué)模型取如下?tīng)顟B(tài)變量：可得到倒擺系統(tǒng)的狀態(tài)方程：2.4狀態(tài)方程的線性化：采用Jacobian 矩陣線性化模型，最終得到系統(tǒng)的線性化狀態(tài)方程為：假定系統(tǒng)的輸出為倒擺的角度和小車的x軸坐標(biāo)，則系統(tǒng)的輸出方程為：三 基于狀態(tài)反饋的倒立擺系統(tǒng)設(shè)計(jì)3.1系統(tǒng)的開(kāi)環(huán)仿真3.1.1開(kāi)環(huán)仿真的系統(tǒng)Simulink結(jié)構(gòu)開(kāi)環(huán)仿真的系統(tǒng)Simulink結(jié)構(gòu)圖如圖3-1所示圖3-1開(kāi)環(huán)仿真Simulink結(jié)構(gòu)圖運(yùn)行后觀察小車位置響應(yīng)曲線如圖3-2所示，小球角度響應(yīng)曲線如圖3-3所示。圖3-2cart pos響應(yīng)曲線圖3-3rod abgle響應(yīng)曲線3.1.2開(kāi)環(huán)系統(tǒng)的分析由圖3-2和圖3-3所示，小球的角度會(huì)隨著小車的位移的增大而增大，并不能自動(dòng)調(diào)整在平衡點(diǎn)附近來(lái)回?cái)[動(dòng)。可見(jiàn)開(kāi)環(huán)系統(tǒng)并不能維持系統(tǒng)的穩(wěn)定性。3.2輸出反饋設(shè)計(jì)方法3.2.1輸出反饋仿真輸出反饋結(jié)構(gòu)Simulink結(jié)構(gòu)圖如圖3-4所示圖3-4輸出反饋Simulink結(jié)構(gòu)圖運(yùn)行后系統(tǒng)波形倒擺的角度的響應(yīng)曲線如圖3-5，小車的位置的響應(yīng)曲線圖3-6所示。圖3-5倒擺的角度響應(yīng)曲線圖3-6小車的位置響應(yīng)曲線3.2.2輸出反饋系統(tǒng)的分析通過(guò)反復(fù)的調(diào)整和研究增益k1、k2對(duì)于系統(tǒng)誤差的敏感性，最終能夠穩(wěn)定系統(tǒng)。然而系統(tǒng)的動(dòng)態(tài)性能遠(yuǎn)不能讓人滿意，對(duì)于k1=-50,k2=-2，系統(tǒng)只是臨界穩(wěn)定，它仍在新的參考點(diǎn)附近反復(fù)震蕩。3.3狀態(tài)反饋設(shè)計(jì)3.3.1基于狀態(tài)反饋控制器的倒立擺設(shè)計(jì)過(guò)程（1） 系統(tǒng)能控性判別，應(yīng)用可控性判別矩陣CM=ctrb（A,B），再判斷該矩陣的秩rank（CM）=4，由開(kāi)環(huán)系統(tǒng)分析部分已經(jīng)得知系統(tǒng)狀態(tài)完全能控。（2） 閉環(huán)系統(tǒng)的極點(diǎn)配置。根據(jù)系統(tǒng)的動(dòng)態(tài)性能，確定閉環(huán)系統(tǒng)的期望幾點(diǎn)clp，clp=-1.5+3.0j -1.5-3.0j -5 -4。（3） 確定反饋增益。應(yīng)用MATLAB的place函數(shù)Ks=place（A,B,clp），確定反饋增益Ks，Ks=-432.6154 -176.2944 -89.5077 -64.1472。（4） 系統(tǒng)設(shè)計(jì)。由狀態(tài)反饋方框圖可得系統(tǒng)的狀態(tài)空間表達(dá)式為 此時(shí)，系統(tǒng)矩陣為，(其中為反饋增益矩陣)，控制矩陣為（其中）,因?yàn)閷?duì)小車的控制要求靜態(tài)終值，所以。此時(shí)的系統(tǒng)設(shè)計(jì)3.3.2狀態(tài)反饋仿真狀態(tài)反饋結(jié)構(gòu)Simulink結(jié)構(gòu)圖如圖3-7所示圖3-7狀態(tài)反饋Simulink結(jié)構(gòu)圖小車位置和狀態(tài)變量的響應(yīng)曲線如圖3-8所示圖3-8小車位置和狀態(tài)變量的響應(yīng)曲線3.3.3狀態(tài)反饋分析從響應(yīng)曲線可以看出，小車開(kāi)始沿x軸正向移動(dòng)，大約3s后靜止在x=1m處。并且此時(shí)所有的狀態(tài)變量都趨于0，x(t)趨于平衡點(diǎn)。3.4全維狀態(tài)觀測(cè)器的倒立擺控制系統(tǒng)設(shè)計(jì)與仿真 3.4.1基于全維狀態(tài)觀測(cè)器的倒立擺系統(tǒng)設(shè)計(jì)步驟（1） 系統(tǒng)能觀性判別。應(yīng)用客觀性判別矩陣N=obsv（A,C），判別該矩陣的秩rank（N）=4，所以系統(tǒng)狀態(tài)完全能觀。（2） 狀態(tài)觀測(cè)器閉環(huán)極點(diǎn)配置。適當(dāng)選擇觀測(cè)器的極點(diǎn)，使觀測(cè)器的動(dòng)態(tài)速度是系統(tǒng)的兩倍以上，所觀測(cè)的極點(diǎn)op=2*clp。（3） 指定極點(diǎn)的觀測(cè)器增益L。同樣應(yīng)用place函數(shù)：G=place（A,B,op），G=G,G=1.0e+00.*-2.882 -9.8401 0.024 0.2382。（4） 系統(tǒng)設(shè)計(jì)。其中3.4.2系統(tǒng)仿真基于狀態(tài)觀測(cè)器的狀態(tài)反饋Simulink結(jié)構(gòu)圖如圖3-9所示圖3-9基于狀態(tài)觀測(cè)器的狀態(tài)反饋Simulink結(jié)構(gòu)圖仿真結(jié)果狀態(tài)曲線圖如圖3-10，圖3-11顯示了系統(tǒng)狀態(tài)與觀測(cè)器得到的估計(jì)狀態(tài)之間的誤差曲線3-10小車位置和倒擺角度響應(yīng)曲線3-11狀態(tài)變量的誤差曲線3.4.3基于狀態(tài)觀測(cè)器的狀態(tài)反饋曲線分析從響應(yīng)曲線可以看出，小車開(kāi)始沿x軸正向移動(dòng)，并且此時(shí)所有的狀態(tài)變量都趨于0，x(t)趨于平衡點(diǎn)。四 鍋爐過(guò)熱汽溫控制系統(tǒng)設(shè)計(jì)及仿真4.1蒸汽溫度控制的任務(wù)鍋爐出口過(guò)熱蒸汽溫度是蒸汽的重要質(zhì)量指標(biāo)，是整個(gè)鍋爐汽水通道中溫度最高的，直接關(guān)系到設(shè)備的安全和系統(tǒng)的生產(chǎn)效率。過(guò)高，使金屬?gòu)?qiáng)度降低，影響設(shè)備安全；過(guò)低，使全廠熱效率顯著下降，每下降 5 oC 使熱效率下降 1%。鍋爐過(guò)熱蒸汽溫度控制的基本任務(wù)就是維持過(guò)熱器出口溫度在允許范圍內(nèi)，保護(hù)設(shè)備安全，并使生產(chǎn)過(guò)程經(jīng)濟(jì)、高效的持續(xù)運(yùn)行。4.2影響蒸汽溫度的因素（1） 減溫水量 QW （控制量）（2）蒸汽流量 D（3）煙氣熱量 QH4.3蒸汽溫度系統(tǒng)開(kāi)環(huán)模型建立4.3.1減溫水量對(duì)蒸汽溫度的影響過(guò)熱器具有多分布參數(shù)的對(duì)象，可以把管內(nèi)蒸汽和金屬管壁看作多個(gè)單容對(duì)象串聯(lián)組成的多容對(duì)象。當(dāng)減溫水流量發(fā)生變化后，需要通過(guò)這些串聯(lián)單容對(duì)象，最終引起出口蒸汽溫度變化。減溫器距離出口越遠(yuǎn)延遲就越大。4.3.2動(dòng)態(tài)特性本實(shí)驗(yàn)采用的動(dòng)態(tài)特性的高階模型為負(fù)荷為100%，動(dòng)態(tài)特性為（1）導(dǎo)前區(qū)： （2）惰性區(qū)4.4蒸汽溫度控制系統(tǒng)設(shè)計(jì)4.4.1開(kāi)環(huán)系統(tǒng)動(dòng)態(tài)特性仿真及分析開(kāi)環(huán)系統(tǒng)動(dòng)態(tài)特性如圖4-1所示圖4-1開(kāi)環(huán)系統(tǒng)動(dòng)態(tài)特性Simulink結(jié)構(gòu)圖運(yùn)行后開(kāi)環(huán)動(dòng)態(tài)特性曲線如圖4-2所示圖4-2開(kāi)環(huán)動(dòng)態(tài)特性曲線4.4.2開(kāi)環(huán)特性曲線分析由圖4-2可知，系統(tǒng)在250秒左右穩(wěn)定在3.8。4.5單回路控制系統(tǒng)4.5.1單回路控制系統(tǒng)仿真及分析單回路控制系統(tǒng)仿真如圖4-3所示圖4-3單回路控制系統(tǒng)Simulink結(jié)構(gòu)圖4.5.2系統(tǒng)PID參數(shù)的整定（1）取Ti=，Td=0。 P較大（Kp較?。┕r穩(wěn)定時(shí)投入自動(dòng)；（2）逐漸減小P（或增大Kp）每改變一次都給系統(tǒng)施加一次定值階躍，觀察輸出曲線，直至出現(xiàn)等幅振蕩（四，五次即可），如圖4-4所示，記錄此時(shí)的Kp=0.61，Pm=1/Kp=1.64，測(cè)出振蕩周期Tm=150；圖4-4等幅震蕩曲線（3）PID參數(shù)整定經(jīng)驗(yàn)公式計(jì)算PTiPI2.2Pm=0.5720.85Tm=10.2根據(jù)整定的參數(shù)，進(jìn)行PID參數(shù)設(shè)置如圖4-5所示圖4-5參數(shù)設(shè)置得到仿真特性曲線如圖4-6所示圖4-6PI調(diào)節(jié)特性曲線（4） 可見(jiàn)振蕩較厲害，響應(yīng)曲線品質(zhì)不夠理想，在此基礎(chǔ)上繼續(xù)調(diào)整，增大積分時(shí)間、減小比例系數(shù)（均為增強(qiáng)穩(wěn)定性）并嘗試加上微分作用。參數(shù)整定如圖4-7所示圖4-7參數(shù)設(shè)置輸出仿真結(jié)果，如圖4-8所示圖4-8單回路控制系統(tǒng)仿真特性曲線4.5.3單回路控制系統(tǒng)仿真曲線分析由圖4-8可見(jiàn)，控制效果大大改善，有效抑制了超調(diào)并增強(qiáng)穩(wěn)定性，快速達(dá)到平衡。4.6串級(jí)控制系統(tǒng)4.6.1串級(jí)控制系統(tǒng)仿真串級(jí)控制系統(tǒng)仿真Simulink結(jié)構(gòu)圖如圖4-9所示圖4-9串級(jí)控制系統(tǒng)仿真Simulink結(jié)構(gòu)圖4.6.2系統(tǒng)PID參數(shù)的整定（1）取Ti=，Td=0。 P較大（Kp較?。┕r穩(wěn)定時(shí)投入自動(dòng)；（2）逐漸減小P（或增大Kp）每改變一次都給系統(tǒng)施加一次定值階躍，觀察輸出曲線，直至出現(xiàn)等幅振蕩（四，五次即可），如圖4-10所示，記錄此時(shí)的Kp=2.4，Pm=1/Kp=0.42，測(cè)出振蕩周期Tm=100；圖4-10等幅震蕩曲線（3）PID參數(shù)整定經(jīng)驗(yàn)公式計(jì)算PTiPI2.2Pm0.85Tm根據(jù)整定的參數(shù)，進(jìn)行PID參數(shù)設(shè)置如圖4-11所示圖4-11參數(shù)設(shè)置得到仿真特性曲線如圖4-12所示圖4-12PI調(diào)節(jié)特性曲線可見(jiàn)振蕩較厲害，響應(yīng)曲線品質(zhì)不夠理想，在此基礎(chǔ)上繼續(xù)調(diào)整，增大積分時(shí)間、減小比例系數(shù)（均為增強(qiáng)穩(wěn)定性）調(diào)整后的PID參數(shù)如圖4-13所示圖4-13調(diào)整后的PID參數(shù)串級(jí)控制系統(tǒng)仿真特性曲線如圖4-14所示圖4-14串級(jí)控制系統(tǒng)仿真特性曲線4.6.3串級(jí)系統(tǒng)的響應(yīng)曲線分析主控制器的輸出即副控制器的給定，而副控制器的輸出直接送往控制閥。主控制器的給定值是由工藝規(guī)定的，是一個(gè)定制，因此，主環(huán)是一個(gè)定值控制系統(tǒng)；而副控制器的給定值是由主控制器的輸出提供的，它隨主控制器輸出變化而變化，因此，副環(huán)是一個(gè)隨動(dòng)控制系統(tǒng)。串級(jí)控制系統(tǒng)中，兩個(gè)控制器串聯(lián)工作，以主控制器為主導(dǎo)，保證主變量穩(wěn)定為目的，兩個(gè)控制器協(xié)調(diào)一致，互相配合。若干擾來(lái)自副環(huán)，副控制器首先進(jìn)行“粗調(diào)”，主控制器再進(jìn)一步進(jìn)行“細(xì)調(diào)”。因此控制質(zhì)量?jī)?yōu)于簡(jiǎn)單控制系統(tǒng)。串級(jí)控制有以下優(yōu)點(diǎn) 由于副回路的存在，減小了對(duì)象的時(shí)間常數(shù)，縮短了控制通道，使控制作用更加及時(shí)； 對(duì)二次干擾具有很強(qiáng)的克服能力，對(duì)客服一次干擾的能力也有一定的提高； 對(duì)負(fù)荷或操作條件的變化有一定的自適應(yīng)能力。五 總結(jié)附 錄倒立擺的.m文件的程序：close all,clear allM=2.0;%小車的質(zhì)量m=0.1;%小球的質(zhì)量l=0.5;%擺桿的長(zhǎng)度g=9.81;%重力加速度%線性化模型的狀態(tài)空間矩陣A=0 1 0 0;(M+m)*g/(M*l) 0 0 0;0 0 0 1;-m*g/M 0 0 0B=0;-1/M/l;0;1/MC=1 0 0 0;0 0 1 0D=0;0ev=eig(A)CM=ctrb(A,B)%輸出秩=4，滿秩，完全能控rank(CM)clp=-1.5+3.0j -1.5-3.0j -5.0 -4.0;Ks=place(A,B,clp)eig(A-B*Ks)%驗(yàn)證閉環(huán)特征值Nr=-1/(C*inv(A-B*Ks)*B)%計(jì)算穩(wěn)態(tài)誤差set(0,showHiddenHandles,on);set(gcf,menubar,figure);set(0,showHiddenHandles,on);set(gcf,menubar,figure);op=2*clp%觀測(cè)器的速度是閉環(huán)系統(tǒng)的2倍G=place(A,C,op)G=Gset(0,showHiddenHandles,on);set(gcf,menubar,figure);