CN113030868A - 一种毫米波雷达角度超分辨方法 - Google Patents
一种毫米波雷达角度超分辨方法 Download PDFInfo
- Publication number
- CN113030868A CN113030868A CN202110330715.3A CN202110330715A CN113030868A CN 113030868 A CN113030868 A CN 113030868A CN 202110330715 A CN202110330715 A CN 202110330715A CN 113030868 A CN113030868 A CN 113030868A
- Authority
- CN
- China
- Prior art keywords
- matrix
- signal
- subspace
- covariance matrix
- millimeter wave
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 24
- 239000011159 matrix material Substances 0.000 claims abstract description 118
- 239000013598 vector Substances 0.000 claims abstract description 39
- 238000012545 processing Methods 0.000 claims abstract description 15
- 238000001228 spectrum Methods 0.000 claims abstract description 13
- 238000000819 phase cycle Methods 0.000 claims abstract description 8
- 238000005070 sampling Methods 0.000 claims abstract description 7
- 238000000354 decomposition reaction Methods 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 6
- 230000003595 spectral effect Effects 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000001514 detection method Methods 0.000 abstract description 5
- 230000035945 sensitivity Effects 0.000 abstract description 3
- 230000009466 transformation Effects 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000012795 verification Methods 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
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/93—Radar or analogous systems specially adapted for specific applications for anti-collision purposes
- G01S13/931—Radar or analogous systems specially adapted for specific applications for anti-collision purposes of land vehicles
-
- 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/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
- G06F17/142—Fast Fourier transforms, e.g. using a Cooley-Tukey type algorithm
-
- 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/15—Correlation function computation including computation of convolution operations
- G06F17/156—Correlation function computation including computation of convolution operations using a domain transform, e.g. Fourier transform, polynomial transform, number theoretic transform
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Algebra (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Computing Systems (AREA)
- Computer Networks & Wireless Communication (AREA)
- Discrete Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种毫米波雷达角度超分辨方法,步骤包括:获取采样数据;对采样数据进行FFT处理,得到距离相序列和第一信号协方差矩阵;对第一信号协方差矩阵相加,重构得到第二信号协方差矩阵;将第二信号协方差矩阵分解为信号和噪声子空间矩阵,并求权值向量;利用权值向量重构噪声子空间矩阵,利用广义逆矩阵处理信号子空间矩阵;利用子空间矩阵构成空间谱密度函数方程;遍历角度寻找峰值以估算角度。本发明的毫米波雷达发射功率低、分辨率与灵敏度高、雷达尺寸小,利用信号子空间与噪声子空间的信息优化MUSIC算法,提高了角分辨率,显著提高雷达远距离模式下的检测能力。
Description
技术领域
本发明属于雷达技术领域,尤其涉及一种毫米波雷达角度超分辨方法。
背景技术
近年来,车载毫米波雷达的应用范围大幅度拓展,毫米波雷达在盲区报警、倒车警告以及辅助驾驶等方面应用较为广泛。然而受限于车载雷达的尺寸、功率、天线设计等问题,车载毫米波雷达在在长距离模式下,无法有效区分出相同径向距离下,切向距离相近、速度相近的多个目标,从而导致了雷达的角分辨率下降,影响雷达性能。常用的MUSIC算法形成的空间谱主瓣较大,容易掩盖角度相近的两个目标,对角分辨率的提高有限。因此需要对算法进行优化,以提高毫米波雷达的角分辨率。
发明内容
有鉴于此,本专利提出一种毫米波雷达角度超分辨方法,对雷达回波ADC信号进行采集与距离FFT处理,得到一维距离相序列,并处理得到噪声协方差矩阵。然后对协方差矩阵进行复共轭变换后,与原协方差矩阵进行相加,得到新的协方差矩阵。对协方差矩阵进行奇异值分解,得到特征值向量。对特征值向量进行区分操作,得到信号和噪声子空间。通过对信号子空间进行广义逆变换、对噪声子空间进行特征向量加权等操作得到空间谱函数。最后结合角度范围与空间谱函数,通过寻求峰值来得到到达角的估计值。
本发明公开的一种毫米波雷达角度超分辨方法,方法包括以下步骤:
对毫米波雷达回波信号进行ADC数据采集,得到采样数据;
对所述采样数据进行一维FFT处理,得到距离相序列,并处理得到第一信号协方差矩阵;
对所述第一信号协方差矩阵进行复共轭运算后与所述第一信号协方差矩阵进行相加,重构得到第二信号协方差矩阵;
对所述第二信号协方差矩阵进行奇异值分解,根据特征值向量将所述第二信号协方差矩阵分解为信号和噪声子空间矩阵,并求得权值向量;
利用所述权值向量重构噪声子空间矩阵,利用广义逆矩阵处理信号子空间矩阵;
利用所述子空间矩阵构成空间谱密度函数方程;
遍历角度寻找峰值以估算角度。
进一步的,所述对所述第一信号协方差矩阵进行复共轭运算后与所述第一信号协方差矩阵进行相加,重构得到第二信号协方差矩阵的步骤包含以下计算过程:
所述第一信号协方差矩阵为R XX ,定义交换矩阵:
进一步的,所述根据特征值向量将所述第二信号协方差矩阵分解为信号和噪声子空间矩阵,并求得权值向量的计算过程为:
已知阵元数为n,信源数为m,取前m大的特征值构成信号子空间特征值对角矩阵:,取后(n-m)大的特征值构成噪声子空间特征值对角矩阵:,信号子空间为对应的前m个特征向量扩展成的矩阵: =span(q 1...q m );噪声子空间为对应的后(n-m)个特征向量扩展成的矩阵: = span(q m+1...q n );
进一步的,所述利用所述权值向量重构噪声子空间矩阵,利用广义逆矩阵处理信号子空间矩阵的计算步骤如下:
进一步的,所述空间谱密度函数的计算步骤如下:
本发明的有益效果:
毫米波雷达通过软件方法进行改进,不增加硬件成本,不需要通过提高发射功率去提高目标检测准确率,具有稳定的探测性能与良好的环境适应性,结构简单、发射功率低、分辨率与灵敏度高、雷达尺寸小。
通过对信号子空间进行特征重构的方式优化了MUSIC算法,构造了新的空间谱函数,极大的提高了远距离下毫米波雷达的雷达角分辨率,使得毫米波雷达可以更好的在远距离下区分目标。
附图说明
图1为本发明改进MUSIC算法的毫米波雷达角度超分辨方法流程图;
图2为本发明改进MUSIC算法的毫米波雷达角度超分辨方法区分结果。
具体实施方式
下面结合附图对本发明作进一步的说明,但不以任何方式对本发明加以限制,基于本发明教导所作的任何变换或替换,均属于本发明的保护范围。
本发明的目的是找到一种方法,使得空间谱主瓣减小,提高角分辨率。
有鉴于此,本发明提出一种毫米波雷达角度超分辨方法,对雷达回波ADC信号进行采集与距离FFT处理,得到一维距离相序列,并处理得到噪声协方差矩阵。然后对协方差矩阵进行复共轭变换后,与原协方差矩阵进行相加,得到新的协方差矩阵。对协方差矩阵进行奇异值分解,得到特征值向量。对特征值向量进行区分操作,得到信号和噪声子空间。通过对信号子空间进行广义逆变换、对噪声子空间进行特征向量加权等操作得到空间谱函数。最后结合角度范围与空间谱函数,通过寻求峰值来得到到达角的估计值。
如图1所示,本发明公开的一种毫米波雷达角度超分辨方法,方法包括以下步骤:
S10:对毫米波雷达回波信号进行ADC数据采集,得到采样数据;
S20:对采样数据进行一维FFT处理,得到距离相序列,并处理得到第一信号协方差矩阵;
S30:对第一信号协方差矩阵进行复共轭运算后与第一信号协方差矩阵进行相加,重构得到第二信号协方差矩阵;
第一信号协方差矩阵为R XX ,定义交换矩阵:
S40:对第二信号协方差矩阵进行奇异值分解,根据特征值向量将第二信号协方差矩阵分解为信号和噪声子空间矩阵,并求得权值向量;
已知阵元数为n,信源数为m,取前m大的特征值构成信号子空间特征值对角矩阵:,取后(n-m)大的特征值构成噪声子空间特征值对角矩阵:,信号子空间为对应的前m个特征向量扩展成的矩阵: =span(q 1...q m );噪声子空间为对应的后(n-m)个特征向量扩展成的矩阵: = span(q m+1...q n );
S50:利用权值向量重构噪声子空间矩阵,利用广义逆矩阵处理信号子空间矩阵;
S60:利用子空间矩阵构成空间谱密度函数方程;
S70:遍历角度寻找峰值以估算角度。
根据上述实施例的方法研制样机进行试验验证,试验结果如图2所示,本发明改进MUSIC算法的毫米波雷达角度超分辨方法相对于传统的MUSIC算法,空间谱主瓣明显减小,显著提高了角分辨率。
本发明的有益效果:
毫米波雷达通过软件方法进行改进,不增加硬件成本,不需要通过提高发射功率去提高目标检测准确率,具有稳定的探测性能与良好的环境适应性,结构简单、发射功率低、分辨率与灵敏度高、雷达尺寸小。
通过对信号子空间进行特征重构的方式优化了MUSIC算法,构造了新的空间谱函数,极大的提高了远距离下毫米波雷达的雷达角分辨率,使得毫米波雷达可以更好的在远距离下区分目标。
上述实施例为本发明的一种实施方式,但本发明的实施方式并不受所述实施例的限制,其他的任何背离本发明的精神实质与原理下所做的改变、修饰、代替、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。
Claims (5)
1.一种毫米波雷达角度超分辨方法,其特征在于,方法包括以下步骤:
对毫米波雷达回波信号进行ADC数据采集,得到采样数据;
对所述采样数据进行一维FFT处理,得到距离相序列,并处理得到第一信号协方差矩阵;
对所述第一信号协方差矩阵进行复共轭运算后与所述第一信号协方差矩阵进行相加,重构得到第二信号协方差矩阵;
对所述第二信号协方差矩阵进行奇异值分解,根据特征值向量将所述第二信号协方差矩阵分解为信号和噪声子空间矩阵,并求得权值向量;
利用所述权值向量重构噪声子空间矩阵,利用广义逆矩阵处理信号子空间矩阵;
利用所述子空间矩阵构成空间谱密度函数方程;
遍历角度寻找峰值以估算角度。
3.根据权利要求1所述的毫米波雷达角度超分辨方法,其特征在于,所述根据特征值向量将所述第二信号协方差矩阵分解为信号和噪声子空间矩阵,并求得权值向量的计算过程为:
已知阵元数为n,信源数为m,取前m大的特征值构成信号子空间特征值对角矩阵: ,取后(n-m)大的特征值构成噪声子空间特征值对角矩阵:,信号子空间为对应的前m个特征向量扩展成的矩阵: =span(q 1...q m );噪声子空间为对应的后(n-m)个特征向量扩展成的矩阵: = span(q m+1...q n );
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110330715.3A CN113030868B (zh) | 2021-03-29 | 2021-03-29 | 一种毫米波雷达角度超分辨方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110330715.3A CN113030868B (zh) | 2021-03-29 | 2021-03-29 | 一种毫米波雷达角度超分辨方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113030868A true CN113030868A (zh) | 2021-06-25 |
CN113030868B CN113030868B (zh) | 2024-06-21 |
Family
ID=76473320
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110330715.3A Active CN113030868B (zh) | 2021-03-29 | 2021-03-29 | 一种毫米波雷达角度超分辨方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113030868B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114428236A (zh) * | 2022-04-06 | 2022-05-03 | 南京楚航科技有限公司 | 一种车载毫米波雷达角度置信度估计方法 |
CN114706056A (zh) * | 2022-04-28 | 2022-07-05 | 珠海微度芯创科技有限责任公司 | 基于协方差拟合的毫米波雷达高分辨测角方法、装置和介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106788653A (zh) * | 2016-11-25 | 2017-05-31 | 西安电子科技大学 | 一种基于协方差矩阵重构的自适应波束形成方法 |
CN109581352A (zh) * | 2018-12-10 | 2019-04-05 | 南京莱斯电子设备有限公司 | 一种基于毫米波雷达的超分辨测角系统 |
CN109975807A (zh) * | 2019-03-27 | 2019-07-05 | 东南大学 | 一种适用于毫米波车载雷达的降维子空间测角方法 |
CN110082741A (zh) * | 2019-03-14 | 2019-08-02 | 哈尔滨工程大学 | 一种基于伪数据重构的超分辨波达角估计算法 |
US20190317191A1 (en) * | 2018-04-11 | 2019-10-17 | Infineon Technologies Ag | Human Detection and Identification in a Setting Using Millimiter-Wave Radar |
CN110673086A (zh) * | 2019-10-31 | 2020-01-10 | 上海无线电设备研究所 | 一种基于数字阵列雷达的二维角度超分辨方法 |
-
2021
- 2021-03-29 CN CN202110330715.3A patent/CN113030868B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106788653A (zh) * | 2016-11-25 | 2017-05-31 | 西安电子科技大学 | 一种基于协方差矩阵重构的自适应波束形成方法 |
US20190317191A1 (en) * | 2018-04-11 | 2019-10-17 | Infineon Technologies Ag | Human Detection and Identification in a Setting Using Millimiter-Wave Radar |
CN109581352A (zh) * | 2018-12-10 | 2019-04-05 | 南京莱斯电子设备有限公司 | 一种基于毫米波雷达的超分辨测角系统 |
CN110082741A (zh) * | 2019-03-14 | 2019-08-02 | 哈尔滨工程大学 | 一种基于伪数据重构的超分辨波达角估计算法 |
CN109975807A (zh) * | 2019-03-27 | 2019-07-05 | 东南大学 | 一种适用于毫米波车载雷达的降维子空间测角方法 |
CN110673086A (zh) * | 2019-10-31 | 2020-01-10 | 上海无线电设备研究所 | 一种基于数字阵列雷达的二维角度超分辨方法 |
Non-Patent Citations (1)
Title |
---|
仝卜匀: "基于时延估计技术的声源定位方法研究", 《中国优秀硕士学位论文全文数据库基础科学辑》, pages 13 - 22 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114428236A (zh) * | 2022-04-06 | 2022-05-03 | 南京楚航科技有限公司 | 一种车载毫米波雷达角度置信度估计方法 |
CN114428236B (zh) * | 2022-04-06 | 2022-07-01 | 南京楚航科技有限公司 | 一种车载毫米波雷达角度置信度估计方法 |
CN114706056A (zh) * | 2022-04-28 | 2022-07-05 | 珠海微度芯创科技有限责任公司 | 基于协方差拟合的毫米波雷达高分辨测角方法、装置和介质 |
Also Published As
Publication number | Publication date |
---|---|
CN113030868B (zh) | 2024-06-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zheng et al. | DOA estimation for coprime linear arrays: An ambiguity-free method involving full DOFs | |
CN107037392B (zh) | 一种基于压缩感知的自由度增加型互质阵列波达方向估计方法 | |
CN110412499B (zh) | 基于压缩感知理论下的rss算法的宽带doa估计算法 | |
CN107544051A (zh) | 嵌套阵列基于k‑r子空间的波达方向估计方法 | |
CN1752771A (zh) | 单通道相控阵接收信号重构及空间信号处理方法 | |
CN113030868B (zh) | 一种毫米波雷达角度超分辨方法 | |
CN109375154A (zh) | 一种冲击噪声环境下基于均匀圆阵的相干信号参数估计方法 | |
CN110082741B (zh) | 一种基于伪数据重构的超分辨波达角估计算法 | |
CN111308416B (zh) | 基于四阶累积量的近场非圆信源参数估计方法 | |
CN106980105B (zh) | 电磁矢量传感器阵列空间旋转解相干测向方法 | |
CN109521393A (zh) | 一种基于信号子空间旋转特性的波达方向估计算法 | |
Xiao et al. | A weighted forward-backward spatial smoothing DOA estimation algorithm based on TLS-ESPRIT | |
CN114779236A (zh) | 一种基于空间平滑music改进的米波雷达低仰角测高方法 | |
CN108594165B (zh) | 一种基于期望最大化算法的窄带信号波达方向估计方法 | |
Ganguly et al. | An efficient DOA estimation and jammer mitigation method by means of a single snapshot compressive sensing based sparse coprime array | |
Li et al. | Direction-of-arrival estimation for coherent signals exploiting moving coprime array | |
CN112327292A (zh) | 一种二维稀疏阵列doa估计方法 | |
He et al. | DOA estimation of wideband signals based on iterative spectral reconstruction | |
Pan et al. | Enhanced doa estimation with co-prime array in the scenario of impulsive noise: A pseudo snapshot augmentation perspective | |
Zhang et al. | Dimension-reduced radio astronomical imaging based on sparse reconstruction | |
WO2020140658A1 (zh) | 波达方向估计方法及装置、雷达、可读存储介质 | |
CN114184999B (zh) | 一种互耦小孔径阵列的生成式模型处理方法 | |
Krieger et al. | Multi-coset sparse imaging arrays | |
Ertan et al. | Spatial and temporal smoothing for covariance estimation in super-resolution angle estimation in automotive radars | |
Wang et al. | Coherent DOA estimation using dynamic metasurface antenna working in space isomeric mode |
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 |