摘" 要： 在雷達(dá)系統(tǒng)中，天線陣列至關(guān)重要。非均勻陣列能減少天線陣元使用，帶來(lái)更高的自由度，但由于其陣元位置不規(guī)則排布在到達(dá)角（DOA）估計(jì)、空間頻譜計(jì)算等方面存在難題，難以滿足雷達(dá)目標(biāo)跟蹤、聲源估計(jì)等場(chǎng)景中高精度和快速DOA估計(jì)的需求。針對(duì)該問(wèn)題，文中在[?0]原子范數(shù)的基礎(chǔ)上提出基于Circle混沌映射的并行無(wú)網(wǎng)格DOA估計(jì)算法。該算法通過(guò)并行處理目標(biāo)函數(shù)的同時(shí)進(jìn)行群體尋優(yōu)，引入自適應(yīng)策略進(jìn)行模型精英學(xué)習(xí)和狀態(tài)評(píng)估，隨后增加變異算子實(shí)現(xiàn)全局最優(yōu)解尋找，從而實(shí)現(xiàn)了非均勻陣列信號(hào)快速DOA估計(jì)。仿真驗(yàn)證表明，相對(duì)于常見(jiàn)DOA估計(jì)算法，該算法可以獲得較小的均方根誤差（RMSE），同時(shí)實(shí)現(xiàn)更精確的角度估計(jì)和更少的計(jì)算復(fù)雜度。

關(guān)鍵詞： DOA估計(jì)； [?0]原子范數(shù)； 非均勻陣列； Circle混沌映射； 自適應(yīng)策略； 變異算子

中圖分類號(hào)： TN911.7?34" " " " " " " " " " " " " " 文獻(xiàn)標(biāo)識(shí)碼： A" " " " " " " " " " "文章編號(hào)： 1004?373X（2024）17?0010?09

Non?uniform array DOA estimation based on atomic norms

LONG Weijun， XU Yizhuo， ZHANG Yulu， ZHAO Qinghua

（School of Electronics and Information Engineering， Nanjing University of Information Science and Technology， Nanjing 210044， China）

Abstract： Antenna arrays are crucial in radar systems. Non?uniform arrays can reduce the use of antenna elements and bring higher degrees of freedom. Due to the irregular arrangement of array elements， non?uniform arrays have difficulties in direction of arrival （DOA） and spatial spectrum calculation， which are difficult to meet the needs of high?precision and fast DOA estimation in radar target tracking and acoustic source estimation. In view of the above， a parallel gridless DOA estimation algorithm based on Circle chaotic map is proposed on the basis of atomic norm [?0 ]. This algorithm performs group optimizing while processing the objective function problem in parallel. An adaptive strategy is introduced for model elite learning and state evaluation. And then， mutation operators are added to find out the global optimal solution， so as to realize fast DOA estimation of non?uniform array signals. Simulation experiments show that the proposed algorithm can obtain a smaller root mean square error （RMSE）， and achieve more accurate angle estimation and less computational complexity in comparison with the existing common DOA estimation algorithms.

Keywords： DOA estimation; atomic norm [?0]; non?uniform array; Circle chaotic mapping; adaptive strategy; mutation operator

0" 引" 言

隨著智能化、信息化的發(fā)展，非均勻陣列的DOA（Direction of Arrival）估計(jì)在軍事、民用航空、汽車(chē)安全和無(wú)人駕駛等場(chǎng)景中得到廣泛應(yīng)用，尤其在雷達(dá)、聲學(xué)、多用戶通信等領(lǐng)域，陣元數(shù)目較多會(huì)增加設(shè)備使用，這將對(duì)信號(hào)處理系統(tǒng)的可靠性帶來(lái)影響。非均勻陣列能減少陣元數(shù)使用，減少資源利用，通過(guò)合成更大孔徑帶來(lái)更高的自由度[1]，可以實(shí)現(xiàn)信號(hào)定位，多目標(biāo)處理和空間譜分析，對(duì)于實(shí)現(xiàn)場(chǎng)景中的目標(biāo)追蹤和干擾消除至關(guān)重要[2?4]。但是非均勻陣列的不規(guī)則排布使得信號(hào)處理變得更加復(fù)雜，可能會(huì)造成角度估計(jì)模糊、信號(hào)區(qū)分困難等難題。

早在20世紀(jì)八九十年代，多重信號(hào)分類（Multiple Signal Classification， MUSIC）[5]、旋轉(zhuǎn)不變技術(shù)信號(hào)參數(shù)估計(jì)（Estimation of Signal Parameters via Rotational Invariance Technique， ESPRIT）[6]、加權(quán)子空間擬合（Weighted Subspace Fitting， WSF）[7]、最大似然函數(shù)估計(jì)法（Maximum Likelihood Estimation， MLE）等DOA估計(jì)算法被相繼提出，這些算法充分考慮了信號(hào)在陣列空間結(jié)構(gòu)上的特性，計(jì)算較為容易，但存在局限性，它們?cè)谠O(shè)計(jì)時(shí)假設(shè)了均勻陣列的特性，非均勻陣列的不規(guī)則排列使得這些算法無(wú)法直接應(yīng)用。2010年基于壓縮感知（Compressive Sensing， CS）的DOA估計(jì)算法被提出，相比傳統(tǒng)DOA算法，可實(shí)現(xiàn)更精確的角度估計(jì)，方法更具優(yōu)勢(shì)。CS算法基于[?1]范數(shù)正則化模型，利用矩陣稀疏結(jié)構(gòu)求解，正交匹配追蹤[8]（Orthogonal Matching Pursuit， OMP）、基于特征值分解的[?1]范數(shù)[9]（[?1]?SVD）和基于陣列協(xié)方差稀疏表示的[?1]范數(shù)[10]（[?1]?SRACV）等算法都是通過(guò)稀疏結(jié)構(gòu)引入DOA估計(jì)。CS算法集中于離散網(wǎng)格結(jié)構(gòu)，通過(guò)劃分角度集實(shí)現(xiàn)角度測(cè)量，但是由于網(wǎng)格點(diǎn)數(shù)量有限，CS算法會(huì)不可避免地產(chǎn)生估計(jì)誤差[11]。無(wú)網(wǎng)格的壓縮感知算法逐漸代替離散網(wǎng)格算法，2012年基于原子范數(shù)最小化[12]（Atomic Norm Minimization， ANM）的無(wú)網(wǎng)格DOA估計(jì)算法、交替方向乘子法（ADMM）等被提出，打破了有網(wǎng)格算法的局限性，使得目標(biāo)的定位估計(jì)或頻譜估計(jì)不需要依賴于網(wǎng)格離散化，可以直接在連續(xù)域中進(jìn)行。目前無(wú)網(wǎng)格DOA估計(jì)是通過(guò)凸松弛[8?10]求解的，對(duì)于非均勻陣列，無(wú)網(wǎng)格DOA仍然可用。

無(wú)網(wǎng)格的壓縮感知算法可以解決半正定規(guī)劃問(wèn)題（SDP），半正定規(guī)劃問(wèn)題通常通過(guò)Matlab的CVX工具箱[13]進(jìn)行求解，求解較慢且精度不高。交替迭代投影[14]（Alternating Projection， AP）、改進(jìn)的模態(tài)分析中主奇異向量法[15]（Enhanced Principal?singular?vector Utilization for Modal Analysis， EPUMA）等算法逐漸應(yīng)用起來(lái)，提高了無(wú)網(wǎng)格的計(jì)算效率。非均勻陣列對(duì)角度的精確性更加敏感，2018年有學(xué)者將ANM技術(shù)推廣到非均勻陣列中，文獻(xiàn)[13]應(yīng)用陣列插值技術(shù)構(gòu)造新的托普利茲矩陣，從虛擬的均勻陣列中得到目標(biāo)的角度值，但是該算法受到計(jì)算復(fù)雜度的影響。文獻(xiàn)[16]在ANM框架的基礎(chǔ)上，使用加速近端梯度算法，減少了空間復(fù)雜度，但也使問(wèn)題求解變得更加復(fù)雜。

然而，不論是均勻陣列還是非均勻陣列，無(wú)網(wǎng)格算法大多數(shù)基于原子[?1]范數(shù)的凸松弛展開(kāi)，在凸松弛問(wèn)題一直存在的情況下，上述算法中的無(wú)網(wǎng)格算法不能很好地保證結(jié)構(gòu)稀疏性，甚至對(duì)角度精確度有很大影響。為了獲得更好的DOA估計(jì)，突破角度精確度的限制，出現(xiàn)了基于非均勻陣列的重新加權(quán)的原子范數(shù)最小化[17]和重新加權(quán)的協(xié)方差擬合準(zhǔn)則的近似[?0]原子范數(shù)[18]等算法。文獻(xiàn)[14]在參考均勻陣列的范德蒙結(jié)構(gòu)的基礎(chǔ)上，使用不規(guī)則范德蒙矩陣分解將[?0]原子范數(shù)最小化的無(wú)網(wǎng)格DOA推廣到非均勻陣列中使用，使用交替迭代投影求解非凸優(yōu)化問(wèn)題，提高角度估計(jì)的精確性。與[?1]原子范數(shù)相比，[?0]原子范數(shù)帶來(lái)了增強(qiáng)的稀疏性和計(jì)算精度，對(duì)于具有廣泛應(yīng)用場(chǎng)景的非均勻陣列來(lái)說(shuō)，研究[?0]原子范數(shù)具有更加重要的意義[18]。

綜上所述，非均勻陣列的DOA估計(jì)在計(jì)算復(fù)雜度、角度精確度等方面存在難題，在此背景下，本文主要進(jìn)行了以下工作：

1） 依據(jù)非均勻陣列進(jìn)行接收信號(hào)模型構(gòu)建，依據(jù)原子[?0]范數(shù)確定目標(biāo)問(wèn)題模型。

2） 提出一種基于Circle映射的并行無(wú)網(wǎng)格DOA估計(jì)方法，根據(jù)目標(biāo)函數(shù)對(duì)初始集合[θ，X]進(jìn)行Circle混沌映射，設(shè)置非支配粒子，通過(guò)自適應(yīng)策略和變異算子迭代實(shí)現(xiàn)非凸計(jì)算。

3） 將2）中算法與各種DOA算法對(duì)比，比較SNR不同時(shí)角度的RMSE（Root Mean Square Error）、算法復(fù)雜度，驗(yàn)證了本文算法的有效性。

1" 模型構(gòu)建

1.1" 接收信號(hào)模型

引入DOA估計(jì)之前，需要構(gòu)造接收信號(hào)的陣列模型?？紤]一個(gè)具有[N]個(gè)陣元的非均勻線性陣列，接收來(lái)自陣列遠(yuǎn)場(chǎng)中的[K]個(gè)窄帶信號(hào)源，假設(shè)信號(hào)源和陣列位于同一平面內(nèi)。信號(hào)源的角度為[θ=θ1，θ2，…，θKT]，陣列的位置位于[r=r1，r2，…，rNT]，其中每一個(gè)[rii=1，2，…，N]以半波長(zhǎng)為單位表示陣元[i]距任意原點(diǎn)的距離，陣列排布是非均勻的。同時(shí)，假設(shè)陣列陣元間滿足各向同性，元件之間沒(méi)有耦合，陣列接收信號(hào)也沒(méi)有額外的誤差。[d]是半波長(zhǎng)，即[d=λ2]，非均勻陣列示意圖如圖1所示。

對(duì)于[L]個(gè)測(cè)量快拍，接收信號(hào)[Y∈CN×L]構(gòu)建為：

[Y=ΑX+N=i=1Kaθixi+N] （1）

式中：[A=aθ1，aθ2，…，aθK∈CN×K]為陣列流形矩陣，陣列流形向量可以表示為：[aθ=e-j2πr1sinθλ，e-j2πr2sinθλ，…，e-j2πrNsinθλT]； [X=x1，x2，…，xKT∈]

[CK×L]包含來(lái)自[K]個(gè)源的信號(hào)；[N∈CN×L]為互不相關(guān)的高斯白噪聲。

1.2" 目標(biāo)問(wèn)題建模

實(shí)現(xiàn)DOA估計(jì)，需要從觀察到的接收信號(hào)模型中恢復(fù)[θk]。早期較為常用的算法是基于原子[?1]范數(shù)的壓縮感知算法，該算法在有網(wǎng)格的基礎(chǔ)上進(jìn)行劃分，它將[θ]的所有可能性取值[θ∈-π，π]離散為[M]個(gè)虛擬網(wǎng)格單元，陣列的流形矩陣變?yōu)椋?/p>

[AM=aθ，?：aθ，?=aθ?：θ∈-π，2-MMπ，…，M-2Mπ，?2=1] （2）

式（2）是離散變化的參數(shù)索引[θ，?]的原子集。由于本文僅討論一維線性陣列，因此僅討論參數(shù)[θ]。定義原子[?0]范數(shù)和原子[?1]范數(shù)分別為：

[YAM，0=infθk，ckK：Y=k=1Ka.（θk）ck，θ∈-π，2-MMπ，…，M-2MπYAM，1=infθk，ckk=1Kck：Y=k=1Ka.（θk）ck，θ∈-π，2-MMπ，…，M-2Mπ] （3）

式中：Inf表示下確界；[ck=xk2]。

由于原子[?0]范數(shù)是NP?hard問(wèn)題，難以計(jì)算。為了便于處理，一般對(duì)原子[?1]范數(shù)凸優(yōu)化求解，通過(guò)使用CVX工具箱解決SDP問(wèn)題或者是構(gòu)造最小化協(xié)方差矩陣實(shí)現(xiàn)擬合。在壓縮感知算法中，可以通過(guò)凸松弛范數(shù)最小化來(lái)恢復(fù)[θk]和[ck]。因此，目標(biāo)的DOA估計(jì)可以通過(guò)求解式（4）基于[?1]范數(shù)的凸優(yōu)化問(wèn)題來(lái)獲得。

[mincY-Ac22+χYAM，1]" （4）

式中：[A]是[N×M]的有網(wǎng)格劃分的陣列流形矩陣；[χ]是與噪聲水平有關(guān)的正則化參數(shù)。然而，求解過(guò)程中會(huì)出現(xiàn)網(wǎng)格精度不匹配的問(wèn)題，使用更加密集的網(wǎng)格劃分會(huì)違反受限等距特性[19]（Restricted Isometry Property， RIP），使得稀疏重建變得更加難以實(shí)施，因此離網(wǎng)網(wǎng)格算法逐步應(yīng)用。

但無(wú)網(wǎng)格算法同樣需要建立[?1]范數(shù)來(lái)模擬最小化問(wèn)題，目前廣泛使用的是核范數(shù)和ANM原子范數(shù)。通過(guò)研究虛擬域原子范數(shù)的性質(zhì)，等效虛擬測(cè)量向量的原子范數(shù)最小化問(wèn)題被公式化，以無(wú)網(wǎng)格方式重構(gòu)內(nèi)插虛擬陣列協(xié)方差矩陣使得能夠進(jìn)行離網(wǎng)DOA估計(jì)。[?1]范數(shù)的凸優(yōu)化問(wèn)題都通過(guò)用其凸松弛代替來(lái)解決。

[minQ，S，Y，uY-Y22+τ2TrQ+TrT（u），s.t." " S=T（u）YYHQ， S?0]" （5）

式中：[Y∈CN×L]是接收信號(hào)的測(cè)量矩陣；[T（u）∈CN×N]是第一列為[u]的托普利茲矩陣；[Q∈CL×L]；[τ]是正則化參數(shù)系數(shù)。

與上述不同，本文不從[?1]范數(shù)考慮，不進(jìn)行非凸轉(zhuǎn)凸規(guī)劃，提出的算法從原子[?0]范數(shù)出發(fā)，受離網(wǎng)非均勻子陣列合成中恢復(fù)重構(gòu)稀疏向量的啟發(fā)[20]，這里引入一種并行處理策略，旨在盡可能多地?cái)M合期望的陣列接收測(cè)量值[Y]，在滿足匹配精度的同時(shí)，最小化[Y]的非零數(shù)量，其可以表示為：

[argminY0s.t." " Y-AX22≤ε]" （6）

簡(jiǎn)化式（6），將原子范數(shù)最小化問(wèn)題簡(jiǎn)化成并行的多個(gè)優(yōu)化最小值問(wèn)題：

[minYA，0minY-AX22]" （7）

式中：[YA，0]和[Y-AX22]分別表示原子[?0]范數(shù)和陣列接收信號(hào)的測(cè)量誤差。并行計(jì)算目標(biāo)函數(shù)兩部分，在引入[?0]原子范數(shù)保證獲得信號(hào)稀疏性的同時(shí)，最小化[Y]的非零數(shù)量，實(shí)現(xiàn)更精確的DOA估計(jì)。

2" 算法描述

本文針對(duì)式（7）目標(biāo)問(wèn)題模型，設(shè)計(jì)了一種基于Circle混沌策略的非支配排序粒子[21]目標(biāo)（Nondominated Sorting Particle Target， NSP）算法。Circle映射是一種混沌映射，它通過(guò)將輸入值映射到一個(gè)新的值來(lái)隨機(jī)產(chǎn)生初始混沌序列[22?23]。

本文從原子[?0]范數(shù)的離網(wǎng)模型構(gòu)造出發(fā)，原子[?0]范數(shù)相較于[?1]原子范數(shù)，增強(qiáng)了稀疏性和計(jì)算精度，通過(guò)最小化[?0]范數(shù)，可以實(shí)現(xiàn)對(duì)信號(hào)中非零元素的選擇，從而達(dá)到信號(hào)稀疏化的目的。針對(duì)非均勻陣列DOA估計(jì)精度問(wèn)題，這里引入混沌映射生成具有高隨機(jī)性的初始序列，減少噪聲的干擾，提高搜索能力以便獲得更優(yōu)的解。對(duì)于求解速度問(wèn)題，針對(duì)多目標(biāo)參數(shù)的優(yōu)化問(wèn)題，計(jì)算Pareto非支配解，通過(guò)非支配解配對(duì)，每一步的迭代中基于天牛算法（BAS）迭代更新步長(zhǎng)，加速更新粒子位置獲得多樣的不同權(quán)衡度的集合，從而得到更多的可行性選擇。

定義種群大小，設(shè)[Nr=30]。用混沌映射生成更高隨機(jī)性[θ]集合，依據(jù)每一個(gè)粒子的[θ]集合構(gòu)造關(guān)于目標(biāo)求解模型合成矩陣，利用天牛特性[24]對(duì)粒子的速度和動(dòng)態(tài)步長(zhǎng)加權(quán)更新粒子位置，降低計(jì)算時(shí)間，借鑒非支配排序獲得多樣性，引入自適應(yīng)策略進(jìn)行模型精英學(xué)習(xí)和狀態(tài)評(píng)估，增加擾動(dòng)算子防止陷入局部最優(yōu)[22]。最后進(jìn)行存檔處理，對(duì)存檔的數(shù)據(jù)依據(jù)擁擠距離進(jìn)行剪裁，獲得最終數(shù)據(jù)進(jìn)行非擬合判別，得到需要的[θ]集，提高估計(jì)精度。算法具體步驟如下：

1） 劃分初始[θ]集合：圖2顯示了混沌策略的優(yōu)勢(shì)。初始定義[θ]的各種集合，將每一個(gè)集合放入每一個(gè)個(gè)體中，引入Circle映射來(lái)增加隨機(jī)性，使種群分布更加多樣性[13]。其方程如下：

[θi+1=θi+0.2-mod0.52πsin（2πxi），1]" （8）

式中，[θi+1]表示映射之后的位置，[θi]表示目標(biāo)原始值，[i]表示維度。

2） 計(jì)算陣列接收信號(hào)的測(cè)量誤差和[A]中合成[Y]的最小原子數(shù)：[minY-AX22]，可以由最小二乘法進(jìn)行求解，將求解得到的值[X]代入誤差函數(shù)項(xiàng)進(jìn)行誤差計(jì)算，獲得誤差向量[FG]；[YA，0]形成最小原子數(shù)向量[TG]；

3） 依據(jù)式（7）將兩個(gè)向量拼接成矩陣，[AG=TGFG]；目標(biāo)矩陣[AG]由[θ]集和輸入源信號(hào)[X]決定，[θ]集合存在于流形矩陣[A]中，合成[Y]的最小原子數(shù)取決于種群中[θ]集合的個(gè)數(shù)。

4） 選擇適配度較高的集合[θ，X，AG]為初始精英集合，并初始化存檔矩陣[AR]，[AR]矩陣通過(guò)目標(biāo)函數(shù)和[θ]的值進(jìn)行合理存取，每經(jīng)過(guò)一次迭代都會(huì)使得存檔矩陣[AR]保留可能性較大的最佳值，避免產(chǎn)生遺漏。

5） 根據(jù)[AG]矩陣計(jì)算帕累托等級(jí)（Pareto）非支配解，實(shí)現(xiàn)非支配排序。給定一個(gè)目標(biāo)矩陣[AG]，根據(jù)目標(biāo)向量對(duì)比來(lái)進(jìn)行解排序，并計(jì)算每個(gè)解的前沿號(hào)和擁擠距離。迭代開(kāi)始時(shí)在確定好Pareto非支配解集后，使用擁擠距離計(jì)算來(lái)選擇最具有更大多樣性的解集， 輸出集合解的前沿號(hào)和擁擠距離。依據(jù)目標(biāo)問(wèn)題算法模型對(duì)同層同Pareto等級(jí)中的個(gè)體進(jìn)行選擇性排序，按照個(gè)體擁擠距離（Crowding Distance）大小排序。

6） 自適應(yīng)策略[21]——精英保留策略+狀態(tài)評(píng)估。依據(jù)前沿號(hào)（每個(gè)解的front rank）和擁擠距離進(jìn)行錦標(biāo)賽選擇輸出個(gè)體的索引，對(duì)精英集合中個(gè)體指導(dǎo)更新，通過(guò)適應(yīng)度值進(jìn)行狀態(tài)評(píng)估，基本思路如下。

① 選擇參與錦標(biāo)賽的個(gè)體。根據(jù)式（9）設(shè)置錦標(biāo)賽大小，表示每次選擇時(shí)比較候選解的數(shù)量。從群體中選擇一定數(shù)量的個(gè)體參與錦標(biāo)賽。每個(gè)個(gè)體以一定概率被選為參與者。找到第一個(gè)適應(yīng)度數(shù)組的所有索引，找出對(duì)應(yīng)的所有非支配解，[Nr]代表種群大小，根據(jù)非支配解的數(shù)量選擇個(gè)體，確保至少選擇[Nr 2]個(gè)個(gè)體。如果非支配解數(shù)量大于[Nr 2]，選擇前[Nr 2]個(gè)解索引；如果非支配解數(shù)量小于[Nr 2]，就從中隨機(jī)選擇一些個(gè)體進(jìn)行填充。

[numindex適應(yīng)度值=1gt;Nr 2，任意1：Nr 2numindex適應(yīng)度值=1lt;Nr 2，隨機(jī)個(gè)體填充] （9）

② 進(jìn)行狀態(tài)估計(jì)。隨機(jī)選擇[Nr 2]個(gè)個(gè)體作為待選擇的個(gè)體，這些個(gè)體將與①中確定的[Nr 2]非支配解配對(duì)。隨機(jī)選擇某個(gè)解與非支配解配對(duì)，將非支配解的索引放在選定個(gè)體索引數(shù)組的奇數(shù)位置，將從隨機(jī)選擇的個(gè)體中獲得的索引放在選定個(gè)體索引數(shù)組的偶數(shù)位置，過(guò)程如圖3所示。

③ 確定勝者。在每一輪的錦標(biāo)賽中，每個(gè)參與個(gè)體都會(huì)進(jìn)行一次競(jìng)爭(zhēng)，最終只有一個(gè)個(gè)體勝出。勝利的個(gè)體被選為當(dāng)前錦標(biāo)賽中的勝者[gbest]。

④ 更新精英。將當(dāng)前錦標(biāo)賽的勝者與當(dāng)前精英解進(jìn)行比較，如果勝者的適應(yīng)度值優(yōu)于當(dāng)前精英解，則將勝者取代為新的精英解。

⑤ 重復(fù)步驟③、步驟④。重復(fù)進(jìn)行錦標(biāo)賽選擇和精英更新的步驟，直到滿足迭代次數(shù)為止。

7） 更新粒子位置。根據(jù)粒子的速度和步長(zhǎng)加權(quán)更新粒子位置，步長(zhǎng)更新由BAS算法[24]獲取，圖4顯示了更新粒子位置的過(guò)程，其中[d]指的是天牛兩個(gè)觸角之間的距離，虛線表示適應(yīng)度函數(shù)曲線，[δk]表示搜索步長(zhǎng)。

8） 變異算子。更新粒子的存檔archive，每一次循環(huán)迭代將會(huì)通過(guò)精英保留策略把每一代的個(gè)體精英粒子[pbest]集合存入其中，由于進(jìn)行步驟6）可能會(huì)丟失一些非支配解或者陷入局部最優(yōu)解，這時(shí)打亂解的順序，對(duì)精英集合的解進(jìn)行變異操作，根據(jù)新生成的集合計(jì)算目標(biāo)函數(shù)值，更新非支配解，同時(shí)find函數(shù)查找排除重復(fù)解，過(guò)程如圖5所示。

9） 當(dāng)存檔遠(yuǎn)遠(yuǎn)超過(guò)自身最大存檔數(shù)時(shí)，可以依據(jù)擁擠距離對(duì)數(shù)據(jù)合理剪裁，然后重復(fù)步驟得到結(jié)果更新種群和存檔。最后將存檔的數(shù)據(jù)進(jìn)行處理，存檔中保存了許多[θ]解對(duì)應(yīng)的非擬合函數(shù)，根據(jù)目標(biāo)函數(shù)的模型，將存檔中保留的可能的精英集合代入式（7）中，找到符合目標(biāo)函數(shù)的轉(zhuǎn)折點(diǎn)，依據(jù)非擬合函數(shù)的斜率法找到符合的[θ]值[25]。

在非擬合函數(shù)的情況下，選取存檔矩陣中的[AG]。本文通過(guò)斜率法實(shí)現(xiàn)轉(zhuǎn)折點(diǎn)尋找：

① 將存檔的數(shù)據(jù)寫(xiě)入，獲得非擬合數(shù)據(jù)集，這個(gè)數(shù)據(jù)集包括了[YA，0]和[Y-AX22]在不同的[θ，X]下的變化。

② 對(duì)于數(shù)據(jù)集中的每一個(gè)數(shù)據(jù)點(diǎn)[x，y]，計(jì)算當(dāng)前點(diǎn)與其前后兩個(gè)點(diǎn)之間的斜率：

[k1=yi-1-yixi-1-xik2=xi+1-xiyi+1-yi]" "（10）

③ 計(jì)算夾角：使用向量之間的夾角公式來(lái)計(jì)算兩個(gè)向量的夾角。對(duì)于數(shù)據(jù)點(diǎn)[x1，y1]和[x2，y2]，可以將它們視為向量，計(jì)算他們之間的夾角。

[θ1=arctank1θ2=arctank2]" （11）

④ 將[θ1]和[θ2]相加，得到當(dāng)前點(diǎn)的總角度，然后將其存儲(chǔ)在[?]數(shù)組中的相應(yīng)位置。

[?=θ1+θ2]" "（12）

找到[?]數(shù)組中具有最大值的索引，這將是轉(zhuǎn)折點(diǎn)的位置。將索引更新到存檔中，得到[θ]。算法整體流程如圖6所示。

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

3.1" 一維DOA估計(jì)結(jié)果

為直觀地顯示DOA估計(jì)結(jié)果，采取空間譜顯示多個(gè)快照DOA性能。通過(guò)常規(guī)波束形成（CBF）算法、改進(jìn)MUSIC算法、ANM算法、加速近端梯度下降（APG）算法以及本文算法對(duì)目標(biāo)進(jìn)行一維DOA估計(jì)。其中快拍數(shù)取500，陣列陣元位置為[0，3，4，6，8，9]，信號(hào)的方向來(lái)自[-60°，-30°，30°，60°]。當(dāng)信噪比較高時(shí)，算法的估計(jì)結(jié)果良好，各個(gè)算法無(wú)明顯差別，但當(dāng)信噪比較低時(shí)，本文算法優(yōu)于其他算法，其他算法得到的DOA估計(jì)結(jié)果與目標(biāo)真實(shí)波達(dá)方位角之間將存在明顯的偏差。

不同SNR條件各算法估計(jì)如圖7所示，角度真值表示入射的4個(gè)期望角度[-60°，-30°，30°，60°]。

3.2" 均方根誤差

在本節(jié)中將評(píng)估所提算法的基本性能。[θ]解的準(zhǔn)確性通過(guò)真實(shí)方位角[θk]和恢復(fù)方位角[θk]之間的均方根誤差（RMSE）來(lái)衡量：

[RMSE=1Kk=1Kθk-θk2]" （13）

式中[K]代表蒙特卡洛次數(shù)。在模擬實(shí)驗(yàn)中，選擇陣列進(jìn)行模擬。本文提出算法（NSP）與已提出的估計(jì)非均勻陣列DOA的MUSIC算法、最大似然估計(jì)（MLE）算法、兩步（TWO?STEP）算法、交替投影迭代法（AP?Gridless）和交替協(xié)方差（AP?COV）算法進(jìn)行比較，最后引入了Cramér?Rao下界（CRLB），CRLB提供了一個(gè)理論上的參數(shù)估計(jì)方差下限，具體表達(dá)形式如下：

[CRLB=12π2NSUsin θ，U=1Nn=1Nrn-r2] （14）

式中：[N]表示非均勻陣列的陣元個(gè)數(shù)；[S]表示線性陣列的信噪比（SNR）；[r]表示陣列位置向量；[rn]表示第[n]個(gè)陣元的位置；[r]表示[r]的平均值。

繪制均方根誤差與SNR的關(guān)系圖，如圖8所示。選取單快拍源，源的角度為-50°和40°。該實(shí)驗(yàn)通過(guò)對(duì)每個(gè)SNR值進(jìn)行1 000次蒙特卡洛模擬實(shí)驗(yàn)獲得。在每次的實(shí)驗(yàn)中，隨機(jī)設(shè)置源的角度從[-50.5°，-49.5°][ ? ][39.5°，40.5°]中均勻隨機(jī)生成。設(shè)置MUSIC和MLE兩種算法的角度搜索范圍為[-70°，70°]，搜索精度為0.01°，TWO?STEP算法首先使用常規(guī)波束形成獲得初始DOA估計(jì)，隨后基于角度偏移的閉合公式精細(xì)化布置離網(wǎng)搜索。AP?GRIDLESS和AP?COV都是在交替投影的基礎(chǔ)上直接解決原子[?1]范數(shù)最小化問(wèn)題。

結(jié)果表明：在低信噪比條件下，SNR取值為-10～10 dB時(shí)，本文算法（NSP）優(yōu)勢(shì)更加明顯。本文算法初始化混沌映射均勻設(shè)置初始范圍，減少噪聲影響，Pareto非支配求解提供了解的多樣性。其他算法也可以提供較高的分辨率，但MUSIC算法計(jì)算復(fù)雜度高，對(duì)噪聲更加敏感。MLE算法需要處理解優(yōu)化的問(wèn)題，在多信號(hào)源場(chǎng)景復(fù)雜度較高，依據(jù)信號(hào)概率處理，估計(jì)的準(zhǔn)確性可能會(huì)受到影響。TWO?STEP算法誤差最明顯，在高信噪比下遠(yuǎn)遠(yuǎn)不如其他算法。通過(guò)對(duì)不同陣元數(shù)進(jìn)行實(shí)驗(yàn)看出，本文提出的算法在低信噪比時(shí)明顯優(yōu)于其他算法，且在信噪比較強(qiáng)的環(huán)境下性能良好且接近CRLB。

3.3" 算法復(fù)雜度比較

DOA快速估計(jì)在雷達(dá)、聲學(xué)、多用戶通信等領(lǐng)域同樣十分重要，為此比較了各種算法在非均勻陣元下的時(shí)間復(fù)雜度和計(jì)算復(fù)雜度，該仿真在具有Intel[?] CoreTM i5?12500 CPU和16 GB RAM的PC上通過(guò)Matlab R2022b實(shí)現(xiàn)。

1） 時(shí)間復(fù)雜度

各算法所需計(jì)算時(shí)間如表1所示。

非均勻陣列陣元數(shù)較少時(shí)，EM?ESPRIT運(yùn)行速度最快，該算法在EM算法的基礎(chǔ)上重構(gòu)陣列協(xié)方差進(jìn)行ESPRIT算法，其次是本文算法和TWO?STEP算法，在非均勻陣列陣元數(shù)較多時(shí)，本文算法運(yùn)行速度最快。這是由于本文算法無(wú)需考慮SDP正半定編程問(wèn)題，并行計(jì)算目標(biāo)函數(shù)，也不需要像MUSIC算法和EM?ESPRIT算法遍歷所有角度，計(jì)算速度相對(duì)較快。TWO?STEP算法與MUSIC算法相比，無(wú)需進(jìn)行空間譜搜索，通過(guò)重投影技術(shù)在信號(hào)子空間直接估計(jì)DOA，隨著DOA維數(shù)增加，TWO?STEP算法計(jì)算優(yōu)勢(shì)會(huì)更加明顯。ADMM算法將問(wèn)題分解成多個(gè)子問(wèn)題進(jìn)行獨(dú)立交替求解，但是ADMM算法收斂速率與問(wèn)題求解相關(guān)，迭代次數(shù)不可知，實(shí)際中需要數(shù)百次迭代達(dá)到精度，AP?GRIDLESS和AP?COV利用[?1]范數(shù)通過(guò)交替迭代半正定規(guī)劃矩陣來(lái)求解，與[?0]范數(shù)相比增加了模型難度，循環(huán)求解次數(shù)較多。

2） 計(jì)算復(fù)雜度

信號(hào)處理中計(jì)算復(fù)雜度的比較也同樣重要，計(jì)算復(fù)雜度可以影響運(yùn)行時(shí)間長(zhǎng)短，對(duì)各種優(yōu)化問(wèn)題直觀反映。[M]代表迭代次數(shù)，[N]代表陣元數(shù)量，[L]表示角度網(wǎng)格搜索數(shù)。各算法的計(jì)算復(fù)雜度如表2所示。

通過(guò)比較表明，本文算法（NSP）相較于其他算法計(jì)算復(fù)雜度較低。本文算法依托Pareto解，非支配前沿排序最壞的情況是有[K]個(gè)前沿，每個(gè)前沿只有一個(gè)解，重復(fù)迭代[M]遍，復(fù)雜度為[O（KM）]。對(duì)于MLE算法來(lái)說(shuō)，復(fù)雜度取決于搜索角度網(wǎng)格的復(fù)雜性。MUSIC算法涉及計(jì)算協(xié)方差矩陣，其大小為[N×N]，以及計(jì)算其特征值和特征向量，這需要[M]次迭代。EM?ESPRIT算法綜合了EM算法和ESPRIT算法的特點(diǎn)，在信號(hào)子空間估計(jì)和信號(hào)參數(shù)估計(jì)步驟中使用了迭代的方法。AP?GRIDLESS實(shí)現(xiàn)了非凸轉(zhuǎn)凸優(yōu)化的解決方案，計(jì)算復(fù)雜度高，在復(fù)雜性上相較于ADMM算法有輕微的損失。

4" 結(jié)" 論

本文提出基于混沌策略的并行無(wú)網(wǎng)格[?0]原子范數(shù)算法，通過(guò)非支配排序粒子目標(biāo)（NSP）算法進(jìn)行自適應(yīng)評(píng)估實(shí)現(xiàn)陣列的DOA估計(jì)。通過(guò)原子[?0]范數(shù)并行處理目標(biāo)模型，運(yùn)用混沌策略初始化解集，采用非支配排序和非擬合函數(shù)進(jìn)行求解，無(wú)需將目標(biāo)函數(shù)轉(zhuǎn)化為凸優(yōu)化進(jìn)行CVX處理，在陣元數(shù)較多的情況下，相對(duì)于MUSIC、AP?GRIDLESS等算法減少了計(jì)算的復(fù)雜度。當(dāng)非均勻陣列陣元數(shù)為24時(shí)，計(jì)算時(shí)間比其他算法平均快2～3 s。其次，本文算法能實(shí)現(xiàn)在信噪比條件下DOA估計(jì)且優(yōu)于其他算法，當(dāng)信噪比為-10～10 dB時(shí)，本文算法的RMSE相較于其他算法更低，當(dāng)信噪比升高時(shí)本文算法也可以接近標(biāo)準(zhǔn)的CRLB。

參考文獻(xiàn)

[1] WU X H， ZHU W P， YAN J. A high?resolution DOA estimation method with a family of nonconvex penalties [J]. IEEE transactions on vehicular technology， 2018， 67（6）： 4925?4938.

[2] GU Y J， GOODMAN N A. Information?theoretic compressive sensing kernel optimization and Bayesian Cramér?Rao bound for time delay estimation [J]. IEEE transactions on signal processing， 2017， 65（17）： 4525?4537.

[3] GU Y J， LESHEM A. Robust adaptive beamforming based on interference covariance matrix reconstruction and steering vector estimation [J]. IEEE transactions on signal processing， 2012， 60（7）： 3881?3885.

[4] PAYAMI S， BALASUBRAMANYA N M， MASOUROS C， et al. Phase shifters versus switches： An energy efficiency perspective on hybrid beamforming [J]. IEEE wireless communications letters， 2019， 8（1）： 13?16.

[5] SCHMIDT R， SCHMIDT R O. Multiple emitter location and signal parameter estimation [J]. IEEE transactions on antennas amp; propagation， 1986， 34（3）： 276?280.

[6] ROY R， KAILATH T. ESPRIT?estimation of signal parameters via rotational invariance techniques [J]. IEEE transactions on acoustics， speech， and signal processing， 1989， 37（7）： 984?995.

[7] VIBERG M， OTTERSTEN B. Detection and estimation in sensor arrays using weighted subspace fitting [J]. IEEE transactions on signal processing， 1991， 39（1）： 2436?2449.

[8] HE Z Q， LIU Q H， JIN L N， et al. Low complexity method for DOA estimation using array covariance matrix sparse representation [J]. Electronics letters， 2013， 49（3）： 228?230.

[9] MALIOUTOV D， CETIN M， WILLSKY A S. A sparse signal reconstruction perspective for source localization with sensor arrays [J]. IEEE transactions on signal processing， 2005， 53（8）： 3010?3022.

[10] YIN J H， CHEN T Q. Direction?of?arrival estimation using a sparse representation of array covariance vectors [J]. IEEE transactions on signal processing， 2011， 59（9）： 4489?4493.

[11] CHI Y， SCHARF L L， PEZESHKI A， et al. Sensitivity to basis mismatch in compressed sensing [J]. IEEE transactions on signal processing， 2011， 59（5）： 2182?2195.

[12] TANG G， BHASKAR B N， SHAH P， et al. Compressed sensing off the grid [J]. IEEE transactions on information theory， 2013， 59（11）： 7465?7490.

[13] ZHOU C， GU Y， FAN X， et al. Direction?of?arrival estimation for coprime array via virtual array interpolation [J]. IEEE transactions on signal processing， 2018， 66（22）： 5956?5971.

[14] WAGNER M， PARK Y， GERSTOFT P. Gridless DOA estimation and root?MUSIC for non?uniform linear arrays [J]. IEEE transactions on signal processing， 2021， 69： 2144?2157.

[15] QIAN C， HUANG L， SIDIROPOULOS N， et al. Enhanced PUMA for direction?of?arrival estimation and its performance analysis [J]. IEEE transactions on signal processing， 2016， 64（16）： 4127?4137.

[16] WANG Y， TIAN Z. IVDST： A fast algorithm for atomic norm minimization in line spectral estimation [J]. IEEE signal processing letters， 2018， 25（11）： 1715?1719.

[17] YANG Z， XIE L H. Enhancing sparsity and resolution via reweighted atomic norm minimization [J]. IEEE transactions on signal processing， 2016， 64（4）： 995?1006.

[18] LIU Q， SO H C， GU Y T. Off?grid DOA estimation with nonconvex regularization via joint sparse representation [J]. Signal processing， 2017， 140： 171?176.

[19] CANDèS E J. The restricted isometry property and its implications for compressed sensing [J]. Comptes rendus. mathematique， 2008， 346（9/10）： 589?592.

[20] YANG S J， Lü W T， ZHANG Z P. An off?grid compressive sensing strategy for the subarray synthesis of non?uniform linear arrays [EB/OL]. [2022?01?05]. https：//arxiv.org/abs/2112.15343.

[21] XUE B， ZHANG M J， BROWNE W N. Particle swarm optimization for feature selection in classification： A multi?objective approach [J]. IEEE transactions on cybernetics， 2013， 43（6）： 1656?1671.

[22] 張?jiān)迄i，左飛，翟正軍.基于雙Logistic變參數(shù)和Chebychev混沌映射的彩色圖像密碼算法[J].西北工業(yè)大學(xué)學(xué)報(bào)，2010，28（4）：628?632.

[23] 張達(dá)敏，徐航，王依柔，等.嵌入Circle映射和逐維小孔成像反向?qū)W習(xí)的鯨魚(yú)優(yōu)化算法[J].控制與決策，2021，36（5）：1173?1180.

[24] WANG T T， YANG L. Beetle swarm optimization algorithm： Theory and application [EB/OL]. [2018?09?02]. http：//arxiv.org/abs/1808.00206.

[25] MIERSWA I， WURST M. Information preserving multi?objective feature selection for unsupervised learning [C]// Genetic and Evolutionary Computation Conference. New York： ACM， 2006： 1545?1552.