帶有落角約束的間接Gauss偽譜最優(yōu)制導(dǎo)律

陳琦，王中原，常思江

（南京理工大學(xué)能源與動力工程學(xué)院，江蘇南京210094）

帶有落角約束的間接Gauss偽譜最優(yōu)制導(dǎo)律

陳琦，王中原，常思江

（南京理工大學(xué)能源與動力工程學(xué)院，江蘇南京210094）

針對帶有落角約束的末制導(dǎo)問題，提出了一種基于極小值原理和Gauss偽譜法的最優(yōu)制導(dǎo)律。以期望落角方向為坐標(biāo)軸定義了落角坐標(biāo)系，并在其中建立了線性化的導(dǎo)引運動關(guān)系方程。將控制系統(tǒng)簡化為1階慣性環(huán)節(jié)，利用極小值原理得到正則方程，然后引入Gauss偽譜法進行離散，將其轉(zhuǎn)化為代數(shù)方程，結(jié)合邊界條件，推導(dǎo)出最優(yōu)制導(dǎo)律的解析表達(dá)式，無需任何積分過程，避免了求解黎卡提微分方程。仿真結(jié)果表明，所提出的算法運算量小，計算效率高，同時也能方便地求解出復(fù)雜加權(quán)矩陣下的最優(yōu)制導(dǎo)律，能夠在滿足落角約束的條件下更快地收斂到落角參考線，并且具有更小的末端需用過載。

兵器科學(xué)與技術(shù)；落角約束；最優(yōu)控制；末制導(dǎo)律；極小值原理；間接Gauss偽譜法

0 引言

在導(dǎo)彈實施末端精確打擊時，為了最大限度提高毀傷效果，往往希望導(dǎo)彈以最小的脫靶量命中目標(biāo)的同時，還能以特定的落角攻擊目標(biāo)要害或薄弱的部位。因此，帶有落角約束的末制導(dǎo)律目前已逐漸成為國內(nèi)外學(xué)者的研究熱點［1-4］。文獻(xiàn)［5］利用圓軌跡導(dǎo)引思想，設(shè)計出了一種帶有落角約束的三維制導(dǎo)律；文獻(xiàn)［6］和文獻(xiàn)［7］采用模型預(yù)測靜態(tài)規(guī)劃技術(shù)，得到了可滿足終端落角要求的非線性次優(yōu)制導(dǎo)律；文獻(xiàn)［8］基于滑?？刂萍夹g(shù)，研究了落角約束導(dǎo)引律設(shè)計問題，并同時考慮了自動駕駛儀1階動力學(xué)延遲的情況；文獻(xiàn)［9］提出了一種非奇異終端滑模制導(dǎo)律，避免了傳統(tǒng)終端滑模制導(dǎo)律在控制量輸入飽和時出現(xiàn)奇異現(xiàn)象。

隨著最優(yōu)控制理論的發(fā)展，越來越多的學(xué)者將最優(yōu)控制理論應(yīng)用于帶有落角約束的末制導(dǎo)律的設(shè)計中。文獻(xiàn)［10］將剩余飛行時間的函數(shù)作為性能指標(biāo)，利用最優(yōu)控制理論得到了帶有落角約束的最優(yōu)制導(dǎo)律。文獻(xiàn)［11］利用線性二次型最優(yōu)控制理論研究了帶落角約束的制導(dǎo)律，并與微分對策理論進行了對比，驗證了最優(yōu)控制理論在求解最優(yōu)制導(dǎo)律方面具有較大的優(yōu)勢。文獻(xiàn)［12］通過引入碰撞三角，將導(dǎo)引運動關(guān)系方程進行了線性化，之后依然采用最優(yōu)控制理論研究了最優(yōu)制導(dǎo)律的問題。在將末制導(dǎo)問題轉(zhuǎn)化為最優(yōu)控制問題后，上述研究均采用極小值原理進行求解，通過積分正則方程得到協(xié)態(tài)變量，然后利用極值條件獲得最優(yōu)解。這種方法雖然精度很高，但推導(dǎo)過程較為繁瑣，特別是對于復(fù)雜的問題，顯示求解正則方程將會變得非常困難。為此，文獻(xiàn)［13］通過求解黎卡提微分方程得到了最優(yōu)制導(dǎo)律，這在一定程度上降低了計算難度，但在一般情況下，要得到黎卡提微分方程的解析解依然是一個不小的挑戰(zhàn)?？紤]到這一點，文獻(xiàn)［14］利用黎卡提微分方程的穩(wěn)態(tài)解作為代替，得到了最優(yōu)制導(dǎo)律，但這只是一種近似的處理方法，局限性較大。與上述研究不同，文獻(xiàn)［15］采用Gauss偽譜法，以實時彈道優(yōu)化的思想研究了帶有落角約束的最優(yōu)制導(dǎo)律，完全避免了積分正則方程和求解黎卡提微分方程，很大程度上簡化了推導(dǎo)過程，并取得了不錯的制導(dǎo)效果。但是該制導(dǎo)律需要不斷地調(diào)用優(yōu)化算法求解非線性規(guī)劃問題，計算量較大，求解效率較低，實際應(yīng)用時有一定的局限性。文獻(xiàn)［16］采用解析方法求解了考慮高階自動駕駛儀情況下的最優(yōu)制導(dǎo)律問題，但其只考慮了簡單的加權(quán)矩陣，在某些任務(wù)需求中，如果需要設(shè)計復(fù)雜的加權(quán)矩陣，該方法求解起來則會變得較為繁瑣。

針對以上問題，本文結(jié)合極小值原理和Gauss偽譜法，提出了一種新型帶有落角約束的最優(yōu)制導(dǎo)律，稱為間接Gauss偽譜最優(yōu)制導(dǎo)律。利用極小值原理推導(dǎo)出正則方程，引入Gauss偽譜法將正則方程離散為代數(shù)方程，結(jié)合邊界條件，得到了最優(yōu)制導(dǎo)律的解析表達(dá)式。與其他文獻(xiàn)的制導(dǎo)律相比，文中所提出的方法不需要任何積分或優(yōu)化迭代過程，避免了求解黎卡提微分方程，同時在推導(dǎo)過程中不對加權(quán)矩陣的具體形式進行限制，可以很方便地處理具有復(fù)雜加權(quán)矩陣的最優(yōu)制導(dǎo)律設(shè)計問題，具有較強的靈活性。

1 問題描述

帶有落角約束的平面末制導(dǎo)幾何關(guān)系如圖1所示。其中:M和T分別表示導(dǎo)彈和目標(biāo)；OxIyI為地面坐標(biāo)系；Txfyf為落角坐標(biāo)系，其原點與目標(biāo)固連，Txf軸與期望落角方向重合；θd為期望落角的大??；vM和aM分別為導(dǎo)彈飛行速度及法向過載；θM為地面坐標(biāo)系下的彈道傾角；R和q分別為彈目距離和視線角；θ為落角坐標(biāo)系下彈道傾角，表示落角誤差；y為導(dǎo)彈和Txf軸的距離偏差；為導(dǎo)彈和Tyf軸的距離偏差；為彈目連線和Txf軸的夾角。

圖1 導(dǎo)彈與目標(biāo)幾何關(guān)系圖Fig.1 Geometrical relationship between missile and target

圖1中用于表示角度的箭頭為逆時針方向時，相應(yīng)的角度定義為正，反之為負(fù)。由此可得到如下的角度關(guān)系:

當(dāng)存在落角約束時，在落角坐標(biāo)系Txfyf下建立相應(yīng)的導(dǎo)引運動關(guān)系可以簡化問題的求解難度。為此，結(jié)合圖1可得如下的導(dǎo)引運動方程:

假設(shè)導(dǎo)彈的飛行速度vM為常數(shù)，落角坐標(biāo)系Txfyf下的彈道傾角θ較小，并定義v=vMθ，則（2）式可線性化為

進一步將導(dǎo)彈控制系統(tǒng)簡化為1階慣性環(huán)節(jié)，

式中:aM，c和Tc分別表示為制導(dǎo)指令和動力學(xué)系統(tǒng)時間常數(shù)。結(jié)合（3）式和（4）式，可得落角坐標(biāo)系下的導(dǎo)引運動方程為

式中:x=［yvaM］T；x0=［y0v00］T；

根據(jù)圖1可知，當(dāng)導(dǎo)彈到達(dá)目標(biāo)時，如果速度vM和Txf軸重合，即y（tf）=0，v（tf）=0，其中tf為末端時刻，那么即可滿足終端零脫靶量和期望的落角約束。因此，本文帶有落角約束的末制導(dǎo)問題可以描述為如下的有限時間最優(yōu)控制問題P:在［t0，tf］時間內(nèi)，確定控制量u，使得導(dǎo)彈M在滿足約束（5）式的條件下，從初始狀態(tài)x0轉(zhuǎn)移到Txf軸上，并使如下的性能指標(biāo)最小。

式中:F、Q和R均為正定的對角矩陣，分別表示對末端狀態(tài)量、狀態(tài)量及控制量的加權(quán)。

2 最優(yōu)制導(dǎo)律設(shè)計

針對第1節(jié)中的最優(yōu)控制問題P，構(gòu)造如下形式的Hamiltonian函數(shù):

式中:λ（t）為協(xié)態(tài)向量。根據(jù)正則方程可得

根據(jù)極值條件

可得最優(yōu)控制量為

將（10）式帶入（8）式，并根據(jù)橫截條件，可得如下的兩點邊值問題:

當(dāng)加權(quán)矩陣Q、R比較簡單時，可以通過解析方法獲得問題P的最優(yōu)解［10，16］，但是對于復(fù)雜的加權(quán)矩陣，解析法需要重新推導(dǎo)演算，計算過程較為繁瑣，實際應(yīng)用有諸多不便。此外，反向積分黎卡提微分方程的方法也可以獲得最優(yōu)控制量，但這種方法耗時較長，實際使用時也存在一定的局限性。為此，本文提出了一種間接Gauss偽譜法，通過全局多項式逼近（11）式中的λ（t）和x（t），引入微分矩陣近似相應(yīng)的導(dǎo)數(shù)項，進而將（11）式轉(zhuǎn)換為一系列的代數(shù)方程，求解這些代數(shù)方程便可得到解析形式的最優(yōu)控制量。采用Gauss偽譜法需要將時間區(qū)域［t0，tf］轉(zhuǎn)換到［-1，1］上，為此引入變量子對時間t進行變換，即

因此，（11）式變?yōu)?/p>

由于（13）式中第1式已知初值，第2式已知末值，因此采用Gauss偽譜法對二者進行離散的方式略有不同。對于x（子），用N個Gauss節(jié)點子1，子2，…，子N和初始點子0=-1上的離散狀態(tài)構(gòu)造Lagrange插值多項式近似x（子），

對（14）式微分得

結(jié)合（16）式和（17）式，（13）式的離散格式為

式中:xk≡x（子k）；λk≡λ（子k）.根據(jù)Gauss求積公式［17］，狀態(tài)量的末值x（1）和協(xié)態(tài)量的初值λ（-1）可通過下式得到:

式中:wk為Gauss求積系數(shù)。

證明 引入常數(shù)變量p（t）=c，利用（16）式對其進行離散可得

接下來引入如下形式的分步積分:

式中:tk（k=1，2，…，N）為Gauss節(jié)點。結(jié)合（16）式和（17）式，將（23）式中的導(dǎo)數(shù)項用微分近似矩陣表示，可得

由于證明過程中并未限制f（t）和g（t）的具體形式，因此（24）式對所有的多項式均成立。取f（t）為N次Lagrange插值基函數(shù)Ll（t），l=0，1，…，N（其所使用的N+1個節(jié)點為初始點t0=-1加上N個Gauss節(jié)點）；g（t）為N次Lagrange插值基函數(shù)，j=1，2，…，N+1（其所使用的N+1個節(jié)點為N個Gauss節(jié)點加上末端點tN+1=1）.以上的取法使f（t）和g（t）各有N+1種形式，如果只使用其中的l=1，2，…，N對應(yīng)的Ll（t）和j=1，2，…，N對應(yīng)的，則有

因此可得

因此有

由定理可知，通過微分插值基函數(shù)Lk（t）得到矩陣D后，由（20）式便可非常方便地計算出其他3個矩陣

式中:

將（25）式寫成簡潔形式為

給定ΛN+1和X0，通過（28）式可以解析地求出N個Gauss節(jié)點上的協(xié)態(tài)量λi和狀態(tài)量xi，i=1，2，…，N，將其帶入到（19）式中便可得到端點上的協(xié)態(tài)量λ（-1）和狀態(tài)量x（1）.至此，所有節(jié)點上的協(xié)態(tài)量和狀態(tài)量均已得出，根據(jù)（10）式便可得到每個節(jié)點上的最優(yōu)控制量，如（29）式所示，

從以上的推導(dǎo)過程可以看出，最優(yōu)控制量的計算中不需要任何的積分或迭代過程，只要已知初始偏差X0，即可解析地確定［t0，tf］時間歷程內(nèi)各節(jié)點上的控制量，相比于傳統(tǒng)的反向積分黎卡提微分方程的方法，這可在很大程度上提高運算速度。由于（29）式中的最優(yōu)控制量只依賴于初始偏差，所以（29）式得到的是開環(huán)最優(yōu)解，無法抑制外界干擾的影響。為了解決這一問題，本文引入了滾動時域方法。流程如下:

步驟1 初始化Gauss節(jié)點和初始時刻t0，給定初始狀態(tài)量x0.

步驟2 根據(jù)當(dāng)前狀態(tài)估計剩余飛行時間tgo，得到tf=t0+tgo.

步驟4 只取第1個控制量，即t0時刻對應(yīng)的控制量u0，將其作用于（5）式中得到下一時刻t′以及狀態(tài)量x.

步驟5 判斷是否命中目標(biāo)。如果沒有命中目標(biāo)，以當(dāng)前時刻t′作為初始時刻t0，當(dāng)前狀態(tài)作為初始狀態(tài)x0，重復(fù)步驟2～步驟4，直至命中目標(biāo)。

值得注意的是，（25）式中每個微分矩陣的計算都需要tf的值，因此，步驟2通過估計剩余飛行時間tgo，對tf進行了近似，結(jié)合圖1，tgo采用下式進行估計:

3 算例仿真與分析

為了驗證所提出算法的準(zhǔn)確性并考察其性能，3.1節(jié)將（28）式的計算結(jié)果和目前成熟的最優(yōu)控制問題求解軟件GPOPS的結(jié)果進行了對比。3.2節(jié)在不同條件下對所提出的末制導(dǎo)算法進行了仿真分析。3.3節(jié)將所提出的末制導(dǎo)算法和其他帶有落角約束的最優(yōu)制導(dǎo)律進行了對比分析。

3.1 準(zhǔn)確性驗證

結(jié)合極小值原理和Gauss偽譜法，（28）式給出了最優(yōu)控制問題P的解析解（包括狀態(tài)量和協(xié)態(tài)量），有必要對其準(zhǔn)確性進行驗證。GPOPS是由A. Patterson和Anil V.Rao等學(xué)者開發(fā)的最優(yōu)控制問題通用求解軟件，該軟件利用Gauss積分偽譜方法將連續(xù)時間的最優(yōu)控制問題直接離散為非線性規(guī)劃問題，然后調(diào)用相關(guān)優(yōu)化算法進行求解。目前該軟件廣泛應(yīng)用于軌跡優(yōu)化、交會對接以及軌道轉(zhuǎn)移等領(lǐng)域，并且其計算精度已被諸多學(xué)者所驗證［18-20］。因此，3.1節(jié)選用GPOPS求解最優(yōu)控制問題P，并與（28）式的計算結(jié)果進行比較。問題P的計算參數(shù)假設(shè)如下:t0=0 s；tf=30 s；x0=［500，50，0］T；離散節(jié)點數(shù)目取為40；時間常數(shù)假設(shè)為Tc=1.0 s；性能指標(biāo)中的加權(quán)矩陣取F=diag（1000，1 000，1 000），Q=diag（1/4002，1/1502，1/502），R=diag（1/502）.對比結(jié)果如圖2和圖3所示，其中x和λ表示（28）式的計算結(jié)果，x*和λ*表示GPOPS的計算結(jié)果。

圖2展示了兩種方法得到的狀態(tài)量和協(xié)態(tài)量的對比效果，從中可以看出，二者的計算結(jié)果吻合得非常好。圖3為狀態(tài)量和協(xié)態(tài)量誤差曲線，從中可以明顯地看出，狀態(tài)量的誤差不超過2×10-5，協(xié)態(tài)量的誤差則更小，低于1×10-9.以上對比結(jié)果表明，解析表達(dá)（28）式具有很高的計算精度。此外，GPOPS在求解過程中必須調(diào)用SNOPT或IPOPT進行迭代尋優(yōu)，這一過程會消耗較多的計算資源并產(chǎn)生較長的計算耗時，而（28）式則避免了這種尋優(yōu)過程。表1對比了GPOPS和（28）式在MATLAB平臺下的計算耗時，可以看出（28）式在計算效率上具有非常明顯的優(yōu)勢，這種優(yōu)勢在滾動時域過程中將更為突出。

圖2 狀態(tài)量和協(xié)態(tài)量對比結(jié)果Fig.2 Comparisons of states and costates

表1 兩種方法的計算耗時對比Tab.1 Comparison of computation times of two methods

圖3 狀態(tài)量和協(xié)態(tài)量誤差曲線Fig.3 Errors of states and costates

3.2 不同條件下末制導(dǎo)性能

為了考察所提出算法的性能，對不同條件下的末制導(dǎo)性能進行了仿真分析。具體仿真參數(shù)見表2，其他參數(shù)（離散節(jié)點數(shù)目、時間常數(shù)及權(quán)重矩陣）與3.1節(jié)一致。仿真結(jié)果如圖4～圖6所示。

表2 末制導(dǎo)仿真參數(shù)Tab.2 Simulation parameters for terminal guidance

圖4給出了在θ0=0°，θd=-60°條件下，不同起始距離偏差y0對應(yīng)的彈道曲線及過載曲線。圖4中的落角參考線與圖1中的Txf軸重合，也表示零過載彈道曲線。由圖4可以看出，當(dāng)導(dǎo)彈沿著落角參考線飛行時，便可滿足落角約束。從圖4（a）可以看出，在不同的起始距離偏差下，導(dǎo)彈均能較好地趨向于落角參考線，以期望的落角去攻擊目標(biāo)。從圖4（b）可以看出，制導(dǎo)初始段過載較大，并且其幅值隨y0的增加而增加，這樣可以使得落角更快地向期望值收斂。同時還可以看出，不同情況下的過載最終都會趨向于0，這從另外一方面驗證了導(dǎo)彈在彈道末端將和落角參考線重合，并且在彈道末端隨著需用過載的減小，導(dǎo)彈的抗干擾能力也會隨之增強。

圖4 不同起始距離偏差情況下仿真結(jié)果Fig.4 Simulation results at different initial miss distances

圖5展示了不同初始偏角θ0情況下的制導(dǎo)效果（y0=500 m，θd=-60°）。由圖5可以看出，導(dǎo)彈能夠很好地趨向于落角參考線，可以在滿足落角約束的條件下攻擊目標(biāo)。圖6給出了在不同期望落角θd條件下制導(dǎo)效果。結(jié)合圖6（a）和圖6（b）可知導(dǎo)彈能夠準(zhǔn)確地命中目標(biāo)，并且可以很好地滿足相應(yīng)的落角約束。從圖6（c）可看出，隨著期望落角的絕對值增加，彈道過載的幅值相應(yīng)地變大，但最終均也逐漸收斂到0.

以上的仿真算例選取的控制量權(quán)函數(shù)為常數(shù)，然而在某些情況下，根據(jù)具體的任務(wù)需求，有時會選取時變的控制量權(quán)函數(shù)，如文獻(xiàn)［10］構(gòu)造的權(quán)函數(shù)為剩余時間的冪函數(shù)，文獻(xiàn)［21］研究了一般加權(quán)函數(shù)情況下的最優(yōu)制導(dǎo)律求解問題。復(fù)雜的權(quán)函數(shù)給問題的求解帶來一定的難度，然而在本文所提方法的推導(dǎo)過程中，并未對權(quán)函數(shù)的具體形式有所限制，因此，本文方法可以非常容易地處理復(fù)雜的加權(quán)函數(shù)。為了驗證所提方法對復(fù)雜加權(quán)函數(shù)的求解效果，參考文獻(xiàn)［21］，選擇指數(shù)函數(shù)形式的加權(quán)函數(shù)，即

式中NG為制導(dǎo)系數(shù)。這種形式的加權(quán)函數(shù)給問題的求解帶來了很大的難度，傳統(tǒng)的方法將會遇到很大的困難。而本文方法只需在程序中更改W值即可，其他求解公式無需作任何改變，仿真結(jié)果如圖7所示。

圖5 不同起始偏角情況下仿真結(jié)果Fig.5 Simulation results in the case of different initial heading errors

圖6 不同期望落角情況下仿真結(jié)果Fig.6 Simulation results at different desired impact angles

從圖7（a）可以看出，隨著系數(shù)NG的增加，導(dǎo)彈將更快地趨向于落角參考線。圖7（b）則展示了導(dǎo)彈過載隨NG的變化情況，從中可以看出，隨著NG的增加，開始階段的過載逐漸增加，但末段的過載則更快地趨向于0，這有助于提高導(dǎo)彈的末端過載裕量。不同的加權(quán)函數(shù)可以產(chǎn)生不同的制導(dǎo)效果，針對特定的任務(wù)需求，可以合理地設(shè)計相應(yīng)的加權(quán)函數(shù)。而本文方法則很容易得到復(fù)雜加權(quán)函數(shù)下的最優(yōu)控制量，相比與傳統(tǒng)的方法，具有很強的靈活性。

3.3 與其他算法的比較

為了進一步驗證所提出算法的有效性，將所提算法分別和文獻(xiàn)［15］及文獻(xiàn)［22］中提出的算法進行對比。

文獻(xiàn)［15］以期望落角作為末端約束，能量最優(yōu)作為性能指標(biāo)，采用實時彈道優(yōu)化的思想求解了帶有落角約束的末制導(dǎo)問題，并取得了不錯的效果。為便于對比描述，本文將這種算法記為彈道優(yōu)化制導(dǎo)律。

圖7 不同過載加權(quán)函數(shù)情況下仿真結(jié)果Fig.7 Simulation results in the case of different acceleration weights

文獻(xiàn)［22］所提出的偏置比例制導(dǎo)律算法可描述如下:

式中:uBPN為制導(dǎo)指令，即為導(dǎo)彈過載；η為制導(dǎo)系數(shù)；

vM為導(dǎo)彈飛行速度；為視線角速度；θd為期望落角。根據(jù)文獻(xiàn)［22］，制導(dǎo)系數(shù)選取為NG=4，η=1.3.

各算法的對比結(jié)果如圖8所示。從圖8（a）可以看出，3種算法均能命中目標(biāo)，但是本文方法能以更快的速度趨向于落角參考線。圖8（b）的過載曲線展示出本文方法在制導(dǎo)初始段需要較大的過載，這樣可以使得導(dǎo)彈更快地收斂于落角參考線，這也和圖8（a）中的現(xiàn)象是一致的。此外在彈道末端，偏置比例制導(dǎo)律產(chǎn)生了很大的過載，彈道優(yōu)化制導(dǎo)律次之，而本文算法則最小。由此可見，本文方法在末端具有較大的過載裕量。由圖8（c）可知，3種算法均能保證落角約束，但顯然本文方法能更快地收斂于期望落角。

圖8 不同制導(dǎo)律情況下仿真結(jié)果Fig.8 Simulation results in the case of different guidance laws

4 結(jié)論

本文提出了一種帶有落角約束具有解析形式的間接Gauss偽譜最優(yōu)制導(dǎo)律，并在不同條件下進行了仿真驗證。

1）在落角坐標(biāo)系中建立了線性化的導(dǎo)引運動方程，簡化了問題的求解，運用最優(yōu)控制理論將末制導(dǎo)問題轉(zhuǎn)換為最優(yōu)控制問題。

2）結(jié)合極小值原理和Gauss偽譜法推出了解析形式的最優(yōu)制導(dǎo)律，避免了間接法中的求解黎卡提微分方程和直接法中的尋優(yōu)過程，在很大程度上提高了末制導(dǎo)效率，同時可以處理一類具有復(fù)雜加權(quán)矩陣的最優(yōu)制導(dǎo)律設(shè)計問題，可為此類問題的求解提供了一個新的思路。

3）通過和GPOPS結(jié)果的對比，驗證了所得解析表達(dá)式的準(zhǔn)確性。不同的仿真算例驗證了所提出的算法的有效性。和已有文獻(xiàn)［15，20］相比，所提出的算法能夠在滿足落角約束的條件下更快地收斂到落角參考線，并且具有更小的末端需用過載。

（

）

［1］ Ashwini R，Debasish G.Impact angle constrained guidance against nonstationary nonmaneuvering targets［J］.Journal of Guidance，Control，and Dynamics，2010，33（1）:269-275.

［2］ Sun M W，Xu Q，Du S Z，et al.Practical solution to impact angle control in vertical plane［J］.Journal of Guidance，Control，and Dynamics，2014，37（3）:1022-1027.

［3］ 張友安，馬培蓓.帶有攻擊角度和攻擊時間控制的三維制導(dǎo)［J］.航空學(xué)報，2008，29（4）:1020-1026. ZHANG You-an，MA Pei-bei.Three-dimensional guidance law with impact angle and impact time constraints［J］.Acta Aeronautica et Astronautica Sinica，2008，29（4）:1020-1026.（in Chinese）

［4］ Harrison G A.Hybrid guidance law for approach angle and timeof-arrival control［J］.Journal of Guidance，Control，and Dynamics，2012，35（4）:1104-1114.

［5］ 胡錫精，黃雪梅.具有碰撞角約束的三維圓軌跡制導(dǎo)律［J］.航空學(xué)報，2012，33（3）:508-519. HU Xi-jing，HUANG Xue-mei.Three-dimensional circular guidance law with impact angle constraints［J］.Acta Aeronautica et Astronautica Sinica，2012，33（3）:508-519.（in Chinese）

［6］ 魏鵬鑫，荊武興，高長生.具有落角約束的彈道導(dǎo)彈再入末制導(dǎo)律設(shè)計［J］.哈爾濱工業(yè)大學(xué)學(xué)報，2013，45（9）:23-30. WEI Peng-xin，JING Wu-xing，GAO Chang-sheng.Design of the reentry terminal guidance law with constraint of impact angle for ballistic missile［J］.Journal of Harbin Institute of Technology，2013，45（9）:23-30.（in Chinese）

［7］ Oza H B，Padhi R.Impact-angle-constrained suboptimal model predictive static programming guidance of air-to-ground missiles［J］. Journal of Guidance，Control，and Dynamics，2012，35（1）: 153-164.

［8］ Lee C H，Kim T H，Tahk M J.Design of impact angle control guidance laws via high-performance sliding mode control［J］.Journal of Aerospace Engineering，2013，227（2）:235-253.

［9］ Kumar S R，Rao S，Ghose D.Nonsingular terminal sliding mode guidance with impact angle constraints［J］.Journal of Guidance， Control，and Dynamics，2014，37（4）:1114-1130.

［10］ Ryoo C K，Cho H，Tahk M J.Time-to-go weighted optimal guidance with impact angle constraints［J］.IEEE Transactions on Control Systems Technology，2006，14（3）:483-492.

［11］ Shaferman V，Shima T.Linear quadratic guidance laws for imposing a terminal intercept angle［J］.Journal of Guidance，Control，and Dynamics，2008，31（5）:1400-1412.

［12］ Cho H，Ryoo C K，Tsourdos A，et al.Optimal impact angle control guidance law based on linearization about collision triangle［J］. Journal of Guidance，Control，and Dynamics，2014，37（3）: 958-964.

［13］ 竇磊，楊新民.大著地角衛(wèi)星制導(dǎo)炸彈最優(yōu)制導(dǎo)律研究［J］.南京理工大學(xué)學(xué)報:自然科學(xué)版，2010，34（3）:314-318. DOU Lei，YANG Xin-min.Optimal guidance law of satellite guided bombs with large landing angle［J］.Journal of Nanjing U-niversity of Science and Technology:Natural Science，2010，34（3）:314-318.（in Chinese）

［14］ 侯明善.精確對地攻擊姿態(tài)約束最優(yōu)末制導(dǎo)設(shè)計［J］.兵工學(xué)報，2008，29（1）:63-67. HOU Ming-shan.Optimum terminal guidance for air-to-ground missile with impact angle constraint［J］.Acta Armamentarii，2008，29（1）:63-67.（in Chinese）

［15］ Zhang L M，Sun M W，Chen Z Q，et al.Receding horizon trajectory optimization with terminal impact specifications［J］.Mathematical Problems in Engineering，2014（5）:604705.

［16］ Rusnak I，Meir L.Modern guidance law for high-order autopilot［J］.Journal of Guidance，Control，and Dynamics，1991，14（5）:1056-1058.

［17］ 關(guān)冶，陸金甫.數(shù)值分析基礎(chǔ)［M］.北京:高等教育出版社，2005. GUAN Ye，LU Jin-pu.Fundamentals of numerical analysis［M］.Beijing:High Education Press，2005.（in Chinese）

［18］ Darby C L，Rao A V.Minimum-fuel low-earth-orbit aeroassisted orbital transfer of small spacecraft［J］.Journal of Spacecraft and Rockets，2011，48（4）:618-628.

［19］ Furfaro R，Cersosimo D，Wibben D R.Asteroid precision landing via multiple sliding surfaces guidance techniques［J］.Journal of Guidance，Control，and Dynamics，2013，36（4）:1075-1092.

［20］ Boyarko G，Romano M，Yakimenko O.Time-optimal reorientation of a spacecraft using an inverse dynamics optimization method［J］. Journal of Guidance，Control，and Dynamics，2011，34（4）: 1197-1208.

［21］ 張友安，黃詰，孫陽平.帶有落角約束的一般加權(quán)最優(yōu)制導(dǎo)律［J］.航空學(xué)報，2014，35（3）:848-856. ZHANG You-an，HUANG Jie，SUN Yang-ping.Generalized weighted optimal guidance laws with impact angle constraints［J］. Acta Aeronautica et Astronautica Sinica，2014，35（3）:848-856.（in Chinese）

［22］ Kim B S，Lee J G，Han H S.Biased PNC law for impact with angular constraint［J］.Aerospace and Electronic Systems，1998，34（1）:277-287.

Optimal Guidance Law with Impact Angle Constraints Based on Indirect Gauss Pseudospectral Method

CHEN Qi，WANG Zhong-yuan，CHANG Si-jiang
（School of Energy and Power Engineering，Nanjing University of Science and Technology，Nanjing 210094，Jiangsu，China）

A novel optimal guidance law is proposed for the terminal guidance with impact angle constraints by using the combination of the minimal principle and Gauss pseudospectral method.An impact angle coordinate system is defined with an coordinate axis in the direction of the desired impact angle，and the linear engagement kinematics is established using this coordinate system.The control system of missile is simplified into a first-order inertial system.The canonical equation is obtained via the minimal principle，and then translated into a set of algebraic equations by employing the Gauss pseudospectral method.According to the boundary conditions，an analytical solution is finally derived for the optimal guidance law with impact angle constraints without any integral process or solving the Riccati differential equation.Numerical simulations show that the proposed guidance law ensures the much fast convergence of impact angle to the reference line，and has smaller required terminal acceleration compared with other guidance laws.In addition，the proposed guidance law can easily tackle with the guidance problem with complex weighting matrices.

ordnance science and technology；impact angle constraint；optimal control；terminal guidance；minimal principle；indirect Gauss pseudospectral method

TJ765.3

A

1000-1093（2015）07-1203-10

10.3969/j.issn.1000-1093.2015.07.008

2014-09-09

國家自然科學(xué)基金項目（11402117）

陳琦（1989—），男，博士研究生。E-mail:qiychan@126.com；王中原（1958—），男，研究員，博士生導(dǎo)師。E-mail:zywang@njust.edu.cn