《電子技術(shù)應用》
您所在的位置:首頁 > 嵌入式技術(shù) > 設計應用 > 偽衛(wèi)星網(wǎng)絡時頻同步系統(tǒng)設計與實現(xiàn)
偽衛(wèi)星網(wǎng)絡時頻同步系統(tǒng)設計與實現(xiàn)
2018年電子技術(shù)應用第6期
紀元法1,2,梁 濤1,2,孫希延1,2,嚴素清1,2,盧偉軍1,2
1.桂林電子科技大學 信息與通信學院,廣西 桂林541004;2.廣西精密導航技術(shù)與應用重點實驗室,廣西 桂林541004
摘要: 為實現(xiàn)偽衛(wèi)星自組織網(wǎng)絡高精度時頻同步,設計了一種應用SCPA結(jié)構(gòu)的偽衛(wèi)星時頻同步系統(tǒng)。該系統(tǒng)以FPGA和DSP作為核心基帶處理芯片,集成了高速數(shù)模轉(zhuǎn)換與模數(shù)轉(zhuǎn)換電路以及上下變頻電路。詳細介紹了偽衛(wèi)星雙向偽距測量技術(shù)和時頻同步設計方案及實現(xiàn)。系統(tǒng)測試結(jié)果表明,偽衛(wèi)星系統(tǒng)星間載波相位同步誤差小于0.1 Hz,時間同步精度優(yōu)于2 ns,可達到偽衛(wèi)星自組織網(wǎng)絡系統(tǒng)定位精度與授時精度要求。
中圖分類號: TN967.1
文獻標識碼: A
DOI:10.16157/j.issn.0258-7998.174680
中文引用格式: 紀元法,梁濤,孫希延,等. 偽衛(wèi)星網(wǎng)絡時頻同步系統(tǒng)設計與實現(xiàn)[J].電子技術(shù)應用,2018,44(6):39-43.
英文引用格式: Ji Yuanfa,Liang Tao,Sun Xiyan,et al. Design and implementation of pseudolites time and frequency synchronization system[J]. Application of Electronic Technique,2018,44(6):39-43.
Design and implementation of pseudolites time and frequency synchronization system
Ji Yuanfa1,2,Liang Tao1,2,Sun Xiyan1,2,Yan Suqing1,2,Lu Weijun1,2
1.School of Information and Communication,Guilin University of Electronic Technology,Guilin 541004,China; 2.Guangxi Key Laboratory of Precision Navigation Technology and Application,Guilin 541004,China
Abstract: In order to achieve high precision of the time and frequency synchronization of pseudolite automatic network,a pseudolite time and frequency synchronization system using the structure of Self-Calibrating Pseudolite Array(SCPA) is designed. The system takes FPGA+DSP as the core processor, and integrates high-speed A/D, D/A conversion circuit and up and down converter circuit. The design and implementation of pseudolite two-way pseudo distance measurement technology and the time and frequency synchronization scheme are introduced in detail. The test results show that the phase synchronization error of carrier phase is less than 0.1 Hz and the time synchronization accuracy is better than 2 ns, which provids a significant guarantee for the positioning and timing accuracy of the pseudolite automatic system.
Key words : pseudolite;SCPA;two-way ranging;time and frequency synchronization

0 引言

    偽衛(wèi)星系統(tǒng)作為一種無線導航發(fā)射設備,可以用作增強GPS星座,也能構(gòu)成獨立的導航定位系統(tǒng)。與GPS衛(wèi)星裝備的原子鐘不同,偽衛(wèi)星的時鐘通常選擇精度不高的低端時鐘,會產(chǎn)生鐘漂誤差[1-2]。依據(jù)衛(wèi)星導航定位原理,為保證用戶接收機的定位精度和授時精度,系統(tǒng)中的偽衛(wèi)星必須保持時間同步。

    偽衛(wèi)星時間同步系統(tǒng)采用無線雙向微波時間同步方案,具有組網(wǎng)靈活、可擴展性好、綜合成本低等優(yōu)點。依據(jù)測量終端間傳播路徑相同特性,基于偽距測量模式的無線雙向時間同步系統(tǒng)可以最大限度地消除由路徑不同引起的誤差,實現(xiàn)精確測距與時間同步[3]。根據(jù)測量終端發(fā)射機與接收機信號閉環(huán)傳輸鏈路,可實現(xiàn)終端間載波頻率同步。

1 系統(tǒng)總體構(gòu)架

    本文設計的偽衛(wèi)星網(wǎng)絡時頻同步系統(tǒng)由1個主站偽衛(wèi)星、4個或以上從站偽衛(wèi)星構(gòu)成,各站均采用具有自校準功能的偽衛(wèi)星收發(fā)器陣列結(jié)構(gòu)(Self-Calibrating Pseudolite Array,SCPA),每個收發(fā)器主要由發(fā)射機和接收機兩個部分組成[2]。采用主從模式的自差收發(fā)器結(jié)構(gòu),每個收發(fā)器使用一個時鐘源,實現(xiàn)發(fā)射部分和接收部分時間基準統(tǒng)一。通過信號分路與合路設計,調(diào)節(jié)發(fā)射信號功率,接收機同時接收天線端與同源發(fā)射端的偽衛(wèi)星信號,形成信號閉環(huán)收發(fā),可減小收發(fā)器系統(tǒng)測量誤差。偽衛(wèi)星系統(tǒng)結(jié)構(gòu)組成如圖1所示。

qrs3-t1.gif

    偽衛(wèi)星系統(tǒng)的收發(fā)器主要包括射頻信號收發(fā)單元和中頻信號處理單元。射頻信號收發(fā)單元主要發(fā)射與接收偽衛(wèi)星頻點的導航信號,并完成中頻信號與射頻信號的轉(zhuǎn)換。射頻上變頻模塊通過對發(fā)射基帶模擬中頻信號與本振混頻生成偽衛(wèi)星射頻信號,通過分路器,一路信號經(jīng)發(fā)射天線播發(fā),另一路信號經(jīng)功率調(diào)節(jié)器發(fā)送到接收端;射頻接收單元下變頻模塊通過對射頻信號與本振混頻產(chǎn)生接收端模擬中頻信號。

    中頻信號處理單元是以FPGA和DSP作為核心基帶處理芯片,主要由發(fā)射信號基帶處理模塊、D/A數(shù)模轉(zhuǎn)換模塊、接收信號基帶處理模塊、A/D模數(shù)轉(zhuǎn)換模塊和時間同步模塊組成。發(fā)射信號基帶處理模塊實現(xiàn)偽衛(wèi)星導航信號的電文編碼與信號調(diào)制,經(jīng)過D/A數(shù)模轉(zhuǎn)換電路產(chǎn)生模擬中頻信號;接收信號基帶處理模塊接收A/D模數(shù)轉(zhuǎn)換電路輸出的數(shù)字中頻信號,實現(xiàn)信號的捕獲、跟蹤與電文解碼,完成與上位機界面的交互;時間同步模塊主要完成從站與主站鐘差的測量與修正,產(chǎn)生同步的時間信號。

2 系統(tǒng)主要硬件電路設計

2.1 下變頻電路設計

    本系統(tǒng)下變頻設計電路選用潤芯公司生產(chǎn)的一款高度集成的射頻芯片RX3007,片上集成了鏡頻抑制混頻器、帶通濾波器、自動增益控制電路、壓控振蕩器、中頻放大器、模數(shù)變換器等電路;支持GPS L1/BD2 B1信號雙通道同時工作;通道噪聲系數(shù)小于2.5 dB,通道增益110 dB,支持有源和無源天線工作模式;參考時鐘輸入范圍為10 MHz~40 MHz,可通過SPI控制接口配置模擬中頻或數(shù)字中頻輸出。本設計A/D模數(shù)轉(zhuǎn)換使用該芯片2 bit模數(shù)轉(zhuǎn)換器,分別對模擬信號進行采樣,采樣時鐘為 16.368 MHz,將頻率為4.092 MHz的模擬中頻信號量化成數(shù)字中頻信號,以SIGN、MAG碼輸出給FPGA基帶芯片。下變頻模塊電路圖如圖2所示。

qrs3-t2.gif

2.2 D/A數(shù)模轉(zhuǎn)換電路設計

    數(shù)模轉(zhuǎn)換電路設計選用Analog Devices公司的AD9744低功耗14 bit數(shù)模轉(zhuǎn)換器,采樣時鐘輸入支持210 MSPS轉(zhuǎn)換速率。輸出端設計采用單電源直流差分耦合電路,模擬差分輸出經(jīng)AD8041放大器實現(xiàn)高速數(shù)據(jù)緩沖。FPGA發(fā)射基帶處理模塊輸出包含B1和L1頻點的偽衛(wèi)星數(shù)字中頻信號,經(jīng)數(shù)模轉(zhuǎn)換器分別生成中頻頻率為11.098 MHz和25.42 MHz的模擬中頻信號。D/A數(shù)模轉(zhuǎn)換電路如圖3所示。

qrs3-t3.gif

2.3 上變頻電路設計

    射頻上變頻單元主要包括頻率合成器、混頻器和濾波器。頻率合成器選用Silicon LABS公司的SI41XX系列芯片,通過MCU單片機配置本振頻率為1 550 MHz;混頻器選用Mini-Circuits公司的無源混頻器JMS-11,中頻信號與本振混頻得到頻點為1 561.098 MHz和1 575.42 MHz的偽衛(wèi)星射頻信號,完成頻譜搬移;濾波器選用臺灣嘉碩科技公司的TA1166A聲表面濾波器(SAW),該濾波器中心頻率為1 575.5 MHz,帶寬為30 MHz,插入損耗不超過3.0 dB,經(jīng)濾波后可剔除不必要的信號及雜訊。射頻上變頻模塊電路圖如圖4所示。

qrs3-t4.gif

3 系統(tǒng)關(guān)鍵技術(shù)

3.1 雙向偽距測量與時間同步技術(shù)設計

    雙向偽距測量通過主站與從站的收發(fā)器設備,接收端利用偽碼和載波相位跟蹤結(jié)果,得到偽距測量值,構(gòu)建雙向測距方程,從而實現(xiàn)兩站間的距離測量與時間同步。DSP接收端采用載波相位平滑偽距算法設計,利用精確的載波相位測量值對粗糙的偽碼測量值進行平滑處理,提高偽距測量值的測量精度[4]。雙向偽距測量的原理圖如圖5所示。其中,Ti為偽距測量值,ti為傳輸時間,qrs3-t4-x1.gif為發(fā)射時延,qrs3-t4-x2.gif為接收時延,Δt為鐘差。

qrs3-t5.gif

    (1)主站發(fā)射端在本地時間0時刻發(fā)射偽衛(wèi)星射頻信號A,主站與從站接收端跟蹤主站信號A,主站通道1偽距測量值TA1

qrs3-gs1.gif

    (2)從站接收端跟蹤信號A且?guī)匠晒螅瑥恼綟PGA時間同步信號處理模塊啟動從站本地時間計數(shù),得到從站通道1的偽距測量值TB1

    qrs3-gs2.gif

    (3)從站FPGA發(fā)射端根據(jù)本地時間同步產(chǎn)生從站信號B,得到鐘差Δt,主站與從站接收端跟蹤從站信號B,主站通道2和從站通道2得到各自的偽距測量值TA2、TB2

     qrs3-gs3-4.gif

    (4)從站DSP發(fā)射端將得到的偽距觀測量寫入電文,主站得到同源發(fā)射端與從站的偽距測量結(jié)果,構(gòu)建雙向測距方程,得到主站偽衛(wèi)星和從站偽衛(wèi)星的時鐘差Δt和兩站間的距離時間tD

qrs3-gs5-6.gif

    時間同步技術(shù)主要在從站偽衛(wèi)星收發(fā)器中實現(xiàn),DSP程序主要完成信號跟蹤環(huán)路設計,并根據(jù)無線雙向偽距測量結(jié)果,通過時延處理模塊實時解算鐘差。FPGA程序設計主要通過時鐘計數(shù)與直接數(shù)字式頻率合成器(Direct Digital Synthesizer,DDS)計數(shù)相結(jié)合的方式實現(xiàn)鐘差修正。具體的實現(xiàn)流程圖如圖6所示。

qrs3-t6.gif

    時間信號同步設計單元一方面完成本地時間與偽衛(wèi)星發(fā)射信號同步,另一方面根據(jù)時延數(shù)據(jù)修正模塊的結(jié)果調(diào)整本地時間與主站時間同步。FPGA發(fā)射端接收到DSP通過解調(diào)電文獲得的主站時間與啟動狀態(tài)標志位后,開始本地時間計數(shù),并在數(shù)據(jù)第一幀時刻啟動偽碼和數(shù)據(jù)碼生成,通過移位器時延修正,實現(xiàn)本地發(fā)射信號與本地時間同步。

    DSP將主站與從站鐘差測量結(jié)果轉(zhuǎn)換成以毫秒計數(shù)與偽碼碼片計數(shù)的粗時延修正值,F(xiàn)PGA通過時鐘計數(shù)與FIFO移位設計實現(xiàn)粗時延修正。將小于一個碼片的時延修正值轉(zhuǎn)換為偽碼數(shù)字控制振蕩器(Numerically Controlled Oscillator,NCO)相位的細時延修正值,F(xiàn)PGA通過DDS技術(shù)調(diào)整NCO相位累加器的方式將時延結(jié)果作為補償值修正從站本地時間與本地偽碼和數(shù)據(jù)碼信號。根據(jù)修正后的本地時間,輸出PPS秒脈沖,同時在整數(shù)秒開始時刻產(chǎn)生同步的從站偽衛(wèi)星發(fā)射信號,實現(xiàn)本地時間與主站時間同步。

    鐘差測量的精確度決定了系統(tǒng)的同步精度,對鐘差的修正誤差會增加系統(tǒng)修正誤差,在保證鐘差精度的同時,應當盡可能地減小系統(tǒng)修正誤差。直接采用時鐘計數(shù)方法對鐘差進行修正,修正的系統(tǒng)誤差為1個時鐘周期,62 MHz系統(tǒng)工作時鐘頻率的修正誤差約為16 ns。本文采用FPGA時鐘計數(shù)與DDS計數(shù)相結(jié)合的方式,可以實現(xiàn)高精度的時延修正,取相位累加器位數(shù)N為32時的時間修正分辨率為:

    qrs3-gs7.gif

    由于DSP接收端跟蹤環(huán)路噪聲的影響,輸出的PPS具有一定的抖動性。采用以最小均方誤差為準則的Kalman濾波算法,通過測量秒脈沖的時間間隔,用前一時刻估計值與當前時刻測量值來估計得到濾波后PPS秒脈沖,實現(xiàn)對秒信號抖動的處理[6]

3.2 載波同步技術(shù)設計

    由于偽衛(wèi)星的晶體振蕩器存在頻率漂移與準確度偏差等原因,在上下變頻模塊實現(xiàn)頻譜搬移過程中會發(fā)生頻偏,會直接影響用戶接收機利用載波相位測量值定位。接收端信號跟蹤環(huán)路以閉環(huán)反饋的形式實現(xiàn)對接收信號的鎖定,本文在從站偽衛(wèi)星DSP接收端跟蹤通道載波環(huán)中加入發(fā)射端載波閉環(huán)修正環(huán)路,實現(xiàn)主站發(fā)射端載波頻率與從站發(fā)射端載波同步。具體的環(huán)路設計如圖7所示。

qrs3-t7.gif

    FPGA混頻模塊將數(shù)字中頻信號sIF(n)與本地正弦載波信號混頻得到同相支路(I支路)結(jié)果,與本地余弦載波信號混頻得到正交支路(Q支路)結(jié)果,兩路結(jié)果通過相關(guān)積分模塊與本地偽碼相關(guān)累加得到1 ms的積分結(jié)果Ip(n)、Qp(n)。

    DSP載波環(huán)路整體設計采用二階鎖頻環(huán)輔助三階鎖相環(huán)結(jié)構(gòu)。鎖頻環(huán)調(diào)整本地復制的載波頻率與接收信號的載波頻率達到一致,鎖相環(huán)調(diào)整環(huán)路輸出信號的相位,使其與輸入信號的相位保持一致,主要區(qū)別在于鑒別器的不同。本文DSP程序設計采用的鑒頻方法為符號函數(shù)sign(·)鑒頻,此方法對數(shù)據(jù)比特跳變不敏感且計算量較小,其計算公式如下:

qrs3-gs8-10.gif

    當鎖相環(huán)鎖定信號后,相位差異基本在零值晃動。FPGA的數(shù)控振蕩器根據(jù)本地載波初始頻率控制字與環(huán)路鑒相結(jié)果通過正弦和余弦函數(shù)查詢表得到與輸入載波同步的本地載波信號。

    FPGA載波積分器累加本地復制的主站與從站載波相位值,通過對載波相位積分結(jié)果作差可以消除本地頻偏,得到與主站相差的載波頻率與相位值,以10 ms的更新速率閉環(huán)修正從站發(fā)射端載波信號,從而實現(xiàn)主站與從站發(fā)射信號載波同步。

    載波同步對比結(jié)果如圖8所示,其中圖8(a)為載波未同步修正下用戶接收機跟蹤主站與從站載波相位作差的結(jié)果;圖8(b)為從站載波同步修正后的載波相位差結(jié)果;圖8(c)為主站與從站同源情況下的載波相位差結(jié)果。

qrs3-t8.gif

4 測試結(jié)果

    本文設計的偽衛(wèi)星系統(tǒng)采用SCPA結(jié)構(gòu),通過雙向測距與時間同步技術(shù)實現(xiàn)偽衛(wèi)星自組織網(wǎng)絡時間同步,通過接收端與從站發(fā)射端載波閉環(huán)修正實現(xiàn)主站載波頻率與從站載波頻率同步。通過系統(tǒng)測試,載波同步結(jié)果如圖8所示,載波同步修正后的結(jié)果與同源情況下的結(jié)果相近,誤差小于0.1 Hz。圖9為主站與從站偽衛(wèi)星雙向測距零時延與50 m距離測試結(jié)果,零時延測試結(jié)果均方根誤差為5.64 cm;無線測試條件下50 m測距結(jié)果均方根誤差為9.26 cm。偽衛(wèi)星主站與從站的時間同步結(jié)果圖10所示(通道2為主站秒脈沖,通道3為從站秒脈沖),時間同步精度優(yōu)于2 ns,可達到系統(tǒng)設計要求。

qrs3-t9.gif

qrs3-t10.gif

5 結(jié)束語

    為解決偽衛(wèi)星自組織網(wǎng)絡時頻同步問題,本文設計了一種基于SCPA結(jié)構(gòu)的偽衛(wèi)星時頻同步系統(tǒng)。在集成DSP+FPGA與上下變頻硬件平臺上完成設計。該系統(tǒng)偽衛(wèi)星星間載波頻率同步優(yōu)于0.1 Hz,時間同步精度優(yōu)于2 ns,結(jié)合載波相位定位終端,可實現(xiàn)厘米級定位精度。該系統(tǒng)可應用于區(qū)域內(nèi)導航與授時、地下定位、室內(nèi)定位等。

參考文獻

[1] 吳剛,劉銀年,王建宇,等.偽衛(wèi)星時鐘同步方法的研究[J].光纖與電纜及其應用技術(shù),2007(2):25-28.

[2] 韓雅蘭,黃智剛,趙昀.基于雙向測距的偽衛(wèi)星網(wǎng)絡高精度時鐘同步方法研究[J].遙測遙控,2010,31(3):8-11.

[3] 李夢,馬紅皎.載波相位平滑偽距算法在雙向測距與時間同步系統(tǒng)中的應用[J].時間頻率學報,2014,37(4):213-220.

[4] MISRA P,ENGE P.Global positioning system signals measurements,and performance[M].Ganga-Jamuna Press,2001.

[5] 鐘興旺,陳豪.雙向單程距離與時差測量系統(tǒng)及零值標定方法[J].電子測量與儀器學報,2009,23(4):13-17.

[6] 范文晶,王召利,周渭.改進卡爾曼濾波算法在時鐘馴服技術(shù)中的應用[J].宇航計測技術(shù),2016,36(3):46-49.



作者信息:

紀元法1,2,梁  濤1,2,孫希延1,2,嚴素清1,2,盧偉軍1,2

(1.桂林電子科技大學 信息與通信學院,廣西 桂林541004;2.廣西精密導航技術(shù)與應用重點實驗室,廣西 桂林541004)

此內(nèi)容為AET網(wǎng)站原創(chuàng),未經(jīng)授權(quán)禁止轉(zhuǎn)載。
亚洲一区二区欧美_亚洲丝袜一区_99re亚洲国产精品_日韩亚洲一区二区
国产精品扒开腿爽爽爽视频| 在线日韩av永久免费观看| 久久成人18免费网站| 亚洲一区美女视频在线观看免费| 91久久精品国产91久久性色tv| 性8sex亚洲区入口| 亚洲欧美日韩中文播放| 亚洲视频在线一区观看| 99re66热这里只有精品3直播| 91久久精品国产| 亚洲电影第三页| 一区二区在线看| 精品999在线播放| 狠狠色香婷婷久久亚洲精品| 韩国v欧美v日本v亚洲v| 国产亚洲精品久久久久婷婷瑜伽| 国产精品专区第二| 国产精品一区二区欧美| 国产伦精品一区二区三区四区免费 | 日韩视频在线免费| 最新精品在线| 亚洲剧情一区二区| 一区二区欧美亚洲| 亚洲永久免费观看| 欧美亚洲日本一区| 欧美在线一二三四区| 久久久久国产精品一区| 麻豆91精品91久久久的内涵| 久久亚洲国产精品一区二区| 噜噜噜久久亚洲精品国产品小说| 美女精品在线观看| 欧美激情视频一区二区三区在线播放| 欧美成人午夜激情| 欧美日本国产| 欧美四级伦理在线| 国产精品入口66mio| 国产欧美一区二区精品仙草咪| 国产性天天综合网| 亚洲大胆在线| 99人久久精品视频最新地址| 亚洲天堂av在线免费| 亚洲欧美经典视频| 亚洲高清不卡在线观看| 99精品国产99久久久久久福利| 亚洲自拍另类| 久久久福利视频| 欧美成人视屏| 国产精品久久激情| 激情久久五月| 亚洲毛片在线| 午夜一区二区三区不卡视频| 久久国产精品一区二区三区四区| 亚洲免费激情| 欧美亚洲尤物久久| 女人色偷偷aa久久天堂| 国产精品国产三级国产普通话三级| 国产一区二区激情| 91久久精品一区| 亚洲影院在线| 亚洲福利视频网站| 中文国产成人精品久久一| 久久精品国产免费观看| 欧美激情中文字幕乱码免费| 国产精品日韩欧美综合 | 国内外成人免费激情在线视频| 91久久精品美女高潮| 亚洲一区二区毛片| 亚洲精品永久免费| 亚洲欧美日韩在线观看a三区| 蜜桃av综合| 国产精品一二一区| 亚洲精品一区在线观看香蕉| 欧美在线高清视频| 亚洲一区自拍| 欧美gay视频| 国产亚洲视频在线| 99精品99| 亚洲精品久久7777| 久久精品国产亚洲高清剧情介绍| 欧美日韩亚洲国产一区| 在线观看视频一区| 亚洲综合视频1区| 9l国产精品久久久久麻豆| 久久精品三级| 国产精品高潮粉嫩av| 亚洲高清自拍| 欧美在线播放高清精品| 亚洲欧美精品suv| 欧美另类变人与禽xxxxx| 国语自产精品视频在线看| 一区二区三区欧美视频| 亚洲欧洲精品一区二区三区| 久久精品国产亚洲高清剧情介绍| 国产精品高潮视频| 999亚洲国产精| 日韩视频免费观看| 美女视频网站黄色亚洲| 国产麻豆视频精品| 亚洲图片欧美一区| 99视频在线精品国自产拍免费观看| 久久久国产精品亚洲一区| 国产精品久久一卡二卡| 一区二区av| 一区二区三区回区在观看免费视频| 另类成人小视频在线| 国产亚洲精品美女| 午夜欧美大尺度福利影院在线看| 中文一区二区在线观看| 欧美理论视频| 亚洲人久久久| 亚洲精品影视| 欧美ed2k| 亚洲电影第1页| 亚洲欧洲日韩综合二区| 美国十次成人| 一区二区在线观看视频| 亚洲电影在线观看| 久久亚洲精品欧美| 激情六月婷婷久久| 亚洲国产精品电影| 免费成人性网站| 136国产福利精品导航| 久久精品天堂| 久久综合中文| 在线观看91精品国产麻豆| 亚洲国产精品日韩| 欧美成人激情视频免费观看| 1000部精品久久久久久久久| 亚洲高清视频一区二区| 久久免费黄色| 激情久久中文字幕| 亚洲人www| 欧美国产视频在线| 亚洲美女av黄| 一区二区三区视频在线看| 欧美日韩一区二区三区四区五区| 夜夜嗨av一区二区三区| 亚洲一区二区三区影院| 国产精品美女在线观看| 午夜精品久久久久久久蜜桃app| 久久精品72免费观看| 在线观看国产一区二区| 亚洲精品一区二区三区婷婷月| 欧美激情二区三区| 一本色道久久综合亚洲精品不| 亚洲女性裸体视频| 国产欧美高清| 亚洲国产精品专区久久| 欧美国产日韩在线观看| 一本一本大道香蕉久在线精品| 午夜精品在线观看| 国产一区二区三区视频在线观看 | 国产毛片一区二区| 欧美中文字幕视频| 欧美成人精品激情在线观看| 日韩视频精品在线观看| 午夜在线a亚洲v天堂网2018| 国内成+人亚洲| 亚洲精品影院在线观看| 国产精品久久久久影院亚瑟| 久久国产婷婷国产香蕉| 欧美精品 国产精品| 在线视频日韩| 久久夜色精品一区| 亚洲伦理一区| 久久精品视频一| 亚洲精品一二区| 欧美一区二区三区在线看| 亚洲成人在线| 午夜宅男久久久| **欧美日韩vr在线| 亚洲欧美另类中文字幕| 激情久久久久久久| 亚洲一区三区电影在线观看| 黄色工厂这里只有精品| 这里是久久伊人| 狠狠入ady亚洲精品经典电影| 中文无字幕一区二区三区| 国产午夜精品一区二区三区欧美| 亚洲精品久久久久久久久久久| 国产精品久久久久9999高清| 久久精品亚洲精品| 欧美日韩一区在线观看视频| 久久成年人视频| 欧美午夜激情在线| 亚洲国产日韩欧美在线图片| 国产精品久久久亚洲一区 | 裸体一区二区| 亚洲影院高清在线| 欧美另类女人| 久久精品国产亚洲精品| 国产精品国产三级国产普通话三级| 亚洲黄色免费| 国产欧美精品日韩| 亚洲视频一起| 亚洲国产精品嫩草影院| 久久精品国产亚洲一区二区三区 | 亚洲视频久久| 欧美国产精品专区| 欧美一级大片在线免费观看|