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

        ?

        空間交會對接多體制接收機設計與實現(xiàn)

        2014-06-24 14:36:18程慶林吳毅杰石云墀
        載人航天 2014年1期
        關鍵詞:空空環(huán)路航天器

        程慶林,洪 亮,吳毅杰,石云墀

        (上海航天電子技術研究所,上海201109)

        空間交會對接多體制接收機設計與實現(xiàn)

        程慶林,洪 亮,吳毅杰,石云墀

        (上海航天電子技術研究所,上海201109)

        空間交會對接時,空空通信系統(tǒng)需在兩個航天器之間建立雙向通信鏈路,完成兩者間的信息交換。在載人航天交會對接任務中,空空通信系統(tǒng)擔負交會對接段以及撤離段航天器間的雙向鏈路數(shù)據(jù)傳輸任務。針對不同的任務模式、分析不同信號的解調(diào)模型和算法特點,基于軟件無線電技術確定本系統(tǒng)的全數(shù)字化解調(diào)算法和硬件平臺設計方案,實現(xiàn)過程運用復用/優(yōu)化設計思路,節(jié)省了邏輯資源、降低了算法復雜度。對關鍵指標的測試表明,該設計可滿足空空通信系統(tǒng)的指標需求。

        多體制;數(shù)字接收機;空空通信系統(tǒng);交會對接

        1 引言

        航天器空間交會對接是空間交會與空間對接的總稱,是兩個航天器之間的在軌運動行為。通常稱在軌運行的航天器為目標航天器,與目標航天器進行交會對接的航天器為追蹤航天器??湛胀ㄐ旁O備在執(zhí)行任務A時,追蹤航天器與目標航天器之間雙向數(shù)據(jù)傳輸速率較低,采用DS?DB?PSK(Direct Sequence Spread Spectrum?Differential Encode Binary Phase?Shift Keying)直接序列擴頻方式進行通信;空空通信設備任務模式為任務B時,追蹤航天器與目標航天器之間雙向數(shù)據(jù)傳輸速率較高,追蹤航天器與目標航天器之間空空通信采用非擴頻的差分二相相移鍵控(Differentially Encoded Binary Phase?Shift Keying,DBPSK)和差分正交相移鍵控(Differentially Encoded Quadra?ture Phase?Shift Keying,DQPSK)進行通信。本文從硬件平臺通用化、軟件功能模塊化的角度出發(fā),在基于“FPGA(Field Programmable Gate Array)、ADC(Analog Digital Converter)”架構的硬件平臺上設計并實現(xiàn)了一種具有三種解調(diào)模式的數(shù)字接收機。在繼承二期空空通信接口和通信指標的基礎上擴展了接收機的DQPSK解調(diào)功能,提高了應用廣度和可靠性。

        2 需求分析

        2.1 工作模式需求

        1)在執(zhí)行任務A時,數(shù)據(jù)通信鏈路繼承二期空空通信的通信模式。整個過程中空空通信鏈路數(shù)據(jù)率較低,綜合考慮發(fā)射功率和捕獲跟蹤作用距離,應采用擴頻方式進行通信。

        擴頻模式下,直接序列擴頻體制將原始數(shù)據(jù)與遠高于其碼速率的偽噪聲序列模二加后生成新序列,可明顯擴展信號頻譜。接收時,利用偽噪聲序列的自相關特性能將信號能量集中在窄帶范圍內(nèi),提高接收信號的信噪比,獲得較高的處理增益[1]。

        2)在執(zhí)行任務B時,數(shù)據(jù)通信鏈路需要在任務A的基礎上增加大量的數(shù)據(jù),雙向通道數(shù)據(jù)量較大,應采用非擴頻的模式進行通信。

        非擴頻模式下,相移鍵控(Phase?Shift Keying,PSK)調(diào)制方式具有恒定的包絡、在給定信號電平中有最低誤碼率、高效率等特點而被廣泛應用在軍事和商業(yè)通信系統(tǒng)中[2]。QPSK調(diào)制體制在相同符號速率情況下碼速率是BPSK調(diào)制的2倍,擁有更高的頻帶利用率[3]。針對雙向通信鏈路需求同時繼承二期的通信需求保留非擴頻DBPSK調(diào)制體制,新增DQPSK調(diào)制體制。

        表1給出了空空通信設備在不同任務模式下的數(shù)據(jù)率和工作模式。

        表1 不同任務模式下工作模式和數(shù)據(jù)率需求Table 1 The modulation mode and bit rate requirementat different stages

        2.2 調(diào)制模式和碼速率

        空空通信設備具有擴頻與非擴頻兩種通信模式,表2給出了兩種通信體制下空空通信機中的傳輸速率及調(diào)制方式。

        表2 空空通信機碼速率及調(diào)制方式Table 2 The bit rate and modulation mode of space to space communicator

        3 多體制接收機的設計

        空空通信機是空間交會對接任務中的核心設備,它的內(nèi)部結構如圖1所示。多體制中頻接收機作為空空通信機中的關鍵部件,需要把天線接收到的信號放大、數(shù)字下變頻、解擴、解調(diào)、RS解碼后送給空空接口[1]。

        圖1 空空通信機內(nèi)部組成Fig.1 The structure of space to space communicator

        多體制接收機中的關鍵技術是中頻信號的解擴、解調(diào),該功能由中頻解擴/解調(diào)模塊實現(xiàn)[1]??紤]到追蹤航天器和目標航天器的通信鏈路中的碼速率、擴頻比不同,若要采用軟件無線電技術實現(xiàn)三種不同的調(diào)制體制的中頻信號的解擴/解調(diào)功能,從平臺通用化的角度出發(fā)可基于相同的硬件平臺加載不同的軟件來實現(xiàn);從軟件功能模塊化的角度出發(fā)可對功能相同的模塊進行復用設計,減少邏輯器件的資源損耗;從工程繼承性角度出發(fā)可涵蓋之前的空空通信指標和接口,增強可靠性。

        多體制中頻解擴/解調(diào)模塊的主要技術指標:中頻7X MHz;采樣率DS_DBPSK:45.X MHz/56. X MHz,DQPSK:45.X MHz,DBPSK:56.X MHz;載波最大多普勒頻偏15 kHz。

        4 多體制接收機的實現(xiàn)

        4.1 硬件設計

        系統(tǒng)的應用目標包含航天產(chǎn)品,對可靠性的要求非常嚴格。減少系統(tǒng)使用的器件種類、降低系統(tǒng)的復雜程度,能夠有效的提高系統(tǒng)的可靠性。綜合以上因素,選擇使用“ADC+FPGA”的方案進行設計,系統(tǒng)的硬件框圖如圖2所示。

        圖2 系統(tǒng)硬件框圖Fig.2 Block diagram of the system hardware

        下面對硬件平臺主要構成進行說明。

        4.1.1 模/數(shù)轉換器

        模/數(shù)轉換器(ADC)使用的是ADI公司生產(chǎn)的AD6645,該芯片采用單電源(+3.3 V)供電,內(nèi)置基準和采樣保持電路,轉換位數(shù)14 bit,最高工作頻率可達200 MHz,最大轉換速率105 MHz,輸入的信號帶寬最大可達3 GHz,輸出的數(shù)據(jù)格式可以設置為二進制補碼或偏移二進制格式。接口形式簡單,只需提供一個采樣時鐘,14位的并行數(shù)據(jù)即可按時間順序依次輸出[4]。ADC的輸入采用差分低壓正發(fā)射極耦合邏輯(Low Voltage Positive Emitter?Couple Logic,LVPECL)電平,AD轉換數(shù)據(jù)可以截取N bit(N≤14)或者全部選用14 bit作為有效數(shù)據(jù),編碼形式為二進制補碼。

        4.1.2 FPGA

        FPGA使用的是Xilinx公司的300萬門FP?GA?XC2V3000,內(nèi)含14436個Slice、1728 Kbit的RAM和96個18 bit×18 bit的乘法器[5],完全能夠滿足應用的需求。解擴/解調(diào)過程中的偽碼捕獲、偽碼跟蹤、載波捕獲、載波跟蹤、數(shù)字變頻、數(shù)字濾波、數(shù)字基帶處理以及后續(xù)的RS解碼等功能均由其實現(xiàn)。

        4.1.3 時鐘源

        由于解擴/解調(diào)系統(tǒng)的工作時鐘需要根據(jù)不同的解調(diào)模式發(fā)生改變。為了提高系硬件平臺的靈活性和通用化程度,配置兩片頻率分別為45.X MHz和56.X MHz晶振,必要時可配合DCM得到不同的時鐘資源以滿足系統(tǒng)需求。

        4.2 軟件設計

        FPGA軟件是接收機的核心,由并行相關器、偽碼發(fā)生器、下變頻器、鑒相器、環(huán)路濾波器、RS解碼解幀等模塊組成。每個功能模塊在數(shù)字域都有其相對應的模型,因此所有的功能部件都工作在數(shù)字域,采用Verilog HDL語言編寫。

        4.2.1 擴頻模式

        擴頻模式下多體制接收機需要對發(fā)送來的數(shù)據(jù)進行解擴、BPSK解調(diào)、RS解碼、解幀后輸出原始信息給空空接口[1]。根據(jù)以上要求,擴頻模式軟件設計框圖如圖3所示。

        1)偽碼的捕獲和跟蹤

        直接擴頻信號解擴的關鍵是保證本地PN序列和與接收信號同步。通過與本地PN序列進行相關運算恢復出原始數(shù)據(jù)。偽碼的同步分捕獲和跟蹤兩個階段:通過擴頻碼的捕獲可使本地偽碼與接收到的碼元基本保持同步,獲得一定的同步精度(如1/2碼元寬)。由于發(fā)射機和接收機的相對運動及時鐘的不穩(wěn)定,特別是在擴頻模式下,須進一步跟蹤接收信號,使本地偽碼盡可能跟隨接收信號的變化,將同步精度控制在更小碼元周期范圍內(nèi),進一步降低數(shù)據(jù)誤碼率。綜合考慮算法復雜度和邏輯資源偽碼捕獲所采用的多通道相關捕獲法,偽碼跟蹤階段采用延遲鎖定環(huán),門限判決階段采用基于瞬時標定功率的自適應算法[1]。

        2)載波的捕獲和跟蹤

        在進行載波捕獲時,為了能夠達到最小的捕獲時間,應該先搜索最大可能性的頻率點,然后再搜索次可能性的頻率點。在不考慮多普勒漂移的情況下,中頻輸入信號的頻率點為18.X MHz,所以這就是可能性最大的頻率點,搜索應該圍繞這個頻率點展開。搜索策略如圖4所示。

        假設頻率搜索帶寬為f,頻率搜索步進為Δf,那么需要在f/Δf個頻率點上進行搜索。由于采用并行相關,因而調(diào)整1次碼相位可以完成多個碼相位位置的相關計算。設在一個頻率點上需要調(diào)整m次碼相位,即m次碼相位調(diào)整可以完成一個完整PN序列所有碼相位的遍歷;并且設一次碼相位調(diào)整需要相位計算時間T才能完成計算,則捕獲時間為mTf/Δf。

        載波跟蹤階段采用經(jīng)典的Costas環(huán),包括下變頻,低通濾波,鑒相器,環(huán)路濾波和壓控振蕩器[6]。由于這一部分與非擴模式下DBPSK解調(diào)的功能相似,考慮對該部分模塊功能進行復用。

        4.2.2 非擴頻模式

        BPSK信號與QPSK信號解調(diào)的主要區(qū)別:

        在載波同步上,兩者的鑒相算法完全不同;在基帶處理上,兩者的編碼特點不同,使得各自的差分解碼在實現(xiàn)時也不同。其次,因為BPSK信號攜帶信息的只有I路,不需要并串轉換模塊,而QPSK信號I、Q兩路都攜帶相應的碼元信息,需要對I、Q路數(shù)據(jù)信息進行串行組合。

        1)鑒相算法

        經(jīng)典的Costas環(huán)的中的鑒相器是一個簡單的乘法器,該環(huán)路被廣泛應用在各種數(shù)字通信解調(diào)接收系統(tǒng)中。但是經(jīng)典Costas環(huán)只能解調(diào)BPSK信號和UQPSK信號,要實現(xiàn)QPSK解調(diào)必須考慮使用其他的環(huán)路。將經(jīng)典的Costas環(huán)進行推廣得到可以用于QPSK載波恢復的改進型硬判決Costas環(huán)[2],其中的鑒相部分就是硬判決鑒相器,其鑒相算法實現(xiàn)框圖如圖5所示。

        圖3 擴頻模式軟件設計框圖Fig.3 The FPGA structure of DSSS demodulation

        圖4 擴頻模式下載波搜索策略Fig.4 The carrier searching strategy of DSSS demodulation

        圖5 硬判決鑒相器算法實現(xiàn)Fig.5 The implementation of hard-decision phase detecting algorithm

        其中sgn()表示符號函數(shù),表3給出了兩種鎖相環(huán)路中鑒相器的鑒相算法。

        表中:Ud表示鑒相誤差,I、Q表示經(jīng)低通濾波后的同相、正交信號。Δω、Δφ分別表示輸入輸出信號的頻差和相差。當環(huán)路鎖定時Δω=0,I(k)2=1、Q(k)2=1,Δφ很小。

        表3 兩種鑒相器的鑒相算法Table 3 The phase detection algorithm of two phase detectors

        2)數(shù)字下變頻器

        下變頻器主要由三個部分組成:本地振蕩器,數(shù)字混頻器,低通濾波器。

        數(shù)控振蕩器(Numerically Controlled Oscillator,NCO)的原理與直接數(shù)字頻率合成的原理相同,也是從相位概念出發(fā)直接合成所需的正/余弦波形[2]。將采樣頻率fclk均分為2N份,輸出頻率是頻率分辨率的s倍,s是頻率控制字(Frequency Control Word,F(xiàn)CW),則有式(1):

        數(shù)字混頻器完成采樣的信號與NCO產(chǎn)生的正/余弦信號相乘;低通濾波器濾除其中的二次倍頻分量。

        環(huán)路濾波器有很多類型,常見的主要包括RC濾波器、無源比例積分濾波器和有源比例積分濾波器[7],這里選用有源比例積分濾波器。典型的一階模擬環(huán)路濾波器的傳遞函數(shù)如式(2)、(3)所示:

        圖6 環(huán)路濾波器基本結構Fig.6 The basic structure of loop filter

        利用雙線性變換,得到環(huán)路濾波器的最終計算公式如式(4)與(5)。

        從式(4)和式(5)可以看出,在一個采樣時間間隔后,可以由給定的環(huán)路自然角頻率ωn和環(huán)路阻尼因子ξ來共同確定環(huán)路濾波器的系數(shù)。

        4.2.3 復用/優(yōu)化設計思想

        非擴頻模式下,考慮到BPSK信號與QPSK信號解調(diào)模型中很多的功能模塊的結構相同,在軟件設計上考慮復用設計。

        1)低通濾波器

        要實現(xiàn)多模數(shù)字解調(diào)功能,最直接的思路是設計多個低通濾波器。通過后續(xù)的工程實踐發(fā)現(xiàn):在所有功能模塊的設計中,濾波器耗費邏輯資源最多。本文采用的方案是結合Matlab FDATool和FIR濾波器IP core對低通濾波器進行復用設計或者融合設計。當兩種解調(diào)模式中的時鐘和采樣率都是碼速率的某一整數(shù)倍時,兩種調(diào)制模式可以考慮完全復用;當采樣率不是碼速率的同一整數(shù)倍時,在設計低通濾波器時可以考慮融合設計。僅設計一個低通濾波器時,先利用Matlab FDATool設計階數(shù)相同的濾波器系數(shù),再通過模式選擇相應的系數(shù)組來實現(xiàn)預期的濾波功能;ISE中提供的用來實現(xiàn)FIR濾波器的IP core(fir_compiler V3.2)最高的采樣時鐘可達250 MHz;如此設計的濾波器還可同時完成整數(shù)倍插值和抽取,內(nèi)部可實現(xiàn)多通道配置,完全能滿足設計需要。

        2)壓控振蕩器

        設計一個壓控振蕩器,通過模式選擇信號對NCO的三種模式下的頻率控制字做選擇,而非單獨設計。ISE中提供的用來實現(xiàn)DDS的IP core(dds_compiler v2.1)最高工作時鐘可達550 M;輸出的正/余弦信號動態(tài)范圍為:18~120 dB,對應的有效位數(shù)為3~21 bit;頻率分辨率范圍為:0.0233 Hz~10 kHz;兩者都能在其范圍內(nèi)實現(xiàn)自主配置以適應不同的需求;內(nèi)部提供頻率偏移和相位偏移的可編程設計,具有較高的通用性。

        (3)環(huán)路濾波器

        將環(huán)路濾波器中的兩個系數(shù)的乘法運算優(yōu)化為數(shù)據(jù)移位運算。圖7說明了環(huán)路濾波器在軟件中的實現(xiàn)過程。在復用過程中,將LF的輸入信號通過上位機選擇不同的移位位數(shù)代替各自的相乘功能,累加器可以通過一個加法器和一個D觸發(fā)器來實現(xiàn)。

        圖7 環(huán)路濾波器的優(yōu)化設計Fig.7 The optimization of loop filter

        其中K1≈log2C1、K2≈log2C2

        圖8表示的是非擴頻模式下兩種PSK信號的解調(diào)軟件的實現(xiàn)過程,其中的相關部件考慮復用或者優(yōu)化處理。

        虛線框內(nèi)的部分表示的數(shù)字下變頻部分,三種模式下該部分功能完全復用。上位機軟件通過模式選擇信號通知NCO、低通濾波器、鑒相器、環(huán)路濾波器等部件選擇對應的參數(shù)組或相對應的功能模式完成非擴模式下的兩種PSK信號解調(diào)。

        4.2.4 設計思路對比

        如果采用傳統(tǒng)的設計思路,即低通濾波器、壓控振蕩器、環(huán)路濾波器等都采用采用獨立設計,并且在算法上不做任何優(yōu)化設計,純源碼實現(xiàn)占據(jù)300萬門FPGA邏輯資源的71%;而采用模塊復用+優(yōu)化設計思路,純源碼實現(xiàn)只占FPGA的46%。

        4.2.5 在線測試與分析

        主要就以下三個指標對多體制中頻接收機進行整機測試:

        圖8 非擴頻模式軟件設計框圖Fig.8 The FPGA structure of DBPSK/DQPSK demodulation

        1)可捕獲跟蹤的中頻信號頻率范圍測試;

        2)解調(diào)系統(tǒng)的抗干擾能力測試:在中頻信號中混入不同信噪比的加性白噪聲,在保證一定誤碼率水平的情況下確定解調(diào)系統(tǒng)所能夠承受的信噪比下限;

        3)接收的中頻動態(tài)范圍即無噪聲中頻信號的功率范圍測試。

        表4、5、6分別給出了三種性能的測試結果。

        表4 不同模式/參數(shù)下的中頻頻率捕獲跟蹤范圍Table 4 The IF frequency capture&trace range under different modes/parameters

        表5 不同模式/參數(shù)下的解調(diào)損耗Table 5 The demodulation loss under differentmode/parameter

        表6 不同模式/參數(shù)下的中頻接收動態(tài)范圍Table 6 The IF receiving dynamic range under different modes/parameters

        分析三表測試結果,可以發(fā)現(xiàn)在擴頻模式下,中頻頻率捕獲跟蹤范圍大于13 kHz,非擴頻模式下,中頻頻率捕獲跟蹤范圍都維持在符號速率的±3%以上。對于空間通信來說,在同時考慮到多普勒頻移和發(fā)射機發(fā)射頻率精度的情況下,該測試結果有足夠大的裕量;抗干擾的測試結果中,沒有引入RS解碼時,解調(diào)損耗都在3.5 dB以下,擴頻模式由于解擴帶來了一定的損耗;如果引入RS解碼模塊,解調(diào)損耗將進一步減小1.5~2 dB。對于中頻接收動態(tài)范圍測試,接收信號功率過大,會引起AD過幅,導致其有效位數(shù)降低;信號功率過小,信號無法被檢測到,動態(tài)范圍就是指這個最大最小的范圍。影響中頻接收動態(tài)范圍的主要有:ADC的量化精度、各功能模塊增益分配等。如果接收機前端引入自動增益控制模塊,動態(tài)范圍能得到進一步的擴展。

        5 結論

        多體制接收機作為空空通信子系統(tǒng)中的核心部件,將在載人航天新一期的交會對接任務中扮演關鍵角色。文章基于軟件無線電技術結合型號的現(xiàn)實需求,研究并實現(xiàn)了一種模式可切換的全數(shù)字化多體制中頻接收機。從功能模塊化、平臺通用化和具有繼承性的角度出發(fā),在常規(guī)設計方法的基礎上提出了復用/優(yōu)化設計方法,減少了FPGA邏輯資源的使用,提高了算法的可移植性,豐富了接收機的解調(diào)功能。還重點對接收機的中頻頻率跟蹤性能、抗噪聲性能和中頻接收動態(tài)范圍等表征接收機重要性能的指標進行了定量評估和分析,驗證了該全數(shù)字多體制中頻接收機的功能與指標,測試結果表明該多體制接收機能夠滿足型號任務需求。

        [1] 石云墀.載人航天空空子系統(tǒng)及其關鍵技術[J].上海航天,2011,38(6):38?42.

        [2] 洪振宏.高碼速率QPSK解調(diào)器載波恢復環(huán)的算法與實現(xiàn)[D].北京:中國科學院空間科學與應用研究中心,2007.

        [3] 樊昌信,曹麗娜.通信原理[M].第6版.北京:國防工業(yè)出版社,2006.

        [4] Analog Device.AD6645 data sheet[EB/OL].2008.ht?tp://www.analog.com/en/analog?to?digital?converters/ad?converter/ad6645/products/product.html.

        [5] Xilinx Inc.Virtex?ii XC2V3000 data sheet[EB/OL].2001. http://www.xilinx.com/support/index.html/cont ent/xil?inx/en/supportNav/silicon_devices/mature_and_discontinued_products/virtex?ii.html.

        [6] 楊小牛,樓才義,徐建良.軟件無線電技術與應用[M].北京:北京理工大學出版社,2010.

        [7] Floyd M.Gardner.鎖相環(huán)技術[M].姚劍清譯.第3版.北京:人民郵電出版社,2007.

        [8] John G.Proakis.數(shù)字通信[M].張力軍等,譯.第四版.北京:電子工業(yè)出版社,2010.

        The Design and Implementation of Multi?mode Receiver for Rendezvous and Docking in Space

        CHENG Qinglin,HONG Liang,WU Yijie,SHI Yunchi
        (Shanghai Institute of Aerospace Electronic Technology,Shanghai201109,China)

        During the manned spaceflight mission,space to space communication system is responsible for the information exchange between two spacecrafts at the rendezvous and docking and evacuation stages.For different mission stages,different communications modes are required.In order to fulfill the different mission requirements in the same hardware,the feature of various demodulation modes and algorithm was analyzed,and full digital demodulation algorithm and hardware design were determined based on software radio technology.A method of restructure and multiplex was used to optimize the design,which leads to less logical resource and decreasing algorithm complexity.The result of the test for the system indicted that the key characteristic of the design could meet the specification of the space to space communication system.

        multi?modes;digital receiver;space to space communication system;rendezvous and docking

        V443+.1;V526

        文章編號:1674?5825(2014)01?0058?07

        2013?07?01;

        2013?12?27

        程慶林(1987?),男,碩士研究生,研究方向為軟件無線電技術。E?mail:cnhbcql@126.com

        猜你喜歡
        空空環(huán)路航天器
        梁潮詩作《了
        ——如空空,空空如……》·源清書法
        藝術品鑒(2023年13期)2023-08-22 10:33:26
        2022 年第二季度航天器發(fā)射統(tǒng)計
        國際太空(2022年7期)2022-08-16 09:52:50
        2019 年第二季度航天器發(fā)射統(tǒng)計
        國際太空(2019年9期)2019-10-23 01:55:34
        2018 年第三季度航天器發(fā)射統(tǒng)計
        國際太空(2018年12期)2019-01-28 12:53:20
        2018年第二季度航天器發(fā)射統(tǒng)計
        國際太空(2018年9期)2018-10-18 08:51:32
        孫空空 捏面人
        上海市中環(huán)路標線調(diào)整研究
        上海公路(2018年4期)2018-03-21 05:57:46
        腦袋空空
        復制人
        Buck-Boost變換器的環(huán)路補償及仿真
        電測與儀表(2014年8期)2014-04-04 09:19:36
        青青草国内视频在线观看| 精品九九人人做人人爱| 久热国产vs视频在线观看| 久久伊人色av天堂九九| 成年女人A级毛片免| av毛片在线播放网址| va精品人妻一区二区三区| 欧美最猛性xxxx| 亚洲啪av永久无码精品放毛片| 熟妇人妻中文av无码| 无码91 亚洲| 日本在线免费不卡一区二区三区| 日韩精品熟女中文字幕| 欧美成人看片一区二区三区尤物| 午夜福利92国语| 国产主播无套内射一区| 久久国产女同一区二区| 黄片视频大全在线免费播放| 亚洲香蕉成人av网站在线观看| 天堂一区人妻无码| 国产精品欧美韩国日本久久| 亚洲第一页在线免费观看| 精品国产亚洲亚洲国产| 久久久www免费人成精品| 国产成人无码aⅴ片在线观看| 精选二区在线观看视频| 国产白色视频在线观看| 粗大的内捧猛烈进出少妇| 亚洲av纯肉无码精品动漫| 国产av天堂亚洲国产av麻豆| 国产在线一区二区三区不卡| 欧美牲交a欧美牲交aⅴ免费下载 | 亚洲熟妇无码久久精品疯| 精品日本韩国一区二区三区| 国产激情视频在线观看大全| 久久久噜噜噜久久| 亚洲日韩欧美国产另类综合| 精品人妻中文字幕一区二区三区| 日本不卡视频一区二区| 九色综合九色综合色鬼| 国产又黄又大又粗视频|