第47卷第3期2021年3月
北京工业大学学报
JOURNAL OF BEIJING UNIVERSITY OF TECHNOLOGY
Vol.47No.3
Mar.2021基于压缩感知理论的二维DOA估计
窦慧晶，梁霄，张文倩
(北京工业大学信息学部，北京100124)
摘要：二维波达方向(direction of arrival,DOA)估计在雷达探测、电子对抗、医学成像等领域有着广泛的应用.针对现有算法估计精度不足、计算量巨大的问题，在基于压缩感知理论的背景下提出一种二维均匀L型阵列信号的DOA估计算法.该算法首先对阵列信号的俯仰角和方位角构建空间合成角，并对空间合成角构建过完备冗余字典；再利用正交化高斯随机矩阵构造观测矩阵；最后通过改进RM-FOCUSS算法和求解三角函数的方法还原出方位角和俯仰角.理论研究表明，该方法在高信噪比、多快拍条件下比传统算法具有更高的估计精度和分辨力，且通过压缩采样降低了运算量.仿真实验验证了上述结论.
关键词：DOA估计；压缩感知；过完备冗余字典；稀疏表示；压缩采样；测量矩阵
中图分类号：TN911文献标志码：A文章编号：0254-0037(2021)03-0231-08
doi：10.11936/bjutxb2019100002
Two-dimensional DOA Estimation Bad on
Compresd Sensing Theory
DOU Huijing,LIANG Xiao,ZHANG Wenqian
(Faculty of Information Technology，Beijing University of Technology，Beijing100124，China)
Abstract：Two-dimensional direction of arrival(DOA)estimation has been widely ud in radar detection，electronic reconnaissance，medical imaging and other fields.Aiming at the problems of inadequate estimation accuracy and enormous computational load of existing algorithms，a DOA estimation algorithm for two-dimensional uniform L-shaped array signals was prented in this paper bad on compresd nsing theory.First，an over-complete redundant dictionary was established by using the space frequency of the azimuth angle and pitch angle.Then the orthogonal Gaussian ra
ndom matrix was ud to construct the measurement matrix.Finally,azimuth and elevation were restored by improving RM -FOCUSS algorithm and solving trigonometric function.The theoretical rearch shows that the propod method has higher estimation accuracy and resolution than the traditional algorithm under the conditions of high SNR and multi-snapshot，and it reduces the computational complexity by compressing sampling.The simulation results verify the effectiveness and correctness.
Key words：direction of arrival(DOA)estimation；compresd nsing(CS)；over-complete redundant dictionary；spare reprentation；compresd sampling；measurement matrix
二维波达方向(direction of arrival,DOA)估计在阵列信号处理领域有着重要的研究意义，与一维
收稿日期：2019-10-11
基金项目：国家自然科学基金资助项目(61171137)；北京市教育委员会科研发展计划资助项目(KM201210005001)
花卉养护作者简介：窦慧晶(1969—)，女，副教授，主要从事数字信号处理、信号参量估计阵列信号处理、语音信号处理方面的研究, E-mail:dhuijing@bjut.edu
噙着泪水
232北京工业大学学报2021年
DOA估计相比，该估计算法能够更精确描述目标的空间特性，因此DOA估计在二维信号领域更具实际应用价值［1-2］.二维多重信号分类(two-dimensional multiple signal classification,2-D MUSIC)算法是目前已有的二维阵列信号DOA估计算法中最为经典的估计算法之一，该算法核心思想是将传统的一维MUSIC估计算法在二维空间进行直接推广，由于该算法需要二维谱峰搜索因而导致计算量巨大，且需要各信源的中心频率已知，因此很难满足实际应用⑶.为了解决上述缺陷，有学者提出一种无须谱峰搜索的二维旋转不变子空间(two-dimensional estimating signal parameter via rotational invariance techniques,2-D ESPRIT)算法以及二维传播算子(two-dimensional propagation method,2-D PM)算法⑷.这些算法的相继问世使阵列信号的处理性能得到一定的提高，但因其在小快拍数及低信噪比情况下估计性能严重下降而无法推广到实际应用中.
在众多阵列结构中，由于L型阵列具有结构简单、实施容易、估计性能佳等优点而被广泛用于工程领域.为解决二维信号角度匹配精度不高且计算复杂的问题，文献［5］提出一种基于L型阵列的无须手动配对的二维DOA估计算法，通过引入新的合成角度计算出新的导向矢量，进而获得原信号的俯仰角和方位角.尽管该方法能够自动完成角度配对，但需要多次谱峰搜索及特征值分解导致计算复杂度过高.文献［6］提出一种新的二维DOA估计方法，该算法首先将方位角和俯仰角分别估计出来，再通过阵列输出的互相关和信号功率对2个角度进行匹配，由于需要大量的采样信号使得该方法不可有效避免
大量的数值计算.为降低运算量有学者提出利用阵列数据的协方差矩阵进行二维角度估计的算法［7-8］.文献［9］提出一种利用多相干信号对方位角和俯仰角进行配对的方法，通过利用协方差矩阵最小化构造的代价函数从而实现角度配对，该算法存在的最大弊端是在构造协方差矩阵的过程中可能会引入外界噪声，从而影响其估计性能.
阅读区观察记录压缩感知(comprehensive nsing,CS)理论的出现为现代信号处理带来一种更高效、更精确的方法，文献［10］提出基于该理论的£-SVD算法，该算法通过对接收信号进行奇异值分解(singular value decomposition,SVD)来降低算法复杂度和对噪声的敏感性，然后利用二阶锥规划的方法求解相应的优化问题，该算法在小快拍数和低信噪比时有很好的性能，并且可以直接用于相干信号［11］.该方法摆脱了传统奈奎斯特采样定理带来超大计算量的束缚.基于此，众多学者将压缩感知理论引入到DOA估计中来，从而达到降低计算量的目的.文献［12］提出一种基于协方差矩阵联合稀疏重构的降维波达方向估计算法，该算法充分利用阵列孔径，无须预先估计目标数目，参数估计性能在低信噪比及小快拍数据长度下优势明显，但在其他方面尚有改进余地.本文在基于压缩感知理论的背景下提出一种二维L 型阵列信号的DOA估计算法.该方法在高信噪比、多快拍条件下相较于传统算法具有更高的估计精度和分辨力，且具有较低的运算量.
1信号模型
本文试验采用L型均匀阵列，该模型中2个子阵互相垂直，成90。夹角，且分别位于x轴和y轴上，如图1所示.该模型的2个子阵中均含有N个阵元，且子阵X和子阵Y共用L型阵列的原点.本文规定该模型存在k个非相关的远场窄带信号(k=1, 2,…,K)，且每个阵元间的距离为半个波长，入射信号与x轴、y轴的夹角分别为琢、茁.定义入射方向角兹、准分别为信号的方位角、俯仰角.在某时刻t接收信号的数据模型为
x(t)=A x(兹，准)s(t)+n”(t)(1)
y(t)=A y(兹，准)s(t)+n(t)(2)式中:t=1,2,…,T,T表示快拍数;s(t)=［S1(t),
S2(t),…,S k(t)］T表示K个信源t时刻信号矢量；导向矢量矩阵为
A X(兹，准)=［a x(兹1，准1),a x(兹2，准2)，…,a x(兹k，准k)］
(3)
A Y(兹，准)=［a y(兹1，准1),®(兹2，准2),…,®(兹k，准k)］
(4)其中
a x(兹k，准k)=［1，浊浊2,…，浊"-1］T(5)
a y(兹k，准k)=I1,"1,"2,…，屛一1］T(6)式(5)(6)分别为子阵X、子阵Y的导向矢量，浊= e-j2^d cos(&k)co s(准»/人滋=e-j2^d cos(兹k)sin(准k)/姿
由图1空间几何位置关系可知
{cos琢=cos兹cos渍
cos茁=sin兹cos渍
(7)
将式(7)分别带入式(5)(6)得到新的导向矢量a(琢)=[1e-j2^d cos(a k)/姿e-j2"d(M-1)ccis(琢k)/A](8) a(茁)=[1e-j2"d c如茁k)/A e-j2仔d(M-1)cos(/3k)/A](9)式中：n x(t)=[n x](t),n x2(t),…,n x N(t)]T 和
第3期窦慧晶，等：基于压缩感知理论的二维DOA估计233
厲(t)=[n y1(t),n y2((t)]T均为方差等于滓n.均值为零的加性高斯白噪声，且该噪声信号与原始接收信号相互独立；姿表示信号的波长;d代表两子阵阵元之间的距离，其值为信号的半波长.
图1均匀L型阵列模型
Fig.1L-type array signal model
2算法原理
2.1压缩感知基本原理
压缩感知理论指出：信号如果满足稀疏的或者可压缩，就可以对信号进行远低于奈奎斯特的采样频率采样，并能实现对原始信号的精确重构[13-14].该理论的核心由以下3个部分组成：
1)对原始信号进行稀疏表示；
2)针对不同情况构造合适的测量矩阵；
3)对上述信号进行稀疏重构[15].
假设信号x为离散信号，长度为N，可以通过一组正交基来表示信号，即
N
X=移鬃.琢.=追琢(10)
i=1
式中：鬃=[鬃]，鬃2,…，鬃/V]表示标准正交基，琢= [的，琢2,…,a N]T表示系数矢量.对于权系数向量琢，可以定义为
琢=追x(11)假如当0<P<2且0<k时，这些系数满足
(N\丄
II琢椰p=(移I琢i I p y臆k(12)
i=1
式中：椰琢椰p表示向量琢的p范数，则表明权系数向量琢在某变换域上可进行稀疏表示，同时表明信号具有可压缩性[16].
在信号X可进行稀疏表示或可压缩的前提下，利用不同的观测矩阵对采样信号进行投影变换，即可将原始高维信号降低维度，从而得到所需的稀疏信号[17鄄18].
观测矢量为
y=椎x=椎鬃琢=专琢(13)式中：专=椎鬃，椎为Q x N维测量矩阵,NAQ且Q逸0(Klg(N/K))，通过式(13)可以实现信号的稀疏采样.
对于压缩感知矩阵专的构造，文献[19]通过数学公式推导的方式证明专在满足限制等距原理(restricte isometry principle，RIP)准则时，可选取最小I。范数法、正交匹配追踪等算法求解式(13)，进
而高精度地完成基于压缩感知的二维DOA估计过程.
2.2基于过完备冗余字典的稀疏表示
在基于压缩感知理论下的DOA估计中，为了提高算法估计精度，文献[20]提出一种对阵列的流形矩阵进行扩展的方式构造二维信号的过完备冗余字典D的方法，该方法可以使字典D中尽可能包含所有存在的空间角度.
D=[a(兹1，准1),a(兹2，准2),…,a(兹，准k)](14)式(14)是由一维DOA估计的经典算法构造字典的经典方式直接将其推广到二维DOA估计领域得到过完备冗余字典，其中a(兹，准k)称为原子.在式(14)中，由于信号维度的扩展，致使字典原子数成倍扩增，为后续处理带来大量不必要的计算.为避免一维经典算法直接推广带来计算困扰，本文采用对第1章所提出空间合成角琢和茁分别构造过完备冗余字典的估计算法，进而得到2个空间合成角的DOA估计，最后通过参数匹配得方式找出其中的关系，从而高精度地还原出二维信号的方位角和俯仰角.
由式(8)(9)可知，定义的空间合成角琢和茁已包含二维信号的方位角和俯仰角等信息，同时已经满足空域稀疏性的要求.接下来分别对琢和茁构建超完备冗余字典
观后感格式范文
D x=[a x(琢1)，a x(琢2)，-，a x(琢N)](15)
D y=[a y(01),a y(茁2)，…，a y(茁)](16)
在构建过完备冗余字典后，由于多数文献采用的空间等角度网格划分算法存在恢复精度不高的缺陷，本文进行一定改进，利用一种空间等余弦网络划分的方法构建过完备冗余字典，及通过取cos琢=-1+(n-1)/N,n=1,2,…，2N+1，将式(1)(2)转换为稀疏表示形式
{X=D X h x+n x
b=D y h y+ U y
(17)式中h x和h y分别代表空间合成角琢和茁的系数
矢量.
234北京工业大学学报2021年
2.3测量矩阵的构造
如何设计测量矩阵椎是信号稀疏表示后的最重要一步.压缩感知理论指出：在信号可以进行稀疏表示的情况下，通常利用一个远小于信号维度的观测矩阵对该信号进行压缩采样，从而利用远低于奈奎斯特采样的方式实现信号重构，但测量矩阵的构造需满足RIP准则.
高斯随机观测矩阵是满足RIP准则最常用测量矩阵构造方法之一，该矩阵元素在(0,1)区间随机独立正态分布性质使它几乎与任何稀疏信号都不相关•但该矩阵的非结构化性质为计算带来较高的复杂度.将高斯矩阵正交化可有效增强它的结构化的特性，为后稀疏分解过程降低计算复杂度，因此本文选择正交化的高斯矩阵作为观测矩阵.具体方案如下:
1)定义高斯测量矩阵椎的RIP参数啄为满足 公式
(1-啄)椰x椰2臆II椎x椰2臆(1+啄)椰x椰2
(18)的最小值啄，式中x为K阶稀疏信号.若啄＜1，则称测量矩阵椎满足K阶RIP.
2)若P T P=I，式中I为单位测量矩阵，则称P 为正交矩阵.设高斯随机矩阵椎为M x N维,P。和P2为N x N维的正交矩阵,P]为M x M的正交矩阵，可以通过公式菲狄亚斯
椎1=椎P。(19)
椎2=P2椎P i(20)得到2个新的观测矩阵.
3)因为椎i=椎P0，所以有
II椎1s||2=W P°S椰2(21)又有
II S I2=椰P-1X||2=||X||2(22)其中椎满足啄的RIP，所以
(1-啄)1X112臆II椎X I2臆(1+啄)椰X I2
(23)由式(19)(20)得
(1-8k)I S I2臆II椎S I2臆(1+啄JI SII2
(24)即椎1也满足RIP准则，同理可得椎2也满足RIP 准则，因此椎1、椎2均可用于构造单位测量矩阵.
2.4稀疏表示的求解
由上述观测矩阵和过完备冗余字典可以得到观测矢量
x cs=椎D x h x+椎n x=专x h x+v x(25)
(26)
Jcs=椎D r
h y+椎n y=&r h y+V y
观测模型为
X cs=专x H x+V x(27)
Y cs=&r H r+V r(28)
式中：专x=椎D x，专-二椎D
r,v x=叽,v r==
椎n y, V x=专X N x,V r=专Y N r.
由式(20)可知，压缩采样可将原始信号数据从N伊1维降低到M伊1维，其中MWN，再利用每次快拍获得的M X1维观测数据求解得到稀疏系数，即得到相应的空间频率，而稀疏系数的求解可以转化成式(27)(28)的优化问题.大多数研究方向在于如何高精度低成本重构原始信号，考虑实际应用中信号的压缩采样过程是在模拟域进行，考虑到模拟器件的非理想性，因此建模时将压缩感知矩阵中含有噪声E的条件加入模型，并假设压缩感知矩阵中存在的噪声是和观测矩阵中存在的噪声相互独立且均为高斯白噪声，因此该问题的含噪模型可表示为
X cs=(专x+E x)H x+V x(29)
Yg=(专r+E r)H r+V『(30)式中：专x、O r为分解后的字典;E x、E r分别为与压缩感知矩阵专x、专r同维度的噪声矩阵.
针对该模型通过改进RM-FOCUSS算法进行求解.利用循环迭代的方法求解最优矩阵H(其中H为H x、H r总称),RM-FOCUSS针对不含噪声的E时优化问题可描述为
表示稀疏度的度量，酌=2A/1p l逸0,p沂[0,1].
该算法的求解过程可以简单描述为
湖南工商大学排名min II V I F+酌椰H I p,q
C(H)=II X cs-0H IF+y.J(p)(H)(32)式中
J(p)(H)
N
=移(II H(i,：)I2)p=
i=1
N
移
i=1
L p/2
(移l H(i,j)l2)
j=1
W;+1=diag(c[i])p/2(33)式中
(L1/2
移1H l(i,j)l2)
j=1
Q i+1=专+1(专1+1专H+1+A I)厂1(34)式中：专l+i=0W l+1(专为专、专r的统称)；姿逸0.
H“1=W i+1Q i+1(35)
第3期窦慧晶，等：基于压缩感知理论的二维DOA估计235
通过J(p)(H)的结果可知，矩阵H x、H y的列在l p5臆1)范数意义下产生了稀疏效果，而行元素是由2范数结合的，所以并没有产生稀疏效果.式(35)中l表示迭代次数，式(31)(32)中的参数酌表示平衡逼近误差与稀疏性的参数，对于不同的信噪比及解的稀疏度采用修正l曲线法可预先选取最优的值.通过查阅相关文献得知0.8臆p臆1.0时，针对信号稀疏性与计算复杂度综合考虑可知：该算法与原算法相比，具有计算量相对较小且估计精度更高的优势.
当字典专中含有噪声E的情况下，本文将RM-FOCUSS算法进行改进，提出一种通过循环迭代从而求解最优矩阵h的方法，其代价函数可表示为
C r(H)=II[V,E]II?+刃(p)(H)
该算法具体过程为：首先初始化噪声E。为与字典专维数相同的零矩阵，在l次迭代时，已知E的估计值为E l，可推导出估计值H i的表达式为
H i=argmin||Y-专H||；+姿椰W-1H||；
E l+1=argmin||E II；+姿椰+||Y-(专+E)H||；
(38)对式(38)求一阶导数后置零，其最优解可表示为&+1=[Y-专H i]H T[/+H z H T]-1(39)多次重复以上
步骤，直到收敛条件可以满足II H-Hl II F/II Hl||F<；啄(本文啄=0.001)时，停止迭代过程，得到满足要求的稀疏解H=H i+1.
通过上述方法即可求出满足条件的稀疏解H x、
h y.
2.5方位角和俯仰角的配对
相较于传统DOA估计算法需要分别求解2个子阵的观测矢量，以及接收数据的协方差矩阵进行特征值分解导致特征矢量的排列顺序产生变换的弊端，本文采用的方法可有效避免特征值分解所带来的冗余计算，最后利用频率或相关性对空间角进行匹配，再根据式
黄山几点可以上山
/cos准\
琢=arctan-------
\cos兹/
0=arcsi^V cos2兹+cos2准
(40)给出的空间角定义反推出所求的方位角和俯仰角. 3仿真结果分析
3.1实验一：3种算法有效性的比较
假设存在3个非相关窄带信号的入射角度分别是(60。，80。)、(100。，100。)、(140。，120。)，设定仿真实验的快拍数等于200，阵列中阵元总个数为16个，信噪比是-5dB，且保持恒定不变.在此条件下，分别对经典2-D MUSIC算法，参考文献[11]算法及本文算法进行仿真试验•得到仿真结果如图2~4所示.
(
。)
磺
旱
S
160
140
120
100怎样克服焦虑心理
80
60
40
20
+<+
n________________________________________________________________
020406080100120140160180
方位角/(°)
2
(
。)
暑
皐
i
图2经典2-D MUSIC算法仿真结果
Traditional2-D MUSIC algorithm simulation results
图3文献[11]算法仿真结果
Fig.3Literature[11]algorithm simulation results
通过图2~4的散点图可以直观看出，在低信噪比及小快拍数下的条件下，传统2-D MUSIC算法，文献[
11]对信号的二维角度估计虽能看出俯仰角与方位角的大致位置，但估计效果与本文算法的估计效果相比而言还是有较明显的差距，本文的估计效果更佳
.
