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

    基于MRF模型的飛機(jī)飛行動(dòng)作識(shí)別劃分算法*

    2022-01-24 02:21:42顏廷龍王鳳芹
    關(guān)鍵詞:飛行數(shù)據(jù)馬爾可夫協(xié)方差

    顏廷龍,李 瑛,王鳳芹

    (1.海軍航空大學(xué)岸防兵學(xué)院,山東 煙臺(tái) 264001;2.海軍航空大學(xué)航空基礎(chǔ)學(xué)院,山東 煙臺(tái) 264001)

    1 引言

    飛行動(dòng)作的識(shí)別和劃分一直是飛機(jī)健康監(jiān)控、飛行模擬和飛行品質(zhì)評(píng)估等應(yīng)用的基礎(chǔ),快速準(zhǔn)確地識(shí)別飛機(jī)的基本飛行動(dòng)作和復(fù)雜飛行動(dòng)作具有重要意義。由于民航飛機(jī)動(dòng)作更加平穩(wěn),識(shí)別劃分難度小,相關(guān)的研究已經(jīng)比較成熟,而軍用戰(zhàn)機(jī)由于其機(jī)動(dòng)性強(qiáng)、速度快的特點(diǎn),其飛行動(dòng)作識(shí)別一直以來都是難點(diǎn)問題。

    目前,飛機(jī)動(dòng)作識(shí)別劃分主要分為以下幾種方法:一是根據(jù)專家的先驗(yàn)知識(shí),預(yù)先設(shè)定規(guī)則,建立知識(shí)庫識(shí)別飛機(jī)的機(jī)動(dòng)動(dòng)作。倪世宏等[1]分析了軍用飛機(jī)的飛行動(dòng)作的變化特征,建立了基于專家先驗(yàn)知識(shí)和飛行動(dòng)作特征的知識(shí)庫來識(shí)別飛機(jī)的飛行動(dòng)作,但這種方法的缺點(diǎn)在于面對(duì)不同的飛機(jī)型號(hào),設(shè)置的規(guī)則不具有通用性,另外,專家的先驗(yàn)知識(shí)可能帶有主觀誤差。二是利用目前的機(jī)器學(xué)習(xí)算法,比如貝葉斯網(wǎng)絡(luò)、支持向量機(jī)和神經(jīng)網(wǎng)絡(luò)監(jiān)督學(xué)習(xí)算法等。沈一超等[2]基于時(shí)間序列的DTW(Dynamic Time Warping)距離和貝葉斯網(wǎng)絡(luò)推理,提出了一種基于貝葉斯網(wǎng)絡(luò)的飛行動(dòng)作識(shí)別算法,可以識(shí)別復(fù)雜飛機(jī)動(dòng)作。這些方法要求預(yù)先將帶有基本飛行動(dòng)作標(biāo)記的飛行時(shí)間片段作為訓(xùn)練集[3 - 6],訓(xùn)練一個(gè)分類器,區(qū)分不同飛行數(shù)據(jù),要求事先人工標(biāo)記訓(xùn)練集,工作量較大,而且模型的通用性不強(qiáng)。

    馬爾可夫隨機(jī)場MRF(Markov Random Field)[7]模型是概率圖模型的一種,它可以表示時(shí)間序列中相鄰2個(gè)隨機(jī)變量的相關(guān)關(guān)系,在飛行數(shù)據(jù)的時(shí)間序列中,可反映相鄰時(shí)間點(diǎn)的參數(shù)的相關(guān)性。比如,飛機(jī)在t時(shí)刻的操縱參數(shù),會(huì)影響t+1時(shí)刻的姿態(tài)參數(shù)。因此,本文在飛行動(dòng)作識(shí)別領(lǐng)域,引入馬爾可夫隨機(jī)場模型,提出了一種基于馬爾可夫隨機(jī)場模型的飛機(jī)飛行動(dòng)作識(shí)別劃分算法,采用基于馬爾可夫隨機(jī)場的時(shí)序數(shù)據(jù)分割聚類算法,將飛行數(shù)據(jù)序列劃分為多個(gè)動(dòng)作類,并且用MRF網(wǎng)絡(luò)描述其每個(gè)動(dòng)作的特征。

    2 飛行動(dòng)作識(shí)別劃分模型

    2.1 馬爾可夫隨機(jī)場模型

    馬爾可夫隨機(jī)場MRF是在隨機(jī)場的基礎(chǔ)上添加馬爾可夫性質(zhì),使得隨機(jī)序列變量的分布僅與前一時(shí)刻有關(guān)。這一性質(zhì)提供了方便且具有一致性的建模方法,可以用來表示時(shí)間序列前后的約束關(guān)系。一維的馬爾可夫隨機(jī)場描述隨機(jī)序列中某一時(shí)刻的狀態(tài)只與前一時(shí)刻有關(guān),二維的馬爾可夫隨機(jī)場常常被用在圖像分割領(lǐng)域,將時(shí)域的馬爾可夫特性轉(zhuǎn)換到空間域,每個(gè)像元對(duì)于除它之外的相鄰的像元組成的鄰近集團(tuán)存在依賴性,通過這一性質(zhì),進(jìn)而有效地描述圖像的局部統(tǒng)計(jì)特征。在多維時(shí)間序列上定義MRF網(wǎng)絡(luò),將多維時(shí)間序列描述為一個(gè)由隨機(jī)變量組成的多層網(wǎng)絡(luò),相鄰層的網(wǎng)絡(luò)通過網(wǎng)絡(luò)節(jié)點(diǎn)的邊相連,網(wǎng)絡(luò)各個(gè)節(jié)點(diǎn)之間的邊描述各個(gè)變量之間的相關(guān)關(guān)系。區(qū)別于傳統(tǒng)飛行動(dòng)作分類方法依賴于飛行數(shù)據(jù)序列距離的度量[8],本文算法應(yīng)用MRF的特性,用MRF網(wǎng)絡(luò)表示飛行數(shù)據(jù)序列前后的依賴關(guān)系,從而得到更精確的飛行動(dòng)作分類結(jié)果。

    2.2 飛行數(shù)據(jù)時(shí)間序列分析

    飛行數(shù)據(jù)是一組離散的時(shí)間序列數(shù)據(jù),飛行動(dòng)作識(shí)別的準(zhǔn)確程度,很大程度取決于各個(gè)飛行動(dòng)作的劃分結(jié)果,本文提出基于馬爾可夫隨機(jī)場的時(shí)間序列分割算法,對(duì)飛行數(shù)據(jù)進(jìn)行分割和聚類,以完成對(duì)飛行動(dòng)作的劃分。設(shè)飛機(jī)的一組長度為T的飛行數(shù)據(jù)為x=[x1,x2,x3,…,xT],其中xi(i≤i≤T)為n維向量,表示同一時(shí)刻采集的n維飛行數(shù)據(jù)。設(shè)需要將待劃分的飛行數(shù)據(jù)分為K類,針對(duì)時(shí)間序列的聚類方式,不能僅對(duì)單獨(dú)數(shù)據(jù)點(diǎn)聚類,還要根據(jù)數(shù)據(jù)前后的相關(guān)變化規(guī)律進(jìn)行聚類,因此通過設(shè)置一個(gè)滑動(dòng)窗口w,其中w?T,w數(shù)值的選取取決于飛行數(shù)據(jù)的粒度和預(yù)期飛行動(dòng)作的長度。以t時(shí)刻為基準(zhǔn),向前截取窗口大小的數(shù)據(jù)片段,記為Xt:

    Xt=[xt-w+1,…,xt-1,xt]

    其中,t=1,2,…,T,Xt為n×w維向量,將原始飛行數(shù)據(jù)截取長度為w的多個(gè)子序列,將此時(shí)針對(duì)飛行數(shù)據(jù)的聚類變?yōu)閷?duì)長度為w的多個(gè)子序列的聚類,這樣在聚類過程中,相鄰的子序列更易被劃分為同一類,實(shí)現(xiàn)時(shí)間一致性的目標(biāo)。

    設(shè)需要將待劃分的飛行數(shù)據(jù)分為K類,屬于第j類的數(shù)據(jù)段集合記為Pj,其中j=1,2,…,K。每個(gè)類用高斯協(xié)方差矩陣定義,類的協(xié)方差逆矩陣Θj反映各個(gè)參數(shù)之間的獨(dú)立性。Θj是一個(gè)nw×nw矩陣,由w×w個(gè)子矩陣組成,每個(gè)子矩陣是大小為n×n的矩陣,位置pq上的子矩陣描述p時(shí)刻和q時(shí)刻之間,n個(gè)維度之間的協(xié)方差逆矩陣。由于飛行數(shù)據(jù)是非時(shí)變的,即每一個(gè)點(diǎn)的參數(shù)值只與相應(yīng)時(shí)間差有關(guān)。

    2.3 基于MRF模型的飛行動(dòng)作識(shí)別劃分算法

    MRF模型需要根據(jù)統(tǒng)計(jì)決策和估計(jì)理論中的最優(yōu)準(zhǔn)則確定時(shí)間序列分割問題的目標(biāo)函數(shù),采用一些優(yōu)化算法可求得滿足這些約束條件的MRF最大似然分布[9]。本文采用多元飛行數(shù)據(jù)時(shí)間序列各個(gè)維度之間的協(xié)方差逆矩陣Θj定義了MRF網(wǎng)絡(luò)的鄰接矩陣,通過估計(jì)稀疏的高斯逆協(xié)方差矩陣來學(xué)習(xí)每個(gè)聚類的MRF,MRF網(wǎng)絡(luò)具有多個(gè)層,層數(shù)對(duì)應(yīng)于定義的短子序列的窗口大小,通過求解帶約束的逆協(xié)方差逆矩陣估計(jì)問題,即可解決飛行動(dòng)作的分類問題。那么,針對(duì)飛行數(shù)據(jù)時(shí)間序列進(jìn)行劃分和識(shí)別變?yōu)榍蠼饷款愶w行動(dòng)作的參數(shù)之間的協(xié)方差逆矩陣Θ={Θ1,Θ2,…,ΘK}和每類飛行動(dòng)作的數(shù)據(jù)段集合P={P1,P2,…,PK}。

    用負(fù)對(duì)數(shù)似然函數(shù)表示飛行數(shù)據(jù)段Xt被歸為基本動(dòng)作j類的代價(jià),如式(1)所示:

    E(Xt∈Pj)=-ll(Xt,Θi)=

    (1)

    其中,ll(Xt,Θi) 表示對(duì)數(shù)似然函數(shù),μi表示矩陣Θi的均值,detΘi表示矩陣Θi的行列式。

    考慮飛行數(shù)據(jù)的連續(xù)性,相鄰時(shí)刻的數(shù)據(jù)段屬于不同類時(shí)施加懲罰項(xiàng)β,β越大,相鄰的飛行數(shù)據(jù)子序列被劃分為同一類動(dòng)作的可能性越大。當(dāng)β=0時(shí),飛行數(shù)據(jù)子序列可以被單獨(dú)劃分,連續(xù)性懲罰項(xiàng)表達(dá)式如式(2)所示:

    (2)

    由于極大似然估計(jì)不能產(chǎn)生稀疏解,導(dǎo)致模型的復(fù)雜度過高,不方便求解,需要加入稀疏性約束,增加正則化懲罰項(xiàng),如式(3)所示:

    λ‖Θi‖

    (3)

    其中,λ是正則化參數(shù)。稀疏化可以極大地簡化網(wǎng)絡(luò)結(jié)構(gòu),一定程度上降低了本文算法的復(fù)雜性,還可以提高模型的泛化能力,解決模型的過擬合問題。

    綜上所述,模型的總體優(yōu)化函數(shù)如式(4)所示:

    (4)

    3 飛行動(dòng)作識(shí)別劃分模型的求解方法

    MRF模型的總體優(yōu)化函數(shù)求解是一個(gè)組合和連續(xù)優(yōu)化問題,飛行動(dòng)作劃分參數(shù)的求解和飛行動(dòng)作類參數(shù)協(xié)方差逆矩陣的求解互相耦合,是高度非凸的優(yōu)化問題。解決這一問題的關(guān)鍵在于采用期望最大化EM(Expectation Maximum)算法將總體優(yōu)化目標(biāo)轉(zhuǎn)化為飛行動(dòng)作識(shí)別和飛行動(dòng)作劃分2個(gè)子問題,交替更新參數(shù),迭代求解,其中飛行動(dòng)作劃分參數(shù)使用動(dòng)態(tài)規(guī)劃算法求解,飛行動(dòng)作識(shí)別參數(shù)采用交替方向乘子法ADMM(Alternating Direction Method of Multipliers)求解[10]。

    3.1 飛行動(dòng)作劃分的參數(shù)求解

    針對(duì)飛行動(dòng)作劃分的參數(shù)Pj的求解問題,首先要給定Θj,此時(shí)Pj的優(yōu)化要考慮2個(gè)方面的問題,一是飛行數(shù)據(jù)段Xi被歸為j類的代價(jià),可以用負(fù)對(duì)數(shù)似然函數(shù)和表示,即式(1),另一個(gè)是飛行數(shù)據(jù)的連續(xù)性約束,應(yīng)用式(3)表示。2個(gè)代價(jià)構(gòu)成典型的流水線問題,可以采用動(dòng)態(tài)規(guī)劃算法進(jìn)行求解[11]。

    動(dòng)態(tài)規(guī)劃算法將T個(gè)子序列X1,…,XT分配到K個(gè)聚類的問題,等效為找到時(shí)間戳1~T的最小代價(jià)路徑,其中節(jié)點(diǎn)代價(jià)是將該數(shù)據(jù)段分配給飛行動(dòng)作的負(fù)對(duì)數(shù)似然函數(shù)和,并且每當(dāng)分類改變時(shí),邊的代價(jià)為β。

    3.2 飛行動(dòng)作識(shí)別參數(shù)矩陣的求解

    首先要給定一類中所有數(shù)據(jù)段的集合Pj,通過求最小化負(fù)對(duì)數(shù)似然總和,可以求解Θi,如式(5)所示:

    E1+E2

    (5)

    其中,E1=-|Pj|log det(Θi),E2可以寫成跡的形式,如式(6)所示:

    (6)

    其中,Si是經(jīng)驗(yàn)協(xié)方差,由當(dāng)前Pj所有數(shù)據(jù)段計(jì)算得到??紤]矩陣的稀疏性,添加一個(gè)正則項(xiàng)λ‖Θi‖,所以逆協(xié)方差逆矩陣的優(yōu)化函數(shù)如式(7)所示:

    (7)

    本文采用ADMM算法,ADMM是一種分布式凸優(yōu)化算法,主要應(yīng)用在大規(guī)模優(yōu)化任務(wù)中[12]。為了使問題更符合ADMM算法的形式,引入變量Z并將原問題重寫為等效問題,如式(8)所示。

    min(-log detΘ)+tr(SiΘ)+‖λ·Z‖1

    subject toΘ=Z,Z∈Γ

    (8)

    其增廣拉格朗日函數(shù)為:

    Lp(Θ,Z,U)=-log detΘ+tr(SiΘ)+

    (9)

    (10)

    (11)

    Uk+1:=Uk+Θk+1-Zk+1

    (12)

    其中k是迭代參數(shù)。

    3.3 EM算法求解整體優(yōu)化問題

    本文利用迭代EM算法[13]解決數(shù)據(jù)段分類和聚類參數(shù)協(xié)方差逆矩陣的問題。隨機(jī)初始化集群,并且交替執(zhí)行E-step和M-step,直到每個(gè)簇的分配已經(jīng)固定,模型已經(jīng)收斂。本文算法的求解步驟如下所示:

    步驟1初始化模型參數(shù)Pj,Θi;

    步驟2執(zhí)行E-step,給定Θi,應(yīng)用動(dòng)態(tài)規(guī)劃算法將飛行數(shù)據(jù)分到各個(gè)飛行動(dòng)作類,求解Pj;

    步驟3執(zhí)行M-step,給定Pj,應(yīng)用ADMM算法更新飛行動(dòng)作識(shí)別參數(shù)Θi;

    步驟4重復(fù)執(zhí)行步驟2和步驟3,不斷地迭代直到算法收斂。

    4 實(shí)驗(yàn)與結(jié)果分析

    4.1 實(shí)驗(yàn)方法

    為了方便驗(yàn)證本文算法的準(zhǔn)確性,采用海軍某型飛機(jī)的50次飛行任務(wù)記錄的飛行數(shù)據(jù)進(jìn)行實(shí)驗(yàn),每次飛行任務(wù)的時(shí)長大約為1 h,將每次飛行任務(wù)的數(shù)據(jù)劃分為10組,共500組。選取其中與飛行動(dòng)作識(shí)別相關(guān)的發(fā)動(dòng)機(jī)轉(zhuǎn)速、飛行速度、飛行高度、航向角、俯仰角和傾斜角6個(gè)參數(shù)。為了驗(yàn)證實(shí)驗(yàn)結(jié)果,飛行數(shù)據(jù)已經(jīng)事先采用人工標(biāo)記各個(gè)階段所屬的基本飛行動(dòng)作,飛行數(shù)據(jù)的采樣間隔均為1 s,算法的參數(shù)設(shè)置為:滑動(dòng)窗口大小w設(shè)置為5,根據(jù)飛行基本動(dòng)作的種類,將聚類簇的個(gè)數(shù)K設(shè)置為5。

    4.2 識(shí)別率對(duì)比

    由于飛機(jī)機(jī)動(dòng)平飛的過程占很大一部分,所以單純使用識(shí)別率不足以來評(píng)價(jià)算法的識(shí)別效果,本文使用Macro-F1 值評(píng)價(jià)各個(gè)算法的識(shí)別準(zhǔn)確率。Macro-F1值計(jì)算方式是先計(jì)算每一個(gè)類的Precison和Recall后,再計(jì)算各個(gè)類的F1值,然后將各個(gè)類的F1值的平均值,其中:

    (13)

    表1是本文算法(記為MRF)、基于高斯混和模型的聚類算法(記為GMM)[14],基于歐氏距離的DTW算法(記為DTW)[15]、基于自組織圖的人工神經(jīng)網(wǎng)絡(luò)聚類方法(記為Neural Gas)[16]和k-means算法針對(duì)測(cè)試樣本聚類的Macro-F1值和部分動(dòng)作識(shí)別準(zhǔn)確率。

    Table 1 Comparison of Macro-F1 values of various algorithms

    從表1可以看出,雖然5種算法對(duì)平飛的識(shí)別準(zhǔn)確率都很高,但DTW、k-means、Neural Gas對(duì)于斤斗類動(dòng)作的識(shí)別準(zhǔn)確率較低。斤斗類的飛行數(shù)據(jù)變化比較復(fù)雜,而這些算法大多使用基于距離的判定規(guī)則,無法識(shí)別飛行數(shù)據(jù)表示的本來特征。所以,本文算法在復(fù)雜動(dòng)作的表現(xiàn)上優(yōu)于其他基于距離的算法。

    接下來測(cè)試各個(gè)算法需要多少樣本才能準(zhǔn)確聚類,識(shí)別飛機(jī)動(dòng)作。圖1所示為5種算法的Macro-F1值與樣本數(shù)量的關(guān)系圖。如圖1所示,本文算法的性能明顯優(yōu)于其他算法,Macro-F1的值為0.91~0.98。與其他基于距離的聚類算法相比,本文算法只需更少的樣本就能達(dá)到相似的準(zhǔn)確率,而基于距離的聚類算法的結(jié)果準(zhǔn)確率不高,難以對(duì)飛機(jī)動(dòng)作進(jìn)行劃分。

    Figure 1 Relationship between Macro-F1 value and data volume圖1 Macro-F1值和數(shù)據(jù)量關(guān)系

    當(dāng)有100個(gè)樣本時(shí),沒有一種算法能夠準(zhǔn)確地識(shí)別飛機(jī)的動(dòng)作。但是,隨著樣本量的增多,基于本文算法的準(zhǔn)確率迅速提高。當(dāng)有200個(gè)樣本時(shí),本文算法的Macro-F1已經(jīng)超過0.9。

    4.3 模型參數(shù)分析

    (1)時(shí)間一致性參數(shù)。

    圖2所示為本文所提算法和去除算法中時(shí)間一致性約束MRF(β=0)的Marco-F1值和數(shù)據(jù)量之間的關(guān)系。由圖2可知,由β定義的時(shí)間一致性約束在數(shù)據(jù)量較小時(shí)僅具有很小的影響,因?yàn)镸RF和MRF(β=0)都獲得了相似的結(jié)果。但是,隨著樣本數(shù)量的增加,MRF(β=0時(shí))的Marco-F1徘徊在0.9附近。這意味著,一旦有足夠的樣本,飛行數(shù)據(jù)連續(xù)性約束這一參數(shù)的設(shè)置是提高算法性能的決定性因素。

    Figure 2 Influence of time consistency parameter 圖2 時(shí)間一致性參數(shù)影響

    (2)滑動(dòng)窗口大小分析。

    圖3所示為本文算法中滑動(dòng)窗口大小和Marco-F1值之間的關(guān)系,可以觀察到算法滑動(dòng)窗口大小設(shè)置為4~15可得到較高M(jìn)arco-F1聚類準(zhǔn)確率評(píng)分。只有在窗口大小降至4以下或升至10以上時(shí),算法準(zhǔn)確率較低。

    Figure 3 Influence of sliding window size 圖3 滑動(dòng)窗口大小影響

    5 結(jié)束語

    本文研究了基于MRF模型的飛機(jī)飛行動(dòng)作劃分和識(shí)別算法,算法在對(duì)飛行動(dòng)作準(zhǔn)確識(shí)別和劃分的同時(shí),還提供了結(jié)果的可解釋性。通過對(duì)比其他算法的Macro-F1值表明,本文算法提高了飛行動(dòng)作識(shí)別的準(zhǔn)確率。

    猜你喜歡
    飛行數(shù)據(jù)馬爾可夫協(xié)方差
    GARMIN1000系統(tǒng)維修平臺(tái)搭建預(yù)研究
    淺談飛行數(shù)據(jù)集中處理流程管理
    保費(fèi)隨機(jī)且?guī)в屑t利支付的復(fù)合馬爾可夫二項(xiàng)模型
    不確定系統(tǒng)改進(jìn)的魯棒協(xié)方差交叉融合穩(wěn)態(tài)Kalman預(yù)報(bào)器
    一種基于廣義協(xié)方差矩陣的欠定盲辨識(shí)方法
    基于廣域信息管理的飛行數(shù)據(jù)服務(wù)設(shè)計(jì)與實(shí)現(xiàn)
    基于SOP的核電廠操縱員監(jiān)視過程馬爾可夫模型
    應(yīng)用馬爾可夫鏈對(duì)品牌手機(jī)市場占有率進(jìn)行預(yù)測(cè)
    認(rèn)知無線網(wǎng)絡(luò)中基于隱馬爾可夫預(yù)測(cè)的P-CSMA協(xié)議
    縱向數(shù)據(jù)分析中使用滑動(dòng)平均Cholesky分解對(duì)回歸均值和協(xié)方差矩陣進(jìn)行同時(shí)半?yún)?shù)建模
    少妇人妻久久综合中文| 亚洲婷婷狠狠爱综合网| 麻豆乱淫一区二区| 亚洲国产精品999| 肉色欧美久久久久久久蜜桃| 1024香蕉在线观看| 国产免费福利视频在线观看| 高清av免费在线| 人妻人人澡人人爽人人| 国产成人系列免费观看| 国产成人精品久久二区二区91 | 欧美国产精品一级二级三级| 超碰成人久久| 丝袜在线中文字幕| 在线天堂中文资源库| 国产精品成人在线| 久久影院123| 丝袜脚勾引网站| 亚洲国产毛片av蜜桃av| 99精国产麻豆久久婷婷| 亚洲欧美成人综合另类久久久| 久久久久精品国产欧美久久久 | 亚洲一级一片aⅴ在线观看| 亚洲成人手机| videos熟女内射| 午夜福利免费观看在线| 菩萨蛮人人尽说江南好唐韦庄| 日日摸夜夜添夜夜爱| 国产乱来视频区| 80岁老熟妇乱子伦牲交| 色精品久久人妻99蜜桃| 国产日韩一区二区三区精品不卡| 亚洲精品国产区一区二| 伦理电影大哥的女人| 午夜日本视频在线| 国产高清国产精品国产三级| 亚洲三区欧美一区| 不卡视频在线观看欧美| 免费观看a级毛片全部| 久久人人爽av亚洲精品天堂| 久久精品久久久久久久性| 另类亚洲欧美激情| 国产精品成人在线| 国产成人免费无遮挡视频| 欧美亚洲日本最大视频资源| 精品酒店卫生间| 99香蕉大伊视频| xxxhd国产人妻xxx| 人人妻,人人澡人人爽秒播 | 成人手机av| 中国国产av一级| 日韩一区二区三区影片| 亚洲精品美女久久久久99蜜臀 | 精品一区二区三区av网在线观看 | 亚洲欧美一区二区三区黑人| 欧美国产精品一级二级三级| 热re99久久精品国产66热6| 亚洲国产精品999| 91国产中文字幕| 欧美久久黑人一区二区| 久久久久久久国产电影| 伦理电影大哥的女人| 亚洲欧美成人综合另类久久久| 精品少妇内射三级| 国产精品麻豆人妻色哟哟久久| 国产av码专区亚洲av| 各种免费的搞黄视频| 无限看片的www在线观看| videos熟女内射| 精品国产乱码久久久久久男人| 18禁裸乳无遮挡动漫免费视频| 欧美人与性动交α欧美软件| av视频免费观看在线观看| 天天躁夜夜躁狠狠久久av| 黑人猛操日本美女一级片| 少妇人妻精品综合一区二区| 久久久久精品人妻al黑| 国产一区二区在线观看av| 99九九在线精品视频| 十八禁高潮呻吟视频| 久久鲁丝午夜福利片| 精品国产一区二区三区久久久樱花| 久久这里只有精品19| 最近最新中文字幕免费大全7| 精品国产国语对白av| 日日啪夜夜爽| 中文字幕人妻丝袜制服| 丰满乱子伦码专区| 晚上一个人看的免费电影| 久久鲁丝午夜福利片| 亚洲一卡2卡3卡4卡5卡精品中文| 久久久国产一区二区| 久久青草综合色| 97人妻天天添夜夜摸| 一区二区三区激情视频| 亚洲成人免费av在线播放| a级片在线免费高清观看视频| 国产成人精品久久久久久| 亚洲人成电影观看| 亚洲av在线观看美女高潮| 熟妇人妻不卡中文字幕| 一边摸一边抽搐一进一出视频| 久久久久国产精品人妻一区二区| 极品少妇高潮喷水抽搐| 欧美激情 高清一区二区三区| 国产伦人伦偷精品视频| 69精品国产乱码久久久| 国产精品三级大全| 精品卡一卡二卡四卡免费| 亚洲三区欧美一区| 日韩 欧美 亚洲 中文字幕| 久久精品人人爽人人爽视色| 一区二区三区精品91| 精品久久蜜臀av无| 国产成人免费无遮挡视频| 国产日韩欧美视频二区| 国产精品麻豆人妻色哟哟久久| 欧美激情极品国产一区二区三区| 18在线观看网站| 国产男人的电影天堂91| 午夜日韩欧美国产| 欧美亚洲日本最大视频资源| 成人免费观看视频高清| 亚洲激情五月婷婷啪啪| av片东京热男人的天堂| 欧美日韩一区二区视频在线观看视频在线| 美女高潮到喷水免费观看| 狠狠婷婷综合久久久久久88av| 只有这里有精品99| 免费不卡黄色视频| 国产1区2区3区精品| 色播在线永久视频| 在线免费观看不下载黄p国产| 欧美在线黄色| 最近中文字幕2019免费版| 亚洲第一区二区三区不卡| 一区二区av电影网| 2018国产大陆天天弄谢| 中文精品一卡2卡3卡4更新| 搡老乐熟女国产| 一级a爱视频在线免费观看| av线在线观看网站| av视频免费观看在线观看| 久久久亚洲精品成人影院| 国产激情久久老熟女| 亚洲综合精品二区| 国产片特级美女逼逼视频| www.精华液| 日日爽夜夜爽网站| 青春草国产在线视频| 欧美黑人精品巨大| 亚洲av日韩在线播放| 一级毛片 在线播放| 国产av国产精品国产| 麻豆乱淫一区二区| 午夜福利一区二区在线看| 看十八女毛片水多多多| 国产熟女欧美一区二区| 80岁老熟妇乱子伦牲交| 夫妻性生交免费视频一级片| 精品国产一区二区三区久久久樱花| 街头女战士在线观看网站| 久久久久视频综合| 中文字幕高清在线视频| 蜜桃国产av成人99| 亚洲精华国产精华液的使用体验| 国产精品久久久久久人妻精品电影 | 成人午夜精彩视频在线观看| 欧美日韩亚洲综合一区二区三区_| 国产极品粉嫩免费观看在线| 成人漫画全彩无遮挡| 97精品久久久久久久久久精品| 精品国产一区二区三区四区第35| 亚洲av日韩在线播放| 精品久久蜜臀av无| 午夜91福利影院| 另类精品久久| 久久亚洲国产成人精品v| 满18在线观看网站| 久久久久视频综合| 美女福利国产在线| 晚上一个人看的免费电影| 亚洲精华国产精华液的使用体验| 美女大奶头黄色视频| xxx大片免费视频| 美国免费a级毛片| 黑人巨大精品欧美一区二区蜜桃| 精品少妇黑人巨大在线播放| 成人毛片60女人毛片免费| 中文精品一卡2卡3卡4更新| 亚洲自偷自拍图片 自拍| 日日啪夜夜爽| 女性生殖器流出的白浆| 99香蕉大伊视频| av国产精品久久久久影院| 性色av一级| 久久毛片免费看一区二区三区| 一级片免费观看大全| 一边亲一边摸免费视频| 色婷婷久久久亚洲欧美| 亚洲精品一二三| 色婷婷av一区二区三区视频| 亚洲,一卡二卡三卡| 日韩中文字幕欧美一区二区 | 亚洲成人一二三区av| 精品少妇黑人巨大在线播放| 国产熟女午夜一区二区三区| 黄色怎么调成土黄色| 97人妻天天添夜夜摸| 日韩一区二区视频免费看| 久久免费观看电影| 777米奇影视久久| 欧美人与性动交α欧美软件| 嫩草影院入口| 亚洲成人手机| 我要看黄色一级片免费的| 乱人伦中国视频| 啦啦啦 在线观看视频| 男女下面插进去视频免费观看| 中文天堂在线官网| 熟女av电影| 永久免费av网站大全| 国产精品亚洲av一区麻豆 | 午夜福利乱码中文字幕| 成人毛片60女人毛片免费| 亚洲av福利一区| 亚洲欧美一区二区三区国产| 欧美黑人精品巨大| 亚洲精品久久久久久婷婷小说| 中文字幕色久视频| 午夜老司机福利片| 夫妻午夜视频| 国产精品一区二区精品视频观看| 十八禁网站网址无遮挡| 久久久久视频综合| 亚洲精品成人av观看孕妇| 国产av一区二区精品久久| 欧美日韩成人在线一区二区| 菩萨蛮人人尽说江南好唐韦庄| 国产1区2区3区精品| 国产野战对白在线观看| 天天操日日干夜夜撸| 中国国产av一级| 亚洲av综合色区一区| 在线 av 中文字幕| 丰满少妇做爰视频| 天天躁夜夜躁狠狠躁躁| 黄色 视频免费看| 777久久人妻少妇嫩草av网站| 日韩视频在线欧美| a级毛片在线看网站| 大片电影免费在线观看免费| 天天躁日日躁夜夜躁夜夜| 自线自在国产av| 成人亚洲精品一区在线观看| 精品国产一区二区久久| 两个人看的免费小视频| 999精品在线视频| 老汉色av国产亚洲站长工具| 少妇被粗大猛烈的视频| 丰满少妇做爰视频| 天天影视国产精品| 人成视频在线观看免费观看| 亚洲国产中文字幕在线视频| 亚洲精品成人av观看孕妇| 亚洲,欧美,日韩| 日韩制服骚丝袜av| 欧美人与性动交α欧美软件| 一区福利在线观看| 亚洲欧洲国产日韩| 亚洲精品中文字幕在线视频| 最新的欧美精品一区二区| 免费在线观看视频国产中文字幕亚洲 | 久久精品久久精品一区二区三区| 人人妻人人爽人人添夜夜欢视频| 久久久精品94久久精品| 一区福利在线观看| 可以免费在线观看a视频的电影网站 | 18在线观看网站| 99久国产av精品国产电影| 久久久久久久久久久免费av| 黄色毛片三级朝国网站| 国产精品二区激情视频| 两个人看的免费小视频| 欧美激情 高清一区二区三区| 又粗又硬又长又爽又黄的视频| 丝瓜视频免费看黄片| 波多野结衣一区麻豆| 亚洲精品第二区| 在线看a的网站| 精品少妇内射三级| 久久久国产精品麻豆| 人人澡人人妻人| 日韩欧美一区视频在线观看| 亚洲综合色网址| 久久久久久久大尺度免费视频| 国产精品国产三级专区第一集| 老司机影院成人| 国产男女超爽视频在线观看| 久久国产亚洲av麻豆专区| 亚洲精品一二三| 国产精品免费大片| 亚洲精品久久成人aⅴ小说| 日日撸夜夜添| 亚洲专区中文字幕在线 | 亚洲成人免费av在线播放| 99热国产这里只有精品6| 啦啦啦啦在线视频资源| 久久精品亚洲av国产电影网| 99热全是精品| 亚洲精品久久久久久婷婷小说| 亚洲国产精品成人久久小说| 亚洲国产精品一区二区三区在线| 国产精品嫩草影院av在线观看| 欧美日韩综合久久久久久| 丁香六月欧美| 伦理电影大哥的女人| 天美传媒精品一区二区| 欧美人与善性xxx| 一级毛片电影观看| 久久毛片免费看一区二区三区| 高清av免费在线| 欧美精品高潮呻吟av久久| 色婷婷av一区二区三区视频| 涩涩av久久男人的天堂| 一级毛片 在线播放| 国产深夜福利视频在线观看| bbb黄色大片| 欧美日韩亚洲高清精品| 久久久久精品久久久久真实原创| 中文字幕制服av| 国产精品一国产av| 亚洲国产欧美日韩在线播放| 大香蕉久久成人网| 国产精品一二三区在线看| 亚洲精品视频女| 一边摸一边抽搐一进一出视频| 肉色欧美久久久久久久蜜桃| 免费日韩欧美在线观看| 亚洲国产精品一区二区三区在线| 美女国产高潮福利片在线看| 国产女主播在线喷水免费视频网站| 日韩免费高清中文字幕av| 日本黄色日本黄色录像| 最新在线观看一区二区三区 | h视频一区二区三区| 亚洲av欧美aⅴ国产| 啦啦啦中文免费视频观看日本| 午夜影院在线不卡| 亚洲天堂av无毛| 欧美国产精品va在线观看不卡| 亚洲人成网站在线观看播放| 日韩不卡一区二区三区视频在线| 曰老女人黄片| 久久亚洲国产成人精品v| 亚洲国产av影院在线观看| 久久久久国产精品人妻一区二区| 黄色视频不卡| www.精华液| 久久久久久人人人人人| 电影成人av| 国产免费视频播放在线视频| 男人添女人高潮全过程视频| 欧美日韩一级在线毛片| 免费日韩欧美在线观看| 日本黄色日本黄色录像| 久久久国产一区二区| 中文字幕制服av| 一级片免费观看大全| 毛片一级片免费看久久久久| 人妻 亚洲 视频| 菩萨蛮人人尽说江南好唐韦庄| 久久人人爽人人片av| 国产有黄有色有爽视频| 精品一区二区三区av网在线观看 | 美女午夜性视频免费| av福利片在线| 最近的中文字幕免费完整| 叶爱在线成人免费视频播放| 中文字幕最新亚洲高清| 国产亚洲av高清不卡| 黄色毛片三级朝国网站| 亚洲av成人不卡在线观看播放网 | 美女福利国产在线| 久久毛片免费看一区二区三区| 精品一区二区三区四区五区乱码 | 美女脱内裤让男人舔精品视频| 中国国产av一级| 久久国产精品男人的天堂亚洲| 在线观看免费午夜福利视频| 亚洲欧美一区二区三区国产| 国产精品av久久久久免费| 日本av手机在线免费观看| 日韩电影二区| 日韩伦理黄色片| 侵犯人妻中文字幕一二三四区| 国产在线视频一区二区| 交换朋友夫妻互换小说| 亚洲第一区二区三区不卡| 国产男女内射视频| 如何舔出高潮| 纵有疾风起免费观看全集完整版| 丁香六月欧美| 久久天堂一区二区三区四区| 女人精品久久久久毛片| 三上悠亚av全集在线观看| 欧美精品一区二区大全| av国产久精品久网站免费入址| 亚洲精品美女久久久久99蜜臀 | 一级片'在线观看视频| 嫩草影视91久久| 午夜av观看不卡| 国产精品免费视频内射| 51午夜福利影视在线观看| 卡戴珊不雅视频在线播放| 亚洲第一青青草原| av网站在线播放免费| 人人妻,人人澡人人爽秒播 | 男女免费视频国产| 久久久久国产精品人妻一区二区| 91aial.com中文字幕在线观看| 亚洲国产成人一精品久久久| 最新在线观看一区二区三区 | 国产一区二区三区综合在线观看| 亚洲婷婷狠狠爱综合网| 满18在线观看网站| 建设人人有责人人尽责人人享有的| 日本色播在线视频| 男女午夜视频在线观看| www.精华液| 尾随美女入室| 嫩草影院入口| 亚洲一码二码三码区别大吗| 秋霞在线观看毛片| 操出白浆在线播放| 一级爰片在线观看| 久久久精品94久久精品| 老汉色av国产亚洲站长工具| 精品久久久精品久久久| 国产精品 国内视频| 精品国产超薄肉色丝袜足j| 99热国产这里只有精品6| 老熟女久久久| 国产一区二区三区av在线| 18禁观看日本| 精品久久蜜臀av无| 亚洲第一青青草原| 色播在线永久视频| 人成视频在线观看免费观看| 另类精品久久| e午夜精品久久久久久久| 亚洲国产看品久久| 两个人看的免费小视频| 夫妻性生交免费视频一级片| 国产精品久久久久久精品电影小说| 亚洲精品久久久久久婷婷小说| 99热全是精品| 亚洲精品日韩在线中文字幕| 中文字幕色久视频| 男人爽女人下面视频在线观看| 国产熟女欧美一区二区| 中文字幕高清在线视频| 女性生殖器流出的白浆| 观看美女的网站| 久久毛片免费看一区二区三区| 欧美日韩亚洲国产一区二区在线观看 | 99香蕉大伊视频| 晚上一个人看的免费电影| 欧美精品一区二区免费开放| 秋霞在线观看毛片| 日韩大码丰满熟妇| 99热全是精品| 免费观看人在逋| 亚洲自偷自拍图片 自拍| 国产av码专区亚洲av| 精品久久久精品久久久| 香蕉丝袜av| 老司机影院成人| 亚洲久久久国产精品| 午夜免费鲁丝| 欧美激情极品国产一区二区三区| 午夜福利视频精品| 欧美最新免费一区二区三区| 久久精品国产综合久久久| 成年人免费黄色播放视频| 亚洲色图综合在线观看| 精品酒店卫生间| 亚洲国产欧美网| 老司机影院毛片| 欧美乱码精品一区二区三区| 精品国产露脸久久av麻豆| av电影中文网址| 亚洲第一区二区三区不卡| av女优亚洲男人天堂| 国产有黄有色有爽视频| 亚洲精品国产av成人精品| 在线天堂最新版资源| 欧美人与性动交α欧美软件| 亚洲欧美一区二区三区黑人| 老司机影院成人| 欧美精品人与动牲交sv欧美| 韩国av在线不卡| 91aial.com中文字幕在线观看| 国产黄色视频一区二区在线观看| 在线看a的网站| 国产精品免费视频内射| 麻豆av在线久日| 一区二区三区激情视频| 日本午夜av视频| 国产一级毛片在线| 波多野结衣一区麻豆| 三上悠亚av全集在线观看| 精品一区二区免费观看| 夜夜骑夜夜射夜夜干| 人成视频在线观看免费观看| 久久久久国产精品人妻一区二区| 国产黄频视频在线观看| 国产精品久久久久久精品古装| 日韩伦理黄色片| 国产亚洲av片在线观看秒播厂| 精品久久蜜臀av无| 精品视频人人做人人爽| 夫妻午夜视频| 亚洲,欧美,日韩| 午夜日韩欧美国产| 欧美日韩一区二区视频在线观看视频在线| 国产精品久久久久久精品电影小说| 久久天躁狠狠躁夜夜2o2o | 18禁动态无遮挡网站| 亚洲美女视频黄频| 一个人免费看片子| 亚洲精品国产一区二区精华液| 久久精品亚洲熟妇少妇任你| 国产有黄有色有爽视频| 亚洲四区av| 纵有疾风起免费观看全集完整版| 免费黄色在线免费观看| 香蕉国产在线看| 国产亚洲av高清不卡| 哪个播放器可以免费观看大片| 超碰97精品在线观看| 精品一区二区免费观看| 国产色婷婷99| www.精华液| 深夜精品福利| av在线播放精品| 午夜福利在线免费观看网站| 成人国语在线视频| 亚洲少妇的诱惑av| 亚洲一级一片aⅴ在线观看| 精品一品国产午夜福利视频| 久久久久人妻精品一区果冻| a级片在线免费高清观看视频| 欧美日韩视频高清一区二区三区二| 久久鲁丝午夜福利片| 国产免费一区二区三区四区乱码| 色94色欧美一区二区| 啦啦啦中文免费视频观看日本| 人妻人人澡人人爽人人| 久久精品人人爽人人爽视色| 色播在线永久视频| 中文字幕人妻丝袜制服| 男女高潮啪啪啪动态图| 亚洲五月色婷婷综合| 免费高清在线观看视频在线观看| 亚洲熟女毛片儿| 国产精品偷伦视频观看了| 免费看av在线观看网站| 国产av码专区亚洲av| 国产精品欧美亚洲77777| 一本一本久久a久久精品综合妖精| 午夜福利视频精品| 国产免费福利视频在线观看| 麻豆乱淫一区二区| 又大又黄又爽视频免费| 男人添女人高潮全过程视频| 欧美日韩视频高清一区二区三区二| 国产精品一区二区在线观看99| 老司机亚洲免费影院| 亚洲精品久久久久久婷婷小说| 狠狠精品人妻久久久久久综合| 18禁裸乳无遮挡动漫免费视频| 久久久精品免费免费高清| 一区二区三区四区激情视频| 国产精品.久久久| 欧美国产精品va在线观看不卡| 这个男人来自地球电影免费观看 | 久久女婷五月综合色啪小说| av电影中文网址| 男女边摸边吃奶| 亚洲国产欧美日韩在线播放| 免费不卡黄色视频| 日韩一区二区视频免费看| 中文字幕最新亚洲高清| 美女高潮到喷水免费观看| 一本大道久久a久久精品| 成人毛片60女人毛片免费| 黄网站色视频无遮挡免费观看| av福利片在线| 无遮挡黄片免费观看| 国产亚洲欧美精品永久| 尾随美女入室| 99久久综合免费| 亚洲精品aⅴ在线观看| 久久精品久久久久久久性| 精品一区二区三区四区五区乱码 | 久久精品久久久久久久性| 日本wwww免费看| 久久精品熟女亚洲av麻豆精品| 久久热在线av| 国产精品一区二区在线观看99| 看非洲黑人一级黄片| 18禁观看日本| 亚洲精品,欧美精品| 波多野结衣av一区二区av|