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

    陀螺隨機誤差補償中一種改進混合降噪法

    2022-12-01 13:23:40閻躍鵬鐘燕清李繼秀
    關(guān)鍵詞:方根陀螺濾波

    田 易,閻躍鵬,鐘燕清,李繼秀,孟 真

    (1.中國科學(xué)院微電子研究所,北京 100029;2.中國科學(xué)院大學(xué) 集成電路學(xué)院,北京 100049)

    微機電系統(tǒng)(Micro ElectroMechanical System,MEMS)陀螺是慣性導(dǎo)航系統(tǒng)中的核心器件之一,在姿態(tài)測量[1-2]領(lǐng)域有著廣泛的應(yīng)用。陀螺的測量精度直接影響姿態(tài)的解算精度,而隨機誤差是影響微機電系統(tǒng)陀螺測量精度的重要因素之一。目前關(guān)于陀螺隨機誤差的補償主要包括基于多尺度小波變換[3-4]的方法、基于卡爾曼濾波(Kalman Filter,KF)[5]的方法、基于前向預(yù)測濾波(Forward Linear Prediction,F(xiàn)LP)的方法[6-8]等?;诮?jīng)驗?zāi)B(tài)分解(Empirical Mode Decomposition,EMD)的方法于1998年由HUANG等提出[9],是一種新型的信號分解方法,適用于對非線性、非平穩(wěn)信號的降噪處理[10],已廣泛應(yīng)用于對微機電系統(tǒng)陀螺隨機誤差的消除[11-15]。但由于經(jīng)驗?zāi)B(tài)分解易產(chǎn)生模態(tài)混疊的問題[16],因此集成經(jīng)驗?zāi)B(tài)分解(Ensemble Empirical Mode Decomposition,EEMD)方法[17]、互補集成經(jīng)驗?zāi)B(tài)分解(Complementary Ensemble Empirical Mode Decomposition,CEEMD)方法[18]、自適應(yīng)噪聲的完整經(jīng)驗?zāi)B(tài)分解(Complete Ensemble Empirical Mode Decomposition with Adaptive Noise,CEEMDAN)方法[19]被相繼提出。但是EEMD和CEEMD方法都存在添加不同隨機噪聲導(dǎo)致分解本征模態(tài)函數(shù)數(shù)量不同的問題,而CEEMDAN通過附加自適應(yīng)白噪聲,不僅能減輕模態(tài)混疊現(xiàn)象,還可以精確重構(gòu)原始信號,可有效提高微機電系統(tǒng)陀螺信號的測量精度[20-21]。

    在常規(guī)CEEMDAN算法中,信號分解為一系列的本征模態(tài)函數(shù)(Intrinsic Mode Functions,IMFs)后,將低階噪聲本征模態(tài)函數(shù)去除且將第一階本征模態(tài)函數(shù)用于混合本征模態(tài)函數(shù)濾波的閾值確定。然而載體運動狀態(tài)突然改變而導(dǎo)致陀螺傳感器數(shù)據(jù)發(fā)生階躍變化時,低階噪聲IMFs中將混入高頻信號,一方面直接去除噪聲本征模態(tài)函數(shù)將引起數(shù)據(jù)高頻信號丟失,另一方面將導(dǎo)致濾波閾值提升引起對混合本征模態(tài)函數(shù)過度濾波問題。因此筆者提出了一種改進的CEEMDAN-FLP混合濾波算法,對噪聲本征模態(tài)函數(shù)和噪聲信號混合本征模態(tài)函數(shù)分別采用軟閾值(Soft Interval Thresholding,SIT)濾波和FLP進行噪聲消除的改進算法,在降噪的同時實現(xiàn)對高頻信號的最大保留,為慣性導(dǎo)航后處理提供高精度的陀螺測量數(shù)據(jù)。

    1 CEEMDAN-FLP基本原理

    1.1 CEEMDAN基本原理

    CEEMDAN是基于EEMD的一種改進算法,可有效解決經(jīng)驗?zāi)B(tài)分解過程中的模態(tài)混疊問題,與EEMD方法相比,有效減少了迭代次數(shù),增加了重構(gòu)精度。CEEMDAN的解算步驟如下:

    (1) 在原始信號中分別添加I組滿足均值為0、標(biāo)準(zhǔn)差為β0的隨機噪聲,構(gòu)成I組序列,其中第i組序列為x(n)+β0ωi(n),ωi(n)為滿足標(biāo)準(zhǔn)正態(tài)分布的隨機白噪聲。對添加噪聲后的序列分別進行EMD分解,則第1階IMF滿足

    (1)

    (2) 殘差滿足

    (2)

    (3) 第2階IMF滿足

    (3)

    (4) 以此遞推,第k階殘差滿足

    (4)

    (5) 第k+1階IMF滿足

    (5)

    (6) 重復(fù)至(4),直到rk滿足終止調(diào)減,因此趨勢項滿足

    (6)

    最終經(jīng)過CEEMDAN分解后,信號x(n)被分解為若干IMFs和一個趨勢項,滿足

    (7)

    1.2 FLP基本原理

    FLP濾波是一種適用于陀螺噪聲消除的濾波方法,其利用之前多個時刻的原始數(shù)據(jù)乘以對應(yīng)的權(quán)重值來預(yù)測當(dāng)前時刻的數(shù)據(jù),作為濾波器的輸出[7],即

    x(n)=WT(n)X(n-1) ,

    (8)

    其中,X(n-1)=[x(n-1)x(n-2) …x(n-N)]T,為x(n)之前的數(shù)據(jù),N為FLP濾波階數(shù);W(n)為權(quán)重系數(shù),滿足

    W(n+1)=W(n)+μ(n)e(n)X(n-1) ,

    (9)

    其中,e(n)為預(yù)測值與實測值的差值,μ(n)則為步長調(diào)整因子,用以控制自適應(yīng)過程的收斂速度和失調(diào)參數(shù)。

    1.3 軟閾值濾波原理

    (10)

    計算第1個IMF的噪聲方差σ1,通過中值估計法得到:

    (11)

    (12)

    各階IMF的噪聲閾值與第1個噪聲閾值的關(guān)系滿足

    (13)

    根據(jù)文獻[13]的結(jié)果,取H=0.5,β和ρ分別為0.719 和2.01。

    2 改進CEEMDAN-FLP混合濾波算法

    改進CEEMDAN-FLP混合濾波算法原理如圖1所示。對陀螺數(shù)據(jù)進行CEEMDAN分解,得到一系列頻率由高到低變化的窄帶IMFs信號和一個趨勢項。由于有用信號中多為低頻信號,噪聲信號中多為高頻信號,因此,隨著IMFs的階數(shù)增加,IMFs中包含的噪聲分量逐漸減小,信號分量逐漸增加。根據(jù)IMFs中包含噪聲分量和信號分量的多少,可以將IMFs區(qū)分為噪聲主導(dǎo)IMFs、信號噪聲混合IMFs和信號主導(dǎo)IMFs,對不同類型的IMFs采取不同的處理方式,達到對信號濾波降噪的目的。

    圖1 改進CEEMDAN-FLP混合濾波算法框圖

    首先,采用連續(xù)均方誤差分析(Consecutive Mean Square Errors,CMSE)來區(qū)分噪聲IMFs與混合IMFs,該方法不需要任何外部參考信息[11]。在噪聲IMFs和混合IMFs分類判別中,將連續(xù)的兩個重構(gòu)信號作比較,來判斷這兩個信號的相似程度。

    (14)

    因此,CMSE值的大小反映了連續(xù)的兩個重構(gòu)信號的相似程度,也反映了各階IMFs中的信號能量多少。隨著IMFs階數(shù)增加,CMSE的值呈先減小后增大的變化規(guī)律。CMSE值先減小,是由于低階IMFs中幾乎不含信號,而噪聲分量隨著階數(shù)的增加而降低;CMSE后增大,是由于IMFs中信號分量的出現(xiàn),并隨階數(shù)增長逐漸增加。因此,將CMSE取值最小時的IMFs階數(shù),作為噪聲IMFs與混合IMFs的分界點M1,滿足[11]

    (15)

    D(i)=dist[PDF(x(t)),PDF(IMFi(t))], 1≤i≤n。

    (16)

    M2=arg max[D(i)]+1 。

    (17)

    將IMFs區(qū)分為噪聲IMFs、混合IMFs和信號IMFs后,提出了一種改進CEEMDAN-FLP混合濾波算法。改進算法對CEEMDAN與FLP混合濾波進行了嘗試,并采用新的濾波方式,對噪聲IMFs進行軟閾值濾波,避免常規(guī)算法將噪聲IMFs直接剔除引起的高頻信號丟失問題;對混合IMFs進行FLP濾波,避免由于閾值提升引起的過度濾波問題。最后,將經(jīng)過軟閾值濾波的噪聲IMFs、經(jīng)過FLP濾波的混合IMFs與信號IMFs和趨勢項進行數(shù)據(jù)重構(gòu),得到最終的降噪信號。改進算法可在消除噪聲的同時,更有效地保留信號中的高頻分量。

    3 仿真驗證與分析

    模擬載體運動狀態(tài)突然改變的情況,在理想的階躍信號上添加均值為0、均方差為2°/s,滿足正態(tài)分布的隨機噪聲作為陀螺信號,對其進行CEEMDAN分解,并根據(jù)CMSE法和IMF與原始信號的2階范數(shù)分析法分別確定IMFs分類,其噪聲IMFs與混合IMFs的分界如圖2所示。由于3階IMFs的CMSE值最小,因此,噪聲IMFs與混合IMFs的分界為M1=3。信號IMFs與混合IMFs的分界如圖3所示,混合IMFs與信號IMFs的分界為M2=6。

    圖2 噪聲IMFs與混合IMFs的分界

    圖3 信號IMFs與混合IMFs的分界

    1~3階IMFs被劃分為噪聲IMFs,4~6階IMFs被劃分為信號噪聲混合IMFs,而7階以上則為信號IMFs。對噪聲IMFs采用軟閾值濾波(以3階IMF為例),濾波結(jié)果及19 s處的細(xì)節(jié)放大圖如圖4所示。由圖可見,噪聲IMFs中明顯包含著有用的高頻信號信息,通過軟閾值濾波,實現(xiàn)了對噪聲的濾除,同時保留了有用的高頻信號。對混合IMFs采用FLP濾波(以4階IMF為例),濾波結(jié)果及19 s處的細(xì)節(jié)放大圖如圖5所示。由圖可見,混合IMFs中的有效信號得到了充分保留,噪聲得到了明顯抑制。

    圖4 噪聲IMF3的軟閾值濾波前后對比

    圖5 混合IMF4的FLP濾波前后對比

    將濾波后的噪聲IMFs和混合IMFs與信號IMFs進行數(shù)據(jù)重構(gòu),得到濾波降噪后的陀螺數(shù)據(jù),將改進CEEMDAN-FLP濾波結(jié)果與原始數(shù)據(jù)、EMD濾波結(jié)果、CEEMDAN濾波結(jié)果及未改進CEEMDAN-FLP濾波結(jié)果進行仿真對比,結(jié)果及29 s和52 s處的細(xì)節(jié)放大圖如圖6所示。從圖6中可以看到,改進算法在信號突變時,具有最優(yōu)的動態(tài)跟蹤能力,CEEMDAN-FLP算法次之,而傳統(tǒng)的EMD數(shù)據(jù)延時最大,信號失真。

    圖6 仿真信號降噪前后對比圖及細(xì)節(jié)放大圖

    與仿真理想數(shù)據(jù)誤差統(tǒng)計對比結(jié)果如表1所示。改進的CEEMDAN-FLP濾波結(jié)果的均方根誤差最小,僅為0.95°/s,而原始信號的均方根誤差為1.96°/s,均方根誤差較濾波前相比降低了約51.53%;基于EMD分解的濾波結(jié)果均方根誤差為1.15°/s,改進算法與EMD濾波算法相比均方根誤差降低了約17.39%;基于CEEMDAN分解的濾波結(jié)果均方根誤差為1.12°/s,改進算法與CEEMDAN算法相比均方根誤差降低了約15.18%。

    表1 濾波前后誤差統(tǒng)計對比 (°)/s

    4 實驗驗證與分析

    采集慣性導(dǎo)航模塊的傳感器數(shù)據(jù),對仿真結(jié)果作進一步驗證。陀螺采用ADXRS 645,傳感器先處于靜止?fàn)顟B(tài),然后繞滾轉(zhuǎn)軸以60°/s的速率旋轉(zhuǎn)一周,靜止20 s左右后,再逆向以60°/s的速率旋轉(zhuǎn)一周,回到初始位置。實采數(shù)據(jù)通過改進CEEMDAN-FLP算法與EMD、CEEMDAN及CEEMDAN-FLP算法的濾波結(jié)果及細(xì)節(jié)放大圖如圖7所示。改進算法在信號突變處的延時最小,基于CEEMDAN-FLP算法和CEEMDAN算法在突變處的延時依次增大,而基于EMD的濾波算法信號延時最大,產(chǎn)生了信號失真。

    圖7 實采信號降噪前后對比圖及細(xì)節(jié)放大圖

    陀螺濾波前后數(shù)據(jù)與轉(zhuǎn)臺運動角速率誤差統(tǒng)計對比結(jié)果如表2所示。改進的CEEMDAN-FLP濾波結(jié)果的均方根誤差最小,僅為1.01°/s?;贑EEMDAN-FLP、CEEMDAN和EMD的濾波結(jié)果均方根誤差依次增大。改進算法與原始信號相比,均方根誤差降低了約57.02%;與基于EMD分解的濾波算法相比,均方根誤差降低了約50.24%;與基于CEEMDAN的濾波算法相比,均方根誤差降低了約38.04%;與基于CEEMDAN- FLP濾波算法相比,均方根誤差降低了約35.67%。

    表2 實采數(shù)據(jù)濾波前后誤差統(tǒng)計對比 (°)/s

    慣導(dǎo)模塊正向旋轉(zhuǎn)一周,在43 s處,理論姿態(tài)角應(yīng)為360°,然后逆向旋轉(zhuǎn)一周,在57 s處,理論姿態(tài)應(yīng)為0°,對比分析正旋轉(zhuǎn)(43 s處)和反旋轉(zhuǎn)時(57 s處)改進CEEMDAN-FLP算法與EMD、CEEMDAN及CEEMDAN-FLP算法的濾波結(jié)果解算后的滾轉(zhuǎn)角解算結(jié)果如圖8所示。改進的CEEMDAN-FLP濾波算法解算的滾轉(zhuǎn)角累積誤差最小。

    圖8 姿態(tài)解算對比圖及細(xì)節(jié)放大圖

    滾轉(zhuǎn)角誤差對比結(jié)果如表3所示。正轉(zhuǎn)和反轉(zhuǎn)結(jié)束后,改進CEEMDAN-FLP算法得到的滾轉(zhuǎn)角誤差為0.44°,誤差最?。换贑EEMDAN-FLP濾波得到的滾轉(zhuǎn)角誤差為1.32°;基于CEEMDAN濾波得到的滾轉(zhuǎn)角誤差為2.14°。而基于EMD分解的濾波算法由于模態(tài)混疊引起信號失真,無法真實反映陀螺數(shù)據(jù)的變化情況,使得最終的滾轉(zhuǎn)角解算結(jié)果甚至不及原始數(shù)據(jù)的滾轉(zhuǎn)角解算結(jié)果。改進算法的姿態(tài)累積誤差僅約為CEEMDAN濾波算法的20.56%,有效降低了姿態(tài)誤差。

    表3 姿態(tài)誤差對比 (°)

    在處理器為Intel Core i7-8550U、操作系統(tǒng)為Windows 10的筆記本電腦上,對實測陀螺數(shù)據(jù)進行處理,各個算法的運行時間對比如表4所示?;贓MD的濾波算法雖然耗時最短,但是模態(tài)混疊的問題限制了其應(yīng)用;而基于CEEMDAN算法、CEEMDAN-FLP濾波算法和改進CEEMDAN-FLP濾波算法的時間消耗相當(dāng),通過仿真和實測數(shù)據(jù)驗證,改進的CEEMDAN-FLP濾波算法的濾波降噪效果最好。

    表4 算法運行時間對比 s

    5 總 結(jié)

    在慣性導(dǎo)航數(shù)據(jù)后處理中,針對當(dāng)載體運動狀態(tài)突然改變,導(dǎo)致MEMS陀螺傳感器數(shù)據(jù)突變的情況,筆者提出了一種改進的CEEMDAN-FLP混合降噪方法。首先通過CEEMDAN算法將信號分解為多個IMFs 和一個趨勢項,并將IMFs劃分為噪聲主導(dǎo)IMFs、信號噪聲混合IMFs及信號主導(dǎo)IMFs。改進算法考慮了噪聲IMFs的影響,對噪聲IMFs和噪聲信號混合IMFs分別采用軟閾值(SIT)濾波和前向線性預(yù)測濾波(FLP)進行噪聲消除,最后對濾波后的IMFs與信號IMFs及趨勢項進行數(shù)據(jù)重構(gòu),在降噪的同時實現(xiàn)對高頻信號的最大保留。將改進的CEEMDAN-FLP算法分別與基于EMD的濾波算法、基于CEEMDAN的濾波算法及基于CEEMDAN-FLP的濾波算法比較。通過仿真驗證,改進算法濾波結(jié)果的均方根誤差與濾波前相比減小了約51.53%,與基于EMD的濾波算法相比減小了約17.39%;通過實測數(shù)據(jù)驗證,改進算法濾波后的陀螺數(shù)據(jù)與基于CEEMDAN的算法濾波后的陀螺數(shù)據(jù)分別用于姿態(tài)解算,在不明顯增加運算負(fù)擔(dān)的同時,改進算法姿態(tài)累積誤差僅約是CEEMDAN算法姿態(tài)累積誤差的20.56%,可見,改進算法可以有效提升姿態(tài)測量精度。

    猜你喜歡
    方根陀螺濾波
    方根拓展探究
    做個紙陀螺
    玩陀螺
    陀螺轉(zhuǎn)轉(zhuǎn)轉(zhuǎn)
    軍事文摘(2018年24期)2018-12-26 00:58:18
    我最喜歡的陀螺
    快樂語文(2018年36期)2018-03-12 00:56:02
    均方根嵌入式容積粒子PHD 多目標(biāo)跟蹤方法
    揭開心算方根之謎
    RTS平滑濾波在事后姿態(tài)確定中的應(yīng)用
    基于線性正則變換的 LMS 自適應(yīng)濾波
    遙測遙控(2015年2期)2015-04-23 08:15:18
    數(shù)學(xué)魔術(shù)
    亚洲午夜理论影院| 日韩 欧美 亚洲 中文字幕| 国产视频一区二区在线看| 99久久九九国产精品国产免费| av专区在线播放| 男女那种视频在线观看| 中文在线观看免费www的网站| 特大巨黑吊av在线直播| 国产成人aa在线观看| 黄片大片在线免费观看| 俺也久久电影网| 母亲3免费完整高清在线观看| 国产精品女同一区二区软件 | 日本黄色片子视频| 成人鲁丝片一二三区免费| 桃色一区二区三区在线观看| 麻豆国产av国片精品| 国产亚洲精品久久久com| 亚洲成av人片免费观看| 久久久国产成人免费| 狂野欧美白嫩少妇大欣赏| 国产欧美日韩一区二区精品| 中亚洲国语对白在线视频| 国产精品久久视频播放| 久久精品91无色码中文字幕| 亚洲精华国产精华精| 国产精品 国内视频| 久久亚洲真实| 中亚洲国语对白在线视频| 天天躁日日操中文字幕| 最近视频中文字幕2019在线8| 欧美性猛交╳xxx乱大交人| 久久久久久久午夜电影| 亚洲狠狠婷婷综合久久图片| 叶爱在线成人免费视频播放| 一夜夜www| 69av精品久久久久久| 国产精品精品国产色婷婷| 琪琪午夜伦伦电影理论片6080| 99精品欧美一区二区三区四区| 长腿黑丝高跟| 听说在线观看完整版免费高清| 国产黄色小视频在线观看| 两个人视频免费观看高清| 一边摸一边抽搐一进一小说| 精品免费久久久久久久清纯| 极品教师在线免费播放| av视频在线观看入口| 狠狠狠狠99中文字幕| 成人午夜高清在线视频| 国产高清有码在线观看视频| 亚洲av熟女| 在线观看午夜福利视频| 国产高清有码在线观看视频| а√天堂www在线а√下载| 亚洲av第一区精品v没综合| 天美传媒精品一区二区| 最近视频中文字幕2019在线8| a级毛片a级免费在线| 黄色视频,在线免费观看| 搡老妇女老女人老熟妇| 啪啪无遮挡十八禁网站| 国产精品 国内视频| 99riav亚洲国产免费| 青草久久国产| 国内精品久久久久久久电影| 99久国产av精品| 99久久综合精品五月天人人| 特大巨黑吊av在线直播| 欧美另类亚洲清纯唯美| 法律面前人人平等表现在哪些方面| 日本在线视频免费播放| 老司机福利观看| 尤物成人国产欧美一区二区三区| 精品久久久久久久久久久久久| 色尼玛亚洲综合影院| 久久精品国产99精品国产亚洲性色| 亚洲av不卡在线观看| 国产高清videossex| 欧美一级a爱片免费观看看| 亚洲成av人片在线播放无| 亚洲一区二区三区不卡视频| 久久久久久久久中文| 国产伦在线观看视频一区| 伊人久久精品亚洲午夜| 99视频精品全部免费 在线| 日韩欧美国产一区二区入口| 精华霜和精华液先用哪个| 最新中文字幕久久久久| 免费电影在线观看免费观看| 国产91精品成人一区二区三区| 亚洲国产欧洲综合997久久,| 亚洲av电影在线进入| 国产伦人伦偷精品视频| 一个人观看的视频www高清免费观看| 国产成人影院久久av| 给我免费播放毛片高清在线观看| 黄色片一级片一级黄色片| 性欧美人与动物交配| 日韩欧美在线二视频| 久久国产精品人妻蜜桃| 亚洲av中文字字幕乱码综合| 亚洲久久久久久中文字幕| 精品日产1卡2卡| 18美女黄网站色大片免费观看| 国产精华一区二区三区| 国产精品美女特级片免费视频播放器| xxxwww97欧美| 欧美中文综合在线视频| 久久久精品大字幕| 免费av观看视频| 国产精品一区二区三区四区免费观看 | 亚洲国产精品成人综合色| 午夜老司机福利剧场| 欧美日韩福利视频一区二区| av片东京热男人的天堂| 亚洲avbb在线观看| 观看免费一级毛片| 亚洲中文日韩欧美视频| 欧美乱码精品一区二区三区| 99国产精品一区二区三区| 黄色日韩在线| 精品人妻偷拍中文字幕| 日韩成人在线观看一区二区三区| avwww免费| 欧美黑人欧美精品刺激| 国产精品自产拍在线观看55亚洲| 偷拍熟女少妇极品色| 两个人视频免费观看高清| 亚洲av免费高清在线观看| x7x7x7水蜜桃| 欧美另类亚洲清纯唯美| 可以在线观看毛片的网站| 成年版毛片免费区| 国产乱人视频| 国产精品香港三级国产av潘金莲| 国产一区二区在线av高清观看| 超碰av人人做人人爽久久 | 中文字幕高清在线视频| 亚洲国产精品合色在线| 国产极品精品免费视频能看的| 免费看a级黄色片| 久久中文看片网| 一个人免费在线观看电影| 日本一二三区视频观看| 亚洲成a人片在线一区二区| 黄色丝袜av网址大全| 成人永久免费在线观看视频| 女同久久另类99精品国产91| 97碰自拍视频| 丰满人妻一区二区三区视频av | 99久久精品国产亚洲精品| 亚洲天堂国产精品一区在线| 亚洲成av人片在线播放无| 女同久久另类99精品国产91| 村上凉子中文字幕在线| 国产91精品成人一区二区三区| 日本免费a在线| 麻豆国产av国片精品| 亚洲不卡免费看| 91在线观看av| 麻豆成人av在线观看| 此物有八面人人有两片| 黄色女人牲交| 免费看美女性在线毛片视频| 成年人黄色毛片网站| 人妻久久中文字幕网| 国产精品电影一区二区三区| 亚洲av二区三区四区| 波多野结衣高清无吗| 亚洲av免费在线观看| 亚洲av熟女| 两人在一起打扑克的视频| 日韩欧美国产在线观看| 久久精品亚洲精品国产色婷小说| 国产一级毛片七仙女欲春2| 91av网一区二区| 哪里可以看免费的av片| 一夜夜www| 在线免费观看不下载黄p国产 | 长腿黑丝高跟| 日韩欧美精品免费久久 | 欧美3d第一页| 国产成人av激情在线播放| 少妇高潮的动态图| 男人舔奶头视频| 黄片小视频在线播放| 我要搜黄色片| 高清毛片免费观看视频网站| 国产中年淑女户外野战色| 在线播放国产精品三级| 色综合婷婷激情| 国产精品亚洲一级av第二区| 国产三级中文精品| 热99在线观看视频| 国产一区二区亚洲精品在线观看| 老鸭窝网址在线观看| 色视频www国产| 成人一区二区视频在线观看| 黄色成人免费大全| 成人国产一区最新在线观看| 天堂√8在线中文| 欧美zozozo另类| 99在线视频只有这里精品首页| 校园春色视频在线观看| 啦啦啦观看免费观看视频高清| 白带黄色成豆腐渣| 在线观看av片永久免费下载| 久久久久久久久久黄片| 观看美女的网站| 精品熟女少妇八av免费久了| 国产 一区 欧美 日韩| 真人做人爱边吃奶动态| 夜夜看夜夜爽夜夜摸| 欧美在线黄色| 女生性感内裤真人,穿戴方法视频| 午夜福利在线观看吧| 色综合欧美亚洲国产小说| 亚洲午夜理论影院| 一进一出抽搐gif免费好疼| 宅男免费午夜| 久久天躁狠狠躁夜夜2o2o| 少妇的丰满在线观看| 啦啦啦韩国在线观看视频| 国产精品久久久久久久电影 | 一级黄片播放器| 黑人欧美特级aaaaaa片| 又黄又爽又免费观看的视频| av中文乱码字幕在线| 九色成人免费人妻av| 性色avwww在线观看| 欧洲精品卡2卡3卡4卡5卡区| 香蕉久久夜色| 三级国产精品欧美在线观看| 99久久精品热视频| 男女之事视频高清在线观看| 免费在线观看亚洲国产| 一级毛片女人18水好多| 国产亚洲精品久久久久久毛片| 一进一出好大好爽视频| 欧美黑人巨大hd| 欧美一级毛片孕妇| 国产熟女xx| 精品人妻1区二区| 少妇的逼好多水| 久久久久九九精品影院| 亚洲不卡免费看| 岛国在线观看网站| 好男人电影高清在线观看| 欧美一区二区国产精品久久精品| 成人性生交大片免费视频hd| 免费人成在线观看视频色| 日韩中文字幕欧美一区二区| avwww免费| 每晚都被弄得嗷嗷叫到高潮| 18禁国产床啪视频网站| 一区二区三区免费毛片| 蜜桃久久精品国产亚洲av| 少妇高潮的动态图| 欧美+亚洲+日韩+国产| 国产一区二区激情短视频| 国产淫片久久久久久久久 | 亚洲av熟女| 亚洲成人久久爱视频| 国产成人av激情在线播放| 一区二区三区高清视频在线| 亚洲18禁久久av| 99在线视频只有这里精品首页| 亚洲人成网站高清观看| 国产黄片美女视频| 欧美bdsm另类| 欧美激情在线99| 精品99又大又爽又粗少妇毛片 | 舔av片在线| 可以在线观看的亚洲视频| 99热这里只有是精品50| 日韩中文字幕欧美一区二区| 精品乱码久久久久久99久播| 午夜福利在线观看吧| 内地一区二区视频在线| 少妇高潮的动态图| 成人三级黄色视频| 亚洲七黄色美女视频| 久久久久性生活片| 真人做人爱边吃奶动态| 日本黄大片高清| 99精品欧美一区二区三区四区| 男人和女人高潮做爰伦理| 国产在视频线在精品| 国产老妇女一区| 午夜免费观看网址| 最后的刺客免费高清国语| 级片在线观看| 日本黄色视频三级网站网址| 国产91精品成人一区二区三区| 岛国在线免费视频观看| 特大巨黑吊av在线直播| 特级一级黄色大片| 久久人人精品亚洲av| 51午夜福利影视在线观看| 少妇的丰满在线观看| 日韩国内少妇激情av| www日本黄色视频网| 天堂动漫精品| 亚洲激情在线av| 女人高潮潮喷娇喘18禁视频| 日本 av在线| 久久久久亚洲av毛片大全| 国产毛片a区久久久久| 男女那种视频在线观看| 成人无遮挡网站| 亚洲国产色片| 怎么达到女性高潮| 欧美+亚洲+日韩+国产| a级毛片a级免费在线| www.www免费av| 亚洲在线自拍视频| 精品电影一区二区在线| 欧美性猛交╳xxx乱大交人| www国产在线视频色| 国语自产精品视频在线第100页| 国产不卡一卡二| 免费看美女性在线毛片视频| 精品日产1卡2卡| 精品不卡国产一区二区三区| 亚洲av美国av| 国产成人aa在线观看| 少妇的逼水好多| 村上凉子中文字幕在线| 亚洲国产精品成人综合色| 中文在线观看免费www的网站| 欧美区成人在线视频| 亚洲国产欧洲综合997久久,| 午夜精品久久久久久毛片777| 美女被艹到高潮喷水动态| 午夜福利高清视频| 欧美性感艳星| 特大巨黑吊av在线直播| 欧美日韩精品网址| 久久国产乱子伦精品免费另类| 搡老妇女老女人老熟妇| 色综合欧美亚洲国产小说| 欧美乱色亚洲激情| 97人妻精品一区二区三区麻豆| 99国产精品一区二区蜜桃av| 熟女少妇亚洲综合色aaa.| 老熟妇乱子伦视频在线观看| 久久婷婷人人爽人人干人人爱| 嫩草影院入口| 亚洲欧美日韩高清在线视频| 91久久精品电影网| 国产三级黄色录像| 日本黄色片子视频| 熟妇人妻久久中文字幕3abv| 91av网一区二区| 精品久久久久久久末码| 国产精品一区二区免费欧美| 亚洲一区二区三区色噜噜| 久久香蕉国产精品| 在线观看免费午夜福利视频| 免费av观看视频| 美女免费视频网站| 夜夜爽天天搞| 亚洲欧美日韩高清在线视频| 国产午夜福利久久久久久| 国产亚洲欧美在线一区二区| 日本撒尿小便嘘嘘汇集6| 国产精品99久久99久久久不卡| 亚洲国产欧洲综合997久久,| 51国产日韩欧美| 99热精品在线国产| 久久国产乱子伦精品免费另类| 欧美成人性av电影在线观看| 精品午夜福利视频在线观看一区| 国产伦一二天堂av在线观看| 他把我摸到了高潮在线观看| 国产成人aa在线观看| 久久久久国产精品人妻aⅴ院| 午夜免费观看网址| 熟女人妻精品中文字幕| 啦啦啦韩国在线观看视频| 一级a爱片免费观看的视频| 免费看十八禁软件| 在线播放国产精品三级| 国产v大片淫在线免费观看| 全区人妻精品视频| 亚洲av日韩精品久久久久久密| 成熟少妇高潮喷水视频| 小说图片视频综合网站| 亚洲狠狠婷婷综合久久图片| 国产又黄又爽又无遮挡在线| 亚洲真实伦在线观看| 嫩草影院入口| 亚洲av中文字字幕乱码综合| 91麻豆精品激情在线观看国产| 首页视频小说图片口味搜索| 91字幕亚洲| 精品欧美国产一区二区三| 欧美不卡视频在线免费观看| 久久久久久久午夜电影| 九色国产91popny在线| 午夜a级毛片| 国产精品久久久久久亚洲av鲁大| 超碰av人人做人人爽久久 | 国产毛片a区久久久久| 国产精品一及| 99久久综合精品五月天人人| 九九在线视频观看精品| 天天躁日日操中文字幕| 日本熟妇午夜| 久久精品影院6| 九色国产91popny在线| 国产蜜桃级精品一区二区三区| 美女高潮的动态| 久久久久亚洲av毛片大全| 99久久无色码亚洲精品果冻| 麻豆一二三区av精品| 亚洲激情在线av| 色av中文字幕| 欧美日本亚洲视频在线播放| 久99久视频精品免费| 亚洲中文字幕日韩| 在线观看美女被高潮喷水网站 | 欧美中文综合在线视频| 日韩精品青青久久久久久| 国产欧美日韩一区二区精品| 亚洲av二区三区四区| 18禁在线播放成人免费| 国产成人av激情在线播放| 国内揄拍国产精品人妻在线| 少妇人妻精品综合一区二区 | 两人在一起打扑克的视频| 成人av一区二区三区在线看| 丝袜美腿在线中文| 丁香六月欧美| 亚洲精品日韩av片在线观看 | 看黄色毛片网站| 欧美黑人欧美精品刺激| 99久国产av精品| 欧美又色又爽又黄视频| 禁无遮挡网站| 国产免费男女视频| www.999成人在线观看| 狂野欧美激情性xxxx| 嫩草影视91久久| 久久久国产成人免费| 又紧又爽又黄一区二区| 69人妻影院| a级一级毛片免费在线观看| 丰满的人妻完整版| 日本撒尿小便嘘嘘汇集6| 日本黄大片高清| www.www免费av| 在线十欧美十亚洲十日本专区| 国产一区二区激情短视频| 神马国产精品三级电影在线观看| 欧美性感艳星| 国产成人av激情在线播放| 久久精品国产自在天天线| 免费在线观看成人毛片| 国内精品美女久久久久久| 国内毛片毛片毛片毛片毛片| 免费在线观看日本一区| 国产精品av视频在线免费观看| 别揉我奶头~嗯~啊~动态视频| 欧美日韩瑟瑟在线播放| 18禁在线播放成人免费| 天堂影院成人在线观看| 成人三级黄色视频| 两人在一起打扑克的视频| 亚洲中文字幕一区二区三区有码在线看| 69人妻影院| 国产综合懂色| 中文字幕人成人乱码亚洲影| 国产真实乱freesex| 国产老妇女一区| 国产黄色小视频在线观看| 亚洲在线观看片| 亚洲国产精品sss在线观看| 91久久精品国产一区二区成人 | av在线蜜桃| 在线观看美女被高潮喷水网站 | 可以在线观看的亚洲视频| av福利片在线观看| 免费大片18禁| 国产亚洲精品综合一区在线观看| 精品乱码久久久久久99久播| 日韩欧美在线二视频| 欧美在线黄色| 国产成人欧美在线观看| 97碰自拍视频| 国产麻豆成人av免费视频| 一区二区三区激情视频| 色播亚洲综合网| 亚洲一区二区三区不卡视频| АⅤ资源中文在线天堂| 天堂影院成人在线观看| 国产午夜精品久久久久久一区二区三区 | 男人舔女人下体高潮全视频| 国产成人av教育| 欧美一区二区精品小视频在线| 中国美女看黄片| 一区二区三区国产精品乱码| 国产欧美日韩一区二区三| 色视频www国产| 国产一区二区三区视频了| 日日干狠狠操夜夜爽| 看免费av毛片| 宅男免费午夜| 亚洲人成网站在线播| 午夜福利18| 三级毛片av免费| 精品一区二区三区视频在线观看免费| 制服丝袜大香蕉在线| 午夜免费男女啪啪视频观看 | 成年版毛片免费区| 国产精品女同一区二区软件 | 色播亚洲综合网| 内射极品少妇av片p| or卡值多少钱| 亚洲av中文字字幕乱码综合| 国产欧美日韩精品亚洲av| 中文字幕久久专区| www.999成人在线观看| 日本一本二区三区精品| 91麻豆精品激情在线观看国产| 高清在线国产一区| 国产精品久久久人人做人人爽| 日本在线视频免费播放| 久久亚洲精品不卡| 欧美国产日韩亚洲一区| 国产蜜桃级精品一区二区三区| 国产av一区在线观看免费| 黄色女人牲交| 国内久久婷婷六月综合欲色啪| 亚洲美女视频黄频| 亚洲avbb在线观看| 搡女人真爽免费视频火全软件 | 久久草成人影院| 久久久久亚洲av毛片大全| 中文字幕精品亚洲无线码一区| 欧美+亚洲+日韩+国产| 三级毛片av免费| 亚洲精品亚洲一区二区| 三级毛片av免费| 欧美+亚洲+日韩+国产| 操出白浆在线播放| 又粗又爽又猛毛片免费看| 中文字幕精品亚洲无线码一区| 久久精品人妻少妇| 国产成人av教育| 波多野结衣高清无吗| 亚洲成a人片在线一区二区| 亚洲自拍偷在线| 91麻豆av在线| 97人妻精品一区二区三区麻豆| 熟女少妇亚洲综合色aaa.| 一a级毛片在线观看| 亚洲天堂国产精品一区在线| 国产91精品成人一区二区三区| 此物有八面人人有两片| 夜夜夜夜夜久久久久| 又粗又爽又猛毛片免费看| 亚洲成人久久爱视频| 老熟妇仑乱视频hdxx| bbb黄色大片| 日韩免费av在线播放| 非洲黑人性xxxx精品又粗又长| 一进一出抽搐gif免费好疼| 日韩有码中文字幕| 一区二区三区高清视频在线| 一本一本综合久久| 国产高潮美女av| 久99久视频精品免费| 亚洲 国产 在线| 女警被强在线播放| 欧美+亚洲+日韩+国产| 欧美黄色淫秽网站| 久久国产乱子伦精品免费另类| 欧美中文日本在线观看视频| 免费看a级黄色片| 国产视频内射| 国产成人系列免费观看| 狂野欧美激情性xxxx| 最新中文字幕久久久久| 丰满人妻熟妇乱又伦精品不卡| 国产伦精品一区二区三区视频9 | 中文在线观看免费www的网站| 亚洲精品日韩av片在线观看 | 国产亚洲精品综合一区在线观看| 成人特级av手机在线观看| 真人做人爱边吃奶动态| 一个人免费在线观看的高清视频| 欧美日本视频| 一个人免费在线观看电影| 亚洲av二区三区四区| 精品久久久久久久末码| 欧美三级亚洲精品| 最近在线观看免费完整版| 12—13女人毛片做爰片一| 18禁黄网站禁片午夜丰满| 国内少妇人妻偷人精品xxx网站| 熟女人妻精品中文字幕| 亚洲自拍偷在线| 精品99又大又爽又粗少妇毛片 | 免费观看精品视频网站| 国产精品三级大全| 又粗又爽又猛毛片免费看| 久久久久性生活片| 日韩成人在线观看一区二区三区| 婷婷丁香在线五月| 国产精品久久久人人做人人爽| 亚洲人成网站高清观看| 变态另类丝袜制服| 亚洲一区二区三区不卡视频|