一種雷達信號處理機的軟件設計*

馬曉靜

(解放軍電子工程學院，合肥 230037)

0 引言

信號處理機是雷達系統的關鍵部件之一，其性能直接影響到雷達作用距離、目標檢測性能和目標的識別。現代雷達信號處理廣泛采用數字脈沖壓縮、動目標顯示/動目標檢測(MTI/MTD）、恒虛警處理等先進的信號處理技術，大大提高了雷達的作戰性能。脈沖壓縮雷達能同時提高雷達的作用距離和距離分辨率，較好地解決了雷達作用距離和距離分辨率之間的矛盾。MTI/MTD 通過對回波相參脈沖串作MTI 濾波和匹配濾波處理來實現雜波抑制，進而提高系統在復雜背景中發現目標的能力。恒虛警處理(CFAR）主要實現當干擾(雜波或噪聲）電平增加時，能夠隨干擾電平的變化自動調整檢測門限值，將虛警率控制在正常的水平。

本文結合了某雷達信號處理機的工程實現。該處理機由線性調頻脈沖壓縮、動目標顯示、動目標檢測和恒虛警處理等工作模塊構成，共處理4800個距離單元，先做8192 點脈沖壓縮，然后取4800個有效單元進行后續處理。本文通過分析該處理機的工作原理，提出了軟件設計的方法。

1 信號處理過程及其工作原理

信號處理機主要實現對前端信號的接收和處理，完成對目標距離以及方位信息的提取，其工作原理如圖1所示。該信號處理機主要包括脈沖壓縮(PC）、數字對消(MTI）、相參積累(MTD）和恒虛警處理(CFAR）4個功能模塊。

圖1 信號處理系統框圖

脈沖壓縮采用頻域脈壓的方法實現。圖2 是頻域脈壓的原理性框圖。在具體工程實現時一般將匹配濾波器權系數存于隨機存取存儲器(RAM）中，當需要加權以降低距離旁瓣時，存于RAM中的權系數應當是匹配濾波器頻率響應與加權函數的乘積。

圖2 頻域數字脈壓原理框圖

MTI 采用兩脈沖對消或三脈沖對消。一次相消器是很簡單的濾波器，其實現無需相乘，每一個輸出采樣只有一個減法器，但其零頻凹口窄，對于雜波抑制效果不理想。二次相消器也是一種常用的MTI 濾波器，可通過級連兩個一次相消器得到。二脈沖對消原理圖見圖3，三脈沖對消原理圖見圖4。MTD 是在MTI后串接一組相鄰且部分重疊的窄帶濾波器組，覆蓋整個重復頻率的范圍，對不同通道進行相參積累處理。

圖3 兩脈沖對消原理圖

圖4 三脈沖對消原理圖

采用單元平均CFAR 實現恒虛警檢測。對于某一個檢測單元，取其兩側一定數目的距離單元作為參考單元，對所有的參考單元相加求和，再乘以適當的門限系數，就可以得到相應于給定的虛警概率的門限電平。此門限電平是隨著目標附近的距離段內的背景雜波或噪聲而連續變化的。大目標的強回波會波及到臨近的距離單元，使平均背景的測量受到影響。因此，檢測單元左右兩側的一些距離單元的輸出應不計入，以免對平均背景值產生影響。這些距離單元稱為保護單元，其原理如圖5所示。

圖5 單元平均CFAR 原理框圖

2 信號處理機的硬件環境

信號處理機采用ADI 公司的新型DSP 芯片ADSP-TS201S為主要器件設計處理板，通過對DSP 軟件編程完成信號處理算法的工程實現。信號處理板在一個標準的6U CPCI 板上集成了ADI 公司的8個600 MHz TS201S處理器，可提供28.8GFLOPS 峰值計算能力，其系統框圖如圖6所示。

圖6 信號處理機系統框圖

每個處理器的64 位總線無縫連接4 片32MX16位SDRAM(MT48LC32M16），4 片SDRAM 配置成32M×64 位，共256MB。整板SDRAM 容量達256MB×8=2048MB。CPCI 主機和DSP 可同時訪問SDRAM。

3 信號處理機的軟件設計

通過對工作模塊中的一片或幾片DSP 編程來完成模塊功能，多個模塊組成的并行處理系統保證了高速實時的數字信號處理能力。在能夠實現預期功能的前提下，軟件編程最重要的考慮是軟件處理的通過率。因此，在設計信號處理軟件時，要對每個模塊的運算通過率作詳細的分析。

3.1 脈沖壓縮的編程原理

脈沖壓縮模塊處理的是一個重復周期的回波數據。以8192 點脈壓為例，開辟一個長度為8192*2的緩沖區buffer1，用于存放一個周期的回波數據，即脈壓的輸入數據;再開辟一個長度為8192*2的緩沖區buffer2，用于存放回波的FFT的結果;另外開辟一個長度為8192*2的緩沖區buffer3，用于存放一個周期匹配函數的FFT的結果。buffer2和buffer3中順序存放的是8192個距離單元的實部和虛部。將同一個距離單元在buffer2和buffer3中分別對應的2個復數進行復數乘法，(a+bj）*(c+dj）=(ac－bd）+(ad+bc）j，在經過4 次乘法運算和2 次加法運算后得到復乘后的復數結果，按照實虛部依次存放在長度為8192*2的緩沖區buffer4中。再對復乘結果作IFFT，將得到的結果放在長度為8192*2的緩沖區buffer5中，其處理流程見圖7。其中，IFFT 是通過取共軛，作FFT，再取共軛，最后乘以1/N來實現的。

圖7 脈沖壓縮的軟件編程流程

3.2 MTI的編程原理

3.2.1 二脈沖對消

由其時域方程y(n）=x(n）－x(n－1）可知，一次相消器(二脈沖對消）就是由相鄰兩個脈沖作減法運算，n個脈沖對消后變為n－1個脈沖，即同一個距離單元的每相鄰2個采樣值就要作一次對消。這樣，5個脈沖經過二脈沖對消后變為4個脈沖，見圖8。

圖8 二脈沖對消的基本原理圖

在DSP 軟件設計中，它的本質就是一次復數減法，也就是由實部和虛部分別作一次實數減法運算。由于采用雙核運算提高運算速度，也就是在X和Y兩個計算塊中進行相同的減法運算，所以只需要一條指令就可以完成。

3.2.2 三脈沖對消

由其時域方程y(n）=x(n）－2x(n－1）+x(n－2）可知，二次相消器(三脈沖對消）就是由相鄰3個脈沖按照一定的規則進行加法運算，n個脈沖對消后變為n－2個脈沖，即同一個距離單元的每相鄰3個采樣值要作二次對消。這樣，10個脈沖經過三脈沖對消后變為8個脈沖，見圖9。

在DSP 軟件設計中，它的本質就是一次復數乘法和兩次復數加法(減法），也就是由實部和虛部分別作一次實數乘法運算和兩次實數加法(減法）運算。由于在X和Y兩個計算塊中進行相同的運算，所以只需要兩條指令就可以完成。

圖9 三脈沖對消的基本原理圖

以4800個距離單元、10個重復周期作三脈沖對消為例，脈壓后的數據按照重復周期存放，MTI后的數據按照距離單元存放。由于作三脈沖對消的3組復數數據在內存區中沒有連續存放，所以需要3個指針j1、j2、j3分別指向這3組數據，通過指針后加4800*2(4800個復數所占的內存單元數）更新到同一個距離單元對應的下一個脈沖(重復周期）的位置上，為下一次對消作準備。將同一個距離單元的8 次三脈沖對消的結果順序存放在一個長度為8*2的緩沖區內。每完成一個距離單元的三脈沖對消后，將3個指針j1、j2、j3 更新到下一個距離單元的位置上，再進行新一輪的三脈沖對消。這樣得到的下一個距離單元的三脈沖對消結果順序存放在MTI后的數據緩沖區內，與上一個距離單元的三脈沖對消的結果連續存放。見圖10。

圖10 MTI 前后的數據存放示意圖

3.3 MTD的編程原理

8 點MTD的基本原理就是對同一個距離單元的8個采樣值作8 點的DFT 運算。由于MTI后的數據已經是按照距離單元存放的，也就是說，同一個距離單元的8個采樣值是放在連續的16個內存單元內。所以，只需要用一個指針j14 指向長度為8*2*4800的輸入緩沖區的首地址，也就是第一個8*2的緩沖區的首地址上，通過調用DFT 這一工作模塊，實現第一個距離單元的8個采樣值的8 點DFT 運算。然后通過指針j14的后加16，將指針指向下一個距離單元的8個采樣值的長度為8*2的緩沖區的首地址，然后再調用DFT 程序，依此通過4800個循環，完成4800個距離單元的8 點MTD的運算。其基本原理見圖11。

圖11 8 點MTD的基本原理

這里利用雙核運算，考慮到要進行四字操作，所以在X和Y 計算塊中分別進行兩個距離單元的平方和運算。這樣一次內循環同時處理4個距離單元，所以8 點MTD循環數=距離單元數*8*2/8=距離單元數*2，4 點MTD 循環數=距離單元數*4*2/8=距離單元數。

3.4 CFAR的編程原理

設當前的參數設置:作兩邊檢測時，只在單側取4個保護單元，8個參考單元;作中間檢測時，在左右兩邊分別取2個保護單元，8個參考單元。

4800個距離單元，左右兩邊同時進行運算，分別放在X和Y 計算塊中，用j 指針和k 指針進行指向，這樣使循環數減半。

在進行中間單元的檢測時，仍然用j2和k2分別指向左邊和右邊的當前檢測單元。j1 是中間檢測時依次指向左邊的檢測單元的左邊的參考單元的指針，其初始位置是其最左邊的參考單元。k1 是中間檢測時依次指向右邊的檢測單元的右邊的參考單元的指針，其初始位置是其最右邊的參考單元。而與之相對應的是，j3 是指向中間檢測時依次指向左邊的檢測單元的右邊的參考單元的指針，其初始位置是其最左邊的參考單元。k3是指向中間檢測時依次指向右邊的檢測單元的左邊的參考單元的指針，其初始位置是其最右邊的參考單元。j4 是指向中間檢測時左邊的檢測單元的左邊的參考單元累加和依次要減去的參考單元，初始位置是其最左邊的參考單元。k4 是指向中間檢測時右邊的檢測單元的右邊的參考單元累加和依次要減去的參考單元，初始位置是最右邊的參考單元。j5 是指向中間檢測時左邊的檢測單元的右邊的參考單元累加和依次要減去的參考單元，初始位置是其最左邊的參考單元。k5 是指向中間檢測時右邊的檢測單元的左邊的參考單元累加和依次要減去的參考單元，初始位置是最右邊的參考單元。其處理流程圖如圖12所示，其參數設置為:M個重復周期，N個距離單元，最左邊N1個距離單元只取其右邊的參考單元求和，最右邊N2個距離單元只取其左邊的參考單元求和。

圖12 單元平均恒虛警的處理流程

4 結束語

本文是結合某雷達信號處理機的設計而展開的。該信號處理機包含脈沖壓縮(PC）、動目標顯示(MTI）、動目標檢測(MTD）、恒虛警檢測(CFAR）等功能模塊，在通用DSP 硬件平臺上通過軟件的方法實現了信號處理功能，具有功能強大、靈活性強的特點，取得了很好的工程應用效果，對軟件化雷達信號處理的實現具有一定的參考價值。

［1］向敬成，張明友.雷達系統［M］.北京:電子工業出版社，2001.

［2］A.V.奧本海姆.數字信號處理［M］.北京:科學出版社，1981.

［3］蘇洪濤.雷達脈沖壓縮系統設計及實現［D］.西安電子科技大學碩士學位論文，1999.

［4］趙樹杰.信號檢測與估計理論［M］.西安電子科技大學出版社，1998.

［5］李世杰.數字信號處理器的匯編程序優化方案［J］.國外電子元器件，2006(4）:48-51.

［6］李偉，牛曉麗，韓松.ADSP_TS101S 匯編語言編程優化［J］.電子測量技術，2006，29(3）:94-95.

［7］李文龍，劉利，湯志忠.軟件流水中的循環展開優化［J］.北京航空航天大學學報，2004，30(11）:1111-1115.