面向GMTI任務(wù)的系繩式InSAR系統(tǒng)展開研究

張錦繡， 張志剛

哈爾濱工業(yè)大學 航天學院， 哈爾濱 150001

面向GMTI任務(wù)的系繩式InSAR系統(tǒng)展開研究

張錦繡*， 張志剛

哈爾濱工業(yè)大學 航天學院， 哈爾濱 150001

研究了短繩系式合成孔徑雷達干涉技術(shù)(InSAR)系統(tǒng)面向地面動目標檢測(GMTI)任務(wù)時的展開問題。使用系統(tǒng)的三維動力學模型，首先分析系統(tǒng)在平衡位置附近的穩(wěn)定性，確定系統(tǒng)執(zhí)行GMTI任務(wù)時面臨的問題。結(jié)合GMTI任務(wù)需求確定系統(tǒng)短系繩情況時的展開方式，使用粒子群算法對系統(tǒng)初始分離速度和方向進行優(yōu)選，使系統(tǒng)在展開結(jié)束時恰好擺動到水平方向。在展開完成后對系統(tǒng)施加阻尼控制和噴氣控制，使系統(tǒng)長期處于水平方向附近，以獲取期望順軌基線，最后通過數(shù)值仿真驗證。結(jié)果表明：考慮展開機構(gòu)靜摩擦力時，在優(yōu)選的初始狀態(tài)下無控展開，系統(tǒng)能夠到達目標位置，且在穩(wěn)定控制下能夠長期保持順軌基線大于99.6 m。

系繩式InSAR系統(tǒng)； GMTI任務(wù)； 穩(wěn)態(tài)性能分析； 粒子群算法； 穩(wěn)定控制

合成孔徑雷達干涉技術(shù)(InSAR)因其能夠全天時、全天候的探測而得到廣泛應(yīng)用[1-2]，其中一項是地面動目標檢測(GMTI)，能夠大面積偵察、檢測、跟蹤、定位地面動目標，有著重要的軍事應(yīng)用價值[3-4]。GMTI系統(tǒng)沿航向編隊時基線性能最優(yōu)[5]，但實際上穩(wěn)定的沿航向編隊難以實現(xiàn)，其他編隊如Cartwheel[6]在運行過程中順軌基線變化很大，切航跡基線亦較大，不適合GMTI任務(wù)。相比于分布式SAR系統(tǒng)[7-8]周期性變化的基線，繩系InSAR系統(tǒng)可在較長時間內(nèi)提供相對穩(wěn)定的順軌基線，且在軌展開時間較短。

在以往研究中，空間系繩被應(yīng)用到不同的太空任務(wù)中，如軌道轉(zhuǎn)移[9]、軌道維持[10]、太空艙返回[11]和衛(wèi)星離軌[12]等。部分任務(wù)對系繩分布方向不作要求，但多數(shù)系統(tǒng)希望系繩能夠沿鉛垂方向分布，以利用重力梯度力矩保持穩(wěn)定。系繩式InSAR系統(tǒng)執(zhí)行GMTI任務(wù)需要沿水平方位分布以使切航跡基線最小化，相比其他任務(wù)，其對系統(tǒng)空間構(gòu)型要求獨特。在系繩式InSAR任務(wù)中，系統(tǒng)快速且穩(wěn)定地成功展開是完成任務(wù)的前提。盡管針對空間繩系SAR任務(wù)的研究較少，但已有諸多學者針對普適性繩系衛(wèi)星系統(tǒng)，尤其是展開控制[13-17]方面開展了大量研究，提出了繩長速率控制、系繩張力反饋控制和最優(yōu)展開控制等方法。但多數(shù)研究中系繩較長，系統(tǒng)展開開始階段展開執(zhí)行機構(gòu)不工作，只有當系繩伸展到一定長度后重力梯度力矩能夠克服執(zhí)行機構(gòu)輸出誤差時才按照各種控制律開始工作。然而GMTI任務(wù)需求順軌基線較短，僅為幾百米，甚至更短，系統(tǒng)重力梯度力矩較小，為避免展開提前停止，短系繩系統(tǒng)在展開時展開執(zhí)行機構(gòu)不工作，且應(yīng)考慮系統(tǒng)展開執(zhí)行機構(gòu)的靜摩擦力。此時系統(tǒng)展開對于初始參數(shù)的依賴增強，而以往關(guān)于空間系繩的展開較少關(guān)注系統(tǒng)的初始參數(shù)選擇問題。

本文針對面向GMTI任務(wù)的短繩系InSAR系統(tǒng)，考慮系繩的彈性和擺動，使用系統(tǒng)三維運動模型，確定系統(tǒng)的平衡位置，分析系統(tǒng)在任務(wù)要求位置附近的穩(wěn)定性。然后結(jié)合系統(tǒng)任務(wù)分析系統(tǒng)展開要求，確定對系統(tǒng)展開初始參數(shù)的要求，使用粒子群算法對初始參數(shù)進行優(yōu)選。展開完成后使用被動阻尼器控制系繩的高頻振蕩，使用噴氣控制系統(tǒng)的面內(nèi)擺動，將其穩(wěn)定在水平位置附近以獲取滿足任務(wù)要求的順軌基線。

1 系繩式InSAR動力學模型

在本文研究的系繩式InSAR系統(tǒng)中，主星和子星的質(zhì)量基本相當，且遠大于系繩質(zhì)量，系統(tǒng)質(zhì)心可近似認為在兩顆衛(wèi)星質(zhì)心連線上。系繩長度與需求基線長度相當，由GMTI任務(wù)需求確定。不過一般GMTI要求順軌基線較短，故系繩長度亦較短，系繩所受大氣阻力較小，產(chǎn)生的大氣阻力力矩也較小，可作為擾動處理。同時為了簡化系統(tǒng)動力學模型，作出如下合理假設(shè)：將衛(wèi)星視為質(zhì)點，認為衛(wèi)星姿態(tài)可控，不會因為衛(wèi)星姿態(tài)失穩(wěn)滾轉(zhuǎn)而導致展開失?。幌到y(tǒng)質(zhì)心軌道為Keplerian軌道；忽略系繩的扭轉(zhuǎn)和彎曲。

系繩式InSAR系統(tǒng)簡化模型如圖1所示，其中c-xyz為系統(tǒng)質(zhì)心軌道坐標系，Mm和Ms分別表示主星和子星，v為系統(tǒng)質(zhì)心速度，θ和φ分別為系統(tǒng)的面內(nèi)擺角和面外擺角。在系統(tǒng)質(zhì)心軌道坐標系c-xyz中，其動力學方程可以表示為[11]

(1)

(2)

(3)

系統(tǒng)質(zhì)心軌道方程為

(4)

式中：μ為地球引力常數(shù)；a為系統(tǒng)質(zhì)心軌道半長軸。

圖1 系繩式InSAR系統(tǒng)簡化模型Fig.1 Simplified model of tethered InSAR system

2 水平方位穩(wěn)定性分析

系繩式InSAR系統(tǒng)執(zhí)行GMTI任務(wù)的理想位置是沿水平方向分布。為獲取系統(tǒng)在水平方向附近的穩(wěn)定性，為后續(xù)展開和控制提供定性分析，首先求解系統(tǒng)的平衡位置，并分析系統(tǒng)在平衡位置附近的穩(wěn)定性。系統(tǒng)處于穩(wěn)定狀態(tài)時，系統(tǒng)狀態(tài)為

(5)

θe=ce1,φe=ce2,Λe=ce3

(6)

式中：θe、φe和Λe分別為系統(tǒng)處于平衡位置時的系統(tǒng)狀態(tài)，式(5)中為與之對應(yīng)狀態(tài)的速度和加速度，下標e表示系統(tǒng)處于平衡位置時的狀態(tài)；cei(i=1,2,3)為任意常數(shù)。

此時系統(tǒng)不受外力作用，系統(tǒng)動力學方程可寫為

2esinf=3sinθecosθe

(7)

sinφecosφe=0

(8)

(9)

式中：系繩張力可以表示為Te=EA(Λe-1)，E為系繩楊氏模量，A為系繩橫截面積，EA為系繩剛度。

由式(8)可知，系統(tǒng)面外平衡角φe與其他狀態(tài)量無關(guān)，故可首先確定系統(tǒng)的面外擺角平衡點為：φe=±iπ/2(i為任意整數(shù))。i為偶數(shù)時，cos2φe=1；i為奇數(shù)時，cos2φe=0。根據(jù)式(9)，為使系繩保持張緊，需Te≥0 N，i只能為偶數(shù)，因此處于平衡位置時系統(tǒng)在質(zhì)心軌道平面內(nèi)。

由式(7)可知，系統(tǒng)的面內(nèi)擺動平衡點與軌道偏心率有關(guān)，且隨著系統(tǒng)質(zhì)心真近點角的變化而變化。當e=0，即系統(tǒng)質(zhì)心軌道為圓軌道時為一特例，此時系統(tǒng)面內(nèi)擺動平衡點為：θe=±iπ/2(i為任意整數(shù))，即平衡位置沿鉛垂或水平方向分布。i為偶數(shù)，系繩沿鉛垂方向分布時，有

當i為奇數(shù)，即系繩沿水平方向分布時，Te=0 N。可見當系統(tǒng)沿水平方向分布時系繩不起作用。

由式(7)可知，系統(tǒng)質(zhì)心軌道為橢圓軌道時，系統(tǒng)平衡位置不再沿鉛垂或水平方向分布，且系統(tǒng)質(zhì)心軌道偏心率越大，整個軌道周期內(nèi)系統(tǒng)平衡位置偏離鉛垂或水平方向的角度越大。系統(tǒng)平衡位置與軌道偏心率及真近點角的關(guān)系見圖2，圖中平衡位置偏離角度與式(7)一致，另外每個時刻系統(tǒng)的平衡位置有4個，圖中僅給出θe在0° 和90° 附近的變化情況，其在180° 和-90° 附近變化規(guī)律亦如此，此處略去。

圖2 系統(tǒng)平衡位置與軌道偏心率及真近點角的關(guān)系Fig.2 Relationship between equilibrium with ecctricity and true anomaly

系統(tǒng)平衡位置若要在水平方向附近，系統(tǒng)軌道偏心率應(yīng)較小。f∈(0,π/2)∪(3π/2,2π)時，Te>0 N；f∈(π/2,3π/2)時,系繩沿鉛垂方向分布時Te>0 N，系繩沿水平方向分布時Te≤0 N。

GMTI任務(wù)需要獲取順軌基線，最理想情況是系統(tǒng)沿水平方向分布。系統(tǒng)處于平衡位置時系統(tǒng)的面外擺角為0°，系統(tǒng)在質(zhì)心軌道面內(nèi)存在4個平衡位置。為簡化分析，只考慮系統(tǒng)面內(nèi)運動，此時系統(tǒng)狀態(tài)變量為

將系統(tǒng)在平衡位置小范圍近似線性化，可得

(10)

式中：

上述方程中包含系繩質(zhì)量的項較多，而系繩質(zhì)量相比于衛(wèi)星質(zhì)量較小，即mt/m?10，忽略系繩質(zhì)量可以簡化分析。此時系數(shù)矩陣M為

(11)

系統(tǒng)特征方程為

[λ(λ-2Q3)+3Q1][λ(λ-2Q3)-Q4]=0

(12)

可得系統(tǒng)特征值為

(13)

式中：

由式(13)可見，運行在橢圓軌道的系統(tǒng)在平衡位置附近的穩(wěn)定性不僅與系統(tǒng)的真近點角相關(guān)，而且還和系統(tǒng)平衡位置相關(guān)，系統(tǒng)的平衡位置又與軌道偏心率和真近點角相關(guān)。當系統(tǒng)平衡位置在水平方向附近時，cos 2θe<0，可得H1>e2sin2f，故系統(tǒng)特征值必有正值，因此系統(tǒng)沿水平方向分布時不穩(wěn)定。然而面向GMTI任務(wù)的系繩式InSAR需要沿水平方向分布，因此最好能夠?qū)⑾到y(tǒng)展開到水平方向并施加控制，使其在水平方向附近運動。

3 系統(tǒng)面向GMTI任務(wù)的展開分析

本文研究對象為攜帶X波段的SAR系統(tǒng)，為高精度測量地面運動目標，需要順軌基線約100 m。因此系繩展開長度設(shè)定為100 m，由于系繩長度較短，系統(tǒng)重力梯度力矩較小，若使用滾輪等展開機構(gòu)容易因為其輸出誤差導致系統(tǒng)展開阻力較大，甚至提前停止，故本文選擇自由展開方式。系統(tǒng)展開流程為：展開前，兩顆衛(wèi)星結(jié)合在一起，通過彈射分離機構(gòu)使系統(tǒng)分離，系統(tǒng)自由展開，當系統(tǒng)展開水平方位時施加控制，使系統(tǒng)停留在水平方位附近。為了使系統(tǒng)快速展開到期望位置，需要對系統(tǒng)彈射分離的速度和方向進行優(yōu)選。系統(tǒng)展開時展開機構(gòu)存在靜摩擦力，初始分離速度需要大于某一個值，以避免在展開完成前停止，但并不是越大越好，較大的初始速度會導致展開完成時剩余速度較大，引發(fā)系繩縱向高頻振蕩，嚴重時系統(tǒng)會發(fā)生碰撞，因此需要選擇一個盡量小的初始分離速度。由于系統(tǒng)的面外運動和面內(nèi)運動耦合較弱，面外初始擺角較小時，后續(xù)的面外擺角也較小且逐漸衰減，故假定系統(tǒng)面外擺角和角速度均零，即φ=0 rad，φ′=0 rad/s，系統(tǒng)的運動方程簡化為面內(nèi)運動方程，即

(14)

(15)

為使系統(tǒng)在無控條件下均能完成展開，系統(tǒng)的展開速度應(yīng)大于某一值，使其能夠在最惡劣情況下展開到目標長度時速度剛好衰減到零，此時的初始展開速度能夠滿足任意初始角度。最惡劣情況下系繩變化方程為

參加會議的費用，一般包括會務(wù)費（少數(shù)會議沒有此費用）、交通費及住宿費等，少則幾百，多則幾千，可以利用個人或團隊的科研經(jīng)費，也可向圖書館領(lǐng)導申請。一般和本職工作密切相關(guān)的會議或有會議征文獲獎的，館領(lǐng)導都會同意出資。例如筆者是采編部主任，所參加的文獻資源建設(shè)方面的相關(guān)會議都是館里資助的，其他會議大多是走個人或團隊的科研經(jīng)費。

(16)

完成系統(tǒng)初始分離速度的優(yōu)選后，還需要對初始分離角度進行優(yōu)選，以使系統(tǒng)展開完成時剛好擺動到水平方向。選擇系統(tǒng)初始分離方向優(yōu)選指標為

J=θf-π/2

(17)

式中：θf為系統(tǒng)展開結(jié)束時的面內(nèi)擺角。

在動力學方程式(14)和式(15)約束下對系統(tǒng)初始分離角度進行優(yōu)選，使性能指標式(17)最小。

4 基于粒子群算法的初值優(yōu)選

粒子群算法[18-19]是一種模擬自然生物群體行為的啟發(fā)式算法，基于群體單位間相互學習行為。群體中每個個體稱為“粒子”，每個粒子代表一個潛在可行解，群體在解空間上搜索最優(yōu)解，且每個粒子每次搜索都有與之對應(yīng)的適應(yīng)度數(shù)值，可以用來調(diào)整搜索方向，另外群體能夠記憶每個粒子的歷史最優(yōu)解和群體歷史最優(yōu)解，每個粒子通過不斷學習自身最優(yōu)解和群體最優(yōu)解而不斷接近最優(yōu)解。

以在D維空間上搜索的群體為例，粒子i的位置信息表示為Pi=[pi1pi2…piD]T，速度為Vi=[vi1vi2…viD]T。第k次搜索粒子i的位置和速度更新公式為

(18)

(19)

式中：Phi為粒子i的歷史最優(yōu)解；Pg為群體歷史最優(yōu)解；w為慣性權(quán)值；c1和c2為學習因子，用來調(diào)節(jié)粒子i向自身歷史最優(yōu)解和群體歷史最優(yōu)解的靠近速度；cr為隨機值。

使用粒子群算法求解最優(yōu)解的步驟如下：

步驟1初始化。在搜索空間上隨機產(chǎn)生N個粒子(N為群體數(shù)量)的初始位置和速度，并計算每個粒子的適應(yīng)度，將當前位置設(shè)為粒子的最優(yōu)解，對比群體中所有粒子的適應(yīng)度值，將適應(yīng)度值最小的粒子位置作為群體歷史最優(yōu)解初值。

步驟2更新粒子的位置和速度。按照式(18)和式(19)更新群體粒子的位置和速度。

步驟3評價更新粒子。計算更新后每個粒子的適應(yīng)度值，如果好于對應(yīng)粒子歷史最優(yōu)解，則將當前粒子位置更新為歷史最優(yōu)解。將本次更新群體中最好粒子對應(yīng)的適應(yīng)度與群體歷史最優(yōu)解對應(yīng)的適應(yīng)度比較，如果效果好于后者，則將對應(yīng)粒子位置更新為群體歷史最優(yōu)解。

步驟4結(jié)果判定。若達到最大迭代次數(shù)或滿足所需適應(yīng)度，則終止迭代，否則返回步驟2。

粒子群算法在解空間上搜索，相比遺傳算法無需編碼、解碼過程，設(shè)置參數(shù)少，調(diào)整方便，易于實現(xiàn)。

5 系統(tǒng)面向GMTI任務(wù)的穩(wěn)態(tài)控制

由于優(yōu)選的最小初始分離速度是按照最壞情況得到的，實際情況中系繩展開到目標長度后依然有殘留速度。雖然這避免了展開末段耗時過長，但也會使系繩產(chǎn)生縱向振蕩。根據(jù)前面分析知，系統(tǒng)面內(nèi)運動耦合較強，為使系統(tǒng)能夠穩(wěn)定在水平方向附近，需對系繩的高頻振動施加阻尼，使其快速衰減。為了節(jié)省能量，此處使用被動阻尼控制方法，通過彈簧阻尼器對系繩振動進行阻尼。為了將系統(tǒng)徑向振動與阻尼器之間的能量轉(zhuǎn)移最大化，使阻尼器的固有頻率與系繩的固有頻率相等[20]。系統(tǒng)展開到目標長度后其徑向振動固有頻率為

(20)

設(shè)計阻尼器參數(shù)為

(21)

式中：ξ為阻尼器阻尼系數(shù)。

通過對系統(tǒng)初始分離角度的優(yōu)選使系統(tǒng)在展開到目標長度的同時擺動到水平位置，然而此時系統(tǒng)的面內(nèi)擺動角速度不為零，需要對系統(tǒng)施加控制力矩使系統(tǒng)的擺動角速度變?yōu)榱?，此處選用噴氣控制提供控制力矩，由于控制時間較短，可近似視為脈沖控制。另外系繩無法提供彎曲力矩，需要兩顆衛(wèi)星同時噴氣來分別控制，每顆衛(wèi)星所需改變動量為

(22)

運行在圓軌道的系統(tǒng)沿水平方向是臨界穩(wěn)定的，系統(tǒng)狀態(tài)稍微偏離平衡點，系統(tǒng)便會偏離水平方位，而實際中無法使系統(tǒng)剛好位于平衡位置。而運行于橢圓軌道的系統(tǒng)沿水平方向分布是不穩(wěn)定的，在后續(xù)運行過程中系統(tǒng)會偏離水平方位，且系統(tǒng)切軌分量會隨著偏離角度的增大而增大，當系統(tǒng)切軌分量較大時會對GMTI任務(wù)造成干擾，因此需要對系統(tǒng)的面內(nèi)擺角施加控制。由于系統(tǒng)位于水平方向附近時角速度較小，若對系統(tǒng)面內(nèi)擺角進行實時控制，所需控制量亦較小，實際中不易執(zhí)行，可選擇離散脈沖控制，當系統(tǒng)偏離水平方位某一角度時施加控制力矩，使系統(tǒng)回到水平方位。

本文中設(shè)定系統(tǒng)最大且航跡基線為10 m，則相對于100 m的系繩，系統(tǒng)容許面內(nèi)擺角偏離值約為5°，故當系統(tǒng)偏離水平方位5° 時對系統(tǒng)施加反方向控制力矩使其獲得大小相等的反向自轉(zhuǎn)角速度，以回到水平位置。

6 仿真驗證

設(shè)定系統(tǒng)攜帶X波段SAR，執(zhí)行GMTI任務(wù)時需要順軌基線100 m，其他系統(tǒng)參數(shù)見表1。

表1 仿真中的系統(tǒng)參數(shù)Table 1 System parameters used in simulation

使用優(yōu)選的初始參數(shù)對系統(tǒng)進行展開，展開到目標長度后對系統(tǒng)進行穩(wěn)定控制。結(jié)合文獻[21]，選擇阻尼器阻尼系數(shù)ξ=0.9。數(shù)值仿真結(jié)果見圖3～圖6。

圖3為系統(tǒng)面內(nèi)擺動圖，其中實線部分為系統(tǒng)展開過程中面內(nèi)擺角變化，可見在該優(yōu)選初始狀態(tài)下系統(tǒng)展開到目標長度時恰好擺動到水平方向，所優(yōu)選的初值能夠滿足要求。隨后使用噴氣控制系統(tǒng)的面內(nèi)擺角，使系統(tǒng)處于水平方向附近，面內(nèi)擺角85° 和95° 為控制點，控制周期約為0.7個軌道周期。

圖4為系統(tǒng)順軌基線變化情況，系統(tǒng)穩(wěn)態(tài)控制階段，系統(tǒng)順軌基線周期性變化，當系統(tǒng)偏離水平方向5° 時系統(tǒng)順軌基線最短，最小值大于99.6 m，其變化規(guī)律與穩(wěn)態(tài)階段的系統(tǒng)面內(nèi)擺角變化情況對應(yīng)。

圖3 系統(tǒng)面內(nèi)擺角變化Fig.3 In-plane angle variation of system

圖4 系統(tǒng)順軌基線變化Fig.4 Length variation of along-track baseline of system

圖5為兩星徑向速度的變化情況，可知系統(tǒng)在展開過程中速度一直減小，這是因為考慮了系統(tǒng)釋放機構(gòu)的靜摩擦力。但展開完成時相對速度不為零，這避免了系統(tǒng)末端展開過慢而導致耗時較長，同時也不會因為殘余速度過大而導致大幅度振蕩。在阻尼器作用下，系統(tǒng)徑向高頻振蕩很快衰減，穩(wěn)態(tài)控制階段系統(tǒng)徑向速度非常小。

圖5 系統(tǒng)徑向速度變化Fig.5 Variation of relative velocity along tether

圖6 展開完成后系繩張力變化Fig.6 Variation of tether tension after deployment

圖6為系繩張力變化情況，展開剛結(jié)束時由于徑向殘余速度導致系繩張力較大，不過隨著徑向高頻振蕩的衰減而迅速減小。在穩(wěn)態(tài)控制階段，系繩張力類周期變化，主要受真近點角影響，與前面的定性分析一致。由圖6的放大圖可知，系統(tǒng)在平衡位置附近時系繩張力會為零，意味著系繩可能出現(xiàn)松弛現(xiàn)象。這是因為在噴氣控制下系統(tǒng)會擺向平衡位置，但是系統(tǒng)接近平衡位置時其相對系統(tǒng)質(zhì)心轉(zhuǎn)動角速度變小，導致系統(tǒng)離心力減小，另外系繩在恢復力作用下收縮，且收縮到標稱長度時速度不為零，導致系繩長度小于標稱長度，系繩張力為零。與此同時系統(tǒng)的轉(zhuǎn)動角速度會因為系繩長度的減小而增大，在無控狀態(tài)下系繩會逐漸張緊，不會出現(xiàn)系繩張力一直為零的情況。另外由動力學方程式(1)～式(3)可知，系繩張力為零時不會對系統(tǒng)姿態(tài)產(chǎn)生不利影響，反而會減弱系繩張力對衛(wèi)星姿態(tài)的影響。

由圖3～圖6可知，使用優(yōu)選得到的初始參數(shù)，系統(tǒng)能夠快速展開到目標長度且恰好擺動到水平方向，衛(wèi)星相對速度亦減小到約0.06 m/s；展開后的阻尼控制能夠使系繩高頻振蕩快速衰減，噴氣控制使系統(tǒng)維持在水平方向附近，且控制頻率約為0.7個軌道周期。系統(tǒng)的順軌基線最小為99.6 m，能夠滿足任務(wù)要求。

7 結(jié) 論

本文研究了面向GMTI任務(wù)的系繩式InSAR系統(tǒng)展開問題，由于系繩較短，系統(tǒng)重力梯度力矩較小，假若展開機構(gòu)工作可能會因控制誤差而導致展開提前停止。故本文在無控狀態(tài)下展開系統(tǒng)。結(jié)合GMTI任務(wù)，使用粒子群算法對系統(tǒng)初值進行優(yōu)選，使其展開到目標長度時恰好擺動到水平方向，并將其控制到水平方向附近。最后通過數(shù)值仿真驗證優(yōu)選的初始參數(shù)和穩(wěn)定控制方法。結(jié)論如下：

1) 使用粒子群優(yōu)選的初值能夠滿足要求，在該初值下無控展開，系統(tǒng)展開速度不會在展開完成前減小到零，且系統(tǒng)在展開到目標長度時恰好擺動到水平位置。

2) 系統(tǒng)展開完成后，使用被動阻尼控制器噴氣控制分別控制系統(tǒng)的徑向高頻振蕩和系統(tǒng)面內(nèi)擺動。系統(tǒng)高頻振蕩很快衰減，而面內(nèi)擺角始終在水平方向附近，獲取的順軌基線大于99.6 m，能夠滿足任務(wù)要求。

[1] OUCHI K. Recent trend and advance of synthetic aperture radar with selected topics[J]. Remote Sensing, 2013, 5(2): 716-727.

[2] MOREIRA A, IRAOLA P P, YOUNIS M, et al. A tutorial on synthetic aperture[J]. IEEE Geoscience and Remote Sensing Magazine, 2013, 1(1): 6-43.

[3] KAHLE R, RUNGE H, ARDAENS J S, et al. Formation flying for along-track interferometric oceanography—First in-flight demonstration with TanDEM-X[J]. Acta Astronautica, 2014, 99: 130-142.

[4] MOCCIA A, RUFINO G. Spaceborne along-track SAR interferometry: Performance analysis and mission scenarios[J]. IEEE Transactions on Aerospace and Electronic Systems, 2001, 37(1): 199-213.

[5] 梁甸農(nóng), 蔡斌, 王敏, 等. 星載SAR-GMTI研究進展[J]. 國防科技大學學報, 2009, 31(4): 87-92.

LIANG D N, CAI B, WANG M, et al. Research process of spaceborne SAR-GMTI systems[J]. Journal of National Univeristy of Defense Technology, 2009, 31(4): 87-92 (in Chinese).

[6] MASSONNET D. Capabilities and limitations of the interferometric Cartwheel[J]. IEEE Transactions on Geoscience and Remote Sensing, 2001, 39(3): 506-520.

[7] MOCCIA A, VETRELLA S. A tethered interferometric synthetic aperture radar (SAR) for a topographic mission[J]. IEEE Transactions on Geoscience and Remote Sensing, 1992, 30(1): 103-109.

[8] BOMBARDELLI C, LORENZINI E C, QUADRELLI M B. Retargeting dynamics of a linear tethered interferometer[J]. Journal of Guidance, Control, and Dynamics, 2004, 27(6): 1061-1067.

[9] 鐘睿, 徐世杰. 基于直接配點法的繩系衛(wèi)星系統(tǒng)變軌控制[J]. 航空學報, 2010, 31(3): 572-578.

ZHONG R, XU S J. Orbit-transfer control for TSS using direct collocation method[J]. Acta Aeronautica et Astronautica Sinica, 2010, 31(3): 572-578 (in Chinese).

[10] 劉剛, 李傳江, 馬廣富. 應(yīng)用非線性模型預測控制的繩系衛(wèi)星Halo軌道保持控制[J]. 航空學報, 2014, 35(9): 2605-2614.

LIU G, LI C J, MA G F. Station-keeping of tethered satellite system around a halo orbit using nonlinear model predictive control[J]. Acta Aeronautica et Astronautica Sinica, 2014, 35(9): 2605-2614 (in Chinese).

[11] WILLIAMS P, HYSLOP A, STELZER M, et al. YES2 optimal trajectories in presence of eccentricity and aerodynamic drag[J]. Acta Astronautica, 2009, 64(7): 745-769.

[12] ZHONG R, ZHU Z H. Optimal control of nanosatellite fast deorbit using electrodynamic tether[J]. Journal of Guidance, Control, and Dynamics, 2014, 37(4): 1182-1194.

[13] YU S H. Range-rate control algorithms and space rendezvous schemes[J]. Journal of Guidance, Control, and Dynamics, 1997, 20(1): 206-208.

[14] GLABEL H, ZIMMERMANN F, BRUCKNER S, et al. Adaptive neural control of the deployment procedure for tether-assisted re-entry[J]. Aerospace Science and Technology, 2004, 8(1): 73-81.

[15] WILLIAMS P. Optimal deployment/retrieval of tethered satellites[J]. Journal of Spacecraft and Rockets, 2008, 45(2): 324-343.

[16] SUN G H, ZHU Z H. Fractional order tension control for stable and fast tethered satellite retrieval[J]. Journal of Guidance, Control, and Dynamics, 2014, 37(6): 2062-2066.

[17] ASLANOV V S, LEDKOV A S. Dynamics of tethered satellite systems[M]. Oxford: Woodhead Publishing, 2012: 65-88.

[18] 徐鶴鳴. 多目標粒子群優(yōu)化算法的研究[D]. 上海: 上海交通大學, 2013: 4-12.

XU H M. Research on multi-objective particle swarm optimization algorithms[D]. Shanghai: Shanghai Jiao Tong University, 2013: 4-12 (in Chinese).

[19] 沈伋, 韓麗川, 沈益斌. 基于粒子群算法的飛機總體參數(shù)優(yōu)化[J]. 航空學報, 2008, 29(6): 1538-1541.

SHEN J, HAN L C, SHEN Y B. Optimization of airplane primary parameters based on particle swarm algorithm[J]. Acta Aeronautica et Astronautica Sinica, 2008, 29(6): 1538-1541 (in Chinese).

[20] LEE N N, ZORN A H, WEST M. Passive vertical stabilization of two tethered nanosatellites with engineered damping[C]//AIAA/AAS Astrodynamics Specialist Conference and Exhibit. Reston: AIAA, 2008: 1-16.

[21] LORENZINI E C. A three-mass tethered system for micro-g/variable-gapplications[J]. Journal of Guidance, Control, and Dynamics, 1987, 10(3): 242-249.

張錦繡男， 博士， 教授， 博士生導師。主要研究方向： 分布式航天器系統(tǒng)設(shè)計與仿真， 系繩式衛(wèi)星系統(tǒng)應(yīng)用研究， 空間垃圾天基主動移除。

Tel.: 0451-86403440-8308

E-mail: jinxiu@hit.edu.cn

張志剛男， 博士研究生。主要研究方向： 航天器動力學與控制， 系繩式InSAR系統(tǒng)。

Tel.: 0451-86402357-8503

E-mail: zhzhghit@126.com

URL：www.cnki.net/kcms/detail/11.1929.V.20160627.0836.002.html

DeploymentresearchoftetheredInSARsystemforGMTImissions

ZHANGJinxiu*,ZHANGZhigang

SchoolofAstronautics,HarbinInstituteofTechnology,Harbin150001,China

Thedeploymentoftetheredinterferometricsyntheticapertureradar(InSAR)systemappliedforgroundmovingtargetindication(GMTI)missionsisresearchedinthepaper.Toascertainthecharacteristicofdeployment,thethree-dimensionaldynamicsmodelisapplied,andthesteadystateanalysisispresented.Accordingtothesystemfeature,thefreedeploymentmethodischosentoavoidhaltcausedbytheoutputerrorofcontrolmechanism.Asaresult,thedeploymentissensitivetotheinitialstateofsystem.Hence,particleswarmalgorithmisusedtooptimizetheinitialstatetodeploythesystemalonghorizontaldirectionwhenthedeploymentfulfills.Thenapassivedamperandjetscontrolareusedtomakethesystemstaynearhorizontaldirectioninlong-term.Finally,thenumericalsimulationisexecutedtoverifytheoptimizedinitialstateandthesteadycontrolmethod.Theresultsshowthatthesystemcoulddeploytothetargetpositionwiththeoptimizedinitialstate,andthealong-trackbaselineofsystemisalwaysbiggerthan99.6mwiththestablecontrol.

tetheredInSARsystem;GMTImission;steadystateanalysis;particleswarmoptimization;stablecontrol

2016-01-11；Revised2016-04-11；Accepted2016-06-14；Publishedonline2016-06-270836

NationalNaturalScienceFoundationofChina(91438202)

.Tel.：0451-86403440-8308E-mailjinxiu@hit.edu.cn

2016-01-11；退修日期2016-04-11；錄用日期2016-06-14； < class="emphasis_bold">網(wǎng)絡(luò)出版時間

時間：2016-06-270836

www.cnki.net/kcms/detail/11.1929.V.20160627.0836.002.html

國家自然科學基金 (91438202)

.Tel.：0451-86403440-8308E-mailjinxiu@hit.edu.cn

張錦繡， 張志剛． 面向GMTI任務(wù)的系繩式InSAR系統(tǒng)展開研究J． 航空學報,2016,37(10):3083-3091.ZHANGJX,ZHANGZG.DeploymentresearchoftetheredInSARsystemforGMTImissionsJ.ActaAeronauticaetAstronauticaSinica,2016,37(10):3083-3091.

http://hkxb.buaa.edu.cnhkxb@buaa.edu.cn

10.7527/S1000-6893.2016.0186

V476.5

A

1000-6893(2016)10-3083-09