CN112630768B - 一种改进调频连续波雷达目标检测的降噪方法 - Google Patents
一种改进调频连续波雷达目标检测的降噪方法 Download PDFInfo
- Publication number
- CN112630768B CN112630768B CN202011051384.1A CN202011051384A CN112630768B CN 112630768 B CN112630768 B CN 112630768B CN 202011051384 A CN202011051384 A CN 202011051384A CN 112630768 B CN112630768 B CN 112630768B
- Authority
- CN
- China
- Prior art keywords
- signal
- filter
- frequency modulation
- doppler
- radial distance
- 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 51
- 238000000034 method Methods 0.000 title claims abstract description 40
- 230000009467 reduction Effects 0.000 title claims abstract description 17
- 239000011159 matrix material Substances 0.000 claims abstract description 77
- 230000003044 adaptive effect Effects 0.000 claims abstract description 45
- 238000001914 filtration Methods 0.000 claims abstract description 44
- 238000005070 sampling Methods 0.000 claims description 30
- 238000002592 echocardiography Methods 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 4
- 230000001131 transforming effect Effects 0.000 claims description 2
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 12
- 230000006978 adaptation Effects 0.000 description 6
- 238000012986 modification Methods 0.000 description 5
- 230000004048 modification Effects 0.000 description 5
- 230000004044 response Effects 0.000 description 4
- 230000009466 transformation Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 230000001427 coherent effect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 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/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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Electromagnetism (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本申请涉及一种改进调频连续波雷达目标检测的降噪方法,应用于汽车电子产品中,所述方法包括:获取并处理雷达的回波混频信号,生成原始数据矩阵;将所述原始数据矩阵通过二维自适应滤波器进行滤波;对滤波后的所述原始数据矩阵进行处理。有益效果是:本申请通过二维自适应滤波器来抑制汽车雷达信号中有效目标周围的噪声;因而在距离多普勒图中大大提高目标幅度与周围噪声幅度的对比度。通过增加这种对比度,目标信号功率与周围噪声功率的比率增加,从而让以幅度或者功率为主要特征的目标检测方法,如恒虚警检测法可以更有效地检测到雷达回波弱目标,从而提高检测率。
Description
技术领域
本申请涉及汽车电子技术领域,更具体地,涉及一种改进调频连续波雷达目标检测的降噪方法。
背景技术
汽车上安装的雷达是高级驾驶辅助系统中的重要传感器之一。雷达用于检测本车附近或者远处的目标,这些目标可以是其他车辆,也可以是行人,或者周围静止目标等。雷达检测这些目标的能力直接与这些目标与其周围背景噪声的信噪比SNR直接相关。高SNR能够有效地降低误检测的概率,从而改善目标检测能力。
通常,各种目标检测算法在低信噪比的情况下性能都会下降。下面这些情况(但不仅限于)都会导致信噪比下降。比如,当目标在远距离时,传输功率由于距离的损耗导致回波信号能量低;或者目标在靠近雷达视场边缘,其与雷达的夹角为大角度,而天线增益在大角度时增益低于天线零度方位(正前方)从而导致回波能量下降;另外通道的噪声或者接收器的噪声本身比较大的情况下,抬高了底噪从而降低了信噪比。因而提高信噪比成为常用的改善目标检测性能的主要手段。改善接收器底噪虽然可以增加信噪比,但是通过提高发射功率的方法可能不适用于功率受限的系统,比如车载雷达。
发明内容
为克服现有技术中如何提高信噪比的问题,本申请提供一种改进调频连续波雷达目标检测的降噪方法。
一种改进调频连续波雷达目标检测的降噪方法,应用于设置有雷达的汽车电子产品中,所述方法包括:
获取并处理雷达的回波混频信号,生成原始数据矩阵;
将所述原始数据矩阵通过二维自适应滤波器进行滤波;
对滤波后的所述原始数据矩阵进行处理;
其中,所述二维自适应滤波器包括对径向距离维滤波的径向距离滤波器、对多普勒维滤波的多普勒滤波器。
可选地,获取所述回波混频信号,包括:
所述雷达发射连续的线性调频信号,所述线性调频信号的回波被雷达接收后与发射线性调频信号混频的信号,经过ADC采样后,生成回波混频信号。
可选地,处理所述回波混频信号,生成原始数据矩阵,包括:
根据所述回波混频信号,建立原始数据矩阵Sn(k,l),其中,n为连续波雷达的第n个接收天线,l为第l个线性调频信号,k为在第l个线性调频信号的回波信号上的第k个采样点。
可选地,所述将所述原始数据矩阵通过二维自适应滤波器进行滤波,包括:
通过所述径向距离滤波器对距离维的径向距离的信号进行滤波;
通过所述多普勒滤波器对多普勒维的多普勒信号进行滤波。
可选地,所述通过所述径向距离滤波器对距离维的径向距离的信号进行滤波,包括:
将径向距离滤波器的第一输入参数a(k)初始化,a(k)是长度为M+1的滤波器系数,即a(k)=[a0(k),a1(k),…,aM(k)];
所述第l个线性调频信号的回波混频信号的原始数据矩阵第k个采样点以及其前M+1个采样点构成中间输入样本它们和第一输入参数a(k)输入至径向距离滤波器中,得到真实信号在第k个采样时刻的估计值y(k),y(k)与x(k,l)之差生成误差信号en;
根据所述误差信号en和中间输入样本通过自适应算法更新在第k+1时刻时的第一输入参数a(k+1);
通过第一输入参数a(k+1)对所述线性调频信号的回波混频信号进行迭代计算,输出初次滤波的线性调频信号矩阵。
可选地,所述通过所述径向距离滤波器对线性调频信号矩阵的同一所述线性调频信号进行滤波之后,还包括:
将所述滤波后的线性调频信号矩阵通过快速傅里叶变换获得径向距离信息。
可选地,所述通过所述多普勒滤波器对线性调频信号矩阵的相邻所述线性调频信号进行滤波,包括:
将多普勒滤波器的第二输入参数b(l)初始化b(l)是长度为P+1的滤波器系数b(l)=[b0(l),b1(l),…,bP(l)];
所述第l个线性调频信号的回波混频信号的原始数据矩阵第k个采样点以及其前P+1个线性调频信号的回波混频信号构成中间输入样本它们和第二输入参数b(l)输入至多普勒滤波器中,得到在第k个采样时刻且第l个线性调频信号的回波混频信号真实信号的估计值y(l),y(l)与X(k,l)之差生成误差信号en(l);
根据所述误差信号en(l)和中间输入样本通过自适应算法更新第二输入参数b(l);
通过输入第二参数b(l)对所述原始数据矩阵进行迭代计算,输出再次滤波的线性调频信号矩阵。
可选地,所述通过所述多普勒滤波器对多普勒维的多普勒信号进行滤波之后,还包括:
将所述线性调频信号矩阵通过快速傅里叶变换计算得到多普勒信息。
可选地,所述自适应算法可为归一化最小均方算法或时变最小均方算法或最小二乘法。
可选地,所述根据所述误差信号en和中间输入样本通过自适应算法更新径向距离滤波器在时刻k+1的参数a(k+1),包括:
通过计算求得,其中,Δ是步进常数。
可选地,所述对滤波后的所述原始数据矩阵进行处理,包括:
将所述各通道的线性调频信号矩阵进行非相干叠加处理;
对所述非相干叠加处理后的结果进行恒虚预警检测;
将恒虚预警检测后的结果进行目标检测。
与现有技术相比,本申请的有益效果是:本申请通过二维自适应滤波器来减少汽车雷达信号中有效目标周围的噪声;特别是在距离多普勒图中(range Doppler map)大大提高目标的幅度与周围噪声幅值的对比度。通过增加这种对比度,目标信号功率与周围噪声的比率增加,从而让基于幅值的目标检测算法,如恒虚警检测法(CFAR)可以更有效地检测到雷达回波较弱的目标,从而提高了雷达的目标检测能力。有了较多的散射点被检测到,通过聚类的方法,雷达便可以进一步感知目标的尺寸,进而提供更多的信息,以帮助目标分类。同时,本申请是在软件中完成,因而无需进行昂贵的硬件修改,从而使其更易于在现有算法之上实施。
附图说明
图1为本申请实施例的方法的示意图。
图2为本申请实施例的线性调频信号的回波混频信号的原始数据矩阵的示意图。
图3为本申请实施例的带有二维自适应滤波器的一般雷达信号处理框图。
图4为本申请实施例的用于过滤具有多个通道的雷达数据的二维自适应滤波器框图。
图5为本申请实施例的范围自适应滤波器结构示意图。
图6为本申请实施例的复值多普勒滤波器结构示意图。
图7为本申请实施例的噪声抑制的模型
图8为本申请实施例的滤波过程的图示。a)带有较强噪声的RD。b)径向距离滤波器幅度响应。c)径向距离滤波之后的RD图。d)多普勒滤波器幅度响应。e)经过径向距离滤波和多普勒滤波器之后的RD图。
具体实施方式
下面结合具体实施方式对本申请作进一步的说明。
本申请实施例的附图中相同或相似的标号对应相同或相似的部件;在本申请的描述中,需要理解的是,若有术语“上”、“下”、“左”、“右”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本申请和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此附图中描述位置关系的用语仅用于示例性说明,不能理解为对本专利的限制。
此外,若有“第一”、“第二”等术语仅用于描述目的,主要是用于区分不同的装置、元件或组成部分(具体的种类和构造可能相同也可能不同),并非用于表明或暗示所指示装置、元件或组成部分的相对重要性和数量,而不能理解为指示或者暗示相对重要性。
在如图1所述的实施例中,本申请提供了一种改进调频连续波雷达目标检测的降噪方法,应用于汽车电子产品中,本方法包括:
100,获取并处理雷达的回波混频信号,生成原始数据矩阵;在步骤100中,获取所述回波混频信号,包括:所述雷达发射连续的线性调频信号,所述线性调频信号的回波被雷达接收后与发射线性调频信号混频的信号,经过ADC采样后,生成回波混频信号。处理所述回波混频信号,生成原始数据矩阵,包括:根据所述回波混频信号,建立原始数据矩阵Sn(k,l),其中,n为连续波雷达的第n个接收天线,l为第l个线性调频信号,k为在第l个线性调频信号的回波信号上的第k个采样点。
200,将所述原始数据矩阵通过二维自适应滤波器进行滤波;在步骤200中,所述二维自适应滤波器包括对径向距离维滤波的径向距离滤波器、对多普勒维滤波的多普勒滤波器。所述将所述原始数据矩阵通过二维自适应滤波器进行滤波,包括:通过所述径向距离滤波器对距离维的径向距离的信号进行滤波;通过所述多普勒滤波器对多普勒维的多普勒信号进行滤波。
300,对滤波后的所述原始数据矩阵进行处理;在步骤300中,对滤波后的所述原始数据矩阵进行处理,包括:将所述各通道的线性调频信号矩阵进行非相干叠加处理;对所述非相干叠加处理后的结果进行恒虚预警检测;将恒虚预警检测后的结果进行目标检测。
在本实施例中,本申请通过二维自适应滤波器来减少汽车雷达信号中有效目标周围的噪声;特别是在距离多普勒图中大大提高目标的幅度与周围噪声幅值的对比度。通过增加这种对比度,目标信号功率与周围噪声的比率增加,从而让基于幅值的目标检测算法,如恒虚警检测法可以更有效地检测到雷达回波较弱的目标,从而提高了雷达的目标检测能力。有了较多的散射点被检测到,通过聚类的方法,雷达便可以进一步感知目标的尺寸,进而提供更多的信息,以帮助目标分类。同时,本申请是在软件中完成,因而无需进行昂贵的硬件修改,从而使其更易于在现有算法之上实施。
在一些实施例中,获取所述回波混频信号,包括:所述雷达发射连续的线性调频信号,所述线性调频信号的回波被雷达接收后与发射线性调频信号混频的信号,经过ADC采样后,生成回波混频信号。在本实施例中,雷达为调频连续波雷达,调频连续波雷达是指发射频率受特定信号调制的连续波雷达。调频连续波雷达通过比较任意时刻回波信号频率与此时刻发射信号的频率的之差方法来得到目标的距离信息,距离正比于两者的频率差。目标的径向速度由获得的多普勒频率呈线性关系。与其他测距测速雷达相比,调频连续波雷达的结构更简单。本申请的回波混频信号为线性调频信号的回波被雷达接收后与发射线性调频信号混频的信号,经过模/数转换器ADC采样后,生成离散的回波混频信号。
在一些实施例中,处理所述回波混频信号,生成原始数据矩阵,包括:根据所述回波混频信号,建立原始数据矩阵Sn(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+I的滤波器系数,即a(k)=[a0(k),a1(k),…,aM(k)];
所述第l个线性调频信号的回波混频信号的原始数据矩阵第k个采样点以及其前M+1个采样点构成中间输入样本它们和第一输入参数a(k)输入至径向距离滤波器中,得到真实信号在第k个采样时刻的估计值y(k),y(k)与x(k,l)之差生成误差信号en;
根据所述误差信号en和中间输入样本通过自适应算法更新在第k+1时刻时的第一输入参数a(k+1);
通过第一输入参数a(k+1)对所述线性调频信号的回波混频信号进行迭代计算,输出初次滤波的线性调频信号矩阵。
图5为径向距离滤波器的结构示意图,线性调频信号矩阵x(k,l)=s(k,l)+n(k,l),其中s(k,l)为真实信号,n(k,l)为噪声,yn(k)是真实信号s(k,l)的估计值。第n个信道的滤波器模块(图5中虚线框内的内容)输出参数en和作为自适应算法的输入。参数en是误差信号,它近似于通道中的噪声。参数/>是中间输入样本。这两个参数是计算将在下一次迭代中使用的滤波器系数所必需的。第一输入参数a(k)是在自适应算法中得到更新的,即适应过程。滤波器系数α在每个线性调频脉冲的开始时被初始化,即使a0(k)的值为1。随着来自线性调频脉冲的回波被滤波,径向距离滤波器的第一输入参数a(k)被更新。移位寄存器包含/>值,这些值是算法在适应中使用的中间输入样本。通过第一输入参数a(k)对线性调频信号矩阵进行迭代计算,输出进行初次滤波的线性调频信号矩阵。所述通过所述径向距离滤波器对距离维的径向距离的信号进行滤波之后,还包括:将所述滤波后的线性调频信号矩阵通过快速傅里叶变换获得径向距离信息。在本实施例中,发射信号是线性调频信号,当遇到目标时其回波是一个延迟的线性调频信号,它与发射信号的混频信号是一个正弦信号,傅里叶变换便得到这个正弦波的频率,该频率与目标的距离有直接对应关系。因此通过傅里叶变换获得目标的径向距离信息。
在上述实施例的一种实施方式中,参见图6,所述通过所述多普勒滤波器对线性调频信号矩阵的相邻所述线性调频信号进行滤波,包括:
将多普勒滤波器的第二输入参数b(l)初始化,b(l)是长度为P+1的滤波器系数b(l)=[b0(l),b1(l),…,bP(l)];
所述第l个线性调频信号的回波混频信号的原始数据矩阵第k个采样点以及其前P+1个线性调频信号的回波混频信号构成中间输入样本它们和第二输入参数b(l)输入至多普勒滤波器中,得到在第k个采样时刻且第l个线性调频信号的回波混频信号真实信号的估计值y(l),y(l)与X(k,l)之差生成误差信号en(l);
根据所述误差信号en(l)和中间输入样本通过自适应算法更新第二输入参数b(l);
通过输入第二参数b(l)对所述原始数据矩阵进行迭代计算,输出再次滤波的线性调频信号矩阵。
在本实施例中,图6为多普勒滤波器的结构示意图,线性调频信号矩阵x(k,l)=S(k,l)+n(k,l),其中S(k,l)为真实信号经过径向距离的傅里叶变换,n(k,l)为噪声,yn(l)是真实信号S(k,l)的估计值。第n个信道的滤波器模块(图6中虚线框内的内容)输出参数en和作为自适应算法的输入。参数en是误差信号,它近似于通道中的噪声。参数/>是中间输入样本。这两个参数是计算将在下一次迭代中使用的滤波器系数所必需的。第二输入参数b(l)是在自适应算法中得到更新的,即适应过程。滤波器系数b在每个线性调频脉冲的开始处被初始化,即使b0(l)的值为1。随着来自线性调频脉冲的回波被滤波,多普勒滤波器的第二输入参数b(l)被更新。移位寄存器包含/>值,这些值是算法在适应中使用的中间输入样本。通过第二输入参数b(l)对线性调频信号矩阵进行迭代计算,输出进行再次滤波的线性调频信号矩阵。所述通过所述多普勒滤波器对多普勒维的多普勒信号进行滤波之后,还包括:将所述线性调频信号矩阵通过快速傅里叶变换计算得到多普勒信息。
在一些实施例中,所述自适应算法为归一化最小均方算法或时变最小均方算法或最小二乘法。所述根据所述误差信号en和中间输入样本通过自适应算法更新径向距离滤波器在时刻k+1的参数a(k+1),包括:
通过计算求得,其中,Δ是步进常数。
图7的框图公开了一个基于最小均方的自适应算法的示例。首先,信号y是真实信号s的估计值。目标是在最小均方意义上最小化误差e:
e=y-(s+n)
对e取平方:
e2=(y-(s+n))2
=y2+s2+n2-2(sy+ny)+2sn
求e2的平均值E(e2):
E(e2)=E(y2)+E(s2)+E(n2)-2E(sy+ny)+2E(sn).
这里假设s和n不相关。当s是正弦信号,而n是宽带噪声时,它们之间时非相关的。因此上面的假设是成立的。E(e2)的最小值为
E(e2)min=E(y2)min+E(s2)+E(n2)-2E(sy+ny)max+0.
由于s和n不相关,所以E(sn)=0。当信号功率大于噪声功率时,且y很接近于真实信号s时,则E(e2)就达到了最小值。
可见,该滤波器实际上是一个带通滤波器,它使得输出y接近真实信号s,因此从另一种意义上讲,噪声的功率被抑制了。如果s是单频正弦信号,那么这个滤波器时一个非常窄的带通滤波器,只允许s通过。而且E(e2)_min等于实际通道的噪声功率。这意味着误差信号e对噪声做到了最理想的估计,因而能被有效的抑制。
对于归一化的LMS,系数更新如下:其中Δ是步进常数。
同理,第二输入参数b(l)可通过进行计算。
在一些实施中,所述对滤波后的所述原始数据矩阵进行处理,包括:
将所述各通道的线性调频信号矩阵进行非相干叠加处理;在本实施例中,将所述N个通道的傅里叶变换后得到的距离一多普勒图(range-doppler map)进行合并获得一个合并后的距离一多普勒图。一种有效快捷的方法是对这N个通道的数据进行非相干叠加处理,即non-coherent sum,因此N个矩阵的数据就合并为一个矩阵。
对所述非相干叠加处理后的结果进行恒虚预警检测;在本实施例中,在雷达信号检测中,当外界干扰强度不断变化时,雷达能自动调整其灵敏度,使雷达的虚警概率保持不变,这种特性称为恒虚警率特性。也就是说,进行检测的门限值不是事先固定的,而是随外界强度变化而做相应的调整,因此它是通过找自适应门限值来进行目标检测的。
将恒虚预警检测后的结果进行目标检测。在本实施例中,在雷达信号检测中,当外界干扰强度不断变化时,雷达能自动调整其灵敏度,使雷达的虚警概率保持不变,这种特性称为恒虚警率特性。也就是说,进行检测的门限值不是事先固定的,而是随外界强度变化而做相应的调整,因此它是通过找自适应门限值来进行目标检测的。
参见图3,本申请公开了一种改进调频连续波雷达目标检测的降噪方法,通过调频连续波获取包含信号和噪音的线性调频信号矩阵,先通过径向距离滤波器对线性调频信号矩阵进行初步滤波,再通过多普勒滤波器对其进行再次滤波后。如果有多个通道的情况,可以将滤波后的多通道线性调频信号矩阵进行非相干叠加,再经过恒虚预警检测和目标检测,获取目标列表。参见图8,a)含有较高噪声的RD图。b)径向距离滤波器幅度响应。c)径向距离过滤器之后的RD图。d)多普勒滤波器幅度响应。e)径向距离和多普勒滤波器之后的RD图。通过图8可清晰看到目标的SNR显著增加,因此提高了这种弱目标的检测能力。本申请通过二维自适应滤波器来减少汽车雷达信号中有效目标周围的噪声;大大提高目标幅度与周围多普勒RD图中周围噪声的对比度。通过增加这种对比度,目标信号功率与周围噪声的比率增加,从而让CFAR可以更有效地检测目标。使雷达可以灵敏地观察到聚类,从而为目标的可能尺寸提供更多信息,以帮助目标分类。同时,本申请在软件中完成,而无需进行昂贵的硬件修改,从而使其更易于在现有算法之上实施。
显然,本申请的上述实施例仅仅是为清楚地说明本申请所作的举例,而并非是对本申请的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本申请的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本申请权利要求的保护范围之内。
Claims (10)
1.一种改进调频连续波雷达目标检测的降噪方法,其特征在于,应用于设置有雷达的汽车电子产品中,所述方法包括:
获取并处理雷达的回波混频信号,生成原始数据矩阵;
将所述原始数据矩阵通过二维自适应滤波器进行滤波;
对滤波后的所述原始数据矩阵进行处理;
其中,所述二维自适应滤波器包括对径向距离维滤波的径向距离滤波器、对多普勒维滤波的多普勒滤波器;
其中,径向距离滤波器的滤波过程包括:
将径向距离滤波器的第一输入参数a(k)初始化,a(k)是长度为M+1的滤波器系数,即a(k)=[a0(k),a1(k),…,aM(k)];
第l个线性调频信号的回波混频信号的原始数据矩阵第k个采样点以及其前M+1个采样点构成中间输入样本它们和第一输入参数a(k)输入至径向距离滤波器中,得到真实信号在第k个采样时刻的估计值y(k),y(k)与x(k,l)之差生成误差信号en(k),其中,线性调频信号矩阵x(k,l)=s(k,l)+n(k,l),s(k,l)为真实信号,n(k,l)为噪声,k为采样点数,l为线性调频信号个数;
根据所述误差信号en和中间输入样本通过自适应算法更新径向距离滤波器在第k+1时刻时的第一输入参数a(k+1);
通过第一输入参数a(k+1)对所述线性调频信号的回波混频信号进行迭代计算,输出初次滤波的线性调频信号矩阵。
2.根据权利要求1所述的一种改进调频连续波雷达目标检测的降噪方法,其特征在于,获取所述回波混频信号,包括:
所述雷达发射连续的线性调频信号,所述线性调频信号的回波被雷达接收后与发射线性调频信号混频的信号,经过ADC采样后,生成回波混频信号。
3.根据权利要求1所述的一种改进调频连续波雷达目标检测的降噪方法,其特征在于,处理所述回波混频信号,生成原始数据矩阵,包括:
根据所述回波混频信号,建立原始数据矩阵Sn(k,l),其中,n为连续波雷达的第n个接收天线,l为第l个线性调频信号,k为在第l个线性调频信号的回波信号上的第k个采样点。
4.根据权利要求1所述的一种改进调频连续波雷达目标检测的降噪方法,其特征在于,所述将所述原始数据矩阵通过二维自适应滤波器进行滤波,包括:
通过所述径向距离滤波器对距离维的径向距离的信号进行滤波;
通过所述多普勒滤波器对多普勒维的多普勒信号进行滤波。
5.根据权利要求1所述的一种改进调频连续波雷达目标检测的降噪方法,其特征在于,通过所述径向距离滤波器对线性调频信号矩阵的同一所述线性调频信号进行滤波之后,还包括:
将所述初次滤波后的线性调频信号矩阵通过快速傅里叶变换获得径向距离信息。
6.根据权利要求4所述的一种改进调频连续波雷达目标检测的降噪方法,其特征在于,通过所述多普勒滤波器对线性调频信号矩阵的相邻所述线性调频信号进行滤波,包括:
将多普勒滤波器的第二输入参数b(l)初始化,b(l)是长度为P+1的滤波器系数,即b(l)=[b0(l),b1(l),…,bP(l)];
第l个线性调频信号的回波混频信号的原始数据矩阵第k个采样点以及其前P+1个采样点构成中间输入样本它们和第二输入参数b(l)输入至多普勒滤波器中,得到在第k个采样时刻且第l个线性调频信号的回波混频信号经过径向距离的傅里叶变换的真实信号S(k,l)的估计值y(l),y(l)与X(k,l)之差生成误差信号en(l),线性调频信号矩阵X(k,l)=S(k,l)+n(k,l);
根据所述误差信号en(l)和中间输入样本通过自适应算法更新第二输入参数b(l);
通过第二输入参数b(l)对所述原始数据矩阵进行迭代计算,输出再次滤波的线性调频信号矩阵。
7.根据权利要求6所述的一种改进调频连续波雷达目标检测的降噪方法,其特征在于,所述通过所述多普勒滤波器对多普勒维的多普勒信号进行滤波之后,还包括:
将再次滤波后的所述线性调频信号矩阵通过快速傅里叶变换计算得到多普勒信息。
8.根据权利要求6所述的一种改进调频连续波雷达目标检测的降噪方法,其特征在于,所述自适应算法可以为归一化最小均方算法或时变最小均方算法或最小二乘法。
9.根据权利要求8所述的一种改进调频连续波雷达目标检测的降噪方法,其特征在于,所述根据所述误差信号en(k)和中间输入样本通过自适应算法更新径向距离滤波器在第k+1时刻时的第一输入参数a(k+1),包括:
通过计算求得,其中,Δ是步进常数,e(k)为第k个采样时刻的误差信号,/>为第k个采样时刻的中间输入样本。
10.根据权利要求1所述的一种改进调频连续波雷达目标检测的降噪方法,其特征在于,所述对滤波后的所述原始数据矩阵进行处理,包括:
将各通道的线性调频信号矩阵进行非相干叠加处理;
对所述非相干叠加处理后的结果进行恒虚预警检测;
将恒虚预警检测后的结果进行目标检测。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011051384.1A CN112630768B (zh) | 2020-09-29 | 2020-09-29 | 一种改进调频连续波雷达目标检测的降噪方法 |
PCT/CN2020/139042 WO2022068097A1 (zh) | 2020-09-29 | 2020-12-24 | 一种改进调频连续波雷达目标检测的降噪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011051384.1A CN112630768B (zh) | 2020-09-29 | 2020-09-29 | 一种改进调频连续波雷达目标检测的降噪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112630768A CN112630768A (zh) | 2021-04-09 |
CN112630768B true CN112630768B (zh) | 2024-04-02 |
Family
ID=75302740
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011051384.1A Active CN112630768B (zh) | 2020-09-29 | 2020-09-29 | 一种改进调频连续波雷达目标检测的降噪方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN112630768B (zh) |
WO (1) | WO2022068097A1 (zh) |
Families Citing this family (8)
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 | 西安瀚博电子科技有限公司 | 一种基于自适应重构滤波的信号特征检测方法 |
CN115166664A (zh) * | 2022-06-12 | 2022-10-11 | 北京理工大学 | 一种基于二维匹配滤波结果的快速自适应旁瓣抑制方法 |
CN116256714B (zh) * | 2022-12-31 | 2023-10-13 | 南京航空航天大学 | 基于雷达成像的电磁散射特性测量方法及装置 |
CN116955444A (zh) * | 2023-06-15 | 2023-10-27 | 共享易付(广州)网络科技有限公司 | 基于大数据分析的采集噪声点挖掘方法及系统 |
CN116577738B (zh) * | 2023-07-12 | 2023-09-12 | 南京隼眼电子科技有限公司 | 雷达同频抗干扰方法、装置、雷达设备及存储介质 |
CN116643250B (zh) * | 2023-07-27 | 2023-10-20 | 西安电子科技大学 | 基于稀疏线性调频信号的宽带雷达高速目标相参积累方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103339525A (zh) * | 2010-12-21 | 2013-10-02 | 塔莱斯公司 | 用于监测地势变化的方法和设备 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2954516B1 (fr) * | 2009-12-18 | 2012-09-21 | Thales Sa | Procede de detection radar, notamment pour des radars aeroportes mettant en oeuvre une fonction de detection et d'evitement d'obstacles |
US9594159B2 (en) * | 2013-07-15 | 2017-03-14 | Texas Instruments Incorporated | 2-D object detection in radar applications |
US9784828B2 (en) * | 2014-08-27 | 2017-10-10 | Texas Insturments Incorporated | FMCW doppler processing algorithm for achieving CW performance |
US10451723B2 (en) * | 2016-12-20 | 2019-10-22 | National Chung-Shan Institute Of Science & Technology | Signal processing apparatus of a continuous-wave (CW) radar sensing system |
DE102018108648B4 (de) * | 2018-04-11 | 2021-11-11 | Infineon Technologies Ag | Fmcw radar mit störsignalunterdrückung |
-
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 (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103339525A (zh) * | 2010-12-21 | 2013-10-02 | 塔莱斯公司 | 用于监测地势变化的方法和设备 |
Also Published As
Publication number | Publication date |
---|---|
CN112630768A (zh) | 2021-04-09 |
WO2022068097A1 (zh) | 2022-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112630768B (zh) | 一种改进调频连续波雷达目标检测的降噪方法 | |
US11016171B2 (en) | Radar sensing with phase correction | |
KR102186191B1 (ko) | 간섭이 억제된 레이더 감지 | |
US10996311B2 (en) | Distributed radar sensor system | |
CN111352102B (zh) | 一种基于调频连续波雷达的多目标个数检测方法及装置 | |
EP1672379B1 (en) | System and method for reducing a radar interference signal | |
CN109444820B (zh) | 杂波和干扰共存时多通道雷达先干扰抑制后目标检测方法 | |
EP2097769B1 (en) | System and method for reducing the effect of a radar interference signal | |
US20200124699A1 (en) | Fmcw radar with interference signal suppression | |
US11592520B2 (en) | FMCW radar with interfering signal suppression in the time domain | |
CN111693946A (zh) | 具有借助于人工神经网络的干扰信号抑制的fmcw雷达 | |
US10656248B2 (en) | Radar post processing for sidelobe suppression | |
US10627481B2 (en) | Multi-resolution doppler processing | |
US11681011B2 (en) | Detection of interference-induced perturbations in FMCW radar systems | |
US20230168367A1 (en) | CFAR Phased Array Pre-Processing Using Noncoherent and Coherent Integration in Automotive Radar Systems | |
CN109507654B (zh) | 一种基于ls的复杂环境下相位信息计算方法 | |
CN109581366B (zh) | 一种基于目标导向矢量失配的离散旁瓣杂波识别方法 | |
CN111812634A (zh) | 警戒线目标监测方法、装置和系统 | |
CN111175714A (zh) | 可抑制雷达近距离谐波的辅助驾驶方法和存储介质 | |
CN110907930B (zh) | 一种基于角度估计的车载雷达目标检测估计方法及装置 | |
US10845475B2 (en) | Method of measuring azimuth of radar target | |
US11300655B2 (en) | Object detection using radar | |
Cha et al. | Implementation of high-resolution angle estimator for an unmanned ground vehicle | |
CN111175715A (zh) | 可抑制雷达近距离谐波的辅助驾驶系统、方法 | |
CN113325421B (zh) | 用于提供所述至少一个目标信息的方法和计算机可读介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |