CN104539340B - 一种基于稀疏表示和协方差拟合的稳健波达角估计方法 - Google Patents
一种基于稀疏表示和协方差拟合的稳健波达角估计方法 Download PDFInfo
- Publication number
- CN104539340B CN104539340B CN201410827592.4A CN201410827592A CN104539340B CN 104539340 B CN104539340 B CN 104539340B CN 201410827592 A CN201410827592 A CN 201410827592A CN 104539340 B CN104539340 B CN 104539340B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msup
- mover
- mtd
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Variable-Direction Aerials And Aerial Arrays (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
本发明提出了一种基于稀疏表示和协方差拟合的稳健波达角估计方法。主要包括根据天线阵列接收信号建立稀疏空间谱表示模型,并将模型误差进行参数化表示,根据协方差拟合准则建立最优化问题。因为所得问题为非凸优化,因此通过等价转换、增加参数、分步求解对问题进行转化和求解。首先不考虑模型误差,原问题可简化为凸优化问题。利用现有方法快速求解该凸优化问题,得到初始解。再迭代求解原问题,估计模型误差参数并更新初始估计。本发明能够以低复杂度获得精确的DOA估计。
Description
技术领域
本发明属于信号处理技术领域,尤其涉及达波角方向估计,可用于无源定位和目标检测等。
背景技术
波达方向(DOA)估计是利用处于空间不同位置的天线阵列接收多个不同方向的信号源发出的信号,运用现代信号处理方法快速准确的获得信号源的方向,在雷达、声纳、无线通信等领域具有重要应用价值。
传统的DOA估计技术包括基于子空间的方法,例如MUSIC方法、最大似然估计。近年来,信号稀疏恢复和稀疏表示的发展为DOA估计技术提供了新思路。信号稀疏表示的主要思路是利用精细网格覆盖待估参数的取值空间,再找出距离参数真实值最近的格点,从而实现参数估计。目前,基于稀疏表示的DOA估计方法中较常用的是L1_SVD方法。L1_SVD方法对阵列接收信号进行稀疏空间谱表示,得到含有L1范数约束的二阶锥规划(SOCP)问题,并利用奇异值分解(SVD)来降低所求问题的维数。但当方向估计精度要求较高时,这种方法仍然面临求解高维SOCP、L1范数罚因子选择等问题。为了避开这些问题,出现了基于稀疏自相关拟合的迭代估计方法(Sparse Iterative Co-variance-based Estimation,SPICE)。此类方法避免了求解SOCP问题,复杂度低。但是该方法要求信号源位于稀疏空间谱表示时所采用的网格上,其估计精确度受到了网格密度的限制。
发明内容
技术问题:本发明针对已有技术的不足,将信号源真实位置和网格格点间的误差参数化、线性化,并结合协方差拟合准则提出一种迭代的方向估计算法。此方法是现有技术的改进和扩展,能够以低复杂度获得DOA在连续取值空间上的估计值,提高了估计精度。
技术方案:为实现上述目标,本发明的主要步骤如下:
一种基于稀疏表示和协方差拟合准则的稳健波达角估计方法,包括以下步骤:
1)采用N个全向天线形成均匀线性阵列,对空间中K个窄带信号源进行方向估计。信号源可能方向集合为Ω,其上覆盖候选方向网格其中θk表示通用方向参数。根据阵列天线的输入信号建立接收信号模型:
其中ak是对应于θk的导向矢量,表示建模误差,为相对于θk的实际位置,sk(t)是未知的源信号,n(t)是噪声,T为总的信号数目。当稀疏空间谱表示采用的网格的密度足够大时ek(θk)可以线性近似为:
其中为网格方向偏差,其均匀的分布在[-lg,lg]中,2lg为网格的步长。
2)根据天线阵列的输出信号计算信号协方差矩阵R:
其中,E(·)表示求数学期望,(·)H表示共轭转置运算,其中,E=[e1,e2,...ek],A=[a1,a2,...ak],P=diag([pT,σ]T),p=[p1,p2,...p2]T。
3)基于参数化的协方差矩阵构造参数估计的最优化问题:
3a)利用步骤2)中协方差基于p,σ,∈的参数化表示得到方向估计的加权协方差拟合准则:
其中σ表示噪声功率,∈表示网格误差,ωe为∈的协方差,∈=[∈1,∈2,...∈k]T。
令并假设和R的逆都存在,并结合参数自身约束,协方差拟合准则可等价转换为如下优化问题:
pK≥0,|∈k|≤lg k=1,2,...K
σ≥0
其中lg为网格步长的1/2,N为天线数,tr(·)表示求矩阵的迹。
3b)由于上一步中建立的最优化问题为非凸优化问题,不易求解。因此对问题进行转换,引入变量C,构造如下最优化问题:
其中B=[A,E,I]H,γ=diag(∈),表示矩阵的伪逆。
3c)求解上述二次凸优化问题,并带回目标函数可得最优值为:此值恰好为3a)中的最优化问题的目标函数。所以可以通过求解如下最优化问题来得到波达角估计的最优值:
pK≥0,|∈k|≤lg,k=1,2,...K
σ≥0
4)上述3c)中的最优化问题为非凸优化问题,求解复杂度极高。因此我们将问题的求解分为两步:首先假设模型误差为0,从而将该问题转化为易于求解的凸优化问题,并利用闭式解迭代的方法求得最优解。然后以该最优解为初始点,求解模型误差参数并迭代更新初始估计。
4a)在不考虑建模误差的情况下利用闭式解迭代的方法求解p,σ,C的初始值。将初始值代入原优化问题得到如下关于∈估计的SOCP问题:
-lg≤∈k≤lg k=1,...K
其中
4b)根据上一步求解得到的∈更新。p,σ将C,∈代入原优化问题得到包含未知数p,σ的优化问题,并求解得:
根据上式更新pk,σ的值。
4c)根据上述步骤中求得的∈,p,σ可进一步更新C:
其中C=[c1,c2,...c2K+N]H
5)由求得的C值重复上述4)中的三个步骤直到最大迭代次数或者满足迭代收敛条件最终求得最优解。
有益效果:本发明采用的技术方案通过将稀疏空间谱表示误差参数化、线性化和协方差拟合准则相结合,得到稳健的基于协方差的连续稀疏迭代估计方法。由于本发明中的方案考虑到了建模误差,能够获得DOA在连续空间上的估计值,以提高估计精度。
附图说明
图1为本发明的流程图。
图2为是本发明与现有波达方向估计方法的性能对比图。
具体实施方式
参照图1,本发明的实现步骤如下:
步骤1:采用N个全向天线接收机形成均匀线性阵列,对空间中K个窄带信号源进行方向估计。信号源可能方向集合为Ω,其上覆盖候选方向网其中θk表示通用方向参数。
步骤2根据天线阵列的输出信号计算信号协方差矩阵R的理论表达式子和实际采样值
步骤3:利用协方差矩阵构造参数估计的优化问题:
利用步骤2)中协方差关于p,σ,∈的参数化表示得到方向估计的加权协方差拟合准则并结合参数自身约束建立优化问题。由于该优化问题为非凸优化,因此引入参数C构造如下最优化问题:
pK≥0,|∈k|≤lg,k=1,2,...K
σ≥0
步骤4:上述步骤3中的最优化问题为非凸优化问题,求解复杂度极高。因此我们将问题的求解分为两步:首先假设模型误差为0,从而将该问题转化为易于求解的凸优化问题,并利用闭式解迭代的方法求得最优解;然后以该最优解为初始点,求解模型误差参数并迭代更新初始估计值。
4a)在不考虑建模误差的情况下利用闭式解迭代的方法求解p,σ,C的初始值。随机产生p,σ的初始值,令
给定∈k=0,迭代求解pk,σ
4b)将4a)中求解的p,σ,C的初始值代入原优化问题,得到如下关于∈估计的SOCP问题:
-lg≤∈k≤lg k=1,...K
4c)根据上一步求得的∈更新p,σ。将C,∈代入原优化问题得到包含未知数p,σ的优化问题,并根据柯西-施瓦茨不等式只有当下式成立时取得等号:
根据上式更新pk,σ的值。
4d)根据上述步骤中求得的∈,p,σ可进一步更新C:
其中C=[c1,c2,...c2K+N]H。
5)由求得的C值重复上述4)中的三个步骤直到最大迭代次数
或者满足迭代收敛条件最终求得最优解。
仿真结果:下面结合仿真分析本发明的性能。
仿真条件与方法:
天线阵列为均匀线性阵列且天线数N为20,各天线接收机的间距d均等于半波长。源信号个数M=13且相互独立系数稳定,其方向均匀分布在-60~60度间。源信号的模型分别为其中分布在[-π,π]。噪声为均值为零方差为σ=10(-SNR/20)的高斯白噪声,SNR定义为s1(t)功率(设a1=1)和σ的比值。采样数T等于200。
2.仿真内容与结果
仿真将本专利所提算法与不考虑稀疏表示误差时方向估计的均方根误差,以及方向估计的克拉美罗界相比,Monte Carlo仿真次数为1000。仿真结果如图2所示。图中,60、180、720分别表示用3、1、0.25度为网格密度覆盖波达方向的取值空间。由图可见,不考虑格点误差时,方向估计精度受到密度的限制,而本专利算法性能逼近克拉美罗界,并且所需的网格密度大大降低。显然,当网络密度降低时覆盖波达方向取值空间的格点数大大降低,因此也大大降低了算法的求解复杂度。
Claims (2)
1.一种基于稀疏表示和协方差拟合的稳健波达角估计方法,其特征在于包括以下步骤:
步骤1.1:通过稀疏空间谱表示对由N个全向天线组成的天线阵列接收的空间中K个窄带信号建立模型,即
<mrow>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<msub>
<mi>e</mi>
<mi>k</mi>
</msub>
<mo>(</mo>
<msub>
<mover>
<mi>&theta;</mi>
<mo>~</mo>
</mover>
<mi>k</mi>
</msub>
<mo>)</mo>
<mo>)</mo>
</mrow>
<msub>
<mi>s</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>3</mn>
<mo>,</mo>
<mn>....</mn>
<mi>T</mi>
</mrow>
其中信号源可能方向集合为Ω,其上覆盖候选方向网格其中θk表示通用方向参数,ak是对应于θk的导向矢量,sk(t)是未知的源信号,n(t)是噪声,T为信号总样本数,为信号源的实际位置,表示建模误差,当稀疏空间谱表示采用的网格的密度足够大时ek(θk)可以线性近似为:
<mrow>
<msub>
<mi>e</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>&theta;</mi>
<mo>~</mo>
</mover>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>e</mi>
<mi>k</mi>
</msub>
<msub>
<mo>&Element;</mo>
<mi>k</mi>
</msub>
</mrow>
其中为网格方向偏差,其均匀的分布在[-lg,lg]中,2lg为网格的步长;
步骤1.2:根据步骤1.1中天线阵列接收信号建立基于协方差拟合准则的稀疏优化问题,具体为:
步骤1.2.1:根据步骤1.1中天线阵列的输出信号计算信号协方差矩阵R
<mrow>
<mi>R</mi>
<mo>=</mo>
<mi>E</mi>
<mo>&lsqb;</mo>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>y</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>H</mi>
</msup>
<mo>&rsqb;</mo>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<msub>
<mi>p</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<msub>
<mi>e</mi>
<mi>k</mi>
</msub>
<msub>
<mo>&Element;</mo>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<msub>
<mi>e</mi>
<mi>k</mi>
</msub>
<msub>
<mo>&Element;</mo>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>H</mi>
</msup>
<mo>+</mo>
<mi>&sigma;</mi>
<mi>I</mi>
</mrow>
其中E(·)表示求数学期望,pk表示候选方向θk上的信号功率,σ表示噪声功率,I为单位矩阵,根据加权的协方差拟合准则得到
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>p</mi>
<mo>,</mo>
<mi>&sigma;</mi>
<mo>,</mo>
<mo>&Element;</mo>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>t</mi>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mover>
<mi>R</mi>
<mo>^</mo>
</mover>
<msup>
<mi>R</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>t</mi>
<mi>r</mi>
<mrow>
<mo>(</mo>
<msup>
<mover>
<mi>R</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>R</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mo>|</mo>
<mo>|</mo>
<mo>&Element;</mo>
<mo>|</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>/</mo>
<msub>
<mi>&omega;</mi>
<mi>e</mi>
</msub>
</mrow>
其中T为信号总样本数,∈=[∈1,∈2,...∈k]T,ωe为∈的协方差,tr(·)表示矩阵的秩,||·||为Frobenius范数;
步骤1.2.2:假设和R的逆都存在,并结合参数自身约束,步骤1.2.1中的准则可等价转换为如下优化问题
<mrow>
<munder>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mrow>
<msub>
<mi>p</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<mi>&sigma;</mi>
<mo>,</mo>
<msub>
<mo>&Element;</mo>
<mi>k</mi>
</msub>
</mrow>
</munder>
<mi>t</mi>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mover>
<mi>R</mi>
<mo>^</mo>
</mover>
<msup>
<mi>R</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mo>|</mo>
<mo>|</mo>
<mo>&Element;</mo>
<mo>|</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>/</mo>
<msub>
<mi>&omega;</mi>
<mi>e</mi>
</msub>
</mrow>
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>t</mi>
<mi>r</mi>
<mrow>
<mo>(</mo>
<msup>
<mover>
<mi>R</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>R</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>N</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
pk≥0,|∈k|≤lg k=1,2,...K
σ≥0
步骤1.2.3:在步骤1.2.2的问题中,引入辅助参量C建立最优化问题
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msup>
<mi>B</mi>
<mi>H</mi>
</msup>
<mi>C</mi>
<mo>=</mo>
<msup>
<mover>
<mi>R</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中γ=diag(∈),表示矩阵的伪逆,(·)H表示共轭转置运算,P=diag([pT,σ]T),p=[p1,p2,...,pK]T;
步骤1.2.4:求解步骤1.2.3中二次凸优化问题,并带回其目标函数可得最优值为此值恰好为步骤1.2.2中的最优化问题的目标函数,所以可以通过求解下列最优化问题来得到波达角估计的最优值
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msup>
<mi>B</mi>
<mi>H</mi>
</msup>
<mi>C</mi>
<mo>=</mo>
<msup>
<mover>
<mi>R</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>,</mo>
<mi>t</mi>
<mi>r</mi>
<mrow>
<mo>(</mo>
<msup>
<mover>
<mi>R</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>R</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>N</mi>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
pk≥0,|∈k|≤lg,k=1,2,...K
σ≥0
其中B=[A,E,I]H,E=[e1,e2,...ek],A=[a1,a2,...ak];
步骤1.3:将步骤1.2.4中得到的非凸优化问题进行转化,求解转化的问题,得到信号的波达角估计,具体为:
步骤1.3.1:首先不考虑模型误差,将步骤1.2.4中得到的最优化问题简化为能够快速求解的凸优化问题,从而求得p,σ,C的初始解,即 其中C=[c1,c2,...c2K+N]H,
步骤1.3.2:然后再利用迭代方法更新初始估计和模型误差,即先固定p,σ,C,求解步骤1.2.4中的优化问题更新模型误差,再固定模型误差,求解步骤1.2.4中的优化问题,更新p,σ,C估计。
2.根据权利要求1所述的一种基于稀疏表示和协方差拟合的稳健波达角估计方法,其特征在于步骤1.3中所述的求解步骤1.2.4中建立的优化问题,具体为:
步骤2.1:首先假设模型误差为零,在不考虑模型误差的情况下利用闭式解迭代的方法求解p,σ,C的初始值,即
步骤2.2:利用求解的p,σ,C代入原优化问题,得到关于∈的SOCP问题:
<mrow>
<munder>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mrow>
<mo>&Element;</mo>
<mo>,</mo>
<mi>t</mi>
<mo>,</mo>
<mi>f</mi>
</mrow>
</munder>
<mi>f</mi>
</mrow>
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>f</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<msub>
<mi>t</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<mi>l</mi>
<mo>,</mo>
<mi>g</mi>
<mo>&GreaterEqual;</mo>
<mo>|</mo>
<mo>|</mo>
<mo>&lsqb;</mo>
<msup>
<mo>&Element;</mo>
<mi>T</mi>
</msup>
<mo>/</mo>
<msqrt>
<msub>
<mi>&omega;</mi>
<mi>e</mi>
</msub>
</msqrt>
<mo>,</mo>
<mi>f</mi>
<mo>&rsqb;</mo>
<mo>|</mo>
<mo>|</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<msub>
<mi>t</mi>
<mi>k</mi>
</msub>
<mo>&GreaterEqual;</mo>
<mfrac>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>c</mi>
<mi>k</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
<mo>|</mo>
<mo>|</mo>
<msub>
<mover>
<mi>a</mi>
<mo>~</mo>
</mover>
<mi>k</mi>
</msub>
<mo>+</mo>
<msub>
<mo>&Element;</mo>
<mi>k</mi>
</msub>
<msub>
<mover>
<mi>e</mi>
<mo>~</mo>
</mover>
<mi>k</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
</mrow>
<msqrt>
<mi>N</mi>
</msqrt>
</mfrac>
<mo>,</mo>
<mo>-</mo>
<msub>
<mi>l</mi>
<mi>g</mi>
</msub>
<mo>&le;</mo>
<msub>
<mo>&Element;</mo>
<mi>k</mi>
</msub>
<mo>&le;</mo>
<msub>
<mi>l</mi>
<mi>g</mi>
</msub>
<mo>,</mo>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mo>...</mo>
<mi>K</mi>
</mrow>
其中根据求解的∈更新p,σ:
<mrow>
<msup>
<mi>&sigma;</mi>
<mo>+</mo>
</msup>
<mo>=</mo>
<msqrt>
<mrow>
<msub>
<mi>Nc</mi>
<mi>&sigma;</mi>
</msub>
<mo>/</mo>
<mi>t</mi>
<mi>r</mi>
<mrow>
<mo>(</mo>
<msup>
<mover>
<mi>R</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
<mo>/</mo>
<mi>&rho;</mi>
<mo>,</mo>
<msubsup>
<mi>p</mi>
<mi>k</mi>
<mo>+</mo>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<msqrt>
<mi>N</mi>
</msqrt>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<msub>
<mover>
<mo>&Element;</mo>
<mo>^</mo>
</mover>
<mi>k</mi>
</msub>
<msub>
<mo>&Element;</mo>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>c</mi>
<mi>k</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
</mrow>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<msubsup>
<mo>&Element;</mo>
<mi>k</mi>
<mn>2</mn>
</msubsup>
<mo>)</mo>
<mo>|</mo>
<mo>|</mo>
<msub>
<mover>
<mi>a</mi>
<mo>~</mo>
</mover>
<mi>k</mi>
</msub>
<mo>+</mo>
<msub>
<mo>&Element;</mo>
<mi>k</mi>
</msub>
<msub>
<mover>
<mi>e</mi>
<mo>~</mo>
</mover>
<mi>k</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
<mi>&rho;</mi>
</mrow>
</mfrac>
<mo>,</mo>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mi>K</mi>
<mo>,</mo>
</mrow>
<mrow>
<msup>
<mi>&rho;</mi>
<mo>+</mo>
</msup>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<mfrac>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<msub>
<mover>
<mo>&Element;</mo>
<mo>^</mo>
</mover>
<mi>k</mi>
</msub>
<msub>
<mo>&Element;</mo>
<mi>k</mi>
</msub>
<mo>)</mo>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>c</mi>
<mi>k</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
<mo>|</mo>
<mo>|</mo>
<msub>
<mover>
<mi>a</mi>
<mo>~</mo>
</mover>
<mi>k</mi>
</msub>
<mo>+</mo>
<msub>
<mo>&Element;</mo>
<mi>k</mi>
</msub>
<msub>
<mover>
<mi>e</mi>
<mo>~</mo>
</mover>
<mi>k</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
</mrow>
<mrow>
<msqrt>
<mi>N</mi>
</msqrt>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<msup>
<msub>
<mo>&Element;</mo>
<mi>k</mi>
</msub>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>+</mo>
<mi>l</mi>
<mo>,</mo>
</mrow>
根据步骤2.2中求解的∈,σ+可进一步更新C:
<mrow>
<msubsup>
<mi>c</mi>
<mi>k</mi>
<mo>+</mo>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>p</mi>
<mi>k</mi>
<mo>+</mo>
</msubsup>
<msup>
<mover>
<mi>R</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<msup>
<mi>R</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<msub>
<mo>&Element;</mo>
<mi>k</mi>
</msub>
<msub>
<mi>e</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mo>...</mo>
<mi>K</mi>
</mrow>
<mrow>
<mo>&lsqb;</mo>
<msubsup>
<mi>c</mi>
<mrow>
<mn>2</mn>
<mi>K</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mo>+</mo>
</msubsup>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msubsup>
<mi>c</mi>
<mrow>
<mn>2</mn>
<mi>K</mi>
<mo>+</mo>
<mi>N</mi>
</mrow>
<mo>+</mo>
</msubsup>
<mo>&rsqb;</mo>
<mo>=</mo>
<msup>
<mi>&sigma;</mi>
<mo>+</mo>
</msup>
<msup>
<mover>
<mi>R</mi>
<mo>^</mo>
</mover>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<msup>
<mi>R</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
步骤2.3:重复步骤2.2直到满足收敛条件求得最优解,其中ε为给定常数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410827592.4A CN104539340B (zh) | 2014-12-26 | 2014-12-26 | 一种基于稀疏表示和协方差拟合的稳健波达角估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410827592.4A CN104539340B (zh) | 2014-12-26 | 2014-12-26 | 一种基于稀疏表示和协方差拟合的稳健波达角估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104539340A CN104539340A (zh) | 2015-04-22 |
CN104539340B true CN104539340B (zh) | 2018-03-13 |
Family
ID=52854820
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410827592.4A Active CN104539340B (zh) | 2014-12-26 | 2014-12-26 | 一种基于稀疏表示和协方差拟合的稳健波达角估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104539340B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104965188B (zh) * | 2015-06-10 | 2018-04-27 | 重庆邮电大学 | 一种阵列误差下的波达方向估计方法 |
CN106814343B (zh) * | 2017-01-22 | 2019-04-30 | 河海大学 | 一种分步估计的空域信号空间谱估计方法 |
CN108802667A (zh) * | 2018-05-25 | 2018-11-13 | 哈尔滨工程大学 | 基于广义正交匹配追踪的波达方向估计方法 |
CN109061554B (zh) * | 2018-06-26 | 2022-07-15 | 哈尔滨工程大学 | 一种基于空间离散网格动态更新的目标到达角度估计方法 |
CN109061556B (zh) * | 2018-09-13 | 2023-03-24 | 电子科技大学 | 一种基于弹性网络的稀疏迭代波达角估计方法 |
CN113050027B (zh) * | 2021-03-08 | 2023-09-19 | 浙江大学 | 一种幅相误差情况下基于稀疏重构的波达方向估计方法 |
CN115825944B (zh) * | 2022-12-22 | 2023-05-16 | 中国人民解放军军事科学院系统工程研究院 | 基于外辐射源雷达的单快拍多目标来波方向估计方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101272168A (zh) * | 2007-03-23 | 2008-09-24 | 中国科学院声学研究所 | 一种信源数估计方法及其波达方向估计方法 |
CN103245941A (zh) * | 2013-04-16 | 2013-08-14 | 哈尔滨工程大学 | 一种基于稳健最小二乘的鲁棒波束形成方法 |
CN103399292A (zh) * | 2013-07-22 | 2013-11-20 | 西安电子科技大学 | 一种基于软稀疏表示的doa估计方法 |
CN103954950A (zh) * | 2014-04-25 | 2014-07-30 | 西安电子科技大学 | 一种基于样本协方差矩阵稀疏性的波达方向估计方法 |
CN103971029A (zh) * | 2014-05-30 | 2014-08-06 | 电子科技大学 | 一种用于网格失配下的doa估计的交替迭代方法 |
CN103983944A (zh) * | 2014-05-30 | 2014-08-13 | 电子科技大学 | 基于协方差矩阵稀疏表示的远场窄带doa估计方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6798380B2 (en) * | 2003-02-05 | 2004-09-28 | University Of Florida Research Foundation, Inc. | Robust capon beamforming |
WO2014066360A1 (en) * | 2012-10-22 | 2014-05-01 | Saab-Sensis Corporation | Sensor system and method for determining target location using sparsity-based processing |
-
2014
- 2014-12-26 CN CN201410827592.4A patent/CN104539340B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101272168A (zh) * | 2007-03-23 | 2008-09-24 | 中国科学院声学研究所 | 一种信源数估计方法及其波达方向估计方法 |
CN103245941A (zh) * | 2013-04-16 | 2013-08-14 | 哈尔滨工程大学 | 一种基于稳健最小二乘的鲁棒波束形成方法 |
CN103399292A (zh) * | 2013-07-22 | 2013-11-20 | 西安电子科技大学 | 一种基于软稀疏表示的doa估计方法 |
CN103954950A (zh) * | 2014-04-25 | 2014-07-30 | 西安电子科技大学 | 一种基于样本协方差矩阵稀疏性的波达方向估计方法 |
CN103971029A (zh) * | 2014-05-30 | 2014-08-06 | 电子科技大学 | 一种用于网格失配下的doa估计的交替迭代方法 |
CN103983944A (zh) * | 2014-05-30 | 2014-08-13 | 电子科技大学 | 基于协方差矩阵稀疏表示的远场窄带doa估计方法 |
Non-Patent Citations (1)
Title |
---|
无线通信系统中的几个典型参数估计问题;蔡曙;《中国博士学位论文全文数据库 信息科技辑》;20140215;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN104539340A (zh) | 2015-04-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104539340B (zh) | 一种基于稀疏表示和协方差拟合的稳健波达角估计方法 | |
CN104749553B (zh) | 基于快速稀疏贝叶斯学习的波达方向角估计方法 | |
CN104537249B (zh) | 基于稀疏贝叶斯学习的波达方向角估计方法 | |
CN104020439B (zh) | 基于空间平滑协方差矩阵稀疏表示的波达方向角估计方法 | |
CN109143152B (zh) | 基于张量建模的极化阵列波达方向和极化参数估计方法 | |
CN107817465A (zh) | 超高斯噪声背景下的基于无网格压缩感知的doa估计方法 | |
CN104408278A (zh) | 一种基于干扰噪声协方差矩阵估计的稳健波束形成方法 | |
CN106405533B (zh) | 基于约束加权最小二乘的雷达目标联合同步与定位方法 | |
CN106569181A (zh) | 基于协方差矩阵重构稳健Capon波束形成的算法 | |
CN103605117B (zh) | 基于干扰抑制的相控阵信号失真实时校正方法 | |
CN105791186B (zh) | 一种大规模mimo系统中的稀疏低秩信道联合估计方法 | |
CN105403874B (zh) | 非均匀阵列欠定波达方向估计方法 | |
Dai et al. | Non-uniform burst-sparsity learning for massive MIMO channel estimation | |
CN102944870A (zh) | 协方差矩阵对角加载的稳健自适应波束形成方法 | |
CN106021637A (zh) | 互质阵列中基于迭代稀疏重构的doa估计方法 | |
CN104181513B (zh) | 一种雷达天线阵元位置的校正方法 | |
CN114338301B (zh) | 一种基于压缩感知的ris辅助毫米波系统的信道估计方法 | |
CN109116293A (zh) | 一种基于离格稀疏贝叶斯的波达方向估计方法 | |
WO2021164282A1 (zh) | 基于低秩矩阵重建的水声宽频散射源的定位方法 | |
CN104076332A (zh) | 一种雷达均匀线性阵列幅度和相位的估计方法 | |
CN105137409A (zh) | 基于幅相约束的目标信号稳健空时自适应处理方法 | |
CN108111208A (zh) | 一种灵活码本配置的毫米波系统角度域信道估计方法 | |
CN106772221A (zh) | 基于机翼形变拟合的共形阵列幅相误差校正方法 | |
CN108933619B (zh) | 一种大规模mimo混合预编码码本设计方法 | |
CN103378888B (zh) | 下行波束赋形方法及设备 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |