基于樽海鞘算法的PID參數(shù)優(yōu)化

劉亞飛，郝玉然，韓超杰

（1.華北水利水電大學 電力學院，河南 鄭州 450000；2.鄭州地鐵集團有限公司，河南 鄭州 450000）

0 引 言

PID控制由于結(jié)構(gòu)簡單、有效性高、控制效果較好、魯棒性強以及可靠性高，廣泛應用在工業(yè)領域。有很多方法可以調(diào)整PID控制器參數(shù)，如手動調(diào)整Ziegler-Nichols法、科恩-庫恩調(diào)整法以及Z-N階躍響應法，但手動調(diào)整存在一定的局限性。隨著智能算法的不斷提出，智能算法優(yōu)化PID參數(shù)取得了較好的結(jié)果。遺傳算法（Genetic Algorithm，GA）[1,2]、粒子群算法（Particle Swarm Optimization，PSO）[3-5]、蟻群算法（Ant Colony Optimization，ACO）[6]、人群搜索算法（Seeker optimization Algorithm，SOA）[7-9]、螢火蟲算法（Firefly Algorithm，F(xiàn)A）[10]以及果蠅優(yōu)化算法（Fruit Fly Optimization Algorithm，F(xiàn)OA）[11,12]等，在優(yōu)化PID控制器參數(shù)中使系統(tǒng)取得了較好的穩(wěn)態(tài)特性和性能指標。

2017年，Seyedali等人通過模擬樽海鞘在海洋航行和覓食時的群集行為，提出了樽海鞘群算法（Salp Swarm Algorithm，SSA）[13]。在多個數(shù)學優(yōu)化函數(shù)上對算法進行測試，觀察并驗證了它們在尋找優(yōu)化問題最優(yōu)解時的有效行為。

鑒于樽海鞘群算法簡單易行、收斂快以及計算量小等特點，本文提出將該算法優(yōu)化不穩(wěn)定系統(tǒng)PID控制器的參數(shù)結(jié)果與粒子群算法（PSO）優(yōu)化結(jié)果進行比較。

1 PID控制器

控制系統(tǒng)中有PID控制、滑模控制以及魯棒控制等控制方法。PID控制是工業(yè)領域最常見的控制，是一種線性控制器。它根據(jù)輸入值r(t)與輸出值y(t)的誤差e(t)控制系統(tǒng)，其中誤差e(t)為：

PID的控制規(guī)律u(t)為：

式中，Kp、Ki以及Kd為系統(tǒng)誤差的比例系數(shù)、積分系數(shù)以及微分系數(shù)，τ為積分變量。

2 樽海鞘算法

樽海鞘是一種生活在海洋深處、身體透明、形狀類似水桶的海洋無脊椎動物。它的運動方式與身體構(gòu)成都與水母十分相似，主要依靠吸收身體周圍的水產(chǎn)生動力。不同于大多數(shù)群居生物，樽海鞘由于其自身結(jié)構(gòu)，分布形式并非以“群”而是以“鏈”的形式首尾相連，后一個追隨前一個的形式移動。

這條鏈分為領導者和追隨者兩組。領導者在鏈的開始處占據(jù)位置，而其余的樽海鞘是追隨者。在SSA的數(shù)學模型中，領導者的位置更新使用方式為：

式中，xi1表示領導者在第i維中的位置；Fi為食物源在第i維中的位置；ubi、lbi為第i維的上、下邊界；c1、c2、c3為隨機數(shù)。

c1是平衡探索與開發(fā)的重要參數(shù)，數(shù)學定義為：

式中，t為當前迭代；T為最大迭代次數(shù)；c2和c3是生成的[0，1]內(nèi)的隨機數(shù)。更新追隨者位置的方式為：

式中，j≥2；xij表示第j個追隨者在第i維中的位置。

運算迭代過程一般情況下都需要設置兩個門限值充當?shù)K止條件。兩個門限值中，一個為最大迭代次數(shù)，另一個是適應度門限。算法對目前的最優(yōu)適應度同提前設置好的門限適應度作比較，當其大于門限適應度時終止循環(huán)，輸出當前最優(yōu)適應度值對應的坐標值用于充當目標位置的估計值。當最優(yōu)適應度始終都未達到其門限值時，看迭代次數(shù)是否達到設置好的最大迭代次數(shù)，如果達到則輸出整個迭代過程中適應度最優(yōu)位置的點坐標作為目標位置的估計值。

3 仿真測試

3.1 實驗平臺

實驗搭建的運行環(huán)境為Windows 7、3.20 GHz處理器、8 GB內(nèi)存，所有算法代碼均用Matlab 2016b編程實現(xiàn)。

3.2 目標函數(shù)的選取和控制對象

本文選取ITAE指標作為算法的目標函數(shù)。ITAE的定義為：

為了驗證提出的算法IMFO的可行性，選擇不穩(wěn)定系統(tǒng)進行控制，并分別與使用PSO、MFO與該系統(tǒng)PID控制器參數(shù)整定結(jié)果對比。

不穩(wěn)定系統(tǒng)拉普拉斯變換（Laplace Transform）為：

式中，s為復頻率。Simulink下的系統(tǒng)模型如圖1所示。

3.3 仿真結(jié)果

利用PSO、SSA算法對該系統(tǒng)進行仿真，得PSO、SSA優(yōu)化適應度函數(shù)曲線和系統(tǒng)階躍響應輸出曲線分別如圖2和圖3所示。獨立運行50次PSO、SSA算法，得出控制系統(tǒng)優(yōu)化結(jié)果，并記錄50次得到的適應度函數(shù)值，分別得到適應度函數(shù)的平均值、標準差、中值、最好值以及最差值，如表1所示。

圖1 PID控制仿真模型

圖2 適應度函數(shù)控制曲線

圖3 算法優(yōu)化后的系統(tǒng)階躍響應

表1 算法評價指標

由圖2適應度函數(shù)控制曲線可知，SSA算法比PSO算法收斂速度快、精度高，易找到全局最優(yōu)解。從階躍響應輸出曲線來判斷，PSO算法平穩(wěn)時間較長，超調(diào)量較大。SSA算法相對于PSO算法，算法控制精度高、收斂快，證明了本文算法的有效性和優(yōu)越性。通過表1中50次運行結(jié)果可以看出，SSA比PSO具有更好的優(yōu)化結(jié)果，算法魯棒性較好。

4 結(jié) 論

為了進一步優(yōu)化PID參數(shù)，本文提出樽海鞘算法（SSA）對PID控制器進行參數(shù)優(yōu)化，通過不穩(wěn)定系統(tǒng)仿真驗證了該算法的有效性和優(yōu)越性。通過與PSO算法對比可知，SSA算法具有更高的控制精度和魯棒性，能夠為系統(tǒng)PID參數(shù)優(yōu)化提供可靠依據(jù)，可使工業(yè)控制領域中系統(tǒng)收斂到最優(yōu)狀態(tài)，具有廣闊的應用前景。