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

    時序分析在制導(dǎo)炸彈測試數(shù)據(jù)判讀中的應(yīng)用

    2018-01-23 04:37張瑞劉懷勛魏志剛
    航空兵器 2018年5期
    關(guān)鍵詞:時間序列測試

    張瑞 劉懷勛 魏志剛

    摘 要:制導(dǎo)炸彈自動化測試是減少測試人員負(fù)擔(dān)、 降低研發(fā)成本、 提高測試效率的重要手段。 測試數(shù)據(jù)的判讀是自動測試系統(tǒng)的重要組成部分。 傳統(tǒng)測試數(shù)據(jù)判讀方法效率低下, 難以應(yīng)對制導(dǎo)炸彈大量測試數(shù)據(jù)的分析判讀需求。 針對這一問題, 本文利用MATLAB軟件, 完成了制導(dǎo)炸彈測試數(shù)據(jù)時間序列判讀方法設(shè)計, 并通過對舵機(jī)角度反饋偏差數(shù)據(jù)的判讀, 對該方法進(jìn)行了工程驗證。 驗證結(jié)果顯示, 設(shè)計的時間序列判讀方法能在先驗數(shù)據(jù)集基礎(chǔ)上, 對測試數(shù)據(jù)進(jìn)行異常檢測, 提高了測試數(shù)據(jù)的判讀效率和準(zhǔn)確性。

    關(guān)鍵詞: 制導(dǎo)炸彈; 時間序列; 數(shù)據(jù)判讀; 測試

    中圖分類號: TJ414文獻(xiàn)標(biāo)識碼: A 文章編號: 1673-5048(2018)05-0073-06[SQ0]

    0 引言

    制導(dǎo)炸彈的測試工作貫穿于其研制、 試驗、 生產(chǎn)、 使用等各個階段, 產(chǎn)生了大量的測試數(shù)據(jù)。 通過對這些測試數(shù)據(jù)的判讀, 測試人員能夠了解到炸彈各個部件的工作狀態(tài), 并能進(jìn)行異常的診斷。

    目前, 對測試數(shù)據(jù)的判讀多停留在依靠測試人員經(jīng)驗比較、 圖像趨勢觀察階段, 效率低下。 采用自動測試技術(shù)和先進(jìn)的數(shù)據(jù)判讀方法, 提高測試數(shù)據(jù)判讀的實時性、 準(zhǔn)確性和工作效率, 減輕測試人員工作負(fù)擔(dān), 是測試工作的迫切需求[1-4]。

    國內(nèi)不少學(xué)者對測試數(shù)據(jù)判讀方法進(jìn)行了研究。 吳限德等運用時間序列分析方法, 實現(xiàn)了對小衛(wèi)星測試數(shù)據(jù)的在線判讀[1]。 曹永梅等針對衛(wèi)星控制系統(tǒng)測試數(shù)據(jù), 基于設(shè)置的判讀規(guī)則方法, 設(shè)計了測試數(shù)據(jù)自動判讀系統(tǒng)[2]。 汪洋等對衛(wèi)星測試系統(tǒng)測試數(shù)據(jù)的變化規(guī)律進(jìn)行了研究, 研制了一種開放的判讀函數(shù)編譯平臺[3]。 蘇振華等利用衛(wèi)星下傳的測試數(shù)據(jù), 提出了一種基于BP神經(jīng)網(wǎng)絡(luò)的衛(wèi)星故障診斷方法[4]。 在這些已應(yīng)用的方法中, 時間序列分析方法可以簡單高效地建立數(shù)據(jù)數(shù)學(xué)模型, 能夠近似描述各種類型的數(shù)據(jù), 達(dá)到預(yù)測數(shù)據(jù)發(fā)展趨勢并診斷異常的目的, 在各個領(lǐng)域都得到了廣泛的應(yīng)用[5]。

    制導(dǎo)炸彈的測試數(shù)據(jù)具有數(shù)據(jù)量大、 數(shù)據(jù)隨環(huán)境影響大、 數(shù)據(jù)類型復(fù)雜多樣等特點。 本文運用時間序列分析理論, 借助MATLAB軟件工具, 對制導(dǎo)炸彈測試數(shù)據(jù)進(jìn)行了判讀, 完成了測試數(shù)據(jù)判讀方法的研究及應(yīng)用驗證。

    1 時間序列分析理論

    1.1 ARMA模型的基本原理

    時間序列模型包括自回歸模型(AR模型)、 滑動平均模型(MA模型)、 自回歸-滑動平均模型(ARMA模型)等。 其中ARMA模型中既包含自回歸部分也包含滑動平均部分, 這使模型在擬合數(shù)據(jù)時具有很大的靈活性[6]。 其表達(dá)式為

    z~t-φ1z~t-1-…-φpz~t-p=at-θ1at-1-…-θqat-q(1)

    式(1)簡記為ARMA(p, q)。 其中: p和q分別為自回歸部分和滑動平均部分的階數(shù);φi(i=1, 2, …, p)和θi(i=1, 2, …, q)分別為自回歸和滑動平均系數(shù)。特殊的, 若p=0, 則模型退化為MA(q); 若q=0, 則模型退化為AR(p)。

    模型類型可以通過數(shù)據(jù)的自相關(guān)函數(shù)(ACF)和偏相關(guān)函數(shù)(PACF)的拖尾性、 截尾性與模型類別的關(guān)系來判讀, 模型識別關(guān)系如表1所示[7]。

    1.2 DDS建模方法

    DDS(Dynamic Data System)方法是Pandit-Wu在1977年研究連續(xù)系統(tǒng)和離散系統(tǒng)的對應(yīng)關(guān)系時提出的[8]。 他們經(jīng)過研究發(fā)現(xiàn), 線性定常系統(tǒng)在白噪聲作用下的響應(yīng), 經(jīng)過均勻間隔采樣所得到的時間序列, 總可以用離散的ARMA(n, n-1)模型來擬合[6]。

    根據(jù)制導(dǎo)炸彈測試數(shù)據(jù)量大、 變化劇烈、 環(huán)境影響大等特點, 本文選用DDS方法來建立測試數(shù)據(jù)的時間序列模型。 DDS建模方法的一般流程如圖1所示[6]。

    1.3 模型階數(shù)的確定

    時間序列模型階數(shù)確定常用的方法有偏相關(guān)定階法、 白度檢驗定階法、 F檢驗法、 AIC、 BIC、 FPE準(zhǔn)則定階法等。 DDS建模方法選用F檢驗法來檢測模型階數(shù)。 F值計算公式為[7]

    F=Q1-Q02Q0N-p-q~F(2, N-p-q)(2)

    式中: Q0為模型ARMA(p, q)的殘差平方和;Q1為模型ARMA(p-1, q-1)的殘差平方和; N為樣本長度;p和q分別為自回歸和滑動平均的階數(shù);2是被檢驗的階數(shù)差數(shù)。 如果F>Fα, 則假設(shè)不成立, 模型階數(shù)仍有上升的可能。 否則ARMA(p-1, q-1)是合適的模型。 Fα從F分布表中查得。

    2 基于MATLAB軟件工具的時間序列建模

    2.1 MATLAB時間序列函數(shù)

    MATLAB時間序列處理函數(shù)主要放在控制類工具箱的系統(tǒng)辨識工具箱中[8]。 涉及時間序列分析的部分函數(shù)如表2所示。

    armaxm=armax(data, [na, nc])對于時間序列data是時間序列, na和nc分別是時間序列自回歸和滑動平均的階數(shù)

    predictyp=predict(m, data)按照模型m, 根據(jù)數(shù)據(jù)data, 進(jìn)行數(shù)據(jù)預(yù)測計算, 預(yù)測數(shù)據(jù)保存在yp中[9]

    residresid(m, data)計算置信度為99%的模型殘差圖[10]

    2.2 MATLAB時間序列建模步驟

    本文根據(jù)DDS建模過程, 運用MATLAB工具, 建立時間序列模型的步驟如下:

    (1) 導(dǎo)入原始數(shù)據(jù), 新建原始數(shù)據(jù)矩陣。

    (2) 使用detrend函數(shù)對原始數(shù)據(jù)進(jìn)行平穩(wěn)化和零均值化處理。

    (3) 運用相關(guān)分析函數(shù)autocorr和parcorr得到自相關(guān)和偏相關(guān)函數(shù)圖, 判斷時間序列屬于哪種典型模型, 進(jìn)而對通過DDS建模方法得到的時間序列模型進(jìn)行驗證。

    (4) 運用armax函數(shù)編寫程序, 建立各個階次的時間序列數(shù)學(xué)模型。

    (5) 運用編寫程序計算F值進(jìn)行F檢驗, 計算預(yù)測數(shù)據(jù)殘差平方和作模型評價參考。

    (6) 使用predict函數(shù)進(jìn)行預(yù)測數(shù)據(jù)計算, 并與原數(shù)據(jù)進(jìn)行比較, 檢查模型擬合情況。

    (7) 運用resid函數(shù)對預(yù)測數(shù)據(jù)進(jìn)行殘差分析, 對殘差序列進(jìn)行白噪聲檢驗計算, 驗證數(shù)學(xué)模型和原數(shù)據(jù)的擬合情況[11]。

    (8) 根據(jù)測試數(shù)據(jù)時間序列的模型階次、 模型參數(shù)、 殘差方差等能表征系統(tǒng)特性的參數(shù)情況, 與正常樣本的數(shù)據(jù)對比, 對測試數(shù)據(jù)異常情況進(jìn)行分析判斷。

    3 舵機(jī)轉(zhuǎn)角反饋偏差數(shù)據(jù)判讀

    3.1 時間序列模型識別

    制導(dǎo)炸彈舵機(jī)轉(zhuǎn)角反饋偏差數(shù)據(jù)是舵機(jī)控制轉(zhuǎn)角和舵機(jī)實際轉(zhuǎn)角反饋值之間的差值, 是隨時間動態(tài)變化的時間序列, 偏差幅度大小受舵機(jī)調(diào)整劇烈程度等因素影響, 是反映舵機(jī)跟隨控制系統(tǒng)命令情況的重要指標(biāo), 對測試舵機(jī)是否存在異常具有重要意義。 某次測試舵機(jī)轉(zhuǎn)角反饋偏差原始數(shù)據(jù)如圖2所示。

    舵機(jī)轉(zhuǎn)角偏差數(shù)據(jù)經(jīng)detrend函數(shù)預(yù)處理后為平穩(wěn)零均值時間序列, 用autocorr, parcorr計算的ACF和PACF圖如圖3~4所示。 由圖可以看出兩者均為拖尾的, 因此選用ARMA模型來擬合。

    3.2 建立數(shù)學(xué)模型

    3.2.1 初步建立模型并進(jìn)行F檢驗

    根據(jù)DDS建模步驟, 用armax函數(shù), 依次建立ARMA(2, 1)及ARMA(2n+2, 2n+1)模型, 并進(jìn)行F檢驗。 計算得到的F值如表3所示, 由表中可以看出在ARMA(28, 27)處時, F<Fα=2.37。

    3.2.2 對經(jīng)典DDS建模流程的修改

    為找到適合模型, 本文實際按照圖5所示的流程圖進(jìn)行數(shù)據(jù)建模。

    圖1中流程圖要求對模型參數(shù)φ2n,θ2n-1的置信區(qū)間進(jìn)行是否包含零的判斷。 實際建模過程中發(fā)現(xiàn)計算φ2n,θ2n-1的置信區(qū)間是困難的, 且DDS方法的流程最后步驟要求設(shè)定“其他模型”進(jìn)行檢驗, 確定最終模型。 而ARMA(2n-1, m), m≤2n-2模型和ARMA(2n, m), m≤2n-1模型就是很好的“其他模型”。

    因此, 首先用相關(guān)性分析判斷模型類型, 以便對DDS建模方法得到的模型進(jìn)行驗證。 當(dāng)使用DDS方法得到ARMA(2n, 2n-1)模型后, 直接對上述兩模型進(jìn)行建模和F檢驗。 這樣雖然增加了計算量, 但可以避免錯漏適合模型。 以現(xiàn)代計算機(jī)的速度和MATLAB算法的高效率, 增加的計算時間實際證明是很少的。

    時間序列建模和F檢驗的部分MATLAB程序如下:

    3.2.3 修改的DDS建模步驟應(yīng)用結(jié)果

    繼續(xù)建立ARMA(27, 26),ARMA(27, 25),ARMA(28, 26),ARMA(28, 25)等模型并與ARMA(28, 27)進(jìn)行F檢驗比較。 計算得到的F值如表4所示。 可見ARMA(28, 27)是合適模型。

    3.3 測試數(shù)據(jù)模型分析

    圖6顯示了殘差序列自相關(guān)系數(shù)圖, 根據(jù)白度檢驗定階法計算得到χ2k,α =34.38, Qk =29.89(顯著水平為α=0.01, k=25), χk,α >Qk即殘差序列是白噪聲。 圖7為反饋差值預(yù)測值與原始值擬合圖,可以看出預(yù)測模型能很好地擬合實際數(shù)據(jù)[12], 驗證了預(yù)測模型和建模方法的有效性。

    利用此方法對四枚制導(dǎo)炸彈實驗測試的16組舵機(jī)轉(zhuǎn)角偏差測試數(shù)據(jù)進(jìn)行時間序列建模, 得到的數(shù)據(jù)計算結(jié)果見表5。

    從表5計算結(jié)果可見, 同一炸彈四個舵機(jī)的轉(zhuǎn)角偏差平方和、 方差、 模型階次均相差不大, 而不同炸彈的數(shù)據(jù)相差較大。 這主要是由于不同炸彈工作工況不同導(dǎo)致的, 符合經(jīng)驗預(yù)期。

    在表5中, B彈第3通道舵機(jī)數(shù)據(jù)的時間序列模型階次(28, 27), 與其他舵機(jī)測試數(shù)據(jù)的階次小于10的情況相差較大, 因此初步判斷此舵機(jī)跟隨性能與其他舵機(jī)相比存在異常。 經(jīng)人工查看舵機(jī)偏差數(shù)據(jù)圖像發(fā)現(xiàn), 此舵機(jī)在剛上電時段跟隨偏差較大, 將此段數(shù)據(jù)刪除后重新建模分析, 計算的模型階次為(7, 6), 繼續(xù)刪除相同長度的數(shù)據(jù)序列繼續(xù)建模, 得到的模型階次依次為(4, 3), (6, 5), (8, 7), 與其他通道相差不大。 這說明偏差數(shù)據(jù)時間序列模型階次的大小與舵機(jī)跟隨偏差較大數(shù)據(jù)存在一定的相關(guān)性。

    表5中的B彈1通道模型階次(3, 2)、 C彈2通道模型階次(5, 4)、 D彈1通道模型階次(9, 8)、 D彈4通道模型階次(6, 4)均是通過驗證模型ARMA(2n-1, m), m≤2n-2和ARMA(2n, m), m≤2n-1得到的擬合模型階次, 這些結(jié)果說明了改進(jìn)DDS建模流程在實際建模過程中的有效性。

    4 總結(jié)

    本文針對制導(dǎo)炸彈測試數(shù)據(jù)自動判讀的需求, 運用時間序列分析理論, 借助MATLAB軟件工具, 建立了測試數(shù)據(jù)的時序數(shù)學(xué)模型, 對測試數(shù)據(jù)進(jìn)行了異常分析, 形成了包括判讀理論、 判讀實施工具和判讀結(jié)果分析的判讀方法。 實際判讀過程表明, 該方法能夠?qū)χ茖?dǎo)炸彈異常數(shù)據(jù)是否存在異常進(jìn)行對比判讀, 提高了制導(dǎo)炸彈測試數(shù)據(jù)判讀的效率和可靠性, 改變了以往僅能依靠數(shù)據(jù)趨勢圖和測試人員經(jīng)驗判讀數(shù)據(jù)的現(xiàn)狀, 為數(shù)據(jù)自動判讀系統(tǒng)的設(shè)計打下基礎(chǔ)。

    本文在實際運用DDS方法進(jìn)行建模過程中, 對DDS建模經(jīng)典流程進(jìn)行改進(jìn), 使其更便于使用計算機(jī)及數(shù)值分析工具軟件建模。 此修改流程可為今后運用DDS方法建模提供參考。

    參考文獻(xiàn):

    [1] 吳限德, 孫兆偉, 陳衛(wèi)東. 小衛(wèi)星測試數(shù)據(jù)在線時間序列判讀方法[J]. 哈爾濱工程大學(xué)學(xué)報, 2012, 33(3):383-388.

    Wu Xiande, Sun Zhaowei, Chen Weidong. An On Line Data Interpretation Method in a Micro Satellite Ground Test System Based on a Time Series Analysis[J].Journal of Harbin Engineering University, 2012, 33(3):383-388.(in Chinese)

    [2] 曹永梅, 周劍敏, 潘立鑫. 衛(wèi)星控制系統(tǒng)測試數(shù)據(jù)自動判讀系統(tǒng)的設(shè)計與實現(xiàn)[J]. 空間控制技術(shù)與應(yīng)用, 2013, 39(5): 58-62.

    Cao Yongmei, Zhou Jianmin, Pan Lixin. Design and Implementation of Automatic Data Discrimination System for Satellite Control System[J]. Aerospace Control and Application, 2013, 39(5): 58-62.(in Chinese)

    [3] 汪洋, 徐犇, 曾齊, 等. 基于函數(shù)解析的衛(wèi)星測試數(shù)據(jù)自動判讀[J]. 遙測遙控, 2016, 37(1):38-44.

    Wang Yang, Xu Ben, Zeng Qi, et al. Automatic Interpretation of Satellite Test Data Based on Analytic Function [J]. Journal of Telemetry, Tracking and Command, 2016, 37(1):38-44.(in Chinese)

    [4] 蘇振華, 陸文高, 奇晶, 等. 基于BP神經(jīng)網(wǎng)絡(luò)的衛(wèi)星故障診斷方法[J]. 計算機(jī)測量與控制, 2016, 24(5):63-66.

    Su Zhenhua, Lu Wengao, Qi Jing, et al. A Method of Satellite Fault Diagnosis Based on BP Neural Network[J]. Computer Measurement & Control, 2016, 24(5):63-66.(in Chinese)

    [5] 吳懷宇. 時間序列分析與綜合[M]. 武漢:武漢大學(xué)出版社, 2004.

    Wu Huaiyu. Time Series Analysis and Applications [M]. Wuhan:Wuhan University Press, 2004.(in Chinese)

    [6] 張樹京, 齊立心. 時間序列分析簡明教程[M]. 北京:清華大學(xué)出版社, 北方交通大學(xué)出版社, 2003.

    Zhang Shujing, Qi Lixin. Time Series Analysis Plain Course[M]. Beijing: Tsinghua University Press,Northern Jiaotong University Press, 2003.(in Chinese)

    [7] 周世瓊, 康龍云, 張彥寧, 等.基于MATLAB系統(tǒng)辨識工具箱的風(fēng)信號預(yù)測[J]. 太陽能學(xué)報, 2008, 29(4):417-421.

    Zhou Shiqiong, Kang Longyun, Zhang Yanning, et al. Wind Signal Forecasting Based on System Identification Toolbox of MATLAB[J]. Acta Energiae Solaris Sinica, 2008, 29(4):417-421.(in Chinese)

    [8] 欒穎. MATLAB R2013a工程分析與仿真[M]. 北京:清華大學(xué)出版社, 2014.

    Luan Ying. MATLAB R2013a Engineering Analysis and Simulation[M]. Beijing:Tsinghua University Press, 2014.(in Chinese)

    [9] 邵璠, 孫育河, 梁嵐珍. 基于時間序列法的風(fēng)電場風(fēng)速預(yù)測研究[J]. 華東電力, 2008, 36(7):26-29.

    Shao Fan, Sun Yuhe, Liang Lanzhen. Study on the Time Series Wind Speed Forecasting of the Wind Farm Based on Time Series[J]. East China Electric Power, 2008, 36(7):26-29.(in Chinese)

    [10] 韓路躍, 杜行檢. 基于MATLAB的時間序列建模與預(yù)測[J]. 計算機(jī)仿真, 2005, 22(4):105-107, 182.

    Han Luyue, Du Xingjian. Modeling and Prediction of Time Series Based on MATLAB[J].Computer Simulation, 2005, 22(4): 105-107, 182.(in Chinese)

    [11] 林卉. 非平穩(wěn)時間序列的建模方法研究[D].武漢:武漢理工大學(xué), 2005.

    Lin Hui. A Modeling Method Research on Nonstationary Time Series[D]. Wuhan:Wuhan University of Techno logy, 2005.(in Chinese)

    [12] 柴嘉薪, 王新龍, 王盾, 等. 一種光纖陀螺隨機(jī)振動誤差高精度建模方法[J]. 航空兵器, 2017(4): 49-54.

    Chai Jiaxin, Wang Xinlong, Wang Dun, et al. A High Precision Modeling Method for Random Vibration Error of Fiber Optic Gyroscope[J]. Aero Weaponry, 2017(4): 49-54.(in Chinese)

    Application of Time Series Analysis in the

    Guided Bomb Test Data Interpretation

    Zhang Rui, Liu Huaixun,Wei Zhigang

    (China Airborne Missile Academy, Luoyang 471009, China)

    Abstract:Automatic test of guided bomb is an important means of reducingoperators workload, reducing development costs, improving test efficiency. Data interpretation is one of the most important partsof automatictest systems. Traditional data interpretations have a low efficiency, therefore, it is difficult to satisfy the requirement of test systems with a large number of test data. To solve this problem, thetime series interpretation method of guided bomb test data is obtained using MATLAB software. Taking advantage of the steering engine angle feedback deviation data, the proposed method is verified. The results show that the interpretation method can detect data abnormalities comprehensively based on a prior data set, thereby improving the accuracy and efficiency of the test data.

    Key words:guided bomb; time series; data interpretation; testing

    猜你喜歡
    時間序列測試
    心理測試
    一年級上冊期末小測試
    逃離體能測試
    上證綜指收益率的影響因素分析
    基于指數(shù)平滑的電站設(shè)備故障時間序列預(yù)測研究
    基于時間序列的我國人均GDP分析與預(yù)測
    基于線性散列索引的時間序列查詢方法研究
    基于組合模型的能源需求預(yù)測
    心理小測試
    測試
    香蕉国产在线看| 久久影院123| 卡戴珊不雅视频在线播放| 日韩中文字幕欧美一区二区 | 亚洲精品久久成人aⅴ小说| 人人妻人人澡人人爽人人夜夜| 一区二区三区精品91| √禁漫天堂资源中文www| 最近手机中文字幕大全| 欧美精品高潮呻吟av久久| 久久久国产精品麻豆| 两性夫妻黄色片| 激情五月婷婷亚洲| 侵犯人妻中文字幕一二三四区| 少妇精品久久久久久久| 国产成人精品久久久久久| 菩萨蛮人人尽说江南好唐韦庄| xxxhd国产人妻xxx| 国产成人欧美| 黄色一级大片看看| 免费在线观看黄色视频的| 两个人看的免费小视频| 丝袜在线中文字幕| 伦理电影大哥的女人| 亚洲美女黄色视频免费看| 亚洲精品国产av蜜桃| 一本一本久久a久久精品综合妖精| 天天躁夜夜躁狠狠久久av| 国产男人的电影天堂91| 久久精品久久精品一区二区三区| 欧美久久黑人一区二区| 欧美日韩综合久久久久久| 国产成人一区二区在线| 国产探花极品一区二区| 美女大奶头黄色视频| 在线观看三级黄色| 女性被躁到高潮视频| av在线老鸭窝| 成人毛片60女人毛片免费| 亚洲国产欧美在线一区| 欧美日韩精品网址| kizo精华| 国产黄色免费在线视频| 人人澡人人妻人| 国产日韩欧美视频二区| 国产在线视频一区二区| bbb黄色大片| 日日摸夜夜添夜夜爱| 欧美 日韩 精品 国产| 99re6热这里在线精品视频| av在线观看视频网站免费| 国产人伦9x9x在线观看| 日韩中文字幕欧美一区二区 | 大陆偷拍与自拍| 一区二区av电影网| 夜夜骑夜夜射夜夜干| av在线app专区| 大陆偷拍与自拍| 色吧在线观看| 丝袜脚勾引网站| 一区二区日韩欧美中文字幕| 中文天堂在线官网| 在线观看国产h片| 在线免费观看不下载黄p国产| 午夜免费观看性视频| 日韩人妻精品一区2区三区| 欧美中文综合在线视频| 一区二区三区激情视频| a级片在线免费高清观看视频| 丝袜人妻中文字幕| 久久精品久久精品一区二区三区| 极品少妇高潮喷水抽搐| 天天操日日干夜夜撸| 午夜激情久久久久久久| 久久人人爽人人片av| 蜜桃在线观看..| 亚洲人成电影观看| 91老司机精品| 亚洲图色成人| 国产有黄有色有爽视频| 国产成人啪精品午夜网站| 久久久久精品性色| 久久毛片免费看一区二区三区| 婷婷色麻豆天堂久久| 丰满迷人的少妇在线观看| 一区二区三区激情视频| 免费高清在线观看日韩| 熟女少妇亚洲综合色aaa.| 欧美精品亚洲一区二区| 国产精品二区激情视频| av片东京热男人的天堂| 男人爽女人下面视频在线观看| 日韩电影二区| 超色免费av| 美女主播在线视频| 久热这里只有精品99| 黄色 视频免费看| 一级爰片在线观看| 亚洲av中文av极速乱| 欧美久久黑人一区二区| 欧美日韩亚洲高清精品| 丝袜美腿诱惑在线| 欧美精品av麻豆av| 伊人久久大香线蕉亚洲五| 欧美日韩综合久久久久久| 久久久久国产一级毛片高清牌| 欧美日韩av久久| 亚洲视频免费观看视频| 99香蕉大伊视频| 99久久综合免费| 在线观看人妻少妇| 中文天堂在线官网| 丁香六月欧美| 国产xxxxx性猛交| 1024视频免费在线观看| 人人澡人人妻人| 久久久久久人妻| 国产男女超爽视频在线观看| 国产精品欧美亚洲77777| 麻豆乱淫一区二区| 国产男女超爽视频在线观看| 免费久久久久久久精品成人欧美视频| 亚洲av综合色区一区| 岛国毛片在线播放| www.精华液| 亚洲精品在线美女| 叶爱在线成人免费视频播放| 久久ye,这里只有精品| 亚洲成人一二三区av| videosex国产| 又粗又硬又长又爽又黄的视频| 97人妻天天添夜夜摸| 国产成人一区二区在线| 精品一区二区三区四区五区乱码 | av卡一久久| tube8黄色片| 亚洲国产欧美一区二区综合| 国产1区2区3区精品| 午夜福利视频在线观看免费| 国产精品国产三级专区第一集| 天美传媒精品一区二区| 精品亚洲成a人片在线观看| 满18在线观看网站| 精品少妇久久久久久888优播| 99久久人妻综合| 男女免费视频国产| 欧美久久黑人一区二区| 少妇精品久久久久久久| 一区在线观看完整版| 9191精品国产免费久久| 久久毛片免费看一区二区三区| 亚洲天堂av无毛| 色播在线永久视频| 巨乳人妻的诱惑在线观看| 婷婷色麻豆天堂久久| 亚洲精品av麻豆狂野| 91成人精品电影| 国产精品久久久久久人妻精品电影 | 日韩成人av中文字幕在线观看| 丰满饥渴人妻一区二区三| 国产一区二区激情短视频 | 国产成人精品无人区| 这个男人来自地球电影免费观看 | 亚洲七黄色美女视频| 亚洲成人av在线免费| 免费高清在线观看视频在线观看| 久久久久精品久久久久真实原创| av网站在线播放免费| av国产精品久久久久影院| 欧美精品亚洲一区二区| av在线老鸭窝| 岛国毛片在线播放| 免费高清在线观看视频在线观看| 一级片'在线观看视频| 亚洲欧美一区二区三区久久| 欧美日韩亚洲高清精品| 亚洲国产精品999| 欧美在线一区亚洲| 久久久久国产一级毛片高清牌| 曰老女人黄片| 久久精品国产亚洲av涩爱| 十分钟在线观看高清视频www| 最近2019中文字幕mv第一页| 最近中文字幕2019免费版| www.精华液| 欧美日韩亚洲高清精品| 久久久亚洲精品成人影院| 久久青草综合色| 国产一区有黄有色的免费视频| 汤姆久久久久久久影院中文字幕| 大片免费播放器 马上看| 精品一区二区三卡| 欧美日韩视频精品一区| 97人妻天天添夜夜摸| 国产成人精品福利久久| 日韩中文字幕视频在线看片| 亚洲精品国产av蜜桃| 男女下面插进去视频免费观看| 天堂8中文在线网| 国产日韩欧美亚洲二区| 国产黄频视频在线观看| 午夜福利视频精品| 高清视频免费观看一区二区| 日韩av免费高清视频| 婷婷色av中文字幕| a级毛片在线看网站| 午夜福利视频精品| 国产成人一区二区在线| 国产一区亚洲一区在线观看| 嫩草影视91久久| av网站免费在线观看视频| 黄色毛片三级朝国网站| kizo精华| 国产97色在线日韩免费| 女的被弄到高潮叫床怎么办| 久久久亚洲精品成人影院| 91精品国产国语对白视频| 青春草视频在线免费观看| 亚洲天堂av无毛| 女性被躁到高潮视频| 亚洲国产日韩一区二区| av天堂久久9| 97在线人人人人妻| 看非洲黑人一级黄片| 亚洲一区中文字幕在线| 亚洲成人av在线免费| 在线观看国产h片| 丝袜喷水一区| 别揉我奶头~嗯~啊~动态视频 | 亚洲一级一片aⅴ在线观看| 黄色怎么调成土黄色| 亚洲美女黄色视频免费看| 亚洲国产欧美在线一区| 又黄又粗又硬又大视频| 久久鲁丝午夜福利片| 美女扒开内裤让男人捅视频| 女的被弄到高潮叫床怎么办| 成人手机av| 最近手机中文字幕大全| 十八禁网站网址无遮挡| 成人18禁高潮啪啪吃奶动态图| tube8黄色片| 日本午夜av视频| 在线 av 中文字幕| 99国产综合亚洲精品| 丰满迷人的少妇在线观看| 在线观看一区二区三区激情| av片东京热男人的天堂| 久久人人97超碰香蕉20202| 日韩一本色道免费dvd| 99九九在线精品视频| a级片在线免费高清观看视频| 中文字幕色久视频| 高清在线视频一区二区三区| 最黄视频免费看| 一级a爱视频在线免费观看| 成人亚洲欧美一区二区av| 久久影院123| 欧美97在线视频| 大香蕉久久成人网| 久热爱精品视频在线9| 最近中文字幕2019免费版| 人人妻人人添人人爽欧美一区卜| 少妇被粗大的猛进出69影院| 午夜激情久久久久久久| 精品亚洲乱码少妇综合久久| 亚洲精品,欧美精品| 国产老妇伦熟女老妇高清| av有码第一页| 亚洲av电影在线观看一区二区三区| 在线观看一区二区三区激情| 色吧在线观看| 97在线人人人人妻| 亚洲av成人精品一二三区| 久久精品久久久久久久性| 成人毛片60女人毛片免费| 国产免费一区二区三区四区乱码| 国产欧美亚洲国产| 最近中文字幕高清免费大全6| 赤兔流量卡办理| 日韩中文字幕视频在线看片| 亚洲美女搞黄在线观看| 别揉我奶头~嗯~啊~动态视频 | 精品一区二区三区av网在线观看 | 极品人妻少妇av视频| 亚洲精品国产区一区二| 亚洲人成电影观看| 亚洲欧美日韩另类电影网站| 狠狠婷婷综合久久久久久88av| 亚洲精品av麻豆狂野| 亚洲精品一二三| 大码成人一级视频| 一个人免费看片子| 成人国产av品久久久| 欧美人与性动交α欧美精品济南到| 9热在线视频观看99| 精品国产一区二区三区久久久樱花| 久久99热这里只频精品6学生| 久热这里只有精品99| 少妇的丰满在线观看| 亚洲专区中文字幕在线 | 亚洲国产精品国产精品| 精品人妻一区二区三区麻豆| 久久性视频一级片| 欧美成人午夜精品| 久久天躁狠狠躁夜夜2o2o | 日韩精品免费视频一区二区三区| 十八禁人妻一区二区| 91老司机精品| 久久免费观看电影| svipshipincom国产片| 国产精品免费视频内射| 亚洲第一青青草原| 99久久精品国产亚洲精品| 热re99久久精品国产66热6| 久久免费观看电影| 日韩 亚洲 欧美在线| 最近最新中文字幕大全免费视频 | 久久久精品区二区三区| av在线观看视频网站免费| 日韩 亚洲 欧美在线| 久久99一区二区三区| 人妻人人澡人人爽人人| 中文字幕另类日韩欧美亚洲嫩草| 女人爽到高潮嗷嗷叫在线视频| 啦啦啦 在线观看视频| 久久久久精品国产欧美久久久 | 老司机在亚洲福利影院| 高清黄色对白视频在线免费看| 欧美日韩国产mv在线观看视频| 岛国毛片在线播放| 少妇猛男粗大的猛烈进出视频| 深夜精品福利| 看免费av毛片| 亚洲国产日韩一区二区| 丁香六月欧美| 性高湖久久久久久久久免费观看| 91精品国产国语对白视频| 国产伦理片在线播放av一区| 黑人猛操日本美女一级片| 精品一品国产午夜福利视频| 国产精品三级大全| 国产又色又爽无遮挡免| 国产人伦9x9x在线观看| 免费在线观看完整版高清| 久久久久人妻精品一区果冻| 蜜桃在线观看..| 一级毛片电影观看| 国产成人午夜福利电影在线观看| 亚洲精品国产色婷婷电影| 午夜福利一区二区在线看| 又粗又硬又长又爽又黄的视频| 午夜影院在线不卡| 一边亲一边摸免费视频| 成人午夜精彩视频在线观看| 国产高清国产精品国产三级| 欧美黑人欧美精品刺激| 免费观看性生交大片5| 高清视频免费观看一区二区| xxxhd国产人妻xxx| 日韩熟女老妇一区二区性免费视频| av国产精品久久久久影院| 精品人妻熟女毛片av久久网站| 国产精品一区二区在线不卡| 深夜精品福利| 国产av国产精品国产| 捣出白浆h1v1| 中文字幕高清在线视频| 久久国产精品大桥未久av| netflix在线观看网站| 女人久久www免费人成看片| 一区二区三区四区激情视频| 国产成人精品在线电影| 丰满少妇做爰视频| 国产免费福利视频在线观看| 国产精品一区二区精品视频观看| 久久天躁狠狠躁夜夜2o2o | 黄片小视频在线播放| 成人18禁高潮啪啪吃奶动态图| 亚洲成人国产一区在线观看 | 国产一级毛片在线| bbb黄色大片| 日韩制服丝袜自拍偷拍| 一级毛片 在线播放| 免费黄频网站在线观看国产| 久久久久国产一级毛片高清牌| 午夜影院在线不卡| 多毛熟女@视频| 午夜福利一区二区在线看| 国产亚洲最大av| 免费看不卡的av| 成人18禁高潮啪啪吃奶动态图| 午夜福利,免费看| 精品第一国产精品| 黄色视频不卡| 欧美日韩视频高清一区二区三区二| 少妇精品久久久久久久| 精品免费久久久久久久清纯 | 女人爽到高潮嗷嗷叫在线视频| 日日啪夜夜爽| 亚洲av中文av极速乱| 人妻人人澡人人爽人人| 久久精品国产亚洲av涩爱| 久久久国产一区二区| 熟妇人妻不卡中文字幕| 国产精品久久久久久人妻精品电影 | 天美传媒精品一区二区| avwww免费| 性少妇av在线| 热re99久久精品国产66热6| 51午夜福利影视在线观看| 校园人妻丝袜中文字幕| 国产一卡二卡三卡精品 | 精品少妇一区二区三区视频日本电影 | 亚洲激情五月婷婷啪啪| 亚洲三区欧美一区| 婷婷色综合大香蕉| 十八禁高潮呻吟视频| 国产97色在线日韩免费| h视频一区二区三区| 国产高清国产精品国产三级| 国产欧美亚洲国产| 国产在视频线精品| 交换朋友夫妻互换小说| 国产亚洲一区二区精品| 操出白浆在线播放| 天天躁狠狠躁夜夜躁狠狠躁| av网站在线播放免费| 另类亚洲欧美激情| xxxhd国产人妻xxx| 最近中文字幕2019免费版| 日本欧美国产在线视频| 免费高清在线观看视频在线观看| 午夜精品国产一区二区电影| 十分钟在线观看高清视频www| 青春草视频在线免费观看| 亚洲精品国产色婷婷电影| 在线精品无人区一区二区三| 精品国产国语对白av| 亚洲国产精品一区二区三区在线| 免费黄色在线免费观看| 欧美日韩一级在线毛片| 国产亚洲av片在线观看秒播厂| 2018国产大陆天天弄谢| 亚洲国产精品999| 丁香六月欧美| 亚洲欧美一区二区三区黑人| 少妇人妻精品综合一区二区| 51午夜福利影视在线观看| 一本久久精品| 男女下面插进去视频免费观看| 在线免费观看不下载黄p国产| 亚洲熟女精品中文字幕| videos熟女内射| 亚洲国产中文字幕在线视频| 国产成人一区二区在线| 午夜91福利影院| 高清黄色对白视频在线免费看| 国产亚洲av高清不卡| 成人黄色视频免费在线看| 嫩草影视91久久| xxx大片免费视频| 国产精品国产av在线观看| 午夜免费鲁丝| 黄色一级大片看看| 欧美国产精品va在线观看不卡| 午夜日本视频在线| 伊人久久大香线蕉亚洲五| 国产一区亚洲一区在线观看| 亚洲成人国产一区在线观看 | av线在线观看网站| 国产片内射在线| 国产探花极品一区二区| 最近手机中文字幕大全| 国产精品国产三级专区第一集| 中文字幕亚洲精品专区| 爱豆传媒免费全集在线观看| 欧美老熟妇乱子伦牲交| 9191精品国产免费久久| 美女脱内裤让男人舔精品视频| 毛片一级片免费看久久久久| 久久国产精品男人的天堂亚洲| 亚洲成av片中文字幕在线观看| 欧美国产精品一级二级三级| 精品国产一区二区三区四区第35| 精品一区二区三区四区五区乱码 | 久久久久精品国产欧美久久久 | 18禁裸乳无遮挡动漫免费视频| 精品亚洲成a人片在线观看| 国产精品秋霞免费鲁丝片| 亚洲专区中文字幕在线 | 国产爽快片一区二区三区| 中文字幕人妻熟女乱码| 婷婷成人精品国产| 亚洲精品,欧美精品| 国产精品国产三级国产专区5o| 一级片免费观看大全| 丝袜脚勾引网站| 欧美精品人与动牲交sv欧美| 色婷婷久久久亚洲欧美| 制服丝袜香蕉在线| 热re99久久精品国产66热6| e午夜精品久久久久久久| 考比视频在线观看| 毛片一级片免费看久久久久| 欧美精品人与动牲交sv欧美| 午夜精品国产一区二区电影| 免费久久久久久久精品成人欧美视频| 国产一区二区三区综合在线观看| 午夜福利影视在线免费观看| 大香蕉久久成人网| 国产亚洲av片在线观看秒播厂| 亚洲欧洲精品一区二区精品久久久 | 男人操女人黄网站| 无遮挡黄片免费观看| 亚洲欧洲精品一区二区精品久久久 | 免费在线观看黄色视频的| 国产xxxxx性猛交| xxxhd国产人妻xxx| 久久99一区二区三区| 国产成人欧美| av网站在线播放免费| 亚洲精品美女久久久久99蜜臀 | av在线老鸭窝| 综合色丁香网| 亚洲第一区二区三区不卡| 韩国高清视频一区二区三区| 亚洲av男天堂| 久久av网站| 无限看片的www在线观看| 国产黄频视频在线观看| 欧美av亚洲av综合av国产av | 精品国产露脸久久av麻豆| 久久人人爽人人片av| 色婷婷av一区二区三区视频| 久久毛片免费看一区二区三区| 精品少妇内射三级| 18在线观看网站| 涩涩av久久男人的天堂| 国产精品女同一区二区软件| 免费观看av网站的网址| 国产一级毛片在线| 高清av免费在线| av有码第一页| 一级a爱视频在线免费观看| 日本wwww免费看| 国产不卡av网站在线观看| 亚洲国产最新在线播放| 亚洲精品,欧美精品| 天天躁狠狠躁夜夜躁狠狠躁| 精品少妇久久久久久888优播| 亚洲av日韩在线播放| 美女视频免费永久观看网站| 欧美日韩亚洲国产一区二区在线观看 | h视频一区二区三区| 在线观看免费高清a一片| 一级,二级,三级黄色视频| 精品亚洲乱码少妇综合久久| 成人漫画全彩无遮挡| 免费观看av网站的网址| 欧美 日韩 精品 国产| 日本猛色少妇xxxxx猛交久久| 日韩熟女老妇一区二区性免费视频| 久热爱精品视频在线9| 99久久精品国产亚洲精品| 18禁国产床啪视频网站| 高清欧美精品videossex| 亚洲天堂av无毛| 在线观看免费视频网站a站| 91精品伊人久久大香线蕉| 日韩av不卡免费在线播放| 十八禁高潮呻吟视频| 久久综合国产亚洲精品| 亚洲七黄色美女视频| 18禁观看日本| 久久久久精品国产欧美久久久 | 在线观看一区二区三区激情| 啦啦啦在线免费观看视频4| 日韩人妻精品一区2区三区| 97在线人人人人妻| 青青草视频在线视频观看| 成人毛片60女人毛片免费| 自线自在国产av| 亚洲av电影在线观看一区二区三区| 男女之事视频高清在线观看 | 欧美在线黄色| 在线观看免费高清a一片| 汤姆久久久久久久影院中文字幕| 国产一区有黄有色的免费视频| 国精品久久久久久国模美| 亚洲精品视频女| av线在线观看网站| 欧美日韩亚洲国产一区二区在线观看 | 亚洲 欧美一区二区三区| 亚洲人成网站在线观看播放| 人人妻人人添人人爽欧美一区卜| 国产成人精品在线电影| 男人爽女人下面视频在线观看| 超色免费av| 亚洲精品国产av蜜桃| 成人三级做爰电影| 亚洲久久久国产精品| 色播在线永久视频| 天天操日日干夜夜撸| 999精品在线视频| 一级毛片 在线播放| 久久久欧美国产精品| 国产1区2区3区精品| 免费女性裸体啪啪无遮挡网站| 一边摸一边做爽爽视频免费| 午夜老司机福利片| 国产xxxxx性猛交| 成人三级做爰电影| 亚洲精品一二三| 国产成人免费无遮挡视频|