CN116050099A - 一种基于压缩感知的面向冲击噪声的嵌套阵doa估计方法 - Google Patents
一种基于压缩感知的面向冲击噪声的嵌套阵doa估计方法 Download PDFInfo
- Publication number
- CN116050099A CN116050099A CN202211682978.1A CN202211682978A CN116050099A CN 116050099 A CN116050099 A CN 116050099A CN 202211682978 A CN202211682978 A CN 202211682978A CN 116050099 A CN116050099 A CN 116050099A
- Authority
- CN
- China
- Prior art keywords
- impact noise
- signal
- matrix
- angle
- compressed sensing
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 32
- 239000011159 matrix material Substances 0.000 claims abstract description 28
- 238000009826 distribution Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 5
- 239000006185 dispersion Substances 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 5
- 230000036039 immunity Effects 0.000 claims description 3
- 238000005315 distribution function Methods 0.000 claims description 2
- 238000012163 sequencing technique Methods 0.000 claims description 2
- 238000004422 calculation algorithm Methods 0.000 description 7
- 238000004088 simulation Methods 0.000 description 7
- 238000003491 array Methods 0.000 description 5
- 238000009499 grossing Methods 0.000 description 4
- 241000284708 Sarcophaga alpha Species 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 238000007635 classification algorithm Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000021615 conjugation Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- 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/10—Noise analysis or noise optimisation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Computing Systems (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Algebra (AREA)
- Geometry (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于压缩感知的面向冲击噪声的嵌套阵DOA估计方法,其特征包括如下步骤:首先,构造嵌套阵模型;其次,构造信号模型添加冲击噪声,计算相控分数低阶矩(PhaseFractionalLowerOrderMoment);接着,将该矩阵虚拟化,并将其排序删除冗余;然后,构造过完备角度集合,定义残差,索引集,重建列集合,计算最佳匹配角标;最后,不断迭代,更新索引集,残差,重建列集合,最后角度集合中所对应角度即为方位角估计。本发明设计的基于压缩感知的面向冲击噪声的嵌套阵DOA估计方法能够将单快拍虚拟信号变换为稀疏信号,同时用于估计冲击噪声环境下的方位角。
Description
技术领域
本发明属于DOA估计领域,尤其涉及一种基于压缩感知的面向冲击噪声的嵌套阵DOA估计方法。
背景技术
高斯分布假设是目前信号处理中的主要问题,它的概率密度函数可以用均值和方差两种形式来表示,从而使信号的理论分析更加容易;在高斯假定的基础上,信号的处理往往是线性的,一般都能得到最优的结果。但是,现实环境中有很多非高斯特征的脉冲干扰,比如通讯线路的瞬间尖峰,大气放电造成的大气噪音,以及部分由雷达产生的杂波。
现由于规则性阵列结构和奈奎斯特采样速率这两方面的约束,均匀阵列是目前研究领域最为常见的天线阵列结构。但是,由于奈奎斯特采样定理,均匀阵列中相邻阵元的间距需小于等于半波长,故阵列孔径相对受到限制;此外,天线阵元个数也限制了均匀线阵的自由度,只有在信号源数目小于天线阵元个数的条件下才能进行有效的信号处理。为了解决上述问题,本发明提出了一种基于压缩感知的面向冲击噪声的嵌套阵DOA估计方法。
发明内容
发明目的:为了解决现有技术存在的问题,本发明提出了一种基于压缩感知的面向冲击噪声的嵌套阵DOA估计方法,在能够有效估计冲击噪声情况下的方位角,同时降低了计算复杂度。
技术方案:本发明所述的一种基于压缩感知的面向冲击噪声的嵌套阵DOA估计方法,包括以下步骤:
一种基于压缩感知的面向冲击噪声的嵌套阵DOA估计方法,包括以下步骤:
3)从观测矩阵AΘ选择与观测信号z内积最大的列,并通过最小二乘计算残差r继续在矩阵AΘ选择与残差r最匹配的列,反复迭代直到获得信源的角度估计。
优选的,步骤1)中所述嵌套阵阵元模型包括串联的第一级均匀线阵和第二级均匀线阵;所述第一级均匀线阵的阵元数为N1,阵元间隔为d1,所述第二级均匀线阵阵元数为N2,阵元间隔为d2,且满足d2=(N1+1)d1;将接收信号X(t)定义为:
X(t)=AS(t)+N(t)
其中A为方向向量,S(t)为发射信号,N(t)为冲击噪声;
将冲击噪声认为服从SαS对称α稳定分布:
获取相控分数低阶矩阵:
其中xi(t)和xj(t)分别表示接收信号x(t)第i个和第j个快拍所代表的元素,*表示复数的共轭。
优选的,步骤2)中,假设Θ={θ1,θ2,…,θD}(D>>K)为过完备冗余字典,其中D为所有可能入射方向数,θi(i=1,2,…,D)为某一方向入射角度,K为待估计信源数,据此将过完备冗余字典Θ扩展成AΘ,即:
其中AΘ为观测方向矩阵,列向量dv为虚拟阵元位置,则有观测信号z:
同时由于l0范数数值上无法实现,且抗噪性能差,故采用l1代替l0,即
优选的,步骤3)中的实现过程为:
3.3)如果t≤K,就回到步骤3.2)继续更新残差,如果t>K,则信源的角度估计就是ΓK在角度集合中所对应的角度。
有益效果:与现有技术相比,本发明的有益效果:本发明所提算法的估计精度优于空间平滑多重信号分类算法(Spatial Smoothing Multiple Signal Classification),空间平滑旋转不变性信号参数估计算法(Spatial Smoothing Estimating SignalParameterviaRotational Invariance Techniques);与SS-MUSIC,SS-ESPRIT算法相比,所提出的方法具有更多的自由度,估计误差更低;本发明可以在保证估计性能的情况下显著降低计算复杂度。
附图说明
图1为本发明的流程图;
图2为嵌套阵模型图;
图3为SαS分布图;
图4为本发明与传统子空间方法计算复杂度对比示意图;
图5为本发明与传统子空间方法在不同信噪比下的求根均方误差性能示意图;
图6为本发明与传统子空间方法在不同快拍数下的求根均方误差性能示意图。
具体实施方式
下面结合附图对本发明做进一步详细说明。
本发明提供一种基于压缩感知的面向冲击噪声的嵌套阵DOA估计方法,如图1所示,具体包括以下步骤:
步骤1:构造如图2所示的嵌套阵阵元模型;得到接收信号信息X(t),并计算相控分数低阶矩阵Rxx。
如图2所示,二级嵌套阵由两个均匀线阵串联构成,第一级均匀线阵的阵元数为N1,阵元间隔为d1,第二级均匀线阵的阵元数为N2阵元间隔为d2,且d2=(N1+1)d1,即该二级嵌套阵是一个阵元位置由集合P1={md1,m=1,2,…,N1}
和集合P2={n(N1+1)d1,n=1,2,....N2}确定的线性阵列。接收信号为
X(t)=AS(t)+N(t)
其中冲击噪声可以认为服从对称α稳定(SαS)分布
其中α是限制于以下值的特征指数(0<α≤2),是位置参数,是分布的离散度。弥散参数决定了分布在其位置参数周围的分布。特征指数决定了分布的形状(见图3)。α越小,alpha稳定密度的尾部越重。这意味着r.v.,s的下列具有小特征指数的阿尔法稳定分布是高度脉冲的。除了高斯(α=2)、柯西(α=1)和Levy(α=1/2)分布之外,一般的SαS概率密度函数没有封闭式。然而,稳定密度函数的幂级数展开是可用的。对于稳定过程,只有p<α阶存在的矩。对于位置参数为零且色散为零的SαS随机变量,这些矩由
E{|x|p}=C(p,α)γp/αfor0<p<α
当α≠2时,则E{|x|p}不存在当p≧α。基于数据的SαS估计方法不能应用。
如果我们把z写成极坐标z=rejθ,那么就很容易表示为z<a>=raejθ,因此PFLOM只作用于其操作数的大小,保留其相位。我们还定义z的αPFLOM的共轭为图3清楚地显示了调谐参数α对数据的影响。设置为α=1时,数据保持不变,设置为α=0时,信号中所有的振幅信息都被消除,只保留其复杂相位。
PFLOM操作对观测数据起到振幅压缩器的作用,以减少其脉冲效应。对空间信号的区别是,所有信号的观测向量除以相同参数矢量的(范数)。PFLOM操作取决于信号本身的振幅下面的定理定义了基于PFLOM的新一类协方差矩阵估计。
其中,列向量dv表明虚拟阵元位置,那么
由于向量z和矩阵AΘ这两项已知,z作为观测信号,AΘ为观测矩阵,pΘ为待解的稀疏信号,因此上式可看作一个压缩感知模型,可通过下式解决:
由于l0范数数值上无法实现,且抗噪性能差,故采用l1代替l0,即
步骤4:最佳匹配列脚标由计算获得,其中AΘ,j表示AΘ的第j列。更新索引集Γt=Γt-1∪{γt}并重建更新残差rt=z-Atpt=z-At(At TAt)-1At Tz,并令t=t+1,如果t≤K,就重复4,如果t>K则进入5;
步骤5:这时候方位角估计就是ΓK在角度集合中对应的角度。
采取多次Monte Carlo实验对算法进行角度估计性能评估。设求根均方误差(RootMean Square Error,RMSE)为
图5为本发明所述方法与传统子空间算法SS-MUSIC和SS-ESPRIT求根均方误差随信噪比变化的性能曲线图。仿真条件为:3个信源,嵌套阵是由两个阵元数为5的均匀线阵构成,快拍数为200,alpha=1,仿真800次。从图5可以看出,本发明实现了更高的估计精度。
图6为本发明所述方法与传统子空间算法SS-MUSIC和SS-ESPRIT求根均方误差随快拍数变化的性能曲线图。仿真条件为:3个信源,嵌套阵是由两个阵元数为5的均匀线阵构成,信噪比为10,alpha=1,仿真800次。从图6可以看出,本发明所提方法相比传统子空间算法SS-MUSIC和SS-ESPRIT具有更好的估计性能。
综上所述,从仿真效果图的分析可知,本发明提出的基于压缩感知的面向冲击噪声的嵌套阵DOA估计方法具有更高的空间自由度和估计精度,能估计更多目标。此外,该方法无需谱峰搜索,显著降低了计算复杂度。
上面结合附图对本发明的实施方式作了详细说明,但是本发明并不限于上述实施方式,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下做出各种变化。
Claims (4)
2.如权利要求1所述的一种基于压缩感知的面向冲击噪声的嵌套阵DOA估计方法,其特征在于,步骤1)中所述嵌套阵阵元模型包括串联的第一级均匀线阵和第二级均匀线阵;所述第一级均匀线阵的阵元数为N1,阵元间隔为d1,所述第二级均匀线阵阵元数为N2,阵元间隔为d2,且满足d2=(N1+1)d1;将接收信号X(t)定义为:
X(t)=AS(t)+N(t)
其中A为方向向量,S(t)为发射信号,N(t)为冲击噪声;
将冲击噪声认为服从SαS对称α稳定分布:
获取相控分数低阶矩阵:
其中xi(t)和xj(t)分别表示接收信号x(t)第i个和第j个快拍所代表的元素,*表示复数的共轭。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211682978.1A CN116050099A (zh) | 2022-12-27 | 2022-12-27 | 一种基于压缩感知的面向冲击噪声的嵌套阵doa估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211682978.1A CN116050099A (zh) | 2022-12-27 | 2022-12-27 | 一种基于压缩感知的面向冲击噪声的嵌套阵doa估计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116050099A true CN116050099A (zh) | 2023-05-02 |
Family
ID=86132512
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211682978.1A Pending CN116050099A (zh) | 2022-12-27 | 2022-12-27 | 一种基于压缩感知的面向冲击噪声的嵌套阵doa估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116050099A (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107450047A (zh) * | 2017-08-28 | 2017-12-08 | 电子科技大学 | 嵌套阵下基于未知互耦信息的压缩感知doa估计方法 |
CN110806557A (zh) * | 2019-11-22 | 2020-02-18 | 中国人民解放军国防科技大学 | 一种冲击噪声背景下单基地mimo雷达doa估计方法 |
WO2021139208A1 (zh) * | 2020-01-08 | 2021-07-15 | 华南理工大学 | 一种基于特定频率组合信号的一维doa估计方法 |
CN115453450A (zh) * | 2022-07-16 | 2022-12-09 | 电子科技大学长三角研究院(湖州) | 一种嵌套阵下对抗互藕的稳健doa估计方法 |
-
2022
- 2022-12-27 CN CN202211682978.1A patent/CN116050099A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107450047A (zh) * | 2017-08-28 | 2017-12-08 | 电子科技大学 | 嵌套阵下基于未知互耦信息的压缩感知doa估计方法 |
CN110806557A (zh) * | 2019-11-22 | 2020-02-18 | 中国人民解放军国防科技大学 | 一种冲击噪声背景下单基地mimo雷达doa估计方法 |
WO2021139208A1 (zh) * | 2020-01-08 | 2021-07-15 | 华南理工大学 | 一种基于特定频率组合信号的一维doa估计方法 |
CN115453450A (zh) * | 2022-07-16 | 2022-12-09 | 电子科技大学长三角研究院(湖州) | 一种嵌套阵下对抗互藕的稳健doa估计方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109655799B (zh) | 基于iaa的协方差矩阵向量化的非均匀稀疏阵列测向方法 | |
CN110045323B (zh) | 一种基于矩阵填充的互质阵稳健自适应波束形成算法 | |
CN107544051A (zh) | 嵌套阵列基于k‑r子空间的波达方向估计方法 | |
AU2009251043A1 (en) | A method and system of data modelling | |
CN112698264A (zh) | 增广互质阵列脉冲噪声环境下相干信源的doa估计方法 | |
CN103885050B (zh) | 基于缩放字典的回波信号参数估计方法 | |
CN109239646B (zh) | 一种冲击噪声环境下连续量子水蒸发的二维动态测向方法 | |
CN105699950B (zh) | 基于自适应迭代前后向平滑共轭梯度的雷达杂波抑制方法 | |
CN113376569B (zh) | 基于最大似然的嵌套阵稀疏表示波达方向估计方法 | |
CN107843875A (zh) | 基于奇异值分解降噪的贝叶斯压缩感知雷达数据融合方法 | |
CN109143154A (zh) | 一种应用于l型阵列的信号二维doa与频率联合估计方法 | |
Chen et al. | A DOA estimation algorithm based on eigenvalues ranking problem | |
CN106501765A (zh) | 一种基于平方和与半定规划的最大似然波达方向估计方法 | |
CN106526529A (zh) | 导向矢量失配情况下基于稀疏表示的波达方向估计方法 | |
CN104155629B (zh) | 一种冲击噪声背景下小快拍数信号波达方向估计方法 | |
CN107290732A (zh) | 一种量子大爆炸的单基地mimo雷达测向方法 | |
CN112444773A (zh) | 基于空域融合的压缩感知二维doa估计方法 | |
CN108387880A (zh) | 一种基于多尺度有向Lyapunov指数的海杂波背景下弱小目标检测方法 | |
CN105572629B (zh) | 一种适用于任意阵列结构的低运算复杂度的二维测向方法 | |
Xiao et al. | A weighted forward-backward spatial smoothing DOA estimation algorithm based on TLS-ESPRIT | |
CN112731280B (zh) | 互质阵列混合噪声环境下的esprit-doa估计方法 | |
Yan et al. | Computationally efficient direction of arrival estimation with unknown number of signals | |
CN106908760B (zh) | 基于阵列自相关矩阵的单站无源定位方法 | |
CN112904298A (zh) | 一种基于局部网格分裂的网格偏离空时自适应处理方法 | |
CN116050099A (zh) | 一种基于压缩感知的面向冲击噪声的嵌套阵doa估计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20230502 |
|
RJ01 | Rejection of invention patent application after publication |