CN101320086A - 一种多普勒测速激光雷达的回波信号处理装置和方法 - Google Patents

一种多普勒测速激光雷达的回波信号处理装置和方法 Download PDF

Info

Publication number
CN101320086A
CN101320086A CNA2008101158144A CN200810115814A CN101320086A CN 101320086 A CN101320086 A CN 101320086A CN A2008101158144 A CNA2008101158144 A CN A2008101158144A CN 200810115814 A CN200810115814 A CN 200810115814A CN 101320086 A CN101320086 A CN 101320086A
Authority
CN
China
Prior art keywords
frequency
echoed signal
signal
sampling rate
module
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
Application number
CNA2008101158144A
Other languages
English (en)
Other versions
CN101320086B (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.)
Beihang University
Original Assignee
Beihang 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 Beihang University filed Critical Beihang University
Priority to CN2008101158144A priority Critical patent/CN101320086B/zh
Publication of CN101320086A publication Critical patent/CN101320086A/zh
Application granted granted Critical
Publication of CN101320086B publication Critical patent/CN101320086B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Optical Radar Systems And Details Thereof (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明提供了一种多普勒测速激光雷达的回波信号处理装置和方法。该装置主要包括:信号存储模块,用于存储被测目标的回波信号,该回波信号为经过A/D转换后的数字信号;DSP模块,用于读取所述信号存储模块中的回波信号并进行FFT,求得信号的功率谱,根据所述功率谱得到粗估计频率值;对所述功率谱进行逆FFT得到所述回波信号的自相关函数,根据该自相关函数得到校正频率值,将所述粗估计频率值和校正频率值相加得到所述回波信号的多普勒频率,根据该多普勒频率计算出所述被测目标的运动速度。本发明将自相关技术与基于频偏校正的频率估计改进算法相结合,可以有效地计算出被测目标的速度。

Description

一种多普勒测速激光雷达的回波信号处理装置和方法
技术领域
本发明涉及光学和电子应用领域,尤其涉及一种多普勒测速激光雷达的回波信号处理装置和方法。
背景技术
激光雷达是以激光波束为信息载体的雷达。激光雷达不仅可精确测距,而且还能精确测速、精确跟踪和高分辨率成像。激光雷达的应用场合非常广泛,可用于军事、航天和工业生产等多个领域。
多普勒测速激光雷达是以光学多普勒效应为基础的一种激光雷达,可以用来测量液体、气体的流动速度或固体目标的运动速度,在科学和工业技术中有非常广泛的应用。
多普勒测速激光雷达采用主动探测方式,系统发射一定波形的激光照射到被测目标上,通过对被测目标反射或散射的回波信号进行检测,获得被测目标的速度信息。多普勒测速激光雷达的探测形式可分为直接探测和相干探测两种。
在外差相干探测多普勒测速激光雷达中,在接收到回波信号后,对该回波信号进行放大、滤波等处理后得到输出的多普勒信号,该多普勒信号频率ωD与被测目标的被测速度
Figure A20081011581400061
之间存在如下的计算公式:
公式1
上述公式中的θ为光束传播方向与被测目标的运动方向之间的夹角。
现有技术中的第一种获取上述多普勒信号频率的方法为:计数法(过零检测法)。该方法的主要工作原理为:从一个过零点开始以非常高的时钟脉冲对多普勒信号进行计数,来获取多普勒信号中的N个信号周期所需的时间(固定周期数计数法);或者,获取多普勒信号中的一定时间长度内的信号周期个数(固定闸门时间计数法)。然后,根据获取的上述时间或信号周期个数信息得到多普勒信号频率。
现有技术中的第一种获取上述多普勒信号频率的方法的缺点为:该方法对噪声比较敏感,不适用于噪声较强的情况。并且测量精度要受计数误差、晶体振荡频率误差的影响。
现有技术中的第二种获取上述多普勒信号频率的方法为:双钳位自相关法。该方法的主要工作原理为:获取回波信号的自相关函数:
R xx ( j ) = 1 R ( 0 ) Σ i = 0 n - 1 x ( i ) x ( i + j )
上述公式中的x(i)和x(i+j)分别为回波信号x(t)的离散N位和M位采样值。R(j)是j个采样延时的自相关函数。若N和M大于1,则为多位自相关。若N大于1,而M等于1,则是单钳位自相关,若N和M都等于1,则是双钳位自相关。一个无量纲的理想自相关系数1/R(0)的值应该为1。
由于信号和噪声是相互独立的过程,根据相关函数的定义,信号只与信号本身相关,而与噪声不相关。因此,该方法计算出双钳位自相关情况下的上述回波信号的自相关函数的频率,根据自相关函数的性质,该频率与上述多普勒信号频率相同。
现有技术中的第二种获取上述多普勒信号频率的方法的缺点为:该方法不适用于噪声较强的情况,对于信噪比较低的应用场合,获得的速度测量精度较差。
发明内容
本发明实施例的目的是提供一种多普勒测速激光雷达的回波信号处理装置和方法,从而可以解决现有技术方案中噪声门限高、测量结果的精度较差的缺点。
本发明实施例的目的是通过以下技术方案实现的:
一种多普勒测速激光雷达的回波信号处理装置,包括:
信号存储模块,用于存储被测目标的回波信号,该回波信号为经过模数A/D转换后的数字信号;
数据信号处理DSP模块,用于读取所述信号存储模块中的回波信号并进行快速傅立叶变换FFT,求得信号的功率谱,根据所述功率谱得到粗估计频率值;对所述功率谱进行逆FFT得到所述回波信号的自相关函数,根据该自相关函数得到校正频率值,将所述粗估计频率值和校正频率值相加得到所述回波信号的多普勒频率,根据该多普勒频率计算出所述被测目标的运动速度。
一种多普勒测速激光雷达的回波信号处理方法,包括:
对被测目标返回的经过A/D转换后的回波信号进行FFT,根据FFT后的回波信号的功率谱获取所述回波信号的粗估计频率值,对所述功率谱进行逆FFT得到所述回波信号的自相关函数,根据该自相关函数得到所述回波信号的校正频率值;
将所述信号频率粗估计值和校正频率值相加得到所述回波信号的多普勒频率,根据该多普勒频率计算出所述被测目标的运动速度。
由上述本发明实施例提供的技术方案可以看出,本发明实施例将自相关技术与基于频偏校正的频率估计改进算法相结合,可以有效地计算出被测目标的速度,并且保证测量结果的精度,通过采取减少算法计算量的措施,提高系统的实时性。
附图说明
图1为本发明实施例所述多普勒测速激光雷达的回波处理装置的结构示意图;
图2为本发明实施例所述多普勒测速激光雷达的回波处理方法的一种处理流程图;
图3为本发明实施例所述频偏校正改进算法的处理流程图。
具体实施方式
在本发明实施例中,预先设定不同的滤波频率范围分别对应不同的采样率,并且相邻的滤波频率范围之间互相重叠,将各个滤波频率范围和对应的采样率进行关联保存。
根据所述放大后的回波信号的初始频率值,查询预先设定的滤波频率范围和采样率之间的对应关系,获取所述初始频率值对应的滤波频率范围和采样率,按照该滤波频率范围对所述回波信号进行带通滤波,按照所述采样率对所述滤波处理后的回波信号进行采样,将采样获取的模拟回波信号转换为数字信号,将该数字信号进行存储。
对所述存储的数字的回波信号进行读取,对读取的的回波信号进行2倍采样点的FFT(快速傅立叶变换,Fast Fourier Transform),求取FFT后的回波信号的功率谱,搜索该功率谱的最大谱峰对应的离散谱线,根据该离散谱线、采样率和采样点数计算出所述回波信号的粗估计频率值;
对所述功率谱进行逆FFT得到所述回波信号的自相关函数,利用所述信号频率粗估计值对该自相关函数进行移频处理,求取移频后信号的自相关函数在指定点的自相关值,根据该自相关值计算出所述回波信号的校正频率值。
将所述信号频率粗估计值和校正频率值相加得到所述回波信号的多普勒频率,根据该多普勒频率计算出所述被测目标的运动速度。
进一步地,判断所述多普勒频率与当前滤波频率范围、采样率是否对应,如果是,则按照所述当前滤波频率范围、采样率继续采集后续回波信号;否则,按照重新选择的所述多普勒频率对应的滤波频率范围和采样率继续采集后续回波信号。
本发明实施例所述多普勒测速激光雷达的回波处理装置的结构示意图如图1所示,包括如下模块:放大处理模块、滤波处理模块、A/D转换模块、信号存储模块、DSP(数据信号处理,Digital Signal Processing)模块、显示模块等。
放大处理模块,用于接收多普勒测速激光雷达中的探测器输出的回波信号,将该回波信号放大到适合A/D转换的幅度,将放大后的回波信号传输给滤波处理模块。探测器输出的回波信号具有较大的带宽,大约为25MHz,且含有较大的噪声,因此,该放大处理模块必须包含大带宽、低噪声的放大电路。
滤波处理模块,用于接收放大处理模块传输的放大后的回波信号,对该回波信号中的噪声、杂波和干扰进行滤波处理,该噪声包括散粒噪声、热噪声和量化噪声等多种噪声。将进行了滤波处理后的回波信号传输给A/D转换模块。滤波处理模块在进行滤波处理时,可以根据所述回波信号的粗测频率值查询预先设定的滤波频率范围和采样率之间的对应关系,获取所述频率值对应的滤波频率范围,按照该滤波频率范围对所述回波信号进行带通滤波。包括:采样率和滤波频率范围分段设置模块。
A/D转换模块,用于接收滤波处理模块输出的进行了滤波处理后的模拟的回波信号,根据该模拟的回波信号的频率值预先设定的滤波频率范围和采样率之间的对应关系,获取所述频率值对应的采样率,按照该采样率对所述模拟的回波信号进行采样,将采样获取的模拟的回波信号转换为数字的回波信号,将该数字的回波信号传输给信号存储模块。包括:采样率和滤波频率范围分段设置模块。
信号存储模块,用于接收A/D转换模块传输的数字的回波信号。该信号存储模块可以通过FPGA(现场可编程门阵列,Field Programmable GateArray)来实现。
DSP模块,用于读取所述信号存储模块中的回波信号并进行FFT,根据FFT后的回波信号的功率谱谱峰值获取粗估计频率值,对所述功率谱进行逆FFT得到所述回波信号的自相关函数,根据该自相关函数得到频率较正值,将所述粗估计频率值和较正频率值相加得到所述回波信号的多普勒频率,根据该多普勒频率计算出所述被测目标的运动速度。包括:信号读取模块、粗估计频率值计算模块、校正频率值计算模块、运动速度获取模块和当前滤波频率范围和采样率控制模块。
显示模块,用于将所述DSP模块所获取的被测目标的运动速度通过液晶等方式进行显示。
上述的滤波处理模块和A/D转换模块中的采样率和滤波频率范围分段设置模块,用于设定不同的滤波频率范围分别对应不同的采样率,并且各段滤波频率范围中的相邻两段之间互相重叠,将各个滤波频率范围和对应的采样率进行关联保存。
上述的DSP模块中的信号读取模块,用于对所述信号存储模块中存储的回波信号进行读取。
上述的DSP模块中的粗估计频率值计算模块,用于对所述信号读取模块所读取的回波信号进行2倍采样点的FFT,求取回波信号的功率谱,搜索该功率谱的最大谱峰对应的离散谱线,根据该离散谱线、采样率和采样点数计算出粗估计频率值;
上述的DSP模块中的校正频率值计算模块,用于对所述功率谱进行逆FFT得到所述回波信号的自相关函数,利用所述信号频率粗估计值对该自相关函数进行移频处理,求取移频后的自相关函数在指定点的自相关值,根据该自相关值计算出校正频率值;
上述的DSP模块中的运动速度获取模块,用于将所述粗估计频率值和校正频率值相加得到所述回波信号的多普勒频率,根据该多普勒频率以及激光束与被测目标的运动方向之间的夹角计算出所述被测目标的运动速度。
上述的DSP模块中的当前滤波频率范围和采样率控制模块,用于判断运动速度计算模块获取的多普勒频率与当前滤波频率范围和采样率是否对应,如果是,则维持当前滤波频率范围和采样率不变;否则,根据所述运动速度计算模块获取的多普勒频率查询预先设定的滤波频率范围和采样率的对应关系,重新确定当前滤波频率范围和采样率。
上述DSP模块可以设置在PC(个人计算机,Personal Computer)机中,信号存储模块通过USB接口或RS-232接口将存储的回波信号传输给PC机,通过PC中DSP模块进行上述计算被测目标的运动速度的处理过程。
本发明实施例所述多普勒测速激光雷达的回波处理方法的一种处理流程如图2所示,包括如下步骤:
步骤2-1、对探测到的回波信号进行放大,根据放大后的回波信号的频率值和预先设定的滤波频率范围和采样率的对应关系,获取放大后的回波信号对应的滤波频率范围,并进行滤波处理。
多普勒测速激光雷达中的探测器探测到目标的回波信号后,将该回波信号传输给放大处理模块。该回波信号具有较大的带宽,大约为25MHz,且含有较大的噪声。
放大处理模块将上述回波信号放大到适合A/D转换的幅度,将放大后的回波信号传输给滤波处理模块。上述放大处理模块必须包含大带宽、低噪声的放大电路。
如果在整个测量频率范围内采用同一种采样率对上述回波信号进行采样,则在对测量频率范围的下限频率进行离散频谱分析时所需的FFT(快速傅立叶变换)点数将非常大。为解决该问题,在本发明实施例中,预先设定不同滤波频率范围对应不同的采样率,对采样率按照不同滤波频率范围分6段设置,该采样率分段设置的情况如下述表1所示,系统上电时,默认采样率为80MSPS,对应的带通滤波频率范围为1kHz-25MHz:
表1:采样率分段设置
  段号   带通滤波频率范围(Hz)   采样率(SPS)
  1   1k~10k   40k
  2   8k~80k   320k
  3   70k~700k   2.8M
  4   600k~6M   24M
  5   5M~12M   48M
  6   10M~25M   80M
  默认   1k~25M   80M
将每个滤波频率范围对应的频率分辨率与频率范围的下限值的比值设定为固定值,这样在后续的信号频谱分析中,就可以采用通用的算法程序在每个滤波频率范围,对采样获取的回波信号分别进行离散频谱分析处理,从而降低频谱分析的难度,并且可以使信号处理算法的相对精度稳定在一个近似相同的水平。相邻的两个滤波频率范围之间存在相互重叠,可以避免在频率分界点附近,采样率的频繁切换。
滤波处理模块接收到放大处理模块传输的放大后的回波信号后,按照上述表1中默认的滤波频率范围对该回波信号中的噪声、杂波和干扰进行滤波处理,该噪声包括散粒噪声、热噪声和量化噪声等多种噪声。将进行了滤波处理后的回波信号传输给A/D转换模块。
滤波处理模块在进行滤波处理时,可以按照预先设定的滤波频率范围通过两个双四选一模拟开关切换进行分段滤波处理,从而进一步地滤除非信号频段的干扰和防止频率混叠,进一步提高系统的信噪比。
步骤2-2、根据滤波处理后的回波信号的频率值和预先设定的滤波频率范围和采样率之间的对应关系,获取滤波处理后的回波信号对应的采样率,并进行采样处理,将采样后的回波信号转换为数字信号,并进行存储。
A/D转换模块接收上述进行了滤波处理后的回波信号,按照上述表1中默认的采样率对该回波信号进行采样处理,将采样后的回波信号转换为数字信号,通过DSP模块进行粗测获取该数字信号的初始频率值。根据所述初始频率值,查询上述表1所示的滤波频率范围和采样率之间的对应关系,获取所述初始频率值对应的滤波频率范围和采样率,并作为当前滤波频率范围和采样率。
按照所述当前滤波频率范围,滤波处理模块重新对上述放大处理模块传输的放大后的回波信号进行带通滤波,A/D转换模块再按照所述当前采样率对所述滤波处理后的回波信号进行采样,将采样获取的模拟回波信号转换为数字信号,将该数字信号传输给FPGA模块,并存储在FPGA模块中的RAM(随机存取存储器,Random Access Memory)中。
由于上述回波信号的频率值是变化的,当该频率值变化到另外一个频率范围内时,则滤波处理模块和A/D转换模块需要根据变化后的频率值查询上述表1所示的滤波频率范围和采样率之间的对应关系,重新获取当前滤波频率范围和采样率。
FPGA模块产生不同频率的方波信号作为A/D转换模块的时钟信号,以控制A/D转换模块在表1所述的六种采样率下完成回波信号的采集。
步骤2-3、对采样后存储的回波信号运用将自相关技术与基于频偏校正的频率估计相结合的频率估计改进算法获取多普勒信号的频率,进而计算出目标速度值。
DSP模块从FPGA模块中的RAM区读取采集数据,对采集数据运用基于频偏校正的频率估计改进算法估计出多普勒信号的频率。该频率估计改进算法的处理流程如图3所示,具体处理过程如下:
DSP模块启动一次采集过程,默认滤波频率范围为1kHz~25MHz,采样率fs为最大值(即为80M),一次从FPGA模块中的RAM区读取N点采样数据。对于较短的采样间隔T=N/fs(N为采样点数,fs为采样频率)内的采样数据,可以看作幅值恒定,则得到的离散采样序列为:
Figure A20081011581400151
上述信号x(n)中的a为采样信号幅度,f0为信号频率,φ0为信号的初始相位,u(n)为零均值的高斯白噪声。
对上述信号x(n),n=0,1,2,…N-1补N个零,并对信号x(n)做2N点FFT,求取信号的功率谱。该2N点FFT需要N(1+log2N)次复数乘法和2N(1+log2N)次复数加法。
在0~fs/2范围内搜索功率谱中的最大谱峰对应的离散谱线k0,并求得对应的频率值 f ^ 0 = k 0 f s / ( 2 N ) , 将该
Figure A20081011581400153
作为多普勒信号的粗估计频率值。上述谱峰搜索过程仅需要很小的运算量。
然后,再对功率谱做逆变换得到x(n)的自相关函数r(n),n=0,1,2,…2N-1。上述逆变换需要与做2N点FFT同样次数的乘法和加法。
由于信号和噪声是相互独立的过程,根据相关函数的定义,信号只与信号本身相关与噪声不相关,而噪声之间一般也是不相关的,r(n)与原信号x(n)具有相同的频率,因此,对采集数据进行自相关运算可以有效的去除干扰,提高信号的信噪比,而且不影响信号的频率特性。
用粗估计得到的频率值
Figure A20081011581400161
创建一个N点的新序列:
w ( n ) = e j ( - 2 π f ^ 0 n ) , n=0,1,2,…2N-1
将该新序列与自相关函数相乘,即对自相关函数进行移频处理,得到 y ( n ) = r ( n ) · w ( n ) = a 2 e j 2 π f Δ nΔ , n=0,1,2,…2N-1
设最终的多普勒信号的频率值为f0,频率偏差 f Δ = f 0 - f ^ 0 . 求取y(n)的自相关函数在m点的自相关值h:
h = ( 2 N - m ) a 4 e j 2 π f Δ nΔ
h的相位φ=2πfΔm为fΔ的一次函数,且相位φ<π,因此,利用h的反正切值可以无模糊的得到:
fΔ=φ/2πm
进而可以得到:
f 0 = f ^ 0 + f Δ
在实际应用中,可以取m=N,则:
h = Na 4 e j 2 π ( f 0 - f ^ 0 ) NΔ = r ( N ) e j ( - 2 π f ^ 0 NΔ )
其中r(N)可以根据自相关函数r(n)得出,且由FFT因子的对称性可知
Figure A20081011581400168
这样,
Figure A20081011581400169
就可以省去N次复数乘法运算。
查询上述表1获得上述f0对应的滤波频率范围和采样率。然后,按照该查询得到的滤波频率范围和采样率继续采集数据,继续运用基于频偏校正的频率估计改进算法重新获取多普勒信号的频率。判断该重新获取的频率是否上述f0对应的频率范围,如果是,则按照当前的滤波频率范围和采样率继续采集后续数据;否则,查询上述表1获得上述重新获取的频率对应的滤波频率范围和采样率,并作为当前滤波频率范围和采样率,按照该当前滤波频率范围和采样率继续采集数据。
最后,按照上述公式1,根据多普勒信号的频率f0计算出被测目标的目标速度值。
DSP模块将得到的目标速度值输出给液晶显示模块进行显示,并且将目标速度值输出给FPGA模块进行保存。
上述DSP模块可以设置在PC机中,此时,PC机通过通信接口从FPGA模块读取数据,运用上述基于频偏校正的频率估计改进算法获取多普勒信号的频率,计算出被测目标的目标速度值。
综上所述,本发明实施例将自相关技术与基于频偏校正的频率估计算法相结合,可以有效地计算出被测目标的运动速度,保证在整个频段范围内的测量结果的误差比较平缓,性能稳定,能够适用于噪声较强的情况。并有效利用进行频率粗估计时得到的信号功率谱,减少整个算法的计算量。
本发明实施例采用分段滤波与变采样率相结合,可以有效的提高信噪比,并且每个分段的频率分辨率对本段的频率下限的比值都近似处于相同的水平,使得在DSP模块中可采用通用的算法进行处理,达到相同的估计精度。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。

Claims (10)

1、一种多普勒测速激光雷达的回波信号处理装置,其特征在于,包括:
信号存储模块,用于存储被测目标的回波信号,该回波信号为经过模数A/D转换后的数字信号;
数据信号处理DSP模块,用于读取所述信号存储模块中的回波信号并进行快速傅立叶变换FFT,求得信号的功率谱,根据所述功率谱得到粗估计频率值;对所述功率谱进行逆FFT得到所述回波信号的自相关函数,根据该自相关函数得到校正频率值,将所述粗估计频率值和校正频率值相加得到所述回波信号的多普勒频率,根据该多普勒频率计算出所述被测目标的运动速度。
2、根据权利要求1所述的多普勒测速激光雷达的回波信号处理装置,其特征在于,所述装置还包括:
滤波处理模块,用于接收放大处理后的回波信号,根据所述回波信号的频率值查询预先设定的滤波频率范围和采样率之间的对应关系,获取所述频率值对应的滤波频率范围,按照该滤波频率范围对所述回波信号进行带通滤波,将滤波处理后的回波信号传输给A/D转换模块;
A/D转换模块,用于根据所述回波信号的频率值查询预先设定的滤波频率范围和采样率之间的对应关系,获取所述频率值对应的采样率,按照该采样率对所述滤波处理模块输出的回波信号进行采样,将采样获取的模拟的回波信号转换为数字的回波信号,将该数字的回波信号传输给信号存储模块。
3、根据权利要求2所述的多普勒测速激光雷达的回波信号处理装置,其特征在于,所述装置还包括:
放大处理模块,用于接收多普勒测速激光雷达中的探测器输出的回波信号,将该回波信号放大到适合模数转换的幅度,将放大后的回波信号传输给滤波处理模块;
显示模块,用于将所述DSP模块所获取的被测目标的运动速度进行显示。
4、根据权利要求1至3任一项所述的多普勒测速激光雷达的回波信号处理装置,其特征在于,所述的滤波处理模块和A/D转换模块包括:
采样率和滤波频率范围分段设置模块,用于设定不同的滤波频率范围分别对应不同的采样率,并且各段滤波频率范围中的相邻两段之间互相重叠,将各个滤波频率范围和对应的采样率进行关联保存。
5、根据权利要求4所述的多普勒测速激光雷达的回波信号处理装置,其特征在于,所述的DSP模块具体包括:
信号读取模块,用于对所述信号存储模块中存储的回波信号进行读取;
粗估计频率值计算模块,用于对所述信号读取模块所读取的回波信号进行2倍采样点的FFT,求取FFT后的回波信号的功率谱,搜索该功率谱的最大谱峰对应的离散谱线,根据该离散谱线、采样率和采样点数计算出所述回波信号的粗估计频率值;
校正频率值计算模块,用于对所述功率谱进行逆FFT得到所述回波信号的自相关函数,利用所述粗估计频率值对该自相关函数进行移频处理,求取移频处理后信号的自相关函数在指定点的自相关值,根据该自相关值计算出所述回波信号的校正频率值;
运动速度获取模块,用于将所述信号频率粗估计值和校正频率值相加得到所述回波信号的多普勒频率,根据该多普勒频率和激光束与被测目标的运动方向之间的夹角计算出所述被测目标的运动速度;
当前滤波频率范围和采样率控制模块,用于判断运动速度计算模块获取的多普勒频率与当前滤波频率范围和采样率是否对应,如果是,则维持当前滤波频率范围和采样率不变;否则,根据所述运动速度计算模块获取的多普勒频率查询预先设定的滤波频率范围和采样率的对应关系,重新确定当前滤波频率范围和采样率。
6、一种多普勒测速激光雷达的回波信号处理方法,其特征在于,包括:
对被测目标返回的经过A/D转换后的回波信号进行FFT,根据FFT后的回波信号的功率谱获取所述回波信号的粗估计频率值,对所述功率谱进行逆FFT得到所述回波信号的自相关函数,根据该自相关函数得到所述回波信号的校正频率值;
将所述信号频率粗估计值和校正频率值相加得到所述回波信号的多普勒频率,根据该多普勒频率计算出所述被测目标的运动速度。
7、根据权利要求6所述的多普勒测速激光雷达的回波信号处理方法,其特征在于,所述方法之前还包括:
接收多普勒测速激光雷达中的探测器输出的回波信号,将该回波信号放大到适合模数转换的幅度,对放大后的回波信号按照默认的滤波频率范围和采样率分别进行滤波处理和采样,对采样的回波信号进行粗测得到回波信号的初始频率值;
根据所述初始频率值,查询预先设定的滤波频率范围和采样率之间的对应关系,获取所述初始频率值对应的滤波频率范围和采样率,按照该滤波频率范围对所述放大后的回波信号进行带通滤波,按照所述采样率对所述滤波处理后的回波信号进行采样,将采样获取的模拟回波信号转换为数字信号,将该数字信号进行存储。
8、根据权利要求7所述的多普勒测速激光雷达的回波信号处理方法,其特征在于,所述预先设定的各段滤波频率范围中的相邻两段之间互相重叠。
9、根据权利要求7或8所述的多普勒测速激光雷达的回波信号处理方法,其特征在于,所述的对被测目标返回的经过A/D转换后的回波信号进行FFT,根据FFT后的回波信号的功率谱获取所述回波信号的粗估计频率值,对所述功率谱进行逆FFT得到所述回波信号的自相关函数,根据该自相关函数得到所述回波信号的校正频率值的过程,具体包括:
对所述存储的数字的回波信号进行读取,对读取的回波信号进行2倍采样点的FFT,求取FFT后的回波信号的功率谱,搜索该功率谱的最大谱峰对应的离散谱线,根据该离散谱线、采样率和采样点数计算出所述回波信号的粗估计频率值;
对所述功率谱进行逆FFT得到所述回波信号的自相关函数,利用所述粗估计频率值对该自相关函数进行移频处理,求取移频后的自相关函数在指定点的自相关值,根据该自相关值计算出所述回波信号的校正频率值。
10、根据权利要求9所述的多普勒测速激光雷达的回波信号处理方法,其特征在于,所述的方法还包括:
判断所述多普勒频率与当前滤波频率范围、采样率是否对应,如果是,则维持当前滤波频率范围和采样率不变;否则,根据所述多普勒频率查询预先设定的滤波频率范围和采样率之间的对应关系,重新确定当前滤波频率范围和采样率。
CN2008101158144A 2008-06-27 2008-06-27 一种多普勒测速激光雷达的回波信号处理装置和方法 Expired - Fee Related CN101320086B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2008101158144A CN101320086B (zh) 2008-06-27 2008-06-27 一种多普勒测速激光雷达的回波信号处理装置和方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2008101158144A CN101320086B (zh) 2008-06-27 2008-06-27 一种多普勒测速激光雷达的回波信号处理装置和方法

Publications (2)

Publication Number Publication Date
CN101320086A true CN101320086A (zh) 2008-12-10
CN101320086B CN101320086B (zh) 2011-03-30

Family

ID=40180245

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2008101158144A Expired - Fee Related CN101320086B (zh) 2008-06-27 2008-06-27 一种多普勒测速激光雷达的回波信号处理装置和方法

Country Status (1)

Country Link
CN (1) CN101320086B (zh)

Cited By (27)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101806889A (zh) * 2010-02-26 2010-08-18 华中科技大学 一种优化调制激光雷达系统参数的装置和方法
CN102200579A (zh) * 2010-03-25 2011-09-28 付陆欣 手持式脉冲激光测距仪
CN102401896A (zh) * 2011-08-26 2012-04-04 湖南湘依铁路机车电器股份有限公司 一种用于计算多普勒信号频谱的方法
CN101839981B (zh) * 2010-01-28 2012-05-23 中国人民解放军电子工程学院 激光成像回波波形和层次特征获取方法及装置
CN101853240B (zh) * 2009-03-31 2012-07-04 华为技术有限公司 一种信号周期的估计方法和装置
CN102736074A (zh) * 2012-06-25 2012-10-17 中国科学院上海技术物理研究所 基于fpga的激光多普勒雷达信号处理器及处理方法
CN102890272A (zh) * 2012-11-05 2013-01-23 中国航天科工集团第二研究院二十三所 一种毫米波云雷达信号处理方法
CN103260780A (zh) * 2010-12-21 2013-08-21 西门子Vai金属科技有限责任公司 用于测量轧件的速度的方法和装置
CN103557753A (zh) * 2013-10-24 2014-02-05 北京理工大学 一种激光回波模拟装置及方法
CN104603636A (zh) * 2012-09-14 2015-05-06 三菱电机株式会社 激光雷达装置以及测定对象物的速度计算方法
CN105277947A (zh) * 2015-10-22 2016-01-27 中国人民解放军国防科学技术大学 一种基于fpga激光测速仪的信号处理系统及方法
CN105372649A (zh) * 2015-11-03 2016-03-02 黄中瑞 基于双路自适应采样滤波的铁路机车测速方法
CN106302292A (zh) * 2016-08-24 2017-01-04 深圳天珑无线科技有限公司 一种载波频率的估算方法和装置
CN107870330A (zh) * 2017-09-27 2018-04-03 西安空间无线电技术研究所 分段变重频模式采样均匀化数字处理方法及装置
CN108226919A (zh) * 2016-12-09 2018-06-29 通用汽车环球科技运作有限责任公司 用于使用goerzel滤波器进行计算上简单的距离多普勒角度跟踪的方法
CN108226918A (zh) * 2016-12-09 2018-06-29 通用汽车环球科技运作有限责任公司 高信噪比下多普勒模糊的解决方法
CN108614257A (zh) * 2011-04-14 2018-10-02 古野电气株式会社 水中探测装置、鱼群探测机以及水中探测方法
CN109283541A (zh) * 2017-07-21 2019-01-29 比亚迪股份有限公司 车载激光雷达测距装置、方法以及汽车
CN109597060A (zh) * 2018-12-07 2019-04-09 北京敏视达雷达有限公司 一种雷达测速方法及装置
CN110319916B (zh) * 2019-06-06 2020-06-30 浙江大学 一种基于水域反射波干扰测量的有限水池低频拓展方法
CN112630746A (zh) * 2020-12-01 2021-04-09 北京遥感设备研究所 一种用于远距目标测量的脉冲多普勒激光雷达
CN112673278A (zh) * 2020-09-10 2021-04-16 华为技术有限公司 基于回波信号的速度探测方法和装置
CN112752986A (zh) * 2020-09-28 2021-05-04 华为技术有限公司 雷达测距方法和装置
CN114324959A (zh) * 2021-12-03 2022-04-12 中国船舶重工集团公司七五0试验场 一种基于fpga和arm的复相关法测频方法及系统
CN114355360A (zh) * 2021-12-14 2022-04-15 中船航海科技有限责任公司 一种多普勒计程仪回波频率估计方法
CN115166667A (zh) * 2022-06-22 2022-10-11 哈尔滨工业大学(威海) 基于谐振区的运动目标识别方法、系统、设备及存储介质
WO2022242463A1 (en) * 2021-05-19 2022-11-24 Huawei Technologies Co.,Ltd. Method and apparatus for a computationally efficient lidar system

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5874914A (en) * 1995-10-09 1999-02-23 Snaptrack, Inc. GPS receiver utilizing a communication link
US6023236A (en) * 1998-11-30 2000-02-08 Kustom Signals, Inc. Speedometer assisted patrol speed search for DSP traffic radar

Cited By (41)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101853240B (zh) * 2009-03-31 2012-07-04 华为技术有限公司 一种信号周期的估计方法和装置
CN101839981B (zh) * 2010-01-28 2012-05-23 中国人民解放军电子工程学院 激光成像回波波形和层次特征获取方法及装置
CN101806889B (zh) * 2010-02-26 2012-07-04 华中科技大学 一种优化调制激光雷达系统参数的装置和方法
CN101806889A (zh) * 2010-02-26 2010-08-18 华中科技大学 一种优化调制激光雷达系统参数的装置和方法
CN102200579A (zh) * 2010-03-25 2011-09-28 付陆欣 手持式脉冲激光测距仪
CN103260780A (zh) * 2010-12-21 2013-08-21 西门子Vai金属科技有限责任公司 用于测量轧件的速度的方法和装置
US10228333B2 (en) 2010-12-21 2019-03-12 Primetals Technologies Austria GmbH Method and device for measuring the speed of a rolling stock
US10753886B2 (en) 2010-12-21 2020-08-25 Primetals Technologies Austria GmbH Method and device for measuring the speed of a rolling stock
CN108614257A (zh) * 2011-04-14 2018-10-02 古野电气株式会社 水中探测装置、鱼群探测机以及水中探测方法
CN102401896A (zh) * 2011-08-26 2012-04-04 湖南湘依铁路机车电器股份有限公司 一种用于计算多普勒信号频谱的方法
CN102736074A (zh) * 2012-06-25 2012-10-17 中国科学院上海技术物理研究所 基于fpga的激光多普勒雷达信号处理器及处理方法
CN104603636B (zh) * 2012-09-14 2016-10-26 三菱电机株式会社 激光雷达装置以及测定对象物的速度计算方法
CN104603636A (zh) * 2012-09-14 2015-05-06 三菱电机株式会社 激光雷达装置以及测定对象物的速度计算方法
CN102890272B (zh) * 2012-11-05 2013-12-04 中国航天科工集团第二研究院二十三所 一种毫米波云雷达信号处理方法
CN102890272A (zh) * 2012-11-05 2013-01-23 中国航天科工集团第二研究院二十三所 一种毫米波云雷达信号处理方法
CN103557753B (zh) * 2013-10-24 2015-05-27 北京理工大学 一种激光回波模拟装置及方法
CN103557753A (zh) * 2013-10-24 2014-02-05 北京理工大学 一种激光回波模拟装置及方法
CN105277947A (zh) * 2015-10-22 2016-01-27 中国人民解放军国防科学技术大学 一种基于fpga激光测速仪的信号处理系统及方法
CN105277947B (zh) * 2015-10-22 2017-12-19 中国人民解放军国防科学技术大学 一种基于fpga激光测速仪的信号处理系统及方法
CN105372649A (zh) * 2015-11-03 2016-03-02 黄中瑞 基于双路自适应采样滤波的铁路机车测速方法
CN106302292A (zh) * 2016-08-24 2017-01-04 深圳天珑无线科技有限公司 一种载波频率的估算方法和装置
CN108226918A (zh) * 2016-12-09 2018-06-29 通用汽车环球科技运作有限责任公司 高信噪比下多普勒模糊的解决方法
CN108226919A (zh) * 2016-12-09 2018-06-29 通用汽车环球科技运作有限责任公司 用于使用goerzel滤波器进行计算上简单的距离多普勒角度跟踪的方法
CN109283541B (zh) * 2017-07-21 2023-10-13 比亚迪股份有限公司 车载激光雷达测距装置、方法以及汽车
CN109283541A (zh) * 2017-07-21 2019-01-29 比亚迪股份有限公司 车载激光雷达测距装置、方法以及汽车
CN107870330B (zh) * 2017-09-27 2019-11-29 西安空间无线电技术研究所 分段变重频模式采样均匀化数字处理方法及装置
CN107870330A (zh) * 2017-09-27 2018-04-03 西安空间无线电技术研究所 分段变重频模式采样均匀化数字处理方法及装置
CN109597060A (zh) * 2018-12-07 2019-04-09 北京敏视达雷达有限公司 一种雷达测速方法及装置
CN109597060B (zh) * 2018-12-07 2021-06-25 北京敏视达雷达有限公司 一种雷达测速方法及装置
CN110319916B (zh) * 2019-06-06 2020-06-30 浙江大学 一种基于水域反射波干扰测量的有限水池低频拓展方法
CN112673278B (zh) * 2020-09-10 2022-01-11 华为技术有限公司 基于回波信号的速度探测方法和装置
CN112673278A (zh) * 2020-09-10 2021-04-16 华为技术有限公司 基于回波信号的速度探测方法和装置
CN112752986A (zh) * 2020-09-28 2021-05-04 华为技术有限公司 雷达测距方法和装置
CN112752986B (zh) * 2020-09-28 2022-04-08 华为技术有限公司 雷达测距方法和装置
CN112630746B (zh) * 2020-12-01 2023-09-15 北京遥感设备研究所 一种用于远距目标测量的脉冲多普勒激光雷达
CN112630746A (zh) * 2020-12-01 2021-04-09 北京遥感设备研究所 一种用于远距目标测量的脉冲多普勒激光雷达
WO2022242463A1 (en) * 2021-05-19 2022-11-24 Huawei Technologies Co.,Ltd. Method and apparatus for a computationally efficient lidar system
CN114324959A (zh) * 2021-12-03 2022-04-12 中国船舶重工集团公司七五0试验场 一种基于fpga和arm的复相关法测频方法及系统
CN114355360A (zh) * 2021-12-14 2022-04-15 中船航海科技有限责任公司 一种多普勒计程仪回波频率估计方法
CN114355360B (zh) * 2021-12-14 2024-05-03 中船航海科技有限责任公司 一种多普勒计程仪回波频率估计方法
CN115166667A (zh) * 2022-06-22 2022-10-11 哈尔滨工业大学(威海) 基于谐振区的运动目标识别方法、系统、设备及存储介质

Also Published As

Publication number Publication date
CN101320086B (zh) 2011-03-30

Similar Documents

Publication Publication Date Title
CN101320086B (zh) 一种多普勒测速激光雷达的回波信号处理装置和方法
CN109061589B (zh) 随机跳频雷达的目标运动参数估计方法
CN107607923B (zh) 基于lfmcw雷达的振动监测系统与信号处理方法
CN107907878B (zh) 高精度获取fmcw雷达距离测量值的方法
CN109964143B (zh) 用于处理由相干激光雷达引起的信号的方法及相关激光雷达系统
US4809703A (en) Ultrasonic doppler blood flow meter
CN101184443B (zh) 速度测量方法及采用该方法的速度测量装置
CN106644030A (zh) 一种基于多普勒雷达的非接触式振动测量方法
CN101595396A (zh) 高分辨率超声位移测量的装置与方法
JPH0268045A (ja) パルスドップラ信号のフーリエプロセッシング方法および装置
Pang et al. Coherent detection algorithm for radar maneuvering targets based on discrete polynomial-phase transform
CN104502898B (zh) 将修正rft和修正mdcft相结合的机动目标参数估计方法
CN101788671B (zh) 应用于外差探测啁啾调幅激光测距装置的多周期调制方法
CN110596671A (zh) 一种lfmcw测速测距雷达的优化处理方法及系统
Yan et al. Through-the-wall human respiration detection using impulse ultra-wide-band radar
Parrish An overview of FMCW systems in MATLAB
RU2018109177A (ru) Способ оценки пространственного размера воздушной цели по частотной протяженности доплеровского портрета
CN109581350B (zh) 基于时频积分插值的雷达测距测速方法及装置
CN111750972A (zh) 一种适用于双频连续激光雷达的目标振动测量方法
JP2013113723A (ja) レーダ装置
JP4937782B2 (ja) レーダ装置
Ciattaglia et al. Assessment of Displacement Measurements by a mmWave Radar
RU2800494C1 (ru) Устройство цифровой обработки сигналов в импульсно-доплеровской рлс с компенсацией миграций целей по дальности и чм доплеровских сигналов за один период излучения и приема пачки радиоимпульсов
RU2807316C1 (ru) Способ определения параметров движения высокоскоростного воздушного объекта
CN114609606B (zh) 目标对象的参数估计方法、装置和电子设备

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C53 Correction of patent for invention or patent application
CB03 Change of inventor or designer information

Inventor after: Ou Pan

Inventor after: Lin Zhili

Inventor after: Jia Yudong

Inventor after: Sun Mingjie

Inventor after: Wang Jiqiang

Inventor after: Zhang Chunxi

Inventor before: Feng Di

Inventor before: Wang Jiqiang

Inventor before: Ou Pan

Inventor before: Zhang Chunxi

Inventor before: Jia Yudong

Inventor before: Tian Lun

COR Change of bibliographic data

Free format text: CORRECT: INVENTOR; FROM: FENG DI WANG JIQIANG OU PAN ZHANG CHUNXI JIA YUDONG TIAN LUN TO: OU PAN LIN ZHILI JIA YUDONG SUN MINGJIE WANG JIQIANG ZHANG CHUNXI

C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20110330

Termination date: 20110627