董智鵬

　　(桂林電子科技大學 機電工程學院，廣西 桂林 541004)

       摘要：微位移精密定位系統中的壓電陶瓷元器件具有非光滑特性，無法直接測得壓電元器件的輸入輸出信號，常規方法難以對其進行有效的辨識和控制。文章采用三明治模型來描繪納米微位移平臺，并提出一種基于最小二乘支持向量機優化的辨識方法解決三明治遲滯模型的辨識問題。最后基于已經辨識的三明治模型，設計一個PID逆補償控制器，解決定位系統的精密軌跡控制問題。

　　關鍵詞：三明治系統；遲滯；最小二乘支持向量機；逆補償前饋控制

　　中圖分類號：TP273.1文獻標識碼：ADOI： 10.19358/j.issn.1674-7720.2017.03.017

　　引用格式：董智鵬.MPT系列微位移定位系統的辨識與控制［J］.微型機與應用，2017,36（3）：55-58.

　0引言

　　微位移定位平臺主要應用于天文射電望遠鏡。因其具有反應速度快和輸出大等特點［13］也適用于工廠自動化中的精確定位。經過柔性鉸鏈裝置將微小的位移量進一步放大，從而得到所需要的實際信號輸出。為了獲得精準、可控的微位移運動軌跡，直接加載在壓電器件上的壓電信號較大。因此，由控制器產生的弱電控制信號必須經過調制解調、放大等處理。圖1描述了PEA的工作過程。遲滯環節的兩端可以直接線性化處理；但具有逆壓電效應的壓電陶瓷是一個非光滑非線性遲滯系統。精密定位系統中的遲滯環節不僅會降低控制精度，更嚴重的還會使得控制系統發散［2 3］。所以對遲滯三明治系統的辨識一直是個熱點。

1含有遲滯三明治系統的描述和分析

　　根據文獻［4］對微位移系統的定義，三明治系統是兩端線性環節中間夾著一個非線性環節的模型結構，如圖2所示。

　　圖2中，L1(·)和L2(·)表示線性動態環節，H(·)表示非線性環節；u(k)為輸入信號,y(k)為輸出信號,ω1(k)為中間環節的輸入信號,ω2(k)是中間環節的輸出信號。本文采用ARMRAX模型來描述L1(·)和L2(·)，則微位移定位系統可表示為：

　　式中，n1a和n1b、a1i和a2i分別為L1(·)和L2(·)的階次和權重，q1和q2為純延時；非線性環節由函數H(·)描述。

2基于LSSVM的三明治系統辨識

　　2.1最小二乘支持向量機

　　與標準的ε支持向量機相比，LSSVM在樣本不是很大的情況下具有較好的擬合性。其基本思想也是通過在高維空間構造決策函數y(x)再引進內積變換和相應的核函數。

　　y(x)=sgn［ω(x)+b］（4）

　　對于最小二乘支持向量機，尋優問題為：

　　為求解上述尋優問題，需要構造拉格朗日等式，把約束決策問題變為無約束超平面問題：

　　式中,αi(i=1,2,…，n)是Lagrange權值。

　　根據KuhnTucker條件,得:

　　

　　可以發現最小二乘支持向量機的訓練實際上是通過奇異值分解法求解一個線性方程組。

　　2.2線性環節辨識

　　在三明治系統中，遲滯環節H(·)不僅具有非光滑、次環、多值映射等復雜特性，而且無法測得其信號，現有的辨識方法很難奏效。從圖3（a）中可知，遲滯具有復雜的特性。首先設計一個退化激勵信號［5］,將遲滯環節限制在單調上升的模態之中，即：

　　ω2(k)=f［ω1(k)］（10）

　　退化后的單模態遲滯環節如圖3（b）所示。

　　這樣就可以將線性環節辨識出來，在具體編寫程序時采用增量斜坡信號代替原始信號。根據文獻［6］的關鍵項分離原則，于是式（1）~式（3）所描述的系統構建了一個具有參數線性化結構的整體模型。

　　這樣就可以采用廣義梯度算法［7］來辨識兩端線性環節。

　　2.3動態遲滯環節的辨識

　　重新設計一個正弦衰減信號形式如下：

　　式中：φ是相位，γ是幅值偏置；Am和fm是信號的峰值和截止頻率，α和β是幅值和信號的退化系數。根據2.2節已經辨識的線性環節將遲滯環節的輸入輸出信號重構出來，即

　　由文獻［8 9］可得“基于最小二乘支持向量機”所構造的微位移模型可以很好地描述遲滯的復雜特性，尤其是參數具有較好的冗余性，可以降低計算機的故障率。從而使得PEA具有更高的運動精度，因此本文采用LSSVM代表遲滯環節。

3基于三明治模型的軌跡控制器設計

　　根據已經辨識的遲滯三明治模型及其逆模型，設計一個帶有逆補償控制器的系統，如圖4所示。

　　圖4中，r(k)、e(k)、u(k)、y(k)分別為閉環控制信號；ω(k)和v(k)為對象及其逆模型的中間信號。控制器由兩部分組成：逆補償控制器和前向控制器，uf(k)和up(k)分別表示系統中兩個子控制器的輸出信號，且有：

　　u(k)=uf(k)+up(k)（15）

　　在所設計的控制器中，逆補償控制器由三明治系統逆模型構成，用于抵消PEA中的遲滯特性。另一方面，當逆模型存在建模誤差或者系統受到外界干擾產生誤差綜上，從性能指標可以看出，本文算法對含有不同程度高斯噪聲或者椒鹽噪聲的巖石顆粒圖像的去噪表現良好，優于基于Contourlet的BayesShrink去噪效果。

4實驗結果

　　采用江蘇匯博公司生產的MPT系列壓電陶瓷執行器來驗證所提出的辨識和控制方案。額定輸入電壓和輸出位移為10 V和10 μm；控制器由研華公司生產的PCL818LS型組成，使用Turbo C編寫控制程序，采用控制頻率為33 kHz。

　　4.1辨識結果

　　(1)線性環節辨識信號采用M序列，幅值mf=3 mV，斜率mA=10>2 mf。這樣就可以使遲滯特性限制在單調上升模態內，具體在相關程序中可以采用3階多項式逼近。接下來采用關鍵項分離原理［6］來辨識兩端線性環節，RGIA的初始值為：P(0)=1.0×106×I(6)，γ(0)=0.001，θ(0)=0.001×J(1,6)。e(k－1)2+e(k)2<ε是判斷線性環節收斂的依據，式中ε=10－8，I是單位陣，J是單位行向量。經過608步迭代后收斂，圖5是線性模型參數的收斂過程。

　　(2)重新設計一個正弦衰減信號：

　　(3)采用最小二乘支持向量機來辨識遲滯環節的神經網絡模型。核函數選取高斯徑向基函數，經過32次迭代后神經網絡收斂，如圖6所示。

整個系統辨識完成后，還需要設計一個正弦調制信號來驗證建模的精度。

　　uv(k)=(3.0×e－10×k)×［sin(2π×(300×e－20×k))+1.0］(17)

　　并且在檢驗信號中隨機加入一個幅值為10 V、寬度為0.1 ms的強干擾噪聲。對應的模型檢驗結果如圖7所示。實驗證明，該建模方法具有良好的動態性能。

　　4.2軌跡控制結果

　　實驗中，控制頻率fs=33 kHz，由于過高的電壓可能會破壞電路，所以需要對輸入信號u(k)進行一定限制，具體限制在［-10,10］之間。依據式(13）、（14)，可以首先將補償環節上的逆模型－11(·)和－12(·)計算出來，再由相應的遲滯環節辨識結果，利用每次給出的反饋來不斷訓練對應的遲滯逆模型。最終的前向控制器C(·)被設計為：C(·):Kp=1.120 4;Ki=0.382 0;Kd=0.152 0;被控對象的輸入輸出如圖8所示。

　　為了更好地說明本建模方法的優越性，本文還采用了普通PID控制器進行軌跡的預測控制，在編程過程中進行了一系列的微調，最終的PID控制器參變量為：Kp=1.756，Ki=0.627和Kd=4.234，軌跡跟蹤對比如圖9所示。

5結論

　　經過一系列的實驗論證，所提出的辨識方案與實際軌跡對比具有很小的誤差；在此基礎上設計的逆補償控制器不僅精度很高，而且動態效果很好，完全滿足了微位移控制的要求。

　　參考文獻

　　［1］ GE P. Modeling and control of hysteresis in piezoceramic actuator ［D］. Kingston: University of Rhode Island, 1996.

　　［2］ HU H.Compensation of hysteresis in piezoceramic actuators and control of nanopositioning system ［D］. Toronto: University of Toronto, 2003.［3］ MAYERGOYZ I D. Mathematical models of hysteresis［M］. New York: SpringerVerlag, 1991.

　　［4］ TAWARE A, TAO G. Control of sandwich nonlinear systems ［M］.New York: Springer, 2003.

　　［5］ 謝揚球,譚永紅.壓電陶瓷執行器的非光滑三明治模型辨識與內模控制［J］.控制理論與應用，2013,30（5）：567-576.

　　［6］ TAN Y H,DONG R L,LI R Y.Recursive identification of sandwich system with deadzone and application［J］.IEEE Transactions on Control System Technology,2009,17(4):945-951.

　　［7］ 方崇智,蕭德云.過程辨識［M］.北京：清華大學出版社,2004.

　　［8］ 閻威武,邵惠鶴.支持向量機和最小二乘支持向量機的比較及應用研究［J］.控制與決策，2003,18（3）：358-360.

　　［9］ 桂衛華，宋海鷹，楊春華.HammersteinWiener模型最小二乘向量機辨識及其應用［J］.控制理論與應用，2008,25（3）：393-397.