基于FPGA的多通道寬帶數(shù)字頻譜儀設(shè)計

劉 欣劉東亮

（①國家無線電監(jiān)測中心 北京100037 ②中國科學(xué)院國家天文臺 北京 100012）

基于FPGA的多通道寬帶數(shù)字頻譜儀設(shè)計

劉 欣①劉東亮②

（①國家無線電監(jiān)測中心 北京100037 ②中國科學(xué)院國家天文臺 北京 100012）

為了進一步提高數(shù)據(jù)處理效率和滿足射電天文后端設(shè)備發(fā)展與類似密集型數(shù)據(jù)處理的需求，文章給出了一種基于FPGA的多通道寬帶數(shù)字頻譜儀系統(tǒng)仿設(shè)計，具有16K通道數(shù)，2GHz的處理帶寬，可以達到128kHz的頻譜分辨率。主要介紹了頻譜儀的設(shè)計思路，給出CAD設(shè)計方案、生成模型驗證和仿真結(jié)果測試。

數(shù)字頻譜儀；FPGA；多通道數(shù)

引言

頻譜儀被廣泛使用在科研工作中，是進行頻譜轉(zhuǎn)換與測量的專用設(shè)備，用來對各種信號進行處理分析。近年來，隨著半導(dǎo)體技術(shù)的發(fā)展，寬帶高速ADC、高性能數(shù)字信號處理芯片（FPGA）以及高性能計算機總線的出現(xiàn)，使基于FFT原理的數(shù)字頻譜技術(shù)迅速發(fā)展起來。

目前國內(nèi)外比較常見的頻譜儀系統(tǒng)是基于Xilinx Virtex系列的FPGA的數(shù)字信號處理系統(tǒng)[1]。例如美國加州大學(xué)伯克利分校近年來也在研制類似的頻譜儀，其功能包括高速采樣量化信號、頻帶調(diào)制、多相濾波、傅里葉變換和對數(shù)字信號的高速傳輸（通過兩個10Gb網(wǎng)絡(luò)端口（CX4）將數(shù)據(jù)以UDP包形式輸出），輸出端可接駁BEE2系統(tǒng)、交換機和高速網(wǎng)卡等，這一系列頻譜儀的主要特點是具有多個Z-DOK接口,可以接駁多種I/O板（包括雙路1GSa/s的采樣卡或者四路250MSa/s的模數(shù)，數(shù)模轉(zhuǎn)換卡)[2]。與此同時，德國馬普研究所近來也在研發(fā)一種擴展帶寬快速傅里葉變換頻譜儀系統(tǒng)[3]，該系統(tǒng)采用賽靈思最新的Xilinx Virtex-6可編程邏輯門電路（FPGA）處理核心，配備實時帶寬單通道4GHz、雙通道2.5GHz的數(shù)字轉(zhuǎn)換和采樣模塊，實現(xiàn)了16K傅里葉變換通道的高分辨率數(shù)字頻譜分析能力。

綜合以上特點，結(jié)合目前Xilinx最新發(fā)展的Virtex系列FPGA技術(shù)發(fā)展，進行了多通道寬帶數(shù)字頻譜儀系統(tǒng)仿真設(shè)計。

1 頻譜儀設(shè)計框架

頻譜儀的主要原理是將接收輸入的基帶混合中頻(IF)信號進行外差式采樣，將輸入的信號進行量化與編碼的數(shù)字測量后，通過數(shù)據(jù)總線傳輸至FPGA芯片核心，進行多相濾波與快速傅里葉變換運算的數(shù)字信號處理，然后計算信號的功率譜數(shù)據(jù)，得到信號測量結(jié)果，并按照X軸為頻率，Y軸為功率譜的數(shù)據(jù)矩陣，將測量結(jié)果存儲在頻譜儀內(nèi)置的存儲芯片中。通過頻譜儀的網(wǎng)絡(luò)接口，可將數(shù)據(jù)輸出至后端接收電腦，實現(xiàn)進一步的測量處理與圖形繪制等工作。頻譜儀機柜將多套頻譜儀整合封裝，并設(shè)置統(tǒng)一的電源與時鐘模塊。通用的頻譜儀設(shè)計框圖如圖1所示。

圖1 頻譜儀設(shè)計框圖

2 多通道頻譜儀的具體設(shè)計

不同于傳統(tǒng)的Verilog語言底層設(shè)計，Xilinx近年來引入系統(tǒng)級模塊化設(shè)計的概念，在這種方式下，不僅可以更直觀、快速和靈活地進行系統(tǒng)級模塊的構(gòu)建，而且可以使設(shè)計者把更多精力集中在核心算法的實現(xiàn)上，而不必拘泥于具體實現(xiàn)方式。

在Matlab Simulink平臺下，調(diào)用Xilinx System Generator函數(shù)庫和Xilinx函數(shù)庫，依次連接ADC、PFB和FFT模塊。ADC模塊的時鐘頻率為2048MHz（1024MHz奈奎斯特采樣率），ADC四路同步輸出，因此全局時鐘為1/4的ADC頻率，512MHz；PFB和FFT模塊中設(shè)置大小為16K通道輸出，之后連接實部虛部的分離模塊和系數(shù)量化模塊，將FFT輸出的18bit按比特位選擇其中的8bit，并將輸出合并為128bit的總數(shù)據(jù)鏈路。設(shè)計如圖2所示。

圖2 Simulink中頻譜儀的ADC與FFT模塊設(shè)計

圖3 基于窗函數(shù)的多相濾波器設(shè)計

傳統(tǒng)的FFT處理過程中會產(chǎn)生一定的頻譜泄露，為了實現(xiàn)更好的相鄰頻帶頻響特性，在設(shè)計中加入了基于窗函數(shù)的相濾波器組，經(jīng)過FFT處理過的每個通道頻譜得到最大信噪比。

如圖3所示，整體設(shè)計調(diào)用Xilinx自帶FIR濾波器模塊進行拼接，F(xiàn)IR采用4抽頭漢寧窗算法。通過多相FIR濾波器對輸入的信號進行預(yù)譜帶處理后，傳輸?shù)胶蠖说?024通道FFT模塊中。

如果直接將128bit總數(shù)據(jù)鏈路輸入到傳輸模塊，則無法控制每個IP數(shù)據(jù)包中對應(yīng)的頻譜信息，因此設(shè)計中Bram寄存器對數(shù)據(jù)流進行重組。此模塊的功能為順序輸入數(shù)據(jù)，按預(yù)存次序輸出。將16K通道對應(yīng)分配到16IP中，每個IP接收對應(yīng)的1024個通道，通過數(shù)據(jù)流調(diào)整模塊和IP數(shù)量的對應(yīng)控制，可以實現(xiàn)對信號頻譜的任意目標分配，非常靈活。

3 性能測試和分析

整體設(shè)計中，由于采用了FIR4抽頭漢寧窗算法的多相濾波器組與FFT相結(jié)合的頻譜儀處理核心，使得其多通道頻譜輸出旁帶得到很好的抑制。比較而言，漢寧窗的主瓣與海明窗相當，比同階的矩形窗寬，旁瓣最大泄漏為-43dB，矩形窗為-21dB，漢寧為-53dB。漢寧窗旁瓣隨其階數(shù)的增加以-5.2dB/個的速度衰減，為最快，海明窗與矩形窗則都以-1.3dB/個的速度在下降。海明窗的旁瓣衰減略比漢寧窗大，但是隨旁瓣增加其衰減很慢。因此,實際信號分析中常常選用漢寧窗，只要選擇恰當?shù)碾A數(shù)，使用漢寧窗截斷的時候，頻譜泄漏可以更小，這樣就可以有效改善FFT輸出的頻率響應(yīng)，抑制頻譜泄露。

本研究在Matlab中對實際設(shè)計進行了測試，如圖4所示，在加入多相濾波器組后的FFT輸出旁瓣抑制近-65dB，較之單獨FFT輸出的-35db旁瓣，旁帶干擾大大降低。

圖4 多相4抽頭FIR濾波器響應(yīng)

對于頻譜儀整體響應(yīng)，在Matlab Simulink平臺下，給系統(tǒng)輸入疊加脈沖信號的白噪聲輸入進行系統(tǒng)測試。頻譜儀輸出測試結(jié)果見圖5，從圖中可以看出在對應(yīng)的通道中得到了疊加信號頻譜響應(yīng)，驗證頻譜儀設(shè)計仿真工作正常。

圖5 白噪聲譜仿真測試結(jié)果

4 結(jié)論

穩(wěn)定高速的實時寬帶數(shù)字頻譜儀系統(tǒng)可以提高整個數(shù)據(jù)接收系統(tǒng)的信號傳輸性能,在此需求下本文重點介紹了頻譜儀的設(shè)計思路，給出CAD設(shè)計方案，生成模型驗證，進行了仿真。從測試結(jié)果來看，達到了在實際應(yīng)用中的設(shè)計要求。在此基礎(chǔ)上，還對比驗證了加窗后的FFT處理，可以有效地降低通道間的頻譜泄露。在后續(xù)工作中，將在現(xiàn)有初步研究的基礎(chǔ)上進行更進一步的系統(tǒng)分析和測試,針對特殊環(huán)境使用的可靠性等方面進行測試。

[1]Parsons,A.,et al.A Scalable Correlator Architecture Based on Modular FPGA Hardware and Data Packetization[J].PASP, 2008,120(873):1207-1221.

[2]A.Parsons,D.Backer,C.Chang,et.al.PetaOp/Second FPGA Signal Processing for SETI and Radio Astronomy[J].Proc. 10th Asilomar Conference on Signals,Systems and Computers,Paci_c Grove,CA,November 2006.

[3]http://www.radiometer-physics.de/rpg/html/Products_Spectrometers.html[EB/OL].

Design the Multi-Channel Wideband Digital Spectrometer based on FPGA

Liu-Xin①Liu Dong-liang②
（①The State Radio-monitoring Center Testing Center,Beijing 100041;②National Astronomical Observatories,CAS,Beijing 100012）

The FPGA based multi-channel wideband digital spectrometer was designed for the further improving the efficiency of data processing and fulfilling the developmental needs on the radio astronomy back-end equipments and similar intensive data processing.The multi-channel wideband digital spectrometer designed has 16K multi-channel with the bandwidth of 2GHz generating a spectral resolution up to 128kHz.We introduce the design frame,give the CAD design scheme and generate the model validation and the simulation results test.

digital spectrometer；FPGA；Multi-channel

TM935.21

A

1005-5738(2014)02-113-04

[責(zé)任編輯：周曉艷]

2014-09-24

2013年度國家自然科學(xué)基金青年項目“基于ROACH高性能射電望遠鏡頻譜分析系統(tǒng)研制”階段性成果，項目號：11203040

劉欣，男，漢族，山東德州人，國家無線電監(jiān)測中心高級工程師，主要研究方向為無線電監(jiān)測、數(shù)字信號處理。