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

    基于特征模型的高超聲速飛行器再入多約束預(yù)測(cè)控制

    2017-03-28 03:07:21李亞輝朱晟楨
    彈道學(xué)報(bào) 2017年1期
    關(guān)鍵詞:超聲速時(shí)變飛行器

    張 軍,李亞輝,朱晟楨

    (1.江蘇大學(xué) 電氣信息工程學(xué)院,江蘇 鎮(zhèn)江 212013;2.中國(guó)運(yùn)載火箭技術(shù)研究院,北京 100191)

    基于特征模型的高超聲速飛行器再入多約束預(yù)測(cè)控制

    張 軍1,李亞輝2,朱晟楨1

    (1.江蘇大學(xué) 電氣信息工程學(xué)院,江蘇 鎮(zhèn)江 212013;2.中國(guó)運(yùn)載火箭技術(shù)研究院,北京 100191)

    針對(duì)存在輸入和輸出多約束的高超聲速飛行器再入控制問(wèn)題,提出了一種基于特征模型的魯棒預(yù)測(cè)控制方法。對(duì)耦合的非線性再入動(dòng)力學(xué)采用低階線性特征模型簡(jiǎn)化,設(shè)計(jì)魯棒預(yù)測(cè)控制器保證滿(mǎn)足再入控制的多約束?;谔卣髂P偷牡碗A等價(jià)性,將再入動(dòng)力學(xué)的輸入輸出動(dòng)態(tài)特性用線性時(shí)變的特征模型描述,并將三通道姿態(tài)耦合和干擾作為廣義復(fù)合干擾。將灰色理論中累加求和的思想引入到特征模型參數(shù)辨識(shí)算法中,弱化了再入隨機(jī)噪聲的影響,通過(guò)灰色辨識(shí)方法在線估計(jì)特征模型的時(shí)變參數(shù)和干擾。設(shè)計(jì)一種基于線性矩陣不等式(LMI)滾動(dòng)優(yōu)化技術(shù)的H2/H∞魯棒預(yù)測(cè)控制律與基于灰色辨識(shí)的補(bǔ)償器共同構(gòu)成復(fù)合控制器,保證了整個(gè)閉環(huán)系統(tǒng)的穩(wěn)定,滿(mǎn)足輸入和輸出約束。數(shù)據(jù)仿真驗(yàn)證了算法的有效性。

    高超聲速飛行器;約束;特征模型;魯棒預(yù)測(cè)控制;灰色辨識(shí)

    高超聲速飛行器具有響應(yīng)快速、機(jī)動(dòng)能力強(qiáng)等特點(diǎn),是未來(lái)國(guó)防裝備的重要發(fā)展方向,其中再入控制系統(tǒng)設(shè)計(jì)最為復(fù)雜。在高超聲速飛行器的再入過(guò)程中,運(yùn)行環(huán)境復(fù)雜,飛行器飛行高度和飛行馬赫數(shù)變化范圍大,飛行器的氣熱特性、氣動(dòng)特性發(fā)生劇烈變化,造成再入動(dòng)力學(xué)模型具有快時(shí)變、非線性耦合、強(qiáng)干擾等特性[1-2]。傳統(tǒng)的PID控制很難使再入系統(tǒng)具有良好的閉環(huán)控制性能,這給再入控制器的設(shè)計(jì)帶來(lái)了困難和挑戰(zhàn)[3-4]。更值得注意的是,高超聲速飛行器在再入過(guò)程中存在再入過(guò)程(動(dòng)壓、過(guò)載、熱流)、執(zhí)行結(jié)構(gòu)等約束,導(dǎo)致姿態(tài)控制系統(tǒng)存在控制和狀態(tài)等多約束,即舵機(jī)的偏角、氣流系角、角速度等,違反了這些約束,再入控制系統(tǒng)性能和穩(wěn)定性將受到嚴(yán)重的惡化[5-6]。

    目前,動(dòng)態(tài)逆控制、變結(jié)構(gòu)控制、神經(jīng)網(wǎng)絡(luò)控制等方法提高了再入系統(tǒng)魯棒性[7-9],但系統(tǒng)的輸入和輸出約束很少明確考慮。雖然采用高增益控制器的非線性控制和魯棒控制能保證系統(tǒng)具有良好的閉環(huán)性能,卻忽略了噪聲放大和需要更大控制能力的問(wèn)題[4],若考慮約束則原有的控制性能和穩(wěn)定性就會(huì)受到破壞,控制律需要重新設(shè)計(jì)。目前,抗飽和控制方法、輸入信號(hào)調(diào)節(jié)方法、預(yù)測(cè)控制等方法是解決約束控制的主要方法,文獻(xiàn)[10]將抗飽和控制和輸入?yún)⒖夹盘?hào)調(diào)節(jié)2種方法組合起來(lái)有效解決輸入約束,但是有一定的保守性,只有發(fā)生飽和行為時(shí),抗飽和補(bǔ)償器才執(zhí)行,影響再入的高精度控制和強(qiáng)機(jī)動(dòng)性能,并且不能處理輸出約束。因此,需要在再入控制系統(tǒng)的高性能要求和滿(mǎn)足約束之間進(jìn)行有效折中。

    預(yù)測(cè)控制具有反饋校正、模型預(yù)測(cè)、滾動(dòng)優(yōu)化等特點(diǎn),能有效處理系統(tǒng)化多約束,逐漸被航天控制所關(guān)注[11-12]。文獻(xiàn)[13]采用精確反饋線性化將再入非線性模型轉(zhuǎn)化為線性模型,通過(guò)滾動(dòng)優(yōu)化機(jī)制解決再入飛行器的輸入和輸出約束,但是控制約束轉(zhuǎn)化成非線性約束,計(jì)算仍然復(fù)雜,并且要求系統(tǒng)模型實(shí)時(shí)可逆。為了簡(jiǎn)化計(jì)算量,文獻(xiàn)[14-15]在高超聲速飛行器的上升段、再入段制導(dǎo)中設(shè)計(jì)了連續(xù)時(shí)變預(yù)測(cè)控制器,但沒(méi)有考慮魯棒性。預(yù)測(cè)模型是影響預(yù)測(cè)控制應(yīng)用的關(guān)鍵因素,若把原有的多耦合、強(qiáng)非線性再入動(dòng)力學(xué)模型直接作為預(yù)測(cè)模型,在線計(jì)算量很大,不能滿(mǎn)足高實(shí)時(shí)性要求。文獻(xiàn)[16]提出的特征模型用低階線性單入單出時(shí)變差分方程來(lái)描述高超聲速飛行器復(fù)雜動(dòng)力學(xué)特性,形式簡(jiǎn)單,工程運(yùn)行方便,對(duì)象特征模型與實(shí)際對(duì)象輸出是等價(jià)的,在工程上得到成功應(yīng)用。文獻(xiàn)[17]設(shè)計(jì)了基于特征模型的自適應(yīng)控制系統(tǒng),采用了黃金分割控制律,分析了閉環(huán)系統(tǒng)的穩(wěn)定性。文獻(xiàn)[18]主要考慮飛行器平穩(wěn)滑翔過(guò)程中滾動(dòng)角保持為零的條件下,基于特征模型設(shè)計(jì)簡(jiǎn)化的自適應(yīng)控制律。但是對(duì)于大角度再入機(jī)動(dòng)的情況,偏航和滾動(dòng)通道強(qiáng)耦合信息和非線性干擾不能簡(jiǎn)單忽略,再入特征模型需要考慮這一因素。

    基于特征模型的低階等價(jià)性,將再入動(dòng)力學(xué)的動(dòng)態(tài)特性用線性時(shí)變的特征模型描述,并把三通道姿態(tài)耦合和外來(lái)干擾作為廣義復(fù)合干擾,對(duì)復(fù)雜的再入動(dòng)力學(xué)進(jìn)行簡(jiǎn)化建模。利用灰色預(yù)測(cè)中“累加生成”的優(yōu)點(diǎn),重新生成辨識(shí)數(shù)據(jù),增強(qiáng)了辨識(shí)數(shù)據(jù)的規(guī)律性,通過(guò)灰色辨識(shí)在線估計(jì)特征模型的時(shí)變參數(shù)和干擾。在此基礎(chǔ)上,設(shè)計(jì)H2/H∞混合預(yù)測(cè)控制指標(biāo),基于LMI滾動(dòng)凸優(yōu)化獲得魯棒預(yù)測(cè)控制律,保證了整個(gè)閉環(huán)系統(tǒng)的穩(wěn)定,滿(mǎn)足系統(tǒng)狀態(tài)和輸入約束。

    1 再入飛行器動(dòng)力學(xué)模型和問(wèn)題描述

    忽略地球自轉(zhuǎn)、曲率、風(fēng)、氣動(dòng)舵面等因素對(duì)飛行器動(dòng)力學(xué)的影響,高超聲速飛行器再入動(dòng)力學(xué)模型為[4]

    (1)

    (2)

    ①舵偏角輸入飽和約束。

    |δa|≤δa,max,|δe|≤δe,max,|δr|≤δr,max。

    ②輸出約束。

    |α(k)|≤αmax,該約束保證了再入飛行器具有良好的氣動(dòng)特性和溫度分布,降低了防熱要求。

    |ωx|≤ωx,max,|ωy|≤ωy,max,|ωz|≤ωz,max,角速度約束保證了飛行器的結(jié)構(gòu)和力學(xué)安全要求。

    再入姿態(tài)控制系統(tǒng)的任務(wù)是:在滿(mǎn)足上述輸入和輸出約束的情況下,控制系統(tǒng)跟蹤制導(dǎo)系統(tǒng)給出的攻角、滾動(dòng)角指令,控制穩(wěn)態(tài)精度小于0.5°,同時(shí)側(cè)滑角保持為0。

    2 再入動(dòng)力學(xué)特征模型

    再入動(dòng)力學(xué)系統(tǒng)(1)和系統(tǒng)(2)是強(qiáng)非線性、時(shí)變、不確定的多變量耦合系統(tǒng),直接設(shè)計(jì)控制器比較困難。基于文獻(xiàn)[16]系統(tǒng)的特征模型輸入輸出等價(jià)原則,對(duì)系統(tǒng)(1)、系統(tǒng)(2)的控制輸出α,β,φ求二階導(dǎo)數(shù),進(jìn)行等價(jià)代換和泰勒級(jí)數(shù)近似展開(kāi),得到輸入輸出時(shí)變的二階攻角離散模型:

    α(k+1)=fα1(k)α(k)+fα2(k)α(k-1)+
    g1(k)δα(k)+Δα(k)
    fα1+fα2≈1,δα≈δe

    (3)

    β(k+1)=fβ1(k)β(k)+fβ2(k)β(k-1)+
    g2(k)δβ(k)+Δβ(k)
    fβ1+fβ2≈1,δβ≈δr

    (4)

    φ(k+1)=fφ1(k)φ(k)+fφ2(k)α(k-1)+
    g3(k)δφ(k)+Δφ
    fφ1+fφ2≈1,δφ≈δa

    (5)

    式中:fα1,fα2,g1為攻角通道的特征模型時(shí)變參數(shù);δα,δβ,δφ為三通道的等效控制輸入;fβ1,fβ2,g2為側(cè)滑角通道的特征模型時(shí)變參數(shù);fφ1,fφ2,g3為滾動(dòng)角通道的特征模型時(shí)變參數(shù)。模型(3)~模型(5)是單入單出的時(shí)變系統(tǒng),明顯降低了原始再入動(dòng)力學(xué)的復(fù)雜性。Δα(k)是攻角通道的廣義擾動(dòng),包含其他特征通道的動(dòng)態(tài)耦合信息、外部干擾,反映了再入飛行器的三通道強(qiáng)耦合、非線性、時(shí)變性,在飛行器大角度機(jī)動(dòng)時(shí)該量不能忽略,這是與文獻(xiàn)[16-17]的不同之處。定義:

    Δα(k)=m11β(k)+m12φ(k)+Δα1(k)

    (6)

    3 再入特征模型的灰色辨識(shí)

    高超聲速飛行器的再入特征模型是時(shí)變的單入單出系統(tǒng),需要線辨識(shí)模型參數(shù)。而再入狀態(tài)中含有近似均值為0的白噪聲,影響辨識(shí)精度。目前,灰色預(yù)測(cè)的數(shù)據(jù)累加生成可有效弱化原始序列中隨機(jī)噪聲因素的影響,增強(qiáng)了數(shù)據(jù)的規(guī)律性,這里采用灰色辨識(shí)估計(jì)特征模型的參數(shù)fα1(k),fα2(k),g1(k),m11(k),m12(k),Δα1(k),具體步驟如下[18]。

    (7)

    進(jìn)行一次累加生成(式中:m為采樣數(shù)據(jù)數(shù)目),得到規(guī)律更強(qiáng)的新序列,消除隨機(jī)噪聲的影響,作為辨識(shí)的輸入數(shù)據(jù):

    {X1(1),X1(2),…,X1(N)}

    (8)式中:N為灰色序列數(shù)據(jù)數(shù)目,N≥n+1。由上述灰色數(shù)據(jù)定義得到再入動(dòng)力學(xué)特征模型(3)的灰色模型為

    (9)

    在辨識(shí)時(shí)間內(nèi),建模誤差和干擾近似認(rèn)為

    令V=(fα1(k)fα2(k)g1(k)m11m12Δα1)T,

    (10)

    由灰色模型可得:

    (11)

    為提高再入特征模型控制的魯棒性,這里設(shè)計(jì)前饋補(bǔ)償控制:

    (12)

    4 H2/H∞魯棒預(yù)測(cè)控制

    令Δfα1(k),Δfα2(k),Δg1(k)為特征模型系數(shù)辨識(shí)誤差;Δ為uc(k)補(bǔ)償Δα(k)后的誤差,它們是有界的。此外,由于α(k),δα(k)也是明確的區(qū)域約束,因此,參數(shù)辨識(shí)誤差、未建模動(dòng)態(tài)、外界干擾等因素引起的綜合擾動(dòng)也是有界的,這里作為廣義擾動(dòng),即

    d(k)=Δfα1(k)α(k)+Δfα2(k)α(k-1)+

    Δg1(k)δα(k)+Δ

    灰色辨識(shí)補(bǔ)償后的特征模型為

    α(k+1)=fα1α(k)+fα2α(k-1)+g1δα1(k)+d(k)
    δα=uc+δα1

    (13)

    x(k+1)=A(k)x(k)+B(k)δα1(k)+Bdd(k)

    (14)

    對(duì)原有系統(tǒng)(1)的約束進(jìn)行轉(zhuǎn)化。

    ①攻角約束轉(zhuǎn)化為系統(tǒng)(14)的輸出約束:

    |α(k)|=|(0 1 0)x(k)|≤αmax

    ②攻角變化率和飛行器最大角速度約束轉(zhuǎn)化為系統(tǒng)(14)的輸出約束:

    (15)

    式中:T為控制周期。

    ③考慮灰色估計(jì)補(bǔ)償器后,新的輸入約束:

    u=|uc(k)+δα1(k)|≤δe,max

    (16)

    |δα1(k)|≤min(|δe,max-uc|,|-δe,max-uc|)=umax

    (17)

    ①H∞性能指標(biāo)。對(duì)于任意給定的γ>0,從d(k)到z(k)的傳遞函數(shù)滿(mǎn)足:

    ‖Gzd‖∞≤γ

    (18)

    ②H2性能指標(biāo)。

    (19)

    這里采用LMI設(shè)計(jì)控制器處理輸入和輸出約束,通過(guò)Lyapunov理論保證系統(tǒng)(14)的魯棒性能和穩(wěn)定性?;贖2/H∞性能指標(biāo)提出如下定理。

    定理1 設(shè)x(k)是系統(tǒng)(14)在采樣k時(shí)刻的狀態(tài)測(cè)量值,系統(tǒng)采用狀態(tài)反饋控制律δα1(k+i)=K(k)x(k+i),F(k)=YQ-1,實(shí)現(xiàn)閉環(huán)系統(tǒng)魯棒漸近穩(wěn)定的充分條件是存在變量:b>0,Q>0,滿(mǎn)足如下LMI多約束的優(yōu)化問(wèn)題:

    C1=(0 1 0)

    證明 選取V(x(k))=xT(k)P(k)x(k),V(x(k+1))=xT(k+1)P(k)x(k+1),定義x(k+i)表示基于x(k)測(cè)量在k+i時(shí)刻的預(yù)測(cè)狀態(tài),δα1(k+i)是k時(shí)刻第i步預(yù)測(cè)控制,則控制系統(tǒng)的二次穩(wěn)定需要滿(mǎn)足:

    (20)

    由灰色辨識(shí)的收斂穩(wěn)定性知道,補(bǔ)償d(k)趨近于0,存在:

    γ2dT(k+i)d(k+i)-zT(k+i)z(k+i)<0。

    當(dāng)M(k)<0時(shí),V(x(k+i+1))-V(x(k+i))<0,當(dāng)i=∞時(shí),V(x(k+i))趨近于0,將式(20)左右兩邊從i=0到i=∞相加,得到:

    (21)

    其中,M(k)<0是實(shí)現(xiàn)設(shè)計(jì)目標(biāo)的充要條件。

    (22)

    利用Schur補(bǔ)性質(zhì),式(22)可以寫(xiě)成下面的矩陣不等式:

    (23)

    令Q=bP-1,K=YQ-1,對(duì)式(23)分別左乘、右乘diag{b1/2P-1,b1/2I,b1/2I,b1/2I,b1/2I},得到如下不等式:

    (24)

    H2性能指標(biāo)式(19)等價(jià)如下LMI:

    (25)

    1)輸出約束分析。

    由M(k)<0,因此V(x(k+i+1))≤V(x(k+i)),推導(dǎo)zT(k)Q-1z(k)≤1 是不變集。

    (26)

    因?yàn)椤瑈(k+i|k)‖2≤ymax,所以

    (27)

    等價(jià)如下不等式:

    (28)

    因此,攻角約束轉(zhuǎn)化為

    (29)

    角速度約束轉(zhuǎn)化為

    (30)

    2)輸入約束分析。

    灰色前饋補(bǔ)償后的舵偏角約束為

    (31)

    令K=YQ-1,則δα1=Kx=YQ-1x,不等式(31)可以寫(xiě)成如下的形式:

    (32)

    利用Schur補(bǔ)性質(zhì),式(32)等效如下的不等式:

    (33)

    3)算法實(shí)現(xiàn)。

    本文提出的算法如下:

    ①基于實(shí)時(shí)的輸入和輸出數(shù)據(jù),利用灰色辨識(shí)攻角通道的特征模型參數(shù)fα1,fα2,g1,m11,m12,Δα1(k)等;

    ②利用灰色辨識(shí)結(jié)果,設(shè)計(jì)補(bǔ)償器uc補(bǔ)償三通耦合運(yùn)動(dòng)信息、未建模信息和干擾Δα1(k);

    ③計(jì)算魯棒預(yù)測(cè)控制律δα1=Kx,將總的控制u=δα1+uc施加到再入飛行器上;

    ④重復(fù)前面過(guò)程 ①。

    其他通道控制算法與此類(lèi)似。

    5 仿真研究

    針對(duì)某面對(duì)稱(chēng)的高超聲速飛行器確定的飛行狀態(tài)(αβφ)=(1000)(單位:(°)),獲得相應(yīng)的特征模型參數(shù)范圍,設(shè)飛行器的高度為40 km,Ma=10,則其氣動(dòng)系數(shù)的范圍如下[4]:

    kxx∈[-0.358 5,-1.434 1],
    kxz∈[-0.010 1,-0.015 2],
    kyy∈[-0.106 0,-0.294 5],
    kzx∈[-0.012 7,-0.050 6],
    kzz∈[-0.025 3,-0.038 0],
    kxa∈[40.590 9,34.339 5],
    kxr∈[3.181 9,3.172 9],
    kye∈[25.322 5,22.790 2],
    kza∈[2.719 6,2.284 5],
    kzr∈[4.536 3,4.133 8]。

    假設(shè)飛行器姿態(tài)角速度的最大值為|ωx|max=|ωy|max=|ωz|max=1.047 rad/s,輸入控制量最大值為30°,采樣時(shí)間為0.01 s。

    攻角控制通道:

    α(k+1)=fα1α(k)+fα2α(k-1)+g1δα(k)+Δα,
    fα1∈[1.996 8,2.003 2],
    fα2∈[-1.003 2,-0.996 8],
    g1∈[2.275 4×10-3,2.536 3×10-3]。

    側(cè)滑角控制通道:

    β(k+1)=fβ1β(k)+fβ2β(k-1)+g2δβ(k)+Δβ,
    fβ1∈[1.996 5,2.003 5],
    fβ2∈[-1.003 5,-0.996 5],
    g2∈[3.376 0×10-4,4.057 7×10-4]。

    滾動(dòng)角控制通道:

    φ(k+1)=fφ1φ(k)+fφ2α(k-1)+g3δφ(k)+Δφ
    fφ1∈[1.982 3,2.018 1],
    fφ2∈[-1.018 1,-0.982 3],
    g3∈[3.305 5×10-3,4.195 7×10-3]。

    控制系統(tǒng)的要求為側(cè)滑角保持在0°附近,攻角跟蹤精度小于0.5°。

    1)有時(shí)變擾動(dòng)特征模型的攻角跟蹤。

    圖1 時(shí)變干擾下的攻角跟蹤曲線

    圖2 時(shí)變干擾下的控制舵角曲線

    圖3 時(shí)變干擾下的角速度曲線

    2)隨機(jī)噪聲和時(shí)變擾動(dòng)下的攻角跟蹤。

    圖4 隨機(jī)噪聲和時(shí)變擾動(dòng)下的攻角跟蹤曲線

    圖5 隨機(jī)噪聲和時(shí)變擾動(dòng)下的控制舵角曲線

    圖6 隨機(jī)噪聲和時(shí)變擾動(dòng)下的角速度曲線

    6 結(jié)論

    針對(duì)有輸入和輸出約束的非線性再入飛行器動(dòng)力學(xué)模型,基于特征模型設(shè)計(jì)穩(wěn)定的魯棒預(yù)測(cè)控制律。采用時(shí)變的特征模型描述再入動(dòng)力學(xué)模型,把三通道姿態(tài)耦合和外來(lái)干擾作為廣義復(fù)合干擾,對(duì)復(fù)雜的再入動(dòng)力學(xué)進(jìn)行簡(jiǎn)化建模,并通過(guò)灰色辨識(shí)估計(jì)特征模型的時(shí)變參數(shù);基于LMI滾動(dòng)凸優(yōu)化設(shè)計(jì)H2/H∞魯棒預(yù)測(cè)控制律,保證了整個(gè)閉環(huán)系統(tǒng)的魯棒穩(wěn)定,滿(mǎn)足系統(tǒng)的輸入和輸出多約束,通過(guò)數(shù)值仿真驗(yàn)證了本文算法的有效性。

    [1] HARPOLD J.Shuttle entry guidance[J].Journal of the Astronautical Sciences,1979,28(3):239-268.

    [2] LU P.Entry guidance:a unified method[J].Journal of Guidance,Control,and Dynamics,2014,37(3):713-728.

    [3] WU H,MENG B.Review of the control research of hypersonic vehicle[J].Advances In Mechanics,2009,39(6):756-765.

    [4] ZHANG Z,HU J.Stability analysis of a hypersonic vehicle controlled by the characteristic model based adaptive controller[J].Science China:Information Sciences,2012,55(10):2 243-2 256.

    [5] SHEN Z J,LU P.On-board generation of three-dimensional constrained entry trajectories[J].Journal of Guidance,Control,and Dynamics,2003,26(1):111-121.

    [6] LU P.Regulation about time-varying trajectories:precision entry guidance illustrated[J].Journal of Guidance,Control,and Dynamics,1999,22(6):784-790.

    [7] XU B,SUN F,LIU H,REN J.Adaptive kriging controller design for hypersonic flight vehicle via back-stepping[J].IET Control Theory and Applications,2012,4(6):487-497.

    [8] CHARLES E H,YURI B S.Sliding mode disturbance observer-based control for a reusable launch vehicle[J].Journal of Guidance,Control,and Dynamics,2006,29(6):1 315-1 328.

    [9] YURI B S,CHARLES E H.Reusable launch vehicle control in multiple-time scale sliding modes[J].Journal of Guidance,Control,and Dynamics,2000,23(6):1 013-1 020.

    [10] ALICIA Z,ANDREA S.Combined reference governor and anti-windup design for constrained hypersonic vehicles models[C]//AIAA Guidance,Navigation,and Control Conference.New York:AIAA,2009:6 283-6 289.

    [11] CHIRISTOPHER P,MORGAN B,ILYA K.Model predictive control guidance with extended command governor inner-loop flight control for hypersonic vehicles[C]//AIAA Guidance,Navigation,and Control Conference and Exhibit.New York:AIAA,2013:5 028-5 048.

    [12] ZHANG J,SUN T R.Disturbance observer-based sliding manifold predictive control for reentry hypersonic vehicles with multi-constraint[J].Journal of Aerospace Engineering,2016,230(3):485-495.

    [13] VAN W R.Combined feedback linearization and constrained model predictive control for entry flight[J].Journal of Guidance,Control,and Dynamics,2006,29(2):427-434.

    [14] BRINDA V,DASGUPTA S.Guidance law for an air-breathing launch vehicle using predictive control concept[J].Journal of Guiance,Control,and Dynamics,2006,29(6):1 460-1464.

    [15] WU H X,HU J,XIE Y C.Characteristic model-based all-coefficient adaptive control method and its applications[J].IEEE Transactions on Systems,Man,and Cybernetics-Prat C:Applications and Reviews,2007,37(2):213-221.

    [16] XU L J.The hypersonic vehicle control based on the all-channel-coupled characteristic model[J].Aerospace Control and Application,2012,38(4):12-18.

    [17] GONG Y L,WU H X.Characteristic model-based adaptive attitude control for hypersonic vehicle[J].Journal of Astronautics,2010,31(9):2 122-2 128.

    [18] LIU L,CHEN R W.Modeling algorithm using accumulated generating operation and weighted product[J].System Engineering and Electronics,2010,32(5):976-979.

    [19] PATIENCE E,ORUKPE I M.Model predictive control based on mixedH2/H∞control approach[C]//Proceedings of the 2007 American Control Conference.New York:IEEE,2007:11-13.

    Predictive Control Based on Characteristic Model for Reentry Hypersonic Vehicle With Multiple Constraints

    ZHANG Jun1,LI Ya-hui2,ZHU Sheng-zhen1

    (1.Electrical and Information Engineering College,Jiangsu University,Zhenjiang 212013,China; 2.China Academy of Launch Vehicle Technology,Beijing 100191,China)

    In order to study the reentry control problem of hypersonic vehicle with multiple constraints,a novel robust predictive-control-method was proposed based on the characteristic model.The low-order linear characteristic model was adopted to simplify the coupled nonlinear reentry dynamics,and the robust-predictive-controller was designed to satisfy the multiple constraints in the reentry control.Based on the low-order equivalence of characteristic model,the input and output dynamic-characteristic of the reentry dynamics was modeled by linear time-varying model,and the three-channel coupling and external disturbance were considered as a generalized composite disturbance.The accumulative sum idea in the grey theory was introduced into the parameter identification of characteristic model,and the effect of random noise factor was reduced,and the regularity of identification data was strengthened.The time-varying parameters and the disturbance were on-line estimated based on grey identification.The mixedH2/H∞r(nóng)obust predictive control-law was proposed by linear matrix inequality(LMI)receding horizon optimization technique,and the compensator of grey identification was feedforward compensated.The closed-loop system is stable,which satisfies the input and state constraints.The validity of the algorithm was verified by simulation results.

    hypersonic vehicle;constraints;characteristic model;robust predictive control;gray identification

    2016-08-19

    國(guó)家自然科學(xué)基金項(xiàng)目(51379044;51405303;61503158)

    張軍(1972- ),男,副教授,博士,研究方向?yàn)轱w行器控制,預(yù)測(cè)控制。E-mail:jzhang@ujs.edu.cn。

    V249

    A

    1004-499X(2017)01-0001-08

    猜你喜歡
    超聲速時(shí)變飛行器
    高超聲速出版工程
    高超聲速飛行器
    超聲速旅行
    復(fù)雜飛行器的容錯(cuò)控制
    電子制作(2018年2期)2018-04-18 07:13:25
    基于時(shí)變Copula的股票市場(chǎng)相關(guān)性分析
    煙氣輪機(jī)復(fù)合故障時(shí)變退化特征提取
    基于MEP法的在役橋梁時(shí)變可靠度研究
    神秘的飛行器
    高超聲速大博弈
    太空探索(2014年5期)2014-07-12 09:53:28
    自適應(yīng)兩級(jí)UKF算法及其在時(shí)變偏差估計(jì)中的應(yīng)用
    免费观看人在逋| 丰满人妻熟妇乱又伦精品不卡| 99久久综合免费| 韩国高清视频一区二区三区| 国产精品二区激情视频| 久久天躁狠狠躁夜夜2o2o | 亚洲精品中文字幕在线视频| 亚洲第一青青草原| 天堂中文最新版在线下载| 极品人妻少妇av视频| 久久热在线av| 精品少妇一区二区三区视频日本电影| 男人添女人高潮全过程视频| 欧美在线一区亚洲| 高清av免费在线| 午夜福利影视在线免费观看| 丝袜喷水一区| 欧美激情高清一区二区三区| 国产一区亚洲一区在线观看| 久久精品亚洲熟妇少妇任你| 日韩精品免费视频一区二区三区| 老鸭窝网址在线观看| 日韩中文字幕视频在线看片| 下体分泌物呈黄色| 午夜免费男女啪啪视频观看| 免费人妻精品一区二区三区视频| 亚洲av电影在线观看一区二区三区| 9热在线视频观看99| 男女午夜视频在线观看| 午夜福利乱码中文字幕| 久久午夜综合久久蜜桃| 久久人人爽av亚洲精品天堂| 啦啦啦中文免费视频观看日本| 日本欧美国产在线视频| 日韩 亚洲 欧美在线| 中文字幕亚洲精品专区| 欧美变态另类bdsm刘玥| 国产av精品麻豆| 国产欧美日韩综合在线一区二区| 国产精品免费大片| 又粗又硬又长又爽又黄的视频| 久久精品久久精品一区二区三区| 亚洲欧美一区二区三区久久| 午夜免费鲁丝| 日日爽夜夜爽网站| 在线观看免费高清a一片| 咕卡用的链子| 国产精品国产三级国产专区5o| 国产一级毛片在线| 久久人人97超碰香蕉20202| 好男人电影高清在线观看| 搡老岳熟女国产| 久久av网站| 90打野战视频偷拍视频| 国产精品av久久久久免费| 最近最新中文字幕大全免费视频 | 国产男人的电影天堂91| 日韩大片免费观看网站| 亚洲国产中文字幕在线视频| 国产99久久九九免费精品| 亚洲中文日韩欧美视频| 欧美人与性动交α欧美软件| 丝袜脚勾引网站| 亚洲欧美日韩另类电影网站| 美女视频免费永久观看网站| 日韩视频在线欧美| 亚洲国产精品999| 永久免费av网站大全| 久久九九热精品免费| 国产成人免费观看mmmm| 亚洲国产欧美网| 乱人伦中国视频| 看十八女毛片水多多多| 久久久国产精品麻豆| 汤姆久久久久久久影院中文字幕| 黄色a级毛片大全视频| 国产精品久久久久久精品古装| 欧美黑人欧美精品刺激| 色婷婷av一区二区三区视频| 热99国产精品久久久久久7| 国产真人三级小视频在线观看| 国产免费福利视频在线观看| 国产高清videossex| 精品少妇一区二区三区视频日本电影| 日韩,欧美,国产一区二区三区| 日本色播在线视频| 精品国产国语对白av| 精品亚洲乱码少妇综合久久| 99国产精品一区二区蜜桃av | 国产精品国产av在线观看| 日韩av不卡免费在线播放| 男女高潮啪啪啪动态图| 亚洲国产欧美日韩在线播放| 你懂的网址亚洲精品在线观看| 一本—道久久a久久精品蜜桃钙片| 老司机深夜福利视频在线观看 | 永久免费av网站大全| 欧美少妇被猛烈插入视频| 一边摸一边做爽爽视频免费| 成人亚洲精品一区在线观看| 国产亚洲欧美精品永久| 建设人人有责人人尽责人人享有的| 国产欧美日韩精品亚洲av| 一区在线观看完整版| 黄色毛片三级朝国网站| 日本五十路高清| 成年美女黄网站色视频大全免费| 2018国产大陆天天弄谢| 性色av一级| 久久精品国产亚洲av涩爱| 在线观看免费午夜福利视频| av不卡在线播放| 久久亚洲精品不卡| 日本色播在线视频| 免费在线观看视频国产中文字幕亚洲 | 久久久久久免费高清国产稀缺| 国产伦理片在线播放av一区| 性少妇av在线| 国产精品99久久99久久久不卡| 赤兔流量卡办理| 亚洲五月婷婷丁香| 99久久99久久久精品蜜桃| 涩涩av久久男人的天堂| 最近最新中文字幕大全免费视频 | 亚洲中文av在线| 国产麻豆69| 热99国产精品久久久久久7| 国产精品香港三级国产av潘金莲 | 亚洲精品国产av蜜桃| 精品国产一区二区三区久久久樱花| 国产一区亚洲一区在线观看| 久久久精品区二区三区| 色网站视频免费| av网站在线播放免费| 又紧又爽又黄一区二区| 色94色欧美一区二区| 在线av久久热| 少妇 在线观看| 国产精品二区激情视频| 中文字幕色久视频| 建设人人有责人人尽责人人享有的| 啦啦啦 在线观看视频| 一区二区av电影网| 一级a爱视频在线免费观看| 99re6热这里在线精品视频| 国产国语露脸激情在线看| 国产在视频线精品| 欧美黄色片欧美黄色片| 久久久久久亚洲精品国产蜜桃av| 国产视频一区二区在线看| 高清视频免费观看一区二区| 国产精品人妻久久久影院| 日韩 亚洲 欧美在线| 精品少妇黑人巨大在线播放| 国产精品人妻久久久影院| 麻豆乱淫一区二区| 国产av精品麻豆| 国产男女内射视频| 国产欧美日韩一区二区三区在线| 纵有疾风起免费观看全集完整版| 日日夜夜操网爽| 丁香六月天网| 国产色视频综合| 丁香六月欧美| 国产精品国产av在线观看| 操出白浆在线播放| 国产精品99久久99久久久不卡| 久久精品国产亚洲av涩爱| 国产成人免费观看mmmm| 超碰97精品在线观看| 天天添夜夜摸| av天堂久久9| 999久久久国产精品视频| 免费观看人在逋| 成人影院久久| tube8黄色片| 免费人妻精品一区二区三区视频| 在现免费观看毛片| av有码第一页| 激情五月婷婷亚洲| 一个人免费看片子| 欧美日韩一级在线毛片| 国产精品免费视频内射| 国产高清不卡午夜福利| 黑人欧美特级aaaaaa片| 亚洲国产精品国产精品| 黄色一级大片看看| 最黄视频免费看| 久久国产亚洲av麻豆专区| 亚洲第一av免费看| bbb黄色大片| 亚洲精品日本国产第一区| 亚洲,欧美精品.| 少妇人妻久久综合中文| 国产精品一二三区在线看| 色网站视频免费| 男女国产视频网站| 一个人免费看片子| 日本一区二区免费在线视频| 精品少妇一区二区三区视频日本电影| 黄频高清免费视频| 国精品久久久久久国模美| 91成人精品电影| 国产精品 国内视频| 真人做人爱边吃奶动态| 亚洲成人国产一区在线观看 | 国产精品国产av在线观看| 黄色片一级片一级黄色片| 一边摸一边抽搐一进一出视频| 亚洲久久久国产精品| 可以免费在线观看a视频的电影网站| 久久国产精品影院| 性色av乱码一区二区三区2| 精品欧美一区二区三区在线| 一本久久精品| www.自偷自拍.com| 亚洲人成77777在线视频| e午夜精品久久久久久久| 少妇 在线观看| 天天添夜夜摸| 免费在线观看日本一区| 一区二区日韩欧美中文字幕| 亚洲欧洲国产日韩| 婷婷成人精品国产| 国产黄色免费在线视频| 91精品伊人久久大香线蕉| 考比视频在线观看| av片东京热男人的天堂| 亚洲国产欧美日韩在线播放| 亚洲精品久久久久久婷婷小说| 欧美日韩国产mv在线观看视频| 亚洲精品一二三| 十分钟在线观看高清视频www| 天堂俺去俺来也www色官网| 欧美激情高清一区二区三区| 热re99久久国产66热| 久久久精品免费免费高清| 一区二区av电影网| 国产高清视频在线播放一区 | 日韩av免费高清视频| 国产高清视频在线播放一区 | 亚洲欧洲精品一区二区精品久久久| 丝瓜视频免费看黄片| 日韩人妻精品一区2区三区| 精品高清国产在线一区| 欧美人与善性xxx| 久久久久国产一级毛片高清牌| 欧美成人精品欧美一级黄| 天天操日日干夜夜撸| 热re99久久国产66热| 欧美激情 高清一区二区三区| 欧美黄色淫秽网站| 亚洲精品国产av蜜桃| 视频在线观看一区二区三区| 老司机亚洲免费影院| 国产欧美亚洲国产| 两个人看的免费小视频| 国产又色又爽无遮挡免| 久久久久国产一级毛片高清牌| 国产欧美日韩一区二区三区在线| 香蕉国产在线看| 亚洲欧美中文字幕日韩二区| 亚洲欧美一区二区三区黑人| 一本综合久久免费| 大香蕉久久网| 亚洲av电影在线观看一区二区三区| 亚洲第一av免费看| 在线av久久热| 精品人妻1区二区| 久久这里只有精品19| 欧美成人午夜精品| 老鸭窝网址在线观看| 女人爽到高潮嗷嗷叫在线视频| 深夜精品福利| 国产视频首页在线观看| 两个人免费观看高清视频| 最新在线观看一区二区三区 | 亚洲av国产av综合av卡| 亚洲国产精品成人久久小说| 国产免费福利视频在线观看| 亚洲欧美清纯卡通| 90打野战视频偷拍视频| 国产免费视频播放在线视频| 九草在线视频观看| 亚洲 欧美一区二区三区| 国产午夜精品一二区理论片| 国产精品亚洲av一区麻豆| 国产高清视频在线播放一区 | 9色porny在线观看| avwww免费| 成人国产av品久久久| 午夜福利视频精品| 美女午夜性视频免费| 天天躁夜夜躁狠狠久久av| 黄色a级毛片大全视频| 亚洲国产看品久久| 精品久久蜜臀av无| 日本91视频免费播放| 99香蕉大伊视频| 亚洲欧美日韩另类电影网站| 中国美女看黄片| 日本欧美国产在线视频| 成在线人永久免费视频| 人人妻人人添人人爽欧美一区卜| 熟女av电影| 国产精品欧美亚洲77777| 午夜福利视频在线观看免费| av片东京热男人的天堂| av在线app专区| 女警被强在线播放| 久久久久久人人人人人| 69精品国产乱码久久久| 欧美亚洲日本最大视频资源| 日本色播在线视频| 每晚都被弄得嗷嗷叫到高潮| 亚洲色图 男人天堂 中文字幕| 国产黄色视频一区二区在线观看| 满18在线观看网站| 不卡av一区二区三区| 成人亚洲欧美一区二区av| 高清不卡的av网站| 亚洲av日韩精品久久久久久密 | 久久久欧美国产精品| 国产亚洲欧美精品永久| 一区福利在线观看| 亚洲国产欧美在线一区| 亚洲第一av免费看| 午夜激情久久久久久久| 亚洲欧美色中文字幕在线| 国产成人免费观看mmmm| 波多野结衣av一区二区av| 少妇人妻久久综合中文| 精品亚洲成a人片在线观看| 亚洲国产av新网站| 精品亚洲乱码少妇综合久久| 精品国产一区二区久久| 嫩草影视91久久| 别揉我奶头~嗯~啊~动态视频 | 亚洲久久久国产精品| 成人18禁高潮啪啪吃奶动态图| 国产精品久久久人人做人人爽| 国产男女内射视频| av视频免费观看在线观看| 一边亲一边摸免费视频| 国产日韩欧美在线精品| 又粗又硬又长又爽又黄的视频| 亚洲欧美色中文字幕在线| 18禁裸乳无遮挡动漫免费视频| 新久久久久国产一级毛片| 18禁黄网站禁片午夜丰满| 国产无遮挡羞羞视频在线观看| 免费高清在线观看视频在线观看| 精品国产超薄肉色丝袜足j| 亚洲精品久久久久久婷婷小说| 狂野欧美激情性bbbbbb| 国产精品国产三级专区第一集| 一级a爱视频在线免费观看| 美女大奶头黄色视频| 下体分泌物呈黄色| 男女边摸边吃奶| 国产精品99久久99久久久不卡| 亚洲黑人精品在线| 亚洲欧美精品综合一区二区三区| 男女边摸边吃奶| 少妇 在线观看| 亚洲精品自拍成人| 亚洲欧美激情在线| 成年人黄色毛片网站| 久久久精品区二区三区| 久9热在线精品视频| 每晚都被弄得嗷嗷叫到高潮| 国产成人av教育| 亚洲 国产 在线| 亚洲欧美成人综合另类久久久| 成年av动漫网址| 国产成人av教育| 国产深夜福利视频在线观看| 日本vs欧美在线观看视频| 新久久久久国产一级毛片| 最新在线观看一区二区三区 | 一区二区av电影网| 久久久国产欧美日韩av| 中文精品一卡2卡3卡4更新| 国产免费又黄又爽又色| 国产一区亚洲一区在线观看| 国产在线观看jvid| 又大又爽又粗| 性色av一级| 制服诱惑二区| 中文字幕制服av| 99热国产这里只有精品6| 999久久久国产精品视频| 亚洲欧洲国产日韩| tube8黄色片| 国产精品人妻久久久影院| 国产日韩欧美在线精品| 国产精品欧美亚洲77777| 日韩视频在线欧美| 一级毛片电影观看| 熟女av电影| 午夜免费鲁丝| 亚洲成色77777| 丰满少妇做爰视频| 丝袜美腿诱惑在线| 亚洲av电影在线观看一区二区三区| 亚洲欧美色中文字幕在线| 国产熟女午夜一区二区三区| 国产在视频线精品| 精品少妇黑人巨大在线播放| 又粗又硬又长又爽又黄的视频| 欧美日韩黄片免| 丝袜美腿诱惑在线| 久久国产精品人妻蜜桃| 午夜影院在线不卡| 波多野结衣一区麻豆| 亚洲欧美清纯卡通| 免费在线观看完整版高清| 99久久精品国产亚洲精品| 又紧又爽又黄一区二区| 日韩欧美一区视频在线观看| 久久精品亚洲av国产电影网| 国产精品三级大全| 丝袜在线中文字幕| 亚洲七黄色美女视频| 国产99久久九九免费精品| av线在线观看网站| 超色免费av| 精品久久久精品久久久| 久久精品国产亚洲av涩爱| xxx大片免费视频| 日韩制服丝袜自拍偷拍| 欧美亚洲 丝袜 人妻 在线| av不卡在线播放| 男的添女的下面高潮视频| 国产亚洲av高清不卡| 成人黄色视频免费在线看| 好男人电影高清在线观看| 桃花免费在线播放| 在线观看免费视频网站a站| 中文字幕最新亚洲高清| www.av在线官网国产| 亚洲欧美一区二区三区久久| 亚洲av日韩在线播放| 久久国产精品影院| 韩国精品一区二区三区| av福利片在线| 韩国精品一区二区三区| 狂野欧美激情性xxxx| 80岁老熟妇乱子伦牲交| 黄片播放在线免费| 亚洲av日韩精品久久久久久密 | av视频免费观看在线观看| 黄片播放在线免费| 成人免费观看视频高清| 无限看片的www在线观看| 国产成人精品无人区| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲av成人精品一二三区| 丁香六月天网| 国产成人啪精品午夜网站| 欧美激情 高清一区二区三区| 亚洲,一卡二卡三卡| 亚洲av电影在线进入| 老司机影院成人| 久久免费观看电影| 免费看av在线观看网站| 精品一区二区三区四区五区乱码 | 别揉我奶头~嗯~啊~动态视频 | 国产av精品麻豆| 国产老妇伦熟女老妇高清| 亚洲精品久久久久久婷婷小说| 中文字幕亚洲精品专区| 女人精品久久久久毛片| 成人午夜精彩视频在线观看| 亚洲欧美日韩高清在线视频 | 久久精品国产亚洲av涩爱| 精品少妇内射三级| 51午夜福利影视在线观看| 久久久久国产精品人妻一区二区| 中文欧美无线码| 欧美+亚洲+日韩+国产| 国产成人欧美在线观看 | 精品一区二区三区av网在线观看 | 亚洲国产精品国产精品| 成人三级做爰电影| 一本一本久久a久久精品综合妖精| 美女福利国产在线| 国产成人精品久久二区二区91| 精品少妇一区二区三区视频日本电影| 后天国语完整版免费观看| 韩国精品一区二区三区| 午夜激情久久久久久久| 99re6热这里在线精品视频| 亚洲国产最新在线播放| 男的添女的下面高潮视频| 老汉色∧v一级毛片| 国产97色在线日韩免费| 久久天堂一区二区三区四区| 久久国产亚洲av麻豆专区| 高清黄色对白视频在线免费看| 国产精品一二三区在线看| 亚洲欧美成人综合另类久久久| 国产真人三级小视频在线观看| 久久国产精品人妻蜜桃| 黑丝袜美女国产一区| 下体分泌物呈黄色| 麻豆国产av国片精品| 波多野结衣av一区二区av| 国产高清videossex| 性高湖久久久久久久久免费观看| 热99国产精品久久久久久7| 久久性视频一级片| av在线app专区| 欧美精品av麻豆av| 国产成人精品在线电影| 亚洲精品一卡2卡三卡4卡5卡 | 18禁裸乳无遮挡动漫免费视频| 人妻 亚洲 视频| 啦啦啦中文免费视频观看日本| 黄频高清免费视频| 另类精品久久| 欧美黄色淫秽网站| 欧美黄色片欧美黄色片| 咕卡用的链子| 久久99精品国语久久久| 新久久久久国产一级毛片| 国产人伦9x9x在线观看| 国产成人精品在线电影| 侵犯人妻中文字幕一二三四区| 国产成人av教育| 精品人妻在线不人妻| 国产免费视频播放在线视频| 欧美精品人与动牲交sv欧美| 国产欧美日韩一区二区三区在线| 2018国产大陆天天弄谢| 一二三四社区在线视频社区8| 国产免费一区二区三区四区乱码| 麻豆av在线久日| 中文字幕制服av| 熟女少妇亚洲综合色aaa.| 亚洲国产日韩一区二区| 亚洲精品日本国产第一区| 麻豆国产av国片精品| 国产真人三级小视频在线观看| 校园人妻丝袜中文字幕| 国产福利在线免费观看视频| 欧美日韩亚洲高清精品| 欧美中文综合在线视频| 日韩熟女老妇一区二区性免费视频| 国产片特级美女逼逼视频| 国产熟女午夜一区二区三区| 搡老岳熟女国产| 久久久久久久精品精品| 欧美精品一区二区大全| 91字幕亚洲| 国产精品欧美亚洲77777| 久久精品人人爽人人爽视色| 久久国产精品影院| 99热网站在线观看| 欧美老熟妇乱子伦牲交| 最近手机中文字幕大全| 麻豆国产av国片精品| 青春草视频在线免费观看| 国产1区2区3区精品| 电影成人av| 精品人妻1区二区| 在线观看国产h片| 少妇裸体淫交视频免费看高清 | 在线亚洲精品国产二区图片欧美| 熟女少妇亚洲综合色aaa.| 亚洲av日韩精品久久久久久密 | 久久人妻熟女aⅴ| 老司机靠b影院| 日韩av免费高清视频| 午夜福利免费观看在线| 欧美成人午夜精品| 亚洲精品乱久久久久久| 国产一区二区三区综合在线观看| 午夜激情av网站| 亚洲伊人久久精品综合| 宅男免费午夜| 不卡av一区二区三区| 日韩伦理黄色片| 国产av一区二区精品久久| 在线观看免费高清a一片| 亚洲精品久久久久久婷婷小说| 欧美av亚洲av综合av国产av| av国产精品久久久久影院| 亚洲国产欧美一区二区综合| 十八禁高潮呻吟视频| 色视频在线一区二区三区| 亚洲欧美精品自产自拍| 老司机靠b影院| 夫妻午夜视频| 欧美日韩亚洲综合一区二区三区_| 建设人人有责人人尽责人人享有的| 国产男人的电影天堂91| 久久精品人人爽人人爽视色| 丰满饥渴人妻一区二区三| 国产精品人妻久久久影院| 日韩熟女老妇一区二区性免费视频| 久久人人爽人人片av| 亚洲精品乱久久久久久| 国产一区二区激情短视频 | 蜜桃在线观看..| 国产精品久久久久成人av| bbb黄色大片| 久久亚洲国产成人精品v| 国产片特级美女逼逼视频| 精品一区二区三区四区五区乱码 | 国产精品欧美亚洲77777| 美女脱内裤让男人舔精品视频| 国产一级毛片在线| 欧美少妇被猛烈插入视频| 男女之事视频高清在线观看 |