孫漢卿，張世杰，王桂芝，連衛(wèi)民

（1.河南牧業(yè)經(jīng)濟學(xué)院 信息工程學(xué)院，鄭州 450001；2.河南工業(yè)大學(xué) 電氣工程學(xué)院，鄭州 450001）

永磁同步電機（permanentmagnet synchronous motor，PMSM）［1］常用的控制系統(tǒng)有轉(zhuǎn)子磁場定向控制（FOC）［2］和直接轉(zhuǎn)矩控制（DTC）［3］。這些系統(tǒng)都是根據(jù)速度控制器輸出速度信號來調(diào)節(jié)磁通和轉(zhuǎn)矩，故速度控制是整個系統(tǒng)的基礎(chǔ)。目前，常用的非線性速度控制方法為比例積分（proportion integration，PI）控制器［4］。然而，傳統(tǒng)PI控制器的控制參數(shù)需要人工設(shè)定，且不會隨著環(huán)境變化而自適應(yīng)改變，不能滿足對環(huán)境的魯棒性要求［5］。為此，有學(xué)者將一些自適應(yīng)機制融入PI控制器中，例如模糊PI參數(shù)自整定控制器［6］，基于遺傳算法［7］、粒子群優(yōu)化［8］等進化算法的PI參數(shù)自整定控制器等。

另外，對時變、非線性控制系統(tǒng)的控制，傳統(tǒng)PI控制器性能不是很好。增強PI控制器性能的較好方法是使用分數(shù)階控制器［9］。分數(shù)階PI控制器PIa中，通常I的階數(shù)為分數(shù)階［10］。因此，除設(shè)置KP、KI外，還需設(shè)置積分的分數(shù)階a。

基于上述分析，本文利用分數(shù)階PI控制器代替PMSM中的PI速度控制器。通過融入粒子群優(yōu)化算法中的信息共享機制和自適應(yīng)步長來改進細菌覓食優(yōu)化（bacteria foraging optimization，BFO）算法，形成改進型細菌覓食優(yōu)化（improved BFO，IBFO）算法，并利用IBFO算法來優(yōu)化KP、KI和a參數(shù)，以獲得最佳的控制器。另外，在適應(yīng)度函數(shù)中考慮了轉(zhuǎn)速控制的時間絕對誤差積分（integrated time absolute error，ITAE）和超調(diào)量，有利于提高控制穩(wěn)定性。通過實驗與傳統(tǒng)PI控制器和傳統(tǒng)BFO算法優(yōu)化參數(shù)的PI控制器進行比較，結(jié)果證明了所提出控制方法具有優(yōu)越性能。

1 PMSM 數(shù)學(xué)及控制模型

1.1 PM SM 數(shù)學(xué)模型

PMSM結(jié)構(gòu)中分為2個軸，分別為d軸和q軸，這2個軸的氣隙長度和電感都不同。d軸定義為勵磁繞組的磁通產(chǎn)生方向［11］，q軸與d軸成90°。PMSM結(jié)構(gòu)如圖1描述。

圖1 PMSM的結(jié)構(gòu)

PMSM模型可以通過轉(zhuǎn)子旋轉(zhuǎn)坐標系（d-q軸）建模為1個電路，如圖2所示。根據(jù)電路電流和電壓關(guān)系式，d軸和q軸的定子電壓分別為：

式中：vd、id表示d軸的電壓和電流，iq、vq表示q軸電壓和電流；λd和λq分別表示2個軸的磁鏈，可通過定子電流、電感和恒定磁鏈λm來計算獲得，如式（3）（4）所示；ωe表示電機轉(zhuǎn)速，rs表示定子電阻。

式中：Ld、Lq分別為d-q軸電感。將式（3）（4）代入式（1）（2）中，定子電壓方程可修改為：

電磁轉(zhuǎn)矩Te為：

式中：Tl為負載轉(zhuǎn)矩；P表示電機極數(shù)；B為阻尼系數(shù)；J為慣性矩。

圖2 PMSM的等效電路

1.2 PMSM 無傳感器控制模型

本文中，選擇磁場定向控制（FOC）方式來控制電極，其系統(tǒng)框圖如圖3所示，其主要由5個模塊組成，分別為速度控制器、電流控制器、逆變器、PWM發(fā)生器和位置估計器［12］。

圖3 PMSM的磁場定向控制系統(tǒng)框圖

1.2.1 FOC控制模型

磁場定向控制的主要思想是將定子電流矢量加在q軸上，把轉(zhuǎn)子磁通加在d軸上，這意味著id（ref）＝0。為此，電磁轉(zhuǎn)矩可修改為

式中：Kt為電機轉(zhuǎn)矩系數(shù)。那么，在狀態(tài)空間中，PMSM方程可建模如下：

1.2.2 基于EKF的速度估計

在FOC無位置傳感器的控制系統(tǒng)中，電機的當前速度是通過一些轉(zhuǎn)速估計算法，根據(jù)電機電流、電壓等信號進行估計，然后將其作為實際轉(zhuǎn)速反饋到PI控制器中，如圖3所示。為此，速度估計的準確性直接影響電機速度控制器的性能。本文將電機三相電流信息變換到αβ軸（兩相靜止坐標系），然后采用基于擴展卡爾曼濾波器（EKF）來估計電機轉(zhuǎn)速。

卡爾曼濾波器（KF）是一種遞歸算法，其使用一種多輸入-輸出系統(tǒng)的隨機狀態(tài)空間模型和噪聲統(tǒng)計來獲得系統(tǒng)狀態(tài)變量的最優(yōu)估計。本文假設(shè)系統(tǒng)中的過程噪聲和測量噪聲均為零均值的高斯白噪聲。通過最小化估計誤差的協(xié)方差來獲得最優(yōu)的狀態(tài)估計。EKF可以實現(xiàn)非線性系統(tǒng)的近似線性化，提高求解精度。PMSM 模型的離散化可表示如下：

式中：xe（k）為k時刻所估計的狀態(tài)向量，即xe＝（isα，isβ，λrα，λrβ，ωr，TL）T；z（k）為k時刻的觀測向量，即為狀態(tài)向量與輸入之間的非線性函數(shù)；ue為輸入向量；He為觀測矩陣；G為噪聲加權(quán)矩陣；v（t）為測量噪聲。

如上所述，EKF將PMSM的離散化非線性模型進行線性化，并在線性化中對狀態(tài)變量進行估計，表示為

2 分數(shù)階PI速度控制器

在PMSM的磁場定向控制系統(tǒng)中，速度控制器是整個系統(tǒng)的核心部分，其決定著后續(xù)轉(zhuǎn)矩和勵磁的控制性能。通常速度控制器采用傳統(tǒng)PI控制，但其性能有限，且參數(shù)不能自適應(yīng)調(diào)整。近年來，越來越多的研究人員在控制領(lǐng)域使用分數(shù)階微積分控制，分數(shù)階微積分是處理分數(shù)冪微分或積分算子的一種技術(shù)［12］。若設(shè)定ωref為電機的參考速度，ωact為驅(qū)動器的實際輸出速度，那么在時域中的分數(shù)階PI控制器PIa的表現(xiàn)形式為

控制器的輸入e（t）為速度差，表示為

式中：為分數(shù)微分算子，表示函數(shù)f（t）的a階導(dǎo)數(shù)，定義如下：

分數(shù)階PI控制器的傳遞函數(shù)在頻域上表示為

式中：KP和Ki為分數(shù)階PI控制器的比例和積分增益值；a為分數(shù)積分器的非整數(shù)階。

對于分數(shù)階PI控制器的設(shè)計，參數(shù)KP、Ki和a決定了其上升時間（tr）、穩(wěn)定時間（ts）、峰值超調(diào)百分比（po%）和穩(wěn)態(tài)誤差（ess）。因此，要獲得良好的瞬態(tài)響應(yīng)，需要獲得合適的控制器參數(shù)組合。由于這些參數(shù)組合的解空間較大，人工調(diào)節(jié)效率較低，所以通常采用智能優(yōu)化算法來尋找最優(yōu)參數(shù)組合。

3 基于IBFO的速度控制器參數(shù)優(yōu)化

3.1 參數(shù)優(yōu)化系統(tǒng)架構(gòu)

本文將分數(shù)階PI速度控制器中的3個參數(shù)KP、Ki和a作為優(yōu)化對象。提出一種改進型BFO算法（IBFO），基于速度差信息來優(yōu)化控制器參數(shù)，以此來提高速度控制的準確性?；贗BFO的速度控制器的系統(tǒng)框架如圖4所示。

圖4 PI速度控制器參數(shù)優(yōu)化框架

通常，PI控制器參數(shù)優(yōu)化的目標函數(shù)為轉(zhuǎn)速控制的時間絕對誤差積分（ITAE）［13］，表示為

為綜合考慮速度超調(diào)情況，本文在目標函數(shù)中融入了超調(diào)量σ，σ由速度峰值與參考速度的比值來計算。這樣有利于PI控制器產(chǎn)生更小的超調(diào)量和振蕩。那么，最終的目標函數(shù)O的表達式為

式中：w為權(quán)重因子，用來權(quán)衡ITAE與超調(diào)量的比重。這里，由于ITAE反映速度控制的穩(wěn)定性，相對而言比超調(diào)量更為重要，所以本文設(shè)置w為0.3。另外，目標函數(shù)值越大說明該參數(shù)組合下控制器的性能越低。

3.2 傳統(tǒng)細菌覓食優(yōu)化算法（BFO）

BFO是由人類大腸中“大腸桿菌”覓食行為而啟發(fā)的一種新型智能優(yōu)化方法。在這種覓食理論中，細菌以單位時間內(nèi)攝入能量最小化的方式來搜索并獲得營養(yǎng)物質(zhì)，細菌會通過趨化性行為聚集到營養(yǎng)豐富的區(qū)域。BFO算法中細菌覓食步驟主要分為3個部分，分別為趨向性操作、繁殖操作和遷徙操作［14］。

趨向性操作：趨向性操作中包含2個基本操作，即旋轉(zhuǎn)和游動［15］，細菌在覓食過程中通過這些操作逐步移動到食物源處。游動操作是細菌根據(jù)上一步的方向進行游動。而旋轉(zhuǎn)操作中，細菌會根據(jù)當前方向上位置和上一步位置的適應(yīng)值旋轉(zhuǎn)方向，朝其他隨機方向游動。

繁殖操作：細菌繁殖過程中，根據(jù)細菌的適應(yīng)度，淘汰一些尋優(yōu)能力較弱的細菌，并復(fù)制相應(yīng)數(shù)量的較優(yōu)細菌進行替換。繁殖過程中產(chǎn)生的較優(yōu)細菌副本具有與本體相同的屬性。

遷徙操作：菌群經(jīng)過若干代繁殖后，將以給定概率執(zhí)行遷徙操作，移動到新的尋優(yōu)空間，以此使細菌進行全局尋優(yōu)，避免陷入局部最優(yōu)。

通常，細菌i在經(jīng)過第j輪趨向性行為、第k輪繁殖行為和第l輪遷徙行為后的位置用θi（j，k，l）來表示。

3.3 提出的改進型BFO算法（IBFO）

在傳統(tǒng)BFO算法中隨機生成旋轉(zhuǎn)操作的方向，隨機方向會導(dǎo)致收斂到全局最優(yōu)解的速度變慢。在IBFO算法中，利用了PSO具有良好的社會信息交換能力，使旋轉(zhuǎn)行為中的方向根據(jù)全局最佳位置和每個細菌的個體最佳位置來決定，而不是隨機確定。在趨化循環(huán)中，游動方向的更新公式為

式中：φ（j）表示細菌在第j次趨向性操作中的方向；pbest為每個細菌的個體最佳位置；gbest為菌群的全局最佳位置；C1、C2為學(xué)習因子，通常設(shè)置為2；wφ為慣性權(quán)重，決定對上一個方向的繼承程度。IBFO算法的流程和具體步驟見圖5所示。

步驟1初始化參數(shù)：p、S、Nc、Ns、Nre、Ned、Ped、C（i）（i＝1，2，3，…，S）、θi。其中：p為搜索空間的維度；S為細菌總數(shù)；Nc為趨化操作的代數(shù)；Ns為最大游動代數(shù)，Nc＞Ns；Nre為繁殖操作的代數(shù)；Ned為遷徙操作的代數(shù)；Ped為遷徙的概率；θi為第i個細菌的位置；C（i）為旋轉(zhuǎn)方向的步長大小。初始化種群，生成隨機向量φ（j），其取值范圍為［-1，1］。

圖5 提出的IBFO算法的流程

步驟2遷徙循環(huán)：l＝l＋1。

步驟3繁殖循環(huán)：k＝k＋1。

步驟4趨向性循環(huán)：j＝j(luò)＋1。

對于每個細菌i（i＝1，2，3，…，S）：

1）計算細菌i的適應(yīng)度函數(shù)O（i，j，k，l）；

2）令當前最佳適應(yīng)度Olast＝O（i，j，k，l）；

3）旋轉(zhuǎn)操作：計算旋轉(zhuǎn)方向

4）根據(jù)以下策略移動到新位置：

5）重新計算O（i，j，k，l）；

6）游動操作：

①設(shè)置游動代數(shù)m＝0，計數(shù)游動代數(shù)；

②當m＜Ns（即沒有達到最大代數(shù)）

· 令m＝m＋1。

· 若適應(yīng)度O（i，j，k，l）＜Olast（即獲得更好的適應(yīng)度）。

· 令Olast＝O（i，j，k，l），并移動到該位置θi（j＋1，k，l）＝θi（j，k，l）＋C（i）φ（j），計算新位置θi（j＋1，k，l）的適應(yīng)度函數(shù)值。

· 否則，令m＝Ns，循環(huán)結(jié)束。

7）如果i＜S，轉(zhuǎn)到下一個細菌（i＋1）（即返回步驟4的1）去處理下一個細菌）。

步驟5如果j＜Nc，轉(zhuǎn)步驟4，繼續(xù)趨化循環(huán)，因為細菌的生命沒有結(jié)束。

步驟6繁殖操作：

1）對于給定的k、l，計算每個細菌i經(jīng)過趨化循環(huán)后的適應(yīng)度值用來表示該細菌i的健康度量。健康度量衡量了其整個壽命中使用了多少營養(yǎng)，以及成功避免有毒物質(zhì)的能力。按照健康度量升序排列細菌Ohealth（更高的適應(yīng)度值意味著更低的健康值）。

2）淘汰具有較低健康值的細菌，并保留具有較高健康值的細菌進行分裂，繁殖產(chǎn)生的副本放置在與其父本相同的位置。

步驟7如果k＜Nre，轉(zhuǎn)至步驟3。由于在這種情況下未達到指定的繁殖代數(shù)，故啟動下一輪趨化循環(huán)。

步驟8對于細菌i，以遷徙概率Ped來移動每個細菌，保持群體中細菌數(shù)量不變。即：若要淘汰一個細菌，則需將另一個細菌分散到該優(yōu)化域上的隨機位置即可。如果l＜Ned，則轉(zhuǎn)到步驟2；否則，結(jié)束算法。

3.4 自適應(yīng)步長

BFO中，趨化性是一種覓食策略，細菌試圖爬到營養(yǎng)濃度高的地方，并避免有害物質(zhì)。在傳統(tǒng)BFO中，趨化循環(huán)中的步長是固定不變的，這使其前期尋優(yōu)速度慢且后期收斂精度低［16］。為此，本文提出了一種自適應(yīng)步長，即趨化步長隨著當前適應(yīng)度值的改變而變化，使其提供更好的收斂能力。對于第i個細菌，其步長C（i）表達式為

式中：ψ為大于0的常數(shù)，ji（θ）為第i個細菌的適應(yīng)度函數(shù)值。如果ji（θ）趨于0，則C（i）→0；如果ji（θ）趨于變大，則C（i）→1。這意味著處于有毒物質(zhì)附近的細菌具有較高的適應(yīng)度值。因此，需要賦予其較大的步長來快速遷移到營養(yǎng)濃度較高的地方。

4 仿真及分析

4.1 仿真設(shè)置

在Simulink環(huán)境里構(gòu)建電機模型和控制系統(tǒng)仿真，表1為構(gòu)建的PMSM模型參數(shù)。本文中，PI速度控制器的基本設(shè)計要求為峰值超調(diào)po%＜3%，上升時間tr＜0.02 s，穩(wěn)定時間ts＜0.8 s，穩(wěn)態(tài)誤差ess＜0.2。

表1 PMSM參數(shù)設(shè)置

對于分數(shù)階PI控制器參數(shù)，設(shè)置KP的變化范圍為0～1，KI的變化范圍為0～10，α是變化范圍為0～1的分數(shù)。

對于BFO和IBFO算法參數(shù)，設(shè)置為：S＝10、

Nc＝5、Ns＝4、Nre＝2、Ned＝2、Ped＝0.25、ψ＝180、C1＝C2＝2、wφ＝0.7。每種算法進行100次迭代，以獲得最優(yōu)參數(shù)組合解。

4.2 速度響應(yīng)

首先，利用傳統(tǒng)梯度下降法、傳統(tǒng)BFO算法和提出的IBFO優(yōu)化算法，根據(jù)電機模型對分數(shù)階PI速度控制系統(tǒng)進行參數(shù)優(yōu)化，所獲得的最優(yōu)參數(shù)值如表2所示。為了公平比較，速度控制器均采用分數(shù)階PI控制器。

表2 各種優(yōu)化方法所獲得的最優(yōu)分數(shù)階PI控制器參數(shù)

然后，將3種參數(shù)優(yōu)化方法所獲得的最優(yōu)參數(shù)配置到分數(shù)階PI速度控制器中。為了驗證這些參數(shù)的有效性，分別進行電機啟動和負載變化仿真，觀察控制器的性能。實驗中，設(shè)置電機啟動時的參考速度為1 000 r／min，并在t＝1.5 s時對電機施加1個負載，負載轉(zhuǎn)矩為TL＝10 N·m。3種PI速度控制器的電機轉(zhuǎn)速響應(yīng)曲線如圖6所示。

圖6 3種參數(shù)優(yōu)化方法下的PI速度控制器速度響應(yīng)曲線

分析圖5發(fā)現(xiàn)：基于傳統(tǒng)梯度下降法的分數(shù)階PI控制器的控制性能最差，在電機啟動時存在較大的超調(diào)量，且穩(wěn)定到參考速度的時間較長。另外，當電機轉(zhuǎn)矩發(fā)生變化時，該方法也不能較好地應(yīng)對。對于經(jīng)過傳統(tǒng)BFO參數(shù)優(yōu)化后的分數(shù)階PI控制器，其在穩(wěn)定性和超調(diào)量方面要明顯優(yōu)于梯度下降法，但其穩(wěn)定速度所需的時間依然較長。而通過本文提出的IBFO優(yōu)化后的分數(shù)階PI控制器，其在啟動時獲得穩(wěn)定速度所需的時間最短，且在負載變化時，對速度的影響最小，具有最好的魯棒性。

為了進一步比較3種方法優(yōu)化分數(shù)階PI速度控制器的性能，進行3次同樣的實驗，分別統(tǒng)計速度控制的平均上升時間、峰值超調(diào)、穩(wěn)定時間和穩(wěn)態(tài)誤差，結(jié)果如表3所示。

通過具體的性能參數(shù)可以看出：提出的IBFO優(yōu)化分數(shù)階PI控制器在各項性能方面都獲得了最佳表現(xiàn)，特別是在超調(diào)量和穩(wěn)態(tài)誤差方面。這是因為，本文在IBFO中的適應(yīng)度函數(shù)中融入了ITAE和超調(diào)量因素，并以此來獲得最佳的控制器參數(shù)。

表3 各種優(yōu)化方法下的速度控制性能指標

5 結(jié)束語

本文提出一種改進型BFO算法用來優(yōu)化PMSM中分數(shù)階PI速度控制器的3個參數(shù)，以此來提高對時變、非線性速度控制系統(tǒng)的控制性能，以及電機負載轉(zhuǎn)矩變化的魯棒性。將所提方法與梯度下降法、BFO方法進行比較，結(jié)果證明了所提方法在上升時間、峰值超調(diào)、穩(wěn)定時間和穩(wěn)態(tài)誤差等方面具有較優(yōu)越的性能，也具有較高的實用價值。

本文中，速度控制器的反饋速度由傳統(tǒng)擴展卡爾曼濾波法估計，而反饋速度的準確性直接影響速度控制的性能。為了進一步提高PMSM控制系統(tǒng)的性能，未來將考慮更先進的速度估計方法，以進一步提高速度控制的準確性。