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

    基于高斯混合模型的無(wú)線傳感器網(wǎng)絡(luò)定位算法

    2022-07-07 01:58:22蔣水賓郝艷妮
    關(guān)鍵詞:三邊測(cè)距高斯

    方 省, 羅 引, 曹 家, 徐 楠, 蔣水賓, 郝艷妮

    (1. 天津大學(xué) 經(jīng)濟(jì)與管理學(xué)院, 天津 300072; 2. 中國(guó)科學(xué)院自動(dòng)化研究所, 北京 100190;3. 北京中科聞歌科技股份有限公司, 北京 100000)

    隨著物聯(lián)網(wǎng)[1-3]和智慧城市建設(shè)的快速發(fā)展, 人們對(duì)位置服務(wù)需求與日俱增. 無(wú)線傳感器網(wǎng)絡(luò)(wireless sensor networks, WSNs)位于物聯(lián)網(wǎng)(internet of things, IoT)的感知層, 節(jié)點(diǎn)的定位信息是物聯(lián)網(wǎng)實(shí)現(xiàn)和普及的基礎(chǔ)和前提[4]. 距離信息可先由TOA[5],AOA[6],TDOA[7-8],RSSI[9],DV-Hop[10-11]和DV-distance[12]等技術(shù)計(jì)算得出, 然后通過三邊測(cè)量法等定位方法計(jì)算未知節(jié)點(diǎn)的位置. 但由于非視距問題[13-14](NLOS)、 定位環(huán)境及路徑傳播模型等因素影響, 測(cè)量的距離通常存在一定的誤差[15].

    針對(duì)距離誤差對(duì)定位的影響, 目前已有許多方法可降低誤差的影響: 1) 根據(jù)定位環(huán)境改善路徑損耗模型參數(shù), 進(jìn)行實(shí)時(shí)估計(jì), 如文獻(xiàn)[16]提出了一種多模型選擇算法, 在分析路徑損耗模型統(tǒng)計(jì)特性的基礎(chǔ)上, 提出了一種進(jìn)行路徑損耗模型的量化方法, 文獻(xiàn)[17]提出了一種根據(jù)環(huán)境變化的RSS室內(nèi)定位方法; 2) 根據(jù)某些方法, 對(duì)于不同的測(cè)距信息給予不同的權(quán)值, 如文獻(xiàn)[18]引入了多個(gè)通信半徑以細(xì)分跳數(shù), 并推導(dǎo)了跳數(shù)隨半徑變化次數(shù)的公式, 通過校正加權(quán)值修正平均跳數(shù)距離對(duì)平均跳距離進(jìn)行加權(quán)處理, 文獻(xiàn)[19]采用最小均方誤差和逆距離加權(quán)法處理平均一跳距離; 3) 采用某種算法篩選出含有誤差的測(cè)距信息, 將其從距離信息中去除. 但在實(shí)際應(yīng)用中, 方法1)過于繁瑣和復(fù)雜, 方法2)效果不佳, 定位結(jié)果仍存在較大誤差.

    基于此, 本文提出一種基于高斯混合模型的無(wú)線傳感器網(wǎng)絡(luò)定位算法, 以解決不準(zhǔn)確距離信息對(duì)定位結(jié)果的影響. 通過構(gòu)建混合高斯模型并分析篩選出誤差較大的距離信息, 將其從距離信息集合中剔除, 再對(duì)剩余的距離信息集合利用三邊測(cè)量法和加權(quán)算法計(jì)算未知節(jié)點(diǎn)的位置. 仿真實(shí)驗(yàn)結(jié)果表明, 利用改進(jìn)算法對(duì)距離信息處理, 剔除大誤差距離信息后, 所得到的定位結(jié)果更精準(zhǔn), 誤差更小.

    1 高斯混合模型及其參數(shù)估計(jì)

    本文使用高斯混合模型(Gaussian mixture model, GMM)擬合由三邊測(cè)量法[20-22]得到的定位結(jié)果集合, 通過模型分析, 找出遠(yuǎn)離中心的定位結(jié)果及定位這些位置所使用的距離信息, 并刪除出現(xiàn)次數(shù)較多的距離信息, 利用剩余的距離信息進(jìn)行位置估計(jì). 假設(shè)p(y|θ)為定位結(jié)果的概率密度函數(shù), 則

    用高斯混合模型擬合定位結(jié)果就是將這些定位結(jié)果進(jìn)行分類, 各類別的概率密度函數(shù)可視為一個(gè)服從高斯分布的函數(shù).常用的高斯混合模型參數(shù)估計(jì)算法是最大似然估計(jì)和最大期望(EM)算法.

    算法1高斯混合模型參數(shù)估計(jì)的EM算法.

    輸入: 觀測(cè)數(shù)據(jù)y1,y2,…,yN, 高斯混合模型;

    輸出: 高斯混合模型參數(shù);

    步驟1) 取參數(shù)的初始值開始迭代;

    步驟2) E步: 根據(jù)當(dāng)前模型參數(shù), 計(jì)算分模型k對(duì)觀測(cè)數(shù)據(jù)yj的響應(yīng)度為

    步驟3) M步: 計(jì)算新一輪迭代的模型參數(shù)為

    步驟4) 重復(fù)步驟2)和步驟3), 直到|p(y|θ)-p(y|θ)′|<ε, 其中p(y|θ)和p(y|θ)′為相鄰兩次迭代的結(jié)果, 即前后兩次迭代的結(jié)果小于一定程度即可終止迭代, 通常ε=10-5.

    2 算法設(shè)計(jì)

    為證明高斯混合模型無(wú)線傳感器網(wǎng)絡(luò)定位算法的合理性和有效性, 本文先介紹三邊測(cè)量法并分析不準(zhǔn)確距離信息對(duì)定位結(jié)果的影響, 以證明提出的算法可提高定位精度.

    由于受環(huán)境和非視距問題(NLOS)等因素的影響, 測(cè)得的距離與真實(shí)距離通常有一定的誤差. 為全面分析不準(zhǔn)確信息對(duì)定位的影響, 分別從不準(zhǔn)距離信息的程度(即不準(zhǔn)距離信息所含誤差大小)及不準(zhǔn)距離信息多少(即不準(zhǔn)距離信息占總距離信息的比例)兩方面進(jìn)行分析和對(duì)比.

    圖1為當(dāng)不準(zhǔn)距離信息個(gè)數(shù)m=2時(shí), 定位結(jié)果的誤差隨不準(zhǔn)距離信息程度變化的情況. 由圖1可見, 當(dāng)ε>10%時(shí), 隨著ε的增大, 定位誤差(RMSE)也隨著增大; 當(dāng)ε>30%時(shí), 由于距離信息所含誤差過大, 導(dǎo)致三邊測(cè)量定位誤差顯著上升.

    圖1 不準(zhǔn)程度對(duì)三邊測(cè)量法的影響Fig.1 Effect of inaccuracy on trilateral measurement

    不準(zhǔn)距離信息多少對(duì)定位結(jié)果的影響是指不準(zhǔn)距離信息占全部距離信息比例的不同會(huì)導(dǎo)致定位結(jié)果的不同.為方便, 本文用η表示不準(zhǔn)信息占全部距離信息的比例, 仿真場(chǎng)景設(shè)置與不準(zhǔn)距離程度的仿真條件相同.仿真條件設(shè)置: 不準(zhǔn)程度ε=5%,10%,20%,30%,m的取值范圍為0≤m≤6, 即η從0變化到40%, 由于當(dāng)η>40%時(shí), 測(cè)距信息集中有近1/2的距離信息不準(zhǔn)確, 此時(shí)定位系統(tǒng)性能非常差, 故對(duì)m>6的情況進(jìn)行仿真無(wú)意義.

    圖2為不同測(cè)距誤差ε下不準(zhǔn)距離信息大小對(duì)三邊測(cè)量結(jié)果的影響.由圖2可見, 隨著不準(zhǔn)程度ε的增高, 定位誤差隨之增加; 在不準(zhǔn)程度ε確定的情況下, 隨著不準(zhǔn)距離信息數(shù)目的提高, 定位誤差也隨著增加; 剔除不準(zhǔn)距離信息后, 利用三邊測(cè)量法對(duì)剩余距離信息進(jìn)行定位, 定位結(jié)果較精確.因此, 不準(zhǔn)確距離信息增加了定位結(jié)果的誤差, 降低了定位的精準(zhǔn)度.從距離集合中剔除誤差大的距離信息, 利用剩余的距離信息進(jìn)行三邊測(cè)量法定位, 將會(huì)使定位結(jié)果更精確.

    圖2 不同測(cè)距誤差ε下不準(zhǔn)距離信息大小對(duì)三邊測(cè)量結(jié)果的影響Fig.2 Effect of inaccurate distance information on trilateral measurement results under different ranging errors ε

    三邊測(cè)量法需要3個(gè)距離信息, 當(dāng)距離信息集合中多數(shù)距離都在合理的誤差范圍內(nèi)時(shí), 通過三邊測(cè)量法得到的一系列定位結(jié)果之間的差距很小, 它們集中在一片密集的區(qū)域.但當(dāng)測(cè)距信息中含有較大誤差, 即通過三邊測(cè)量所對(duì)應(yīng)的3個(gè)距離信息只要有一個(gè)含有較大的誤差時(shí), 則所得到的位置將偏離該區(qū)域, 本文據(jù)此通過高斯混合模型分析對(duì)偏離密集區(qū)域的定位結(jié)果進(jìn)行分析, 剔出大誤差距離信息.

    高斯混合模型中K值的選擇對(duì)定位精度和運(yùn)行時(shí)間有影響, 子高斯模型過多會(huì)延長(zhǎng)算法的運(yùn)行時(shí)間, 對(duì)位置樣本點(diǎn)A中樣本的分析將更詳細(xì), 定位精度有所提高; 子高斯模型過少會(huì)縮短算法的運(yùn)行時(shí)間, 但模型過于簡(jiǎn)單, 對(duì)于位置樣本點(diǎn)A中樣本的分析會(huì)粗糙, 定位精度低.通過仿真實(shí)驗(yàn)可知選擇K=2最合適.

    定位結(jié)果高斯混合模型建立的過程: 概率密度函數(shù)為

    基于上述分析, 本文利用高斯混合模型對(duì)三邊測(cè)量法得到的一系列定位結(jié)果建模, 并對(duì)建模結(jié)果進(jìn)行分析, 可知集合B和C所含的定位結(jié)果偏離真實(shí)定位結(jié)果, 根據(jù)一個(gè)定位結(jié)果中含有3個(gè)距離信息, 分析可得B和C所含有的距離信息有較大誤差, 再通過反向分析可得出含有誤差較大的距離信息.將誤差較大測(cè)距信息從距離信息集合中剔除, 利用剩余距離信息集合進(jìn)行定位分析和計(jì)算.該方法必須滿足已知位置節(jié)點(diǎn)數(shù)目必須超過6個(gè), 這樣通過三邊測(cè)量方法將會(huì)產(chǎn)生足夠多的定位結(jié)果樣本, 而錨節(jié)點(diǎn)多于6個(gè)的情況在應(yīng)用中是普遍現(xiàn)象[19].

    加權(quán)定位算法表示為對(duì)定位集合A, 假設(shè)由任意3個(gè)節(jié)點(diǎn)通過三邊測(cè)量法得到的未知節(jié)點(diǎn)坐標(biāo)為(xn,yn), 其由di,dj,dk三個(gè)測(cè)距信息獲取, 針對(duì)定位集合A的加權(quán)定位算法[18]求得節(jié)點(diǎn)(x,y)為

    基于高斯混合模型的定位算法步驟如下:

    2) 根據(jù)高斯混合模型輸出的最終結(jié)果找出元素個(gè)數(shù)最少的子高斯模型φ2, 其含有的定位結(jié)果集合為B和距離φ1中心最遠(yuǎn)的定位結(jié)果集合C;

    3) 找出集合B和C中定位結(jié)果所包含的點(diǎn), 由于每個(gè)點(diǎn)對(duì)應(yīng)著3個(gè)距離信息, 通過分析找出所對(duì)應(yīng)距離信息集合中出現(xiàn)頻率最高的前m個(gè)距離信息, 假設(shè)d1,d2,…,dm即為誤差較大的距離信息;

    4) 對(duì)剩余的距離信息進(jìn)行三邊測(cè)量法, 得到一系列定位結(jié)果, 對(duì)其進(jìn)行高斯模型擬合, 得高斯模型的均值μ=(μx,μy), 對(duì)定位結(jié)果進(jìn)行加權(quán)定位, 用(x,y)表示;

    本文取δ=0.5 m.閾值的設(shè)置使算法適應(yīng)性更強(qiáng), 若m設(shè)置過小, 則不能將不準(zhǔn)確距離信息全部剔除, 需通過估計(jì)位置與高斯模型中心的歐氏距離判斷是否需要繼續(xù)剔除不準(zhǔn)確信息, 使得出的定位精度更高, 與真實(shí)值的位置更接近.

    3 仿真實(shí)驗(yàn)

    (1)

    3.1 不同測(cè)距誤差分布下算法的有效性

    仿真環(huán)境設(shè)置: 將錨節(jié)點(diǎn)隨機(jī)分布在100 m×100 m的監(jiān)測(cè)區(qū)域中, 未知節(jié)點(diǎn)U位于坐標(biāo)為(50,50)位置處, 假設(shè)未知節(jié)點(diǎn)能接收與15個(gè)錨節(jié)點(diǎn)的距離信息.

    在本文仿真場(chǎng)景中, 未知節(jié)點(diǎn)接收到的15個(gè)距離信息{d1,d2,…,d15}所包含的誤差εi(i=1,2,…,15)(當(dāng)ε>50%時(shí), 測(cè)距誤差較大, 對(duì)定位結(jié)果影響過大, 故仿真不予考慮)分別服從均勻分布和正態(tài)分布, 仿真結(jié)果如圖3所示.由圖3可見, 無(wú)論測(cè)距誤差服從均勻分布還是正態(tài)分布, 所得定位結(jié)果均較相近.從而驗(yàn)證了在不同測(cè)距誤差分布下本文算法都有效.

    圖3 不同測(cè)距誤差分布下定位結(jié)果的對(duì)比Fig.3 Comparison of positioning results under different ranging error distributions

    表1為計(jì)算中m分別取1,2,3, 高斯混合建模分析剔除誤差較大的1,2或3個(gè)信息(該仿真算法不包含基于高斯混合模型定位算法中的步驟5)進(jìn)行自適應(yīng)調(diào)整m值), 然后使用三邊測(cè)量法和加權(quán)算法求其定位結(jié)果.

    表1 當(dāng)測(cè)距誤差服從不同分布時(shí), 兩種算法定位誤差的仿真結(jié)果

    由表1可見, 無(wú)論測(cè)距誤差ε服從均勻分布還是正態(tài)分布, 改進(jìn)算法對(duì)于測(cè)距信息中不存在不準(zhǔn)距離信息的情況也適用.當(dāng)m增大時(shí), 去除誤差大的前m個(gè)測(cè)距信息, 隨著m的增大, 改進(jìn)算法所得的定位結(jié)果與真實(shí)定位結(jié)果更接近.

    3.2m值對(duì)定位結(jié)果的影響

    當(dāng)m設(shè)置過小時(shí), 算法會(huì)根據(jù)定位結(jié)果與聚類中心的歐氏距離進(jìn)行調(diào)節(jié), 其對(duì)定位結(jié)果的影響較小, 只需多次迭代花費(fèi)更多的時(shí)間; 如果m設(shè)置過大, 大于含有不準(zhǔn)距離信息的個(gè)數(shù)時(shí), 考慮是否會(huì)對(duì)定位結(jié)果產(chǎn)生影響.不失一般性, 仿真環(huán)境與3.1相同, 不準(zhǔn)確距離信息數(shù)目n=3,ε=20%, 分別取m=1,2,3,4,5, 然后使用三邊測(cè)量法和加權(quán)算法求其定位結(jié)果, 定位誤差仿真結(jié)果列于表2.

    表2 當(dāng)m取不同值時(shí)不同算法定位誤差的仿真結(jié)果

    由表2可見: 當(dāng)m=1,2,3時(shí), 定位結(jié)果非常接近, 驗(yàn)證了算法對(duì)于閾值δ設(shè)置的合理性; 當(dāng)m=4,5時(shí), 雖然大于不準(zhǔn)確距離信息的個(gè)數(shù), 但對(duì)于定位結(jié)果的影響較小, 誤差波動(dòng)不大.驗(yàn)證了在合理的取值范圍內(nèi),m的取值不會(huì)對(duì)定位結(jié)果有較大影響.

    3.3K值對(duì)定位結(jié)果的影響

    K值的大小將對(duì)定位精度和運(yùn)行時(shí)間產(chǎn)生影響, 為得出合適的K值進(jìn)行仿真對(duì)比. 仿真環(huán)境與3.1相同, 不準(zhǔn)確距離信息數(shù)目n=3,ε=20%.當(dāng)K值分別取2,3時(shí), 計(jì)算其RMSE, 結(jié)果列于表3. 由表3可見, 當(dāng)K=3時(shí), 定位誤差比K=2時(shí)小, 但其波動(dòng)浮動(dòng)較小, 而運(yùn)行時(shí)間延長(zhǎng)浮動(dòng)較大, 權(quán)衡兩者對(duì)定位的影響, 本文選擇K=2.

    表3 當(dāng)K=2,3時(shí)不同算法定位誤差的仿真結(jié)果

    3.4 改進(jìn)算法與三邊測(cè)量法的定位精度對(duì)比

    針對(duì)在不同誤差ε和不同大誤差測(cè)距信息數(shù)目的情況下, 為驗(yàn)證改進(jìn)算法的有效性和實(shí)用性, 進(jìn)行改進(jìn)定位算法和三邊測(cè)量法所得定位誤差的仿真對(duì)比, 結(jié)果如圖4所示. 仿真環(huán)境與3.1相同.

    圖4 定位誤差隨不準(zhǔn)距離信息個(gè)數(shù)的變化Fig.4 Positioning errors vary with number of inaccurate distance information

    由圖4可見: 隨著測(cè)距誤差增大, 大誤差測(cè)距信息數(shù)目增多, 三邊測(cè)量定位法受到的影響越大, 性能下降, 導(dǎo)致定位誤差越大, 精度更低; 改進(jìn)定位算法所得定位結(jié)果基本不隨大誤差測(cè)距信息的增大而增大, 在測(cè)距誤差逐漸增大和不準(zhǔn)測(cè)距信息數(shù)目逐漸增多的情況下, 改進(jìn)定位算法剔除大誤差測(cè)距信息后, 定位誤差有明顯減少, 提高了定位精度.

    3.5 DV-Hop算法的仿真分析

    為驗(yàn)證改進(jìn)定位算法應(yīng)用的廣泛性和合理性, 將改進(jìn)定位算法應(yīng)用到DV-Hop算法中, 對(duì)改進(jìn)DV-Hop算法與原始DV-Hop算法所得定位誤差進(jìn)行對(duì)比, 結(jié)果如圖5所示. 仿真環(huán)境參數(shù)設(shè)置為: 在80 m×80 m的區(qū)域上布置80個(gè)節(jié)點(diǎn), 通信半徑分別為20 m和30 m, 其中K=2, 考慮到錨節(jié)點(diǎn)數(shù)目的不同, 當(dāng)錨節(jié)點(diǎn)數(shù)目小于10時(shí)m設(shè)置為2, 當(dāng)錨節(jié)點(diǎn)數(shù)目大于10小于20時(shí)m設(shè)置為3, 當(dāng)錨節(jié)點(diǎn)數(shù)目大于20小于30時(shí)m設(shè)置為4, 當(dāng)錨節(jié)點(diǎn)數(shù)目大于30小于40時(shí)m設(shè)置為4.

    圖5為當(dāng)通信半徑R=20,30 m時(shí)改進(jìn)的DV-Hop算法與原DV-Hop算法的定位精度隨錨節(jié)點(diǎn)數(shù)目變化的情況. 由圖5可知: 通信半徑R=30 m大于R=20 m時(shí)的定位精度, 這是因?yàn)楫?dāng)通信半徑增大時(shí), 未知節(jié)點(diǎn)可通信的錨節(jié)點(diǎn)數(shù)目增多, 即得到的測(cè)距信息更多, 當(dāng)測(cè)距信息更多時(shí)可獲得的定位結(jié)果更精確; 改進(jìn)DV-Hop定位算法的定位誤差遠(yuǎn)小于原DV-Hop算法的定位誤差, 這是因?yàn)楦倪M(jìn)算法提高了定位精度, 使定位結(jié)果與真實(shí)定位結(jié)果更接近.

    圖5 當(dāng)R=20,30 m時(shí)不同算法定位精度隨錨節(jié)點(diǎn)個(gè)數(shù)的變化Fig.5 Positioning accuracy of different algorithms changes with number of anchor nodes when R=20,30 m

    綜上所述, 本文在分析不準(zhǔn)確距離信息對(duì)定位影響的基礎(chǔ)上, 結(jié)合高斯混合模型提出了一種提高定位精度的定位算法. 該算法通過高斯混合模型擬合定位結(jié)果, 剔除大誤差測(cè)距信息, 結(jié)合三邊測(cè)量定位法對(duì)剩余的距離信息集合得到定位估計(jì)集合, 再對(duì)定位估計(jì)集合進(jìn)行加權(quán)處理, 得到離未知節(jié)點(diǎn)真實(shí)位置較近的坐標(biāo). 仿真結(jié)果表明, 改進(jìn)算法削弱了不準(zhǔn)測(cè)距信息對(duì)定位結(jié)果的影響, 定位精度有顯著提高, 與真實(shí)位置更接近.

    猜你喜歡
    三邊測(cè)距高斯
    小高斯的大發(fā)現(xiàn)
    三角形中線與高之間的三個(gè)幾何不等式
    九點(diǎn)圓圓心關(guān)于三邊的對(duì)稱點(diǎn)的性質(zhì)
    走三邊
    類星體的精準(zhǔn)測(cè)距
    科學(xué)(2020年3期)2020-01-06 04:02:51
    天才數(shù)學(xué)家——高斯
    淺談超聲波測(cè)距
    電子制作(2017年7期)2017-06-05 09:36:13
    三 邊 柳
    基于PSOC超聲測(cè)距系統(tǒng)設(shè)計(jì)
    有限域上高斯正規(guī)基的一個(gè)注記
    欧美高清性xxxxhd video| 黄色女人牲交| 中出人妻视频一区二区| 久久久精品欧美日韩精品| 亚洲成人精品中文字幕电影| 婷婷丁香在线五月| 在现免费观看毛片| 精品久久久久久久末码| 日本 欧美在线| 国产精品精品国产色婷婷| 国产精品久久久久久精品电影| 中文字幕人成人乱码亚洲影| 亚洲电影在线观看av| 国产精品98久久久久久宅男小说| 亚洲中文日韩欧美视频| 窝窝影院91人妻| 精品久久久久久久久av| 亚洲av二区三区四区| 麻豆成人午夜福利视频| 久久精品国产亚洲av香蕉五月| 91九色精品人成在线观看| 免费一级毛片在线播放高清视频| 90打野战视频偷拍视频| 精品午夜福利在线看| www日本黄色视频网| 日韩精品青青久久久久久| 国产一区二区三区视频了| 亚洲国产精品合色在线| 日韩国内少妇激情av| 中文字幕人妻熟人妻熟丝袜美| 一个人免费在线观看的高清视频| 丁香欧美五月| 国产亚洲精品综合一区在线观看| 免费搜索国产男女视频| 国产野战对白在线观看| 老司机深夜福利视频在线观看| 真实男女啪啪啪动态图| 免费在线观看亚洲国产| 中国美女看黄片| 亚洲狠狠婷婷综合久久图片| 亚洲av成人不卡在线观看播放网| 国产成+人综合+亚洲专区| 特大巨黑吊av在线直播| 99国产精品一区二区蜜桃av| 亚洲精品亚洲一区二区| av专区在线播放| 99久国产av精品| 熟女人妻精品中文字幕| 中文字幕熟女人妻在线| 搡老岳熟女国产| 国产精品不卡视频一区二区 | 露出奶头的视频| 久久欧美精品欧美久久欧美| av天堂中文字幕网| 91九色精品人成在线观看| 国产aⅴ精品一区二区三区波| 免费在线观看成人毛片| 色精品久久人妻99蜜桃| 欧美日韩中文字幕国产精品一区二区三区| 成人国产一区最新在线观看| 中出人妻视频一区二区| 色精品久久人妻99蜜桃| 深夜精品福利| 亚洲精品粉嫩美女一区| 国产精品,欧美在线| 91av网一区二区| 亚洲精品在线美女| 亚洲经典国产精华液单 | 久久人人爽人人爽人人片va | 丰满乱子伦码专区| 久久精品夜夜夜夜夜久久蜜豆| 成人高潮视频无遮挡免费网站| 又爽又黄无遮挡网站| 99久国产av精品| 亚洲第一欧美日韩一区二区三区| 色综合欧美亚洲国产小说| a级一级毛片免费在线观看| 国产熟女xx| 亚洲国产日韩欧美精品在线观看| 99久久精品一区二区三区| 久久久久久久久久成人| 夜夜看夜夜爽夜夜摸| 无遮挡黄片免费观看| 九九久久精品国产亚洲av麻豆| 中文字幕熟女人妻在线| 欧美成人一区二区免费高清观看| 欧美xxxx性猛交bbbb| 性色av乱码一区二区三区2| 高清日韩中文字幕在线| 国内少妇人妻偷人精品xxx网站| 国产老妇女一区| 别揉我奶头 嗯啊视频| 最近在线观看免费完整版| 91在线精品国自产拍蜜月| 国产一区二区激情短视频| 男女下面进入的视频免费午夜| 又粗又爽又猛毛片免费看| 一级黄色大片毛片| 99久久99久久久精品蜜桃| 亚洲avbb在线观看| 国内精品久久久久久久电影| 桃红色精品国产亚洲av| 噜噜噜噜噜久久久久久91| 久久中文看片网| 3wmmmm亚洲av在线观看| 久久欧美精品欧美久久欧美| 一区二区三区激情视频| 欧美日韩中文字幕国产精品一区二区三区| 日韩有码中文字幕| 最近最新免费中文字幕在线| 最近中文字幕高清免费大全6 | 熟妇人妻久久中文字幕3abv| 欧美一级a爱片免费观看看| 校园春色视频在线观看| 热99re8久久精品国产| 色播亚洲综合网| 免费大片18禁| 国产国拍精品亚洲av在线观看| 国产精品久久久久久精品电影| 91午夜精品亚洲一区二区三区 | 白带黄色成豆腐渣| 国产免费男女视频| 国产精品乱码一区二三区的特点| 国产在视频线在精品| 婷婷丁香在线五月| xxxwww97欧美| 久99久视频精品免费| 亚洲人成伊人成综合网2020| 久9热在线精品视频| 99久久无色码亚洲精品果冻| 免费大片18禁| 变态另类丝袜制服| 午夜a级毛片| 国产高潮美女av| 国产精品98久久久久久宅男小说| 在线观看舔阴道视频| 窝窝影院91人妻| 国产欧美日韩一区二区精品| 欧美乱妇无乱码| 欧美xxxx黑人xx丫x性爽| 欧美又色又爽又黄视频| 搞女人的毛片| 日本黄色片子视频| 99在线视频只有这里精品首页| 丝袜美腿在线中文| 嫩草影院精品99| 国产黄a三级三级三级人| 欧美午夜高清在线| eeuss影院久久| 九九久久精品国产亚洲av麻豆| av天堂在线播放| 日韩欧美国产一区二区入口| 午夜福利欧美成人| 综合色av麻豆| 国产精品日韩av在线免费观看| 国产免费av片在线观看野外av| 精品久久久久久,| 亚洲一区二区三区色噜噜| 九九在线视频观看精品| 韩国av一区二区三区四区| 欧美潮喷喷水| 久久久精品大字幕| 国产69精品久久久久777片| 日韩有码中文字幕| 亚洲av免费高清在线观看| 国产亚洲精品久久久com| 午夜福利视频1000在线观看| 最近最新免费中文字幕在线| 精品久久久久久久末码| 久久久久国内视频| 一边摸一边抽搐一进一小说| 日韩高清综合在线| 18禁裸乳无遮挡免费网站照片| 嫁个100分男人电影在线观看| 夜夜爽天天搞| 精品久久久久久久末码| 俺也久久电影网| 能在线免费观看的黄片| 别揉我奶头~嗯~啊~动态视频| 能在线免费观看的黄片| 国产成人av教育| 夜夜夜夜夜久久久久| 久久草成人影院| 夜夜夜夜夜久久久久| 99久久精品国产亚洲精品| 热99re8久久精品国产| x7x7x7水蜜桃| 又爽又黄无遮挡网站| 国内久久婷婷六月综合欲色啪| 日韩欧美在线二视频| 天堂影院成人在线观看| 国产成人啪精品午夜网站| 日本黄色片子视频| 精品99又大又爽又粗少妇毛片 | 丰满人妻熟妇乱又伦精品不卡| 老鸭窝网址在线观看| 性色av乱码一区二区三区2| 婷婷色综合大香蕉| 五月玫瑰六月丁香| 韩国av一区二区三区四区| 亚洲天堂国产精品一区在线| 亚洲成a人片在线一区二区| 日韩高清综合在线| 精品久久国产蜜桃| 草草在线视频免费看| 三级国产精品欧美在线观看| 一本综合久久免费| 亚洲av不卡在线观看| 久久天躁狠狠躁夜夜2o2o| 长腿黑丝高跟| eeuss影院久久| 国产乱人视频| 欧美激情久久久久久爽电影| 国产精品久久电影中文字幕| 久99久视频精品免费| 国产一区二区亚洲精品在线观看| 18禁黄网站禁片免费观看直播| 成人特级黄色片久久久久久久| 欧美日韩国产亚洲二区| 级片在线观看| 99久久精品国产亚洲精品| 男插女下体视频免费在线播放| 在线观看美女被高潮喷水网站 | 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 男女下面进入的视频免费午夜| 中文字幕av成人在线电影| 白带黄色成豆腐渣| 欧美高清性xxxxhd video| 色视频www国产| 欧美zozozo另类| 亚洲国产欧洲综合997久久,| 国内揄拍国产精品人妻在线| 老熟妇仑乱视频hdxx| 成年版毛片免费区| 国产精品不卡视频一区二区 | 亚洲无线观看免费| 亚洲欧美日韩高清在线视频| 欧美日韩亚洲国产一区二区在线观看| 国产成人欧美在线观看| 午夜精品久久久久久毛片777| 欧美成人a在线观看| 免费观看人在逋| www.999成人在线观看| 精华霜和精华液先用哪个| 久久人人精品亚洲av| 日韩免费av在线播放| 亚洲国产高清在线一区二区三| 国产单亲对白刺激| 久久久国产成人免费| 国产亚洲精品久久久com| 亚洲av熟女| 有码 亚洲区| 一区福利在线观看| 18禁黄网站禁片免费观看直播| 国产白丝娇喘喷水9色精品| 在线播放无遮挡| 日韩大尺度精品在线看网址| 国产成人欧美在线观看| 两人在一起打扑克的视频| 真人一进一出gif抽搐免费| 亚洲av电影在线进入| 亚洲18禁久久av| 丰满人妻一区二区三区视频av| 日韩欧美 国产精品| 熟女电影av网| 网址你懂的国产日韩在线| 色5月婷婷丁香| 天堂网av新在线| 一卡2卡三卡四卡精品乱码亚洲| 在线免费观看不下载黄p国产 | 久久久久亚洲av毛片大全| 精品不卡国产一区二区三区| netflix在线观看网站| 国内精品久久久久精免费| 亚洲美女搞黄在线观看 | 91久久精品国产一区二区成人| 琪琪午夜伦伦电影理论片6080| 久久久久久国产a免费观看| 成人三级黄色视频| 村上凉子中文字幕在线| 91久久精品电影网| 国产精品乱码一区二三区的特点| 琪琪午夜伦伦电影理论片6080| 九九久久精品国产亚洲av麻豆| 最新中文字幕久久久久| 99久久久亚洲精品蜜臀av| 三级男女做爰猛烈吃奶摸视频| 波野结衣二区三区在线| 日日夜夜操网爽| 中文字幕人妻熟人妻熟丝袜美| 桃红色精品国产亚洲av| 69av精品久久久久久| 欧美黑人巨大hd| 我要看日韩黄色一级片| 我要搜黄色片| 精品免费久久久久久久清纯| 好男人在线观看高清免费视频| 麻豆av噜噜一区二区三区| 精品一区二区三区av网在线观看| 在线国产一区二区在线| 国产伦精品一区二区三区四那| 白带黄色成豆腐渣| 又粗又爽又猛毛片免费看| 丰满的人妻完整版| 欧美一区二区亚洲| 国产主播在线观看一区二区| 村上凉子中文字幕在线| 97超视频在线观看视频| 欧美性猛交╳xxx乱大交人| 国产亚洲精品久久久com| 婷婷丁香在线五月| a级一级毛片免费在线观看| 亚洲专区国产一区二区| 亚洲成人精品中文字幕电影| 高清日韩中文字幕在线| a级毛片免费高清观看在线播放| 婷婷六月久久综合丁香| 一级毛片久久久久久久久女| 亚洲欧美日韩高清专用| 国产单亲对白刺激| 波多野结衣巨乳人妻| 日本黄色视频三级网站网址| 一区福利在线观看| 亚洲av一区综合| 99久国产av精品| 国产亚洲精品综合一区在线观看| 悠悠久久av| 在线观看一区二区三区| 我要搜黄色片| 非洲黑人性xxxx精品又粗又长| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 内射极品少妇av片p| 动漫黄色视频在线观看| 国产高清激情床上av| 一级毛片久久久久久久久女| 亚洲欧美清纯卡通| 欧洲精品卡2卡3卡4卡5卡区| 欧美+日韩+精品| 欧美又色又爽又黄视频| 久久人人爽人人爽人人片va | 精品久久久久久成人av| 久久久久久久久久成人| 人人妻,人人澡人人爽秒播| 亚洲国产欧洲综合997久久,| 青草久久国产| 国产aⅴ精品一区二区三区波| 欧美+日韩+精品| 国产免费一级a男人的天堂| 国产精品免费一区二区三区在线| 久久天躁狠狠躁夜夜2o2o| 欧洲精品卡2卡3卡4卡5卡区| 91九色精品人成在线观看| 日韩免费av在线播放| 欧美黑人欧美精品刺激| 啪啪无遮挡十八禁网站| 99久国产av精品| 韩国av一区二区三区四区| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品色激情综合| 日韩av在线大香蕉| 国产精品亚洲av一区麻豆| 日日夜夜操网爽| 国内精品美女久久久久久| 亚洲激情在线av| 给我免费播放毛片高清在线观看| 老司机午夜福利在线观看视频| 国产精品不卡视频一区二区 | 日韩中文字幕欧美一区二区| 黄色女人牲交| 性色avwww在线观看| 国产成人欧美在线观看| 好看av亚洲va欧美ⅴa在| 熟女电影av网| 色噜噜av男人的天堂激情| 欧美黑人巨大hd| 婷婷丁香在线五月| a级一级毛片免费在线观看| 免费av不卡在线播放| 美女免费视频网站| 色综合站精品国产| 禁无遮挡网站| 国产av麻豆久久久久久久| 99视频精品全部免费 在线| 美女大奶头视频| 听说在线观看完整版免费高清| 好男人在线观看高清免费视频| 如何舔出高潮| 国产熟女xx| 日韩免费av在线播放| 国产精品免费一区二区三区在线| 禁无遮挡网站| 午夜精品一区二区三区免费看| 白带黄色成豆腐渣| 久久99热这里只有精品18| 亚洲av电影在线进入| 亚洲第一欧美日韩一区二区三区| 中文字幕高清在线视频| 麻豆久久精品国产亚洲av| 哪里可以看免费的av片| 18禁裸乳无遮挡免费网站照片| 老司机福利观看| 757午夜福利合集在线观看| 国产成+人综合+亚洲专区| 久久久久久九九精品二区国产| 国产三级黄色录像| 很黄的视频免费| 窝窝影院91人妻| 国产高潮美女av| 亚洲无线在线观看| 国产精品美女特级片免费视频播放器| 午夜久久久久精精品| 欧美激情久久久久久爽电影| 国产免费av片在线观看野外av| 免费观看人在逋| 精品熟女少妇八av免费久了| www.熟女人妻精品国产| 亚洲电影在线观看av| 国产野战对白在线观看| 露出奶头的视频| 欧美黄色淫秽网站| 国产亚洲精品av在线| 久久人人精品亚洲av| 亚洲欧美日韩卡通动漫| 久久草成人影院| 中文资源天堂在线| 97碰自拍视频| 久久精品综合一区二区三区| 桃色一区二区三区在线观看| 搡老熟女国产l中国老女人| 亚洲精品亚洲一区二区| 免费观看人在逋| 欧美3d第一页| 国产熟女xx| 热99在线观看视频| 欧美性感艳星| 国产精品精品国产色婷婷| 亚洲国产色片| av专区在线播放| 国产乱人伦免费视频| 在线观看66精品国产| 亚洲18禁久久av| 国产一区二区三区在线臀色熟女| 国产精品嫩草影院av在线观看 | 国产精品,欧美在线| 一本一本综合久久| 男女那种视频在线观看| 欧美区成人在线视频| 亚洲自偷自拍三级| 午夜精品在线福利| av专区在线播放| 日本免费一区二区三区高清不卡| 91九色精品人成在线观看| 两人在一起打扑克的视频| 看片在线看免费视频| 国产视频一区二区在线看| av天堂在线播放| 亚洲无线在线观看| 亚洲av不卡在线观看| 亚洲欧美日韩东京热| 国产不卡一卡二| a在线观看视频网站| 高清在线国产一区| 夜夜看夜夜爽夜夜摸| 欧美日韩中文字幕国产精品一区二区三区| 精品人妻熟女av久视频| 日日摸夜夜添夜夜添小说| 丁香欧美五月| 一个人观看的视频www高清免费观看| 日韩欧美在线二视频| 色哟哟哟哟哟哟| 十八禁网站免费在线| 日韩有码中文字幕| 国内精品一区二区在线观看| 久久这里只有精品中国| 亚洲成av人片免费观看| 亚洲成av人片在线播放无| 人人妻人人澡欧美一区二区| 精品午夜福利在线看| 国产美女午夜福利| 99久久99久久久精品蜜桃| 村上凉子中文字幕在线| 99riav亚洲国产免费| 国产色婷婷99| 国产精品久久久久久久电影| 国产精品精品国产色婷婷| 免费看日本二区| 国产精品1区2区在线观看.| 日韩欧美精品v在线| 国产淫片久久久久久久久 | 免费av观看视频| 久久国产乱子伦精品免费另类| 久久九九热精品免费| 亚洲av成人精品一区久久| 亚洲成人久久性| 亚洲性夜色夜夜综合| 高清在线国产一区| 我的老师免费观看完整版| 狠狠狠狠99中文字幕| 乱码一卡2卡4卡精品| av中文乱码字幕在线| 精品久久国产蜜桃| 亚洲不卡免费看| 日韩精品青青久久久久久| 我要搜黄色片| 日韩大尺度精品在线看网址| 大型黄色视频在线免费观看| 亚洲性夜色夜夜综合| 亚洲av不卡在线观看| 女同久久另类99精品国产91| 深夜精品福利| 国产精品美女特级片免费视频播放器| 亚洲午夜理论影院| 91av网一区二区| 最近最新免费中文字幕在线| 美女被艹到高潮喷水动态| 久久精品久久久久久噜噜老黄 | 中国美女看黄片| 国产大屁股一区二区在线视频| 日本免费一区二区三区高清不卡| 我要搜黄色片| 精品久久久久久久久久久久久| 亚洲中文字幕日韩| 少妇的逼好多水| 国产成人欧美在线观看| 99视频精品全部免费 在线| 亚洲无线在线观看| 国产久久久一区二区三区| 不卡一级毛片| 搡老妇女老女人老熟妇| 国产av一区在线观看免费| 成人毛片a级毛片在线播放| 最近中文字幕高清免费大全6 | 757午夜福利合集在线观看| 亚洲美女视频黄频| 久久天躁狠狠躁夜夜2o2o| 国产激情偷乱视频一区二区| 婷婷色综合大香蕉| 好男人在线观看高清免费视频| 99久久精品一区二区三区| 亚洲精品成人久久久久久| 欧美成狂野欧美在线观看| 在线观看舔阴道视频| 丰满的人妻完整版| 别揉我奶头~嗯~啊~动态视频| 18+在线观看网站| 国产午夜精品久久久久久一区二区三区 | 国产私拍福利视频在线观看| 国产精品嫩草影院av在线观看 | 欧美日韩国产亚洲二区| 日韩欧美国产一区二区入口| 色精品久久人妻99蜜桃| 美女被艹到高潮喷水动态| 亚洲国产精品sss在线观看| 禁无遮挡网站| 欧美+日韩+精品| 我的女老师完整版在线观看| 99久久精品国产亚洲精品| 欧美日韩国产亚洲二区| 高清日韩中文字幕在线| aaaaa片日本免费| 91狼人影院| 宅男免费午夜| 免费看日本二区| 人妻夜夜爽99麻豆av| 搡老妇女老女人老熟妇| 日本在线视频免费播放| 久久欧美精品欧美久久欧美| 国模一区二区三区四区视频| 好看av亚洲va欧美ⅴa在| 高清日韩中文字幕在线| 国产69精品久久久久777片| 一区二区三区激情视频| 超碰av人人做人人爽久久| 亚洲在线观看片| 美女免费视频网站| 日韩 亚洲 欧美在线| 国产精品1区2区在线观看.| 综合色av麻豆| 国产在视频线在精品| 此物有八面人人有两片| 午夜福利在线观看吧| 天堂av国产一区二区熟女人妻| 日韩中字成人| 中文字幕精品亚洲无线码一区| 九九在线视频观看精品| 深爱激情五月婷婷| 久久人妻av系列| 午夜a级毛片| 亚洲精品456在线播放app | 婷婷丁香在线五月| 国产爱豆传媒在线观看| 午夜老司机福利剧场| 一个人看视频在线观看www免费| 一进一出好大好爽视频| 99热这里只有是精品在线观看 | .国产精品久久| 国产高清视频在线观看网站| 欧美丝袜亚洲另类 | 嫩草影院新地址| 高清在线国产一区| 日韩欧美国产一区二区入口| 给我免费播放毛片高清在线观看| 亚洲最大成人中文| 亚洲精品粉嫩美女一区| 欧美日韩亚洲国产一区二区在线观看| 午夜福利在线观看免费完整高清在 | 日韩人妻高清精品专区| 偷拍熟女少妇极品色| 精品国内亚洲2022精品成人| 永久网站在线| 亚洲黑人精品在线| 少妇熟女aⅴ在线视频| 日韩欧美免费精品| 欧美高清成人免费视频www| 九九久久精品国产亚洲av麻豆|