WO2022068097A1 - 一种改进调频连续波雷达目标检测的降噪方法 - Google Patents
一种改进调频连续波雷达目标检测的降噪方法 Download PDFInfo
- Publication number
- WO2022068097A1 WO2022068097A1 PCT/CN2020/139042 CN2020139042W WO2022068097A1 WO 2022068097 A1 WO2022068097 A1 WO 2022068097A1 CN 2020139042 W CN2020139042 W CN 2020139042W WO 2022068097 A1 WO2022068097 A1 WO 2022068097A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- signal
- filter
- doppler
- echo
- continuous wave
- Prior art date
Links
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
- 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
- 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/35—Details of non-pulse systems
- G01S7/352—Receivers
-
- 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/418—Theoretical aspects
-
- 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
Abstract
一种改进调频连续波雷达目标检测的降噪方法,应用于汽车电子产品中,该方法包括:获取并处理雷达的回波混频信号,生成原始数据矩阵(100);将原始数据矩阵通过二维自适应滤波器进行滤波(200);对滤波后的原始数据矩阵进行处理(300)。该方法通过二维自适应滤波器来抑制汽车雷达信号中有效目标周围的噪声,因而在距离多普勒图中提高了目标幅度与周围噪声幅度的对比度,通过增加这种对比度,目标信号功率与周围噪声功率的比率增加,从而使以幅度或者功率为主要特征的目标检测方法,如恒虚警检测法,可以更有效地检测到雷达回波弱目标,从而提高检测率。
Description
本申请涉及汽车电子技术领域,更具体地,涉及一种改进调频连续波雷达目标检测的降噪方法。
汽车上安装的雷达是高级驾驶辅助系统中的重要传感器之一。雷达用于检测本车附近或者远处的目标,这些目标可以是其他车辆,也可以是行人,或者周围静止目标等。雷达检测这些目标的能力直接与这些目标与其周围背景噪声的信噪比SNR直接相关。高SNR能够有效地降低误检测的概率,从而改善目标检测能力。
通常,各种目标检测算法在低信噪比的情况下性能都会下降。下面这些情况(但不仅限于)都会导致信噪比下降。比如,当目标在远距离时,传输功率由于距离的损耗导致回波信号能量低;或者目标在靠近雷达视场边缘,其与雷达的夹角为大角度,而天线增益在大角度时增益低于天线零度方位(正前方)从而导致回波能量下降;另外通道的噪声或者接收器的噪声本身比较大的情况下,抬高了底噪从而降低了信噪比。因而提高信噪比成为常用的改善目标检测性能的主要手段。改善接收器底噪虽然可以增加信噪比,但是通过提高发射功率的方法可能不适用于功率受限的系统,比如车载雷达。
发明内容
为克服现有技术中如何提高信噪比的问题,本申请提供一种改进调频连续波雷达目标检测的降噪方法。
一种改进调频连续波雷达目标检测的降噪方法,应用于设置有雷达的汽车电子产品中,所述方法包括:
获取并处理雷达的回波混频信号,生成原始数据矩阵;
将所述原始数据矩阵通过二维自适应滤波器进行滤波;
对滤波后的所述原始数据矩阵进行处理;
其中,所述二维自适应滤波器包括对径向距离维滤波的径向距离滤波器、对多普勒维滤波的多普勒滤波器。
可选地,获取所述回波混频信号,包括:
所述雷达发射连续的线性调频信号,所述线性调频信号的回波被雷达接收后与发射线性调频信号混频的信号,经过ADC采样后,生成回波混频信号。
可选地,处理所述回波混频信号,生成原始数据矩阵,包括:
根据所述回波混频信号,建立原始数据矩阵S
n(k,l),其中,n为连续波雷达的第n个接收天线,l为第l个线性调频信号,k为在第l个线性调频信号的回波信号上的第k个采样点。
可选地,所述将所述原始数据矩阵通过二维自适应滤波器进行滤波,包括:
通过所述径向距离滤波器对距离维的径向距离的信号进行滤波;
通过所述多普勒滤波器对多普勒维的多普勒信号进行滤波。
可选地,所述通过所述径向距离滤波器对距离维的径向距离的信号进行滤波,包括:
将径向距离滤波器的第一输入参数a(k)初始化,a(k)是长度为M+1的滤波器系数,即a(k)=[a
0(k),a
1(k),…,a
M(k)];
所述第l个线性调频信号的回波混频信号的原始数据矩阵第k个采样点以及其前M+1个采样点构成中间输入样本
它们和第一输入参数a(k)输入至径向距离滤波器中,得到真实信号在第k个采样时刻的估计值y(k),y(k)与x(k,l)之差生成误差信号e
n;
通过第一输入参数a(k+1)对所述线性调频信号的回波混频信号进行迭代计算,输出初次滤波的线性调频信号矩阵。
可选地,所述通过所述径向距离滤波器对线性调频信号矩阵的同一所述线性调频信号进行滤波之后,还包括:
将所述滤波后的线性调频信号矩阵通过快速傅里叶变换获得径向距离信息。
可选地,所述通过所述多普勒滤波器对线性调频信号矩阵的相邻所述线性调频信号进行滤波,包括:
将多普勒滤波器的第二输入参数b(l)初始化,b(l)是长度为P+1的滤波器系数b(l)=[b
0(l),b
1(l),…,b
P(l)];
所述第l个线性调频信号的回波混频信号的原始数据矩阵第k个采样点以及其前P+1个线性调频信号的回波混频信号构成中间输入样本
它们和第二输入参数b(l)输入至多普勒滤波器中,得到在第k个采样时刻且第l个线性调频信号的回波混频信号真实信号的估计值y(l),y(l)与X(k,l)之差生成误差信号e
n(l);
通过输入第二参数b(l)对所述原始数据矩阵进行迭代计算,输出再次滤波的线性调频信号矩阵。
可选地,所述通过所述多普勒滤波器对多普勒维的多普勒信号进行滤波之后,还包括:
将所述线性调频信号矩阵通过快速傅里叶变换计算得到多普勒信息。
可选地,所述自适应算法可为归一化最小均方算法或时变最小均方算法或最小二乘法。
可选地,所述对滤波后的所述原始数据矩阵进行处理,包括:
将所述各通道的线性调频信号矩阵进行非相干叠加处理;
对所述非相干叠加处理后的结果进行恒虚预警检测;
将恒虚预警检测后的结果进行目标检测。
与现有技术相比,本申请的有益效果是:本申请通过二维自适应滤波器来减少汽车雷达信号中有效目标周围的噪声;特别是在距离多普勒图中(range Doppler map)大大提高目标的幅度与周围噪声幅值的对比度。通过增加这种对比度,目标信号功率与周围噪声的比率增加,从而让基于幅值的目标检测算法,如恒虚警检测法(CFAR)可以更有效地检测到雷达回波较弱的目标,从而提高了雷达的目标检测能力。有了较多的散射点被检测到,通过聚类的方法,雷达便可以进一步感知目标的尺寸,进而提供更多的信息,以帮助目标分类。同时,本申请是在软 件中完成,因而无需进行昂贵的硬件修改,从而使其更易于在现有算法之上实施。
图1为本申请实施例的方法的示意图。
图2为本申请实施例的线性调频信号的回波混频信号的原始数据矩阵的示意图。
图3为本申请实施例的带有二维自适应滤波器的一般雷达信号处理框图。
图4为本申请实施例的用于过滤具有多个通道的雷达数据的二维自适应滤波器框图。
图5为本申请实施例的范围自适应滤波器结构示意图。
图6为本申请实施例的复值多普勒滤波器结构示意图。
图7为本申请实施例的噪声抑制的模型
图8为本申请实施例的滤波过程的图示。a)带有较强噪声的RD。b)径向距离滤波器幅度响应。c)径向距离滤波之后的RD图。d)多普勒滤波器幅度响应。e)经过径向距离滤波和多普勒滤波器之后的RD图。
下面结合具体实施方式对本申请作进一步的说明。
本申请实施例的附图中相同或相似的标号对应相同或相似的部件;在本申请的描述中,需要理解的是,若有术语“上”、“下”、“左”、“右”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本申请和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此附图中描述位置关系的用语仅用于示例性说明,不能理解为对本专利的限制。
此外,若有“第一”、“第二”等术语仅用于描述目的,主要是用于区分不同的装置、元件或组成部分(具体的种类和构造可能相同也可能不同),并非用于表明或暗示所指示装置、元件或组成部分的相对重要性和数量,而不能理解为指示或者暗示相对重要性。
在如图1所述的实施例中,本申请提供了一种改进调频连续波雷达目标检测的降噪方法,应用于汽车电子产品中,本方法包括:
100,获取并处理雷达的回波混频信号,生成原始数据矩阵;在步骤100中,获取所述回波混频信号,包括:所述雷达发射连续的线性调频信号,所述线性调频信号的回波被雷达接收后 与发射线性调频信号混频的信号,经过ADC采样后,生成回波混频信号。处理所述回波混频信号,生成原始数据矩阵,包括:根据所述回波混频信号,建立原始数据矩阵S
n(k,l),其中,n为连续波雷达的第n个接收天线,l为第l个线性调频信号,k为在第l个线性调频信号的回波信号上的第k个采样点。
200,将所述原始数据矩阵通过二维自适应滤波器进行滤波;在步骤200中,所述二维自适应滤波器包括对径向距离维滤波的径向距离滤波器、对多普勒维滤波的多普勒滤波器。所述将所述原始数据矩阵通过二维自适应滤波器进行滤波,包括:通过所述径向距离滤波器对距离维的径向距离的信号进行滤波;通过所述多普勒滤波器对多普勒维的多普勒信号进行滤波。
300,对滤波后的所述原始数据矩阵进行处理;在步骤300中,对滤波后的所述原始数据矩阵进行处理,包括:将所述各通道的线性调频信号矩阵进行非相干叠加处理;对所述非相干叠加处理后的结果进行恒虚预警检测;将恒虚预警检测后的结果进行目标检测。
在本实施例中,本申请通过二维自适应滤波器来减少汽车雷达信号中有效目标周围的噪声;特别是在距离多普勒图中大大提高目标的幅度与周围噪声幅值的对比度。通过增加这种对比度,目标信号功率与周围噪声的比率增加,从而让基于幅值的目标检测算法,如恒虚警检测法可以更有效地检测到雷达回波较弱的目标,从而提高了雷达的目标检测能力。有了较多的散射点被检测到,通过聚类的方法,雷达便可以进一步感知目标的尺寸,进而提供更多的信息,以帮助目标分类。同时,本申请是在软件中完成,因而无需进行昂贵的硬件修改,从而使其更易于在现有算法之上实施。
在一些实施例中,获取所述回波混频信号,包括:所述雷达发射连续的线性调频信号,所述线性调频信号的回波被雷达接收后与发射线性调频信号混频的信号,经过ADC采样后,生成回波混频信号。在本实施例中,雷达为调频连续波雷达,调频连续波雷达是指发射频率受特定信号调制的连续波雷达。调频连续波雷达通过比较任意时刻回波信号频率与此时刻发射信号的频率的之差方法来得到目标的距离信息,距离正比于两者的频率差。目标的径向速度由获得的多普勒频率呈线性关系。与其他测距测速雷达相比,调频连续波雷达的结构更简单。本申请的回波混频信号为线性调频信号的回波被雷达接收后与发射线性调频信号混频的信号,经过模/数转换器ADC采样后,生成离散的回波混频信号。
在一些实施例中,处理所述回波混频信号,生成原始数据矩阵,包括:根据所述回波混频信号,建立原始数据矩阵S
n(k,l),其中,n为连续波雷达的第n个接收天线,l为第l个线性调频信号,k为在第l个线性调频信号chirp的回波信号上的第k个采样点。参见图2,根据线性调频信号建立原始数据矩阵。在本实施例中,每一个发射的chirp信号遇到目标后返回的回波存放在这个矩阵的每一列上,不同的chirp的回波就放在不同列上。对于每一个chirp的回波信号共有K个采样点,因此自适应滤波器对这K个采样点进行噪声抑制,此操作应用于l个chirp。由于这K个采样点里带有跟径向距离相关的信息,因此这第一步的操作是优化了在距离域上的信噪比。这里将在距离维上使用的自适应滤波称为距离滤波器。
在一些实施例中,所述将所述原始数据矩阵通过二维自适应滤波器进行滤波,包括:
通过所述径向距离滤波器对距离维的径向距离的信号进行滤波;
通过所述多普勒滤波器对多普勒维的多普勒信号进行滤波。
在本实施例中,二维自适应滤波器被应用于雷达的每个信道,雷达采集的线性调频信号被排放为如图2所示的矩阵。每一个发射的线性调频信号遇到目标后的回波经混频采样后存放在这个矩阵的每一列上;发射的不同的线性调频信号,其回波经混频采样后排放在该矩阵的不同列上。对每一列的信号进行傅里叶变换便得到跟径向距离相关的信息,因此当本申请里提到的自适应噪声滤波器用在这一列的信号上,就称为径向距离滤波器。通过连续的线性调频信号的各自接收信号的相互关系,便得到运动目标的相对运动信息。因此,本专利提出的自适应滤波器应用在线性调频信号之间时,就称为多普勒滤波器。参见图4,图4为本申请的二维自适应滤波器框图;二维自适应滤波器用于信道1、信道2、…、信道N;每个信道采用相同的二维自适应滤波器,其中,二维自适应滤波器包括径向距离滤波器Range filter和多普勒滤波器Doppler filter,且径向距离滤波器和多普勒滤波器输出端均通过快速傅里叶变换模块而得到滤波后的距离维信息以及多普勒维的信息。经过径向距离滤波后,沿距离维,即图2中的行,进行快速傅里叶变换FFT。FFT的结果是一个复数矩阵信号。因此,要将第二个滤波器应用于速度维,即多普勒维,滤波器必须是复数的。与距离滤波器不同的是,多普勒滤波器应用于多普勒维,即图2中的列,另外其运算需要复数运算。同样,该多普勒滤波器会进一步抑制在多普勒域上的噪声。本专利正是利用这两个滤波器构造了一个二维自适应滤波器,从而在两个维度 上有效地抑制噪声。
在上述实施例的一种实施方式中,参见图5,所述通过所述径向距离滤波器对距离维的径向距离的信号进行滤波,包括:
将距离径向滤波器的第一输入参数a(k)初始化,a(k)是长度为M+1的滤波器系数,即a(k)=[a
0(k),a
1(k),…,a
M(k)];
所述第l个线性调频信号的回波混频信号的原始数据矩阵第k个采样点以及其前M+1个采样点构成中间输入样本
它们和第一输入参数a(k)输入至径向距离滤波器中,得到真实信号在第k个采样时刻的估计值y(k),y(k)与x(k,l)之差生成误差信号e
n;
通过第一输入参数a(k+1)对所述线性调频信号的回波混频信号进行迭代计算,输出初次滤波的线性调频信号矩阵。
图5为径向距离滤波器的结构示意图,线性调频信号矩阵x(k,l)=s(k,l)+n(k,l),其中s(k,l)为真实信号,n(k,l)为噪声,y
n(k)是真实信号s(k,l)的估计值。第n个信道的滤波器模块(图5中虚线框内的内容)输出参数e
n和
作为自适应算法的输入。参数e
n是误差信号,它近似于通道中的噪声。参数
是中间输入样本。这两个参数是计算将在下一次迭代中使用的滤波器系数所必需的。第一输入参数a(k)是在自适应算法中得到更新的,即适应过程。滤波器系数α在每个线性调频脉冲的开始时被初始化,即使a
0(k)的值为1。随着来自线性调频脉冲的回波被滤波,径向距离滤波器的第一输入参数a(k)被更新。移位寄存器包含
值,这些值是算法在适应中使用的中间输入样本。通过第一输入参数a(k)对线性调频信号矩阵进行迭代计算,输出进行初次滤波的线性调频信号矩阵。所述通过所述径向距离滤波器对距离维的径向距离的信号进行滤波之后,还包括:将所述滤波后的线性调频信号矩阵通过快速傅里叶变换获得径向距离信息。在本实施例中,发射信号是线性调频信号,当遇到目标时其回波是一个延迟的线性调频信号,它与发射信号的混频信号是一个正弦信号,傅里叶变换便得到这个正弦波的频率,该频率与目标的距离有直接对应关系。因此通过傅里叶变换获得目标的径向距离信息。
在上述实施例的一种实施方式中,参见图6,所述通过所述多普勒滤波器对线性调频信号 矩阵的相邻所述线性调频信号进行滤波,包括:
将多普勒滤波器的第二输入参数b(l)初始化,b(l)是长度为P+1的滤波器系数b(l)=[b
0(l),b
1(l),…,b
P(l)];
所述第l个线性调频信号的回波混频信号的原始数据矩阵第k个采样点以及其前P+1个线性调频信号的回波混频信号构成中间输入样本
它们和第二输入参数b(l)输入至多普勒滤波器中,得到在第k个采样时刻且第l个线性调频信号的回波混频信号真实信号的估计值y(l),y(l)与X(k,l)之差生成误差信号e
n(l);
通过输入第二参数b(l)对所述原始数据矩阵进行迭代计算,输出再次滤波的线性调频信号矩阵。
在本实施例中,图6为多普勒滤波器的结构示意图,线性调频信号矩阵X(k,l)=S(k,l)+n(k,l),其中S(k,l)为真实信号经过径向距离的傅里叶变换,n(k,l)为噪声,y
n(l)是真实信号S(k,l)的估计值。第n个信道的滤波器模块(图6中虚线框内的内容)输出参数e
n和
作为自适应算法的输入。参数e
n是误差信号,它近似于通道中的噪声。参数
是中间输入样本。这两个参数是计算将在下一次迭代中使用的滤波器系数所必需的。第二输入参数b(l)是在自适应算法中得到更新的,即适应过程。滤波器系数b在每个线性调频脉冲的开始处被初始化,即使b
0(l)的值为1。随着来自线性调频脉冲的回波被滤波,多普勒滤波器的第二输入参数b(l)被更新。移位寄存器包含
值,这些值是算法在适应中使用的中间输入样本。通过第二输入参数b(l)对线性调频信号矩阵进行迭代计算,输出进行再次滤波的线性调频信号矩阵。所述通过所述多普勒滤波器对多普勒维的多普勒信号进行滤波之后,还包括:将所述线性调频信号矩阵通过快速傅里叶变换计算得到多普勒信息。
图7的框图公开了一个基于最小均方的自适应算法的示例。首先,信号y是真实信号s的估计值。目标是在最小均方意义上最小化误差e:
e=y-(s+n)
对e取平方:
e
2=(y-(s+n))
2
=y
2+s
2+n
2-2(sy+ny)+2sn
求e
2的平均值E(e
2):
E(e
2)=E(y
2)+E(s
2)+E(n
2)-2E(sy+ny)+2E(sn).
这里假设s和n不相关。当s是正弦信号,而n是宽带噪声时,它们之间时非相关的。因此上面的假设是成立的。E(e
2)的最小值为
E(e
2)
min=E(y
2)
min+E(s
2)+E(n
2)-2E(sy+ny)
max+0.
由于s和n不相关,所以E(sn)=0。当信号功率大于噪声功率时,且y很接近于真实信号s时,则E(e
2)就达到了最小值。
可见,该滤波器实际上是一个带通滤波器,它使得输出y接近真实信号s,因此从另一种意义上讲,噪声的功率被抑制了。如果s是单频正弦信号,那么这个滤波器时一个非常窄的带通滤波器,只允许s通过。而且E(e
2)_min等于实际通道的噪声功率。这意味着误差信号e对噪声做到了最理想的估计,因而能被有效的抑制。
在一些实施中,所述对滤波后的所述原始数据矩阵进行处理,包括:
将所述各通道的线性调频信号矩阵进行非相干叠加处理;在本实施例中,将所述N个通道的傅里叶变换后得到的距离-多普勒图(range-doppler map)进行合并获得一个合并后的距离-多普勒图。一种有效快捷的方法是对这N个通道的数据进行非相干叠加处理,即non-coherent sum,因此N个矩阵的数据就合并为一个矩阵。
对所述非相干叠加处理后的结果进行恒虚预警检测;在本实施例中,在雷达信号检测中,当外界干扰强度不断变化时,雷达能自动调整其灵敏度,使雷达的虚警概率保持不变,这种特性称为恒虚警率特性。也就是说,进行检测的门限值不是事先固定的,而是随外界强度变化而做相应的调整,因此它是通过找自适应门限值来进行目标检测的。
将恒虚预警检测后的结果进行目标检测。在本实施例中,在雷达信号检测中,当外界干扰强度不断变化时,雷达能自动调整其灵敏度,使雷达的虚警概率保持不变,这种特性称为恒虚警率特性。也就是说,进行检测的门限值不是事先固定的,而是随外界强度变化而做相应的调整,因此它是通过找自适应门限值来进行目标检测的。
参见图3,本申请公开了一种改进调频连续波雷达目标检测的降噪方法,通过调频连续波获取包含信号和噪音的线性调频信号矩阵,先通过径向距离滤波器对线性调频信号矩阵进行初步滤波,再通过多普勒滤波器对其进行再次滤波后。如果有多个通道的情况,可以将滤波后的多通道线性调频信号矩阵进行非相干叠加,再经过恒虚预警检测和目标检测,获取目标列表。参见图8,a)含有较高噪声的RD图。b)径向距离滤波器幅度响应。c)径向距离过滤器之后的RD图。d)多普勒滤波器幅度响应。e)径向距离和多普勒滤波器之后的RD图。通过图8可清晰看到目标的SNR显著增加,因此提高了这种弱目标的检测能力。本申请通过二维自适应滤波器来减少汽车雷达信号中有效目标周围的噪声;大大提高目标幅度与周围多普勒RD图中周围噪声的对比度。通过增加这种对比度,目标信号功率与周围噪声的比率增加,从而让CFAR可以更有效地检测目标。使雷达可以灵敏地观察到聚类,从而为目标的可能尺寸提供更多信息,以帮助目标分类。同时,本申请在软件中完成,而无需进行昂贵的硬件修改,从而使其更易于在现有算法之上实施。
显然,本申请的上述实施例仅仅是为清楚地说明本申请所作的举例,而并非是对本申请的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本申请的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本申请权利要求的保护范围之内。
Claims (11)
- 一种改进调频连续波雷达目标检测的降噪方法,其特征在于,应用于设置有雷达的汽车电子产品中,所述方法包括:获取并处理雷达的回波混频信号,生成原始数据矩阵;将所述原始数据矩阵通过二维自适应滤波器进行滤波;对滤波后的所述原始数据矩阵进行处理;其中,所述二维自适应滤波器包括对径向距离维滤波的径向距离滤波器、对多普勒维滤波的多普勒滤波器。
- 根据权利要求1所述的一种改进调频连续波雷达目标检测的降噪方法,其特征在于,获取所述回波混频信号,包括:所述雷达发射连续的线性调频信号,所述线性调频信号的回波被雷达接收后与发射线性调频信号混频的信号,经过ADC采样后,生成回波混频信号。
- 根据权利要求1所述的一种改进调频连续波雷达目标检测的降噪方法,其特征在于,处理所述回波混频信号,生成原始数据矩阵,包括:根据所述回波混频信号,建立原始数据矩阵S n(k,l),其中,n为连续波雷达的第n个接收天线,l为第l个线性调频信号,k为在第l个线性调频信号的回波信号上的第k个采样点。
- 根据权利要求1所述的一种改进调频连续波雷达目标检测的降噪方法,其特征在于,所述将所述原始数据矩阵通过二维自适应滤波器进行滤波,包括:通过所述径向距离滤波器对距离维的径向距离的信号进行滤波;通过所述多普勒滤波器对多普勒维的多普勒信号进行滤波。
- 根据权利要求4所述的一种改进调频连续波雷达目标检测的降噪方法,其特征在于,所述通过所述径向距离滤波器对距离维的径向距离的信号进行滤波,包括:将径向距离滤波器的第一输入参数a(k)初始化,a(k)是长度为M+1的滤波器系数,即a(k)=[a 0(k),a 1(k),…,a M(k)];所述第l个线性调频信号的回波混频信号的原始数据矩阵第k个采样点以及其前M+1个采样点构成中间输入样本 它们和第一输入参数a(k)输入至径向距离滤波器中,得到真实信号在第k个采样时刻的估计值y(k),y(k)与x(k,l) 之差生成误差信号e n;通过第一输入参数a(k+1)对所述线性调频信号的回波混频信号进行迭代计算,输出初次滤波的线性调频信号矩阵。
- 根据权利要求5所述的一种改进调频连续波雷达目标检测的降噪方法,其特征在于,所述通过所述径向距离滤波器对线性调频信号矩阵的同一所述线性调频信号进行滤波之后,还包括:将所述滤波后的线性调频信号矩阵通过快速傅里叶变换获得径向距离信息。
- 根据权利要求4所述的一种改进调频连续波雷达目标检测的降噪方法,其特征在于,所述通过所述多普勒滤波器对线性调频信号矩阵的相邻所述线性调频信号进行滤波,包括:将多普勒滤波器的第二输入参数b(l)初始化,b(l)是长度为P+1的滤波器系数b(l)=[b 0(l),b 1(l),…,b P(l)];所述第l个线性调频信号的回波混频信号的原始数据矩阵第k个采样点以及其前P+1个线性调频信号的回波混频信号构成中间输入样本 它们和第二输入参数b(l)输入至多普勒滤波器中,得到在第k个采样时刻且第l个线性调频信号的回波混频信号真实信号的估计值y(l),y(l)与X(k,l)之差生成误差信号e n(l);通过输入第二参数b(l)对所述原始数据矩阵进行迭代计算,输出再次滤波的线性调频信号矩阵。
- 根据权利要求7所述的一种改进调频连续波雷达目标检测的降噪方法,其特征在于,所述通过所述多普勒滤波器对多普勒维的多普勒信号进行滤波之后,还包括:将所述线性调频信号矩阵通过快速傅里叶变换计算得到多普勒信息。
- 根据权利要求5或7任一项所述的一种改进调频连续波雷达目标检测的降噪方法,其特征在于,所述自适应算法可以为归一化最小均方算法或时变最小均方算法或最小二乘法。
- 根据权利要求1所述的一种改进调频连续波雷达目标检测的降噪方法,其特征在于,所述对滤波后的所述原始数据矩阵进行处理,包括:将所述各通道的线性调频信号矩阵进行非相干叠加处理;对所述非相干叠加处理后的结果进行恒虚预警检测;将恒虚预警检测后的结果进行目标检测。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011051384.1A CN112630768B (zh) | 2020-09-29 | 2020-09-29 | 一种改进调频连续波雷达目标检测的降噪方法 |
CN202011051384.1 | 2020-09-29 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2022068097A1 true WO2022068097A1 (zh) | 2022-04-07 |
Family
ID=75302740
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/CN2020/139042 WO2022068097A1 (zh) | 2020-09-29 | 2020-12-24 | 一种改进调频连续波雷达目标检测的降噪方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN112630768B (zh) |
WO (1) | WO2022068097A1 (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115166664A (zh) * | 2022-06-12 | 2022-10-11 | 北京理工大学 | 一种基于二维匹配滤波结果的快速自适应旁瓣抑制方法 |
CN116256714A (zh) * | 2022-12-31 | 2023-06-13 | 南京航空航天大学 | 基于雷达成像的电磁散射特性测量方法及装置 |
CN116577738A (zh) * | 2023-07-12 | 2023-08-11 | 南京隼眼电子科技有限公司 | 雷达同频抗干扰方法、装置、雷达设备及存储介质 |
CN116643250A (zh) * | 2023-07-27 | 2023-08-25 | 西安电子科技大学 | 基于稀疏线性调频信号的宽带雷达高速目标相参积累方法 |
CN116955444A (zh) * | 2023-06-15 | 2023-10-27 | 共享易付(广州)网络科技有限公司 | 基于大数据分析的采集噪声点挖掘方法及系统 |
CN115166664B (zh) * | 2022-06-12 | 2024-05-10 | 北京理工大学 | 一种基于二维匹配滤波结果的快速自适应旁瓣抑制方法 |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113325410A (zh) * | 2021-05-28 | 2021-08-31 | 浙江大华技术股份有限公司 | 一种雷达天线信号处理方法、装置、控制设备及存储介质 |
CN113534125A (zh) * | 2021-06-04 | 2021-10-22 | 惠州市德赛西威汽车电子股份有限公司 | 一种估算目标模糊速度的方法 |
CN113627398B (zh) * | 2021-10-11 | 2022-01-04 | 西安瀚博电子科技有限公司 | 一种基于自适应重构滤波的信号特征检测方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110187586A1 (en) * | 2009-12-18 | 2011-08-04 | Thales | Radar detection method, notably for airborne radars implementing an obstacle detection and avoidance function |
CN103339525A (zh) * | 2010-12-21 | 2013-10-02 | 塔莱斯公司 | 用于监测地势变化的方法和设备 |
US20160061947A1 (en) * | 2014-08-27 | 2016-03-03 | Texas Instruments Incorporated | Fmcw doppler processing algorithm for achieving cw performance |
US20160245911A1 (en) * | 2013-07-15 | 2016-08-25 | Texas Instruments Incorporated | 2-D Object Detection in Radar Applications |
US20180172816A1 (en) * | 2016-12-20 | 2018-06-21 | National Chung-Shan Institute Of Science & Technology | Signal processing apparatus of a continuous-wave (cw) radar sensing system |
US20190317187A1 (en) * | 2018-04-11 | 2019-10-17 | Infineon Technologies Ag | Fmcw radar with interference signal rejection |
-
2020
- 2020-09-29 CN CN202011051384.1A patent/CN112630768B/zh active Active
- 2020-12-24 WO PCT/CN2020/139042 patent/WO2022068097A1/zh active Application Filing
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110187586A1 (en) * | 2009-12-18 | 2011-08-04 | Thales | Radar detection method, notably for airborne radars implementing an obstacle detection and avoidance function |
CN103339525A (zh) * | 2010-12-21 | 2013-10-02 | 塔莱斯公司 | 用于监测地势变化的方法和设备 |
US20160245911A1 (en) * | 2013-07-15 | 2016-08-25 | Texas Instruments Incorporated | 2-D Object Detection in Radar Applications |
US20160061947A1 (en) * | 2014-08-27 | 2016-03-03 | Texas Instruments Incorporated | Fmcw doppler processing algorithm for achieving cw performance |
US20180172816A1 (en) * | 2016-12-20 | 2018-06-21 | National Chung-Shan Institute Of Science & Technology | Signal processing apparatus of a continuous-wave (cw) radar sensing system |
US20190317187A1 (en) * | 2018-04-11 | 2019-10-17 | Infineon Technologies Ag | Fmcw radar with interference signal rejection |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115166664A (zh) * | 2022-06-12 | 2022-10-11 | 北京理工大学 | 一种基于二维匹配滤波结果的快速自适应旁瓣抑制方法 |
CN115166664B (zh) * | 2022-06-12 | 2024-05-10 | 北京理工大学 | 一种基于二维匹配滤波结果的快速自适应旁瓣抑制方法 |
CN116256714A (zh) * | 2022-12-31 | 2023-06-13 | 南京航空航天大学 | 基于雷达成像的电磁散射特性测量方法及装置 |
CN116256714B (zh) * | 2022-12-31 | 2023-10-13 | 南京航空航天大学 | 基于雷达成像的电磁散射特性测量方法及装置 |
CN116955444A (zh) * | 2023-06-15 | 2023-10-27 | 共享易付(广州)网络科技有限公司 | 基于大数据分析的采集噪声点挖掘方法及系统 |
CN116577738A (zh) * | 2023-07-12 | 2023-08-11 | 南京隼眼电子科技有限公司 | 雷达同频抗干扰方法、装置、雷达设备及存储介质 |
CN116577738B (zh) * | 2023-07-12 | 2023-09-12 | 南京隼眼电子科技有限公司 | 雷达同频抗干扰方法、装置、雷达设备及存储介质 |
CN116643250A (zh) * | 2023-07-27 | 2023-08-25 | 西安电子科技大学 | 基于稀疏线性调频信号的宽带雷达高速目标相参积累方法 |
CN116643250B (zh) * | 2023-07-27 | 2023-10-20 | 西安电子科技大学 | 基于稀疏线性调频信号的宽带雷达高速目标相参积累方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112630768B (zh) | 2024-04-02 |
CN112630768A (zh) | 2021-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2022068097A1 (zh) | 一种改进调频连续波雷达目标检测的降噪方法 | |
CN111352102B (zh) | 一种基于调频连续波雷达的多目标个数检测方法及装置 | |
KR102186191B1 (ko) | 간섭이 억제된 레이더 감지 | |
CN107683422B (zh) | 车辆雷达系统 | |
JP2019090824A (ja) | レーダーにおける角度分解能 | |
US7403153B2 (en) | System and method for reducing a radar interference signal | |
US10031221B2 (en) | System and method for estimating number and range of a plurality of moving targets | |
CN109444820B (zh) | 杂波和干扰共存时多通道雷达先干扰抑制后目标检测方法 | |
US9140783B2 (en) | Radar device | |
JP2020067455A (ja) | 妨害信号抑圧を行うfmcwレーダー | |
CN109425854B (zh) | 多分辨率多普勒处理 | |
US10656248B2 (en) | Radar post processing for sidelobe suppression | |
US4916452A (en) | Process for range measurement with a pulse radar of high pulse repetition frequency | |
US20050179579A1 (en) | Radar receiver motion compensation system and method | |
CN111896926B (zh) | 一种基于强杂波抑制的低空目标检测方法及系统 | |
US8760340B2 (en) | Processing radar return signals to detect targets | |
CN110907929B (zh) | 一种基于双门限检测的车载雷达目标检测方法及装置 | |
CN111220955B (zh) | 基于垂直阵列协方差矩阵特征分解的机载气象雷达地杂波抑制方法 | |
US20230168367A1 (en) | CFAR Phased Array Pre-Processing Using Noncoherent and Coherent Integration in Automotive Radar Systems | |
CN115877344B (zh) | 一种宽带检测跟踪识别一体化处理的雷达探测方法及系统 | |
CN111856406B (zh) | 一种基于fmcw雷达回波的目标检测方法及装置 | |
CN112014806A (zh) | 一种复杂干扰场景下的机载雷达无意干扰抑制方法 | |
CN109061626A (zh) | 一种步进频相参处理检测低信杂比动目标的方法 | |
CN107976663A (zh) | 一种基于子空间投影的外辐射源雷达对目标压缩检测方法 | |
CN109581366B (zh) | 一种基于目标导向矢量失配的离散旁瓣杂波识别方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 20956113 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 20956113 Country of ref document: EP Kind code of ref document: A1 |