(公安部第一研究所，北京 100048)

針對(duì)近年國際上出現(xiàn)的新型人體炸彈威脅事件及由此產(chǎn)生的緊張航空安全形式，毫米波成像設(shè)備越來越受到機(jī)場(chǎng)和旅客的關(guān)注，多家安檢研究機(jī)構(gòu)和設(shè)備制造商都在開發(fā)基于毫米波成像技術(shù)的人體安全檢查系統(tǒng)[1，2]。毫米波人體成像設(shè)備是指通過測(cè)量人體與物體，利用其在時(shí)域、頻域、空域、極化域和功率域物體的散射信息差異，檢測(cè)隱匿在人體衣物下的各種危險(xiǎn)物品。數(shù)字信號(hào)采集與處理技術(shù)是毫米波成像設(shè)備關(guān)鍵技術(shù)之一，設(shè)備采集到的信號(hào)經(jīng)過AD9643后將模擬信號(hào)轉(zhuǎn)變成FS高頻數(shù)字信號(hào)，數(shù)字采集系統(tǒng)再將采集的回波數(shù)據(jù)通過網(wǎng)線上傳給上位機(jī)進(jìn)行信號(hào)處理。直接使用高采樣率數(shù)字信號(hào)造成采集的數(shù)據(jù)量大，以太網(wǎng)傳輸帶寬不夠，導(dǎo)致數(shù)據(jù)傳輸不暢，需要通過抽取的方式，降低采樣數(shù)據(jù)量，且降頻抽取后的采樣頻率滿足奈奎斯特采樣定理，完整地保留中頻回波信號(hào)中的信息。濾除帶外系統(tǒng)噪聲信號(hào)，避免對(duì)目標(biāo)物形成干擾，因此需要對(duì)AD轉(zhuǎn)換后的數(shù)字信號(hào)在抽樣前進(jìn)行數(shù)字濾波操作。

對(duì)于常規(guī)采樣數(shù)據(jù)的大倍數(shù)抽取多采用逐級(jí)濾波逐級(jí)抽取的級(jí)聯(lián)形式[3，4]，最常見是使用多個(gè)FIR濾波器[5-7]串聯(lián)級(jí)聯(lián)的方式，但是這種方式只能滿足采集系統(tǒng)采集固定中頻信息，一旦輸出頻帶范圍改變，就必須重新設(shè)計(jì)濾波器參數(shù)。實(shí)際經(jīng)過多代射頻收發(fā)機(jī)設(shè)計(jì)，收發(fā)機(jī)最終輸出中頻信號(hào)頻帶范圍一直在隨著設(shè)計(jì)不同有所變化，為了滿足數(shù)字采集與處理系統(tǒng)能夠采集到不同頻段回波信息，將濾波器設(shè)計(jì)為并聯(lián)型多級(jí)級(jí)聯(lián)結(jié)構(gòu)，根據(jù)回波信號(hào)的頻帶選擇降頻抽取倍數(shù)，避免反復(fù)修改濾波器參數(shù)。并聯(lián)結(jié)構(gòu)的每一個(gè)分支滿足采集不同中頻信息，每一個(gè)分支且都采用了多級(jí)濾波抽取的串聯(lián)級(jí)聯(lián)結(jié)構(gòu)。相比傳統(tǒng)的濾波器設(shè)計(jì)，并聯(lián)型多級(jí)濾波器體現(xiàn)出對(duì)多頻段采樣的靈活性和便利性。在本領(lǐng)域研究中，目前還沒有對(duì)多頻段采樣的濾波器設(shè)計(jì)的相關(guān)報(bào)道，設(shè)計(jì)基于產(chǎn)品的實(shí)際應(yīng)用，提出了并聯(lián)型多級(jí)濾波器的新設(shè)計(jì)方法。

1 基于MATLAB并聯(lián)型多級(jí)濾波器設(shè)計(jì)

高頻FS=250MHz的AD直接采集回波數(shù)字信息的數(shù)據(jù)量大，不滿足以太網(wǎng)傳輸要求，需要通過抽取的方式，降低采樣數(shù)據(jù)量。同時(shí)經(jīng)過多代優(yōu)化設(shè)計(jì)的射頻收發(fā)機(jī)輸出中頻信號(hào)頻帶會(huì)隨著設(shè)計(jì)不同有所變化，為了滿足數(shù)字采集與處理系統(tǒng)能夠采集到不同頻段回波信息，低通濾波器設(shè)計(jì)為并聯(lián)型多級(jí)結(jié)構(gòu)見圖1，將AD采樣率FS經(jīng)低通濾波抽取后選擇性設(shè)計(jì)輸出FS/4、FS/8、FS/16，濾波器可實(shí)現(xiàn)4/8/16倍降頻濾波，從而通過降低采集的數(shù)據(jù)量達(dá)到減緩中頻數(shù)據(jù)傳輸速率的目的，且避免反復(fù)修改濾波器參數(shù)濾除無用高頻信號(hào)，得到特定頻率的有用信號(hào)。

1.1 并聯(lián)型多級(jí)濾波器設(shè)計(jì)

本研究設(shè)計(jì)的并聯(lián)型多級(jí)濾波器，基于最大化降低濾波器設(shè)計(jì)與計(jì)算復(fù)雜度，利用半帶濾波器的優(yōu)越性，將前級(jí)采用了FIR半帶濾波器[8]，最后一級(jí)采用一般FIR濾波器。如果最后一級(jí)也采用FIR半帶濾波器，濾波器系統(tǒng)的濾波效果不能滿足要求，濾波功能不能實(shí)現(xiàn)，并且出現(xiàn)噪聲信號(hào)。并聯(lián)型多級(jí)低通濾波器框架如圖1。

圖1 并聯(lián)型多級(jí)濾波器結(jié)構(gòu)

1.2 基于MATLAB濾波器仿真設(shè)計(jì)

使用MATLAB設(shè)計(jì)時(shí)[9-11]，AD采樣頻率為FS=250MHz，要分別實(shí)現(xiàn)3種抽取頻率，分別為FS/4、FS/8、FS/16，這里選取并聯(lián)型多級(jí)濾波器中抽樣頻率是FS/8即8抽取做為示例，既抽樣后的采樣頻率為31.25MHZ，這次射頻收發(fā)機(jī)輸出的中頻回波信號(hào)為6MHz左右的帶寬信號(hào)，濾波器帶寬為6MHz，帶寬為6 MHz，帶外頻率信號(hào)都會(huì)被濾除。MATLAB頻率設(shè)計(jì)時(shí)，先設(shè)定采樣頻率、有用信號(hào)頻率信號(hào)，再濾除噪聲頻率，能夠?qū)崿F(xiàn)8抽取且濾除噪聲的功能,前兩級(jí)FIR半帶濾波器只需要考慮通帶截止頻率即FBC2，最后一級(jí)FIR低通濾波器見圖2。

使用設(shè)計(jì)的多級(jí)抽取濾波器處理毫米波攜帶目標(biāo)物信息并夾雜系統(tǒng)噪聲的時(shí)域回波信號(hào)見圖3，A圖為抽取濾波前采集的回波信號(hào)，B為抽取濾波后回波信號(hào)，可看出初始信號(hào)經(jīng)過8倍抽取的濾波器后，采樣點(diǎn)由4096減少到4096/8=512，剛好滿足8抽取的結(jié)果，并且有效濾除帶外噪聲。

圖2 設(shè)計(jì)的低通濾波器

圖3 抽取濾波前后時(shí)域信號(hào)比較A .抽取濾波前時(shí)域信號(hào) ；B .抽取濾波后時(shí)域信號(hào)

圖3的頻域信號(hào)如圖4所示，A代表濾波前頻域信號(hào)，B代表濾波后頻域信號(hào)，此時(shí)橫軸代表信號(hào)頻率。濾波器的通帶范圍為3MHz到9MHz，圖4有4個(gè)頻率點(diǎn)分別為6MHz、12MHz、20MHz和50MHz，其中6MHz是回波信號(hào)頻率，經(jīng)過濾波器后，將12MHz、20MHz和50MHz 3個(gè)頻率濾除。濾波前有多個(gè)信號(hào)頻率，包括有用信號(hào)6MHz與其他信號(hào)的噪聲頻率，經(jīng)過濾波器后，此時(shí)只有1個(gè)波峰，即有用的信號(hào)頻率，即6MHz，與理論推導(dǎo)結(jié)果一致，其余信號(hào)都被濾除，滿足設(shè)計(jì)要求。

圖4 抽取濾波前后頻域信號(hào)比較A.抽取濾波前頻域信號(hào)；B.抽取濾波后頻域信號(hào)

3 并聯(lián)型多級(jí)濾波器在毫米波成像系統(tǒng)中實(shí)現(xiàn)

3.1并聯(lián)型多級(jí)濾波器在數(shù)字處理系統(tǒng)中實(shí)現(xiàn)

數(shù)字采集及處理系統(tǒng)是毫米波成像系統(tǒng)的重要部分，作用為采集射頻接收機(jī)輸出的中頻回波信號(hào)，信號(hào)經(jīng)AD轉(zhuǎn)化后分兩路進(jìn)入FPGA，等待下一步的信號(hào)處理。下一步信號(hào)處理主要是針對(duì)中頻信號(hào)的正交混頻、抽取和濾波，并聯(lián)型多級(jí)濾波器就是其中的抽取和濾波，它的目標(biāo)是根據(jù)需要選定相應(yīng)的抽取倍數(shù)和相對(duì)應(yīng)的濾波器，這種并聯(lián)型多級(jí)濾波器滿足多種收發(fā)機(jī)輸出不同頻段回波信息的信號(hào)處理要求，增加了信號(hào)處理的多樣性和靈活性。

正交混頻[12]是將攜帶目標(biāo)物信息的中頻信號(hào)變換為基帶信號(hào)，便于后面三維毫米波全息成像算法的計(jì)算，圖5中濾波器就是本研究設(shè)計(jì)的并聯(lián)型多級(jí)濾波器，可以通過降頻選擇的方式選擇抽取倍數(shù)。系統(tǒng)FPGA硬件處理系統(tǒng)框架為實(shí)現(xiàn)中頻信號(hào)的下變頻濾波處理。

圖5 FPGA下變頻濾波預(yù)處理系統(tǒng)框架

正交混頻如下：

Iout=Iin×cos(wt)+Qin×sin(wt)，

(1)

Qout=Iin×sin(wt)+Qin×cos(wt)，

(2)

其中w是通過UDP報(bào)文給出的中頻頻率。采樣間隔設(shè)置為固定采樣率FS；中頻頻率w在設(shè)計(jì)中要求按相參的FS采樣頻率折算。其中正弦函數(shù)按查表獲取，數(shù)據(jù)表精度為16位整型。輸出為17位精度，進(jìn)入濾波器處理.

FPGA設(shè)計(jì)遵循MATLAB設(shè)計(jì)思路，經(jīng)過兩級(jí)半帶濾波器實(shí)現(xiàn)降采樣，通過FIR濾波器實(shí)現(xiàn)最終濾波功能。FPGA芯片采用Xilinx廠家的XC7K325T，芯片資源較多，首先是板卡芯片的全局時(shí)鐘管理以及復(fù)位管理，然后為AD芯片使用FPGA模擬SPI協(xié)議實(shí)現(xiàn)AD芯片的寄存器配置。濾波器輸入以及三級(jí)濾波輸出信號(hào)如圖6方框所示，可看到經(jīng)過濾波器后的信號(hào)保留了輸入信號(hào)中的低頻信息，因此證明濾波器架構(gòu)無誤，完美濾除頻帶外高頻信息，保留低頻信息。

圖6 FPGA 測(cè)試抽取濾波前后時(shí)域信號(hào)

3.2 并聯(lián)型多級(jí)濾波器對(duì)毫米波全息圖像的影響

近程毫米波三維全息成像原理是射頻收發(fā)機(jī)通過天線陣列向目標(biāo)物發(fā)射具有一定帶寬的毫米波信號(hào)，經(jīng)反射后信號(hào)與本振信號(hào)進(jìn)行相干處理，后經(jīng)一系列近似處理及相位補(bǔ)償，最后反演出清晰的三維全息像[1，2]。本研究設(shè)計(jì)的并聯(lián)型多級(jí)濾波器主要應(yīng)用于毫米波三維探測(cè)系統(tǒng)中的數(shù)字采集與處理系統(tǒng)。當(dāng)攜帶目標(biāo)物信息的中頻信號(hào)進(jìn)入數(shù)字采集與處理系統(tǒng)后，經(jīng)過正交混頻的、抽取、濾波處理后將數(shù)據(jù)信息上傳給上位機(jī)，通過三維全息成像算法運(yùn)算，最后反演出探測(cè)物全息像。

毫米波圖像的噪聲來源于系統(tǒng)噪聲，所以有效的濾除這些干擾因素，濾波器起到至關(guān)重要的作用。圖7展示了距離向信息使用并聯(lián)型多級(jí)濾波器的去噪過程，A圖是濾波前距離向信息，B圖是濾波后距離向信息，可以從B中看到多級(jí)濾波器有效濾除了帶外噪聲。圖8是三維全息像在距離方向最大投影后的二維像，可以清晰地看到圖像中測(cè)試平臺(tái)上放置的刀具、剪刀和鉗子。A圖是沒有經(jīng)過并聯(lián)型多級(jí)濾波器處理的圖像，背景噪聲明顯，B圖是經(jīng)過并聯(lián)型多級(jí)濾波器處理的圖像。

為更清晰比對(duì)濾波前后對(duì)毫米波圖像的影響，計(jì)算了濾波器前后圖像的信噪比SNR(SIGNAL-NOISE RATIO)，SNR是衡量圖像性能一項(xiàng)重要指標(biāo)，是系統(tǒng)中信號(hào)與噪聲的比例，信噪比越大，基底噪聲越小。首先選取圖8中A、B兩圖同一目標(biāo)物區(qū)域面積，計(jì)算該區(qū)域內(nèi)像元均值，然后選中圖8中A、B兩圖同一位置均勻背景區(qū)域，計(jì)算區(qū)域內(nèi)標(biāo)準(zhǔn)差，圖像信噪比SNR為均值和標(biāo)準(zhǔn)差的比值[13]。經(jīng)計(jì)算可知，濾波前(圖8A)SNR為4.7204，濾波后(圖8B)SNR為12.6933(見表1)，濾波后信噪比提高3倍左右，基底噪聲變小，圖像背景更加干凈透亮，從而進(jìn)一步證明了并聯(lián)型多級(jí)濾波器設(shè)計(jì)的合理性，驗(yàn)證了其結(jié)構(gòu)的正確性。

圖7 濾波前后距離向信息比較 A.濾波前距離向信息；B.濾波后距離向信息

圖8 濾波前后毫米波圖像比較A.濾波前毫米波圖像；B.濾波后毫米波圖像

表1 濾波前后毫米波圖像信噪比比較

4 結(jié)論

針對(duì)毫米波成像設(shè)備數(shù)字采集與處理系統(tǒng)，研究提出一種基于多頻帶選擇的并聯(lián)型多級(jí)濾波器，該濾波器前兩級(jí)采用了FIR半帶濾波器，最后一級(jí)采用一般FIR濾波器，這種設(shè)計(jì)結(jié)構(gòu)最大化降低計(jì)算復(fù)雜度，有效消除帶外噪聲。通過分析具體頻帶信號(hào)實(shí)際特性，使用MATLAB設(shè)計(jì)確定濾波器參數(shù)，并將其加載在FPGA。通過測(cè)試驗(yàn)證，該濾波器可針對(duì)回波信息頻段特性進(jìn)行選頻采集，靈活便捷、高效實(shí)用；可對(duì)帶有雜波信息的回波信號(hào)進(jìn)行抽取濾波，回波數(shù)據(jù)的數(shù)據(jù)量呈倍數(shù)減少，降低數(shù)據(jù)速率，滿足以太網(wǎng)傳輸帶寬傳輸速率的要求，實(shí)現(xiàn)回波數(shù)據(jù)實(shí)時(shí)傳輸；可有效去除毫米波圖像背景噪聲，顯著提升毫米波圖像質(zhì)量。