鄧 嵐 康懷祺 李宏波 鄢 林

(電子科技大學(xué) 成都 611731)

1 引言

針對(duì)數(shù)字化雷達(dá)接收機(jī)采樣前端，目前市面上已經(jīng)出現(xiàn)了很多超高速多比特的A/D采樣芯片，其采樣速度可達(dá)到幾個(gè)Gsps，而由此帶來的問題就是數(shù)字信號(hào)處理器很難實(shí)時(shí)處理如此高速的數(shù)據(jù)流，數(shù)字下變頻技術(shù)就顯得至關(guān)重要。數(shù)字下變頻的基本功能是將輸入的寬帶中頻信號(hào)下變頻為數(shù)字基帶信號(hào)，并轉(zhuǎn)換成較低的數(shù)據(jù)流。常規(guī)的數(shù)字下變頻結(jié)構(gòu)［1］，采樣后的數(shù)據(jù)通過與數(shù)控振蕩NCO產(chǎn)生的解調(diào)信號(hào)相乘，得到兩路正交的信號(hào)，同時(shí)實(shí)現(xiàn)頻譜搬移，再由低通濾波器濾波后，進(jìn)行多倍抽取，達(dá)到降速的效果。對(duì)于寬帶信號(hào)，經(jīng)過高速A/D采樣后，數(shù)據(jù)速率可達(dá)幾Gsps，由于常規(guī)的數(shù)字下變頻方法，工作速率與ADC輸出數(shù)據(jù)的速率一致，而FPGA器件無法達(dá)到如此高的運(yùn)行速度，因此采用常規(guī)的數(shù)字下變頻方法無法實(shí)現(xiàn)對(duì)寬帶中頻信號(hào)的處理。

為克服常規(guī)的數(shù)字下變頻方法的缺陷，國(guó)外Jeffrey O.Coleman，James J.Alter和Dan Scholnik 等人提出了利用希爾伯特濾波器實(shí)現(xiàn)頻譜搬移的方法［2］。信號(hào)通過帶通濾波器濾去復(fù)信號(hào)，再通過Hilbert濾波器，將得到的信號(hào)下變頻為零中頻，最后進(jìn)行D倍抽取降速。這種下變頻方法雖然可以硬件實(shí)現(xiàn)，但是其資源占用量大，系統(tǒng)內(nèi)部設(shè)計(jì)復(fù)雜。

本文為解決上述兩個(gè)數(shù)字下變頻結(jié)構(gòu)的缺陷，提出了一種寬帶中頻信號(hào)的數(shù)字下變頻方法，不僅可以實(shí)現(xiàn)對(duì)超高速率中頻采樣信號(hào)的處理，而且設(shè)計(jì)簡(jiǎn)單，其資源占用量低。

2 系統(tǒng)結(jié)構(gòu)設(shè)計(jì)

本文設(shè)計(jì)的采集系統(tǒng)主要應(yīng)用于某機(jī)載雷達(dá)的數(shù)據(jù)記錄儀中，完成實(shí)時(shí)采集帶寬為300MHz，載頻為450MHz的雷達(dá)回波信號(hào)。根據(jù)帶通采樣定理:

式中，fs為采樣速率;fH為信號(hào)的上限頻率;fL為信號(hào)的下限頻率;f0為信號(hào)的中心頻率。只要n取滿足fs≥2(fH－fL)的最大正整數(shù)(0，1，2，…)，則采樣所得到的信號(hào)采樣值x(nTs)能準(zhǔn)確確定原信號(hào)x(t)。為簡(jiǎn)化系統(tǒng)設(shè)計(jì)，在設(shè)計(jì)中選擇n=0時(shí)的fs，即fs=1.8GHz。

圖1 數(shù)字下變頻I路結(jié)構(gòu)圖

回波信號(hào)經(jīng)過采樣芯片處理，產(chǎn)生4路速率為450Msps的8bit信號(hào)，且依時(shí)序排列。經(jīng)過串并轉(zhuǎn)換模塊，形成16路速率為112.5Msps的并行信號(hào)，送入到DDC模塊中進(jìn)行數(shù)字下變頻處理。首先將16路信號(hào)通過基于多相濾波的數(shù)字正交變換，實(shí)現(xiàn)寬帶中頻信號(hào)下變頻為基帶信號(hào)，產(chǎn)生正交的I、Q兩路數(shù)據(jù)［3］;其次利用多相濾波結(jié)構(gòu)和分布式算法實(shí)現(xiàn)4倍抽取，產(chǎn)生8路相互正交的I信號(hào)和Q信號(hào)，其數(shù)據(jù)率為112.5Msps存入到存儲(chǔ)芯片中。其數(shù)字下變頻I路結(jié)構(gòu)框圖如圖1所示，Q路結(jié)構(gòu)框圖如圖2所示。

圖2 數(shù)字下變頻Q路結(jié)構(gòu)圖

2.1 并行化模塊

在并行化處理過程中，需將4路速率為450Msps的信號(hào)降速為16路速率為112.5Msps的信號(hào)。因此每一路數(shù)據(jù)按照1路串行輸入，4路并行輸出，即1:4進(jìn)行轉(zhuǎn)換。在此模塊中，采用FPGA中串轉(zhuǎn)并原語(yǔ)ISERDES［4］，對(duì)每路信號(hào)進(jìn)行1:4的串并轉(zhuǎn)換。由于在ISERDES原語(yǔ)中，需要4分頻的時(shí)鐘信號(hào)，因此還需利用DCM模塊，對(duì)輸入的時(shí)鐘信號(hào)進(jìn)行4分頻處理。

2.2 正交化模塊

4路信號(hào)經(jīng)過并行化模塊處理，得到的16路112.5MHz的8bit信號(hào)并依時(shí)序排列。

設(shè)輸入信號(hào)為:

其中a(t)為輸入信號(hào)的包絡(luò)，φ(t)為信號(hào)的相位，載頻f0=450MHz，對(duì)其進(jìn)行采樣頻率為fs=4f0帶通采樣，得到采樣序列為:

分別為x［n］的同相分量和正交分量(零中頻信號(hào))，根據(jù)基于多相濾波的正交變換［5］，得到:

為簡(jiǎn)化系統(tǒng)設(shè)計(jì)，對(duì)多相濾波正交化結(jié)構(gòu)進(jìn)行改進(jìn)。將2倍抽取置于乘法器之后，與抽取濾波模塊中的D倍抽取結(jié)構(gòu)相結(jié)合，形成2D倍抽取。其改進(jìn)的多相濾波正交變換為:將xI［2n+1］和xQ［2n］都置零，則x［n］的同相分量和正交分量可以轉(zhuǎn)換為下式:

2.3 抽取濾波模塊

16路數(shù)據(jù)經(jīng)過正交化處理后，形成的16路I和16路 Q數(shù)據(jù)，每路數(shù)據(jù)的傳輸速率仍為112.5Msps。由于傳輸給存儲(chǔ)芯片時(shí)，只有4路I通道和4路Q通道，即每路通道數(shù)據(jù)率為450Msps，其速率過高，后續(xù)模塊無法處理，則需要對(duì)形成的正交16路I和16路Q數(shù)據(jù)進(jìn)行降速處理。在本文中選擇抽取率D=2，即實(shí)現(xiàn)2D=4倍抽取，就可以滿足系統(tǒng)要求，最后產(chǎn)生需要存儲(chǔ)的4路I數(shù)據(jù)和4路Q數(shù)據(jù)，其每路數(shù)據(jù)速率為112.5MHz。

在抽取濾波過程中，首先應(yīng)對(duì)信號(hào)進(jìn)行低通濾波，再進(jìn)行抽取。根據(jù)Noble恒等式知，可先進(jìn)行抽取，再進(jìn)行濾波，更簡(jiǎn)化濾波器的設(shè)計(jì)［6］。本文采用32階的FIR濾波器實(shí)現(xiàn)4倍抽取，即h［n］=h［－n］。而為了節(jié)省資源消耗，提高運(yùn)算速率，采用了多相濾波和改進(jìn)型的DA算法［7］進(jìn)行實(shí)現(xiàn)。

利用Matlab軟件的FDAtool工具，生成一個(gè)32階的第I類FIR等紋波低通濾波器，將其系數(shù)導(dǎo)出。由于 N=32(即 n=0，1，2，…，31)，在進(jìn)行抽取濾波時(shí)，只需要利用 h［n］，n=0，1，2，…，15 這 16 個(gè)系數(shù)進(jìn)行處理。

信號(hào)經(jīng)過多相濾波正交化后，得到:

利用DA算法實(shí)現(xiàn)抽取濾波過程，其DA算法的數(shù)學(xué)表達(dá)式為:

為了節(jié)省資源占用量，采用圖3所示結(jié)構(gòu)，同時(shí)采用塊ROM和分布式ROM，利用位移寄存器的方式，產(chǎn)生8個(gè)地址，進(jìn)入到LUT中查表，得到相應(yīng)結(jié)果。將得到結(jié)果通過流水線寄存器相加，得到最后濾波結(jié)果。

圖3 DA算法結(jié)構(gòu)圖

在塊ROM中，組成LUT的系數(shù)為:

在分布式ROM中，組成LUT的系數(shù)為:

此四組系數(shù)，分別由Matlab生成.coe文件導(dǎo)入到對(duì)應(yīng)ROM中，生成所需查找表。由于LUT中，略去了 xBI［n］=0 和xBQ［n］=0 的處理，且 xBI［n］和xBQ［n］的零點(diǎn)都是間隔出現(xiàn)，則此處便完成了一次2倍抽取的過程。

根據(jù)上述LUT的設(shè)計(jì)，本系統(tǒng)的濾波器等效于一個(gè)16階濾波器，因此4倍抽取多相濾波結(jié)構(gòu)中，每一個(gè)子結(jié)構(gòu)只有16個(gè)輸入。而在LUT設(shè)計(jì)時(shí)，已實(shí)現(xiàn)了2倍抽取，所以在多相濾波抽取過程中，只需完成2倍抽取即可。

為更容易理解，將 xI［n］和xQ［n］中為0的序列去掉，得到新的序列 x'I［n］和x'Q［n］。

在實(shí)現(xiàn)2倍抽取過程中，只需要選x'I［n］和x'Q［n］作為濾波器輸入便可實(shí)現(xiàn)2倍抽取。

3 實(shí)驗(yàn)結(jié)果及分析

選擇Chirp信號(hào)模擬雷達(dá)中頻信號(hào)，利用Matlab對(duì)上述設(shè)計(jì)進(jìn)行仿真，得到如圖4－a所示的chirp信號(hào)頻譜圖和如圖4－b所示的基帶信號(hào)頻譜圖。

圖4 a為chirp信號(hào)頻譜圖，b為基帶信號(hào)頻譜圖

在FPGA實(shí)現(xiàn)過程中，首先將chirp信號(hào)數(shù)據(jù)轉(zhuǎn)換成.txt的數(shù)據(jù)文本文件，以此作為數(shù)據(jù)激勵(lì)源，送入到FPGA中，利用Xilinx公司的xc4vlx25ff668－12 FPGA來實(shí)現(xiàn)數(shù)字下變頻，采用ISE10.1綜合軟件對(duì)其進(jìn)行綜合，采用Modelsim6.3c仿真軟件對(duì)其進(jìn)行仿真［7］。數(shù)字下變頻的時(shí)序仿真圖如圖5所示。FPGA系統(tǒng)資源占用率如表1所示。由此可見，采用FPGA模塊化設(shè)計(jì)能夠較好地對(duì)寬帶中頻信號(hào)進(jìn)行下變頻處理，采用分布式算法，節(jié)省了資源消耗，更利于系統(tǒng)設(shè)計(jì)。

表1 數(shù)字下變頻資源利用表

圖5 數(shù)字下變頻時(shí)序仿真圖

4 結(jié)論

本設(shè)計(jì)通過對(duì)超高速中頻采樣信號(hào)進(jìn)行并行處理，形成多路并行數(shù)據(jù)，降低了每路數(shù)據(jù)的運(yùn)算速度，滿足下變頻處理中FPGA對(duì)工作頻率的要求。通過基于多相濾波的正交化處理，簡(jiǎn)化了下變頻結(jié)構(gòu)，節(jié)省了很大的資源占用量，提高了運(yùn)算速度。除此之外，利用多相濾波結(jié)構(gòu)和改進(jìn)后的DA算法，完成了對(duì)中頻信號(hào)的多倍抽取，實(shí)現(xiàn)了下變頻結(jié)構(gòu)中降低基帶信號(hào)傳輸速率的要求。仿真及實(shí)驗(yàn)結(jié)果表明:本設(shè)計(jì)不僅在軟件上得到了實(shí)現(xiàn)，在硬件上，采用FPGA實(shí)現(xiàn)數(shù)字下變頻延時(shí)短，最高系統(tǒng)主頻可達(dá)到206MHz，輸入輸出時(shí)間延遲為7個(gè)時(shí)鐘周期，滿足實(shí)時(shí)處理要求。

［1］劉紅文.基于FPGA的軟件無線電數(shù)字接收機(jī)的研究［D］.武漢:武漢理工大學(xué)，2009.

［2］Jeffrey O.Coleman，James J.Alter，Dan Scholnik.FPGA Architecture for gigahertz－sampling wideband if－to－baseband conversion［C］.2000 Int’l Conf.on Signal Processing Applications and Technology Dallas TX.2000，10.

［3］Wenhao Zhang，Jun Wang，Yuxi Zhang，Wei Li.Design of RF digital receiver based on FPGA［C］.International Communication Conference on Wireless Mobile and Computing，2009:699－702.

［4］Vitex－4 User Guide.Xilinx.2007，www.xilinx.com .

［5］楊小牛，婁才義，徐建良.軟件無線電技術(shù)與應(yīng)用［M］.北京:北京理工大學(xué)出版社，2010:87－91.

［6］周浩.寬帶數(shù)字接收機(jī)關(guān)鍵技術(shù)研究［D］.南京:南京理工大學(xué)，2008.

［7］Hua-Ming Liu，Guang-Jun Li，Bo Yan，Qiang Li.A 100MHz digital down converter with modified filter for wideband software-defined radios［C］.2010 International Conference on Electronics and Information Engineering，2010:V2540-V2544.