楊希祥，楊慧欣，王 鵬

偽譜法及其在飛行器軌跡優(yōu)化設(shè)計領(lǐng)域的應(yīng)用綜述*

楊希祥，楊慧欣，王 鵬

(國防科技大學(xué) 航天科學(xué)與工程學(xué)院， 湖南 長沙 410073)

采用偽譜法進行飛行器軌跡優(yōu)化設(shè)計是近年來的熱點研究方向，然而較全面地對各種方法進行綜合分析的文獻卻很少。在對國內(nèi)外相關(guān)文獻進行系統(tǒng)研究的基礎(chǔ)上，闡述了航空航天領(lǐng)域應(yīng)用較為廣泛的幾種偽譜法的基本原理；歸納了偽譜法將連續(xù)最優(yōu)控制問題轉(zhuǎn)化為非線性規(guī)劃問題的思路和具體步驟；總結(jié)了偽譜法在飛行器軌跡優(yōu)化設(shè)計領(lǐng)域的應(yīng)用情況；對偽譜法及其在飛行器軌跡優(yōu)化設(shè)計領(lǐng)域應(yīng)用的未來研究方向進行了分析。

偽譜法；飛行器；軌跡；優(yōu)化設(shè)計；綜述

(College of Aerospace Science and Engineering, National University of Defense Technology, Changsha 410073, China)

飛行器軌跡優(yōu)化問題一般為非線性、帶有狀態(tài)約束和控制約束的最優(yōu)控制問題。求解最優(yōu)控制問題的數(shù)值方法很多，一般分為間接法和直接法[1]。間接法基于Pontryagin極小值原理推導(dǎo)最優(yōu)控制的一階必要條件，將最優(yōu)控制問題轉(zhuǎn)換為Hamilton邊值問題，然后通過打靶法等數(shù)值方法求解。間接法在早期的飛行器軌跡優(yōu)化研究中應(yīng)用廣泛[2-3]。直接法無須求解最優(yōu)性必要條件，而是采用參數(shù)化方法將連續(xù)空間的最優(yōu)控制問題求解轉(zhuǎn)化為非線性規(guī)劃(NonLinear Programming，NLP) 問題，通過數(shù)值求解NLP問題獲得最優(yōu)解。根據(jù)參數(shù)化方法不同，直接法分為兩種基本類型：僅離散控制變量的打靶法與同時離散控制變量和狀態(tài)變量的配點法。由于直接法彌補了間接法的諸多不足，近30年來得到迅速發(fā)展，并被更為廣泛地應(yīng)用于飛行器軌跡優(yōu)化領(lǐng)域[1,4]。

近年來，一類直接配點法——偽譜法 (pseudospectral method)備受關(guān)注[5-6]。該方法采用全局插值多項式在一系列離散點上近似狀態(tài)變量和控制變量，通過引入類似有限差分矩陣的偽譜差分矩陣將微分方程約束轉(zhuǎn)換為代數(shù)約束。從插值多項式的角度，偽譜法可稱為全局方法，而用分段多項式近似的傳統(tǒng)配點方法可稱為局部方法。由于偽譜法的配點一般為正交多項式的根，因此又稱為正交配點(orthogonal collocation)法。與局部方法相比，偽譜法能夠以較小的計算代價獲得較高的求解精度，因此在航空航天領(lǐng)域得到廣泛應(yīng)用，包括空間站姿態(tài)控制[7]、小推力下航天器軌道轉(zhuǎn)移[8]、月球軟著陸軌跡優(yōu)化、運載火箭軌跡優(yōu)化[9-10]等，并已在實際飛行中得到應(yīng)用[7]。

本文對航空航天領(lǐng)域常見偽譜法及其在飛行器軌跡快速優(yōu)化領(lǐng)域的應(yīng)用情況進行綜述。國外對各種偽譜法進行綜合研究的文獻并不多見，最近的一篇由Fahroo和Ross發(fā)表于2008年[11]，而國內(nèi)尚未見公開發(fā)表的相關(guān)文獻。

1 飛行器軌跡優(yōu)化問題描述

對于運載火箭、導(dǎo)彈和可重復(fù)使用運載器等不同類型的飛行器，軌跡優(yōu)化問題的目標(biāo)函數(shù)、約束條件和控制變量等各有特點，但都可以描述為一般的最優(yōu)控制問題：尋找控制變量u(t)∈m，最小化Bolza型性能指標(biāo)函數(shù)[12]

(1)

其中，狀態(tài)變量x(t)∈n，初始時間t0和終端時間tf(自由或固定)滿足：

動力學(xué)微分方程約束

(2)

邊界條件

φ[x(t0),t0,x(tf),tf]=0

(3)

以及路徑約束

C[x(t),u(t),t]≤0

(4)

式(1)～(4)中，函數(shù)Φ,L,f,φ,C定義為

Φ :n××n×→

L :n×m×→

f :n×m×→n

φ :n××n×→q

C :n×m×→c

應(yīng)用偽譜法時，都將最優(yōu)控制問題的時間區(qū)間由t∈[t0,tf]轉(zhuǎn)換到τ∈[-1,1]，對時間t作變換

(5)

將式(5)應(yīng)用于最優(yōu)控制問題的表達式(1)～(4)， 最小化性能指標(biāo)

(6)

動力學(xué)微分方程約束

(7)

邊界條件

φ[x(-1),t0,x(1),tf]=0

(8)

路徑約束

C[x(τ),u(τ),τ]≤0

(9)

2 軌跡快速優(yōu)化領(lǐng)域常用偽譜法

偽譜法源于求解微分方程數(shù)值解的譜方法。常見的譜方法有三種：Galerkin法、Tau法和配點法，其中譜配點法又稱偽譜法[13]。由于偽譜法的配點一般為正交多項式(Chebyshev多項式或Legendre多項式)的根，因此又被稱為正交配點法。

飛行器軌跡快速優(yōu)化領(lǐng)域常見的偽譜法包括：Legendre偽譜法、Radau 偽譜法、Gauss 偽譜法以及Chebyshev 偽譜法。各種偽譜法的不同在于所選取的插值基函數(shù)和配點類型不同，如表1所示。

表1 各種偽譜法插值基函數(shù)及配點類型

2.1 Legendre偽譜法

τ0=-1,τ1,τ2,…,τN-1,τN=1

(10)

(11)

對于k=0,1,…,N，

(12)

為N階Lagrange多項式。由Lagrange多項式的性質(zhì)可得，在配點處，利用插值多項式得到的值和利用動力學(xué)微分方程計算所得到的值相等。在第i個配點處，狀態(tài)變量的導(dǎo)數(shù)為

(13)

{Dik}為(N+1)×(N+1)的差分矩陣，且有

(14)

這樣，在N+1個配點處，原來連續(xù)時間最優(yōu)控制問題的動力學(xué)微分方程約束轉(zhuǎn)換為代數(shù)約束，即

i=0,1,…,N

(15)

類似的，性能指標(biāo)函數(shù)轉(zhuǎn)換為

(16)

式中，wi為權(quán)重。

邊界條件約束轉(zhuǎn)換為

φ[x(-1),t0,x(1),tf]=0

(17)

路徑約束轉(zhuǎn)換為

C[x(τi),u(τi)]≤0i=0,1,…,N

(18)

2.2 Radau偽譜法

Radau偽譜法[16-17]的配點采用(-1,1]上的N個Legendre-Gauss-Radau (LGR)點，插值節(jié)點由N個配點與初始點τ0=-1構(gòu)成，即配點數(shù)比插值節(jié)點數(shù)少1個，這也是Radau偽譜法和Legendre偽譜法的不同之處。

狀態(tài)變量近似表達式為

(19)

這樣，在N個配點處，原來連續(xù)時間最優(yōu)控制問題的動力學(xué)微分方程約束轉(zhuǎn)換為代數(shù)約束，即

i=1,…,N

(20)

{Dik}為N×(N+1)的差分矩陣，且有

(21)

g(τk)=(1+τk)[LN(τk)-LN-1(τk)]

(22)

控制變量近似表示為

(23)

邊界控制點可采用樣條插值外推得到[18]。性能指標(biāo)函數(shù)轉(zhuǎn)換為

(24)

邊界條件約束轉(zhuǎn)換為

φ[x(-1),t0,x(1),tf]=0

(25)

路徑約束轉(zhuǎn)換為

C[x(τi),u(τi)]≤0i=1,…,N

(26)

2.3 Gauss偽譜法

Gauss偽譜法[13]的配點采用(-1,1)上的N-1個Legendre-Gauss (LG)點，插值節(jié)點由N-1個配點與初始點τ0=-1構(gòu)成，即配點數(shù)比插值節(jié)點數(shù)少1個。

狀態(tài)變量近似表達式為

(27)

這樣，在N-1個配點處，原來連續(xù)時間最優(yōu)控制問題的動力學(xué)微分方程約束轉(zhuǎn)換為代數(shù)約束，即

(28)

{Dik}為(N-1)×N的差分矩陣，且有

控制變量近似表示為

(29)

性能指標(biāo)函數(shù)轉(zhuǎn)換為

(30)

邊界條件約束轉(zhuǎn)換為

φ[x(-1),t0,x(1),tf]=0

(31)

路徑約束轉(zhuǎn)換為

C[x(τi),u(τi)]≤0i=1,…,N-1

(32)

式(28)未定義終端狀態(tài)xf，終端狀態(tài)也應(yīng)滿足動力學(xué)方程約束，將終端狀態(tài)約束條件離散并用Gauss積分來近似，可得

(33)

文獻[19]從近似精度和計算效率等方面對上述三種偽譜法進行了比較。Radau偽譜法和Gauss偽譜法在狀態(tài)變量、控制變量和共軛變量的近似精度上均優(yōu)于Legendre偽譜法，同時，Gauss偽譜法對共軛變量邊界值的估計精度高于Radau偽譜法，且在處理含初始和終端約束的問題上具有優(yōu)勢。在計算效率方面，求解相同規(guī)模問題時三種方法耗時差別不大。

2.4 Chebyshev偽譜法

和上述三種偽譜法的配點基于Legendre多項式不同，Chebyshev 偽譜法[20-22]的配點基于Chebyshev多項式，采用[-1,1]上的N+1個Chebyshev-Gauss-Lobatto (CLG)點，即

(34)

狀態(tài)變量和控制變量的近似表達式為

(35)

(36)

對于k=0,1,…,N，

(37)

(38)

這樣，在N+1個配點處，連續(xù)時間最優(yōu)控制問題的動力學(xué)微分方程約束轉(zhuǎn)換為代數(shù)約束，即

i=0,1,…,N

(39)

{Dik}為(N+1)×(N+1)的差分矩陣，且有

性能指標(biāo)函數(shù)轉(zhuǎn)換為

邊界條件約束轉(zhuǎn)換為

φ[x(-1),t0,x(1),tf]=0

(41)

路徑約束轉(zhuǎn)換為

C[x(τi),u(τi)]≤0i=0,1,…,N

(42)

3 軌跡快速優(yōu)化領(lǐng)域偽譜法應(yīng)用綜述

早在1822年，F(xiàn)ourier就已將譜理論引入微分方程的研究。利用譜方法求解微分方程數(shù)值解的思想，可以追溯到1938年Lanczos所做的工作，他建立了選擇譜配點法的基函數(shù)與配點的理論[23]。20世紀五六十年代，隨著Chebyshev多項式應(yīng)用于初值問題，譜方法再次復(fù)蘇。20世紀70年代，Gottlieb等首次建立了譜方法的統(tǒng)一數(shù)學(xué)理論；譜方法的成名源于其在流體動力學(xué)領(lǐng)域的應(yīng)用[24]；1972年，Orszag首次將譜配點法稱為偽譜法[25]。80年代后期，Vlassenbroeck等首次利用偽譜法求解最優(yōu)控制問題[20]，當(dāng)時采用的是Chebyshev偽譜法；緊接著，Elnagar等提出了Legendre偽譜法[14]；美國海軍研究院的Fahroo和Ross等對偽譜法的研究作出了重要貢獻，并將其引入飛行器軌跡優(yōu)化領(lǐng)域[9,22]。

2000年以后，尤其是近幾年來，伴隨著利用偽譜法求解最優(yōu)控制問題研究的逐漸增多[6,26]，其在航空航天領(lǐng)域的應(yīng)用研究開始受到重視，并由理論計算發(fā)展到地面試驗，進而發(fā)展到飛行試驗驗證[7]，而飛行器軌跡優(yōu)化是偽譜法應(yīng)用最廣泛、最成功的領(lǐng)域。Rao等采用Legendre偽譜法研究了高升阻比機動再入飛行器再入軌跡優(yōu)化問題，并對得到的最優(yōu)解進行了驗證和分析[27]；Bollino等將Legendre偽譜法應(yīng)用于無人飛行器軌跡在線快速生成[28]；Fahroo等研究了Legendre偽譜法在可重復(fù)使用運載器可達區(qū)快速生成問題中的應(yīng)用，研究結(jié)果表明，采用偽譜法可以避免傳統(tǒng)求解方法(如能態(tài)近似法) 的很多不足，提高求解速度和精度[29]；Williams研究了Legendre偽譜法在飛行器地形跟蹤軌跡實時優(yōu)化問題中的應(yīng)用[30]；文獻[31]和文獻[32]分別采用Gauss偽譜法實現(xiàn)飛機參數(shù)近實時辨識和飛機軌跡優(yōu)化; Ogawa等將進化算法與Gauss偽譜法相結(jié)合，解決火箭基組合動力兩級入軌可重復(fù)使用運載器軌跡優(yōu)化問題[33]。

2001年以來，麻省理工學(xué)院航空航天系對偽譜法進行了系統(tǒng)深入研究，取得了眾多研究成果，出版了多篇學(xué)位論文。Rea將Legendre偽譜法應(yīng)用于運載火箭上升段軌跡快速優(yōu)化問題，分別研究了二維平面內(nèi)的單級運載火箭軌跡優(yōu)化問題和三維空間兩級運載火箭軌跡優(yōu)化問題，并用一個預(yù)測制導(dǎo)算法研究了采用偽譜法實現(xiàn)在線控制指令生成的可行性，文章同時提出了一種引入二階差分矩陣減小問題求解規(guī)模的方法[34]；Clarke將Legendre偽譜法成功應(yīng)用于通用航空飛行器性能優(yōu)化[35]；Taylor研究了Legendre偽譜法在航天飛機應(yīng)急下降軌跡優(yōu)化問題中的應(yīng)用[36]；Stanton采用Legendre偽譜法求解航天器最優(yōu)轉(zhuǎn)移軌道問題，包括共面轉(zhuǎn)移問題、簡單軌道面改變問題和復(fù)雜的軌道尺寸/外形/軌道面改變問題，目標(biāo)函數(shù)均選擇為燃料消耗最省，仿真結(jié)果證明，求解精度和求解效率都能夠得到很好保證[37]；Benson較系統(tǒng)地研究了Gauss偽譜法，并將其應(yīng)用于小推力航天器軌跡優(yōu)化問題和Delta運載火箭上升段軌跡優(yōu)化問題[38]，研究結(jié)果表明，與Legendre偽譜法相比，Gauss偽譜法都具有更快的收斂速度和更精確的共軛變量估算結(jié)果；Hawkins研究了Legendre偽譜法在月球軟著陸軌跡快速優(yōu)化問題中的應(yīng)用[39]。

從2007年開始，有關(guān)偽譜法在航空航天領(lǐng)域中應(yīng)用的文章陸續(xù)出現(xiàn)在《宇航學(xué)報》等國內(nèi)中文核心期刊中。袁建平教授等將Legendre偽譜法應(yīng)用于月球軟著陸軌跡快速優(yōu)化問題[40]和空間飛行器有限推力軌道轉(zhuǎn)移優(yōu)化設(shè)計問題[41]，研究結(jié)果證明，算法具有收斂速度快、對初值不敏感、魯棒性強等優(yōu)點；曹喜濱教授則采用Gauss偽譜法解決日-火Halo轉(zhuǎn)移軌道快速優(yōu)化設(shè)計問題[42]；楊希祥等研究了Gauss偽譜法在多級固體運載火箭上升段軌跡快速優(yōu)化和空空導(dǎo)彈最優(yōu)中制導(dǎo)律設(shè)計中的應(yīng)用，提出采用串行策略進一步提高偽譜法計算效率[43-44]；雍恩米等將Gauss偽譜法應(yīng)用于高超聲速飛行器再入軌跡快速優(yōu)化問題，研究了采用偽譜法在線生成飛行器軌跡時的初值生成問題[45-46]；李海陽教授所在課題組對Gauss偽譜法在月球軟著陸軌道優(yōu)化中的應(yīng)用進行了大量研究[47]；此外，楊希祥及其所在課題組對Gauss偽譜法在超空泡水下航行體和平流層飛艇軌跡優(yōu)化中的應(yīng)用進行了研究[48-49]?？傮w來看，國內(nèi)在偽譜法理論研究方面還不夠系統(tǒng)深入，在將其應(yīng)用于飛行器軌跡快速優(yōu)化設(shè)計領(lǐng)域研究方面正處于快速發(fā)展和范圍延伸階段。

4 結(jié)論與展望

當(dāng)前，偽譜法在飛行器軌跡快速優(yōu)化領(lǐng)域的應(yīng)用正日益受到重視，根據(jù)對國內(nèi)外研究現(xiàn)狀的研究和本文作者們所在課題組的研究工作體會，今后一段時期內(nèi)，偽譜法及其在飛行器軌跡優(yōu)化領(lǐng)域應(yīng)用研究將重點集中在兩個方面，一是偽譜法理論的完善研究，二是針對飛行器軌跡快速優(yōu)化領(lǐng)域的現(xiàn)實需求，深化和拓展偽譜法的應(yīng)用。

4.1 偽譜法理論深化研究

盡管譜方法用于求解微分方程已有70余年的歷史，利用偽譜法求解最優(yōu)控制問題也已有20多年的歷史，但偽譜法理論還有很多方面亟待完善。

1)NLP問題的 Karush-Kuhn-Tucker (KKT)條件與Hamilton邊值問題一階最優(yōu)必要條件的等價性問題。此問題可轉(zhuǎn)化成NLP問題的最優(yōu)解是否是原最優(yōu)控制問題最優(yōu)解問題，這一問題也是直接法普遍面臨的關(guān)鍵問題，很多學(xué)者對這一問題的證明做了研究，例如，Benson證明了采用Gauss偽譜法轉(zhuǎn)化得到NLP問題的KKT條件與Hamilton邊值問題一階最優(yōu)必要條件是等價的[13]，同時，證明了節(jié)點上滿足路徑約束即可保證整個過程滿足路徑約束；Huntington在他的博士學(xué)位論文中深化了這一工作[50]，但對于Radau和Legendre等偽譜法，這一關(guān)鍵問題的證明尚未進行或完成。

2)配點的選取問題。配點是決定偽譜法求解精度和求解效率的關(guān)鍵因素之一，很多學(xué)者仍在對配點的選取問題進行積極研究，例如，文獻[51]提出了Legendre 偽譜法的改進形式Jacobi偽譜法，它使用Jacobi多項式尋找配點，Legendre 多項式是Jacobi多項式的子集；文獻[52]提出了Hermite-LGL方法，它不使用Legendre 多項式，而是使用分段三次多項式，在LGL點的子集進行配點。隨著偽譜法理論的不斷完善和數(shù)值計算方法的發(fā)展，這一方面的研究仍將受到關(guān)注。

3)偽譜法中共軛變量的估計問題。共軛變量的精確估計對驗證最優(yōu)控制問題解的最優(yōu)性、分析狀態(tài)變量對優(yōu)化性能指標(biāo)的敏感度和進行計算網(wǎng)格加細(mesh refinement)是十分重要的。間接法的求解表達式中都包含了共軛變量，因此其解也一般都給出共軛變量估計值。由于直接法不需要構(gòu)建Hamilton邊值問題方程，因此不能對共軛變量進行顯式近似，很多學(xué)者都一直致力于估計方法的研究[15,53-54]，雖然取得了一定成果，但仍未成熟。

4)計算網(wǎng)格加細問題。偽譜法處理最優(yōu)控制問題的高精度和高效率是在假定問題完全光滑的前提下得到的。但實際的復(fù)雜系統(tǒng)狀態(tài)變量往往是不連續(xù)的，或者控制變量是不連續(xù)的，或者兩者都是不連續(xù)的，而且間斷點的數(shù)目和位置往往是未知的。計算網(wǎng)格加細算法能夠檢測到間斷點的位置和數(shù)目，從而大大提高偽譜法解決問題的能力。

偽譜法理論中需要深入研究的問題還包括邊界控制量的估計、收斂性證明[55]等。

4.2 軌跡優(yōu)化領(lǐng)域偽譜法應(yīng)用深化研究

偽譜法在飛行器軌跡快速優(yōu)化領(lǐng)域的進一步研究，一是要拓寬算法的應(yīng)用對象和應(yīng)用范圍，二是要針對軌跡優(yōu)化中的具體問題應(yīng)用和改進偽譜法。

1)離線軌跡優(yōu)化問題。快速、機動發(fā)射是現(xiàn)代作戰(zhàn)飛行器發(fā)展的必然方向。對于采用攝動制導(dǎo)的快速機動發(fā)射運載火箭或彈道導(dǎo)彈，需要在射前預(yù)先生成一條最優(yōu)標(biāo)準彈道，從而完成射擊諸元裝訂。射擊諸元的快速裝訂對縮短發(fā)射準備時間具有重要作用，而發(fā)射準備時間關(guān)系著運載工具的射前生存能力和所運送有效載荷執(zhí)行任務(wù)的快速性，機動發(fā)射地點的時變性，對彈道快速優(yōu)化設(shè)計提出了新要求，偽譜法的高計算效率可應(yīng)對這一挑戰(zhàn)，它在這一領(lǐng)域的應(yīng)用值得關(guān)注。高效率和高精度是優(yōu)化計算不斷追求的目標(biāo)。偽譜法的節(jié)點個數(shù)是影響計算效率和精度的因素之一，將節(jié)點個數(shù)同時作為優(yōu)化問題的設(shè)計變量，提高優(yōu)化效率和精度，是偽譜法求解最優(yōu)控制問題的研究方向之一[46]。

2)在線軌跡優(yōu)化問題。軌跡在線生成(快速重規(guī)劃)是當(dāng)前飛行器軌跡優(yōu)化領(lǐng)域的熱門研究方向。傳統(tǒng)的飛行器制導(dǎo)算法是射前預(yù)先計算好一條最優(yōu)標(biāo)準彈道，利用一定的制導(dǎo)控制算法，使飛行器克服外界干擾，盡量沿設(shè)計的彈道飛行。這樣的飛行彈道實際上是次最優(yōu)彈道。軌跡在線生成技術(shù)的引入，使飛行器在因外界干擾作用而偏離設(shè)計的最優(yōu)彈道時，能夠快速規(guī)劃出一條新的最優(yōu)彈道。簡化的、能轉(zhuǎn)化為線性規(guī)劃的最優(yōu)控制問題，已經(jīng)實現(xiàn)了軌跡在線生成，對于只能轉(zhuǎn)化為NLP問題的復(fù)雜最優(yōu)控制問題，實現(xiàn)軌跡在線生成，還有很多問題需要研究[51]，偽譜法計算高精度和高效率特點，使其恰恰具有應(yīng)用于這一領(lǐng)域的潛力。很多軌跡在線生成問題可行域的非凸性，使得任意給定初值條件下，利用偽譜法很難得到最優(yōu)軌跡，因此初值生成方法成為亟待研究的問題，良好的初值生成方法，將大大減小偽譜法的計算時間[45]，從而提高軌跡在線生成速度。

在偽譜法研究方面，國外經(jīng)歷了由理論研究到開展地面試驗[56]，再到應(yīng)用于實際飛行的歷程[7,57]。在理論研究和數(shù)值仿真成熟的條件下，我們也可以有條件地適時設(shè)計開展地面試驗，以驗證理論研究的正確性，并逐步實現(xiàn)其在實際型號中的應(yīng)用。

References)

[1] Betts J T. Survey numerical methods for trajectory optimization[J].Journal of Guidance, Control and Dynamics, 1998, 21(2):193-207.

[2] Vinh N X.Optimal multiple-pass aeroassisted plane change[J]. Acta Astronautica, 1990, 21(11/12): 749-758.

[3] Vinh N X, Lu P. Necessary conditions for maximax problems with application to aeroglide of hypervelocity vehicles[J]. Acta Astronautica, 1987, 13(6/7): 413-420.

[4] 雍恩米, 陳磊, 唐國金. 飛行器軌跡優(yōu)化數(shù)值方法綜述[J].宇航學(xué)報, 2008, 29(2):397-406.YONG Enmi, CHEN Lei, TANG Guojin. A survey of numerical methods for trajectory optimization of spacecraft[J]. Journal of Astronautics, 2008, 29(2):397-406.(in Chinese)[5] Fahroo F, Ross I M. Convergence of pseudospectral discretizations of optimal control problems[C]// Proceedings of the 40th IEEE Conference on Decision and Control, 2001, 4: 3175-3177.

[6] Gong Q, Fahroo F, Ross I M. Spectral algorithm for pseudospectral methods in optimal control[J]. Journal of Guidance, Control and Dynamics, 2008, 31(3): 460-471.

[7] Kang W, Bedrossian N. Pseudospectral optimal control theory makes debut flight, saves NASA $1M in under three hours [J]. SIAM News, 2007, 40(7).[8] Riehl J P, Paris S W, Sjaw W K. Comparison of implicit integration methods for solving aerospace trajectory optimization problems[C]//Proceedings of AIAA/AAS Astodynamics Specialist Conference and Exhibit, 2006.

[9] Fahroo F, Ross I M, Souza C D, et al. A fast approach to multi-stage launch vehicle trajectory optimizaiton[C]//Proceedings of AIAA Guidance, Navigation, and Control Conference and Exhibit, 2003.

[10] Rea J. Launch vehicle trajectory optimizaiton using a Legendre pseudospectral method[C]//Proceedings of AIAA Guidance, Navigation, and Control Conference and Exhibit, 2003.

[11] Fahroo F, Ross I M. Advances in pseudospectral methods for optimal control[C] //Proceedings of AIAA Guidance, Navigation, and Control Conference and Exhibit, 2008.

[12] Pietz J A. Pseudospectral collocation methods for the direct transcription of optimal control problems[D]. USA: Rice University, 2003.

[13] Benson D. A gauss pseudospectral transcription for optimal control[D].USA: Massachusetts Institute of Technology, 2005.

[14] Elnagar G, Kazemi M A, Razzaghi M. The pseudospectral Legendre method for discretizintg optimal control problems[J].IEEE Transactions on Automatic Control, 1995, 40(10): 1793-1796.

[15] Fahroo F, Ross I M. Costate estimation by a Legendre pseudospectral method[J].Journal of Guidance, Control and Dynamics, 2001, 24(2): 270-277.

[16] Fahroo F, Ross I M. Pseudospectral method for infinite horizon nonlinear optimal control problems[C]//Proceedings of AIAA Guidance, Navigation, and Control Conference and Exhibit, 2005.

[17] Kameswaran S, Biegler L T. Convergence rates for direct transcription of optimal control problems at Radau points[J]. Computational Optimization and Applications, 2008, 41(1) : 81-126.

[18] Betts J T. Practical methods for optimal control using nonlinear programming[M].USA:Society for Industrial and Applied Mathematics Press, 2001.

[19] Fahroo F, Ross I M. On discrete-time optimality conditions for pseudospectral methods[C]//Proceedings of AIAA/ AAS Astrodynamics Specialist Conference and Exhibit, 2006.

[20] Vlassenbroeck J, Van Dooren R. A Chebyshev technique for solving nonlinear optimal control problems[J]. IEEE Transactions on Automatic Control, 1988, 33(4): 333-340.[21] Elnagar G N, Kazemi M A. Pseudospectral Chebyshev optimal control of constrained nonlinear dynamical systems[J].Computational Optimization and Applications, 1998, 11(2) : 195-217.

[22] Fahroo F, Ross I M. Direct trajectory optimization by a Chebyshev pseudospectral method[J]. Journal of Guidance, Control and Dynamics, 2002, 25(1):160-166.

[23] Lanczos C. Trigonometric interpolation of empirical and analytical functions [J]. Journal of Mathmatics and Physics, 1938(103): 123-199.

[24] Canuto C, Hussaini M Y,Quarteroni A, et al. Spectral methods in fluid dynamics[M].USA: Springer-Verlag, 1988.

[25] Orszag S A. Comparison of pseudospectral and spectral approximations[J]. Study of Application Mathmatics, 1972(51): 689-703.[26] Gong Q, Kang W, Ross I M. A pseudospectral methods for the optimal control of constrained feedback linearizable[J]. IEEE Transactions on Automatic Control, 2006, 51(7):1115-1129.

[27] Rao A V, Clarke K A. Performance optimization of a maneuvering re-entry vehicle using a Legendre pseudospectral method[C]//Proceedings of AIAA Atmospheric Flight Mechanics Conference and Exhibit, 2002.

[28] Bollino K P, Lewis L R, Sekhavat P,et al. Pseudospectral optimal control：a clear road for autonomous intelligent path planning[C]//Proceedings of AIAA Infotech and Aerospace Conference and Exhibit, 2007.

[29] Fahroo F, Doman D B, Ngo A D.Footprint generation for reusable launch vehicle using a direct pseudospectral method[C]//Proceedings of American Control Conference, 2003, 3: 2163-2168.

[30] Williams P. Three-dimensional aircraft terrain-following via real-time optimal control[J]. Journal of Guidance, Control and Dynamics, 2007, 30(4): 1201-1206.

[31] Jorris T R, Mccracken B. Aircraft system identification using pseudospectral parameter optimization wtih adaptive nodes[C]//Proceedings of AIAA Atmospheric Flight Mechanics Conference, 2011.

[32] Bittner M, Fisch F, Holzapfel F. A multi-model Gauss pseudospectral optimization method for aircraft trajecto-ries[C]//Proceedings of AIAA Atmospheric Flight Mechanics Conference, 2012.

[33] Ogawa H, Kodera M, Tomioka S. Multi-phase trajectory optimization for access-to-space with rbcc-powered tsto via surrogated-assisted hybrid evolutionary algorithms incorporating pseudo-spectral methods[C]//Proceedings of 19th AIAA International Space Planes and Hypersonic Systems and Technologies Conference, 2014.

[34] Rea J R. A Legendre pseudospectral method for rapid optimization of launch vehicle trajectories [D].USA: Massachusetts Institute of Technology, 2001.

[35] Clarke K A. Performance optimization study of common aero vehicle using a Legendre pseudospectral method [D].USA: Massachusetts Institute of Technology, 2003.

[36] Taylor C P. Optimization study of a trans-atlantic abort for the U.S. space shuttle using a Legendre pseudospectral method [D].USA: Massachusetts Institute of Technology, 2003. [37] Stanton S A. Optimal orbital transfer using a Legendre pseudospectral method[D].USA: Massachusetts Institute of Technology, 2003.

[38] Benson D. A Gauss pseudospectral transcription for optimal control[D].USA: Massachusetts Institute of Technology, 2005.

[39] Hawkins A M. Constrained trajectory optimization of a soft lunar landing from a parking orbit[D].USA: Massachusetts Institute of Technology, 2005.

[40] 羅建軍, 王明光, 袁建平. 基于偽光譜方法的月球軟著陸軌道快速優(yōu)化[J]. 宇航學(xué)報, 2007, 28(5) : 54-57.

LUO Jianjun, WANG Mingguang, YUAN Jianping.The Legendre pseudospectral method for rapid lunar soft-landing trajectory optimization for hypersonic reentry vehicle via Gauss pseudospectral method [J]. Journal of Astronautics, 2007, 28(5) : 54-57.(in Chinese)

[41] 涂良輝, 袁建平, 羅建軍. 基于偽光譜方法的有限推力軌道轉(zhuǎn)移優(yōu)化設(shè)計 [J]. 宇航學(xué)報, 2008, 29(4) :1189-1193.

TU Lianghui, YUAN Jianping, LUO Jianjun.Optimal design of orbital transfer with finite thrust based on Legendre pseudospectral method[J]. Journal of Astronautics, 2008, 29(4): 1189-1193.(in Chinese)

[42] 曹喜濱, 張相宇，王峰. 基于Gauss偽譜法的小推力日-火Halo軌道轉(zhuǎn)移優(yōu)化設(shè)計[J]. 宇航學(xué)報, 2013, 34(8): 1047-1054.

CAO Xibin, ZHANG Xiangyu, WANG Feng.Optimization of low-thrust transfer trajectory for the sun-mars Halo orbit based on Gauss pseudospectral method [J]. Journal of Astronautics, 2013, 34(8): 1047-1054.(in Chinese)

[43] 楊希祥, 張為華. 基于Gauss偽譜法的固體運載火箭上升段軌跡快速優(yōu)化研究[J]. 宇航學(xué)報, 2011, 32(1) :15-21.

YANG Xixiang, ZHANG Weihua.Rapid optimization of ascent trajectory for solid launch vehicle based on Gauss pseudospectral method [J]. Journal of Astronautics, 2011, 32(1) : 15-21.(in Chinese)

[44] 楊希祥, 張為華. 基于Gauss偽譜法的空空導(dǎo)彈最優(yōu)中制導(dǎo)律設(shè)計[J]. 國防科技大學(xué)學(xué)報, 2013, 35(1):28-32.

YANG Xixiang, ZHANG Weihua.Midcourse guidance law optimal design for air-to-air missiles based on Gauss pseudospectral method [J]. Journal of National University of Defense Technology, 2013, 35(1): 28-32.(in Chinese)

[45] 雍恩米, 唐國金, 陳磊. 基于Gauss偽譜方法的高超聲速飛行器再入軌跡快速優(yōu)化[J]. 宇航學(xué)報, 2008, 29(6) :1766-1772.

YONG Enmi, TANG Guojin, CHEN Lei.Rapid trajectory optimization for hypersonic reentry vehicle via Gauss pseudospectral method [J]. Journal of Astronautics, 2008, 29(6): 1766-1772.(in Chinese)

[46] 雍恩米.高超聲速滑翔式再入飛行器軌跡優(yōu)化與制導(dǎo)方法研究[D]. 長沙：國防科學(xué)技術(shù)大學(xué), 2008.

YONG Enmi. Study on trajectory optimization and guidance approach for hypersonic glide-reentry vehicle[D]. Changsha: National University of Defense Technology, 2008. (in Chinese)

[47] 彭祺擘, 李海陽, 沈紅新,等. 基于Gauss偽譜法和直接打靶法結(jié)合的月球定點著陸軌道優(yōu)化[J]. 國防科技大學(xué)學(xué)報, 2012, 34(2): 119-124.

PENG Qibo, LI Haiyang, SHEN Hongxin, et al. Lunar exact-landing trajectory optimization via the method combining gpm with direct shooting method[J]. Journal of National University of Defense Technology, 2012, 34(2): 119-124.(in Chinese)

[48] 林明東, 楊希祥, 張為華,等. 無尾舵布局超空泡航行體彈道優(yōu)化設(shè)計[J]. 國防科技大學(xué)學(xué)報, 2012, 34(6): 116-120.

LIN Mingdong, YANG Xixiang, ZHANG Weihua, et al. Trajectory optimization of planning supported supercavitating flight[J]. Journal of National University of Defense Technology, 2012, 34(6): 116-120.(in Chinese)

[49] 張禮學(xué)， 王中偉, 楊希祥, 等. 基于Gauss偽譜法的平流層飛艇上升段航跡規(guī)劃[J]. 上海交通大學(xué)學(xué)報, 2013, 47(8): 1205-1209.

ZHANG Lixue, WANG Zhongwei, YANG Xixiang,et al. Ascent trajectory planning for stratospheric airship based on Gauss pseudospectral method [J]. Journal of ShangHai JiaoTong University, 2013, 47(8): 1205-1209.(in Chinese)

[50] Huntington G T. Advancement and analysis of a Gauss pseudospectral transcription for optimal control problems[D].USA: Massachusetts Institute of Technology, 2007.

[51] Williams P. Jacobi pseudospectral method for solving optimal control problems[J]. Journal of Guidance, Control and Dynamics, 2004, 27(2):293-297.

[52] Williams P. Hermite-Legendre-Gauss-Lobatto direct transcription methods in trajectory optimization[J]. Journal of Guidance, Control, and Dynamics, 2009, 32(4):1392-1395.

[53] Herman A L, Conway B A. Direct optimization using collocation based on high order Gauss-Lobatto quadrature rules[J]. Journal of Guidance, Control and Dynamics, 1996, 19(3): 592-599.

[54] Seywald H, Kumar R R. Method for automatic costate calculation[J]. Journal of Guidance, Control and Dynamics, 1996, 19(6):1252-1261.

[55] 黃詰, 張友安, 王麗英. 基于Radau偽譜法的非線性最優(yōu)控制問題的收斂性[J].控制理論與應(yīng)用, 2014, 31(2):263-267.

HUANG Jie, ZHANG Youan, WANG Liying. Convergence of nonlinear optimal control problem using Radau pseudospectral method[J]. Journal of Astronautics, 2014, 31(2):263-267.(in Chinese)

[56] Ross I M, Sekhavat P, Fleming A, et al. Pseudospectral feedback control: foundations, examples and experimental results[C]///Proceedings of AIAA Guidance, Navigation, and Control Conference and Exhibit, 2006.

[57] Bedrossian N, Bhatt S, Lammers M, et al. First ever flight demonstration of zero propellant maneuver attitude control concept[C] //Proceedings of AIAA Guidance, Navigation, and Control Conference and Exhibit, 2007.

Overview of pseudo-spectral method and its application intrajectory optimum design for flight vehicles

YANG Xixiang， YANG Huixin， WANG Peng

Trajectory optimum design for flight vehicles using pseudo-spectral methods has been a hot research direction in recent years, but studies about the overview of this field are very few. Based on the systematic research of domestic and overseas related literatures, the basic principles of several pseudo-spectral methods which have a wide range of applications in aerospace field were expounded; the thoughts and formulas for transforming continuous optimal control problems into nonlinear programming problems were concluded; the applications of pseudo-spectral method in the trajectory optimum design for flight vehicles were summarized; the development directions of pseudo-spectral method and its application in trajectory optimization design for flight vehicles were discussed.

pseudo-spectral method; flight vehicles; trajectory; optimum design; overview

2015-04-01

國家自然科學(xué)基金資助項目(11102229)

楊希祥(1982—)，男，河北阜城人，副教授，博士，碩士生導(dǎo)師，E-mail：nkyangxixiang@163.com

10.11887/j.cn.201504001

http://journal.nudt.edu.cn

V221; V421.1

A

1001-2486(2015)04-001-08