文思奇 林榮文 黃李威

（福州大學(xué)電氣工程與自動(dòng)化學(xué)院，福州 350108）

在交流伺服系統(tǒng)中，多數(shù)采用三相永磁同步電動(dòng)機(jī)（PMSM）。這類電動(dòng)機(jī)與無刷直流電動(dòng)機(jī)（BLDC）和感應(yīng)異步電動(dòng)機(jī)（IM）相比，其效率高，運(yùn)行可靠。永磁同步電動(dòng)機(jī)有兩種類型，一種是內(nèi)裝式，另一種是外裝式。外裝式常見的隱極式同步電動(dòng)機(jī)，其特點(diǎn)是在d軸和q軸上，磁路的磁阻基本相同。內(nèi)裝式常見的凸極式同步電動(dòng)機(jī)[1]，其特點(diǎn)是在q軸上的磁阻小于d軸，且結(jié)構(gòu)簡單，制造成本低，易對永磁磁極進(jìn)行優(yōu)化設(shè)計(jì)[2]。

目前在工業(yè)上大量使用常規(guī)的PID參數(shù)調(diào)節(jié)來實(shí)現(xiàn)控制伺服系統(tǒng)，這種固定參數(shù)的調(diào)節(jié)方式由于系統(tǒng)的非線性和被控對象特性變化時(shí)，控制性能比較差。于是可以借助智能控制算法來達(dá)到更好的控制效果。

本文主要介紹基于遺傳算法的伺服系統(tǒng)PID參數(shù)整定，通過Matlab仿真Kp、Ki、Kd參數(shù)自整定得到最優(yōu)解。

1 數(shù)控機(jī)床進(jìn)給伺服系統(tǒng)的數(shù)學(xué)模型

1.1 PMSM的數(shù)學(xué)模型

通過坐標(biāo)變換將三相坐標(biāo)系下的定子交流電流變換為兩相旋轉(zhuǎn)坐標(biāo)系下的直軸電流和交軸電流[3]。

為了得到理想情況：①假設(shè)磁路不飽和；②定子三相繞組對稱，忽略空間高次諧波；③忽略鐵心損耗（渦流損耗和磁滯損耗）；④定子空載電動(dòng)勢是比較標(biāo)準(zhǔn)的正弦波；⑥系統(tǒng)中頻率和溫度對繞組的電阻假設(shè)無影響，一直保持恒定值。圖1為PMSM等效結(jié)構(gòu)坐標(biāo)圖。

圖1 PMSM等效結(jié)構(gòu)坐標(biāo)圖

圖1 中，Oa、Ob、Oc為三相定子繞組的軸線，取轉(zhuǎn)子的軸線與定子a相繞組的軸線的電氣角為θ，可以建立永磁同步電動(dòng)機(jī)的物理方程[4]為

式（1）、式（2）分別是感應(yīng)電動(dòng)機(jī)的電壓方程和磁鏈方程，ua、ub、uc是三相定子繞組的電壓，ia、ib、ic是三相定子繞組的電流，?a、?b、?c是三相定子繞組的磁鏈，定子電阻Ra＝Rb＝Rc＝R，定子電感La＝Lb＝Lc=L；?f為轉(zhuǎn)子磁場的等效磁鏈。

三相定子交流電的主要作用是產(chǎn)生一個(gè)旋轉(zhuǎn)的磁場[5]，當(dāng)兩相相位正交的對稱繞組通以兩相相位相差90°的交流電時(shí)，也能產(chǎn)生旋轉(zhuǎn)磁場。根據(jù)矢量變換原理，建立一個(gè)與永磁同步電動(dòng)機(jī)轉(zhuǎn)子同步旋轉(zhuǎn)的d-q坐標(biāo)系，讓d軸與轉(zhuǎn)子磁極重合，順著旋轉(zhuǎn)方向超前90°電角度為q軸，把永磁同步電動(dòng)機(jī)定子的各參量都轉(zhuǎn)化到d-q旋轉(zhuǎn)坐標(biāo)系下。

1.2 矢量變換原理

矢量變換分為Clarke和Park兩種，變換的基本的原則是變換前后磁動(dòng)勢完全等效。矢量變換時(shí)，先進(jìn)行 Clarke變換，式（3）中將三相靜止坐標(biāo)系的電流值iA、iB、iC轉(zhuǎn)換為兩相靜止坐標(biāo)系中的電流值 iα、iβ。然后再經(jīng)過 Park變換，兩相靜止坐標(biāo)系就變成了兩相旋轉(zhuǎn)坐標(biāo)系[6]。式（5）中將兩相旋轉(zhuǎn)坐標(biāo)系中的電流值 iα、iβ轉(zhuǎn)換為兩相靜止坐標(biāo)系中的電流值id、iq。

圖2為OABC到Oab的轉(zhuǎn)換。

圖2 OABC到Oαβ的轉(zhuǎn)換

圖3 Oαβ到OMT的轉(zhuǎn)換

可以得到Clarke正變換為

Clarke和Park變換其用意就是將交流電動(dòng)機(jī)電流解耦成勵(lì)磁電流跟電驅(qū)電流來控制，模擬直流電動(dòng)機(jī)一樣的控制特性。

1.3 同步旋轉(zhuǎn)坐標(biāo)上的PMSM數(shù)學(xué)模型

凸極式電動(dòng)機(jī)的自感互感都是余弦函數(shù)，將電壓方程變成含有時(shí)變系數(shù)的微分方程。圖4為永磁同步電動(dòng)機(jī)矢量圖。

圖4 永磁同步電動(dòng)機(jī)矢量圖

設(shè)定Lmd、Lmq為d、q軸的勵(lì)磁電感；Lsσ為d、q軸線圈的漏感；ud、uq為d、q軸定子電壓；id、iq為 d、q軸定子電流；φd、φq為 d、q軸定子磁鏈；Lq=Lmq+Lsσ、Ld=Lmd+Lsσ為 d、q 軸的總電感；φf是轉(zhuǎn)子磁場的等效磁鏈。

得到PMSM電流方程為

2 遺傳算法在PID控制中的應(yīng)用

遺傳算法（genetic algorithms, GA）是以自然選擇和基因遺傳理論為基礎(chǔ)[7]，模擬大自然生物進(jìn)化過程，通過選擇，交叉，變異等操作[8]，在解空間中通過指定的迭代次數(shù)進(jìn)行尋優(yōu)，向著全局最優(yōu)的方向收斂。

在遺傳算法中，根據(jù)實(shí)際需求設(shè)定一個(gè)適應(yīng)度函數(shù)，篩選出適應(yīng)度較高個(gè)體，以提高遺傳算法的搜索效率和全局搜索能力[9]。遺傳算法過程用到了適應(yīng)度準(zhǔn)則，適應(yīng)度高的會(huì)擁有較高概率遺傳到下一代中，適應(yīng)度低的則會(huì)擁有較高的概率被淘汰[10]。

選用誤差絕對值時(shí)間積分性能指標(biāo)作為參數(shù)選擇的最小目標(biāo)函數(shù)，即

式中，e(t)為系統(tǒng)誤差；u(t)為控制器輸出；tu為上升時(shí)間；w1、w2、w3、w4為權(quán)值，而且 w4＞＞w1。

本文取適應(yīng)度函數(shù)為目標(biāo)函數(shù)J的倒數(shù)，即

式中，遺傳算法在進(jìn)化搜索中基本不利用外部信息，僅以適應(yīng)度函數(shù)為依據(jù)，利用種群每個(gè)個(gè)體的適應(yīng)度來進(jìn)行搜索。因?yàn)檫m應(yīng)度函數(shù)的復(fù)雜度是遺傳算法復(fù)雜度的主要組成部分，所以適應(yīng)度函數(shù)的設(shè)計(jì)應(yīng)盡可能簡單，使計(jì)算的時(shí)間復(fù)雜度最小。

在伺服系統(tǒng)中應(yīng)用遺傳算法時(shí)，首先需要對參數(shù)進(jìn)行編碼和譯碼操作。編碼一般有二進(jìn)制編碼和浮點(diǎn)數(shù)編碼兩種方式，這里選用二進(jìn)制編碼。由于二進(jìn)制編碼不能直接在控制器算式中使用，因此必須解碼后使用。編碼時(shí)，將3個(gè)參數(shù)分別代表10位二進(jìn)制碼，同樣在解碼過程中，將這個(gè)數(shù)字串分為三部分，每十位為一個(gè)二進(jìn)制編碼。然后轉(zhuǎn)換成十進(jìn)制數(shù)。例如在工業(yè)過程控制中，Kp的變化范圍是從0.1到100。Kp譯碼公式如下：

式中，m為十位無符號(hào)二進(jìn)制編碼。

再確定一個(gè)解碼方案，相當(dāng)于基因型到表現(xiàn)型的轉(zhuǎn)換。然后進(jìn)行解碼操作。再通過已經(jīng)確定的轉(zhuǎn)換規(guī)則可得到個(gè)體的適應(yīng)度。適應(yīng)度值非負(fù)，適應(yīng)度值增加的方向就是目標(biāo)函數(shù)優(yōu)化的方向。圖5為基于遺傳算法的PID優(yōu)化控制方案原理圖。

圖5 基于遺傳算法的PID優(yōu)化控制方案原理圖

本文選用比例選擇算子。因?yàn)楸壤x擇算子嚴(yán)格的遵循個(gè)體被選中的概率和適應(yīng)度呈線性的正比關(guān)系。

式中，Pi代表被選中的概率。在這個(gè)群體中適應(yīng)度高的個(gè)體要比適應(yīng)度低的個(gè)體更容易遺傳下去，因?yàn)楸贿x中進(jìn)行復(fù)制的概率大。這樣就使得這種遺傳關(guān)系變得特別容易從適應(yīng)度中看出來，并且比較標(biāo)準(zhǔn)。

根據(jù)遺傳算法的流程，可以得到基于遺傳算法的PID參數(shù)尋優(yōu)的流程圖。與單純的GA算法類似，也是一個(gè)迭代循環(huán)。

圖6是基于遺傳算法的PID整定流程圖，可以看出在經(jīng)過主要的選擇變異交叉操作后，進(jìn)行判斷，若群體不穩(wěn)定，則繼續(xù)迭代；若達(dá)到穩(wěn)定條件，則停止操作，輸出尋優(yōu)結(jié)果。

這種算法方法簡單，僅僅需要對數(shù)字串進(jìn)行復(fù)制交換或者改變某一位。由于復(fù)制、交叉和變異功能以及群點(diǎn)尋優(yōu)的方式避免了陷入局部最優(yōu)解問題，所以這種方法不受模型初始條件的影響，可以在完全不知道系統(tǒng)初始狀態(tài)的情況下，在一個(gè)解空間中進(jìn)行迭代尋優(yōu)。一般來說，先根據(jù)傳統(tǒng)工業(yè)控制中的實(shí)際操作經(jīng)驗(yàn)，選取一組參數(shù)，然后進(jìn)行遺傳算法的迭代，直到達(dá)到控制要求為止。這樣就可以節(jié)省很多時(shí)間以及避免不必要的時(shí)間浪費(fèi)以及計(jì)算浪費(fèi)，提高了控制效率。

把遺傳算法運(yùn)用于PID的參數(shù)整定，不但克服了常規(guī)PID整定方法的缺點(diǎn)，而且實(shí)現(xiàn)了PID參數(shù)離線整定。這種方法把遺傳學(xué)的思想與控制工程相結(jié)合，模擬自然進(jìn)化過程，搜尋最優(yōu)解。

圖6 遺傳算法（GA）流程

3 進(jìn)行Matlab仿真

圖7 階躍響應(yīng)輸出曲線

圖8 階躍響應(yīng)誤差曲線

圖9 最優(yōu)個(gè)體適應(yīng)值

圖 10 Kp、Ki、Kd優(yōu)化曲線

使用經(jīng)驗(yàn)法選定Kp、Ki、Kd參數(shù)時(shí)，假設(shè)選定的Kp=0.5、Ki=60、Kd=0，當(dāng)輸入為一階躍函數(shù)時(shí)，得到輸出曲線和誤差曲線，分別如圖7和圖8所示。

圖7中調(diào)節(jié)時(shí)間為0.071s，很明顯能看出此次調(diào)節(jié)時(shí)間較長，且系統(tǒng)波動(dòng)較大。若應(yīng)用在實(shí)際的工業(yè)控制中，控制效果不太好。使用遺傳算法時(shí)，進(jìn)行Matlab仿真。

其中，Kp=0.1656、Ki=65.8216、Kd=0，適應(yīng)度函數(shù)的最優(yōu)值是Jbest=11.64。可以看出在100次迭代中，最優(yōu)個(gè)體適應(yīng)度值在 20次循環(huán)后逐漸趨于穩(wěn)定。在圖10中，采樣點(diǎn)取每次迭代的最后一次結(jié)果，也就是最優(yōu)解，采樣100次后，繪制出優(yōu)化曲線圖?？梢钥闯龅玫終p、Ki、Kd都已經(jīng)穩(wěn)定，取出這個(gè)穩(wěn)定值就是最優(yōu)解。Kp、Ki、Kd得到最優(yōu)解后繪制輸出曲線和誤差曲線，分別如圖11和圖12所示。

圖11 GA優(yōu)化階躍響應(yīng)輸出曲線

圖12 GA優(yōu)化階躍響應(yīng)誤差曲線

調(diào)節(jié)時(shí)間為0.0481s，相比于使用經(jīng)驗(yàn)法，調(diào)節(jié)時(shí)間已經(jīng)明顯變短，圖形波動(dòng)也沒有經(jīng)驗(yàn)法變化的劇烈，波形更加平滑。因此，可以知道基于遺傳算法的PID參數(shù)調(diào)節(jié)比僅僅依靠經(jīng)驗(yàn)值進(jìn)行估算Kp、Ki、Kd的方式收斂速度更快。

4 結(jié)論

通過本文的研究，傳統(tǒng)的PID控制依靠經(jīng)驗(yàn)選取參數(shù)，曲線波動(dòng)較大，調(diào)節(jié)時(shí)間也比較長。相比于傳統(tǒng)的PID控制，基于遺傳算法的PID控制效果更好，實(shí)現(xiàn)了PID參數(shù)調(diào)節(jié)的全局最優(yōu)化。經(jīng)過仿真實(shí)驗(yàn)，可以得出基于遺傳算法的PID參數(shù)整定與傳統(tǒng)的PID參數(shù)選定相對比，具有更加精確的優(yōu)點(diǎn)，多目標(biāo)尋優(yōu)，運(yùn)算速度更快，魯棒性強(qiáng)，是理想快捷的優(yōu)化算法。

[1] 陳堅(jiān). 交流電機(jī)數(shù)學(xué)模型及調(diào)速系統(tǒng)[M]. 北京: 國防工業(yè)出版社, 1989.

[2] 陳修亮, 車倍凱. 永磁同步電機(jī)矢量控制解耦方法的研究[J]. 電氣技術(shù), 2013, 14(4): 37-40, 43.

[3] 廖效果. 數(shù)控技術(shù)[M]. 武漢: 湖北科學(xué)技術(shù)出版社,2000.

[4] 邱麗, 曾貴娥, 朱學(xué)峰. 幾種PID控制器參數(shù)整定方法的比較研究[J]. 自動(dòng)化技術(shù)與應(yīng)用, 2005, 24(11):28-31.

[5] 夏加寬, 王桐. 一種新穎的不對稱六相 SVPWM 算法[J]. 電氣技術(shù), 2014, 15(2): 9-13, 19.

[6] 胡金高, 程國揚(yáng). 魯棒近似時(shí)間最優(yōu)控制及其在電機(jī)伺服系統(tǒng)的應(yīng)用[J]. 電工技術(shù)學(xué)報(bào), 2014, 29(7):163-172.

[7] 左月飛, 張捷, 劉闖, 等. 基于自抗擾控制的永磁同步電機(jī)位置伺服系統(tǒng)一體化設(shè)計(jì)[J]. 電工技術(shù)學(xué)報(bào),2016, 31(11): 51-58.

[8] 肖曦, 許青松, 王雅婷, 等. 基于遺傳算法的內(nèi)埋式永磁同步電機(jī)參數(shù)辨識(shí)方法[J]. 電工技術(shù)學(xué)報(bào),2014, 29(3): 21-26.

[9] 張曉楠, 方浩, 戴冠. 遺傳算法的編碼機(jī)制研究[J].信息與控制, 1997, 26(2): 134-139.

[10] 郭慶鼎, 李蒙, 郭威. PID控制器參數(shù)的遺傳算法優(yōu)化設(shè)計(jì)[J]. 沈陽工業(yè)大學(xué)學(xué)報(bào), 2000, 22(1): 31-33.