PID控制算法的MATLAB仿真研究課程設(shè)計(jì)報(bào)告

1、鄧M上學(xué)電氣工程學(xué)院課程設(shè)計(jì)報(bào)告課程名稱：計(jì)算機(jī)控制技術(shù)設(shè)計(jì)題目：PID控制算法的MATLA防真研究專(zhuān)業(yè)班級(jí)：學(xué)號(hào)：姓名：時(shí)間： 2015 年12月24-31日以下由指導(dǎo)教師填寫(xiě)分項(xiàng)成績(jī)：出勤成品答辯及考核 一總成績(jī)：總分 成績(jī)指導(dǎo)教師（簽名）：目錄一、任務(wù)提出 1二、具體內(nèi)容及實(shí)現(xiàn)功能 12.1 PID參數(shù)整定 12.2 動(dòng)態(tài)特性參數(shù)法（Ziegler-Nichols 整定公式）簡(jiǎn)介12.3 按照Z(yǔ)-N公式設(shè)計(jì)PID校正器 22.4 模型失配時(shí)PID控制器的控制效果仿真 42.4.1 被控對(duì)象比例系數(shù)改變時(shí)的控制效果 52.4.2 被控對(duì)象時(shí)間常數(shù)改變時(shí)的控制效果 62.4.3 純滯后時(shí)間改

2、變時(shí)的控制效果 72.5 執(zhí)行機(jī)構(gòu)非線性對(duì)PID控制器控制效果的影響 82.6 擾動(dòng)作用對(duì)PID控制器控制效果的影響 92.6.1 反饋回路中的擾動(dòng)對(duì)控制效果的影響 92.6.2 控制器輸出位置的擾動(dòng)信號(hào)對(duì)控制效果的影響 10三、總結(jié) 10四、參考文獻(xiàn) 12五、附錄 13PID控制算法的 MATLA昉真研究任務(wù)提出采用帶純滯后的一階慣性環(huán)節(jié)的系統(tǒng)作為被控對(duì)象模型，傳遞函數(shù)為：G（s）=需？ ° ?= 633?7?平?？ .MATLA斷真框圖如圖1所示:圖1-1仿真控制系統(tǒng)二、具體內(nèi)容及實(shí)現(xiàn)功能2.1 PID參數(shù)整定PID控制器調(diào)節(jié)參數(shù)Kp Ki、Kd的整定對(duì)控制系統(tǒng)能否得到較 好的控

3、制效果是至關(guān)重要的。常用的PID參數(shù)整定方法有理論整定法（如Z-N法）和實(shí)驗(yàn)整定法（如擴(kuò)充臨界比例度法、試湊法等），也可以采用模糊自適應(yīng)參數(shù)整定、遺傳算法參數(shù)整定等新型的PID參數(shù)整定方法。在本仿真中采用了動(dòng)態(tài)參數(shù)整定法（ZN法）進(jìn)行參數(shù)整定。2.2 動(dòng)態(tài)特性參數(shù)法（Ziegler-Nichols整定公式）簡(jiǎn)介對(duì)于被控廣義對(duì)象為帶延遲的一階慣性環(huán)節(jié)的系統(tǒng)， 即傳遞函數(shù) 表達(dá)式為式（2-1）的系統(tǒng)，具PID控制的參數(shù)值可以用一組經(jīng)驗(yàn)公 式來(lái)計(jì)算。已知被控對(duì)象的傳遞函數(shù)模型K、T、p三個(gè)參數(shù)時(shí)，整定PID調(diào)節(jié)器參數(shù)的計(jì)算公式如表2-1所示G(s) =(2-1)KpTiTdPT/(Kr )PI0.9

4、T/(Kr )3.3 rPID1.2T/(Kr )2.2 r0.5 TK ? T ?T?+1 一表2-1調(diào)節(jié)器ZN整定公式2.3按照Z(yǔ)-N公式設(shè)計(jì)PID校正器將被控對(duì)象中的參數(shù)帶入表2-1中的公式可得參數(shù)如下:1.2?630Kp=30?60=0.42Ti=2.2*60=132Td=0.5*60=30選擇采樣時(shí)間Ts=1s,則可計(jì)算相應(yīng)的控制系數(shù)如下:比例系數(shù)為Kp=0.42比例系數(shù)為 Ki=?0.0032?圖2-1單位階躍響應(yīng)曲線1由該響應(yīng)曲線可以看出: 單位階躍響應(yīng)曲線上有滯后；此響應(yīng)曲線上毛刺較多、過(guò)渡過(guò)程不平滑；響應(yīng)的超調(diào)量過(guò)大，振蕩次數(shù)較多。針對(duì)以上問(wèn)題，分析如下：階躍響應(yīng)曲線的滯后是

5、由于對(duì)象所具有的純滯后特性所引起的， 采用單純地PID控制無(wú)法克服這一滯后特性；過(guò)渡過(guò)程不平滑，是由于采樣時(shí)間選擇的不合適。由相應(yīng)的控制 理論知識(shí)，可通過(guò)減小采樣周期以保證過(guò)渡過(guò)程足夠平滑；由于采樣時(shí)間的減小，勢(shì)必減弱系統(tǒng)的穩(wěn)定性。因此可增大積分 時(shí)間常數(shù)Ti ,以減弱這種影響；可通過(guò)減小比例系數(shù)Kp和微分時(shí)間常數(shù)Td減小系統(tǒng)的超調(diào)量綜合了上述的分析，在原有系數(shù)的基礎(chǔ)上多次試湊，選取采樣時(shí)間Ts=0.2 比例系數(shù)Kp=0.25積分時(shí)間常數(shù)Ti=125 微分時(shí)間常數(shù)Td=2.4這時(shí)，相應(yīng)的控制系數(shù)為:比例系數(shù)為Kp=0.25比例系數(shù)為Ki=? =0.0004?微分系數(shù)KdF=3.00上述參數(shù)在圖

6、1的仿真模型中進(jìn)行設(shè)置，仿真并記錄記錄單位階躍響應(yīng)曲線如圖2-2所示:0200400600900100012001400160018002000圖2-2單位階躍響應(yīng)曲線2結(jié)合曲線可得，在該組參數(shù)下系統(tǒng)的性能指標(biāo)如下:超調(diào)量Mp=9.2%上升時(shí)間tr=147.2s調(diào)整時(shí)間ts=262.4s穩(wěn)態(tài)誤差e (oo) =0。2.4 模型失配時(shí)PID控制器的控制效果仿真由于在實(shí)際生產(chǎn)的過(guò)程控制中，用模型表示被控對(duì)象時(shí)往往存在一定誤差，且參數(shù)也不可能是固定不變的。在已確定控制器最優(yōu) PID 調(diào)節(jié)參數(shù)下，仿真驗(yàn)證對(duì)象模型的三個(gè)參數(shù)(K T、中某一個(gè)參 數(shù)變化(不超過(guò)原值的士 5%時(shí)，系統(tǒng)出現(xiàn)模型失配時(shí)控制效果

7、改變的現(xiàn)象并分析原因。為驗(yàn)證在2.3中確定的控制器的控制效果，本設(shè)計(jì)分別對(duì)K T、p進(jìn)行增加5咐口減小5痢次變動(dòng)(共6次仿真)。然后將獲得的控制 曲線上的點(diǎn)保存到“ Workspace"，調(diào)用“plot ()”函數(shù)，對(duì)比不同 對(duì)象參數(shù)下同一控制器的控制效果，仿真記錄如下。2.4.1 被控對(duì)象比例系數(shù)改變時(shí)的控制效果如下圖2-3所示是K變化士 5%寸系統(tǒng)的單位階躍響應(yīng)曲線，其性能指標(biāo)記錄如表2-1所示：表2-1 K變化時(shí)的性能指標(biāo)超調(diào)量Mp上升時(shí)間 tr調(diào)整時(shí)間 ts穩(wěn)態(tài)誤差 e(°°)K減少5%8.4%192.6s268.2s0K增加5%14.1%177s256.

8、2s0圖2-3 K變化時(shí)的響應(yīng)曲線I11初始曲線K顏小5%K增加5% 'rjL-'1-=1iLi10100200300400500600700800900 WOO與初始被控對(duì)象相比，比例系數(shù)K的增加使系統(tǒng)超調(diào)量增加，調(diào)節(jié)時(shí)間減??；而K的減小使超調(diào)量減小，調(diào)節(jié)時(shí)間增加；但相應(yīng)的變化量都不是很大。這是因?yàn)椋豢貙?duì)象比例系數(shù)的增加（或減?。┦沟孟到y(tǒng)的開(kāi)環(huán)增益變大（變?。?，故而系統(tǒng)響應(yīng)的快速性得到提高（減 弱），但超調(diào)量也隨之增大（減?。Ｓ纱丝芍?，當(dāng)對(duì)象的比例系數(shù)在一定的范圍內(nèi)變化時(shí)，所設(shè)計(jì)的 PID控制器仍有較好的控制效果！2.4.2 被控對(duì)象時(shí)間常數(shù)改變時(shí)的控制效果如下圖2-4所

9、示是T變化士 5%寸系統(tǒng)的單位階躍響應(yīng)曲線，其性能指標(biāo)記錄如表2-2所示：表2-2 T變化時(shí)的性能指標(biāo)超調(diào)量Mp上升時(shí)間 tr調(diào)整時(shí)間 ts穩(wěn)態(tài)誤差 e（°°）T減少5%13.7%174.8s251.4s0T增加5%9.1%193.2s275.6s0圖2-4 T變化時(shí)的響應(yīng)曲線4 J«inii)ii|:初始曲線1!1!1!T摑小 £%4 T-1 -Z；iT 增加 5%1. f；：；n r_qj!13!31 A C n a UNA1:U fni i i i i i i i iu01002003004C0500 GOO 700 SOO 903 1000與初始

10、被控對(duì)象相比，對(duì)象時(shí)間常數(shù)T的增加使系統(tǒng)超調(diào)量減小, 調(diào)節(jié)時(shí)間增大；而T的減小使超調(diào)量增加，調(diào)節(jié)時(shí)間減??；但相應(yīng)的 變化量都不是很大，整體性能指標(biāo)仍較為優(yōu)良。由此可知，當(dāng)被控對(duì)象的時(shí)間常數(shù)T在一定的范圍內(nèi)變化時(shí)，所設(shè)計(jì)的PID控制器的控制效果仍較為優(yōu)良！2.4.3 純滯后時(shí)間改變時(shí)的控制效果如下圖2-5所示是p變化士 5%寸系統(tǒng)的單位階躍響應(yīng)曲線，其性能指標(biāo)記錄如表2-3所示：表2-3 r變化時(shí)的性能指標(biāo)超調(diào)量Mp上升時(shí)間 tr調(diào)整時(shí)間 ts穩(wěn)態(tài)誤差 e(°°)T減少5%7.9%182.2s251s0T增加5%14.6%186.6s271.8s0圖2-5 r變化時(shí)的響應(yīng)曲線

11、與初始被控對(duì)象相比，對(duì)象純滯后時(shí)間T的增加使系統(tǒng)超調(diào)量增 加，調(diào)節(jié)時(shí)間增加；而P的減小使超調(diào)量減小，調(diào)節(jié)時(shí)間減??；但相 應(yīng)的變化量都不是很大，整體性能指標(biāo)較為優(yōu)良。由此可知，當(dāng)被控對(duì)象的時(shí)間常數(shù)T在一定的范圍內(nèi)變化時(shí)，所設(shè)計(jì)的PID控制器的應(yīng)仍能滿足控制要求。2.5 執(zhí)行機(jī)構(gòu)非線性對(duì)PID控制器控制效果的影響實(shí)際的控制系統(tǒng)中往往會(huì)存在非線性環(huán)節(jié)，如執(zhí)行機(jī)構(gòu)的非線性等。系統(tǒng)的非線性會(huì)對(duì)控制器的控制效果造成一定的影響?？紤]到這一點(diǎn)，本設(shè)計(jì)在原控制系統(tǒng)仿真框圖控制器的輸出后加入了飽和非線 性環(huán)節(jié)，以此來(lái)研究非線性對(duì)原控制器控制效果的影響。仿真框圖如圖2-6所示。圖2-6加入了非線性環(huán)節(jié)的仿真框圖下圖

12、2-7是加入了飽和非線性后控制系統(tǒng)的響應(yīng)曲線， 與原系統(tǒng) 響應(yīng)相比，超調(diào)量、調(diào)節(jié)時(shí)間等都有所增加。因?yàn)榭刂破鬏敵鲞_(dá)到飽 和值時(shí)，控制量的增加對(duì)控制效果的影響不再增加， 因而達(dá)不到期望 的控制效果。理論上非線性環(huán)節(jié)的加入會(huì)使控制效果變差，仿真結(jié)果與此是一致的！圖2-7飽和非線性環(huán)節(jié)響應(yīng)曲線與原響應(yīng)對(duì)比1 dIT1 210.80 6040 20 (初始曲線加入非荻性環(huán)節(jié)后的響應(yīng)曲線kKal-1-1，-UIFr1fi)100200300400500600700800900 WOO2.6 擾動(dòng)作用對(duì)PID控制器控制效果的影響在本小節(jié)中主要分析給系統(tǒng)施加小的擾動(dòng)信號(hào)時(shí)，控制器的控制 效果。這里選擇階躍信

13、號(hào)為擾動(dòng)信號(hào)，分別在反饋回路和控制器輸出 兩個(gè)位置施加這一擾動(dòng)作用。相應(yīng)的響應(yīng)曲線如下所示初始曲甘1-反憬回番加階法擾市/飛，J% ; - 11 »! V- : !' ' BmaBirniamaBnnf105 DOlogo 1EQQ 2QQQ 25002.6.1 反饋回路中的階 躍擾動(dòng)對(duì)控制效果的影 響如下圖2-8所示，是 在反饋回路中施加階躍 信號(hào)時(shí)的響應(yīng)曲線。圖2-8反饋回路加擾動(dòng)由此可知，系統(tǒng)仍然穩(wěn)定。但由于在反饋回路上 t=500s時(shí)加了 一個(gè)幅度為0.1的階躍擾動(dòng)，相當(dāng)于在t=500s時(shí)把系統(tǒng)的輸入值減 小了 0.1。這樣系統(tǒng)穩(wěn)定時(shí)穩(wěn)態(tài)值就變?yōu)榱?0.9 ,

14、亦即在反饋回路上 的擾動(dòng)沒(méi)有被克服。2.6.2 控制器輸出位置的擾動(dòng)信號(hào)對(duì)控制效果的影響下圖2-9所示是在控制器輸出位置加階躍擾動(dòng)信號(hào)時(shí)的響應(yīng)曲線。系統(tǒng)仍然按照原先的給定值進(jìn)行調(diào)節(jié)，最終穩(wěn)定在 1.0處；只是 在施加擾動(dòng)以后系統(tǒng)重新調(diào)平，有一個(gè)更大的超調(diào)，調(diào)節(jié)時(shí)間也略有 增加。圖2-9控制器輸出位置加擾動(dòng)初始曲線控制器愉出位況加擾動(dòng)金一- -+_-/ /L I,_111050010001500 20CO 2500三、總結(jié)本次課程設(shè)計(jì)通過(guò)設(shè)計(jì) PID控制器，并對(duì)控制系統(tǒng)在改變被控對(duì)象參數(shù)、施加干擾等情況下進(jìn)行 MATLA的真，研究分析了所設(shè)計(jì)控制器的控制效果。將以前所學(xué)的控制理論知識(shí)與仿真實(shí)驗(yàn)相

15、結(jié)合，鞏固了原有知識(shí)，也提高了對(duì) Simulink工具箱的運(yùn)用能力。當(dāng)然，課設(shè)中也遇到了一些問(wèn)題。其一是理論計(jì)算的控制參數(shù)控 制性能較差；其二是將多組曲線進(jìn)行合并比較。針對(duì)第一個(gè)問(wèn)題，本 設(shè)計(jì)主要以理論值為參考，單獨(dú)改變某個(gè)參數(shù)以獲得相應(yīng)響應(yīng)曲線； 針對(duì)新的響應(yīng)曲線上所顯現(xiàn)的問(wèn)題， 進(jìn)一步調(diào)整其它參數(shù)，最終獲得 較為滿意的控制參數(shù)。而對(duì)于第二個(gè)問(wèn)題，作者最終將響應(yīng)曲線中的 數(shù)據(jù)以結(jié)構(gòu)體的方式保存到“工作空間”中，將其視為數(shù)組，調(diào)用“plot()”函數(shù)完成相應(yīng)的觀察任務(wù)。在附錄中也添加了大部分程序段， 它們其實(shí)僅改變了調(diào)用的數(shù)據(jù)，程序其實(shí)是類(lèi)似的。上述兩個(gè)問(wèn)題的解決可以說(shuō)是本次課設(shè)的一大收獲，由

16、此可見(jiàn)， 解決問(wèn)題應(yīng)該總攬全局，活學(xué)活用。就像調(diào)整控制參數(shù)一樣，不能按 理論值“一根筋”地鉆進(jìn)去，換個(gè)角度結(jié)合試湊法，問(wèn)題就能迎刃而 解。最后我要感謝辛勤教育我們的每一個(gè)老師， 正是老師們的付出教 會(huì)了我們相關(guān)的知識(shí)和技能，讓我們能順利完成任務(wù)！四、參考文獻(xiàn)1 .王書(shū)鋒，譚建豪.計(jì)算機(jī)控制技術(shù)M.武漢：華中科技大學(xué)出版社2 .楊佳，許強(qiáng)，徐鵬，余成波.控制系統(tǒng)MATLAB©真與設(shè)計(jì)M.北京：清華大 學(xué)出版社3 .薛定宇.反饋控制系統(tǒng)設(shè)計(jì)與分析-Matlab 語(yǔ)言應(yīng)用M.北京：清華大學(xué)出 版社，2000.五、附錄在合并不同參數(shù)下的圖像時(shí)，對(duì)保存在工作空間中的數(shù)據(jù)進(jìn)行了 操作。本附錄是部

17、分程序。可對(duì)應(yīng)相應(yīng)注釋進(jìn)行查看，此處不再具體 劃分對(duì)照。clc;%改變對(duì)象的比例系數(shù)Kplot(h1.time,h1.signals.values,'k');grid on;hold on;%原始參數(shù)下的響應(yīng)曲線plot(h2.time,h2.signals.values,'r-');hold on;%K減小5%后的響應(yīng)曲線plot(h3.time,h3.signals.values,'b-.');hold on;%K增加5%后的響應(yīng)曲線legend(初始曲線','K減小5%','K增加5%');clc;%

18、改變對(duì)象的時(shí)間常數(shù)Tplot(hT1.time,hT1.signals.values,'k');grid on;hold on;%原始參數(shù)下的響應(yīng)曲線plot(hT2.time,hT2.signals.values,'r-');hold on;%T減小5%后的響應(yīng)曲線plot(hT3.time,hT3.signals.values,'b-.');hold on;%T增加5%后的響應(yīng)曲線legend(初始曲線,'T減小5%','T增加5%');clc;%改變對(duì)象的純滯后時(shí)間tplot(ht1.time,ht1.sign