亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        面向大中型機房的RSSI定位算法研究

        2017-10-24 11:31:45林杰華樂志偉
        關(guān)鍵詞:信號

        林 航, 李 葵, 林杰華, 樂志偉, 王 浩, 周 立

        (1.國網(wǎng)安徽省電力公司 信息通信分公司, 安徽 合肥 230061; 2.合肥工業(yè)大學(xué) 計算機與信息學(xué)院, 安徽 合肥 230009)

        面向大中型機房的RSSI定位算法研究

        林 航1, 李 葵1, 林杰華2, 樂志偉2, 王 浩2, 周 立2

        (1.國網(wǎng)安徽省電力公司 信息通信分公司, 安徽 合肥 230061; 2.合肥工業(yè)大學(xué) 計算機與信息學(xué)院, 安徽 合肥 230009)

        針對大中型機房中信號閱讀器只能部署在較高安全位置,不能直接采用基于測距的接收的信號強度指標(biāo)(

        signal strength indication,RSSI)定位算法進行三維定位的問題,文章通過分析機房內(nèi)特有的室內(nèi)環(huán)境,提出了在原定位算法計算結(jié)果的基礎(chǔ)上,采用逼近定位算法的方案。考慮障礙物對信號傳播的遮擋,修正測距算法中閱讀器到有源標(biāo)簽的距離矢量;引入?yún)⒖脊?jié)點方向矢量權(quán)重,改善加權(quán)質(zhì)心法,迭代矯正定位,并對異常結(jié)果進行矯正處理。相比基于RSSI原定的定位算法,其結(jié)果不僅給出了三維坐標(biāo),在平面上的定位精度也提高了約60%。

        接收的信號強度指示(RSSI)定位;三維定位;有源標(biāo)簽;障礙物;加權(quán)質(zhì)心法;定位精度

        大中型機房內(nèi)部具有人員流動較小、機柜形狀規(guī)則、高度相同、位置固定且噪音較大等環(huán)境特點。出于安全要求,閱讀器通常部署在機柜上方,以防止人員觸碰。閱讀器均部署在同一高度的平面上,因此如果采用文獻(xiàn)[1-10]提出的定位算法,只能進行平面定位,本文在文獻(xiàn)[4]基礎(chǔ)上進一步考慮機房障礙物的部署特點及其對信號傳播的影響,修正信號傳播的理論衰減距離,并通過改進文獻(xiàn)[5]質(zhì)心法的加權(quán)因子重新定位,通過迭代算法,使得定位結(jié)果逐漸逼近有源標(biāo)簽真實的三維坐標(biāo)。其定位結(jié)果在平面上的精度也提高了60%。

        1 基于機柜障礙物的信號衰減模型

        1.1 信號傳播的理論衰減模型

        無線射頻信號能量傳播隨著距離的增加呈對數(shù)衰減,其理想信號衰減模型如下:

        (1)

        其中,PL(d)為在距離有源標(biāo)簽的d處信號的衰減量;PL(d0)為參考距離d0處的信號衰減量;X為標(biāo)準(zhǔn)方差;n為信號衰減系數(shù),該系數(shù)因環(huán)境而異。由該信號量衰減模型可得距離信號源d處的信號強度Pr(d)為:

        (2)

        其中,P0為信號源的額定功率強度;Pr(d)為距離標(biāo)簽d處閱讀器的信號強度接收值。令d0為1,X為0,通過接收的Pr(d)可求得衰減距離d,即

        PLA=P0-PL(d0),

        Pr(d)=PLA-10nlgd-X

        (3)

        1.2 基于路徑分段的信號傳播衰減模型

        機房存在大量整齊固定的機柜箱障礙物,而信號傳播穿過障礙物時會發(fā)生反射、折射及透射等現(xiàn)象,給基于測距的定位算法帶來誤差。信號在不同障礙物內(nèi)傳播時,能量的衰減系數(shù)不同,針對障礙物影響基于RSSI測距準(zhǔn)確性的問題,本文提出路徑分段衰減模型,修正各閱讀器與標(biāo)簽對應(yīng)的距離矢量。根據(jù)文獻(xiàn)[4]定位結(jié)果坐標(biāo)(x,y,z),建立閱讀器Rj到其的路徑距離矢量。信號源距離對應(yīng)Rj的路徑距離矢量進行分段處理,信號衰減路徑分段如圖1所示。

        圖1 信號傳播路徑分段衰減

        圖1中,P(i-1)為信號穿過障礙物前的強度;di-1為信號在障礙物內(nèi)部的傳播路徑;P(i)為傳出障礙物界面的信號強度。由于信號強度是對數(shù)衰減的,則有:

        (4)

        其中,ni-1為信號在障礙物中的損耗系數(shù);l為信號源功率在障礙物中傳播衰減到P(i)時對應(yīng)的理論距離;l-di-1為信號源功率在障礙物中傳播衰減到P(i-1)時對應(yīng)的理論路徑長度。通過衰減模型求得l,當(dāng)l

        (5)

        1.3 障礙節(jié)點信號衰減系數(shù)的可靠性分析

        本文采用網(wǎng)格采點方式采集信息,在障礙物表面建立虛擬網(wǎng)格采集點,網(wǎng)格的邊長為L(實驗長度定為30 cm),每個采集點與每個閱讀器都有一條固定的物理路徑,并包含信號在該條路徑上貫穿障礙物時的衰減系數(shù),這些采集點的信息為經(jīng)驗數(shù)據(jù)。

        通過文獻(xiàn)[4]所提的算法定位結(jié)果及幾何原理可求得信號傳播到閱讀器Rj的路徑穿過障礙物表面時的障礙節(jié)點p1點,在距離p1點的L范圍內(nèi)選取K個距離最近的采集點作為信息參考點集,以確定在該障礙物內(nèi),信號在該節(jié)點與對應(yīng)閱讀器路徑上的衰減系數(shù)。

        選取距離待測點最近的K個信號采集點為TKi,i∈(1,K)。根據(jù)信號傳播理論衰減模型,如果采用其中a個采集點中與閱讀器Rj對應(yīng)衰減系數(shù)和信號接收強度計算出的衰減距離l>di-1,說明信號在這些點對應(yīng)的路徑上可以貫穿障礙物。

        (1) 當(dāng)a≠0時,這些點的可信度更高(保留了上一次定位的準(zhǔn)確性),設(shè)p2點為這些點的中心,則信號源在p2點和Rj坐標(biāo)的連接直線上,而其對應(yīng)的衰減系數(shù)為這些采集點對應(yīng)衰減系數(shù)的平均數(shù),相對于閱讀器Rj來說,信號源的方向由Rj到p1,變成了由Rj到p2,從而修正閱讀器到標(biāo)簽的方向Vj,其對應(yīng)的角度修正偏移量為θj,這個變量是本文核心算法中的一個重要參數(shù)。信號在p2到Rj連接直線上穿過障礙物時的衰減系數(shù)為這些點的平均值。

        (6)

        (7)

        求得Rj到p1的距離d1;Rj到p2的距離d2;p2到p1的距離d3。由余弦定理可知θj為:

        (8)

        (2) 當(dāng)a=0時, 則θj=0;求得Tdi為TKi與p1點的距離,距離越近參考衰減系數(shù)ni的權(quán)值wi越大,本文權(quán)重因子公式如下:

        (9)

        目標(biāo)點對應(yīng)的信號衰減系數(shù)n為:

        (10)

        2 逼近定位

        2.1LCV原定位算法

        出于安全考慮,機房內(nèi)接收信號的閱讀器部署在同一平面,通過其接收信號強度計算出理論距離,因此可采用文獻(xiàn)[4]提出的LCV算法,即基于加權(quán)質(zhì)心法進行定位,并作為本文的原定位算法,求出平面坐標(biāo),定位公式如下:

        (11)

        其中,A、B、C為信號接收最強的3個參考節(jié)點; (dA,dB,dC)為相應(yīng)的理論衰減距離。該算法加強了信號強度較強參考節(jié)點的權(quán)重,弱化了信號能量較弱節(jié)點的參考價值。

        2.2 改進加權(quán)質(zhì)心法逼近定位

        由于上述定位算法無法給出目標(biāo)Z軸坐標(biāo),本文考慮了機房內(nèi)部的環(huán)境特征,采用逼近矯正定位算法的思路,即在LCV算法結(jié)果的基礎(chǔ)上選定Z軸上一個初始值,并以此為迭代初始參考點,進行不斷矯正。通過信號衰減理論模型求得各個閱讀器到信號源的理論距離dj,dj∈(d1,d2,d3,…,dn),而根據(jù)LCV算法結(jié)果求得其與各個閱讀器在X、Y軸平面投影的距離vdj∈(vd1, vd2, vd3,…,vdn), 機柜上閱讀器高度為H,根據(jù)勾股定理即可求得每個閱讀器到信號源的高度zj∈(z1,z2,z3,…,zn),即

        (12)

        對于z0的選取,以信號接收強度作為確定Z軸初始坐標(biāo)zj的權(quán)重依據(jù),Rj接收的信號越強,dj越小,zj權(quán)重越大,即

        (13)

        得出LCV算法結(jié)果確定為 (x0,y0,z0),閱讀器Rj到定位結(jié)果的方向矢量Vj(x0-xj,y0-yj,z0-zj),由閱讀器Rj信號接收強度,采用信號傳播路徑分段衰減求距,修正Dj并修正角度偏移量θj與Vj,求得在閱讀器Rj的Vj方向上,距離Dj處的坐標(biāo)VRj,如圖2所示。

        圖2 信號傳播路徑分段衰減

        (14)

        對于閱讀器Rj來說,VRj便是其對應(yīng)的有源標(biāo)簽的坐標(biāo),新的定位矯正結(jié)果轉(zhuǎn)化為對VRj質(zhì)心的求解,再以所求質(zhì)心結(jié)果為新的目標(biāo)點進行迭代矯正,逐漸逼近目標(biāo)的真實坐標(biāo)。

        文獻(xiàn)[1-10]研究的基于距離質(zhì)心加權(quán)算法均強調(diào)以距離作為分析參考節(jié)點權(quán)重的因素,即閱讀器Rj信號接收強度越高,到標(biāo)簽距離Dj越小,VRj的坐標(biāo)加權(quán)因子越大。由于通過路徑分段衰減模型修正了Dj,減少了理論距離與實際距離的誤差,信號較弱的VRj對應(yīng)的Dj反而可能更小。 LCV算法已不適用求VRj的質(zhì)心,本文提出新的加權(quán)質(zhì)心法:首先引入閱讀器Rj接收信號強度作為參考權(quán)重,即信號Rj接收強度越大,VRj坐標(biāo)參考權(quán)重越大。

        因為VRj本身就是帶有方向性的,θj變化越大,說明其修正后對定位結(jié)果的影響越大,本文引用Vj的修正偏移量θj作為參考權(quán)重,考慮到θj可能為0的情況,將其以線性方程形式轉(zhuǎn)換成φj,即

        φj=βθj+1

        (15)

        其中,β為常數(shù)項,調(diào)節(jié)整體權(quán)重,本文實驗中θj以弧度制計算,β取1。

        各閱讀器接收待測信號源能量為(pr1,pr2,pr3,…,prn),根據(jù)文獻(xiàn)[5]所提方案,選取VRj,j∈(1,n)中任意4個節(jié)點組成若干個四方體, (r1,r2,r3,r4)為所選取的4個點信號能量的映射,都是正數(shù),每個四方體進行加權(quán)質(zhì)心定位,即

        (16)

        質(zhì)心處的角度偏移量φ和信號強度r為:

        (17)

        (18)

        再以新生成的點組成新的若干個四方體,直至最后只剩小于4個的節(jié)點,進行加權(quán)平均求值。

        3 Knn異常定位結(jié)果矯正

        由于迭代矯正后的結(jié)果可能分布在機房中機柜箱的陰影區(qū)域,所求結(jié)果坐標(biāo)需要進一步矯正。結(jié)合文獻(xiàn)[11-14]分析關(guān)于設(shè)定虛擬標(biāo)簽算法的思想,本文提出在障礙物表面設(shè)置虛擬參考標(biāo)簽的觀點,記錄各個閱讀器信號的接收強度,每個標(biāo)簽點包含一組信號強度的接收量表示其特征,即VPr(vpr1,vpr2,vpr3, …,vprn),vprj表示Rj的信號接收強度經(jīng)驗值。

        根據(jù)異常定位結(jié)果,選取空間幾何距離該點最近的一個表面PN作為參考標(biāo)簽選取區(qū)域,求得異常結(jié)果點在該表面的投影坐標(biāo)PN-P。并在該平面上選擇距離PN-P最近的選擇范圍Dist-thr(實驗選取對應(yīng)障礙物的寬度值)內(nèi)的所有虛擬標(biāo)簽作為參考,得到選定的虛擬參考標(biāo)簽集合VPrl,l∈(1,2,3,4, …,n)。

        出現(xiàn)異常定位結(jié)果的情況時,根據(jù)各個閱讀器接收待測信號源的實際能量TPr(pr1,pr2,pr3, …,prn)選取其中m個最強的信號值作為可信度較高的目標(biāo)計算矢量, TPr(pra,prb,…,prc),其中a,b,c∈(1,n)。

        從每個選定的虛擬參考標(biāo)簽信息中選取相應(yīng)的信號能量TPr(pra,prb,…,prc)作為參考,計算每個虛擬標(biāo)簽VPrl到目標(biāo)實際TPr的歐氏距離δl,歐式距離越小,虛擬標(biāo)簽的參考價值越大。

        設(shè)ζi=pri-vpri,引入加權(quán)歐式距離因子λi,本文采用的歐式距離δl計算公式如下:

        (19)

        (20)

        (21)

        (19)式與(20)式引入信號理論衰減距離,根據(jù)對應(yīng)距離的差異,賦予了pri-vpri不同的權(quán)重,強調(diào)對比信號越強,對應(yīng)的參考價值越大。由表達(dá)式可知,信號能量偏差量相等的情況下,距離越遠(yuǎn),信號越弱,歸一化后的λi影響因子越大。通過計算選取K個歐式距離最小的虛擬參考標(biāo)簽,繼而計算出矯正后的坐標(biāo)位置。虛擬標(biāo)簽信息參考的歐氏距離越小,參考因數(shù)越大,選取其中歐式距離最小的K個虛擬標(biāo)簽作為參考標(biāo)簽。設(shè)選定的K個虛擬標(biāo)簽中,虛擬坐標(biāo)posu、wu為對應(yīng)的歐式距離參考因子,糾正后的目標(biāo)坐標(biāo)為correct-pos,即

        (22)

        (23)

        其中,posu為選定的虛擬標(biāo)簽對應(yīng)的坐標(biāo);wu為對應(yīng)的權(quán)重。

        4 定位流程

        本文定位算法思想是結(jié)合機房內(nèi)部機柜分布及閱讀器部署的環(huán)境特征,矯正基于測距的加權(quán)質(zhì)心算法,并將質(zhì)心算法矢量化,逐步逼近目標(biāo)真實坐標(biāo)值,定位流程如下:

        (1) LCV原定位算法定位,通過勾股定理及質(zhì)心算法,在Z軸上選擇一個適當(dāng)?shù)某跏贾?得到定位結(jié)果Tagi其中i=0。

        (2) 求閱讀器Rj到Tagi的方向Vj,結(jié)合室內(nèi)障礙物坐標(biāo)及其信號衰減信息,基于路徑衰減模型修正Dj與Vj,并得到θj。

        (3) 根據(jù)修正后的Dj與Vj,推算出每個閱讀器認(rèn)為對應(yīng)的有源標(biāo)簽坐標(biāo)VRj,j∈(1,2,3,4, …,n)。

        (4) 引入VRj對應(yīng)的角度矯正因子θj,加權(quán)定位,得出Tagi+1。

        (5) 當(dāng)i

        5 實驗效果

        本文實驗機房大小為16 m×16 m×4 m,如圖3所示。

        圖3 機房俯視圖

        經(jīng)測算信號源在機房內(nèi)傳播1 m 處的信號接收強度PLA為-51 dbm,信號能量衰減系數(shù)PLN為1.9。本實驗根據(jù)閱讀器與目標(biāo)源的位置中的障礙物遮擋數(shù)量,分為3種類型,即障礙物為1、2、3個的3種情況,每種情況隨機采集20個點, 計算用基于路徑分段衰減模型測距算法所得結(jié)果,并與不考慮障礙物衰減模型算法所求結(jié)果作比較,前者相比于后者減少的距離誤差量如圖4所示。通過對比分析,在標(biāo)簽與閱讀器障礙物比較少的情況下,采用路徑分段求距,修正的誤差不是很明顯,但隨著障礙物多時,誤差的減少比較明顯,修正后的距離更可靠,尤其是在閱讀器接收信號較弱的情況下,誤差減小比較明顯。

        圖4 路徑分段求距

        為了分析本文所提算法對LCV原定算法定位精度的影響,本文選取均勻分布在機房內(nèi)較有代表性的50個測試點,每個測試點選取5個不同高度的相應(yīng)坐標(biāo)放置有源標(biāo)簽,求得總體誤差的平均值。并通過添加閱讀器的個數(shù),分析2種算法結(jié)果在X、Y軸平面上的定位誤差,其中本文算法迭代1次,數(shù)據(jù)結(jié)果如圖5所示。

        圖5 X、Y軸平面定位誤差

        不難看出,多數(shù)情況下,本文所提算法對原定算法的重定位后,誤差減小比較明顯,但是當(dāng)閱讀器大于13個時,矯正重定位效果不明顯,考慮到閱讀器部署上的性價比,最終選擇閱讀器的數(shù)量為10。為了進一步測試算法迭代的有效性,以下是在機房內(nèi)隨機選取的10個待測點,閱讀器個數(shù)為10,分別進行10次迭代運算法,運行后的結(jié)果在X、Y、Z軸上立體定位誤差如圖6所示。

        圖6 算法迭代誤差分析

        從圖6中可以看出,通過迭代算法,定位誤差得到了進一步的減少,并且大致在迭代到第3次以后,誤差范圍在一定范圍內(nèi)波動,總體上誤差減小的幅度不大。當(dāng)開始誤差比較小的時候,通過迭代算法,減小誤差效果不明顯。考慮到系統(tǒng)的開銷和定位的實時性,室內(nèi)機房定位的算法迭代次數(shù)N定為3,而對應(yīng)的誤差進一步減小,相比于原定算法而言,給出了在Z軸上0.45 m誤差范圍內(nèi)的坐標(biāo),而在平面其定位精度提高了約60%。

        6 結(jié) 論

        本文針對大中型機房環(huán)境特點,修正基于測距定位算法的距離矢量,并采用逼近矯正定位算法的新思路,不僅在平面上較原定位算法(LCV)精度提高了約60%,而且給出了目標(biāo)的三維坐標(biāo),但本文未涉及對部署的閱讀器拓?fù)浣Y(jié)構(gòu)的研究。

        [1] 陳順明,李平.基于RSSI權(quán)值的環(huán)境適應(yīng)型室內(nèi)定位算法研究[J].計算機工程與應(yīng)用,2015,51(22):99-103.

        [2] 葛斌,鄭建寶,韓江洪.RSSI輔助的三維空間坐標(biāo)四面體質(zhì)心定位算法[J].計算機科學(xué),2015,42(4):81-83.

        [3] 吳杰,馮鋒,丁志義.基于RFID與WSN融合技術(shù)的井下定位算法研究[J].計算機工程與設(shè)計,2014,35(7):2280-2282.

        [4] 王振朝,張琦,張峰.基于RSSI測距的改進加權(quán)質(zhì)心定位算法[J].電測與儀表,2014,51(21):63-65.

        [5] 陳均,李穩(wěn)操.基于RSSI的三維加權(quán)質(zhì)心定位研究[J].軟件,2014,35(3):76-79.

        [6] 尚少鋒,張雪英,王峰.一種基于RSSI矯正的加權(quán)質(zhì)心定位算法[J].科學(xué)技術(shù)與工程,2013,13(10):2833-2835.

        [7] 路標(biāo).質(zhì)心定位中的加權(quán)坐標(biāo)誤差修正方法[J].電視技術(shù),2013,37(23):147-149.

        [8] ZHANG W G,ZHANG G H,CAO J,et al.The wireless sensor network node localization algorithm based on multitudinous beacon[C]//6th International Conference on Electronics and Information Engineering.Dalian,China:ICEIE,2015:97940M.

        [9] KNAUTH S, STORZ M, DASTAGEERI H,et al.Fingerprint calibrated centroid and scalar product correlation RSSI positioning in large environments[C]//International Conference on Indoor Positioning and Indoor Navigation.Banff, AB, Canada:IPIN,2015:1-6.

        [10] LIU Yong, ZHANG Jinlong, ZHANG Yanbo,et al.Improved weighted centroid localization algorithm in narrow space[J].Journal of Computer Applications,2015,35(5):1274-1276.

        [11] 曾碧,毛勤.改進的室內(nèi)三維模糊位置指紋定位算法[J].山東大學(xué)學(xué)報(工學(xué)版),2015,45(3):23-26.

        [12] 董永峰,周艷聰,孫陸楠.基于RFID的虛擬參考標(biāo)簽三維室內(nèi)定位算法[J].計算機工程與設(shè)計,2015,36(6):1536-1538.

        [13] 李衛(wèi)平,楊杰,王鋼.比例逆權(quán)重kNN算法及其流處理應(yīng)用[J].計算機工程與設(shè)計,2015,36(12): 3556-3557.

        [14] 潘登.KNN算法的相近似度研究[D].沈陽:東北師范大學(xué),2014:12-15.

        ResearchonRSSIlocationalgorithmforlargeandmediumgeneratorroom

        LIN Hang1, LI Kui1, LIN Jiehua2, LE Zhiwei2, WANG Hao2, ZHOU Li2

        (1.Information and Telecommunication Branch, State Grid Anhui Electric Power Company, Hefei 230061, China; 2.School of Computer and Information, Hefei University of Technology, Hefei 230009, China)

        As the signal reader in large and medium generator room can only be deployed in high security positions, three-dimensional positioning can not be made directly by the received signal strength indication(RSSI) location based on the ranging algorithm. A localization algorithm using approximation scheme is proposed on the basis of the results of original positioning method by considering the unique features of generator room. The distance vector of each signal reader to the active tag in the ranging algorithm is corrected by considering the blocking of obstacles to the signal propagation. The improved weighted centroid algorithm is presented by introducing the weights of the reference node directional vector, the approach of iterative correction is adopted and the correction processing about abnormal results is made. The presented method not only gives the three-dimensional coordinate but also improves the accuracy of planar positioning by nearly 60% compared to the original method of RSSI location.

        received signal strength indication(RSSI) location; three-dimensional positioning; active tag; obstacle; weighted centroid algorithm; positioning accuracy

        signal strength indication,RSSI)定位技術(shù)由于其低成本、底功耗等特點得到廣泛運用和關(guān)注,其由發(fā)射無線信號的有源標(biāo)簽以及可以接收識別信號的閱讀器組成。文獻(xiàn)[1-10]研究并提出了基于RSSI測距定位的各類算法,但由于無線信號在復(fù)雜的室內(nèi)環(huán)境中會發(fā)生反射、衍射、投射等現(xiàn)象,特別是大中型機房由于障礙物較多,這些算法運用于大中型機房定位時,定位精度不穩(wěn)定。

        2016-03-23;

        2016-10-14

        國家自然科學(xué)基金資助項目(61175033)

        林 航(1975-),男,江蘇海安人,國網(wǎng)安徽省電力公司高級工程師;

        王 浩(1962-),男,江蘇泰州人,博士,合肥工業(yè)大學(xué)教授,博士生導(dǎo)師;

        林杰華(1965-),男,安徽望江人,合肥工業(yè)大學(xué)講師,通訊作者,E-mail:gis@hfut.edu.cn.

        10.3969/j.issn.1003-5060.2017.09.009

        TP399

        A

        1003-5060(2017)09-1193-06

        (責(zé)任編輯 張 镅)

        猜你喜歡
        信號
        信號
        鴨綠江(2021年35期)2021-04-19 12:24:18
        完形填空二則
        7個信號,警惕寶寶要感冒
        媽媽寶寶(2019年10期)2019-10-26 02:45:34
        孩子停止長個的信號
        《鐵道通信信號》訂閱單
        基于FPGA的多功能信號發(fā)生器的設(shè)計
        電子制作(2018年11期)2018-08-04 03:25:42
        基于Arduino的聯(lián)鎖信號控制接口研究
        《鐵道通信信號》訂閱單
        基于LabVIEW的力加載信號采集與PID控制
        Kisspeptin/GPR54信號通路促使性早熟形成的作用觀察
        91精品国产高清久久福利| 亚洲区福利视频免费看| 精品丝袜一区二区三区性色| 亚洲精品国产二区在线观看| 国产精品久久av色婷婷网站| 一本色道久久综合狠狠躁中文| 精品午夜中文字幕熟女| 手机看片久久第一人妻| 精品福利一区二区三区蜜桃| 久久精品欧美日韩精品| 国产成人亚洲精品无码h在线| 国产视频毛片| 亚洲欧美另类日本久久影院| 日本一区二区精品色超碰| 国产精品毛片极品久久| 国产大屁股视频免费区| 日本最大色倩网站www| 精品少妇人妻av免费久久久| 亚洲AV无码专区国产H小说| 成年女人18毛片毛片免费| 在线观看亚洲视频一区二区| 不卡一区二区黄色av| 国产麻传媒精品国产av| 国产精品爽爽va在线观看无码| 亚洲成在人线电影天堂色| 国产成人夜色在线视频观看| 国产精品成人av大片| 天天躁日日躁aaaaxxxx| 怡红院免费的全部视频| 国产亚洲欧美成人久久片| 欧美人与物videos另类| 熟女高潮av一区二区| 免费国产自拍在线观看| 国产精品人妻一码二码| 久久亚洲精品ab无码播放| 日韩女人毛片在线播放| 白色白在线观看免费2| 成年性生交大片免费看| 成人看片黄a免费看那个网址| 国产精品无码专区综合网| 美女一区二区三区在线观看视频|