陳林杰，顏毅華，劉 飛，王 威

(1.中國(guó)科學(xué)院國(guó)家天文臺(tái)，北京 100012；2.中國(guó)科學(xué)院研究生院，北京 100049)

射電頻譜儀是分析天體射電頻譜信息的重要觀測(cè)設(shè)備。隨著科技的發(fā)展，為了滿足新的觀測(cè)要求和發(fā)現(xiàn)更多新的射電頻譜現(xiàn)象，建造同時(shí)具有寬頻帶、高時(shí)間分辨率、高頻率分辨率的射電頻譜儀顯得越來(lái)越重要。

傳統(tǒng)的射電頻譜儀采用混頻、濾波，得到基帶信號(hào)后進(jìn)行功率檢波，從而得到信號(hào)的頻譜強(qiáng)度。這種傳統(tǒng)頻譜儀的頻譜通道數(shù)取決于硬件的資源量，因此帶寬不可能很寬，通道數(shù)也不可能做到很多。隨著FPGA(現(xiàn)場(chǎng)可編程門陣列)的運(yùn)行速度越來(lái)越快，在FPGA內(nèi)運(yùn)用高速FFT變換對(duì)射電信號(hào)進(jìn)行頻譜變換，為寬帶射電頻譜儀的設(shè)計(jì)提供了另一種選擇。但是受限于當(dāng)前的FPGA內(nèi)部處理時(shí)鐘最高只能到500MHz，為了實(shí)現(xiàn)帶寬更寬、速率更高的射電頻譜分析，文中提出了一種基于多相濾波器的新型寬帶射電頻譜儀的設(shè)計(jì)方案。通過(guò)多相濾波器組對(duì)寬帶射電信號(hào)進(jìn)行濾波，濾波后的各頻率通道信號(hào)通過(guò)數(shù)字功率檢波得到信號(hào)的功率，再通過(guò)積分控制，靈活地進(jìn)行信號(hào)的積分，最終得到射電信號(hào)的頻譜強(qiáng)度信息。

1 實(shí)信號(hào)的多信道模型

為建立實(shí)信號(hào)多信道接收機(jī)的數(shù)學(xué)模型，首先，對(duì)實(shí)信號(hào)的數(shù)字頻譜做如下信道劃分[1]：

(1)

式中ωk為第k信道的歸一化中心角頻率；D為數(shù)據(jù)抽取率。圖1為實(shí)信號(hào)4個(gè)信道的頻譜分配圖。

圖1 實(shí)信號(hào)的信道劃分示意圖Fig.1 Illustration of channel assignment of real signal

如圖1所示，由(1)式得到的實(shí)信號(hào)信道存在著對(duì)應(yīng)的鏡頻ωk’，并且信道總數(shù)受數(shù)據(jù)抽取率D的限制。

根據(jù)低通抽取原理，把寬帶信號(hào)與復(fù)本振ejωkn,k=0,1,2,…,D-1，相乘可實(shí)現(xiàn)將第k個(gè)信道的數(shù)字譜移到基帶的目的，實(shí)現(xiàn)結(jié)構(gòu)[2-3]如圖2。

圖2 基于低通濾波器組的實(shí)信號(hào)信道化接收機(jī)Fig.2 Channelized receiver based on the low-pass filters

圖中，每個(gè)低通濾波器hLP(n)的帶寬為π/D，對(duì)應(yīng)的原型理想低通濾波器的頻率響應(yīng)為HLP(ejω)[2]：

(2)

由于經(jīng)過(guò)復(fù)本振的信號(hào)已為復(fù)信號(hào)，故可以對(duì)低通濾波后的信號(hào)進(jìn)行2D倍抽取。

2 基于多相濾波器的信道化接收機(jī)

基于低通濾波器組的信道化接收機(jī)結(jié)構(gòu)，抽取器位于濾波器之后，故當(dāng)抽取率D很大時(shí)，低通濾波器所需的階數(shù)可能會(huì)變得非常大，而且每一信道都要分配這樣一個(gè)濾波器，實(shí)現(xiàn)的效率非常低。為了獲得頻率分辨率一致的濾波器，這里引入多相濾波器結(jié)構(gòu)，有效地實(shí)現(xiàn)了多通道信號(hào)的濾波，并提取帶寬為基帶帶寬1/D的信號(hào)。

基于多相濾波器的實(shí)信號(hào)信道化接收機(jī)數(shù)學(xué)模型[4]如圖3。

圖3 實(shí)信號(hào)多相濾波信道化接收機(jī)模型Fig.3 Model of a channelized receiver with polyphase filtering

其中，多相濾波器的低通濾波器原型采用Parks-McClellan窗。這種窗在一致意義上對(duì)低通濾波器作最佳逼近，可以提供理想的頻率響應(yīng)[5]。確定的多相濾波器系數(shù)如下：

hk(m)=h(mD+k)m=0,1,2……Q-1;k=0,1,2……D-1.

(3)

3 基于多相結(jié)構(gòu)的寬帶射電頻譜儀

3.1 數(shù)字功率檢波

數(shù)字功率檢波的過(guò)程就是數(shù)字信號(hào)平方后經(jīng)過(guò)積分器得到信號(hào)的功率信息，如圖4。原理比較簡(jiǎn)單[6],但卻可以很好地記錄信號(hào)的功率強(qiáng)度信息。

圖4 數(shù)字功率檢波Fig.4 Principle of the digital power detector

3.2 基于多相結(jié)構(gòu)的射電頻譜儀

射電頻譜儀的目的是分析射電信號(hào)頻譜強(qiáng)度信息，為了得到信號(hào)的頻譜，通過(guò)多相濾波器可以高效地實(shí)現(xiàn)多通道信號(hào)的并行濾波，同時(shí)又可以大量減少計(jì)算量，因此可以考慮在信號(hào)經(jīng)過(guò)多相濾波器后，對(duì)每一通道的信號(hào)進(jìn)行數(shù)字功率檢波，通過(guò)控制積分序列的長(zhǎng)度來(lái)控制信號(hào)的積分時(shí)間，從而實(shí)現(xiàn)一定時(shí)間分辨率的射電頻譜分析。

本文提出基于多相結(jié)構(gòu)的新型寬帶射電頻譜儀的設(shè)計(jì)，如圖5。

圖5 多相濾波結(jié)構(gòu)的寬帶射電頻譜儀Fig.5 Wideband spectrum analyzer(for radio waves)with polyphase filtering

寬帶信號(hào)經(jīng)過(guò)A/D采樣，完成對(duì)輸入IF模擬信號(hào)的量化處理，再通過(guò)多相濾波器實(shí)現(xiàn)數(shù)字信號(hào)的基帶轉(zhuǎn)換，然后對(duì)各通道基帶信號(hào)進(jìn)行上述數(shù)字功率檢波，由于是數(shù)字信號(hào)，因此可以靈活地控制對(duì)其積分的序列長(zhǎng)度，從而調(diào)整了積分時(shí)間，實(shí)現(xiàn)了頻譜儀靈活的時(shí)間分辨率。

3.3 新型射電頻譜儀的實(shí)現(xiàn)

假設(shè)頻譜儀模擬接收機(jī)的中頻輸出帶寬為512MHz，頻譜分辨率為1MHz，時(shí)間分辨率為1ms。此時(shí)，采用1024MHz Sps進(jìn)行采樣，抽取率D=512，多相濾波器的原型采用Parks-McClellan窗的6144階低通濾波器，那么信號(hào)將被分為512個(gè)時(shí)鐘僅為1MHz的基帶信號(hào)，每個(gè)通道的采樣率為2MHz。FFT輸出的每通道信號(hào)為窄帶復(fù)信號(hào)，對(duì)其求模取平方，然后通過(guò)計(jì)數(shù)器控制積分器對(duì)每長(zhǎng)度為2000的序列進(jìn)行積分，即可得到1ms的時(shí)間分辨率下的功率信息。這里可以調(diào)整積分序列的長(zhǎng)度以實(shí)現(xiàn)幾乎任意長(zhǎng)度的積分，從而得到靈活的時(shí)間分辨率。由于采用硬件積分，可以極大地減少在后續(xù)軟件處理中對(duì)信號(hào)重新積分的運(yùn)算量。積分后的輸出即為基帶信號(hào)的功率信息，至此完成對(duì)射電信號(hào)的頻譜分析。

在頻譜儀進(jìn)行FPGA的實(shí)現(xiàn)中，對(duì)于多相濾波器，由于其階數(shù)比較高，即在FPGA中乘法器的消耗會(huì)比較多，即使現(xiàn)今最新的FPGA自身所帶的乘法器也無(wú)法滿足要求，因此這里考慮通過(guò)其他方法來(lái)實(shí)現(xiàn)乘法器。通過(guò)改進(jìn)的Booth算法產(chǎn)生部分積，用一種Wallace樹(shù)結(jié)構(gòu)壓縮部分積，并使用減少符號(hào)位填充和減少尾部0填充的方法來(lái)有效減小部分積壓縮器的面積，從而實(shí)現(xiàn)高效靈活的乘法運(yùn)算。另外，經(jīng)過(guò)抽取后的信號(hào)數(shù)據(jù)率很低，降低為2MB/s，多相濾波器的濾波器系數(shù)也是固定的，且系數(shù)關(guān)于中心對(duì)稱，因此可以通過(guò)查表法來(lái)實(shí)現(xiàn)乘法。通過(guò)對(duì)濾波器系數(shù)進(jìn)行量化，假定為8位，如上例，則僅需要256×256的16位查數(shù)表即可滿足多相濾波器對(duì)乘法器的要求。

4 仿真結(jié)果分析

為了驗(yàn)證以上的設(shè)計(jì)，文中對(duì)其進(jìn)行了仿真。輸入為50～450MHz的中頻白噪聲信號(hào)，采樣時(shí)鐘為1024MHz，8bit量化，時(shí)間序列長(zhǎng)度為12ms，在7ms時(shí)加入了4個(gè)信號(hào)，如下：

其中f1=63.5MHz；f2=64.5MHz；f3=255.5MHz；f4=384.5MHz。采樣率fs=1024MHz，L序列長(zhǎng)度為12ms。抽取率D為512，低通濾波器原型采用Parks-McClellan窗進(jìn)行設(shè)計(jì)，6144階。對(duì)以上輸入信號(hào)進(jìn)行了基于多相濾波器的頻譜儀輸出仿真,分別得到頻譜儀通帶內(nèi)的整個(gè)頻譜輸出(圖6)以及頻譜儀的時(shí)間和頻譜區(qū)線(圖7)。

圖6 多相濾波結(jié)構(gòu)的寬帶頻譜儀的輸出Fig.6 Output of the wideband spectrum analyzer with polyphase filtering

分析以上仿真曲線可以看出，輸出的信號(hào)功率在7ms時(shí)階躍增大，這和仿真輸入一致。第64、65、256和385通道的中心頻率分別為64.5MHz、65.5MHz、255.5MHz和384.5MHz，輸出中的4個(gè)信號(hào)也分別出現(xiàn)在這4個(gè)通道，證明多相濾波器數(shù)學(xué)模型的正確性。分析信號(hào)的頻譜功率可得表1。

由表1可以看出，信號(hào)的功率輸出和理論值保持著高度一致，誤差在1.5%以內(nèi)，這說(shuō)明信號(hào)的輸出與輸入成很好的線性關(guān)系，而這正是射電頻譜儀的一個(gè)最重要的性能指標(biāo)，從而也進(jìn)一步證明了該設(shè)計(jì)方案的準(zhǔn)確性和有效性。

圖7 頻譜儀的時(shí)間和頻譜曲線Fig.7 Power-versus-time and power-versus-frequency curves of the spectrum analyzer with polyphase filtering

信號(hào)功率之差理論值(dB)仿真值(dB)相對(duì)誤差P1-P23521834861101%P1-P31583615642122%P1-P46020659839061%P2-P42498824978004%P3-P21938219546085%P3-P44437044523034%

頻譜曲線中除了幾個(gè)點(diǎn)頻信號(hào)外，整個(gè)信號(hào)的頻帶為50～450MHz，這和輸入的白噪聲信號(hào)帶寬一致，仔細(xì)分析發(fā)現(xiàn)4個(gè)點(diǎn)頻信號(hào)也分別和輸入一致，輸入中f1和f2頻率之間僅相差1MHz，由此也可證明文中頻譜儀的頻譜分辨率達(dá)到了512MHz/D=1MHz。

5 結(jié) 論

文中提出了一種基于多相濾波結(jié)構(gòu)的寬帶射電頻譜儀設(shè)計(jì)，通過(guò)多相濾波器對(duì)寬帶中頻信號(hào)進(jìn)行濾波完成基帶轉(zhuǎn)換，然后通過(guò)數(shù)字功率檢波，再進(jìn)行可控的時(shí)間積分得到信號(hào)的基帶功率信息，從而完成射電信號(hào)的頻譜分析。仿真結(jié)果顯示，此設(shè)計(jì)得到了很好的效果，降低了信號(hào)的處理速率，極大地提高了硬件的利用效率，為寬帶射電頻譜儀的實(shí)現(xiàn)提供了一種切實(shí)有效的方案。

[1] 周欣,吳瑛.基于多相濾波的寬帶接收機(jī)信道化算法研究[J].現(xiàn)代雷達(dá)，2006，28(11)：71-74.

ZHOU Xin,WU Ying.A Study on Broad-band Channelized Receiver Algorithm Based on Polyphase Filter[J].Modern Rada,2006,28(11):71-74.

[2] JOHN D BUNTON.SKA CORRELATOR ADVANCES[J].Experimental Astronomy，2004，17:251-259.

[3] John Bunton.ALMA Memo 342-An Improved FX Correlator[C].CSIRO,Telecommunications and Industrial Physics,Australia，2000.

[4] 付永慶,李裕.基于多相濾波器的信道化接收機(jī)及其應(yīng)用研究[J].信號(hào)處理,2004,20(5):517-520.

FU Yong-qing,LI Yu.The realization of channelized receivers based on polyphase filters[J].Signal Processing,2004,20(5):417-520.

[5] 陳嵐，張秀忠.用于VLBI數(shù)字基帶轉(zhuǎn)換的多相濾波技術(shù)研究[J].天文學(xué)進(jìn)展,2008,26(1):87-94.

CHEN Lan,ZHANG Xiu-zhong.The Study of DBBC Based on Poly-phase Filter Banks and FFT in VLBI[J].Progress in Astronomy,2008,26(1):87-94.

[6] 宋勇輝,張志明,唐莽.一種基于多相濾波器的信道化數(shù)字檢波方法[J].航天電子對(duì)抗,2008,24(5):59-60.

SONG Yong-hui,ZHANG Zhi-ming,TANG Mang.A method of channelized digital demodulation based on poly-phase filter[J].Aerospace Electronic Warfare,2008,24(5):59-60.