基于改進雙態粒子群算法的BLDCM分數階控制器

金 鵬

(遼寧工程職業學院 電氣工程系，鐵嶺 112008)

0 引 言

無刷直流電機(以下簡稱BLDCM)具有體積小、質量輕、出力大、響應快、速度高、慣量小、轉動平滑、力矩穩定的優點[1]，并且容易實現智能化控制，現在越來越多地被作為精密控制電機使用，如機械手臂關節電機等[1]。但BLDCM是一個典型的多變量、非線性、強耦合、時變的復雜系統[1]，難以建立精確模型。普通PID控制必須建立精確的數學模型，并且參數調節難度大，控制精度不高[1]。而一些智能PID算法往往編程復雜，并且需要大量的參數調節經驗，在實際的工程應用中可行度不高。

分數階PIλDμ(以下簡稱FOPID)控制器是近年來提出的一種新型控制器，具有kp、ki、kd、λ、μ5個調節參數，因此，控制精度和靈活度更高[1]。目前，有學者提出通過模糊算法或復雜智能算法調節5個參數[2-4]，但這些調節方法需要大量的參數經驗，并且迭代計算復雜。

本文結合混沌算法和雙態粒子群算法(以下簡稱BPSO)，形成CBPSO算法，利用CBPSO在線調節FOPID控制器的參數。仿真實驗及電機轉速測試結果表明：基于CBPSO-FOPID算法，BLDCM伺服系統獲得更好的動態性能和魯棒性能。

1 BLDCM數學模型

在忽略電機損耗及電樞反應對氣隙磁場影響的情況下，將整個BLDCM看作一個整體，則BLDCM的電壓平衡方程[5]：

(1)

式中：u為電樞電壓；i為電樞電流；r0為相電阻,r0=2r；L0為相電感，L0=2(L-M)；ke為反電動勢系數；ω為電機轉速[5]。

轉矩平衡方程[5]：

(2)

Te=kti

(3)

式中：Te為電磁轉矩；TL為負載轉矩；J為轉動慣量；B為粘滯摩擦系數；kt為轉矩系數[5]。

2 FOPID控制器設計

BLDCM控制系統采用雙閉環控制系統。內環為電流環，采用PI控制器；外環為轉速環，是本文研究重點，采用FOPID控制器。

2.1 FOPID控制器整數化

Podlubony在分數階微積分的基礎上提出了FOPID微分方程[1-6]：

u(t)=Kpe(t)+KiDi-λe(t)+KdDdμe(t)

(4)

式中：λ>0，μ>0。

(5)

得出FOPID的傳遞函數：

G(s)=Kp+Kis-λ+Kdsμ

(6)

λ和μ可以取大于0的任意實數[7]，因此，FOPID具更高的控制精度、靈活度及無限維特性[1-6]，但不能直接應用于過程控制中?，F廣泛采用的方法是利用Oustaloup濾波算法在有限頻域段(wb,wh)內[7]，通過頻域擬合的方法將sα近似成整數階傳遞函數的形式[7]。但此方法在高頻區域近似效果不好[2]，這里采用改進Oustaloup濾波算法[5]：

(7)

2N+1為濾波器階次，wh及wb分別為(wb,wh)的上、下限，且b=10,d=9。

為證明FOPID控制器近似為PID控制器的近似度，設FOPID控制器傳遞函數：

Oustaloup算法近似整數階傳遞函數：

改進Oustaloup算法近似整數階傳遞函數：

在頻段(101,104)，圖1為G1(s)與G(s)的頻率特性對比圖，圖2為G2(s)與G(s)的頻率特性對比圖。可以看出，改進Oustaloup濾波算法在高頻段[7](103,104)的近似度要明顯高于Oustaloup濾波算法。

圖1 Oustaloup頻率特性圖2 改進Oustaloup頻率特性

2.2 FOPID控制器參數尋優

FOPID共有5個參數需要尋優求解。針對PSO在迭代計算的后期，收斂速度變慢、易于陷入局部極值的問題，本文結合混沌算法和BPSO，形成CBPSO算法。首先，通過Tent混沌映射初始化粒子群的位置，使粒子群隨機、均勻分布在搜索空間內[1-7]，然后利用雙態粒子群算法對5個參數同時進行搜索。Tent混沌映射如下：

xk+1=1-2|xk-0.5| 0≤xk≤1

(8)

(1)BPSO

捕食狀態：雙態粒子群算法將粒子群分為兩個不同行為狀態的子群——探索群體和捕食群體[8]。在粒子群迭代初期，所有粒子處于捕食狀態。捕食狀態的粒子群行為狀態與普通粒子群算法一樣。設粒子規模為M，搜索空間維數為D，種群的最大迭代次數為T，粒子群中粒子i(1≤i≤M)的位置為Xi=(xi1,…,xid)，速度為Vi=(vi1,…,vid)，它根據自己的個體最優值pid和全局最優值pgd在D維空間中搜索得出最優解，粒子i的速度和位置按式(9)和式(10)更新下一代的速度和位置[1]。

vid(t+1)=ω(t)·vid(t)+r1·c1·[pid(t)-

xid(t)]+r2·c2·[pgd(t)-xid(t)]

(9)

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

(10)

式中：t為粒子i的迭代次數；c1，c2為粒子i的學習因子；r1，r2分別為(0,1)間的隨機數。ω(t)為粒子i的慣性權重，本文采用自適應慣性權重，如下式：

(11)

式中：ωmax，ωmin分別為慣性權重的最大值和最小值；fi為粒子i的適應度，favg為當前粒子群適應度平均值；fg為當前粒子群最優適應度。本文選取BLDCM系統指標函數JITAE為粒子群適應度函數，即：

(12)

探索狀態：在粒子群尋優過程中，若算法陷入局部最優值，保留一部分優良的粒子繼續以捕食狀態搜索，其余大部分粒子轉為探索狀態，在搜索空間內重新初始化，繼續尋優，這部分粒子位置：

(13)

速度更新公式和位置更新公式：

xkd(t+1)=xkd(t)+vkd(t+1)

(14)

vkd(t+1)=ω·sign(r)·vkd+c1·r1·

[pkd(t)-xkd(t)]+c2·r2·

[pgd(t)-xkd(t)]+c3·u(0,1)·

[xkd(t)-hpgd(t)]

(15)

(16)

式中：sign(r)為符號函數；c3為(0,1)之間隨機數；u(0,1)為高斯分布函數；hpgd為搜索狀態的粒子在解集空間內的全局最優值。

在迭代計算中，根據式(17)判斷粒子群是否陷入局部最優值，若連續5次δ(t)≤0.1，即粒子群適應度基本無變化，則認為算法陷入局部最優值。

(17)

另外，從上述分析過程也可看出，BPSO是兩個階段的迭代算法，算法的編程難度不高，便于應用于實際工程領域。

(2)算法流程：

步驟1：粒子群參數初始化。

步驟2：按式(8)初始化粒子位置和速度，生成FOPID的5參數的初始位置序列。

步驟3：按式(9)、式(10)、式(12)更新粒子位置、速度和適應度。

步驟4：按式(17)計算δ(t)。若算法陷入局部極值，跳步驟5，否則執行步驟7。

步驟5：按適應度排序粒子，性能較好的30%粒子以捕食狀態行動，按式(9)、式(10)更新速度和位置，其余粒子轉為探索狀態，按式(14)、式(15)、式(12)更新位置、速度和適應度。

步驟6：分別計算捕食狀態的和探索狀態的全局最優值，比較后更新全局最優值。

步驟7：算法迭代次數達到T，尋優結束，求得Kp、Ki、Kd、λ、μ全局最優解。

3 仿真實驗

在MATLAB平臺上對本文算法進行仿真。電機參數設置如下：額定電壓UN=36 V，額定轉矩TN=0.23 N·m，定子相繞組電阻R=0.8 Ω，定子相繞組有效電感L-M=1.5 mH，轉動慣量J=0.001 05 kg·m2，極對數p=4。粒子種群參數設置如下：M=100，T=100，D=5，c1=2、c2=2，搜索空間上限G1=[50 50 50 2 2]、下限G2=[0 0 0 0 0]。

(1)CBPSO與PSO收斂性能對比

當速度環采用FOPID算法時，分別用CBPSO、BPSO搜索FOPID的5個參數，適應度收斂曲線如圖3所示。從圖3可以看出，CBPSO收斂速度更快，ITAE性能指標更優。

圖3 適應度收斂曲線

(2)系統性能對比

速度環分別采用普通PID、普通FOPID、PSO-FOPID、CBPSO-FOPID算法，設定轉速為2 000 r/min，在t=2.5 s時，加入負載TL=0.2 N·m。對比4種算法下BLDCM系統的動態性能及抗干擾性能。圖4為速度響應曲線，圖5為電磁轉矩響應曲線，圖6為相電流曲線。

圖4 電機轉速響應曲線

圖5 電機電磁轉矩響應曲線

圖6 電機相電流響應曲線

從圖4可以看出：普通PID算法下，BLDCM產生了較大的超調；普通FOPID算法下，BLDCM也產生了較小幅度的超調；PSO-FOPID和CBPSO-FOPID算法下，系統無超調，但CBPSO-FOPID的響應速度更快，在系統遇到擾動時，CBPSO-FOPID的抗干擾能力更強。

從圖5和圖6可以看出，當系統加入負載時，CBPSO-FOPID算法下，BLDCM的電磁轉矩及相電流的波動最小，再次證明CBPSO-FOPID的抗干擾能力更強。

4 實驗測試

實驗測試采用TMS320F28335電機專用DSP處理器[7]作為無刷直流伺服電機的主控制器。電機參數與仿真參數一致。實驗時，設置轉速為700 r/min。在2.5 s施加負載轉矩0.2 N·m。基于PID、FOPID、PSO-FOPID、CBPSO-FOPID算法的電機轉速如圖7所示，性能指標對比如表1所示。

表1 算法性能指標

圖7 電機轉速曲線

通過圖7及表1可以看出：PID和FOPID都產生了不同程度的超調，但PID的穩態誤差明顯較大，擾動量也比PID大；PSO-FOPID和CBPSO-FOPID無超調，調節時間短、擾動量也很?。?種算法中，CBPSO-FOPID的性能指標最好。實驗結果再次證明，本文算法具有更優越的動態性能及魯棒性能。

5 結 語

通過前文分析可以看出，改進的分數階控制器可以有效提高伺服系統的控制精度、靈活性和抗干擾性。并且通過仿真及實測結果可以得出以下結論：(1)利用CBPSO尋優FOPID控制器參數時，尋優速度快，控制精度較高，魯棒性較好。CBPSO-FOPID控制器與PID控制器、FOPID控制器、PSO-FOPID控制器相比，具有更高的性能指標。(2)CBPSO算法收斂速度快，避免了普通PSO算法陷入局部極小值的缺陷，搜索的解值更優。