• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    復雜約束條件下再入高超聲速滑翔飛行器軌跡快速優(yōu)化①

    2017-01-05 09:40:20邢清華毛藝帆
    固體火箭技術 2016年6期
    關鍵詞:偽譜禁飛區(qū)約束條件

    王 路,邢清華,毛藝帆

    (空軍工程大學 防空反導學院,西安 710051)

    復雜約束條件下再入高超聲速滑翔飛行器軌跡快速優(yōu)化①

    王 路,邢清華,毛藝帆

    (空軍工程大學 防空反導學院,西安 710051)

    作為直接法的典型代表,高斯偽譜方法在處理復雜約束條件(含路徑點或禁飛區(qū)約束)下再入高超聲速滑翔飛行器軌跡優(yōu)化問題時,僅能保證所得軌跡在各高斯節(jié)點處嚴格滿足各項約束,而節(jié)點之間軌跡的可行性無法保證,為解決這一問題,文章提出改進多階段高斯偽譜方法。該方法首先使用新定義的兩類節(jié)點(固定節(jié)點和虛擬節(jié)點)將軌跡分段,其中固定節(jié)點是用來保證路徑點與高斯節(jié)點重合,虛擬節(jié)點則是用來保證禁飛區(qū)附近分布更多的高斯節(jié)點,此分段方式能保證所得軌跡在任意位置可行;然后,向各分段軌跡插入指定數(shù)目的高斯節(jié)點;最后,使用序列二次規(guī)劃方法對此多階段軌跡優(yōu)化模型進行求解。仿真結果表明,改進多階段高斯偽譜方法在處理復雜約束條件下再入高超聲速滑翔飛行器軌跡優(yōu)化問題時快速有效。

    再入高超聲速滑翔飛行器;軌跡優(yōu)化;路徑點;禁飛區(qū);高斯偽譜方法

    0 引言

    高超聲速滑翔飛行器因其高升阻比外形設計具有良好的氣動性能,能憑借空氣動力的作用,實現(xiàn)遠距離無動力滑翔飛行。由于飛行器運動方程形式復雜、再入滑翔軌跡對控制量敏感,且飛行過程中受過載、熱流密度等約束影響較大,其軌跡優(yōu)化問題一直是研究的熱點。

    通常軌跡優(yōu)化的數(shù)值解法從本質(zhì)上可分為直接法和間接法[1]。間接法根據(jù)一階最優(yōu)必要條件和龐特里亞金最小化原理,將最優(yōu)化問題轉化為兩點邊值問題[2],缺點是對初值高度敏感,很難收斂,時效性不高。而直接法將原始最優(yōu)控制問題轉化為非線性規(guī)劃問題,在求解過程中具有較大的收斂半徑,且不受一階必要條件的限制[3]。

    作為直接法中的典型代表,高斯偽譜方法使用少量的節(jié)點,并具有較快的收斂速度,目前已被廣泛應用于普通約束條件下高超聲速飛行器上升段[4]或再入段[5-11]的軌跡優(yōu)化問題,這里的普通約束主要包括動壓、過載、熱流密度等過程約束、控制約束及終端狀態(tài)約束。然而,在飛行器實際再入過程中,不得不考慮某些真實約束,如路徑點和禁飛區(qū)約束。路徑點是指偵察或載荷投送時,飛行器必須經(jīng)過的空間位置,禁飛區(qū)是指受地理政治或敵方武器系統(tǒng)威脅,飛行器應嚴格避讓的某些特殊區(qū)域。

    高斯偽譜法將狀態(tài)變量和控制變量在一系列高斯節(jié)點上離散,所有約束條件在高斯節(jié)點上嚴格滿足,但高斯節(jié)點的分布呈中間疏松、兩端密集的特點。此時,考慮路徑點和禁飛區(qū)約束條件,可能會帶來以下問題:

    (1)相鄰高斯節(jié)點之間的軌跡不一定滿足禁飛區(qū)約束或過程約束;

    (2)路徑點約束條件不一定能夠滿足高精度。

    針對問題1和問題2,Jiang[12]提出了基于高斯偽譜方法的多階段軌跡優(yōu)化技術,通過所提出的3類插入點,將整個再入軌跡劃分為不同階段,由于其所設計的3類插入點都是固定的,所以最終得到的軌跡必須嚴格經(jīng)過各插入點。但通過研究發(fā)現(xiàn),在軌跡生成之前,只有路徑點是確定的,而禁飛區(qū)附近的飛行軌跡是未知的,所以只有路徑點可提前設置為固定插入點,Jorris[13-14]和Hu[15]所研究的軌跡優(yōu)化問題考慮了路徑點和禁飛區(qū)等復雜約束條件,但沒有具體討論高斯偽譜方法存在的上述兩類問題。針對以上問題,還有一類解決方法,就是使用大量的高斯節(jié)點對問題進行求解,但代價是十分耗時,實用價值不強。

    為克服傳統(tǒng)高斯偽譜方法在處理復雜約束條件下軌跡優(yōu)化問題時存在的不足,本文提出改進多階段高斯偽譜方法,通過使用新定義的兩類節(jié)點(固定節(jié)點和虛擬節(jié)點)將軌跡分段,在不增加總節(jié)點數(shù)目的同時,改變節(jié)點的分布。一方面,使路徑點與高斯節(jié)點重合,從而保證路徑點約束被嚴格滿足;另一方面,是禁飛區(qū)附近分布高密度的高斯節(jié)點,從而保證禁飛區(qū)約束滿足要求,又由于節(jié)點數(shù)目并沒有增加。所以,可保證所提算法的快速性。仿真結果驗證了其有效性。

    1 問題描述

    本節(jié)主要介紹高超聲速飛行器再入軌跡優(yōu)化問題模型,包括再入動態(tài)方程、約束條件和目標函數(shù)。

    1.1 再入動態(tài)方程

    在高超聲速飛行器的再入階段,假設目標為質(zhì)點,并將地球視為非旋轉勻質(zhì)球體,則飛行器再入運動方程如下[16]:

    (1)

    式中V為飛行器速度,m/s;θ為航跡傾角,rad;ψ為航跡偏角,rad;r為飛行器質(zhì)心至地心的距離,m;λ為地理經(jīng)度,rad;φ為地理緯度,rad;γc為滾轉角,rad;m為飛行器質(zhì)量,kg;D為阻力,D=CDρSV2/2;L為升力,L=CLρSV2/2;CD為阻力系數(shù);CL為升力系數(shù);S為氣動參考面積。

    在高超聲速條件下,升力、阻力系數(shù)CL、CD可近似表示為攻角α(°)與速度V(Ma)的函數(shù)[17]。

    (2)

    式中CLi和CDi(i=0,…,3)為常系數(shù)。

    1.2 約束條件

    1.2.1 過程約束

    高超聲速飛行器再入滑翔階段是一個高動態(tài)的過程,為滿足系統(tǒng)熱防護和結構穩(wěn)定性的需要,再入過程需滿足一系列苛刻約束條件的限制,如動壓限制、熱流密度限制、過載限制。

    (3)

    1.2.2 控制約束

    為了保證系統(tǒng)的穩(wěn)定性與可控性,攻角與滾轉角應滿足如下約束:

    (4)

    式中 下標“min”與“max”為對應變量的下限與上限。

    1.2.3 終端約束

    為使飛行器具備一定的末段打擊能力,對終端約束條件也應進行一定的限制。

    (5)

    式中 下標“f”為終端狀態(tài)。

    1.2.4 路徑點與禁飛區(qū)約束

    路徑點約束:為了完成諸如偵察或載荷投擲等特定任務,飛行器需精準的飛過指定路徑點上空,對應路徑點約束表示如下:

    (6)

    式中Npoint為路徑點的總個數(shù)。

    禁飛區(qū)約束:受地理政治或敵方武器系統(tǒng)威脅的限制,飛行器應嚴格避讓某些特殊區(qū)域,在此假設禁飛區(qū)為不限高度的圓柱體[13],第j個禁飛區(qū)的中心點坐標用(λj,φj)表示,對應的禁飛區(qū)半徑用Rj表示,則禁飛區(qū)約束表示為

    (7)

    式中Nzone為禁飛區(qū)的個數(shù)。

    1.3 目標函數(shù)

    通常情況下,不同的任務對應不同的目標函數(shù)。如果是為了維持飛行器較好的熱防護,目標函數(shù)可設置為最小熱流密度或最小峰值熱流密度;如果打擊或終點目標是固定的,目標函數(shù)可設置為達到時間最短;而如果是為了測試飛行器的性能,目標函數(shù)可設置為橫向或縱向航程最大等等。

    1.4 軌跡優(yōu)化

    對于高超聲速飛行器,其再入軌跡優(yōu)化問題就是尋找控制量α和γ,在滿足各項約束條件的同時,使得目標函數(shù)最優(yōu)。

    2 改進多階段高斯偽譜方法

    2.1 高斯偽譜法及其分析

    2.1.1 高斯偽譜法

    可用一般最優(yōu)控制問題描述高超聲速飛行器再入軌跡優(yōu)化問題,這里采用Bolza問題形式[18]進行表述。連續(xù)Bolza問題可描述為:求解控制變量u(t)∈Rm和狀態(tài)變量x(t)∈Rn,使如下目標函數(shù)[7]最小化。

    (8)

    同時滿足如下運動方程、過程約束以及終端約束:

    (9)

    式(8)、式(9)被稱為連續(xù)Bolza 問題的數(shù)學描述。高斯偽譜法在高斯節(jié)點上將狀態(tài)變量和控制變量進行離散,并利用這些離散點構造全局拉格朗日插值多項式來近似系統(tǒng)的動力學方程,從而將連續(xù)的無限維最優(yōu)控制問題轉換為非線性規(guī)劃問題進行求解[19]。對于高超聲速飛行器,飛行持續(xù)時間為從t0到tf,而高斯節(jié)點的分布區(qū)間為[-1,1]。因此,需將時間區(qū)域從t∈[t0,tf]轉化至τ∈[-1,1],轉化公式如式(10):

    (10)

    需要注意的是,N個離散的節(jié)點包含了初始節(jié)點τ0=-1、終端節(jié)點τf=1以及K=N-2個內(nèi)部節(jié)點,且這K個節(jié)點為K階拉格朗日插值多項式的根,則狀態(tài)和控制變量可近似表示如下:

    (11)

    其中,Li(τ)和Li*(τ)分別為N+1和N階拉格朗日插值多項式的基函數(shù),計算公式如下:

    (12)

    為獲得動態(tài)運動方程的約束,對式(11)求導可得

    (13)

    高斯偽譜方法中每個節(jié)點處目標狀態(tài)變量的倒數(shù)都應滿足式(13)。微分矩陣D∈RN×N+1可離線生成,其各元素計算公式如下[20]:

    式中PN(τ)為N階拉格朗日插值多項式;i=1,2,…,N;k=1,2,…,N。

    然后,將式(9)所述動態(tài)方程約束轉化為如下等式約束:

    其中,k=1,2,…,N。由于式(12)只表示高斯節(jié)點的目標狀態(tài),不包括初始時刻和終端時刻的目標狀態(tài),因此應增加2個端點處的等式約束如下[21]:

    (16)

    式中ωk為高斯積分的權重。

    然后,目標函數(shù)中的積分項通過高斯積分近似,則目標函數(shù)可表示為

    (17)

    同時,可將終端約束和過程約束表示為

    至此,將最優(yōu)控制問題轉化為非線性規(guī)劃問題,描述如下:求解非線性規(guī)劃問題的變量:

    (19)

    在滿足式(15)、式(16)和(式18)所述約束條件的同時,使得式(17)所述目標函數(shù)達到最小。

    2.1.2 分析

    如圖1所示,高斯節(jié)點的分布呈中間疏松兩端密集的特點。

    圖1 拉格朗日高斯節(jié)點分布示意圖Fig.1 Distribution of LG nodes

    高斯節(jié)點的這種分布特性會帶來以下幾方面的問題:

    (1)當高斯節(jié)點數(shù)目確定時,節(jié)點的分布位置在區(qū)間[-1,1]是確定的,如圖2所示,由于約束條件只在節(jié)點處嚴格滿足,所以很難保證所生成軌跡滿足路徑點約束,除非路徑點和某一高斯節(jié)點恰好重合。

    (2)如圖3所示,盡管禁飛區(qū)約束在所有高斯節(jié)點均嚴格滿足,但高斯節(jié)點之間的飛行器軌跡穿過了禁飛區(qū)。

    (3)過程約束條件不能保證嚴格滿足,如圖4所示。盡管某項過程約束指標在所有高斯節(jié)點均嚴格滿足,但該指標在多處超出了允許范圍,這會為飛行器維持其結構的穩(wěn)定性帶來嚴重挑戰(zhàn)。

    圖2 路徑點附近軌跡示意圖Fig.2 Example of trajectory near the waypoints

    圖3 禁飛區(qū)附近軌跡示意圖Fig.3 Example of trajectory near the no-fly zone

    圖4 過程約束變化示意圖Fig.4 Change flow of path constraint

    2.2 改進多階段高斯偽譜算法

    為克服前文所述傳統(tǒng)高斯偽譜方法在處理復雜約束條件下的軌跡優(yōu)化問題時存在的問題,本節(jié)提出改進多階段高斯偽譜算法。

    2.2.1 相關定義

    定義1 固定節(jié)點:只包含位置信息,不包含時間信息的節(jié)點,如路徑點。

    定義2 虛擬節(jié)點:只含時間信息,不包含位置信息的節(jié)點。

    2.2.2 算法流程

    (1)步驟1 使用固定節(jié)點將軌跡劃分為多階段

    如圖5所示,將路徑點設置為固定節(jié)點,使用固定節(jié)點將軌跡劃分為多階段,如此以來,路徑點約束可嚴格滿足。轉步驟2。

    圖5 軌跡分段示意圖(使用固定節(jié)點)Fig.5 Trajectory divided by fixed nodes

    (2)步驟2 模型求解

    假設軌跡被劃分為Z段,每階段的高斯節(jié)點數(shù)目為N,則第z階段的目標函數(shù)可表示為

    (20)

    則總目標函數(shù)可表示為

    (21)

    對應的終端和過程約束表示為

    (22)

    式(15)轉變?yōu)?/p>

    (23)

    式(16)轉變?yōu)?/p>

    (24)

    為保證各階段連接點處的連續(xù)性,增加如下約束:

    (25)

    至此,式(20)~式(25)將傳統(tǒng)高斯偽譜法轉化為多階段高斯偽譜法。模型的處理方式同傳統(tǒng)高斯偽譜法一致,本文選用適用于求解大規(guī)模非線性規(guī)劃問題的序列二次規(guī)劃方法。

    模型求解后,得到當前最優(yōu)軌跡,轉步驟3。

    (3)步驟3 軌跡可行性判別

    (26)

    經(jīng)檢測如果各約束條件在所有檢測點處嚴格滿足,則停止,所得軌跡即為最優(yōu)軌跡;否則,轉步驟4。

    (4)步驟4 增加虛擬節(jié)點將軌跡劃分為多階段

    步驟3中不滿足約束條件的檢測點由于只包含時間信息,因此將其視為虛擬節(jié)點,如圖7所示。軌跡被固定節(jié)點和新增虛擬節(jié)點劃分為多階段,由于高斯節(jié)點的分布特性,在虛擬節(jié)點的周圍會有較多的高斯節(jié)點分布,通過這種方式,可保證過程約束或路徑約束得到滿足。然后,轉步驟2。

    圖6 高斯節(jié)點和檢測點分布示意圖Fig.6 Distribution of LG nodes and checkpoints

    圖7 虛擬節(jié)點插入后禁飛區(qū)約束得到滿足示意圖Fig.7 No-fly zone constraint being satisfied by inserting virtual node

    2.2.3 初值生成算法

    本文使用傳統(tǒng)的初始化策略:先采用較少的高斯節(jié)點,求解出一條低精度的軌跡;然后,通過對此低精度軌跡進行插值,得到理想數(shù)目高斯節(jié)點所對應的各節(jié)點初始值。

    2.2.4 改進多階段高斯偽譜法分析

    改進多階段高斯偽譜法主要思想為:首先將路徑點視為固定節(jié)點,并用固定節(jié)點將軌跡分段,進行求解,從而使得到的軌跡嚴格滿足路徑點約束條件;為使全局軌跡嚴格滿足禁飛區(qū)和路經(jīng)約束,需要對所生成軌跡進行可行性驗證,若不滿足要求,則需增加虛擬節(jié)點,將軌跡重新分段,并進行求解,直至所生成軌跡嚴格滿足所有約束條件后即可。

    相比傳統(tǒng)高斯偽譜方法,改進多階段高斯偽譜方法在處理復雜約束條件下的軌跡優(yōu)化問題時,具有如下優(yōu)勢:

    (1)通過使用固定節(jié)點,將軌跡進行分段,所得軌跡嚴格滿足路徑點約束條件。

    (2)通過使用虛擬節(jié)點,將軌跡進行分段,所得軌跡不僅在高斯節(jié)點處嚴格滿足禁飛區(qū)以及過程約束,高斯節(jié)點之間的可行性同樣可保證。

    (3)固定節(jié)點和虛擬節(jié)點的加入,沒有改變高斯節(jié)點的數(shù)目,而是改變了高斯節(jié)點的分布。所以,算法的速度沒有受到影響。

    3 仿真及分析

    表1 初始時刻及終端時刻飛行器狀態(tài)Table 1 Initial and terminal condition of vehicle

    表2 路徑點與禁飛區(qū)參數(shù)設置Table 2 Parameters of waypoint and no-fly zone constraints

    3.1 粒子群算法仿真結果

    粒子的數(shù)目為20,其他參數(shù)設置與文獻[22]一致,所得最優(yōu)軌跡如圖8所示。

    (a) 三維軌跡 (b) 俯視軌跡 (c) 禁飛區(qū)1附近軌跡

    (d) 禁飛區(qū)2附近軌跡 (e) 路徑點附近軌跡圖8 粒子群算法仿真結果Fig.8 Results of optimal trajectory by PSO

    3.2 傳統(tǒng)高斯偽譜方法仿真結果

    高斯偽譜方法選用的高斯節(jié)點數(shù)目為30,仿真結果如圖9所示。

    (a) 三維軌跡

    (b) 俯視軌跡

    (d) 禁飛區(qū)2附近軌跡圖9 高斯偽譜方法仿真結果Fig.9 Results of optimal trajectory by GPM

    3.3 改進多階段高斯偽譜方法仿真結果

    為與傳統(tǒng)高斯偽譜方法進行比較,改進多階段高斯偽譜方法所選取的高斯節(jié)點數(shù)目同樣為30,仿真結果如圖10所示。

    (a) 三維軌跡

    (b) 俯視軌跡

    (c) 禁飛區(qū)1附近軌跡

    (d) 禁飛區(qū)2附近軌跡圖10 改進多階段高斯偽譜方法仿真結果Fig.10 Results of optimal trajectory by improved GPM

    3.4 仿真分析

    由圖8(a)、(b)可看出,所得軌跡較為平滑,圖8(c)、(d)說明禁飛區(qū)約束能夠嚴格滿足,圖8(e)說明路徑點約束不能滿足。

    圖9(a)給出了飛行器的三維運動軌跡。可見,所得軌跡較為平滑。圖9(b)~(d)顯示出了傳統(tǒng)高斯偽譜方法的不足,即只能保證所得軌跡在高斯節(jié)點處嚴格滿足約束條件,而無法保證其他部分的軌跡是否滿足禁飛區(qū)等復雜約束。

    由圖10給出的仿真結果可見,通過使用固定節(jié)點和虛擬節(jié)點將軌跡分為了4段,而總的高斯節(jié)點的數(shù)目沒有增加,只是改變了高斯節(jié)點的分布,所得軌跡嚴格滿足路徑點與禁飛區(qū)約束。

    此外,粒子群算法求得最終飛行器飛行時間為1 592 s,傳統(tǒng)高斯偽譜方法和改進高斯偽譜方法的最終結果分別為1 540 s和1 563 s。由此可見,高斯偽譜方法和改進高斯偽譜方法的精度都優(yōu)于粒子群算法,相比改進高斯偽譜方法,傳統(tǒng)高斯偽譜法雖然求得的目標函數(shù)值更小,但由于所得軌跡穿越了禁飛區(qū),因此軌跡無法使用。

    最后,分析不同仿真方法的耗時。粒子群算法的平均耗時約為2 min,但所得軌跡也只是接近全局最優(yōu);高斯偽譜方法的耗時約為1 min,但30個高斯節(jié)點不能保證節(jié)點之間的軌跡滿足約束條件,如果將高斯節(jié)點的數(shù)目增加,雖然所得軌跡精度會提高,但仿真耗時同樣會提高;相比之下,同樣是30個高斯節(jié)點,但本文所提算法的耗時為30 s。這是由于固定節(jié)點的插入將路徑點約束轉化成了終端約束,而虛擬節(jié)點的插入將禁飛區(qū)約束只在某分段模型中進行求解,而傳統(tǒng)高斯偽譜方法是將路徑點和禁飛區(qū)約束在整個模型中進行求解。因此,本文所提算法耗時較小,且精度更高。

    從以上仿真結果,可得出以下結論:

    (1)使用固定節(jié)點將軌跡進行分段,所生成軌跡能夠嚴格滿足路徑點約束;

    (2)使用虛擬節(jié)點將軌跡進行分段,會有更多的高斯節(jié)點分布于虛擬節(jié)點附近,從而實現(xiàn)所生成軌跡能嚴格滿足禁飛區(qū)約束和過程約束;

    (3)所提出的改進多階段高斯偽譜算法只是改變了高斯節(jié)點的分布,而沒有增加高斯節(jié)點的數(shù)目。因此,相比其他類型多階段偽譜方法,本文所提算法復雜度明顯降低。

    4 結論

    為解決復雜約束條件下再入高超聲速滑翔飛行器軌跡優(yōu)化問題,本文提出改進多階段高斯偽譜方法。通過使用固定節(jié)點和虛擬節(jié)點,將軌跡進行分段,在不增加高斯節(jié)點數(shù)目的前提下,改變了高斯節(jié)點的分布,以實現(xiàn)路徑點和禁飛區(qū)附近有較多的高斯節(jié)點分布。仿真結果表明,所得軌跡具有較高的精度,各項約束條件能夠嚴格滿足。

    [1] Huntington G T, Benson D, Rao A V. A comparison of accuracy and computational efficiency of three pseudospectral methods[C]//Proc. of the AIAA Guidance, Navigation, and Control Conference and Exhibit.2007, AIAA 20-23.

    [2] Guo Xiao, Zhu Ming. Direct trajectory optimization based on a mapped Chebyshev pseudospectral method[J]. Chinese Journal of Aeronautics, 2013, 26(2):401-412.

    [3] Ross I M, Fahroo F B. A perspective on methods for trajectory optimization[R].AIAA 2002-4727.

    [4] Bailing Tian. 3DOF ascent phase trajectory optimization for aircraft based on adaptive gauss pseudospectral method[C]//2012 Third International Conference on Intelligent Control and Information Processing, 2012: 431-453.

    [5] 孫勇,張卯瑞,梁曉玲. 求解含負責約束非線性最優(yōu)控制問題的改進Gauss偽譜法[J]. 自動化學報,2013,39(5):672-678.

    [6] 陳琦,王中原,常思江. 基于Gauss偽譜法的滑翔彈道快速優(yōu)化[J]. 彈道學報,2013,21(15):105-108.

    [7] 李柯,聶萬勝,馮必鳴. 助推-滑翔飛行器彈道分段優(yōu)化研究[J]. 指揮控制與仿真,2012,34(5):21-25.

    [8] 湯亮,楊建民,陳風雨. 多約束條件下升力滑翔式再入軌跡優(yōu)化[J].導彈與航天運載技術,2013(1):1-5.

    [9] 曹紅錦,葛致磊. 滑翔增程彈彈道優(yōu)化設計研究[J].四川兵工學報,2013,34(11):8-10.

    [10] 張鼎逆,劉毅. 基于改進遺傳算法和序列二次 規(guī)劃的再入軌跡優(yōu)化[J]. 浙江大學學報,2014,48(1):161-167.

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

    [12] Jiang Z, Rui Z. Reentry trajectory optimization for hypersonic vehicle satisfying complex constraints[J]. Chinese Journal of Aeronautics, 2013, 26(6): 1544-1553.

    [13] Jorris T R, Cobb R G. Three-dimensional trajectory optimization satisfying waypoint and no-fly zone constraints[J]. Journal of Guidance, Control, and Dynamics, 2009, 32(2): 551-572.

    [14] Jorris T R. Common aero vehicle autonomous reentry trajectory optimization satisfying waypoint and no-fly zone constraints[D]. Ohio, USA: Air Force Institute of Technology, 2007.

    [15] Hu Z D. Research on trajectory planning and guidance for space-based strike weapon[D]. Changsha, China: Graduate School of National University of Defense Technology, 2009.

    [16] Vinh N X. Optimal trajectories in atmospheric flight[D]. New York:Elsevier Scientific Publishing Company,1981.

    [17] 孫勇,段廣仁,張卯瑞,等. 高超聲速飛行器再入過程改進氣動系數(shù)模型[J]. 系統(tǒng)工程與電子技術,2011,33(1):134-137.

    [18] Benson D A, Huntington G T, Thorvaldsen T P, et al. Direct trajectory optimization and costate estimation via an orthogonal collocation method[J]. Journal of Guidance, Control, and Dynamics, 2006,29(6): 1435-1440.

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

    [20] Liu Y, Qian Y J, Li J Q, et al. Mars exploring trajectory optimization using gauss pseudo-spectral method[C]//Proc. of the International Conference on Mechatronics and Automation.2012: 2371-2377.

    [21] Yong E M, Chen L, Tang G J. Trajectory optimization of hypersonic gliding reentry vehicle based on the physical programming[J]. Acta Aeronautica et Astronautica Sinica, 2008, 29(5):1092-1097.

    [22] 謝富強,吳浩,唐靈靈. 基于粒子群算法的飛行器再入軌跡優(yōu)化[J].計算技術與自動化,2008,27(4):72-75.

    (編輯:薛永利)

    Trajectory rapid optimization for reentry hypersonic glide vehicle satisfying complex constraints

    WANG Lu,XING Qing-hua,MAO Yi-fan

    (School of Air and Missile Defense, Air Force Engineering University, Xi'an 710051, China)

    As a typical representation of direct method, in dealing with reentry hypersonic glide vehicle trajectory optimization problem under complex constraints including waypoint or no-fly zone constraints, gauss pseudo-spectral method (GPM) can merely ensure that all the constraints can be satisfied strictly at each Legendre-Gauss (LG) node, but cannot ensure the feasibility of the trajectory between each two neighbor nodes. To overcome this problem, an improved multi-segment gauss pseudo-spectral method was proposed. In this method, first of all, the trajectory was divided into multiple segments by two kinds of proposed nodes, fixed node and virtual node, the first one is used to ensure that each waypoint locates at one of LG nodes, and the second one is used to ensure that there will be more LG nodes locating near the no-fly zones, such that all constraints for whole trajectory can be satisfied; then the number of LG nodes for each segment was given; finally, the sequential quadratic programming method was used to solve the model. Simulation results demonstrate that the proposed optimization method can rapidly generate a reentry trajectory satisfying complex constraints.

    reentry hypersonic glide vehicle;trajectory optimization;waypoint;no-fly zone;gauss pseudospectral method

    2015-07-18;

    2015-09-02。

    國家自然科學基金(61272011);全軍軍事類研究生資助課題(2014JY526)。

    王路(1987—),男,博士生,研究方向為地基反臨近空間飛行器指揮決策問題。E-mail:408191081@qq.com

    V448.2

    A

    1006-2793(2016)06-0839-08

    10.7673/j.issn.1006-2793.2016.06.018

    猜你喜歡
    偽譜禁飛區(qū)約束條件
    基于一種改進AZSVPWM的滿調(diào)制度死區(qū)約束條件分析
    矩陣偽譜的新定位集及其在土壤生態(tài)系統(tǒng)的應用
    A literature review of research exploring the experiences of overseas nurses in the United Kingdom (2002–2017)
    大疆更新多邊形禁飛區(qū)策略
    航空模型(2017年7期)2017-08-15 10:24:44
    紊流環(huán)境下四維軌跡優(yōu)化的偽譜方法研究
    飛行力學(2017年1期)2017-02-15 03:00:58
    線性規(guī)劃的八大妙用
    偽譜法及其在飛行器軌跡優(yōu)化設計領域的應用綜述*
    基于Radau偽譜方法的軌跡優(yōu)化
    人臉禁飛區(qū)
    自然與科技(2013年6期)2013-04-29 00:44:03
    一類禁飛區(qū)后方安全撤離軌跡的設計方法研究
    99热全是精品| 日韩av不卡免费在线播放| 成人午夜精彩视频在线观看| 国产日韩一区二区三区精品不卡 | 国产成人freesex在线| 纯流量卡能插随身wifi吗| 好男人视频免费观看在线| 国产伦精品一区二区三区视频9| 成人免费观看视频高清| 自线自在国产av| 婷婷色综合www| 国产一区二区在线观看av| 日本-黄色视频高清免费观看| 久久婷婷青草| 国产精品国产三级国产专区5o| 国产日韩欧美视频二区| 嘟嘟电影网在线观看| 特大巨黑吊av在线直播| 欧美日韩av久久| 久久久久精品性色| 美女大奶头黄色视频| 亚洲国产精品999| 99热这里只有是精品在线观看| 久久97久久精品| 国产一区二区在线观看av| 人妻夜夜爽99麻豆av| 97超碰精品成人国产| 特大巨黑吊av在线直播| 最近最新中文字幕免费大全7| 五月开心婷婷网| 妹子高潮喷水视频| 一区二区三区免费毛片| 国产国拍精品亚洲av在线观看| 国产精品国产三级专区第一集| 国产一区二区三区综合在线观看 | 欧美激情国产日韩精品一区| 丝袜在线中文字幕| 91久久精品国产一区二区成人| 亚洲av成人精品一二三区| 永久网站在线| 中国美白少妇内射xxxbb| 亚洲国产欧美日韩在线播放| 99九九在线精品视频| 国产精品免费大片| 久久99蜜桃精品久久| 色94色欧美一区二区| av在线老鸭窝| 午夜激情av网站| 9色porny在线观看| 男的添女的下面高潮视频| 国产精品人妻久久久久久| 国产爽快片一区二区三区| 亚洲精品日本国产第一区| 青春草亚洲视频在线观看| 最新的欧美精品一区二区| 色94色欧美一区二区| 精品人妻在线不人妻| 久久久久久人妻| 国产 精品1| 午夜福利视频精品| 亚洲欧洲日产国产| 看非洲黑人一级黄片| 交换朋友夫妻互换小说| 人人妻人人爽人人添夜夜欢视频| 高清在线视频一区二区三区| 国产成人aa在线观看| 黄色毛片三级朝国网站| 又粗又硬又长又爽又黄的视频| 国产午夜精品久久久久久一区二区三区| 精品酒店卫生间| 久久精品熟女亚洲av麻豆精品| 嘟嘟电影网在线观看| 中文乱码字字幕精品一区二区三区| 中国国产av一级| 亚洲伊人久久精品综合| 午夜日本视频在线| 全区人妻精品视频| 91精品三级在线观看| 午夜福利视频在线观看免费| 在线观看www视频免费| 国内精品宾馆在线| 欧美日韩视频高清一区二区三区二| 日韩大片免费观看网站| 国产一区二区三区av在线| 亚洲国产精品一区二区三区在线| 99国产精品免费福利视频| 免费少妇av软件| 熟妇人妻不卡中文字幕| 免费黄色在线免费观看| 热re99久久国产66热| 多毛熟女@视频| 亚洲伊人久久精品综合| 国产一区二区三区av在线| 边亲边吃奶的免费视频| 国产精品.久久久| 少妇被粗大的猛进出69影院 | 黄片播放在线免费| 中文精品一卡2卡3卡4更新| 精品久久久久久久久亚洲| 欧美精品人与动牲交sv欧美| .国产精品久久| 亚洲国产精品999| 男女无遮挡免费网站观看| av线在线观看网站| 免费看av在线观看网站| 丝袜脚勾引网站| 一级毛片我不卡| 一边亲一边摸免费视频| 久久99精品国语久久久| 亚洲精品aⅴ在线观看| 看免费成人av毛片| 亚洲人成网站在线观看播放| av在线老鸭窝| av线在线观看网站| 如何舔出高潮| 久久精品夜色国产| 我要看黄色一级片免费的| 中国国产av一级| 精品久久久久久久久av| 成人国产av品久久久| 天天操日日干夜夜撸| 大香蕉久久网| 日本免费在线观看一区| 成人毛片a级毛片在线播放| 人妻夜夜爽99麻豆av| 在线天堂最新版资源| 久久精品国产a三级三级三级| 国产精品一区www在线观看| 亚洲精品国产av成人精品| 国产av国产精品国产| 国产男女内射视频| 在线观看国产h片| 一本一本综合久久| 日本爱情动作片www.在线观看| 黑人猛操日本美女一级片| 日本vs欧美在线观看视频| 黄片无遮挡物在线观看| 国产成人精品久久久久久| 老司机影院毛片| 岛国毛片在线播放| 欧美精品高潮呻吟av久久| 精品99又大又爽又粗少妇毛片| 青春草国产在线视频| 免费黄色在线免费观看| 国产成人aa在线观看| 在现免费观看毛片| 国产精品久久久久久av不卡| 欧美日韩视频高清一区二区三区二| 免费av中文字幕在线| 中文精品一卡2卡3卡4更新| 一本—道久久a久久精品蜜桃钙片| 婷婷色综合大香蕉| 成人手机av| 91久久精品电影网| 成人毛片60女人毛片免费| 国产熟女欧美一区二区| 国产成人一区二区在线| 国产成人午夜福利电影在线观看| 久久婷婷青草| 日韩电影二区| 国产色婷婷99| 中文字幕久久专区| 日韩电影二区| 久久婷婷青草| 黄色配什么色好看| 久久99一区二区三区| 日产精品乱码卡一卡2卡三| 多毛熟女@视频| 日本欧美国产在线视频| 日本91视频免费播放| 国产白丝娇喘喷水9色精品| 18在线观看网站| 汤姆久久久久久久影院中文字幕| 日韩av在线免费看完整版不卡| 国产永久视频网站| 成人毛片a级毛片在线播放| av网站免费在线观看视频| 天天操日日干夜夜撸| 中文天堂在线官网| 久久久亚洲精品成人影院| 视频在线观看一区二区三区| 亚洲精品视频女| 日韩人妻高清精品专区| 男女啪啪激烈高潮av片| 久久久亚洲精品成人影院| 日韩av在线免费看完整版不卡| 亚洲精品亚洲一区二区| 日韩av在线免费看完整版不卡| 婷婷色麻豆天堂久久| 久久99精品国语久久久| 欧美日韩在线观看h| 精品久久久久久久久亚洲| 免费播放大片免费观看视频在线观看| 日韩伦理黄色片| 九九久久精品国产亚洲av麻豆| 极品人妻少妇av视频| 国产亚洲午夜精品一区二区久久| 大片免费播放器 马上看| 亚洲高清免费不卡视频| 一本久久精品| 国产精品熟女久久久久浪| 精品人妻熟女av久视频| 午夜激情福利司机影院| 在线 av 中文字幕| 少妇 在线观看| 青春草国产在线视频| 18禁动态无遮挡网站| 亚洲图色成人| 亚洲欧美色中文字幕在线| 国产淫语在线视频| av视频免费观看在线观看| 色5月婷婷丁香| 五月开心婷婷网| 国产有黄有色有爽视频| 国产免费福利视频在线观看| 男的添女的下面高潮视频| 欧美亚洲 丝袜 人妻 在线| 亚洲美女黄色视频免费看| 99热国产这里只有精品6| 亚洲伊人久久精品综合| av又黄又爽大尺度在线免费看| 少妇猛男粗大的猛烈进出视频| 日韩成人av中文字幕在线观看| 亚洲,一卡二卡三卡| 精品国产国语对白av| 美女cb高潮喷水在线观看| 久久久久网色| 国产日韩欧美视频二区| 欧美 日韩 精品 国产| 涩涩av久久男人的天堂| 日韩精品免费视频一区二区三区 | 国产黄片视频在线免费观看| 丝袜美足系列| 美女国产高潮福利片在线看| 久久精品久久久久久久性| 中文字幕免费在线视频6| 乱码一卡2卡4卡精品| 人妻 亚洲 视频| 成人午夜精彩视频在线观看| 国产探花极品一区二区| 少妇高潮的动态图| 你懂的网址亚洲精品在线观看| 日韩一本色道免费dvd| 免费观看a级毛片全部| 国产成人91sexporn| 亚洲欧洲日产国产| 亚洲色图综合在线观看| 一区二区三区乱码不卡18| 欧美 亚洲 国产 日韩一| 麻豆乱淫一区二区| 99热这里只有精品一区| 日本欧美视频一区| 高清视频免费观看一区二区| 国产精品一区二区在线观看99| 在线看a的网站| 久久99精品国语久久久| 精品久久久久久电影网| 夫妻午夜视频| 美女福利国产在线| 国产 精品1| 乱码一卡2卡4卡精品| 热99久久久久精品小说推荐| 精品人妻一区二区三区麻豆| 国产一级毛片在线| .国产精品久久| 人妻夜夜爽99麻豆av| 美女主播在线视频| 爱豆传媒免费全集在线观看| 午夜久久久在线观看| 成人影院久久| 久久99热这里只频精品6学生| 久久 成人 亚洲| 国产 精品1| 九九爱精品视频在线观看| 男女无遮挡免费网站观看| 国产探花极品一区二区| 亚洲国产精品一区二区三区在线| 国产精品无大码| 自拍欧美九色日韩亚洲蝌蚪91| 精品少妇久久久久久888优播| 国产精品熟女久久久久浪| 欧美成人精品欧美一级黄| 91成人精品电影| 在线精品无人区一区二区三| 日本黄大片高清| 最黄视频免费看| 一级a做视频免费观看| av卡一久久| 日韩一区二区三区影片| 亚洲少妇的诱惑av| 十八禁网站网址无遮挡| 久久久亚洲精品成人影院| 久久精品久久久久久噜噜老黄| av又黄又爽大尺度在线免费看| 一级a做视频免费观看| 王馨瑶露胸无遮挡在线观看| 久久久久国产精品人妻一区二区| 伊人久久国产一区二区| 日韩成人伦理影院| 亚洲国产精品一区二区三区在线| 婷婷色综合大香蕉| 一级毛片电影观看| 日韩欧美精品免费久久| 韩国高清视频一区二区三区| 亚洲欧美色中文字幕在线| 精品酒店卫生间| 2022亚洲国产成人精品| 国产欧美另类精品又又久久亚洲欧美| 久久99热这里只频精品6学生| 人妻少妇偷人精品九色| 91精品国产九色| 日韩视频在线欧美| 亚洲成色77777| 欧美日韩成人在线一区二区| 欧美精品国产亚洲| 99久久综合免费| 超碰97精品在线观看| 亚洲av在线观看美女高潮| 亚洲性久久影院| 亚洲国产精品一区二区三区在线| 久久久久精品久久久久真实原创| 国产男女内射视频| 一本色道久久久久久精品综合| 国产亚洲欧美精品永久| 我的女老师完整版在线观看| 日本av手机在线免费观看| 51国产日韩欧美| 丝袜喷水一区| 色94色欧美一区二区| 午夜av观看不卡| 国产亚洲午夜精品一区二区久久| 夜夜骑夜夜射夜夜干| 亚洲av欧美aⅴ国产| 一区在线观看完整版| 中文字幕最新亚洲高清| 在现免费观看毛片| 伊人久久国产一区二区| 黄色配什么色好看| 黄片无遮挡物在线观看| 男女国产视频网站| 精品一区二区三区视频在线| 97在线人人人人妻| 国产一区二区三区综合在线观看 | av.在线天堂| 国产一区二区在线观看av| 亚洲精品久久午夜乱码| 三级国产精品片| 一区二区日韩欧美中文字幕 | 精品人妻一区二区三区麻豆| 黄色配什么色好看| 亚洲美女黄色视频免费看| 国产视频首页在线观看| 熟女人妻精品中文字幕| 夜夜爽夜夜爽视频| 免费观看在线日韩| 成人国产av品久久久| 熟女av电影| 亚洲精品日韩在线中文字幕| 午夜久久久在线观看| 成年人午夜在线观看视频| 一个人免费看片子| 欧美日韩综合久久久久久| 亚洲欧洲国产日韩| av在线播放精品| 如何舔出高潮| 久久av网站| 最新的欧美精品一区二区| 国国产精品蜜臀av免费| 中文乱码字字幕精品一区二区三区| 亚洲国产毛片av蜜桃av| 色网站视频免费| 99热网站在线观看| 一区二区三区免费毛片| 国产深夜福利视频在线观看| 国产精品.久久久| 精品久久久久久久久亚洲| 青春草亚洲视频在线观看| 色婷婷久久久亚洲欧美| 久久久a久久爽久久v久久| 亚洲成人av在线免费| 精品久久久久久电影网| 亚洲av电影在线观看一区二区三区| 亚洲精品第二区| 久久婷婷青草| 久久热精品热| 高清av免费在线| 精品人妻熟女av久视频| 人妻 亚洲 视频| 久久免费观看电影| 一级毛片电影观看| 国产精品99久久99久久久不卡 | 久久久久久久亚洲中文字幕| 久久久久久久久久人人人人人人| 成人国产麻豆网| 99九九在线精品视频| 黑人巨大精品欧美一区二区蜜桃 | 人妻一区二区av| 美女国产视频在线观看| 水蜜桃什么品种好| 波野结衣二区三区在线| 一边亲一边摸免费视频| 在线观看www视频免费| 亚洲久久久国产精品| 精品国产乱码久久久久久小说| 中文精品一卡2卡3卡4更新| 国语对白做爰xxxⅹ性视频网站| 乱人伦中国视频| 美女脱内裤让男人舔精品视频| 69精品国产乱码久久久| 午夜福利视频精品| 久久久久久久久久久丰满| 久久国产精品大桥未久av| 极品人妻少妇av视频| 亚洲欧美清纯卡通| 丝袜喷水一区| 亚洲国产精品一区二区三区在线| 啦啦啦在线观看免费高清www| 精品一区二区三区视频在线| 久久久精品免费免费高清| 男的添女的下面高潮视频| 欧美国产精品一级二级三级| 亚洲精品第二区| 人人妻人人爽人人添夜夜欢视频| 欧美日韩在线观看h| 在线观看免费高清a一片| 亚洲国产精品国产精品| 99视频精品全部免费 在线| 日韩中文字幕视频在线看片| 午夜激情av网站| 亚洲精华国产精华液的使用体验| 2022亚洲国产成人精品| 欧美日韩视频高清一区二区三区二| 国产一级毛片在线| 成人免费观看视频高清| 色吧在线观看| 精品一区二区三区视频在线| 久久久精品免费免费高清| 免费高清在线观看视频在线观看| 亚洲国产精品一区三区| 女的被弄到高潮叫床怎么办| 亚洲av日韩在线播放| 在线观看国产h片| 国产精品熟女久久久久浪| 人妻制服诱惑在线中文字幕| 亚洲精品久久午夜乱码| 欧美激情 高清一区二区三区| 日韩av免费高清视频| 午夜视频国产福利| 色94色欧美一区二区| 国产一区亚洲一区在线观看| 一个人免费看片子| av免费在线看不卡| 亚洲不卡免费看| 免费人妻精品一区二区三区视频| 91在线精品国自产拍蜜月| 高清在线视频一区二区三区| 母亲3免费完整高清在线观看 | av电影中文网址| 日韩欧美一区视频在线观看| av在线观看视频网站免费| 国精品久久久久久国模美| 97在线人人人人妻| 国产精品 国内视频| 国产亚洲最大av| 精品亚洲成国产av| 亚洲精品乱码久久久v下载方式| 天堂中文最新版在线下载| 91久久精品国产一区二区三区| 国产一级毛片在线| 欧美三级亚洲精品| 激情五月婷婷亚洲| 国产亚洲一区二区精品| 亚洲av二区三区四区| 夫妻午夜视频| 日本黄色片子视频| 久久国产精品大桥未久av| 一边亲一边摸免费视频| 国产无遮挡羞羞视频在线观看| 日韩熟女老妇一区二区性免费视频| 中国三级夫妇交换| 亚洲精品第二区| av有码第一页| 男女边吃奶边做爰视频| 欧美成人精品欧美一级黄| 你懂的网址亚洲精品在线观看| 久久久精品94久久精品| 纯流量卡能插随身wifi吗| 欧美三级亚洲精品| 18禁在线播放成人免费| 亚洲av在线观看美女高潮| 伊人亚洲综合成人网| av国产精品久久久久影院| 国产高清三级在线| 亚洲精品自拍成人| 黄色毛片三级朝国网站| 天天影视国产精品| 国产极品天堂在线| 丝袜美足系列| 久久久精品免费免费高清| 欧美人与性动交α欧美精品济南到 | 婷婷色麻豆天堂久久| 中文字幕av电影在线播放| 国产高清三级在线| 久久99精品国语久久久| 99九九在线精品视频| 蜜桃在线观看..| 女人精品久久久久毛片| 99热6这里只有精品| 在线天堂最新版资源| 青青草视频在线视频观看| 国产亚洲av片在线观看秒播厂| 蜜臀久久99精品久久宅男| 丰满少妇做爰视频| 亚洲国产精品999| 国产成人精品婷婷| 国产欧美日韩一区二区三区在线 | av卡一久久| 国产又色又爽无遮挡免| 亚洲av.av天堂| 欧美3d第一页| 日本黄大片高清| 亚洲,一卡二卡三卡| 午夜免费男女啪啪视频观看| 亚洲国产精品专区欧美| 国产日韩欧美视频二区| 你懂的网址亚洲精品在线观看| 视频中文字幕在线观看| 亚洲三级黄色毛片| 久久精品久久久久久噜噜老黄| 日韩 亚洲 欧美在线| 老司机影院毛片| 久久人人爽av亚洲精品天堂| 亚洲第一av免费看| 日韩伦理黄色片| 秋霞在线观看毛片| 亚洲欧美精品自产自拍| 高清不卡的av网站| 热99久久久久精品小说推荐| 精品午夜福利在线看| 中文字幕亚洲精品专区| 91午夜精品亚洲一区二区三区| 一级毛片我不卡| 人妻制服诱惑在线中文字幕| 激情五月婷婷亚洲| 国产精品麻豆人妻色哟哟久久| 制服诱惑二区| 久久久久精品性色| 亚洲人成77777在线视频| 一级二级三级毛片免费看| 国产日韩欧美在线精品| 少妇丰满av| 一边摸一边做爽爽视频免费| 永久网站在线| 午夜免费鲁丝| 热99久久久久精品小说推荐| 大码成人一级视频| 特大巨黑吊av在线直播| 亚洲精品久久久久久婷婷小说| 高清欧美精品videossex| 成年人免费黄色播放视频| 多毛熟女@视频| 免费黄频网站在线观看国产| av福利片在线| 国产日韩欧美视频二区| 亚洲欧美精品自产自拍| videossex国产| 黄色一级大片看看| 国产视频首页在线观看| 国产一区有黄有色的免费视频| 69精品国产乱码久久久| 最近手机中文字幕大全| 日本黄色日本黄色录像| 日本av手机在线免费观看| 两个人免费观看高清视频| 91成人精品电影| 日本午夜av视频| 少妇人妻精品综合一区二区| 人妻夜夜爽99麻豆av| 国模一区二区三区四区视频| 伊人久久精品亚洲午夜| 激情五月婷婷亚洲| 亚洲av成人精品一区久久| 国产伦精品一区二区三区视频9| 亚洲av综合色区一区| 蜜臀久久99精品久久宅男| 精品一区二区三区视频在线| 欧美xxⅹ黑人| 丰满少妇做爰视频| 国产成人精品在线电影| 国产黄片视频在线免费观看| 80岁老熟妇乱子伦牲交| xxx大片免费视频| 国产亚洲av片在线观看秒播厂| 观看av在线不卡| 日本黄色片子视频| 欧美精品一区二区大全| 丝袜脚勾引网站| 久久人妻熟女aⅴ| 一级毛片我不卡| 亚洲成人av在线免费| 国产免费视频播放在线视频| 99re6热这里在线精品视频| 亚洲情色 制服丝袜| 久久av网站| 欧美性感艳星| 国产精品人妻久久久久久| 国产成人精品婷婷| 国产成人精品一,二区| 99热网站在线观看| videos熟女内射| 最近最新中文字幕免费大全7| 乱人伦中国视频| 日韩不卡一区二区三区视频在线| 女的被弄到高潮叫床怎么办| 如日韩欧美国产精品一区二区三区 |