• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    CPU?GPU平臺上的高速MPSK并行解調(diào)算法

    2019-07-08 05:33:51吳濤閆迪劉燕都趙江
    現(xiàn)代電子技術(shù) 2019年13期

    吳濤 閆迪 劉燕都 趙江

    摘 ?要: 針對多進(jìn)制數(shù)字相位調(diào)制信號特點,采用數(shù)據(jù)并行方法,將連續(xù)信號流分解為多個信號塊以實現(xiàn)并行處理。研究一種基于三維迭代搜索的載波相位?碼元相位聯(lián)合估計算法,在本地構(gòu)造載波信號的搜索庫,利用分段碼元相關(guān)法從搜索庫中挑選出與待處理信號相似度最高的載波信號作為同步載波。通過迭代搜索方法逐步縮小搜索庫,降低計算量。仿真結(jié)果表明,與克拉美?羅下限相比,解調(diào)損失小于0.1 dB。在惠普工作站與英偉達(dá)K20平臺上搭建驗證系統(tǒng),系統(tǒng)的處理速率可達(dá)625.9 MB/s。

    關(guān)鍵詞: CPU?GPU平臺; 高速MPSK; 并行解調(diào); 相位模糊; 三維迭代搜索; 數(shù)據(jù)并行

    中圖分類號: TN911.72?34;V443+.1 ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)識碼: A ? ? ? ? ? ? ? ? 文章編號: 1004?373X(2019)13?0009?06

    A parallel demodulation algorithm based on CPU?GPU platform

    for high?rate MPSK signal

    WU Tao1, YAN Di1, 2, LIU Yandu1, ZHAO Jiang3

    (1. Department of Electronic and Optical Engineering, Space Engineering University, Beijing 101416, China;

    2. The Unit 61646 of PLA, Beijing 100192, China; 3. China Satellite Maritime Tracking and Control Department, Jiangyin 214431, China)

    Abstract: According the characteristics of multiple phase shift keying (MPSK) signal, the data parallel method is used to decompose the continuous signal stream into multiple signal blocks to realize parallel processing. A carrier phase?code element phase joint estimation algorithm based on three?dimensional iterative search is studied. The search library of carrier signal is constructed locally, and then the carrier signal with the highest similarity with the signal under processing is chosen from the search library by means of segmentation code element correlation method to act as the synchronization carrier. The iterative search method is used to gradually deflate the search library to reduce the computational quantity. The simulation results show that, in comparison with the Cramer?Rao lower bound, the demodulation loss of the proposed algorithm is less than 0.1 dB. The verification system was built on HP workstation and NVIDIA K20 GPU, which realized the processing rate of 625.9 MB/s.

    Keywords: CPU?GPU platform; high?rate MPSK; parallel demodulation; phase ambiguity; three?dimensional iterative search; data parallelism

    0 ?引 ?言

    衛(wèi)星通信數(shù)傳正朝著更高碼率、更大帶寬的方向發(fā)展,而高速數(shù)傳解調(diào)是一個重要問題。CPU?GPU平臺計算能力飛速發(fā)展,在其上進(jìn)行數(shù)字信號處理已成為一種發(fā)展趨勢[1?3]。本文針對CPU?GPU異構(gòu)平臺上的高速MPSK信號解調(diào)算法進(jìn)行研究。

    目前文獻(xiàn)公開的三種高速并行解調(diào)架構(gòu)分別是IEC(Interstate Electronics Corporation)[4],PRX(Parallel Receiver)[5]和APRX(All?digital Parallel Receiver)[6]。其中,APRX實現(xiàn)簡單,資源消耗少,被廣泛使用。國內(nèi)在APRX基礎(chǔ)上進(jìn)行了一些數(shù)傳解調(diào)器的研究[7?8]。盡管如此,APRX仍有兩點不足:一是APRX要求采樣率必須為數(shù)據(jù)率整數(shù)倍,為避免采樣混疊,信號中心頻率需隨數(shù)據(jù)率選取[9?10];二是基于APRX的并行解調(diào)器不適合在CPU?GPU平臺上運行,基于APRX的各種解調(diào)算法均存在反饋調(diào)整,用來糾正載波頻差和定時誤差等,在CPU,GPU中實現(xiàn)這種反饋需要頻繁的線程同步,帶來較大的同步開銷,且反饋操作增加了算法內(nèi)部的數(shù)據(jù)依賴,降低了算法的并行性。因此,APRX并不適合于CPU?GPU平臺。

    考慮到CPU?GPU平臺的特殊性,本文采用數(shù)據(jù)并行思想,將高速信號流分解為多個信號塊,通過同時處理多個信號塊實現(xiàn)并行解調(diào)。這要求信號塊的處理必須是獨立的、互不相關(guān)的。一般而言,信號的相干解調(diào)使用反饋環(huán)路方法,如果用于本文的信號塊處理,必然不能滿足要求。當(dāng)信號塊長度足夠短時,可以認(rèn)為在這段時間內(nèi)信號的載波參數(shù)是恒定不變的,因此,本文將信號塊的相干解調(diào)問題轉(zhuǎn)化為載波參數(shù)估計問題,各信號塊的載波參數(shù)估計可以獨立進(jìn)行。首先提出一種基于三維迭代搜索的載波相位?碼元相位聯(lián)合估計算法,迭代搜索空間是<載波頻率,載波初相,碼元相位>。通過不斷縮小搜索空間的每個維度的范圍,大幅減小搜索計算量,因此稱之為迭代搜索。然后研究了算法在CPU?GPU平臺上的具體實現(xiàn)方法。最后在惠普Z820+NVIDIA TESLA K20 硬件平臺上搭建驗證系統(tǒng)。

    1 ?并行解調(diào)算法

    1.1 ?算法原理

    算法利用互相關(guān)的原理,在本地構(gòu)造出大量搜索信號,將搜索信號逐一與接收信號做相關(guān),當(dāng)相關(guān)值最大時,就得到了接收信號的參數(shù)估值。由于有3個待估參數(shù),所以搜索空間是三維的,用[]表示,[fd]是載波多普勒頻率,單位是Hz;[φc]和[φb]分別是載波初相和碼元初相,單位是rad。

    互相關(guān)是反映兩個信號相似程度的常用方法,兩個信號越接近,相關(guān)值越大。以QPSK信號為例,可表示為:

    碼元[a]是隨機分布的,故[J]為0??梢姡捎谡{(diào)制信息影響,即使[f=fc+fd],[φc=φ],[J]依然為0。若按式(5)的方法進(jìn)行相關(guān),[J]就能反映接收信號和本地搜索信號的相似程度。通過尋找最大互相關(guān)值,就能獲得載波頻率和載波初相的估計值。本文稱之為分段碼元相關(guān)法。

    1.2 ?算法流程

    算法流程如圖1所示,基本步驟如下:

    步驟1:初始化搜索空間。搜索空間初始值[0],其中,[i∈1,N1,j∈1,N2,][l∈1,N3],[N1]是[fd]的個數(shù),[N2]是[φc]的個數(shù),[N3]是[φb]的個數(shù),下標(biāo)0表示當(dāng)前迭代次數(shù)為0。

    圖1 ?算法流程圖

    步驟3:求相關(guān)值。利用分段碼元相關(guān)法,求接收信號與搜索信號的相關(guān)值[CORRi,j,l],對于每一種搜索信號有[N3]個不同的碼元初相,故有[N3]個相關(guān)值。

    1) 本次迭代的相關(guān)值小于上次迭代;

    2) 載波頻率的步長小于50 Hz;

    3) 載波初相的步長小于0.01。

    搜索空間的更新方法是:從[p]中選取第[p]次迭代得到的參數(shù)估值的左右相鄰值,并作為[p+1]次迭代搜索空間的邊界值,步長縮小為第[p]次步長的[13]。

    1.3 ?分段碼元相關(guān)法

    分段碼元相關(guān)法在計算互相關(guān)值時,消除調(diào)制信息影響。其原理是在每個碼元周期內(nèi)做相關(guān),得到一系列相關(guān)值,將這些相關(guān)值的絕對值求和得到總相關(guān)值,共分為兩個步驟:

    步驟1:確定碼元分布情況。

    碼元分布由碼元初相和碼元周期決定,[φb]分布在[0,2π]內(nèi),碼元周期由實際碼速率[Rs]決定,[Rs]由碼多普勒[Rsd]和標(biāo)準(zhǔn)碼速率[Rs0]組成:

    若虛擬出一個頻率為[Rs2],初相為[φb2]的正弦波,那么正弦波的過零點就是碼元變化時刻,介于兩個相鄰過零點之間的采樣點是一個碼元周期,如圖2所示。

    圖2 ?碼元分布情況確定方法

    步驟2:計算積分。

    接收信號和本地搜索信號的混頻信號為[mix(k)]。假設(shè)碼元分布情況如下:第一個不完整碼元的采樣點數(shù)量是[x],完整碼元的個數(shù)是[y],每個完整碼元內(nèi)的采樣點數(shù)是[N(i),i∈[1,y]],最后一個不完整碼元的采樣點數(shù)量是[z]。那么,積分值[CORR]如下:

    1.4 ?相位模糊匹配方法

    用鎖相環(huán)方法對QPSK信號進(jìn)行解調(diào)時,由于本地恢復(fù)載波與QPSK載波相位存在多種可能的相位差,從而導(dǎo)致解調(diào)結(jié)果存在模糊[11],本文同樣存在這種模糊現(xiàn)象。假設(shè)本地載波與QPSK載波同頻,只存在相位差[Δφ],正交下變頻后[I],[Q]路如式(9)所示:

    表1 ?不同[Δφ]下QPSK解調(diào)的相位模糊

    每個信號塊獨立解調(diào),其[Δφ]不一致,導(dǎo)致相位模糊也不一致,需要對不同信號塊的相位模糊進(jìn)行匹配。本文提出一種基于數(shù)據(jù)冗余的相位模糊匹配方法,基本原理是利用相鄰信號塊的重疊,調(diào)整后一個信號塊的模糊使之與前一個模塊一致,最終使得所有信號塊的模糊都一致。如圖3所示為信號塊冗余分割方法,[L]為信號塊長度,[L0]為信號塊重疊區(qū)長度。

    圖3 ?信號塊冗余分割方法

    圖4 ?相位模糊匹配方法流程圖

    2 ?仿真與測試

    為了驗證算法解調(diào)性能和實際執(zhí)行效率,本文進(jìn)行Matlab仿真,并在CPU?GPU平臺上對算法進(jìn)行測試。統(tǒng)一使用以下調(diào)制參數(shù):調(diào)制方式為QPSK,符號率[Rs=] 1 Gb/s,載頻頻率[fc=1] GHz,采樣頻率[fs=4] GHz,多普勒頻率[fd=]12 345 Hz,載波初相[φc=0.45] rad,碼元初相[φb=0.5] rad。

    2.1 ?迭代搜索過程的仿真

    設(shè)置信噪比[EbN0=]10 dB,載波多普勒的初始搜索空間是(-20 000,-10 000,0,10 000,20 000),載波相位的初始搜索空間是[π8,π4,3π8,π2],碼元相位的初始搜索空間是[π4,3π4,5π4,7π4]。仿真結(jié)果如圖5所示,圖5a)是載波多普勒估計結(jié)果,圖5b)是載波初相估計結(jié)果,圖5c)是碼元初相估計結(jié)果,圖5d)是相關(guān)值變化??梢?,參數(shù)估值隨著迭代過程逐漸逼近真值,相關(guān)值增加到一定程度后變化很小。

    大量仿真表明,經(jīng)過4次迭代后能完成參數(shù)估計,共進(jìn)行80+64+64+64=272次相關(guān)運算。如果直接進(jìn)行三維搜索,為了達(dá)到與迭代搜索同樣的估計精度,搜索空間步長需要設(shè)置為當(dāng)前的[127],需進(jìn)行(5×27)×(4×27)×(4×27)=1 574 640次相關(guān)運算。可見,迭代搜索方法將相關(guān)運算次數(shù)減少為0.017%,大大減少了計算量。

    圖5 ?仿真結(jié)果

    2.2 ?算法解調(diào)性能仿真

    為了驗證算法的解調(diào)誤碼性能,對[EbN0]在[0,10] dB區(qū)間內(nèi)的信號進(jìn)行解調(diào),統(tǒng)計誤碼率結(jié)果與克拉美?羅下界(Cramer?Rao Lower Bound,CRLB)如圖6所示。從仿真結(jié)果可以看到,在[0,10] dB區(qū)間內(nèi),解調(diào)損失在0.1 dB以內(nèi)。

    2.3 ?CPU?GPU平臺上的測試

    為進(jìn)一步驗證算法的執(zhí)行效率,在CPU?GPU平臺上搭建驗證系統(tǒng),系統(tǒng)硬件配置如表2所示。

    圖6 ?誤碼率仿真結(jié)果比較

    驗證系統(tǒng)由信號調(diào)制機和信號解調(diào)機兩部分組成,二者通過萬兆網(wǎng)連接,信號調(diào)制機產(chǎn)生8 bit量化的QPSK信號后,傳給信號解調(diào)機進(jìn)行處理。信號調(diào)制機采用“多線程計算,輪流輸出”的方法,每個線程產(chǎn)生時間上間隔的QPSK信號塊,按順序輪流讀取各線程產(chǎn)生的信號塊,就得到了時間上連續(xù)的QPSK信號流。

    信號解調(diào)機采用如圖7所示的框架,在CPU中完成信號流分割,將信號塊填充到信號池。進(jìn)程管理器創(chuàng)建多個解調(diào)進(jìn)程,將信號塊傳入GPU進(jìn)行處理。結(jié)果綜合使用相位模糊匹配方法,將不同信號塊的解調(diào)結(jié)果再加工合并為最終解調(diào)結(jié)果。

    表2 ?系統(tǒng)硬件配置參數(shù)

    如圖8所示為運行狀態(tài)的系統(tǒng),100 s內(nèi)信號解調(diào)機的網(wǎng)絡(luò)接收速度如圖9所示,平均處理速率是625.9 MB/s,換算為碼速率后是313.0 Mb/s。

    圖8 ?系統(tǒng)硬件組成

    圖9 ?信號解調(diào)機的網(wǎng)絡(luò)接收速度

    3 ?結(jié) ?論

    本文提出一種CPU?GPU平臺上的高速MPSK解調(diào)算法,給出了算法的具體步驟,并在惠普Z820工作站上進(jìn)行了編程實現(xiàn)。仿真結(jié)果表明,[0,10] dB區(qū)間內(nèi)的解調(diào)損失在0.1 dB以內(nèi),驗證系統(tǒng)實現(xiàn)了313.0 Mb/s的QPSK實時解調(diào)。313.0 Mb/s解調(diào)速度是在一塊NVIDIA TESLA K20上實現(xiàn)的,與FPGA平臺相比并沒有明顯優(yōu)勢。但本文方法的最大優(yōu)勢是其具有可擴(kuò)展性,通過擴(kuò)展GPU計算資源就能方便地增加解調(diào)進(jìn)程數(shù)目,進(jìn)而提高系統(tǒng)解調(diào)速率。

    參考文獻(xiàn)

    [1] 牛炯,申晉晉,陳倩.高頻地波雷達(dá)實時信號并行處理方案及實現(xiàn)[J].中國海洋大學(xué)學(xué)報,2017,47(2):56?61.

    NIU Jiong, SHEN Jinjin, CHEN Qian. Real?time signal parallel processing scheme and realization of HF ground wave radar based on CPU+GPU architecture [J]. Periodical of Ocean University of China, 2017, 47(2): 56?61.

    [2] 王子瀟.基于CUDA GPU平臺的脈沖多普勒雷達(dá)信號模擬軟件[D].北京:北京理工大學(xué),2016:1?3.

    WANG Zixiao. A PD radar simulation software based on CUDA GPU platform [D]. Beijing: Beijing Institute of Technology, 2016: 1?3.

    [3] 武新波.基于CUDA的GPS軟件接收機研究[D].北京:北京理工大學(xué),2015:3?4.

    WU Xinbo. Research on GPS software receiver based on CUDA [D]. Beijing: Beijing Institute of Technology, 2015: 3?4.

    [4] 林長星.2 Gb/s高速通信解調(diào)技術(shù)及其實現(xiàn)研究[D].北京:清華大學(xué),2012:31?33.

    LIN Changxing. Research on demodulation technique and its implementation for 2 Gb/s high speed communication [D]. Beijing: Tsinghua University, 2012: 31?33.

    [5] SRINIVASAN M, CHEN C C, GREBOWSKY G, et al. An all?digital high data?rate parallel receiver [J]. Surface science, 1997, 604(15/16): 1294?1299.

    [6] SADR R, CAIDYANATHAN P, RAPHAELI D, et al. Parallel digital modem using multirate digital filter banks [M]. California: JPL Publication, 1994: 1?112.

    [7] 梁俠,任海根,徐先超,等.800 Mb/s高速解調(diào)器的定時恢復(fù)算法及實現(xiàn)研究[J].現(xiàn)代電子技術(shù),2007,30(23):1?3.

    LIANG Xia, REN Haigen, X? Xianchao, et al. 800 Mb/s timing recovery algorithm and implementation for high speed demodulator [J]. Modern electronics technique, 2007, 30(23): 1?3.

    [8] 郭曉峰,鄭雪峰,盧滿宏,等.高速數(shù)傳QAM解調(diào)器設(shè)計及應(yīng)用[J].遙測遙控,2011,32(3):21?25.

    GUO Xiaofeng, ZHENG Xuefeng, LU Manhong, et al. Design and application of high?data?rate QAM demodulator [J]. Journal of telemetry, tracking and command, 2011, 32(3): 21?25.

    [9] 唐婷,杜瑜.全數(shù)字化高速數(shù)傳解調(diào)器的設(shè)計與實現(xiàn)[J].科學(xué)技術(shù)與工程,2016,16(14):32?35.

    TANG Ting, DU Yu. The design and realization of fully?digitalized high?rate data demodulator [J]. Science technology and engineering, 2016, 16(14): 32?35.

    [10] 曾輝.基于全數(shù)字高速并行接收結(jié)構(gòu)APRX的解調(diào)技術(shù)研究[D].成都:電子科技大學(xué),2013:3?5.

    ZENG Hui. Research on demodulation techniques based on a parallel all?digital receiving structure named as APRX [D]. Chengdu: University of Electronic Science and Technology of China, 2013: 3?5.

    [11] 保駿.QPSK、SQPSK信號解調(diào)相位模糊及其對Viterbi譯碼的影響[J].兵器裝備工程學(xué)報,2011,32(3):53?55.

    BAO Jun. Demodulation phase ambiguity of QPSK and SQPSK and its influence on Viterbi decoding [J]. Journal of ordnance equipment engineering, 2011, 32(3): 53?55.

    大码成人一级视频| 99久久综合免费| 国产激情久久老熟女| 久久精品亚洲av国产电影网| 一级,二级,三级黄色视频| 国产在线免费精品| 蜜桃在线观看..| 久久久国产欧美日韩av| 亚洲精品成人av观看孕妇| 又黄又粗又硬又大视频| 亚洲色图综合在线观看| 男女高潮啪啪啪动态图| 欧美中文综合在线视频| 中国三级夫妇交换| 国产熟女午夜一区二区三区| 天天躁夜夜躁狠狠久久av| 亚洲成色77777| 亚洲国产日韩一区二区| 亚洲国产看品久久| 国产精品久久久久久精品古装| 97精品久久久久久久久久精品| av女优亚洲男人天堂| 曰老女人黄片| 国产片特级美女逼逼视频| 春色校园在线视频观看| 卡戴珊不雅视频在线播放| 我的亚洲天堂| 欧美国产精品一级二级三级| 国产精品欧美亚洲77777| 大片电影免费在线观看免费| 欧美精品国产亚洲| av免费在线看不卡| 色吧在线观看| 一级毛片 在线播放| 久久久久久伊人网av| 免费播放大片免费观看视频在线观看| 制服诱惑二区| 男女边摸边吃奶| 最近最新中文字幕大全免费视频 | 国产熟女午夜一区二区三区| 这个男人来自地球电影免费观看 | 国产精品 国内视频| 亚洲欧美精品自产自拍| 狠狠婷婷综合久久久久久88av| 国产精品.久久久| 久久久久精品人妻al黑| 人人妻人人爽人人添夜夜欢视频| 高清视频免费观看一区二区| 精品少妇黑人巨大在线播放| 亚洲国产日韩一区二区| 大香蕉久久网| 国产成人精品福利久久| 欧美激情极品国产一区二区三区| 亚洲,欧美,日韩| 蜜桃在线观看..| 国产精品二区激情视频| 国产伦理片在线播放av一区| 亚洲精品一二三| 日韩在线高清观看一区二区三区| 欧美97在线视频| 日本91视频免费播放| 精品亚洲成国产av| 91成人精品电影| 精品人妻一区二区三区麻豆| 亚洲国产毛片av蜜桃av| 精品亚洲成国产av| 国产精品麻豆人妻色哟哟久久| 久久精品国产a三级三级三级| 亚洲图色成人| 久久精品国产a三级三级三级| 欧美av亚洲av综合av国产av | 婷婷成人精品国产| 免费观看无遮挡的男女| 中文乱码字字幕精品一区二区三区| 午夜日韩欧美国产| 色播在线永久视频| 国语对白做爰xxxⅹ性视频网站| 亚洲国产最新在线播放| 99九九在线精品视频| 精品卡一卡二卡四卡免费| 大片免费播放器 马上看| 免费在线观看视频国产中文字幕亚洲 | 久久久a久久爽久久v久久| 成人午夜精彩视频在线观看| 18+在线观看网站| 不卡视频在线观看欧美| 捣出白浆h1v1| 亚洲中文av在线| 国产高清不卡午夜福利| 在线天堂最新版资源| 叶爱在线成人免费视频播放| 90打野战视频偷拍视频| 久久女婷五月综合色啪小说| 免费观看在线日韩| 天天躁夜夜躁狠狠久久av| 亚洲熟女精品中文字幕| 69精品国产乱码久久久| 制服诱惑二区| 国产麻豆69| 国产精品亚洲av一区麻豆 | 久久99蜜桃精品久久| 青草久久国产| 少妇被粗大猛烈的视频| 最新中文字幕久久久久| av片东京热男人的天堂| 午夜久久久在线观看| 亚洲情色 制服丝袜| 91在线精品国自产拍蜜月| 国产成人a∨麻豆精品| 超碰成人久久| 欧美精品亚洲一区二区| 亚洲精品中文字幕在线视频| 丝瓜视频免费看黄片| 巨乳人妻的诱惑在线观看| 亚洲欧美成人精品一区二区| 日本黄色日本黄色录像| 精品酒店卫生间| 男女下面插进去视频免费观看| 边亲边吃奶的免费视频| 成人手机av| 国产精品偷伦视频观看了| 一区二区av电影网| 高清av免费在线| 丰满迷人的少妇在线观看| 久久精品久久精品一区二区三区| videossex国产| 亚洲国产精品成人久久小说| 国产乱来视频区| 18禁观看日本| 日本91视频免费播放| 精品人妻在线不人妻| 国产一区二区三区综合在线观看| 亚洲精品美女久久久久99蜜臀 | 精品少妇内射三级| 一级毛片我不卡| 亚洲欧美精品综合一区二区三区 | 欧美日本中文国产一区发布| 熟女电影av网| a级毛片在线看网站| 免费在线观看黄色视频的| av免费在线看不卡| 久久综合国产亚洲精品| 亚洲婷婷狠狠爱综合网| 国产精品99久久99久久久不卡 | www.自偷自拍.com| 男女边吃奶边做爰视频| 卡戴珊不雅视频在线播放| 欧美日韩精品网址| 蜜桃国产av成人99| 在线观看一区二区三区激情| 日本欧美视频一区| 看免费成人av毛片| 边亲边吃奶的免费视频| 欧美日韩精品成人综合77777| 极品人妻少妇av视频| 亚洲精品国产一区二区精华液| 天美传媒精品一区二区| 午夜福利影视在线免费观看| 国产日韩欧美视频二区| 男人舔女人的私密视频| 一个人免费看片子| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 91aial.com中文字幕在线观看| 亚洲成国产人片在线观看| 久久精品久久精品一区二区三区| 亚洲四区av| 久久久久久人人人人人| 久久国产精品男人的天堂亚洲| 深夜精品福利| 久久人妻熟女aⅴ| 日本色播在线视频| 啦啦啦在线观看免费高清www| 丰满迷人的少妇在线观看| a级毛片黄视频| 午夜精品国产一区二区电影| 欧美变态另类bdsm刘玥| 伦理电影大哥的女人| 男的添女的下面高潮视频| 精品一区二区免费观看| 久久ye,这里只有精品| 亚洲国产最新在线播放| 久久久精品区二区三区| 国产xxxxx性猛交| 成年人免费黄色播放视频| 男的添女的下面高潮视频| 日韩三级伦理在线观看| 亚洲,一卡二卡三卡| www日本在线高清视频| 国产精品蜜桃在线观看| 久久精品熟女亚洲av麻豆精品| 性高湖久久久久久久久免费观看| 成年美女黄网站色视频大全免费| 美女主播在线视频| 一边亲一边摸免费视频| 天美传媒精品一区二区| 天美传媒精品一区二区| 日日爽夜夜爽网站| 久久 成人 亚洲| 最近手机中文字幕大全| 妹子高潮喷水视频| 欧美av亚洲av综合av国产av | 久久婷婷青草| 精品人妻熟女毛片av久久网站| 精品国产乱码久久久久久男人| 国产精品偷伦视频观看了| av免费在线看不卡| 综合色丁香网| 人妻人人澡人人爽人人| 韩国av在线不卡| 久久久久久人妻| 91精品三级在线观看| 中文字幕人妻丝袜制服| 国产精品成人在线| 国产男人的电影天堂91| 大片免费播放器 马上看| 人妻系列 视频| 国产又爽黄色视频| 在线精品无人区一区二区三| 制服丝袜香蕉在线| 国产精品秋霞免费鲁丝片| av网站在线播放免费| 91在线精品国自产拍蜜月| 国产精品免费视频内射| www.自偷自拍.com| 日韩在线高清观看一区二区三区| 成人亚洲欧美一区二区av| 亚洲精品久久久久久婷婷小说| 午夜激情久久久久久久| 亚洲 欧美一区二区三区| 欧美日韩av久久| 久久人人97超碰香蕉20202| 亚洲人成网站在线观看播放| 亚洲中文av在线| 大香蕉久久成人网| 777米奇影视久久| 日韩欧美精品免费久久| 亚洲国产成人一精品久久久| 18+在线观看网站| 永久网站在线| 日韩精品有码人妻一区| 国产激情久久老熟女| 热99国产精品久久久久久7| 女人被躁到高潮嗷嗷叫费观| 久热这里只有精品99| 日本免费在线观看一区| 青草久久国产| 自拍欧美九色日韩亚洲蝌蚪91| 少妇猛男粗大的猛烈进出视频| 日韩中文字幕欧美一区二区 | 日本av免费视频播放| 欧美另类一区| 国产老妇伦熟女老妇高清| 亚洲伊人色综图| 91精品三级在线观看| 亚洲av国产av综合av卡| 午夜久久久在线观看| 三级国产精品片| 欧美成人午夜免费资源| kizo精华| 边亲边吃奶的免费视频| 久久久久久伊人网av| 看免费av毛片| 亚洲五月色婷婷综合| 美女中出高潮动态图| 欧美xxⅹ黑人| 成人国产av品久久久| 亚洲精品一区蜜桃| 国产在线视频一区二区| 国产 一区精品| 国产人伦9x9x在线观看 | 久久久久久久久久久免费av| 大香蕉久久成人网| 久久久久视频综合| 我要看黄色一级片免费的| 亚洲av在线观看美女高潮| 免费大片黄手机在线观看| 男女边摸边吃奶| 黑人猛操日本美女一级片| av不卡在线播放| 国产亚洲av片在线观看秒播厂| 欧美日韩综合久久久久久| 亚洲美女视频黄频| 午夜日韩欧美国产| 日韩制服丝袜自拍偷拍| 岛国毛片在线播放| 亚洲人成电影观看| 国产亚洲最大av| 久久久久精品性色| 国产又爽黄色视频| 午夜福利在线观看免费完整高清在| 美女国产视频在线观看| 丰满少妇做爰视频| 日日啪夜夜爽| 美女脱内裤让男人舔精品视频| 亚洲国产av影院在线观看| 色播在线永久视频| 精品亚洲成a人片在线观看| 国产成人精品福利久久| 黄片无遮挡物在线观看| 观看美女的网站| 国产野战对白在线观看| www.av在线官网国产| 只有这里有精品99| 亚洲国产精品一区三区| 欧美成人午夜精品| 国产成人一区二区在线| 成人黄色视频免费在线看| 中文字幕av电影在线播放| 欧美 日韩 精品 国产| 国产精品.久久久| 精品久久久久久电影网| 最近手机中文字幕大全| 久久精品熟女亚洲av麻豆精品| 天天操日日干夜夜撸| 国产精品欧美亚洲77777| 另类精品久久| 伦理电影大哥的女人| 成人国产麻豆网| 欧美日韩国产mv在线观看视频| 亚洲欧美色中文字幕在线| 天天操日日干夜夜撸| 黑人猛操日本美女一级片| 少妇 在线观看| 91精品三级在线观看| 亚洲av成人精品一二三区| 精品卡一卡二卡四卡免费| 国产欧美日韩一区二区三区在线| 女性生殖器流出的白浆| 日韩制服骚丝袜av| 久久精品国产综合久久久| 精品一区二区三卡| 婷婷色麻豆天堂久久| 国产精品女同一区二区软件| 成人二区视频| 国产一区二区激情短视频 | 少妇的丰满在线观看| 最近的中文字幕免费完整| 欧美黄色片欧美黄色片| 最新中文字幕久久久久| 天天躁日日躁夜夜躁夜夜| 老女人水多毛片| 欧美精品一区二区大全| 欧美另类一区| 欧美成人精品欧美一级黄| 亚洲av电影在线进入| 国产一区亚洲一区在线观看| 最近的中文字幕免费完整| 欧美人与善性xxx| 久久精品久久久久久噜噜老黄| 久久精品国产综合久久久| 日韩视频在线欧美| 青草久久国产| av国产久精品久网站免费入址| 亚洲三级黄色毛片| 香蕉国产在线看| 精品卡一卡二卡四卡免费| 久久精品熟女亚洲av麻豆精品| 国产 一区精品| 一本久久精品| 久久婷婷青草| 中文字幕人妻熟女乱码| 一区二区三区乱码不卡18| 精品国产一区二区三区四区第35| 亚洲激情五月婷婷啪啪| 国产成人av激情在线播放| 精品亚洲乱码少妇综合久久| 亚洲av日韩在线播放| 久久ye,这里只有精品| 性色av一级| 久久精品人人爽人人爽视色| 国产成人91sexporn| 国产精品嫩草影院av在线观看| 亚洲精品一二三| 免费在线观看黄色视频的| 婷婷色综合www| 国产午夜精品一二区理论片| 午夜av观看不卡| 大片免费播放器 马上看| 成年女人在线观看亚洲视频| 超色免费av| 国产免费一区二区三区四区乱码| 欧美精品高潮呻吟av久久| 丝袜人妻中文字幕| 日日摸夜夜添夜夜爱| 久久久国产欧美日韩av| 精品福利永久在线观看| 色视频在线一区二区三区| 国产野战对白在线观看| 国产成人午夜福利电影在线观看| 美女主播在线视频| 亚洲综合精品二区| 亚洲精华国产精华液的使用体验| 亚洲欧美成人综合另类久久久| 欧美人与性动交α欧美软件| 欧美老熟妇乱子伦牲交| 亚洲av.av天堂| 麻豆乱淫一区二区| 波野结衣二区三区在线| 国产亚洲欧美精品永久| 麻豆乱淫一区二区| 99精国产麻豆久久婷婷| 色哟哟·www| 999久久久国产精品视频| 寂寞人妻少妇视频99o| 久久午夜综合久久蜜桃| 97在线人人人人妻| 亚洲av中文av极速乱| 卡戴珊不雅视频在线播放| 精品一区二区三卡| 一区在线观看完整版| a级毛片黄视频| 欧美日本中文国产一区发布| 成人国产av品久久久| 国产成人精品在线电影| 各种免费的搞黄视频| 精品国产超薄肉色丝袜足j| 国产熟女欧美一区二区| 伦理电影免费视频| 日日啪夜夜爽| 亚洲激情五月婷婷啪啪| 国产1区2区3区精品| 国产免费又黄又爽又色| 成人手机av| 五月天丁香电影| 久久av网站| av又黄又爽大尺度在线免费看| 国产无遮挡羞羞视频在线观看| 日本免费在线观看一区| 国产高清国产精品国产三级| 黑人巨大精品欧美一区二区蜜桃| 欧美激情高清一区二区三区 | 欧美日韩视频高清一区二区三区二| tube8黄色片| 国产精品国产三级专区第一集| 一级毛片我不卡| 深夜精品福利| 又粗又硬又长又爽又黄的视频| 欧美成人精品欧美一级黄| 熟女电影av网| 国产精品久久久久久av不卡| 久久精品aⅴ一区二区三区四区 | 精品一区在线观看国产| 亚洲三区欧美一区| 最近最新中文字幕大全免费视频 | 赤兔流量卡办理| 亚洲欧洲国产日韩| 欧美精品亚洲一区二区| 国产白丝娇喘喷水9色精品| 国产片内射在线| 老女人水多毛片| 亚洲av中文av极速乱| av卡一久久| 丝袜在线中文字幕| 久久久久国产网址| 777米奇影视久久| 久久久久精品性色| 国产精品不卡视频一区二区| 各种免费的搞黄视频| 18禁动态无遮挡网站| a级毛片在线看网站| av国产久精品久网站免费入址| 宅男免费午夜| 日本猛色少妇xxxxx猛交久久| 国产乱人偷精品视频| 国产一区二区三区av在线| 不卡av一区二区三区| 黄色一级大片看看| 多毛熟女@视频| 国产在线一区二区三区精| freevideosex欧美| 日本欧美国产在线视频| 久久久久久久久久久免费av| √禁漫天堂资源中文www| 少妇的逼水好多| 在线天堂中文资源库| 中文精品一卡2卡3卡4更新| 久热久热在线精品观看| 天天操日日干夜夜撸| 国产精品久久久久久av不卡| 亚洲,欧美精品.| 又粗又硬又长又爽又黄的视频| 丝瓜视频免费看黄片| 免费人妻精品一区二区三区视频| 国产在视频线精品| 国产在线视频一区二区| 在线观看三级黄色| 欧美成人午夜免费资源| 一区二区三区精品91| 在线 av 中文字幕| 人体艺术视频欧美日本| 黑人猛操日本美女一级片| 免费黄频网站在线观看国产| av免费观看日本| 亚洲精品成人av观看孕妇| 男的添女的下面高潮视频| 久久精品国产综合久久久| 久久久欧美国产精品| 可以免费在线观看a视频的电影网站 | 1024香蕉在线观看| 成人国产av品久久久| 91精品国产国语对白视频| 日本黄色日本黄色录像| 日韩成人av中文字幕在线观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 女性被躁到高潮视频| 中文天堂在线官网| 麻豆乱淫一区二区| 如日韩欧美国产精品一区二区三区| 久热久热在线精品观看| 亚洲人成电影观看| 国产片内射在线| 精品亚洲成a人片在线观看| 免费高清在线观看视频在线观看| 在线 av 中文字幕| 中文字幕另类日韩欧美亚洲嫩草| 久久 成人 亚洲| 黄片小视频在线播放| 亚洲av福利一区| 新久久久久国产一级毛片| 又大又黄又爽视频免费| 少妇 在线观看| 欧美精品一区二区大全| 啦啦啦中文免费视频观看日本| 日韩在线高清观看一区二区三区| 亚洲成国产人片在线观看| 水蜜桃什么品种好| 国产在线视频一区二区| 亚洲精品国产一区二区精华液| 国产日韩欧美亚洲二区| 免费观看a级毛片全部| 曰老女人黄片| 久久精品国产亚洲av高清一级| 一本久久精品| 熟女电影av网| 国精品久久久久久国模美| 欧美黄色片欧美黄色片| 制服人妻中文乱码| 亚洲国产看品久久| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲少妇的诱惑av| 亚洲精品视频女| 热re99久久国产66热| 秋霞伦理黄片| 色婷婷av一区二区三区视频| 一级a爱视频在线免费观看| 久久久久久久久久人人人人人人| 国产人伦9x9x在线观看 | 国产极品粉嫩免费观看在线| 久久精品夜色国产| 国产成人精品久久久久久| 午夜老司机福利剧场| av有码第一页| 夫妻午夜视频| 天天躁狠狠躁夜夜躁狠狠躁| 看非洲黑人一级黄片| 性色av一级| 黄片无遮挡物在线观看| 三级国产精品片| 日产精品乱码卡一卡2卡三| 精品少妇一区二区三区视频日本电影 | 色哟哟·www| 美女午夜性视频免费| 国产麻豆69| 欧美 日韩 精品 国产| 午夜av观看不卡| 18禁观看日本| 日韩一卡2卡3卡4卡2021年| 国产亚洲精品第一综合不卡| 亚洲av男天堂| 亚洲第一av免费看| 久久久精品94久久精品| 91成人精品电影| 一本色道久久久久久精品综合| 免费黄频网站在线观看国产| 成年女人毛片免费观看观看9 | 欧美成人午夜精品| 国产男女超爽视频在线观看| 久久精品久久精品一区二区三区| 深夜精品福利| 青青草视频在线视频观看| 尾随美女入室| 免费观看性生交大片5| 精品人妻熟女毛片av久久网站| 免费观看性生交大片5| 美女高潮到喷水免费观看| videossex国产| 亚洲欧美精品自产自拍| 色婷婷久久久亚洲欧美| 亚洲国产精品一区二区三区在线| a级毛片黄视频| 日韩av在线免费看完整版不卡| 少妇人妻精品综合一区二区| 汤姆久久久久久久影院中文字幕| 韩国av在线不卡| 国产黄频视频在线观看| 国产精品一二三区在线看| 久久99精品国语久久久| 亚洲精品中文字幕在线视频| 免费在线观看黄色视频的| 老熟女久久久| 日韩大片免费观看网站| 天天躁狠狠躁夜夜躁狠狠躁| 一本—道久久a久久精品蜜桃钙片| 精品亚洲成a人片在线观看| 亚洲av福利一区| 少妇的逼水好多| 两个人看的免费小视频| 免费在线观看完整版高清| 国精品久久久久久国模美| 90打野战视频偷拍视频| 超碰成人久久| 亚洲一区二区三区欧美精品| 国产精品国产三级国产专区5o| 欧美亚洲 丝袜 人妻 在线| 另类精品久久|