CN111431575B - 基于常规波束形成的来波方向稀疏重构方法 - Google Patents
基于常规波束形成的来波方向稀疏重构方法 Download PDFInfo
- Publication number
- CN111431575B CN111431575B CN202010144833.0A CN202010144833A CN111431575B CN 111431575 B CN111431575 B CN 111431575B CN 202010144833 A CN202010144833 A CN 202010144833A CN 111431575 B CN111431575 B CN 111431575B
- Authority
- CN
- China
- Prior art keywords
- array
- signal
- matrix
- algorithm
- sparse reconstruction
- 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
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000001228 spectrum Methods 0.000 claims abstract description 23
- 239000011159 matrix material Substances 0.000 claims description 29
- 230000005540 biological transmission Effects 0.000 claims 1
- 238000005070 sampling Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/08—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
- H04B7/0837—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station using pre-detection combining
- H04B7/0842—Weighted combining
- H04B7/086—Weighted combining using weights depending on external parameters, e.g. direction of arrival [DOA], predetermined weights or beamforming
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
- G01S3/14—Systems for determining direction or deviation from predetermined direction
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
Description
技术领域
本发明涉及波束形成的来波方向稀疏重构方法。属于DOA估计技术领域。
背景技术
信号的波达方向(direction of arrival,DOA)估计技术在通信、雷达、声纳、地震等领域有着广泛的应用,在国防军事和国民经济领域发挥着越来越重要的作用。目前常见的方位估计算法有波束形成法、子空间法、稀疏重构法。
波束形成法的基本原理为利用设计的波束在空间内进行扫描,从而得到空间谱,即信号在空间方向上的能量分布,根据空间谱即可判断信号的入射方向。,但是常规波束形成算法受限于“瑞利限,即无法分辨处于一个波束宽度内的信号,其角度分辨率取决于阵列尺寸,阵列尺寸越大,波束宽度越窄,则角度分辨率越好。然而实际应用中,阵列尺寸不可能无限增大。为了突破瑞利限的限制,进一步提高DOA估计的角度分辨率,有学者提出了子空间类算法,最为经典的是多重信号分类(MUSIC)算法和旋转不变子空间(ESPRIT)算法,子空间类算法利用了阵列信号的二阶统计量即协方差矩阵,先对其进行分解,从特征向量中重构出信号子空间和噪声子空间,然后利用信号子空间和噪声子空间的特性估计信号方位,极大的提高了角度分辨率。子空间类算法的缺点为由于其利用了信号的二阶统计量,因此计算量较大。近年来,稀疏重构算法在DOA估计领域得到了广泛的应用,相比传统的DOA估计算法,稀疏重构算法能够在少快拍数、低信噪比和相关信源等不利条件下具有良好的性能。稀疏重构算法利用阵列信号的空间稀疏特性,设计一个完备的冗余字典,字典中的每一列称为‘原子’,通过字典中的原子的线性组合重构信号,从而估计信号的方位。冗余字典的设计依赖阵列的几何结构因此稀疏重构算法对阵列误差的鲁棒性较差。
发明内容
本发明是为了解决DOA估计中现有的稀疏重构算法存在冗余字典的设计依赖阵列的几何结构而导致稀疏重构算法对阵列误差的鲁棒性较差的问题。进而提出一种基于常规波束形成的来波方向稀疏重构方法。
基于常规波束形成的来波方向稀疏重构方法,包括以下步骤:
将角度空间进行均匀划分,入射信号是其中的一个角度值,基于接收信号的空间谱,通过正交匹配跟踪算法估计目标方位。
进一步地,所述利用基阵接收目标发射的信号得到的信号为:
X(t)=[x1(t),x2(t),...,xN(t)]
其中,接收信号X(t)是一个N×L的矩阵,L为接收信号的快拍数;xk(t)表示每个基元对应的接收信号,k=1,2,...,N。
进一步地,所述接收目标发射信号的基阵是由N个基元组成均匀直线阵列。
进一步地,所述基阵的相邻阵元间距为接收信号波长的一半。
进一步地,所述基阵的阵元个数大于待测目标的个数。
进一步地,所述基阵的输出信号为:
进一步地,通过正交匹配跟踪算法估计目标方位的过程中需要设计冗余字典:
其中,S为稀疏矩阵,只在目标方位处有值,其余位置均为零。
P=AS
进一步地,通过正交匹配跟踪算法估计目标方位的过程如下:
1.初始化:被选择的原子所构成的字典子集D=空集,残差r=P,S全部元素初始化为0;
2.找出冗余字典A中与残差r内积最大的原子ai,记录下这个原子的索引i;把ai添加入字典子集D中;
3.由最小二乘可得系数向量的解,即下列方程的最小二乘解:
5.判断是否达到迭代终止条件;若是,结束算法,输出稀疏矩阵S;若否,转至第2步;所述迭代终止条件为迭代次数达到迭代次数阈值或残差的范数小于残差阈值;
通过正交匹配追踪算法能够确定S的值,根据S中的值与方位角度的一一对应关系,即可求得目标方位。
有益效果:
本发明将常规波束形成算法和稀疏重构算法结合,充分利用常规波束形成算法的鲁棒性和稀疏重构算法在低快拍、低信噪比环境下的良好性能,本发明基于常规波束形成的来波方向稀疏重构方法,能够在低快拍、低信噪比、高阵元位置差的条件下精确估计来波方位。本发明解决了DOA估计中现有的稀疏重构算法存在冗余字典的设计依赖阵列的几何结构而导致稀疏重构算法对阵列误差的鲁棒性较差的问题。
附图说明
图1为接收信号基本模型示意图;
图2为空间角度划分方式示意图。
具体实施方式
本发明的目的在于提出一种基于常规波束形成的来波方向稀疏重构方法,能够在低快拍、低信噪比、高阵元位置差的条件下精确估计来波方位。
具体实施方式一:
本实施方式所述基于常规波束形成的来波方向稀疏重构方法,包括以下步骤:
接收信号模型如附图1所示,假设接收基阵是由N个基元组成均匀直线阵列,阵元间距为d,相邻阵元间距为接收信号波长的一半,阵元个数大于待测目标的个数;待测目标发射的信号为窄带信号;一号基元为参考基元,有一个目标声源位于远场,因此入射到基阵表面的声波可视为远场平面波且声波与基阵法线的夹角为θ,则基阵接收到的信号可以表示为:
X(t)=[x1(t),x2(t),...,xN(t)]
其中,接收信号X(t)是一个N×L的矩阵,L为接收信号的快拍数;xk(t)表示每个基元对应的接收信号,k=1,2,...,N;
通过常规波束形成算法计算阵列接收信号的空间谱:所述的常规波束形成算法是波束形成算法的一种,根据波束指向对接收信号进行时延补偿,然后将各个阵元的信号相加获得阵列的输出信号,输出信号的功率即为该方向上谱值,改变波束指向即可获得整个空间的方位谱。
改变波束指向进行空间波束扫描即可获得接收信号的空间谱:
其中,M为空间角度扫描的角度个数,角度个数越多,精度越高。
本实施方式使用的是常规波束形成算法,其在方位估计中的优点在于该方法对阵列误差有最佳的鲁棒性;将角度空间进行均匀划分,得到接收信号的空间谱后,即可通过正交匹配跟踪算法估计目标方位:
将角度空间进行均匀划分的间隔需要根据所需的方位估计精度进行设定,精度要求越高则间隔越小。在本实施方式中角度空间为[-90°,90°],将角度空间等分为M份(如附图2所示),入射信号是其中的一个角度值,则入射信号的方位在整个角度空间中是稀疏的,且入射信号的空间谱在整个角度空间的谱中也是稀疏的;
根据正交匹配追踪的基本原理,首先设计冗余字典,即:
则:
P=AS
正交匹配追踪算法的基本原理为给定一个有所有方位组成的冗余字典矩阵,字典中的每一列称为原子,信号可以通过字典中原子的线性组合表示。正交匹配追踪算法能够迭代地找到具有剩余最大绝对内积的原子,将该原子添加到一个新的集合中并更新残差(实际观察值与估计值或拟合值之间的差),直到迭代次数达到一定值或残差的范数足够小。所述正交匹配追踪算法的流程如下:
1.初始化:被选择的原子所构成的字典子集D=空集(每次迭代被选中的原子添加进这个子集作为新的一列),残差r=P,系数S全部元素初始化为0;
2.找出冗余字典A中与残差r内积最大的那个原子ai,记录下这个原子的索引i(即该原子在字典D中的位置);把ai添加入字典子集D中(每迭代一次D中就会增加一个新的原子,这个加入的新原子就是此次迭代中发现的与残差内积最大的那个原子);
3.由最小二乘可得系数向量的解,即下列方程的最小二乘解:
5.判断是否达到迭代终止条件;若是,结束算法,输出稀疏矩阵S;若否,转至第2步;所述迭代终止条件为迭代次数达到迭代次数阈值或残差的范数小于残差阈值;
通过正交匹配追踪算法可以求得S的值,根据S中的值与方位角度的一一对应关系,即可求得目标方位。
本发明的核心内容在于将常规波束形成方法与正交匹配追踪算法结合,充分利用常规波束形成算法的鲁棒性和稀疏重构算法在低快拍、低信噪比环境下的良好性能,从而在低快拍、低信噪比、高阵元位置差的条件下精确估计来波方位。鲁棒性指对与基阵本身的误差如基元位置误差、通道幅相误差等具有更好的适应性。快拍指一次接受数据的采样点个数,数值等于采样时间和采样频率的乘积。
最后应该说明的是,以上实例仅用以描述本发明的技术方案而不是对本技术方法进行限制,本发明在应用上可以延伸为其他的修改、变化、应用和实施例,并且因此认为所有这样的修改、变化、应用、实施例都在本发明的精神和教导范围内。
Claims (5)
1.基于常规波束形成的来波方向稀疏重构方法,其特征在于,包括以下步骤:
X(t)=[x1(t),x2(t),...,xN(t)]
其中,接收信号X(t)是一个N×L的矩阵,L为接收信号的快拍数;xk(t)表示每个基元对应的接收信号,k=1,2,...,N,N接收基阵中基元的数量;
将角度空间进行均匀划分,入射信号是其中的一个角度值,基于接收信号的空间谱,通过正交匹配跟踪算法估计目标方位;
通过正交匹配跟踪算法估计目标方位的过程中需要设计冗余字典:
其中,只在目标方位处有值,其余位置均为零;M为空间角度扫描的角度个数,m表示空间角度扫描的角度序号;
P=AS
通过正交匹配跟踪算法估计目标方位的过程如下:
1.初始化:被选择的原子所构成的字典子集D=空集,残差r=P,S全部元素初始化为0;
2.找出冗余字典A中与残差r内积最大的原子ai,记录下这个原子的索引i;把ai添加入字典子集D中;
3.由最小二乘可得系数向量的解,即下列方程的最小二乘解:
5.判断是否达到迭代终止条件;若是,结束算法,输出稀疏矩阵S;若否,转至第2步;所述迭代终止条件为迭代次数达到迭代次数阈值或残差的范数小于残差阈值;
通过正交匹配追踪算法能够确定S的值,根据S中的值与方位角度的一一对应关系,即可求得目标方位。
2.根据权利要求1所述的基于常规波束形成的来波方向稀疏重构方法,其特征在于,所述接收目标发射信号的基阵是由N个基元组成均匀直线阵列。
3.根据权利要求2所述的基于常规波束形成的来波方向稀疏重构方法,其特征在于,所述基阵的相邻阵元间距为接收信号波长的一半。
4.根据权利要求3所述的基于常规波束形成的来波方向稀疏重构方法,其特征在于,所述基阵的阵元个数大于待测目标的个数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010144833.0A CN111431575B (zh) | 2020-03-04 | 2020-03-04 | 基于常规波束形成的来波方向稀疏重构方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010144833.0A CN111431575B (zh) | 2020-03-04 | 2020-03-04 | 基于常规波束形成的来波方向稀疏重构方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111431575A CN111431575A (zh) | 2020-07-17 |
CN111431575B true CN111431575B (zh) | 2022-02-11 |
Family
ID=71547485
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010144833.0A Active CN111431575B (zh) | 2020-03-04 | 2020-03-04 | 基于常规波束形成的来波方向稀疏重构方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111431575B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116973900B (zh) * | 2023-07-17 | 2024-08-16 | 哈尔滨工程大学 | 冰下声纳目标方位估计方法、系统、终端及计算机存储介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009059279A1 (en) * | 2007-11-01 | 2009-05-07 | University Of Maryland | Compressive sensing system and method for bearing estimation of sparse sources in the angle domain |
CN105974358A (zh) * | 2016-05-25 | 2016-09-28 | 天津商业大学 | 基于压缩感知的智能天线doa估计方法 |
CN106443621A (zh) * | 2016-09-13 | 2017-02-22 | 哈尔滨工程大学 | 冲击噪声条件的基于正交匹配稀疏重构的相干信源动态doa追踪方法 |
CN104937856B (zh) * | 2013-12-26 | 2018-06-05 | 华为技术有限公司 | 一种估计到达角的方法、装置及电子设备 |
CN109782246A (zh) * | 2018-12-31 | 2019-05-21 | 深圳市华讯方舟太赫兹科技有限公司 | 波达方向估计方法及装置、雷达、可读存储介质 |
CN110824415A (zh) * | 2019-11-19 | 2020-02-21 | 中国人民解放军国防科技大学 | 一种基于多发多收阵列的稀疏波达方向角度估计方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104459606A (zh) * | 2014-12-25 | 2015-03-25 | 武汉大学 | 一种阵列空间信号稀疏构建及其重建方法 |
CN107566305B (zh) * | 2017-08-15 | 2020-03-31 | 南京邮电大学 | 一种低复杂度的毫米波系统信道估计方法 |
-
2020
- 2020-03-04 CN CN202010144833.0A patent/CN111431575B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009059279A1 (en) * | 2007-11-01 | 2009-05-07 | University Of Maryland | Compressive sensing system and method for bearing estimation of sparse sources in the angle domain |
CN104937856B (zh) * | 2013-12-26 | 2018-06-05 | 华为技术有限公司 | 一种估计到达角的方法、装置及电子设备 |
CN105974358A (zh) * | 2016-05-25 | 2016-09-28 | 天津商业大学 | 基于压缩感知的智能天线doa估计方法 |
CN106443621A (zh) * | 2016-09-13 | 2017-02-22 | 哈尔滨工程大学 | 冲击噪声条件的基于正交匹配稀疏重构的相干信源动态doa追踪方法 |
CN109782246A (zh) * | 2018-12-31 | 2019-05-21 | 深圳市华讯方舟太赫兹科技有限公司 | 波达方向估计方法及装置、雷达、可读存储介质 |
CN110824415A (zh) * | 2019-11-19 | 2020-02-21 | 中国人民解放军国防科技大学 | 一种基于多发多收阵列的稀疏波达方向角度估计方法 |
Non-Patent Citations (1)
Title |
---|
阵列信号DOA估计算法研究;陈松;《中国优秀硕士学位论文全文数据库》;20170215;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111431575A (zh) | 2020-07-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107290732B (zh) | 一种量子大爆炸的单基地mimo雷达测向方法 | |
CN113189592B (zh) | 考虑幅相互耦误差的车载毫米波mimo雷达测角方法 | |
CN114814830B (zh) | 一种基于鲁棒主成分分析降噪的米波雷达低仰角测高方法 | |
CN113111304B (zh) | 强冲击噪声下基于量子射线机理的相干分布源测向方法 | |
Li et al. | Parameter estimation based on fractional power spectrum density in bistatic MIMO radar system under impulsive noise environment | |
CN112285647A (zh) | 一种基于稀疏表示与重构的信号方位高分辨估计方法 | |
CN113820653B (zh) | 基于动态和差波束的米波雷达低仰角目标doa估计方法 | |
CN111431575B (zh) | 基于常规波束形成的来波方向稀疏重构方法 | |
CN108594165B (zh) | 一种基于期望最大化算法的窄带信号波达方向估计方法 | |
CN111722178B (zh) | 一种基于指向性模型数值求解的远场窄带信号来波方向估计方法 | |
CN111812607A (zh) | 基于波束空间的米波mimo雷达低仰角估计方法 | |
CN111381212A (zh) | 基于子阵划分的虚拟超短基线定位方法 | |
CN105572631B (zh) | 基于多波位联合处理的最大似然目标doa估计方法 | |
CN114167346B (zh) | 基于协方差矩阵拟合阵元扩展的doa估计方法及系统 | |
CN114563760B (zh) | 一种基于sca阵型的二阶超波束形成方法、设备及介质 | |
CN113093098B (zh) | 基于lp范数补偿的轴向不一致矢量水听器阵列测向方法 | |
CN115575941A (zh) | 面向稀疏阵列的频控阵mimo雷达目标参数估计方法 | |
CN114184999B (zh) | 一种互耦小孔径阵列的生成式模型处理方法 | |
CN113791379B (zh) | 嵌套阵列非高斯环境下的正交匹配追踪doa估计方法 | |
CN109298381A (zh) | 一种基于变分贝叶斯推断的互质阵相干信号方位角估计方法 | |
CN114755628A (zh) | 非均匀噪声下声矢量传感器阵列波达方向估计方法 | |
CN114325560A (zh) | 波束扫描雷达超分辨目标测向方法 | |
CN113820654A (zh) | 基于波束域降维的s波段雷达目标低仰角doa估计方法 | |
CN110231590B (zh) | 一种基于dft变换的阵列目标角度测量方法 | |
Li et al. | DOA Estimation Based on Sparse Reconstruction via Acoustic Vector Sensor Array under Non-uniform Noise |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |