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

    一種基于DSPSR和GPU的脈沖星相干消色散基帶數(shù)據(jù)處理技術(shù)

    2021-10-30 02:17:20李浩羅近濤潘之辰孫鵬飛姜金辰
    時間頻率學(xué)報 2021年3期
    關(guān)鍵詞:共享內(nèi)存脈沖星色散

    李浩,羅近濤,潘之辰,孫鵬飛,姜金辰,5

    一種基于DSPSR和GPU的脈沖星相干消色散基帶數(shù)據(jù)處理技術(shù)

    李浩1,2,羅近濤1,3,潘之辰4,孫鵬飛1,姜金辰4,5

    (1. 中國科學(xué)院 國家授時中心,西安 710600;2. 中國科學(xué)院大學(xué),北京 100049; 3. 中國科學(xué)院 精密導(dǎo)航定位與定時技術(shù)重點實驗室,西安 710600;4. 中國科學(xué)院 國家天文臺,北京 100101;5. 北京大學(xué),北京 100871)

    相比于非相干消色散,相干消色散具有色散消除更徹底、可保留原始數(shù)據(jù)相位、時間分辨率信息等優(yōu)勢,但同時具有運算量大、數(shù)據(jù)處理耗時長等問題。為實現(xiàn)脈沖星基帶數(shù)據(jù)相干消色散處理,對數(shù)據(jù)處理流程進行加速,開發(fā)一種以DSPSR(digital signal processing software for pulsar)和GPU(graphics processing unit)為數(shù)據(jù)處理核心的脈沖星數(shù)據(jù)處理技術(shù)。利用高速共享內(nèi)存實現(xiàn)對環(huán)形緩存區(qū)兩路極化基帶數(shù)據(jù)的讀取,結(jié)合圖形處理器高性能計算技術(shù),使用DSPSR對基帶數(shù)據(jù)進行相干消色散處理。數(shù)據(jù)處理結(jié)果表明,相比于多線程非相干消色散系統(tǒng),此數(shù)據(jù)處理技術(shù)在數(shù)據(jù)處理速度方面具有更好的優(yōu)越性,并可獲得更精細的脈沖輪廓。

    脈沖星;相干消色散;脈沖星數(shù)字信號處理程序包;圖形處理器;脈沖星分布式采集和數(shù)據(jù)分析處理軟件

    0 引言

    脈沖星信號在傳播過程中由于受到星際介質(zhì)(interstellar medium,ISM)影響產(chǎn)生色散現(xiàn)象,導(dǎo)致觀測到的脈沖輪廓展寬,需采用消色散技術(shù)對脈沖星數(shù)據(jù)進行處理。脈沖星消色散方法主要分為非相干消色散和相干消色散兩種。非相干消色散在時域進行,多通道濾波器將脈沖星信號等分為若干通道的窄帶信號,再將每個通道的信號經(jīng)過平方檢波后在時域上對齊相加,補償因色散而產(chǎn)生的時間延遲,但是同一通道內(nèi)色散無法消除[1]。相干消色散核心算法為快速傅里葉變換,將時域數(shù)字信號變換到頻域,再使用星際介質(zhì)傳輸函數(shù)進行處理,消除整個通帶內(nèi)的色散。相比于非相干消色散,相干消色散可消除帶內(nèi)色散,效果更徹底,并能夠保留原始數(shù)據(jù)的相位關(guān)系、時間分辨率。已有學(xué)者研究表明,在一定的觀測頻率之內(nèi),相干消色散的效果優(yōu)于非相干消色散[2]。

    目前主流脈沖星數(shù)字終端系統(tǒng)由現(xiàn)場可編程門陣列(field programmable gate array,F(xiàn)PGA)、圖形處理器(graphics processing unit,GPU)、中央處理器(central processing unit,CPU)三部分組成,F(xiàn)PGA平臺主要用于ADC(analog to digital converter)控制、數(shù)據(jù)流控制和信號預(yù)處理,GPU負(fù)責(zé)消色散等核心算法的完成,CPU用于數(shù)據(jù)傳輸、邏輯控制和文件讀寫[3]。FPGA端輸出的基帶數(shù)據(jù)可保留脈沖星信號原始信息,結(jié)合相干消色散技術(shù)可開展超高時間分辨率脈沖星輻射特性研究,對于持續(xù)時間極短的快速射電暴(fast radio burst,F(xiàn)RB)搜尋以及脈沖星超短時標(biāo)脈沖等研究有著重要意義。以中國FAST 500 m望遠鏡①http://www.cas.cn/zt/kjzt/fast/、澳大利亞Parkes 64 m望遠鏡[4]、德國100 m口徑Effelsberg望遠鏡[5]等為代表的大口徑射電望遠鏡以脈沖星觀測為主要科學(xué)目標(biāo),寬帶寬、超高時間分辨率和高頻率分辨率數(shù)字終端技術(shù)成為發(fā)展趨勢。

    國內(nèi)脈沖星數(shù)字化終端系統(tǒng)的研制已取得一系列進展。中國科學(xué)院新疆天文臺依托25 m射電望遠鏡和VLBI記錄終端MK5A系統(tǒng),開展脈沖星相干消色散觀測系統(tǒng)研究[6]。中國科學(xué)院云南天文臺以ROACH2平臺和脈沖星數(shù)字信號處理程序包(digital pulsar signal processing for pulsars,DSPSR)為數(shù)據(jù)處理核心構(gòu)建相干消色散觀測系統(tǒng)[7]。中國科學(xué)院上海天文臺65 m射電望遠鏡配備脈沖星數(shù)字化終端(digital backend system,DIBAS),同時具備相干與非相干消色散功能[3]。中國科學(xué)院國家天文臺面向FAST開發(fā)了基于ROACH2平臺的非相干數(shù)字化脈沖星終端[8]。中國科學(xué)院國家授時中心40 m射電望遠鏡HRT脈沖星觀測系統(tǒng)使用FPGA+GPU架構(gòu),可實現(xiàn)非相干消色散與相干消色散兩種觀測模式[9]。國際上,美國NRAO(National Radio Astronomy Observatory)在GUPPI(green bank ultimate pulsar processing instrument)系統(tǒng)的基礎(chǔ)上,研制硬件平臺和性能指標(biāo)更為先進的VEGAS(versatile GBT astronomical spectrometer),基本系統(tǒng)由ADC,ROACH系列FPGA平臺與GPU處理集群組成②http://www.cas.cn/zt/kjzt/fast/。德國Effelsberg射電望遠鏡PSRIX后端采用ROACH系列平臺實現(xiàn),GPU集群進行數(shù)據(jù)處理[10]。澳大利亞Parkes射電望遠鏡的BPSR(Berkeley Parkes Swinburne Recorder)脈沖星數(shù)字終端系統(tǒng)采用FPGA平臺和CPU+GPU架構(gòu),數(shù)據(jù)處理后端采用開源軟件PSRDADA(processing software for pulsar distributed acquisition and data analysis)實現(xiàn)數(shù)據(jù)緩沖區(qū)、DSPSR實現(xiàn)實時相干消色散和RFI處理[11]。

    結(jié)合相干消色散脈沖星觀測系統(tǒng)發(fā)展趨勢,本文基于CPU+GPU架構(gòu),開發(fā)一種以GPU和DSPSR為數(shù)據(jù)處理核心的相干消色散基帶數(shù)據(jù)處理技術(shù)。利用高速共享內(nèi)存原理實現(xiàn)最高兩路極化基帶數(shù)據(jù)輸入,后續(xù)解碼、消色散、折疊等數(shù)據(jù)處理流程采用脈沖星數(shù)字信號處理程序包DSPSR結(jié)合GPU加速技術(shù)進行處理。數(shù)據(jù)處理結(jié)果以PSRFITS格式存儲。面向FAST脈沖星觀測系統(tǒng)海量數(shù)據(jù)相干消色散處理需求,將此數(shù)據(jù)處理流程在FAST基帶數(shù)據(jù)服務(wù)器集群上進行部署測試,驗證流程的正確性與可行性,為后續(xù)寬帶實時相干消色散脈沖星觀測系統(tǒng)的開發(fā)提供技術(shù)儲備。

    1 相干消色散原理

    2 DSPSR+GPU數(shù)據(jù)處理技術(shù)

    脈沖星數(shù)據(jù)進行相干消色散處理需進行大量的浮點運算,相比于非相干消色散算法,對計算機的性能和運算能力要求更高,隨著圖形處理器高性能計算的發(fā)展,數(shù)字化相干消色散脈沖星終端技術(shù)研究成為新的趨勢。本文構(gòu)建的相干消色散數(shù)據(jù)處理流程如圖1所示,利用高速共享內(nèi)存,即環(huán)形緩存區(qū)(ringbuffer)實現(xiàn)基帶數(shù)據(jù)的讀寫,包括數(shù)據(jù)輸入、數(shù)據(jù)處理、數(shù)據(jù)輸出3個部分。

    圖1 數(shù)據(jù)處理流程

    首先從磁盤中讀取基帶數(shù)據(jù)、header文件并寫入環(huán)形緩存區(qū),環(huán)形緩存區(qū)由兩部分組成:Data buffer和Header buffer,Data buffer用于存放基帶數(shù)據(jù),Header buffer用于存放觀測信息以及基帶數(shù)據(jù)基本信息。然后DSPSR從環(huán)形緩存區(qū)讀取基帶數(shù)據(jù)和對應(yīng)的信息,進行相干消色散、折疊等后續(xù)處理,并將處理結(jié)果存入磁盤。

    2.1 共享內(nèi)存技術(shù)

    共享內(nèi)存是Unix系統(tǒng)下一種高效的多進程間通信(inter-process communication,IPC)方式,多用于同一程序不同進程間的數(shù)據(jù)傳輸。進程可直接對內(nèi)存進行訪問、讀寫,不需要進行數(shù)據(jù)拷貝,是連接到共享內(nèi)存的多個進程間共享和傳遞數(shù)據(jù)的有效方式。

    共享內(nèi)存中的數(shù)據(jù)可被連接到此共享內(nèi)存的所有進程讀取,此時進程間數(shù)據(jù)傳輸不再通過內(nèi)核,可有效減少系統(tǒng)調(diào)用時間,提高數(shù)據(jù)傳輸效率。主流的共享內(nèi)存實現(xiàn)方式有跨進程間通訊的基于SYS V的共享內(nèi)存、基于POSIX mmap文件映射共享內(nèi)存和基于DMABUF框架,可跨設(shè)備間通訊的共享內(nèi)存兩大類[14]。

    脈沖星終端更寬帶寬、更高時間分辨率的發(fā)展趨勢,面臨海量數(shù)據(jù)高速傳輸技術(shù)難題。FAST 19波束接收機在8 bits采樣下基帶數(shù)據(jù)數(shù)據(jù)量為38 Gbit/s,PSRFITS格式數(shù)據(jù)量為3 040 Mbit/s(4 k通道, 50 μs,8 bits)[3],海量數(shù)據(jù)快速處理成為基帶數(shù)據(jù)處理難題,共享內(nèi)存高效數(shù)據(jù)傳輸?shù)膬?yōu)勢,是解決海量基帶數(shù)據(jù)處理的途徑。

    2.2 PSRDADA簡介

    PSRDADA是用于脈沖星基帶記錄的開源軟件項目,提供靈活且易于管理的環(huán)形緩存區(qū)(ringbuffer)實現(xiàn)數(shù)據(jù)讀寫。可通過基于文本的套接字接口來控制讀、寫程序,以實現(xiàn)從ringbuffer中將數(shù)據(jù)寫入磁盤陣列、將數(shù)據(jù)發(fā)送到節(jié)點陣列以及直接在ringbuffer中處理數(shù)據(jù)③http://psrdada.sourceforge.net/。

    在PSRDADA的邏輯中,環(huán)形緩沖區(qū)通常有一個指向ringbuffer中可讀數(shù)據(jù)的讀指針和一個指向可寫緩存區(qū)的寫指針,通過控制讀、寫指針實現(xiàn)緩沖區(qū)數(shù)據(jù)的輸入和輸出。在PSRDADA的具體實現(xiàn)中,寫程序創(chuàng)建并鎖定環(huán)形緩存區(qū)后,未被程序讀取的數(shù)據(jù)區(qū)域定義為使用區(qū)域,通過將已讀取數(shù)據(jù)區(qū)域標(biāo)記為已清除,成為新的未使用區(qū)域且可被寫程序識別,來構(gòu)成邏輯上的循環(huán)。

    2.3 脈沖星數(shù)字信號處理程序包

    DSPSR是一款開源、面向?qū)ο蟮臄?shù)字信號處理軟件包,可以選擇性地利用多核處理器和通用圖形處理單元進行射電脈沖星數(shù)字信號處理,旨在提高算法和數(shù)據(jù)的可移植性[13],支持CPSR,F(xiàn)ITS,Mark5等多種數(shù)據(jù)格式。DSPSR還可實現(xiàn)多種模塊化算法:如相干消色散、生成濾波器組、全極化檢測、格式轉(zhuǎn)換等。數(shù)據(jù)結(jié)果以PSRFITS格式輸出,可用于PSRCHIVE等軟件進行后續(xù)數(shù)據(jù)處理。

    DSPSR通過共享內(nèi)存區(qū),可直接從RAM中讀取脈沖星數(shù)據(jù)。并可結(jié)合GPU加速特性進行并行運算,提升消色散等數(shù)據(jù)處理效率。

    2.4 GPU加速技術(shù)

    GPU即圖形處理器,配備有高速內(nèi)存以及大量算術(shù)邏輯單元(arithmetic and logic unit,ALU),是專為執(zhí)行大量幾何計算而設(shè)計的一種高度并行化多核處理器。相比于CPU采用串行架構(gòu),功能模塊較多,能適應(yīng)復(fù)雜運算環(huán)境的特點,GPU可執(zhí)行大規(guī)模并發(fā)計算,擁有更強大的浮點運算能力。

    CUDA(compute unified device architecture)是NVIDIA公司推出的軟硬件通用計算設(shè)備架構(gòu)[15],采用標(biāo)準(zhǔn)C語言的拓展作為編程語言,提供大量高性能計算指令。利用CUDA進行編程可采用新的計算架構(gòu)來使用GPU提供的硬件資源,使線程的創(chuàng)建和調(diào)度所占用的開銷大大減少,從而提升程序的運算效率。

    GPU的邏輯架構(gòu)包括thread,block和grid三個部分,如圖2所示,GPU執(zhí)行計算可調(diào)用的最小單位為thread,數(shù)個thread共同組成一個block,一個block中的thread共用同一塊共享內(nèi)存,并可快速實現(xiàn)數(shù)據(jù)同步。在物理層,GPU硬件的核心組件為流式多處理器(streaming multiprocessor,SM)。SM資源主要包括CUDA核、共享內(nèi)存、寄存器、乘加器等,可以并行執(zhí)行數(shù)百個線程,資源數(shù)量越多則GPU的并行計算能力越高。

    圖2 GPU線程結(jié)構(gòu)

    脈沖星數(shù)據(jù)處理流程中,消色散所需時間占比最大,是影響數(shù)據(jù)處理效率最主要的因素[16]。隨著GPU的可編程能力及并行處理能力不斷提高,應(yīng)用范圍迅速擴展,可利用GPU高性能并行計算特性來滿足海量基帶數(shù)據(jù)的實時處理需求。本文利用DSPSR結(jié)合GPU加速技術(shù)對相干消色散算法進行加速,提高數(shù)據(jù)處理效率。

    3 數(shù)據(jù)處理

    本流程對基帶數(shù)據(jù)的處理可分為基帶數(shù)據(jù)的讀取與寫入、調(diào)用DSPSR讀取數(shù)據(jù)、相干消色散處理、折疊、PSRFITS格式輸出等5個步驟。基帶數(shù)據(jù)的讀寫通過編寫程序?qū)崿F(xiàn),通過在編譯時添加相應(yīng)的數(shù)據(jù)格式來實現(xiàn)DSPSR對ringbuffer中數(shù)據(jù)的讀取,相干消色散、折疊等處理通過調(diào)用DSPSR實現(xiàn),數(shù)據(jù)處理結(jié)果由DSPSR以PSRFITS格式保存。為提高數(shù)據(jù)處理速度,使用GPU對DSPSR的數(shù)據(jù)處理進行加速。

    3.1 數(shù)據(jù)處理原理

    在本數(shù)據(jù)處理流程中,環(huán)形緩存區(qū)為基帶數(shù)據(jù)讀寫的接口,包含指針、唯一標(biāo)識符(KEY)、內(nèi)存塊大小、內(nèi)存塊數(shù)量4個參數(shù)。這4個參數(shù)中,指針指向未使用的環(huán)形緩沖區(qū)句柄,通過此指針來實現(xiàn)對共享內(nèi)存的訪問,ringbuffer的大小由內(nèi)存塊大小和內(nèi)存塊數(shù)量共同決定。程序通過指定命令及以上參數(shù)來對ringbuffer進行創(chuàng)建、連接、鎖定、讀取、重置等操作。如圖3所示,對ringbuffer的使用以及數(shù)據(jù)處理流程步驟如下:

    ① 創(chuàng)建ringbuffer,通過設(shè)置內(nèi)存塊大小以及內(nèi)存塊數(shù)量來控制ringbuffer大小,并為其指定唯一標(biāo)識符‘KEY’。

    ② 寫程序通過唯一標(biāo)識符‘KEY’連接并鎖定ringbuffer,此時ringbuffer僅允許當(dāng)前寫程序進行訪問。

    ③ 寫程序通過命令行參數(shù)所指定的路徑與文件名,讀取header文件與基帶數(shù)據(jù)文件,并對ringbuffer進行寫操作。

    ④ 寫程序通過判斷是否讀取到文件末尾以保證文件的連續(xù)讀取,直至所有文件讀取完畢。

    ⑤ DSPSR通過唯一標(biāo)識符‘KEY’連接并鎖定到ringbuffer,對基帶數(shù)據(jù)進行讀取并進行后續(xù)相干消色散、折疊等處理,并輸出結(jié)果文件。此時ringbuffer僅允許當(dāng)前讀程序(DSPSR)進行訪問。

    ⑥ 全部數(shù)據(jù)處理完成后,寫程序會返回寫入ringbuffer中基帶數(shù)據(jù)文件總數(shù)、總大小以及平均傳輸速度到屏幕,然后通過對ringbuffer進行刷新,以釋放內(nèi)存。

    圖3 數(shù)據(jù)讀取原理

    處理兩路基帶數(shù)據(jù)時,需保證兩路數(shù)據(jù)文件的數(shù)量、每個文件大小分別對應(yīng)相等,目前實現(xiàn)兩路極化各2 048個文件輸入,根據(jù)服務(wù)器內(nèi)存配置,可調(diào)整文件輸入數(shù)量上限。在數(shù)據(jù)處理服務(wù)器內(nèi)存容量允許的情況下,可通過指定不同的標(biāo)識符‘KEY’創(chuàng)建多個ringbuffer來同時進行不同數(shù)據(jù)的讀寫操作。

    2001-2012年海南省國際旅游外匯收入除2003年和2009年下降外,總體呈波動上升趨勢。期間,2003年受亞洲“非典”公共危機事件影響,2008-2009年受全球金融危機影響出現(xiàn)下滑。其國際旅游外匯收入從2001年的1.06億美元上升到2012年的3.48億美元。12年間的平均國際旅游外匯收入為2.21億美元,約占全國的0.65%,排名在二十一位名上下波動。2001-2012年海南省國際旅游外匯收入一直低于全國平均值,差距越來越大。

    3.2 數(shù)據(jù)處理方式

    根據(jù)pipeline處理數(shù)據(jù)要求,使用ringbuffer前,需創(chuàng)建脈沖星觀測信息文本文件,文件中包括觀測頻率、射電源名稱、觀測模式、極化、通道、帶寬、中心頻率、采樣比特數(shù)、臺站信息、數(shù)據(jù)大小等參數(shù)信息,并對參數(shù)信息根據(jù)實際觀測情況進行修改,再將基帶數(shù)據(jù)文件名按照指定格式保存到文本文件中。

    依據(jù)3.1節(jié)數(shù)據(jù)處理原理,首先創(chuàng)建指定大小的ringbuffer,寫程序根據(jù)‘KEY’值使用指定命令連接到ringbuffer,將header文件及數(shù)據(jù)分別寫入header buffer和data buffer中,并輸出header文件內(nèi)容到屏幕。然后使用DSPSR相關(guān)命令讀取ringbuffer中數(shù)據(jù),進行相干消色散、分通道、折疊以及PSRFITS格式輸出等處理。通過修改DSPSR命令參數(shù)調(diào)整圖形處理器以及中央處理器線程數(shù)的使用,并對通道數(shù)、輸出文件名及路徑、積分時間等數(shù)據(jù)處理信息進行調(diào)整。所使用的脈沖星歷文件由Psrcat軟件生成,采用Tempo2生成預(yù)報多項式進行折疊。最后刷新ringbuffer,清空共享內(nèi)存區(qū)中數(shù)據(jù)。

    為測試、驗證數(shù)據(jù)處理流程的可行性與準(zhǔn)確性,在FAST基帶數(shù)據(jù)服務(wù)器集群上安裝配置此流程,并使用FAST對脈沖星PSR J0534+2200、PSR J1939+2134觀測所得的基帶數(shù)據(jù)進行測試處理。

    4 結(jié)果分析

    本文所處理的基帶數(shù)據(jù)觀測信息為:中心頻率1 250 MHz、帶寬500 MHz、雙極化、8比特采樣。首先,分別采用3.2節(jié)的數(shù)據(jù)處理方式和DSPSR CPU多線程方式對數(shù)據(jù)進行處理,記錄處理1 s觀測時長的數(shù)據(jù)量所需時間,結(jié)果如表1所示。根據(jù)FAST基帶數(shù)據(jù)服務(wù)器集群配置信息,本文介紹的數(shù)據(jù)處理流程處理8比特采樣、觀測時長1 s的雙極化基帶數(shù)據(jù),所需時間大約為0.8 s,與CPU模式下調(diào)用32線程處理時間5 s相比,速度提升約6.25倍。FAST基帶數(shù)據(jù)服務(wù)器配置信息如表2所示。

    表1 數(shù)據(jù)處理信息和所需時間

    表2 FAST基帶數(shù)據(jù)GPU配置信息

    其次,使用3.2節(jié)的數(shù)據(jù)處理方式對PSR J0534+2200和PSR J1939+2134的總觀測時長基帶數(shù)據(jù)進行相干消色散處理,處理結(jié)果如圖4至圖7所示。同時對PSR J1939+2134基帶數(shù)據(jù)進行非相干消色散處理,結(jié)果如圖8所示。PSR J0534+2200和PSR J1939+2134經(jīng)過相干消色散處理后的數(shù)據(jù)信噪比分別為16.677和1 270.176,可以利用Tempo2和PSRCHIVE對數(shù)據(jù)進行計時等研究工作。

    圖4 PSR J0534+2200相干消色散脈沖輪廓圖(中心頻率1 250 MHz,帶寬500 MHz,雙極化)

    從圖4中可清晰看到PSR J0534+2200脈沖輪廓的兩個主成分:分別位于相位0.3和0.7附近的主脈沖和中間脈沖。從圖5的頻率—相位圖中可以看到經(jīng)過相干消色散處理不同的通道已經(jīng)對齊,消除了色散效應(yīng)的影響。圖6為PSR J1939+2134經(jīng)過相干消色散處理后得到的脈沖輪廓圖,可看到位于相位0.2附近主脈沖的精細結(jié)構(gòu),在頻域圖中也可看到該結(jié)構(gòu)。圖7顯示各頻率通道也已經(jīng)對齊。

    圖6 PSR J1939+2134相干消色散脈沖輪廓圖(中心頻率1 250 MHz,帶寬500 MHz,雙極化)

    圖7 PSR J1939+2134頻率—相位圖

    圖8 PSR J1939+2134非相干消色散脈沖輪廓圖(中心頻率1 250 MHz,帶寬500 MHz,雙極化)

    最后,為驗證本文所論述的數(shù)據(jù)處理流程的數(shù)據(jù)處理質(zhì)量,根據(jù)Australia Telescope National Facility④https://data.csiro.au/dap/public/atnf/pulsarSearch.zul的資料,與Parkes 對PSR J1939+2134在L波段觀測所得到的輪廓進行比較,結(jié)果如圖9所示。圖9所示為Parkes的PDFB終端系統(tǒng)對PSR J1939+2134在1378 MHz觀測數(shù)據(jù)的相干消色散處理結(jié)果,可清晰看到位于相位0.2附近的主脈沖、相位0.7附近的中間脈沖以及主脈沖的精細結(jié)構(gòu),與圖6中FAST觀測數(shù)據(jù)處理結(jié)果一致,證明本數(shù)據(jù)處理流程的數(shù)據(jù)處理結(jié)果是符合要求的。

    圖9 Parkes觀測得到的PSR J1939+2134相干消色散脈沖輪廓圖

    5 結(jié)論

    本文提出的基于DSPSR和GPU的脈沖星數(shù)據(jù)處理流程,可實現(xiàn)最高兩路極化基帶數(shù)據(jù)相干消色散處理。數(shù)據(jù)讀取及處理方面結(jié)合高速共享內(nèi)存技術(shù)和GPU加速技術(shù),實現(xiàn)基于CPU+GPU異構(gòu)系統(tǒng)的相干消色散基帶數(shù)據(jù)處理技術(shù),相比于非相干消色散,色散消除更徹底,可以獲得更精確的脈沖輪廓。數(shù)據(jù)處理效率方面:在TITAN XP圖形處理器配置下,對8 bit采樣、觀測時長1 s的雙極化PSR J0534+2200基帶數(shù)據(jù)進行處理,所需時間約為0.8 s,與CPU模式下,DSPSR采用32線程處理同樣數(shù)據(jù)量的基帶數(shù)據(jù)所需時間約5 s相比,本數(shù)據(jù)處理流程速度提升約6.25倍。在觀測頻率、帶寬相近的條件下,將PSR J1939+2134數(shù)據(jù)處理結(jié)果與Parkes在L波段觀測所得的輪廓圖進行對比,本流程處理所得結(jié)果與Parkes處理結(jié)果相同,驗證了此數(shù)據(jù)處理流程在脈沖星數(shù)據(jù)處理方面的可行性與準(zhǔn)確性。本數(shù)據(jù)處理技術(shù)的開發(fā),在實現(xiàn)相干消色散基帶數(shù)據(jù)處理的基礎(chǔ)上,提升了數(shù)據(jù)處理效率,為海量脈沖星基帶數(shù)據(jù)實時處理以及國內(nèi)開發(fā)相干消色散脈沖星終端提供數(shù)據(jù)處理方面的技術(shù)儲備和技術(shù)參考。

    在今后的工作中將繼續(xù)研究DSPSR在GPU并行計算場景下的應(yīng)用,進一步優(yōu)化其數(shù)據(jù)處理速度,為實時寬帶相干消色散的接收機做技術(shù)儲備。并拓展基于PSRDADA共享內(nèi)存技術(shù)的數(shù)據(jù)讀寫功能,如增加支撐更多數(shù)據(jù)格式,實現(xiàn)特定格式文件輸出等。

    致謝:感謝中國科學(xué)院國家天文臺FAST現(xiàn)場工作人員對數(shù)據(jù)處理流程部署的大力支持和熱情幫 助,在此深表謝意。

    [1] 劉建偉, 金乘進. 基于計算機集群的脈沖星相干消色散系統(tǒng)研究[J]. 科學(xué)技術(shù)與工程, 2009, 9(13): 3612-3615+3620.

    [2] 黃玉祥, 汪敏, 郝龍飛, 等. 脈沖星信號相干消色散與非相干消色散的比較研究[J]. 天文研究與技術(shù), 2019, 16(1): 16-24.

    [3] 張海龍, 張萌, 聶俊, 等. 脈沖星數(shù)字終端技術(shù)綜述[J].中國科學(xué)(物理學(xué) 力學(xué) 天文學(xué)), 2019, 49(9): 15-28.

    [4] BAILES M. Precision timing at the Parkes 64-m radio telescope[J]. Radio Pulsars, ASP Conference Proceedings, San Francisco: Astronomical Society of the Pacific, 2003: 57-60.

    [5] HACHENBERG O, GRAHL B, WIELEBINSKI R. The 100-meter radio telescope at Effelsberg[J]. Proceedings of the IEEE, 1973, 61(9): 1288-1295.

    [6] 劉立勇, 艾力·伊沙木丁, 張晉. 烏魯木齊天文站建立脈沖星相干消色散觀測系統(tǒng)[J]. 天文研究與技術(shù), 2007, 4(1): 72-78.

    [7] 徐永華, 李紀(jì)云, 張穎倩, 等. 相干消色散脈沖星觀測系統(tǒng)的研究[J]. 天文研究與技術(shù), 2015, 12(4): 480-486.

    [8] 南仁東, 李會賢.FAST的進展—科學(xué)、技術(shù)與設(shè)備[J]. 中國科學(xué)(物理學(xué) 力學(xué) 天文學(xué)), 2014, 44(10): 1063-1074.

    [9] LUO J T, GGAO Y P, YANG T G, et al.Pulsar timing observations with Haoping radio telescope[J]. Research in Astronomy and Astrophysics, 2020, 20(7): 187-192.

    [10] LAZARUS P, KARUPPUSAMY R, GRAIKOU E, et al.Prospects for high-precision pulsar timing with the new Effelsberg PSRIX backend[J]. Monthly Notices of the Royal Astronomical Society, 2016, 458(1): 868-880.

    [11] SARKISSIAN J M, CARRETTI E, STRATEN W V.The Parkes pulsar backends[J]. AIP Conference Proceedings, 2011, 1357(1): 351-352.

    [12] LORMER D R, KRAMER M.Handbook of pulsar astronomy[K]. 2012.

    [13] STRATEN W V, BAILES M.DSPSR: digital signal processing software for pulsar astronomy[J]. Publications of the Astronomical Society of Australia, 201l, 28(1): 1-14.

    [14] Microsystems S.Shared Memory[C] // System Interface Guide, [S.l.]: Iuniverse Inc, 2001: 106-107.

    [15] CHENG J, GROSSMAN M, MCKERCHER T.Professional CUDA C programming[M]. John Wiley & Sons, [S.l.]: Wrox, 2014.

    [16] 李佳功, 徐永華, 李志玄, 等. 基于Mark5B+GPU脈沖星觀測系統(tǒng)[J]. 天文研究與技術(shù), 2014, 11(4): 335-342.

    A coherent de-dispersion processing technology on pulsar baseband data based on DSPSR and GPU

    LI Hao1,2, LUO Jin-tao1,3, PAN Zhi-chen4, SUN Peng-fei1, JIANG Jin-chen4,5

    (1. National Time Service Center, Chinese Academy of Sciences, Xi’an 710600, China;2.University of Chinese Academy of Sciences, Beijing 101048, China; 3. Key Laboratory of Precise Navigation and Timing Technology, Chinese Academy of Sciences, Xi’an 710600, China;4. National Astronomical Observatory, Chinese Academy of Sciences, Beijing 100101, China; 5. Peking University, Beijing 100871, China)

    Compared with incoherent de-dispersion, coherent de-dispersion has the advantages of more thorough dispersion elimination, and retaining the original phase information and the time resolution of the original data. However, it comes with challenges such as the massive computations and the longer data processing time. In order to realize a coherent de-dispersion processing of baseband data and accelerate the data processing, a pulsar data processing pipeline using DSPSR (the digital signal processing for pulsars) and GPU (graphics processing unit) as the data processing core is developed. Firstly, through the high-speed shared memory, the circular buffer is used to read the two-polarization baseband data. Combined with the high-performance computing technologies based on the graphics processor, the DSPSR is used to perform coherent de-dispersion, and folding process on the baseband data. Results from this pipeline show that compared with the multi-threaded incoherent de-dispersion system, this data processing technology has advantages in data processing speed and obtaining finer pulse profiles.

    pulsar; coherent achromatic dispersion; digital signal processing for pulsars (DSPSR); graphics processing unit (GPU); processing software for pulsar distributed acquisition and data analysis (PSRDADA)

    10.13875/j.issn.1674-0637.2021-03-0204-11

    李浩, 羅近濤, 潘之辰, 等. 一種基于DSPSR和GPU的脈沖星相干消色散基帶數(shù)據(jù)處理技術(shù)[J]. 時間頻率學(xué)報, 2021, 44(3): 204-214.

    2021-04-21;

    2021-05-22

    國家自然科學(xué)基金資助項目(U1931128;11973046;91736207;42030105);國家重點研發(fā)計劃資助項目(2017YFA0402602)

    猜你喜歡
    共享內(nèi)存脈沖星色散
    “光的折射”“光的色散”知識鞏固
    “中國天眼”已發(fā)現(xiàn)740余顆新脈沖星
    軍事文摘(2023年12期)2023-06-12 07:51:00
    “光的折射”“光的色散”知識鞏固
    發(fā)現(xiàn)脈沖星的女天文學(xué)家——貝爾
    科學(xué)(2022年4期)2022-10-25 02:43:42
    “光的折射”“光的色散”知識鞏固
    『光的折射』『光的色散』隨堂練
    通過QT實現(xiàn)進程間的通信
    基于PCI總線的多處理器協(xié)同機制研究
    基于虛擬觀測值的X射線單脈沖星星光組合導(dǎo)航
    QNX下PEX8311多路實時數(shù)據(jù)采集的驅(qū)動設(shè)計
    電子世界(2014年21期)2014-04-29 06:41:36
    99视频精品全部免费 在线| 国产男靠女视频免费网站| 亚洲国产精品sss在线观看| 国产精品三级大全| 国产亚洲av嫩草精品影院| 国产精品久久久久久久久免 | 久久久久免费精品人妻一区二区| 亚洲一区二区三区色噜噜| 亚洲乱码一区二区免费版| 久久精品综合一区二区三区| 一本一本综合久久| 免费人成在线观看视频色| 99视频精品全部免费 在线| 久久这里只有精品中国| 少妇裸体淫交视频免费看高清| 国产麻豆成人av免费视频| 免费观看精品视频网站| 老熟妇仑乱视频hdxx| 午夜免费观看网址| 俺也久久电影网| 国产精品 欧美亚洲| 一本综合久久免费| 精品久久久久久,| 噜噜噜噜噜久久久久久91| 可以在线观看的亚洲视频| 亚洲人成网站高清观看| 国产99白浆流出| 色吧在线观看| a在线观看视频网站| 亚洲不卡免费看| 久久久久国产精品人妻aⅴ院| 久久国产乱子伦精品免费另类| 色综合亚洲欧美另类图片| 日韩av在线大香蕉| 亚洲久久久久久中文字幕| 国产精品乱码一区二三区的特点| 一区福利在线观看| 一本久久中文字幕| 免费观看精品视频网站| 久久久久久人人人人人| 黄色视频,在线免费观看| 美女 人体艺术 gogo| 天堂av国产一区二区熟女人妻| 两个人的视频大全免费| a级毛片a级免费在线| 亚洲人成电影免费在线| 啦啦啦韩国在线观看视频| 国产一区二区三区视频了| 精品乱码久久久久久99久播| 欧美色视频一区免费| 欧美一级a爱片免费观看看| 成人特级黄色片久久久久久久| 99久久精品热视频| 一本一本综合久久| 人妻丰满熟妇av一区二区三区| 亚洲中文字幕一区二区三区有码在线看| 18禁国产床啪视频网站| 好男人在线观看高清免费视频| 欧美一级a爱片免费观看看| 亚洲一区二区三区色噜噜| 亚洲av成人av| 操出白浆在线播放| 国产亚洲精品一区二区www| 精品人妻1区二区| 噜噜噜噜噜久久久久久91| 久久久久性生活片| 国产精品国产高清国产av| 亚洲精品一区av在线观看| 亚洲18禁久久av| 九色成人免费人妻av| 麻豆国产97在线/欧美| 夜夜爽天天搞| 日韩欧美精品免费久久 | 日韩欧美免费精品| 麻豆成人午夜福利视频| 少妇的丰满在线观看| 国产在线精品亚洲第一网站| av专区在线播放| 黑人欧美特级aaaaaa片| 亚洲av不卡在线观看| 日韩大尺度精品在线看网址| 欧美绝顶高潮抽搐喷水| 国产99白浆流出| 无遮挡黄片免费观看| 精品人妻1区二区| 国产成人系列免费观看| 欧美午夜高清在线| 亚洲性夜色夜夜综合| 男插女下体视频免费在线播放| 国产精品亚洲av一区麻豆| 亚洲av成人精品一区久久| 亚洲va日本ⅴa欧美va伊人久久| av在线蜜桃| 久久精品国产自在天天线| АⅤ资源中文在线天堂| 亚洲国产精品成人综合色| 51国产日韩欧美| 午夜久久久久精精品| 女同久久另类99精品国产91| 欧美性猛交黑人性爽| 非洲黑人性xxxx精品又粗又长| 在线播放无遮挡| 亚洲国产精品999在线| 成人午夜高清在线视频| 精品一区二区三区av网在线观看| 国产三级中文精品| av视频在线观看入口| 国内久久婷婷六月综合欲色啪| 欧美日韩综合久久久久久 | 亚洲专区国产一区二区| 少妇的逼好多水| 99在线视频只有这里精品首页| av片东京热男人的天堂| 国产不卡一卡二| 美女高潮喷水抽搐中文字幕| 欧美黄色片欧美黄色片| 国产色婷婷99| 国产午夜福利久久久久久| 国产精品久久久久久精品电影| 国产精品女同一区二区软件 | 久久欧美精品欧美久久欧美| 日韩欧美三级三区| 天天躁日日操中文字幕| 美女高潮的动态| a级一级毛片免费在线观看| 熟女电影av网| 亚洲片人在线观看| 香蕉丝袜av| 很黄的视频免费| 中文字幕人妻熟人妻熟丝袜美 | 噜噜噜噜噜久久久久久91| 国模一区二区三区四区视频| 国产69精品久久久久777片| 91麻豆av在线| 看片在线看免费视频| av天堂中文字幕网| 国产三级在线视频| www.999成人在线观看| 日日干狠狠操夜夜爽| 一进一出抽搐动态| 最近最新中文字幕大全电影3| 一个人看视频在线观看www免费 | 精品欧美国产一区二区三| 亚洲男人的天堂狠狠| 免费av观看视频| 欧美三级亚洲精品| 18禁在线播放成人免费| 久久人妻av系列| 亚洲专区中文字幕在线| 又紧又爽又黄一区二区| 国产一区二区亚洲精品在线观看| 欧美成人免费av一区二区三区| 中文亚洲av片在线观看爽| 国产一区二区在线av高清观看| 国产亚洲欧美98| 欧美不卡视频在线免费观看| 99热只有精品国产| 女人被狂操c到高潮| 又紧又爽又黄一区二区| 欧美性猛交╳xxx乱大交人| 99热这里只有精品一区| 久久精品人妻少妇| 日本 欧美在线| 一级毛片女人18水好多| 国产野战对白在线观看| 亚洲成av人片免费观看| 欧美日韩国产亚洲二区| 美女高潮喷水抽搐中文字幕| 色综合欧美亚洲国产小说| 国产亚洲精品久久久久久毛片| 国产精品三级大全| 法律面前人人平等表现在哪些方面| 亚洲精品在线美女| 日日干狠狠操夜夜爽| 99热这里只有是精品50| 午夜福利在线在线| 黄片大片在线免费观看| 一个人看视频在线观看www免费 | 一本综合久久免费| 99热只有精品国产| 欧美极品一区二区三区四区| 大型黄色视频在线免费观看| 啪啪无遮挡十八禁网站| 韩国av一区二区三区四区| 草草在线视频免费看| 国产真实乱freesex| 成人永久免费在线观看视频| 国产精品免费一区二区三区在线| e午夜精品久久久久久久| 国产精品亚洲一级av第二区| 中文字幕熟女人妻在线| 国产aⅴ精品一区二区三区波| 又粗又爽又猛毛片免费看| 少妇裸体淫交视频免费看高清| 尤物成人国产欧美一区二区三区| 日本与韩国留学比较| 一区二区三区免费毛片| 国产精品亚洲美女久久久| 国产精品,欧美在线| 狂野欧美白嫩少妇大欣赏| 亚洲av五月六月丁香网| 九色成人免费人妻av| 午夜福利成人在线免费观看| 亚洲专区中文字幕在线| 国产精华一区二区三区| 国产伦在线观看视频一区| 男女午夜视频在线观看| 亚洲狠狠婷婷综合久久图片| 亚洲黑人精品在线| 日韩欧美免费精品| 欧美性猛交╳xxx乱大交人| 黄色片一级片一级黄色片| 久久国产精品人妻蜜桃| 小说图片视频综合网站| 999久久久精品免费观看国产| 国产欧美日韩一区二区三| 男女之事视频高清在线观看| 色视频www国产| 少妇的逼水好多| 大型黄色视频在线免费观看| 亚洲熟妇熟女久久| 日韩亚洲欧美综合| 99热这里只有是精品50| 熟女电影av网| 免费av毛片视频| 亚洲精品影视一区二区三区av| 精品久久久久久成人av| 国产视频一区二区在线看| 极品教师在线免费播放| 国产激情偷乱视频一区二区| 黄片小视频在线播放| 欧美3d第一页| 一区二区三区高清视频在线| 日本一二三区视频观看| 白带黄色成豆腐渣| 18禁在线播放成人免费| 国产成人av激情在线播放| eeuss影院久久| 亚洲成人久久爱视频| 白带黄色成豆腐渣| 国产伦人伦偷精品视频| 亚洲av中文字字幕乱码综合| 国产成年人精品一区二区| 欧美成人免费av一区二区三区| 久久久久久九九精品二区国产| 欧美黑人巨大hd| 欧美乱码精品一区二区三区| 欧美一区二区国产精品久久精品| 制服人妻中文乱码| 搡老岳熟女国产| 亚洲人成电影免费在线| 亚洲人成网站在线播放欧美日韩| 亚洲aⅴ乱码一区二区在线播放| 精品久久久久久久人妻蜜臀av| 国产精品久久久人人做人人爽| 最新在线观看一区二区三区| 中文字幕av成人在线电影| 日韩欧美精品免费久久 | 国内精品美女久久久久久| 日韩av在线大香蕉| 日本熟妇午夜| 综合色av麻豆| 99精品在免费线老司机午夜| 亚洲,欧美精品.| 老司机福利观看| 久久欧美精品欧美久久欧美| 好男人在线观看高清免费视频| 国产精品99久久99久久久不卡| 中文字幕熟女人妻在线| 精品久久久久久久人妻蜜臀av| 久久中文看片网| 婷婷丁香在线五月| 日韩成人在线观看一区二区三区| 黄色片一级片一级黄色片| 国产亚洲精品久久久久久毛片| 亚洲七黄色美女视频| 国产精品野战在线观看| 久久精品影院6| av黄色大香蕉| avwww免费| 亚洲欧美激情综合另类| 此物有八面人人有两片| 亚洲成av人片免费观看| 成年女人看的毛片在线观看| 日韩欧美精品免费久久 | 国产主播在线观看一区二区| 精品午夜福利视频在线观看一区| 成人三级黄色视频| 国产成人欧美在线观看| 亚洲真实伦在线观看| 午夜精品久久久久久毛片777| 欧美激情久久久久久爽电影| 观看美女的网站| 丁香六月欧美| 哪里可以看免费的av片| 午夜精品一区二区三区免费看| 国内揄拍国产精品人妻在线| 欧美极品一区二区三区四区| 中文字幕人成人乱码亚洲影| 亚洲第一欧美日韩一区二区三区| 亚洲美女黄片视频| 首页视频小说图片口味搜索| 色在线成人网| 成年免费大片在线观看| 婷婷精品国产亚洲av在线| 精品欧美国产一区二区三| 国产老妇女一区| 色播亚洲综合网| 日韩国内少妇激情av| 日本免费a在线| 真人做人爱边吃奶动态| 久久久久久久午夜电影| 天堂av国产一区二区熟女人妻| 亚洲精品色激情综合| 久久伊人香网站| 搡老岳熟女国产| 日韩有码中文字幕| 两个人的视频大全免费| 一级作爱视频免费观看| 麻豆久久精品国产亚洲av| a在线观看视频网站| 国语自产精品视频在线第100页| 免费看美女性在线毛片视频| 国产成+人综合+亚洲专区| 久久久国产精品麻豆| 看免费av毛片| 九色成人免费人妻av| 精品一区二区三区av网在线观看| 精品国产美女av久久久久小说| 色av中文字幕| 久久精品国产自在天天线| 成年女人毛片免费观看观看9| 女人高潮潮喷娇喘18禁视频| 色视频www国产| 神马国产精品三级电影在线观看| 黄色片一级片一级黄色片| 老熟妇仑乱视频hdxx| 1024手机看黄色片| 日韩欧美 国产精品| 给我免费播放毛片高清在线观看| 成年女人看的毛片在线观看| 好男人在线观看高清免费视频| 亚洲电影在线观看av| 18禁裸乳无遮挡免费网站照片| 激情在线观看视频在线高清| 国产精品久久久久久人妻精品电影| 老司机深夜福利视频在线观看| 亚洲性夜色夜夜综合| 亚洲国产欧美网| 国产一区二区在线av高清观看| 国产探花在线观看一区二区| 床上黄色一级片| 国产探花在线观看一区二区| 国产一区二区亚洲精品在线观看| 日韩欧美在线乱码| 超碰av人人做人人爽久久 | 国产黄a三级三级三级人| 精品一区二区三区av网在线观看| 亚洲熟妇熟女久久| 看黄色毛片网站| 九九热线精品视视频播放| 精品久久久久久久人妻蜜臀av| 女生性感内裤真人,穿戴方法视频| 成人一区二区视频在线观看| 国产一区二区在线av高清观看| 美女被艹到高潮喷水动态| 床上黄色一级片| 国产激情欧美一区二区| 亚洲无线在线观看| 国产精品亚洲美女久久久| 在线观看美女被高潮喷水网站 | 夜夜看夜夜爽夜夜摸| 男插女下体视频免费在线播放| 日韩中文字幕欧美一区二区| 嫁个100分男人电影在线观看| 国产精品久久久久久久电影 | 伊人久久大香线蕉亚洲五| bbb黄色大片| 日韩精品青青久久久久久| 国产美女午夜福利| 欧美日本亚洲视频在线播放| 色综合欧美亚洲国产小说| 国产不卡一卡二| 欧美中文综合在线视频| 搡女人真爽免费视频火全软件 | 搡老岳熟女国产| 99久国产av精品| av国产免费在线观看| 熟女少妇亚洲综合色aaa.| 热99在线观看视频| 国产综合懂色| 99在线人妻在线中文字幕| 成人一区二区视频在线观看| 亚洲狠狠婷婷综合久久图片| 伊人久久精品亚洲午夜| 亚洲精华国产精华精| 女同久久另类99精品国产91| 日韩亚洲欧美综合| 久久久久亚洲av毛片大全| 免费在线观看亚洲国产| 欧美在线黄色| 国产精品永久免费网站| 色吧在线观看| 国产欧美日韩精品一区二区| 欧美性猛交╳xxx乱大交人| 搞女人的毛片| 狠狠狠狠99中文字幕| 欧美激情久久久久久爽电影| 日韩高清综合在线| 91久久精品电影网| av片东京热男人的天堂| 色噜噜av男人的天堂激情| 久久伊人香网站| 国产精品,欧美在线| 国内精品久久久久精免费| 亚洲国产欧美人成| 好看av亚洲va欧美ⅴa在| 国产高清有码在线观看视频| 国产精品98久久久久久宅男小说| 嫁个100分男人电影在线观看| 精品一区二区三区人妻视频| 黄色片一级片一级黄色片| 2021天堂中文幕一二区在线观| 狂野欧美激情性xxxx| 国内揄拍国产精品人妻在线| 久久人妻av系列| 国产黄a三级三级三级人| 18+在线观看网站| 免费大片18禁| 亚洲国产精品合色在线| 桃红色精品国产亚洲av| 成人国产一区最新在线观看| 国产精品亚洲av一区麻豆| 午夜a级毛片| 国产 一区 欧美 日韩| 一二三四社区在线视频社区8| 亚洲aⅴ乱码一区二区在线播放| 成年版毛片免费区| 成年人黄色毛片网站| 亚洲成人久久性| 国产中年淑女户外野战色| 久久国产乱子伦精品免费另类| 校园春色视频在线观看| 深夜精品福利| 9191精品国产免费久久| 亚洲aⅴ乱码一区二区在线播放| 国产老妇女一区| 一级作爱视频免费观看| 日韩成人在线观看一区二区三区| 亚洲avbb在线观看| 国产精品一区二区三区四区久久| 9191精品国产免费久久| 一个人观看的视频www高清免费观看| 在线观看66精品国产| 久久午夜亚洲精品久久| 亚洲真实伦在线观看| 人妻久久中文字幕网| 宅男免费午夜| 午夜日韩欧美国产| 99久久无色码亚洲精品果冻| 欧美日韩黄片免| 国内精品久久久久久久电影| 一a级毛片在线观看| 中文字幕人妻熟人妻熟丝袜美 | 欧美日韩瑟瑟在线播放| 国模一区二区三区四区视频| 精品电影一区二区在线| 老司机福利观看| 精品乱码久久久久久99久播| 久久久久久久亚洲中文字幕 | av专区在线播放| 久久人人精品亚洲av| 国产成人aa在线观看| 成人高潮视频无遮挡免费网站| 一夜夜www| 又黄又粗又硬又大视频| 99国产精品一区二区蜜桃av| 波野结衣二区三区在线 | 久久精品国产99精品国产亚洲性色| 国产精品影院久久| 最近在线观看免费完整版| www日本在线高清视频| 极品教师在线免费播放| 美女被艹到高潮喷水动态| 男女床上黄色一级片免费看| 久久草成人影院| 国产av在哪里看| 色视频www国产| 人人妻,人人澡人人爽秒播| 51国产日韩欧美| 亚洲成人久久性| 手机成人av网站| 天天添夜夜摸| 可以在线观看毛片的网站| 国产精品久久久人人做人人爽| 亚洲成a人片在线一区二区| 男人舔奶头视频| 一进一出抽搐动态| 嫩草影视91久久| 亚洲人成电影免费在线| 久久精品91蜜桃| 给我免费播放毛片高清在线观看| 欧美日韩综合久久久久久 | 日本黄色视频三级网站网址| 中亚洲国语对白在线视频| 国产精品国产高清国产av| 亚洲国产精品成人综合色| 久久久久免费精品人妻一区二区| 国产精品女同一区二区软件 | aaaaa片日本免费| 国内少妇人妻偷人精品xxx网站| 亚洲欧美日韩高清专用| av天堂中文字幕网| 亚洲五月婷婷丁香| 老汉色∧v一级毛片| 日韩高清综合在线| 国内精品美女久久久久久| 免费观看的影片在线观看| 国产精品亚洲一级av第二区| 天天躁日日操中文字幕| 变态另类成人亚洲欧美熟女| 两个人视频免费观看高清| 欧美日韩综合久久久久久 | 日韩免费av在线播放| bbb黄色大片| 神马国产精品三级电影在线观看| 桃红色精品国产亚洲av| 欧美激情久久久久久爽电影| 夜夜爽天天搞| 成人永久免费在线观看视频| 俄罗斯特黄特色一大片| 中文在线观看免费www的网站| 3wmmmm亚洲av在线观看| 热99re8久久精品国产| av天堂中文字幕网| 国产高清有码在线观看视频| 韩国av一区二区三区四区| 国产精品三级大全| 亚洲在线观看片| 少妇的逼好多水| 最近最新中文字幕大全电影3| 一a级毛片在线观看| www国产在线视频色| 一区福利在线观看| 精品午夜福利视频在线观看一区| 精品人妻1区二区| 美女高潮的动态| 亚洲精品亚洲一区二区| 99久久精品国产亚洲精品| 女人高潮潮喷娇喘18禁视频| 欧美日韩一级在线毛片| 亚洲精品一区av在线观看| 中出人妻视频一区二区| avwww免费| 两个人看的免费小视频| 久久精品影院6| 国产成人a区在线观看| 国产精品久久久人人做人人爽| 午夜福利成人在线免费观看| 在线视频色国产色| 亚洲第一欧美日韩一区二区三区| 少妇丰满av| 亚洲国产中文字幕在线视频| 97人妻精品一区二区三区麻豆| 一卡2卡三卡四卡精品乱码亚洲| 久久中文看片网| 看免费av毛片| 欧美日韩中文字幕国产精品一区二区三区| 美女被艹到高潮喷水动态| 中文字幕熟女人妻在线| 一本综合久久免费| 岛国在线免费视频观看| 国语自产精品视频在线第100页| 男插女下体视频免费在线播放| 1000部很黄的大片| 国内精品久久久久久久电影| 97人妻精品一区二区三区麻豆| 国内久久婷婷六月综合欲色啪| 欧美黑人欧美精品刺激| 亚洲av中文字字幕乱码综合| 国产精品女同一区二区软件 | 很黄的视频免费| 精华霜和精华液先用哪个| 精品久久久久久久久久久久久| 国产国拍精品亚洲av在线观看 | 国产成人系列免费观看| 日本三级黄在线观看| 88av欧美| 在线国产一区二区在线| 色精品久久人妻99蜜桃| 无遮挡黄片免费观看| 国产欧美日韩精品一区二区| 国模一区二区三区四区视频| 少妇人妻一区二区三区视频| 国产伦精品一区二区三区四那| 18美女黄网站色大片免费观看| av中文乱码字幕在线| 人妻丰满熟妇av一区二区三区| 天天躁日日操中文字幕| 村上凉子中文字幕在线| 男女那种视频在线观看| 欧美日韩中文字幕国产精品一区二区三区| 老司机深夜福利视频在线观看| 乱人视频在线观看| 国模一区二区三区四区视频| 国产一区二区三区视频了| 桃红色精品国产亚洲av| 国产爱豆传媒在线观看| 在线观看免费午夜福利视频| 男人舔奶头视频| 中出人妻视频一区二区| 18禁国产床啪视频网站| 免费一级毛片在线播放高清视频| 99riav亚洲国产免费| 久久精品国产自在天天线|