【正文】 理一，重力轉(zhuǎn)彎軟著陸的最優(yōu)制導(dǎo)律是一種開關(guān)（BangBang）控制，只須控制發(fā)動機(jī)開關(guān)，不需要調(diào)節(jié)推力的大小。又因為不與此時由（6b）式有反證假設(shè)矛盾。根據(jù)式及性質(zhì)2）可知，由性質(zhì)3）必有根據(jù)是時間t的斜率非零的線性函數(shù)，）若定，根據(jù)橫截條件有在區(qū)間內(nèi)為常數(shù)。用反證法，假設(shè)存在奇異條件，則在某個閉區(qū)間設(shè)，并由（5）式得。月球重力轉(zhuǎn)彎軟著陸系統(tǒng)（2）的燃耗最優(yōu)制導(dǎo)或時間最優(yōu)制導(dǎo)問題不存在奇異條件。又由（9）式可得T（t）=0，4）根據(jù)極大值原理，系統(tǒng)的狀態(tài)變量和共軛變量都是時間的連續(xù)可微函數(shù)，將切換函數(shù)對時間求導(dǎo)，利用（2），（6）式和性質(zhì)2）得 軟著陸最優(yōu)控制中奇異條件的分析對于月球重力轉(zhuǎn)彎軟著陸問題，最優(yōu)制導(dǎo)律具有兩個很好的性質(zhì)。此時如果最優(yōu)解存在，則稱為奇異解，（8）式稱為奇異條件。根據(jù)pontryagin極大值原理，系統(tǒng)的哈密頓函數(shù)及其對u的偏導(dǎo)數(shù)為使哈密頓函數(shù)（5）式達(dá)到極大地控制輸入u就是最優(yōu)控制，科表示為。軟著陸燃耗最優(yōu)問題的描述 對于最終著陸段，可假設(shè)為一小角度。重力轉(zhuǎn)彎過程中，探測器的高度、速度和姿態(tài)角度可由雷達(dá)高度表、多普勒雷達(dá)及慣性儀表測得。重力轉(zhuǎn)彎軟著陸過程中探測器質(zhì)心動力學(xué)方程可表示為上式中各變量的物理意義如圖1中所示，其中m0為探測器質(zhì)量；k0為制動發(fā)動機(jī)比沖；u表示制動發(fā)動機(jī)的秒耗量可通過一定的機(jī)構(gòu)加以調(diào)節(jié)，故作為軟著陸問題的控制變量。上述結(jié)論對上注探測器關(guān)鍵點(diǎn)的選取有著較強(qiáng)的指導(dǎo)意義，比如基于最優(yōu)軌線的斜率對路徑點(diǎn)合并、基于最優(yōu)軌線簇的對稱性對上注軌線進(jìn)行等效延伸、或者嘗試僅將 S 型和 C 型的轉(zhuǎn)折點(diǎn)作為路徑點(diǎn)等，這樣可以大大降低探測器自主存儲與計算需求，進(jìn)而有效提升任務(wù)的可靠性。3)當(dāng)探測器初始水平速度為零時，圓錐體軸線垂直于火星地表，所有最優(yōu)軌線關(guān)于該軸線中心對稱。圖 2 各種不同初始速度對應(yīng)的火星著陸器動力下降段燃料最優(yōu)軌跡簇1)對任意探測器初始位置，特定初始速度對應(yīng)的燃料最優(yōu)著陸軌跡在末端必然收斂到一個固定的近似圓錐體內(nèi)。因此，可行的方案是通過在地面計算大量的燃料最優(yōu)軌跡，并尋找規(guī)律，選取關(guān)鍵路徑點(diǎn)狀態(tài)存儲到著陸器計算機(jī)中，通過在線查表或者在利用對計算量要求較小的反饋制導(dǎo)律完成安全著陸任務(wù)。此外，通過利用如 TOMLAB 等商業(yè)最優(yōu)控制軟件進(jìn)行復(fù)核計算，也驗證了此計算結(jié)果的燃料最優(yōu)性能。由優(yōu)化結(jié)果可以看出，探測器在給定時間飛行并軟著陸到指定位置，且在整個下降過程始終與火星地表保持一定的安全距離，驗證了下降傾角約束的有效性。二階錐優(yōu)化問題可以通過大量免費(fèi)的優(yōu)化工具求解，如 CSDP、DSDP、OpenOpt、SeDuMi、SDPA、SDPLR等。其中探測器各參數(shù)分別取為：m0=2000kg,g=[]ms2,c=2kms,T1=,T2=13kN.。0T1vr)p+1vTz(Φky0+Akg4)+z0,z179。rkp 控制上限:(vzΨk+TT[TTv0]T163。R，二階錐約束參數(shù)維數(shù)n(Ai,bi,ci,di)由相應(yīng)約束確定則式(17)～式(23)可最終轉(zhuǎn)換為如下最優(yōu)化問題： 指標(biāo)函數(shù)：min(vpp)滿足：初值約束:MxΨ0p+Mx(Ψ0y0)+A0g4r0末值約束:MxΨ0p+Mx(Ψ0y0)+A0g4控制約束:Murkp163。R，線性約束參數(shù)D206。b+dinTiTi（k=1,L，n）n*pp其中x206。000000Y=Fy0+Yp+Lg4分別定義如下常值矩陣：最終可得離散化后的燃料最優(yōu)化問題如下： 指標(biāo)函數(shù)：式(9)可表示為邊界條件：式(3)可表示為控制約束：式(10)和式(11)分別可表示為狀態(tài)約束：式(5)和式(12)分別可表示為含有 p個線性約束和 q個二階錐約束的最優(yōu)化問題的標(biāo)準(zhǔn)形式為 指標(biāo)函數(shù)min(lTx)滿足約束DTx+f179。4(n+1)234。234。234。234。MMA+AB+BL3ABB002B000234。234。0234。234。233。233。235。A234。234。M234。M234。2234。2234。234。1234。0234。Y0249。7(n+1)180。235。1235。1235。235。234。n234。234。234。MM234。234。234。234。234。 Y=234。=234。F=234。 ,p=234。234。234。234。234。234。234。234。234。233。233。233。有系統(tǒng)性質(zhì)可知，整個控制時域內(nèi)系統(tǒng)狀態(tài)滿足 y3=Ay2+B(p2+g4)=A3y0+A2B(p0+g4)+AB(p1+g4)+B(p2+g4)Myn=Ayn1+B(pn1+g4)=Any0+An1B(p0+g4)+L+AB(pn2+g4)+B(pn1+g4)y1=Ay0+B(p0+g4)y2=Ay1+B(p1+g4)=A2y0+AB(p0+g4)+B(pn2+g4)+B(p1+g4)為表達(dá)方便，令233。e(Dts)AcBcds=242。4分別為離散系統(tǒng)的系統(tǒng)矩陣和輸入矩陣12A=eDtAc187。7,B206。c=[tanqaltT000000] 等效燃料最優(yōu)精確著陸問題的離散化首先將整個飛行時間均分成 n 段（對應(yīng) n +1 個點(diǎn)），每段步長為Dt，離散化后的著陸器運(yùn)動方程為：yk+1=Ayk+B(pk+g4)其中A206。0010000235。S=234。0（12）其中233。d163。控制約束：由文獻(xiàn)[10]可知，控制約束(4)可等效表示為u163。d],g4=[gTTD0]Tt指標(biāo)函數(shù)：min242。a=Acy+Bc(p+g4)（8）d235。234。233。7*7234。235。235。234。234。I0v233。233。v0（5）進(jìn)一步地，若著陸區(qū)域附近表面崎嶇不平，僅僅確保地表約束不能滿足需求時，可以考慮下降傾角約束，即將著陸器下降軌線約束到以著陸點(diǎn)為頂點(diǎn)的圓錐體內(nèi) 等效后燃料最優(yōu)精確著陸問題 定義等效變換變量Ttrx2+ry2rh163。min242。190。m(t)dtt0f設(shè)計主減速段制導(dǎo)控制律 2動力下降段燃料最優(yōu)精確著陸問題描述 燃料最優(yōu)精確著陸問題著陸器運(yùn)動方程：考慮采用變推力發(fā)動機(jī)情況，有r=v.v=g+a(1)a=Tmm=aT..其中r=[rhrxry]T，v=[vhvxvy]T分別表示著陸器相對期望著陸點(diǎn)的位置和速度矢量；T為推力器提供的推力矢量，幅值為 T，對應(yīng)控制加速度矢量 a；g為火星的重力加速度矢量，此處認(rèn)為是常值；m為著陸器質(zhì)量，對應(yīng)推力器質(zhì)量排除系數(shù)a。優(yōu)化變量為制動發(fā)動機(jī)推力方向角y(t)。根據(jù)動力下降段的起點(diǎn)位置可以確定動力學(xué)方程初始條件,由于起點(diǎn)處于霍曼轉(zhuǎn)移軌道的近地點(diǎn),故其初始條件為: r0=rpq0=0v0=0 w0=1rpmrp(2ra)ra+rp其中rp和ra分別為霍曼轉(zhuǎn)移段的近地點(diǎn)半徑和遠(yuǎn)地點(diǎn)半徑。其示意圖如圖1所示,其中o為月球質(zhì)心,x軸方向為由月心指向著陸器的初始位置,y軸方向為初始位置著陸器速度方向。由于月球表面附近沒有大氣,所以在飛行器的動力學(xué)模型中沒有大氣阻力項。當(dāng)著陸器運(yùn)行到近月點(diǎn)時,制動發(fā)動機(jī)開始工作,其主要任務(wù)是抵消著陸器的初始動能和勢能,使著陸器接觸地面時,相對月面速度為零,即實(shí)現(xiàn)所謂的軟著陸,這一階段稱為動力下降段。162。=43by根據(jù)曲率半徑公式有 r=(1+y162。)=0 (9) 將(8)、(9)兩式聯(lián)立得a2b2y2+a4x2 (10) y162。+yy162。=2 (8)by將(7)式再次對x求導(dǎo)得2a2+2b2(y162。如圖1所示,橢圓的軌跡方程為x2y2+2=1 (5) 2ba將(5)式變形為a2x2+b2y2=a2b2 (6)根據(jù)隱函數(shù)的求導(dǎo)法則將(6)式對x求導(dǎo)有2a2x+2b2yy162。眾所周知,太陽系中的八大行星都在按照各自的橢圓軌道繞太陽進(jìn)行公轉(zhuǎn),太陽位于橢圓的一個焦點(diǎn)上,行星的運(yùn)動遵循開普勒三定律,筆者發(fā)現(xiàn),在各類物理競賽中,常會涉及到天體運(yùn)動速度的計算,本文擬從能量和行星運(yùn)動的軌跡方程兩個不同的角度來探索行星在近日點(diǎn)和遠(yuǎn)日點(diǎn)的速度。即近月點(diǎn)位置坐標(biāo)為(,)海拔15km,遠(yuǎn)月點(diǎn)位置坐標(biāo)為(,)海拔100km。t(T22x+Ty2)=7500Nv2=2at180。0232。m242。0231。a247。49012（）=（沿切線方向）v2=，徑向速度vk=0。R0== R0+R1利用能量平衡式求得近地點(diǎn)速度為2180。R0+R