CN111736118B - 一种线列阵阵列扩展方法 - Google Patents
一种线列阵阵列扩展方法 Download PDFInfo
- Publication number
- CN111736118B CN111736118B CN202010708433.8A CN202010708433A CN111736118B CN 111736118 B CN111736118 B CN 111736118B CN 202010708433 A CN202010708433 A CN 202010708433A CN 111736118 B CN111736118 B CN 111736118B
- Authority
- CN
- China
- Prior art keywords
- array
- signal
- sub
- odd
- linear
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 239000011159 matrix material Substances 0.000 claims abstract description 37
- 238000003491 array Methods 0.000 claims abstract description 19
- 238000001514 detection method Methods 0.000 claims abstract description 12
- 230000004044 response Effects 0.000 claims abstract description 5
- 239000013598 vector Substances 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 description 4
- 238000011426 transformation method Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
-
- 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
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computer Networks & Wireless Communication (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Signal Processing (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种线列阵阵列扩展方法,属于阵列波束形成技术领域。所述方法包括:获取线列阵的接收信号;将线列阵按奇、偶阵元分成两个子阵;计算奇、偶阵元构成的两个子阵的接收信号协方差矩阵;根据接收信号协方差矩阵构造扩展接收阵列信号;利用常规波束形成或最小方差无畸变响应波束形成方法进行波束形成和目标检测。本发明解决了现有线列阵扩展方法低信噪比下性能下降,鲁棒性差;要求信源信号具有非圆对称性,应用有局限性;阵列扩展方法复杂,计算量大等问题。
Description
技术领域
本发明涉及阵列信号处理方法,具体涉及一种线列阵阵列扩展方法。
背景技术
阵列信号处理技术广泛应用于雷达、通信、声呐等众多军用和民用领域,是这些领域中的一个焦点问题。信噪比是影响阵列信号处理性能的一个关键因素。对于给定的阵列配置,物理阵元数和阵列孔径是确定的,如何提高低信噪比下弱目标的探测性能是工程应用中亟待解决的问题。利用阵列扩展技术,通过虚拟增加阵元数目,实现阵列孔径扩展,提高低信噪比下弱目标的探测性能不失为一种有效途径。目前常用的阵列扩展技术主要有:基于高阶累积量法、基于内插变换法、基于信源和接收信号特性的重构数据法和基于宽带信号时延特性和信号非圆对称性的线列阵扩展方法等。
高阶累积量法有良好稳定的阵列扩展特性,且根据高斯信号的高阶累积量为零的性质,能够对***中的高斯噪声起到抑制作用,使得其能够在不同的高斯噪声环境下具有良好的估计性能,但该方法计算量庞大,存在很多的冗余信息,当阵元数增加时,这些冗余信息可能使阵元之间发生耦合;内插变换法通过增加阵元数实现阵列扩展,但子区域的长度和内插变换角度的步长难以确定,如何兼顾算法计算量和计算精度是内插变换法的一大难点;基于信源和接收信号特性重构数据法的提出旨在解决内插变换法中存在的“角度敏感”、“插值步长敏感”以及运算量大等问题,但主要应用于L型阵列等特殊的二维阵列,如何将该方法运用到一维阵列模型有待研究;目前常用的一维线列阵扩展方法有基于线列阵时延特性的扩展方法和基于接收信号非圆特性的线列阵扩展方法,前者可以在阵列扩展的同时,实现相干信源的解相干,但在低信噪比条件下鲁棒性较差,后者要求入射信号具有非圆对称性,应用中有局限性。
因此,有必要提出一种新的线列阵的阵列扩展方法。
发明内容
发明目的:为克服上述现有技术的不足,本发明提供一种低信噪比条件下鲁棒性好、运算量小的线列阵阵列扩展方法。
技术方案:为实现上述目的,本发明采用如下技术方案:
一种线列阵阵列扩展方法,包括以下步骤:
进一步的,所述步骤(1)中接收阵为均匀线列阵,阵元数为2M,K个远场信号源以平面波入射到接收阵,阵列接收信号为:
进一步的,所述步骤(2)中,线列阵按奇、偶阵元分成两个子阵,两个子阵的接收信号分别为:
进一步的,所述步骤(3)中,奇、偶阵元构成的两个子阵的接收信号协方差矩阵为:
进一步的,所述步骤(4)中,构造的扩展接收阵列信号为:
进一步的,所述步骤(5)中,扩展接收阵列的信号协方差矩阵为:
更进一步的,所述步骤(6)中,利用常规波束形成CBF或最小方差无畸变响应MVDR波束形成方法对扩展阵列进行波束形成和目标检测,CBF和MVDR波束形成的波束输出功率为:
有益效果:本发明利用各阵元噪声的不相关性,以及线列阵的旋转不变性重构阵列数据,实现了一维线列阵的扩展。通过阵列扩展提高弱目标的探测性能,可以大大降低探测***的硬件成本。与现有技术相比,本发明在阵列扩展的同时,进一步减小了噪声的影响,降低了波束输出旁瓣,有利于低信噪比下弱目标的检测,提高了鲁棒性;阵列扩展方法简单,计算量小,不要求信源信号具有非圆对称性,提高了实用性。
附图说明
图1是本发明实施例提供的线列阵阵列扩展方法流程图;
图2是根据本发明实施例得到的扩展阵常规波束形成波束输出功率随方位变化曲线示意图;
图3是根据本发明实施例得到的扩展阵最小方差无畸变响应波束形成波束输出功率随方位变化曲线示意图;
图4是根据本发明实施例得到的扩展阵检测概率随信噪比变化的曲线示意图。
具体实施方式
下面结合附图对本发明的技术方案作进一步说明。
参照图1,在一实施例中,针对一维线列阵的阵列扩展方法包括以下步骤:
步骤(1),获取线列阵接收信号:
给定目标源信号为带限噪声,频带为2000Hz-2400Hz,采样频率16000Hz,采用阵元数为6的均匀线列阵,声速1500m/s,阵元间距为中心频率波长的一半,单次快拍数为1024,信源数为1,目标方位为60°,信噪比为-10dB。
阵列接收信号为:
步骤(2),按奇、偶阵元将线列阵划分为两个子阵:
奇、偶阵元构成的两个子阵的接收信号为:
步骤(3),计算奇、偶阵元构成的两个子阵的接收信号协方差矩阵:
其中,和分别为奇、偶阵元构成的两个子阵的阵列流型矩阵,可由阵列流型矩阵A得到,如上所述,A是维矩阵,提取A中的奇数行K维得到,将A中的偶数行K维得到,为源信号协方差矩阵。上标T表示转置,上标H表示共轭。由于各阵元的噪声不相关,奇、偶阵元构成的子阵的接收信号协方差矩阵去除了噪声的影响。
步骤(4),构造扩展接收阵列信号:
其中,,,表示Y1和Y2是从中取出元素构成的矩阵,括号中逗号前面表示所取的行,逗号后面表示所取的列,:表示取所有的,则和分别为协方差矩阵的前M-1列和后M-1列构成的矩阵,前M-1列指的是1到M-1列,后M-1列指的是2到M列,矩阵Y3、Y4由线列阵的旋转不变性得到,分别为矩阵Y2、Y1的共轭和行、列倒置;具体地,Y3是由Y2先共轭再行列倒置得到,Y4是由Y1先共轭再行列倒置得到。
步骤(5),计算扩展接收阵列的信号协方差矩阵:
步骤(6),计算常规波束形成(Conventional Beamforming, CBF)和最小方差无畸变响应波束形成(Minimum Variance Distortionless Response, MVDR)的波束输出功率:
图2示出了扩展线列阵与未扩展线列阵常规波束形成波束输出功率随方位的变化,图3为扩展线列阵与未扩展线列阵MVDR波束形成波束输出功率随方位的变化,由图2和图3可见,扩展线列阵具有更低的旁瓣,表明对噪声的抑制能力增强;图4为扩展线列阵与未扩展线列阵检测概率随信噪比的变化,由100次蒙特卡洛实验得到,由图4可见,低信噪比下扩展线列阵的检测概率明显高于未扩展线列阵,有利于弱目标的检测。
Claims (4)
1.一种线列阵阵列扩展方法,其特征在于,包括以下步骤:
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010708433.8A CN111736118B (zh) | 2020-07-22 | 2020-07-22 | 一种线列阵阵列扩展方法 |
AU2020356795A AU2020356795B2 (en) | 2020-07-22 | 2020-12-21 | Method of array extension for linear arrays |
PCT/CN2020/137877 WO2022016796A1 (zh) | 2020-07-22 | 2020-12-21 | 一种线列阵阵列扩展方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010708433.8A CN111736118B (zh) | 2020-07-22 | 2020-07-22 | 一种线列阵阵列扩展方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111736118A CN111736118A (zh) | 2020-10-02 |
CN111736118B true CN111736118B (zh) | 2020-11-17 |
Family
ID=72657023
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010708433.8A Expired - Fee Related CN111736118B (zh) | 2020-07-22 | 2020-07-22 | 一种线列阵阵列扩展方法 |
Country Status (3)
Country | Link |
---|---|
CN (1) | CN111736118B (zh) |
AU (1) | AU2020356795B2 (zh) |
WO (1) | WO2022016796A1 (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111736118B (zh) * | 2020-07-22 | 2020-11-17 | 东南大学 | 一种线列阵阵列扩展方法 |
CN113447887B (zh) * | 2021-09-02 | 2022-01-21 | 广东省新一代通信与网络创新研究院 | 全空间定位方法、装置、设备与计算机可读存储介质 |
CN114858271B (zh) * | 2022-07-05 | 2022-09-23 | 杭州兆华电子股份有限公司 | 一种声音探测用的阵列放大方法 |
CN116244561A (zh) * | 2023-03-14 | 2023-06-09 | 哈尔滨工程大学 | 一种基于组合子阵协方差矩阵的线列阵扩展方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103091671A (zh) * | 2013-01-15 | 2013-05-08 | 西安电子科技大学 | 基于非同心电磁矢量阵列雷达的两维波达方向估计方法 |
CN103308877A (zh) * | 2013-05-21 | 2013-09-18 | 西安电子科技大学 | 基于分离式偶极子对阵列的多参数估计方法 |
CN104698433A (zh) * | 2015-03-16 | 2015-06-10 | 电子科技大学 | 基于单快拍数据的相干信号doa估计方法 |
CN105824002A (zh) * | 2016-04-15 | 2016-08-03 | 西安电子科技大学 | 基于嵌套式子阵阵列的波达方向估计方法 |
CN109283486A (zh) * | 2018-10-26 | 2019-01-29 | 大连理工大学 | 基于广义复数相关熵的相干分布式非圆信号角度参数联合估计方法 |
JP2020091263A (ja) * | 2018-12-07 | 2020-06-11 | 株式会社東芝 | レーダシステム及び信号処理方法 |
CN111273219A (zh) * | 2020-02-17 | 2020-06-12 | 华南理工大学 | 一种基于圆与非圆混合信号的一维水下波达方向估计方法 |
CN111580039A (zh) * | 2020-03-29 | 2020-08-25 | 重庆邮电大学 | 基于非圆信号的单基地展开互质阵列mimo雷达doa估计方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103546221B (zh) * | 2013-10-25 | 2015-09-16 | 东南大学 | 一种宽带相干信号波达角估计方法 |
CN111736118B (zh) * | 2020-07-22 | 2020-11-17 | 东南大学 | 一种线列阵阵列扩展方法 |
-
2020
- 2020-07-22 CN CN202010708433.8A patent/CN111736118B/zh not_active Expired - Fee Related
- 2020-12-21 WO PCT/CN2020/137877 patent/WO2022016796A1/zh active Application Filing
- 2020-12-21 AU AU2020356795A patent/AU2020356795B2/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103091671A (zh) * | 2013-01-15 | 2013-05-08 | 西安电子科技大学 | 基于非同心电磁矢量阵列雷达的两维波达方向估计方法 |
CN103308877A (zh) * | 2013-05-21 | 2013-09-18 | 西安电子科技大学 | 基于分离式偶极子对阵列的多参数估计方法 |
CN104698433A (zh) * | 2015-03-16 | 2015-06-10 | 电子科技大学 | 基于单快拍数据的相干信号doa估计方法 |
CN105824002A (zh) * | 2016-04-15 | 2016-08-03 | 西安电子科技大学 | 基于嵌套式子阵阵列的波达方向估计方法 |
CN109283486A (zh) * | 2018-10-26 | 2019-01-29 | 大连理工大学 | 基于广义复数相关熵的相干分布式非圆信号角度参数联合估计方法 |
JP2020091263A (ja) * | 2018-12-07 | 2020-06-11 | 株式会社東芝 | レーダシステム及び信号処理方法 |
CN111273219A (zh) * | 2020-02-17 | 2020-06-12 | 华南理工大学 | 一种基于圆与非圆混合信号的一维水下波达方向估计方法 |
CN111580039A (zh) * | 2020-03-29 | 2020-08-25 | 重庆邮电大学 | 基于非圆信号的单基地展开互质阵列mimo雷达doa估计方法 |
Non-Patent Citations (1)
Title |
---|
联合互协方差矩阵的快速波达方向估计;闫锋刚 等;《***工程与电子技术》;20180430;第40卷(第4期);1阵列结构及数据模型以及2联合互协方差矩阵的DOA估计 * |
Also Published As
Publication number | Publication date |
---|---|
AU2020356795B2 (en) | 2022-12-15 |
CN111736118A (zh) | 2020-10-02 |
WO2022016796A1 (zh) | 2022-01-27 |
AU2020356795A1 (en) | 2022-02-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111736118B (zh) | 一种线列阵阵列扩展方法 | |
CN108344967B (zh) | 基于互质面阵的二维波达方向快速估计方法 | |
CN108896954B (zh) | 互质阵中一种基于联合实值子空间的波达角估计方法 | |
CN110927661A (zh) | 基于music算法的单基地展开互质阵列mimo雷达doa估计方法 | |
CN107544051A (zh) | 嵌套阵列基于k‑r子空间的波达方向估计方法 | |
CN111239678A (zh) | 一种基于l型阵列的二维doa估计方法 | |
CN104020440B (zh) | 基于l型干涉式线性阵列的二维波达角估计方法 | |
CN106872934B (zh) | L型电磁矢量传感器阵列解相干esprit参数估计方法 | |
CN108872930B (zh) | 扩展孔径二维联合对角化doa估计方法 | |
CN109491009B (zh) | 一种光纤组合阵及基于光纤组合阵的栅瓣抑制方法 | |
CN109270483B (zh) | 三维阵虚拟扩展相干源二维波达方向估计方法 | |
CN113820654B (zh) | 基于波束域降维的s波段雷达目标低仰角doa估计方法 | |
CN111830458B (zh) | 一种平行线阵单快拍二维测向方法 | |
CN115795249A (zh) | 一种基于l型互质阵列的二维doa估计方法 | |
CN111431575B (zh) | 基于常规波束形成的来波方向稀疏重构方法 | |
CN114563760A (zh) | 一种基于sca阵型的二阶超波束形成方法、设备及介质 | |
CN114371441A (zh) | 虚拟阵列波达方向估计方法、装置、产品及存储介质 | |
CN106950530B (zh) | 基于互协方差的声矢量传感器阵列解相干参数估计方法 | |
CN112698263A (zh) | 一种基于正交传播算子的单基地互质mimo阵列doa估计算法 | |
Zhai et al. | DOA estimation of noncircular signals for coprime linear array via polynomial root-finding technique | |
CN114841024B (zh) | 一种基于矢量扩展互质阵的波束形成方法及装置 | |
Yan et al. | Underdetermined direction of arrival estimation with coprime array constrained by approximated zero norm | |
CN114397620B (zh) | 改进和差非均匀阵列的高精度波达方向估计方法 | |
CN114755627A (zh) | 基于压缩感知和最小化处理相结合的互质面阵二维doa估计方法 | |
CN116244561A (zh) | 一种基于组合子阵协方差矩阵的线列阵扩展方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20201117 |