计算机频谱仪仿真,ROACH构架中的多相频谱仪仿真测试

刘东亮++李建玲

摘 要：针对射电天文观测应用需求，结合最新工业界可编程逻辑门阵列(FPGA)技术的发展，研究了基于天文信号处理与电子技术研究集合库(ROACH)构架中的多相宽带频谱仪仿真设计。文章重点介绍了频谱仪的设计流程、Simulink中的设计构架、数据流模块，并进行了仿真测试结果分析。实现并验证了具有16～32K通道数，2～2.5 GHz的处理带宽，可以达到68.5 kHz的频谱分辨率的频谱仪设计，从而进一步提高了相关数据处理效率和满足了频谱仪设备发展与类似密集型数据处理的需求。

关键词：ROACH Simulink 多相 频谱分辨率

中图分类号：P161 文献标识码：A 文章编号：1674-098X(2015)01(b)-0018-02

频谱仪通过测量特定设备在特定频率范围内的响应，是进行频谱转换与测量的专用设备，用来对各种信号进行处理分析。近年来，频谱仪不但被更加广泛的使用在科研工作中，而且伴随着半导体技术的发展，宽带高速ADC、高性能数字信号处理芯片/算法，以及高性能计算机总线的出现，基于FFT原理的数字频谱技术近年来迅速发展起来。

目前国内外比较常见的频谱仪系统是基于Xilinx平台系列的FPGA的数字信号处理系统[1]。例如美国加州大学伯克利分校近年来也在研制类似的频谱仪，其ROACH构架平台包括具有2Gsps的高速采样率，进行对频带调制，多相滤波，傅里叶变换，和对数字信号的高速传输(通过CX4网络端口将数据以UDP包形式输出)，输出端可接驳GPU系统，交换机和高速网卡等，这一系列频谱仪的主要特点是具有多个Z-DOK接口，可以接驳多种I/O板(包括双路2GSa/s的采样卡或者四路500MSa/s的模数，数模转换卡)[2]。与此同时德国马普研究所近来也在研发一种扩展带宽快速傅里叶变换频谱仪系统[3]，该系统采用赛灵思最新的Xilinx Virtex-6可编程逻辑门电路(FPGA)处理核心，配备实时带宽单通道4 GHz、双通道2.5 GHz的数字转换和采样模块，实现了32 K傅里叶变换通的高分辨率数字频谱分析能力。

综合以上特点，结合目前Xilinx最新发展的Virtex系列FPGA技术发展，我们在ROACH平台中进行多相宽带数字频谱仪系统仿设计。

1 频谱仪设计框架

频谱仪的主要原理是将接收输入的基带混合中频(IF)信号进行外差式采样，将输入的信号进行量化与编码的数字测量后，通过数据总线传输至FPGA芯片核心，进行多相滤波与快速傅里叶变换运算的数字信号处理，然后计算信号的功率谱数据，得到信号测量结果，并按照X轴为频率，Y轴为功率谱的数据矩阵，将测量结果存储在频谱仪内置的存储芯片中。通过频谱仪的网络接口，可将数据输出至后端接收电脑，实现进一步的测量处理与图形绘制等工作。频谱仪机柜将多套频谱仪整合封装，并设置统一的电源与时钟模块。

2 多通道频谱仪的具体设计

不同与传统的Verilog语言底层设计，Xilinx近年来引入系统级模块化设计的概念，在这种方式下，不仅可以更直观、快速和灵活的进行系统级模块的构建，而且可以使设计者把更多精力集中在核心算法的实现，而不必拘泥于具体实现方式。

在Matlab Simulink平台下，我们调用Xilinx System Generator函数库和Xilinx函数库，依次连接ADC，PFB和FFT模块，ADC模块的时钟频率为2048MHz(1024Mhz的奈奎斯特采样率)，ADC四路同步输出，因此全局时钟为1/4的ADC频率，512 MHz；PFB和FFT模块中设置大小为2^14=16K通道输出，之后连接实部虚部的分离模块和系数量化模块，将FFT输出的18bit按比特位选择其中的8bit，并将输出合并为128bit的总数据链路。设计如图1所示。

传统的FFT处理过程中会产生一定的频谱泄露，为了实现更好的相邻频带频响特性，我们在设计中加入了多基于窗函数的相滤波器组，使得经过FFT处理过的每通道频谱得到最大信噪比。

如图2所示，整体设计调用Xilinx自带FIR滤波器模块进行拼接，FIR采用4抽头汉宁窗算法。通过多相FIR滤波器对输入的信号进行预谱带处理后，传输到后端的1024通道FFT模块中。

如果直接将128bit总数据链路输入到传输模块，则无法控制每个IP数据包中对应的频谱信息，因此设计中Bram寄存器对数据流进行重组。此模块的功能为顺序输入数据，按预存次序输出。将16K通道对应分配到16IP中，每个IP接收对应的1024个通道，具体的数据流信息如表1所示。通过数据流调整模块和IP数量的对应控制，可以实现对信号频谱的任意目标分配，非常灵活。

3 性能测试和分析

整体设计中，由于采用了FIR4抽头汉宁窗算法的多相滤波器组与FFT相结合的频谱仪处理核心，使得其多通道谱输出旁带得到很好的抑制。比较而言，汉宁窗的主瓣与海明窗相当，比同阶的矩形窗宽，旁瓣最大泄漏为-42 dB，矩形窗为-23 db，汉宁为-52db。汉宁窗旁瓣随其阶数的增加以- 5.5dB/个的速度衰减，为最快，海明窗与矩形窗则都以-1.5dB/个的速度下降。海明窗的旁瓣衰减略比汉宁窗大，但是随旁瓣增加其衰减很慢。因此实际信号分析中常常选用汉宁窗，只要选择恰当的阶数，使用汉宁窗截断的时候，频谱泄漏可以更小，这样就可以有效改善FFT输出的频率响应，抑制频谱泄露。

我们在Matlab中对实际设计进行了测试，如图3所示，在加入多相滤波器组后的FFT输出旁瓣抑制达到近-70db较之单独FFT输出的-30db旁瓣，旁带干扰大大降低。

对于频谱仪整体响应，在Matlab Simulink平台下，给系统输入叠加脉冲信号的白噪声输入进行系统测试。频谱仪输出测试结果如图4所示，从图中可以看出在对应的通道中得到了叠加信号频谱响应，验证频谱仪设计仿真工作正常。

4 结论

通过对数字频谱仪系统的分析和测试，达到了我们在实际应用中需要的设计要求；对比验证了加窗后的FFT处理有效地降低了通道间的频谱泄露，初步讨论和测试了数字频谱仪对后端系统的影响，稳定高速的实时宽带数字频谱仪系统可以提高整个数据接收系统的信号传输性能；最后对系统进行了上电测试，验证了设计方案的可行性和可靠性。要在射电望远镜系统中实际应用，还需要在现有初步研究的基础上进行更进一步的系统分析和测试，尤其是特殊环境的可靠性测试。

参考文献

[1]Aaron Parsons， Donald Backer，Andrew Siemion，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，Proc. 10th Asilomar Conference on Signals， Systems and Computers， Paci_c Grove， CA，2006.

[3]http：//www.radiometer-physics.de/rpg/html/Products_Spectrometers.html.endprint