現代電磁信號環境越來越復雜密集，要求電子戰接收機必須具有很寬的處理帶寬、高靈敏度、大動態范圍、多信號并行處理和大量信息實時處理的能力。而數字信道化接收機不僅可以較好地滿足上述要求，還可實現監視信道內信號的全概率截獲。
數字信道化過程是寬帶數字接收機的核心，目前廣泛采用基于多相濾波的數字信道化結構。這種結構先用高速的模數轉換器（A／D）進行數據采樣，得到的高速數據流經抽取降低數據速率后進入多相濾波器組，該濾波器組是由一個原型濾波器調制到多個支路。現場可編程門陣列（FPGA）中豐富的乘法器、鎖存器及數字信號處理算法IP核等資源，可以非常靈活地實現寬帶數字信道化接收處理算法。本文采用基于多相濾波器的結構實現了一種高效高速的寬帶數字信道化接收機，并在Altera公司的EP3SE110F1152C4上綜合實現，輸出載頻、相位信息。
1 信道化接收機的基礎理論
1.1 信道劃分
為建立實信號多信道接收機的數學模型，首先，對實信號的數字譜作如下信道劃分：
式（1）中，ωk為第k信道的歸一化中心角頻率；K為劃分信道數。圖1給出對應k=8時，實信道的頻譜分配情況。需要指出的是由于實信號的頻譜是對稱的，所以只有4個獨立的信道。
圖1 實信號的信道劃分示意圖
采用上述方法進行信道劃分有些頻點無法識別，為確保整個覆蓋帶寬內無盲區，信道的劃分選擇相鄰信道50％交疊，即擴大每個信道的處理帶寬，如圖2所示。
圖2 50％交疊的信道劃分示意圖
1. 2 多相濾波器組結構
本文采用高速高效數字信道化結構完成信道化接收，其數學模型如圖3所示。
圖3 高速高效數字信道化接收機結構圖
圖3中K為系統劃分的信道數，M為每個信道的抽取倍數且K=FM。H（n）為原型低通濾波器的單位沖擊響應，K個帶通濾波器都是通過該原型濾波器調制生成的，即均勻濾波器組的多相濾波分量。
第k個信道輸出為：
引入多相概念可得：
將wk=2πk/K帶入可得：
即為圖3的結構，該結構的信道化過程是在1／M的信號輸入速率下進行的，可以降低整個過程的運算量，使系統的復雜度和數據速率大大降低，實時處理能力得到提高。
要實現480～960 MHz的16信道劃分，所以選取K=16。選擇無盲區相鄰信道50％交疊的信道化分形式，F應該為2。根據上述原理M=8，信號需進行8倍抽取。
2 基于FPGA的信道化接收機實現
2.1 主要芯片介紹
ADC10D1000是NS最新推出的一款超高速低功耗10位模擬／數字轉換器，單通道最高采樣頻率可達到2.0 GHz，全功率帶寬為2.8 GHz。該芯片采用單電源1.9 V供電，總功耗只有2.8 W，比同級的A／D低33％，被NS列為Power Wise系列的高能源效率產品之一。該芯片采用292個球體的BGA封裝，令產品更小巧輕盈，而且散熱能力更強，即使沒有散熱器，系統也可在攝氏-40°～85°的工業級溫度范圍內工作。該芯片的無雜散信號動態范圍（SFDR）可高達66 dBc，達到業界最高水平，而且有效位數（ENOB）高達9.1位，為提高寬帶數字接收機的動態范圍提供了有力的條件。
ADC10D1000與8位高速A／D相比，在許多性能上有了提高，但輸入的最大模擬電壓的峰峰值為860 mV，相較于8位高速A／D較低，使得輸入信號的功率應在3 dB以下，建議使用時功率在2 dB以下。
選用的StratixIII系列EP3SE110F1152C4型號的FPGA。該系列的FPGA是世界上結合了最佳性能、最大密度和最低功耗的65-nm器件。具有最低的靜態和動態功耗，比上一代器件快了25％。Stratix III FPGA系列有33.8萬的邏輯單元（LE）和27萬的寄存器、擁有17.2Mb的600MHz內存和896個18x18的乘法器。Stratix III FPCA支持40多個I／O接口標準，支持高速內核以及高速I／O，已實現400 MHz DDR3，并且具有業界最佳的信號完整性。
2.2 系統硬件設計
系統的硬件設計框圖如圖4所示。
圖4 系統的硬件設計框圖
A／D轉換是進行數字化處理的前提，其性能直接影響接收機的整體性能。其性能指標主要有采樣速率和分辨率。射頻前端輸出信號的中心頻率為720 MHz，帶寬為480 MHz，根據帶通采樣定理，所需A／D器件的采樣速率應為960 MHz。要想得到大動態范圍的接收機，所需A／D器件的分辨率應越小越好，即輸出數據位數越多越好。綜合以上兩點選用了ADC10D1000。
為了給ADC10D1000提供更穩定相噪更好的時鐘信號，該時鐘信號由外部晶振和鎖相環（LMX2312和VCO190-964）產生。VCO190-964的頻率范圍為951-977 MHz，單端輸出。由FPGA控制LMX2312的工作方式及工作頻率，設計選用200 kHz為相位監測比較頻率，LMX2312通過比較自身時鐘信號與VCO反饋信號產生控制電壓，鎖定VCO的輸出頻率為960 MHz。
ADC10D1000輸入的時鐘信號要求為差分形式，因此要通過變壓器ADTL2-18對VCO輸出信號進行轉換，且變壓器輸出端應接100 Ω差分阻抗匹配A／D的輸入阻抗。A／D的輸出為LVDS信號，所以在與FPGA連接時要注意100 Ω匹配電阻要靠近FPGA管腳。為保證ADC10D1000的輸出不減
少數據吞吐率，設計采用內部1：2Demux增加數據寬度的方法，即同時并行輸出2組10位采樣數據，及DDR模式在時鐘上升沿和下降沿均輸出數據的方法降低了時鐘速率，使輸出時鐘頻率降為時鐘信號960 MHz的1／4，即240 MHz。
為了在調試時，可以很方便地修改FPGA內部判定信號的幅度閾值，不用等待FPGA長時間的編譯過程，在設計中加入DSP。DSP還可以校正相位差編碼，確立相位差的零點。
2.3 FPGA應用
軟件模塊主要包括數字信道化過程和后續數據輸出整理兩部分。數字信道化部分主要由上述數學模型構建，包括：數據抽取、符號轉換、多相濾波和IFFT。后續數據處理主要根據信道化結果得到頻率和相位信息及對同時到達三路信號的判斷。FPGA內部處理模塊框圖如圖5所示。
圖5 FPGA內部處理模塊框圖
2.3.1 數據抽取變換
A／D輸出的240 MHz高速采樣差分信號進入FPGA。根據上述模型，若信道數K=16，則抽取倍數M=8，FPGA通過LVDS接口的串并轉換實現8倍抽取。設計采用QuartusⅡ的LNDS模塊完成信號的串并轉換，降低信號及時鐘速率。設計中將模塊設置成INDS receiver形式，選擇8倍轉換因子，可得輸出信號16組以及和A／D的隨路時鐘同步的60 MHz時鐘。用該時鐘作為全局時鐘驅動后續所有處理模塊。由于A／D采樣輸出信號為偏移二進制類型，需經符號變換模塊后變為二進制補碼類型。
2.3.2 多相濾波器的設計與仿真
多相濾波器組采用具有穩定系統、可以實現線性相位的FIR型數字濾波器。FIR原型濾波器的設計主要考慮采樣頻率fs、通帶波紋rp、阻帶衰減rs以及過渡帶寬。例如采用fs=960 MHz，rp=0.1 dB，rs=63 dB，過渡帶起始頻率15 MHz，截止頻率30 MHz，得到原型濾波器幅頻特性曲線如圖6所示。該低通原型濾波器的階數為192階，將原型濾波器分為32相，每相濾波器為6階。由于采用50％交疊的結構需間隔插零，每相濾波器階數增至12階。設計使用程序編寫乘法累加運算實現FIR濾波。FIR原型濾波器的系數通過MATLAB生成導出，量化后寫入FPGA的濾波器程序中。
圖6 原型濾波器幅頻特性曲線
2.3.3 IFFT運算
IFFT運算采用按時間抽選的基-2算法。為了加快信號處理的速度，IFFT模塊采用多級流水線設計，并且運算模塊利用Quartus的宏產生。例如IFFT運算的核心蝶形運算可由Altmult_complex宏和lpm_add_sub宏實現。每次復數乘法會占用4個18x18 DSP乘法器資源，所以單路信道化的IFFT共需占用136個乘法器資源。
2.3.4 信道輸出
因為輸入的是實信號，經IFFT得到16個信道的子帶信號。對每個信道采用旋轉數字計算機算法（CORDIC）計算每個信道信號的幅度及瞬時相位。根據CORDIC輸出的信號幅度判斷信號是否存在以及信號的起始點和結束點，給出對應的包絡脈沖。同時利用CORDIC輸出相位根據瞬時相位差法計算頻率。為了提升測頻的準確度，用脈沖上升沿平穩后的連續4個無模糊的相位差平均值測頻，輸出載頻編碼。用兩通道信號的CORDIC輸出相位測算兩通道信號的相位差，輸出相位差編碼。
為了節省對外接口資源，最多只輸出三路信號即同時處理三路不同信號，當某路信道上出現包絡脈沖時才將該信道的頻率碼和相位差碼輸出，否則不輸出。16個信道都要進行判斷，確定是否輸出。具體流程如圖7所示，當判斷不成立或者語句執行結束時，結束程序。
圖7 輸出的判斷邏輯流程圖
3 系統硬件仿真與結果分析
本設計在EP3SE110F1152C4上完成了兩通道的信道化過程、信號包絡脈沖輸出及對載頻、相位差信息的編碼輸出。在硬件驗證仿真時，用到了內嵌式邏輯分析儀——SignalTapⅡLogic Analyzer。它是一種調試工具，能捕獲和顯示FPGA中的實時信號特性，通過JTAG接口下載FPGA配置數據和上載捕獲的信號數據，并在計算機中觀察FPGA內部節點信號，使用戶可以在整個設計工作過程中以系統級的速度觀察硬件和軟件的交互作用。FPGA芯片各項資源消耗情況如表1所示，共占用82％的資源，其中包括SignalTapⅡLogicAnalyzer所占用的資源。
表1 FPGA芯片各項資源消耗情況
A／D的采樣精度直接影響后面的精度，因此首先對A／D進行性能測試。信噪比RSN定義為信號峰值點的功率與去掉零頻以及前五階諧波分量后的所有噪聲的功率比值。信號噪聲失真比SINAD定義為信號峰值點的功率與去掉零頻后的所有諧波及噪聲的功率比值，其值較信噪比小。無雜散動態范圍SFDR定義為單信號輸入時信號與最大的諧波或雜散的功率比值。
實驗一：輸入信號頻率為由信號源Agilent 83752A產生的正弦波，頻率為720 MHz，幅度為-1 dBFS，采樣頻率為960MHz，從FPGA中導出采樣數據作8 k點的FFT，得信號頻譜如圖8所示。
圖8 A/D輸出720MHZ信號頻譜圖
經計算得，信噪比RSN為47.5 dB，信號噪聲失真比SINAD為46.3 dB，有效位數ENOB為7.4 bits，無雜散動態范圍SFDR為59 dBc。
實驗二：用Agilent的E4438C矢量信號發生器作為中頻輸入，輸入載頻為725 MHz，PRI=10μs，PW=2μs的脈沖信號測試結果如圖9所示。圖9中第一行表示輸入信號經過LVDS降速后的輸出波形，中間15行表示15個信道包絡脈沖輸出，倒數第二行表示有包絡脈沖輸出的那一路輸出載頻碼，最后一行表示有包絡脈沖輸出的那一路輸出的相位差碼。
圖9 輸出界面
由以上的分析可知，載頻為725 MHz信號應該出現在705～735 MHz的第9信道上，輸出載頻碼為725-480=245，DSP寫入校正編碼使輸出的相位差碼為0。由圖9可看出，只有第9信道有包絡脈沖輸出，輸出載頻碼為245，輸出的相位差碼為0，這與理論結果一致。
實驗三：用Agilent的E4438C矢量信號發生器作為中頻輸入，輸入載頻為725 MHz，PRI=10μs，PW=2μs的脈沖信號。用示波器同時采集輸入中頻脈沖信號和輸出的信號包絡脈沖，可得信號載頻碼和相位差碼輸出延遲時間，即整個系統延遲時間測試結果如圖10所示。上邊的一條線為輸入的中頻脈沖信號，下邊的一條線為輸出的信號包絡脈沖，由圖10可以看出系統延遲時間小于1.3μs，保證了系統的實時處理。
圖10 系統延遲時間
實驗四：用一臺Agilent的E4438C矢量信號發生器和兩臺Agilent的83752A作為中頻輸入，分別輸入載頻510MHz，PRI=100μs，PW=10μs；載頻為720MHz，PRI=90μs，PW=8μs；載頻為930 MHz，PRI=80μs，PW=20μs的三路脈沖信號。用示波器采集三路信號包絡脈沖輸出接口信號，可得系統對多信號處理結果如圖11所示。最上邊的線為第一路包絡脈沖輸出接口，中間的線為第二路包絡脈沖輸出接口，下邊的線為第三路包絡脈沖輸出接口。當信號在時域交疊時，由不同的輸出接口輸出包絡脈沖；否則在第一路輸出接口輸出。由圖11可以看出系統完成了對同時到達多信號的處理。
圖11 系統的多信號處理結果
4 結論
本文結合工程實際，完成了960MHz的16通道數字信道化接收機的FPGA實現。采用多相濾波器的高速高效數字信道化結構實現的數字信道化接收機，既能保證寬瞬時帶寬要求，又能達到實時處理的目的；與傳統的數字信道化結構節省硬件資源，提高系統的整體工作性能。FPGA仿真結果表明該模型在FPGA上實現的可行性以及實用性，并且實現了預期的指標要求。