CN105242272B - 基于自回归时间序列模型的车载毫米波防撞雷达恒虚警检测方法 - Google Patents
基于自回归时间序列模型的车载毫米波防撞雷达恒虚警检测方法 Download PDFInfo
- Publication number
- CN105242272B CN105242272B CN201510634038.9A CN201510634038A CN105242272B CN 105242272 B CN105242272 B CN 105242272B CN 201510634038 A CN201510634038 A CN 201510634038A CN 105242272 B CN105242272 B CN 105242272B
- Authority
- CN
- China
- Prior art keywords
- time series
- sequence
- time
- cfar detection
- auto
- 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
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
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Abstract
基于自回归时间序列模型的车载毫米波防撞雷达恒虚警检测方法,属于恒虚警检测领域,用于解决现有恒虚警检测方法无法保留毫米波雷达系统的实时性的问题,技术要点是:S1.对毫米波雷达信号收到的回波进行快速傅里叶变换;S2.对当前门限时间序列及该序列的前三个时间周期的时间序列分别进行数据更新;S3.求取当前门限时间序列的前三个时间周期的时间序列所对应的新的时间序列在对应点上的均值;S4.应用自回归模型,得到当前周期的阈值门限。有益效果:保证了恒虚警检测过程中系统的实时性。
Description
技术领域
本发明属于恒虚警检测领域,涉及一种车载毫米波防撞雷达恒虚警检测方法。
背景技术
随着现代化水平的高度发展,交通问题已经成为各国亟需解决的问题,毫米波汽车防撞雷达就是解决这类问题最好的途径.毫米波汽车防撞雷达最大优点就是可以全天候工作,它广泛应用于汽车和智能交通领域.毫米波汽车防撞雷达工作于杂波环境时,检测器的门限设置必须自适应于杂波功率水平的变化,使杂波引起的目标检测即虚警保持在一个较低的可接受水平上,从而解决传统雷达门限是一个常数存在的缺点。
雷达恒虚警(constant false alarm rate,CFAR)检测的实际检测性能与背景杂噪的分布和强度有很大关系。雷达自动检测和跟踪中的虚警问题是每个雷达系统和设计人员不可避免的重要问题之一。二十多年来,雷达自动检测与CFAR处理技术逐渐发展成为国际雷达信号处理界的一大热门研究领域和关键性课题。CFAR技术是雷达自动检测系统中控制虚警率的最重要手段,它是雷达自动检测过程中起到极其重要的作用。现在,CFAR研究已经出现了多个研究方向。根据模拟杂波背景所使用的咋波分布模型分为:瑞利分布、韦布尔分布、对数正态分布、K分布和莱斯分布模型中的CFAR研究;按照数据处理方法分为:参量和非参量CFAR技术;按处理所在的数域分为:时域和频域CFAR研究方向;根据数据的形式分为:标量和向量CFAR技术;根据信号的相关程度分为:相关和不相关信号几部分相关信号的CFAR方法。此外,还可分为单参数和多参数CFAR技术,单传感器和多传感器分布式CFAR技术,以及其他的一些研究方法。
CFAR处理的主要过程就是形成杂波平均功率水平估计,相应的估计方法有两大类:一类是空域方法,即利用与检测单元在空间上相邻的参考单元采样值,这类方法适用于在空间上平稳的杂波;另一类是时域的方法,即杂波图,它利用检测单元以往多次扫描测量值形成杂波平均功率水平估计,并依靠新的测量值迭代更行。
杂波图存储每一个单元背景的估计值,每一个值依靠新的和以前若干次的扫描测量值来进行迭代更行,并把它作为当前的杂波背景强度估计值。当空域杂波强烈变化剧烈时,若采用传统的均值类恒虚警方案只能采用很少的参考单元,因而恒虚警损失很大,而且虚警率不易保持恒定。一般情况下早搏虽然在距离和方位上的变化十分剧烈,但同一距离单元的杂波强度随时间的变化使缓慢的,因而可以采用“时间单元”恒虚警处理方法,在时间上对以往各次雷达回波的测量值进行迭代处理,依旧是所谓的CFAR方法。
Nitzberg等人首先提出和分析了这种方法,在文献中成为Nitzberg杂波图方法,这一处理过程利用了雷达杂波环境时域的相对平稳性。因此,杂波图处理不受注入地杂波在空域非平稳性的影响。Levanon教授给出了另外一种计算Nitzberg杂波图虚警率的公式,可以较快的速度收敛。为了增强杂波图的鲁棒性,Lops M提出混合杂波图技术(CM/L-CFAR),即将空域处理与时域处理结合起来。它将几个雷达分辨率单元的回波信号组合成一个杂波图单元,对杂波图单元内的回波样本基于L滤波器进行空域处理作为当前早不单元的输入信号;在对每个杂波图单元的输入信号对以往各次的结果进行迭代,得到杂波图对检测位置处杂波强度的估计以检测目标。Hamadouche M研究了Nitzberg杂波图技术在韦布尔背景中的检测性能,分析指出Nitzberg杂波图处理在拖尾严重的非高斯杂波时其目标检测能力很差。孟祥伟采用排序数据变量(ODV)方法来替换CM/L-CFAR方法中进行固定删除处理的L滤波器,删除了CM/ODV-CFAR方法,由于它对杂波图单元内的杂波样本进行自适应删除而获得了一定程度上性能的改善。
以上很多雷达恒虚警检测方法,在理论研究中具有重要的意义,但是却牺牲了实际车载毫米波雷达系统的实时性这一系统的重要指标。
发明内容
为了解决上述技术问题,本发明旨在得到一种保留毫米波雷达系统的实时性而检测恒虚警的方法,以提高现有恒虚警检测方法的准确性和实时性。
为了实现上述目的,本发明采用的技术方案是:一种基于自回归时间序列模型的车载毫米波防撞雷达恒虚警检测方法,包括如下步骤:
S1.对毫米波雷达信号收到的回波,在NFFT点进行快速傅里叶变换,得到各个时间周期内NFFT点的快速傅里叶变换后的值,作为各个时间周期的时间序列;
S2.对当前门限时间序列及该序列的前三个时间周期的时间序列分别进行数据更新,分别获得当前门限时间序列及该序列的前三个时间周期的时间序列所对应的新的时间序列;
S3.求取当前门限时间序列的前三个时间周期的时间序列所对应的新的时间序列在对应点上的均值;
S4.使用步骤S3中所述的均值,应用自回归模型,得到当前周期的阈值门限。
作为技术方案的补充,所述方法还包括下述步骤:得到最终的门限值,其为步骤S4中得到的当前周期的阈值门限与门限调整因子的乘积。
作为技术方案的进一步补充,步骤S2中对时间序列更新的方法为:
S2.1.将当前时间序列进行升序排序,删除百分之M个序列的最大值Ymax_M%(t)或者直接删除M个最大值,同时删除百分之N个序列的最小值Ymin_N%(t)或者直接删除N个最小指,M表示排序后序列后面的M个最大值,N表示排序后序列前面的N个最小值;
S2.2.对余下的时间序列的所有值求取平均值YQ(t);
S2.3.用求取的平均值YQ(t)替换掉所删除的最大值和最小值,完成当前时间序列的更新。
有益效果:本发明所述的车载毫米波防撞雷达恒虚警检测方法,采用了时间域上的恒虚警检测方法,该方法对于处理高速等路面或空域非平稳的地杂波有很好的特性,同时具有较佳的准确性和实时性等特性,使得该算法在工程应用中具有较好的适用价值,且使用该方法,保证了恒虚警检测过程中系统的实时性。
本发明所述的自回归模型算法,将当前周期的前三个周期的均值作为历史数据,而不是采用上一个周期数据作为历史数据,增加了系统的稳定性和可靠性。
附图说明
图1为本发明实施例中所述方法的流程示意图。
具体实施方式
为了更为准确的阐述本发明创造,结合具体实施方式对本发明进一步说明。
实施例1:
一种基于自回归时间序列模型的车载毫米波防撞雷达恒虚警检测方法,包括如下步骤:
S1.对毫米波雷达信号收到的回波,在NFFT点进行快速傅里叶变换,得到各个时间周期内NFFT点的快速傅里叶变换后的值,作为各个时间周期的时间序列;
S2.对当前门限时间序列Yn(t)(t≥4)及该序列的前三个时间周期的时间序列Yn(t-1),Yn(t-2),Yn(t-3)分别进行数据更新,分别获得当前门限时间序列及该序列的前三个时间周期的时间序列所对应的新的时间序列Y′n(t),Y′n(t-1),Y′n(t-2),Y′n(t-3)。
S3.求取当前门限时间序列的前三个时间周期的时间序列所对应的新的时间序列在对应点上的均值:
(n∈[1,NFFT]),NFFT表示快速傅里叶变换后的点数;
S4.使用步骤S3中所述的均值,应用自回归模型,得到当前周期的阈值门限,获得阈值门限是一种自适应门限,具有一定的恒虚警检测概率,阈值门限:
Yn_yz=w·Y′n(t)+(1-w)·Y′n_Q(t),w为遗忘因子。
实施例2:
与实施例1的区别在于,所述方法还包括下述步骤:S5.得到最终的门限值Yn_mx:
Yn_mx=Yn_yz×μ,μ为门限调整因子。
由于在实施例1中,门限计算出来后可能会只比原始数据大不了多少,很容易存在一定的虚警率,所以将门限通过门限调整因子提高一些,这样可以有效的降低虚警率。
实施例3:
上述两个实施例中的步骤S2中对时间序列更新的方法为:
S2.1.将当前时间序列进行升序排序,删除百分之M个序列的最大值Ymax_M%(t)或者直接删除M个最大值,同时删除百分之N个序列的最小值Ymin_N%(t)或者直接删除N个最小指,M表示排序后序列后面的M个最大值,N表示排序后序列前面的N个最小值;
S2.2.对余下的时间序列的所有值求取平均值YQ(t);
S2.3.用求取的平均值YQ(t)替换掉所删除的最大值和最小值,完成当前时间序列的更新。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明披露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。
Claims (3)
1.一种基于自回归时间序列模型的车载毫米波防撞雷达恒虚警检测方法,其特征在于,包括如下步骤:
S1.对毫米波雷达信号收到的回波,在NFFT点进行快速傅里叶变换,得到各个时间周期内NFFT点的快速傅里叶变换后的值,作为各个时间周期的时间序列;
S2.对当前门限时间序列及该序列的前三个时间周期的时间序列分别进行数据更新,分别获得当前门限时间序列及该序列的前三个时间周期的时间序列所对应的新的时间序列;
S3.求取当前门限时间序列的前三个时间周期的时间序列所对应的新的时间序列在对应点上的均值;
S4.使用步骤S3中所述的均值,应用自回归模型,得到当前周期的阈值门限。
2.如权利要求1所述的基于自回归时间序列模型的车载毫米波防撞雷达恒虚警检测方法,其特征在于,还包括下述步骤:
S5.得到最终的门限值,其为步骤S4中得到的当前周期的阈值门限与门限调整因子的乘积。
3.如权利要求1所述的基于自回归时间序列模型的车载毫米波防撞雷达恒虚警检测方法,其特征在于,步骤S2中对时间序列更新的方法为:
S2.1.将当前时间序列进行升序排序,删除百分之M个序列的最大值Ymax_M%(t)或者直接删除M个最大值,同时删除百分之N个序列的最小值Ymin_N%(t)或者直接删除N个最小值,M表示排序后序列后面的M个最大值,N表示排序后序列前面的N个最小值;
S2.2.对余下的时间序列的所有值求取平均值YQ(t);
S2.3.用求取的平均值YQ(t)替换掉所删除的最大值和最小值,完成当前时间序列的更新。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510634038.9A CN105242272B (zh) | 2015-09-29 | 2015-09-29 | 基于自回归时间序列模型的车载毫米波防撞雷达恒虚警检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510634038.9A CN105242272B (zh) | 2015-09-29 | 2015-09-29 | 基于自回归时间序列模型的车载毫米波防撞雷达恒虚警检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105242272A CN105242272A (zh) | 2016-01-13 |
CN105242272B true CN105242272B (zh) | 2017-10-27 |
Family
ID=55039975
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510634038.9A Active CN105242272B (zh) | 2015-09-29 | 2015-09-29 | 基于自回归时间序列模型的车载毫米波防撞雷达恒虚警检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105242272B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107783112A (zh) * | 2016-08-25 | 2018-03-09 | 大连楼兰科技股份有限公司 | 基于组合波形的旋翼无人机复杂环境防碰撞系统信号处理方法 |
CN107783127B (zh) * | 2016-08-25 | 2021-04-09 | 大连楼兰科技股份有限公司 | 旋翼无人机防撞毫米波雷达信号处理方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104406792A (zh) * | 2014-09-18 | 2015-03-11 | 内蒙航天动力机械测试所 | 一种固体火箭发动机故障诊断方法 |
CN104914433A (zh) * | 2015-06-15 | 2015-09-16 | 南京航空航天大学 | 一种基于链表排序的os-cfar多目标提取的实现方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4881209B2 (ja) * | 2007-04-02 | 2012-02-22 | 株式会社東芝 | 目標検出装置 |
KR101109150B1 (ko) * | 2010-07-26 | 2012-02-24 | 재단법인대구경북과학기술원 | 계산량을 향상시킨 순차통계 일정 오경보율 검파방법 |
US8976059B2 (en) * | 2012-12-21 | 2015-03-10 | Raytheon Canada Limited | Identification and removal of a false detection in a radar system |
-
2015
- 2015-09-29 CN CN201510634038.9A patent/CN105242272B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104406792A (zh) * | 2014-09-18 | 2015-03-11 | 内蒙航天动力机械测试所 | 一种固体火箭发动机故障诊断方法 |
CN104914433A (zh) * | 2015-06-15 | 2015-09-16 | 南京航空航天大学 | 一种基于链表排序的os-cfar多目标提取的实现方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105242272A (zh) | 2016-01-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105137498B (zh) | 一种基于特征融合的地下目标探测识别系统及方法 | |
CN108802722B (zh) | 一种基于虚拟谱的弱目标检测前跟踪方法 | |
CN111045008B (zh) | 基于展宽计算的车载毫米波雷达目标识别方法 | |
CN104931934B (zh) | 一种基于pam聚类分析的雷达点迹凝聚方法 | |
CN103217673B (zh) | 一种非均匀Weibull杂波背景下的CFAR检测方法 | |
CN107765228B (zh) | 一种基于区域相似性的在线雷达目标检测方法 | |
CN109100696B (zh) | 基于点状杂波图的慢速运动目标消除方法 | |
CN110837079B (zh) | 一种基于雷达的目标检测方法及装置 | |
CN104200471A (zh) | 基于自适应权值图像融合的sar图像变化检测方法 | |
CN101482969B (zh) | 基于同质点计算的sar图像去斑方法 | |
CN108196241B (zh) | 一种基于Hough变换的高速动目标速度估计方法 | |
CN111008585A (zh) | 基于自适应分层高分辨sar图像的舰船目标检测方法 | |
CN106468771A (zh) | 一种低可观测高杂波条件下的多目标检测与跟踪方法 | |
CN105372636A (zh) | 一种基于自适应遗忘因子的杂波图更新方法及装置 | |
CN105044686A (zh) | 一种雷达密集假目标干扰抑制方法 | |
CN112083390B (zh) | 一种时空关联下雷达杂波的动态抑制方法及装置 | |
CN109932695A (zh) | 一种提高目标物识别速度的方法及装置 | |
CN111999726B (zh) | 一种基于毫米波雷达的人员定位方法 | |
CN105913074A (zh) | 基于幅度与径向速度联合的sar图像动目标聚类方法 | |
CN111562570A (zh) | 基于毫米波雷达的面向自动驾驶的车辆感知方法 | |
CN105242272B (zh) | 基于自回归时间序列模型的车载毫米波防撞雷达恒虚警检测方法 | |
CN106056115B (zh) | 一种非均匀背景下的红外小目标检测方法 | |
CN106772357B (zh) | 信噪比未知条件下的ai-phd滤波器多目标跟踪方法 | |
CN103576131A (zh) | 距离扩展目标内禀模态特征能量智能融合检测方法 | |
CN104931951A (zh) | 一种基于Hough变换域杂波图的重杂波区小目标的检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C41 | Transfer of patent application or patent right or utility model | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20160122 Address after: 116023 Dalian province high tech park, Hui Park, No. 7, take off the garden of the 2 phase of the 11 floor of the Applicant after: Dalian Loulan Technology Co., Ltd. Address before: 116023 Dalian province high tech park, Hui Park, No. 7, take off the garden of the 2 phase of the 11 floor of the Applicant before: XI'AN ZHIJI ANTENNA TECHNOLOGY CO., LTD. |
|
GR01 | Patent grant | ||
GR01 | Patent grant |