【正文】 加實現(xiàn)高分辨成像，其基本原理可以用點目標結(jié)合線性陣列模型來解釋。 假定 SAR 平臺沿直線勻速運動采集數(shù)據(jù)，由此采樣點可以構造成一空間直線陣列，SAR 的雙程傳輸模型可以等效為目標向天線陣列的輻射模型，具體參見圖43。 圖43 SAR的回波歷程圖44 BP算法的聚焦過程根據(jù)圖43中SAR的幾何關系，可以得到目標到天線陣列的時延曲線(如圖44所示)，具體的時延表示式為： （436）為了簡便起見，發(fā)射信號直接表示為沖激信號 ，那么SAR的時域回波數(shù)據(jù)為： （437）BP算法的基本思想是通過計算目標點到各個天線孔徑點之間的雙程時延，而后沿相應的時延曲線進行相干疊加(見圖 (b))，求疊加結(jié)果的幅度得到所需目標點的后向散射強度，即： （438）對成像區(qū)域中的各個點遍歷上述“延遲—求和”過程，即可得到每個點的散射強度值。當成像區(qū)域的像素點與目標位置重合時，式(438)就會沿著目標的時延曲線積分，實現(xiàn)相干疊加。當像素點與目標位置不重合時，式(438)將無法進行相干疊加，從而就可得到目標的精確成像結(jié)果。 標準BP算法描述MIMO雷達的回波數(shù)據(jù)表示式是關于發(fā)射孔徑位置變量、接收孔徑位置變量和時間t的函數(shù)形式，在方位向上無論是信號包絡還是載頻相位項都與變量和有關，在距離時間向已經(jīng)通過匹配濾波進行了壓縮，常規(guī)的頻域SAR成像算法難以同時完成對發(fā)射孔徑和接收孔徑的方位聚焦。因此，MIMO雷達成像的首選處理方法就是不受陣列形式限制的BP算法。應用標準BP算法對MIMO雷達回波數(shù)據(jù)進行成像處理之前，首先對MIMO雷達回波數(shù)據(jù)的時延曲線進行分析。目標散射點 q對應的傳播時延可表示為： （439）MIMO 雷達方位向成像處理時需要同時聚焦收發(fā)孔徑的時延曲線，而傳統(tǒng)的RD算法、RM算法及CS算法等一般只能聚焦發(fā)射孔徑或是接收孔徑的時延曲線，所以它們難以實現(xiàn)MIMO雷達成像，但 MIMO雷達復雜的時延曲線特征對于標準BP算法卻不會產(chǎn)生影響。按照 BP 算法的基本原理，建立如圖45所示的成像處理幾何結(jié)構。圖中對成像區(qū)域進行了網(wǎng)格劃分，即通過均勻采樣將連續(xù)圖像離散化。區(qū)域劃分后像素點之間的間隔一般與分辨率相當，因為間隔太大會造成圖像的欠采樣，影響成像質(zhì)量，而間隔過小又會加大算法的運算量。圖 ，這只是為了方便繪制 MIMO 雷達成像處理的幾何結(jié)構，并不代表它們之間的實際分布關系。 假定成像區(qū)域劃分為K個像素點，以(k=1,2,…,K)和(l=1,2,…,L)分別表示距離向和方位向上像素點的坐標值，和則表示像素點（）到發(fā)射陣元 m和接收陣元n的距離，如圖45所示。圖45 MIMO雷達標準BP算法的成像處理幾何結(jié)構在球面波輻射接收模型下，像素點（）與收發(fā)天線位置組合（）之間的雙程時延為： （4310）距離壓縮后的MIMO雷達回波數(shù)據(jù)表示為： （4311）標準 BP算法需要對它完成方位向的聚焦處理。成像區(qū)域像素點（）的聚焦結(jié)果可以表示為： （4312）上述BP算法的成像處理包含三個過程：基頻包絡的時延調(diào)整， 載頻項的相位校正和數(shù)據(jù)相干疊加，具體的算法處理流程如圖46所示。從圖中可見，標準BP算法的處理流程簡單，相位補償容易，同時在成像過程中不存在任何假設條件的情況下對MIMO雷達成像的陣列形式?jīng)]有限制，是一種魯棒性很強的時域成像算法。但標準BP算法的處理過程需要逐點精確聚焦，這使得算法的運算量很大，這是標準BP算法的主要缺點，所以MIMO雷達成像使用標準BP算法時需要特別注意這個問題。 MIMO雷達回波數(shù)據(jù)劃分成像區(qū)域K選擇初始像素點信號分選(k,l)e（，t）確定像素點的時延曲線時間延遲e(,)(k+1,l)或(k,l+1)改變像素點位置相位校正e(，)否判斷是否遍歷整個成像區(qū)域沿發(fā)射孔徑和接收孔徑相干疊加是成像結(jié)果圖46 MIMO雷達標準BP算法流程圖 標準BP算法仿真實驗場地設在一棟大樓的樓頂，因此收發(fā)孔徑及目標場景距離的取值范圍都會受到一定的限制，但參數(shù)的選擇都滿足誤差分析的要求。為了便于對實驗系統(tǒng)參數(shù)進行描述，以天線導軌方向為x軸方向，以導軌的水平垂直方向為y軸方向建立二維直角坐標系。接收陣列的孔徑長度為1m，分布于x軸上的[ 0. 5 m， 0. 5 m] 區(qū)間；發(fā)射陣列包含3個發(fā)射陣元，它們的坐標分別為( 0. 7 m， 0 m) 、( 0 m，0 m) 和( 0. 7 m，0 m) 。目標的參考中心到天線導軌的距離為10m，相應的坐標為( 0 m，10 m)。多散射點目標包括4個散射點，它們的坐標分別為( 0. 25 m， 9. 5 m) 、 ( 0. 25 m，10. 5 m)、 ( 0. 25 m，9. 5 m)、 ( 0. 25 m，10. 5 m) 。發(fā)射信號的中心波長為0. 03 m，帶寬為200 MHz。利用標準BP算法對上述MIMO雷達實驗系統(tǒng)所錄取的實測數(shù)據(jù)進行成像處理，處理結(jié)果如圖47所示。圖47 多散射點標準BP算法成像結(jié)果第五章 總結(jié)與展望本文考慮到MIMO雷達成像技術的應用背景，針對MIMO雷達的成像基礎問題、成像算法等進行了深入的研究，并且對雷達正演做了模擬成像，主要內(nèi)容如下： 研究了MIMO雷達成像所涉及的數(shù)據(jù)獲取方式、陣列特性等基礎性問題。通過MIMO雷達成像的數(shù)據(jù)獲取方式研究，證明多觀測通道效應是MIMO雷達實時成像的基本保障條件。模擬了探底雷達正演的成像過程。設定目標模型的對應參數(shù)，利用GprMax軟件對目標體進行二維建模。研究了MIMO雷達數(shù)字波束形成技術。如果能夠?qū)ο嗫仃嚴走_的波束形成及自適應處理技術進行改進應用于 MIMO 雷達中，則能夠提高目標定位精度和多目標分辨能力，并能夠更有效地抑制雜波和干擾。研究了MIMO雷達成像BP算法?；赟AR BP算法原理，導出了MIMO雷達標準BP算法?？紤]到標準BP算法的實現(xiàn)效率問題，提出了MIMO雷達TCCBP算法，它是基于時延曲線校正的。結(jié)合SAR RD算法，提出了MIMO雷達RDBP算法，它具有高于標準BP算法和TCCBP算法的運算效率。從MIMO雷達概念提出至今，人們對MIMO雷達各方面關鍵技術的研究日趨深入，MIMO 雷達的多方面性能優(yōu)勢也逐步得到證實。由于MIMO雷達成像技術尚處于起步階段，其技術本身比較復雜，尚有許多理論和實際問題需要進一步深入研究： (1) 正交發(fā)射信號及其分選方法設計。設計適用于MIMO 雷達成像技術的發(fā)射信號波形及其分選方法需要進行研究。(2) 雷達正演成像的效果圖有待提高。(3) 文中提出了很多種算法技術，但是，各種算法都有自己的優(yōu)缺點，在今后的研究中，應該結(jié)合各種算法綜合考慮，逐步完善成像算法。 致謝四年的本科生學習即將結(jié)束，在這段時光里充滿了艱辛與快樂、挑戰(zhàn)與收獲，我得到了諸多老師、同學和朋友的關心、支持與幫助。值此論文完成之際，謹向他們表示衷心的感謝。參考文獻[1] 王懷軍，MIMO雷達成像算法研究，博士學位論文，國防科學技術大學，201009。[2] 雷文太、黃春琳、王懷軍等，合成孔徑與實孔徑雷達譜域成像算法對比分析，數(shù)據(jù)采集與處理，2010。[3] 王懷軍、許紅波、陸珉等，MIMO雷達技術與應用分析，雷達科學與技術，2009。[4] 王懷軍、黃春琳、陸珉等，MIMO雷達反向投影成像算法，系統(tǒng)工程與電子技術，2010。[5] 姜偉，MIMO雷達信號處理關鍵技術研究，博士學位論文，北京理工大學，2009年11月30日。[6] 黃玲，多輸入多輸出探地雷達方法研究，博士學位論文，吉林大學，2010年4月。[7] 黃韜、袁超偉、楊睿哲等，MIMO相關技術與應用，北京: 機械工業(yè)出版社，2006年3月4日。[8] 夏威，MIMO雷達模型與信號處理研究，博士學位論文，西安: 西安電子科技大學, 2007年3月8日。[9] 孔令講，淺地層探地雷達信號處理算法的研究，碩士學位論文，西安：西安電子科技大學，2003年。[10] 粟毅、黃春琳、雷文太，探地雷達理論與應用，北京: 科學出版社，2006年。[11] 曾昭發(fā)、劉四新、王者江、薛建，探地雷達方法原理與應用，北京: 科學出版社，2006年。[12] 王青，MIMO雷達陣列設計方法研究，碩士學位論文，國防科學技術大學研究生院，2011年7月。[13] 何子述、韓春林、劉波，MIMO雷達概念及其技術特點分析，電子學報，2005年。[14] Nikolaus Lehmann, Some Contributions on MIMO radar,Newark,New jersey Institute of Technology, Doctor thesis of Philosophy,2007.[15] James Irving，Rosemary Knight，Numerical modeling of groundpenetrating radar in 2D using MATLAB，Computers amp。 Geosciences 32，2006.28