基于TMS320C6678的ISAR實時成像方法研究

2015-04-14 03:05:50袁子喬

火控雷達技術 2015年2期

關鍵詞：方法

楊剛向聰袁子喬趙博

(西安電子工程研究所西安 710100)

0 引言

伴隨著寬帶雷達技術的發展，逆合成孔徑雷達ISAR(Inverse Synthetic Aperture Radar)的在國防和民用領域的應用不斷增加。ISAR 具有全天候、全天時的對飛機、艦船等遠程目標進行高分辨成像的優點，所以一直是雷達領域的研究熱點之一。經過眾多科研工作者這些年來的努力，ISAR 技術的逐漸發展成熟。但ISAR 成像算法計算量極大，在實時性要求高的場合，ISAR 的實現是一項繁重的任務。距離-多普勒(R-D)方法是一種研究比較成熟且應用廣泛的ISAR 成像方法，這種方法是基于轉臺模型的，該模型假設雷達與目標之間不存在平動分量，也就是進行了非常理想的運動補償［1］，這種方法相對簡單，便于工程實現，本文將采用R-D 方法進行成像。

TMS320C6678(后文中簡稱C6678)是TI公司于2010年推出的一款具有業界領先處理能力的多核DSP 處理器，其采用的是KeyStone 多內核架構，片內集成了8 個C66x 核，每個核都具有定點和浮點運算能力，同時也集成了SRIO(Serial RapidIO)，PCIe(Peripheral Component Interconnect Express)等接口，為方便數據的搬移，配置了EDMA(Enhanced Direct Memory Access)。TI公司也為開發人員提供了C6678 的MathLib，DSPlib 等函數庫。

本文在分析了ISAR 成像流程的基礎上，對成像流程進行了合理的劃分，在C6678 上進行了并行處理。最后通過仿真數據，驗證了并行處理的高效性和成像結果的正確性。

1 步進頻ISAR 成像原理

ISAR 成像的目的是從回波數據中求出目標散射點的幅度及其在二維平面上的分布。ISAR 成像的方式是，雷達不動，目標相對于雷達運動。步進頻ISAR 成像的第一步是對每個CPI的數據生成一維距離像。目標和雷達的相對運動可以分解為目標相對于雷達的平移運動(平動)、目標相對于雷達的旋轉運動(轉動)和目標自身的旋轉運動。平動分量會破壞回波信號之間的相干性，造成ISAR 成像結果的失真［2］，運動補償的目標就是補償掉目標相對于雷達的平動部分。

圖1 ISAR 成像流程

1.1 一維距離像原理

ISAR 成像的第一步是利用每個CPI內的回波數據得到一維距離像。線性調頻步進信號合成高分辨距離像的方法中，首先對子脈沖回波進行脈沖壓縮得到粗分辨距離像，然后對一個CPI內的粗分辨距離像進行IFFT 得到高分辨距離像。為降低采樣損失，一般參數的選擇使得IFFT 處理結果有部分重疊，按照一定的準則，從IFFT 結果中選取一些點，組成完備一維距離像的過程就是目標抽取算法［3］。

目標抽取算法已經非常成熟，應用較多的方法有同距離舍棄法、同距離選大法和同距離累加法，這三種方法是分別對IFFT 結果中重疊部分進行舍棄、選大和累加取平均操作實現的。另外，文獻［4］提到了移位抽取法一種降低抽取損失，該方法是通過將抽取窗與目標同步移動實現的。綜合考慮，本文中采用了移位抽取法來實現一維距離像。

文獻［5］指出當目標速度比較高時，在脈沖持續時間內目標和雷達之間距離不變的假設不再成立，要想得到高分辨一維距離像，需要進行高速運動補償，假設通過窄帶計算的目標速度為v1，第k 個脈沖對應的回波數據為s(k)，脈沖重復時間為Tr，則一維距離像的高速運動補償公式見式(1)。ISAR成像流程如圖1所示。

1.2 包絡對齊

包絡對齊的目的是消除由于目標相對于雷達的平動帶來的相鄰一維距離像的錯位問題。目前應用較多的方法是相鄰互相關法［6］，該方法假設目標平動造成的相鄰一維距離像之間的包絡變化較小，有較強的相關性，一般情況下，ISAR 成像可以滿足該條件。互相關的計算可以采用復包絡，也可以采用實包絡，本文中采用實包絡，公式為:

根據互相關最大值的位置可以確定相鄰距離像偏移的距離單元數(第i個距離像偏移的距離單元數記為ri)。為減少累加誤差，本文中是將每個距離像與前面所有距離像移位后的加權值進行互相關運算來求偏移的距離單元數。

由式(2)所得的偏移量為采樣間隔的整數倍，為提高包絡對齊的精度，一般是先對距離像進行插值操作，然后再做相鄰互相關。考慮到距離像的間隔時間相同，成像時間很短，假設在成像期間目標速度不變，那么，求得的偏移距離單元數與距離像編號應為線性關系，通過對所求距離像偏移距離單元進行最小二乘擬合，得到偏移距離單元r和距離像i之間的線性關系:

利用最后一個距離像的偏移單元rN和第一個距離像的偏移單元ri可以求得目標速度v，進而對每個距離像進行包絡對齊。

1.3 相位補償

相位補償中應用比較廣泛的方法是多普勒中心跟蹤法［6］，其實現方法是，將包絡對齊后的距離像中，相鄰距離像的各距離單元相位差進行加權平均，所得結果作為相鄰距離像的多普勒中心相位差，計算公式為:

2 實時成像方案

本系統中所采用的CPI數為256，脈沖數為16，距離單元數為128，目標抽取時在每個IFFT 結果中抽取4 個點。本系統采用了FPGA+DSP 的處理方式，二者通過RapidIO 接口傳輸數據。FPGA 對每個CPI內的回波數據做完距離向壓縮后，連同通過窄帶信號計算出的速度v1一起發給DSP，DSP 繼續后續流程，本文主要討論DSP 處理流程。并行處理模型有兩種，主從模式和數據流模式［7］，主從模式比較適合此系統的處理流程，本系統中選取核0 作為主核，核1-7 作為從核，八個核進行并行處理。

圖2 主核流程圖

2.1 數據存儲

C6678 中每個核都包含512KB 的二級存儲器(LL2)，八個核共享4MB 的共享存儲區(MSM)，同時DSP 還支持8GB 的DDR3 存儲器。本系統中采用的存儲原則是，數據量大的和八個核共享的數據存儲在MSM 中，并行運算時每個核獨立處理的、數據量小的和臨時的數據存儲在LL2 中。

在計算時為避免多核對MSM 的訪問沖突消耗時間，在計算前把每個核需要的數據通過EDMA 搬進LL2，在計算完成后將計算結果同樣通過EDMA搬至MSM。C6678 中擁有3 個CC(Channel Controller)，10 個TC(Transfer controller)，對應關系為CC0(TC0，TC1)，CC1(TC0，TC1，TC2，TC3)和CC2(TC0，TC1，TC2，TC3)。10 個TC 可以被配置起來同時執行數據傳輸，本系統中為每個核配置了一個固定的TC 用來傳輸數據，這樣在多核同時傳輸數據時可以提高效率。

2.2 核間同步

主核和從核可以通過IPC 和Semaphore 進行同步的，當需要開始執行一項任務時，主核根據任務類型給從核發送帶有不同中斷信息的IPC 中斷，然后主核開始進行自身的計算。從核X 接收到中斷后，首先請求占有Semaphore X，然后根據不同的中斷信息執行不同的任務，當任務執行結束后，釋放Semaphore X。主核執行完自身計算后，通過檢查Semaphore 1-7 的狀態來判斷從核的計算是否完成，待所有核計算完成一個任務后，主核開始進行下一任務的分發計算。

2.3 任務分配

本系統中任務分配的原則是計算量小的，不利于并行計算的以及時間能夠滿足要求的任務由一個核進行計算，對于計算量比較大的且可以并行計算的任務，分發給八個核同時進行計算。分發的原則是根據CPI號或者距離單元號進行的，具體來說，如果是根據CPI號進行分發，那么核X 執行CPI(32 ×X)至CPI(32 ×X+31)之間的數據;如果是根據一維距離像中距離單元號進行的，那么核X 執行距離單元(64 × X)至距離單元(64 × X+63)之間的數據。

本系統詳細的DPS 實現流程圖如圖2 和圖3所示，其中圖2 為主核流程圖，圖3 為從核流程圖(虛線框內的內容只有核1 有，其他從核沒有，因為虛線框中流程耗時很少，在一個CPI時間內即可完成，沒有必要進行并行處理)。圖中DNUM 代表核的編號。