999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

空時分組編碼的無人機中繼通信航跡規劃方法

2017-11-20 01:20:39劉海濤趙文強李春鳴李冬霞
航空學報 2017年9期
關鍵詞:規劃方法系統

劉海濤, 趙文強, 李春鳴, 李冬霞

中國民航大學 天津市智能信號與圖像處理重點實驗室, 天津 300300

空時分組編碼的無人機中繼通信航跡規劃方法

劉海濤*, 趙文強, 李春鳴, 李冬霞

中國民航大學 天津市智能信號與圖像處理重點實驗室, 天津 300300

無人機(UAV)中繼通信是實現遠距離點對點無線通信的一種重要技術手段。為提高無人機中繼通信系統鏈路傳輸的可靠性,提出基于空時分組編碼的無人機放大轉發中繼通信傳輸方案,并基于雙跳鏈路遍歷容量最大化的準則給出了無人機最佳航跡規劃方法,并進一步利用FM-EM算法給出了基于空時分組編碼無人機中繼通信系統的中斷概率及遍歷信道容量計算公式。計算機仿真表明:提出的方法顯著優于傳統的單發單收(SISO)無人機中繼通信系統。

無人機中繼; 航跡規劃; 空時分組編碼; FM-EM算法; 性能分析

相對于固定中繼通信,無人機(UAV)中繼通信具有通信距離遠、部署方便、中繼位置靈活可控、系統構建迅捷、維護成本低廉等諸多優點,因此無人機中繼通信系統在軍用與民用領域獲得廣泛的應用[1-2]。與此同時,由于無人機的高速機動性以及無人機能量受限等特性,無人機中繼通信也帶來一系列新問題[3-4],例如無人機中繼網絡拓撲結構快速變化,中繼無人機的航跡規劃等,因此針對無人機中繼通信帶來的特殊問題開展關鍵技術的研究具有重要的意義。

目前,無人機中繼通信的研究主要集中在以下5個領域:無人機空地信道的建模[5-6]、無人機中繼網絡優化[7-9]、無人機中繼通信傳輸方法[10-12]、無人機中繼通信性能分析[13-15]、中繼無人機的航跡規劃[16-20]。圍繞著中繼通信無人機的航跡規劃問題,相關研究如下:針對多用戶接入的無人機中繼通信系統的性能優化問題,文獻[16]基于各態歷經歸一化傳輸速率最大化的準則給出了無人機航跡規劃方法,然而在優化無人機航跡時僅考慮了用戶節點-無人機的單跳鏈路,沒有考慮無人機-基站鏈路對整體系統性能的影響。同樣針對多用戶接入的無人機中繼通信最優化問題,文獻[17]基于平均和速率最大化的準則及用戶最小速率最大化的準則提出兩種無人機航跡規劃方法,然而該文章也僅考慮了用戶節點-無人機單跳鏈路,沒有考慮無人機-基站鏈路對性能的影響。針對無人機中繼廣播通信系統的航跡規劃問題,文獻[18]基于最小各態歷經鏈路容量最大化的準則提出了一種低復雜度的無人機航跡規劃方法,該方法主要特點是無需精確知曉各個用戶的位置信息。針對點對點無人機中繼通信的航跡規劃問題,文獻[19]提出基于發射波束成形與接收波束成形的無人機中繼傳輸方法,并基于信噪比最大化的準則給出無人機航跡規劃方法,然而該方法要求接入節點發射機與基站接收機均需精確知曉信道的衰落信息,實際應用中發射機難以獲取信道的衰落信息,因此該方法應用于實際系統非常困難。

為提高點對點無人機中繼通信系統鏈路傳輸的可靠性,同時降低無人機中繼通信系統實現的復雜度,本文在文獻[7]的基礎上提出了基于空時分組編碼的無人機放大轉發中繼傳輸方法,給出相應的無人機中繼通信系統模型,以該模型為基礎,推導給出無人機中繼通信系統遍歷容量計算公式,基于系統遍歷容量最大化的準則給出了中繼無人機的最佳航跡規劃方法,并利用FM-EM(Finite Mixture-Expectation Maximum)算法給出了系統的中斷概率與遍歷容量計算公式,最后通過仿真驗證了論文所提方法的正確性及有效性。本文存在3個方面的貢獻:① 提出了基于空時分組編碼(STBC)的無人機中繼通信系統無人機航跡規劃方法;② 理論分析給出所提出傳輸方法的中斷概率及遍歷容量計算公式;③ 證明了空時分組編碼的無人機中繼通信系統優于傳統單發單收無人機中繼傳輸方法。

1 空時分組編碼的無人機中繼系統模型

1.1 系統模型

圖1給出了空時分組編碼的無人機中繼通信系統的模型。系統由可移動式接入節點(Access Point, AP)、高速運動的固定翼無人機及固定地面基站(Base Station, BS)組成。假設AP與BS的距離很遠,系統不存在AP與BS之間的直達通信鏈路,為實現AP與BS節點間的相互通信,必須通過UAV中繼來實現AP與BS之間的通信,為敘述方便,本文僅考慮AP至BS節點的通信。為提高UAV中繼通信系統的鏈路傳輸可靠性,同時降低系統實現的復雜性,本文采用空時分組編碼的UAV中繼通信方案,AP采用兩天線空時分組編碼發射分集方式傳輸信號,UAV采用單天線放大轉發協議中繼信號,BS采用多天線最大比值合并(Maximal-Ratio Combining, MRC)方法合并信號以提高鏈路傳輸的可靠性。

圖1 空時分組編碼的無人機中繼通信系統Fig.1 UAV relay communication system with space-time block coding

1.2 信號模型

(1)

(2)

(3)

然后UAV以一定發射功率將信號轉發至BS。在第3個時隙,BS兩個接收天線(為敘述方便,BS使用兩個接收天線,本方法容易推廣到BS使用多個接收天線的場合)接收到的信號表示為

(4)

在第4個時隙,基站兩個天線接收到的信號表示為

(5)

(6)

(7)

(8)

假設BS接收機通過信道估計可精確獲得所有信道的衰落系數,則BS接收機進一步按式(9)給出方法對BS兩個時隙接收到的信號進行相關合并處理:

(9)

將式(4)與式(5)分別代入式(9),化簡后得到

(10)

式中:

(11)

(12)

進一步假設無人機的飛行高度始終保持為h,且無人機以恒定速度v飛行,則UAV在t時刻的位置可根據其t-Δt時刻位置信息,并利用以下方程更新獲得[17,21]:

(13)

式中:δt為t時刻UAV的航向角,其滿足δt-Δt-δmax≤δt≤δt-Δt+δmax,其中δmax為UAV的最大轉彎角;Δt為UAV位置更新的時間間隔。將式(13)代入式(12),化簡后得到

(14)

式(14)表明:在t-Δt時刻UAV的位置給定后,t時刻AP-UAV及UAV-BS的距離僅決定于t時刻UAV的航向角δt。

2 無人機航跡規劃方法與性能分析

2.1 UAV最佳航跡規劃方法

在無人機中繼通信系統中,由于AP與UAV均處于運動狀態,AP-UAV及UAV-BS的距離處于不斷地變化中,為保證無人機中繼通信系統鏈路傳輸可靠性最優,需要根據AP移動路徑的變化優化UAV的飛行航跡。考慮到遍歷容量是衡量無線通信系統鏈路傳輸可靠性的主要技術指標,因此論文基于系統遍歷容量最大化準則來優化無人機的航跡。t時刻無人機中繼通信系統的遍歷容量表示為[22]

(15)

式中:由于直接對式(15)進行統計平均運算非常困難,論文利用泰勒級數的方法展開式(15),進一步考慮到展開式中高階項的取值較小,且對高階項進行統計平均運算的復雜度較高,因此僅保留展開式中的線性項得到

(16)

式中:E(rout,t)為t時刻接收機解調器輸入的平均信噪比,其表達式為

(17)

考慮到式(17)計算仍比較困難,再次利用泰勒級數對式(17)展開,同理取展開式中的第1項,式(17)可簡化為

(18)

根據式(18)可得出,影響t時刻UAV中繼通信系統BS接收機解調器輸入平均信噪比的影響因素為大尺度衰耗(dAU,t,dUB,t)。因此,將式(14)代入式(18),可觀測得到以下結論:在t-Δt時刻UAV位置給定情況下,t時刻UAV中繼通信系統BS接收機解調器輸入平均信噪比僅決定于t時刻UAV的航向角δt。為使t時刻中繼通信系統鏈路容量最大化,本文基于最大化遍歷容量準則來優化UAV的航向角,以上優化準則可表述為

s.t. |δt-δt -Δt|≤δmax

(19)

此外進一步考慮到log2(·)為單調遞增函數,因此式(19)的優化問題最后表述為

s.t. |δt-δt -Δt|≤δmax

(20)

2.2 性能分析

(21)

式中:b(opt)為b的最佳值。

2.2.1 系統的中斷概率

中斷概率是衡量無線通信系統鏈路傳輸可靠性的一項重要技術指標,其定義為接收信噪比低于某一門限值rth的概率[22],即

(22)

(23)

利用式(23)得到的t時刻BS接收機解調器輸入瞬時信噪比的概率密度函數,式(22)可進一步化簡為

(24)

式中:Q(η)為高斯Q函數,具體可表示為

(25)

2.2.2 系統的遍歷容量

為了計算得到系統的遍歷容量,本文首先通過FM-EM算法來估計得到t時刻系統遍歷容量Ct的瞬時概率密度函數fCt(c),隨后利用估計得到的概率密度函數fCt(c),可計算得到UAV中繼系統的遍歷容量為

(26)

3 數值仿真

3.1 仿真參數設置

3.2 仿真結果與分析

圖2為UAV最佳航跡隨著AP運動路線變化的曲線,圖2橫坐標與縱坐標分別代表直角坐標系的x軸與y軸,藍色的曲線代表AP的運動路線,粉色曲線代表UAV的最佳航跡。曲線比較表明:① 當AP運動時,UAV通過調整其航向角使UAV跟隨AP路徑的變化飛行;② 當 UAV跟隨AP運動飛行時,UAV的航跡呈現為圓形,原因為由于UAV的飛行速度大于AP的移動速度,且UAV飛行受到最大轉彎角的限制,為了保證UAV中繼通信系統的性能達到最佳,UAV在某些時刻必須以繞圓方式飛行。

圖3給出了不同最大轉彎角對UAV飛行航跡的影響(最大轉彎角分別為10° 與15°),其中,黑色的曲線是最大轉彎角為15° 時UAV最佳飛行航跡,粉色的曲線表示的是最大轉彎角為10° 時UAV的最佳航跡。曲線比較表明:當最大轉彎角取值較小時,無人機繞圓飛行的半徑較大;當最大轉彎角取值較大時,無人機繞圓飛行的半徑較小。產生以上現象的原因是:無人機飛行速度給定情況下,最大轉彎角的增大意味著無人機可選擇最佳航向角的范圍增大, 導致無人機繞圓飛行的半徑變小。

圖2 UAV最佳航跡隨AP運動路線變化的曲線Fig.2 Optimal trajectories of UAV changing with the mobility of AP

圖3 不同最大轉彎角對UAV的飛行航跡的影響 Fig.3 Trajectories of UAV with different maximum heading

圖4和圖5分別給出了無人機中繼通信系統的中斷概率及遍歷容量隨時間變化的曲線。圖4和圖5包含兩類曲線,標識為實線的曲線代表理論計算得到的曲線,標識為“o”型的曲線代表蒙特卡羅仿真得到的曲線;其中紅色曲線代表的是采用空時分組編碼的無人機中繼系統的中斷概率和遍歷容量曲線,黑色曲線代表的是采用單發單收(Single-Input Single-Output, SISO)無人機中繼系統的中斷概率與遍歷容量曲線。曲線比較表明:① 理論計算得到曲線與蒙特卡羅仿真得到的曲線完全一致,表明本文提出的航跡規劃方法的正確性;② 本文提出的空時分組編碼的無人機中繼通信方案顯著優于單發單收無人機中繼通信系統方案,表明論文提出方法的有效性。

圖4 中繼系統中斷概率隨時間的變化曲線(δmax=10°) Fig.4 Outage probability of relay system vs time(δmax=10°)

圖5 中繼系統遍歷容量隨時間的變化曲線(δmax=10°) Fig.5 Ergodic capacity of relay system vs time(δmax=10°)

圖6與圖7分別顯示最大轉彎角對UAV中繼通信系統中斷概率及遍歷容量性能的影響(最大轉彎角分別為10° 與15°)。曲線比較表明:① 不同最大轉彎角情況下,空時分組編碼的無人機中繼通信方案仍優于單發單收無人機中繼通信方案;② UAV的最大轉彎角對UAV中繼通信系統中斷概率及信道遍歷容量的影響不大。

圖6 最大轉彎角對中繼系統中斷概率性能的影響 Fig.6 Outage probability of relay system withdifferent maximum turning angle

圖7 最大轉彎角對中繼系統遍歷容量性能的影響Fig.7 Ergodic capacity of relay system with different maximum turning angle

4 結 論

1) 本文提出了空時分組編碼的無人機中繼通信系統無人機最佳航跡規劃方法。

2) 空時分組編碼的無人機中繼通信系統優于單發單收無人機中繼通信系統。

3)無人機的最大轉彎角對無人機中繼通信系統的中斷概率及遍歷容量影響不大。

[1] FREW E W, BROWN T X. Airborne communication networks for small unmanned aircraft systems[J]. Proceedings of the IEEE, 2008, 96(12): 2008-2027.

[2] HAYAT S, YANMAZ E, MUZAFFAR R. Survey on unmanned aerial vehicle networks for civil applications: A communications viewpoint[J]. IEEE Communications Surveys & Tutorials, 2016, 18(4): 2624-2661.

[3] ZENG Y, ZHANG R, TENG J L. Wireless communications with unmanned aerial vehicles: Opportunities and challenges[J]. IEEE Communications Magazine, 2016, 54(5): 36-42.

[4] GUPTA L, JAIN R, VASZKUN G. Survey of important issues in UAV communication networks[J]. IEEE Communications Surveys & Tutorials, 2016, 18(2): 1123-1152.

[5] SIMUNEK M, FONTAN F P, PECHAC P. The UAV low elevation propagation channel in urban areas: Statistical analysis and time-series generator[J]. IEEE Transactions on Antennas and Propagation, 2013, 61(7): 3850-3858.

[6] MATOLAK D W, SUN R. Unmanned aircraft systems: Air-ground channel characterization for future applications[J]. IEEE Vehicular Technology Magazine, 2015, 10(2): 79-85.

[7] PALAT R C, ANNAMALAU A, REED J R. Cooperative relaying for ad-hoc ground networks using swarm UAVs[C]//2005 Military Communications Conference. Piscataway, NJ: IEEE Press, 2005: 1588-1594.

[8] HAN Z, SWINDLEHURST A L, LIU K. Optimization of MANET connectivity via smart deployment/movement of unmanned air vehicles[J]. IEEE Transactions on Vehicular Technology, 2009, 58(7): 3533-3546.

[9] 徐贊新, 袁堅, 王鉞, 等. 一種支持移動自組網通信的多無人機中繼網絡[J]. 清華大學學報: 自然科學版, 2011, 51(2): 150-155.

XU Z X, YUAN J, WANG Y, et al. UAV relay network to provide communications in mobile ad hoc networks[J]. Journal of Tsinghua University: Science and Technology, 2011, 51(2): 150-155 (in Chinese).

[10] PALAT R C, ANNAMALAU A, REED J H. Cooperative relaying for ad-hoc ground networks using swarm UAVs[C]//2005 IEEE Military Communications Conference. Piscataway, NJ: IEEE Press, 2005: 1588-1594.

[11] CHENG C M, HSIAO P H, KUNG H T, et al. Maximizing throughput of UAV-relaying networks with the load-carry-and-deliver paradigm[C]//2007 IEEE Wireless Communications and Networking Conference. Piscataway, NJ: IEEE Press, 2007: 4417-4424.

[12] ONO F, OCHIAI H, MIURA R. A wireless relay network based on unmanned aircraft system with rate optimization[J]. IEEE Transactions on Wireless Communications, 2016, 15(11): 7699-7708.

[13] ABUALHAOL I Y, MATALGAH M M. Performance analysis of multi-carrier relay-based UAV network over fading channels[C]//2010 IEEE GLOBECOM Workshops. Piscataway, NJ: IEEE Press, 2010: 1811-1815.

[14] 歐陽鍵, 莊毅, 薛羽, 等. 非對稱衰落信道下無人機中繼傳輸方案及性能分析[J]. 航空學報, 2013, 34(1): 130-140.

OUYANG J, ZHUANG Y, XUE Y, et al. UAV relay transmission scheme and its performance analysis over asymmetric fading channels[J]. Acta Aeronautica et Astronautica Sinica, 2013, 34(1): 130-140 (in Chinese).

[15] 林敏, 魏恒, 歐陽鍵, 等. 無人飛行器中繼雙跳無線鏈路中的優化設計及性能分析[J]. 系統工程與電子技術, 2015, 37(6): 1391-1398 (in Chinese).

LIN M, WEI H, OUYANG J, et al. Optimal design and performance analysis for two-hop wireless links with UAV relaying[J]. Systems Engineering and Electronics, 2015, 37(6): 1391-1398 (in Chinese).

[16] ZHAN P, YU K, SWINDLEHURST A L. Wireless relay communications with unmanned aerial vehicles: Performance and optimization[J]. IEEE Transactions on Aerospace and Electronics Systems, 2011, 47(3): 2068-2085.

[17] FENG J, SWINDLEHURST A L. Optimization of UAV heading for the ground-to-air uplink[J]. IEEE Journal on Selected Areas in Communications, 2012, 30(5): 993-1005.

[18] CHOI D H, JUNG B H, DAN K S. Low-complexity maneuvering control of a UAV-based relay without location information of mobile ground nodes[C]//2014 IEEE Symposium on Computers and Communications (ISCC). Piscataway, NJ: IEEE Press, 2014: 1-6.

[19] OUYANG J, ZHUANG Y, LIN M, et al. Optimization of beamforming and path planning for UAV-assisted wireless relay networks[J]. Chinese Journal of Aeronautics, 2014, 27(2): 313-320.

[20] 符小衛, 程思敏, 高曉光. 無人機協同中繼過程中的路徑規劃與通信優化[J]. 系統工程與電子技術, 2014, 36(5): 890-894.

FU X W, CHENG S M, GAO X G. Path planning and communication optimizing in multi-UAVs cooperative communication relay[J]. Systems Engineering and Electronics, 2014, 36(5): 890-894 (in Chinese).

[21] DUBINS L E. On curves of minimal length with a constraint on average curvature, and with prescribed initial and terminal positions and tangents[J]. American Journal of Mathematics, 1957, 79(3): 497-516.

[22] SIMON M K, ALOUINI M S. Digital communication over fading channels[M]. 2nd ed. Hoboken, NJ: John Wiley and Sons, Inc., 2005: 20-24.

[23] ZAIDI Z R, MARK B L, THOMAS R K. A two-tier representation of node mobility in ad hoc networks[C]//2004 First IEEE Communications Society Conference on Sensor and Ad Hoc Communications and Networks. Piscataway, NJ: IEEE Press, 2004: 153-161.

[24] MCLACHLAN G, PEEL D. Finite mixture models[M]. New York: Wiley, 2000: 47-50.

附錄A: FM-EM算法

假設在FM模型中,被估計的概率密度函數fZ(z)由N個高斯分布的概率密度函數加權求和組成。因此對于單變量的FM模型,隨機變量Z的概率密度函數可表示為[24]

(A1)

(A2)

對于上述FM模型中混合項的參數的估計可利用經典算法——EM算法。利用EM算法估計FM模型中混合項的參數,首先,①要確定混合項的個數N;②被估計參數的初始值Θ0;③ EM算法迭代的終止條件ε。其次,可利用EM算法的迭代公式不斷更新這些參數值直到滿足EM算法的終止條件。下面將給出上述參數的具體迭代公式:

(A3)

k=1,2,…,N;p=1,2,…,M

(A4)

(A5)

再根據式(A3)對μk求偏導數可得μk的更新公式為

(A6)

(A7)

現在把應用EM算法進行FM模型估計的實施步驟歸納如下:

步驟1 確定FM模型中混合項的項數N。

步驟3 利用式(A4)計算每一個觀測數據zp所對應的后驗概率。

步驟4 利用式(A5)~式(A7)更新參數集Θ。

步驟5 重復第3步與第4步直到EM算法收斂。

附錄B: AP運動模型

AP運動模型采用一階自回歸模型[23]。假定AP在t-Δt時刻的狀態用sA,t-Δt表示,即

(B1)

sA,t=TAsA,t -Δt+wA,t

(B2)

(B3)

(責任編輯: 蘇磊)

*Corresponding author. E-mail: hitliucauc@qq.com

Path planning method for UAV relay communication system with space-time block coding

LIU Haitao*, ZHAO Wenqiang, LI Chunming, LI Dongxia

TianjinKeyLaboratoryforAdvancedSignalProcessing,CivilAviationUniversityofChina,Tianjin300300,China

Unmanned Aerial Vehicle (UAV) relay communication is an important technical scheme for long range point-to-point wireless communications. To improve the reliability of transmission link for the UAV relay communication system, a novel UAV relay transmission scheme is proposed based on space-time block coding. According to the criterion of maximizing the average mutual information of the dual-hop link of the UAV relaying communication system, an optimal path planning method for the UAV relay communication system is presented. The outage probability and ergodic capacity of the UAV relay communication system are derived using Finite-Mixture with Expectation-Maximization (FM-EM) algorithm. The computer simulation results indicate that the proposed scheme is superior to the traditional Single-Input Single-Output (SISO) UAV relay communication scheme.

Unmanned Aerial Vehicle (UAV) relay; path planning; space-time block coding; FM-EM algorithm; performance analysis

2016-12-13; Revised: 2017-03-10; Accepted: 2017-03-25; Published online: 2017-04-11 10:21

URL: www.cnki.net/kcms/detail/11.1929.V.20170411.1021.004.html

National Key Research and Development Plan (2016YFB0502402)

2016-12-13; 退修日期: 2017-03-10; 錄用日期: 2017-03-25; 網絡出版時間: 2017-04-11 10:21

www.cnki.net/kcms/detail/11.1929.V.20170411.1021.004.html

國家重點研發計劃(2016YFB0502402)

*通訊作者.E-mail: hitliucauc@qq.com

劉海濤, 趙文強, 李春鳴, 等. 空時分組編碼的無人機中繼通信航跡規劃方法[J]. 航空學報, 2017, 38(9): 321048. LIU H T, ZHAO W Q, LI C M, et al. Path planning method for UAV relay communication system with space-time block coding[J]. Acta Aeronautica et Astronautica Sinica, 2017, 38(9): 321048.

http://hkxb.buaa.edu.cn hkxb@buaa.edu.cn

10.7527/S1000-6893.2017.321048

猜你喜歡
規劃方法系統
Smartflower POP 一體式光伏系統
工業設計(2022年8期)2022-09-09 07:43:20
WJ-700無人機系統
ZC系列無人機遙感系統
北京測繪(2020年12期)2020-12-29 01:33:58
規劃引領把握未來
連通與提升系統的最后一塊拼圖 Audiolab 傲立 M-DAC mini
快遞業十三五規劃發布
商周刊(2017年5期)2017-08-22 03:35:26
多管齊下落實規劃
中國衛生(2016年2期)2016-11-12 13:22:16
用對方法才能瘦
Coco薇(2016年2期)2016-03-22 02:42:52
迎接“十三五”規劃
四大方法 教你不再“坐以待病”!
Coco薇(2015年1期)2015-08-13 02:47:34
主站蜘蛛池模板: 超碰色了色| 天堂在线www网亚洲| 国产一区亚洲一区| 亚洲成人网在线观看| 国产成人av一区二区三区| 在线观看网站国产| 国产sm重味一区二区三区| 99偷拍视频精品一区二区| 国产丝袜无码精品| 色亚洲成人| 国产区精品高清在线观看| 女人av社区男人的天堂| 久久一日本道色综合久久| 精品国产91爱| 日本www在线视频| 天堂岛国av无码免费无禁网站| 91最新精品视频发布页| 日韩精品一区二区深田咏美| 国产一级妓女av网站| 日韩a在线观看免费观看| 香蕉色综合| 中文字幕啪啪| 亚洲人成影院在线观看| 在线精品自拍| 真实国产乱子伦视频| 婷婷色中文| 亚洲无码在线午夜电影| 极品私人尤物在线精品首页| 日本欧美成人免费| 国产国拍精品视频免费看 | 成人午夜网址| 国产福利免费视频| 热九九精品| 国产丝袜第一页| 久久久久九九精品影院| 欧美性爱精品一区二区三区 | 久久久噜噜噜久久中文字幕色伊伊 | 免费人成网站在线观看欧美| 国产区91| 欧美a网站| 亚洲一级毛片在线观| 中文字幕乱码中文乱码51精品| 欧美精品v| 亚洲小视频网站| 久久综合丝袜长腿丝袜| 国产精品色婷婷在线观看| 免费在线国产一区二区三区精品| 女人av社区男人的天堂| 亚洲一区毛片| 99热免费在线| 中文字幕乱码二三区免费| 在线观看国产网址你懂的| 欧美国产视频| 免费在线成人网| 91在线无码精品秘九色APP| 91在线播放免费不卡无毒| 欧美精品另类| 色噜噜在线观看| 免费国产高清精品一区在线| www成人国产在线观看网站| 色窝窝免费一区二区三区| 亚洲资源站av无码网址| 国产男女免费视频| 国产最新无码专区在线| 国产成人艳妇AA视频在线| 国产自在线播放| 国产一区在线观看无码| 日韩精品一区二区三区免费| 国产在线91在线电影| 91亚洲视频下载| 久久人午夜亚洲精品无码区| 午夜精品久久久久久久无码软件 | 97在线免费| 欧美区一区| 拍国产真实乱人偷精品| 国产爽妇精品| 亚洲男人天堂2020| 成人福利在线视频| 福利在线不卡一区| 精品国产一二三区| 亚洲an第二区国产精品| 欧洲日本亚洲中文字幕|