艾志偉，嵇建波，李 靜，黃書童

快速反射鏡狀態(tài)模型構(gòu)建方法及其控制系統(tǒng)設(shè)計

艾志偉，嵇建波，李 靜，黃書童

（桂林航天工業(yè)學院，廣西 桂林 541004）

為減少快速反射鏡狀態(tài)空間建模過程中所需的結(jié)構(gòu)參數(shù)數(shù)量，提出了一種基于系統(tǒng)辨識的狀態(tài)模型構(gòu)建方法，采用該方法建立狀態(tài)模型時只需使用音圈電阻和電感兩個結(jié)構(gòu)參數(shù)?；跔顟B(tài)模型，設(shè)計了一套由降階觀測器、狀態(tài)反饋、內(nèi)模和鎮(zhèn)定補償器構(gòu)成的組合控制系統(tǒng)，利用狀態(tài)反饋完成對內(nèi)模和鎮(zhèn)定補償器的設(shè)計，通過設(shè)計降階觀測器實現(xiàn)對電流和角速度的獲取，組合系統(tǒng)可同時實現(xiàn)對輸入信號的漸進跟蹤和干擾的抑制。在SIMULINK中建立仿真模型，仿真結(jié)果顯示，不考慮干擾作用時，相較于不完全微分PID（Proportion Integration Differentiation）控制系統(tǒng)，組合系統(tǒng)的調(diào)節(jié)時間下降了53.6%，超調(diào)量上升了131.2%；加入干擾信號后，不完全微分PID控制系統(tǒng)的動穩(wěn)態(tài)性能有明顯下降，而組合系統(tǒng)的輸出性能基本不受影響。仿真結(jié)果驗證了理論分析的正確性。

快速反射鏡；狀態(tài)模型；漸進跟蹤；干擾抑制

0 引言

快速反射鏡是一種精密調(diào)整光束傳播方向的裝置，具有結(jié)構(gòu)緊湊、調(diào)整精度高、調(diào)節(jié)速度快的優(yōu)勢[1]，與大慣量的機架共同使用時，可組成粗-精跟蹤系統(tǒng)，達到穩(wěn)定視軸和消除相移的目的，因而被廣泛地應(yīng)用在激光通信[2]、天文望遠鏡[3]和自適應(yīng)光學系統(tǒng)的波前校正[4]中。

目前，對快速反射鏡系統(tǒng)的研究主要集中在支承結(jié)構(gòu)的設(shè)計優(yōu)化[5]以及控制系統(tǒng)的高精度高帶寬設(shè)計、擾動抑制[6]等方面。為了提高音圈電機驅(qū)動快速反射鏡的動態(tài)響應(yīng)特性，中科院長春光機所的王帥[7]提出了一種基于電流環(huán)的控制方法，通過為每個音圈電機配置一個高帶寬的電流環(huán)，實現(xiàn)同一運動軸上的兩個音圈電機的同步推拉作用；中科院光電所的鄧超、田竟[8-9]在分析建立音圈電機驅(qū)動的快速反射鏡傳遞函數(shù)模型基礎(chǔ)上，設(shè)計了基于多種多環(huán)路控制的干擾觀測器，實現(xiàn)對擾動的抑制；北京航空航天大學的桂瀟怡[10]通過在快速反射鏡底座中設(shè)置兩組MEMS加速度傳感器測量干擾形成的控制前饋，并采用自適應(yīng)方法調(diào)節(jié)前饋環(huán)節(jié)系數(shù)，提高了振動環(huán)境下系統(tǒng)的穩(wěn)定性能。

這些方法在進行系統(tǒng)設(shè)計時，使用的被控對象數(shù)學模型都是頻域傳遞函數(shù)模型，提高控制精度或者進行干擾抑制時，均使用了額外的傳感器，多傳感器的使用必然增加硬件系統(tǒng)的復(fù)雜性。為了在不增加控制系統(tǒng)物理硬件復(fù)雜性的基礎(chǔ)上得到期望的系統(tǒng)穩(wěn)態(tài)性能和瞬態(tài)性能，同時實現(xiàn)對干擾信號的抑制，提出了一種基于系統(tǒng)辨識的狀態(tài)模型建立方法，該方法在建立狀態(tài)模型的過程中除辨識出的系統(tǒng)參數(shù)外僅用到了音圈電機的線圈電阻和電感，大大減少了建模過程中結(jié)構(gòu)參數(shù)的使用，基于狀態(tài)模型，設(shè)計了由降階觀測器、狀態(tài)反饋、內(nèi)模和鎮(zhèn)定補償器構(gòu)成的組合系統(tǒng)，組合系統(tǒng)在不增加額外傳感器的基礎(chǔ)上實現(xiàn)了多環(huán)控制的目的，同時實現(xiàn)了漸進跟蹤和干擾抑制的目標。最后，對控制系統(tǒng)進行仿真驗證，驗證所設(shè)計的控制系統(tǒng)的跟蹤性能和擾動抑制能力。

1 快速反射鏡的狀態(tài)模型

音圈電機驅(qū)動的二自由度快速反射鏡具有良好的對稱性，為了便于分析，對其中的一個運動軸進行動力學分析，轉(zhuǎn)矩平衡方程為：

式中：為鏡片的轉(zhuǎn)動慣量；為反射鏡的輸出轉(zhuǎn)角；c為音圈電機線圈的質(zhì)量；為等效的阻尼系數(shù)；為柔性支承結(jié)構(gòu)的轉(zhuǎn)動剛度；為音圈電機作用點到轉(zhuǎn)動中心的距離；是音圈電機作用在運動方向上的力矩，大小為：

＝2f(2)

式中：為通過音圈電機的電流；f是音圈電機力矩系數(shù)。

音圈電機為反射鏡運動機構(gòu)提供驅(qū)動力矩，根據(jù)基爾霍夫電壓定律有

式中：是音圈電機的工作電壓；是音圈電機線圈的電感；e是音圈電機的反電動勢系數(shù)；為音圈電機運動的位移?？焖俜瓷溏R的運動范圍較小，近似有：

≈(4)

聯(lián)立以上各式，采用古典控制理論的建模方法可知快速反射鏡系統(tǒng)數(shù)學模型為一個振蕩環(huán)節(jié)和慣性環(huán)節(jié)組成的三階線性定常系統(tǒng)，其傳遞函數(shù)可表示為：

式中：1、2、3、4是待辨識的快速反射鏡參數(shù)；0、1、2、0是常系數(shù)。

根據(jù)式(1)、式(2)、式(3)和式(4)，選取負載反射鏡的轉(zhuǎn)動角度、角速度、通過音圈電機的電流為系統(tǒng)狀態(tài)變量，可以得到快速反射鏡的狀態(tài)空間方程：

聯(lián)立式(5)、式(6)有：

建立快速反射鏡的狀態(tài)空間模型，方法之一是獲取到快速反射鏡的各個參數(shù)，這樣狀態(tài)模型自然可以隨之得到，但是由于影響快速反射鏡結(jié)構(gòu)參數(shù)的量有20個之多，其中很多參數(shù)是很難甚至是無法得到，因此采用這種方法建立快速反射鏡的狀態(tài)模型也變得很困難。

通過對式(7)進行分析可知，如果可以通過系統(tǒng)辨識得到快速反射鏡傳遞函數(shù)的參數(shù)0,1,2,0，那么確定狀態(tài)模型時只需要得到音圈電機的電阻和電感，使用這6個參量即可完成狀態(tài)空間模型的搭建。

2 控制系統(tǒng)設(shè)計

在經(jīng)典控制理論中，提高快速反射鏡系統(tǒng)的控制精度的普遍的做法是將負載反射鏡的位置、速度和加速度或者位置、速度和電流等信號同時進行檢測并組成多環(huán)控制系統(tǒng)，再為每個環(huán)路設(shè)計合適的控制器達到提高系統(tǒng)輸出精度的目的。工程上，速度、加速度以及電流這些物理量雖然是可以測量的，但是檢測裝置價格也都比較昂貴，這必然會增加控制系統(tǒng)的成本與增加系統(tǒng)的復(fù)雜性。基于快速反射鏡的狀態(tài)模型，可以有效地解決上述問題，通過設(shè)計內(nèi)模、補償器以及進行狀態(tài)反饋，同樣可以實現(xiàn)多環(huán)控制的效果，達到提高控制精度的目的。針對狀態(tài)反饋過程中難以獲取的物理量，可以設(shè)計降階觀測器，對狀態(tài)變量進行估計，然后再將其應(yīng)用到狀態(tài)反饋環(huán)節(jié)中?；跔顟B(tài)模型的組合系統(tǒng)控制框圖如圖1所示，下面分別對各環(huán)節(jié)進行參數(shù)設(shè)計。

圖1 快速反射鏡組合控制系統(tǒng)

先將快速反射鏡系統(tǒng)方程寫成如下形式：

假定系統(tǒng)能控能觀，、、分別是對應(yīng)式(6)的矩陣，()是干擾信號且滿足狀態(tài)方程：

輸入信號()由下列方程產(chǎn)生：

假定上述干擾信號和輸入信號均可測量，()和()分別是干擾信號()和輸入信號()的特征多項式，即：

令()為()和()在右半閉平面零點的最小公倍式，則其展開式可以表示如下：

根據(jù)內(nèi)模原理，內(nèi)模－1()的狀態(tài)空間模型可實現(xiàn)為：

在不考慮干擾信號時，誤差信號＝－，可得內(nèi)模狀態(tài)方程為：

快速反射鏡系統(tǒng)與內(nèi)模串聯(lián)組成的系統(tǒng)狀態(tài)方程為：

假設(shè)組合系統(tǒng)能控，引入狀態(tài)反饋控制c和鎮(zhèn)定補償器，控制輸出為：

引入狀態(tài)反饋后的組合系統(tǒng)的狀態(tài)方程為：

下面繼續(xù)對受控快速反射鏡系統(tǒng)進行研究。系統(tǒng)輸出是反射鏡的轉(zhuǎn)角，將這個狀態(tài)變量用1表示。組合系統(tǒng)需要使用的是轉(zhuǎn)角，角速度和電流3個狀態(tài)變量，因為轉(zhuǎn)角已知，要得到全部3個狀態(tài)變量，只需要構(gòu)造一個低階的觀測器估計剩余的2個狀態(tài)即可，這2個狀態(tài)變量用2表示，據(jù)此，可將受控快速反射鏡系統(tǒng)的狀態(tài)空間模型寫成如下形式：

對上式進行展開，有：

下面對不能測量得到的電流和角速度兩個狀態(tài)變量設(shè)計觀測器進行觀測，令：

則上式可以變成：

如果將視為系統(tǒng)輸入量，視為輸出量，22和12分別視作系統(tǒng)矩陣和輸出矩陣，則上式可以看成是受控快速反射鏡系統(tǒng)的一個子系統(tǒng)，2個狀態(tài)變量均不能直接量測得到，為此可以通過構(gòu)建一個2維的狀態(tài)觀測器，實現(xiàn)對這2個狀態(tài)的重構(gòu)，狀態(tài)不可量測子系統(tǒng)的狀態(tài)觀測器為：

為了在原系統(tǒng)中可以實現(xiàn)該狀態(tài)觀測器，需要消去中間變量，，將式(20)帶入式(22)，可得：

工程中的信號通常包含有大量的噪聲，在系統(tǒng)中直接對輸出信號進行微分會在控制系統(tǒng)引入噪聲干擾，為避免微分運算，引入中間變量，令：

對變量進行求導(dǎo)，經(jīng)整理可得：

降階觀測器得到系統(tǒng)狀態(tài)變量為：

根據(jù)快速反射鏡控制系統(tǒng)的性能指標要求，確定出控制系統(tǒng)閉環(huán)主導(dǎo)極點的位置，保證另一極點的實部比主導(dǎo)極點的實部大3～6倍以上時可忽略該極點對控制系統(tǒng)輸出性能的影響。降階觀測器的特征根應(yīng)該有比狀態(tài)反饋系統(tǒng)特征根更大的實部模。

由閉環(huán)系統(tǒng)主導(dǎo)極點確定的特征方程式和組合系統(tǒng)的特征方程式相等，即可確定反饋環(huán)節(jié)和鎮(zhèn)定補償器的參數(shù)。

3 仿真實驗與結(jié)果分析

對一款工程中使用的快速反射鏡進行了系統(tǒng)辨識，通過擬合得到的快速反射鏡開環(huán)模型為：

音圈電機的標稱電感和電阻分別為1.9mH和3.6W，得到快速反射鏡的狀態(tài)方程為：

采用掃頻法辨識出的快速反射鏡開環(huán)系統(tǒng)模型、擬合得到的系統(tǒng)模型以及狀態(tài)模型幅頻特性曲線對比圖如圖2所示，從圖中可以看到狀態(tài)模型與辨識出的系統(tǒng)曲線基本重合。

圖2 幅頻特性曲線對比圖

對于系統(tǒng)來說，階躍信號是最具挑戰(zhàn)性的，故選擇單位階躍信號()＝1()作為系統(tǒng)測試輸入，擾動信號選為正弦信號()＝sin()，可得()＝3＋，內(nèi)?？蓪崿F(xiàn)為：

給定快速反射鏡閉環(huán)控制系統(tǒng)的設(shè)計指標是：超調(diào)量小于10%，調(diào)節(jié)時間小于10ms。通過配置狀態(tài)反饋環(huán)節(jié)的參數(shù)將系統(tǒng)的極點配置在復(fù)平面希望位置上，結(jié)合主導(dǎo)極點的定義，可將系統(tǒng)的預(yù)期的主導(dǎo)極點取為1,2＝－440±j448，其他極點為：

3,4＝－2000±j1000，4,5＝－3000±j1500

計算并對參數(shù)修正后得狀態(tài)反饋環(huán)節(jié)和鎮(zhèn)定補償器分別為：

＝[1057.9, 0.4552, 16.9829]

c＝[2.2303, 7.3806×106, 1.2892×106]

由于采用改進根軌跡的方法設(shè)計不完全微分PID控制系統(tǒng)的參數(shù)時，也是以時域性能指標要求為設(shè)計出發(fā)點，因此選擇以同樣性能指標要求時采用該種方法確定的控制系統(tǒng)作為對照實驗系統(tǒng)。不考慮擾動信號作用時，兩種方法的系統(tǒng)輸出曲線如圖3所示。從圖中可以得到，不完全微分PID控制系統(tǒng)的峰值時間為0.0070s，調(diào)節(jié)時間為0.0049s，超調(diào)量為4.0195%，而組合系統(tǒng)的峰值時間為0.0019s，調(diào)節(jié)時間為0.0023s，超調(diào)量為9.2950%，相較于不完全微分PID控制系統(tǒng)，組合系統(tǒng)的峰值時間下降了72.85%，調(diào)節(jié)時間下降了53.06%，超調(diào)量則上升了131.2%，組合系統(tǒng)的超調(diào)量雖有所上升，但仍然滿足性能指標的要求。

圖3 不考慮干擾時兩種控制方法輸出曲線

在系統(tǒng)中施加幅值為1，圓頻率為100rad/s的正弦干擾信號前后不完全微分PID控制系統(tǒng)的輸出曲線如圖4所示。從圖中可以看到，施加擾動作用后系統(tǒng)輸出會在平衡點上下震蕩，影響了控制系統(tǒng)的穩(wěn)態(tài)精度。

加入同樣的擾動后，組合系統(tǒng)的輸出曲線如圖5所示，從圖中可以看出，加入干擾信號作用后組合系統(tǒng)的峰值時間為0.0019s，調(diào)節(jié)時間為0.0023s，超調(diào)量為9.2950%，加入干擾信號前后，組合系統(tǒng)的輸出曲線基本重合，峰值時間、調(diào)節(jié)時間和超調(diào)量等時域性能指標也完全一樣。仿真結(jié)果證明了基于狀態(tài)反饋的組合系統(tǒng)可以同時實現(xiàn)對輸入的漸進跟蹤和對干擾信號的抑制。

圖4 幅值為1頻率為100rad/s的干擾作用前后不完全微分系統(tǒng)輸出曲線對比

圖5 幅值為1頻率為100rad/s的干擾作用前后組合系統(tǒng)輸出曲線對比

加大干擾信號作用，在系統(tǒng)中施加幅值為10，圓頻率為1000rad/s的正弦干擾信號，加入干擾信號前后，不完全微分PID控制系統(tǒng)的輸出曲線和組合系統(tǒng)的輸出曲線分別如圖6和圖7所示。從圖6中可以看到，加大干擾信號作用后，不完全微分PID控制系統(tǒng)的輸出完全失真；從圖7中可以看到，加大干擾信號作用后，組合系統(tǒng)的峰值時間為0.0019s，調(diào)節(jié)時間為0.0023s，超調(diào)量為9.4326%，可以發(fā)現(xiàn)組合系統(tǒng)受干擾信號的影響很小。通過對比圖5和圖7還可以發(fā)現(xiàn)組合系統(tǒng)對擾動信號的頻率也具有較低的靈敏度，這是由于在進行組合系統(tǒng)設(shè)計過程中用內(nèi)模抵消了干擾信號的作用，實現(xiàn)了組合系統(tǒng)對擾動的全頻段抑制，與理論期待相符。

圖6 幅值為10頻率為1000rad/s的干擾作用前后不完全微分系統(tǒng)輸出曲線對比

圖7 幅值為10頻率為1000rad/s的干擾作用前后組合系統(tǒng)輸出曲線對比

4 結(jié)論

為使得設(shè)計的控制系統(tǒng)可以同時實現(xiàn)對輸入信號的漸進跟蹤和對可測干擾信號的抑制，亦不增加控制系統(tǒng)硬件復(fù)雜程度，本文從狀態(tài)模型的快速建立為出發(fā)點，提出一種基于系統(tǒng)辨識的狀態(tài)模型建立方法，該方法只需要使用6個參量即可完成狀態(tài)模型的建立，大大減少了建模過程中結(jié)構(gòu)參數(shù)的使用?；诳焖俜瓷溏R的狀態(tài)模型，設(shè)計了一套組合控制系統(tǒng)，并著重分析了組合系統(tǒng)中內(nèi)模、狀態(tài)反饋和鎮(zhèn)定補償器的設(shè)計方法，給出了組合系統(tǒng)中降階觀測器實現(xiàn)方法，分析計算出組合系統(tǒng)狀態(tài)方程。對工程上使用的一款快速反射鏡應(yīng)用了上述組合控制系統(tǒng)，給出了控制系統(tǒng)的具體參數(shù)，并進行了仿真驗證，仿真結(jié)果表明理論分析的正確性。

[1] Kluk D J, Boulet M T, Trumper D L. A high-bandwidth, high-precision, two-axis steering mirror with moving iron actuator[J]., 2012, 22: 257-270.

[2] 趙磊, 紀明, 趙振海, 等. 艦載激光武器穩(wěn)定平臺粗精復(fù)合控制[J]. 激光與紅外, 2019, 49(1): 86-92.

ZHAO Lei, JI Ming, ZHAO Zhenhai, et, al. Primary-precise compounded control for stabilized platform in shipborne laser weapon[J]., 2019, 49(1): 86-92.

[3] 呂世良, 劉金國, 周懷得, 等. 星載紅外探測器快速反射鏡控制系統(tǒng)設(shè)計[J]. 紅外與激光工程, 2017, 46(9): 102-107.

LYU Shiliang, LIU Jinguo, ZHOU Huaide, et al. Design of control system for fast steering mirror of infrared detector based on satellite[J]., 2017, 46(9): 102-107.

[4] WANG Kaidi, SU Xiuqin, LI Zhe, et al. Time-frequency characteristics optimal control of fast steering mirror for image motion compensation[J]., 2018, 47(S1): 120-126.

[5] 魯亞飛, 范大鵬, 范世珣, 等. 快速反射鏡兩軸柔性支承設(shè)計[J]. 光學精密工程, 2010, 18(12): 2574-2582.

LU Yafei, FAN Dapeng, FAN Shixun, et, al. Design of two-axis elastic support for fast steering mirror[J]., 2010, 18(12): 2574-2582.

[6] 彭樹萍, 于洪君, 王偉國, 等. 新型快速反射鏡伺服系統(tǒng)設(shè)計[J]. 紅外與激光工程, 2014, 43(5): 1610-1615.

PENG Shuping, YU Hongjun, WANG Weiguo, et al. Design of servo system for novel fast-steering mirror[J]., 2014, 43(5): 1610-1615.

[7] 王帥, 曹玉巖. 音圈電機型快速反射鏡的驅(qū)動控制系統(tǒng)設(shè)計[J]. 電子測量與儀器學報, 2017, 31(7): 1115-1121.

WANG Shuai, CAO Yuyan. Design of control system for fast steering mirror driven by voice coil actuator[J]., 2017, 31(7): 1115 - 1121.

[8] DENG Chao, TANG Tao, MAO Yao, et al. Enhanced disturbance observer based on acceleration measurement for fast steering mirror system[J]., 2017, 9(3): 1-11.

[9] TIAN Jing, YANG Wenshu, PENG Zhenming. Application of MEMS Accelerometers and Gyroscopes in Fast Steering Mirror Control Systems[J]., 2016, 440: 1-13.

[10] 桂瀟怡, 李偉鵬, 鄭新濤, 等. 基于MEMS加速度計的快速反射鏡復(fù)合控制[J]. 航天控制, 2018, 36(2): 30-36, 41.

GUI Xiaoyi, LI Weipeng, ZHENG Xintao, et al. Composite Control of Fast Steering Mirror Based on MEMS Accelerometer[J]., 2018, 36(2): 30-36, 41.

State Model Construction Method for Fast Steering Mirror and Its Control System Design

AI Zhiwei，JI Jianbo，LI Jing，HUANG Shutong

(Guilin University of Aerospace Technology, Guilin 541004, China)

To reduce structural parameters needed in the process of the state space modeling of a fast steering mirror (FSM), a modeling method based on system identification is proposed. Only two structural parameters, namely, coil resistance and inductance, are needed when using this modeling method. Based on the state model, a set of combined control systems consisting of a reduced-order observer, state feedback, an internal model, and a stabilization compensator were designed; the internal model and stabilization compensator were designed by using state feedback. Furthermore, current and angular velocity can be obtained using a reduced-order observer. The combined system can progressively track of input signals and suppress disturbances; the simulation model was built using SIMULINK. The simulation results show that, compared with the incomplete differential proportional integral derivative (PID) control system, the setting time decreases by 53.6% and the overshoot increases by 131.2% without considering the disturbance in the composite system. The dynamic and steady-state performance of the incomplete differential PID control system decreases significantly after adding the disturbance signal, while the composite system is maintained. The simulation results verify the correctness of the theoretical analysis.

fast steering mirror, state model, progressive tracking, disturbance suppression

TP273

A

1001-8891(2020)01-0040-06

2019-07-09；

2019-12-31.

艾志偉（1992-），男，碩士，助教，主要研究方向為光電跟蹤控制。E-mail：aizhiwei752@163.com。

廣西高校中青年教師科研基礎(chǔ)能力提升項目（2019KY0793）。