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

    高頻地波雷達(dá)實(shí)時(shí)信號并行處理方案及實(shí)現(xiàn)*

    2017-02-15 05:26:09申晉晉
    關(guān)鍵詞:信號處理線程預(yù)處理

    牛 炯,申晉晉,陳 倩

    (中國海洋大學(xué)工程學(xué)院,山東 青島 266100)

    高頻地波雷達(dá)實(shí)時(shí)信號并行處理方案及實(shí)現(xiàn)*

    牛 炯,申晉晉,陳 倩

    (中國海洋大學(xué)工程學(xué)院,山東 青島 266100)

    近些年,隨著高頻地波雷達(dá)技術(shù)的發(fā)展,大量的遙感數(shù)據(jù)給實(shí)時(shí)信號處理帶來了巨大的挑戰(zhàn)。高頻地波雷達(dá)信號處理機(jī)的主流設(shè)計(jì)方案通常采用DSP+FPGA技術(shù)或基于高性能CPU技術(shù)來實(shí)現(xiàn)信號的實(shí)時(shí)處理,但DSP+FPGA方案中系統(tǒng)的軟硬件耦合較為緊密,開發(fā)周期長且不易維護(hù),而高性能CPU技術(shù)對計(jì)算機(jī)性能要求較高,大大增加了硬件成本。在保證信號實(shí)時(shí)處理的前提下,為了縮短開發(fā)周期,提高系統(tǒng)的可維護(hù)性和降低開發(fā)成本,本文基于CPU+GPU架構(gòu)設(shè)計(jì)了的信號處理方案,實(shí)現(xiàn)了高頻地波雷達(dá)實(shí)時(shí)信號處理的基本流程,獲取了目標(biāo)的距離信息和速度信息。實(shí)驗(yàn)結(jié)果表明,本文提出的CPU+GPU架構(gòu)下的信號處理方案與同期中央處理器(CPU)運(yùn)算平臺相比,具有10倍以上的加速比,驗(yàn)證了方案的可行性。

    CPU+GPU; 高頻地波雷達(dá); 實(shí)時(shí)信號處理; 協(xié)同工作;

    高頻地波雷達(dá),又稱為高頻地波超視距雷達(dá),是海上目標(biāo)大范圍、連續(xù)探測的主要手段,它利用高頻電磁波(3~30 MHz)沿海面繞射傳播的特性,實(shí)現(xiàn)海上船只和低空飛機(jī)等移動目標(biāo)和海洋動力參數(shù)的大范圍連續(xù)監(jiān)測[1-3]。因此,高頻地波雷達(dá)受到了全球研究機(jī)構(gòu)的關(guān)注,其主要研究內(nèi)容涵蓋系統(tǒng)設(shè)計(jì)、信號處理、雜波抑制、目標(biāo)提取以及海態(tài)反演等方面。隨著數(shù)字技術(shù)的飛速發(fā)展,雷達(dá)正向數(shù)字化和軟件化的方向發(fā)展。一直以來,軟件雷達(dá)實(shí)現(xiàn)的瓶頸問題之一是實(shí)時(shí)性問題[4]。其中信號處理的實(shí)時(shí)性尤為突出,因而雷達(dá)系統(tǒng)研發(fā)過程中,相當(dāng)一部分工作集中在信號處理的軟件化工作上,通常對于高頻地波雷達(dá)的信號處理機(jī),研究的主流方案有DSP+FPGA板卡方案[5-7]和高性能計(jì)算機(jī)方案[8]。對于DSP+FPGA方案來說,該方案的軟件和硬件耦合較為緊密,雖然計(jì)算效率高、可擴(kuò)展性好,但開發(fā)周期通常較長,并且不易維護(hù);高性能計(jì)算機(jī)方案雖然沒有軟硬件耦合的缺點(diǎn),但其為了滿足雷達(dá)實(shí)時(shí)信號處理的需求,對硬件性能有較高的要求通常需要花費(fèi)較高的成本。

    GPU(Graphic Processing Unit)在傳統(tǒng)上的應(yīng)用局限在處理圖形渲染計(jì)算任務(wù),受軍事視景仿真與圖像處理需求的牽引,GPU的數(shù)據(jù)處理能力取得了飛速的發(fā)展。目前大數(shù)據(jù)量進(jìn)行FFT運(yùn)算時(shí)GPU已經(jīng)是當(dāng)前CPU處理能力的幾十倍,而高頻地波雷達(dá)數(shù)字信號處理的核心就是FFT,在高頻地波雷達(dá)數(shù)據(jù)的前端處理過程中,需要進(jìn)行大量的FFT運(yùn)算,經(jīng)實(shí)測數(shù)據(jù)處理可知其耗時(shí)最長。目前GPU單精度浮點(diǎn)處理能力超過11TFlops,同時(shí)其雙精度浮點(diǎn)處理能力也已經(jīng)超過1.4TFlops。因此將GPU應(yīng)用于雷達(dá)數(shù)字信號處理是一種趨勢。近些年已經(jīng)有大量的研究機(jī)構(gòu)將GPU技術(shù)用于數(shù)字信號處理。文獻(xiàn)[9]中將GPU用于合成孔徑雷達(dá)(SAR)的信號處理,利用GPU通用并行計(jì)算,使用CUDA(Compute Unified Device Architecture)實(shí)現(xiàn)了SAR成像算法。文獻(xiàn)[10]采用CPU-GPU技術(shù)設(shè)計(jì)了微波雷達(dá)終端處理系統(tǒng),數(shù)據(jù)處理加速比達(dá)到5.56 。文獻(xiàn)[11]提出一種基于GPU的外輻射源雷達(dá)信號處理實(shí)時(shí)方法,利用GPU的并行處理能力實(shí)現(xiàn)了外輻射源雷達(dá)信號處理實(shí)時(shí)處理。GPU在微波雷達(dá)領(lǐng)域的應(yīng)用已經(jīng)較為廣泛,但目前尚未見到將CPU+GPU技術(shù)應(yīng)用于高頻地波雷達(dá)數(shù)字信號處理的相關(guān)文獻(xiàn)。本文結(jié)合高頻地波雷達(dá)數(shù)字信號處理的流程及實(shí)時(shí)處理要求,提出采用CPU+GPU的并行異構(gòu)模式對高頻地波雷達(dá)實(shí)測數(shù)據(jù)進(jìn)行并行處理。

    1 理論方法及數(shù)據(jù)處理流程

    高頻地波雷達(dá)回波信號的數(shù)字化處理主要在信號處理機(jī)內(nèi)完成,前端數(shù)據(jù)處理流程見圖1。獲取指定長度的原始數(shù)據(jù)后,首先需要對原始數(shù)據(jù)進(jìn)行預(yù)處理,包括去直達(dá)波處理、去斜處理和加窗處理,然后對信號進(jìn)行解距離處理后獲得目標(biāo)的距離信息,最后對解距離后的數(shù)據(jù)進(jìn)行解速度處理得到目標(biāo)的速度信息。

    假設(shè)高頻地波雷達(dá)的發(fā)射波形采用線性調(diào)頻斷續(xù)波(FMICW),發(fā)射信號表示為:

    (1)

    其中:

    為了方便分析,不考慮信號在傳輸過程中幅度的衰減,則接收的回波信號可以表示為:

    (2)

    (3)

    式中:K為單個(gè)脈沖重復(fù)周期內(nèi)的采樣點(diǎn)數(shù);P為一個(gè)調(diào)頻周期的脈沖個(gè)數(shù);m為雷達(dá)最大探測距離單元。當(dāng)解距離后數(shù)據(jù)的時(shí)間長度滿足設(shè)定的積累周期時(shí),對Y(m)進(jìn)行第二次FFT處理,即解速度處理,可得到目標(biāo)的速度信息。高頻地波雷達(dá)前端信號處理詳細(xì)的公式推導(dǎo)見文獻(xiàn)[5]。

    圖1 前端數(shù)據(jù)處理流程圖Fig.1 Front-end data processing flowchart

    2 基于CPU+GPU架構(gòu)的協(xié)同并行計(jì)算方案

    本文提出了基于CPU+GPU架構(gòu)的信號處理方案,本方案采用CPU+GPU協(xié)同并行的計(jì)算來實(shí)現(xiàn)高頻地波雷達(dá)前端數(shù)據(jù)的處理。

    2.1 基于CPU+GPU架構(gòu)的系統(tǒng)方案

    高頻地波雷達(dá)的接收站包括接收天線陣列、接收機(jī)、信號處理機(jī)以及顯控終端等。接收天線陣列主要用于接收來自空間的電磁波,并將雷達(dá)回波信號通過饋線傳輸給接收機(jī);接收機(jī)主要用于接收來自接收天線陣列的回波信號,并對回波信號進(jìn)行預(yù)處理,同時(shí)監(jiān)測空間的電磁環(huán)境。接收機(jī)接收到回波數(shù)據(jù)后,需要對原始回波數(shù)據(jù)進(jìn)行放大、濾波、A/D采樣以及I/Q信號分離等處理,然后將處理后的數(shù)據(jù)傳輸?shù)叫盘柼幚頇C(jī),信號處理機(jī)接收到來自接收機(jī)的數(shù)字化正交信號(I/Q信號)后,在內(nèi)部對前端數(shù)據(jù)進(jìn)行預(yù)處理(包括去直達(dá)波、去斜和加窗)、解距離、解速度、點(diǎn)跡檢測以及雜波抑制等(干擾包括射頻干擾、電離層干擾、海雜波以及瞬態(tài)干擾)相關(guān)處理,并對原始數(shù)據(jù)、脈沖壓縮后數(shù)據(jù)、頻域數(shù)據(jù)以及點(diǎn)跡數(shù)據(jù)進(jìn)行實(shí)時(shí)存儲,最后將點(diǎn)跡數(shù)據(jù)實(shí)時(shí)的傳送到顯控終端,在顯控終端對目標(biāo)點(diǎn)跡進(jìn)行航跡關(guān)聯(lián),并實(shí)時(shí)顯示目標(biāo)運(yùn)動的事態(tài)?;贑PU+GPU架構(gòu)的信號處理機(jī)方案框圖見圖2。

    圖2 基于CPU+GPU架構(gòu)的信號處理機(jī)方案框圖

    在整個(gè)系統(tǒng)中,信號處理機(jī)采用CPU+GPU架構(gòu),CPU與GPU在信號處理機(jī)內(nèi)協(xié)同工作,CPU作為主機(jī)(Host),GPU作為協(xié)處理器。CPU、GPU各自擁有相互獨(dú)立的存儲器地址空間:主機(jī)端的內(nèi)存和設(shè)備端的顯存。

    2.2 基于CPU+GPU架構(gòu)的并行處理算法

    高頻地波雷達(dá)前端數(shù)據(jù)在CPU+GPU架構(gòu)下協(xié)同并行處理的方案見圖3,其中CPU主要負(fù)責(zé)邏輯性強(qiáng)的事件處理和串行計(jì)算,包括計(jì)算需要開辟內(nèi)存的大小、獲取雷達(dá)系統(tǒng)參數(shù)、截取指定時(shí)長的原始數(shù)據(jù)、清理上一個(gè)內(nèi)核函數(shù)、啟動下一個(gè)內(nèi)核函數(shù)、數(shù)據(jù)存儲以及RD譜顯示等;GPU則專注于執(zhí)行高度線程化的并行處理事務(wù),包括原始數(shù)據(jù)的預(yù)處理、解距離和解速度等處理。

    為了保證方案的通用性,系統(tǒng)應(yīng)該采用模塊化和參數(shù)化設(shè)計(jì)。CPU根據(jù)系統(tǒng)參數(shù)設(shè)置,獲取了指定長度的高頻雷達(dá)原始數(shù)據(jù)后,通過PCIE總線將該數(shù)據(jù)傳輸?shù)紾PU,假設(shè)該批數(shù)據(jù)的大小為N個(gè)字節(jié),則:

    N=nfre×nch×nT×TSweep×fs×2×4。

    (4)

    其中:nfre為雷達(dá)載頻數(shù);nch為接收天線陣元數(shù);nT為該批數(shù)據(jù)包含的掃頻周期數(shù);TSweep為掃頻周期;fs為采樣率。

    (4)式中乘以2是因?yàn)榛夭〝?shù)據(jù)為復(fù)數(shù),包含實(shí)部和虛部,乘以4是由于獲取的原始數(shù)據(jù)為float型。在GPU內(nèi)對該批數(shù)據(jù)進(jìn)行預(yù)處理和解距離后,獲得的數(shù)據(jù)大小為N1個(gè)字節(jié),則:

    N1=nfre×nch×m×nT×2×4。

    (5)

    其中m表示雷達(dá)探測的最大距離單元,根據(jù)公式4和公式5可以看出,原始數(shù)據(jù)進(jìn)行預(yù)處理和解距離處理后,數(shù)據(jù)壓縮比為TSweep*fs/m。根據(jù)表1中雷達(dá)實(shí)際的相關(guān)參數(shù),做完解距離處理后,數(shù)據(jù)壓縮比約為100∶1,得到的脈壓后數(shù)據(jù)存儲在GPU顯存中相應(yīng)的位置,當(dāng)一個(gè)積累周期對應(yīng)的時(shí)域數(shù)據(jù)全部處理完畢后,對其進(jìn)行解速度處理,得到頻域數(shù)據(jù)??梢?,大量原始數(shù)據(jù)的復(fù)雜處理都在GPU內(nèi)完成,利用GPU并行處理的特點(diǎn),極大的提高了數(shù)據(jù)處理速度,節(jié)省了CPU資源,保證了系統(tǒng)的實(shí)時(shí)性、穩(wěn)定性和可靠性。

    圖3 CPU+GPU協(xié)同并行處理方案

    為了更好的理解并行優(yōu)化策略,下面針對CPU+GPU架構(gòu)下數(shù)據(jù)處理流程的關(guān)鍵技術(shù)進(jìn)行介紹,包括CPU/GPU之間的內(nèi)存分配和數(shù)據(jù)傳輸、GPU中的數(shù)據(jù)預(yù)處理以及GPU中的FFT處理。

    2.2.1 CPU+GPU的內(nèi)存分配與數(shù)據(jù)傳輸

    CPU與GPU之間通過PCI Express接口進(jìn)行數(shù)據(jù)通信。系統(tǒng)所要處理的數(shù)據(jù)首先要傳到Host端的內(nèi)存中即CPU的內(nèi)存中,然后通過PCIE總線傳輸?shù)斤@存,在GPU內(nèi)進(jìn)行數(shù)據(jù)處理。

    首先根據(jù)雷達(dá)系統(tǒng)相關(guān)參數(shù)計(jì)算需要開辟內(nèi)存大小,在Host端開辟相應(yīng)的內(nèi)存,用來存儲原始數(shù)據(jù)、本振信號數(shù)據(jù)以及窗函數(shù)數(shù)據(jù)等。系統(tǒng)完成一個(gè)積累周期的數(shù)據(jù)處理需要調(diào)用多次數(shù)據(jù)預(yù)處理函數(shù)和解距離函數(shù),考慮到GPU顯存資源有限,采用在解距離之前一次性開辟出固定大小的內(nèi)存,供解距離時(shí)重復(fù)使用,在每次解距離之前需要開辟相應(yīng)大小的內(nèi)存,然后在完成最后一次解距離處理時(shí)統(tǒng)一釋放。在GPU 內(nèi)開辟內(nèi)存可以用CudaMalloc函數(shù)來實(shí)現(xiàn)。GPU每執(zhí)行一次解距離運(yùn)算,都需要將nch個(gè)接收通道的CIT_base_Num個(gè)掃頻周期的數(shù)據(jù)從CPU傳輸?shù)紾PU,我們使用CUDAMemcpy函數(shù)來實(shí)現(xiàn)CPU與GPU間的數(shù)據(jù)傳輸。其中CudaMalloc函數(shù)與CudaMemcpy函數(shù)均為CUDA的庫函數(shù)。2.2.2 GPU中的數(shù)據(jù)預(yù)處理

    在解距離之前,首先需要對回波數(shù)據(jù)進(jìn)行預(yù)處理。

    對于采用收發(fā)共址的模式高頻地波雷達(dá)系統(tǒng),發(fā)射線性調(diào)頻斷續(xù)波,在發(fā)射機(jī)工作的時(shí)候,接收機(jī)關(guān)閉接收通道,不接收回波信號,但由于發(fā)射機(jī)發(fā)射功率較大,即使接收機(jī)關(guān)閉了接收通道,發(fā)射信號也會泄露進(jìn)入接收機(jī),而這部分能量通常高于接收機(jī)正常工作時(shí)接收到的回波信號的能量一個(gè)數(shù)量級以上,因此必須對直達(dá)波進(jìn)行剔除。去直達(dá)波操作就是對每個(gè)掃頻周期內(nèi)的PulseNum個(gè)脈沖的直達(dá)波對應(yīng)位置的采樣點(diǎn)置零,即將每個(gè)脈沖前nTpp*fs個(gè)采樣點(diǎn)置零。實(shí)際操作中GPU啟動Nthread個(gè)線程,如公式(6)所示:

    Nthread=nTpp×fs×PulseNum×nch×nfre×BaseNum。

    (6)

    其中,包含Setzero_block個(gè)線程塊,每個(gè)塊包含Setzero_thread個(gè)線程,表達(dá)式如下:

    Setzero_block=nTpp×fs×PulseNum×nch,

    (7)

    Setzero_thread=nfre×Base_Num,

    (8)

    其中Base_Num表示用戶設(shè)置的每次處理數(shù)據(jù)包含的掃頻周期數(shù)。

    去斜處理是將去除直達(dá)波后的回波信號與本振信號進(jìn)行復(fù)乘。去斜處理時(shí),一個(gè)線程內(nèi)執(zhí)行一個(gè)復(fù)乘,每個(gè)線程執(zhí)行復(fù)乘時(shí)需要占用較多的寄存器,而GPU硬件上寄存器的數(shù)量是一定的,啟動內(nèi)核時(shí),系統(tǒng)會將寄存器平均分配到每個(gè)線程中,操作時(shí)應(yīng)開辟較少的線程以保證每個(gè)線程可以分配到足夠的寄存器。實(shí)際操作中GPU一次啟動開辟nTs個(gè)線程。

    加窗操作是將去斜處理后數(shù)據(jù)的實(shí)部和虛部都與相應(yīng)的窗加權(quán)系數(shù)相乘,并行處理方法方法同去斜處理相同,GPU一次啟動nTs個(gè)線程執(zhí)行計(jì)算。

    2.2.3 GPU中的FFT處理

    在高頻地波雷達(dá)信號前端數(shù)據(jù)的處理過程中,通過解距離和解速度分別得到目標(biāo)的距離信息和速度信息,而這2個(gè)模塊的關(guān)鍵技術(shù)就是FFT。NVIDIA為開發(fā)人員提供了多種庫函數(shù),每一類庫針對某一特定的應(yīng)用,CUFFT(CUDA Fast Fourier Transform)則是CUDA專門用于進(jìn)行快速傅里葉變換的函數(shù)庫,該函數(shù)庫是經(jīng)過充分優(yōu)化的。該函數(shù)庫為開發(fā)人員提供了完善的用戶接口,開發(fā)人員可以任意設(shè)置FFT的點(diǎn)數(shù)、次數(shù)以及數(shù)據(jù)類型。

    3 實(shí)測數(shù)據(jù)驗(yàn)證

    為驗(yàn)證提出的基于CPU+GPU架構(gòu)的高頻地波雷達(dá)信號并行處理技術(shù)的加速性能,分別采用CPU+GPU架構(gòu)和純CPU架構(gòu)兩種方案對青島高頻地波雷達(dá)系統(tǒng)于2015年1月8日錄取的實(shí)測數(shù)據(jù)采用脫機(jī)處理的方式進(jìn)行處理,該高頻地波雷達(dá)為收發(fā)共址,發(fā)射線性調(diào)頻斷續(xù)波,信號處理機(jī)內(nèi)CPU為Intel Corei5-3470,主頻3.2 GHz,GPU計(jì)算卡采用英偉達(dá)的NVIDIA GTX960,其內(nèi)部有1 024個(gè)CUDA core,內(nèi)置4 096 MB GDDR5顯存,單精度峰值處理能力和峰值顯存帶寬可達(dá)到2.6 T Flops和232 GB/s。按照圖1構(gòu)建完整的地波雷達(dá)信號處理流程,將獲取的地波雷達(dá)數(shù)據(jù)分別在不同的開發(fā)平臺下進(jìn)行處理,表1給出了兩種方案處理一個(gè)積累周期數(shù)據(jù)時(shí)各個(gè)模塊的耗時(shí)結(jié)果。采用CPU+GPU架構(gòu)方案時(shí),處理一個(gè)積累周期的數(shù)據(jù)耗時(shí)約8 s,完全能夠滿足高頻地波雷達(dá)實(shí)時(shí)信號處理系統(tǒng)的滑窗需求,同時(shí)極大的節(jié)省了CPU的資源,使CPU能夠可靠的進(jìn)行指令的上傳和下達(dá);而純CPU方案處理相同量級的數(shù)據(jù)則耗時(shí)較長,無法滿足高頻雷達(dá)實(shí)時(shí)性處理需求,若要提高其處理速度,則需要提升硬件性能,增加成本。

    表1 一個(gè)積累周期下兩種方案耗時(shí)對比

    Note:①M(fèi)emory development;②Data preprocessing and resolving Distance (16 Times);③Time consuming for the first call program;④Time consuming for each subsequent call program;⑤Total time-consuming;⑥Resolving speed;⑦Total time-consuming

    圖4給出了兩種方案連續(xù)處理10個(gè)積累周期耗時(shí)的曲線圖。從表1和圖4中可以看出對高頻地波雷達(dá)數(shù)據(jù)進(jìn)行預(yù)處理、解距離和解速度處理時(shí),CPU+GPU方案較純CPU方案的處理速度提高了10倍以上。

    圖4 兩種方案處理10個(gè)積累周期數(shù)據(jù)耗時(shí)對比

    圖5 GPU+CPU方案處理第一個(gè)陣元回波數(shù)據(jù)的RD譜Fig.5 The RD spectrum of the first element echo data processed by the GPU+CPU scheme

    圖5所示為采用CPU+GPU方案處理4.7 MHz載頻下第一個(gè)接收通道回波數(shù)據(jù)后得到的RD譜,積累時(shí)間260 s。對兩種方案獲得的頻域數(shù)據(jù)的實(shí)部和虛部分別做差求模,兩種方案頻域數(shù)據(jù)做差后各點(diǎn)數(shù)據(jù)的相對偏差小于10-7,該偏差不會影響雷達(dá)的探測性能,驗(yàn)證了方案的可行性。兩種方案處理同批數(shù)據(jù)產(chǎn)生的數(shù)據(jù)偏差是由于不同編譯軟件中數(shù)據(jù)的精度誤差導(dǎo)致。圖6表示兩種方案獲得的第一個(gè)陣元的頻域數(shù)據(jù)做差后的結(jié)果。

    圖6 兩種方案的頻域數(shù)據(jù)做差的結(jié)果Fig.6 The results of the subtraction of the frequency domain data of the two schemes

    4 結(jié)語

    本文設(shè)計(jì)了基于CPU+GPU架構(gòu)的高頻地波雷達(dá)信號處理系統(tǒng),并在此架構(gòu)上實(shí)現(xiàn)了基于GPU的信號并行處理。實(shí)測數(shù)據(jù)的處理結(jié)果表明,基于GPU加速的雷達(dá)信號處理系統(tǒng)能夠?qū)崟r(shí)完成信號處理,且該方案與純CPU方案相比,數(shù)據(jù)加速比達(dá)到了10,驗(yàn)證了該方案的可行性,能夠很好的滿足軟件化雷達(dá)信號實(shí)時(shí)處理要求,為軟件化雷達(dá)的數(shù)據(jù)快速處理奠定了基礎(chǔ)。另外,由于雜波抑制算法涉及到較為復(fù)雜的矩陣運(yùn)算,使用純CPU處理將無法滿足實(shí)時(shí)性的要求,下一步將在GPU內(nèi)實(shí)現(xiàn)雜波抑制算法。

    [1] Headrick J M,Skolnik M I.Over-the-Horizon radar in the HF band[J].Proceedings of the IEEE,1974,62(6):664-673.

    [2] Ponsford A M.Surveillance of the 200 nautical mile Exclusive Economic Zone (EEZ) using High Frequency Surface Wave Radar (HFSWR)[J].Canadian Journal of Remote Sensing,2001,27(4):354-360.

    [3] 韓樹宗,王海龍,郭佩芳,等.高頻地波雷達(dá)反演海浪的海上對比驗(yàn)證方法研究[J].青島海洋大學(xué)學(xué)報(bào):自然科學(xué)版,2002,32(5):687-694.Han ShuZong,Wang HaiLong,Guo Peifang,et al.A technological study of sensing ocean hydrological parameters with High Frequency Ground Wave Radar[J].Journal of Ocean University of Qingdao,2002,32(5):687-694.

    [4] 秦華,周沫,察豪,等.基于GPU加速的雷達(dá)信號處理并行技術(shù)[J].艦船科學(xué)技術(shù),2013,35(7):77-82.Qin Hua,Zhou Mo,Cha Hao,Sheng Kuo.Research on the Parallel Technology of GPU Acceleration on Radar Signal Processing[J].Ship Science and Technology,2013,35(7):77-82.

    [5] 周廣超.高頻雷達(dá)信號處理機(jī)系統(tǒng)設(shè)計(jì)及基本信號處理軟件研制[D].哈爾濱工業(yè)大學(xué),2006.Zhou Guang Chao.System Design of HFRadar Signal Processor and Developing of Basal Signal Processing Software[D].Harbin:Harbin Institute of Technology,2006.

    [6] 王報(bào)華.基于高頻雷達(dá)的信號處理板設(shè)計(jì)[D].哈爾濱工業(yè)大學(xué),2011.Wang BaoHua.Designof Signal Processing Board Based on HF Radar[D].Harbin:Harbin Institute of Technology,2011.

    [7] 劉霜.基于FPGA的軟件雷達(dá)脈壓處理器的實(shí)現(xiàn)[D].哈爾濱工業(yè)大學(xué),2006.Liu Shuang.Fpga-Based Implementation of Pulse Compression Processor of Software Radar[D].Harbin:Harbin Institute of Technology,2006.

    [8] 曾美艷.高頻地波雷達(dá)實(shí)時(shí)信號處理軟件的設(shè)計(jì)與實(shí)現(xiàn)[D].哈爾濱工業(yè)大學(xué),2008.Zeng MeiYan.Design And Realization of Real-Time Signal Processing Software in High Frequency Ground Wave Radar[D].Harbin Institute of Technology,2008.

    [9] 俞驚雷,柳彬,王開志,等.一種基于GPU的高效合成孔徑雷達(dá)信號處理器[J].信息與電子工程,2010,08(4):415-419.Yu JingLei,LiuBin,WangKaiZhi,LiuXinZhao,Yu WenXian.A high Efficient GPU-Based Signal Processor of Synthetic Aperture Radar[J].Information And Electronic Engineering,2010,08(4):415-419.

    [10] 曹成忠.基于GPU+CPU軟件化雷達(dá)終端系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn)[D].武漢理工大學(xué),2014.Cao ChengZhong.The Design And Implementation of Radar System[D].Wuhan University of Technology,2014.

    [11] 李曉波,關(guān)欣,仲利華,等.基于GPU的外輻射源雷達(dá)信號處理實(shí)時(shí)實(shí)現(xiàn)方法[J].系統(tǒng)工程與電子技術(shù),2014,36(11):2192-2198.Li Xiao Bo,Guan Xin,Zhong Li Hua,et al.Real-time implementation of signal processing for passive radars Based on GPU[J].Systems Engineering and Electronics,2014,36(11):2192-2198.

    責(zé)任編輯 陳呈超

    Real-Time Signal Parallel Processing Scheme and Realization of HF Ground Wave Radar Based on CPU+GPU Architecture

    NIU Jiong,SHEN Jin-Jin,CHEN Qian

    (College of Engineering,Ocean University of China,Qingdao 266100,China)

    With the development of High Frequency Surface Wave Radar(HFSWR) technologies in recent years,the huge amount of remote sensing data brings challenges for real-time imaging processing.The mainstream design scheme of High frequency ground wave signal processor is usually DSP+FPGA technology or based on high performance CPU technology to achieve real-time signal processing,but the hardware and software of system using DSP+FPGA is more compact ,it is difficult to maintain because of long development cycle.High performance CPU technology requires high performance on the hardware,it greatly increasing the cost of hardware.In order to shorten the development cycle,improve the maintainability of the system and reduce the development cost under the premise of guaranteeing signal real-time processing,this paper designs a signal processing scheme based on CPU+GPU architecture,realizes the basic flow of real-time signal processing of HF Ground wave Radar,including data pre-processing,the distance information and speed information of the aim are acquired.Experiment results shows that the proposed signal processing scheme under CPU+GPU architecture has a speedup of 10 times compared with the CPU platform of the same period,which verifies the feasibility of the scheme.

    CPU+GPU; HFSWR; real-time signal processing; cooperative work

    海洋公益性行業(yè)科研專項(xiàng)項(xiàng)目(201505002);國家自然科學(xué)基金重點(diǎn)項(xiàng)目(61132005)資助

    Supported by the Public Science and Technology Research Funds Projects of Ocean(201505002); National Natural Science Foundation of China(61132005)

    2016-09-15;

    2016-11-20

    牛炯(1982-),男,工程師。E-mail:niu_buaa@163.com

    TN957.74

    A

    1672-5174(2017)02-056-06

    10.16441/j.cnki.hdxb.20160238

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

    NIU Jiong,SHEN Jin-Jin,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.

    猜你喜歡
    信號處理線程預(yù)處理
    《信號處理》征稿簡則
    信號處理(2018年5期)2018-08-20 06:16:02
    《信號處理》第九屆編委會
    信號處理(2018年5期)2018-08-20 06:16:00
    《信號處理》征稿簡則
    信號處理(2018年8期)2018-07-25 12:25:42
    《信號處理》第九屆編委會
    信號處理(2018年8期)2018-07-25 12:24:56
    基于預(yù)處理MUSIC算法的分布式陣列DOA估計(jì)
    淺談linux多線程協(xié)作
    淺談PLC在預(yù)處理生產(chǎn)線自動化改造中的應(yīng)用
    絡(luò)合萃取法預(yù)處理H酸廢水
    基于自適應(yīng)預(yù)處理的改進(jìn)CPF-GMRES算法
    Linux線程實(shí)現(xiàn)技術(shù)研究
    亚洲一区高清亚洲精品| 深爱激情五月婷婷| 国产女主播在线喷水免费视频网站 | 日日撸夜夜添| 在线免费观看的www视频| 国产亚洲午夜精品一区二区久久 | av国产免费在线观看| 精品一区二区三区人妻视频| 成人毛片a级毛片在线播放| 亚洲精品久久久久久婷婷小说 | 91在线精品国自产拍蜜月| 毛片一级片免费看久久久久| 亚洲国产最新在线播放| 久99久视频精品免费| 欧美丝袜亚洲另类| 国产91av在线免费观看| 国产精品麻豆人妻色哟哟久久 | 亚洲精华国产精华液的使用体验| 亚洲av成人精品一区久久| 免费在线观看成人毛片| 国产一级毛片在线| 午夜精品国产一区二区电影 | 精品国产三级普通话版| 国产精品人妻久久久影院| 九九热线精品视视频播放| 免费观看性生交大片5| 最后的刺客免费高清国语| 亚洲欧美一区二区三区国产| 99久国产av精品| 十八禁国产超污无遮挡网站| 老师上课跳d突然被开到最大视频| 精品久久久久久久久久久久久| 成人鲁丝片一二三区免费| 国产久久久一区二区三区| 日韩欧美国产在线观看| 可以在线观看毛片的网站| 国产探花极品一区二区| 插逼视频在线观看| 国产三级中文精品| 久久这里只有精品中国| 亚洲欧洲日产国产| 卡戴珊不雅视频在线播放| 日韩亚洲欧美综合| 又爽又黄无遮挡网站| 真实男女啪啪啪动态图| 中国美白少妇内射xxxbb| 在线a可以看的网站| 久久久亚洲精品成人影院| 日本色播在线视频| 99在线人妻在线中文字幕| 99国产精品一区二区蜜桃av| 亚洲精品影视一区二区三区av| 午夜福利视频1000在线观看| 久久久色成人| 亚洲精品456在线播放app| 日本黄色片子视频| 欧美xxxx性猛交bbbb| 亚洲国产成人一精品久久久| 99热这里只有是精品50| 91狼人影院| 免费观看在线日韩| 国产精品无大码| 亚洲精品色激情综合| 偷拍熟女少妇极品色| 一本久久精品| 国产亚洲一区二区精品| 久久久精品大字幕| 白带黄色成豆腐渣| 亚洲精品自拍成人| 日韩大片免费观看网站 | 韩国av在线不卡| 国产午夜精品一二区理论片| 深夜a级毛片| 日日干狠狠操夜夜爽| 亚洲国产成人一精品久久久| 日本黄色视频三级网站网址| 欧美性猛交黑人性爽| 韩国av在线不卡| 欧美日韩国产亚洲二区| 久久这里只有精品中国| 97在线视频观看| 日韩成人av中文字幕在线观看| 国产精品国产三级专区第一集| 久久国产乱子免费精品| 婷婷六月久久综合丁香| 三级毛片av免费| 1024手机看黄色片| 国产久久久一区二区三区| 久久久欧美国产精品| 一个人看的www免费观看视频| 日产精品乱码卡一卡2卡三| 精品久久久久久久久亚洲| 久久久久精品久久久久真实原创| 国产精品,欧美在线| 在线观看美女被高潮喷水网站| 日本欧美国产在线视频| 女的被弄到高潮叫床怎么办| 大又大粗又爽又黄少妇毛片口| 三级经典国产精品| 久久久午夜欧美精品| 亚洲一级一片aⅴ在线观看| 69人妻影院| 亚洲人与动物交配视频| 中文字幕人妻熟人妻熟丝袜美| 久久久久免费精品人妻一区二区| 午夜福利网站1000一区二区三区| 亚洲精品色激情综合| 免费观看的影片在线观看| 久久久久久久国产电影| 午夜视频国产福利| 联通29元200g的流量卡| 听说在线观看完整版免费高清| 国产亚洲最大av| 禁无遮挡网站| 如何舔出高潮| 亚洲欧美成人精品一区二区| 亚洲国产高清在线一区二区三| 哪个播放器可以免费观看大片| 国产精品人妻久久久影院| 我要搜黄色片| 99在线视频只有这里精品首页| 99热这里只有精品一区| 国产成人午夜福利电影在线观看| 九九在线视频观看精品| eeuss影院久久| 久久久精品94久久精品| 男人的好看免费观看在线视频| 中文字幕亚洲精品专区| 免费av观看视频| 日韩欧美 国产精品| 久久久成人免费电影| av视频在线观看入口| 99热这里只有是精品在线观看| 啦啦啦韩国在线观看视频| 亚洲欧美精品自产自拍| 精品久久久久久电影网 | 18+在线观看网站| 国产一区二区三区av在线| 日韩av在线大香蕉| 国产精品电影一区二区三区| 自拍偷自拍亚洲精品老妇| 两个人视频免费观看高清| 久久人妻av系列| 久久热精品热| 男女视频在线观看网站免费| 国产成人午夜福利电影在线观看| 少妇被粗大猛烈的视频| 一区二区三区乱码不卡18| 成人一区二区视频在线观看| 日韩精品有码人妻一区| 久久久久网色| 白带黄色成豆腐渣| 欧美3d第一页| 国产av码专区亚洲av| 久久人人爽人人爽人人片va| 青春草国产在线视频| www.av在线官网国产| 日本wwww免费看| 日本欧美国产在线视频| 久久久久免费精品人妻一区二区| 国产精品一及| 蜜桃亚洲精品一区二区三区| 国产女主播在线喷水免费视频网站 | 舔av片在线| 18禁动态无遮挡网站| 在线天堂最新版资源| 国产一区二区亚洲精品在线观看| 亚洲av不卡在线观看| 亚洲欧美精品自产自拍| 国产免费视频播放在线视频 | 成人毛片60女人毛片免费| av在线天堂中文字幕| 日本三级黄在线观看| 1000部很黄的大片| 丝袜美腿在线中文| 欧美变态另类bdsm刘玥| 欧美bdsm另类| 狂野欧美白嫩少妇大欣赏| 国产毛片a区久久久久| av又黄又爽大尺度在线免费看 | 两个人的视频大全免费| 日日摸夜夜添夜夜爱| 日韩视频在线欧美| 成人国产麻豆网| 麻豆成人午夜福利视频| 黄色日韩在线| 18禁在线播放成人免费| 99热精品在线国产| 久久久成人免费电影| 最近2019中文字幕mv第一页| 免费av不卡在线播放| 麻豆国产97在线/欧美| 国产精品蜜桃在线观看| ponron亚洲| 成人午夜高清在线视频| 啦啦啦韩国在线观看视频| 亚洲精品成人久久久久久| 国产亚洲精品久久久com| 欧美高清成人免费视频www| 久久99热这里只有精品18| 亚洲成人久久爱视频| 黄色一级大片看看| 国产高清国产精品国产三级 | 精品午夜福利在线看| 51国产日韩欧美| 免费观看精品视频网站| av在线蜜桃| 91在线精品国自产拍蜜月| 夜夜爽夜夜爽视频| 色综合亚洲欧美另类图片| 最近2019中文字幕mv第一页| 亚洲图色成人| 精品人妻偷拍中文字幕| 一边摸一边抽搐一进一小说| 搞女人的毛片| 日本猛色少妇xxxxx猛交久久| 91精品国产九色| 两个人的视频大全免费| 亚洲va在线va天堂va国产| 尾随美女入室| 亚洲欧洲日产国产| 久久久久久久亚洲中文字幕| 最新中文字幕久久久久| 99热全是精品| 成人毛片a级毛片在线播放| 久久韩国三级中文字幕| 搞女人的毛片| a级毛片免费高清观看在线播放| 亚洲精品乱码久久久v下载方式| 又粗又爽又猛毛片免费看| 日韩av在线大香蕉| 日韩成人伦理影院| 麻豆国产97在线/欧美| 老司机福利观看| 国产黄片美女视频| 国产精品99久久久久久久久| 国产精品一区二区三区四区免费观看| 国产精品久久久久久久电影| 天堂√8在线中文| 午夜免费激情av| 极品教师在线视频| 欧美人与善性xxx| 国产精品综合久久久久久久免费| 淫秽高清视频在线观看| 韩国高清视频一区二区三区| 亚洲在线自拍视频| 国产高清不卡午夜福利| 国产女主播在线喷水免费视频网站 | 直男gayav资源| 寂寞人妻少妇视频99o| 亚洲18禁久久av| 能在线免费看毛片的网站| 国产精品久久电影中文字幕| 国产伦一二天堂av在线观看| av免费观看日本| 久久99热这里只频精品6学生 | 欧美日韩综合久久久久久| 色吧在线观看| 国产精品一区www在线观看| 看黄色毛片网站| 久久99精品国语久久久| 少妇裸体淫交视频免费看高清| 免费看日本二区| 夫妻性生交免费视频一级片| 久久久久久伊人网av| 1024手机看黄色片| 麻豆成人av视频| 国产精品av视频在线免费观看| 18禁在线无遮挡免费观看视频| 日本一本二区三区精品| 日韩中字成人| 男人狂女人下面高潮的视频| 日韩一本色道免费dvd| 一区二区三区免费毛片| 免费无遮挡裸体视频| 我要看日韩黄色一级片| a级一级毛片免费在线观看| 亚洲av男天堂| 别揉我奶头 嗯啊视频| 国产真实伦视频高清在线观看| 麻豆成人午夜福利视频| 国内精品美女久久久久久| 欧美成人a在线观看| 国产精品99久久久久久久久| 青春草亚洲视频在线观看| 国产又黄又爽又无遮挡在线| 婷婷色综合大香蕉| 国产一级毛片在线| 深夜a级毛片| 热99re8久久精品国产| 亚洲电影在线观看av| 六月丁香七月| 中国美白少妇内射xxxbb| 日本猛色少妇xxxxx猛交久久| 老女人水多毛片| eeuss影院久久| 美女高潮的动态| 91午夜精品亚洲一区二区三区| 免费观看的影片在线观看| 国产黄片美女视频| 91精品伊人久久大香线蕉| 高清视频免费观看一区二区 | 非洲黑人性xxxx精品又粗又长| 亚洲av日韩在线播放| 欧美成人一区二区免费高清观看| 岛国毛片在线播放| 日日摸夜夜添夜夜爱| 高清视频免费观看一区二区 | 久久99热6这里只有精品| 日本免费一区二区三区高清不卡| 国产淫片久久久久久久久| 熟女电影av网| 一区二区三区乱码不卡18| 三级经典国产精品| 六月丁香七月| 在线a可以看的网站| 国产免费视频播放在线视频 | 在线观看一区二区三区| 欧美极品一区二区三区四区| 26uuu在线亚洲综合色| 又爽又黄无遮挡网站| 蜜臀久久99精品久久宅男| 能在线免费观看的黄片| 男女啪啪激烈高潮av片| 精品无人区乱码1区二区| 国产高潮美女av| 精品欧美国产一区二区三| av.在线天堂| 成人无遮挡网站| 亚洲色图av天堂| 最近中文字幕高清免费大全6| 中国美白少妇内射xxxbb| 久久婷婷人人爽人人干人人爱| 午夜爱爱视频在线播放| 日韩欧美精品免费久久| 国产黄色视频一区二区在线观看 | 日本色播在线视频| 99久久精品国产国产毛片| 看黄色毛片网站| 国产伦精品一区二区三区视频9| 水蜜桃什么品种好| h日本视频在线播放| 免费电影在线观看免费观看| 乱系列少妇在线播放| 亚洲,欧美,日韩| 日日撸夜夜添| 亚洲国产色片| 亚洲最大成人中文| 国产一区亚洲一区在线观看| 国产亚洲一区二区精品| 中文亚洲av片在线观看爽| or卡值多少钱| 级片在线观看| 日产精品乱码卡一卡2卡三| 久久午夜福利片| kizo精华| 国产精品久久久久久精品电影| 三级国产精品片| 全区人妻精品视频| 美女xxoo啪啪120秒动态图| 国产又黄又爽又无遮挡在线| 免费av毛片视频| 欧美成人午夜免费资源| 色综合亚洲欧美另类图片| 国产精品永久免费网站| 亚洲一级一片aⅴ在线观看| 99热6这里只有精品| 我的老师免费观看完整版| 亚洲最大成人手机在线| 边亲边吃奶的免费视频| 精品人妻熟女av久视频| 永久免费av网站大全| 亚洲四区av| 亚洲中文字幕一区二区三区有码在线看| 99久久精品国产国产毛片| 亚洲第一区二区三区不卡| 亚洲四区av| 亚洲av日韩在线播放| 久久久久网色| 亚洲av中文av极速乱| 人妻系列 视频| 日产精品乱码卡一卡2卡三| 国产视频内射| 日韩中字成人| 老司机福利观看| 少妇丰满av| 国产久久久一区二区三区| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲av中文av极速乱| 春色校园在线视频观看| 国产高清三级在线| 99久久人妻综合| 亚洲电影在线观看av| 久久久久久久午夜电影| 日本黄色片子视频| 欧美三级亚洲精品| or卡值多少钱| 又黄又爽又刺激的免费视频.| 亚洲欧美日韩无卡精品| av女优亚洲男人天堂| 国产黄片美女视频| 国产精品精品国产色婷婷| 国产一区有黄有色的免费视频 | 国内精品一区二区在线观看| 久久久久久久午夜电影| 国产大屁股一区二区在线视频| av国产免费在线观看| 最近的中文字幕免费完整| 熟女电影av网| 女人十人毛片免费观看3o分钟| 免费av观看视频| 久久人人爽人人爽人人片va| 久久久久久久久久久丰满| 午夜福利在线在线| 国产亚洲一区二区精品| 久久精品久久久久久噜噜老黄 | 欧美一区二区国产精品久久精品| av在线亚洲专区| 国产黄片视频在线免费观看| 成人亚洲欧美一区二区av| 亚洲欧美清纯卡通| 久久久亚洲精品成人影院| 久久99热6这里只有精品| 日本午夜av视频| 高清毛片免费看| 国产私拍福利视频在线观看| 一级毛片aaaaaa免费看小| 淫秽高清视频在线观看| 黄色配什么色好看| 亚洲婷婷狠狠爱综合网| 一级黄色大片毛片| 免费观看人在逋| 美女被艹到高潮喷水动态| 亚洲av成人精品一二三区| 婷婷色av中文字幕| 日日啪夜夜撸| 只有这里有精品99| 狠狠狠狠99中文字幕| 黄色欧美视频在线观看| 欧美日韩精品成人综合77777| 国产视频内射| 日本免费一区二区三区高清不卡| 亚洲婷婷狠狠爱综合网| 久久久久久久久久成人| 激情 狠狠 欧美| 亚洲人成网站在线观看播放| 一边亲一边摸免费视频| 久久草成人影院| 麻豆久久精品国产亚洲av| 最后的刺客免费高清国语| 伦理电影大哥的女人| or卡值多少钱| 国产一区有黄有色的免费视频 | 日本免费一区二区三区高清不卡| 男人狂女人下面高潮的视频| 久久久久久伊人网av| 啦啦啦观看免费观看视频高清| 激情 狠狠 欧美| 国产精品国产高清国产av| 国产在视频线精品| 国产三级中文精品| 亚洲人成网站在线播| 国产精品一及| 18+在线观看网站| 亚洲av熟女| 亚洲电影在线观看av| 狂野欧美白嫩少妇大欣赏| 免费观看a级毛片全部| 中文字幕人妻熟人妻熟丝袜美| 有码 亚洲区| 内射极品少妇av片p| 性插视频无遮挡在线免费观看| 国产成人精品久久久久久| 99热全是精品| 男女啪啪激烈高潮av片| 日本黄色片子视频| 亚洲av.av天堂| 免费一级毛片在线播放高清视频| 最近中文字幕高清免费大全6| 欧美区成人在线视频| 人妻少妇偷人精品九色| 97人妻精品一区二区三区麻豆| 午夜久久久久精精品| 久久久久久国产a免费观看| 日韩,欧美,国产一区二区三区 | 国产av在哪里看| 九九爱精品视频在线观看| 亚洲av成人av| 天天一区二区日本电影三级| 国产精华一区二区三区| 2022亚洲国产成人精品| 久久久久久久久久久丰满| 精品久久久久久久久av| www.av在线官网国产| 国产精品熟女久久久久浪| 久久99蜜桃精品久久| 精品久久久久久久久久久久久| 国产精品三级大全| 岛国毛片在线播放| av女优亚洲男人天堂| 日本色播在线视频| 国产成人精品一,二区| 深爱激情五月婷婷| 中文亚洲av片在线观看爽| 中文字幕精品亚洲无线码一区| 黄片wwwwww| 中文字幕久久专区| 一夜夜www| 欧美三级亚洲精品| 免费看av在线观看网站| 免费不卡的大黄色大毛片视频在线观看 | 久久久久久久久久成人| 最近手机中文字幕大全| 亚洲乱码一区二区免费版| 国产精品永久免费网站| 国产精品爽爽va在线观看网站| 国产黄片美女视频| 51国产日韩欧美| 亚洲欧美成人综合另类久久久 | 一夜夜www| 欧美三级亚洲精品| 免费看日本二区| 五月玫瑰六月丁香| 亚洲最大成人中文| 国产真实乱freesex| 26uuu在线亚洲综合色| 久久精品熟女亚洲av麻豆精品 | 一本一本综合久久| av国产免费在线观看| 啦啦啦观看免费观看视频高清| 最近中文字幕高清免费大全6| 国产免费福利视频在线观看| 久热久热在线精品观看| 一边亲一边摸免费视频| 99久久九九国产精品国产免费| 午夜爱爱视频在线播放| 亚洲美女视频黄频| 国产在线男女| 97热精品久久久久久| 亚洲真实伦在线观看| 老女人水多毛片| 超碰av人人做人人爽久久| 国产午夜福利久久久久久| 国产一级毛片在线| 九色成人免费人妻av| 成人高潮视频无遮挡免费网站| 国产片特级美女逼逼视频| 欧美bdsm另类| 亚洲人成网站在线观看播放| 成人无遮挡网站| 亚洲精品,欧美精品| 十八禁国产超污无遮挡网站| 一个人免费在线观看电影| 永久网站在线| 亚洲欧美精品综合久久99| 亚洲av福利一区| www.av在线官网国产| 国产成人aa在线观看| 草草在线视频免费看| 亚洲成av人片在线播放无| 五月玫瑰六月丁香| 日韩亚洲欧美综合| av在线蜜桃| 成人av在线播放网站| 国产精品国产三级专区第一集| 国产单亲对白刺激| 哪个播放器可以免费观看大片| 蜜桃亚洲精品一区二区三区| 国产伦精品一区二区三区视频9| 色尼玛亚洲综合影院| 精品一区二区免费观看| 观看免费一级毛片| 亚洲18禁久久av| 欧美高清成人免费视频www| 一级毛片久久久久久久久女| 不卡视频在线观看欧美| 午夜亚洲福利在线播放| 久久鲁丝午夜福利片| 国产精品久久久久久av不卡| 高清午夜精品一区二区三区| 成人性生交大片免费视频hd| 亚洲精品久久久久久婷婷小说 | 热99re8久久精品国产| 亚洲欧美日韩无卡精品| 联通29元200g的流量卡| 18禁在线无遮挡免费观看视频| 久久99精品国语久久久| 国产免费又黄又爽又色| 性插视频无遮挡在线免费观看| av在线观看视频网站免费| 亚洲在线观看片| 男女视频在线观看网站免费| 自拍偷自拍亚洲精品老妇| 啦啦啦观看免费观看视频高清| 天天躁日日操中文字幕| 久久99蜜桃精品久久| 免费观看精品视频网站| 一区二区三区免费毛片| 深夜a级毛片| 免费一级毛片在线播放高清视频| 日韩高清综合在线| 日本午夜av视频| 久久久欧美国产精品| 午夜精品在线福利| 国内揄拍国产精品人妻在线| 中国美白少妇内射xxxbb| 久久久国产成人免费| 久久久久久伊人网av| 少妇熟女aⅴ在线视频| 99热6这里只有精品| 国产男人的电影天堂91| 免费无遮挡裸体视频| 欧美3d第一页| 国产精品日韩av在线免费观看| 国产熟女欧美一区二区|