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

    應(yīng)用偽譜法的運(yùn)載火箭在線制導(dǎo)方法研究

    2017-03-31 02:21:05張志國余夢倫耿光有
    宇航學(xué)報(bào) 2017年3期
    關(guān)鍵詞:偽譜制導(dǎo)基點(diǎn)

    張志國,余夢倫,耿光有,宋 強(qiáng)

    (北京宇航系統(tǒng)工程研究所,北京100076)

    應(yīng)用偽譜法的運(yùn)載火箭在線制導(dǎo)方法研究

    張志國,余夢倫,耿光有,宋 強(qiáng)

    (北京宇航系統(tǒng)工程研究所,北京100076)

    研究Gauss偽譜法(GPM)在液體運(yùn)載火箭拋罩結(jié)束到入軌飛行段制導(dǎo)律設(shè)計(jì)中的應(yīng)用性。在每一個(gè)制導(dǎo)周期內(nèi),采用高效高精度數(shù)值軌跡優(yōu)化方法計(jì)算當(dāng)前制導(dǎo)周期內(nèi)的制導(dǎo)律。通過合理選擇非線性規(guī)劃問題的基點(diǎn)數(shù)量和制導(dǎo)周期,節(jié)省制導(dǎo)方法計(jì)算時(shí)間。將基于偽譜法的制導(dǎo)方法與運(yùn)載火箭中使用的迭代制導(dǎo)方法進(jìn)行對比,在保證同等入軌精度的條件下,該方法對于復(fù)雜約束問題處理方法更為便捷,滿足在線制導(dǎo)的需求。同時(shí)仿真表明,該方法能夠有效應(yīng)對各種偏差,是比較接近工程應(yīng)用的一種方法。

    偽譜法;實(shí)時(shí)在線制導(dǎo);軌跡優(yōu)化;非線性規(guī)劃(NLP)

    0 引 言

    制導(dǎo)技術(shù)是影響飛行器性能的關(guān)鍵技術(shù),從國內(nèi)外運(yùn)載火箭的發(fā)展歷程來看,制導(dǎo)方法經(jīng)歷了從開環(huán)制導(dǎo)到閉環(huán)制導(dǎo),從攝動制導(dǎo)、迭代制導(dǎo)到實(shí)時(shí)在線預(yù)測制導(dǎo)的發(fā)展。制導(dǎo)方法的發(fā)展主要圍繞兩個(gè)關(guān)鍵問題,即制導(dǎo)精度和制導(dǎo)效率。傳統(tǒng)制導(dǎo)方法受硬件性能限制,為保證制導(dǎo)效率主要采用解析方法,犧牲了部分制導(dǎo)精度,如攝動制導(dǎo),射前裝訂好制導(dǎo)參數(shù),可以在線實(shí)時(shí)制導(dǎo),但精度不高。隨后迭代制導(dǎo)方法的出現(xiàn),提升了火箭的制導(dǎo)精度,但由于迭代制導(dǎo)在方法中對引力場采用了近似處理,方法誤差無法完全消除,同時(shí)迭代制導(dǎo)對末端程序角不加限制,無法控制火箭入軌的姿態(tài),對于有姿態(tài)需求的任務(wù)需要單獨(dú)增加調(diào)姿段。隨著計(jì)算機(jī)技術(shù)的發(fā)展,使得高精度的制導(dǎo)數(shù)值計(jì)算方法成為可能,尤其是高效高精度軌跡優(yōu)化計(jì)算方法的發(fā)展,衍生出了一系列基于軌跡優(yōu)化方法的制導(dǎo)新方法。如基于數(shù)值間接法、直接法[1]、混合法[2]軌跡優(yōu)化的制導(dǎo)方法,本文采用的偽譜法即為直接法中的一種。

    每一次制導(dǎo)運(yùn)算實(shí)際上是進(jìn)行一次軌跡優(yōu)化計(jì)算,只是不需要獲得全局優(yōu)化和控制參數(shù),而只需要獲得當(dāng)前制導(dǎo)周期內(nèi)的控制變量?;跀?shù)值間接法的制導(dǎo)方法受到間接法初值敏感,計(jì)算時(shí)間長而在當(dāng)前計(jì)算機(jī)條件下實(shí)時(shí)制導(dǎo)應(yīng)用受到限制?;谥苯臃ǖ闹茖?dǎo)方法通過合理選擇問題的階數(shù),一方面相對間接法能夠節(jié)省制導(dǎo)方法計(jì)算時(shí)間,另一方面針對復(fù)雜的約束條件仍能獲得較高的制導(dǎo)精度。是比較接近工程的一種方法。

    偽譜法[3-4]將高度非線性微分方程的解用插值多項(xiàng)式表示出來,即通過有限個(gè)離散基點(diǎn)來近似狀態(tài)變量和控制變量,相應(yīng)基點(diǎn)處的微分方程、約束條件和系統(tǒng)積分指標(biāo)都可以表示成基點(diǎn)的代數(shù)方程。原微分問題轉(zhuǎn)化為標(biāo)準(zhǔn)非線性規(guī)劃(Nonlinear programming, NLP)問題,對于大規(guī)模NLP問題的求解方法主要有序列二次規(guī)劃(Sequential quadratic programming, SQP)法和內(nèi)點(diǎn)法,基于這兩種方法已有一些成熟的軟件,如基于SQP方法的SNOPT[5]、NPSOL等軟件,基于內(nèi)點(diǎn)法的IPOPT軟件[6],近年來,這些方法均得到了快速發(fā)展,使得偽譜法的應(yīng)用越來越廣泛。

    偽譜法在應(yīng)用過程中,若在Matlab環(huán)境下編程(如:采用GPOPS工具箱),計(jì)算精度高但效率偏低,適合用于離線計(jì)算,如文獻(xiàn)[7]中求解月球定點(diǎn)著陸問題和文獻(xiàn)[8]求解固體火箭上升段軌跡優(yōu)化問題,計(jì)算耗時(shí)為分鐘的量級。若要提高數(shù)值方法計(jì)算效率,從而滿足在線軌跡規(guī)劃和制導(dǎo)律設(shè)計(jì)的實(shí)時(shí)性需求,解決思路包括采用更高效的編程環(huán)境、簡化問題模型、改進(jìn)NLP問題求解方法的效率等,如文獻(xiàn)[9]將軌跡跟蹤制導(dǎo)問題轉(zhuǎn)化為線性時(shí)變系統(tǒng)調(diào)節(jié)器問題來進(jìn)行求解,文獻(xiàn)[10]中采用狀態(tài)量縮減的方法改進(jìn)閉環(huán)制導(dǎo)方法的計(jì)算效率。

    本文采用C語言版本SNOPT求解器,通過選取不同配點(diǎn)和制導(dǎo)周期,比較得出合適的制導(dǎo)參數(shù)選取方案,在保證入軌精度的同時(shí),進(jìn)一步提高制導(dǎo)律的求解效率,滿足在線設(shè)計(jì)需求。該方法能夠有效處理帶末端姿態(tài)約束火箭制導(dǎo)問題,處理過程簡便,并能很好地適應(yīng)一定范圍內(nèi)的雙向偏差。

    1 偽譜法軌跡優(yōu)化方法

    已知一般的非線性系統(tǒng)最優(yōu)控制問題,包括狀態(tài)方程、等式(不等式)約束、性能指標(biāo)三部分,整個(gè)最優(yōu)控制問題要求在滿足狀態(tài)方程和等式(不等式)約束的條件下,使得性能指標(biāo)達(dá)到極小值,統(tǒng)一形式可以寫為:

    (1)

    式中:x(t)、u(t)分別表示狀態(tài)變量和控制變量,t0、tf分別表示最優(yōu)控制問題的起始時(shí)間和末端時(shí)間,E和C分別表示等式約束和不等式約束。偽譜法插值多項(xiàng)式的選擇主要有Legendre多項(xiàng)式和Chebyshev多項(xiàng)式兩種,根據(jù)基點(diǎn)位置、插值多項(xiàng)式種類、積分形式的不同選取方法,偽譜法分為Legendre偽譜法、Gauss偽譜法、Radau偽譜法和Chebyshev偽譜法,幾種方法的計(jì)算效率和精度相當(dāng)[11],其中Gauss偽譜法的解滿足庫恩-塔克(Karush-Kuhn-Tucker,KKT)條件,利用余向量映射定理,可以證明Gauss偽譜法的解同間接法的解具有一致完備性,因此在軌跡優(yōu)化問題中被廣泛采用[12]。本文采用Gauss偽譜法,配點(diǎn)選擇正交Legendre多項(xiàng)式的零點(diǎn),積分采用Gauss積分公式,正交Legendre多項(xiàng)式的基點(diǎn)在(-1,1)之間取值,需要首先將時(shí)間變量[t0,tf]投影變換到[-1,1]內(nèi)。

    (2)

    將連續(xù)最優(yōu)控制問題轉(zhuǎn)換為以Legendre-Gauss基點(diǎn)處的狀態(tài)變量和控制變量為未知系數(shù)的標(biāo)準(zhǔn)NLP問題[4]。

    (3)

    式中:k為有限的離散基點(diǎn),Xk、Uk分別為離散時(shí)刻τk處的狀態(tài)變量和控制變量。采用Lagrange插值,Dk,i為Lagrange插值多項(xiàng)式的導(dǎo)數(shù),用于計(jì)算狀態(tài)量微分方程的等式約束。wk為Gauss積分公式中的積分權(quán)重,當(dāng)基點(diǎn)數(shù)給定時(shí)為常值。其中Dk,i的表達(dá)式為

    (4)

    對于式(3)構(gòu)成的適當(dāng)規(guī)模的NLP問題,目前有多種數(shù)值方法可以高效地給出求解結(jié)果。

    2 基于偽譜法的實(shí)時(shí)在線制導(dǎo)方法

    火箭與導(dǎo)彈的制導(dǎo)目標(biāo)不同,火箭需要瞄準(zhǔn)目標(biāo)軌道根數(shù)。迭代制導(dǎo)開始作用在拋罩后飛行段[13]。本文基于偽譜法的實(shí)時(shí)制導(dǎo)方法選擇同樣飛行狀態(tài)作為初始狀態(tài),制導(dǎo)開始時(shí),火箭已經(jīng)拋掉整流罩,在真空飛行,沒有空氣動力,主要受到重力和發(fā)動機(jī)推力影響,發(fā)動機(jī)參數(shù)偏差,主要包括比沖、推進(jìn)劑秒耗量偏差,推力作用線橫移和偏斜等。飛行過程中做瞬時(shí)平衡假設(shè),制導(dǎo)參數(shù)的選取僅包括俯仰角和偏航角兩個(gè)變量。因此火箭動力學(xué)微分方程可以寫為:

    (5)

    式中:狀態(tài)量x、y、z,Vx、Vy、Vz,m分別表示火箭的位置,速度和質(zhì)量。控制量φ、ψ分別表示火箭的俯仰程序角和偏航程序角。F(t) 表示火箭發(fā)動機(jī)的推力。性能指標(biāo)函數(shù)為:

    (6)

    對于推力大小不可調(diào)的火箭來說,時(shí)間最優(yōu)性能指標(biāo)等價(jià)于燃料最優(yōu)。考慮一級飛行結(jié)束時(shí)程序角已經(jīng)固定,因此制導(dǎo)開始時(shí),程序角初值是問題的一個(gè)約束,如果考慮火箭末端姿態(tài),過程中角加速度限制,最優(yōu)問題的約束方程可以寫為:

    (7)

    (8)

    (9)

    在每個(gè)制導(dǎo)周期內(nèi),根據(jù)火箭當(dāng)前的狀態(tài)量,進(jìn)行一次偽譜法軌跡優(yōu)化計(jì)算,獲得剩余飛行時(shí)間內(nèi)的控制變量基點(diǎn)值,通過Lagrange插值得到當(dāng)前狀態(tài)下的控制量,僅取當(dāng)前控制量作為火箭制導(dǎo)控制參數(shù)

    (10)

    式中:U0(τ=-1)即為當(dāng)前時(shí)刻的制導(dǎo)律。飛行全過程按照此方法進(jìn)行滾動時(shí)域優(yōu)化控制,即可實(shí)現(xiàn)火箭實(shí)時(shí)在線制導(dǎo),基于偽譜法的實(shí)時(shí)在線制導(dǎo)流程如圖1所示。

    偽譜法最大的優(yōu)勢在于能夠用更少的參數(shù)得到較高的計(jì)算精度。從理論上講,插值點(diǎn)選取的越多,計(jì)算結(jié)果精度將越高,但是隨著基點(diǎn)的增多將帶來計(jì)算量和計(jì)算時(shí)間的大大增加,這將使直接法計(jì)算效率的優(yōu)勢不再明顯。下面分別選擇基點(diǎn)數(shù)=4、8、12三種情況進(jìn)行單次偽譜法軌跡優(yōu)化方法計(jì)算,選擇同樣的收斂控制精度,統(tǒng)計(jì)單步偽譜法的計(jì)算時(shí)間,并獲得俯仰角和偏航角變化曲線圖2和圖3(數(shù)值仿真過程采用臺式機(jī)CPU3.19GHZ,內(nèi)存1.93GB)。

    表1 不同基點(diǎn)數(shù)量條件下偽譜制導(dǎo)法計(jì)算效率

    由表1可知,選擇同等要求的數(shù)值算法收斂精度,隨著離散基點(diǎn)數(shù)量的增加,單步偽譜法計(jì)算時(shí)間不斷增加,對于很少的基點(diǎn),偽譜法仍然能夠得到較高的收斂精度,但注意該精度是插值意義下的收斂精度,即保證通過配點(diǎn)進(jìn)行Lagrange插值得到的末端狀態(tài)收斂精度。

    實(shí)時(shí)制導(dǎo)周期的選擇至少要大于單步計(jì)算周期,考慮到選擇8~12個(gè)基點(diǎn)的單步偽譜法計(jì)算周期在100~200 ms的量級,本文分別采用10 s、5 s和1 s三種制導(dǎo)周期進(jìn)行實(shí)時(shí)制導(dǎo)方法計(jì)算,仿真結(jié)果如圖4和圖5所示??疾觳煌茖?dǎo)周期對入軌精度的影響,同時(shí)用制導(dǎo)周期為1 s的迭代制導(dǎo)[15]入軌精度作對比。

    從位置速度和軌道根數(shù)同標(biāo)稱入軌條件的偏差(見表2和表3)可以看出,偽譜制導(dǎo)法選取的制導(dǎo)周期越短,入軌精度越高,并且和迭代制導(dǎo)能夠獲得同樣量級的入軌精度。圖5的偏航程序角歷程和圖3有一定的差異,是由于圖3為單次全局軌跡優(yōu)化曲線;圖5為周期制導(dǎo)曲線,控制量通過偽譜法的有限個(gè)基點(diǎn)插值獲得,并且在每個(gè)制導(dǎo)周期內(nèi)才更新,導(dǎo)致積分獲得的狀態(tài)量和標(biāo)準(zhǔn)彈道略有差異,求解最優(yōu)控制問題進(jìn)而引起偏航程序角微小變化。偏航程序角設(shè)計(jì)是在0°附近,火箭基本在射面內(nèi)飛行,雖然圖上顯示趨勢不同,但是相對量級很小(在1°以內(nèi)),符合工程設(shè)計(jì)要求。而俯仰程序角由于基準(zhǔn)值較大,圖4和圖2曲線一致。

    在制導(dǎo)末段,為了防止制導(dǎo)程序角出現(xiàn)不穩(wěn)定情況,在最后兩到三個(gè)制導(dǎo)周期內(nèi)不進(jìn)行偽譜法更新計(jì)算,而是采用前一個(gè)制導(dǎo)周期已經(jīng)求得的末端最優(yōu)制導(dǎo)程序角進(jìn)行插值。因此和末端近似的迭代制導(dǎo)方法獲得同等量級的入軌精度,如果想進(jìn)一步提升入軌精度并且工程上可用,無論何種方法都需要首先解決末端程序角跳變不穩(wěn)定問題。但其實(shí)當(dāng)前方法的入軌精度相對于導(dǎo)航精度已經(jīng)足夠。

    表2 不同制導(dǎo)周期下位置速度入軌精度

    表3 不同制導(dǎo)周期下軌道根數(shù)入軌精度

    表4 不同末端條件下軌道根數(shù)的入軌精度

    3 末端帶約束制導(dǎo)方法研究

    迭代制導(dǎo)對末端程序角不加限制,無法控制火箭入軌的姿態(tài),但對于基于偽譜法的制導(dǎo)方法不存在這一問題,因?yàn)閭巫V法的優(yōu)勢就在于處理約束復(fù)雜的問題,無論是過程約束還是末端約束,不需要解析推導(dǎo)新的一階必要條件,僅需通過將約束施加在原問題的離散點(diǎn)上,即可求得有約束問題的制導(dǎo)律控制??紤]末端有姿態(tài)角約束的問題,分別選擇俯仰角-10°、-20°、-30°三種狀態(tài),偏航角0°、-2°、+2°三種狀態(tài),用偽譜制導(dǎo)法進(jìn)行制導(dǎo)控制得到軌道根數(shù)的入軌精度(見表4)。

    從圖6和圖7可以看出,對于末端有姿態(tài)角約束的問題,基于偽譜法的制導(dǎo)算法能夠獲得與末端無約束問題同等量級入軌精度的制導(dǎo)率,入軌軌道根數(shù)中半長軸、偏心率、軌道傾角、升交點(diǎn)赤經(jīng)均在1.0×10-4(°)量級以下,唯一變化的是影響入軌位置的真近點(diǎn)角,這是由于對于火箭來說,不考慮發(fā)動機(jī)搖擺,火箭的姿態(tài)決定了推力的方向,同時(shí)假設(shè)火箭一直滿推的情況,當(dāng)末端增加姿態(tài)約束時(shí),如果目標(biāo)仍然瞄準(zhǔn)初始設(shè)計(jì)的全部六個(gè)軌道根數(shù),新的最優(yōu)軌跡可能已不再是運(yùn)載能力最優(yōu)的軌跡[14],因此需要放松一個(gè)瞄準(zhǔn)參數(shù),即真近點(diǎn)角,調(diào)整合適的入軌位置,保證運(yùn)載能力最大,如果采用緯度幅角參數(shù)(u=ω+f),則其入軌精度也可達(dá)到1.0×10-4(°)量級。

    從不同末端約束條件下的控制程序角曲線(圖6和圖7)可以看出,俯仰和偏航程序角的斜率都發(fā)生了變化,也可以從另一個(gè)側(cè)面說明新的制導(dǎo)控制飛行路線已經(jīng)不再是原來的標(biāo)準(zhǔn)彈道,如果初始程序角沒有變化,真正決定施加末端約束條件下能否可控的條件是角加速度的控制能力,尤其是初始程序角能否快速過渡到最優(yōu)控制程序角曲線上來。另一方面,該方法在控制末端程序角雖然采用插值控制,但基本是線性變化,即無論前期程序角如何變化,末端程序角都表現(xiàn)為線性控制,這樣利于工程實(shí)施的可行性和制導(dǎo)控制的可靠性。

    4 制導(dǎo)方法應(yīng)對偏差的收斂性

    火箭實(shí)際飛行過程中會受到各種偏差干擾的影響,因此需要考察制導(dǎo)方法對于各種偏差的適應(yīng)性。理論上只要偏差沒有超出火箭能力所能到達(dá)的范圍時(shí),制導(dǎo)方法應(yīng)該能夠應(yīng)對各種偏差并且控制火箭到達(dá)既定軌道。這里參考茹家欣[13]在迭代制導(dǎo)方法中加偏差測試的方法,將火箭各種狀態(tài)偏差折合到出發(fā)時(shí)刻的偏差。文中考慮火箭二級飛行段制導(dǎo)律設(shè)計(jì),因?yàn)槎夛w行初始偏差主要來自于起飛和一級飛行段,火箭制導(dǎo)的誤差來源主要包括工具誤差和方法誤差:工具誤差的主要來源是慣性器件的測量偏差,根據(jù)現(xiàn)有一級飛行段工程中使用的慣性器件所能達(dá)到的測量精度進(jìn)行估算,如速度、位置的偏差;方法誤差包括質(zhì)量、發(fā)動機(jī)偏差、環(huán)境等因素,如文中考慮質(zhì)量約±1%的偏差。分別考慮正負(fù)極偏差的情況,選擇同樣的入軌條件,得到基于偽譜法的制導(dǎo)方法應(yīng)對偏差制導(dǎo)結(jié)果如表5所示。

    表5 正負(fù)極初始偏差條件下軌道根數(shù)的入軌精度

    對于含偏差問題,如果瞄準(zhǔn)初始設(shè)計(jì)的六個(gè)軌道根數(shù),同樣有運(yùn)載能力變化的問題,因此仍然考慮瞄準(zhǔn)五個(gè)入軌軌道根數(shù),調(diào)整入軌位置,將真近點(diǎn)角約束放松作為f0±10°的約束,通過施加幾組正負(fù)偏差,數(shù)值仿真表明,偽譜制導(dǎo)法都能夠很好的適應(yīng),并且保證較高的入軌精度。

    5 結(jié)束語

    本文設(shè)計(jì)了基于偽譜法的運(yùn)載火箭在線制導(dǎo)方法,針對火箭拋罩結(jié)束到入軌動力飛行段制導(dǎo)律進(jìn)行數(shù)值仿真,結(jié)果表明基于偽譜法的制導(dǎo)方法應(yīng)用于實(shí)時(shí)制導(dǎo)律設(shè)計(jì)時(shí),確有較高的制導(dǎo)效率和制導(dǎo)精度,可以和迭代制導(dǎo)有同樣甚至更高的入軌精度。選擇合適的基點(diǎn)數(shù)量和制導(dǎo)周期可以平衡制導(dǎo)效率和制導(dǎo)精度的關(guān)系。該方法在處理多約束制導(dǎo)問題時(shí)處理簡單,效率高,并且能很好的適應(yīng)一定范圍內(nèi)的雙向偏差,為偽譜制導(dǎo)方法的工程應(yīng)用提供一定參考。

    [1] Lu P, Pan B. Highly constrained optimal launch ascent guidance[J]. Journal of Guidance Control & Dynamics, 2012, 33(2):404-414.

    [2] 崔乃剛, 黃盤興, 韋常柱, 等. 基于混合優(yōu)化的運(yùn)載器大氣層內(nèi)閉環(huán)制導(dǎo)方法[J]. 中國慣性技術(shù)學(xué)報(bào), 2015,23(3):328-333. [Cui Nai-gang, Huang Pan-xing, Wei Chang-zhu, et al.Colsed-loop endo-atmospheric guidance of launch vehicle based on hybrid optimization approach [J]. Journal of Chinese Inertial Technology, 2015,23(3):328-333.]

    [3] Fahroo F, Michael I. Advances in Pseudo-spectral methods for optimal control[C]. AIAA Guidance, Navigation and Control Conference, Honolulu, Hawaii, USA, August 18-21, 2008.

    [4] 孫勇. 基于改進(jìn)Gauss偽譜法的高超聲速飛行器軌跡優(yōu)化與制導(dǎo)[D]. 哈爾濱: 哈爾濱工業(yè)大學(xué), 2012. [Sun Yong. Trajectory optimization and guidance of hypersonic vehicle based on improved Gauss Pseudo-spectral method [D]. Harbin: Harbin Institute of Technology, 2012.]

    [5] Gill P E, Murray W, Saunders M A. SNOPT: An SQP algorithm for large scale constrained optimization [J]. SIAM Journal on Optimization, 2002, 12(4):979-1006.

    [6] W?chter A, Biegler L T. On the implementation of an interior-point filter line-search algorithm for large-scale nonlinear programming [J]. Mathematical Programming, 2006, 106(1):25-57.

    [7] 彭祺擘, 李海陽, 沈紅新. 基于高斯—偽譜法的月球定點(diǎn)著陸軌道快速優(yōu)化設(shè)計(jì)[J]. 宇航學(xué)報(bào), 2010, 31(4):1012-1016. [Peng Qi-bo, Li Hai-yang, Shen Hong-xin. Rapid lunar exact-landing trajectory optimization via gauss Pseudo-spectral method [J]. Journal of Astronautics, 2010, 31(4):1012-1016.]

    [8] 楊希祥, 張為華. 基于Gauss偽譜法的固體運(yùn)載火箭上升段軌跡快速優(yōu)化研究[J]. 宇航學(xué)報(bào), 2011, 32(1):15-21. [Yang Xi-xiang, Zhang Wei-hua. Rapid optimization of ascent trajectory for solid launch vehicles based on Gauss Pseudo-spectral method [J]. Journal of Astronautics, 2011, 32(1):15-21.]

    [9] 廖宇新, 李惠峰, 包為民. 基于間接Radau偽譜法的滑翔段軌跡跟蹤制導(dǎo)律[J]. 宇航學(xué)報(bào), 2015, 36(12):1398-1405. [Liao Yu-xin, Li Hui-feng, Bao Wei-min. Gliding trajectory tracking guidance law based on indirect Radau Pseudo-spectral method [J]. Journal of Astronautics, 2015, 36(12):1398-1405.]

    [10] 閆循良, 廖守億, 何泰,等. 基于改進(jìn)偽譜反饋控制的遠(yuǎn)程變軌閉環(huán)制導(dǎo)[J]. 宇航學(xué)報(bào), 2016, 37(1):84-93. [Yan Xun-liang, Liao Shou-yi, He Tai, et al. Modified Pseudo-spectral feedback control based closed-loop guidance for long-range orbit maneuver [J]. Journal of Astronautics, 2016, 37(1):84-93.]

    [11] Huntington G T, Benson D, Rao A V. A comparison of accuracy and computational efficiency of three pseudo-spectral methods[C]. Proceedings of the AIAA Guidance, Navigation, and Control Conference, Washington, DC, USA, 2007.

    [12] Drake D, Xin M, Balakrishnan S N. New nonlinear control technique for ascent phase of reusable launch vehicles [J]. Journal of Guidance, Control, and Dynamics, 2004, 27(6):930-937.

    [13] 茹家欣. 液體運(yùn)載火箭的一種迭代制導(dǎo)方法[J], 中國科學(xué), 2009, 39(4):696-706. [Ru Jia-xin. An iterative guidance method of liquid rocket [J]. Science in China, 2009, 39(4):696-706.]

    [14] 韓祝齋. 用于大型運(yùn)載火箭的迭代制導(dǎo)方法[J]. 宇航學(xué)報(bào), 1983(1). [Han Zhu-zhai. An iterative guidance method for the large launch vehicle [J]. Journal of Astronautics, 1983(1).]

    [15] 陳新民, 余夢倫. 迭代制導(dǎo)在運(yùn)載火箭上的應(yīng)用研究[J]. 宇航學(xué)報(bào), 2003, 24(5):484-489. [Chen Xin-min, Yu Meng-lun. Study of iterative guidance application to launch vehicles [J]. Journal of Astronautics, 2003, 24(5):484-489.]

    通信地址:北京市豐臺區(qū)南大紅門路1號(100076)

    電話:(010)68758250

    E-mail:zhangzhiguo08@yeah.net

    (編輯:牛苗苗)

    Research on Application of Pseudo-Spectral Method in Online Guidance Method for a Launch Vehicle

    ZHANG Zhi-guo, YU Meng-lun, GENG Guang-you, SONG Qiang

    (Beijing Institute of Aerospace System Engineering, Beijing 100076, China)

    The study conducts an application research on the guidance law design of the ascent trajectory for a launch vehicle based on the Gauss pseudo-spectral method (GPM). During each period, high efficiency and high precision trajectory optimization numerical methods are used to design the current guidance law. Through making appropriate choices about the number of nodes and the period of guidance, this method is not only to reduce the calculation time, but also to solve the constrained problem rapidly. Compared to the iterative guidance method, the GPM guidance method is much more easily to deal with the constrained problem under the same orbit injection precision. Simulations indicate that this method can also effectively deal with a variety of deviations, so it is very closer to the engineering application.

    Pseudo-spectral method; Real-time online guidance; Trajectory optimization; Nonlinear programming (NLP)

    2016-09-05;

    2017-01-03

    V448.13

    A

    1000-1328(2017)03-0262-08

    10.3873/j.issn.1000-1328.2017.03.006

    張志國(1991-),男,博士生,主要從事飛行器軌跡優(yōu)化、制導(dǎo)與控制方法研究。

    猜你喜歡
    偽譜制導(dǎo)基點(diǎn)
    矩陣偽譜的新定位集及其在土壤生態(tài)系統(tǒng)的應(yīng)用
    紊流環(huán)境下四維軌跡優(yōu)化的偽譜方法研究
    基于MPSC和CPN制導(dǎo)方法的協(xié)同制導(dǎo)律
    基于在線軌跡迭代的自適應(yīng)再入制導(dǎo)
    帶有攻擊角約束的無抖振滑模制導(dǎo)律設(shè)計(jì)
    偽譜法及其在飛行器軌跡優(yōu)化設(shè)計(jì)領(lǐng)域的應(yīng)用綜述*
    基于Radau偽譜方法的軌跡優(yōu)化
    復(fù)合制導(dǎo)方式確保精確入軌
    太空探索(2014年1期)2014-07-10 13:41:49
    亚洲第一区二区三区不卡| 亚洲精华国产精华液的使用体验| 成人av在线播放网站| 人妻少妇偷人精品九色| 亚洲成人久久爱视频| 国内精品宾馆在线| 午夜精品在线福利| 国产人妻一区二区三区在| 男人爽女人下面视频在线观看| 国内精品宾馆在线| 免费黄色在线免费观看| 人妻少妇偷人精品九色| 在线免费观看不下载黄p国产| 久久人人爽人人片av| 国产单亲对白刺激| 国产精品久久久久久久久免| 非洲黑人性xxxx精品又粗又长| 麻豆av噜噜一区二区三区| 亚洲精品中文字幕在线视频 | 欧美成人a在线观看| 欧美成人a在线观看| 精品久久久精品久久久| 好男人视频免费观看在线| 性插视频无遮挡在线免费观看| 一夜夜www| 在线观看美女被高潮喷水网站| 最近视频中文字幕2019在线8| 日日干狠狠操夜夜爽| 国产精品久久视频播放| 高清欧美精品videossex| 日本欧美国产在线视频| 精品欧美国产一区二区三| 2021少妇久久久久久久久久久| 国产老妇伦熟女老妇高清| av国产久精品久网站免费入址| 春色校园在线视频观看| 在线 av 中文字幕| 午夜爱爱视频在线播放| 免费黄色在线免费观看| 麻豆av噜噜一区二区三区| 亚洲图色成人| 老女人水多毛片| 国产成人午夜福利电影在线观看| 伊人久久精品亚洲午夜| 美女脱内裤让男人舔精品视频| 成人毛片a级毛片在线播放| 国产亚洲5aaaaa淫片| 亚洲精品久久久久久婷婷小说| 免费观看a级毛片全部| 色综合亚洲欧美另类图片| 国产一区亚洲一区在线观看| 亚洲欧洲日产国产| 青青草视频在线视频观看| 少妇猛男粗大的猛烈进出视频 | 国产国拍精品亚洲av在线观看| 亚州av有码| 亚洲激情五月婷婷啪啪| 永久免费av网站大全| 一区二区三区高清视频在线| 肉色欧美久久久久久久蜜桃 | 国产老妇女一区| 久久久久久久久久人人人人人人| 国产有黄有色有爽视频| 两个人的视频大全免费| 我的老师免费观看完整版| 日本黄大片高清| 夫妻性生交免费视频一级片| 成年版毛片免费区| 国产黄片美女视频| 色综合色国产| 免费看光身美女| 国产精品99久久久久久久久| 国产伦理片在线播放av一区| 亚洲性久久影院| 黄色配什么色好看| 国产免费一级a男人的天堂| av在线蜜桃| 国产一区二区三区av在线| 日日撸夜夜添| 国产成人精品婷婷| 国产女主播在线喷水免费视频网站 | av.在线天堂| 久久久亚洲精品成人影院| 亚洲精品aⅴ在线观看| 精品人妻偷拍中文字幕| 亚洲欧美日韩无卡精品| 菩萨蛮人人尽说江南好唐韦庄| 亚洲国产欧美在线一区| 国产一区二区亚洲精品在线观看| 少妇裸体淫交视频免费看高清| 两个人视频免费观看高清| 日日干狠狠操夜夜爽| 免费观看的影片在线观看| 午夜爱爱视频在线播放| 人妻系列 视频| 99热6这里只有精品| 国产亚洲一区二区精品| 久久久久久久国产电影| 国产日韩欧美在线精品| 日韩在线高清观看一区二区三区| 日本熟妇午夜| 全区人妻精品视频| 人人妻人人澡欧美一区二区| 亚洲av成人精品一二三区| 日本猛色少妇xxxxx猛交久久| 青青草视频在线视频观看| 插阴视频在线观看视频| 在线免费十八禁| 青春草国产在线视频| 精品99又大又爽又粗少妇毛片| 国产综合懂色| 午夜免费男女啪啪视频观看| 伊人久久精品亚洲午夜| 国产有黄有色有爽视频| 91午夜精品亚洲一区二区三区| 91久久精品电影网| 男女下面进入的视频免费午夜| 丰满人妻一区二区三区视频av| 久久久久久久午夜电影| 久久6这里有精品| 久久久色成人| 亚洲av不卡在线观看| 波多野结衣巨乳人妻| 少妇熟女aⅴ在线视频| 黄色日韩在线| 久久久久久九九精品二区国产| 久久久午夜欧美精品| 欧美bdsm另类| av又黄又爽大尺度在线免费看| a级毛片免费高清观看在线播放| 久久韩国三级中文字幕| 99热6这里只有精品| 熟妇人妻不卡中文字幕| av在线老鸭窝| 国产探花极品一区二区| 亚洲电影在线观看av| 成人毛片a级毛片在线播放| 欧美bdsm另类| 久久久久久国产a免费观看| 91精品伊人久久大香线蕉| 成人av在线播放网站| 18禁动态无遮挡网站| 精品国内亚洲2022精品成人| 亚洲,欧美,日韩| 国产精品一区二区性色av| 日韩精品有码人妻一区| www.色视频.com| 综合色av麻豆| 在现免费观看毛片| 久久久久久久国产电影| 免费看不卡的av| 日本免费在线观看一区| 91av网一区二区| 久久久国产一区二区| 亚洲婷婷狠狠爱综合网| 日韩一区二区三区影片| 亚洲av成人精品一区久久| 看十八女毛片水多多多| 亚洲精品,欧美精品| 国产一区亚洲一区在线观看| 国产高清有码在线观看视频| av女优亚洲男人天堂| 亚洲人成网站在线播| 亚洲国产日韩欧美精品在线观看| 精品熟女少妇av免费看| 2021天堂中文幕一二区在线观| 免费不卡的大黄色大毛片视频在线观看 | 欧美xxxx黑人xx丫x性爽| 久久精品夜色国产| 国产在视频线在精品| 日韩,欧美,国产一区二区三区| 久久热精品热| 爱豆传媒免费全集在线观看| 免费观看无遮挡的男女| 久久精品夜色国产| 国产 一区 欧美 日韩| 国产黄片美女视频| 亚洲精品乱码久久久久久按摩| 亚州av有码| 六月丁香七月| 日韩 亚洲 欧美在线| 精品99又大又爽又粗少妇毛片| 国产伦理片在线播放av一区| 一级毛片黄色毛片免费观看视频| 少妇被粗大猛烈的视频| 久久久久久国产a免费观看| 美女cb高潮喷水在线观看| 亚洲国产成人一精品久久久| 国产成人精品一,二区| 成年版毛片免费区| 亚洲精品日韩在线中文字幕| 中文资源天堂在线| 亚洲经典国产精华液单| 日韩av在线免费看完整版不卡| 黑人高潮一二区| 天天躁日日操中文字幕| 成人二区视频| 赤兔流量卡办理| 国产免费福利视频在线观看| 久久99热这里只有精品18| 色播亚洲综合网| 亚洲av电影不卡..在线观看| 日韩成人伦理影院| 国产午夜福利久久久久久| 国产一区二区在线观看日韩| 最近中文字幕高清免费大全6| 久久精品夜色国产| 日韩视频在线欧美| 高清午夜精品一区二区三区| 看非洲黑人一级黄片| 最近中文字幕高清免费大全6| 麻豆乱淫一区二区| 国产淫语在线视频| 亚洲经典国产精华液单| av播播在线观看一区| 青春草亚洲视频在线观看| 国产黄色视频一区二区在线观看| 欧美高清性xxxxhd video| 午夜视频国产福利| 国产女主播在线喷水免费视频网站 | 国产视频首页在线观看| 久久午夜福利片| 久久精品久久久久久久性| a级一级毛片免费在线观看| 欧美日韩一区二区视频在线观看视频在线 | 亚洲av电影不卡..在线观看| 午夜精品一区二区三区免费看| 欧美不卡视频在线免费观看| freevideosex欧美| 精品久久国产蜜桃| 日本欧美国产在线视频| 男女那种视频在线观看| 99久久九九国产精品国产免费| 国产高清不卡午夜福利| 国产 一区精品| 七月丁香在线播放| 欧美zozozo另类| 国产精品久久久久久久久免| 国产精品人妻久久久久久| 午夜福利在线观看免费完整高清在| 亚洲欧洲国产日韩| 精品一区在线观看国产| 最近的中文字幕免费完整| 天堂网av新在线| 亚洲国产高清在线一区二区三| 少妇被粗大猛烈的视频| 欧美日韩视频高清一区二区三区二| 国产乱来视频区| 久久久a久久爽久久v久久| 国产熟女欧美一区二区| 少妇人妻一区二区三区视频| 亚洲精品国产成人久久av| 亚洲欧美一区二区三区国产| 精品一区二区三卡| 色网站视频免费| 日韩中字成人| 国产综合精华液| 国产精品一区www在线观看| 免费黄网站久久成人精品| 十八禁网站网址无遮挡 | 蜜臀久久99精品久久宅男| 97在线视频观看| 亚洲国产色片| 国产成人a区在线观看| 国产有黄有色有爽视频| 亚洲伊人久久精品综合| videossex国产| 老师上课跳d突然被开到最大视频| 亚洲经典国产精华液单| 亚洲欧洲国产日韩| 亚洲aⅴ乱码一区二区在线播放| 天堂俺去俺来也www色官网 | 寂寞人妻少妇视频99o| 七月丁香在线播放| 国产av不卡久久| 69av精品久久久久久| 建设人人有责人人尽责人人享有的 | 小蜜桃在线观看免费完整版高清| 永久网站在线| 亚洲国产av新网站| 啦啦啦中文免费视频观看日本| 中文字幕亚洲精品专区| 亚洲最大成人中文| 汤姆久久久久久久影院中文字幕 | 少妇熟女欧美另类| 人妻一区二区av| 只有这里有精品99| 激情 狠狠 欧美| 51国产日韩欧美| 国产成人午夜福利电影在线观看| 麻豆久久精品国产亚洲av| 亚洲av男天堂| 国产精品久久视频播放| 亚洲精品一区蜜桃| 亚洲经典国产精华液单| 日韩不卡一区二区三区视频在线| 69人妻影院| 中文字幕亚洲精品专区| 成人亚洲精品一区在线观看 | 欧美zozozo另类| 99热这里只有精品一区| 精品国产露脸久久av麻豆 | 国产视频内射| 日本三级黄在线观看| 亚洲经典国产精华液单| 久久久成人免费电影| av国产久精品久网站免费入址| 熟妇人妻不卡中文字幕| 国产成人freesex在线| 亚洲精品,欧美精品| 日韩成人av中文字幕在线观看| 男人狂女人下面高潮的视频| 80岁老熟妇乱子伦牲交| 日韩人妻高清精品专区| 有码 亚洲区| 久久久久精品性色| 国产成人精品一,二区| 成人亚洲欧美一区二区av| 国产成人91sexporn| 一区二区三区免费毛片| 国产在线一区二区三区精| 色综合色国产| 亚洲欧洲国产日韩| 日韩欧美精品免费久久| 亚洲人成网站高清观看| 成人毛片a级毛片在线播放| 青春草亚洲视频在线观看| 亚洲成人一二三区av| 中国美白少妇内射xxxbb| 男女下面进入的视频免费午夜| 亚洲国产最新在线播放| 别揉我奶头 嗯啊视频| 三级经典国产精品| 国产 亚洲一区二区三区 | 久久午夜福利片| 国产成年人精品一区二区| 国产精品久久久久久av不卡| 亚洲av成人精品一区久久| 久久精品人妻少妇| 最近中文字幕高清免费大全6| av黄色大香蕉| av.在线天堂| 夫妻午夜视频| 亚洲成人久久爱视频| 高清视频免费观看一区二区 | 51国产日韩欧美| 啦啦啦韩国在线观看视频| 久久热精品热| 乱人视频在线观看| 婷婷色av中文字幕| 内射极品少妇av片p| 少妇被粗大猛烈的视频| videossex国产| 日韩电影二区| 午夜福利在线在线| 一夜夜www| 九九爱精品视频在线观看| 国内精品美女久久久久久| 丰满乱子伦码专区| 亚洲国产成人一精品久久久| 亚洲经典国产精华液单| 国产男人的电影天堂91| 国产亚洲av嫩草精品影院| 国产伦精品一区二区三区四那| 汤姆久久久久久久影院中文字幕 | 国产精品一二三区在线看| 少妇高潮的动态图| 国产成人91sexporn| 久久久色成人| 尾随美女入室| 少妇丰满av| 亚洲四区av| 国产在线男女| 看非洲黑人一级黄片| 一本久久精品| 中国美白少妇内射xxxbb| 少妇被粗大猛烈的视频| 国语对白做爰xxxⅹ性视频网站| 中文字幕久久专区| 大话2 男鬼变身卡| 搡老妇女老女人老熟妇| 高清在线视频一区二区三区| 日本免费在线观看一区| 亚洲av福利一区| 成年免费大片在线观看| 伊人久久精品亚洲午夜| 久久这里只有精品中国| 91在线精品国自产拍蜜月| 日日啪夜夜爽| 亚洲精品日韩在线中文字幕| 成人午夜高清在线视频| 欧美成人精品欧美一级黄| 九草在线视频观看| 亚洲自偷自拍三级| 美女高潮的动态| 亚洲成人av在线免费| 国产成人免费观看mmmm| 禁无遮挡网站| 久久这里有精品视频免费| 国产免费一级a男人的天堂| 久久精品综合一区二区三区| 亚洲精品456在线播放app| 黄色一级大片看看| 亚洲av福利一区| 天堂俺去俺来也www色官网 | 国产一区二区三区av在线| 美女内射精品一级片tv| 美女被艹到高潮喷水动态| 欧美日韩精品成人综合77777| 亚洲最大成人手机在线| 欧美+日韩+精品| 高清av免费在线| 如何舔出高潮| 成人毛片60女人毛片免费| 超碰av人人做人人爽久久| 国产精品麻豆人妻色哟哟久久 | av女优亚洲男人天堂| 亚洲,欧美,日韩| 麻豆国产97在线/欧美| 午夜精品国产一区二区电影 | 久久久久久久久中文| 又大又黄又爽视频免费| 婷婷色综合www| 搞女人的毛片| 亚洲自偷自拍三级| 欧美极品一区二区三区四区| 亚洲av日韩在线播放| 91在线精品国自产拍蜜月| 精品国内亚洲2022精品成人| 亚洲真实伦在线观看| 日韩av在线大香蕉| 黄色欧美视频在线观看| 国产成人一区二区在线| 国产乱人视频| 一二三四中文在线观看免费高清| 丝袜喷水一区| 国产午夜精品久久久久久一区二区三区| 欧美xxxx黑人xx丫x性爽| 亚洲精华国产精华液的使用体验| 精品一区二区免费观看| 在线免费观看不下载黄p国产| 三级经典国产精品| 久久精品综合一区二区三区| 国产91av在线免费观看| 欧美精品一区二区大全| 一个人观看的视频www高清免费观看| 亚洲人成网站高清观看| 少妇人妻一区二区三区视频| 麻豆成人午夜福利视频| 久久鲁丝午夜福利片| 国产精品三级大全| av线在线观看网站| 国产黄色免费在线视频| 亚洲综合精品二区| 高清视频免费观看一区二区 | 亚洲天堂国产精品一区在线| 久久久久久久久久成人| 丰满人妻一区二区三区视频av| 国产精品一二三区在线看| 大陆偷拍与自拍| 亚洲精品久久午夜乱码| 精品一区二区三卡| 在线观看一区二区三区| 国产美女午夜福利| 欧美最新免费一区二区三区| 欧美一区二区亚洲| 中文字幕久久专区| 亚洲av成人精品一二三区| 午夜激情福利司机影院| 黄色一级大片看看| 国产亚洲午夜精品一区二区久久 | 极品教师在线视频| 寂寞人妻少妇视频99o| 成人午夜高清在线视频| 久久久精品免费免费高清| av免费在线看不卡| 啦啦啦啦在线视频资源| 午夜福利在线观看免费完整高清在| 嫩草影院新地址| 秋霞伦理黄片| 日韩大片免费观看网站| 99热6这里只有精品| 免费观看av网站的网址| 国产黄色小视频在线观看| 成人一区二区视频在线观看| 亚洲综合色惰| 成年免费大片在线观看| 亚洲电影在线观看av| 亚洲欧美日韩卡通动漫| 国产成人午夜福利电影在线观看| 日韩欧美三级三区| 精品久久久久久成人av| 午夜福利在线在线| 国产老妇伦熟女老妇高清| 国产男女超爽视频在线观看| 免费看美女性在线毛片视频| 国产免费又黄又爽又色| 久久久成人免费电影| 晚上一个人看的免费电影| 亚洲精品亚洲一区二区| 美女cb高潮喷水在线观看| .国产精品久久| 嫩草影院新地址| 国产午夜精品一二区理论片| 亚洲自偷自拍三级| 狠狠精品人妻久久久久久综合| 亚洲真实伦在线观看| 久久久久精品性色| 永久网站在线| 精品久久久久久久末码| 人妻制服诱惑在线中文字幕| 91久久精品国产一区二区三区| 黄色欧美视频在线观看| 99久久中文字幕三级久久日本| 国产久久久一区二区三区| 免费看a级黄色片| 狂野欧美白嫩少妇大欣赏| 高清午夜精品一区二区三区| 一二三四中文在线观看免费高清| 熟妇人妻久久中文字幕3abv| 91久久精品国产一区二区三区| 老女人水多毛片| 日本wwww免费看| 中文字幕免费在线视频6| 日韩精品青青久久久久久| av女优亚洲男人天堂| 欧美日韩一区二区视频在线观看视频在线 | 伦精品一区二区三区| 国产一区二区亚洲精品在线观看| 亚洲人成网站高清观看| 亚州av有码| 国产乱人偷精品视频| 精品久久久久久久久亚洲| 成年免费大片在线观看| 亚洲av中文av极速乱| 久久午夜福利片| 搞女人的毛片| 偷拍熟女少妇极品色| 国产在线男女| 国产一区亚洲一区在线观看| 国产色爽女视频免费观看| 久久久久久伊人网av| 中文字幕免费在线视频6| 黄片无遮挡物在线观看| 99热6这里只有精品| 久久久亚洲精品成人影院| 久久久久网色| 肉色欧美久久久久久久蜜桃 | 一本久久精品| 人人妻人人澡人人爽人人夜夜 | 熟妇人妻不卡中文字幕| 自拍偷自拍亚洲精品老妇| 非洲黑人性xxxx精品又粗又长| 久久精品人妻少妇| 成年人午夜在线观看视频 | 高清av免费在线| 十八禁国产超污无遮挡网站| 久久久久久伊人网av| 又爽又黄a免费视频| 菩萨蛮人人尽说江南好唐韦庄| 丝瓜视频免费看黄片| 日韩伦理黄色片| 精品99又大又爽又粗少妇毛片| 在线观看一区二区三区| 亚洲人成网站在线播| 3wmmmm亚洲av在线观看| 成人二区视频| 街头女战士在线观看网站| 韩国高清视频一区二区三区| 亚洲真实伦在线观看| 久久99热6这里只有精品| 只有这里有精品99| 天堂俺去俺来也www色官网 | 乱人视频在线观看| 午夜免费观看性视频| 久久人人爽人人爽人人片va| 大又大粗又爽又黄少妇毛片口| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 成人漫画全彩无遮挡| 80岁老熟妇乱子伦牲交| 一边亲一边摸免费视频| 亚洲丝袜综合中文字幕| 一级片'在线观看视频| 免费少妇av软件| 国产成年人精品一区二区| 成人二区视频| 亚洲乱码一区二区免费版| 亚洲精品第二区| av在线观看视频网站免费| 韩国高清视频一区二区三区| 日日啪夜夜撸| 亚洲精品成人av观看孕妇| 亚洲精品日韩在线中文字幕| 水蜜桃什么品种好| 日韩精品青青久久久久久| 人妻系列 视频| 精品人妻视频免费看| 欧美 日韩 精品 国产| 久久精品夜夜夜夜夜久久蜜豆| 亚洲精品乱码久久久v下载方式| 男女下面进入的视频免费午夜| 亚洲美女搞黄在线观看| 99久久精品热视频| 最近最新中文字幕大全电影3| 777米奇影视久久| 免费人成在线观看视频色| 成人特级av手机在线观看| 91精品一卡2卡3卡4卡| 青春草亚洲视频在线观看| 少妇熟女aⅴ在线视频| 六月丁香七月| 亚洲经典国产精华液单| videossex国产| 日本黄大片高清| 秋霞在线观看毛片|