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

    基于最小均方自適應(yīng)算法的光纖陀螺信號實時濾波方法

    2020-04-06 08:27:14杜建邦何金陽
    中國慣性技術(shù)學(xué)報 2020年6期
    關(guān)鍵詞:均方階數(shù)陀螺

    杜建邦,何金陽,2,卓 超

    (1. 北京航天自動控制研究所,北京 100854;2. 宇航智能控制技術(shù)國家級重點實驗室,北京 100854)

    高速、高機動飛行器具有突出的快時變、強耦合、非線性和不確定性等特點,這給其上安裝的旋轉(zhuǎn)調(diào)制式慣導(dǎo)系統(tǒng)的旋轉(zhuǎn)軸穩(wěn)定控制帶來了嚴峻的挑戰(zhàn)。光纖陀螺具有工作原理先進、結(jié)構(gòu)簡單、無活動部件及抗沖擊能力強等優(yōu)點,是高精度穩(wěn)定控制系統(tǒng)的理想角運動敏感元件。但光纖陀螺相比于傳統(tǒng)的機電陀螺,其輸出信號易受光路、電路隨機噪聲的干擾,產(chǎn)生的角速率漂移會嚴重影響系統(tǒng)控制精度[1]。為了實現(xiàn)系統(tǒng)的高精度穩(wěn)定控制,需要對控制系統(tǒng)中的反饋信號—光纖陀螺的輸出進行降噪濾波處理[2,3]。

    目前,工程中廣泛采用低通濾波器對陀螺信號進行濾波,濾波器的截止頻率、通帶等參數(shù)主要根據(jù)工程要求及經(jīng)驗確定。由于傳統(tǒng)低通濾波器采用固定的濾波器系數(shù),難以滿足光纖陀螺非平穩(wěn)隨機信號條件下的精確濾波。自適應(yīng)濾波器則能夠?qū)崟r地調(diào)節(jié)當前時刻的濾波系數(shù)以適應(yīng)隨機信號的時變統(tǒng)計特性,可以實現(xiàn)一定意義下的最優(yōu)濾波,因此更加適用于光纖陀螺輸出信號的濾波處理。在自適應(yīng)濾波器設(shè)計中,最小均方(Least Mean Square,LMS)算法使用隨機梯度下降的方法實現(xiàn)代價函數(shù)的最小化,具有計算復(fù)雜度低、無需統(tǒng)計數(shù)據(jù)的先驗知識和均值無偏地收斂到維納解等優(yōu)點,成為自適應(yīng)算法中應(yīng)用最廣泛的一種[4]。為解決各類實際問題,相關(guān)學(xué)者提出了諸多變步長LMS 算法及變換域的LMS 算法[5,6]。文獻[7]提出了級聯(lián)的LMS 自適應(yīng)的降噪方法。該算法采用了變步長和不同濾波器初始權(quán)值的自適應(yīng)算法,通過仿真結(jié)果表明,該方法在濾波精度和實時性等綜合性能方面具有一定的優(yōu)越性。文獻[8]針對光纖陀螺輸出噪聲的問題,提出了一種基于最小均方法與二代小波變換相結(jié)合的降噪方法,該方法提高了捷聯(lián)慣性導(dǎo)航初始對準姿態(tài)精度。文獻[9]提出了基于歸一化LMS 算法的前向線性預(yù)測濾波器的降噪方法,該方法顯著降低了角隨機游走噪聲,實現(xiàn)了單個光纖陀螺對地球自轉(zhuǎn)速度的精確測量。然而,目前大多LMS 自適應(yīng)濾波器的相關(guān)研究均以靜態(tài)濾波精度作為主要設(shè)計指標,沒有考慮時延特性,適用于陀螺尋北和初始對準等動態(tài)性要求較低的場景,無法滿足高動態(tài)條件下旋轉(zhuǎn)調(diào)制式慣導(dǎo)系統(tǒng)旋轉(zhuǎn)軸穩(wěn)定控制系統(tǒng)的濾波需求。

    本文針對旋轉(zhuǎn)慣導(dǎo)穩(wěn)定控制對濾波器高頻消噪與時延等要求,在綜合考慮了期望信號、濾波器階數(shù)和信號采樣率等影響因素后,設(shè)計了一種優(yōu)化的LMS自適應(yīng)濾波器。該濾波器兼顧了動靜態(tài)性能要求,適用于大動態(tài)、高精度穩(wěn)定控制系統(tǒng),具有結(jié)構(gòu)簡單、工程實現(xiàn)方便等優(yōu)點。

    1 LMS 自適應(yīng)濾波器原理

    最小均方(LMS)自適應(yīng)濾波器的結(jié)構(gòu)如圖1(a)所示,包括橫向濾波器與LMS 算法兩部分[9]。橫向濾波器(圖1(b)所示)本質(zhì)為FIR 結(jié)構(gòu)的維納濾波器,其輸出信號是輸入信號u(n- 1 ),u(n- 2)…u(n-M)的線性組合。LMS 算法根據(jù)期望信號d與的誤差調(diào)整橫向濾波器的權(quán)系數(shù)向量,從而適應(yīng)隨機信號的時變統(tǒng)計特性。

    圖1 LMS 自適應(yīng)濾波器原理Fig.1 Principle of LMS adaptive filter

    LMS 自適應(yīng)濾波器的輸出是當前采樣值u(n)的估計值,的計算公式為:

    式(1)中,BT= [b1,b2,b3,…bm]是橫向濾波器系數(shù)組成的向量;U(n- 1) = [u(n- 1 ),u(n- 2),u(n- 3),…u(n-M)]T是濾波器輸入向量。

    定義估計誤差fM為:

    將估計誤差的最小均方誤差作為最優(yōu)的統(tǒng)計準則,定義代價函數(shù)J為:

    式(3)中,E表示估計誤差平方的數(shù)學(xué)期望;J隨時間的變化表征了濾波器的收斂速度與收斂精度。

    需要解決的問題是如何選擇濾波器的權(quán)值使得的均方誤差最小,或者說當J最小時濾波器的權(quán)值最優(yōu)。最優(yōu)權(quán)向量Bopt可表示為:

    式(4)中,R為輸入信號的自相關(guān)矩陣,rxd為輸入信號與期望響應(yīng)的互相關(guān)矢量。在實際系統(tǒng)中,特別是環(huán)境未知時,一般很難獲得濾波器輸入信號與期望響應(yīng)的統(tǒng)計特性(R和rxd),因此,采用對非平穩(wěn)環(huán)境具有學(xué)習(xí)能力的自適應(yīng)LMS 算法對權(quán)系數(shù)向量進行實時估計,逐步逼近最優(yōu)權(quán)值。根據(jù)LMS 算法,通過以下公式反復(fù)迭代,權(quán)值逐漸逼近其最優(yōu)解:

    式(5)中,μ是權(quán)值更新步長,μ的收斂范圍為:

    λmax為輸入信號自相關(guān)矩陣R的最大特征值。

    相關(guān)學(xué)者在傳統(tǒng)LMS 算法上提出了許多優(yōu)化算法,其中歸一化LMS 算法采用了可變更新步長,具有收斂速度更快、濾波器時延小的優(yōu)點,因此適用于穩(wěn)定控制系統(tǒng)的濾波需求,其權(quán)值迭代公式為:

    式(7)中,γ是一個很小的常數(shù)以保證分母不為零,μ的收斂范圍為:

    2 LMS 自適應(yīng)濾波器優(yōu)化設(shè)計

    針對高精度穩(wěn)定控制系統(tǒng)對濾波器的性能要求,為了實現(xiàn)LMS 自適應(yīng)濾波器的優(yōu)化設(shè)計,需綜合考慮各種因素對濾波器性能的影響。根據(jù)第一節(jié),歸一化LMS 自適應(yīng)濾波器的性能主要取決于權(quán)值更新步長μ、未知環(huán)境下的期望信號d和濾波器階數(shù)M三個參數(shù)。本節(jié)結(jié)合實際控制系統(tǒng)對濾波器的性能要求以及光纖陀螺實驗數(shù)據(jù),分析相關(guān)參數(shù)對濾波性能的影響并進行濾波器的優(yōu)化設(shè)計。

    2.1 權(quán)值更新步長μ

    由式(7)可得,權(quán)值更新步長μ決定了權(quán)值大小。為了進一步分析μ對濾波性能的影響,以光纖陀螺靜態(tài)數(shù)據(jù)作為輸入信號,其他濾波器參數(shù)一定,對比μ=0.001,μ=0.1 和μ=1(0<μ≤1)三種不同步長濾波器的濾波效果。圖2 展示了濾波結(jié)果的均方誤差J隨時間變化的曲線。

    圖2 采用不同步長濾波的均方誤差曲線Fig.2 MSE curves with different update step length

    由圖2 可得,步長μ影響著收斂速度和濾波精度。隨著μ的增大,收斂速度變快,但收斂后在最優(yōu)解附近波動幅值變大,穩(wěn)態(tài)精度相應(yīng)降低;μ過小,雖可保證濾波精度,但收斂時間大大增加,無法保證系統(tǒng)實時性。因此需要根據(jù)收斂速度和濾波精度折中選取μ。旋轉(zhuǎn)慣導(dǎo)穩(wěn)定控制系統(tǒng)要求濾波器在動態(tài)突變環(huán)境下能夠快速收斂,相比于收斂后較低量級的波動幅值,μ的選取應(yīng)側(cè)重考慮收斂速度。

    2.2 期望信號d

    目前已發(fā)表的相關(guān)文獻中,自適應(yīng)濾波器的期望信號通常采用當前時刻的信號采樣值,即d(n) =u(n),但是當前時刻采樣值u(n)包含了較大的噪聲,導(dǎo)致期望信號的波動,影響了濾波精度。為了減小期望信號這種不穩(wěn)定性的影響,對期望信號的選取進行優(yōu)化,選取過去時刻采樣值的平均值作為期望信號,即

    式(9)中Δ表示選取的采樣點數(shù)。Δ選取的原則為既要有效地減小當前時刻采樣值的噪聲,也不能對當前時刻采樣值的真值造成影響以及產(chǎn)生太大的時延。隨著Δ的增大,當前時刻信號的噪聲強度隨之降低,但增加了濾波器輸出信號的時延。另外,隨著Δ的增大,與當前時刻信號的相關(guān)性逐漸降低,也會影響當前時刻期望信號真值的確定。

    通過以上分析,利用光纖陀螺靜態(tài)數(shù)據(jù)作為輸入信號,其他濾波器參數(shù)一定,選取Δ=0,Δ=5 和Δ=M(M為濾波器階數(shù),此處選定M=10)三種典型的期望值進行濾波器性能的比對分析。不同期望信號濾波后的結(jié)果及其均值、標準差統(tǒng)計如圖3 和表1 所示。

    不同期望信號濾波后的結(jié)果及其均值、標準差統(tǒng)計如圖3 和表1 所示。

    圖3 采用不同期望信號的濾波結(jié)果Fig.3 Filtering results with different expected signal

    表1 不同期望信號下濾波后信號濾波均值與標準差統(tǒng)計Tab.1 Mean value and standard deviation of filtering results with different expected signal

    由圖3 和表1 可得,當Δ由0 增大至M時,濾波器降噪能力明顯提高。

    進一步驗證Δ增大對動態(tài)響應(yīng)的影響,令濾波器跟蹤正弦信號,濾波后產(chǎn)生的相位延遲可用式(10)估算:

    其中,N為濾波器輸出信號滯后的采樣點個數(shù),ω為信號的角頻率,fs為信號采樣頻率。對于角頻率為20 Hz 的正弦信號濾波后相位延遲如表2 所示。

    表2 不同Δ 的相位延遲Tab.2 Phase delay with different Δ

    由表2 可知,隨著Δ增大,期望信號包含了更多過去時刻的采樣信息,濾波器的動態(tài)響應(yīng)相應(yīng)降低。因此,Δ影響濾波精度和時延,需綜合考慮信號的噪聲特性、信號的相關(guān)性以及控制帶寬的要求選取Δ值。

    2.3 濾波器階數(shù)M

    濾波器階數(shù)M影響著濾波器收斂速度和濾波精度。為了探究M對濾波性能的影響,利用光纖陀螺靜態(tài)數(shù)據(jù)作為輸入信號,其他濾波器參數(shù)一定,對比M=10,M=50 和M=100 三種階數(shù)濾波器的濾波效果。濾波結(jié)果信號的均值、標準差統(tǒng)計如表3 所示,濾波結(jié)果的均方誤差J隨時間變化的曲線如圖4 所示。

    表3 采用不同階數(shù)濾波后信號的均值與標準差統(tǒng)計Tab.3 Mean value and standard deviation of filtering results with different filter orders

    圖4 采用不同階數(shù)濾波的均方誤差曲線Fig.4 MSE curves with different filter orders

    由圖4 和表3 可得,濾波器階數(shù)M增大,濾波精度有所提高,但濾波器收斂速度明顯降低。針對高精度穩(wěn)定控制系統(tǒng)對濾波器短時延的需求,濾波器階數(shù)M的選取也必須優(yōu)先考慮收斂速度。

    2.4 采樣率fs 的影響

    目前的相關(guān)文獻大多針對某一特定采樣率進行分析,而采樣率對濾波器的設(shè)計非常重要。高頻采樣得到的濾波器輸入信號通常噪聲較強,通過對高頻采樣數(shù)據(jù)進行適當?shù)仄交幚恚梢杂行Ы档透哳l噪聲,但會引起濾波器時延的增加。在設(shè)計濾波器時要綜合考慮采樣頻率的影響,合理選取濾波器參數(shù)。

    3 實驗驗證與分析

    根據(jù)2.1~2.4 節(jié)分析可得,LMS 自適應(yīng)濾波器的濾波精度主要取決于期望信號d表達式中的Δ值,收斂速度主要取決于權(quán)值更新步長μ和濾波器階數(shù)M。在綜合考慮濾波精度和收斂速度的需求、光纖陀螺輸出信號特性與采樣率后,優(yōu)化選取的LMS 自適應(yīng)濾波器參數(shù)為:

    為了驗證優(yōu)化后LMS 自適應(yīng)濾波器的濾波效果,設(shè)計了靜態(tài)與動態(tài)兩類實驗并與傳統(tǒng)低通濾波器進行性能對比。采用具有線性相位的10 階FIR 濾波器作為對比,采樣率為1000 Hz,濾波器系數(shù)見表4,F(xiàn)IR 濾波器截止頻率為250 Hz,通帶衰減小于10-2dB,遵循最小相角設(shè)計,20 Hz 處的相角為-0.1988 rad。

    表4 FIR 濾波器系數(shù)Tab.4 Coefficients of FIR filter

    3.1 靜態(tài)實驗

    將光纖陀螺靜置于轉(zhuǎn)臺,采集陀螺輸出數(shù)據(jù)并利用LMS 與FIR 濾波器進行實時濾波處理,圖5 展示了兩種濾波器的濾波效果,表5 統(tǒng)計了濾波后信號的均值與標準差。

    由表5 和圖5 可見,優(yōu)化的LMS 自適應(yīng)濾波器有效地降低了高頻噪聲,高頻噪聲降低了約42.35%,降噪能力比低通FIR 濾波器提高了約9.16%。

    圖5 靜態(tài)條件下LMS 和FIR 濾波器的濾波結(jié)果對比Fig.5 Comparison of filtering results of LMS and FIR filters in static experiment

    表5 經(jīng)兩種濾波器濾波后信號的均值與標準差統(tǒng)計Tab.5 Mean value and standard deviation of the output signal processed by two filters

    3.2 動態(tài)實驗

    將光纖陀螺置于轉(zhuǎn)臺進行搖擺試驗,受轉(zhuǎn)臺性能所限,只控制轉(zhuǎn)臺做1 Hz 搖擺運動,采集陀螺輸出數(shù)據(jù)并做實時濾波處理,兩種濾波器的動態(tài)跟蹤情況如圖6 所示。

    由圖6 可知,在1 Hz 搖擺運動條件下,優(yōu)化的LMS 自適應(yīng)濾波器比FIR 濾波器的相位延遲小0.108 °,提高了約18.96%。

    圖6 搖擺狀態(tài)下LMS 和FIR 濾波器的濾波效果對比Fig.6 Filtering effects of LMS and FIR filters in swaying experiment

    3.2 動態(tài)突變實驗

    仿真產(chǎn)生一組動態(tài)突變信號,選取不連續(xù)角頻率的正弦運動,角頻率由1 Hz→20 Hz→40 Hz→100 Hz突變,兩種濾波器對動態(tài)突變信號的跟蹤情況如圖7(a)所示,LMS 自適應(yīng)濾波器收斂過程的均方誤差曲線如圖7(b)所示。由圖7 可得,LMS 自適應(yīng)濾波器可快速收斂,能夠較好地跟蹤動態(tài)突變信號,且相位延遲小于低通FIR 濾波器。

    圖7 動態(tài)突變情況下LMS 和FIR 濾波器的濾波結(jié)果對比Fig.7 Comparison of filtering results of LMS and FIR filters in sudden-change dynamic experiment

    4 結(jié) 論

    本文提出了一種適用于旋轉(zhuǎn)調(diào)制式慣導(dǎo)的旋轉(zhuǎn)軸穩(wěn)定控制的優(yōu)化LMS 自適應(yīng)濾波方法。通過對濾波器主要參數(shù)進行優(yōu)化選取,濾波性能可以較好地滿足高頻消噪與時延的要求。相比于工程中常用的傳統(tǒng)FIR 濾波器,其降噪能力提高了約9.16%,相位延遲降低了約18.96%。該濾波器結(jié)構(gòu)簡單、易于實現(xiàn),具有工程應(yīng)用參考價值。后續(xù)將根據(jù)飛行試驗結(jié)果,進一步完善算法設(shè)計。

    猜你喜歡
    均方階數(shù)陀螺
    一類隨機積分微分方程的均方漸近概周期解
    關(guān)于無窮小階數(shù)的幾點注記
    確定有限級數(shù)解的階數(shù)上界的一種n階展開方法
    Beidou, le système de navigation par satellite compatible et interopérable
    做個紙陀螺
    玩陀螺
    陀螺轉(zhuǎn)轉(zhuǎn)轉(zhuǎn)
    軍事文摘(2018年24期)2018-12-26 00:58:18
    我最喜歡的陀螺
    快樂語文(2018年36期)2018-03-12 00:56:02
    基于抗差最小均方估計的輸電線路參數(shù)辨識
    基于隨機牽制控制的復(fù)雜網(wǎng)絡(luò)均方簇同步
    一区二区三区高清视频在线| xxx96com| 可以在线观看毛片的网站| 精品国内亚洲2022精品成人| 91九色精品人成在线观看| 久久久久国产精品人妻aⅴ院| 老司机在亚洲福利影院| 国产高清激情床上av| 黄片小视频在线播放| 一个人观看的视频www高清免费观看 | 亚洲狠狠婷婷综合久久图片| 精品欧美国产一区二区三| 国产成人av教育| 夜夜躁狠狠躁天天躁| 人人妻人人看人人澡| 国产精品综合久久久久久久免费| 欧美一区二区国产精品久久精品 | 亚洲av五月六月丁香网| 国产精品av视频在线免费观看| 日本免费a在线| 国产成人影院久久av| 1024手机看黄色片| 亚洲人成网站高清观看| 一区二区三区高清视频在线| 淫秽高清视频在线观看| 97超级碰碰碰精品色视频在线观看| 国产激情偷乱视频一区二区| 日本五十路高清| 国产一区二区激情短视频| 好男人在线观看高清免费视频| 亚洲第一电影网av| 国产视频一区二区在线看| 两个人免费观看高清视频| 亚洲七黄色美女视频| 又爽又黄无遮挡网站| 欧美乱码精品一区二区三区| 日韩成人在线观看一区二区三区| 人成视频在线观看免费观看| 99久久综合精品五月天人人| 熟妇人妻久久中文字幕3abv| 人人妻人人澡欧美一区二区| а√天堂www在线а√下载| 国产精品亚洲av一区麻豆| 搡老熟女国产l中国老女人| 日韩欧美国产一区二区入口| 日韩欧美免费精品| 欧美中文日本在线观看视频| 国产黄片美女视频| 欧美黑人欧美精品刺激| 成人av在线播放网站| 亚洲欧美日韩高清在线视频| 天天躁夜夜躁狠狠躁躁| 91av网站免费观看| 亚洲国产中文字幕在线视频| 亚洲国产高清在线一区二区三| 国产精品一及| 欧美日韩一级在线毛片| 久久精品国产综合久久久| 欧美又色又爽又黄视频| 久久久精品欧美日韩精品| 国产麻豆成人av免费视频| 国产激情欧美一区二区| 人人妻,人人澡人人爽秒播| 一个人免费在线观看的高清视频| 欧美中文日本在线观看视频| 亚洲国产精品合色在线| 丰满人妻一区二区三区视频av | 一进一出好大好爽视频| 1024手机看黄色片| 欧美乱码精品一区二区三区| 亚洲精品久久成人aⅴ小说| 69av精品久久久久久| 色av中文字幕| 三级男女做爰猛烈吃奶摸视频| 久久精品aⅴ一区二区三区四区| 人成视频在线观看免费观看| www.熟女人妻精品国产| 我的老师免费观看完整版| www.999成人在线观看| 成年人黄色毛片网站| 精品一区二区三区av网在线观看| 成人国产一区最新在线观看| 男女之事视频高清在线观看| www.熟女人妻精品国产| 亚洲专区中文字幕在线| 18禁美女被吸乳视频| 婷婷丁香在线五月| 在线观看66精品国产| 欧美午夜高清在线| 在线观看日韩欧美| 在线国产一区二区在线| 男女下面进入的视频免费午夜| 一个人观看的视频www高清免费观看 | 久久久国产成人精品二区| 日韩国内少妇激情av| 美女大奶头视频| 久久久久久久精品吃奶| 成人亚洲精品av一区二区| 国产乱人伦免费视频| 国产精品久久久久久精品电影| a级毛片在线看网站| 欧美大码av| 国产精品亚洲一级av第二区| 欧美日韩中文字幕国产精品一区二区三区| 特大巨黑吊av在线直播| 午夜两性在线视频| 色综合站精品国产| 亚洲男人的天堂狠狠| 熟女电影av网| 九色成人免费人妻av| 午夜激情福利司机影院| 又黄又爽又免费观看的视频| 99在线人妻在线中文字幕| 日本一区二区免费在线视频| 久久久久久大精品| www国产在线视频色| 亚洲午夜精品一区,二区,三区| 两个人免费观看高清视频| 国产精华一区二区三区| 好看av亚洲va欧美ⅴa在| 亚洲成人中文字幕在线播放| 91国产中文字幕| 人妻丰满熟妇av一区二区三区| 欧美一区二区精品小视频在线| 舔av片在线| 国产单亲对白刺激| 亚洲电影在线观看av| 19禁男女啪啪无遮挡网站| 99热6这里只有精品| 国产伦一二天堂av在线观看| 国产私拍福利视频在线观看| 久久这里只有精品19| 中文资源天堂在线| 91麻豆精品激情在线观看国产| 精品欧美国产一区二区三| 久久久久久久久中文| 97人妻精品一区二区三区麻豆| 一级片免费观看大全| 久久午夜综合久久蜜桃| 久久欧美精品欧美久久欧美| 国产精品 国内视频| 观看免费一级毛片| 村上凉子中文字幕在线| 久久久久九九精品影院| 好看av亚洲va欧美ⅴa在| 久久久国产欧美日韩av| 午夜福利视频1000在线观看| 中文字幕av在线有码专区| 男女午夜视频在线观看| 这个男人来自地球电影免费观看| 久久精品国产综合久久久| 日韩精品免费视频一区二区三区| 成人av在线播放网站| 国产成+人综合+亚洲专区| 久久精品夜夜夜夜夜久久蜜豆 | 国产乱人伦免费视频| 天堂动漫精品| 成熟少妇高潮喷水视频| 午夜福利高清视频| 亚洲黑人精品在线| 久久 成人 亚洲| 国产野战对白在线观看| 精品欧美国产一区二区三| 黑人操中国人逼视频| 国产精品香港三级国产av潘金莲| 一进一出抽搐动态| av福利片在线观看| 久久久精品欧美日韩精品| 法律面前人人平等表现在哪些方面| 精品国产亚洲在线| 欧美成人性av电影在线观看| 成人欧美大片| 欧美午夜高清在线| 无人区码免费观看不卡| av福利片在线观看| 久久中文字幕一级| 国产探花在线观看一区二区| 一区福利在线观看| 99精品久久久久人妻精品| 国产伦在线观看视频一区| 少妇被粗大的猛进出69影院| 国内精品久久久久久久电影| 日本黄色视频三级网站网址| а√天堂www在线а√下载| 欧美久久黑人一区二区| 亚洲成人中文字幕在线播放| 一进一出抽搐动态| 亚洲国产高清在线一区二区三| 国产一区二区三区在线臀色熟女| 中文在线观看免费www的网站 | 久久草成人影院| 黑人巨大精品欧美一区二区mp4| 国产一区二区激情短视频| 欧美一级毛片孕妇| 亚洲全国av大片| 免费搜索国产男女视频| 搡老熟女国产l中国老女人| 在线永久观看黄色视频| 最新在线观看一区二区三区| 麻豆久久精品国产亚洲av| 51午夜福利影视在线观看| 男女那种视频在线观看| 1024手机看黄色片| 麻豆一二三区av精品| 国产一区二区在线观看日韩 | 亚洲av日韩精品久久久久久密| 亚洲avbb在线观看| av天堂在线播放| 1024视频免费在线观看| 超碰成人久久| 国产黄a三级三级三级人| 欧美另类亚洲清纯唯美| 亚洲乱码一区二区免费版| 免费看十八禁软件| 特级一级黄色大片| 午夜福利免费观看在线| 色老头精品视频在线观看| 99在线人妻在线中文字幕| 欧美成人免费av一区二区三区| 妹子高潮喷水视频| 色精品久久人妻99蜜桃| 在线国产一区二区在线| 亚洲精品美女久久av网站| 成年女人毛片免费观看观看9| 夜夜看夜夜爽夜夜摸| 少妇裸体淫交视频免费看高清 | 欧美黑人精品巨大| 久久草成人影院| 两个人的视频大全免费| 99久久无色码亚洲精品果冻| 久久久久免费精品人妻一区二区| 在线免费观看的www视频| 久久草成人影院| 亚洲成人久久性| 99久久无色码亚洲精品果冻| 欧美日韩国产亚洲二区| 久久久久久久午夜电影| 国产在线精品亚洲第一网站| 国产精品久久久av美女十八| 久久久久性生活片| 亚洲专区中文字幕在线| 精品欧美国产一区二区三| 成人国产综合亚洲| 午夜福利欧美成人| 国模一区二区三区四区视频 | av福利片在线观看| 色在线成人网| 最近最新中文字幕大全免费视频| 国产久久久一区二区三区| 欧美日韩黄片免| www.999成人在线观看| 久久久久国内视频| 首页视频小说图片口味搜索| 一本大道久久a久久精品| www.自偷自拍.com| 日韩av在线大香蕉| 亚洲av中文字字幕乱码综合| 美女高潮喷水抽搐中文字幕| 精品日产1卡2卡| 日韩精品青青久久久久久| 高清毛片免费观看视频网站| 男人舔女人的私密视频| 亚洲无线在线观看| 脱女人内裤的视频| 五月玫瑰六月丁香| 最近最新中文字幕大全电影3| 国内少妇人妻偷人精品xxx网站 | 91大片在线观看| 久久午夜综合久久蜜桃| 99国产精品一区二区三区| 日本免费一区二区三区高清不卡| 国产私拍福利视频在线观看| 99riav亚洲国产免费| 成人精品一区二区免费| 我的老师免费观看完整版| 亚洲精品粉嫩美女一区| 久久久久久亚洲精品国产蜜桃av| 9191精品国产免费久久| 精品电影一区二区在线| 听说在线观看完整版免费高清| 国产熟女午夜一区二区三区| 一区二区三区国产精品乱码| 99国产极品粉嫩在线观看| 叶爱在线成人免费视频播放| 少妇的丰满在线观看| 精品福利观看| 天堂动漫精品| 亚洲成a人片在线一区二区| 国产99白浆流出| 亚洲国产欧美一区二区综合| www日本在线高清视频| videosex国产| 男男h啪啪无遮挡| 十八禁网站免费在线| 国产成人啪精品午夜网站| 日韩欧美精品v在线| 国产麻豆成人av免费视频| 国产午夜福利久久久久久| 人成视频在线观看免费观看| 欧美激情久久久久久爽电影| 午夜亚洲福利在线播放| 亚洲一区二区三区不卡视频| 国产黄a三级三级三级人| 久久国产精品人妻蜜桃| 免费av毛片视频| 亚洲av电影不卡..在线观看| 欧美在线黄色| 精品欧美一区二区三区在线| 久久久国产欧美日韩av| 欧美中文日本在线观看视频| 在线视频色国产色| 老司机午夜福利在线观看视频| 午夜福利在线在线| 免费在线观看成人毛片| 亚洲真实伦在线观看| 国产熟女午夜一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 一级作爱视频免费观看| 国内少妇人妻偷人精品xxx网站 | 欧美久久黑人一区二区| 1024香蕉在线观看| 久久午夜综合久久蜜桃| 久久天堂一区二区三区四区| 国内久久婷婷六月综合欲色啪| 国产一区在线观看成人免费| 亚洲欧美精品综合久久99| 欧美又色又爽又黄视频| 国产视频内射| 精品久久久久久,| 老司机午夜十八禁免费视频| 女人高潮潮喷娇喘18禁视频| 久久久久性生活片| 国产亚洲精品第一综合不卡| 成人国语在线视频| 搡老岳熟女国产| 最新在线观看一区二区三区| 欧美av亚洲av综合av国产av| 久久精品国产99精品国产亚洲性色| 亚洲精品在线观看二区| 亚洲免费av在线视频| 日韩欧美国产一区二区入口| 国产成人av教育| 午夜福利高清视频| 非洲黑人性xxxx精品又粗又长| 免费在线观看影片大全网站| 国产真实乱freesex| 久久久久久久久免费视频了| 亚洲男人天堂网一区| 午夜久久久久精精品| 久久国产乱子伦精品免费另类| 一本久久中文字幕| 免费在线观看日本一区| 久久精品91蜜桃| 国产免费av片在线观看野外av| 在线观看免费日韩欧美大片| 日日爽夜夜爽网站| 亚洲成人中文字幕在线播放| 久久午夜亚洲精品久久| 欧美性长视频在线观看| 国产精品一区二区三区四区久久| 很黄的视频免费| 一级黄色大片毛片| 床上黄色一级片| 日本三级黄在线观看| 欧美+亚洲+日韩+国产| 午夜福利在线观看吧| 久久久久久久久中文| 91麻豆精品激情在线观看国产| 性欧美人与动物交配| 两个人的视频大全免费| 欧美黑人巨大hd| netflix在线观看网站| 免费电影在线观看免费观看| 久久精品国产99精品国产亚洲性色| 99久久久亚洲精品蜜臀av| 国产精品久久久久久人妻精品电影| netflix在线观看网站| 嫩草影院精品99| 韩国av一区二区三区四区| 亚洲人成77777在线视频| 国产1区2区3区精品| 国产av不卡久久| 精品久久久久久久久久免费视频| 亚洲熟妇熟女久久| 99国产精品99久久久久| 极品教师在线免费播放| 欧美日韩瑟瑟在线播放| 国模一区二区三区四区视频 | 久久亚洲精品不卡| 91大片在线观看| 美女高潮喷水抽搐中文字幕| 国产精品野战在线观看| 叶爱在线成人免费视频播放| 久久久久国产精品人妻aⅴ院| 99久久无色码亚洲精品果冻| 夜夜看夜夜爽夜夜摸| 久久午夜亚洲精品久久| 日韩精品青青久久久久久| 97超级碰碰碰精品色视频在线观看| 人人妻人人澡欧美一区二区| 免费看日本二区| 男人舔女人下体高潮全视频| 国产aⅴ精品一区二区三区波| 精品国产美女av久久久久小说| 国产又黄又爽又无遮挡在线| 久久久久久久久久黄片| 亚洲精品av麻豆狂野| 女生性感内裤真人,穿戴方法视频| 中文字幕av在线有码专区| av在线天堂中文字幕| 日韩国内少妇激情av| 免费搜索国产男女视频| 久久精品91蜜桃| 一级作爱视频免费观看| 男女午夜视频在线观看| 美女扒开内裤让男人捅视频| 国产精品,欧美在线| 夜夜躁狠狠躁天天躁| 国产一区二区激情短视频| 亚洲午夜理论影院| 99国产精品一区二区蜜桃av| 又紧又爽又黄一区二区| 一本一本综合久久| 99国产极品粉嫩在线观看| 精品高清国产在线一区| 三级毛片av免费| 全区人妻精品视频| cao死你这个sao货| 国产亚洲精品av在线| 日韩三级视频一区二区三区| 国产人伦9x9x在线观看| 亚洲自拍偷在线| 女警被强在线播放| 999久久久精品免费观看国产| 欧美zozozo另类| 女人爽到高潮嗷嗷叫在线视频| 变态另类丝袜制服| 黑人欧美特级aaaaaa片| 国产精品久久久人人做人人爽| 法律面前人人平等表现在哪些方面| 国产精品久久久av美女十八| aaaaa片日本免费| 男女午夜视频在线观看| 亚洲国产精品久久男人天堂| 十八禁网站免费在线| 在线国产一区二区在线| 69av精品久久久久久| 国产亚洲精品久久久久久毛片| 午夜精品久久久久久毛片777| 97碰自拍视频| 午夜福利成人在线免费观看| 深夜精品福利| 午夜日韩欧美国产| 啪啪无遮挡十八禁网站| 色哟哟哟哟哟哟| 亚洲欧美激情综合另类| 在线观看一区二区三区| 亚洲av成人精品一区久久| 欧美性猛交黑人性爽| 亚洲va日本ⅴa欧美va伊人久久| 18美女黄网站色大片免费观看| 三级国产精品欧美在线观看 | √禁漫天堂资源中文www| 深夜精品福利| 91九色精品人成在线观看| 亚洲黑人精品在线| 两个人的视频大全免费| 国产成人欧美在线观看| 91麻豆精品激情在线观看国产| 国产激情欧美一区二区| а√天堂www在线а√下载| 在线观看免费视频日本深夜| 三级毛片av免费| 中文字幕人成人乱码亚洲影| xxx96com| 俄罗斯特黄特色一大片| 在线a可以看的网站| av福利片在线| 中文字幕人成人乱码亚洲影| 国产成人av激情在线播放| 舔av片在线| 777久久人妻少妇嫩草av网站| 好男人在线观看高清免费视频| 人人妻人人澡欧美一区二区| 国产av一区在线观看免费| 19禁男女啪啪无遮挡网站| 欧美 亚洲 国产 日韩一| 夜夜爽天天搞| ponron亚洲| 婷婷丁香在线五月| 十八禁网站免费在线| 国产亚洲精品一区二区www| av超薄肉色丝袜交足视频| 悠悠久久av| 可以免费在线观看a视频的电影网站| 他把我摸到了高潮在线观看| 国产成人精品久久二区二区免费| 一级片免费观看大全| 一进一出抽搐gif免费好疼| 欧美色欧美亚洲另类二区| 久久婷婷人人爽人人干人人爱| 成人欧美大片| 亚洲aⅴ乱码一区二区在线播放 | 久久这里只有精品19| 久久精品夜夜夜夜夜久久蜜豆 | 18禁国产床啪视频网站| 欧美高清成人免费视频www| 国产精品久久久久久人妻精品电影| 老汉色∧v一级毛片| 狂野欧美激情性xxxx| 国产精品免费一区二区三区在线| 午夜日韩欧美国产| 国产精品久久视频播放| 波多野结衣高清无吗| 欧美另类亚洲清纯唯美| 国产一区二区在线av高清观看| 日本熟妇午夜| 两性夫妻黄色片| 精品久久久久久,| 最近最新中文字幕大全电影3| 舔av片在线| 麻豆成人av在线观看| 亚洲成人久久爱视频| 午夜a级毛片| 久久 成人 亚洲| 色播亚洲综合网| 午夜福利高清视频| √禁漫天堂资源中文www| 亚洲中文av在线| 国产乱人伦免费视频| 国产麻豆成人av免费视频| 国产高清激情床上av| 精品久久久久久久人妻蜜臀av| 母亲3免费完整高清在线观看| 国内少妇人妻偷人精品xxx网站 | 国产视频一区二区在线看| 久久久久性生活片| 精品乱码久久久久久99久播| 欧美日韩国产亚洲二区| 成年人黄色毛片网站| 97人妻精品一区二区三区麻豆| 高潮久久久久久久久久久不卡| 亚洲欧美一区二区三区黑人| 少妇粗大呻吟视频| 法律面前人人平等表现在哪些方面| 久久精品综合一区二区三区| 国产乱人伦免费视频| 99久久国产精品久久久| 国产欧美日韩一区二区三| 香蕉久久夜色| 制服诱惑二区| 亚洲成人中文字幕在线播放| 99精品在免费线老司机午夜| 一个人观看的视频www高清免费观看 | 亚洲欧美精品综合一区二区三区| 一级作爱视频免费观看| 国产精品av视频在线免费观看| 精品高清国产在线一区| www.www免费av| 亚洲成人精品中文字幕电影| 国产成人av激情在线播放| 欧美成人性av电影在线观看| 999久久久国产精品视频| 亚洲第一电影网av| 久久中文字幕一级| 我的老师免费观看完整版| 美女扒开内裤让男人捅视频| 九九热线精品视视频播放| 亚洲人成网站在线播放欧美日韩| e午夜精品久久久久久久| 99精品欧美一区二区三区四区| 久久精品亚洲精品国产色婷小说| 最近在线观看免费完整版| 国内揄拍国产精品人妻在线| 真人做人爱边吃奶动态| 日韩高清综合在线| 丁香欧美五月| 成人精品一区二区免费| 嫁个100分男人电影在线观看| 国产在线观看jvid| 久久 成人 亚洲| 日韩高清综合在线| 97人妻精品一区二区三区麻豆| 亚洲精品av麻豆狂野| 1024视频免费在线观看| 久久国产乱子伦精品免费另类| 免费观看精品视频网站| 日本 欧美在线| 国产野战对白在线观看| 亚洲 欧美一区二区三区| 国产精品久久视频播放| 天天躁狠狠躁夜夜躁狠狠躁| 成人欧美大片| 99riav亚洲国产免费| 亚洲成人中文字幕在线播放| 成人精品一区二区免费| 91国产中文字幕| av欧美777| 亚洲国产欧美网| tocl精华| 黄色女人牲交| 超碰成人久久| 午夜激情福利司机影院| 99热只有精品国产| 亚洲色图av天堂| e午夜精品久久久久久久| 亚洲精品美女久久av网站| 久久国产精品影院| 成人国语在线视频| 久久久久久国产a免费观看| 亚洲熟女毛片儿| 成人亚洲精品av一区二区| 久久久精品欧美日韩精品| 我的老师免费观看完整版| 久久久久精品国产欧美久久久|