Joumal of Computer Applications
计算机应用,
2012
,
32( 8) : 2106 - 2108, 2127
ISSN 1001-9081
CODENJYIIDU
2012-08-01
http://www.joca.cn
文章编号
:1001
-9081(2012)08
-2106
-03
doi:10.3724/SP.J.1087.2012.02106
基于稀疏表示和约束优化的波达方向估计方法
郭莹,孟彩云
(沈阳工业大学信息科学与工程学院,沈阳
110870)
(
*通信作者电子邮箱
lovelygy2002@
yahoo.
com.
cn)
摘
要:对于噪声环境中信号源的波达方向
(DOA)
估计,传统的多信号分类(
MUSIC)
算法只对不相干信号有效,
且所需较多样本。针对此问题,将进行
DOA
估计的搜索范围看作冗余字典,从而待估计的
DOA
成为该冗余字典中的
某些元素,可以由冗余字典对其进行稀疏表示;其次,利用单次快才自数据,应用二阶锥
(SOC)
约束优化的方法对该稀疏
表示问题进行描述,并进而转化为标准的二阶锥形式,采用有效的优化工具
SeDuMi
末实现
DOA
的估计。仿真结果表
明,与现有的子空间方法相比,该方法只需羊拍数据即可得到较好的估计结果,且无需对信源个数有先验知识,同时
适用于相干和非非目干信号。
关键词:波达方向估计;羊次快拍;二阶锥;子空间;稀疏表示
中图分类号:
TN91
1.
23
文献标志码
:A
DOA
estimation method based on sparse representation and constrained optimization
GUO
Ying
瑜,
MENG
Cai-yun
(School
01
11
伪
mα
tion
Science and Engineering, Shenyang University
ollnd
旧
try
,
Shenyang
Li
α
oning
110870
,
Chin
α)
Abstract:
For Direction-Of-Arrival
(DOA)
estimation of signal in additive noise, the traditional Muhiple Signal
Classification (MUSI
C)
algorithm cannot process the coherent signal with fewer snapshots. The searching scope of estimated
DOA
was
considered
as
redundant dictionary. Consequently, the estimated
DOA
was
taken as some elements in the dictionary,
and could be represented sparsely
by
the dictionary. Then, this problem
was
thrown into the Second Order Cone (SOC)
constraints and an efficient estimation algorithm using a single snapshot
was
developed. This constrained problem could be
depicted as a standard
SOC
form
and be solved
by
the SeDuMi, an optimization toolbox.
Th
e simulation results show that the
proposed algorithm has a
few
advantages over the existing
s
由
space
method including one single snapshot to be needed,
no
requirement for the number of source signals, ability to work with coherent and non-coherent signals.
Key
words:
Direction-Of-Arrival (DOA) estimation; single snapshot; Second Order Cone
(SOC);
subspace; sparse
representat
lO
n
。
引言
噪声中信号源的波达方向
(Direction-Of-Arrival
,
DOA)
估
计是传感器网络和阵列信号处理领域所研究的主要问题之
一,在水下声呐、故障诊断、电子侦察等众多场合均具有广泛
应用口气基于参数化的最大似然口气加权子空间拟合[
4J
等
方法能逼近克拉美
罗下界,但却需要全局极值的多维搜索,
多目标情况下计算量较大,且初始值的选择比较敏感。旋转
不变
(Estimating
Signal Parameters via Rotational Invariance ,
ESPRIT)
[5
人多信号分类
(M
Ul
tiple
Signal Classification ,
MUSIC)
[6-7J
等特征子空间类方法是一种次优估计方法,可采
用尽量少的数据量来计算获得尽量高的分辨率,但是这类算
法对信噪比
(Signal-to-Noise
Ratio ,
SNR)
门限和采样快拍数
有较高要求。对于相关或多径信号,需预先进行解相干处理,
例如前后向平滑(
F orward -Backward Spatial Smoothing , FBSS)
算法及其多种变形
[8
-飞从而增加了计算量,而且上述方法在
选取信号或噪声子空间时,通常需要假设已知信源个数。
若将进行
DOA
估计的搜索范围看作冗余字典,则待估计
的
DOA
只是该冗余字典中的某些元素。因此对于冗余字典
而言,待估计的
DOA
可以由冗余字典进行稀疏表示,并可由
特定函数对其进行约束。所谓信号的稀疏表示就是在构造的
冗余字典中,找到信号矢量
b
的最稀疏表示,使其满足
y
=
Ab
。其中
:y
被称为观测向量
;
A E C
MxN
(M
<
的,被称为冗
余字典或冗余矩阵
[IOJ
。若要得到该式的最优解,必须根据信
号矢量
b
的稀疏性对其进行约束,否则该问题是病态的,具有
无穷个解。最直接的约束函数是
l
。范数函数(即对矢量中非
零元素的个数进行限制)
,但是
l
。范数约束是个
NP
难解问
题,更常用的约束函数是
II
范数。另外一个不可避免的问题
是在实际应用场合中的背景噪声,因此具有实际意义的稀疏
信号表示问题应该描述为:
y=Ab+z
(1)
其中
z
是噪声向量。即该问题要解决的是由带噪观测向量
y
和冗余矩阵
A
得到
b
的最稀疏表示。
基于以上稀疏表示的基本概念和理论,本文针对信号源
的
DOA
估计问题,对信号进行
II
范数约束,对噪声进行
l2
范
数约束,并将该约束优化问题变形为凸优化问题,采用二阶锥
优化理论,借助有效的
SeDuMi
工具箱进行参数估计。与
MUSIC
、
FBSS
等子空间算法相比,本文的方法具有所需采样
收稿日期
:2012-01-
06;
修回日期:
2012- 03-
04
。
基金项目:国家自然科学基金资助项目(
61102123)
;辽宁省教育厅一般项目
(口
011019);
沈阳工业大学中青年骨干教师资助项目;沈阳工业大学博士启动项目。
作者简介:郭莹(1
975
- )
,女,辽宁铁岭人,讲师,博士,主要研究方向:非高斯信号处理、参数估计;
孟彩云(1
986
- )
,女,安徽淮北人,硕
士研究生,主要研究方向:非高斯信号处理、自适应滤波。