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

    基于卡爾曼濾波的Ad Hoc網(wǎng)絡(luò)節(jié)點(diǎn)間時(shí)間同步算法研究

    2021-01-15 14:27:50孫力鄒德財(cái)韓濤
    時(shí)間頻率學(xué)報(bào) 2020年4期
    關(guān)鍵詞:卡爾曼濾波時(shí)鐘偏差

    孫力,鄒德財(cái),韓濤

    基于卡爾曼濾波的Ad Hoc網(wǎng)絡(luò)節(jié)點(diǎn)間時(shí)間同步算法研究

    孫力1,2,3,鄒德財(cái)1,2,4,韓濤5

    (1. 中國(guó)科學(xué)院 國(guó)家授時(shí)中心,西安 710600;2. 中國(guó)科學(xué)院 精密導(dǎo)航定位與定時(shí)技術(shù)重點(diǎn)實(shí)驗(yàn)室,西安 710600;3. 中國(guó)科學(xué)院大學(xué),北京 100049;4. 中國(guó)科學(xué)院大學(xué) 天文與空間科學(xué)學(xué)院,北京 101048;5.西安理工大學(xué) 理學(xué)院,西安 713300)

    本文建立了Ad Hoc網(wǎng)絡(luò)節(jié)點(diǎn)的時(shí)鐘模型,設(shè)計(jì)了適用于網(wǎng)絡(luò)節(jié)點(diǎn)時(shí)鐘模型的Kalman濾波算法,利用時(shí)差噪聲的分形特征估計(jì)Kalman濾波的協(xié)方差矩陣,并對(duì)時(shí)鐘噪聲序列特征進(jìn)行分析。仿真結(jié)果表明,該算法不僅提高了濾波算法效率,同時(shí)能明顯降低網(wǎng)絡(luò)節(jié)點(diǎn)間的時(shí)鐘偏差,實(shí)現(xiàn)了Ad Hoc網(wǎng)絡(luò)節(jié)點(diǎn)間的高精度時(shí)間同步,同步精度優(yōu)于25ns。

    Ad Hoc;網(wǎng)絡(luò)節(jié)點(diǎn);分形特性;Kalman濾波

    0 引言

    近年來,無線通信技術(shù)取得了快速發(fā)展,傳統(tǒng)的無線移動(dòng)通信網(wǎng)絡(luò)已不能滿足某些特殊場(chǎng)合,如軍事、救災(zāi)搶險(xiǎn)等對(duì)自主定位的需求。在這一背景下,Ad Hoc網(wǎng)絡(luò)[1]應(yīng)運(yùn)而生,在網(wǎng)絡(luò)部署方面比較靈活,可以在一些沒有通信基礎(chǔ)設(shè)施的地區(qū)自由組網(wǎng)。Ad Hoc網(wǎng)絡(luò)所有節(jié)點(diǎn)都可以隨意進(jìn)入或退出,各節(jié)點(diǎn)之間的地位和功能是平等的[2]。但由于缺少標(biāo)準(zhǔn)時(shí)鐘為網(wǎng)絡(luò)中的節(jié)點(diǎn)提供時(shí)間參考,各個(gè)網(wǎng)絡(luò)節(jié)點(diǎn)時(shí)鐘不同步,其時(shí)間差異很有可能導(dǎo)致整個(gè)網(wǎng)絡(luò)架構(gòu)不能正常運(yùn)行,自主定位功能無法實(shí)現(xiàn)。因此,網(wǎng)絡(luò)節(jié)點(diǎn)間的時(shí)間同步是整個(gè)網(wǎng)絡(luò)系統(tǒng)有效運(yùn)行的重要保證。

    Ad Hoc網(wǎng)絡(luò)架構(gòu)如圖1所示。Ad Hoc網(wǎng)絡(luò)不依賴基礎(chǔ)網(wǎng)絡(luò)設(shè)施的建立,避免了基礎(chǔ)設(shè)施網(wǎng)絡(luò)的故障影響到整個(gè)網(wǎng)絡(luò)系統(tǒng)通信的問題,它是一個(gè)沒有中心的網(wǎng)絡(luò)節(jié)點(diǎn),所有的網(wǎng)絡(luò)節(jié)點(diǎn)都是平等的,沒有優(yōu)先級(jí)之分。在網(wǎng)絡(luò)系統(tǒng)結(jié)構(gòu)中,采用分布式的組網(wǎng)特征,網(wǎng)絡(luò)節(jié)點(diǎn)之間自由組網(wǎng),具有自主對(duì)等的特性,具有良好的抗毀性和健壯性。

    圖1 Ad Hoc網(wǎng)絡(luò)架構(gòu)圖

    1 Ad Hoc網(wǎng)絡(luò)時(shí)間同步算法

    當(dāng)前,許多學(xué)者在無線自組織網(wǎng)絡(luò)時(shí)間同步方面開展了研究并提出了各種算法。在這些不同的網(wǎng)絡(luò)節(jié)點(diǎn)時(shí)間同步算法中,S. Ganeriwal等人[3]提出的傳感器網(wǎng)絡(luò)定時(shí)同步算法(TPSN)、J. Elson等專家[4-5]提出的參考廣播同步算法(RBS)、B. Kesy等人[6]提出的防洪時(shí)間同步算法(FTSP)和Pings等人提出的延時(shí)測(cè)量時(shí)間同步算法(DMTS)[7]是比較典型的幾個(gè)方法。傳統(tǒng)網(wǎng)絡(luò)節(jié)點(diǎn)時(shí)間同步算法的特點(diǎn)如表1所示。

    表1 網(wǎng)絡(luò)時(shí)間同步算法

    在Ad Hoc網(wǎng)絡(luò)中節(jié)點(diǎn)具有多跳的特征。據(jù)此FTSP算法采用層級(jí)的網(wǎng)絡(luò)結(jié)構(gòu),這種網(wǎng)絡(luò)結(jié)構(gòu)在大規(guī)模的網(wǎng)絡(luò)節(jié)點(diǎn)下是比較適用的,多個(gè)接收者與根節(jié)點(diǎn)發(fā)送的廣播信息進(jìn)行同步。對(duì)于網(wǎng)絡(luò)結(jié)構(gòu)的動(dòng)態(tài)變化以及一些網(wǎng)絡(luò)節(jié)點(diǎn)的異常情況,F(xiàn)TSP算法的容錯(cuò)性都是比較好的。在時(shí)間同步的過程中利用線性回歸的方式不斷精確地計(jì)算偏移量,補(bǔ)償了網(wǎng)絡(luò)節(jié)點(diǎn)的時(shí)間偏差,從而獲得更高的同步精度。對(duì)比之下,F(xiàn)TSP算法更適合Ad Hoc網(wǎng)絡(luò)節(jié)點(diǎn)間的時(shí)間同步。

    2 Ad Hoc網(wǎng)絡(luò)時(shí)鐘模型

    式(5)中

    那么,離散的時(shí)鐘模型可以表示為

    3 Ad Hoc網(wǎng)絡(luò)時(shí)鐘模型的Kalman濾波方法

    Kalman濾波廣泛應(yīng)用于導(dǎo)航、傳感器數(shù)據(jù)融合、雷達(dá)系統(tǒng)以及導(dǎo)彈追蹤等。近來更被應(yīng)用于人臉識(shí)別、圖像分割等計(jì)算機(jī)圖像的處理。它也是一種經(jīng)常用來描述網(wǎng)絡(luò)節(jié)點(diǎn)時(shí)鐘差的方法,可以用來描述節(jié)點(diǎn)時(shí)鐘的噪聲狀態(tài)。預(yù)測(cè)與更新是在使用卡爾曼濾波進(jìn)行的兩個(gè)過程[15]。首先想要對(duì)當(dāng)前時(shí)刻的狀態(tài)做出估計(jì),就需要使用上一時(shí)刻的狀態(tài)對(duì)它進(jìn)行預(yù)測(cè)估計(jì),這個(gè)就是預(yù)測(cè)過程。然后想要優(yōu)化預(yù)測(cè)操作中預(yù)測(cè)的這個(gè)估計(jì)值,那就需要當(dāng)前時(shí)刻的觀測(cè)狀態(tài)量來對(duì)其進(jìn)行優(yōu)化,這個(gè)就是更新過程。整個(gè)卡爾曼濾波過程即:先預(yù)測(cè),做出先驗(yàn)估計(jì)。之后計(jì)算卡爾曼增益,再得到估計(jì)值,這個(gè)后驗(yàn)估計(jì)值就與真值更為接近[16]。

    假設(shè)在Ad Hoc網(wǎng)絡(luò)中,兩個(gè)節(jié)點(diǎn)的頻率是相同的。這個(gè)假設(shè)與整個(gè)網(wǎng)絡(luò)架構(gòu)的實(shí)際情況是不相符的,但是通過時(shí)鐘模型的Kalman濾波方法可以改進(jìn)頻率不同的情況。則基于上述的假設(shè)條件,網(wǎng)絡(luò)節(jié)點(diǎn)間時(shí)差的Kalman濾波的數(shù)學(xué)表達(dá)式為

    狀態(tài)轉(zhuǎn)移矩陣和噪聲的協(xié)方差矩陣分別表示為:

    則:

    由上述可以知道,對(duì)協(xié)方差矩陣的估計(jì)、節(jié)點(diǎn)時(shí)鐘譜密度的估計(jì)以及時(shí)鐘漂移率譜密度的估計(jì)是使用Kalman濾波算法處理節(jié)點(diǎn)時(shí)鐘偏差的關(guān)鍵,可以參考系統(tǒng)時(shí)鐘噪聲序列特征對(duì)這些數(shù)值進(jìn)行估計(jì)。

    3.1 時(shí)鐘噪聲序列特征分析

    對(duì)時(shí)鐘噪聲序列進(jìn)行特征分析,可以使得節(jié)點(diǎn)時(shí)鐘模型能夠通過數(shù)學(xué)方法更準(zhǔn)確的描述,是建立更精確的時(shí)鐘模型、提高網(wǎng)絡(luò)間時(shí)間同步精度和效率的前提。

    3.1.1 時(shí)差序列獲取

    為了更加準(zhǔn)確地利用數(shù)學(xué)方法描述節(jié)點(diǎn)時(shí)鐘模型,必須對(duì)其數(shù)據(jù)特征進(jìn)行描述。為此,本文基于中國(guó)科學(xué)院國(guó)家授時(shí)中心的1 PPS時(shí)間基準(zhǔn)對(duì)3個(gè)網(wǎng)絡(luò)節(jié)點(diǎn)晶振的時(shí)間求取時(shí)差序列,分析其噪聲特性。由式(8)可知,對(duì)這3個(gè)觀測(cè)量進(jìn)行線性擬合,其擬合殘差描述了它們各自的噪聲情況。殘差序列的數(shù)字特征如表2所示。

    表2 節(jié)點(diǎn)時(shí)鐘偏差殘差序列數(shù)字特征 單位:s

    3.1.2 噪聲序列的非正態(tài)性

    對(duì)噪聲序列的分析是建立更精確的時(shí)鐘模型的前提,是提高網(wǎng)絡(luò)間時(shí)間同步精度和效率的必要前提。從它們的頻數(shù)直方圖(圖2)中很難看出它們是否服從正態(tài)分布。

    圖2 3個(gè)網(wǎng)絡(luò)節(jié)點(diǎn)晶振時(shí)差序列頻數(shù)統(tǒng)計(jì)直方圖

    本節(jié)將通過雅各-貝拉檢驗(yàn)(Jarque-Bera)來驗(yàn)證本文所采用的網(wǎng)絡(luò)節(jié)點(diǎn)時(shí)鐘數(shù)據(jù)是否具有符合正態(tài)分布的偏度和峰度的擬合優(yōu)度的檢驗(yàn)。

    構(gòu)造JB統(tǒng)計(jì)量[13]:

    本文的3個(gè)網(wǎng)絡(luò)時(shí)間節(jié)點(diǎn)時(shí)差噪聲序列的JB統(tǒng)計(jì)結(jié)果如表3所示。其中是判斷該檢驗(yàn)是否接受原假設(shè)的值,當(dāng)=0時(shí),認(rèn)為數(shù)據(jù)服從正態(tài)分布,接收原假設(shè);當(dāng)=1時(shí),則拒絕原假設(shè),認(rèn)為數(shù)據(jù)不服從正態(tài)分布。,JBSTAT和CV分別指檢驗(yàn)的概率值、檢驗(yàn)統(tǒng)計(jì)量的值和判斷是否拒絕原假設(shè)的關(guān)鍵值。由于各個(gè)殘差序列的JB統(tǒng)計(jì)結(jié)果均在顯著性水平0.05時(shí)拒絕正態(tài)分布假設(shè),因此,網(wǎng)絡(luò)節(jié)點(diǎn)時(shí)鐘的噪聲是非正態(tài)的,需要進(jìn)一步討論和分析。

    表3 節(jié)點(diǎn)時(shí)鐘偏差殘差JB統(tǒng)計(jì)量檢驗(yàn)結(jié)果(顯著性水平α = 0.05)

    3.2 時(shí)差噪聲的分形特征

    ① 計(jì)算均值:

    圖3 網(wǎng)絡(luò)節(jié)點(diǎn)時(shí)鐘偏差的分析結(jié)果

    其中,3個(gè)網(wǎng)絡(luò)節(jié)點(diǎn)時(shí)鐘晶振噪聲的Hurst指數(shù)線性擬合表達(dá)式為

    式(27)的計(jì)算結(jié)果表明,3個(gè)網(wǎng)絡(luò)節(jié)點(diǎn)時(shí)鐘晶振噪聲序列的Hurst指數(shù)均大于0.5,具有分形特性,是分形序列。

    由式(10)可知,這3個(gè)噪聲序列的傅里葉頻譜指數(shù)分別為:2.984 0、2.997 2和2.984 6。將這些計(jì)算結(jié)果代入式(15),即可對(duì)Kalman濾波的協(xié)方差矩陣進(jìn)行估計(jì)。

    通過上述分析,得到了網(wǎng)絡(luò)節(jié)點(diǎn)時(shí)鐘模型、Kalman濾波算法以及分形特征參數(shù)。其中,分形特性可以用來估計(jì)Kalman濾波的協(xié)方差矩陣,從而提高濾波算法效率。

    4 仿真分析

    本文通過軟件Matlab來進(jìn)行移動(dòng)自組織網(wǎng)絡(luò)的時(shí)間同步算法的仿真驗(yàn)證,網(wǎng)絡(luò)節(jié)點(diǎn)時(shí)鐘的本振頻率為10 MHz,仿真時(shí)間為180 s,同步周期設(shè)置為1 s[19]。通過仿真可以得到卡爾曼濾波前的時(shí)鐘偏差情況,其結(jié)果如圖4所示。

    圖4 卡爾曼濾波前的時(shí)鐘偏差

    從圖4可以看到卡爾曼濾波前的時(shí)鐘偏差是隨機(jī)抖動(dòng)的,通過分析可知,時(shí)鐘偏差的均值和方差分別為5.051 4×10-6s和1.229 6×10-13s??柭鼮V波后的時(shí)鐘偏差情況,其結(jié)果如圖5所示。

    圖5 卡爾曼濾波后的時(shí)鐘偏差

    從圖5可以看到卡爾曼濾波后的時(shí)鐘偏差也是隨機(jī)的,對(duì)仿真數(shù)據(jù)進(jìn)行分析可知,時(shí)鐘偏差的均值和方差分別為2.375 7×10-8s和1.363 8×10-13s。

    通過仿真分析與計(jì)算,可以看出使用卡爾曼濾波使得節(jié)點(diǎn)的時(shí)鐘偏差得到了明顯的改善,對(duì)時(shí)鐘系統(tǒng)狀態(tài)模型進(jìn)行調(diào)整,使得節(jié)點(diǎn)間實(shí)現(xiàn)時(shí)間同步。

    5 結(jié)語

    本文對(duì)建立的Ad Hoc網(wǎng)絡(luò)的節(jié)點(diǎn)時(shí)鐘模型進(jìn)行了數(shù)學(xué)描述,利用分形方法對(duì)網(wǎng)絡(luò)節(jié)點(diǎn)時(shí)鐘噪聲進(jìn)行分析和研究。為了對(duì)節(jié)點(diǎn)時(shí)鐘的噪聲狀態(tài)進(jìn)行描述,采用Kalman濾波方法對(duì)時(shí)鐘噪聲序列特征進(jìn)行分析,其中包括了時(shí)差序列獲取和噪聲序列的非正態(tài)性分析?;诜中翁卣鞣治鼋Y(jié)果,估計(jì)Kalman濾波的協(xié)方差矩陣,從而提高了濾波算法效率。通過Matlab對(duì)文中的基于Kalman濾波算法對(duì)Ad Hoc網(wǎng)絡(luò)節(jié)點(diǎn)間同步精度進(jìn)行仿真,結(jié)果表明該算法能明顯降低網(wǎng)絡(luò)節(jié)點(diǎn)間的時(shí)鐘偏差,基于該模型實(shí)現(xiàn)了Ad Hoc網(wǎng)絡(luò)節(jié)點(diǎn)間的高精度時(shí)間同步,同步精度優(yōu)于25 ns。

    [1] ZHENG S R, WANG H T, ZHAO Z F, et al. Ad Hoc networks technology[M]. Beijing: Posts & Telecom Press, 2005.

    [2] 楊海東, 鄧勇.Ad Hoc網(wǎng)絡(luò)時(shí)間同步技術(shù)研究[J]. 計(jì)算機(jī)工程與設(shè)計(jì), 2008, 29(21): 5447-5450.

    [3] GANERIWAL S, KUMAR R, SRIVASTAVA M B. Timing-syncprotocol for sensor networks[C] // SenSys’03: Proceedings of the1st International Conference on Embedded Networked Sensor Systems, New York: ACM, 2003: 138-149.

    [4] 李文鋒, 王汝傳, 孫力娟. 基于RBS的無線傳感器網(wǎng)絡(luò)時(shí)間同步算法[J]. 通信學(xué)報(bào), 2008, 29(6): 82-86.

    [5] 葉玲, 李太華, 代學(xué)武. 無線傳感器網(wǎng)絡(luò)環(huán)境下基于卡爾曼濾波的PTP協(xié)議[J]. 智能系統(tǒng)學(xué)報(bào), 2012, 7(6): 518-524.

    [6] ELSON J, GIROD L, ESTRIN D. Fine-grained network time synchronization using reference broadcasts[J]. ACM SIGOPS Operating Systems Review, 2002, 36(1): 147-163.

    [7] CHEN L, LENEUTRE J. Toward secure and scalable time synchronization in Ad Hoc networks[M]. Elsevier Science Publishers B: Computer Communications, 2007.

    [8] 宋曉雪. 無線自組織網(wǎng)絡(luò)時(shí)間同步技術(shù)研究[D]. 成都: 電子科技大學(xué), 2018.

    [9] 袁海波, 武福平, 董紹武. 原子鐘信號(hào)主要噪聲類別的判定與分析[J]. 儀器儀表學(xué)報(bào), 2008, 29(8): 187-191.

    [10] 郭海榮, 郭樹人, 焦文海, 等. 原子鐘時(shí)域頻率穩(wěn)定性計(jì)算方法[J]. 飛行器測(cè)控學(xué)報(bào), 2007, 26(2): 47-52.

    [11] 張敏. 原子鐘噪聲類型和頻率穩(wěn)定度估計(jì)的自由度分析與探討[D]. 西安: 中國(guó)科學(xué)院國(guó)家授時(shí)中心, 2008.

    [12] 陳伊卿. 無線傳感器網(wǎng)絡(luò)時(shí)間同步算法研究[D]. 西安: 西安電子科技大學(xué), 2011.

    [13] HUANG G, ZOMAYA A Y, DELICATO F C. An accurate on-demand time synchronization protocol for wireless sensor networks[J]. Journal of Parallel & Distributed Computing, 2012, 72(10): 1332-1346.

    [14] MAKARAVA N, HOLSCHNEIDER M. Estimation of the Hurst exponent from noisy data: a Bayesian approach[J]. The European Physical Journal B, 2012, 85(8): 272-471.

    [15] 顧仁財(cái), 劉飛. 自適應(yīng)Kalman濾波的戰(zhàn)術(shù)數(shù)據(jù)鏈自主時(shí)間同步算法[J]. 火力與指揮控制, 2019, 44(2): 76-79.

    [16] 莊曉燕, 王厚軍. 基于卡爾曼濾波器的IEEE 1588時(shí)鐘同步算法[J]. 電子測(cè)量與儀器學(xué)報(bào), 2012, 26(9): 747-751.

    [17] SARANN L, CATHERINE C, AURORE D. Different methods for spatial interpolation of rainfall data for operational hydrology and hydrological modeling at watershed scale. A review[J]. Biotechnology Agronomy Society & Environment, 2013, 17(2): 392-406.

    [18] HASAN M K, SAEED R A, HASHIM A H A, et al. Femtocell network time synchronization protocols and schemes[J]. Research Journal of Applied Science, Engineering and Technology, 2012, 4(23): 5136-5143.

    [19] 胡愛華, 鄧中亮, 張耀. 基于改進(jìn)TPSN和卡爾曼濾波的時(shí)間同步算法[J]. 現(xiàn)代電子技術(shù), 2018, 41(13): 5-9.

    Study on time synchronization algorithm of Ad Hoc networks based on Kalman filter

    SUN Li1,2,3, ZOU De-cai1,2,4, HAN Tao5

    (1. National Time Service Center, Chinese Academy of Sciences, Xi’an 710600, China;2. Key Laboratory of Precise Positioning and Timing Technology, Chinese Academy of Sciences, Xi’an 710600, China;3. University of Chinese Academy of Sciences, Beijing 100049, China;4. School of Astronomy and Space Science, University of Chinese Academy of Sciences, Beijing 101048, China;5.Faculty of Sciences, Xi’an University of Technology, Xi’an 713300, China)

    In this study, the clock model of the Ad Hoc network node is established, the Kalman filter algorithm which suitable for the network node clock model is designed, the covariance matrix of the Kalman filter is estimated using the fractal characteristics of the time difference noise, and the characteristics of the clock noise sequence are analyzed. The simulation results shown that the filtering algorithm not only can improve the efficiency, but also can significantly reduce the clock deviation between network nodes, and achieveshigh-precision time synchronization between Ad Hoc network nodes. The accuracy of synchronization is better than 25 ns.

    Ad Hoc; network node; fractal characteristics; Kalman filter

    10.13875/j.issn.1674-0637.2020-04-0279-09

    孫力, 鄒德財(cái), 韓濤. 基于卡爾曼濾波的Ad Hoc網(wǎng)絡(luò)節(jié)點(diǎn)間時(shí)間同步算法研究[J]. 時(shí)間頻率學(xué)報(bào), 2020, 43(4): 279-287.

    2020-04-17;

    2020-05-06

    載人航天創(chuàng)新性預(yù)研資助項(xiàng)目(17680601)

    猜你喜歡
    卡爾曼濾波時(shí)鐘偏差
    別樣的“時(shí)鐘”
    古代的時(shí)鐘
    如何走出文章立意偏差的誤區(qū)
    兩矩形上的全偏差
    基于遞推更新卡爾曼濾波的磁偶極子目標(biāo)跟蹤
    有趣的時(shí)鐘
    時(shí)鐘會(huì)開“花”
    基于模糊卡爾曼濾波算法的動(dòng)力電池SOC估計(jì)
    關(guān)于均數(shù)與偏差
    基于擴(kuò)展卡爾曼濾波的PMSM無位置傳感器控制
    成人18禁高潮啪啪吃奶动态图| 欧美日韩亚洲综合一区二区三区_| 99国产综合亚洲精品| 国产色视频综合| 大码成人一级视频| 国产高清国产精品国产三级| 日韩 亚洲 欧美在线| 18禁观看日本| 女性生殖器流出的白浆| 久久香蕉激情| 欧美久久黑人一区二区| 日韩中文字幕欧美一区二区| 美女中出高潮动态图| 人人妻人人澡人人爽人人夜夜| 欧美日韩视频精品一区| 精品少妇久久久久久888优播| 一级a爱视频在线免费观看| 国产在线一区二区三区精| 精品少妇黑人巨大在线播放| 久久精品亚洲av国产电影网| 淫妇啪啪啪对白视频 | 国产97色在线日韩免费| 精品卡一卡二卡四卡免费| 欧美 日韩 精品 国产| 欧美老熟妇乱子伦牲交| 久久精品亚洲av国产电影网| 黄色怎么调成土黄色| 久久人人爽人人片av| 日韩精品免费视频一区二区三区| 啪啪无遮挡十八禁网站| 多毛熟女@视频| 99国产精品一区二区三区| 97在线人人人人妻| 免费观看av网站的网址| 午夜激情久久久久久久| 丝袜脚勾引网站| 国产日韩欧美在线精品| 电影成人av| 色综合欧美亚洲国产小说| 91九色精品人成在线观看| 久久人人97超碰香蕉20202| 日本黄色日本黄色录像| 国产成人av教育| 亚洲精品久久午夜乱码| 国产一级毛片在线| 日韩三级视频一区二区三区| 欧美日韩一级在线毛片| 亚洲成人手机| 午夜福利视频在线观看免费| 狠狠狠狠99中文字幕| 十分钟在线观看高清视频www| 搡老熟女国产l中国老女人| 黑人巨大精品欧美一区二区蜜桃| 黑丝袜美女国产一区| 国产成人av教育| 热re99久久国产66热| 18禁观看日本| 免费观看人在逋| 桃花免费在线播放| 精品久久久精品久久久| 久久狼人影院| 一级片免费观看大全| 国产激情久久老熟女| 国产老妇伦熟女老妇高清| 啦啦啦免费观看视频1| 欧美午夜高清在线| 啦啦啦 在线观看视频| 19禁男女啪啪无遮挡网站| 久久久水蜜桃国产精品网| 十八禁网站免费在线| 国产免费现黄频在线看| 久久久国产精品麻豆| 狠狠狠狠99中文字幕| 亚洲熟女精品中文字幕| 久久久久久亚洲精品国产蜜桃av| 国产片内射在线| 制服人妻中文乱码| 国产男人的电影天堂91| www.精华液| 老司机午夜福利在线观看视频 | 国产成人欧美在线观看 | 制服人妻中文乱码| 每晚都被弄得嗷嗷叫到高潮| 精品久久久久久久毛片微露脸 | 妹子高潮喷水视频| 女性生殖器流出的白浆| 首页视频小说图片口味搜索| 丝袜脚勾引网站| 另类亚洲欧美激情| 丝袜人妻中文字幕| 久久青草综合色| 午夜精品国产一区二区电影| 亚洲av日韩精品久久久久久密| 国产成人一区二区三区免费视频网站| 国产成人精品久久二区二区免费| 国产又爽黄色视频| 欧美精品人与动牲交sv欧美| 久久人人爽av亚洲精品天堂| 美女视频免费永久观看网站| 国产在线免费精品| 狂野欧美激情性xxxx| 精品少妇内射三级| av天堂在线播放| 国产不卡av网站在线观看| 久久精品人人爽人人爽视色| 欧美日韩成人在线一区二区| 黑人猛操日本美女一级片| 少妇 在线观看| 欧美午夜高清在线| 99热国产这里只有精品6| 日日爽夜夜爽网站| 国产在线观看jvid| 国产一区二区三区av在线| 国产99久久九九免费精品| 日韩欧美国产一区二区入口| 久久99一区二区三区| 色综合欧美亚洲国产小说| 久久精品亚洲熟妇少妇任你| 18禁国产床啪视频网站| 99热国产这里只有精品6| 91九色精品人成在线观看| 在线观看一区二区三区激情| 最新的欧美精品一区二区| 两个人看的免费小视频| 精品人妻1区二区| 成年美女黄网站色视频大全免费| 法律面前人人平等表现在哪些方面 | 久久人人爽人人片av| 亚洲精品一二三| 亚洲av日韩在线播放| 国产97色在线日韩免费| 亚洲专区国产一区二区| 岛国毛片在线播放| 爱豆传媒免费全集在线观看| 777久久人妻少妇嫩草av网站| 久久久精品免费免费高清| 久久久久久久久久久久大奶| 精品少妇久久久久久888优播| 久久久久久免费高清国产稀缺| 亚洲国产日韩一区二区| 夜夜骑夜夜射夜夜干| 国产成人系列免费观看| 高清黄色对白视频在线免费看| 成人影院久久| 女人被躁到高潮嗷嗷叫费观| 国产免费福利视频在线观看| 久久久久久久精品精品| 日韩制服丝袜自拍偷拍| 最近最新免费中文字幕在线| 国产免费一区二区三区四区乱码| 欧美变态另类bdsm刘玥| 国产精品秋霞免费鲁丝片| 国产精品久久久久成人av| 欧美黑人精品巨大| 久久国产精品影院| 免费人妻精品一区二区三区视频| 免费人妻精品一区二区三区视频| 欧美日韩福利视频一区二区| 99热网站在线观看| 久久久久久久久免费视频了| 在线十欧美十亚洲十日本专区| 国产高清视频在线播放一区 | 国产色视频综合| 多毛熟女@视频| 老司机深夜福利视频在线观看 | 精品一区二区三区四区五区乱码| 国产精品免费大片| 国产成人精品久久二区二区免费| 亚洲精品中文字幕在线视频| 亚洲男人天堂网一区| 久久天躁狠狠躁夜夜2o2o| av又黄又爽大尺度在线免费看| 亚洲av美国av| 色婷婷av一区二区三区视频| 亚洲欧美日韩另类电影网站| 国产精品九九99| 亚洲综合色网址| 午夜两性在线视频| 亚洲精品国产精品久久久不卡| 欧美黄色淫秽网站| 另类精品久久| 国产男女超爽视频在线观看| 欧美精品av麻豆av| 亚洲少妇的诱惑av| 亚洲激情五月婷婷啪啪| 精品少妇内射三级| 高清视频免费观看一区二区| 婷婷丁香在线五月| 韩国高清视频一区二区三区| 最近最新中文字幕大全免费视频| 一级毛片精品| 国产免费av片在线观看野外av| 亚洲一卡2卡3卡4卡5卡精品中文| a级毛片黄视频| 高清视频免费观看一区二区| 免费久久久久久久精品成人欧美视频| 999久久久国产精品视频| 欧美xxⅹ黑人| 中文字幕av电影在线播放| 波多野结衣av一区二区av| av超薄肉色丝袜交足视频| 日日爽夜夜爽网站| 国产91精品成人一区二区三区 | 亚洲国产欧美网| 男女无遮挡免费网站观看| 欧美成狂野欧美在线观看| 精品久久久精品久久久| 亚洲,欧美精品.| 极品人妻少妇av视频| 午夜激情av网站| 国产欧美亚洲国产| 久久午夜综合久久蜜桃| 精品亚洲成国产av| 国产精品成人在线| 老司机亚洲免费影院| 高潮久久久久久久久久久不卡| 亚洲av片天天在线观看| 看免费av毛片| 97在线人人人人妻| 欧美日韩视频精品一区| 五月天丁香电影| 在线精品无人区一区二区三| 美女扒开内裤让男人捅视频| 亚洲精品一区蜜桃| 国产av又大| 亚洲午夜精品一区,二区,三区| 日韩,欧美,国产一区二区三区| 久久久久国产一级毛片高清牌| 亚洲精品在线美女| 亚洲黑人精品在线| netflix在线观看网站| 久久中文看片网| 亚洲成国产人片在线观看| 999久久久精品免费观看国产| 日韩欧美免费精品| 99久久综合免费| 在线永久观看黄色视频| av又黄又爽大尺度在线免费看| 悠悠久久av| 狠狠狠狠99中文字幕| 一级毛片精品| 91精品三级在线观看| 国产亚洲欧美精品永久| 99国产极品粉嫩在线观看| 久久精品亚洲av国产电影网| 国产精品99久久99久久久不卡| 一个人免费在线观看的高清视频 | 欧美97在线视频| 精品国产乱码久久久久久小说| 曰老女人黄片| 精品免费久久久久久久清纯 | 国产日韩一区二区三区精品不卡| 国产在线观看jvid| 色婷婷av一区二区三区视频| 12—13女人毛片做爰片一| 亚洲一码二码三码区别大吗| 丝袜美足系列| 欧美激情极品国产一区二区三区| 精品一区二区三区av网在线观看 | 亚洲国产精品999| av又黄又爽大尺度在线免费看| 天天躁夜夜躁狠狠躁躁| 国产成人影院久久av| 亚洲va日本ⅴa欧美va伊人久久 | 国产欧美日韩综合在线一区二区| 国产精品欧美亚洲77777| 啦啦啦在线免费观看视频4| 99热全是精品| 深夜精品福利| 精品国产乱码久久久久久小说| 免费一级毛片在线播放高清视频 | 黄网站色视频无遮挡免费观看| 91精品三级在线观看| 国产区一区二久久| 曰老女人黄片| 亚洲精品粉嫩美女一区| 女人精品久久久久毛片| 精品久久久精品久久久| 老熟妇仑乱视频hdxx| 亚洲精品国产精品久久久不卡| 国产成人av教育| 欧美激情久久久久久爽电影 | 国产精品影院久久| 90打野战视频偷拍视频| 黄色片一级片一级黄色片| 美女福利国产在线| 国产欧美日韩一区二区三区在线| 99国产精品免费福利视频| 99精品久久久久人妻精品| 中文字幕高清在线视频| 黑人巨大精品欧美一区二区蜜桃| 国产av精品麻豆| 国产男人的电影天堂91| 国产精品99久久99久久久不卡| 精品人妻熟女毛片av久久网站| e午夜精品久久久久久久| 91老司机精品| 中文字幕制服av| 黑丝袜美女国产一区| 老司机福利观看| 久久精品亚洲熟妇少妇任你| 欧美97在线视频| 国产精品久久久久久人妻精品电影 | 亚洲精品一二三| 男女边摸边吃奶| 欧美大码av| 亚洲精品在线美女| 精品亚洲成a人片在线观看| 久久国产精品大桥未久av| 久久久国产一区二区| 国产无遮挡羞羞视频在线观看| 亚洲精品国产精品久久久不卡| av片东京热男人的天堂| 在线 av 中文字幕| 日韩欧美一区视频在线观看| 中亚洲国语对白在线视频| 人妻久久中文字幕网| 免费在线观看视频国产中文字幕亚洲 | 久久毛片免费看一区二区三区| 欧美精品人与动牲交sv欧美| 大码成人一级视频| 亚洲欧美日韩另类电影网站| 捣出白浆h1v1| bbb黄色大片| www日本在线高清视频| 深夜精品福利| 一本—道久久a久久精品蜜桃钙片| 岛国在线观看网站| 日韩精品免费视频一区二区三区| 日韩视频一区二区在线观看| 十分钟在线观看高清视频www| 9热在线视频观看99| av视频免费观看在线观看| 久久精品亚洲av国产电影网| 精品一区在线观看国产| 久久人人爽人人片av| 久久久欧美国产精品| 国产精品久久久久久人妻精品电影 | 国产精品九九99| 三级毛片av免费| av不卡在线播放| 亚洲国产精品999| 亚洲av男天堂| 国产亚洲av高清不卡| 啦啦啦 在线观看视频| 十分钟在线观看高清视频www| 9热在线视频观看99| kizo精华| 18禁黄网站禁片午夜丰满| www日本在线高清视频| 亚洲久久久国产精品| 一本—道久久a久久精品蜜桃钙片| 亚洲欧美日韩另类电影网站| av有码第一页| 欧美一级毛片孕妇| 亚洲精品中文字幕在线视频| 免费观看人在逋| 国产一区二区 视频在线| 嫩草影视91久久| 人成视频在线观看免费观看| 淫妇啪啪啪对白视频 | 51午夜福利影视在线观看| 国产精品一区二区免费欧美 | 久久久久久亚洲精品国产蜜桃av| 亚洲五月色婷婷综合| av有码第一页| 国产精品免费大片| 三上悠亚av全集在线观看| 亚洲精品中文字幕一二三四区 | 中文字幕高清在线视频| 国产主播在线观看一区二区| 亚洲精品久久久久久婷婷小说| 丝袜美足系列| 久久久精品免费免费高清| 亚洲欧美成人综合另类久久久| 日韩 亚洲 欧美在线| 国产三级黄色录像| 欧美国产精品一级二级三级| 天天躁夜夜躁狠狠躁躁| 久久人妻熟女aⅴ| 亚洲欧美日韩高清在线视频 | 久久久欧美国产精品| 极品人妻少妇av视频| 久久国产精品影院| 午夜福利在线免费观看网站| 国产激情久久老熟女| 纯流量卡能插随身wifi吗| 国产成人精品久久二区二区91| 精品一区二区三区av网在线观看 | 悠悠久久av| 人成视频在线观看免费观看| 国产1区2区3区精品| 国产精品久久久人人做人人爽| 黄色片一级片一级黄色片| 一二三四在线观看免费中文在| 国产一区二区三区在线臀色熟女 | 老司机在亚洲福利影院| 国产免费现黄频在线看| 国产精品欧美亚洲77777| 97在线人人人人妻| 欧美成狂野欧美在线观看| 女性被躁到高潮视频| 亚洲精品美女久久久久99蜜臀| 老司机福利观看| 国产深夜福利视频在线观看| 午夜福利在线观看吧| 一本综合久久免费| 亚洲成av片中文字幕在线观看| 人妻久久中文字幕网| 伊人久久大香线蕉亚洲五| 久久久久久久久久久久大奶| 亚洲国产欧美网| 99精品欧美一区二区三区四区| 中文字幕色久视频| 12—13女人毛片做爰片一| 欧美在线一区亚洲| netflix在线观看网站| 80岁老熟妇乱子伦牲交| 日韩制服骚丝袜av| 精品国产超薄肉色丝袜足j| 成在线人永久免费视频| 亚洲欧洲精品一区二区精品久久久| 成年av动漫网址| 视频在线观看一区二区三区| 91九色精品人成在线观看| 一二三四社区在线视频社区8| √禁漫天堂资源中文www| 乱人伦中国视频| 久久九九热精品免费| 手机成人av网站| 岛国毛片在线播放| 国产精品一区二区在线不卡| 国产麻豆69| 国产男女内射视频| 女人高潮潮喷娇喘18禁视频| 99精品久久久久人妻精品| 老司机深夜福利视频在线观看 | 亚洲 欧美一区二区三区| 亚洲视频免费观看视频| 一级片'在线观看视频| 国产日韩欧美亚洲二区| 视频区图区小说| 亚洲专区中文字幕在线| 99精品久久久久人妻精品| 老汉色∧v一级毛片| 中亚洲国语对白在线视频| 手机成人av网站| 男女免费视频国产| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品久久久久成人av| 岛国毛片在线播放| 久久久久久久久免费视频了| 欧美黑人精品巨大| 日韩中文字幕视频在线看片| 99国产极品粉嫩在线观看| 午夜精品久久久久久毛片777| 欧美成人午夜精品| 欧美在线一区亚洲| 中文字幕精品免费在线观看视频| 性高湖久久久久久久久免费观看| 丰满人妻熟妇乱又伦精品不卡| 午夜影院在线不卡| 亚洲国产av新网站| 亚洲国产成人一精品久久久| 久久性视频一级片| 亚洲专区字幕在线| 国产欧美日韩精品亚洲av| 淫妇啪啪啪对白视频 | 香蕉丝袜av| 色播在线永久视频| 视频区欧美日本亚洲| 日日夜夜操网爽| 久久天堂一区二区三区四区| 午夜久久久在线观看| 久久国产精品男人的天堂亚洲| 两人在一起打扑克的视频| 成人av一区二区三区在线看 | 亚洲成人免费电影在线观看| 最新的欧美精品一区二区| 久热爱精品视频在线9| 最近最新中文字幕大全免费视频| 中文字幕高清在线视频| 久久av网站| 国产高清国产精品国产三级| 日韩精品免费视频一区二区三区| 成人国产av品久久久| 久久国产精品影院| 亚洲国产欧美一区二区综合| 老司机亚洲免费影院| 黄色视频,在线免费观看| 99久久99久久久精品蜜桃| 久久精品人人爽人人爽视色| 国精品久久久久久国模美| 国产一区二区三区综合在线观看| 精品国产乱子伦一区二区三区 | 成人18禁高潮啪啪吃奶动态图| 亚洲视频免费观看视频| 中文字幕最新亚洲高清| 精品人妻在线不人妻| 色综合欧美亚洲国产小说| 久久国产精品男人的天堂亚洲| 成年女人毛片免费观看观看9 | 免费在线观看黄色视频的| 99国产极品粉嫩在线观看| 亚洲avbb在线观看| 国产一区二区三区综合在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 两性午夜刺激爽爽歪歪视频在线观看 | 99国产精品一区二区三区| 欧美日韩成人在线一区二区| 国产日韩欧美亚洲二区| 欧美大码av| 黑丝袜美女国产一区| 久久久久国产精品人妻一区二区| 欧美日韩一级在线毛片| 巨乳人妻的诱惑在线观看| 丰满迷人的少妇在线观看| 精品国产一区二区三区久久久樱花| 97精品久久久久久久久久精品| 亚洲欧美成人综合另类久久久| 一区福利在线观看| 国产99久久九九免费精品| 国产免费福利视频在线观看| 妹子高潮喷水视频| 成人国产av品久久久| 国产成人精品久久二区二区91| 国产精品久久久久久精品古装| 操美女的视频在线观看| 亚洲第一青青草原| 欧美日韩国产mv在线观看视频| 久久免费观看电影| 久久精品国产a三级三级三级| 免费av中文字幕在线| 美女高潮到喷水免费观看| 9热在线视频观看99| 精品熟女少妇八av免费久了| h视频一区二区三区| 免费在线观看日本一区| 高清黄色对白视频在线免费看| 12—13女人毛片做爰片一| 午夜激情久久久久久久| 国产一区二区激情短视频 | 波多野结衣一区麻豆| 啦啦啦免费观看视频1| 久久青草综合色| 最近最新中文字幕大全免费视频| 女人被躁到高潮嗷嗷叫费观| 亚洲国产中文字幕在线视频| 日韩熟女老妇一区二区性免费视频| 国产精品99久久99久久久不卡| 黄色片一级片一级黄色片| 久久精品国产亚洲av香蕉五月 | 免费在线观看日本一区| 69av精品久久久久久 | 久久久国产精品麻豆| 欧美日韩精品网址| 日韩制服骚丝袜av| 精品乱码久久久久久99久播| 欧美 亚洲 国产 日韩一| 成人黄色视频免费在线看| 亚洲国产看品久久| 国产成人免费无遮挡视频| 欧美精品一区二区大全| 成年av动漫网址| av一本久久久久| 叶爱在线成人免费视频播放| 另类精品久久| 高清av免费在线| 午夜福利免费观看在线| 大片免费播放器 马上看| 免费女性裸体啪啪无遮挡网站| 男人爽女人下面视频在线观看| 男女边摸边吃奶| 建设人人有责人人尽责人人享有的| 中国美女看黄片| 99热全是精品| 久久天躁狠狠躁夜夜2o2o| 啦啦啦啦在线视频资源| 精品熟女少妇八av免费久了| av天堂在线播放| 久久久国产一区二区| 18禁国产床啪视频网站| 窝窝影院91人妻| 精品国产乱码久久久久久男人| 手机成人av网站| 制服诱惑二区| 十八禁高潮呻吟视频| 亚洲国产欧美在线一区| 欧美黑人精品巨大| 精品少妇久久久久久888优播| a级毛片黄视频| 满18在线观看网站| 久久久久久免费高清国产稀缺| 久久这里只有精品19| 久热这里只有精品99| 这个男人来自地球电影免费观看| 热re99久久精品国产66热6| 女警被强在线播放| 巨乳人妻的诱惑在线观看| 久久人人爽av亚洲精品天堂| 久久久久精品人妻al黑| 一本大道久久a久久精品| 亚洲精品自拍成人| 欧美日韩亚洲高清精品| 在线av久久热| 精品乱码久久久久久99久播| 美女福利国产在线| 日本av免费视频播放| 999久久久国产精品视频| 91麻豆av在线| kizo精华| 国产成人精品久久二区二区免费| 国内毛片毛片毛片毛片毛片| 国产精品国产三级国产专区5o| 美女扒开内裤让男人捅视频|