• <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
    欧美日韩黄片免| 老司机靠b影院| 美女福利国产在线| 亚洲精品国产一区二区精华液| 18禁黄网站禁片午夜丰满| a在线观看视频网站| 露出奶头的视频| 一二三四社区在线视频社区8| 欧美av亚洲av综合av国产av| 久久香蕉精品热| 高清欧美精品videossex| 精品国内亚洲2022精品成人 | 国产亚洲精品久久久久5区| 黄片小视频在线播放| 国产真人三级小视频在线观看| 每晚都被弄得嗷嗷叫到高潮| 成年版毛片免费区| 欧美性长视频在线观看| 午夜日韩欧美国产| 久久 成人 亚洲| 欧美精品av麻豆av| 法律面前人人平等表现在哪些方面| 精品国产一区二区久久| 国产av又大| 国产黄色免费在线视频| cao死你这个sao货| 丝袜美腿诱惑在线| 黄色视频,在线免费观看| 丰满的人妻完整版| 天天影视国产精品| 国产在线一区二区三区精| 免费在线观看影片大全网站| 亚洲精品成人av观看孕妇| 久久久久国产精品人妻aⅴ院 | 黄片小视频在线播放| 精品国产一区二区三区四区第35| av中文乱码字幕在线| 一区二区三区激情视频| 手机成人av网站| a级毛片黄视频| 国产亚洲精品一区二区www | 1024香蕉在线观看| 精品人妻在线不人妻| 亚洲国产精品sss在线观看 | 黑人欧美特级aaaaaa片| 国产高清国产精品国产三级| 美女扒开内裤让男人捅视频| 亚洲精品乱久久久久久| 久久久国产成人精品二区 | 一级作爱视频免费观看| 免费看十八禁软件| 亚洲五月婷婷丁香| 久久精品成人免费网站| 中亚洲国语对白在线视频| 亚洲aⅴ乱码一区二区在线播放 | 久久国产精品人妻蜜桃| 男女高潮啪啪啪动态图| 国产激情久久老熟女| 欧洲精品卡2卡3卡4卡5卡区| 老司机午夜福利在线观看视频| 国产区一区二久久| 性少妇av在线| 国产色视频综合| 亚洲av电影在线进入| 19禁男女啪啪无遮挡网站| tube8黄色片| 国产视频一区二区在线看| 国产日韩一区二区三区精品不卡| av网站在线播放免费| 狂野欧美激情性xxxx| 精品欧美一区二区三区在线| 99久久综合精品五月天人人| 欧美一级毛片孕妇| 久久久国产成人精品二区 | 亚洲人成伊人成综合网2020| 日韩有码中文字幕| 免费观看a级毛片全部| 欧美精品一区二区免费开放| 青草久久国产| 成人黄色视频免费在线看| 91国产中文字幕| 日本wwww免费看| 欧美日韩视频精品一区| 男人的好看免费观看在线视频 | 高清欧美精品videossex| 高清欧美精品videossex| 亚洲自偷自拍图片 自拍| 无限看片的www在线观看| 一本综合久久免费| 老司机亚洲免费影院| 欧美日韩国产mv在线观看视频| 国产亚洲精品一区二区www | 在线观看66精品国产| 成年女人毛片免费观看观看9 | 脱女人内裤的视频| 免费在线观看日本一区| 亚洲色图综合在线观看| 精品高清国产在线一区| 热re99久久国产66热| 黑人操中国人逼视频| 亚洲国产欧美网| 咕卡用的链子| 不卡av一区二区三区| 美女国产高潮福利片在线看| 国产乱人伦免费视频| 午夜两性在线视频| 老汉色∧v一级毛片| 天天躁夜夜躁狠狠躁躁| 99国产极品粉嫩在线观看| 久久久久久久国产电影| 亚洲伊人色综图| 国产在线观看jvid| av欧美777| 日韩欧美一区二区三区在线观看 | 国产精品一区二区在线观看99| 国产一区二区三区在线臀色熟女 | 人妻一区二区av| 在线天堂中文资源库| 国产一区二区激情短视频| 大陆偷拍与自拍| 精品国产乱码久久久久久男人| 中文字幕另类日韩欧美亚洲嫩草| 亚洲国产精品一区二区三区在线| 51午夜福利影视在线观看| 亚洲欧洲精品一区二区精品久久久| 亚洲av熟女| 久久久久久久精品吃奶| 久久久久国内视频| 波多野结衣一区麻豆| 国产成人av教育| 中出人妻视频一区二区| 在线播放国产精品三级| 搡老乐熟女国产| 免费在线观看日本一区| 咕卡用的链子| 欧美日韩av久久| 久久香蕉精品热| 飞空精品影院首页| 国产亚洲精品第一综合不卡| 久久精品熟女亚洲av麻豆精品| 亚洲第一青青草原| 捣出白浆h1v1| 久久中文字幕人妻熟女| 久久久久久久精品吃奶| 在线观看免费日韩欧美大片| 别揉我奶头~嗯~啊~动态视频| 很黄的视频免费| 99久久综合精品五月天人人| 在线播放国产精品三级| 欧美色视频一区免费| 亚洲一卡2卡3卡4卡5卡精品中文| 无人区码免费观看不卡| 国产黄色免费在线视频| 国产精品久久久久久人妻精品电影| 久久久久精品人妻al黑| 精品乱码久久久久久99久播| 19禁男女啪啪无遮挡网站| 高清黄色对白视频在线免费看| 国产精品久久视频播放| 欧美色视频一区免费| 操美女的视频在线观看| 欧美成人免费av一区二区三区 | 亚洲中文日韩欧美视频| 欧美日韩瑟瑟在线播放| 精品一区二区三区视频在线观看免费 | 激情视频va一区二区三区| 大型黄色视频在线免费观看| 国产不卡一卡二| 国产精品国产av在线观看| av视频免费观看在线观看| 99国产精品免费福利视频| 国产男女超爽视频在线观看| 满18在线观看网站| 亚洲人成电影观看| 老司机午夜十八禁免费视频| 午夜久久久在线观看| 久久午夜综合久久蜜桃| 看片在线看免费视频| 国产欧美日韩一区二区三| 日本黄色日本黄色录像| 亚洲欧美一区二区三区黑人| 国产黄色免费在线视频| 三级毛片av免费| 欧洲精品卡2卡3卡4卡5卡区| 日日夜夜操网爽| 国产精品乱码一区二三区的特点 | 国产男女超爽视频在线观看| 亚洲国产精品一区二区三区在线| 99riav亚洲国产免费| 日韩一卡2卡3卡4卡2021年| 国产精品永久免费网站| 亚洲 欧美一区二区三区| 中文字幕精品免费在线观看视频| 一二三四在线观看免费中文在| 黄网站色视频无遮挡免费观看| 精品人妻1区二区| 少妇被粗大的猛进出69影院| 久久久久久久午夜电影 | 久久精品人人爽人人爽视色| 高清黄色对白视频在线免费看| 制服诱惑二区| 女人被狂操c到高潮| 黄色a级毛片大全视频| netflix在线观看网站| 欧美精品亚洲一区二区| av线在线观看网站| 欧美午夜高清在线| 国产免费现黄频在线看| 嫁个100分男人电影在线观看| 波多野结衣一区麻豆| 国产成人欧美在线观看 | 日韩 欧美 亚洲 中文字幕| 99热只有精品国产| 亚洲av第一区精品v没综合| 成人黄色视频免费在线看| 18禁裸乳无遮挡动漫免费视频| 黑人猛操日本美女一级片| 久久精品国产综合久久久| 亚洲三区欧美一区| 国产xxxxx性猛交| 中文字幕人妻丝袜一区二区| 高清黄色对白视频在线免费看| 如日韩欧美国产精品一区二区三区| 99久久精品国产亚洲精品| 国产高清国产精品国产三级| 国产亚洲精品一区二区www | 交换朋友夫妻互换小说| 91国产中文字幕| 99久久综合精品五月天人人| 国产激情欧美一区二区| 亚洲专区国产一区二区| 十八禁高潮呻吟视频| 宅男免费午夜| 国产av精品麻豆| 成人手机av| 日韩欧美在线二视频 | 高清视频免费观看一区二区| 91av网站免费观看| 69精品国产乱码久久久| 国产精品 欧美亚洲| 中文字幕av电影在线播放| 亚洲精品成人av观看孕妇| 在线视频色国产色| 成人免费观看视频高清| 中国美女看黄片| 久久草成人影院| 国产1区2区3区精品| 黄色 视频免费看| 巨乳人妻的诱惑在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 欧美日本中文国产一区发布| 香蕉丝袜av| 黄色片一级片一级黄色片| √禁漫天堂资源中文www| 一本大道久久a久久精品| 欧美亚洲日本最大视频资源| 久久精品熟女亚洲av麻豆精品| 国产精品.久久久| 嫁个100分男人电影在线观看| 91麻豆av在线| 午夜福利影视在线免费观看| 婷婷成人精品国产| 宅男免费午夜| 91av网站免费观看| 精品无人区乱码1区二区| 精品一区二区三区av网在线观看| 欧美精品啪啪一区二区三区| 91字幕亚洲| 欧美黑人欧美精品刺激| 99国产精品一区二区蜜桃av | 欧美成人午夜精品| 日韩欧美国产一区二区入口| 中文字幕人妻丝袜一区二区| 国产精品久久久久久人妻精品电影| 99香蕉大伊视频| 成人18禁高潮啪啪吃奶动态图| 亚洲国产精品合色在线| 涩涩av久久男人的天堂| 深夜精品福利| 最近最新免费中文字幕在线| av网站在线播放免费| 啦啦啦在线免费观看视频4| 国产无遮挡羞羞视频在线观看| 亚洲美女黄片视频| 色94色欧美一区二区| 麻豆国产av国片精品| 国产片内射在线| 色婷婷久久久亚洲欧美| av线在线观看网站| 91成年电影在线观看| 亚洲人成电影观看| 欧美久久黑人一区二区| 精品少妇一区二区三区视频日本电影| 日本欧美视频一区| 国产欧美亚洲国产| 嫁个100分男人电影在线观看| www.精华液| 又黄又粗又硬又大视频| 国产亚洲精品第一综合不卡| av有码第一页| 欧美日韩视频精品一区| 亚洲精品国产色婷婷电影| 在线十欧美十亚洲十日本专区| 99国产极品粉嫩在线观看| 啦啦啦 在线观看视频| 国产有黄有色有爽视频| 99久久国产精品久久久| 窝窝影院91人妻| 妹子高潮喷水视频| 国产精品国产av在线观看| 高清黄色对白视频在线免费看| 亚洲欧美激情综合另类| 色播在线永久视频| 午夜视频精品福利| 深夜精品福利| 夜夜爽天天搞| 国产单亲对白刺激| 成人永久免费在线观看视频| 成人国产一区最新在线观看| 亚洲久久久国产精品| tocl精华| 高潮久久久久久久久久久不卡| 亚洲专区国产一区二区| 他把我摸到了高潮在线观看| 在线观看午夜福利视频| 多毛熟女@视频| 欧美日韩精品网址| 国产淫语在线视频| 交换朋友夫妻互换小说| 午夜久久久在线观看| 女性被躁到高潮视频| av中文乱码字幕在线| 久久中文字幕人妻熟女| 午夜老司机福利片| 亚洲 欧美一区二区三区| 欧美日韩黄片免| 国产精品九九99| 国产精品 国内视频| 黑人猛操日本美女一级片| 少妇 在线观看| 亚洲,欧美精品.| 这个男人来自地球电影免费观看| 91麻豆av在线| 亚洲欧美精品综合一区二区三区| 精品视频人人做人人爽| 一级片免费观看大全| 免费久久久久久久精品成人欧美视频| 久久久久国产精品人妻aⅴ院 | 美女 人体艺术 gogo| 真人做人爱边吃奶动态| 18禁观看日本| 久久国产精品大桥未久av| 最新在线观看一区二区三区| 欧美 亚洲 国产 日韩一| 极品少妇高潮喷水抽搐| 高清视频免费观看一区二区| 69精品国产乱码久久久| 中文字幕av电影在线播放| 日韩欧美国产一区二区入口| 国产男女超爽视频在线观看| 亚洲免费av在线视频| 久久精品亚洲精品国产色婷小说| 成人特级黄色片久久久久久久| 99热只有精品国产| 亚洲熟女精品中文字幕| 精品国内亚洲2022精品成人 | 欧美午夜高清在线| 成人手机av| 国产精品电影一区二区三区 | 亚洲熟女精品中文字幕| 丁香六月欧美| 午夜影院日韩av| 天堂中文最新版在线下载| 777米奇影视久久| 亚洲专区国产一区二区| 亚洲男人天堂网一区| 在线观看免费午夜福利视频| 91九色精品人成在线观看| 国产一区有黄有色的免费视频| 亚洲七黄色美女视频| 国产亚洲欧美98| a级毛片黄视频| 一进一出好大好爽视频| 久久精品亚洲av国产电影网| 亚洲成国产人片在线观看| 色尼玛亚洲综合影院| 亚洲熟女毛片儿| 国产一区二区三区综合在线观看| 日本一区二区免费在线视频| 精品久久久精品久久久| 久久狼人影院| 大香蕉久久成人网| 看片在线看免费视频| 国产精品久久久av美女十八| 亚洲精华国产精华精| 黄色a级毛片大全视频| 久久久久久久国产电影| 国产男女超爽视频在线观看| 狂野欧美激情性xxxx| 97人妻天天添夜夜摸| 欧美亚洲日本最大视频资源| 岛国毛片在线播放| 操美女的视频在线观看| 亚洲国产精品合色在线| 久久天堂一区二区三区四区| 成年人黄色毛片网站| 脱女人内裤的视频| 十八禁高潮呻吟视频| 夜夜爽天天搞| 91老司机精品| 日本五十路高清| 最近最新免费中文字幕在线| 亚洲情色 制服丝袜| 国产高清视频在线播放一区| 91成年电影在线观看| 久久久久国产一级毛片高清牌| 在线观看舔阴道视频| 久久精品熟女亚洲av麻豆精品| a在线观看视频网站| 99精品在免费线老司机午夜| 国产男靠女视频免费网站| 极品教师在线免费播放| 久久久国产一区二区| 欧美精品亚洲一区二区| 天堂动漫精品| 一进一出抽搐动态| 欧美性长视频在线观看| 亚洲成人国产一区在线观看| 欧美日韩福利视频一区二区| 无限看片的www在线观看| 另类亚洲欧美激情| 麻豆国产av国片精品| 久久精品亚洲精品国产色婷小说| 国产又爽黄色视频| 欧美成狂野欧美在线观看| 国产精品久久久久成人av| 热99久久久久精品小说推荐| 国产av又大| 看黄色毛片网站| 悠悠久久av| 成人亚洲精品一区在线观看| 色综合欧美亚洲国产小说| 国产视频一区二区在线看| 99在线人妻在线中文字幕 | 欧美日韩黄片免| 亚洲av日韩精品久久久久久密| 成熟少妇高潮喷水视频| 欧美日韩亚洲高清精品| 免费观看人在逋| 操美女的视频在线观看| 国产精品美女特级片免费视频播放器 | 日韩欧美三级三区| bbb黄色大片| 免费日韩欧美在线观看| 天堂动漫精品| 自拍欧美九色日韩亚洲蝌蚪91| 成人特级黄色片久久久久久久| 欧美日韩亚洲综合一区二区三区_| 成人影院久久| 老熟妇仑乱视频hdxx| 久久人妻熟女aⅴ| 50天的宝宝边吃奶边哭怎么回事| 十八禁高潮呻吟视频| 一级片免费观看大全| 色综合婷婷激情| 国产一卡二卡三卡精品| 日韩有码中文字幕| 国产精品1区2区在线观看. | 亚洲五月天丁香| 在线国产一区二区在线| av在线播放免费不卡| 制服诱惑二区| 国产在视频线精品| 99re6热这里在线精品视频| 十分钟在线观看高清视频www| 欧美成人午夜精品| 国内毛片毛片毛片毛片毛片| 国产片内射在线| 国产在线一区二区三区精| 超色免费av| av电影中文网址| 亚洲精品美女久久av网站| 91麻豆av在线| 三上悠亚av全集在线观看| 一级片免费观看大全| 色综合婷婷激情| 免费在线观看完整版高清| 久久精品国产99精品国产亚洲性色 | 9191精品国产免费久久| 在线观看免费视频网站a站| 三级毛片av免费| 天堂俺去俺来也www色官网| 91av网站免费观看| 啦啦啦在线免费观看视频4| 一a级毛片在线观看| aaaaa片日本免费| 亚洲精品成人av观看孕妇| 免费观看人在逋| 一级毛片女人18水好多| 亚洲片人在线观看| 国产麻豆69| 亚洲精品中文字幕在线视频| 久久99一区二区三区| 丝袜美足系列| 亚洲av美国av| 国产成人啪精品午夜网站| 国产精品影院久久| 国产av又大| 老司机在亚洲福利影院| 亚洲成av片中文字幕在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 好男人电影高清在线观看| 最新的欧美精品一区二区| 欧美国产精品va在线观看不卡| 亚洲熟妇熟女久久| 成人18禁在线播放| 国产日韩欧美亚洲二区| 国产熟女午夜一区二区三区| a级毛片黄视频| 日韩欧美免费精品| 啦啦啦在线免费观看视频4| 80岁老熟妇乱子伦牲交| 亚洲综合色网址| 婷婷精品国产亚洲av在线 | 水蜜桃什么品种好| 黑丝袜美女国产一区| 黄色怎么调成土黄色| 国产高清激情床上av| 国产成人欧美在线观看 | 日韩精品免费视频一区二区三区| 人人妻人人澡人人看| 不卡一级毛片| 亚洲精华国产精华精| 久久精品亚洲精品国产色婷小说| 黑丝袜美女国产一区| 午夜免费鲁丝| 国产精品免费一区二区三区在线 | 99国产精品99久久久久| 亚洲,欧美精品.| 最新在线观看一区二区三区| 三上悠亚av全集在线观看| 日韩大码丰满熟妇| 黄色成人免费大全| 丝袜美腿诱惑在线| 9191精品国产免费久久| 捣出白浆h1v1| 最新美女视频免费是黄的| 精品一品国产午夜福利视频| 国产亚洲欧美98| 欧美+亚洲+日韩+国产| 欧美日韩成人在线一区二区| 午夜激情av网站| 男女之事视频高清在线观看| 精品卡一卡二卡四卡免费| a级片在线免费高清观看视频| 亚洲avbb在线观看| 国产xxxxx性猛交| 精品久久蜜臀av无| 亚洲国产精品sss在线观看 | 人妻一区二区av| 中文字幕最新亚洲高清| 日韩 欧美 亚洲 中文字幕| 亚洲第一欧美日韩一区二区三区| 大型黄色视频在线免费观看| 99精品久久久久人妻精品| 妹子高潮喷水视频| a在线观看视频网站| 国产精品一区二区精品视频观看| 人成视频在线观看免费观看| 99国产精品一区二区三区| 国产成人精品在线电影| 真人做人爱边吃奶动态| 国产三级黄色录像| 亚洲 欧美一区二区三区| 高清视频免费观看一区二区| 中文字幕最新亚洲高清| 久久天堂一区二区三区四区| 黑人猛操日本美女一级片| 午夜福利乱码中文字幕| 中文字幕色久视频| 国产成人系列免费观看| 成年女人毛片免费观看观看9 | 欧美日本中文国产一区发布| 午夜免费鲁丝| 亚洲久久久国产精品| 国产精品一区二区免费欧美| 国产91精品成人一区二区三区| 男女高潮啪啪啪动态图| 少妇被粗大的猛进出69影院| 日韩欧美在线二视频 | 操美女的视频在线观看| 极品少妇高潮喷水抽搐| 女人精品久久久久毛片| 人人妻,人人澡人人爽秒播| 精品久久久久久电影网| 99久久人妻综合| 日本vs欧美在线观看视频| 99国产极品粉嫩在线观看| 性少妇av在线| av中文乱码字幕在线| av片东京热男人的天堂| av线在线观看网站| 精品久久久久久电影网| 波多野结衣av一区二区av| 国产成人免费观看mmmm| 国产男靠女视频免费网站| 亚洲免费av在线视频| 欧美精品高潮呻吟av久久| 777米奇影视久久| 日韩精品免费视频一区二区三区| 法律面前人人平等表现在哪些方面| av欧美777| 中文字幕最新亚洲高清| 99热只有精品国产| 一级毛片女人18水好多|