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

    基于非線性最優(yōu)反饋控制的遠(yuǎn)程攔截閉環(huán)制導(dǎo)*

    2016-10-18 10:43:08閆循良廖守億王仕成
    火力與指揮控制 2016年9期
    關(guān)鍵詞:偽譜開(kāi)環(huán)制導(dǎo)

    閆循良,廖守億,王仕成

    (火箭軍工程大學(xué),西安710025)

    基于非線性最優(yōu)反饋控制的遠(yuǎn)程攔截閉環(huán)制導(dǎo)*

    閆循良,廖守億,王仕成

    (火箭軍工程大學(xué),西安710025)

    利用非線性最優(yōu)反饋控制和偽譜軌跡快速重構(gòu)技術(shù),設(shè)計(jì)了一種有限推力空間遠(yuǎn)程攔截自適應(yīng)閉環(huán)制導(dǎo)方法。首先建立了有限推力遠(yuǎn)程攔截最優(yōu)制導(dǎo)問(wèn)題模型,并給出了適用于該問(wèn)題的非線性最優(yōu)反饋控制求解原理。然后,將空間變軌動(dòng)力學(xué)模型特點(diǎn)和偽譜法相結(jié)合,設(shè)計(jì)了基于狀態(tài)量縮減的計(jì)算效率改進(jìn)策略以提高軌跡優(yōu)化的實(shí)時(shí)性?;诟倪M(jìn)偽譜法進(jìn)行連續(xù)軌跡快速重構(gòu),利用開(kāi)環(huán)最優(yōu)解形成閉環(huán)反饋,從而保證制導(dǎo)指令的實(shí)時(shí)更新,并通過(guò)引入控制邏輯對(duì)制導(dǎo)算法進(jìn)行改進(jìn)。時(shí)間最短遠(yuǎn)程攔截仿真表明,該閉環(huán)制導(dǎo)方法在保證任務(wù)指標(biāo)最優(yōu)性的同時(shí),可以有效抑制J2攝動(dòng)和計(jì)算誤差的影響,具有較高的制導(dǎo)精度、自適應(yīng)性和魯棒性。

    遠(yuǎn)程攔截,閉環(huán)制導(dǎo),非線性最優(yōu)反饋控制,改進(jìn)偽譜法,軌跡重構(gòu)

    0 引言

    空間動(dòng)能攔截遠(yuǎn)程變軌具有空間包絡(luò)大、耗時(shí)長(zhǎng)、限制條件和干擾多等典型特點(diǎn),該段制導(dǎo)規(guī)律是實(shí)現(xiàn)精確制導(dǎo)交班并保證攔截任務(wù)順利完成的關(guān)鍵技術(shù)之一。因此,進(jìn)行空間遠(yuǎn)程攔截制導(dǎo)理論和算法的研究至關(guān)重要。

    目前研究較為成熟的遠(yuǎn)程攔截制導(dǎo)律多為非最優(yōu)制導(dǎo)律,如比例導(dǎo)引、速度增益制導(dǎo)、基于零控?cái)r截流形的制導(dǎo)等[2-6]。上述制導(dǎo)律往往基于簡(jiǎn)化的變軌動(dòng)力學(xué)模型,考慮最優(yōu)性要求較少,同時(shí)對(duì)發(fā)動(dòng)機(jī)性能及變軌時(shí)間往往具有一定的限制。為了提高飛行任務(wù)的可行性、最優(yōu)性和安全性,最優(yōu)控制理論被引入到遠(yuǎn)程軌道設(shè)計(jì)領(lǐng)域并得到廣泛應(yīng)用。一般而言,均將空間軌道設(shè)計(jì)問(wèn)題轉(zhuǎn)化為最優(yōu)控制問(wèn)題進(jìn)行處理,并采用間接法或直接法進(jìn)行數(shù)值求解[1,7-8]。然而,此類制導(dǎo)規(guī)律為開(kāi)環(huán)形式,考慮實(shí)際飛行中的不確定性因素和干擾時(shí),存在較大的誤差。近年來(lái),隨著航空航天技術(shù)、先進(jìn)控制理論和計(jì)算技術(shù)的發(fā)展,尤其是高速微型計(jì)算機(jī)的發(fā)展,最優(yōu)控制領(lǐng)域涌現(xiàn)出了新的非傳統(tǒng)、非線性反饋控制技術(shù),如Ross等[9-10]提出的偽譜最優(yōu)反饋控制算法。該算法借鑒采樣保持反饋控制的基本原理,利用偽譜法的高精度和快速計(jì)算優(yōu)勢(shì),實(shí)現(xiàn)了利用低頻自由采樣完成高精度實(shí)時(shí)反饋控制的目的,同時(shí)避免了復(fù)雜偏微分方程的求解。鑒于上述優(yōu)點(diǎn),該方法在諸多制導(dǎo)、控制問(wèn)題中得到了應(yīng)用[11-13]。

    本文基于非線性最優(yōu)反饋控制和偽譜軌跡快速重構(gòu)技術(shù),設(shè)計(jì)了一種空間遠(yuǎn)程攔截閉環(huán)制導(dǎo)方法。利用連續(xù)軌跡快速重構(gòu)實(shí)時(shí)反饋更新制導(dǎo)指令,并引入制導(dǎo)邏輯對(duì)算法進(jìn)行改進(jìn)。通過(guò)有限推力時(shí)間最短遠(yuǎn)程攔截?cái)?shù)值仿真驗(yàn)證了其可行性和有效性。

    1 有限推力遠(yuǎn)程攔截最優(yōu)制導(dǎo)問(wèn)題

    1.1遠(yuǎn)程變軌動(dòng)力學(xué)模型

    在平方反比引力場(chǎng)和有限推力假設(shè)下,以赤道慣性系中的位置、速度為狀態(tài)量的航天器變軌攝動(dòng)動(dòng)力學(xué)方程可描述為

    將式(1)進(jìn)行無(wú)量綱化處理,得到

    1.2最優(yōu)制導(dǎo)問(wèn)題描述

    遠(yuǎn)程攔截最優(yōu)制導(dǎo)規(guī)律設(shè)計(jì)問(wèn)題可描述為:確定最優(yōu)推力大小T、方向u的變化規(guī)律,以及其他系統(tǒng)參數(shù),使得航天器在給定的端點(diǎn)條件和約束下,以一定的最優(yōu)性指標(biāo)和精度完成預(yù)定攔截任務(wù),同時(shí)保證對(duì)不確定性因素和干擾具有較好的魯棒性和自適應(yīng)性。

    對(duì)于空間快速打擊等任務(wù)而言,時(shí)間最短是首要考慮的性能指標(biāo)??紤]時(shí)間最短攔截問(wèn)題,有

    除需要滿足動(dòng)力學(xué)約束、狀態(tài)及控制約束及端點(diǎn)約束條件外,考慮到推力方向矢量u為一單位矢量,因此,需要滿足以下內(nèi)點(diǎn)路徑約束

    若考慮變軌高度限制,還需滿足不等式約束

    顯然,該制導(dǎo)問(wèn)題為一典型非線性、多約束、包含干擾和不確定因素的復(fù)雜最優(yōu)控制問(wèn)題。

    基于最優(yōu)反饋控制原理,可以對(duì)上述問(wèn)題進(jìn)行數(shù)值求解。即利用采樣時(shí)間間隔對(duì)控制時(shí)間區(qū)間[t0,tf]分段,并順序求?。踭i,ti+1]上的最優(yōu)開(kāi)環(huán)控制解(即最優(yōu)反饋解)u(t)=k(t,x(ti)),依次執(zhí)行開(kāi)環(huán)控制即等價(jià)于閉環(huán)反饋控制。關(guān)于解的穩(wěn)定性、連續(xù)性等細(xì)節(jié)討論,可參考相關(guān)文獻(xiàn)[9-10,13]。

    2 基于改進(jìn)偽譜反饋控制的閉環(huán)制導(dǎo)策略設(shè)計(jì)

    考慮到最優(yōu)反饋控制的實(shí)現(xiàn)取決于開(kāi)環(huán)最優(yōu)控制的實(shí)時(shí)獲取能力,以下結(jié)合空間變軌動(dòng)力學(xué)模型特點(diǎn)和偽譜法,設(shè)計(jì)了基于狀態(tài)量縮減的計(jì)算效率改進(jìn)策略以提高開(kāi)環(huán)控制解算的實(shí)時(shí)性。

    2.1基于狀態(tài)量縮減的改進(jìn)偽譜法

    近年來(lái),偽譜法(PseudospectralMethod,PSM)的大量研究表明,PSM具有提供實(shí)時(shí)或近實(shí)時(shí)最優(yōu)控制解的巨大潛力[8-9,14]。Legendre-Gauss-Lobatto(LGL)PSM是目前應(yīng)用最為廣泛的偽譜方法之一,具有諸多性能優(yōu)勢(shì),在軌跡優(yōu)化研究中得到了廣泛應(yīng)用[14-15]。為提高計(jì)算效率,本節(jié)結(jié)合空間變軌的動(dòng)力學(xué)特性,提出利用二階微分方程技術(shù)消去部分狀態(tài)量,實(shí)現(xiàn)優(yōu)化變量數(shù)目的縮減,以進(jìn)一步提高傳統(tǒng)LGLPSM數(shù)值解算的實(shí)時(shí)性。

    將式(3)改寫(xiě)為以下微分方程組的形式

    采用LGLPSM對(duì)上述方程組進(jìn)行離散化,得到包含微分矩陣D的代數(shù)方程組,為

    式中,D=[Dki]為(N+1)×(N+1)維的微分矩陣。R、V、A、M、T分別為相應(yīng)變量在LGL點(diǎn)的離散值。故優(yōu)化變量為X=[RTVTMT]T。利用成熟的非線性規(guī)劃算法即可實(shí)現(xiàn)數(shù)值求解,并可對(duì)解的可行性和最優(yōu)性進(jìn)行驗(yàn)證。

    考慮式(9)的前兩式,將其變形為

    進(jìn)而得到新的離散化代數(shù)方程為

    其中

    結(jié)合離散化的質(zhì)量方程,有新的優(yōu)化變量為X'=[RTMT]T,同時(shí)得到離散化的動(dòng)力學(xué)方程約束為

    并有

    經(jīng)過(guò)變換處理后,優(yōu)化變量由X縮減為X',即消去了狀態(tài)分量V;動(dòng)力學(xué)約束由式(9)變?yōu)槭剑?3),消去了與V有關(guān)的矢量方程。考慮到V為速度矢量v在LGL點(diǎn)的離散值,其維數(shù)為3(N+1),故優(yōu)化變量和約束數(shù)目均縮減了3(N+1)。因此,狀態(tài)量縮減后,數(shù)值解算效率將會(huì)明顯提高,算法實(shí)時(shí)性會(huì)得到進(jìn)一步保證。

    2.2自適應(yīng)閉環(huán)制導(dǎo)策略設(shè)計(jì)

    圖1給出了基于軌跡快速重構(gòu)的遠(yuǎn)程變軌閉環(huán)制導(dǎo)策略圖。該策略主要包括離線開(kāi)環(huán)軌跡生成、數(shù)據(jù)存儲(chǔ)及處理、遠(yuǎn)程變軌動(dòng)力學(xué)方程、制導(dǎo)邏輯、軌跡快速重構(gòu)等5個(gè)模塊。

    圖1 基于軌跡快速重構(gòu)的遠(yuǎn)程變軌閉環(huán)制導(dǎo)策略

    制導(dǎo)算法流程具體可描述為:

    ①給定初、末端條件x0、xf,指標(biāo)J及約束C,利用PSM離線求解開(kāi)環(huán)最優(yōu)控制解(u0*(t),x0*(t))及tf*等,并將其作為啟動(dòng)閉環(huán)制導(dǎo)的初始猜測(cè);初始化參數(shù),,積分步長(zhǎng)tp;

    ②將ui*映射至積分步長(zhǎng)對(duì)應(yīng)的網(wǎng)格,從而獲得積分所需的制導(dǎo)指令,由ti時(shí)刻起積分運(yùn)動(dòng)方程,通過(guò)模擬飛行以獲得實(shí)時(shí)飛行軌跡xiR;

    ③執(zhí)行制導(dǎo)算法控制邏輯模塊:

    設(shè)當(dāng)前時(shí)刻狀態(tài)偏差

    終端狀態(tài)偏差

    其中,a,b,m,n均為權(quán)重系數(shù),可根據(jù)任務(wù)要求選擇;tc為當(dāng)前時(shí)刻,分別為積分得到的實(shí)際地心距和速度矢量,分別為前一次軌跡重構(gòu)得到的理論最優(yōu)狀態(tài)和速度矢量。

    令ε1和ε2為與設(shè)計(jì)要求相關(guān)的控制門限值(可根據(jù)不同的任務(wù)要求進(jìn)行選取,本文不作詳細(xì)討論)。若εf≤ε1,表明實(shí)際終端狀態(tài)與期望終端狀態(tài)滿足偏差要求,故終止軌跡重構(gòu),采用當(dāng)前最優(yōu)控制解完成后續(xù)遠(yuǎn)程變軌機(jī)動(dòng),制導(dǎo)過(guò)程結(jié)束;εf>ε1時(shí),若εc≤ε2,表明當(dāng)前實(shí)際狀態(tài)與前一次最優(yōu)解狀態(tài)滿足偏差要求,則以一定步長(zhǎng)繼續(xù)積分運(yùn)動(dòng)方程而不進(jìn)行軌跡重構(gòu),即轉(zhuǎn)到②;否則,執(zhí)行開(kāi)環(huán)最優(yōu)解的計(jì)算,轉(zhuǎn)至④;

    ④積分運(yùn)動(dòng)方程的同時(shí),以ti時(shí)刻的狀態(tài)(ti)作為初始條件,以等參數(shù)作為初始猜測(cè),重新計(jì)算該問(wèn)題的最優(yōu)開(kāi)環(huán)控制解;根據(jù)前述定義,知優(yōu)化計(jì)算耗時(shí),故得優(yōu)化完成時(shí)刻新的最優(yōu)反饋解及參數(shù)

    ⑤在優(yōu)化完成時(shí)刻ti+1,實(shí)現(xiàn)控制量的反饋更新,即令;令i=i+1,返回②。

    針對(duì)上述過(guò)程作以下討論:①此處,反饋更新時(shí)間即為制導(dǎo)指令更新周期,該周期可隨開(kāi)環(huán)控制的解算自適應(yīng)變化。由于采用前一步的優(yōu)化結(jié)果作為軌跡重構(gòu)猜測(cè)值,原則上,若外界干擾引起的狀態(tài)偏差較小,軌跡重構(gòu)耗時(shí)會(huì)較小,反饋更新時(shí)間較小,反之,則較大;②若當(dāng)前軌跡重構(gòu)未獲得可行解,則可繼續(xù)使用上一步得到的制導(dǎo)指令,以一定周期積分運(yùn)動(dòng)方程,并采用新的飛行狀態(tài)作為初值重新計(jì)算。

    3 仿真算例及分析

    3.1仿真條件

    以近地軌道時(shí)間最短遠(yuǎn)程攔截為例,驗(yàn)證所提算法的性能。追蹤器和目標(biāo)的初始軌道要素見(jiàn)表1。

    表1 追蹤器與目標(biāo)初始軌道要素

    追蹤器初始質(zhì)量m0=1 000 kg,燃料質(zhì)量mp= 400 kg,發(fā)動(dòng)機(jī)參數(shù)為:比沖Is=300 s,推力上限Tmax=500N。仿真計(jì)算在CPU為E7500/2.93GHz,2G內(nèi)存的臺(tái)式計(jì)算機(jī)上進(jìn)行,基于MATLAB環(huán)境進(jìn)行編程。為了提高計(jì)算效率,離線偽譜優(yōu)化采用初值生成技術(shù)和串行優(yōu)化策略[8]。初值生成器的節(jié)點(diǎn)數(shù)目取N1=5,離線最優(yōu)解的多節(jié)點(diǎn)數(shù)目及在線重構(gòu)節(jié)點(diǎn)數(shù)目均取N2=25。

    針對(duì)不考慮攝動(dòng)影響的時(shí)間最短遠(yuǎn)程攔截問(wèn)題,利用改進(jìn)的LGLPSM分別進(jìn)行多次離線軌跡優(yōu)化仿真,平均CPU耗時(shí)為2.61 s,優(yōu)化指標(biāo)J=tf= 2 241.72 s,發(fā)動(dòng)機(jī)全程保持最大推力工作,攔截成功時(shí)刻消耗燃料380.46 kg。部分離線優(yōu)化結(jié)果如圖2所示??梢?jiàn),該方法能以較少的節(jié)點(diǎn)和計(jì)算時(shí)間獲得優(yōu)化軌跡以及光滑控制量。圖2(b)、圖2(c)給出了解的最優(yōu)性驗(yàn)證結(jié)果。圖2(b)中,OPT代表數(shù)值優(yōu)化得到的控制解,PMP代表基于龐特里亞金極值原理(PMP)和協(xié)狀態(tài)映射原理(CMP)獲得的控制解,二者吻合良好,數(shù)值結(jié)果滿足PMP的一階最優(yōu)必要條件。圖2(c)給出了由PMP和CMP得到的Hamiltonian函數(shù),其值在0附近震蕩。事實(shí)上,由PMP可知,該問(wèn)題的Hamiltonian函數(shù)沿最優(yōu)軌跡恒為0,因此,數(shù)值解滿足該一階最優(yōu)必要條件。

    圖2 軌跡離線優(yōu)化結(jié)果

    為了驗(yàn)證本文制導(dǎo)算法的精度、魯棒性等性能,考慮地球扁率J2攝動(dòng),分別進(jìn)行開(kāi)環(huán)制導(dǎo)(Open-loop optimal guidance,OLOG)及閉環(huán)制導(dǎo)(Closed-loop optimalguidance,CLOG)仿真。仿真中,閉環(huán)制導(dǎo)邏輯中的權(quán)重和門限值選取為:a,b,m均取1,n取0,ε1=50/rref,ε2=1×10-6,其中,參考距離rref= re為地球平均半徑,即分別考慮當(dāng)前時(shí)刻狀態(tài)偏差和攔截終端位置偏差要求,終端位置偏差不大于50m。表2給出了不同情況下的制導(dǎo)仿真部分?jǐn)?shù)據(jù)結(jié)果,其中,OLOG和CLOG分別代表不考慮J2攝動(dòng)的開(kāi)環(huán)和閉環(huán)制導(dǎo)結(jié)果,“OLOG+J2”和“CLOG+J2”分別代表考慮J2攝動(dòng)的開(kāi)環(huán)和閉環(huán)制導(dǎo)仿真結(jié)果,Tg為閉環(huán)制導(dǎo)連續(xù)軌跡重構(gòu)結(jié)束的對(duì)應(yīng)時(shí)刻,該時(shí)刻表征了制導(dǎo)計(jì)算結(jié)束的時(shí)間,Δrf、Δvf、mf分別為飛行結(jié)束tf時(shí)刻的位置、速度偏差及剩余質(zhì)量。

    表2 不同情況的仿真結(jié)果

    分析表2可知:①考慮J2攝動(dòng)影響時(shí),開(kāi)環(huán)制導(dǎo)具有較大的制導(dǎo)偏差,而閉環(huán)制導(dǎo)策略能夠顯著提高攔截精度Δrf;各種情況的終端相對(duì)攔截速度Δvf均較大,達(dá)到了近1 400m/s,可以較好地實(shí)現(xiàn)直接碰撞攔截;②閉環(huán)制導(dǎo)結(jié)束時(shí)刻均小于實(shí)際飛行時(shí)間,即不需要全程進(jìn)行軌跡重構(gòu),便可達(dá)到預(yù)定的制導(dǎo)精度;CLOG+J2對(duì)應(yīng)的Tg較大,表明J2攝動(dòng)的存在使得飛行誤差增大,延長(zhǎng)了飛行器軌跡重構(gòu)和反饋制導(dǎo)時(shí)間;③與其他3種情況相比,CLOG+J2對(duì)應(yīng)的mf和tf分別增加了0.77 kg和70.87 s;事實(shí)上,OLOG和OLOG+J2只是利用了離線軌跡優(yōu)化得到的最優(yōu)控制指令,故消耗的燃料和攔截時(shí)間是一致的,CLOG雖然對(duì)制導(dǎo)指令進(jìn)行了調(diào)整,但仍滿足時(shí)間最短優(yōu)化指標(biāo)的結(jié)果;而為了消除J2攝動(dòng)的影響,OLOG+J2需要進(jìn)行連續(xù)軌跡重構(gòu)以生成新的制導(dǎo)指令,從而增加了飛行時(shí)間,但由于發(fā)動(dòng)機(jī)并不是全程工作,故剩余燃料有一定的增加;④OLOG對(duì)應(yīng)的攔截精度偏差為98.01m,這在一定意義上驗(yàn)證了離線軌跡優(yōu)化解的可行性;事實(shí)上,該偏差是由于離線優(yōu)化節(jié)點(diǎn)數(shù)目較少和計(jì)算誤差所引起的,可以通過(guò)增大節(jié)點(diǎn)數(shù)目和改進(jìn)積分策略減小該誤差,或通過(guò)本文的連續(xù)軌跡重構(gòu)和制導(dǎo)指令更新來(lái)減小該誤差,如對(duì)于CLOG而言,其攔截精度偏差即減小為49.47m。

    圖3給出了不同情況對(duì)應(yīng)的狀態(tài)軌跡、制導(dǎo)指令與離線優(yōu)化解(PSM)的部分對(duì)比結(jié)果??梢?jiàn),由于計(jì)算誤差和J2攝動(dòng)等影響因素的存在,閉環(huán)制導(dǎo)對(duì)應(yīng)的飛行軌跡、制導(dǎo)信號(hào)與離線最優(yōu)解并不相同,但最終逐漸收斂于預(yù)定的指標(biāo)和約束條件,滿足攔截精度要求。閉環(huán)制導(dǎo)信號(hào)在絕大部分時(shí)間內(nèi)仍能夠保持較好的光滑性,對(duì)于不光滑的部分,則可以通過(guò)引入控制量變化率約束加以改進(jìn)。

    對(duì)于以上仿真而言,CLOG僅需要5次軌跡重構(gòu)便可滿足終端攔截要求,重構(gòu)耗時(shí)最大值為0.17 s,均值為0.08 s;由于J2攝動(dòng)是連續(xù)作用的,CLOG+J2進(jìn)行了13 033次軌跡重構(gòu),重構(gòu)耗時(shí)最大值為2.05 s,均值為0.12 s。可見(jiàn),本文的計(jì)算效率改進(jìn)策略可以顯著降低優(yōu)化的CPU耗時(shí)(即制導(dǎo)指令更新周期),保證制導(dǎo)策略在線應(yīng)用的可行性,同時(shí)實(shí)現(xiàn)了制導(dǎo)更新周期隨外界干擾的自適應(yīng)調(diào)整。

    圖4給出了CLOG+J2情況對(duì)應(yīng)的狀態(tài)偏差量變化曲線,其橫坐標(biāo)Iter表示迭代次數(shù)??梢?jiàn),當(dāng)前狀態(tài)偏差εc的變化沒(méi)有具體規(guī)律可循,但是在整個(gè)過(guò)程中多次出現(xiàn)滿足控制門限值的情況,此時(shí),飛行器以一定步長(zhǎng)繼續(xù)飛行而不進(jìn)行制導(dǎo)指令的計(jì)算和更新,故而可以減少軌跡重構(gòu)次數(shù);在迭代計(jì)算的最后階段,εc達(dá)到并保持控制門限值要求。而終端狀態(tài)偏差εf總體呈現(xiàn)先增大后減小的變化趨勢(shì),一方面是因?yàn)棣與存在滿足控制門限的情況,故制導(dǎo)指令停止更新導(dǎo)致終端偏差增大,另一方面是由制導(dǎo)指令連續(xù)更新和計(jì)算偏差等因素造成的。但隨著迭代計(jì)算的進(jìn)行,εf最終滿足終端控制門限(該時(shí)間對(duì)應(yīng)制導(dǎo)計(jì)算結(jié)束時(shí)間,此處為2 154.57 s),此時(shí),飛行器不再反饋更新制導(dǎo)信號(hào),即僅采用上一次在線優(yōu)化得到的制導(dǎo)指令,便可完成交會(huì)任務(wù)。

    圖3 不同情況的制導(dǎo)仿真結(jié)果對(duì)比

    圖4 CLOG+J2對(duì)應(yīng)的狀態(tài)偏差曲線

    4 結(jié)論

    本文基于最優(yōu)反饋控制原理和軌跡快速重構(gòu)技術(shù),設(shè)計(jì)了一種空間遠(yuǎn)程攔截自適應(yīng)最優(yōu)制導(dǎo)策略,并進(jìn)行了數(shù)值仿真驗(yàn)證。結(jié)果表明,利用軌跡重構(gòu)實(shí)時(shí)反饋更新制導(dǎo)指令,無(wú)需傳統(tǒng)的參考軌跡跟蹤和控制增益調(diào)整過(guò)程,即可對(duì)J2攝動(dòng)影響和計(jì)算誤差進(jìn)行修正,制導(dǎo)策略具有較好的自主性和魯棒性;基于最優(yōu)實(shí)時(shí)反饋和制導(dǎo)邏輯,可實(shí)現(xiàn)制導(dǎo)指令更新周期的自適應(yīng)調(diào)整,且不需要更新周期和在線計(jì)算時(shí)間的先驗(yàn)信息,因此,算法具有較強(qiáng)的自適應(yīng)能力;同時(shí),偽譜優(yōu)化和基于狀態(tài)量縮減的計(jì)算效率改進(jìn)策略提高了軌跡重構(gòu)的實(shí)時(shí)性能,為算法的在線應(yīng)用提供了一定的保證。

    [1]湯善同.飛行器的軌道優(yōu)化與制導(dǎo)規(guī)律研究綜述[J].現(xiàn)代防御技術(shù),2002,30(1):32-37.

    [2]ZARCHAN P.Midcouse guidance strategies for exoatmospheric intercept[R].AD-A355820,1998.

    [3]陳峰,肖業(yè)倫,陳萬(wàn)春.基于需用速度增益曲面的大氣層外超遠(yuǎn)程攔截導(dǎo)引方法[J].航空學(xué)報(bào),2010,31(2):342-349.

    [4]NEWMAN B.Strategic midcourse guidance using modified zero effortmiss steering[J].Journal of Guidance,Control,and Dynamics,1996,19(1):107-112.

    [5]尚海濱,崔平遠(yuǎn),欒恩杰.基于最優(yōu)狀態(tài)反饋的小推力轉(zhuǎn)移軌道制導(dǎo)策略[J].吉林大學(xué)學(xué)報(bào),2007,37(4):949-954.

    [6]朱戰(zhàn)霞,金偉,韓沛.一種基于引力場(chǎng)線性化的制導(dǎo)方法[J].西北工業(yè)大學(xué)學(xué)報(bào),2012,30(2):181-186.

    [7]GAO Y,LIW Q.Systematic direct approach for optimizing continuous-thrustearth-orbit transfers[J].Chinese Journal ofAeronautics,2009,22(5):56-69.

    [8]閆循良,廖守億,張金生,等.基于節(jié)點(diǎn)改進(jìn)策略的偽譜軌跡優(yōu)化[J].宇航學(xué)報(bào),2013,34(7):891-900.

    [9]ROSS IM,SEKHAVATP,F(xiàn)LEMINGA,etal.Optimal feedback control:foundations,examples,and experimental results for a new approach[J].Journal of Guidance,Control,and Dynamics,2008,31(2):307-321.

    [10]ROSS IM,GONGQ,F(xiàn)AHROO F,etal.Practical stabilization through real-time optimal control[C]//Proceedings of the 2006 American Control Conference,Minneapolis,MN,June14-16,2006.

    [11]YAN H,ROSS IM,ALFRIEND K T.Pseudospectral feedback control for three-axismagnetic attitude stabilization in elliptic orbits[J].Journal of Guidance,Control,and Dynamics,2007,30(4):1107-1115.

    [12]BOLLINO K P,ROSS I M.A pseudospectral feedback method for real-time optimal guidance of reentry vehicles[C]//Proceedings of the 2007 American Control Conference,New York:IEEE,2007,7:3861-3867.

    [13]張友安,王麗英,趙國(guó)榮.基于偽譜法的自由采樣實(shí)時(shí)最優(yōu)反饋控制及應(yīng)用[J].控制理論與應(yīng)用,2012,29(9):1151-1156.

    [14]GARG D,PATTERSON M A,HAGERW W.An overview of three pseudospectralmethods for thenumericalsolution of optimalproblems[C]//AAS09-332.

    [15]ROSS IM,GONGQ,SEKHAVATP.Low-thrust,high-accuracy trajectory optimization[J].Journal of Guidance,Control,and Dynamics,2007,30(4):921-933.

    Nonlinear Optimal Feedback ControlBased Closed-Loop Guidanceof Long-range Intercept

    YANXun-liang,LIAOShou-yi,WANGShi-cheng
    (Rocket Force University of Engineering,Xi’an 710025,China)

    An adaptive closed-loop guidance algorithm,which is based on nonlinear optimal feedback and rapid pseudospectral-based trajectory reconfiguration,is studied and proposed for longrange orbital interceptwith finite thrust.Themathematic model is established for the optimal guidance problem,which can be solved through nonlinear optimal feedback control theory.Then,combining the characteristics of the dynamics model with the pseudospectral method,a computational efficiency refinementmethod with the reduced state parameters,is proposed to improve the efficiency of trajectory optimization.The modified pseudospectral method is used for the successive rapid trajectory reconfiguration such that the open-loop optimal controls are fed back to the system to update the guidance commands in real-time.A control logic algorithm is introduced to heighten the performance of the algorithm.The simulation results of long-range intercept indicate that the algorithm can not only supply excellent optimal solution for practical mission with high accuracy,but can manage J2 perturbation and computation error autonomously,adaptively and robustly.

    long-range intercept,closed-loop guidance,nonlinear optimal feedback control,modified pseudospectralmethod,trajectory reconfiguration

    V421.4

    A

    1002-0640(2016)09-0084-06

    2015-07-09

    2015-08-17

    中國(guó)博士后科學(xué)基金資助項(xiàng)目(2014T70974)

    閆循良(1984-),男,山東濟(jì)寧人,博士后。研究方向:飛行動(dòng)力學(xué)與控制、軌跡設(shè)計(jì)與優(yōu)化。

    猜你喜歡
    偽譜開(kāi)環(huán)制導(dǎo)
    矩陣偽譜的新定位集及其在土壤生態(tài)系統(tǒng)的應(yīng)用
    轉(zhuǎn)速開(kāi)環(huán)恒壓頻比的交流調(diào)速系統(tǒng)的分析與仿真研究
    紊流環(huán)境下四維軌跡優(yōu)化的偽譜方法研究
    基于MPSC和CPN制導(dǎo)方法的協(xié)同制導(dǎo)律
    基于在線軌跡迭代的自適應(yīng)再入制導(dǎo)
    帶有攻擊角約束的無(wú)抖振滑模制導(dǎo)律設(shè)計(jì)
    一種溫和環(huán)醚開(kāi)環(huán)成雙酯的新方法
    基于開(kāi)環(huán)補(bǔ)償?shù)娘w機(jī)偏航角控制系統(tǒng)設(shè)計(jì)及仿真
    偽譜法及其在飛行器軌跡優(yōu)化設(shè)計(jì)領(lǐng)域的應(yīng)用綜述*
    基于Radau偽譜方法的軌跡優(yōu)化
    小说图片视频综合网站| 人人妻,人人澡人人爽秒播| 国产乱人偷精品视频| 一级毛片电影观看 | 久久久久久久亚洲中文字幕| 亚洲人成网站在线观看播放| 全区人妻精品视频| 国产 一区精品| 淫秽高清视频在线观看| 国产中年淑女户外野战色| 亚洲在线观看片| 日韩精品中文字幕看吧| 看非洲黑人一级黄片| 18禁在线播放成人免费| 性欧美人与动物交配| 小蜜桃在线观看免费完整版高清| 国产精品久久久久久亚洲av鲁大| a级毛片免费高清观看在线播放| 最后的刺客免费高清国语| 亚洲av电影不卡..在线观看| 久久久久久久久久成人| 国产精品无大码| 中文字幕熟女人妻在线| 精品一区二区三区视频在线| 99精品在免费线老司机午夜| ponron亚洲| 97超视频在线观看视频| 中国美女看黄片| 精品日产1卡2卡| 久久久久免费精品人妻一区二区| 亚洲乱码一区二区免费版| 国产欧美日韩精品亚洲av| 插阴视频在线观看视频| 又黄又爽又刺激的免费视频.| 小说图片视频综合网站| 国产精品人妻久久久影院| 国模一区二区三区四区视频| 久久久a久久爽久久v久久| 免费av不卡在线播放| 色噜噜av男人的天堂激情| 日韩成人av中文字幕在线观看 | 精品久久久久久久久亚洲| 俄罗斯特黄特色一大片| 国产日本99.免费观看| 午夜视频国产福利| 亚洲真实伦在线观看| 色哟哟·www| 国产淫片久久久久久久久| 精品久久久久久久久久免费视频| 悠悠久久av| 国产高清视频在线观看网站| 能在线免费观看的黄片| a级一级毛片免费在线观看| 午夜福利在线观看免费完整高清在 | 午夜视频国产福利| 国产一区二区三区av在线 | 麻豆av噜噜一区二区三区| 亚洲人成网站在线观看播放| 日韩欧美精品v在线| 国产精品精品国产色婷婷| 精品无人区乱码1区二区| av专区在线播放| 成人av一区二区三区在线看| 女生性感内裤真人,穿戴方法视频| 美女cb高潮喷水在线观看| 成年免费大片在线观看| 又黄又爽又刺激的免费视频.| 精品熟女少妇av免费看| 亚洲内射少妇av| 麻豆国产av国片精品| 成人综合一区亚洲| 欧美又色又爽又黄视频| 嫩草影院新地址| 久久精品国产99精品国产亚洲性色| 青春草视频在线免费观看| 麻豆久久精品国产亚洲av| 十八禁网站免费在线| 国内精品宾馆在线| 赤兔流量卡办理| 亚洲av中文av极速乱| a级毛片a级免费在线| 国产欧美日韩一区二区精品| 小蜜桃在线观看免费完整版高清| 美女免费视频网站| а√天堂www在线а√下载| 国产v大片淫在线免费观看| 伦精品一区二区三区| 无遮挡黄片免费观看| 天美传媒精品一区二区| 99热这里只有精品一区| 精品一区二区三区视频在线观看免费| 久久久精品94久久精品| 最近手机中文字幕大全| 可以在线观看的亚洲视频| av在线播放精品| 亚洲av中文av极速乱| 国产成人freesex在线 | 亚洲,欧美,日韩| 中国美白少妇内射xxxbb| 国产大屁股一区二区在线视频| 内地一区二区视频在线| av卡一久久| 少妇人妻精品综合一区二区 | 看免费成人av毛片| 成人国产麻豆网| 一区二区三区免费毛片| 女人被狂操c到高潮| 狂野欧美白嫩少妇大欣赏| 国产免费一级a男人的天堂| 中文字幕久久专区| 中文字幕人妻熟人妻熟丝袜美| 国产高清三级在线| 国产欧美日韩精品一区二区| 久久久久久久久中文| 国产亚洲精品久久久久久毛片| 我的女老师完整版在线观看| 免费在线观看成人毛片| 国产亚洲精品av在线| 91久久精品国产一区二区成人| 午夜老司机福利剧场| 毛片女人毛片| 精品午夜福利视频在线观看一区| 一级毛片aaaaaa免费看小| 97超碰精品成人国产| 在线国产一区二区在线| 国产aⅴ精品一区二区三区波| 久久久久精品国产欧美久久久| 99热这里只有精品一区| 99精品在免费线老司机午夜| 国产精品国产高清国产av| av在线观看视频网站免费| 变态另类丝袜制服| 午夜亚洲福利在线播放| 亚洲最大成人中文| 国产三级中文精品| 国产精品福利在线免费观看| АⅤ资源中文在线天堂| 一级毛片我不卡| 成人av一区二区三区在线看| 国产av麻豆久久久久久久| 国产激情偷乱视频一区二区| 色视频www国产| 色哟哟哟哟哟哟| 赤兔流量卡办理| 大香蕉久久网| 97碰自拍视频| 亚洲精品色激情综合| 深爱激情五月婷婷| 欧美色欧美亚洲另类二区| 干丝袜人妻中文字幕| 变态另类成人亚洲欧美熟女| 亚洲一区高清亚洲精品| 免费看a级黄色片| h日本视频在线播放| 男女之事视频高清在线观看| 亚洲精品国产av成人精品 | 天天一区二区日本电影三级| 欧美绝顶高潮抽搐喷水| 欧美不卡视频在线免费观看| 国产精品国产高清国产av| 成年免费大片在线观看| 麻豆乱淫一区二区| 欧美最黄视频在线播放免费| 最近最新中文字幕大全电影3| 日日摸夜夜添夜夜添av毛片| 国产一级毛片七仙女欲春2| 日本与韩国留学比较| 一级a爱片免费观看的视频| 久久精品影院6| 久久久久久大精品| 午夜福利视频1000在线观看| 亚洲图色成人| 91在线观看av| 亚洲精品日韩在线中文字幕 | 亚洲av免费在线观看| 午夜老司机福利剧场| 直男gayav资源| 国产精品国产高清国产av| 亚洲国产精品sss在线观看| 听说在线观看完整版免费高清| 精品一区二区三区视频在线| 亚洲欧美日韩高清在线视频| 精品一区二区三区视频在线观看免费| 男女做爰动态图高潮gif福利片| 色噜噜av男人的天堂激情| .国产精品久久| 亚洲欧美日韩无卡精品| 免费观看人在逋| 男人狂女人下面高潮的视频| 午夜视频国产福利| 97碰自拍视频| 久久精品综合一区二区三区| 午夜福利在线观看吧| 高清毛片免费看| 伦理电影大哥的女人| 搡老岳熟女国产| 免费搜索国产男女视频| 亚洲精品乱码久久久v下载方式| 精品不卡国产一区二区三区| 久久精品久久久久久噜噜老黄 | 国产精品国产高清国产av| 变态另类成人亚洲欧美熟女| 久久精品91蜜桃| 国产 一区精品| 久久久久久大精品| 午夜精品在线福利| 精品久久久久久久久久久久久| 成人综合一区亚洲| 少妇熟女aⅴ在线视频| 最新中文字幕久久久久| 一夜夜www| 国产日本99.免费观看| 欧美最黄视频在线播放免费| 九九在线视频观看精品| 国产91av在线免费观看| 国产精品伦人一区二区| 全区人妻精品视频| 亚洲精品日韩在线中文字幕 | 成人美女网站在线观看视频| 午夜视频国产福利| 国产aⅴ精品一区二区三区波| 22中文网久久字幕| 日韩强制内射视频| 在线看三级毛片| 国产成人一区二区在线| 国产高清视频在线播放一区| 国产av不卡久久| 亚洲人与动物交配视频| 韩国av在线不卡| 成人永久免费在线观看视频| 国产爱豆传媒在线观看| 国语自产精品视频在线第100页| 久久久久国产精品人妻aⅴ院| 在线观看66精品国产| 国产老妇女一区| 亚洲美女黄片视频| 国产一区二区三区在线臀色熟女| 国产日本99.免费观看| 国产精品99久久久久久久久| 久久精品国产亚洲av涩爱 | 免费观看人在逋| 国产v大片淫在线免费观看| 亚洲av中文字字幕乱码综合| 亚洲成人av在线免费| 国产蜜桃级精品一区二区三区| 美女cb高潮喷水在线观看| 99热这里只有是精品50| av.在线天堂| 日韩制服骚丝袜av| av女优亚洲男人天堂| 熟妇人妻久久中文字幕3abv| 少妇人妻一区二区三区视频| 日日干狠狠操夜夜爽| 久久99热6这里只有精品| 热99在线观看视频| 一进一出好大好爽视频| 最近中文字幕高清免费大全6| 狂野欧美白嫩少妇大欣赏| 中文字幕av在线有码专区| 大型黄色视频在线免费观看| 菩萨蛮人人尽说江南好唐韦庄 | 中文字幕av在线有码专区| 精品99又大又爽又粗少妇毛片| 亚洲国产色片| 人妻丰满熟妇av一区二区三区| 久久草成人影院| 久久久色成人| 国产一区二区三区av在线 | 国产精品不卡视频一区二区| 国产真实乱freesex| 两个人视频免费观看高清| 亚洲熟妇中文字幕五十中出| 免费看美女性在线毛片视频| 国产真实乱freesex| 免费观看精品视频网站| 赤兔流量卡办理| 国产精品1区2区在线观看.| 欧美bdsm另类| 一区二区三区免费毛片| 亚洲久久久久久中文字幕| 三级经典国产精品| www.色视频.com| 午夜亚洲福利在线播放| 别揉我奶头~嗯~啊~动态视频| 亚洲最大成人手机在线| 最近的中文字幕免费完整| 99久久成人亚洲精品观看| 内射极品少妇av片p| 日日啪夜夜撸| 日本-黄色视频高清免费观看| avwww免费| 国产精品人妻久久久影院| 成人无遮挡网站| 国产v大片淫在线免费观看| 天天一区二区日本电影三级| 久久人人精品亚洲av| 亚洲熟妇熟女久久| 亚洲国产精品成人久久小说 | 午夜福利成人在线免费观看| 伦精品一区二区三区| 日本与韩国留学比较| 亚洲性夜色夜夜综合| 成人性生交大片免费视频hd| av专区在线播放| 熟女电影av网| av国产免费在线观看| 日日摸夜夜添夜夜爱| 国产久久久一区二区三区| 国内精品美女久久久久久| 久久久久久久午夜电影| 99热6这里只有精品| 真人做人爱边吃奶动态| 精品一区二区三区视频在线观看免费| 日韩三级伦理在线观看| 久99久视频精品免费| 免费av观看视频| 一级a爱片免费观看的视频| 一夜夜www| 国产精品久久电影中文字幕| 小说图片视频综合网站| 精品乱码久久久久久99久播| 一本一本综合久久| 欧美潮喷喷水| 亚洲精品在线观看二区| 露出奶头的视频| 国产精品伦人一区二区| 内射极品少妇av片p| 日本 av在线| 日韩高清综合在线| 国产色爽女视频免费观看| 熟女电影av网| 有码 亚洲区| 久久精品国产亚洲网站| av在线亚洲专区| 神马国产精品三级电影在线观看| 久久人人爽人人爽人人片va| 性色avwww在线观看| 国产aⅴ精品一区二区三区波| 国产一级毛片七仙女欲春2| 男女边吃奶边做爰视频| 非洲黑人性xxxx精品又粗又长| 国产aⅴ精品一区二区三区波| 一区二区三区高清视频在线| 国产av不卡久久| 亚洲欧美精品综合久久99| 国内少妇人妻偷人精品xxx网站| 亚洲四区av| 变态另类成人亚洲欧美熟女| 嫩草影院入口| 婷婷精品国产亚洲av| 一夜夜www| 亚洲,欧美,日韩| 露出奶头的视频| 天堂av国产一区二区熟女人妻| 国产精品久久久久久精品电影| 色av中文字幕| 精品久久久久久久久亚洲| 韩国av在线不卡| 国内少妇人妻偷人精品xxx网站| 欧美色视频一区免费| 99视频精品全部免费 在线| 国产精品久久视频播放| 波野结衣二区三区在线| 久久午夜福利片| 欧美不卡视频在线免费观看| 69人妻影院| 日韩成人伦理影院| 中文字幕熟女人妻在线| 身体一侧抽搐| 久久久久久国产a免费观看| 99精品在免费线老司机午夜| 九色成人免费人妻av| 国产精品永久免费网站| 久久久久久久久大av| 欧美高清性xxxxhd video| 亚洲欧美成人综合另类久久久 | 九九热线精品视视频播放| 国产探花在线观看一区二区| 久久精品综合一区二区三区| 少妇的逼水好多| 综合色av麻豆| 中文字幕av在线有码专区| 亚洲精品国产成人久久av| 午夜视频国产福利| 欧美性猛交╳xxx乱大交人| 国产av麻豆久久久久久久| 九色成人免费人妻av| 国产精品电影一区二区三区| 国产黄a三级三级三级人| 成人综合一区亚洲| 亚州av有码| 一级av片app| 国产高清视频在线播放一区| 国产美女午夜福利| 亚洲乱码一区二区免费版| 免费在线观看成人毛片| 无遮挡黄片免费观看| 99久国产av精品| 97人妻精品一区二区三区麻豆| 欧美3d第一页| 午夜爱爱视频在线播放| 国产精品综合久久久久久久免费| 欧美色视频一区免费| 国产高潮美女av| 成人av在线播放网站| 一进一出抽搐gif免费好疼| 国产成人aa在线观看| 色哟哟哟哟哟哟| 能在线免费观看的黄片| 最好的美女福利视频网| 免费看日本二区| 日韩欧美精品v在线| 欧美丝袜亚洲另类| 欧美+亚洲+日韩+国产| 国产黄色视频一区二区在线观看 | 亚洲性久久影院| 免费大片18禁| 日韩精品中文字幕看吧| 一区二区三区免费毛片| 麻豆久久精品国产亚洲av| 乱系列少妇在线播放| 亚洲精品亚洲一区二区| 日本撒尿小便嘘嘘汇集6| 亚洲精品在线观看二区| aaaaa片日本免费| 少妇人妻一区二区三区视频| 午夜精品国产一区二区电影 | 国产黄a三级三级三级人| 小蜜桃在线观看免费完整版高清| 美女xxoo啪啪120秒动态图| 特级一级黄色大片| 久久韩国三级中文字幕| 亚洲电影在线观看av| 男人和女人高潮做爰伦理| 亚洲成人久久性| 可以在线观看的亚洲视频| 97碰自拍视频| 国产单亲对白刺激| 午夜a级毛片| 国产av一区在线观看免费| 深夜精品福利| 天堂av国产一区二区熟女人妻| 听说在线观看完整版免费高清| 日日摸夜夜添夜夜添av毛片| 成人午夜高清在线视频| 成人高潮视频无遮挡免费网站| 色播亚洲综合网| 日韩国内少妇激情av| 精品久久久久久久末码| 国产午夜福利久久久久久| 日本爱情动作片www.在线观看 | 亚洲精品国产成人久久av| 亚洲性久久影院| 色播亚洲综合网| 全区人妻精品视频| 午夜爱爱视频在线播放| 成人三级黄色视频| 色尼玛亚洲综合影院| 美女黄网站色视频| 看黄色毛片网站| 欧美一区二区国产精品久久精品| 久久久久久久久久成人| av卡一久久| av.在线天堂| 亚洲欧美中文字幕日韩二区| 精品人妻视频免费看| 成年av动漫网址| av免费在线看不卡| 午夜久久久久精精品| 99热只有精品国产| 人人妻人人看人人澡| 露出奶头的视频| 久久精品国产亚洲av涩爱 | 日本爱情动作片www.在线观看 | 日日干狠狠操夜夜爽| 少妇熟女aⅴ在线视频| 久久韩国三级中文字幕| 在线观看66精品国产| 欧美xxxx黑人xx丫x性爽| 亚洲国产精品合色在线| 卡戴珊不雅视频在线播放| 夜夜看夜夜爽夜夜摸| 久久九九热精品免费| 国产熟女欧美一区二区| 精品欧美国产一区二区三| 日本免费一区二区三区高清不卡| 亚洲成a人片在线一区二区| 人人妻人人澡欧美一区二区| 国产麻豆成人av免费视频| 精品一区二区三区视频在线观看免费| 毛片女人毛片| 国产高清视频在线播放一区| 男人和女人高潮做爰伦理| 欧美又色又爽又黄视频| 性欧美人与动物交配| www日本黄色视频网| 中文字幕人妻熟人妻熟丝袜美| 国产高清视频在线观看网站| 小蜜桃在线观看免费完整版高清| 久久精品影院6| 啦啦啦韩国在线观看视频| 色综合站精品国产| 欧美潮喷喷水| 久久久精品欧美日韩精品| 国产免费男女视频| 亚洲av第一区精品v没综合| 白带黄色成豆腐渣| 校园人妻丝袜中文字幕| 国产精品乱码一区二三区的特点| 男女之事视频高清在线观看| 免费无遮挡裸体视频| 男女边吃奶边做爰视频| 91狼人影院| 亚洲精品亚洲一区二区| 最新在线观看一区二区三区| 婷婷色综合大香蕉| 最近最新中文字幕大全电影3| 亚洲精品色激情综合| 一级毛片aaaaaa免费看小| 亚洲精品国产成人久久av| 免费看光身美女| 亚洲高清免费不卡视频| 国产免费一级a男人的天堂| 最近的中文字幕免费完整| 18禁在线播放成人免费| 在线观看午夜福利视频| 成人综合一区亚洲| 少妇人妻一区二区三区视频| 人妻少妇偷人精品九色| 淫妇啪啪啪对白视频| 久久精品国产亚洲av香蕉五月| 99热这里只有是精品在线观看| 欧美区成人在线视频| 亚洲av免费在线观看| 老女人水多毛片| 国产精品人妻久久久影院| 日日摸夜夜添夜夜添小说| 淫秽高清视频在线观看| 亚洲三级黄色毛片| 男女下面进入的视频免费午夜| av女优亚洲男人天堂| 悠悠久久av| 波多野结衣高清作品| 亚洲不卡免费看| 菩萨蛮人人尽说江南好唐韦庄 | 一卡2卡三卡四卡精品乱码亚洲| av天堂在线播放| 可以在线观看的亚洲视频| 国产女主播在线喷水免费视频网站 | 99热这里只有是精品在线观看| 久久久午夜欧美精品| 少妇人妻精品综合一区二区 | 亚洲成a人片在线一区二区| 男女那种视频在线观看| 国产熟女欧美一区二区| 狠狠狠狠99中文字幕| 一个人看视频在线观看www免费| 亚洲av熟女| 亚洲专区国产一区二区| 国产亚洲av嫩草精品影院| 两个人的视频大全免费| 国产中年淑女户外野战色| 婷婷精品国产亚洲av在线| 久久99热6这里只有精品| 欧美日韩在线观看h| 免费看日本二区| 男女下面进入的视频免费午夜| 91精品国产九色| 色噜噜av男人的天堂激情| 国产久久久一区二区三区| 菩萨蛮人人尽说江南好唐韦庄 | 国产免费一级a男人的天堂| 国产精品美女特级片免费视频播放器| 国产爱豆传媒在线观看| 国产精品久久视频播放| 中文亚洲av片在线观看爽| 亚洲欧美清纯卡通| 成人av一区二区三区在线看| 亚洲欧美中文字幕日韩二区| 久久久久久九九精品二区国产| 国产成人aa在线观看| 村上凉子中文字幕在线| 99热这里只有是精品在线观看| 黄片wwwwww| 亚洲久久久久久中文字幕| 精品一区二区免费观看| 亚洲欧美日韩高清在线视频| 久久久久久久午夜电影| 久久久久久久久久成人| 国产精品久久久久久av不卡| 一个人看的www免费观看视频| 国产精品久久视频播放| 亚洲第一区二区三区不卡| 可以在线观看的亚洲视频| 国产男靠女视频免费网站| 啦啦啦啦在线视频资源| 五月伊人婷婷丁香| 国产av一区在线观看免费| 日韩欧美三级三区| 国产欧美日韩精品一区二区| 嫩草影院精品99| 久久久久久久久久成人| 色播亚洲综合网| 亚洲电影在线观看av| 精品人妻熟女av久视频| 人人妻人人澡欧美一区二区| 精品99又大又爽又粗少妇毛片| 日日撸夜夜添| 成人美女网站在线观看视频| 欧美一区二区国产精品久久精品| 亚洲欧美日韩东京热| 禁无遮挡网站| 亚洲精品日韩在线中文字幕 | 久久久a久久爽久久v久久| a级毛片免费高清观看在线播放|