《電子技術應用》
您所在的位置:首頁 > 通信與網絡 > 設計應用 > 基于改進的粒子群算法在WSN節點定位中的研究
基于改進的粒子群算法在WSN節點定位中的研究
2016年微型機與應用第24期
陳玲君
紹興職業技術學院,浙江 紹興 312000
摘要: 如何能夠減小無線傳感中的節點定位誤差一直都是研究的熱點。提出一種基于改進的粒子群優化算法以修正DVHop誤差的傳感器節點定位方法,通過分析粒子間距離、雙變異因子和權重設置改進了粒子群算法,改進后的粒子群算法減少了未知節點與錨節點間距離的估計誤差。仿真實驗表明,相對于DVHOP算法,本文的算法可以有效地提高傳感器節點定位精度。
Abstract:
Key words :

  陳玲君

  (紹興職業技術學院,浙江 紹興 312000)

       摘要:如何能夠減小無線傳感中的節點定位誤差一直都是研究的熱點。提出一種基于改進的粒子群優化算法以修正DVHop誤差的傳感器節點定位方法,通過分析粒子間距離、雙變異因子和權重設置改進了粒子群算法,改進后的粒子群算法減少了未知節點與錨節點間距離的估計誤差。仿真實驗表明,相對于DVHOP算法,本文的算法可以有效地提高傳感器節點定位精度。

  關鍵詞DV-Hop算法;定位精度;估計誤差

  中圖分類號:TP393文獻標識碼:ADOI: 10.19358/j.issn.1674-7720.2016.24.020

  引用格式:陳玲君. 基于改進的粒子群算法在WSN節點定位中的研究[J].微型機與應用,2016,35(24):70-72,76.

0引言

  如何能夠在無線傳感網中進行節點定位一直都是研究的熱點,目前大多數的研究都是基于通過錨節點來計算相關未知節點的定位研究[12]。國內外學者對此進行了不斷深入的研究。DVHop是一種應用最廣泛的免測距方法,針對其定位精度較低的問題,許多學者對其進行了改進。文獻[34]利用改進的粒子群優化算法對DVHop算法的位置估計進行校正,并利用該算法進行求解。文獻[5]從兩個方面提出改進:一是基于節點的通信半徑對節點間的跳數進行修正;二是借助信標節點間的估計距離與實際距離的偏差對平均每跳的跳距進行修正,仿真實驗取得了比較好的效果。文獻[67]提出將加權運用到無線傳感網的節點定位中,取得了比較好的效果。文獻[89]提出采用其他的人工智能算法求解無線傳感網中的節點定位算法,取得了一定的效果。

  為了能夠更好地提高節點定位的效果,本文首先描述了基本的定位算法,其次對粒子群算法進行了改進,分析了節點定位誤差的原因,采用改進的粒子群算法對DVHop算法定位誤差進行修正,最后通過仿真說明本文改進算法取得的效果。

1基本算法描述

  1.1DV-Hop算法

  DVHop節點定位算法只需要包含少量的錨節點,通過定位算法來確定未知節點的位置,它具有不需要通過具體測量距離、硬件要求低等優點,特別適合在硬件條件有限的無線傳感網中的應用,具體的算法步驟如下:

  (1)錨節點在通信網絡范圍內向鄰居節點發送自己的位置信息。接收節點記錄該節點到每個錨節點之間的最小跳數,同時刪除同一個錨節點發送的最大跳數信息,然后跳數值加1,繼續轉發個下一個鄰居節點。

  (2)每一個錨節點會根據其他錨節點的坐標信息和跳數估算網絡平均跳距距離,采用公式(1)表示。

  U26G$GE]7]I8CICZG$EG36U.png

  式中,hopSij表示錨節點i和j之間的跳數。

  錨節點將平均跳距發送到整個網絡后,未知節點僅記錄所收到的第一個平均跳距,通過公式(2)估算未知節點與錨節點的距離:

  Li=Si×HopSize(2)

  (3)設P1(x1,y1),P2(x2,y2),…,Pn(xn,yn)表示n個錨節點的坐標位置,待定位節點D的位置為(x,y),其與錨節點的估計距離分別為d1,d2,…,dn,可以建立式(3)所示的方程。

  G~N$3]XXO3ZU3_NBIP7${SY.png

  在實際的測距過程中,會產生一些不確定誤差,因此線性方程組應該將誤差考慮進去,變成AL+N=b,N為隨機誤差向量,通過最小二乘法得到的解為:

  L=(ATA)-1ATb(5)

  1.2粒子群算法

  粒子群算法是一種模擬鳥群在覓食過程中的遷徙的一種群集行為的算法,假設在d維的空間搜索,粒子i的速度和位置分別為vi=(vi1,vi2,…,vid)和xi=(xi1,xi2,…,xid)粒子i自身的歷史最優位置為pi=(pi1,pi2,…,pid),種群歷史的最優位置為pg=(pg1,pg2,…,pgd),其速度和位置更新方式為:

  vk+1id=ωkid+c1r1(pid-xk+1id)+c2r2(pgd-xkid)(6)

  xk+1id=xkid+vk+1id(7)

  式中,c1和c2分別為加速系數,r1和r2為[0,1]之間的隨機數,k為迭代次數,ω為慣性權值。

  粒子群算法具有算法結構簡單、容易實現等優點,缺點是算法粒子容易發生早熟,易陷入局部最優。為了能夠更好地實現粒子群算法在DVHop中的定位,本文對粒子群算法進行優化,使得優化后的粒子群算法能夠解決粒子群算法陷入局部最優的缺陷,為節點定位提供理論支持。

2改進的粒子群算法

  2.1粒子之間的平均距離

  根據文獻[10]研究發現,粒子群算法的早熟收斂問題與種群的多樣性密切相關,粒子群算法在后期的執行過程中,整個過種群多樣性會逐步收斂,種群的多樣性也在不斷地喪失。本文采用粒子間的平均距離來衡量種群的多樣性,如下:

  FUXLUBLO2_@]0~S}[%S_W]I.png

  式中,m表示種群中的粒子的數目,D表示空間的維數,Lmax表示粒子群空間的半徑,pij表示第i個粒子在第j維空間變量,pj是所有粒子在第j維空間變量中的平均值。

  2.2雙變異因子

  本文基于對文獻[11]的研究,結合遺傳算法的變異概念,提出了雙變異因子來對目標粒子進行跟隨,雙變異因子由最優因子Ybest和惰性因子Yworst組成,前者主要用來跟隨每次迭代中的適應度函數值最優的粒子,后者主要是用來跟隨每次迭代中適應度函數值最差的粒子。當進行迭代的時候對兩個因子進行高速變異,對最優因子采用公式(9)進行變異,這樣有利于在局部最優解的附近的范圍內提高搜索精度,對惰性因子跟蹤的粒子按照公式(10)進行變異,這樣有利于擴大種群的范圍,持續更新種群的“生命力”。

  Yworsti+1=Yworsti(1+0.621*random)(9)

  Ybesti+1=Ybesti(1+0.5*random)(10)

  其中,random∈Gauss(0,1)。

  2.3權重因子的設置

  為了保證粒子群算法具有很好的收斂性,本文提出了在粒子速度上引入權重因子w,通過權重因子的調整來降低粒子在全局和局部最優中調節粒子的速度,從而保證粒子能夠獲得最優解。

  (1)全局解-線性微分策略

  在粒子群算法中,伴隨著迭代次數的不斷增加,粒子的速度會呈現不同的變化,本文采用典型的線性遞減策略,權重系數根據迭代次數的更新如下:

  w(t)=wmin+wT×t(11)

  式中,wmin表示權值系數的最小值,T為最大迭代次數,t為當前迭代次數。通過使用比較大的慣性權重能夠快速定位最優解,伴隨著算法的迭代次數不斷增大,粒子的速度逐漸減慢,收斂速度加快,算法性能提高,但在算法的初始階段,會隨著慣性權值的減少,搜索能力逐漸變弱,因此局部搜索能力變強,促使算法陷入局部最優。為了避免這種情況的發生,使用如下微分遞減的策略來進行慣性權重的更新:

  }HO{@1VIDAIRXV~V[5U%_[6.png

  (2)局部解非線性策略

  權重系數通過線性遞減策略改進后,算法的性能已經有了很大的提高,但由于線性遞減策略自身的特性導致算法容易陷入局部最優,因此采用非線性策略來解決局部最優解:

  EY${O}]NHAQ$18(OLQTU3W4.png

3基于改進的粒子群算法在節點定位中的研究

  3.1節點定位誤差分析

  設錨節點(xi,yi)(i=1,2,…,n)與未知節點(x,y)的實際距離為ri,i=1,2,…,n,測距誤差為εi,那么|ri-di|<εi,i=1,2,…,n。根據式(2)可知,(x,y)應該滿足如下約束條件:

  4_A{%F35C)M8APM[N05U%SQ.png

  式中如果f(x,y)的值最小,則未知節點的解與真實值之間的偏差就最小,而此時的坐標(x,y)為最優解,即滿足下式的未知節點坐標。

  8NV%2}IN`Q0M`J08]MI6K~O.png

  式(16)是一個非線性最優化問題,因此將它作為粒子群的目標函數,通過改進的粒子群算法來提高整個粒子群之間的信息交流能力,求出最優解,從而實現未知節點坐標的計算。

  3.2粒子群算法的修正誤差的步驟

  (1)通過分析,將無線傳感網中的節點定位問題轉換為非線性優化問題。

  (2)設置粒子群算法的相關參數 。

  (3)計算每一個粒子的速度和位置,并根據式(15)計算其適應度值,將適應度最大的對應的粒子的位置保存pi中,將子群體的最優位置保存到Pkg(表示第 k個子群的最優位置)中。

  (4)根據式(8)計算當前粒子群的粒子之間的平均距離averagedistance,根據公式(9)和(10)對粒子進行變異。

  (5)根據公式(10)和公式(11)計算權重,從而代入式(6)和(7)進行計算。

  (6)每個粒子位置與自身的歷史最優位置pi對比,選擇出最優的粒子位置。

  (7)每個粒子位置與所在子群的歷史最位置 Pkg進行比較,選擇出種群的最優位置Pkg。

  (8)當次數達到最大迭代次數,算法循環截止,計算每個子群的最優位置對應的適應度值,選取整個子群對應的最優位置為pg=min(P1g,P2g,…,Pkg,…,PNg);否則轉至步驟(4)繼續迭代。

  (9)輸出pg對應粒子的坐標作為待定位的傳感器定點坐標。

4仿真實驗

  為了說明本文算法在節點定位中具有的有效性,本文選擇在MATLAB2010平臺上進行仿真實驗。選擇計算機硬件配置為:酷睿i3 CPU,4GB DDR3內存,500 GB硬盤;軟件環境為Windows Xp。選取250個節點隨機分布在邊長為100 m的區域中,錨節點和未知節點的坐標位置隨機產生。將DVHop算法和本文算法在錨節點數量和測距誤差兩個方面進行對比。

  4.1錨節點數量下的比較

  設定錨節點所占的比例不超過10%,即不超過25, DVHop算法和本文算法錨節點數量與定位誤差變化曲線如圖1所示。從圖1可知,伴隨著錨節點個數的不斷增大,兩種算法的平均定位誤差均逐漸減小,本文算法的定位誤差明顯小于 DVHop算法的定位誤差,減少了19.17%。

  

001.jpg

  4.2測距誤差下的定位性能比較

  在實際的定位中測距誤差是無法避免的,DVHop 算法和本文算法定位誤差變化和測距誤差的關系曲線如圖2所示。從圖2可知,隨著測距誤差的不斷增多,兩種算法的定位誤差逐漸增大,相對于DvHop算法而言,本文算法的定位誤差分別減少了9.25%。

002.jpg

5結論

  本文分析WSN中基于DVHop定位算法精度低的原因,提出一種基于粒子群的DVHop算法,通過對粒子群算法的改進,使得算法的性能得到提高。將改進后的算法運用到了節點定位中,仿真實驗結果表明,相比于DVHop 算法,本文算法在定位精度方面具有明顯優勢,是一種簡單實用的無線傳感器網絡節點定位方法。

參考文獻

 ?。?] Li Mo,Liu Yunhao. Rendeved path: rangefree localization in anisotropic sensornetworks with holes[J].IEEE/ACM Transactions on Networking,2010,18(1):320 332.

  [2] LAZOS L,POOVENDRAN R.Highresolution robust localization for wireless sensor networks[J].IEEE Journal on Selected Areas in Communications,2006,24(2):233 246.

 ?。?] 魏玉宏,田杰,高志強.基于混合粒子群優化的DVHop定位算法[J].計算機測量與控制,2015,23(12):42144216.

 ?。?] 歐陽丹彤,何金勝,白洪濤.一種約束粒子群優化的無線傳感網絡節點定位算法[J].計算機科學,2011,38(7):46 50.

 ?。?] 夏少波.無線傳感器網絡DVHop定位算法的改進[J].計算機應用,2015,35(2):340 344.

 ?。?] 周杭霞,崔晨,葉佳駿.一種基于加權處理和誤差修的DVHop定位算法研究[J].傳感技術學報,2014,27(12):16991703.

 ?。?] 周彥,文寶,李建勛.無線傳感器網絡節點近點加權質心定位方法[J].計算機工程與應用,2012,48(1):87 89.

 ?。?] 程超,錢志鴻,付彩欣.一種基于誤差距離加權與跳段算法選擇的遺傳優化DV-Hop定位算法[J].電子與信息學報,2015,37(10):24182422.

  [9] 江濤.基于混合人工蜂群策略的改進DVHop定位算法[J],電子器件,2014,37(5):912 915.

 ?。?0] LIU F, LIU G Z.Markov chain analysis and the convergence speed of genetic algorithms[J].Systems Engineering Journal,1998, 13(4): 79 85

 ?。?1] Xue Wentao, Wu Xiaobei, Xu Zhiliang. An immune network algorithm based on double mutation operators[J]. Controland Decision,2008,23(12):1417 1422.


此內容為AET網站原創,未經授權禁止轉載。
热re99久久精品国产66热_欧美小视频在线观看_日韩成人激情影院_庆余年2免费日韩剧观看大牛_91久久久久久国产精品_国产原创欧美精品_美女999久久久精品视频_欧美大成色www永久网站婷_国产色婷婷国产综合在线理论片a_国产精品电影在线观看_日韩精品视频在线观看网址_97在线观看免费_性欧美亚洲xxxx乳在线观看_久久精品美女视频网站_777国产偷窥盗摄精品视频_在线日韩第一页
  • <strike id="ygamy"></strike>
  • 
    
      • <del id="ygamy"></del>
        <tfoot id="ygamy"></tfoot>
          <strike id="ygamy"></strike>
          永久555www成人免费| 国产日韩欧美综合精品| 欧美专区日韩视频| 亚洲区一区二| 韩日在线一区| 久久精品一区二区国产| 国产一区二区三区在线观看免费| 欧美激情中文字幕在线| 国产一区自拍视频| 美女脱光内衣内裤视频久久影院| 国产人成一区二区三区影院| 国产精品视频免费一区| 欧美亚洲免费高清在线观看| 久久精品99无色码中文字幕| 欧美日本不卡高清| 亚洲一区精彩视频| 欧美大尺度在线观看| 香蕉亚洲视频| 中文在线一区| 国产精品视频最多的网站| | 狠狠色狠狠色综合系列| 欧美日韩在线一二三| 亚洲欧美成人精品| 136国产福利精品导航网址应用| 能在线观看的日韩av| 亚洲午夜激情网页| 国产午夜一区二区三区| 国产精品毛片高清在线完整版| 一区二区三区免费在线观看| 国产精品国产自产拍高清av| 国产九色精品成人porny| 一本久道综合久久精品| 欧美另类综合| 在线亚洲欧美视频| 欧美影院午夜播放| 欧美色欧美亚洲另类七区| 亚洲欧美高清| 欧美精品www在线观看| 国产三级欧美三级| 欧美成人dvd在线视频| 国产日韩在线亚洲字幕中文| 欧美激情精品久久久久久黑人| 国内精品久久久| 中文日韩在线视频| 久久久www免费人成黑人精品| 美女露胸一区二区三区| 亚洲精品久久久久久久久| 久久黄色网页| 亚洲欧美国产日韩天堂区| 亚洲欧洲日产国产网站| 欧美日韩精品久久| 亚洲自拍啪啪| 欧美中文字幕视频在线观看| 欧美日韩免费一区二区三区视频| 欧美日韩国产一区精品一区| 久久男人av资源网站| 国产精品美女久久| 午夜精品在线看| 亚洲一级二级在线| 久久精品99国产精品酒店日本| 久久久久久精| 国产精品老牛| 欧美在线不卡视频| 欧美日本亚洲| 欧美日本高清视频| 久久久青草婷婷精品综合日韩| 国产精品永久| 亚洲精品国产系列| 欧美高清影院| 亚洲乱码视频| 亚洲一区二区三区乱码aⅴ蜜桃女| 在线亚洲免费视频| 欧美日韩三级视频| 国产精品一区二区久久精品| 香蕉成人啪国产精品视频综合网| 亚洲在线日韩| 亚洲综合大片69999| 久久一区免费| 欧美日韩精品一区二区三区| 国产精品自在在线| 亚洲影院在线观看| 亚洲国产岛国毛片在线| 欧美日韩高清在线| 尤物yw午夜国产精品视频明星| 在线播放国产一区中文字幕剧情欧美| 国产精品v欧美精品v日本精品动漫| 99精品免费| 狠狠久久综合婷婷不卡| 亚洲欧美日韩爽爽影院| 亚洲男女自偷自拍| 亚洲一线二线三线久久久| 国产精品久久久久一区| 国产精品乱子乱xxxx| 91久久午夜| 国产日韩精品视频一区二区三区| 午夜免费电影一区在线观看| 久久精品伊人| 国产日韩欧美在线播放不卡| 欧美日本二区| 久久亚洲午夜电影| 欧美精品在线一区二区三区| 免费看成人av| 久久精品一本| 国产精品女主播在线观看| 蜜臀av性久久久久蜜臀aⅴ| 国产精品久久福利| 午夜在线不卡| 国产一区二区久久久| 国产人久久人人人人爽| 欧美专区亚洲专区| 久久噜噜噜精品国产亚洲综合| 激情亚洲成人| 国产精品大片免费观看| 亚洲一区二区三区乱码aⅴ| 欧美午夜精品久久久| 久久亚洲午夜电影| 久久综合久久综合这里只有精品| 午夜一区二区三区在线观看| 国产精品福利在线观看| 欧美日韩精品免费在线观看视频| 欧美四级电影网站| 中国成人亚色综合网站| 亚洲欧美一级二级三级| 韩国av一区| 欧美福利视频在线观看| 亚洲欧美在线播放| 免费精品99久久国产综合精品| 久久人人爽人人| 另类人畜视频在线| 国产综合在线视频| 亚洲欧美一区二区三区极速播放| 欧美黄污视频| 国产精品你懂的在线| 夜夜嗨av一区二区三区| 国产欧美日韩免费看aⅴ视频| 国产精品视区| 久久久久久久网站| 免费成人在线视频网站| 欧美国产综合视频| 国产日韩在线看| 中国亚洲黄色| 国产亚洲人成网站在线观看| 一本色道久久综合亚洲精品小说| 欧美伊久线香蕉线新在线| 久久蜜桃精品| 久久久99免费视频| 好吊一区二区三区| 欧美成人精品在线视频| 亚洲男人天堂2024| 亚洲人精品午夜在线观看| 欧美日本精品在线| 一区二区高清在线观看| 国产亚洲高清视频| 在线欧美日韩国产| 欧美天天综合网| 美女露胸一区二区三区| 亚洲综合日本| 一本色道久久99精品综合| 日韩一二三在线视频播| 午夜伦欧美伦电影理论片| 欧美日韩一卡| 免费在线观看一区二区| 亚洲综合国产激情另类一区| 亚洲一区影院| 国产精品嫩草99av在线| 亚洲人成网站在线观看播放| 欧美一区中文字幕| 米奇777超碰欧美日韩亚洲| 亚洲欧美电影在线观看| 美女91精品| 国产综合第一页| 国产一区av在线| 午夜国产欧美理论在线播放| 欧美日韩专区| 午夜精品久久99蜜桃的功能介绍| 黄色av日韩| 久久精视频免费在线久久完整在线看| 久久字幕精品一区| 亚洲黄色av一区| 午夜精品久久久久99热蜜桃导演| 国产老女人精品毛片久久| 久久亚洲视频| 国产精品一区免费在线观看| 久久aⅴ国产欧美74aaa| 欧美激情第8页| 久久久久综合| 91久久久久久久久| 在线视频国产日韩| 欧美日韩妖精视频| 久久蜜桃av一区精品变态类天堂| 激情六月婷婷久久| 99在线视频精品| 在线播放一区| 怡红院精品视频在线观看极品| 久久国产精品网站| 在线一区二区日韩| 欧美一区二区三区电影在线观看| 国产精品久久久久永久免费观看| 欧美精品网站| 欧美日韩视频专区在线播放| 久久国产色av| 亚洲欧美日韩中文播放| 久久成人精品| 一区二区日韩欧美| 国产精品激情电影| 亚洲婷婷综合色高清在线| 欧美麻豆久久久久久中文| 欧美人牲a欧美精品| 欧美在线国产| 快she精品国产999| 国产久一道中文一区| 欧美色精品天天在线观看视频| 欧美日韩亚洲国产精品| 日韩视频中文字幕| 欧美成人在线免费观看| 久久天堂成人| 激情久久久久| 国产麻豆精品theporn| 亚洲欧美一区二区精品久久久| 欧美二区视频| 红桃视频国产一区| 国产精品区免费视频| 欧美大片免费| 国产精品人人爽人人做我的可爱| 欧美特黄a级高清免费大片a级| 国产日韩一级二级三级| 亚洲一区二区精品在线观看| 欧美韩国日本一区| 一区三区视频| 亚洲激情网站| 国产色产综合产在线视频| 亚洲免费影视第一页| 欧美精品一区在线发布| 亚洲视频中文| 国产精品亚洲片夜色在线| 欧美一区二区三区日韩| 亚洲黄色在线看| 国产啪精品视频| 欧美久久综合| 久久婷婷影院| 久久精品国产久精国产一老狼| 亚洲欧美日韩在线综合| 欧美成人四级电影| 另类国产ts人妖高潮视频| 久久精品视频网| 夜夜嗨av一区二区三区网站四季av| 红杏aⅴ成人免费视频| 亚洲一区二区三区中文字幕| 另类尿喷潮videofree| 欧美性感一类影片在线播放| 国产人成一区二区三区影院| 亚洲欧美激情一区二区| 亚洲激情一区二区三区| 一区二区在线视频| 久久亚洲春色中文字幕久久久| 欧美久久久久久久久| 亚洲另类春色国产| 欧美人成在线视频| 在线亚洲+欧美+日本专区| 国产精品久久久久久久久久免费| 激情久久久久久| 亚洲精品乱码久久久久久黑人| 国产区二精品视| 国内精品久久久久影院薰衣草| 亚洲电影激情视频网站| 久久久亚洲欧洲日产国码αv| 亚洲激情不卡| 久久一区免费| 在线观看视频欧美| 久久激情婷婷| 午夜免费日韩视频| 亚洲精品一区久久久久久| 亚洲欧美另类在线观看| 午夜日韩激情| 黄色成人小视频| 激情久久五月天| 性欧美超级视频| 欧美中文在线观看| 在线欧美电影| 欧美mv日韩mv国产网站| 欧美日本一区二区视频在线观看| 国产欧美日韩另类一区| 你懂的亚洲视频| 久久夜色精品国产| 在线精品视频在线观看高清| 久久网站热最新地址| 久久亚洲私人国产精品va媚药| 欧美日韩国产美女| 免费观看在线综合色| 久久精品日产第一区二区三区| 欧美成在线视频| 亚洲在线中文字幕| 欧美日韩亚洲一区三区| 国产精品jizz在线观看美国| 久久国产精品久久精品国产| 免费视频最近日韩| 欧美成va人片在线观看| 国产丝袜美腿一区二区三区| 国产婷婷色一区二区三区四区| 国产精品伊人日日| 亚洲高清在线视频| 久久久久久欧美| 亚洲黄色av| 国产亚洲一区在线播放| 亚洲午夜精品久久久久久浪潮| 国产伦精品一区二区三区高清版| 一区二区亚洲| 国产精品久久久一区二区| 欧美日韩性视频在线| 欧美日韩精品是欧美日韩精品| 亚洲永久免费视频| 欧美日韩国产综合一区二区| 久久蜜桃香蕉精品一区二区三区| 伊人久久综合97精品| 在线国产日韩| 久久av资源网| 久久久一二三| 国产欧美韩日| 羞羞视频在线观看欧美| 欧美一区在线看| 久久亚洲国产精品一区二区| 亚洲一区欧美| 国产视频一区在线| 国产欧美一区二区精品秋霞影院| 亚洲精品乱码久久久久久黑人| 欧美一激情一区二区三区| 久久午夜电影网|