CN102798863B - 基于汽车防撞雷达的道路中央隔离带检测方法 - Google Patents

基于汽车防撞雷达的道路中央隔离带检测方法 Download PDF

Info

Publication number
CN102798863B
CN102798863B CN201210244436.6A CN201210244436A CN102798863B CN 102798863 B CN102798863 B CN 102798863B CN 201210244436 A CN201210244436 A CN 201210244436A CN 102798863 B CN102798863 B CN 102798863B
Authority
CN
China
Prior art keywords
dimension
target vehicle
spectrum
car
unit
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
Application number
CN201210244436.6A
Other languages
English (en)
Other versions
CN102798863A (zh
Inventor
宋骊平
张昱
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xidian University
Original Assignee
Xidian University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xidian University filed Critical Xidian University
Priority to CN201210244436.6A priority Critical patent/CN102798863B/zh
Publication of CN102798863A publication Critical patent/CN102798863A/zh
Application granted granted Critical
Publication of CN102798863B publication Critical patent/CN102798863B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种基于汽车防撞雷达的道路中央隔离带检测方法,其实现步骤为:(1)发射信号;(2)接收回波信号;(3)产生中频信号;(4)二维信号处理;(5)二维恒虚警检测;(6)确定动目标车辆位置和速度;(7)一维恒虚警检测;(8)确定静目标车辆和中央隔离带位置;(9)比较目标车辆与中央隔离带之间的位置关系;(10)防撞预警分析处理。本发明既可以较好的检测动目标车辆的位置和速度,又可以较好的检测静目标车辆和道路中央隔离带的位置,从而鉴别往返车道车辆,兼顾检测结果中的动目标漏检信息和预警***中的虚警信息,具有较好的实时性和较高的检测精度,更好的辅助驾驶员,保障行驶安全、防止汽车碰撞,可应用于高速公路环境中的汽车防撞、导航以及自动驾驶等领域。

Description

基于汽车防撞雷达的道路中央隔离带检测方法
技术领域
本发明属于通信技术领域,更进一步涉及雷达信号处理技术领域中的基于汽车防撞雷达的道路中央隔离带检测方法。本发明融入于传统汽车防撞雷达***中,能够在设置中央隔离带的高速公路上更好的辅助驾驶员,保障行驶安全、防止汽车碰撞。
背景技术
在汽车防撞雷达预警方法的研究中,常见的汽车防撞雷达是利用回波信号和发射信号进行相干混频,得到包含目标信息的中频信号,然后对中频信号进行频谱分析得到目标的距离、速度和方位角,进行简单的预警分析,但该方法存在的不足是,无法对高速公路中央隔离带进行检测,所以往往会将反向行驶于隔离带另一侧的汽车识别为危险目标,这就容易在雷达目标检测中产生较高的虚警率。
阮树成在其专利申请“毫米波时分线性调频多目标检测汽车防撞雷达”(专利申请号:200810120605.9,公开号:CN101354438)中提出了一种在道路环境中同时测定多个目标车辆的距离、速度和方位角的毫米波时分线性调频多目标检测汽车防撞雷达方法。该方法应用于毫米波段,运用时分线性调频的方法,能够同时测定多个目标车辆,然后通过防撞预警***对距离自车最近的目标车辆进行防撞预警分析,若判断此目标车辆为危险目标,则对驾驶员发出报警,保证车辆安全行驶,但是该方法存在的不足是,防撞预警***只对距离自车最近的目标车辆进行预警分析,而没有计算分析自车与其它目标车辆之间的理论安全距离,忽视了其它车辆对自车的潜在威胁,并且会将中央隔离带另一侧反向车道的车辆检测为预警目标,导致预警***发出错误报警。
侯志在文献“复杂探测背景下的LFMCW雷达动目标二维检测”(西安电子科技大学学报,2011,38(4):1001-2400.)中提出了一种LFMCW雷达动目标二维检测方法。该方法使用LFMCW雷达二维信号处理方法和二维恒虚警检测方法,对复杂探测背景中的动目标车辆和静止大目标角反射架进行检测识别,能够检测淹没在一维距离谱中的运动目标,降低了复杂探测背景对动目标检测的影响,但是该方法存在的不足是,它只对静止大目标进行检测,没有对道路环境中的连续静止目标即中央隔离带进行检测识别,将这种方法应用于防撞雷达***中,不能完整的对道路环境中存在的目标进行检测,导致预警***出现较高的虚警率。
发明内容
本发明针对上述现有技术存在的不足,提出了一种基于汽车防撞雷达的道路中央隔离带检测方法。本发明既可以较好的检测高速公路上车辆的行驶状态,同时又可以兼顾中央隔离带的检测识别,既保证了车辆的安全行驶,又降低了预警分析***的虚警率。
本发明实现上述目的的思路是:将雷达接收的回波信号和发射信号相干混频得到中频信号后,先对中频信号进行二维信号处理得到距离速度二维频谱,然后对多个雷达扫描周期的二维频谱幅值积累结果作二维恒虚警检测,分离动目标车辆频谱与静止目标频谱,确定动目标车辆的位置和速度,再对静止目标频谱做一维恒虚警检测,进一步将静目标车辆和中央隔离带分离开,确定静目标车辆和中央隔离带的位置,实现对中央隔离带的检测识别,最后通过目标车辆与中央隔离带之间的位置关系,判断目标车辆是否处于同向车道,对处于同向车道的目标车辆作防撞预警分析处理,提高了防撞预警***的预警效率,有效的保证了车辆的安全行驶。
本发明的步骤包括如下:
(1)由自车防撞雷达阵列发射天线,连续发射多周期对称三角波调制的线性调频连续波信号;
(2)接收回波信号
2a)在一个雷达扫描周期内,防撞雷达的接收天线,持续接收每一时刻连续波信号经由前方道路环境中车辆、中央隔离带以及其它障碍物反射回来的多周期回波信号;
2b)对雷达扫描周期内每一时刻由相邻阵列天线单元接收的不同相位回波信号进行比相,测出它们之间的相位差;
2c)根据来波方向与相位差之间的换算公式得到每一时刻回波信号的来波方向;
(3)产生中频信号
3a)将步骤(1)发射的连续波信号输入中频信号处理***中,作为混频器的本振信号;
3b)将步骤2a)接收的回波信号输入混频器中与本振信号相乘,得到一个混频信号;
3c)将混频信号通过低通滤波器产生一个中频信号;
(4)二维信号处理
4a)将这个中频信号作为中频信号处理单元;
4b)根据一维快速傅里叶变换公式,对中频信号处理单元中的中频信号作一维快速傅里叶变换,得到包含距离信息的一维频谱,此一维频谱包含上扫频段一维频谱和下扫频段一维频谱;
4c)根据离散快速傅里叶变换公式,对一维频谱作离散快速傅里叶变换,得到包含距离和相对速度信息的二维频谱,此二维频谱包含上扫频段二维频谱和下扫频段二维频谱;
4d)根据多普勒频率配对公式,对上扫频段二维频谱和下扫频段二维频谱,作消除配对处理,得到距离-相对速度二维频谱;
4e)根据速度与相对速度之间的变换公式,对距离-相对速度二维频谱进行变换,得到距离-速度二维频谱;
4f)对步骤3b)至4e)重复执行10次,将每次得到的距离-速度二维频谱幅值相加积累,得到高信噪比二维频谱,将其幅值以二维矩阵的形式记录,此二维矩阵由M1行N1列元素组成;
(5)二维恒虚警检测
5a)将此M1行N1列元素组成的二维矩阵输入二维恒虚警检测***;
5b)选取二维矩阵中的第i行第j列矩阵元素,令i等于1,j等于1;
5c)以第i行第j列矩阵元素作为二维被检测单元,分别提取以二维被检测单元为中心的相邻行和列中的4个矩阵元素和6个矩阵元素,组成二维检测参考单元,根据二维恒虚警检测参考单元的均值公式,求得二维检测参考单元的均值;
5d)将二维检测参考单元的均值乘以阈值因子得到二维恒虚警检测门限;
5e)将二维被检测单元的幅值与二维恒虚警检测门限作比较,当幅值大于检测门限时,记录此幅值,否则记为0;
5f)若i的值不等于M1,则将i的值加1后,转到步骤5c);否则,转到步骤5g);
5g)若j的值不等于N1,则将j的值加1后,转到步骤5c);否则,转到步骤5h);
5h)将步骤5e)记录的全部幅值组成静止频谱和动目标频谱,将静止频谱幅值以一维矩阵的形式记录,此一维矩阵由L1个元素组成;
(6)确定动目标车辆位置和速度
6a)将动目标频谱最大幅值对应的横坐标值,作为自车与动目标车辆之间的距离,对应的纵坐标值作为动目标车辆的速度;
6b)将与动目标频谱最大幅值对应的步骤2c)中同一时刻接收到的回波信号的来波方向作为动目标车辆的方位角;
6c)以自车为坐标原点,根据三角变换公式,计算得到动目标车辆的位置;
(7)一维恒虚警检测
7a)将步骤5h)中L1个元素组成的一维矩阵输入一维恒虚警检测***;
7b)选取一维矩阵中的第k个矩阵元素,令k等于1;
7c)以第k个矩阵元素作为一维被检测单元,分别提取以一维被检测单元为中心的相邻的4个矩阵元素,组成一维检测参考单元,根据一维恒虚警检测参考单元的均值公式,求得一维检测参考单元的均值;
7d)将一维检测参考单元的均值乘以阈值因子得到一维恒虚警检测门限;
7e)将一维被检测单元的幅值与一维恒虚警检测门限作比较,当幅值大于检测门限时,记录此幅值,否则记为0;
7f)若k的值不等于L1,则将k的值加1后,转到步骤7c);否则,转到步骤7g);
7g)将步骤7e)记录的全部幅值组成静目标频谱;
7h)将步骤7e)记为0的每一个被检测单元的幅值作为一个中央隔离带单元频谱幅值,将所有的中央隔离带单元频谱幅值组成中央隔离带频谱;
(8)确定静目标车辆和中央隔离带位置
8a)将静目标频谱中最大幅值对应的横坐标值,作为自车与静目标车辆之间的距离;
8b)将每一个中央隔离带单元频谱的横坐标值,作为自车与每一个中央隔离带单元之间的距离;
8c)将与静目标频谱最大幅值对应的步骤2c)中同一时刻接收到的回波信号的来波方向作为静目标车辆的方位角;
8d)将与频谱幅值对应的步骤2c)中同一时刻接收到的的回波信号的来波方向作为每一个中央隔离带单元的方位角;
8e)以自车为坐标原点,根据三角变换公式,计算得到静目标车辆和每一个中央隔离带单元的位置;
(9)比较目标车辆与中央隔离带之间的位置关系,若目标车辆的位置在中央隔离带位置的左方,则认为此目标车辆处于反向车道,视为虚警;否则,认为此目标车辆处于同向车道;
(10)防撞预警分析处理
10a)由防撞雷达预警***对处于同向车道的目标车辆进行防撞预警分析,根据安全距离公式,计算得到自车与每个目标车辆之间的安全距离;
10b)当自车与目标车辆之间的距离小于对应的安全距离时,自车的防撞雷达预警***对驾驶员发出声光报警或对车辆采取紧急制动。
本发明与现有技术相比具有以下优点:
第一,本发明利用二维信号处理方法和恒虚警检测方法,将动目标车辆、静目标车辆与中央隔离带分别检测出来,克服了现有技术不能对连续静止目标即中央隔离带检测识别的缺点,提高了本发明对道路目标的检测能力。
第二,本发明能够精确测量目标车辆的距离、速度和方位角,同时能够对中央隔离带进行定位,从而鉴别往返车道车辆,克服了现有技术存在较多虚警信息的缺点,提高了本发明中防撞预警***的预警效率。
第三,本发明利用预警***安全距离计算公式对道路中探测雷达范围以内的所有目标车辆进行预警分析,克服了现有技术仅对与自车距离最近的目标车辆进行预警分析的缺点,使得本发明提高了预警***的安全系数。
附图说明
图1是本发明的流程图;
图2是目标车辆与中央隔离带二维频谱图;
图3是目标车辆与中央隔离带检测位置示意图。
具体实施方式
下面结合附图1对本发明的步骤做进一步的详细描述。
步骤1,发射信号。
由自车防撞雷达阵列发射天线,连续发射多周期对称三角波调制的线性调频连续波信号。
上扫频段发射信号为:
S t _ up ( t 1 ) = A t cos 2 π [ ( f 0 - Δf 2 ) t 1 + μt 1 2 2 ]
下扫频段发射信号为:
S t _ down ( t 2 ) = A t cos 2 π [ ( f 0 + 3 Δf 2 ) t 2 + μ 2 2 2 ]
其中, kT < t 1 < kT + T 2 , kT + T 2 < t 2 < ( k + 1 ) T , k=0,1,…,N-1,N表示信号积累周期数,At表示发射信号幅度,f0表示载频,三角波调制斜率为:
μ=2Δf/T
其中,T表示调制周期,Δf表示调制带宽。
步骤2,接收回波信号。
2a)在一个雷达扫描周期内,防撞雷达的接收天线,持续接收每一时刻连续波信号经由前方道路环境中车辆、中央隔离带以及其它障碍物反射回来的多周期回波信号。
上扫频段回波信号为:
S r _ up ( t 1 ) = A r cos 2 &pi; [ ( f 0 - &Delta;f 2 ) ( t 1 - &tau; ) + &mu; 2 ( t 1 - &tau; ) 2 ]
下扫频段回波信号为:
S r _ down ( t 2 ) = A r cos 2 &pi; [ ( f 0 + 3 &Delta;f 2 ) ( t 2 - &tau; ) + &mu; 2 ( t 2 - &tau; ) 2 ]
其中,Ar=KAt,表示回波信号幅度,K表示衰减因子,τ=2r/c,表示回波延时,r表示目标与自车之间的距离,c表示光速,其它因子同上。
2b)对雷达扫描周期内每一时刻由相邻阵列天线单元接收的不同相位回波信号进行比相,测出它们之间的相位差。
2c)根据下面的来波方向与相位差之间的换算公式得到每一时刻回波信号的来波方向。
Figure BSA00000751551500067
其中,θ(t)表示t时刻回波信号的来波方向,λ表示连续波信号的波长,d表示相邻阵列天线单元的间距,
Figure BSA00000751551500071
表示t时刻相邻阵列天线单元所接收的不同回波信号之间的相位差。
步骤3,产生中频信号。
3a)将步骤1发射的连续波信号输入中频信号处理***中,作为混频器的本振信号。
3b)将步骤2a)接收的回波信号输入混频器中与本振信号相乘,得到一个混频信号。
3c)将混频信号通过低通滤波器产生一个中频信号。
上扫频段中频信号为:
S b _ up ( t 1 ) = 1 2 A t A r cos 2 &pi; [ f b t 1 + ( f 0 - &Delta;f 2 ) &tau; - &mu;&tau; 2 2 ]
下扫频段中频信号为:
S b _ down ( t 2 ) = 1 2 A t A r cos 2 &pi; [ f b t 1 + ( f 0 + 3 &Delta;f 2 ) &tau; - &mu;&tau; 2 2 ]
其中,fb=μτ,表示多普勒频移,其它因子同上。
步骤4,二维信号处理。
4a)将这个中频信号作为中频信号处理单元。
4b)根据下面一维快速傅里叶变换公式,对中频信号处理单元中的中频信号作一维快速傅里叶变换,得到包含距离信息的一维频谱,此一维频谱包含上扫频段一维频谱和下扫频段一维频谱。
上扫频段一维频谱:
X up ( &omega; ) = &Integral; kT + &tau; kT + T / 2 S b _ up ( t 1 ) exp ( - j&omega; t 1 ) dt 1
下扫频段一维频谱:
X down ( &omega; ) = &Integral; kT + T / 2 ( k + 1 ) T S b _ down ( t 2 ) exp ( - j&omega; t 2 ) dt 2
其中,k表示周期序号,T表示三角波调制周期,Sb_up(t1)表示上扫频段中频信号,Sb_down(t2)表示下扫频段中频信号。
4c)对一维频谱按照下式作离散快速傅里叶变换,得到包含距离和相对速度信息的二维频谱,此二维频谱包含上扫频段二维频谱和下扫频段二维频谱。
上扫频段二维频谱为:
V up ( &omega; ) = &Sigma; k = 0 N - 1 X up ( &omega; ) exp ( - jk&theta; )
下扫频段二维频谱为:
V down ( &omega; ) = &Sigma; k = 0 N - 1 X down ( &omega; ) exp ( - jk&theta; )
其中,N表示中频信号的周期数,k表示周期序号,Xup(ω)表示上扫频段二维频谱,Xdown(ω)表示下扫频段二维频谱。
4d)根据下面的多普勒频率配对公式,对上扫频段二维频谱和下扫频段二维频谱,作消除配对处理,得到距离-相对速度二维频谱。
r = cT 4 &Delta;F ( f - + f + ) &Delta;v = c 2 f 0 ( f - - f + )
其中,r表示道路前方目标与自车之间的距离,c表示光速,T表示三角波调制周期,ΔF表示调制带宽,f-表示下扫频段二维频谱频率,f+表示上扫频段二维频谱频率,Δv表示目标与自车之间的相对速度,f0表示中心频率。
4e)根据下面的速度与相对速度之间的变换公式,对距离-相对速度二维频谱进行变换,得到距离-速度二维频谱。
Vt=v0-ΔV/cosθ(t)
其中,Vt表示目标的速度,v0表示自车的速度,Δv表示目标与自车之间的相对速度,θ(t)表示t时刻回波信号的来波方向。
4f)对步骤3b)至4e)重复执行10次,将每次得到的距离-速度二维频谱幅值相加积累,得到高信噪比二维频谱,将其幅值以二维矩阵的形式记录,此二维矩阵由M1行N1列元素组成。
步骤5,二维恒虚警检测。
5a)将此M1行N1列元素组成的二维矩阵输入二维恒虚警检测***。
5b)选取二维矩阵中的第i行第j列矩阵元素,令i等于1,j等于1。
5c)以第i行第j列矩阵元素作为二维被检测单元,分别提取以二维被检测单元为中心的相邻行和列中的4个矩阵元素和6个矩阵元素,组成二维检测参考单元,根据下面的二维恒虚警检测参考单元的均值公式,求得二维检测参考单元的均值。
Z i , j = 1 MN - 1 ( &Sigma; m = i - M / 2 i + M / 2 &Sigma; n = j - N / 2 j + N / 2 x m , n - x i , j )
其中,Zi,j表示二维恒虚警检测参考单元的均值,M和N分别表示提取的以二维被检测单元为中心相邻的行和列矩阵元素的个数,xm,n表示二维检测参考单元幅值,xi,j表示二维被检测单元的幅值。
5d)将二维检测参考单元的均值乘以阈值因子得到二维恒虚警检测门限。
U = P 1 M + N + K
其中,U表示阈值因子,P表示虚警概率,一般取10-6,M和N分别表示提取的以二维被检测单元为中心相邻的行和列矩阵元素的个数,K表示线性因子,一般取0.5。
5e)将二维被检测单元的幅值与二维恒虚警检测门限作比较,当幅值大于检测门限时,记录此幅值,否则记为0。
5f)若i的值不等于M1,则将i的值加1后,转到步骤5c);否则,转到步骤5g)。
5g)若j的值不等于N1,则将j的值加1后,转到步骤5c);否则,转到步骤5h)。
5h)将步骤5e)记录的全部幅值组成静止频谱和动目标频谱,将静止频谱幅值以一维矩阵的形式记录,此一维矩阵由L1个元素组成。
步骤6,确定动目标车辆位置和速度。
6a)将动目标频谱最大幅值对应的横坐标值,作为自车与动目标车辆之间的距离Ri,对应的纵坐标值作为动目标车辆的速度vi
6b)将与动目标频谱最大幅值对应的步骤2c)中同一时刻接收到的回波信号的来波方向作为第i个动目标车辆的方位角θi
6c)以自车为坐标原点,根据下面的三角变换公式,计算得到动目标车辆的位置。
Rxi=Ri×sinθi   Ryi=Ri×cosθi
其中,Rxi表示第i个动目标车辆的横向位置,Ri表示自车与第i个动目标车辆之间的距离,Ryi表示第i个动目标车辆的纵向位置,θi表示第i个动目标车辆的方位角。
步骤7,一维恒虚警检测。
7a)将步骤5h)中L1个元素组成的一维矩阵输入一维恒虚警检测***。
7b)选取一维矩阵中的第k个矩阵元素,令k等于1。
7c)以第k个矩阵元素作为一维被检测单元,分别提取以一维被检测单元为中心的相邻的4个矩阵元素,组成一维检测参考单元,根据下面的一维恒虚警检测参考单元的均值公式,求得一维检测参考单元的均值。
Z k = 1 L - 1 ( &Sigma; l = k - L / 2 k + L / 2 y l - y k )
其中,Zk表示一维恒虚警检测参考单元的均值,L表示提取的与一维被检测单元为中心相邻元素的个数。,y1表示一维检测参考单元幅值,yk表示一维被检测单元的幅值。
7d)将一维检测参考单元的均值乘以阈值因子得到一维恒虚警检测门限。
U = P 1 L + K
其中,U表示阈值因子,P表示虚警概率,一般取10-6,L表示提取的与一维被检测单元为中心相邻元素的个数,K表示线性因子,一般取0.5。
7e)将一维被检测单元的幅值与一维恒虚警检测门限作比较,当幅值大于检测门限时,记录此幅值,否则记为0。
7f)若k的值不等于L1,则将k的值加1后,转到步骤7c);否则,转到步骤7g);
7g)将步骤7e)记录的全部幅值组成静目标频谱。
7h)将步骤7e)记为0的每一个被检测单元的幅值作为一个中央隔离带单元频谱幅值,将所有的中央隔离带单元频谱幅值组成中央隔离带频谱。
步骤8,确定静目标车辆和中央隔离带位置。
8a)将静目标频谱中最大幅值对应的横坐标值,作为自车与静目标车辆之间的距离。
8b)将每一个中央隔离带单元频谱的横坐标值,作为自车与每一个中央隔离带单元之间的距离。
8c)将与静目标频谱最大幅值对应的步骤2c)中同一时刻接收到的回波信号的来波方向作为静目标车辆的方位角θi′。
8d)将与频谱幅值对应的步骤2c)中同一时刻接收到的的回波信号的来波方向作为第i个中央隔离带单元的方位角θj
8e)以自车为坐标原点,根据下面的三角变换公式,计算得到静目标车辆和每一个中央隔离带单元的位置。
静目标车辆三角变换公式:
Rxi′=Ri′×sinθi′   Ryi′=Ri′×cosθi
中央隔离带单元三角变换公式:
Rxj=Rj×sinθj   Ryj=Rj×cosθj
其中,Rxi′表示第i个静目标车辆的横向位置,Ri′表示自车与第i个静目标车辆的距离,θi′表示第i个静目标车辆的方位角,Ryi′表示第i个静目标车辆的纵向位置,Rxj表示第j个中央隔离带单元的横向位置,Rj表示自车与第j个中央隔离带单元之间的距离,Ryj表示第i个中央隔离带单元的纵向位置,θj表示第j个中央隔离带单元的方位角。
步骤9,比较目标车辆与中央隔离带之间的位置关系,若目标车辆的位置在中央隔离带位置的左方,则认为此目标车辆处于反向车道,视为虚警;否则,认为此目标车辆处于同向车道。
步骤10,防撞预警分析处理。
10a)由防撞雷达预警***对处于同向车道的目标车辆进行防撞预警分析,根据下面的安全距离公式,计算得到自车与每个目标车辆之间的安全距离。
S f = v t &times; t 0 + ( &Delta;v ) 2 2 &times; a + S 0
其中,Sf表示自车与目标车辆之间的安全距离,vt表示目标车辆速度,t0表示驾驶员的反应时间,一般取1~2s,Δv表示自车与目标车辆之间的相对速度,a表示自车的最大刹车减速度,由自车刹车***决定,S0表示道路安全法中规定的两车静止时的最小间距,一般取2~3m。
10b)当自车与目标车辆之间的距离小于对应的安全距离时,自车的防撞雷达预警***对驾驶员发出声光报警或对车辆采取紧急制动。
下面结合附图2、附图3对本发明的仿真效果做进一步的描述。
1.仿真条件
本发明的仿真是在主频3.1GHZ的Intel(R) Core CPU(TM)i3-2100、内存1.91GB的硬件环境和MATLAB R2008a的软件环境下进行的。
2.仿真内容
本实验采用配置有毫米波汽车防撞雷达的车辆,实验环境设置为有中央隔离带的三车道高速公路环境下,以中央隔离带为基准,从内至外依次为内车道、中车道、外车道,自车的防撞雷达***对前方车辆目标和中央隔离带进行探测识别。已知自车行驶于内车道,速度为15m/s,高速公路中央隔离带位于自车左侧,与自车之间的距离为2m左右,目标车辆甲行驶于反向车道,速度为15m/s,目标车辆乙静止于中车道,目标车辆丙行驶于外车道,速度为12m/s。
本发明实施例中采用的调制信号的中心频率为77GHz,调频带宽为300MHz,三角波调制周期为2ms。
3.仿真效果分析
仿真实验中通过运用二维信号处理和恒虚警检测,能够分离出车辆目标和中央隔离带,并确定各自在道路中的分布。仿真结果分别如图2和图3所示,其中:
图2表示防撞雷达信号处理***通过对中频信号进行二维信号处理后的距离-速度频谱图,从图中可以看出,两个离散频谱峰值之间存在一段连续频谱,在这段连续频谱中存在一个频谱突变,根据坐标能够确定前方道路探测环境中存在两个动目标车辆、一个静目标车辆以及道路中央隔离带,通过对此距离-速度频谱作二维恒虚警检测,能够确定两个动目标车辆甲和丙的位置和速度,通过对二维恒虚警检测出的连续频谱作一维恒虚警检测,能够确定静目标车辆乙的位置和中央隔离带的位置,将动目标车辆、静目标车辆和中央隔离带的位置分布输出在显示屏上,如图3,图3中◇表示自车位置,□表示动目标车辆甲,○表示静目标车辆乙,Δ表示动目标车辆丙,·表示中央隔离带。从图3中可以看出此发明能够精确测量目标的距离、速度和方位角,同时能够对中央隔离带进行定位,有效的检测出高速公路中央隔离带并确定出其位置分布,克服了现有技术不能对连续静止目标即中央隔离带检测识别的缺点,提高了本发明对道路目标的检测能力。

Claims (11)

1.基于汽车防撞雷达的道路中央隔离带检测方法,包括如下步骤: 
(1)由自车防撞雷达阵列发射天线,连续发射多周期对称三角波调制的线性调频连续波信号; 
(2)接收回波信号 
2a)在一个雷达扫描周期内,防撞雷达的接收天线,持续接收每一时刻连续波信号经由前方道路环境中车辆、中央隔离带以及其它障碍物反射回来的多周期回波信号; 
2b)对雷达扫描周期内每一时刻由相邻阵列天线单元接收的不同相位回波信号进行比相,测出它们之间的相位差; 
2c)根据来波方向与相位差之间的换算公式得到每一时刻回波信号的来波方向; 
(3)产生中频信号 
3a)将步骤(1)发射的连续波信号输入中频信号处理***中,作为混频器的本振信号; 
3b)将步骤2a)接收的回波信号输入混频器中与本振信号相乘,得到一个混频信号; 
3c)将混频信号通过低通滤波器产生一个中频信号; 
其特征在于,还包括如下步骤: 
(4)二维信号处理 
4a)将中频信号作为中频信号处理单元; 
4b)根据一维快速傅里叶变换公式,对中频信号处理单元中的中频信号作一维快速傅里叶变换,得到包含距离信息的一维频谱,此一维频谱包含上扫频段一维频谱和下扫频段一维频谱; 
4c)根据离散快速傅里叶变换公式,对一维频谱作离散快速傅里叶变换,得到包含距离和相对速度信息的二维频谱,此二维频谱包含上扫频段二维频谱和下扫频段二维频谱; 
4d)根据多普勒频率配对公式,对上扫频段二维频谱和下扫频段二维频谱,作消除配对处理,得到距离-相对速度二维频谱; 
4e)根据速度与相对速度之间的变换公式,对距离-相对速度二维频谱进行变换,得到距离-速度二维频谱; 
4f)对步骤3b)至4e)重复执行10次,将每次得到的距离-速度二维频谱幅值相加积累,得到高信噪比二维频谱,将其幅值以二维矩阵的形式记录,此二维矩阵由M1行N1列元素组成; 
(5)二维恒虚警检测 
5a)将此M1行N1列元素组成的二维矩阵输入二维恒虚警检测***; 
5b)选取二维矩阵中的第i行第j列矩阵元素,其中,i等于3,j等于4; 
5c)以第i行第j列矩阵元素作为二维被检测单元,分别提取以二维被检测单元为中心的相邻行和列中的4个矩阵元素和6个矩阵元素,组成二维检测参考单元,根据二维恒虚警检测参考单元的均值公式,求得二维检测参考单元的均值; 
5d)将二维检测参考单元的均值乘以阈值因子得到二维恒虚警检测门限; 
5e)将二维被检测单元的幅值与二维恒虚警检测门限作比较,当幅值大于检测门限时,记录此幅值,否则记为0; 
5f)若i的值不等于M1-2,则将i的值加1后,转到步骤5c);否则,转到步骤5g); 
5g)若j的值不等于N1-3,则将j的值加1后,转到步骤5c);否则,转到步骤5h); 
5h)将步骤5e)记录的全部幅值组成静止频谱和动目标频谱,将静止频谱幅值以一维矩阵的形式记录,此一维矩阵由L1个元素组成; 
(6)确定动目标车辆位置和速度 
6a)将动目标频谱最大幅值对应的横坐标值,作为自车与动目标车辆之间的距离,对应的纵坐标值作为动目标车辆的速度; 
6b)将与动目标频谱最大幅值对应的步骤2c)中同一时刻接收到的回波信号的来波方向作为动目标车辆的方位角; 
6c)以自车为坐标原点,根据三角变换公式,计算得到动目标车辆的位置; 
(7)一维恒虚警检测 
7a)将步骤5h)中L1个元素组成的一维矩阵输入一维恒虚警检测***; 
7b)选取一维矩阵中的第k个矩阵元素,其中,k等于3; 
7c)以第k个矩阵元素作为一维被检测单元,分别提取以一维被检测单元为中心的相邻的4个矩阵元素,组成一维检测参考单元,根据一维恒虚警检测参考单元的均值公式,求得一维检测参考单元的均值; 
7d)将一维检测参考单元的均值乘以阈值因子得到一维恒虚警检测门限; 
7e)将一维被检测单元的幅值与一维恒虚警检测门限作比较,当幅值大于检测门限时,记录此幅值,否则记为0; 
7f)若k的值不等于L1-2,则将k的值加1后,转到步骤7c);否则,转到步骤7g); 
7g)将步骤7e)记录的全部幅值组成静目标频谱; 
7h)将步骤7e)记为0的每一个被检测单元的幅值作为一个中央隔离带单元频谱幅值,将所有的中央隔离带单元频谱幅值组成中央隔离带频谱; 
(8)确定静目标车辆和中央隔离带位置 
8a)将静目标频谱中最大幅值对应的横坐标值,作为自车与静目标车辆之间的距离; 
8b)将每一个中央隔离带单元频谱的横坐标值,作为自车与每一个中央隔离带单元之间的距离; 
8c)将与静目标频谱最大幅值对应的步骤2c)中同一时刻接收到的回波信号的来波方向作为静目标车辆的方位角; 
8d)将与频谱幅值对应的步骤2c)中同一时刻接收到的的回波信号的来波方向作为每一个中央隔离带单元的方位角; 
8e)以自车为坐标原点,根据三角变换公式,计算得到静目标车辆和每一个中央隔离带单元的位置; 
(9)比较目标车辆与中央隔离带之间的位置关系,若目标车辆的位置在中央隔离带位置的左方,则认为此目标车辆处于反向车道,视为虚警;否则,认为此目标车辆处于同向车道; 
(10)防撞预警分析处理 
10a)由防撞雷达预警***对处于同向车道的目标车辆进行防撞预警分析,根据安全距离公式,计算得到自车与每个目标车辆之间的安全距离; 
10b)当自车与目标车辆之间的距离小于对应的安全距离时,自车的防撞雷达预 警***对驾驶员发出声光报警或对车辆采取紧急制动。 
2.根据权利要求1所述的基于汽车防撞雷达的道路中央隔离带检测方法,其特征在于:步骤2c)所述的来波方向与相位差之间的换算公式为: 
Figure FSB0000124060350000041
其中,θ表示回波信号的来波方向,λ表示连续波信号的波长,d表示相邻阵列天线单元的间距,
Figure FSB0000124060350000042
表示相邻阵列天线单元所接收的不同回波信号之间的相位差。 
3.根据权利要求1所述的基于汽车防撞雷达的道路中央隔离带检测方法,其特征在于:步骤4b)所述的一维快速傅里叶变换公式为: 
其中,X(ω)表示一维快速傅里叶变换结果,k表示周期序号,T表示三角波调制周期,τ表示回波信号最大时间延迟,S(t)表示中频信号。 
4.根据权利要求1所述的基于汽车防撞雷达的道路中央隔离带检测方法,其特征在于:步骤4c)所述的离散快速傅里叶变换公式为: 
Figure FSB0000124060350000044
其中,V(ω)表示离散快速傅里叶变换结果,N表示中频信号的周期数,k表示周期序号,X(ω)表示一维快速傅里叶变换结果。 
5.根据权利要求1所述的基于汽车防撞雷达的道路中央隔离带检测方法,其特征在于:步骤4d)所述的多普勒频率配对公式为: 
Figure FSB0000124060350000045
其中,r表示道路前方目标与自车之间的距离,c表示光速,T表示三角波调制周期,ΔF表示调制带宽,f-表示下扫频段二维频谱频率,f+表示上扫频段二维频谱频率,Δv表示目标与自车之间的相对速度,f0表示中心频率。 
6.根据权利要求1所述的基于汽车防撞雷达的道路中央隔离带检测方法,其特征在于:步骤4e)所述的速度与相对速度之间的变换公式为; 
v0=v0-Δv/cosθ 
其中,vt表示目标的速度,v0表示自车的速度,Δv表示目标与自车之间的相对速度,θ表示回波信号的来波方向。 
7.根据权利要求1所述的基于汽车防撞雷达的道路中央隔离带检测方法,其特征在于:步骤5c)所述的二维恒虚警检测参考单元的均值公式为: 
Figure FSB0000124060350000051
其中,Z表示二维检测参考单元的均值,M和N分别表示提取的以二维被检测单元为中心相邻的行和列矩阵元素的个数,S表示全部二维检测参考单元幅值的和,E表示二维被检测单元的幅值。 
8.根据权利要求1所述的基于汽车防撞雷达的道路中央隔离带检测方法,其特征在于:步骤7c)所述的一维恒虚警检测参考单元的均值公式为: 
Figure FSB0000124060350000052
其中,Z表示一维恒虚警检测参考单元的均值,L表示提取的与一维被检测单元为中心相邻元素的个数,S表示全部一维检测参考单元幅值的和,E表示一维被检测单元的幅值。 
9.根据权利要求1所述的基于汽车防撞雷达的道路中央隔离带检测方法,其特征在于:步骤5d)和7d)所述的阈值因子由阈值因子公式得到: 
U=PW+K 
其中,U表示阈值因子,P表示虚警概率,P取10-6,W表示矩阵因子,步骤5d)中
Figure FSB0000124060350000053
M和N分别表示提取的以二维被检测单元为中心相邻的行和列矩阵元素的个数,步骤7d)中
Figure FSB0000124060350000054
L表示提取的与一维被检测单元为中心相邻元素的个数,K表示线性因子,K取0.5。 
10.根据权利要求1所述的基于汽车防撞雷达的道路中央隔离带检测方法,其特征在于:步骤6c)和8e)所述的三角变换公式为: 
Rx=R×sinθ    Ry=R×cosθ 
其中,Rx表示目标车辆或每一个中央隔离带单元的横向位置,R表示自车与目标车辆或每一个中央隔离带单元之间的距离,Ry表示目标车辆或每一个中央隔离带 单元的纵向位置,θ表示目标车辆或每一个中央隔离带单元的方位角。 
11.根据权利要求1所述的基于汽车防撞雷达的道路中央隔离带检测方法,其特征在于:步骤10a)所述的安全距离公式为: 
Figure FSB0000124060350000061
其中,Sf表示自车与目标车辆之间的安全距离,vt表示目标车辆速度,t0表示驾驶员的反应时间,t0取1~2s,Δv表示自车与目标车辆之间的相对速度,a表示自车的最大刹车减速度,由自车刹车***决定,S0表示道路安全法中规定的两车静止时的最小间距,S0取2~3m。 
CN201210244436.6A 2012-07-04 2012-07-04 基于汽车防撞雷达的道路中央隔离带检测方法 Expired - Fee Related CN102798863B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210244436.6A CN102798863B (zh) 2012-07-04 2012-07-04 基于汽车防撞雷达的道路中央隔离带检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210244436.6A CN102798863B (zh) 2012-07-04 2012-07-04 基于汽车防撞雷达的道路中央隔离带检测方法

Publications (2)

Publication Number Publication Date
CN102798863A CN102798863A (zh) 2012-11-28
CN102798863B true CN102798863B (zh) 2014-06-18

Family

ID=47198019

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210244436.6A Expired - Fee Related CN102798863B (zh) 2012-07-04 2012-07-04 基于汽车防撞雷达的道路中央隔离带检测方法

Country Status (1)

Country Link
CN (1) CN102798863B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9869765B2 (en) 2014-03-14 2018-01-16 Wistron Neweb Corporation Alarm system and method for vehicle

Families Citing this family (32)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104597452A (zh) * 2013-11-01 2015-05-06 中国科学院上海技术物理研究所 对称三角线性调频连续波激光雷达探测目标的方法
CN104943605B (zh) * 2014-03-27 2017-10-13 启碁科技股份有限公司 用于交通工具的警示***及警示方法
EP3139194B1 (en) * 2015-09-02 2022-04-20 Veoneer Sweden AB A close range filtering fmcw radar
CN106054193A (zh) * 2016-05-24 2016-10-26 深圳市雷博泰克科技有限公司 一种车辆周围多目标检测方法、处理器及毫米波雷达***
CN107783125B (zh) * 2016-08-25 2021-04-20 大连楼兰科技股份有限公司 旋翼无人机防撞毫米波雷达***及信号处理方法
CN107783091B (zh) * 2016-08-25 2021-08-03 大连楼兰科技股份有限公司 自动驾驶汽车防撞毫米波雷达信号处理方法
CN107783088B (zh) * 2016-08-25 2021-02-26 大连楼兰科技股份有限公司 应用在恒虚警检测中的目标配对及数据处理方法
CN107783124B (zh) * 2016-08-25 2021-07-06 大连楼兰科技股份有限公司 基于组合波形的旋翼无人机复杂环境防碰撞雷达***及信号处理方法
CN106597411B (zh) * 2016-12-30 2019-06-04 无锡市雷华科技有限公司 雷达信号处理方法
CN109917390A (zh) * 2017-12-12 2019-06-21 比亚迪股份有限公司 基于雷达的车辆检测方法和***
CN108398684A (zh) * 2018-01-12 2018-08-14 米传科技(上海)有限公司 一种减少地面及隔离带对汽车毫米波雷达干扰的方法
CN109358324B (zh) * 2018-10-31 2020-07-14 安徽四创电子股份有限公司 基于地面监视雷达谱数据的目标检测方法
CN109559525B (zh) * 2018-11-26 2020-08-04 厦门精益远达智能科技有限公司 一种基于毫米波雷达的超速监控方法、装置和设备
CN109669182B (zh) * 2018-12-05 2022-05-17 南京邮电大学 无源双基地sar动/静目标联合稀疏成像方法
CN109856626A (zh) * 2018-12-10 2019-06-07 北京瑞达速通科技有限公司 基于微波雷达的高铁站台端部防穿越报警***和方法
CN109594496A (zh) * 2018-12-10 2019-04-09 北京瑞达速通科技有限公司 一种分布式微波雷达复杂车型防砸***和方法
CN111452723B (zh) * 2019-01-18 2023-06-13 长沙智能驾驶研究院有限公司 智能车辆及其预警控制方法、装置、***及存储介质
CN109814102B (zh) * 2019-01-31 2020-10-27 厦门精益远达智能科技有限公司 一种单车道的超高监控方法、装置、设备和存储介质
CN110060468B (zh) * 2019-01-31 2021-04-30 厦门精益远达智能科技有限公司 基于毫米波雷达的异常路况的监测方法、装置和设备
CN109814103B (zh) * 2019-01-31 2020-09-04 厦门精益远达智能科技有限公司 一种基于雷达的超高监控方法、装置、设备和存储介质
CN109872566B (zh) * 2019-01-31 2021-01-26 厦门精益远达智能科技有限公司 一种车辆信息提示方法、装置、设备和存储介质
CN109615882A (zh) * 2019-01-31 2019-04-12 厦门精益远达智能科技有限公司 一种单车道的车流量统计方法、装置、设备和存储介质
CN109814104B (zh) * 2019-01-31 2021-01-22 厦门精益远达智能科技有限公司 一种基于雷达的障碍物检测方法、装置、设备和存储介质
CN109658715B (zh) * 2019-01-31 2021-01-22 厦门精益远达智能科技有限公司 多车道的车流量统计方法、装置、设备和存储介质
CN109919999B (zh) * 2019-01-31 2021-06-11 深兰科技(上海)有限公司 一种目标位置检测的方法及装置
CN110045376B (zh) * 2019-04-28 2021-06-01 森思泰克河北科技有限公司 可行驶区域获取方法、计算机可读存储介质及终端设备
CN112526498B (zh) * 2019-08-29 2023-12-12 比亚迪股份有限公司 车辆以及车载雷达的目标检测方法和装置
CN110531336A (zh) * 2019-09-20 2019-12-03 山东大学 一种物体检测识别方法及***
EP4053584A4 (en) * 2019-10-29 2023-11-29 Kyocera Corporation ELECTRONIC DEVICE, ELECTRONIC DEVICE CONTROL METHOD AND PROGRAM
KR102653404B1 (ko) * 2020-05-08 2024-04-02 주식회사 에이치엘클레무브 차량의 제어 장치 및 제어 방법
CN111751798A (zh) * 2020-07-22 2020-10-09 上海英恒电子有限公司 一种雷达测角方法
CN112014836B (zh) * 2020-09-21 2022-03-04 四川长虹电器股份有限公司 一种基于毫米波雷达的短距人员目标跟踪方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6085151A (en) * 1998-01-20 2000-07-04 Automotive Systems Laboratory, Inc. Predictive collision sensing system
JP3489514B2 (ja) * 1999-12-09 2004-01-19 株式会社デンソー Fmcwレーダ装置
DE102008038731A1 (de) * 2008-08-12 2010-02-18 Continental Automotive Gmbh Verfahren zur Erkennung ausgedehnter statischer Objekte

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9869765B2 (en) 2014-03-14 2018-01-16 Wistron Neweb Corporation Alarm system and method for vehicle

Also Published As

Publication number Publication date
CN102798863A (zh) 2012-11-28

Similar Documents

Publication Publication Date Title
CN102798863B (zh) 基于汽车防撞雷达的道路中央隔离带检测方法
Saponara et al. Radar-on-chip/in-package in autonomous driving vehicles and intelligent transport systems: Opportunities and challenges
US6518916B1 (en) Object recognition apparatus
US10627481B2 (en) Multi-resolution doppler processing
CN113296098B (zh) 用于雷达目标检测的侧向块监测
CN105197108B (zh) 一种基于汽车辅助驾驶***的多目标测向***及其方法
CN115575922B (zh) 基于车载fmcw毫米波雷达的动目标检测方法和装置
EP4036600A1 (en) Radar system to detect angles
US20240027611A1 (en) Obstacle detection method and system of mmwave radar and vehicle
US20210323560A1 (en) Vehicle speed calculation method, system, device, and storage medium
CN104977566A (zh) 用于一调频连续波雷达感测***的信号处理方法及装置
Pandey et al. Automated parking test using isar images from automotive radar
CN206649156U (zh) 汽车防撞雷达***
Yeh et al. Simultaneous range and radial velocity estimation with a single narrowband LFM pulse
WO2021085083A1 (ja) 電子機器、電子機器の制御方法、及びプログラム
Kruse et al. Target classification based on near-distance radar sensors
Zhang et al. A Novel target classification and identification algorithm for 77G LFMCW automotive collision avoidance radar
US11300660B2 (en) Determining relative velocity in a vehicle radar system
Schiementz Postprocessing architecture for an automotive radar network
Nathaniel et al. Modeling, Simulation and Analysis of Automotive Radar Signal Using Wavelet Transform Technique
CN111175717A (zh) 可抑制雷达近距离谐波的辅助驾驶方法和场景应用
EP4239364A1 (en) Electronic device, electronic device control method, and program
EP4174523A1 (en) Electronic device, method for controlling electronic device, and program
EP4390456A1 (en) Beamvector filtering and disambiguation using velocity profile pattern for low level radar ego-motion
EP4177638A1 (en) Detection and localization of non-line-of-sight objects using multipath radar reflections and map data

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20140618

Termination date: 20200704

CF01 Termination of patent right due to non-payment of annual fee