CN106021637B - 互质阵列中基于迭代稀疏重构的doa估计方法 - Google Patents
互质阵列中基于迭代稀疏重构的doa估计方法 Download PDFInfo
- Publication number
- CN106021637B CN106021637B CN201610237414.5A CN201610237414A CN106021637B CN 106021637 B CN106021637 B CN 106021637B CN 201610237414 A CN201610237414 A CN 201610237414A CN 106021637 B CN106021637 B CN 106021637B
- Authority
- CN
- China
- Prior art keywords
- array
- relatively prime
- sparse
- function
- iteration
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Abstract
本发明公开了互质阵列中一种基于迭代稀疏重构的DOA估计方法,接收天线阵列采用非线性互质阵列,通过对接收信号的二阶统计特性协‑方差矩阵的向量化处理,可确定孔径长度更大的差分阵列,提高检测能力。对目标所在的角度域进行离散化处理,目标可视为稀疏分布于栅格点之上或附近,并构建出对数和形式的稀疏信号重建问题。利用对数和函数的凸紧上界,重建该原始稀疏问题,以迭代的方式动态调整更新角度域的离散点以逼近目标真实的到达角度。
Description
技术领域
本发明涉及通信信号处理技术领域,尤其涉及一种使用互质阵列的基于迭代稀疏重构的DOA估计方法。
背景技术
基于天线阵列的信号处理方法被广泛应用于无线通信、电磁场、雷达、声呐等诸多领域。波达方向(DOA)估计是阵列信号处理领域中的一个重要问题。
传统的估计方法通常研究阵元间隔为半波长的均匀线性阵列,且适用于检测目标数少于阵元个数的场合,例如对N天线均匀线性阵列而言,传统估计方法(如基于子空间的估计方法等)至多可检测N-1个目标,非线性互质阵列通过利用协方差矩阵的特性,将原始阵列构建成具有更多虚拟天线和更大孔径长度的差分阵列,进而可显著提高其自由度,即检测能力,传统估计方法通常需要目标个数的先验信息和大样本才能估计出目标的波达方向,在小样本及目标个数未知条件下不适用,此外,传统估计方法也难以适用于信噪比低的情况,即可能无法检测到发射功率较低(例如小于噪声功率)的弱目标。
近年来,随着稀疏重构理论的发展,稀疏表示已逐渐开始应用于小波去噪、雷达成像等领域,基于稀疏信号重构的DOA估计方法能够充分发掘互质阵列自由度高的优势,可检测目标数目较传统方法显著提高,然而,常规的稀疏重构方法需要将角度域进行静态的栅格化处理,建立以l0范数(通常以lp范数近似,p≤1)最小化为目标的稀疏优化问题,这类方法的问题在于估计精度严重依赖于初始的栅格化过程,若目标大部分位于或极其靠近栅格点,则估计性能良好;反之,若大部分目标流离于栅格点之外,则估计性能难以保证,因此,由于角度域的静态栅格化处理造成的栅格失配问题会严重影响重构的效果。
发明内容
本发明针对上述现有技术存在的不足,提出了一种适用于互质阵列的基于迭代稀疏重构的DOA估计方法,本发明的方法具有更高的自由度,可检测比实际阵元数更多的目标,本发明的方法的估计性能不依赖于初始的栅格化处理,具有更强的灵活性,适用于小样本、低信噪比及目标个数未知等复杂环境,其中所述DOA全称为Direction of Arrival,即信号的波达方向。
为实现上述发明目的,本发明所采取的技术方案为:
一种互质阵列中基于迭代稀疏重构的DOA估计方法,它包括以下步骤:
(1)建立接收阵列模型
用非线性的互质阵列作为接收阵列采集信号,其中互质阵列可分解为两个均匀线性子阵列,每个子阵列的阵元间隔均大于半波长;
(2)构建虚拟差分阵列
对互质阵列接收信号的协方差矩阵进行向量化处理,在不同阵元的差分位置形成虚拟阵元,可将原始阵列等效成具有更多虚拟阵元及更大孔径尺度的差分阵列;
(3)栅格化处理,建立稀疏优化问题
将目标区域进行栅格化处理,分析差分阵列中的接收信号,建立以对数和函数为优化目标的稀疏优化问题,兼顾估计结果的稀疏性和准确性;
(4)DOA估计迭代实现
利用对数和函数的凸紧上界重构原始稀疏问题,以迭代方式逐步更新栅格点位置来使目标函数最小化,逐渐逼近真实信号源的方向,直至满足终止条件,所述方法可动态调整栅格点的位置,能有效克服角度域的静态离散化造成的栅格失配问题。
本发明的有益效果为:
1、本发明利用互质阵列构建出具有更多阵元数目和更大孔径尺寸的虚拟差分阵列,显著提高了自由度,可实现多于阵元数目的目标检测;
2、相对于传统角度域的静态栅格化,本发明以迭代方式动态调整栅格点来逐渐逼近目标的真实位置,可避免估计性能对栅格初始位置的依赖,改善由于静态的栅格化造成的失配问题,提高估计精度和分辨率;
3、本发明无需目标个数的先验信息及对协方差矩阵的满秩性要求,可适用于目标个数未知及小样本环境,具有很强的灵活性。在单样本条件下,该方法可实现目标的动态追踪;
4、本发明鲁棒性好,对信噪比要求低,可实现低发射功率(小于噪声功率)目标的检测。
附图说明
图1为本发明方法的流程图;
图2为本发明互质阵列及其子阵列示意图;
图3为本发明的仿真实验中归一化功率谱图;
图4为本发明的仿真实验中估计精度图a;
图5为本发明的仿真实验中估计精度图b。
具体实施方式
下面结合附图对本发明作进一步的详细描述:
本发明方法的流程图如图1所示,具体实现过程如下:
(1)构建非线性互质阵列模型,获取接收信号;
(2)计算接收信号的协方差矩阵并进行向量化处理,建立虚拟差分阵列;
(3)对角度域进行栅格化处理,建立稀疏优化问题;
(4)以迭代方式动态调整栅格点位置,直至满足终止条件。分析稀疏解,确定最终波达方向。
本发明的DOA估计方法过程如下:
1.互质阵列及其接收信号
本发明所涉及的互质阵列如图2所示,该阵列可分解为两个均匀线性子阵列,其中子阵列1包含M1个天线,相邻天线间隔为M2λ/2,子阵列2包含2M2个天线,相邻天线间隔为M1λ/2,这里M1和M2为互质的正整数,λ表示载波的波长,子阵列1和2的全体构成非线性互质阵列,由于子阵列1和2共享第一个天线,因此互质阵列的天线数为M=M1+2M2-1。
有未知个数(假定为K)的不相关目标从不同方向Θ=[θ1,θ2,…,θK]到达互质阵列,该阵列在时刻t(1≤t≤T)的接收信号为
其中,A=[a(θ1),a(θ2),…,a(θK)]表示由阵列位置确定的已知的阵列流形矩阵,s(t)=[s1(t),s2(t),…,sK(t)]T表示K个目标的发射信号矢量,n(t)为独立同分布加性高斯白噪声矢量,上标T表示转置。
2.协方差矩阵及虚拟差分阵列
接收信号x(t)的协方差矩阵可表征为
其中,和σ2分别表示第k个信号的功率和噪声功率,E表示期望,上标H表示共轭转置, 表示M1+2M2-1维单位阵,矩阵Rxx的第(m,n)项为可视为在lm-ln处存在的虚拟天线的接收信号,在Rxx中,由第m个和第n个天线的差值产生的虚拟天线在位置lm-ln(1≤m,n≤M1+2M2-1),lm和ln分别表示第m和第n个天线的实际位置。
对Rxx进行向量化处理,则有
z=vec(Rxx)=Φ(θ1,θ2,…,θK)p+σ21n, (3)
其中 和vec表示向量化处理,表示Kronecker乘积,z为虚拟差分阵列的接收信号,Φ(θ1,θ2,…,θK)表示虚拟差分阵列的阵列流行矩阵。
3.栅格化处理,建立稀疏优化问题
为利用稀疏方法进行DOA估计,需要将目标角度域进行栅格化,因此可建立稀疏优化问题为
其中||·||0表示0-范数,该优化问题表示的含义为在给定静态栅格点,即在给定虚拟阵列流行矩阵的前提下,用尽可能少的信号功率p去重构接收信号z,该问题是NP-难问题,计算量极大,
为此,利用对数和函数去近似(4)中的0-范数,并构建无约束优化问题为
其中ε>0用来确定对数函数的存在性,||·||2表示最小二乘代价函数,η>0衡量稀疏性和最小二乘代价之间的均衡,由于对数函数的非凸性,优化问题(5)容易陷入局部最优解,进一步,利用对数和函数的凸上界函数
来代替优化问题中的对数和函数,式中,pi为p中的第i个元素,为pi在第t次迭代中的估计,替代中忽略掉与未知变量{p,θg,σ2}无关的项之后,优化问题(5)可转化为
其中
4.DOA估计迭代实现
该迭代方法的具体实施步骤如下:
步骤1:初始化离散角度集合θg,(0),对应的信号功率集合和噪声功率σ2 , (0),并令t=1,
步骤2:根据当前角度集合θg,(t-1)和噪声功率σ2,(t-1),优化问题(7)对p求导并置零,计算得
步骤3:根据当前角度集合θg , (t-1)和信号功率估计计算噪声功率σ2 , (t)为
步骤4:根据当前信号功率估计构建对数和函数的凸上界函数更新
步骤5:将估计值σ2 , (t)和D(t)代入优化问题(7),优化问题变为
由于Φ(θg)是关于θg的非线性函数,直接获取最优θg难以实现,可采用迭代方式逐步靠近最优的θg,为此寻找新的估计θg,(t)满足下式
f(θg,(t))≤f(θg,(t-1))
θg , (t)可根据梯度下降法估计为
其中μ为较小的正数,
令t=t+1,
步骤6:若满足终止条件,算法结束,否则跳至步骤2。
以下通过仿真对比本发明方法和其他传统方法(如空间平滑方法)说明本发明的优越性能:
本发明仿真实验采用M1=5和M2=3的互质阵列模型,其天线总个数为10,存在K=11个等功率信号源,其角度分别为[-49.3,-37.2,-26.8,-17.3,-8.3,0.45,9.2,18.3,27.8,38.3,50.6]度,将角度域按照3度的间隔进行等间隔栅格化,信噪比定义为输入功率与噪声功率之比。
图3为本发明的仿真实验中归一化功率谱图,其中图2中的横坐标表示波达方向,纵坐标表示归一化能量,虚线表示真实角度,上图实线表示空间平滑算法估计的角度,下图实线表示本发明方法估计的角度,如图3所示,本发明方法能成功检测所有的目标,而空间平滑方法漏掉了其中一个,因此本发明方法具有更强的检测能力。
图4为本发明的仿真实验中估计精度图,定量分析了本发明的估计精度,图4中横坐标表示信噪比,纵坐标表示估计均方误差,图5中横坐标表示样本数目,纵坐标表示估计的均方误差,由图4可见,在不同信噪比条件下,本发明的估计均方误差均小于空间平滑算法的估计误差,在低信噪比环境下更加明显,由图5可见,在不同样本数目条件下,本发明的估计均方误差也均小于空间平滑算法的估计误差,在小样本条件下更加明显,显然,本发明的DOA估计精度高于现有的空间平滑算法。
Claims (4)
1.互质阵列中基于迭代稀疏重构的DOA估计方法,其特征在于:所述方法按以下步骤进行:
(1)建立接收阵列模型
用非线性的互质阵列作为接收阵列采集信号,其中互质阵列可分解为两个均匀线性子阵列,每个子阵列的阵元间隔均大于半波长;
(2)构建虚拟差分阵列
对互质阵列接收信号的协方差矩阵进行向量化处理,在不同阵元的差分位置形成虚拟阵元,可将原始阵列等效成具有更多虚拟阵元及更大孔径尺度的差分阵列;
(3)栅格化处理,建立稀疏优化问题
将目标区域进行栅格化处理,分析差分阵列中的接收信号,建立以对数和函数为优化目标的稀疏优化问题,兼顾估计结果的稀疏性和准确性:
为利用稀疏方法进行DOA估计,需要将目标角度域进行栅格化,D>>K,因此可建立稀疏优化问题为
其中||·||0表示0-范数,该优化问题表示的含义为在给定静态栅格点,即在给定虚拟阵列流行矩阵的前提下,用尽可能少的信号功率p去重构接收信号z,该问题是NP-难问题,计算量极大,
为此,利用对数和函数去近似(4)中的0-范数,并构建无约束优化问题为
其中ε>0用来确定对数函数的存在性,||·||2表示最小二乘代价函数,η>0衡量稀疏性和最小二乘代价之间的均衡,由于对数函数的非凸性,优化问题(5)容易陷入局部最优解,进一步,利用对数和函数的凸上界函数
来代替优化问题中的对数和函数,式中,pi为p中的第i个元素,为pi在第t次迭代中的估计,替代中忽略掉与未知变量{p,θg,σ2}无关的项之后,优化问题(5)可转化为
其中
(4)DOA估计迭代实现
利用对数和函数的凸紧上界重构原始稀疏问题,以迭代方式逐步更新栅格点位置来使目标函数最小化,逐渐逼近真实信号源的方向,直至满足终止条件,所述方法可动态调整栅格点的位置,能有效克服角度域的静态离散化造成的栅格失配问题。
2.根据权利要求1所述的互质阵列中基于迭代稀疏重构的DOA估计方法,其特征在于:所述的非线性互质阵列可分解为两个包含M1和2M2个天线的均匀子阵列,其对应的阵元间隔为M2和M1个半波长,这里M1和M2为互质的正整数。
3.根据权利要求1所述的互质阵列中基于迭代稀疏重构的DOA估计方法,其特征在于:所述差分阵列中的虚拟天线的位置由互质阵列中实际天线位置的差值确定,该差分阵列具有更大的孔径尺寸,可以用来检测比天线数目更多的目标。
4.根据权利要求1所述的互质阵列中基于迭代稀疏重构的DOA估计方法,其特征在于:所述的迭代稀疏重构的DOA估计方法中,利用对数和函数近似传统的0-范数来确保稀疏性,原问题可转化为迭代最小化对数和函数的凸上界函数问题,由此逐步更新角度域的栅格点的位置以逼近真实值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610237414.5A CN106021637B (zh) | 2016-04-15 | 2016-04-15 | 互质阵列中基于迭代稀疏重构的doa估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610237414.5A CN106021637B (zh) | 2016-04-15 | 2016-04-15 | 互质阵列中基于迭代稀疏重构的doa估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106021637A CN106021637A (zh) | 2016-10-12 |
CN106021637B true CN106021637B (zh) | 2019-02-19 |
Family
ID=57082030
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610237414.5A Active CN106021637B (zh) | 2016-04-15 | 2016-04-15 | 互质阵列中基于迭代稀疏重构的doa估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106021637B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106501773B (zh) * | 2016-12-23 | 2018-12-11 | 云知声(上海)智能科技有限公司 | 基于差分阵列的声源方向定位方法 |
CN106803124B (zh) * | 2017-01-21 | 2020-04-07 | 中国海洋大学 | 基于流形正则和范数正则的领域迁移极限学习机方法 |
CN107015190A (zh) * | 2017-03-01 | 2017-08-04 | 浙江大学 | 基于虚拟阵列协方差矩阵稀疏重建的互质阵列波达方向估计方法 |
CN107329110B (zh) * | 2017-08-24 | 2019-08-30 | 浙江大学 | 基于稀疏阵列直接内插的波达方向估计方法 |
CN108594164B (zh) * | 2017-11-30 | 2020-09-15 | 山东农业大学 | 一种平面阵列doa估计方法及设备 |
CN111007457B (zh) * | 2018-10-08 | 2020-10-16 | 哈尔滨工业大学 | 一种基于块稀疏贝叶斯模型的辐射源直接定位方法 |
CN109917330B (zh) * | 2019-03-04 | 2022-11-01 | 东南大学 | 一种存在相位误差时基于稀疏正交匹配追踪理论的到达角估计方法 |
CN110133656B (zh) * | 2019-06-06 | 2022-05-03 | 电子科技大学 | 一种基于互质阵列的分解与融合的三维sar稀疏成像方法 |
CN110749856B (zh) * | 2019-09-10 | 2021-07-27 | 杭州电子科技大学 | 一种基于零化去噪技术的互质阵欠定测向方法 |
CN111077522B (zh) * | 2019-12-25 | 2023-07-04 | 中国科学院上海微系统与信息技术研究所 | 一种目标检测方法、装置、稀疏阵列和成像设备 |
CN113030853B (zh) * | 2021-03-07 | 2022-07-08 | 中国人民解放军陆军工程大学 | 基于rss和aoa联合测量的多辐射源无源定位方法 |
CN117434511B (zh) * | 2023-12-13 | 2024-03-01 | 广东大湾区空天信息研究院 | 一种基于毫米波雷达的多目标角度解模糊方法及相关设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105182293A (zh) * | 2015-08-25 | 2015-12-23 | 西安电子科技大学 | 基于互质阵列mimo雷达doa与dod估计方法 |
CN105388462A (zh) * | 2015-11-20 | 2016-03-09 | 成都理工大学 | 一种基于互质阵列的虚拟波束形成方法 |
CN105403874A (zh) * | 2015-12-25 | 2016-03-16 | 西安电子科技大学 | 非均匀阵列欠定波达方向估计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9562968B2 (en) * | 2012-10-22 | 2017-02-07 | Saab-Sensis Corporation | Sensor system and method for determining target location using sparsity-based processing |
-
2016
- 2016-04-15 CN CN201610237414.5A patent/CN106021637B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105182293A (zh) * | 2015-08-25 | 2015-12-23 | 西安电子科技大学 | 基于互质阵列mimo雷达doa与dod估计方法 |
CN105388462A (zh) * | 2015-11-20 | 2016-03-09 | 成都理工大学 | 一种基于互质阵列的虚拟波束形成方法 |
CN105403874A (zh) * | 2015-12-25 | 2016-03-16 | 西安电子科技大学 | 非均匀阵列欠定波达方向估计方法 |
Non-Patent Citations (2)
Title |
---|
Partial spectral search-based DOA estimation method for co-prime linear arrays;Fenggang Sun;《Electronics Letters》;20151130;第51卷(第24期);第2053-2055页 |
基于部分校准极化敏感阵列的信号DOA和极化参数迭代估计;常文秀 等;《电子与信息学报》;20080831;第30卷(第8期);第1893-1896页 |
Also Published As
Publication number | Publication date |
---|---|
CN106021637A (zh) | 2016-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106021637B (zh) | 互质阵列中基于迭代稀疏重构的doa估计方法 | |
CN106788653B (zh) | 一种基于协方差矩阵重构的自适应波束形成方法 | |
CN105652273B (zh) | 一种基于混合匹配追踪算法的mimo雷达稀疏成像算法 | |
CN107450047B (zh) | 嵌套阵下基于未知互耦信息的压缩感知doa估计方法 | |
Chen et al. | Hybrid spherical-and planar-wave modeling and DCNN-powered estimation of terahertz ultra-massive MIMO channels | |
CN103744076B (zh) | 基于非凸优化的mimo雷达动目标检测方法 | |
CN110045323B (zh) | 一种基于矩阵填充的互质阵稳健自适应波束形成算法 | |
Jeng | Near-field scattering by physical theory of diffraction and shooting and bouncing rays | |
CN106501765B (zh) | 一种基于平方和与半定规划的最大似然波达方向估计方法 | |
CN107576931B (zh) | 一种基于协方差低维度迭代稀疏重构的相关/相干信号波达方向估计方法 | |
Wang et al. | Adaptive beamforming with compressed sensing for sparse receiving array | |
CN108089147B (zh) | 改进的短波单站定位方法 | |
CN107121665B (zh) | 一种基于稀疏阵的近场相干源的无源定位方法 | |
CN111273229B (zh) | 基于低秩矩阵重建的水声宽频散射源的定位方法 | |
Cai et al. | Trajectory-aided maximum-likelihood algorithm for channel parameter estimation in ultrawideband large-scale arrays | |
CN112379327A (zh) | 一种基于秩损估计的二维doa估计与互耦校正方法 | |
CN107290732A (zh) | 一种量子大爆炸的单基地mimo雷达测向方法 | |
Poursheikhali et al. | Source localization in inhomogeneous underwater medium using sensor arrays: Received signal strength approach | |
Guo et al. | Robust DOA estimation for burst impulsive noise | |
Qi et al. | Time-frequency DOA estimation of chirp signals based on multi-subarray | |
Reaz et al. | A comprehensive analysis and performance evaluation of different direction of arrival estimation algorithms | |
CN104360338A (zh) | 一种基于对角加载的阵列天线自适应波束形成方法 | |
Li et al. | Robust DOA Estimation in Satellite Systems in Presence of Coherent Signals Subject to Low SNR | |
Zoubir et al. | Efficient DSPE algorithm for estimating the angular parameters of coherently distributed sources | |
Ma et al. | Computation-efficient 2-D DOA estimation algorithm with array motion strategy |
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 |