CN104215952A - 基于微动特性的车载目标识别系统及其识别方法 - Google Patents
基于微动特性的车载目标识别系统及其识别方法 Download PDFInfo
- Publication number
- CN104215952A CN104215952A CN201410428475.0A CN201410428475A CN104215952A CN 104215952 A CN104215952 A CN 104215952A CN 201410428475 A CN201410428475 A CN 201410428475A CN 104215952 A CN104215952 A CN 104215952A
- Authority
- CN
- China
- Prior art keywords
- frequency
- signal
- target
- time
- input end
- 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
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/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
- G01S13/08—Systems for measuring distance only
- G01S13/32—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated
- G01S13/34—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated using transmission of continuous, frequency-modulated waves while heterodyning the received signal, or a signal derived therefrom, with a locally-generated signal related to the contemporaneously transmitted signal
- G01S13/345—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated using transmission of continuous, frequency-modulated waves while heterodyning the received signal, or a signal derived therefrom, with a locally-generated signal related to the contemporaneously transmitted signal using triangular modulation
-
- 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/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
- G01S13/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
- G01S13/583—Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of continuous unmodulated waves, amplitude-, frequency-, or phase-modulated waves and based upon the Doppler effect resulting from movement of targets
- G01S13/584—Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of continuous unmodulated waves, amplitude-, frequency-, or phase-modulated waves and based upon the Doppler effect resulting from movement of targets adapted for simultaneous range and velocity measurements
-
- 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
- G01S7/415—Identification of targets based on measurements of movement associated with the target
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明属于车辆目标识别技术领域,特别涉及基于微动特性的车载目标识别系统及其识别方法。该基于微动特性的车载目标识别系统包括:安装在车辆前方的信号发射模块、安装在车辆前方的信号接收模块、设置在车辆上的模数转换器、设置在车辆上的数字信号处理器、以及安装在车辆操控台上的显示器;所述模数转换器的输入端电连接所述信号接收模块的输出端,输出端电连接所述数字信号处理器的输入端;所述数字信号处理器的输出端电连接显示器。
Description
技术领域
本发明属于车辆目标识别技术领域,特别涉及基于微动特性的车载目标识别系统及其识别方法,可用于识别行车环境中的车辆和行人。
背景技术
交通安全主动探测系统可以有效避免驾驶员疏忽、异常天气以及夜间行驶等因素导致的交通安全问题。车载探测系统获取目标信息,利用探测结果可以提前发现危险目标并进行报警,同时使驾驶人员有充分的反应时间尽量避免事故发生。现有车载探测系统通常选用毫米波雷达。毫米波雷达首先发射调频连续波信号,然后通过接收的回波提取目标的距离和速度信息。驾驶人员需要对行车环境和目标信息有更好的判断,因此,需要有效识别不同的目标。
目前,目标识别分为基于结构特征的目标识别和基于高分辨距离像的目标识别。基于结构特征的目标识别方法主要是利用目标自身的结构特性,如长度、宽度和厚度等具有明显几何特性的特征。基于高分辨距离像的目标识别方法主要利用目标回波特性的距离像。上述方法主要适用于脉冲体制雷达,对车辆和行人识别率低。
发明内容
本发明的目的在于提出基于微动特性的车载目标识别系统及其识别方法。本发明适用于车载系统,可以向驾驶人员提供更详细的行车环境和目标信息,可以有效地识别行车环境中的车辆和行人。
本发明的思路是:微多普勒频率是目标或者目标上微运动部件运动形式在雷达目标回波中的体现。因此,对雷达目标回波中的微多普勒信息进行分析处理,能够得到目标或者目标上微运动部件运动形式的相关信息。由于不同类型的目标决定了其微运动形式的差异,对微多普勒信号的时频分析便有助于提取关于目标类别的信息,从而实现目标类别的判断。本发明通过微多普勒时频分析,获取特征频率有效地识别车辆和行人。
为实现上述技术目的,本发明采用如下技术方案予以实现。
技术方案一:
基于微动特性的车载目标识别系统,其特征在于,包括:安装在车辆前方的信号发射模块、安装在车辆前方的信号接收模块、设置在车辆上的模数转换器、设置在车辆上的数字信号处理器、以及安装在车辆操控台上的显示器;所述模数转换器的输入端电连接所述信号接收模块的输出端,输出端电连接所述数字信号处理器的输入端;所述数字信号处理器的输出端电连接显示器。
该技术方案的特点和进一步改进在于:
所述基于微动特性的车载目标识别系统,还包括用于产生周期性的三角波的三角波发生器,所述信号发射模块包括发射天线、载波信号源和混频器;所述混频器的两个输入端对应电连接三角波发生器和载波信号源;
所述信号接收模块包括接收天线、信号分配器、90度移相器、第一变频器、第二变频器;所述信号分配器的输入端电连接接收天线,所述信号分配器的两个输出端对应电连接第一变频器的第一输入端和第二变频器的第一输入端,所述混频器的输出端分别电连接发射天线、第一变频器的第二输入端、以及90度移相器的输入端,所述90度移相器的输出端电连接第二变频器的第二输入端,所述模数转换器的输入端分别电连接第一变频器的输出端和第二变频器的输出端。
在第一变频器和模数转换器之间依次串接有第一抗泄露滤波器和第一低通滤波器,在第二变频器和模数转换器之间依次串接有第二抗泄露滤波器和第二低通滤波器;所述第一变频器的输出端电连接第一抗泄露滤波器的输入端,所述第二变频器的输出端电连接第二抗泄露滤波器的输入端;所述第一抗泄露滤波器的输出端电连接第一低通滤波器的输入端,所述第二抗泄露滤波器的输出端电连接第二低通滤波器的输入端;所述模数转换器的输入端分别电连接第一低通滤波器的输出端和第二低通滤波器的输出端。
技术方案二:
基于微动特性的目标识别方法,基于上述基于微动特性的车载目标识别系统,其特征在于,包括以下步骤:
S1:信号发射模块向车辆前方发射射频信号,并将射频信号传输至信号接收模块;
S2:所述信号接收模块接收回波信号,然后,根据射频信号和回波信号,生成相互正交的I通道差频模拟信号和Q通道差频模拟信号;
S3:所述模数转换器分别对I通道差频模拟信号和Q通道差频模拟信号进行模数转换,将模数转换得出的I通道差频数字信号和Q通道差频数字信号发送至数字信号处理器;
S4:所述数字信号处理器将I通道差频数字信号和Q通道差频数字信号合并为回波基频信号;通过对所述回波基频信号进行恒虚警检测,来判断对应时刻是否出现目标;所述目标为车辆或行人;
S5:对应时刻有目标出现时,所述数字信号处理器得出对应目标的速度;然后根据对应目标的速度,得出对应目标的类别。
该技术方案的特点和进一步改进在于:
所述基于微动特性的车载目标识别系统还包括用于产生周期性的三角波的三角波发生器,所述信号发射模块包括发射天线、载波信号源和混频器;所述混频器的两个输入端对应电连接三角波发生器和载波信号源;所述信号接收模块包括接收天线、信号分配器、90度移相器、第一变频器、第二变频器;
所述信号分配器的输入端电连接接收天线,所述信号分配器的两个输出端对应电连接第一变频器的第一输入端和第二变频器的第一输入端,所述混频器的输出端分别电连接发射天线、第一变频器的第二输入端、以及90度移相器的输入端,所述90度移相器的输出端电连接第二变频器的第二输入端;所述模数转换器的输入端分别电连接第一变频器的输出端和第二变频器的输出端;
在步骤S1中,所述三角波发生器产生周期性的三角波,根据三角波产生周期性的线性调频连续波,所述线性调频连续波的周期为所述三角波的周期,所述三角波的幅度用于控制线性调频连续波的频率;然后,将线性调频连续波传输至所述混频器;所述载波信号源生成载波信号,并将所述载波信号传输至所述混频器;所述混频器对线性调频连续波和载波信号进行混频,得到射频信号,然后将射频信号分别发送至发射天线、第一变频器、以及90度移相器;所述发射天线向车辆前方发射所述射频信号。
在步骤S2中,所述接收天线接收回波信号,并将回波信号分别发送至第一变频器和第二变频器,所述90度移相器对所述射频信号进行90度移相处理,将移相处理后信号发送至第二变频器;所述第一变频器对所述回波信号和所述射频信号进行混频处理,得出I通道差频模拟信号;所述第二变频器对所述回波信号和所述移相处理后信号进行混频处理,得出Q通道差频模拟信号。
在步骤S3中,所述I通道差频模拟信号表示为I(t),t为离散时间变量,所述Q通道差频模拟信号表示为Q(t);
在步骤S4中,所述数字信号处理器将I通道差频数字信号I(t)和Q通道差频数字信号Q(t)合并为回波基频信号s(t);
在步骤S4中,得出所述目标的距离R和所述目标的速度v:
其中,c为光速,Δf为线性调频连续波的最大频偏,f0为线性调频连续波的中心频率,正向调频时差拍信号频率,为负向调频时差拍信号频率,和分别为:
其中,ft为射频信号的频率,fr为回波信号的频率;
在步骤S5中,所述目标的速度v,判断目标的类别,当目标的速度v超过设定速度阈值时,所述目标的类别为车辆;反之当目标的速度v小于或等于设定速度阈值时,所述目标的类别为行人。
在步骤S5中,对应时刻有目标出现时,所述数字信号处理器对对应回波基频信号进行时频分析,得出对应回波基频信号的时频谱S(t,f),S(t,f)为:
S(t,f)=∫s(τ)·w(τ-t)·exp(-j2πfτ)dτ
其中,f表示离散频率变量,τ时间积分变量,w(t)表示窗函数,
然后得出对应回波基频信号的时频图,所述对应回波基频信号的时频图为S(t,f)离散化后时频图;
然后,根据S(t,f)离散化后时频图,得出对应目标的时频特征,所述对应目标的时频特征为时频熵f1或对应目标的平均瞬时多普勒谱的熵f2,f1为:
其中,i、j分别表示对应回波基频信号的时频图中第i行第j列位置处的像素点,M和N分别为对应回波基频信号的时频图的行数和列数,S(ti,fj)表示对应回波基频信号的时频图中第i行第j列位置处的像素点的能量;
f2为:
其中,F(fk)为F(f)的第k个傅里叶变换值,K为g(t)的快速傅里叶变换点数,S(t,fk)为对应回波基频信号的时频图中第k列的所有能量值构成的列向量,|S(t,fk)|表示对S(t,fk)取模。
在步骤S5中,在累积获得l个目标的时频熵之后,根据每个目标的时频熵,得出对应的训练样本;当对应目标的时频特征为时频熵f1时,l=M×N,当对应目标的时频特征为对应目标的平均瞬时多普勒谱的熵f2时,l=K;第i个训练样本包括:第i个目标的时频特征xi,以及第i个目标的分类标记yi,i取1至l;当第i个目标的速度超过设定速度阈值时,yi=1,否则yi=-1;
然后求解以下优化问题:
其中,w和b为待求解的两个参数(标量),||w||表示w的绝对值;
如果(w·xi)+b)>0,则第i个目标为车辆目标,否则第i个目标为行人目标。
本发明的有益效果为:第一,本发明适用于车载探测系统,能够向驾驶人员提供更详细的行车环境和目标信息。第二,本发明基于目标的微多普勒特性,通过时频分析,可以有效地识别行车环境中的车辆和行人。
附图说明
图1为本发明的基于微动特性的车载目标识别系统的结构示意图;
图2为本发明的基于微动特性的车载目标识别方法的流程图;
图3为本发明中三角波发生器产生的三角波和线性调频连续波的关系示意图;
图4为本发明实施例中恒虚警检测的流程图;
图5为本发明进行目标类别判断的流程图;
图6为仿真实验得出的多散射点恒虚警检测的结果示意图。
具体实施方式
下面结合附图对本发明作进一步说明:
参照图1,为本发明的基于微动特性的车载目标识别系统的结构示意图。该基于微动特性的车载目标识别系统包括:安装在车辆前方的信号发射模块、安装在车辆前方的信号接收模块、设置在车辆上的模数转换器、设置在车辆上的数字信号处理器(DSP芯片)、以及安装在车辆操控台上的显示器。所述模数转换器的输入端电连接所述信号接收模块的输出端,输出端电连接所述数字信号处理器的输入端;所述数字信号处理器的输出端电连接显示器。
本发明实施例中,还设置有用于产生周期性的三角波的三角波发生器,三角波发生器根据产生的周期性的三角波的幅度,产生线性调频连续波,其中三角波的周期即为调频连续波的周期。上述信号发射模块包括发射天线、载波信号源和混频器;所述混频器的两个输入端对应电连接三角波发生器和载波信号源;混频器用于对载波信号源产生的载波信号和调频连续波进行混频处理,产生两路相同的信号。所述混频器的输出端分别电连接发射天线、第一变频器的第二输入端、以及90度移相器的输入端,对于混频器产生的两路相同的信号,一路信号作为发射信号向外发射(通过发射天线向外发射),另一路传输给信号接收部分(即90度移相器和第一变频器)。
上述信号接收模块包括接收天线、信号分配器、90度移相器、第一变频器、第二变频器;所述信号分配器的输入端电连接接收天线,所述信号分配器的两个输出端对应电连接第一变频器的第一输入端和第二变频器的第一输入端,信号分配器将接收到的回波信号分别发送至第一变频器和第二变频器。上述90度移相器的输出端电连接第二变频器的第二输入端,90度移相器将接收到的信号进行90度移相处理,然后将经90度移相处理的信号传输至第二变频器。每个变频器对接收的两组信号进行下变频处理,其中,第一变频器产生I通道差频模拟信号,第二变频器产生Q通道差频模拟信号。所述模数转换器的输入端分别电连接第一变频器的输出端和第二变频器的输出端,模数转换器用于对接收到的信号进行模数转换(分别对I通道差频模拟信号、生Q通道差频模拟信号进行采样),由于模数转换器的输出端电连接所述数字信号处理器的输入端,所以数字信号处理器接收到经模数转换后的信号。
本发明实施例中,数字信号处理器将接收的信号合并为回波基频信号,然后通过恒虚警检测获取目标的距离、速度信息;同时,通过对回波基频信号进行时频分析,并提取对应目标的时频特征,将对应目标的时频特征传输给支持向量机(SVM)分类器进行判决,得出目标识别结果。数字信号处理器将目标的距离、速度信息以及目标识别结果呈献给驾驶员,保证其有效识别车辆和行人,进行预判。
本发明实施例中,在第一变频器和模数转换器之间依次串接有第一抗泄露滤波器和第一低通滤波器,在第二变频器和模数转换器之间依次串接有第二抗泄露滤波器和第二低通滤波器;所述第一变频器的输出端电连接第一抗泄露滤波器的输入端,所述第二变频器的输出端电连接第二抗泄露滤波器的输入端;所述第一抗泄露滤波器的输出端电连接第一低通滤波器的输入端,所述第二抗泄露滤波器的输出端电连接第二低通滤波器的输入端;所述模数转换器的输入端分别电连接第一低通滤波器的输出端和第二低通滤波器的输出端。
本发明实施例还提出了一种基于微动特性的车载目标识别方法,参照图2,为本发明的基于微动特性的车载目标识别方法的流程图。该基于微动特性的车载目标识别方法,以上述基于微动特性的车载目标识别系统为基础,包括以下步骤:
S1:信号发射模块向车辆前方发射射频信号,并将射频信号传输至信号接收模块。
具体地说,三角波发生器产生周期性的三角波,根据三角波产生周期性的线性调频连续波,所述线性调频连续波的周期为所述三角波的周期,所述三角波的幅度用于控制线性调频连续波的频率。参照图3,为本发明中三角波发生器产生的三角波和线性调频连续波的关系示意图。图3中,横轴表示本发明中三角波的电压,纵轴表示线性调频连续波的频率,T+为波形正半周持续时间,T-为波形负半周持续时间。由图3看出,三角波形式表现为正半周时间点电压从0线性增加,负半周时间点电压从最高值线性降低。对于本发明的三角波发生器,主要通过压控信号源实现,压控信号源输出的信号频率随输入信号的电压而改变,并且频率调谐率是线性的。
然后,将线性调频连续波传输至所述混频器;所述载波信号源生成载波信号,并将所述载波信号传输至所述混频器;所述混频器对线性调频连续波和载波信号进行混频,得到射频信号,然后将射频信号分别发送至发射天线、第一变频器、以及90度移相器;所述发射天线向车辆前方发射所述射频信号。
本发明实施例中,射频信号由正斜率调频信号和负斜率调频信号构建而成,其中,正斜率调频信号s(t0)和负斜率调频信号s(t1)分别为:
s(t0)=cos(2π(f0t0+Kt0 2/2))
s(t1)=cos(2π(f0t1-Kt1 2/2))
其中,K为线性调频连续波的调频斜率,f0为中心频率,t0为正周期时间点,t1为负周期时间点。
S2:所述信号接收模块接收回波信号,然后,根据射频信号和回波信号,生成相互正交的I通道差频模拟信号和Q通道差频模拟信号。
具体地说,接收天线接收回波信号(前方目标或散射点的回波信号),并将回波信号分别发送至第一变频器和第二变频器,所述90度移相器对所述射频信号进行90度移相处理,将移相处理后信号发送至第二变频器;所述第一变频器对所述回波信号和所述射频信号进行混频处理,得出I通道差频模拟信号;所述第二变频器对所述回波信号和所述移相处理后信号进行混频处理,得出Q通道差频模拟信号。
本发明实施例中,接收天线接收的回波信号的多普勒频移fd为:
其中,v表示目标速度,c为光速。
接收天线接收的回波信号的延迟效应fdelay为:
其中,R表示目标与车辆的距离,T为调频连续波的周期,Δf为调频连续波的最大频偏。
S3:所述模数转换器分别对I通道差频模拟信号和Q通道差频模拟信号进行模数转换,将模数转换得出的I通道差频数字信号和Q通道差频数字信号发送至数字信号处理器。
具体地说,在步骤S3中,所述I通道差频模拟信号表示为I(t),t为离散时间变量,所述Q通道差频模拟信号表示为Q(t);当t为I通道差频模拟信号的正半周波形持续时间(Q通道差频模拟信号的正半周波形持续时间)时,I(t)和Q(t)的具体形式分别为:
I(t)=Acos(πKt2)
Q(t)=Asin(πKt2)
当t为I通道差频模拟信号的负半周波形持续时间(Q通道差频模拟信号的负半周波形持续时间)时,I(t)和Q(t)的具体形式分别为:
I(t)=Acos(-πKt2)
Q(t)=Asin(-πKt2)
S4:所述数字信号处理器将I通道差频数字信号和Q通道差频数字信号合并为回波基频信号;通过对所述回波基频信号进行恒虚警检测,来判断对应时刻是否出现目标;所述目标为车辆或行人;
具体地说,回波基频信号表示为s(t),当t为I通道差频模拟信号的正半周波形持续时间(Q通道差频模拟信号的正半周波形持续时间)时,s(t)的具体形式为:
当t为I通道差频模拟信号的负半周波形持续时间(Q通道差频模拟信号的负半周波形持续时间)时,s(t)的具体形式为:
其中,K为线性调频连续波的调频斜率,A为常数(与目标散射率有关的常数)。
在步骤4中,对所述回波基频信号进行恒虚警检测,其具体过程为:
参照图4,为本发明实施例中恒虚警检测的流程图。采用单元平均恒虚警检测器,将回波基频信号s(t)送到由2L+1个延迟距离单元构成的延迟线上,检测距离单元D位于中心位置,两侧各取L个距离单元为参考单元,分别对应不同离散时间的输入信号。将所有参考单元中的s值求和后再除以2L就可以得到被检测单元处杂波背景的均值估计
其中,检测门限K'为设定的门限乘子。当调整门限乘子K'的大小时,可以改变门限U0的大小,从而保证了虚警率大小的恒定。当输入信号(回波基频信号s(t))超过检测门限时,判定为检测距离单元D所对应时刻出现目标,传输到下一步;未达到回波基频信号s(t)门限时,则判定为杂波或噪声,不予考虑。
在步骤S4中,得出所述目标的距离R和所述目标的速度v:
其中,c为光速,Δf为线性调频连续波的最大频偏,f0为线性调频连续波的中心频率,正向调频时差拍信号频率,为负向调频时差拍信号频率,和分别为:
其中,ft为射频信号的频率,fr为回波信号的频率。
S5:对应时刻有目标出现时,所述数字信号处理器得出对应目标的速度;然后根据对应目标的速度,得出对应目标的类别。
参照图5,为本发明进行目标类别判断的流程图。在步骤S5中,所述目标的速度v,判断目标的类别,当目标的速度v超过设定速度阈值时,所述目标的类别为车辆;反之当目标的速度v小于或等于设定速度阈值时,所述目标的类别为行人。
对应时刻有目标出现时,所述数字信号处理器对对应回波基频信号进行时频分析,得出对应回波基频信号的时频谱S(t,f),S(t,f)为:
S(t,f)=∫s(τ)·w(τ-t)·exp(-j2πfτ)dτ
其中,f表示离散频率变量,τ时间积分变量,w(t)表示窗函数,
然后得出对应回波基频信号的时频图,所述对应回波基频信号的时频图为S(t,f)离散化后时频图;
然后,根据S(t,f)离散化后时频图,得出对应目标的时频特征,所述对应目标的时频特征为时频熵f1、对应目标的平均瞬时多普勒谱的熵f2或对应目标的平均瞬时多普勒谱最大值所占比重f3,f1为:
其中,i、j分别表示对应回波基频信号的时频图中第i行第j列位置处的像素点,M和N分别为对应回波基频信号的时频图的行数和列数,S(ti,fj)表示对应回波基频信号的时频图中第i行第j列位置处的像素点的能量;
f2为:
其中,F(fk)为F(f)的第k个傅里叶变换值,F(f)=FFT(g(t)),F(f)为g(t)的快速傅里叶变换值,K为g(t)的快速傅里叶变换点数,S(t,fk)为对应回波基频信号的时频图中第k列的所有能量值构成的列向量,|S(t,fk)|表示对S(t,fk)取模。F(f)为平均瞬时多普勒谱,对于任意的时刻t′,时频谱S(t′,f)反映了该时刻目标所有散射点的瞬时频率。g(t)是散射点瞬时频率的线性组合,因此可以通过对g(t)进行快速傅里叶变换得到散射点的瞬时频率即F(f)。行走的人体目标时频谱中包含的瞬时频率分量比车辆要多,因此车辆平均瞬时频谱的熵比行走的人体目标要小。
f3为:f3=max(ρ),其中,f3描述了平均瞬时多普勒谱的最大值所占总能量的比例,归一化后最大值越大,特征三就越大。
步骤S5中,在累积获得l个目标的时频熵之后,根据每个目标的时频熵,得出对应的训练样本;当对应目标的时频特征为时频熵f1时,l=M×N,当对应目标的时频特征为对应目标的平均瞬时多普勒谱的熵f2时,l=K;第i个训练样本包括:第i个目标的时频特征xi,以及第i个目标的分类标记yi,i取1至l;当第i个目标的速度超过设定速度阈值时,yi=1,否则yi=-1;
然后求解以下优化问题:
其中,w和b为待求解的两个参数(标量),||w||表示w的绝对值;在得出w和b之后,根据w和b,建立分类面(w·x)+b=0;然后根据分类面(w·x)+b=0,对目标类别进行识别;
如果(w·xi)+b)>0,则第i个目标为车辆目标,否则第i个目标为行人目标。
本发明的效果可以通过以下仿真实验进行说明:
1)仿真条件
信号发射模块的参数如下:载频f0=3GHz,Nyquist采样率30MHz,信号接收模块的观测窗口为[10,2000]m。目标参数为:观测区间内有3个散射点,分别位于[600,780,1200]m,3个散射点的后向散射系数分别为1、0.5和0.8,3个散射点的散射点速度分别为100m/s、60m/s和15m/s,在进行恒虚警检测时,参数L=24,理论恒虚警率设为为1×10-6。背景噪声为零均值的高斯白噪声,方差为σ2,定义信噪比
2)仿真内容
利用以上仿真参数进行恒虚警检测过程,参照图6,为仿真实验得出的多散射点恒虚警检测的结果示意图。图6中,横轴表示时间,单位为s,纵轴表示经恒虚警处理后的回波基频信号的幅度,单位为v。
3)仿真结果分析
由图6可以看出,通过恒虚警检测,能够在正确的位置上得到3个散射点,即3个检测目标。如果恒虚警概率选的过低,则图中检测门限会上移,导致散射系数较小的第二个目标被误认为杂波或噪声,淹没在噪声里,造成漏警;反之,检测门限会下移,导致过多虚警现象的产生。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (9)
1.基于微动特性的车载目标识别系统,其特征在于,包括:安装在车辆前方的信号发射模块、安装在车辆前方的信号接收模块、设置在车辆上的模数转换器、设置在车辆上的数字信号处理器、以及安装在车辆操控台上的显示器;所述模数转换器的输入端电连接所述信号接收模块的输出端,输出端电连接所述数字信号处理器的输入端;所述数字信号处理器的输出端电连接显示器。
2.如权利要求1所述的基于微动特性的车载目标识别系统,其特征在于,还包括用于产生周期性的三角波的三角波发生器,所述信号发射模块包括发射天线、载波信号源和混频器;所述混频器的两个输入端对应电连接三角波发生器和载波信号源;
所述信号接收模块包括接收天线、信号分配器、90度移相器、第一变频器、第二变频器;所述信号分配器的输入端电连接接收天线,所述信号分配器的两个输出端对应电连接第一变频器的第一输入端和第二变频器的第一输入端,所述混频器的输出端分别电连接发射天线、第一变频器的第二输入端、以及90度移相器的输入端,所述90度移相器的输出端电连接第二变频器的第二输入端,所述模数转换器的输入端分别电连接第一变频器的输出端和第二变频器的输出端。
3.如权利要求1所述的基于微动特性的车载目标识别系统,其特征在于,在第一变频器和模数转换器之间依次串接有第一抗泄露滤波器和第一低通滤波器,在第二变频器和模数转换器之间依次串接有第二抗泄露滤波器和第二低通滤波器;所述第一变频器的输出端电连接第一抗泄露滤波器的输入端,所述第二变频器的输出端电连接第二抗泄露滤波器的输入端;所述第一抗泄露滤波器的输出端电连接第一低通滤波器的输入端,所述第二抗泄露滤波器的输出端电连接第二低通滤波器的输入端;所述模数转换器的输入端分别电连接第一低通滤波器的输出端和第二低通滤波器的输出端。
4.基于微动特性的目标识别方法,基于权利要求1所述的基于微动特性的车载目标识别系统,其特征在于,包括以下步骤:
S1:信号发射模块向车辆前方发射射频信号,并将射频信号传输至信号接收模块;
S2:所述信号接收模块接收回波信号,然后,根据射频信号和回波信号,生成相互正交的I通道差频模拟信号和Q通道差频模拟信号;
S3:所述模数转换器分别对I通道差频模拟信号和Q通道差频模拟信号进行模数转换,将模数转换得出的I通道差频数字信号和Q通道差频数字信号发送至数字信号处理器;
S4:所述数字信号处理器将I通道差频数字信号和Q通道差频数字信号合并为回波基频信号;通过对所述回波基频信号进行恒虚警检测,来判断对应时刻是否出现目标;所述目标为车辆或行人;
S5:对应时刻有目标出现时,所述数字信号处理器得出对应目标的速度;然后根据对应目标的速度,得出对应目标的类别。
5.如权利要求4所述的基于微动特性的目标识别方法,其特征在于,所述基于微动特性的车载目标识别系统还包括用于产生周期性的三角波的三角波发生器,所述信号发射模块包括发射天线、载波信号源和混频器;所述混频器的两个输入端对应电连接三角波发生器和载波信号源;所述信号接收模块包括接收天线、信号分配器、90度移相器、第一变频器、第二变频器;
所述信号分配器的输入端电连接接收天线,所述信号分配器的两个输出端对应电连接第一变频器的第一输入端和第二变频器的第一输入端,所述混频器的输出端分别电连接发射天线、第一变频器的第二输入端、以及90度移相器的输入端,所述90度移相器的输出端电连接第二变频器的第二输入端;所述模数转换器的输入端分别电连接第一变频器的输出端和第二变频器的输出端;
在步骤S1中,所述三角波发生器产生周期性的三角波,根据三角波产生周期性的线性调频连续波,所述线性调频连续波的周期为所述三角波的周期,所述三角波的幅度用于控制线性调频连续波的频率;然后,将线性调频连续波传输至所述混频器;所述载波信号源生成载波信号,并将所述载波信号传输至所述混频器;所述混频器对线性调频连续波和载波信号进行混频,得到射频信号,然后将射频信号分别发送至发射天线、第一变频器、以及90度移相器;所述发射天线向车辆前方发射所述射频信号。
6.如权利要求5所述的基于微动特性的目标识别方法,其特征在于,在步骤S2中,所述接收天线接收回波信号,并将回波信号分别发送至第一变频器和第二变频器,所述90度移相器对所述射频信号进行90度移相处理,将移相处理后信号发送至第二变频器;所述第一变频器对所述回波信号和所述射频信号进行混频处理,得出I通道差频模拟信号;所述第二变频器对所述回波信号和所述移相处理后信号进行混频处理,得出Q通道差频模拟信号。
7.如权利要求4所述的基于微动特性的目标识别方法,其特征在于,在步骤S3中,所述I通道差频模拟信号表示为I(t),t为离散时间变量,所述Q通道差频模拟信号表示为Q(t);
在步骤S4中,所述数字信号处理器将I通道差频数字信号I(t)和Q通道差频数字信号Q(t)合并为回波基频信号s(t);
在步骤S4中,得出所述目标的距离R和所述目标的速度v:
其中,c为光速,Δf为线性调频连续波的最大频偏,f0为线性调频连续波的中心频率,正向调频时差拍信号频率,为负向调频时差拍信号频率,和分别为:
其中,ft为射频信号的频率,fr为回波信号的频率;
在步骤S5中,所述目标的速度v,判断目标的类别,当目标的速度v超过设定速度阈值时,所述目标的类别为车辆;反之当目标的速度v小于或等于设定速度阈值时,所述目标的类别为行人。
8.如权利要求7所述的基于微动特性的目标识别方法,其特征在于,在步骤S5中,对应时刻有目标出现时,所述数字信号处理器对对应回波基频信号进行时频分析,得出对应回波基频信号的时频谱S(t,f),S(t,f)为:
S(t,f)=∫s(τ)·w(τ-t)·exp(-j2πfτ)dτ
其中,f表示离散频率变量,τ时间积分变量,w(t)表示窗函数,
然后得出对应回波基频信号的时频图,所述对应回波基频信号的时频图为S(t,f)离散化后时频图;
然后,根据S(t,f)离散化后时频图,得出对应目标的时频特征,所述对应目标的时频特征为时频熵f1或对应目标的平均瞬时多普勒谱的熵f2,f1为:
其中,i、j分别表示对应回波基频信号的时频图中第i行第j列位置处的像素点,M和N分别为对应回波基频信号的时频图的行数和列数,S(ti,fj)表示对应回波基频信号的时频图中第i行第j列位置处的像素点的能量;
f2为:
其中,F(fk)为F(f)的第k个傅里叶变换值,K为g(t)的快速傅里叶变换点数,S(t,fk)为对应回波基频信号的时频图中第k列的所有能量值构成的列向量,|S(t,fk)|表示对S(t,fk)取模。
9.如权利要求8所述的基于微动特性的目标识别方法,其特征在于,在步骤S5中,在累积获得l个目标的时频熵之后,根据每个目标的时频熵,得出对应的训练样本;当对应目标的时频特征为时频熵f1时,l=M×N,当对应目标的时频特征为对应目标的平均瞬时多普勒谱的熵f2时,l=K;第i个训练样本包括:第i个目标的时频特征xi,以及第i个目标的分类标记yi,i取1至l;当第i个目标的速度超过设定速度阈值时,yi=1,否则yi=-1;
然后求解以下优化问题:
其中,w和b为待求解的两个参数(标量),||w||表示w的绝对值;
如果(w·xi)+b)>0,则第i个目标为车辆目标,否则第i个目标为行人目标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410428475.0A CN104215952B (zh) | 2014-08-27 | 2014-08-27 | 基于微动特性的车载目标识别系统及其识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410428475.0A CN104215952B (zh) | 2014-08-27 | 2014-08-27 | 基于微动特性的车载目标识别系统及其识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104215952A true CN104215952A (zh) | 2014-12-17 |
CN104215952B CN104215952B (zh) | 2017-06-06 |
Family
ID=52097651
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410428475.0A Active CN104215952B (zh) | 2014-08-27 | 2014-08-27 | 基于微动特性的车载目标识别系统及其识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104215952B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106646455A (zh) * | 2016-09-09 | 2017-05-10 | 电子科技大学 | 一种路面散射特性测量方法及雷达装置 |
CN107291451A (zh) * | 2017-05-25 | 2017-10-24 | 深圳市冠旭电子股份有限公司 | 语音唤醒方法及装置 |
CN111406224A (zh) * | 2019-01-30 | 2020-07-10 | 深圳市大疆创新科技有限公司 | 一种目标可信度确定方法、一种目标识别方法、系统、车辆及存储介质 |
CN113219462A (zh) * | 2021-04-29 | 2021-08-06 | 森思泰克河北科技有限公司 | 基于时频图的目标识别方法、装置和终端设备 |
CN113341405A (zh) * | 2020-02-18 | 2021-09-03 | 南京大学 | 一种消除近距离处车轮微多普勒影响的雷达跟踪方法 |
CN114355337A (zh) * | 2021-12-31 | 2022-04-15 | 北京太格时代自动化系统设备有限公司 | 基于24GHz毫米波雷达技术的铁路施工侵限预警系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060066473A1 (en) * | 2004-09-29 | 2006-03-30 | Tdk Corporation | Pulse wave radar device |
CN102707266A (zh) * | 2012-05-24 | 2012-10-03 | 北京理工大学 | 一种具有抗干扰和多目标识别功能的雷达及其检测方法 |
CN103257346A (zh) * | 2013-05-15 | 2013-08-21 | 桂林电子科技大学 | 一种汽车防撞雷达多目标探测方法与系统 |
CN103914703A (zh) * | 2014-05-12 | 2014-07-09 | 西安电子科技大学 | 一种行人与车辆微动目标的分类识别方法 |
-
2014
- 2014-08-27 CN CN201410428475.0A patent/CN104215952B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060066473A1 (en) * | 2004-09-29 | 2006-03-30 | Tdk Corporation | Pulse wave radar device |
CN102707266A (zh) * | 2012-05-24 | 2012-10-03 | 北京理工大学 | 一种具有抗干扰和多目标识别功能的雷达及其检测方法 |
CN103257346A (zh) * | 2013-05-15 | 2013-08-21 | 桂林电子科技大学 | 一种汽车防撞雷达多目标探测方法与系统 |
CN103914703A (zh) * | 2014-05-12 | 2014-07-09 | 西安电子科技大学 | 一种行人与车辆微动目标的分类识别方法 |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106646455A (zh) * | 2016-09-09 | 2017-05-10 | 电子科技大学 | 一种路面散射特性测量方法及雷达装置 |
CN106646455B (zh) * | 2016-09-09 | 2019-02-12 | 电子科技大学 | 一种路面散射特性测量方法及雷达装置 |
CN107291451A (zh) * | 2017-05-25 | 2017-10-24 | 深圳市冠旭电子股份有限公司 | 语音唤醒方法及装置 |
CN111406224A (zh) * | 2019-01-30 | 2020-07-10 | 深圳市大疆创新科技有限公司 | 一种目标可信度确定方法、一种目标识别方法、系统、车辆及存储介质 |
CN113341405A (zh) * | 2020-02-18 | 2021-09-03 | 南京大学 | 一种消除近距离处车轮微多普勒影响的雷达跟踪方法 |
CN113341405B (zh) * | 2020-02-18 | 2023-10-13 | 南京大学 | 一种消除近距离处车轮微多普勒影响的雷达跟踪方法 |
CN113219462A (zh) * | 2021-04-29 | 2021-08-06 | 森思泰克河北科技有限公司 | 基于时频图的目标识别方法、装置和终端设备 |
CN114355337A (zh) * | 2021-12-31 | 2022-04-15 | 北京太格时代自动化系统设备有限公司 | 基于24GHz毫米波雷达技术的铁路施工侵限预警系统 |
Also Published As
Publication number | Publication date |
---|---|
CN104215952B (zh) | 2017-06-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108008391B (zh) | 一种基于fmcw的车载雷达多目标识别方法 | |
CN104215952A (zh) | 基于微动特性的车载目标识别系统及其识别方法 | |
CN100507601C (zh) | 双基地合成孔径雷达的双门限恒虚警运动目标检测方法 | |
CN102707285A (zh) | 车载毫米波防撞雷达系统的频域恒虚警检测方法 | |
CN102707266B (zh) | 一种具有抗干扰和多目标识别功能的雷达及其检测方法 | |
CN105842685B (zh) | 一种多目标雷达探测方法 | |
CN102866391B (zh) | 基于短时傅里叶变换和分数阶傅里叶变换的多目标检测方法 | |
CN102162848B (zh) | 基于伪随机编码与线性调频连续波的汽车主动防撞雷达 | |
US10802130B2 (en) | Marine target detection in cluttered environments | |
CN104898103B (zh) | 基于多通道杂波图的低速目标检测方法 | |
CN109324317A (zh) | 毫米波雷达系统及其定位测速方法 | |
CN102798863A (zh) | 基于汽车防撞雷达的道路中央隔离带检测方法 | |
CN105044712A (zh) | 一种微波栅栏雷达装置与目标检测方法 | |
CN103235310A (zh) | 一种车载毫米波列车防撞雷达系统 | |
DE112015003685B4 (de) | Niederschlagsbestimmungsvorrichtung | |
CN104160296A (zh) | 目标检测方法 | |
CN110058239A (zh) | 一种车载毫米波雷达装置及目标探测方法 | |
EP2942641A1 (en) | Radar device | |
CN102288942A (zh) | 一种毫米波雷达信号波形的设计方法 | |
CN106054193A (zh) | 一种车辆周围多目标检测方法、处理器及毫米波雷达系统 | |
CN104730507A (zh) | 基于预调制aic雷达距离成像的车载路障告警方法 | |
CN102253370A (zh) | 一种77GHz毫米波汽车防碰撞雷达信号处理装置 | |
CN203012135U (zh) | 一种调频连续波雷达系统 | |
CN102540163A (zh) | 一种基于常规相参雷达的直升机分类识别方法 | |
Wan et al. | Moving target detection using the 2D-FFT algorithm for automotive FMCW radars |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20220506 Address after: Qiuzhen building, Southeast campus, Changshu Institute of technology, No.99 Hushan Road, Changshu Southeast Economic Development Zone, Suzhou City, Jiangsu Province Patentee after: Changshu southeast high tech Venture Service Co.,Ltd. Address before: 215500 No. 68 Southeast Avenue, Changshu High-tech Zone, Suzhou City, Jiangsu Province Patentee before: SUZHOU WENJIE SENSING TECHNOLOGY Co.,Ltd. |
|
TR01 | Transfer of patent right |