互質陣列信號處理研究進展：波達方向估計與自適應波束成形

周成偉 鄭 航 顧宇杰 王 勇 史治國*

①(浙江大學控制科學與工程學院 杭州 310027)

②(浙江大學信息與電子工程學院 杭州 310027)

③(天普大學電子與計算機工程系 費城 19122)

1 引言

陣列信號處理通過布設的傳感器陣列對空域信號進行采樣和處理，以提取信號特征及其信息，在雷達、聲吶、語音、天文成像、無線通信等領域均有著廣泛的應用。波達方向估計和波束成形是陣列信號處理領域的基本問題，在理論、技術和應用層面均受到研究人員的廣泛關注。

以雷達領域為例，陣列信號處理在相控陣雷達、多輸入多輸出(Multiple Input Multiple Output,MIMO)雷達、機載雷達等各類應用中發(fā)揮著重要的作用，雷達技術的發(fā)展與陣列信號處理理論與技術的研究相輔相成、互相促進。舉例而言，近幾十年間，相控陣在雷達系統(tǒng)中逐步推廣應用，它通過大量可獨立控制的天線單元排列成天線陣面，并利用各單元形成不同的相位波束，在空間中輻射出具有不同方向性的波束[1]。相控陣具備電子掃描特性和相位可控的陣列天線結構，具有快速波束掃描的能力，能夠同時對多個目標進行搜尋和跟蹤[2]。通過在相控陣雷達體制的基礎上結合數(shù)字波束成形(digital beamforming)技術[3]，數(shù)字陣列雷達(digital array radar)實現(xiàn)了雷達系統(tǒng)的高效化運作[4]；與此同時，為了獲取相控陣雷達信號的角度和多普勒頻率聯(lián)合估計[5,6]，空域參數(shù)估計方面形成了一系列高精度算法，包括2維多重信號分類(MUltiple SIgnal Classification,MUSIC)方法[7]、2維子空間旋轉不變(Estimation of Signal Parameters via Rotational Invariance Techniques,ESPRIT)方法[8]以及應用于機載相控陣脈沖多普勒雷達中的基于分數(shù)低階矩方法[9]。相較于通過多個天線陣元發(fā)射相同波形的相控陣雷達，空間分集思想的引入使得能夠通過多天線陣元發(fā)射多種信號波形的MIMO雷達開始興起[10]。MIMO雷達利用多天線陣元結構同步發(fā)射分集波形，并同時使用多個天線接收回波信號進行收發(fā)信號的集中處理[11]，在目標檢測、參數(shù)估計和識別、分辨能力等方面均具有明顯優(yōu)勢[12–14]。為了在MIMO雷達中實現(xiàn)精確的目標方位識別，相關研究形成了基于ESPRIT方法的雙基地MIMO雷達系統(tǒng)到達角和離開角聯(lián)合估計[15]，基于張量模型的雙基地MIMO雷達角度估計和陣列互耦自校準方法[16,17]，基于車載MIMO雷達的低復雜度高分辨參數(shù)估計[18]，基于分布式MIMO雷達的移動目標和天線聯(lián)動定位[19]，面向MIMO雷達中波形非正交化情況下的新型波達方向估計[20,21]，以及MIMO雷達在存在未知空間有色噪聲及陣列誤差等非理想場景下的魯棒參數(shù)估計[22–25]等一系列理論與技術應用。相控陣雷達和MIMO雷達作為經典的雷達體制，在包括車載雷達、機載雷達、星載雷達等在內的諸多場景中得到了廣泛的推廣應用。其中，機載雷達通常采用空時自適應處理(Space-Time Adaptive Processing,STAP)技術檢測低空運動目標并抑制環(huán)境雜波[26–28]。STAP方法通過空域和時域二維聯(lián)合自適應濾波以實現(xiàn)雜波的有效過濾，典型研究包括基于STAP的無人機載雷達運動物體檢測[29]，基于稀疏字典原子選擇方法[30]，基于子陣列柵瓣雜波抑制方法[31]的機載MIMO雷達STAP算法，以及基于STAP的離網(wǎng)型角度/多普勒頻率聯(lián)合估計方法[32]。

無論是相控陣雷達、MIMO雷達還是機載雷達，其共同特點是基于傳感器陣列天線，并利用陣列信號處理方法實現(xiàn)目標檢測、方位估計及跟蹤等一系列功能。出于規(guī)則性陣列結構和奈奎斯特采樣速率的限制，均勻陣列是傳統(tǒng)雷達系統(tǒng)中最為常用的天線陣列結構。然而，為了滿足奈奎斯特采樣定理，均勻陣列中相鄰天線陣元的間距需不大于半波長，故陣列孔徑相對受限；此外，采用均勻陣列的方法自由度受限于天線陣元個數(shù)，無法在信號源數(shù)目大于天線陣元個數(shù)的條件下進行有效的信號處理。若通過增加天線陣元個數(shù)的方式來獲取更大的陣列孔徑和更高的自由度，實際系統(tǒng)的部署成本和算法復雜度將隨之大幅增加。因此，隨著應用背景的革新和雷達技術的迅猛發(fā)展，均勻陣列已逐漸無法滿足當前雷達領域應用日益增長的高效性和準確性等需求，傳統(tǒng)方法在性能和成本方面的矛盾亟待解決。

為了克服上述挑戰(zhàn)，美國加州理工學院課題組于2011年在文獻[33]中首次提出了互質采樣的構想和互質陣列的結構，奠定了互質陣列信號處理的理論基礎?；ベ|陣列是一種具有系統(tǒng)化結構的稀疏陣列，由一對陣元數(shù)滿足互質條件的稀疏均勻線性陣列構成?；ベ|陣列相較于傳統(tǒng)均勻陣列，主要具備3方面優(yōu)勢：一是互質陣列的稀疏陣元排布能夠實現(xiàn)入射信號的欠采樣，從而突破奈奎斯特定理對天線陣元間距的限制；二是陣列孔徑的擴展能夠有效提升分辨率性能；三是互質陣列能夠獲得遠超其物理陣元個數(shù)的自由度[33]，使得算法所能識別的信源數(shù)突破天線陣元數(shù)目的限制，從而節(jié)約了系統(tǒng)軟硬件成本開銷。為了充分利用互質陣列的上述優(yōu)勢以深入推進其在雷達系統(tǒng)中的應用，互質陣列的非均勻性及其信號模型的匹配問題亟待解決；為此，面向互質陣列的陣列信號處理理論與技術在近年來得到了廣泛的關注和研究。

本文介紹了當前互質陣列信號處理領域的研究進展，分別從波達方向估計和自適應波束成形角度回顧了最新研究工作。具體而言，在互質陣列波達方向估計方面，本文介紹了基于互質子陣分解和虛擬陣列信號處理等兩大類典型技術路線。其中，互質子陣分解方法通過探索稀疏陣列相位模糊的規(guī)則性，并利用質數(shù)的性質實現(xiàn)信號源的唯一性波達方向估計；虛擬陣列信號處理方法則利用增廣虛擬陣列所對應的2階等價信號統(tǒng)計量進行處理，相比于互質子陣分解方法具有更大的自由度。進一步地，本文分別從低復雜度和超分辨估計的角度介紹了互質陣列波達方向估計的最新研究：考慮到計算復雜度對于系統(tǒng)的實時性需求，本文介紹了基于壓縮感知技術的互質陣列波達方向估計算法；另一方面，由于早期的虛擬域奈奎斯特匹配方法未充分利用增廣虛擬陣列的信息，且存在基不匹配所導致的估計準確度受限問題，本文介紹了基于虛擬域陣元內插和無網(wǎng)格化技術的互質陣列波達方向估計算法。在互質陣列波束成形方面，本文在介紹互質陣列波束成形信號建模與工作原理的同時，指出了互質陣列信號處理框架下波束成形問題與波達方向估計問題的本質區(qū)別，并提出了面向非均勻互質陣列的波束成形器設計框架?；谏鲜隹蚣?，基于互質子陣分解的波達方向估計和基于互質子陣協(xié)方差矩陣聯(lián)合優(yōu)化的功率估計方法被提出，用于重建干擾加噪聲協(xié)方差矩陣和期望信號導引向量，并以此構造互質陣列波束成形器的權重向量。最后，本文對互質陣列信號處理方向的現(xiàn)有研究工作進行了總結，并從互質陣列結構優(yōu)化設計、互質陣列MIMO雷達以及非理想場景下的互質陣列魯棒參數(shù)估計等方向提出了研究展望。

2 互質陣列信號建模

在介紹互質陣列波達方向估計和互質陣列波束成形的研究之前，本文首先介紹互質陣列的結構及其信號模型，作為后續(xù)算法設計的模型基礎?；ベ|陣列由如圖1(a)所示的一對滿足互質條件的稀疏均勻線性子陣列構成，子陣列的陣元數(shù)分別為M和N，且陣元間距分別為Nd和Md。其中，M和N為互質整數(shù)，單位間隔d取為半波長。將上述兩個互質子陣以首個陣元疊加的方式進行線性組合，如圖1(b)所示，除首個陣元外的其余陣元均不重疊[33]，故互質陣列共包含M+N-1個陣元，各陣元位置的數(shù)集形式可表示為

考慮空間中有K個方向為θk，k=1,2,···,K的非相關信號源入射至互質陣列 S上，則t時刻的互質陣列接收信號可建模為

圖1 互質陣列結構示意圖Fig.1 Illustration of the coprime array structure

其中，AS=[a(θ1),a(θ2),···,a(θK)]∈C(M+N-1)×K為互質陣列導引矩陣，s(t)=[s1(t),s2(t),···,sK(t)]T為信號波形向量，n(t)為加性高斯白噪聲分量，[·]T表示轉置操作。對應于第k個入射信號的互質陣列導引向量a(θk)可表示為

其中，λ為信號波長，u? ∈S,?=1,2,···,M+N-1為互質陣列中各物理陣元的實際位置，?；ベ|陣列接收信號的協(xié)方差矩陣定義為

與均勻陣列相比，互質陣列中的天線陣元被稀疏放置，因而極大地增大了陣列孔徑，減少了陣元之間的互耦效應，進而為估計準確度和分辨率性能的提升奠定了基礎?；ベ|陣列具有系統(tǒng)化的稀疏陣列結構，在陣列設計的過程中，只需給定互質整數(shù)M和N就可以得到陣列架構，不需要通過復雜的優(yōu)化問題求解或遍歷搜索來確定陣列結構。同時，通過對陣列接收信號x(t)的2階統(tǒng)計量RS進行向量化處理與建模，可將互質陣列信號模型擴展至虛擬域實現(xiàn)陣列的增廣，并利用對應的2階等價接收信號進行虛擬域信號處理。虛擬域信號處理能夠有效擺脫傳統(tǒng)方法自由度受物理陣元個數(shù)的限制，僅通過M+N-1個天線陣元即可獲得高達O(MN)的自由度[33]。而與此同時，互質陣列本身的非均勻性不容忽視，故基于互質陣列的波達方向估計和波束成形方法需要充分考慮陣列本身的物理結構特征及其接收信號模型。

3 互質陣列波達方向估計

波達方向估計通過對陣列接收信號的統(tǒng)計處理實現(xiàn)入射信源相對方位角的估計，是雷達領域各類應用中的重要任務之一，也一直是國內外學術界、工業(yè)界研究的熱點問題[34–36]。當前，互質陣列波達方向估計的算法設計可分為互質子陣分解和虛擬陣列信號處理等兩類方法。其中，互質子陣分解方法具有計算復雜度低、估計精度高等優(yōu)勢，而虛擬陣列信號處理方法能夠在充分利用互質陣列大孔徑所帶來的分辨率優(yōu)勢的同時，突破奈奎斯特采樣定理的限制，實現(xiàn)波達方向估計自由度的增加。

3.1 基于互質子陣分解的波達方向估計

互質子陣分解方法的思路是將互質陣列分解為兩個滿足互質條件的稀疏均勻子陣列以分別進行波達方向估計，并利用質數(shù)的特性分析相位模糊的規(guī)律性，以獲得對應于每一個信號源的唯一波達方向估計結果[37]。在文獻[37]提出的“DECOM”算法中，傳統(tǒng)的MUSIC思想被分別應用到一對互質稀疏均勻子陣列所對應的信號模型中，并通過譜峰搜索對生成的兩個MUSIC譜分別進行搜索。由于稀疏子陣列的陣元間距大于半波長，波達方向估計存在相位模糊，圖2給出了M=7和N=5條件下兩個子陣MUSIC空間譜的相位模糊示意圖。

圖2 互質子陣列的MUSIC空間譜相位模糊示意圖Fig.2 Phase ambiguity of the pair of coprime subarray in the MUSIC spatial spectrum

由于這對稀疏均勻子陣列滿足互質的條件，其所對應的相位模糊具有一定的規(guī)律性。具體而言，對陣元數(shù)目為M的稀疏均勻子陣列，通過譜峰搜索可得到N個估計結果，從第k個信源得到的真實角度θk和相位模糊角度之間存在如式(5)所示關系

其中，D為非0整數(shù)。N個波達方向估計結果中除包含一個對應于真實波達方向的估計結果之外，還存在其他N-1個相位模糊角度；相應地，對于陣元數(shù)目為N的稀疏均勻子陣列，其相位模糊角度的規(guī)律性與上述結論類似。文獻[37]指出，在兩個稀疏均勻子陣列分別獲得的波達方向估計結果中，有且僅有一個相同的結果，即真實的互質陣列波達方向估計結果。上述結論可利用質數(shù)的性質予以證明：假設除以外，還存在一個在兩個子陣列分別獲得的波達方向估計結果中是相同的，則陣元數(shù)目為M的子陣列滿足關系為

其中，。類似地，對陣元數(shù)目為N的子陣列，存在關系DM=-(N-1),-(N-2),···,-1,1,···,N-2,N-1

其中，DN=-(M-1),-(M-2),···,-1,1,···,M-2,M-1，則DM/N=DN/M。由于M與N為互質整數(shù)，無法找到滿足上述等價條件的參數(shù)對{DM,DN}，故不存在，為唯一的波達方向估計結果。在實際應用中，由于存在噪聲等外部因素的干擾，在兩個子陣所對應MUSIC空間譜中所對應的值并非一致；為此，可通過尋找兩個空間譜中最為接近的兩個峰值，并將其所對應的角度求平均以獲得唯一的互質陣列波達方向估計結果。

上述互質子陣分解方法雖然克服了相位模糊問題并有效降低了軟硬件開銷，但是，DECOM算法作為傳統(tǒng)MUSIC方法在互質陣列上的衍生與應用，最終的角度估計精度仍然取決于譜峰搜索的間距，且子陣列分解帶來的互信息損失也會在一定程度上降低波達方向估計的效果；此外，由于該方法將互質陣列分解成兩個子陣單獨進行計算，可辨識的目標數(shù)與傳統(tǒng)均勻陣列相比將減少至少1/2以上?；谏鲜龌ベ|子陣分解的技術思想，文獻[38]提出了一種通過局部譜搜索代替全局譜搜索的方法，進一步降低了計算復雜度；與此同時，應用Root-MUSIC和ESPRIT等方法形成的閉式解可有效避免譜峰搜索過程帶來的精度受限問題[38–40]；文獻[41]提出的快速波達方向估計方法在保證自由度性能的前提下，進一步減少互質陣列的物理陣元個數(shù)，擴大了陣列孔徑，提高了估計精度。此外，為了提升子陣分解方法的自由度性能，通過對子陣列空間譜進行聯(lián)合計算，同時利用了兩個子陣的自信息和互信息，可以有效解決互信息損失問題，從而實現(xiàn)子陣列分解技術框架下自由度無損的波達方向估計[42,43]。

互質子陣分解算法遵循將互質稀疏均勻子陣和傳統(tǒng)均勻陣列波達方向估計方法相結合的思路，通過質數(shù)性質的應用實現(xiàn)波達方向的有效估計。得益于互質陣列的稀疏結構特性，該類方法在分辨率和計算復雜度方面優(yōu)于傳統(tǒng)采用均勻陣列的波達方向估計算法，實現(xiàn)方法簡單、易操作，適用于目標數(shù)量有限且對估計精確度要求較高的應用場景。

3.2 基于虛擬陣列信號處理的互質陣列波達方向估計

虛擬陣列信號處理方法通過計算互質陣列差集數(shù)組(difference coarray)在虛擬域上形成一個增廣虛擬陣列，并利用其2階等價虛擬域信號的處理，實現(xiàn)面向互質陣列的有效波達方向估計。由于虛擬陣元個數(shù)大于物理陣元個數(shù)，基于虛擬陣列信號處理的方法在算法自由度性能方面有顯著的提升。其中，虛擬域奈奎斯特匹配方法通過在虛擬域中引入傳統(tǒng)用于均勻陣列的空間平滑和多重信號分類技術，在提升自由度的同時實現(xiàn)超分辨波達方向估計；基于稀疏恢復(sparse recovery)的方法通過對虛擬域統(tǒng)計量進行空間域過完備表示，并通過稀疏條件約束下的優(yōu)化重建，實現(xiàn)自由度增加的波達方向估計；基于虛擬陣列協(xié)方差矩陣稀疏重建的方法在前述方法的基礎上，對入射信號的波達方向、功率和信源數(shù)等重要參數(shù)實現(xiàn)全面而精確的估計。

3.2.1 虛擬域奈奎斯特匹配和稀疏恢復方法

文獻[44]提出的虛擬域奈奎斯特匹配方法，通過推導增廣虛擬均勻陣列結構和對應的等價虛擬域信號模型，引入空間平滑技術和虛擬域多重信號分類處理，得到超分辨的波達方向估計結果。相較于上述互質子陣列分解方法，該方法不需要對互質陣列進行分解，故估計結果不存在因稀疏子陣列而帶來的相位模糊問題；同時，由于利用了對應于增廣虛擬陣列的虛擬域信號處理，算法的自由度性能得到了顯著的提升。

在該方法中，互質陣列 S通過差集數(shù)組計算的方式被推導至一個增廣虛擬陣列 V，該虛擬陣列可表示為

為了克服非連續(xù)虛擬陣列對于傳統(tǒng)奈奎斯特方法造成的信號模型失配問題，文獻[44]采用包含2M+N-1個物理陣元的擴展互質陣列，該擴展互質陣列的差集數(shù)組可表示為

該數(shù)組中包含了一個虛擬陣元位置由-MNd到MNd的連續(xù)子集，意味著物理陣元個數(shù)為2M+N-1的擴展互質陣列可以獲得2MN+1的自由度。根據(jù)式(2)和式(4)的定義，對于空間中K個來自θ=[θ1,θ2,···,θK]T方向的遠場非相關窄帶信號，得到互質陣列在t時刻的接收信號x(t)及協(xié)方差矩陣RS。通過向量化陣列接收信號的協(xié)方差矩陣RS，推導出實際接收信號x(t)與構造的虛擬陣列等價信號z之間的數(shù)學映射關系

其中，Ri為第i個虛擬子陣列所對應的協(xié)方差矩陣。最后，將傳統(tǒng)多重信號分類方法應用于空間平滑后得到虛擬域協(xié)方差矩陣，通過在虛擬域匹配奈奎斯特方法實現(xiàn)精確的波達方向估計。在類似的框架之下，通過改用結合Root-MUSIC,ESPRIT等其他子空間類算法，可以避免獲取角度結果時的譜峰搜索過程[45,46]。

虛擬域奈奎斯特匹配方法通過建立互質陣列接收信號1階統(tǒng)計量和增廣虛擬陣列2階等價信號之間的關系，利用空間平滑方法得到增廣虛擬陣列所對應的滿秩協(xié)方差矩陣，并結合多重信號分類處理實現(xiàn)互質陣列超分辨波達方向估計，大幅度地提升了算法自由度。然而，在該方法中，空間平滑方法將導致虛擬陣列孔徑的減小，并且非連續(xù)部分的虛擬陣元在平滑過程中會被忽略，增廣虛擬陣列的信息沒有得到充分利用。相對應地，基于稀疏恢復(sparse recovery)思想的虛擬域陣列信號稀疏重建(Sparse Signal Reconstruction,SSR)算法[47–52]在稀疏約束條件下，對過完備基表示的空間功率譜進行優(yōu)化重建以實現(xiàn)波達方向的估計，因該方法不需要一個對應于均勻線性陣列的虛擬域信號模型，且不需要經過空間平滑步驟，即可得到優(yōu)化的空間功率譜，故能夠充分利用非均勻虛擬陣列的全部信息。在稀疏恢復思想的基礎上，廣義化互質陣列結構[53]和Toeplitz先驗矩陣結構信息[54]的引入可進一步提升互質陣列波達方向估計的精確性和高效性。進一步地，基于稀疏恢復的方法成功實現(xiàn)了在平面互質陣列結構下的2維空間角度估計[55]，在應用上，稀疏恢復方法也被用于進行MIMO雷達的連續(xù)非相關目標定位檢測[56]。

對應于用預定義空間網(wǎng)格點過完備表示的擴展互質陣列導引矩陣，對角矩陣包含預定義空間網(wǎng)格點中個波達方向上所對應的信號功率，為(2M+N-1)×(2M+N-1)維單位矩陣。從而，信號稀疏恢復方法轉換為以下的優(yōu)化問題

3.2.2 基于虛擬陣列協(xié)方差矩陣稀疏重建的互質陣列波達方向估計算法

前述的虛擬域奈奎斯特匹配和信號稀疏重建方法實現(xiàn)了信號源數(shù)目大于物理陣元數(shù)目條件下波達方向的有效估計，但上述方法生成的空間譜普遍存在額外的虛峰。由于真實信源的數(shù)目在大多數(shù)實際應用中是未知的，故虛峰的辨識較為困難，容易造成信源的辨識偏差進而影響波達方向估計的精確度。此外，基于子空間類方法所生成的偽空間譜無法實現(xiàn)各個信號源的功率估計。為了解決這些問題，文獻[57,58]提出了基于虛擬陣列協(xié)方差矩陣稀疏重建的互質陣列波達方向估計算法，其思想是在虛擬域增廣陣列模型的框架下，利用預定義網(wǎng)格點對虛擬陣列協(xié)方差矩陣中的波達方向進行過完備表示，并在虛擬域通過協(xié)方差矩陣稀疏重建的方式設計優(yōu)化問題，最后實現(xiàn)信源個數(shù)、波達方向和信號功率等重要參數(shù)的精確估計。

與3.2.1小節(jié)所介紹的信號稀疏重建方法思路類似，文獻[57]中引入包含個波達方向的過完備預定義網(wǎng)格點對互質陣列接收信號x(t)中的波達方向參數(shù)進行過完備表示。在虛擬域上，上述操作則等效于對式(11)所定義的增廣虛擬陣列滿秩協(xié)方差矩陣Rss進行過完備稀疏表示，得到對應的稀疏化虛擬陣列協(xié)方差矩陣可表示為

對應一個陣列孔徑為MNd的虛擬均勻線性陣列導引矩陣，為(MN+1)×(MN+1)維單位矩陣，稀疏對角矩陣包含預定義空間網(wǎng)格點中個波達方向上對應的信號功率。

基于上述定義，虛擬陣列協(xié)方差矩陣稀疏重建優(yōu)化問題的基本思想為最小化空間平滑矩陣Rss與其對應的過完備表示虛擬陣列協(xié)方差矩陣之間的擬合誤差，從而尋得最優(yōu)化表示的稀疏空間功率譜，該優(yōu)化問題表示為

在此基礎上，文獻[58]設計了基于滑動窗口思路的信源枚舉方法，對稀疏空間譜中存在的虛峰響應進行移除，以獲得對應于個入射信號的波達方向估計。與此同時，以為先驗信息，可將稀疏重建的框架作用于互質陣列接收信號的協(xié)方差矩陣，形成增強型的功率估計優(yōu)化問題

基于虛擬陣列協(xié)方差矩陣稀疏重建的互質陣列波達方向估計算法在虛擬域信號處理的框架下實現(xiàn)了波達方向的精確估計，并大幅提升了算法的自由度性能；其次，通過信源枚舉步驟排除了現(xiàn)有方法中存在的空間譜虛峰問題，有效避免了因虛峰擾動所造成的估計精度受限問題；同時，通過增強型功率估計優(yōu)化問題設計，進一步保證了信源功率估計的準確性。圖3利用最優(yōu)子模式分配(Optimal Sub-Pattern Assignment,OSPA)準則[59]來對比前述算法與SSR[52]算法在信號源個數(shù)不確定情況下的波達方向估計性能，OSPA指數(shù)越低，表明估計結果越理想。仿真中設置一對互質線性子陣列分別包含2M=2×3=6和N=5個物理陣元，擴展互質陣列包含2M+N-1=10個物理陣元，假設空間中有12個均勻分布于–60°～60°方向的入射信號源，預定義空間網(wǎng)格點均勻分布于[–90,90°]，且采樣間距0.1°。圖3(a)表示采樣快拍數(shù)設置為500條件下各算法OSPA與輸入信噪比之間的性能對比，圖3(b)表示信噪比設置為0 dB條件下各算法OSPA與采樣快拍數(shù)之間的性能對比。從圖3的對比結果可知，所提算法在不同信噪比(Signal Noise Ratio,SNR)和不同采樣快拍數(shù)(snapshots)條件下的OSPA均低于SSR算法，且隨信噪比和采樣快拍數(shù)的增加呈現(xiàn)下降趨勢；而SSR算法由于在空間譜中存在若干隨機且不規(guī)則的虛峰，故OSPA較高，且在不同的采樣快拍數(shù)和信噪比條件下性能類似。上述結果表明，基于虛擬陣列協(xié)方差矩陣稀疏重建的互質陣列波達方向估計算法在估計性能上具有顯著的優(yōu)越性。

圖3 波達方向估計性能對比Fig.3 Comparison of DOA estimation performance

3.3 互質陣列波達方向估計算法的延伸發(fā)展

3.2節(jié)所述的兩大類方法都是在虛擬域信號處理的框架下實現(xiàn)面向互質陣列的高效波達方向估計，在自由度性能方面表現(xiàn)出較大的優(yōu)勢。本節(jié)將在上述研究的基礎上，分別從算法的高效性和估計的精確度等兩方面介紹互質陣列波達方向估計算法的深入延伸研究。為了進一步提升算法的高效性，利用壓縮感知技術對互質陣列接收信號進行壓縮，并基于壓縮信號進行處理，可在保留互質陣列性能優(yōu)勢的同時大幅降低計算復雜度；為了進一步提高算法的估計精度，利用無網(wǎng)格化方法可以實現(xiàn)波達方向估計參量的連續(xù)化表示，進而從根本上杜絕傳統(tǒng)稀疏類方法中因過完備網(wǎng)格點所造成的估計精度受限問題。

3.3.1 基于壓縮感知的低復雜度互質陣列波達方向估計算法

為了進一步降低波達方向估計算法的計算復雜度，壓縮感知技術被應用到陣列信號處理算法中，通過對互質陣列接收信號中的冗余信息進行壓縮，實現(xiàn)欠采樣條件下的高效信號處理[60–64]。文獻[64]提出采用一個隨機壓縮感知核Φ ∈CW×(M+N-1)對互質陣列接收信號x(t)∈C(M+N-1)×1進行降維處理，通過隨機投影的方式將(M+N-1)×1維的互質陣列接收信號x(t)壓縮成一個W×1維的壓縮信號y(t)

其中，Φ中的元素可通過高斯分布、伯努利分布等隨機產生，W <

壓縮信號y(t)的理想?yún)f(xié)方差矩陣Ry為

利用協(xié)方差矩陣Ry計算Capon空間功率譜，可以得到式(19)的空間譜函數(shù)

其中，θ為譜峰搜索網(wǎng)格點上假定的波達方向，為壓縮信號的采樣協(xié)方差矩陣，壓縮信號對應的W×1維導引向量ac(θ)表示為ac(θ)=Φa(θ)。進而，通過搜索Capon空間功率譜pCapon(θ)中最大的K個譜峰，即可得到波達方向角度估計值。

進一步地，為了保留互質陣列在自由度方面的優(yōu)勢，將壓縮感知技術推廣至虛擬域信號處理方法中，通過稀疏信號重建的方式實現(xiàn)自由度提升的波達方向估計。具體而言，將壓縮信號的協(xié)方差矩陣進行向量化操作，以構造等價虛擬擴展陣列及其2階信號統(tǒng)計量

上述算法通過引入壓縮感知技術，通過降維壓縮信號的Capon功率譜計算以及虛擬域信號稀疏重建，在大幅度降低計算復雜度的同時，分別實現(xiàn)了高分辨率和自由度增加的波達方向估計。此外，通過優(yōu)化壓縮感知信號的結構[65,66]，可以有效降低信號壓縮過程造成的原始信息損失，例如，通過對壓縮感知核的結構進行優(yōu)化，使其元素基于一定的先驗信息概率分布進行生成[66]，而不是采用隨機生成的方式[64]，可以減少信號降維壓縮過程中的信息損失，從而使算法具備更高的精確性和更低的計算復雜度。

3.3.2 基于虛擬域陣元內插的無網(wǎng)格化互質陣列波達方向估計算法

現(xiàn)有方法在估計精度性能上主要受以下兩方面的因素限制：一是由于虛擬陣列的非均勻性所導致的信息損失；二是由于預定義空間網(wǎng)格點在稀疏過完備表示或譜峰搜索過程中所造成的固有估計誤差。

為了避免信息損失問題，基于虛擬域信號稀疏重建的互質陣列波達方向估計算法對采樣協(xié)方差矩陣向量化后得到的2階等價虛擬信號進行了整體重建[52]，但是構造的大量重復虛擬陣元使得計算負擔變重，且虛擬陣列本身的非均勻性將導致空間功率譜中出現(xiàn)不規(guī)則的虛峰響應?；诙囝l率的互質陣列波達方向估計算法利用等價頻率填充差集數(shù)組中的缺失元素[67]，保證了在充分利用非均勻虛擬陣列提供的全部自由度基礎上實現(xiàn)高精度波達方向估計。進一步地，基于核范數(shù)最小化的互質陣列波達方向估計算法通過利用矩陣填充準則恢復虛擬陣列協(xié)方差矩陣中孔洞位置所對應的缺失元素，且不需要引入自定義用戶參數(shù)即可將非均勻虛擬陣列中所包含的全部信息用于波達方向估計[68]。文獻[69]改進了用于矩陣填充的最小已知虛擬陣元個數(shù)。

另一方面，基于稀疏恢復的互質陣列波達方向估計算法通過預定義網(wǎng)格對波達方向進行過完備參數(shù)化表示，但是實際入射信號的波達方向不一定會正好落在這些預定義空間網(wǎng)格點上，從而導致基不匹配問題，進而影響波達方向的估計精度。盡管通過將預定義空間網(wǎng)格點密集化可以在一定程度上減小稀疏網(wǎng)格帶來的匹配誤差[70]，但是過高的網(wǎng)格密度將導致計算復雜度的指數(shù)型提升；文獻[71]提出了一種利用凸優(yōu)化方法進行線性頻譜稀疏恢復的數(shù)學理論，并基于該理論建立起基于無網(wǎng)格化稀疏恢復的波達方向估計算法框架[72]；文獻[73,74]總結了當前幾種有效的離網(wǎng)型(off-grid)和無網(wǎng)格化(gridless)波達方向估計方法，例如基于稀疏貝葉斯學習的離網(wǎng)型波達方向估計[75]，基于稀疏恢復思想的互質陣列離網(wǎng)型2維波達方向估計[76]和結合原子范數(shù)最小化方法的無網(wǎng)格化稀疏恢復波達方向估計算法[77]。

為了同時解決信息損失和基不匹配問題，文獻[78–81]結合上述虛擬陣列填充和無網(wǎng)格化方法的思路，提出了基于內插虛擬陣列協(xié)方差矩陣無網(wǎng)格化重建的互質陣列波達方向估計算法。該方法在虛擬域擴展陣列的框架下，通過引入陣元內插的思想對增廣非均勻虛擬陣列中的非連續(xù)部分進行填充，將非均勻虛擬陣列轉化為一個虛擬均勻線性陣列Vc，如圖4所示。內插虛擬陣列 Vc既包含了非均勻虛擬陣列中全部的虛擬陣元，對應的等價信號又滿足奈奎斯特采樣速率的要求。將內插虛擬陣列的2階等價信號表示為通過設計優(yōu)化問題恢復中各內插虛擬陣元上的未知信號，以實現(xiàn)波達方向的有效估計。

在內插虛擬陣列信號模型的基礎上，原子范數(shù)的引入能夠實現(xiàn)波達方向這一參數(shù)的連續(xù)化表示，進而避免傳統(tǒng)稀疏方法帶來的基不匹配問題。為了實現(xiàn)虛擬域信號的原子范數(shù)表示，本文首先將理想條件下內插虛擬陣列2階等價接收信號建模為

隨后，為了對波達方向進行連續(xù)化參數(shù)表示，定義ZV的原子集合為

其中，e(θ)=[1,e-jπsin(θ),···,e-jπ(E-1)sin(θ)]T為E個虛擬子陣列和參考虛擬陣列之間的相位偏移。相應地，將ZV的原子范數(shù)表示為

圖4 互質陣列及其對應的各種虛擬域陣列結構示意圖Fig.4 Illustration of the coprime array and its corresponding virtual array structures

其中，conv(?)表示原子集合?的凸包(Convex hull)。式(23)定義的多采樣虛擬信號ZV的原子范數(shù)提供了一個2階等價虛擬信號的無網(wǎng)格化表示模型，實現(xiàn)了波達方向這一參數(shù)的連續(xù)化表示。進一步地，文獻[79]根據(jù)其等價半正定規(guī)劃表現(xiàn)形式及其對應的推論揭示了內插虛擬陣列協(xié)方差矩陣與2階等價虛擬信號之間的關系，由此通過原子范數(shù)的最小化實現(xiàn)了內插虛擬陣列協(xié)方差矩陣T(l)的無網(wǎng)格化重建。具體而言，T(l)的首列向量l與首個虛擬陣列對應的2階等價虛擬信號相同，其原子范數(shù)可表示為

基于虛擬均勻線性陣列理想?yún)f(xié)方差矩陣具有的低秩特性和厄米特Toeplitz結構[82–84]，該方法將對應于非均勻虛擬陣列的2階相關統(tǒng)計量作為參考值，結合矩陣重建的思想，構建了基于原子范數(shù)最小化的內插虛擬陣列協(xié)方差矩陣Toeplitz化重建的優(yōu)化問題，即通過最小化式(24)中向量l的原子范數(shù)，重建內插虛擬陣列協(xié)方差矩陣T(l)，具體表示為

基于內插虛擬陣列協(xié)方差矩陣無網(wǎng)格化重建的互質陣列波達方向估計算法具備以下優(yōu)點：一是算法中重建的內插均勻虛擬陣列包含了非均勻虛擬陣列的全部信號信息，從而能在充分利用互質陣列特性的基礎上對波達方向作出高精度估計；二是重建的內插虛擬陣列協(xié)方差矩陣滿足半正定厄米特Toeplitz結構，該結構與均勻線性陣列接收信號所對應的理想?yún)f(xié)方差矩陣一致；三是算法采用了虛擬域原子范數(shù)的定義和波達方向參數(shù)連續(xù)化表示方法，從而能夠無網(wǎng)格化地重建內插虛擬陣列協(xié)方差矩陣，有效避免了基不匹配問題。圖5采用均方根誤差(Root Mean Square Error,RMSE)作為衡量波達方向估計性能的參數(shù)，通過仿真對比結果以說明上述結論，相對比的算法包括：基于協(xié)方差矩陣稀疏重建的互質陣列波達方向估計算法[57,58](Covariance Matrix Sparse Reconstruction,CMSR)，基于稀疏信號重建的互質陣列波達方向估計算法[52](Sparse Signal Reconstruction,SSR)，克拉美羅界(Cramer-Rao Bound,CRB)也同時在仿真圖中給出。仿真設置中，構成互質陣列的一對互質整數(shù)選取為M=3和N=5，單個信源入射角度由隨機方式生成，當以信噪比為變量進行均方根誤差對比時，采樣快拍數(shù)設置為500，當以采樣快拍數(shù)作為變量進行均方根誤差對比時，信噪比固定為20 dB。從圖5所示對比結果可知，在單個入射信號源情況下，CMSR算法和SSR的RMSE在信噪比大于10 dB的情況下趨向平緩，相比之下，所提基于內插虛擬陣列協(xié)方差矩陣無網(wǎng)格化重建算法由于在波達方向估計過程中無需使用預定義的空間網(wǎng)格點，故波達方向估計性能不受空間網(wǎng)格點間距的限制，RMSE曲線隨著信噪比的增加而不斷降低，且與CRB的趨勢保持一致。在算法復雜度方面，由圖6的對比結果可知，基于內插虛擬陣列協(xié)方差矩陣無網(wǎng)格化重建算法在不同采樣間隔設置下均保持相對穩(wěn)定，且計算復雜度低于其他對比算法。

圖5 單個隨機信號源情況下的波達方向估計性能對比Fig.5 Comparison of DOA estimation performance under the single random source scenario

圖6 算法復雜度性能對比Fig.6 Comparison of algorithm complexity

4 互質陣列自適應波束成形

波束成形通過抑制干擾信號并增強期望信號(Signal of Interest,SOI)，以實現(xiàn)信號的定向傳送，從而避免不必要的開銷，保證整體系統(tǒng)效率。作為陣列信號處理領域的又一個基本問題，魯棒性和高效性是波束成形算法設計中的關鍵考量因素。其中，魯棒性要求在外部環(huán)境擾動、陣列結構失準、信號模型失配等非理想條件下保證波束成形器輸出性能的穩(wěn)定，而高效性則要求算法具備快速精確的處理能力，并同時能夠節(jié)省系統(tǒng)的軟硬件成本開銷。因此，設計兼?zhèn)漪敯粜院透咝缘牟ㄊ尚嗡惴ǎ瑢τ诋斍芭畈l(fā)展的雷達技術應用具有重大的意義和應用價值。

4.1 互質陣列波束成形設計準則

均勻陣列是傳統(tǒng)波束成形算法中最為常用的陣列結構，由于波束主瓣與陣列孔徑成反比，奈奎斯特采樣定理的限制使得傳統(tǒng)算法需要通過增加物理陣元數(shù)目來實現(xiàn)波束的精尖化，進而導致了硬件成本和計算復雜度的大幅增加。相對地，互質陣列的稀疏陣列架構、大孔徑特性和高自由度等特性能夠在相同的開銷下獲得更好的波束控制效果，從而實現(xiàn)更加高效的能量集中、更加精確的波束對準和更加強勁的干擾抑制能力。因此，采用互質陣列進行波束成形在魯棒性和高效性方面均具有較大潛在優(yōu)勢。

與前述互質陣列波達方向估計問題相比，互質陣列波束成形問題的設計準則與其在原理上存在本質區(qū)別。具體而言，由于波達方向估計屬于參數(shù)估計問題，故互質陣列波達方向估計可直接通過虛擬域信號處理以實現(xiàn)自由度的增加。值得注意的是，雖然2階等價虛擬陣列信號z與1階陣列接收信號x(t)在信號建模上存在類似的結構，但式(2)中的1階統(tǒng)計量包含信號波形s(t)，而式(10)中的虛擬域2階統(tǒng)計量則包含信號功率p。由于波束成形的原理是陣列接收信號的空域濾波，即對物理陣元接收信號的加權求和，故互質陣列波束成形無法直接在虛擬域中完成，其波束成形器的權重設計以及信號的權值累加必須在實際的物理陣元上實現(xiàn)。因此，互質陣列波束成形算法在設計準則上與前述互質陣列波達方向估計存在本質區(qū)別，如何有效利用包括虛擬域信號處理在內的一系列互質陣列優(yōu)勢特性，并同時克服因其物理陣列的非均勻性所帶來的一系列難題，是互質陣列自適應波束成形算法研究的核心問題。

相較于近年來被廣泛研究的互質陣列波達方向估計問題，互質陣列自適應波束成形的研究仍處于初步階段，尚有較多富有前景的研究課題有待進一步深入發(fā)掘。在介紹具體的互質陣列波束成形算法之前，本文首先回顧了傳統(tǒng)的波束成形技術，以說明下文所采用波束成形器設計思路的原因。最小方差無失真響應(Minimum Variance Distortionless Response,MVDR)準則[5]是理想化波束成形設計的經典準則之一，其目標是最大化陣列輸出的信干噪比(Signal to Interference plus Noise Ratio,SINR)，并同時保證期望信號分量的無畸變通過。但是，由于理想的干擾加噪聲協(xié)方差矩陣及期望信號導引向量在實際應用中未知，波束成形算法的性能易受SOI方向失配、陣列校準誤差、以及有限觀測樣本下采樣協(xié)方差矩陣估計誤差等問題帶來的性能衰減。為了緩解這一問題，一系列魯棒自適應波束成形算法被提出。文獻[85]總結了兩類典型魯棒自適應波束成形算法設計，它們分別通過對采樣協(xié)方差矩陣和假定的期望信號導引向量進行統(tǒng)計信號處理和優(yōu)化，以實現(xiàn)輸出SINR性能的提升。但是，由于干擾加噪聲協(xié)方差矩陣中的期望信號分量未被有效消除，信號自相消現(xiàn)象將隨SINR的增大而逐步凸顯。子空間方法[86,87]和對角加載方法[88]可以提高導向矢量失配情況下MVDR波束成形的穩(wěn)健性，但是基于子空間的波束成形算法在低信噪比或大信干噪比的情況下效果較差，而對角加載方法的缺點是最優(yōu)對角加載因子難以確定。此外，貝葉斯方法[89]和基于最壞情況性能優(yōu)化算法[90]也都沒有從本質上解決信號的自相消問題。

為了解決上述經典算法中存在的問題，文獻[91]提出了一種基于干擾加噪聲協(xié)方差矩陣重建和期望信號導引向量估計的魯棒自適應波束成形算法，首次實現(xiàn)了逼近理論值的輸出SINR性能。該算法的核心思想是通過干擾加噪聲協(xié)方差矩陣的重建來移除其中的期望信號分量，并同時設計了一個2次約束2次規(guī)劃問題以優(yōu)化期望信號導引向量，從而大幅度提升了自適應波束成形的魯棒性和輸出SINR性能。鑒于干擾加噪聲協(xié)方差矩陣重建和期望信號導引向量估計方法的性能優(yōu)勢，互質陣列波束成形擬采用該技術框架進行面向非均勻物理陣元的自適應波束成形權重設計。

4.2 基于協(xié)方差矩陣和導引向量重建的互質陣列自適應波束成形

基于前述自適應波束成形設計準則，文獻[92,93]在文獻[91]的思想上，充分考慮了互質陣列的結構特征和信號建模特性，利用互質陣列的大孔徑優(yōu)勢，提出了基于干擾加噪聲協(xié)方差矩陣和期望信號導引向量重建的互質陣列自適應波束成形算法。其中，文獻[93]給出了互質陣列自適應波束成形算法設計的完整框架。具體而言，對于一個包含M+N-1個陣元的互質陣列，假定有K+1個遠場窄帶信號入射至互質陣列，其中期望信號的角度方向為θs，其余K個來自{θk,k=1,2,···,K}方向的信號為干擾信號，w ∈CM+N-1為波束成形的權重向量。根據(jù)MVDR準則，波束成形權重向量可通過求解如下優(yōu)化問題獲得

上述優(yōu)化問題中權重向量的最優(yōu)解為

其中，a(θs)是對應于期望信號θs的導引向量，Ri+n為干擾加噪聲協(xié)方差矩陣，可以表示為

其中，a(θk)為對應各干擾信號入射方向θk的導引向量。

由式(27)可知，基于MVDR準則的波束成形器權重向量是由期望信號導引向量a(θs)和干擾加噪聲協(xié)方差矩陣Ri+n共同決定的，因此所提互質陣列波束成形算法的思路是充分利用互質陣列的特性以精確重建a(θs)和Ri+n這兩個組成波束成形器的核心參量，從而使所設計的互質陣列波束成形權重向量逼近理論值考慮到期望信號入射角度是期望信號導引向量重建的核心參數(shù)，且由式(28)可知各干擾信號的入射角度為干擾加噪聲協(xié)方差矩陣重建的必要參數(shù)，互質陣列自適應波束成形算法將首先對上述參數(shù)進行估計，以作為波束成形權重向量重建的先驗信息。

為了精確地獲得期望信號和各干擾信號的波達方向估計結果，文獻[93]采用3.1小節(jié)介紹的互質子陣分解方法，將互質陣列分解成包含M和N個陣元的稀疏均勻子陣列，并分別對子陣列的MUSIC空間譜函數(shù)進行譜峰搜索，得到包含相位模糊的波達方向估計結果，可表示為：。為了去除相位模糊，期望信號的波達方向可通過選取子陣估計結果中最接近的一組候選估計值實現(xiàn)，即

其中，φ為期望信號的角度范圍。期望信號的波達方向估計值由式(29)求得的最優(yōu)值取平均獲得

類似地，干擾信號的波達方向可通過求解式(31)得到

由式(28)關于干擾加噪聲協(xié)方差矩陣的定義可知，除了干擾信號的波達方向之外，其對應的信號功率k=1,2,···,以及噪聲功率也是實現(xiàn)該矩陣重建的必要參數(shù)，而前述所采用的子陣分解方法得到的MUSIC空間偽譜在各入射信號方向上的響應值無法表征功率。為此，文獻[93]將前述得到的期望信號和干擾信號波達方向估計值作為先驗信息，設計了基于互質稀疏子陣協(xié)方差矩陣聯(lián)合優(yōu)化問題以進行功率估計，其核心思想是考慮兩個互質稀疏子陣所對應的采樣協(xié)方差矩陣中所包含的信號功率信息一致，并構造如下不等式約束最小二乘問題

其中，λmin(·)表示矩陣的最小特征值。通過式(32)求得優(yōu)化變量Γ的最優(yōu)值包含了期望信號和個干擾信號的功率估計值。

其中，d1,d2,···,dM+N-1表示互質陣列中各物理陣元的實際位置，d1=0。同時，根據(jù)式(28)的定義，利用求得的干擾信號波達方向估計值和功率估計值，以及噪聲功率估計值，干擾加噪聲協(xié)方差矩陣可重建為

最后，結合式(27)關于最優(yōu)權重向量的定義和上述重建的參量，互質陣列波束成形的權重向量可表示為

上述互質陣列波束成形算法的設計方案充分利用了互質陣列的大孔徑特性，通過高精度的參數(shù)估計實現(xiàn)了參量的精確重建；與此同時，重建思想的引入有效克服了信號自相消問題，避免了輸出SINR在高信噪比情況下的性能衰減；其次，互質稀疏子陣協(xié)方差矩陣聯(lián)合優(yōu)化問題可估計所有先驗波達方向所對應的功率，以分辨真實信號源所對應的波達方向估計值，故所提算法無需將信號源個數(shù)作為先驗信息，對于實際應用中信號源個數(shù)未知情形具有明顯的優(yōu)勢；再者，互質陣列的應用使得所提波束成形算法在硬件開銷和計算復雜度方面均具有優(yōu)于傳統(tǒng)基于均勻陣列的波束成形算法。因此，上述互質陣列自適應波束成形算法兼?zhèn)漪敯粜院透咝浴?/p>

圖7給出了上述互質陣列自適應波束成形算法與傳統(tǒng)的采樣矩陣求逆(Sample Matrix Inversion,SMI)波束成形算法[94]、對角加載波束成形算法[88]、基于子空間的波束成形算法[87]、基于最壞情況性能優(yōu)化的波束成形算法[90]和基于協(xié)方差矩陣重建的波束成形算法[91]在輸出SINR特性上的性能對比。在圖7中，這些算法分別采用“PROPOSED”,“SMI”,”DLSMI”,“EIGENSPACE”,“WORST-CASE”和“RECONSTRUCTION”作為圖例說明。在實驗中，互質參數(shù)M和N選用M=5,N=6；當以期望信號的輸入信噪比為變量進行輸出SINR對比時，采樣快拍數(shù)設置為30；當以采樣快拍數(shù)為變量進行輸出SINR對比時，期望信號的輸入信噪比設置為20 dB。

圖7 觀測方向存在隨機誤差情況下的輸出性能對比Fig.7 Comparison of output SINR performance under the random signal look direction mismatch scenario

由圖7(a)所示結果可以看出，隨著輸入信噪比的增加，SMI波束成形算法、對角加載波束成形算法以及基于子空間的波束成形算法的輸出SINR曲線隨著信噪比的增加而趨向平緩。基于最壞情況性能優(yōu)化的波束成形算法盡管在信噪比較高的情況下性能優(yōu)于其他幾類傳統(tǒng)算法，但性能仍受信號自相消現(xiàn)象的限制。相較之下，基于重建思想的波束成形算法由于移除了干擾加噪聲協(xié)方差矩陣中的期望信號分量，從而有效地克服了信號自相消導致的輸出SINR性能衰減問題，并且具有與最優(yōu)SINR趨勢一致的輸出SINR性能，且所提互質陣列自適應波束成形算法性能優(yōu)于傳統(tǒng)采用均勻陣列的重建思想波束成形算法。從圖7(b)所示的對比結果可知，互質陣列自適應波束成形算法在輸出SINR性能上均優(yōu)于其余對比算法。

基于上述算法框架，文獻[95]在保證互質陣列輸出SINR性能的條件下，通過結合虛擬域信號處理實現(xiàn)了算法自由度的進一步提升；除此之外，考慮到波束成形算法在實際應用中的復雜度，壓縮感知等理論也被應用于自適應波束成形方法[96]中，例如基于3.1.1節(jié)的壓縮感知模型與MVDR準則的互質陣列波束成形算法，在虛擬域壓縮信號模型上實現(xiàn)了高效的波束成形設計[97]，相比傳統(tǒng)波束成形技術大幅降低了計算復雜度。

5 結束語

陣列信號處理是雷達領域各類應用中的核心技術之一，確保了雷達在不同體制、不同場景中實現(xiàn)目標檢測、定位、追蹤等功能。近幾十年來，傳統(tǒng)陣列信號處理算法的發(fā)展趨近成熟，為雷達領域提供了大量高效的算法模型。然而，隨著大規(guī)模檢測和超分辨定位等更高需求的提出，傳統(tǒng)均勻陣列帶來的性能受限問題開始凸顯，稀疏非均勻陣列的應用開始受到關注。其中，互質陣列由于其稀疏結構帶來的優(yōu)越特性，能為基于互質陣列結構的算法帶來在自由度、分辨率及復雜度等性能上的顯著優(yōu)勢，以突破傳統(tǒng)均勻陣列的技術瓶頸。為此，本文總結了大量基于互質陣列結構的陣列信號處理算法，主要包含波達方向估計和波束成形等兩大類方向。

在互質陣列波達方向估計方面，本文總結了基于互質子陣分解的方法和基于虛擬陣列信號處理的方法。其中，互質子陣分解方法通過子空間類方法的引入和互質稀疏子陣相位模糊規(guī)律性的探索，對互質稀疏子陣的估計結果進行分辨后得到唯一有效的波達方向估計結果，在有效利用大陣列孔徑帶來的高分辨率的同時，具有操作簡便，估計準確度高等優(yōu)勢；而基于虛擬陣列信號處理的方法則通過構造增廣虛擬陣列的2階等價虛擬信號，在虛擬域上進行信號處理以實現(xiàn)波達方向估計，能夠突破可分辨信號源個數(shù)受限于物理陣元數(shù)的瓶頸，在自由度性能方面具有明顯的優(yōu)勢。在虛擬陣列信號處理的框架下，主要介紹了奈奎斯特匹配、信號稀疏恢復、以及基于虛擬域協(xié)方差矩陣稀疏重建等3類代表性方法。其中基于虛擬域奈奎斯特匹配的方法將傳統(tǒng)采用均勻陣列的波達方向估計的思路擴展到了虛擬域，通過建立虛擬域擴展陣列信號與原始信號統(tǒng)計量之間的類比關系，實現(xiàn)了傳統(tǒng)奈奎斯特信號處理方法在虛擬域中的應用；基于稀疏恢復的方法則在虛擬域中引入壓縮感知的思想，通過重建虛擬陣列等價信號統(tǒng)計量直接獲得波達方向估計結果；與此同時，基于虛擬域協(xié)方差矩陣稀疏重建的互質陣列波達方向算法在保證自由度提升的條件下，能夠同時實現(xiàn)信號源波達方向及其功率的精確估計。上述算法充分利用了互質陣列的優(yōu)勢，實現(xiàn)了波達方向估計算法在自由度、分辨率、軟硬件開銷等層面的性能提升。進一步地，本文總結了兩種優(yōu)化算法思路，一種是針對計算復雜度優(yōu)化的基于壓縮感知技術的互質陣列波達方向估計算法，通過引入壓縮感知核對互質陣列接收信號進行核心信息提取與處理，在保證優(yōu)勢性能的同時大幅降低了計算復雜度；另一種是基于內插虛擬陣列協(xié)方差矩陣無網(wǎng)格化重建的互質陣列波達方向估計算法，它有效解決了由于增廣虛擬陣列非連續(xù)所導致的虛擬域信號處理信息損失問題，以及稀疏恢復類方法中因預定義網(wǎng)格點不夠密而導致的估計精度受限問題。

在互質陣列波束成形方面，本文回顧了波束成形技術的研究背景與發(fā)展歷程，分析了互質陣列波達方向估計與互質陣列波束成形問題在作用原理和統(tǒng)計量建模層面的本質區(qū)別，并構建了互質陣列波束成形方法的設計準則?；谏鲜鼋Y論，本文重點介紹了基于干擾加噪聲協(xié)方差矩陣和期望信號導引向量重建的互質陣列波束成形算法。該算法通過剖析MVDR準則獲得用于重建最優(yōu)權重向量的核心參數(shù)，即干擾源來波方向信息與信號源功率信息，并分別通過基于子陣分解的互質陣列波達方向估計、基于互質稀疏子陣協(xié)方差矩陣聯(lián)合優(yōu)化的功率估計實現(xiàn)上述核心參數(shù)的高效估計。其中，重建思想的引入消除了干擾加噪聲協(xié)方差矩陣中的期望信號分量，能夠有效克服信號的自相消問題；而信源功率的估計使得傳統(tǒng)波達方向估計問題中信源個數(shù)這一先驗條件得以去除；與此同時，面向非均勻互質陣列的權重設計有效利用了互質陣列的結構特性，最終實現(xiàn)了高效性和魯棒性兼?zhèn)涞母咝阅懿ㄊ尚巍?/p>

除波達方向估計和自適應波束成形這兩類問題之外，互質陣列信號處理仍有較多值得研究的話題及應用，例如：互質陣列結構優(yōu)化設計，非理想場景下的魯棒性算法，以及互質MIMO雷達等。在互質陣列結構優(yōu)化方面，廣義化的互質陣列結構設計[98,99]，互質采樣技術框架下的優(yōu)化陣列設計方案[100]，以及考慮保證相同虛擬陣列孔徑和自由度的精簡型互質陣列(thinned coprime array)方案[101]，都在互質采樣的框架下進行改進，以優(yōu)化互質陣列的稀疏結構，從而達到算法在自由度、計算復雜度等性能指標上的強化。與此同時，由于實際應用的需要，如何在非理想場景下保證互質陣列參數(shù)估計方法的有效性，是當前亟需解決的問題之一，這其中包含了克服陣列誤差影響[102]和排除未知噪聲干擾[103]等方面。除此之外，由于互質陣列在結構特性上的優(yōu)越性能，采用互質陣列結構的互質MIMO雷達[104,105]在空間信號分辨能力、識別精確度性能上都有大幅提升，面向互質MIMO雷達相關的結構設計和參數(shù)估計算法等方面的研究正在受到關注。針對互質MIMO雷達的結構設計，當前的研究包括廣義化的互質MIMO雷達模型[106]，以及在互質采樣技術的框架下對發(fā)射和接收端進行優(yōu)化的MIMO雷達陣列結構[107]等。面向互質MIMO雷達的算法研究不僅有結合硬件加速的參數(shù)估計方法[108–111]，也有大量高精度的測向算法研究[112–116]。因此，如何將互質陣列信號處理方法在雷達領域深入推廣應用，不斷優(yōu)化相關核心技術的算法性能并增強其高效性和實用性，是未來陣列信號處理領域的重要課題，也是推動雷達技術發(fā)展的核心工具之一。