要亞斌，王利杰，王佳榮

(1.河北省電力勘測設(shè)計研究院，石家莊 050031；2.國網(wǎng)河北省電力公司電力科學(xué)研究院，石家莊 050021；3.青島鴻瑞電力工程咨詢有限公司，山東 青島 266100)

?

基于PSO的PID參數(shù)優(yōu)化及其工程實現(xiàn)

要亞斌1，王利杰2，王佳榮3

(1.河北省電力勘測設(shè)計研究院，石家莊 050031；2.國網(wǎng)河北省電力公司電力科學(xué)研究院，石家莊 050021；3.青島鴻瑞電力工程咨詢有限公司，山東 青島 266100)

介紹運用粒子群優(yōu)化算法(PSO)對單回路和雙回路系統(tǒng)PID控制器參數(shù)的優(yōu)化，并進行數(shù)字仿真，得出控制品質(zhì)參數(shù)。分析在被控對象的表達形式為n階慣性環(huán)節(jié)加純遲延和一階慣性環(huán)節(jié)串聯(lián)2種不同情況下，采用粒子群優(yōu)化算法尋優(yōu)的方法，說明根據(jù)給定被控對象參數(shù)計算出PID參數(shù)以及相應(yīng)的仿真曲線和控制品質(zhì)數(shù)據(jù)的工程實現(xiàn)過程。

PID控制器；參數(shù)優(yōu)化；參數(shù)整定；粒子群優(yōu)化算法

1 概述

PID控制器是工程實際中應(yīng)用最為廣泛的控制器[1]，其結(jié)構(gòu)簡單、容易實現(xiàn)、并且具有較強的魯棒性，成為工業(yè)控制的主要技術(shù)之一。自PID控制器提出以來，人們便致力于PID整定方法的研究，涌現(xiàn)出了各種各樣的方法。隨著計算機技術(shù)的發(fā)展，一些新的智能算法得到了迅速發(fā)展和廣泛應(yīng)用[2-4]，并在理論和應(yīng)用方面都有重要的意義。在這些算法中，粒子群優(yōu)化算法(PSO)是近年來發(fā)展起來的一種較為實用的進化算法[5-7]，它基于群體智能理論，通過群體中粒子間的合作與競爭產(chǎn)生的群體智能指導(dǎo)優(yōu)化搜索，具有很強的搜索能力。PSO作為一種強大的優(yōu)化工具已被廣泛應(yīng)用于優(yōu)化問題的求解[3]。以下利用PSO算法，實現(xiàn)了單回路和雙回路系統(tǒng)的PID控制器優(yōu)化。在此基礎(chǔ)上，開發(fā)出了一套PID控制器優(yōu)化程序，為PID控制系統(tǒng)的參數(shù)整定提供了一個簡單、方便的平臺。

2 PSO對PID參數(shù)優(yōu)化的方法

PSO 算法屬于進化算法的一種，和遺傳算法相似，它也是從隨機解出發(fā)，通過迭代尋找最優(yōu)解，通過適應(yīng)度來評價解的品質(zhì)，但它比遺傳算法規(guī)則更為簡單，沒有遺傳算法的“交叉”和“變異”操作，它通過追隨當(dāng)前搜索到的最優(yōu)值來尋找全局最優(yōu)。

在PSO算法中, 每個粒子代表解空間的一個候選解,粒子在搜索空間以一定的速度飛行,飛行速度根據(jù)飛行經(jīng)驗進行動態(tài)調(diào)整。該算法基于群智能的并行全局搜索策略,采用速度位搜索模型實現(xiàn)對整個空間的尋優(yōu)操作。PSO算法是模仿生物社會性行為而得出的一種全局優(yōu)化算法,是一種高效、簡單的并行搜索算法,其優(yōu)點在于概念簡單、實現(xiàn)容易、魯棒性好,并且能以較大概率收斂到全局最優(yōu),而且它對所優(yōu)化目標(biāo)的先驗知識要求甚少,一般只需知道其數(shù)值關(guān)系即可。但是該算法的慣性權(quán)重對算法性能具有很大的影響,另外,在初始群體的生成上,它是根據(jù)經(jīng)驗估計出PID 3個參數(shù)的取值范圍,并在此范圍內(nèi)采用隨機生成的方式,對其可行解空間進行搜索的,因此需要合理估計PID 3個參數(shù)的取值范圍。

假設(shè)在一個D維的目標(biāo)搜索空間中有m個粒子組成一個群體，其中第i個粒子的位置表示為向量xi=(xi1,xi2,…xiD)，(i=1,2,…，m),其速度也是一個D維的向量vi=(vi1,vi2,…viD)。第i個粒子群搜索到的最優(yōu)位置為Pg=(pg1,pg2,…,pgD)粒子群公式如下：

vid(t+1)=vid+c1r1(pid-xid(t))+c2r2(pgd-xid(t))

(1)

xid(t+1)=xid(t)+vid(t+1)

(2)

式中：i=1,2,…m，d=1,2,…,D；加速常數(shù)c1和c2為非負常數(shù)；r1和r2服從[0,1]上的均勻分布隨機數(shù)；vid(t)是第i個粒子當(dāng)前位置；Pid是第i個粒子的當(dāng)前速度；vid∈[-Vmax,Vmax],Vmax為最大限制速度，是非負數(shù)，當(dāng)vid>Vmax時，取vid=Vmax，當(dāng)vid<-Vmax時，取vid<-Vmax。

鑒于蟻群算法具有分布式并行搜索能力，且易于與其它算法結(jié)合，是一種全局優(yōu)化的算法。因此，可以利用蟻群算法的全局優(yōu)化能力優(yōu)化PID參數(shù),得到PID參數(shù)的初始值,然后根據(jù)系統(tǒng)當(dāng)前的誤差和誤差變化率,用模糊推理方法在線優(yōu)化調(diào)整PID參數(shù)的權(quán)值來動態(tài)地調(diào)整參數(shù)?；诹Ｗ尤旱腜ID參數(shù)優(yōu)化的基本思想：將PID 3個參數(shù)組合在一起設(shè)為粒子群的粒子，在搜索空間中每一個粒子以一定的速度飛行，飛行過程中每個微粒飛行速度將依據(jù)微粒的飛行經(jīng)驗和群體的飛行經(jīng)驗進行動態(tài)調(diào)整，以一定的速度向目標(biāo)逼近,最后找到最優(yōu)目標(biāo),即PID控制器的最優(yōu)參數(shù)。

3 工程實現(xiàn)

3.1 單回路熱工系統(tǒng)PID整定及其工程實現(xiàn) 單回路控制系統(tǒng)是最簡單、最基本、也是使用最廣泛的一種形式，單回路系統(tǒng)框圖見圖1。對于控制系統(tǒng)來說，在設(shè)定值發(fā)生變化或系統(tǒng)受到擾動作用后，系統(tǒng)將從原來的穩(wěn)態(tài)經(jīng)歷一個過程進入另一個新的穩(wěn)態(tài)。一個自動控制系統(tǒng)的好壞在穩(wěn)態(tài)下是難以判別，只有在過渡過程中才充分反映出來。一個良好的控制系統(tǒng)，在經(jīng)歷擾動后，一般應(yīng)平穩(wěn)、迅速和準(zhǔn)確地趨近或恢復(fù)到設(shè)定值。這就需要對調(diào)節(jié)器的控制參數(shù)進行準(zhǔn)確的設(shè)定才能達到理想的效果，而PID 調(diào)節(jié)器在出廠時一般都不會提供現(xiàn)成的PID 參數(shù)，它們只能由技術(shù)人員在生產(chǎn)過程中根據(jù)實際情況自己摸索出來。

圖1 單回路系統(tǒng)框圖

文章所提供的方法就是用最優(yōu)算法，可以通過被控對象的參數(shù)計算出一種理論上的PID調(diào)節(jié)參數(shù)，并將計算出的PID參數(shù)帶入控制系統(tǒng)進行仿真，得到這組PID參數(shù)的調(diào)節(jié)效果。使得系統(tǒng)在調(diào)試時，可以方便順利的找到最優(yōu)的參數(shù)值。

3.1.1 n階慣性環(huán)節(jié)加純遲延

在單回路控制系統(tǒng)中，當(dāng)被控對象為n階慣性環(huán)節(jié)加純遲延，即被控對象的傳遞函數(shù)表達形式為：

(3)

通過給定的K、T、n及τ的值，運用粒子群優(yōu)化算法，反算出PID控制器的優(yōu)化參數(shù)：Tp、KI、Kd以及階躍響應(yīng)曲線的超調(diào)量、調(diào)節(jié)時間和衰減率。

在熱工控制系統(tǒng)中被控對象為n階慣性環(huán)節(jié)加純遲延的時候，傳遞函數(shù)為式(3)時，這只是自平衡過程的傳遞函數(shù)。非自平衡對象應(yīng)加一個積分環(huán)節(jié)。

表1 有自平衡并帶有純遲延系統(tǒng)優(yōu)化計算結(jié)果

KpKiKd超調(diào)量/%調(diào)節(jié)時間/s衰減率/%0.495270.0067608.83452.6490.18

當(dāng)被控對象傳遞函數(shù)為：

圖2 有自平衡并帶有純遲延系統(tǒng)優(yōu)化后的單位階躍響應(yīng)曲線

圖3 有自平衡優(yōu)化后的單位階躍響應(yīng)曲線

表2 有自平衡系統(tǒng)優(yōu)化計算結(jié)果

KpKiKd超調(diào)量/%調(diào)節(jié)時間/s衰減率/%0.521680.0044407.34717.0790.92

3.1.2 一階慣性環(huán)節(jié)串聯(lián)

如果被控對象的傳遞函數(shù)形式不是n階慣性環(huán)節(jié)加純遲延的形式，而是由一階慣性環(huán)節(jié)串聯(lián)形式，對程序中做了相應(yīng)的功能，可以計算這種形式。

(4)

例如被控對象傳遞函數(shù)如為以下表示時

(5)

輸出的結(jié)果見圖4和表3。

圖4 一階慣性環(huán)節(jié)串聯(lián)系統(tǒng)PID優(yōu)化后的單位階躍響應(yīng)曲線

表3 一階慣性環(huán)節(jié)串聯(lián)系統(tǒng)優(yōu)化計算結(jié)果

KpKiKd超調(diào)量/%調(diào)節(jié)時間/s衰減率/%0.591740.0107204.38320.6497.48

當(dāng)被控對象傳遞函數(shù)表達式為：

圖5 3個環(huán)節(jié)串聯(lián)系統(tǒng)PID優(yōu)化后的單位階躍響應(yīng)曲線

表4 3個環(huán)節(jié)串聯(lián)系統(tǒng)優(yōu)化計算結(jié)果

KpKiKd超調(diào)量/%調(diào)節(jié)時間/s衰減率/%10.1176714.4327.698.86

從圖2-5及表1-4可以看出，通過PSO優(yōu)化方法，可以獲取單回路系統(tǒng)的PID優(yōu)化參數(shù)，并且取得了較好的優(yōu)化效果。證明了軟件可以有效解決單回路系統(tǒng)的PID控制器優(yōu)化問題。

3.2 雙回路(串級)熱工控制系統(tǒng)PID整定及其工程實現(xiàn)

串級系統(tǒng)是由2只調(diào)節(jié)器串聯(lián)起來工作，其中一個調(diào)節(jié)器的輸出作為另一個調(diào)節(jié)器的給定值的系統(tǒng)。前一個調(diào)節(jié)器稱為主調(diào)節(jié)器，它所檢測和控制的變量稱主變量(主被控參數(shù))，即工藝控制指標(biāo)；后一個調(diào)節(jié)器稱為副調(diào)節(jié)器，它所檢測和控制的變量稱副變量(副被控參數(shù))，是為了穩(wěn)定主變量而引入的輔助變量。

整個系統(tǒng)包括2個控制回路，主回路和副回路。副回路由副變量檢測變送、副調(diào)節(jié)器、調(diào)節(jié)閥和副過程構(gòu)成；主回路由主變量檢測變送、主調(diào)節(jié)器、副調(diào)節(jié)器、調(diào)節(jié)閥、副過程和主過程構(gòu)成。當(dāng)擾動發(fā)生時，破壞了穩(wěn)定狀態(tài)，調(diào)節(jié)器進行工作，見圖6。根據(jù)擾動施加點的位置不同，分3種情況進行分析：擾動作用于副回路；擾動作用于主過程；擾動同時作用于副回路和主過程。

圖6 雙回路系統(tǒng)示意

3.2.1 n階慣性環(huán)節(jié)加純遲延

當(dāng)串級系統(tǒng)內(nèi)外回路的被控對象也就是主副被控對象的傳遞函數(shù)表達式都為n階慣性環(huán)節(jié)加純遲延的時候。

測試的被控對象1為：

(6)

(7)

程序輸出結(jié)果見圖7和表5。

圖7 雙回路系統(tǒng)PID優(yōu)化后的單位階躍響應(yīng)曲線

表5 雙回路系統(tǒng)優(yōu)化計算結(jié)果

PID控制器參數(shù)KpKiKd超調(diào)量/%調(diào)節(jié)時間/s衰減率/%主回路0.51280.00350副回路-0.6706-0.013408.66905.690.96

測試的被控對象2為：

(8)

(9)

程序輸出結(jié)果見圖8和表6。

圖8 具有負調(diào)特性雙回路系統(tǒng)PID優(yōu)化后的單位階躍響應(yīng)曲線

表6 具有負調(diào)特性雙回路系統(tǒng)優(yōu)化計算結(jié)果

PID控制器參數(shù)KpKiKd超調(diào)量/%調(diào)節(jié)時間/s衰減率/%主回路0.55220.00530.1534副回路-1.3633-0.027808.15567.3693.61

3.2.2 一階慣性環(huán)節(jié)串聯(lián)

當(dāng)被控對象的傳遞函數(shù)不是n階慣性環(huán)節(jié)純遲延表達形式的時候，可以轉(zhuǎn)換成一階慣性環(huán)節(jié)串聯(lián)的形式，這里還用在單回路中用到的傳遞函數(shù)。

副回路傳遞函數(shù)為：

(10)

(11)

程序結(jié)果見圖9和表7。

圖9 混合形式傳遞函數(shù)雙回路系統(tǒng)PID優(yōu)化后的單位階躍響應(yīng)曲線表7 混合形式傳遞函數(shù)雙回路系統(tǒng)優(yōu)化計算結(jié)果

PID控制器參數(shù)KpKiKd超調(diào)量/%調(diào)節(jié)時間/s衰減率/%主回路0.52510.00280.5096副回路0.55460.007807.181073.8194.06

從圖7-9以及表5-7可以看出，利用該文設(shè)計的優(yōu)化軟件，可以獲取雙回路系統(tǒng)的PID優(yōu)化參數(shù)，并且取得了較好的優(yōu)化效果。證明了軟件可以有效解決雙回路系統(tǒng)的PID控制器優(yōu)化問題。

4 結(jié)束語

在簡單的熱工控制系統(tǒng)中，當(dāng)被控對象的動態(tài)特性已知，PID控制器參數(shù)就會直接影響控制系統(tǒng)的控制品質(zhì)。所以，通過對熱工控制系統(tǒng)中的被控對象動態(tài)特性參數(shù)計算出一組最優(yōu)的PID控制器參數(shù)是一項非常有實際意義的研究。針對PID單回路熱工控制系統(tǒng)和串級熱工控制系統(tǒng)，運用面向?qū)ο蟮木幊陶Z言環(huán)境，進行了這一問題的工程實現(xiàn)。在設(shè)計的程序中運用了粒子群優(yōu)化算法(PSO)計算最優(yōu)的PID控制器參數(shù)，進行了相應(yīng)的控制系統(tǒng)仿真，并求出了控制品質(zhì)參數(shù)和系統(tǒng)響應(yīng)曲線。使用這個程序進行單回路PID控制系統(tǒng)和串級PID控制系統(tǒng)控制器的參數(shù)整定會十分方便。只要被控對象的動態(tài)特性符合n階慣性環(huán)節(jié)加純遲延或一階慣性環(huán)節(jié)串聯(lián)的形式，在程序界面中輸入相應(yīng)的數(shù)據(jù)就可以得到比較理想的PID控制器參數(shù)和仿真結(jié)果。

為了更好的解決PID參數(shù)優(yōu)化這一問題，還有許多功能可以拓展，比如加入多變量的PID控制系統(tǒng)，這將需要解決多變量的解耦控制問題，也非常具有實際意義。另外，為了將程序作為一個整定PID控制器參數(shù)的工具使用，應(yīng)該能夠在運算完某一種系統(tǒng)的參數(shù)后，仍可以繼續(xù)返回主界面進行其他系統(tǒng)的選擇和運算。而且由于程序中只考慮到了被控對象傳遞函數(shù)為n階慣性環(huán)節(jié)加純遲延或一階慣性環(huán)節(jié)串聯(lián)這2種形式，仍有一定的局限性，可以在以后考慮添加運用有理分式表達的傳遞函數(shù)。

[1] 田育奇.時滯對象PID參數(shù)優(yōu)化設(shè)計方法的研究[D].北京：北京化工大學(xué)，2014.

[2] 宋莉莉，朱 超，孫萬麟.基于微粒群算法的PID控制器優(yōu)化研究[J].昌吉學(xué)院學(xué)報，2014(6):76-80.

[3] 許駿馬，侯石超，何偉光，等.遺傳算法與蟻群算法在PID優(yōu)化中的應(yīng)用比較[J].機電信息，2013(15):114-116.

[4] 劉 明，華 亮.基于PSO算法的PID控制器參數(shù)整定[J].控制工程，2016,23(1)：64-68.

[5] 朱小明，張慧斌.PSO算法的穩(wěn)定性分析及算法改進[J].計算機科學(xué)，2013，40(3)：275-278.

[6] 陳 亮，江 明，林園勝.基于PSO算法的PID控制器參數(shù)優(yōu)化研究分析[J].重慶工商大學(xué)學(xué)報(自然科學(xué)版)，2014，31(5):77-80.

[7] 曾文權(quán)，余愛民.基于改進PSO算法的任務(wù)分配研究[J].計算機工程與應(yīng)用，2013,49(13):51-55.

本文責(zé)任編輯：王洪娟

PID Parameters Optimization and Its Engineering Realization Based on PSO

Yao Yabin1,Wang Lijie2,Wang Jiarong3

(1. Hebei Electric Power Design&Reseach Institute,Shijiazhuang 050031,China; 2.State Grid Hebei Electric Power Research Institute,Shijiazhuang 050021,China;3.Qingdao Hongrui Power Engineering Consulting Co.,Ltd.,Qingdao 266100,China)

This paper introduces optimization method of PID controller parameters in single-loop and double-loop system by particle swarm optimization algorithm and obtains control quality parameters on the basis of digital simulation,analyze control object in form of series connection composed by n order intertial link with pure delay and first order inertia and use particle swarm optimization to find out PID parameters as well as the engineering process of corresponding simulation curve and control quality data.

PID controller;parameters optimization;parameter adjustment;PSO

2016-06-17

要亞斌(1979-)，男，高級工程師，主要從事電廠熱控專業(yè)設(shè)計研究工作。

TM621

B

1001-9898(2017)03-0003-04