摘" 要：

傳統(tǒng)波束成形算法在應(yīng)對(duì)維數(shù)較高的多維陣列時(shí)，訓(xùn)練快拍數(shù)難以滿足遠(yuǎn)大于信號(hào)維數(shù)這一要求，從而導(dǎo)致波束成形器性能急劇下降。針對(duì)這一問(wèn)題，本文給出了多維陣列輸出數(shù)據(jù)的張量模型，基于多維陣列子維度的可分離性，引入張量波束成形方法，分析了其在訓(xùn)練快拍數(shù)需求上具有的優(yōu)勢(shì)。然后，基于張量波束成形中子維度干擾協(xié)方差矩陣的模型，通過(guò)直接估計(jì)干擾協(xié)方差矩陣實(shí)現(xiàn)了一種魯棒張量波束成形方法。分析表明：該方法可更好地應(yīng)對(duì)非均勻雜波環(huán)境以及相干干擾，并得到更高的輸出信干噪比，以給出的二維極化敏感陣列為例，所提方法的訓(xùn)練快拍數(shù)需求為傳統(tǒng)方法的1/3，相干干擾下輸出信干噪比提高了約2.5 dB，仿真驗(yàn)證了分析的有效性。

關(guān)鍵詞：

魯棒波束成形; 張量波束成形; 多維陣列; 相干干擾

中圖分類(lèi)號(hào)：

V 243.2

文獻(xiàn)標(biāo)志碼： A""" DOI：10.12305/j.issn.1001-506X.2024.07.01

Robust beamforming method for multi-dimensional array tensor models

BI Quanyang， LI Dan*， ZHANG Jianqiu

（School of Information Science and Engineering， Fudan University， Shanghai 200433， China）

Abstract：

Traditional beamforming performance will significantly degrade with insufficient training samples on multi-dimensional arrays with much higher dimensions. In this paper， a tensor beamforming method that has advantages in the number of training snapshots required is introduced based on the proposed tensor model and the separability of sub-dimensions for multi-dimensional arrays. Then， for coherent interference， based on the model of the sub-dimension’s interference covariance matrix in tensor beamforming， a robust tensor beamforming method is proposed by directly estimating the interference covariance matrix. Analysis shows that the proposed method could overcome the non-homogeneous clutter environment and coherent interference and obtain a higher output signal-to-interference-noise ratio. Taking the given two-dimensional polarization-sensitive array as an example， the number of training snaps required by the proposed method is 1/3 of the traditional method， and the output signal-to-interference-to-noise ratio increases by about 2.5 dB under coherent interference scenario. The simulation results verify the effectiveness of the analysis.

Keywords：

robust beamforming; tensor beamforming; multi-dimensional array; coherent interference

＼=

0" 引" 言

波束成形器是天線/傳感器陣列技術(shù)與數(shù)字信號(hào)處理技術(shù)結(jié)合的產(chǎn)物，其目的在于選擇性地傳輸或接收信號(hào)，即讓期望信號(hào)無(wú)失真地通過(guò)而盡可能地抑制干擾。它作為一種經(jīng)典而傳統(tǒng)的陣列處理技術(shù)，已廣泛應(yīng)用于通信、雷達(dá)、聲納、麥克風(fēng)陣列、醫(yī)學(xué)影像等多個(gè)領(lǐng)域［13］。其實(shí)現(xiàn)過(guò)程就是：對(duì)陣列中每一個(gè)陣元接收到的信號(hào)進(jìn)行加權(quán)合成以獲得所需的期望信號(hào)。最小方差無(wú)失真響應(yīng)（minimum variance distortionless response， MVDR）算法，是目前應(yīng)用最為廣泛的波束成形算法之一。其可在固定目標(biāo)接收信號(hào)方向響應(yīng)的基礎(chǔ)上，通過(guò)最小化陣列輸出響應(yīng)的功率來(lái)獲得自適應(yīng)波束成形器的權(quán)矢量。MVDR算法最大化了陣列輸出的信干噪比（signal-to-interference-plus-noise ratio， SINR）［45］。

在MVDR波束形成器的設(shè)計(jì)中，權(quán)矢量的設(shè)計(jì)需對(duì)干擾信號(hào)與噪聲協(xié)方差矩陣（interference-plus-noise covariance matrix， INCM）求逆，而INCM在實(shí)際應(yīng)用中一般是未知的，因此就必須根據(jù)陣列輸出數(shù)據(jù)對(duì)其進(jìn)行估計(jì)。通常的方法是利用樣本協(xié)方差矩陣（sample covariance matrix， SCM），該方法也稱(chēng)為樣本矩陣求逆（sample matrix inversion， SMI）法。在陣列應(yīng)用中，如均勻線性陣列（uniform linear array， ULA），其信號(hào)的維數(shù)低，陣列輸出的樣本數(shù)，即快拍數(shù)，可以遠(yuǎn)遠(yuǎn)大于信號(hào)的維數(shù)，因此基于SCM的MVDR算法性能良好。但隨著應(yīng)用需求的增加，陣列天線已發(fā)展出了各種結(jié)構(gòu)更為復(fù)雜的結(jié)構(gòu)，如多輸入多輸出（multiple input multiple output， MIMO）陣列、包含兩個(gè)維度的面陣和T型陣列、含有3個(gè)維度的立體型陣列、能夠感知信號(hào)極化的極化敏感天線陣列等。這些陣列天線結(jié)構(gòu)相對(duì)于ULA，包含了更多的維度，此時(shí)信號(hào)的維度也隨之增加。針對(duì)這種多維輸出的數(shù)據(jù)，通常的做法是將多維數(shù)據(jù)矢量化，即將各維度天線的輸出排列成一個(gè)長(zhǎng)矢量數(shù)據(jù)。此時(shí)，由于信號(hào)的維數(shù)劇增，就導(dǎo)致樣本數(shù)難以滿足遠(yuǎn)遠(yuǎn)大于信號(hào)維數(shù)這一波束形成器設(shè)計(jì)的要求，尤其在非均勻雜波環(huán)境下，SCM的估計(jì)誤差會(huì)進(jìn)一步增大。此外，當(dāng)用于估計(jì)INCM的訓(xùn)練快拍中包含干擾信號(hào)時(shí)，干擾信號(hào)之間或者與期望信號(hào)之間不滿足獨(dú)立同分布的要求時(shí)以及陣列存在多種誤差時(shí)，均會(huì)導(dǎo)致SCM估計(jì)誤差的增加，就使得基于SCM設(shè)計(jì)的波束成形器性能急劇下降［68］。因此，研究者們提出了許多魯棒算法，例如子空間投影理論，不確定集約束，INCM重構(gòu)等［915］，這些魯棒算法提高了傳統(tǒng)MVDR算法的魯棒性，但也有其各自的局限性，例如子空間投影需要已知信源數(shù)，非確定集約束計(jì)算復(fù)雜度高等。

·" ·

同時(shí)，對(duì)于多維陣列，傳統(tǒng)波束成形方法的矢量化處理將多維陣列不同維度感知的信息完全混合在一起，不能充分表達(dá)和利用多維陣列中存在的結(jié)構(gòu)信息［1625］。在過(guò)去的幾年中，張量運(yùn)算由于其固有的多維特性在多維信號(hào)處理領(lǐng)域引起了廣泛的關(guān)注，隨著張量代數(shù)不斷完善和發(fā)展，張量分析方法具有的一些獨(dú)特優(yōu)勢(shì)，正被越來(lái)越多的研究者所發(fā)現(xiàn)。例如，張量理論越來(lái)越多地應(yīng)用于自適應(yīng)多維波束成形中：文獻(xiàn)［17］考慮到傳統(tǒng)時(shí)空自適應(yīng)處理（space time adaptive processing， STAP）算法對(duì)數(shù)據(jù)的矢量化，沒(méi)有充分利用多維數(shù)據(jù)的結(jié)構(gòu)信息;將文獻(xiàn)［16］中的張量波束成形器引入到了STAP之中，并提出了一種交替展開(kāi)高階奇異值分解（alternative unfolding high order singular valve decomposition， AU-HOSVD）算法，提升了STAP算法的性能;文獻(xiàn)［21］基于張量正則多元分解（canonical polyadic decomposition， CPD）分解，在MVDR準(zhǔn)則下，提出了一種新的波束形成器;利用張量表示的多維數(shù)據(jù)，文獻(xiàn)［22］提出了一種多維廣義旁瓣抵消器（generalized sidelobe canceller， GSC）的波束成形器，相對(duì)于傳統(tǒng)的GSC，它能實(shí)現(xiàn)更好的角度波束成形并降低計(jì)算復(fù)雜度。

一些研究者還利用張量理論，研究了多維平移不變性陣列的可分離性，然后基于子陣列級(jí)而不是全局陣列設(shè)計(jì)張量波束形成器。例如，文獻(xiàn)［23］在假設(shè)所設(shè)計(jì)的權(quán)張量具有張量外積形式的基礎(chǔ)上，基于最小均方誤差（minimum mean square error， MMSE）準(zhǔn)則與張量分解提出的張量分解MMSE （tensor decomposition MMSE， TD-MMSE）波束成形器。文獻(xiàn)［24］在文獻(xiàn)［23］的基礎(chǔ)上，將MVDR準(zhǔn)則應(yīng)用于TD-MMSE，提出了自適應(yīng)的TD-MVDR算法。文獻(xiàn)［25］則通過(guò)構(gòu)建多個(gè)虛擬陣列，對(duì)這些虛擬陣列分別進(jìn)行波束成形設(shè)計(jì)，利用張量中的克羅內(nèi)克乘積（Kronecker product， KP）合并所得到的多個(gè)波束成形器，提出了克羅內(nèi)克乘積最小方差無(wú)失真響應(yīng)（KP-MVDR）算法。

而TD-MMSE/TD-MVDR以及KP-MVDR算法均需假設(shè)所設(shè)計(jì)的權(quán)張量具有張量外積形式，然后基于多維陣列具有的多維平移不變性，討論其子陣列的可分離性或是構(gòu)建虛擬陣列。不同于此，文獻(xiàn)［18］則在STAP背景下提出了張量子波束合成（tensor sub-beam synthesis， TSS）法，并證明了多維陣列張量波束成形器設(shè)計(jì)時(shí)，其子維度的可分離性，進(jìn)一步完善并拓展了基于子維度或子陣列級(jí)設(shè)計(jì)波束成形器的理論基礎(chǔ)。與經(jīng)典的全局波束形成器相比，這些從張量理論出發(fā)，基于子陣列級(jí)或子維度級(jí)分別設(shè)計(jì)子波束成形器，然后合成多個(gè)子波束成形器的方法，在獲得性能提升的同時(shí)，也使得計(jì)算復(fù)雜度大大降低。

上述基于張量理論的波束成形方法，與傳統(tǒng)的MVDR算法相比，具有很好的性能優(yōu)勢(shì)，但其對(duì)INCM的估計(jì)與傳統(tǒng)的MVDR算法一樣，仍是基于SCM方法。因此，在實(shí)際應(yīng)用中，會(huì)由于各種誤差的存在導(dǎo)致模型適配，進(jìn)而使得波束成形器性能下降。針對(duì)該局限性，本文將進(jìn)一步完善在文獻(xiàn)［18］中提出的TSS法并基于此提出一種魯棒的張量波束合成法，以提高在實(shí)際應(yīng)用中的性能。首先，本文在給出的泛化多維陣列張量信號(hào)模型基礎(chǔ)上提出TSS法，進(jìn)一步擴(kuò)大TSS法的應(yīng)用范圍，同時(shí)給出TSS算法針對(duì)非均勻雜波環(huán)境時(shí)的魯棒性分析。其次，針對(duì)傳統(tǒng)SCM估計(jì)INCM方法存在的問(wèn)題，通過(guò)對(duì)TSS算法中子維度干擾協(xié)方差矩陣的分析，直接對(duì)INCM進(jìn)行估計(jì)，從而提出魯棒的TSS算法。

本文首先介紹了矢量MVDR算法，然后給出和分析了張量波束成形算法以及魯棒張量波束成形方法。接著給出了算法復(fù)雜度分析，最后給出了仿真實(shí)驗(yàn)結(jié)果。

1" 矢量波束成形算法

給定一個(gè)天線陣元數(shù)為N的ULA，在一個(gè)觀測(cè)周期內(nèi)有K+1個(gè)入射信號(hào)，共有P個(gè)觀測(cè)快拍，對(duì)于第P∈{1，2，…，P}個(gè)快拍x（j）∈CN，其可描述［10］為

x（P）=s0（P）a（θ0）+∑Ki=1si（P）a（θi）+n（P）（1）

式中：a（θ）=［1，ej2πfs，…，ej2π（N－1）fs］T為信號(hào)的導(dǎo)向矢量;fs=dsin（θ）/λ為空間頻率，λ為信號(hào)波長(zhǎng)，d為天線陣元間距，θ∈（－π，π）為入射信號(hào)的波達(dá)角;s0（P）為期望信號(hào)，{si（P）}Ki=1為干擾信號(hào);n（P）∈CN為獨(dú)立、零均值且方差為δ2n的復(fù)高斯白噪聲，即cov（n（P））=δ2nIN，IN表示維度為N的單位陣，并假設(shè)其與信號(hào)si（P）在統(tǒng)計(jì)上相互獨(dú)立。

根據(jù)式（1）描述的x（P），干擾加噪聲的協(xié)方差矩陣可表示為

Ci+n=E∑Ki=1sia（θi）+n∑Ki=1sia（θi）+nH=

∑Ki=1|si|2a（θi）aH（θi）+δ2nIN（2）

式中：E［·］代表期望運(yùn)算;H表示共軛轉(zhuǎn)置。

設(shè)待設(shè)計(jì)波束成形器的權(quán)矢量為w∈CN，其在期望信號(hào)無(wú)失真通過(guò)的約束下，滿足輸出功率MVDR的設(shè)計(jì)準(zhǔn)則［4］：

minw wHCi+nw

s.t. wHa（θ0）=1（3）

對(duì)于期望方向θ0，此時(shí)MVDR波束成形器的權(quán)矢量［4］為

wMVDR=C－1i+na（θ0）aH（θ0）C－1i+na（θ0）（4）

在應(yīng)用中，由于式（2）定義的INCM Ci+n一般是未知的，需要根據(jù)陣列輸出數(shù)據(jù)對(duì)其進(jìn)行估計(jì)。如果估計(jì)得到的協(xié)方差矩陣為C^N，那么MVDR權(quán)矢量則可表示為

w^MVDR=C^－1Na（θ0）aH（θ0）C^－1Na（θ0）（5）

2" 張量波束合成法

2.1" 多維陣列信號(hào)的張量模型

為了使討論具有一般性，本文首先以一個(gè)構(gòu)想的泛化多維陣列為例，來(lái)描述其輸出數(shù)據(jù)的張量模型，即考慮一個(gè)維度為Z的多維陣列，該多維陣列可包含的維度有：空域的全部3個(gè)維度或部分維度（方位角或俯仰角）、矢量天線中三維電場(chǎng)和磁場(chǎng)的6個(gè)維度，或它們的部分維度，極化維和子陣列維等。假設(shè)在t時(shí)刻，該Z維陣列接收到了期待的目標(biāo)信號(hào)和K個(gè)干擾信號(hào)，那么理想的陣列輸出數(shù)據(jù)就構(gòu)成了一個(gè)Z維的張量數(shù)據(jù)，本文稱(chēng)之為張量快拍X（t）∈CI1×I2×…×IZ，如此其張量模型就可表示為

X（t）=S0s0（t）+∑Kk=1Sksk（t）+N（t）（6）

式中：si（t）（i=0，1，2，…，K）為目標(biāo)信號(hào)和干擾信號(hào)的復(fù)幅度;S0∈CI1×I2×…×IZ為目標(biāo)信號(hào)的導(dǎo)向張量;Sk（k=1，2，…，K）為第k個(gè)干擾信號(hào)的導(dǎo)向張量;N（t）表示噪聲張量，并假設(shè)其與信號(hào)s（t）統(tǒng)計(jì)獨(dú)立?？梢钥吹?，式（6）中的期望信號(hào)S0s0（t）是一個(gè)多維信號(hào)，張量波束成形器的目的就是設(shè)計(jì)某一權(quán)張量W∈CI1×I2×…×Iz，它可使期望張量信號(hào)無(wú)失真通過(guò)同時(shí)，盡可能地抑制式（6）中的干擾。

導(dǎo)向張量可直接描述多維陣列的多維結(jié)構(gòu)信息，其也可寫(xiě)成多維陣列各子維度導(dǎo)向矢量的外積［17］：

S=s1s2…sZ（7）

式中：表示張量外積運(yùn)算;si表示第i個(gè)子維度的導(dǎo)向矢量，在應(yīng)用中，可根據(jù)給定的陣列結(jié)構(gòu)，得到相應(yīng)的導(dǎo)向矢量，然后替換式（7）中的si。此外，若將多維陣列矢量化，則有a0=vec（S0），vec（）為依列將一個(gè)多維數(shù)組轉(zhuǎn)為矢量的運(yùn)算符。

下面將以一個(gè)由Mx×My個(gè)電磁矢量傳感器構(gòu)成的均勻平面陣列作為例子來(lái)說(shuō)明本文上述多維張量信號(hào)模型的應(yīng)用。設(shè)該陣列的矢量陣元分別沿x和y方向均勻分布，如圖1所示。

那么該陣列的一個(gè)快拍，就包含了3個(gè)維度，它們分別為極化維對(duì)應(yīng)的導(dǎo)向矢量sp（γ，η）、x方向的空間導(dǎo)向矢量sx（θ，）和y方向的空間導(dǎo)向矢量sy（θ，），它們可分別定義［26］為

sp（γ，η）=－sin cos θcos

cos cos θsin

0－sin θ

cos θcos sin

cos θsin -cos

-sin θ0

cos γ

sin γejη

sx（θ，）=［1，q1x，q2x，…，qMx－1x］T， qx=e－jπsin θcos

sy（θ，）=［1，q1y，q2y，…，qMy－1y］T， qy=e－jπsin θsin （8）

式中：γ∈［0，π/2］為極化輔助角;η∈［－π，π）為極化相位差異角;θ∈［0，π］表示俯仰角;∈［0，2π）表示方位角。此時(shí)，如果令式（7）中s1=sp（γ，η），s2=sx（θ，），s3=sy（θ，），那么再根據(jù)式（6），就可得極化敏感面陣信號(hào)的張量模型。

類(lèi)似于協(xié)方差矩陣的估計(jì)，通過(guò)有限觀測(cè)樣本，協(xié)方差張量R^∈CI1×I2×…×IZ×I1×I2×…×IZ的估計(jì)［18］可寫(xiě)為

R^=1P∑Pt=1X（t）X（t）*（9）

式中：*表示共軛運(yùn)算。

2.2" TSS法

對(duì)于一個(gè)多維陣列，文獻(xiàn)［18］中推導(dǎo)證明了多維陣列子維度的可分離性，并基于此提出了一種利用TSS來(lái)設(shè)計(jì)張量波束成形器的方法，即設(shè)計(jì)權(quán)張量W的方法，其可以描述為

WTSS=w1w2…wZ

（10）

即，期待目標(biāo)多維信號(hào)的波束成形器，可通過(guò)設(shè)計(jì)各子維度波束成形器的權(quán)矢量w1=（U10s01），w2=（U20s02），…，wZ=（UZ0s0Z），再通過(guò)外積運(yùn)算而合成的方法進(jìn)行設(shè)計(jì)，這就為多維波束成形器的設(shè)計(jì)提供了一個(gè)新方法。也就是說(shuō)，所需設(shè)計(jì)的多維波束成形器，即張量波束成形器的權(quán)張量W，可先從式（6）定義的張量快拍X（t）的各子維度出發(fā)，分別設(shè)計(jì)各子維度對(duì)應(yīng)的子波束成形器{wi，i=1，2，…，Z}，然后再通過(guò)張量的外積運(yùn)算，合成各子波束成形器，獲得權(quán)張量W，該波束成形方法稱(chēng)為T(mén)SS法。

設(shè)計(jì)第i個(gè)子維度的最優(yōu)權(quán)矢量wi［18］：

wi=C－1is0isH0iC－1is0i（11）

式中：Ci=E［［X（t）］i［X（t）］Hi］為i維度信號(hào)的協(xié)方差矩陣，而［X］i表示將張量X 沿子維度i展開(kāi)的一個(gè)矩陣?？墒窃趹?yīng)用中，它一般是不可知的，但通?？衫肞個(gè)快拍近似得到它的估計(jì)C" ^i［18］：

C^i=1P∑Pt=1［X（t）］i［X（t）］Hi（12）

就本文考慮的極化敏感面陣而言，也就是可根據(jù)式（11）分別設(shè)計(jì)極化維，以及x和y方向空域維的子波束權(quán)矢量，再通過(guò)式（10）的外積運(yùn)算，合成設(shè)計(jì)出張量波束成形器所需的權(quán)張量：

WTSS=wpwxwy（13）

文獻(xiàn)［18］從平滑去相關(guān)的角度，分析了TSS法相對(duì)于傳統(tǒng)矢量MVDR在應(yīng)對(duì)相干干擾信號(hào)時(shí)所具有的魯棒性，本文則將從應(yīng)對(duì)非均勻雜波干擾的角度，繼續(xù)對(duì)TSS算法的魯棒特性進(jìn)行分析。

對(duì)于式（5）的矢量波束成形器，其歸一化輸出功率［27］可表示為

ρ（w^MVDR）=1σ2sH（θ1）C^－1NCi+nC^－1Ns（θ1）（sH（θ1）C^－1Ns（θ1））2（14）

觀察式（14）可發(fā)現(xiàn)：當(dāng)且僅當(dāng)C^－1N=C－1i+n時(shí)，式（14）才能取得最小值ρmin=ρ（wMVDR），即達(dá)到理想MVDR波束成形器的性能。否則，ρ（w^MVDR）將會(huì)變大，這也就反映了由于非理想的協(xié)方差可能導(dǎo)致的性能損失。

文獻(xiàn)［27］的分析表明：非理想?yún)f(xié)方差的估計(jì)誤差造成的性能損失會(huì)更加嚴(yán)重，一般而言，傳統(tǒng)矢量波束成形器協(xié)方差矩陣的估計(jì)均為SCM，其可表示為

C^SCM=1P∑Pj=1［x（j）］［x（j）］H（15）

將式（15）代入式（5）中，即得到了基于樣本采樣協(xié)方差矩陣估計(jì)的MVDR實(shí)現(xiàn)w^SCM。文獻(xiàn)［21］的分析表明：當(dāng)用于估計(jì)C^SCM的快拍數(shù)P與信號(hào)維數(shù)N相近時(shí)，基于SCM的波束成形器的輸出功率將顯著大于理論最小值ρmin，即波束成形器的性能將顯著下降，其滿足如下關(guān)系：當(dāng)N，P→∞ 且c=N/P∈（0，1），

則

ρ（w^SCM）ρmin=11－c （16）

式（16）表明，即使訓(xùn)練快拍數(shù)量是信號(hào)維數(shù)的兩倍（且兩者都足夠大），也會(huì)導(dǎo)致所得到的波束成形器的總輸出功率損失3 dB。此外，當(dāng)訓(xùn)練樣本數(shù)P接近N時(shí)，損失甚至可以變得無(wú)限大

根據(jù)式（16）的描述，當(dāng)c=N/P趨向于0時(shí)，MVDR的歸一化輸出功率將達(dá)到ρmin，那么對(duì)于式（14）描述的多維陣列信號(hào)而言，當(dāng)用于估計(jì)SCM的樣本數(shù)P遠(yuǎn)大于信號(hào)維數(shù)N=I1×I2×…×IZ時(shí)，傳統(tǒng)的矢量MVDR才會(huì)達(dá)到最優(yōu)性能。

對(duì)式（10）定義的TSS而言，當(dāng)構(gòu)建權(quán)張量的各子維度權(quán)矢量{wi，i=1，2，…，Z}都達(dá)到最優(yōu)解時(shí)，TSS將達(dá)到最優(yōu)性能。而據(jù)式（11）和式（16）知：當(dāng)估計(jì)子維度協(xié)方差矩陣R^i的樣本數(shù)P遠(yuǎn)遠(yuǎn)大于子維度i的維數(shù)Ii時(shí)，wi將趨近于最優(yōu)性能。這是因?yàn)榫妥泳S度Ii而言，原多維信號(hào)的其余信號(hào)維度經(jīng)模n展開(kāi)后，將并入到回波信號(hào)中，這將使子維度中信號(hào)的維數(shù)大大降低，且各子維度在設(shè)計(jì)各子波束成形器時(shí)，都容易滿足樣本數(shù)P遠(yuǎn)遠(yuǎn)大于信號(hào)維數(shù)的這一要求，即對(duì)TSS而言，當(dāng)樣本數(shù)P遠(yuǎn)遠(yuǎn)大于{Ii}Zi=1中的最大值時(shí)，其就可達(dá)到漸近的最優(yōu)性能。

上述分析表明：MVDR要達(dá)到最優(yōu)性能的條件比TSS法更加嚴(yán)苛。這就使得在相同樣本數(shù)的條件下，或者說(shuō)是非均勻環(huán)境下，TSS波束成形器將具有更優(yōu)的性能。

2.3" 魯棒TSS法

在應(yīng)用中，INCM通常是不可知的，需利用有限快拍進(jìn)行估計(jì)，即SCM估計(jì)，如式（12）和式（15）所示，SCM估計(jì)方法除了第2.2節(jié)分析的快拍數(shù)不足或非均勻環(huán)境導(dǎo)致的性能下降外，還會(huì)由于各種非理想因素的存在導(dǎo)致模型適配，使得波束成形器性能下降，進(jìn)而導(dǎo)致基于SCM估計(jì)的算法性能明顯下降。具體分析如下：假設(shè)一個(gè)ULA接收的信源，包含一個(gè)目標(biāo)信號(hào)以及一個(gè)干擾信號(hào)，且其導(dǎo)向矢量分別為a0和a1，則第t個(gè)快拍可表示為

x（t）=s0（t）a0+s1（t）a1+n（t）（17）

那么，根據(jù)式（15），其協(xié)方差矩陣可表示為

C^=1P∑Pt=1［x（t）］［x（t）］H=C^s+C^i+n+C^cross（18）

其中

C^s=1P∑Pt=1s20（t）a0aH0

C^i+n=1P∑Pt=1s21（t）a1aH1+δ2nI

C^cross=1P∑Pt=1s0（t）sH1（t）a0aH1（19）

由式（19）可知，SCM估計(jì)得到的協(xié)方差矩陣由3部分構(gòu)成，分別為目標(biāo)信號(hào)協(xié)方差矩陣C^s，干擾加噪聲的協(xié)方差矩陣C^i+n，以及目標(biāo)信號(hào)與干擾信號(hào)的互協(xié)方差矩陣C^cross，而根據(jù)式（4）定義的最優(yōu)MVDR波束成形器，協(xié)方差矩陣應(yīng)當(dāng)僅包含C^i+n部分。那么，為了獲得最優(yōu)性能，第一種方法是移除包含目標(biāo)信號(hào)的C^s和C^cross;第二種方法則是根據(jù)式（19）中定義的C^i+n，對(duì)其直接進(jìn)行估計(jì)?？紤]第一種方法，對(duì)于C^s，可利用目標(biāo)信號(hào)的先驗(yàn)信息直接移除。而對(duì)于C^cross，若目標(biāo)信號(hào)與干擾信號(hào)之間不相關(guān)，C^cross將隨著快拍數(shù)P的增加而趨向于0。當(dāng)目標(biāo)信號(hào)與干擾信號(hào)之間出現(xiàn)相干時(shí)，即使快拍數(shù)無(wú)限多，基于SCM也無(wú)法移除C^cross。顯然，第2.2節(jié)給出的基于式（12）的TSS法也存在上述魯棒問(wèn)題，為了盡可能降低C^s和C^cross對(duì)算法性能的影響，本文期望對(duì)C^i+n直接進(jìn)行估計(jì)。

為了對(duì)C^i+n進(jìn)行估計(jì)，首先需要對(duì)子維度信號(hào)協(xié)方差矩陣的數(shù)學(xué)形式進(jìn)行推導(dǎo)，考慮第t個(gè)快拍，對(duì)第t維度經(jīng)模n展開(kāi)后可表示為

［X（t）］i=∑Kk=0［Sk］isk（t）+［N（t）］i=

∑Kk=0siksHSksk（t）+［N（t）］i（20）

式中：集合S={1，2，…，i－1，i+1，…，Z}，即為除i維度外其他維度的下標(biāo)集合，根據(jù)式（7），則sHSk=sZks（Z－1）k…s（i+1）ks（i－1）k…s1k是一個(gè)維度為D=∏s∈SIs的行矢量?？紤]將sHSksk（t）視為經(jīng)模n展開(kāi)后新的源信號(hào)gSk（t），即gSk（t）=［s1k（t），s2k（t），…，sDk（t）］包含D個(gè)源信號(hào)，則經(jīng)模n展開(kāi)后的［X（t）］i可視為D個(gè)快拍，即

［X（t）］i=∑Kk=0sikgSk（t）+［N（t）］i=

［xi1（t），xi2（t），…，xiD（t）］∈CIi×D

xid（t）=∑Kk=0siksdk（t）+nid（21）

觀察式（21）可知，［X（t）］i經(jīng)模n展開(kāi)后可得到D個(gè)快拍，但這D個(gè)快拍在i維度上包含的信息是相同的，均為sik，而其他維度的信息則并入到了源信號(hào)之中，即原來(lái)的源信號(hào)s（t）變成了gSk（t），這一過(guò)程稱(chēng)為去相關(guān)［18］，或稱(chēng)為平滑。根據(jù)式（21）、式（12）和式（19），可知待估計(jì)的i維度INCM可表示為

Ci=1D∑Dd=1∑Kk=1σ2dksiksHik+σ2niIIi（22）

式中：σ2dk=E{sdk（t）sHdk（t）}。接下來(lái)，考慮利用訓(xùn)練快拍估計(jì)出所有干擾信號(hào)的波達(dá)角（direction of arrival， DOA）位置以及對(duì)應(yīng)的功率，進(jìn)而得到INCM，對(duì)于式（22）描述的協(xié)方差矩陣，其估計(jì)可表示為

C～i=∑Kk=1σ^2ks^iks^Hik+σ^2niIIi（23）

觀察式（23），待估計(jì)量為干擾信號(hào)的功率σ^2k以及對(duì)應(yīng)的導(dǎo)向矢量s^ik和噪聲功率σ^2ni。文獻(xiàn)［14］提出導(dǎo)向矢量s^ik可通過(guò)譜峰搜索獲得，本文以Capon譜為例，其定義為

PCapon（θ）=1sHi（θ）C^－1isi（θ）（24）

通過(guò)式（24）描述的Capon譜，根據(jù)譜峰位置對(duì)應(yīng)的θ即可得到待估計(jì)量s^ik，需要特別指出的是，用式（12）估計(jì)得到的協(xié)方差矩陣C^i來(lái)計(jì)算Capon譜，對(duì)相干干擾信號(hào)具有魯棒性，即出現(xiàn)相干干擾時(shí)，矩陣C^i不會(huì)出現(xiàn)虧秩問(wèn)題，這有效避免了Capon譜在相干干擾下給出偽峰或干擾信號(hào)漏估計(jì)的問(wèn)題，這點(diǎn)相對(duì)于傳統(tǒng)矢量SCM估計(jì)方法具有很大的優(yōu)勢(shì)，這也是TSS與傳統(tǒng)Capon譜結(jié)合帶來(lái)的優(yōu)勢(shì)之一。而TSS法的解相干原理，本文在式（21）中，從模n展開(kāi)給出了簡(jiǎn)單分析，文獻(xiàn)［18］則從協(xié)方差C^i的推導(dǎo)展開(kāi)了詳細(xì)的分析，本文不再贅述。

關(guān)于估計(jì)干擾信號(hào)功率σ^2k，文獻(xiàn)［15］提出，干擾信號(hào)的功率σ^2k可通過(guò)對(duì)協(xié)方差矩陣C^i做奇異值分解近似獲得，假設(shè)λ1≥λ2…≥λK≥λK+1=λK+2…=λIi為C^i從大到小依次排列的特征值，則σ^2k（t）以及σ^2ni可表示為

σ^2k=λkIi，" k=1，2，…，K

σ^2ni=λIi（25）

針對(duì)式（25）中給出的σ^2k，文獻(xiàn)［15］指出，干擾信號(hào)功率的欠估計(jì)會(huì)導(dǎo)致波束成形器性能的顯著下降，即無(wú)法有效抑制干擾信號(hào)，而干擾信號(hào)功率的過(guò)估計(jì)則可以有效提升干擾抑制性能?；诖?，文獻(xiàn)［15］給出的σ^2k估計(jì)為

σ^2k=λ1Ii（26）

即所有干擾信號(hào)的功率均通過(guò)最大特征值來(lái)進(jìn)行估計(jì)，文獻(xiàn)［15］從輸出信雜比角度詳細(xì)分析了干擾信號(hào)功率的過(guò)估計(jì)所帶來(lái)的好處，但其并未對(duì)干擾信號(hào)功率過(guò)估計(jì)的局限性作出具體分析，即干擾信號(hào)功率的過(guò)高估計(jì)會(huì)使得噪聲抑制能力下降，下面本文對(duì)此做出具體分析。

波束成形器w的輸出信噪比［9］（signal to noise ratio， SNR）可表示為

SNR=δ21|wHa1|2wHCnw（27）

式中：δ21為期望信號(hào)的功率。噪聲協(xié)方差矩陣Cn=δ2nI，并將式（4）代入式（27）中，輸出SNR可表示為

SNR=δ21δ2n·（aH1C－1i+na1）2aH1C－1i+nC－1i+na1（28）

為了便于分析，考慮一個(gè)干擾信號(hào)，那么式（2）中估計(jì)協(xié)方差矩陣的逆可表示為

C－1i+n=1δ^2nI－δ^2iδ^2n（δ^2n+Nδ^2i）aiaHi=

V－QaiaHi（29）

式中：V=1/δ^2n=Q（N+δ^2n/δ^2i）， 將式（29）代入式（28），簡(jiǎn)化可得

SNR=δ21δ2n·（VN－QL2）2NV2－2VQL2+NQ2L2（30）

式中：L=a1aHi，且滿足LN。通過(guò)分析可得，式（30）是一個(gè)隨δ^2n/δ^2i單調(diào)遞增的函數(shù)，因此對(duì)干擾信號(hào)功率的過(guò)高估計(jì)會(huì)降低輸出SNR，使得噪聲抑制能力下降。那么，為了補(bǔ)償對(duì)噪聲的抑制能力，引入對(duì)角加載算法，通過(guò)增加δ^2n以提高輸出SNR，進(jìn)一步提高波束成形器的性能，則噪聲功率σ^2ni的估計(jì)可表示為

σ^2ni=λIi+Δni（31）

式中：Δni為對(duì)角加載因子。本文采用文獻(xiàn)［28］的加載因子選擇方法，即Δni=Std（diag（C^i）），其中Std表示標(biāo)準(zhǔn)差，diag表示一個(gè)矩陣的對(duì)角元素。

綜上所述，待估計(jì)的i維度干擾信號(hào)協(xié)方差矩陣可表示為

C～i=λ1Ii·∑Kk=1s^i（θk）s^Hi（θk）+（λIi+Δni）IIi（32）

將式（32）代入式（11）即可獲得子維度wi，最后通過(guò)式（10）得到權(quán)張量W，具體算法流程如下。

步驟 1" 根據(jù)式（12）計(jì)算各子維度的INCM；

步驟 2" 根據(jù)式（24）進(jìn)行譜峰搜索；

步驟 3" 根據(jù)式（26）和式（32）重構(gòu)INCM；

步驟 4" 將重構(gòu)的INCM代入式（11）獲得子維度權(quán)矢量wi；

步驟 5" 依據(jù)式（10）計(jì)算得到權(quán)張量W。

3" 算法復(fù)雜度分析

本節(jié)將對(duì)本文方法、文獻(xiàn)［15］以及文獻(xiàn)［17］所提算法的計(jì)算復(fù)雜度進(jìn)行分析比較。以第2.1節(jié)描述極化敏感陣列為例，本文方法的計(jì)算代價(jià)主要為計(jì)算Capon譜和特征值分解O（M3x+M3y），計(jì)算Capon譜過(guò)程計(jì)算代價(jià)主要為矩陣求逆，因此其復(fù)雜度為O（M3x+M3y），那么本文方法的計(jì)算復(fù)雜度可近似表示為O（2（M3x+M3y））。文獻(xiàn)［15］中的INCM重構(gòu)算法計(jì)算復(fù)雜度約為O（2（MxMy）3）。而文獻(xiàn)［17］中基于張量AU-HOSVD的MVDR算法，其計(jì)算復(fù)雜度約為O（MyM3x+MxM3y）。綜上所述，本文方法計(jì)算復(fù)雜度顯著低于文獻(xiàn)［15］的INCM重構(gòu)算法， 同時(shí)也比文獻(xiàn)［17］的方法效率高。

4" 仿真實(shí)驗(yàn)

本節(jié)將通過(guò)仿真實(shí)驗(yàn)來(lái)對(duì)比本文方法、文獻(xiàn)［15］的INCM重構(gòu)算法、文獻(xiàn)［17］和文獻(xiàn)［18］的張量波束成形算法，以及傳統(tǒng)長(zhǎng)矢量MVDR算法的性能。主要分為3部分：第1部分是在不考慮誤差，期望信號(hào)與干擾信號(hào)不相干的條件下對(duì)比不同算法的性能;第2部分則考慮期望信號(hào)與干擾信號(hào)相干的條件下，驗(yàn)證各算法的性能;第3部分為同時(shí)考慮相干干擾和隨機(jī)的陣列視向誤差情況下，各算法性能的對(duì)比。

本文仿真實(shí)驗(yàn)條件如下：考慮圖1所示的電磁矢量均勻平面陣列，分別沿x軸均勻分布6個(gè)陣元，沿y軸均勻分布8個(gè)陣元，x軸與y軸交點(diǎn)為共用的參考陣元，陣元間距為1/2波長(zhǎng)，此時(shí)該陣列的輸出快拍包含3個(gè)維度：電磁矢量傳感器具有的極化維以及陣元維的兩個(gè)空間維，分別用1、2和3表示這3個(gè)維度，則3個(gè)維度的維數(shù)分別為I1=6，I2=6，I3=8。給定的期望信源為ej2πφ1（t），干擾源為ej2πφ2（t），相位φ1（t）與φ2（t）是［0，2π］內(nèi)的隨機(jī)值，期望信號(hào)與干擾信號(hào)的方位角均為40°，極化參數(shù)相同，俯仰角不同，期望信號(hào)的俯仰角為40°，干擾信號(hào)的俯仰角為-45°，SNR與干噪比（interference to noise ratio， INR）一致，除SNR變化的實(shí)驗(yàn)外，均為20 dB，除快拍數(shù)變化的實(shí)驗(yàn)外，快拍數(shù)固定為500。

首先，為了驗(yàn)證第2.2節(jié)的分析結(jié)果，圖2給出了傳統(tǒng)矢量MVDR算法和本文波束成形器的輸出功率對(duì)比結(jié)果。如圖2所示，基于采樣實(shí)現(xiàn)的矢量MVDR算法的歸一化輸出功率ρ會(huì)隨著參數(shù)c=N/P的減小而顯著降低，與式（16）給出的結(jié)果相符，這也驗(yàn)證了當(dāng)樣本數(shù)難以滿足遠(yuǎn)遠(yuǎn)大于信號(hào)維數(shù)的要求時(shí)，矢量MVDR無(wú)法達(dá)到理想的最小化陣列輸出功率，從而導(dǎo)致性能?chē)?yán)重下降。而隨著快拍數(shù)的增加，矢量MVDR的性能得到提高。本文方法對(duì)快拍數(shù)的要求顯著低于矢量MVDR的要求，其在圖2中的表現(xiàn)是：較低快拍數(shù)下也可達(dá)到比較理想的性能，這也驗(yàn)證了第2.2節(jié)分析得到的結(jié)論，即矢量MVDR要達(dá)到最優(yōu)性能的條件比本文提出的TSS法要更加嚴(yán)苛，在相同樣本數(shù)下，TSS法具有更優(yōu)的性能。

圖3給出了在不考慮陣列誤差時(shí)，信號(hào)非相干條件（即φ1（t）≠φ2（t））下各算法的輸出SINR與SNR的關(guān)系，可以看出，在SNR較低時(shí)（小于0 dB），除文獻(xiàn)［17］方法外，其余算法性能相似，這是因?yàn)槲墨I(xiàn)［17］方法是子空間投影，而SNR較低時(shí)，噪聲子空間與信號(hào)子空間可區(qū)分度較差，導(dǎo)致了其性能變差。而當(dāng)SNR大于0 dB小于15 dB時(shí)，本文方法與文獻(xiàn)［15］均優(yōu)于文獻(xiàn)［17］與文獻(xiàn)［18］的張量方法，說(shuō)明通過(guò)INCM重構(gòu)可有效增加對(duì)干擾信號(hào)的抑制能力，獲得更好的輸出SINR。同時(shí)，本文通過(guò)引入對(duì)角加載算法，進(jìn)一步提高對(duì)噪聲的抑制能力，使得本文方法的輸出SINR優(yōu)于文獻(xiàn)［15］方法。隨著SNR的進(jìn)一步增加，除傳統(tǒng)長(zhǎng)矢量MVDR外，其余算法的輸出SINR趨向收斂于一個(gè)值，說(shuō)明各算法的最優(yōu)性能相似，在理想高SNR條件下，均遠(yuǎn)遠(yuǎn)優(yōu)于傳統(tǒng)的MVDR算法。

圖4給出了理想條件下，各算法的輸出SINR與快拍數(shù)的關(guān)系。可以看出，本文方法、文獻(xiàn)［17］和文獻(xiàn)［18］方法收斂最快，50快拍即可收斂。而基于長(zhǎng)矢量SCM估計(jì)的文獻(xiàn)［15］和MVDR方法，在150個(gè)快拍時(shí)，其輸出SINR收斂，這說(shuō)明在非均勻雜波環(huán)境下，本文方法、文獻(xiàn)［17］和文獻(xiàn)［18］方法具有更優(yōu)的性能，進(jìn)一步驗(yàn)證了第2.2節(jié)分析得到的結(jié)論，從子維度出發(fā)的本文方法和文獻(xiàn)［17］方法，在設(shè)計(jì)各子波束成形器時(shí)，所需要的訓(xùn)練快拍數(shù)低于傳統(tǒng)的矢量方法。

圖5給出了干擾信號(hào)與期望信號(hào)相干時(shí)（即φ1（t）=aφ2（t）， a為常數(shù)）各算法的輸出SINR與SNR的關(guān)系?？梢钥闯?，文獻(xiàn)［15］方法與長(zhǎng)矢量MVDR性能相較于圖3，均出現(xiàn)了性能下降，因?yàn)槠谕盘?hào)與干擾信號(hào)相干時(shí)，傳統(tǒng)SCM估計(jì)得到的SCM中包含期望信號(hào)成分，這導(dǎo)致了基于MVDR設(shè)計(jì)的波束成形器在干擾方向和期望方向的增益相似，無(wú)法在接收期望信號(hào)的同時(shí)，有效抑制干擾信號(hào)，使得性能?chē)?yán)重下降。而對(duì)于文獻(xiàn)［15］方法，SCM中包含期望信號(hào)成分則會(huì)導(dǎo)致Capon譜估計(jì)誤差較大，無(wú)法準(zhǔn)確地重構(gòu)干擾協(xié)方差，導(dǎo)致性能?chē)?yán)重下降。而本文方法、文獻(xiàn)［17］和文獻(xiàn)［18］方法，性能與非相干干擾情況下基本相同，這是因?yàn)楦鶕?jù)式（18），設(shè)計(jì)各子波束成形器時(shí)，通過(guò)模n展開(kāi)將其他維度的信息并入源信號(hào)中，實(shí)現(xiàn)了去相關(guān)。而文獻(xiàn)［17］則是通過(guò)對(duì)協(xié)方差張量的模n展開(kāi)實(shí)現(xiàn)了去相關(guān)，因此本文方法、文獻(xiàn)［17］和文獻(xiàn)［18］方法對(duì)于相干干擾具有很好的魯棒性。

圖6給出了干擾信號(hào)與期望信號(hào)相干時(shí)各算法的輸出SINR與快拍數(shù)的關(guān)系圖。相干干擾對(duì)各算法的收斂速度沒(méi)有明顯影響，本文方法在低快拍數(shù)下，依舊有著最高的輸出SINR。而在足夠訓(xùn)練快拍下，本文方法、文獻(xiàn)［17］和文獻(xiàn)［18］方法相比非相干干擾沒(méi)有明顯性能下降，再次驗(yàn)證了這些方法對(duì)于相干干擾具有很好的魯棒性。而文獻(xiàn)［15］則因?yàn)椴荒茉谙喔筛蓴_下準(zhǔn)確地重構(gòu)干擾協(xié)方差，導(dǎo)致性能?chē)?yán)重下降。同時(shí)，長(zhǎng)矢量MVDR在相干干擾下，最優(yōu)性能也有所下降。

為了進(jìn)一步驗(yàn)證各算法的魯棒性，除了干擾信號(hào)與期望信號(hào)相干外，假設(shè)陣列存在隨機(jī)的視向誤差，即待估計(jì)信號(hào)與干擾信號(hào)產(chǎn)生隨機(jī)的DOA失配，本文假設(shè)其方位角和俯仰角失配誤差范圍為4°，即期望信號(hào)的方位角與俯仰角以及干擾信號(hào)的方位角均勻分布于［36°， 44°］之間，干擾信號(hào)的俯仰角則均勻分布于［－41°， 49°］之間。

圖7給出了同時(shí)考慮相干干擾和隨機(jī)陣列視向誤差時(shí)，各算法輸出SINR與SNR的關(guān)系，可以看出，相比較圖5的結(jié)果， 本文方法和文獻(xiàn)［15］方法沒(méi)有出現(xiàn)明顯性能下降。這得益于協(xié)方差重構(gòu)中，基于Capon譜對(duì)真實(shí)DOA的重估計(jì)，因此其對(duì)于DOA失配具有魯棒性，而文獻(xiàn)［17］方法相較于圖5也沒(méi)有出現(xiàn)性能下降，這是因?yàn)槠淅昧烁髯泳S度信號(hào)數(shù)的先驗(yàn)信息，能較準(zhǔn)確地估計(jì)信號(hào)子空間，通過(guò)子空間投影實(shí)現(xiàn)了對(duì)于DOA失配的魯棒性。而基于SCM估計(jì)的文獻(xiàn)［18］與長(zhǎng)矢量MVDR方法，由于DOA失配導(dǎo)致SCM估計(jì)誤差變大，相較于圖5性能?chē)?yán)重下降。綜上，本文提出方法在相干干擾和隨機(jī)陣列視向誤差下依然具有最優(yōu)性能。圖8則比較了相干干擾和隨機(jī)陣列視向誤差下，各算法輸出SINR與快拍數(shù)的關(guān)系。如圖8所示，本文方法的輸出SINR最高，僅有利用了各子維度信號(hào)數(shù)先驗(yàn)信息的文獻(xiàn)［17］方法與本文方法性能相近，而擁有相近收斂速度的文獻(xiàn)［18］方法由于DOA失配，導(dǎo)致最終收斂的性能較差。

5" 結(jié)" 論

本文在給出多維陣列信號(hào)張量模型的基礎(chǔ)上，提出了一種新的魯棒張量波束成形方法。通過(guò)分析，所提方法有以下優(yōu)勢(shì)：

（1） 本文在給出的泛化多維陣列張量信號(hào)模型的基礎(chǔ)上，引入張量波束成形器設(shè)計(jì)方法，利用多維陣列子維度的可分離性，通過(guò)不同的張量展開(kāi)分層進(jìn)行子波束成形器的設(shè)計(jì)，實(shí)現(xiàn)了維度的降低，進(jìn)而有效降低了算法對(duì)訓(xùn)練快拍數(shù)的需求，也有效降低了其魯棒算法的計(jì)算復(fù)雜度。

（2） 本文方法在給出的子維度干擾協(xié)方差矩陣的模型基礎(chǔ)上，直接對(duì)其進(jìn)行估計(jì)，避免了傳統(tǒng)SCM估計(jì)INCM存在的局限性，實(shí)現(xiàn)了對(duì)相干干擾信號(hào)的魯棒。

（3） 通過(guò)在協(xié)方差估計(jì)中引入對(duì)角加載算法，提高對(duì)噪聲的抑制能力，進(jìn)一步提升了算法的魯棒性能。

最后以二維極化敏感陣列為例，仿真驗(yàn)證了上述結(jié)論。

參考文獻(xiàn)

［1］ KROLIK J L. The performance of matched-field beamformers with Mediterranean vertical array data［J］. IEEE Trans.on Signal Processing， 1996， 44（10）： 26052611.

［2］ GERSHMAN A B， TURCHIN V I， ZVEREV V A. Experimental results of localization of moving underwater signal by adaptive beamforming［J］. IEEE Trans.on Signal Processing， 1995， 43（10）： 22492257.

［3］ GODARA L C. Application of antenna arrays to mobile communications. II. Beamforming and direction-of-arrival considerations［J］. Proceedings of the IEEE， 1997， 85（8）： 11951245.

［4］ CAPON J. High-resolution frequency-wavenumber spectrum analysis［J］. Proceedings of the IEEE， 1969， 57（8）： 14081418.

［5］ VAN TREES H L. Optimum array processing： part IV of detection， estimation， and modulation theory［M］. New Jersey： John Wiley amp; Sons， 2004.

［6］ BOROSON D M. Sample size considerations for adaptive arrays［J］. IEEE Trans.on Aerospace and Electronic Systems， 1980， 16（4）： 446451.

［7］ GUO J Y， YANG H C， YE Z F. A novel robust adaptive beamforming algorithm based on subspace orthogonality and projection［J］. IEEE Sensors Journal， 2023， 11： 1207612083.

［8］ HUANG Y W， FU H， VOROBYOV S A， et al. Robust adaptive beamforming via worst-case SINR maximization with nonconvex uncertainty sets［J］. IEEE Trans.on Signal Processing， 2023， 71： 218232.

［9］ LIU C， ZHEN J. Diagonal loading beamforming based on aquila optimizer［J］. IEEE Access， 2023， 11： 6909169100.

［10］ MUHAMMAD M， LI M， ABBASI Q H， et al. Adaptive diago-nal loading technique to improve direction of arrival estimation accuracy for linear antenna array sensors［J］. IEEE Sensors Journal， 2022， 22（11）： 1098610994.

［11］ LI Z Z， CAO S L， ZENG W G， et al. Robust beamforming algorithm based on steering vector estimation and interference plus noise covariance matrix reconstruction［C］∥Proc.of the IEEE Signal Processing： Algorithms， Architectures， Arrangements， and Applications， 2022： 2427.

［12］ LIAO B， CHAN S C， TSUI K M. Recursive steering vector estimation and adaptive beamforming under uncertainties［J］. IEEE Trans.on Aerospace and Electronic Systems， 2013， 49（1）： 489501.

［13］ LIAO B， TSUI K M， CHAN S C. Robust beamforming with magnitude response constraints using iterative second-order cone programming［J］. IEEE Trans.on Antennas and Propagation， 2011， 59（9）： 34773482.

［14］ ZHANG P， YANG Z， LIAO G， et al. An RCB-like steering vector estimation method based on interference matrix reduction［J］. IEEE Trans.on Aerospace and Electronic Systems， 2020， 57（1）： 636646.

［15］ ZHENG Z， YANG T， WANG W Q， et al. Robust adaptive beamforming via simplified interference power estimation［J］. IEEE Trans.on Aerospace and Electronic Systems， 2019， 55（6）： 31393152.

［16］ HAARDT M， ROEMER F， DEL G G. Higher-order SVD-based subspace estimation to improve the parameter estimation accuracy in multidimensional harmonic retrieval problems［J］. IEEE Trans.on Signal Processing， 2008， 56（7）： 31983213.

［17］ BRIGUI F， BOIZARD M， GINOLHAC G， et al. New low-rank filters for MIMO-STAP based on an orthogonal tensorial decomposition［J］. IEEE Trans.on Aerospace and Electronic Systems， 2018， 54（3）： 12081220.

［18］ 畢權(quán)楊， 李旦， 張建秋. 空時(shí)自適應(yīng)處理張量波束成形器的外積合成法［J］. 航空學(xué)報(bào)， 2019， 40（10）： 158170.

BI Q Y， LI D， ZHANG J Q. An outer product synthesis approach to tensor beamformer for space-time adaptive processing［J］. Acta Aeronautica et Astronautica Sinica， 2019， 40（10）： 158170.

［19］ XIE D， PENG H P， LI L X， et al. Semi-tensor compressed sensing［J］. Digital Signal Processing， 2016， 58： 8592.

［20］ QIAN C， FU X， SIDIROPOULOS N D， et al. Tensor-based channel estimation for dual-polarized massive MIMO systems［J］. IEEE Trans.on Signal Processing， 2018， 66（24）： 63906403.

［21］ GONG X F， LIN Q H. Spatially constrained parallel factor analysis for semi-blind beamforming［C］∥Proc.of the IEEE 7th International Conference on Natural Computation， 2011， 1： 416420.

［22］ MIRANDA R K， DA COSTA J P C L， ROEMER F， et al. Generalized sidelobe cancellers for multidimensional separable arrays［C］∥Proc.of the IEEE 6th International Workshop on Computational Advances in Multi-Sensor Adaptive Processing， 2015： 193196.

［23］ RIBEIRO L N， DE ALMEIDA A L F， MOTA J C M. Tensor beamforming for multilinear translation invariant arrays［C］∥Proc.of the IEEE International Conference on Acoustics， Speech and Signal Processing， 2016： 29662970.

［24］ LIU L， XIE J， WANG L， et al. Robust tensor beamforming for polarization sensitive arrays［J］. Multidimensional Systems and Signal Processing， 2019， 30（2）： 727748.

［25］ SHARMA R， COHEN I， BENESTY J. Adaptive and hybrid Kronecker product beamforming for far-field speech signals［J］. Speech Communication， 2020， 120： 4252.

［26］ 虞翔， 李旦， 張建秋. 魯棒成形極化敏感陣列波束的方法及極化估計(jì)［J］. 航空學(xué)報(bào)， 2017， 38（6）： 234241.

YU X， LI D， ZHANG J Q. A robust beamformer with a pola-rization sensitive array and polarization estimation［J］. Acta Aeronautica et Astronautica Sinica， 2017， 38（6）： 234241.

［27］ YANG L， MCKAY M R， COUILLET R. High-dimensional MVDR beamforming： optimized solutions based on spiked random matrix models［J］. IEEE Trans.on Signal Processing， 2018， 66（7）： 19331947.

［28］ MA N， GOH J T. Efficient method to determine diagonal loading value［C］∥Proc.of the IEEE International Conference on Acoustics， Speech， and Signal Processing， 2003， 5： V341.

作者簡(jiǎn)介

畢權(quán)楊（1995—），男，博士研究生，主要研究方向?yàn)樾盘?hào)處理。

李" 旦（1982—），男，副教授，博士，主要研究方向?yàn)樾盘?hào)處理及其應(yīng)用。

張建秋（1962—），男，教授，博士研究生導(dǎo)師，博士，主要研究方向?yàn)樾盘?hào)處理及其應(yīng)用。