CN105629235B - 多目标检测的组合波形汽车变道辅助系统的信号处理装置 - Google Patents
多目标检测的组合波形汽车变道辅助系统的信号处理装置 Download PDFInfo
- Publication number
- CN105629235B CN105629235B CN201511005838.0A CN201511005838A CN105629235B CN 105629235 B CN105629235 B CN 105629235B CN 201511005838 A CN201511005838 A CN 201511005838A CN 105629235 B CN105629235 B CN 105629235B
- Authority
- CN
- China
- Prior art keywords
- matrix
- value
- frequency
- wave
- fmcw1
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 32
- 230000008859 change Effects 0.000 title claims abstract description 24
- 238000012545 processing Methods 0.000 title claims abstract description 18
- 239000011159 matrix material Substances 0.000 claims description 118
- 230000035559 beat frequency Effects 0.000 claims description 35
- 238000004364 calculation method Methods 0.000 claims description 25
- 238000000034 method Methods 0.000 claims description 15
- 230000009471 action Effects 0.000 claims description 2
- 230000009466 transformation Effects 0.000 claims 1
- 238000013461 design Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 241001269238 Data Species 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 206010039203 Road traffic accident Diseases 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
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
- 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
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
多目标检测的组合波形汽车变道辅助系统的信号处理装置,属于信号处理领域,用于解决现有汽车变道辅助系统,很难实现对多目标的解算的问题,技术要点是:所述波形包括第一段恒频波CW1、第二段锯齿波FMCW1、第三段恒频波CW2、第四段锯齿波FMCW2。效果是:由于采用四段波形,可以实现多目标检测功能,并且可以检测出真实目标,去除虚假目标。
Description
技术领域
本发明属于信号处理领域,尤其涉及一种组合波形汽车变道辅助系统的信号处理装置。
背景技术
随着现代社会的进步,汽车的使用量以及使用率呈现出快速的增长,但是伴随着汽车引发的问题也越来越多,甚至越来越严重,例如,汽车对于能源的消耗、环境危害,以及交通安全等。其中交通安全问题,越来越引起人们的广泛关注,相应的汽车防撞功能,汽车变道辅助功能,自动泊车功能,自动巡航功能等被更过的应用到了汽车上。
汽车变道辅助系统,是为驾驶员在进行变更道路时提供有效地依据,可以有效地降低由于汽车变道而发生的交通事故问题。随着传感器技术的不断的发展,毫米波雷达传感器,由于其体较小、质量轻。可以在相对恶劣的雨雪天气使用的优点,被越来越多的应用于汽车中。同时线性调频连续波(LFMCW)雷达具有很高的速度分辨率与距离分辨率,因此线性调频连续波被更多的应用到毫米波雷达系统中。但是,线性调频连续波雷达在进行运动目标检测时,如果只采用单一的波形设计,比如单独的恒频波,三角波、锯齿波等波形,很难实现对多目标的解算问题。
发明内容
为了更好的解决现有的车变道辅助系统对多目标的解算的问题,本发明提供了一种多目标检测的组合波形汽车变道辅助系统的信号处理装置,以实现对多目标解算。
为了实现上述目的,本发明的技术方案是:一种多目标检测的组合波形汽车变道辅助系统的信号处理装置,所述波形包括第一段恒频波CW1、第二段锯齿波FMCW1、第三段恒频波CW2、第四段锯齿波FMCW2,该装置包括:
FFT计算模块.对各段波形,A/D采集到的IQ数据,进行FFT计算;
门限检测模块,将各段波形FFT变换后的复数模值做门限检测,输出过门限点位置。
多普勒频率值计算模块,计算得到恒频段的多普勒频率值;
差频频率值计算模块,计算得到锯齿波段的差频频率值。
有益效果:由于采用四段波形,可以实现多目标检测功能,并且可以检测出真实目标,去除虚假目标。
附图说明
图1恒频波CW与锯齿波FMCW在一个扫频周期范围内的频率变化图;
图2单目标的(R,V)空间图;
图3多目标的(R,V)空间图;
图4基于组和波形的汽车变道辅助系统信号处理流程图。
具体实施方式
实施例1:一种多目标检测的组合波形汽车变道辅助系统的信号处理装置,所述波形包括第一段恒频波CW1、第二段锯齿波FMCW1、第三段恒频波CW2、第四段锯齿波FMCW2,其中,第一段恒频波CW1和第三段恒频波CW2为恒频段,第二段锯齿波FMCW1和第四段锯齿波FMCW2为锯齿波段,该装置包括:FFT计算模块、门限检测模块、多普勒频率值计算模块、差频频率值计算模块。
FFT计算模块,对各段波形,A/D采集到的IQ数据,进行FFT计算;作为技术方案的优选,FFT计算模块对通道1中的第一段恒频波CW1、第二段锯齿波FMCW1、第三段CW2和第四段锯齿波FMCW2,A/D采集到的IQ数据,选取各段线性度高的256点数据,分别进行256点FFT,对通道2中的第一段恒频波CW1和第二段锯齿波FMCW1,A/D采集到的IQ数据,选取各段线性度高的256点数据,分别进行256点FFT。其中,线性度高应该用指线性度好来表示会更好,线性度好意思就是采集的数据存在一定的线性关系,线性度不好,就是该段的数据没有呈现出一定的线性关系,表现出的是非线性关系等特点。为了更准确的分析数据,所以在采集到的数据中只选取呈现出线性关系的部分做FFT分析。
门限检测模块,将各段波形FFT变换后的复数模值做门限检测,输出过门限点位置。作为技术方案的优选,其检测方法是:设通道1中,第一段恒频波CW1,有n1点过门限,则对应的第一段恒频波过门限点的位置矩阵ACW1=[a1,a2,…an1],第二段锯齿波FMCW1,有n2点过门限,则对应的第二段锯齿波FMCW1过门限点的位置矩阵为BFMCW1=[b1,b2,…bn2],第三段恒频波CW2,有n3点过门限,则对应的第三段恒频波CW2过门限点的位置矩阵为CCW2=[c1,c2,…cn3],第四段锯齿波FMCW2,有n4点过门限,则对应的第四段锯齿波FMCW2过门限点的位置矩阵为DFMCW2=[d1,d2,…dn4],设通道2中,第一段恒频波CW1,有n5点过门限,则对应的第一段恒频波CW1过门限点的位置矩阵为ECW1=[e1,e2,…en5],第二段锯齿波FMCW1,有n6点过门限,则对应的第二段锯齿波FMCW1过门限点的位置矩阵为GFMCW1=[g1,g2,…gn6]。若过门限的位置点等于1,则认为其是直流分量,不作为目标判定,直接剔除该位置点。
多普勒频率值计算模块,计算得到恒频段的多普勒频率值;优选的,计算得到横频段的多普勒频率值的方法是:通道1中,第一段恒频波CW过门限点的位置矩阵ACW1=[a1,a2,…an1],根据如下规则计算对应点上的多普勒频率值,得到多普勒频率矩阵为FACW1=[fda1,fda2,…fdan1];对于通道1中,第三段恒频波CW2过门限点的位置矩阵CCW2=[c1,c2,…cn3],根据如下规则计算对应点上的多普勒频率值,得到多普勒频率矩阵为FCCW1=[fdc1,fdc2,…fdcn3];
该规则为,若点数为1≤xi≤128(1≤i≤n),判断目标靠近,其对应点上的多普勒频率(1≤i≤n);若点数为128<xi≤256(1≤i≤n),判断目标远离,其对应点上的多普勒频率(1<i≤n)。
差频频率值计算模块,计算得到锯齿波段的差频频率值。优选的,计算得到锯齿波段的差频频率值的方法是:通道1中,第二段锯齿波FMCW1过门限点的位置矩阵BFMCW1=[b1,b2,…bn2],根据如下规则计算对应点上的差频频率矩阵FBFMCW1=[fob1,fob2,…fobn2],对于通道1中,第四段锯齿波FMCW2过门限点的位置矩阵DFMCW2=[d1,d2,…dn4],根据如下规则计算对应点上的差频频率矩阵FDFMCW1=[fod1,fod2,…fodn4];
该规则为,即若点数为1≤yj≤128(1≤j≤n),其对应点上的差频频率值(1≤j≤n);若点数为128<yj≤256(1≤j≤n),其对应点上的差频频率值(1<j≤n)。
本实施例是记载了一种中心频率在24GHz或77GHz,且基于恒频波调制的CW信号以及锯齿波调制的FMCW信号组合而成的波形(作为上述波形的优选形式),根据该调制波形实现一种汽车变道辅助系统的信号处理装置。由该装置设计的汽车变道辅助系统可以实现对汽车行驶后方,毫米波雷达波束覆盖范围内的目标,实现相对距离、相对速度以及方向角的解算,同时由于采用恒频波CW与锯齿波FMCW设计的组合波形,可以实现多目标的检测问题,使得系统对于多目标检测具有更好的准确性、快速性,保障了驾驶员在驾驶中变换道路时的安全性。
波形设计以及波形分析:
本实施例给出在中心频率f为24.128GHz工作频率下的具体波形图,第一段波形为恒频波CW1,工作频率为24.128GHz,第二段波形为上升的锯齿波FMCW1,工作频率变化范围为从24.128GHz变化到24.278GHz,带宽为150MHz,第三段为恒频波CW2,工作频率为24278GHz,第四段为锯齿波FMCW2,工作频率变化范围为从24.278GHz到24.128GHz。每一段的周期T为5ms。恒频波CW与锯齿波FMCW在一个扫频周期范围内的频率变化图如图1所示。
本实施例选择此设计波形的原因有:
(1)提高相对速度与相对距离的解算精度。
在第一段波形——恒频波CW1,根据恒频波的特性,可以获得由速度引起的多普勒频率值fd1,根据第二段波形——锯齿波FMCW1,可以计算出目标的差频频率值fo1。通过第一段波形获得的多普勒频率值fd1,可以计算出目标的相对速度值v1,通过第一段波形获得的多普勒频率值以及第二段波形获得的差频频率值fo1,可以计算出目标的相对距离R1。
同理,第三段波形——恒频波CW2,同样可以获得目标的多普勒频率值fd2,根据第四段波形——锯齿波FMCW2,可以计算出目标的差频频率值fo2。通过第三段波形获得的多普勒频率值fd2,可以计算出目标的相对速度值v2,通过第四段波形获得的多普勒频率值以及第二段波形获得的差频频率值fo2,可以计算出目标的相对距离R2。
利用通过四段组合波形,根据后期多目标匹配算法等处理最终获得真实目标的相对距离以及相对速度。
通过四段组合波形,获得相对速度v1和相对速度v2,相对距离R1和相对距离R2。理论上,v1=v2,R1=R2,但是由于设备使用环境的复杂性等因素,可能造成干扰等因素,造成测量的v1值与v2值存在偏差,R1与R2也同样存在偏差,如果单一的使用v1或是v2作为目标的相对速度,R1或是R2作为目标的相对速度,都会造成获取目标相对速度以及距离的不准确性。如图2,为单目标的(R,V)空间图,从图中可以看出,四条直线确定一个交点,所获得的相对距离值与相对速度值的准确性,要高于两条直线确定的相对距离值与相对速度值,同时,即使有一个波段计算得值相对不准确,其他三条直线确定的相对距离值和相对速度值,也会强过两条直线确定的相对距离值与相对速度值,使得系统具有一定的抗干扰能力以及鲁棒特性。因此,通过四个波段,可以有效的提高对目标速度以及距离解算的准确度,同时使系统具有抗干扰的能力以及鲁棒特性;
(2)由于采用四段波形,可以实现多目标检测功能,并且可以检测出真实目标,去除虚假目标。
由于真实目标的相对速度和相对距离与调频周期T无关,而虚假目标的距离和速度计算值与调频周期T有关系。所以,对于真实目标来说,在不同周期中,其真实值是一直存在,而在不同周期虚假值是随周期T变化的。如图3,为多目标的(R,V)空间图。从图2中可以看出,四个波段会更加准确的获得一个目标的相对距离和相对速度值,从图3中可以看出,采用四段波形可以实现对多目标的检测,同时可以有效地将真实目标检测出来,去除掉虚假目标。
通过如上两条证明,本实施例所设计的波形,可以有效地提高对目标相对距离值与相对速度值的解算的准确性,同时使得系统具有抗干扰特性以及具有鲁棒特性,同时可以有效地检测多目标,并且可以有效地获得真实目标,去除掉虚假目标。
实施例2:具有与实施例1相同的技术方案,更为具体的是:所述装置还包括:相对速度值计算模块、相对距离值计算模块、方向角计算模块。
相对速度值计算模块,根据得到的多普勒频率值,计算相对速度矩阵;
相对距离值计算模块,根据计算得到的多普勒频率矩阵与差频频率矩阵,计算得到相对距离矩阵;
方向角计算模块,进行多目标的方位角的计算。
实施例3:具有与实施例1-2任一项相同的技术方案,更为具体的是:根据各段计算得到的多普勒频率值矩阵,计算目标的速度矩阵的具体方法为:计算速度公式为其中,c为光速,c=3×108,f为中心频率f=24.128GHz。根据第一段恒频波CW1得到的多普勒频率矩阵FACW1=[fda1,fda2,…fdan1],得到其速度矩阵为VACW1=[va1,va2,…van1],根据第二段恒频波CW2得到的多普勒频率矩阵FCCW1=[fdc1,fdc2,…fdcn3],得到其速度矩阵为VCCW1=[vc1,vc2,…vcn3]。
实施例4:具有实施例1-3任一项相同的技术方案,更为具体的是:根据计算得到的多普勒频率矩阵与差频频率矩阵,计算得到相对距离矩阵的具体步骤是:计算距离公式为其中,T为每一段波形的作用时间,T=5ms,B为调频带宽,B=150MHz,fd为多普勒频率值,fo为差频频率值。
根据通道1中,第一段恒频波CW1得到的多普勒频率矩阵FACW1=[fda1,fda2,…fdan1]和第二段锯齿波FMCW1得到的差频频率矩阵FBFMCW1=[fob1,fob2,…fobn2],将多普勒矩阵中所有元素和差频频率矩阵中的所有元素进行一一配对计算相对距离矩阵,计算得到相对距离矩阵为其中raibj(1≤i≤n1,1≤j≤n2),表示是由第一段恒频波CW1得到的多普勒频率矩阵中第i个元素与第二段锯齿波FMCW1得到的差频频率矩阵中第j个元素进行计算得到的距离值。同理,对于第二段恒频波CW2得到的多普勒频率矩阵FCCW1=[fdc1,fdc2,…fdcn3]和第四段锯齿波FMCW2得到的差频频率矩阵FDFMCW2=[fod1,fod2,…fodn4],同样进行上述处理,最后得到相对距离矩阵为其中rcicj(1≤i≤n3,1≤j≤n4),表示是由第三段恒频波CW2得到的多普勒频率矩阵中第i个元素与第四段锯齿波FMCW2得到的差频频率矩阵中第j个元素进行计算得到的距离值。
实施例5:具有实施例1-4任一项相同的技术方案,更为具体的是:进行多目标的方位角的计算的具体步骤是:
(1)首先计算各过门限点对应的相位值。
通道1中,第一段恒频波CW1过门限点的位置矩阵ACW1=[a1,a2,…an1],根据如下计算方法计算对应点上的相位值,得到相位矩阵为ψACW1=[ψa1,ψa2,…ψan1],通道1中,第二段锯齿波FMCW1过门限点的位置矩阵BFMCW1=[b1,b2,…bn2],根据如下计算方法计算对应点上的差频频率值ψBFMCW1=[ψb1,ψb2,…ψbn2]。通道2中,第一段恒频波CW2过门限点的位置矩阵ECW1=[e1,e2,…en5],根据如下计算方法计算对应点上的相位值,得到相位矩阵为ψECW1=[ψe1,ψe2,…ψen5],通道2中,第二段锯齿波FMCW1过门限点的位置矩阵GFMCW1=[g1,g2,…gn6],根据如下计算方法计算对应点上的差频频率值ψGFMCW1=[ψg1,ψg2,…ψgn6]。
其中,计算相位的方法为,在各段波形计算FFT后,得到各自复数值的矩阵。根据复
数的特性求相位值的方法为,假设复数为c=a+j*b=cosθ+j*sinθ,则得
到复数的相位值
(2)计算相位差。
计算通道1的第一段恒频波CW1和通道2的第一段恒频波CW2之间的相位差,得到相位差矩阵为计算通道1的第二段锯齿波FMCW1和通道2的第二段锯齿波FMCW1之间的相位差,得到相位差矩阵为
具体相位差的计算公式为:1≤i≤
n1,1≤j≤5
(1≤i≤n1,1≤j≤n5)。
(3)计算方位角。
得到相位差矩阵后,根据公式方位角公式,其中,d=7.5mm为天线间距,λ=12.4mm。
计算通道1的第一段恒频波CW1和通道2的第一段恒频波CW2之间的方位角,得到方位角矩阵为计算通道1的第二段锯齿波FMCW1和通道2的第二段锯齿波FMCW1之间的方位角,得到相位差矩阵为
实施例6:具有与实施例1-5任一项相同的技术方案,更为具体的是:还包括查找真实距离模块,查找多目标的真实距离。优选的,查找多目标的真实距离的具体步骤是:
(1)由于第一段恒频波CW1获得的多普勒值与第二段锯齿波FMCW1获得的差频频率是一一对应的,也就是在相对距离矩阵的第一行中,仅有一个值是目标的距离真值,其他的距离值都是由于多普勒值与错误的差频频率值相匹配而产生的虚假目标的距离值。同理,对于第三段恒频波CW2获得的多普勒值和第四段锯齿波FMCW2获得的差频频率值也是一一对应的。
(2)由于真实目标的值是与周期无关,而虚假目标的值是与周期有关系的,也就是在相对距离矩阵RAB和RCD的矩阵中,真实目标的距离值是相等的,虚假目标的距离值是完全不同的。只要找到两个矩阵中距离相等的值,则为真实目标的距离值。
(3)由于相对距离矩阵RAB和相对距离矩阵RCD的维数可能存在不同,这是由于,在进行过门限检测的时候,各个波段过门限后所检测的目标点数存在一定的差异造成的。处理的原则是,以相对距离矩阵RAB和相对距离矩阵RCD的行数和列数之和为最小的矩阵为基准,进行距离真值的查找。
(4)在相对距离矩阵RAB和相对距离矩阵RCD中找到一个距离真值的时候,记录其真值所在矩阵的行值与列值,同时,将获得的真值所在矩阵中的行和列的值全部去除掉,依次进行全部真值的查找,直到作为基准的矩阵维数为0为止。
实施例7:具有与实施例1-6任一项相同的技术方案,更为具体的是:还包括速度值匹配模块,进行多目标的真实距离值与其对应的速度值的匹配和/或方位角匹配模块,根据多目标的真实距离值,进行多目标方位角的匹配。
实施例8:具有与实施例1-7任一项相同的技术方案,更为具体的是:
速度值匹配的具体步骤是:利用找到的真实距离值所在的行数值,作为相对速度矩阵VACW1=[va1,va2,…van1]或是VCCW1=[vc1,vc2,…vcn3]矩阵的列数值,找到真实目标距离值所对应的速度值,这样则完成真实目标的距离与速度的匹配;
方位角匹配的具体步骤是:在相对距离矩阵RAB中找到距离真值的时候,同时记录目标真值所在相对距离矩阵的行值与列值,在相对速度矩阵VACW1中找到速度真值的时候,记录相对速度矩阵的列值。利用找到的真实距离值所在的行值,对方位角矩阵θAE相同行中的所有方位角值和真实速度值所在的列值,对方位角矩阵θBG相同数值的行中的所有方位角值,查找相同方位角的值,查找到目标所对应的方向角后,记录该方向角所在矩阵的行值和列值,同时,将θAE和θBG获得的方向角所在矩阵中,对应的行和列的值全部去除掉,根据相对距离矩阵RAB中,真实目标所在的行数,依次进行目标对应方向角的查找,直到找出所有真实目标对应的方向角。
本实施例给出了一种可实现多目标检测的组合波形设计方案,同时给出了,可实现多目标检测的理论分析,对于在设计其他波形实现多目标检测时,提供了一种波形设计思路;
本实施例所设计变道辅助系统,可以实现对多目标的相对距离以及相对速度的检测,同时可以实现目标方向角的检测功能,实现对多目标空间方位上的区分。
本实施例给出了详细的信号处理过程,包括多目标相对速度的解算、相对距离解算,真实目标的相对距离以及相对速度的匹配方法、相位差的计算,方向角的解算,以及真实目标方向角的匹配方法等处理过程以及相关公式。该部分对于设计变道辅助系统人员提供了一种信号处理方法。
实施例9:一种多目标检测的组合波形汽车变道辅助系统的信号处理装置,所述波形包括第一段恒频波CW1、第二段锯齿波FMCW1、第三段恒频波CW2、第四段锯齿波FMCW2,该装置包括:
FFT计算模块,对各段波形,A/D采集到的IQ数据,进行FFT计算。
门限检测模块,将各段波形FFT变换后的复数模值做门限检测,输出过门限点位置。
多普勒频率值计算模块,计算得到恒频段的多普勒频率值。
差频频率值计算模块,计算得到锯齿波段的差频频率值。
相对速度值计算模块,根据得到的多普勒频率值,计算相对速度矩阵。
相对距离值计算模块,根据计算得到的多普勒频率矩阵与差频频率矩阵,计算得到相对距离矩阵。
方向角计算模块,进行多目标的方位角的计算。
查找真实距离模块,查找多目标的真实距离。
速度值匹配模块,进行多目标的真实距离值与其对应的速度值的匹配。
方位角匹配模块,根据多目标的真实距离值,进行多目标方位角的匹配。
实施例10:一种多目标检测的组合波形汽车变道辅助系统的信号处理方法,对应于实施例1-9中任一项装置技术方案的方法,所述波形包括第一段恒频波CW1、第二段锯齿波FMCW1、第三段恒频波CW2、第四段锯齿波FMCW2,该方法包括以下步骤:
S1.对各段波形,A/D采集到的IQ数据,进行FFT计算;
S2.将各段波形FFT变换后的复数模值做门限检测,输出过门限点位置;
S3.计算得到恒频段的多普勒频率值。
S4.计算得到锯齿波段的差频频率值。
S5.根据得到的多普勒频率值,计算相对速度矩阵。
S6.根据计算得到的多普勒频率矩阵与差频频率矩阵,计算得到相对距离矩阵。
S7.查找多目标的真实距离.
S8.进行多目标的真实距离值与其对应的速度值的匹配和/或S10.根据多目标的真实距离值,进行多目标方位角的匹配。
S9.进行多目标的方位角的计算。
实施例11:一种汽车变道辅助系统,安装有实施例1-10任一项技术方案所述的多目标检测的组合波形汽车变道辅助系统的信号处理装置。
本实施例所设计恒频波与锯齿波的组和波形,以及波形所涉及的参数选取不局限于本实施例公开的参数,本领域技术人员可以根据具体应用场景,选取不同的设计参数,或是对波形进行改进等;在本发明披露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。
Claims (6)
1.一种多目标检测的组合波形汽车变道辅助系统的信号处理装置,其特征在于:所述波形包括第一段恒频波CW1、第二段锯齿波FMCW1、第三段恒频波CW2、第四段锯齿波FMCW2,该装置包括:
FFT计算模块,对各段波形,A/D采集到的IQ数据,进行FFT计算;
门限检测模块,将各段波形FFT变换后的复数模值做门限检测,输出过门限点位置;
多普勒频率值计算模块,计算得到恒频段的多普勒频率值;
差频频率值计算模块,计算得到锯齿波段的差频频率值;相对速度值计算模块,根据得到的多普勒频率值,计算相对速度矩阵;
相对距离值计算模块,根据计算得到的多普勒频率矩阵与差频频率矩阵,计算得到相对距离矩阵;
方向角计算模块,进行多目标的方位角的计算;
相对距离值计算模块基于如下方式计算得到相对距离矩阵:计算距离公式为其中,T为每一段波形的作用时间,B为调频带宽,fd为多普勒频率值,fo为差频频率值;
根据通道1中,第一段恒频波CW1得到的多普勒频率矩阵FACW1=[fda1,fda2,…fdan1]和第二段锯齿波FMCW1得到的差频频率矩阵FBFMCW1=[fob1,fob2,…fobn2],将多普勒频率矩阵中所有元素和差频频率矩阵中的所有元素进行一一配对计算相对距离矩阵,计算得到相对距离矩阵为其中raibj,1≤i≤n1,1≤j≤n2,表示是由第一段恒频波CW1得到的多普勒频率矩阵中第i个元素与第二段锯齿波FMCW1得到的差频频率矩阵中第j个元素进行计算得到的距离值;对于第三段恒频波CW2得到的多普勒频率矩阵FCCW2=[fdc1,fdc2,…fdcn3]和第四段锯齿波FMCW2得到的差频频率矩阵FDFMCW2=[fod1,fod2,…fodn4],同样进行上述处理,最后得到相对距离矩阵为其中rcicj,1≤i≤n3,1≤j≤n4,表示是由第三段恒频波CW2得到的多普勒频率矩阵中第i个元素与第四段锯齿波FMCW2得到的差频频率矩阵中第j个元素进行计算得到的距离值。
2.如权利要求1所述的装置,其特征在于,还包括:查找真实距离模块,查找多目标的真实距离。
3.如权利要求2所述的装置,其特征在于,还包括:速度值匹配模块,进行多目标的真实距离值与其对应的速度值的匹配和/或方位角匹配模块,根据多目标的真实距离值,进行多目标方位角的匹配。
4.如权利要求1所述的装置,其特征在于,多普勒频率值计算模块基于如下方式计算得到恒频段的多普勒频率值:通道1中,第一段恒频波CW1过门限点的位置矩阵ACW1=[a1,a2,…an1],根据如下规则计算对应点上的多普勒频率值,得到多普勒频率矩阵为FACW1=[fda1,fda2,…fdan1];对于通道1中,第三段恒频波CW2过门限点的位置矩阵CCW2=[c1,c2,…cn3],根据如下规则计算对应点上的多普勒频率值,得到多普勒频率矩阵为FCCW2=[fdc1,fdc2,…fdcn3];
该规则为,若点数为1≤xi≤128,1≤i≤n,判断目标靠近,其对应点上的多普勒频率1≤i≤n;若点数为128<xi≤256,1≤i≤n,判断目标远离,其对应点上的多普勒频率1<i≤n。
5.如权利要求1所述的装置,其特征在于,差频频率值计算模块基于如下方式计算得到锯齿波段的差频频率值:通道1中,第二段锯齿波FMCW1过门限点的位置矩阵BFMCW1=[b1,b2,…bn2],根据如下规则计算对应点上的差频频率矩阵FBFMCW1=[fob1,fob2,…fobn2],对于通道1中,第四段锯齿波FMCW2过门限点的位置矩阵DFMCW2=[d1,d2,…dn4],根据如下规则计算对应点上的差频频率矩阵FDFMCW2=[fod1,fod2,…fodn4];
该规则为,即若点数为1≤yj≤128,1≤j≤n,其对应点上的差频频率值1≤j≤n;若点数为128<yj≤256,1≤j≤n,其对应点上的差频频率值1<j≤n。
6.如权利要求1所述的装置,其特征在于,方向角计算模块基于如下方式计算多目标的方位角:
(1)首先计算各过门限点对应的相位值
通道1中,第一段恒频波CW1过门限点的位置矩阵ACW1=[a1,a2,…an1],根据如下计算方法计算对应点上的相位值,得到相位矩阵为ψACW1=[ψa1,ψa2,…ψan1];
通道1中,第二段锯齿波FMCW1过门限点的位置矩阵BFMCW1=[b1,b2,…bn2],根据如下计算方法计算对应点上的差频频率值ψBFMCW1=[ψb1,ψb2,…ψbn2];
通道2中,第一段恒频波CW1过门限点的位置矩阵ECW1=[e1,e2,…en5],根据如下计算方法计算对应点上的相位值,得到相位矩阵为ψECW1=[ψe1,ψe2,…ψen5];通道2中,第二段锯齿波FMCW1过门限点的位置矩阵GFMCW1=[g1,g2,…gn6],根据如下计算方法计算对应点上的差频频率值ψGFMCW1=[ψg1,ψg2,…ψgn6];
其中,计算相位值的方法为,在各段波形计算FFT后,得到各自复数值的矩阵,根据复数的特性求相位值的方法为,假设:复数为c=a+j*b=cosθ+j*sinθ,则得到复数的相位值
(2)计算相位差
计算通道1的第一段恒频波CW1和通道2的第一段恒频波CW2之间的相位差,得到相位差矩阵为计算通道1的第二段锯齿波FMCW1和通道2的第二段锯齿波FMCW1之间的相位差,得到相位差矩阵为
具体相位差的计算公式为: 1≤i≤n1,1≤j≤n5;
(3)计算方位角
得到相位差矩阵后,根据公式方位角公式,其中,d为天线间距,λ为波长;
计算通道1的第一段恒频波CW1和通道2的第一段恒频波CW1之间的方位角,得到方位角矩阵为计算通道1的第二段锯齿波FMCW1和通道2的第二段锯齿波FMCW1之间的方位角,得到相位差矩阵为
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511005838.0A CN105629235B (zh) | 2015-12-29 | 2015-12-29 | 多目标检测的组合波形汽车变道辅助系统的信号处理装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511005838.0A CN105629235B (zh) | 2015-12-29 | 2015-12-29 | 多目标检测的组合波形汽车变道辅助系统的信号处理装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105629235A CN105629235A (zh) | 2016-06-01 |
CN105629235B true CN105629235B (zh) | 2019-03-15 |
Family
ID=56044351
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201511005838.0A Active CN105629235B (zh) | 2015-12-29 | 2015-12-29 | 多目标检测的组合波形汽车变道辅助系统的信号处理装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105629235B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107783128B (zh) * | 2016-08-25 | 2021-02-19 | 大连楼兰科技股份有限公司 | 基于毫米波雷达的固定翼无人机多目标防撞系统 |
CN107783101B (zh) * | 2016-08-25 | 2021-04-09 | 大连楼兰科技股份有限公司 | 基于组合波形的旋翼无人机短距离防撞系统信号处理方法及装置 |
CN107783126B (zh) * | 2016-08-25 | 2021-03-30 | 大连楼兰科技股份有限公司 | 基于组合波形的自动驾驶汽车防撞雷达系统信号处理方法及装置 |
CN107783127B (zh) * | 2016-08-25 | 2021-04-09 | 大连楼兰科技股份有限公司 | 旋翼无人机防撞毫米波雷达信号处理方法 |
CN108415010B (zh) * | 2018-04-28 | 2022-04-19 | 西安电子科技大学 | 一种基于梯形lfmcw调制的雷达多目标检测方法 |
CN111537973A (zh) * | 2020-05-15 | 2020-08-14 | 山东汇贸电子口岸有限公司 | 对雷达接收到的回波进行处理的方法及装置 |
CN112381136B (zh) * | 2020-11-12 | 2022-08-19 | 深兰智能科技(上海)有限公司 | 目标检测方法和装置 |
CN112319505B (zh) * | 2020-11-26 | 2022-05-03 | 东风商用车有限公司 | 自动驾驶换道判断方法、装置、设备及存储介质 |
CN112882006B (zh) * | 2021-01-20 | 2024-02-13 | 南京理工大学 | 基于复合调频连续波的毫米波近程目标探测方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090315761A1 (en) * | 2006-07-13 | 2009-12-24 | Thomas Walter | FMCW Radar Sensor |
CN102944877A (zh) * | 2012-11-23 | 2013-02-27 | 北京理工大学 | 一种利用梯形调制波提高调频连续波测距精度的方法 |
CN102998671A (zh) * | 2012-11-23 | 2013-03-27 | 北京理工大学 | 一种利用梯形调制波判断调频连续波测距时有无强信号干扰的方法 |
CN103257346A (zh) * | 2013-05-15 | 2013-08-21 | 桂林电子科技大学 | 一种汽车防撞雷达多目标探测方法与系统 |
US20140049420A1 (en) * | 2012-08-14 | 2014-02-20 | Jenoptik Robot Gmbh | Method for Classifying Moving Vehicles |
CN103852756A (zh) * | 2012-11-30 | 2014-06-11 | 中国科学院沈阳自动化研究所 | 一种连续波雷达目标检测跟踪方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3460453B2 (ja) * | 1995-12-11 | 2003-10-27 | 株式会社デンソー | Fmcwレーダ装置 |
-
2015
- 2015-12-29 CN CN201511005838.0A patent/CN105629235B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090315761A1 (en) * | 2006-07-13 | 2009-12-24 | Thomas Walter | FMCW Radar Sensor |
US20140049420A1 (en) * | 2012-08-14 | 2014-02-20 | Jenoptik Robot Gmbh | Method for Classifying Moving Vehicles |
CN102944877A (zh) * | 2012-11-23 | 2013-02-27 | 北京理工大学 | 一种利用梯形调制波提高调频连续波测距精度的方法 |
CN102998671A (zh) * | 2012-11-23 | 2013-03-27 | 北京理工大学 | 一种利用梯形调制波判断调频连续波测距时有无强信号干扰的方法 |
CN103852756A (zh) * | 2012-11-30 | 2014-06-11 | 中国科学院沈阳自动化研究所 | 一种连续波雷达目标检测跟踪方法 |
CN103257346A (zh) * | 2013-05-15 | 2013-08-21 | 桂林电子科技大学 | 一种汽车防撞雷达多目标探测方法与系统 |
Non-Patent Citations (3)
Title |
---|
"基于梯形波FMCW雷达差拍信号配对的研究";薛效龙等;《全国微波毫米波会议》;20131231;第1474-1476页 * |
"汽车防撞毫米波雷达信号处理";周刚等;《计算机测量与控制》;20111231;第19卷(第7期);第1656-1658页 * |
薛效龙等."基于梯形波FMCW雷达差拍信号配对的研究".《全国微波毫米波会议》.2013,第1474-1476页. * |
Also Published As
Publication number | Publication date |
---|---|
CN105629235A (zh) | 2016-06-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105629235B (zh) | 多目标检测的组合波形汽车变道辅助系统的信号处理装置 | |
CN105629211B (zh) | 多目标检测的组合波形汽车变道辅助系统的信号处理方法和汽车变道辅助系统 | |
US10180494B2 (en) | Method for angle estimation and radar sensor for motor vehicles | |
JP7081046B2 (ja) | レーダセンサにおける角度測定誤差を検出する方法 | |
CN103389493B (zh) | 基于波束选择方法的多波束单脉冲测角方法 | |
CN107209253A (zh) | 具有水平数字波束形成和通过在彼此错开的发射机的情况下进行相位比较引起的垂直对象测量的成像雷达传感器 | |
CN102298142B (zh) | 雷达测距模糊条件下多微弱目标ht-stc-tbd检测方法 | |
CN106405556A (zh) | 车辆目标信息探测识别系统及其信号处理方法 | |
CN105487060B (zh) | 一种双通道四斜率调制的多目标提取方法 | |
CN100365429C (zh) | 一种合成孔径雷达运动目标成像方法 | |
EP2437079A1 (en) | Radar system | |
CN102798863A (zh) | 基于汽车防撞雷达的道路中央隔离带检测方法 | |
CN109598946B (zh) | 一种基于雷达体制的多车道测速方法 | |
CN104898119B (zh) | 一种基于相关函数的动目标参数估计方法 | |
CN107831489B (zh) | 一种连续波4d雷达及其测量多目标方法 | |
CN105005049A (zh) | 基于微波探测技术的车载防撞报警器及车距测量示警方法 | |
CN102608598A (zh) | 基于子空间投影的实孔径前视成像方法 | |
CN111929685A (zh) | 基于线性调频连续波的雷达检测方法及装置、雷达装置和机动车 | |
Askeland et al. | Tracking with a high-resolution 2D spectral estimation based automotive radar | |
Roldan et al. | See further than cfar: a data-driven radar detector trained by lidar | |
Greco | Automotive radar | |
Rizik et al. | Single target recognition using a low-cost FMCW radar based on spectrum analysis | |
CN102062851A (zh) | 基于改进l阵的星载宽带多目标的测向方法 | |
CN107783130B (zh) | 基于组合波形的无人驾驶汽车复杂环境防碰撞系统信号处理方法 | |
CN108549063B (zh) | 一种星载合成孔径雷达图像中运动目标航迹向速度快速估计方法 |
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 |