CN111337902B - 多通道高重频大动态范围测距测速激光雷达方法及装置 - Google Patents

多通道高重频大动态范围测距测速激光雷达方法及装置 Download PDF

Info

Publication number
CN111337902B
CN111337902B CN202010354911.XA CN202010354911A CN111337902B CN 111337902 B CN111337902 B CN 111337902B CN 202010354911 A CN202010354911 A CN 202010354911A CN 111337902 B CN111337902 B CN 111337902B
Authority
CN
China
Prior art keywords
frequency
target
signal
phase
light beam
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
Application number
CN202010354911.XA
Other languages
English (en)
Other versions
CN111337902A (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.)
Sun Jianfeng
Original Assignee
Hangzhou Ilidar Technology Co ltd
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 Hangzhou Ilidar Technology Co ltd filed Critical Hangzhou Ilidar Technology Co ltd
Priority to CN202010354911.XA priority Critical patent/CN111337902B/zh
Publication of CN111337902A publication Critical patent/CN111337902A/zh
Application granted granted Critical
Publication of CN111337902B publication Critical patent/CN111337902B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/48Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
    • G01S7/481Constructional features, e.g. arrangements of optical elements
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S17/00Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/02Systems using the reflection of electromagnetic waves other than radio waves
    • G01S17/50Systems of measurement based on relative movement of target
    • G01S17/58Velocity or trajectory determination systems; Sense-of-movement determination systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Electromagnetism (AREA)
  • Optical Radar Systems And Details Thereof (AREA)

Abstract

本发明公开了一种多通道高重频大动态范围测距测速激光雷达方法及装置,激光光源产生光束通过分束器分为N个通道的输出光束,每通道输出光束进行宽带线性调频,然后经过不同幅度的移相后驱动相位调制器,使每通道输出光束具有不同的时间延迟;再经分束器分为本振光束和发射光束;将N个发射光束并行发射至目标并接收目标的回波光束,再与对应的本振光束进行相干光混频,获得包含目标距离和速度信息的中频信号,对中频信号的采样数据进行并行快速傅里叶变换和互谱处理,实现目标距离和速度的并行同步测量,最后由主控计算机实现N个通道点云图像的合并输出。本发明能够有效克服测距模糊,具有高重频、大动态测距范围、高分辨率、高精度、可靠性高、体积小和质量轻等优点。

Description

多通道高重频大动态范围测距测速激光雷达方法及装置
技术领域
本发明涉及激光雷达技术领域,特别涉及一种多通道高重频大动态范围测距测速激光雷达方法及装置。
背景技术
调频连续波激光雷达将现代雷达技术中的调频连续波测距与激光探测技术相结合。该技术在时间上线性改变发射信号的频率,并测量回波信号与本振信号的差拍信号瞬时频率的方法来获得目标距离信息和径向速度信息。相较于脉冲激光雷达,调频连续波激光雷达具有同步测距测速、测距/测速范围大、分辨率高、多目标探测及目标成像等优势,在自动驾驶、高精度三维成像和遥感测绘等领域得到了广泛应用。
但是目前啁啾调频连续波激光雷达依然存在以下问题:为了取得较大调谐范围从而获得较高距离分辨率,脉冲重复频率(PRF)严重受限;大范围扫频带来的频率调制非线性依然是一个未解决的难题,严重影响测速测距的精度;线性调频去斜解调要求回波信号与本振信号在时域上有重叠,决定了其不适用于较大测距动态范围内的探测。
如图1所示的对称三角线性调频连续波相干激光雷达测距测速信号示意图。注意这里考虑多普勒频移远低于扫频带宽,图中仅仅是为了显示方便。这里定义进行有效频谱分析的时间Teff,时域信号有效重叠比例η:
Figure GDA0003311751660000021
当回波光束与本振光束的相对时间差τSL较小时,回波光束可以和本振光束在线性调频阶段可以做到较好的重叠,在探测器上较长一段时间内回波信号和本振信号存在固定频差,频谱分析时,可得到较高的信噪比及谱分辨率。有效重叠比例η(τ)概念可以用于分析距离模糊:
当τSL≈0时,回波信号和本振信号在相同的变频阶段几乎完全重叠,两者在变频阶段的频率差一直是一个固定的值,有效重叠比例为η(0)=100%;
Figure GDA0003311751660000022
时,有效重叠比例很低,这时获得的信号信噪比很低。
Figure GDA0003311751660000023
时,即回波信号和本振信号在同一时刻变频方向完全相反,在相同的变频阶段完全没有重叠,则两者的频率差不固定,将无法分析出任何有效的频率,有效重叠比例为η(±T/2)=0%。
Figure GDA0003311751660000024
时,回波信号和本振信号不能和同一个对称三角调频信号在同一个变频方向的线性调频部分重叠,无法测量,这里称为测距模糊。
延迟造成的回波-本振频率差和有效重叠比例均具有周期性。对称三角FMCW激光雷达***探测的飞行距离动态范围小于二分之一个周期所允许的距离。
吴军等提出了双频双调制双本振对称三角线性调频相干测距测速激光雷达体制,扩展了对称三角线性调频测距测速体制的距离作用范围,且有高的探测重复频率(吴军,洪光烈,何志平,舒嵘,一种大测距动态范围高重频相干测距测速激光雷达(Ⅰ):体制及性能,红外与毫米波学报,2014,Vol.33,No.6,680-690);理论上,通过延长信号的调制周期可以增加***的测距动态范围,但该方式一方面探测时需要大量采集相干拍频的数据点作频谱分析,另一方面也要求对称三角线性调频的调制时间带宽积增大,而长时间高线性度的宽带线性调频难以实现,大大增加了***复杂性,影响了应用前景。
发明内容
本发明的目的在于,提供一种多通道高重频大动态范围测距测速激光雷达方法及装置。本发明能够有效克服测距模糊,具有高重频、大动态测距范围、高分辨率、高精度、体积小和质量轻等优点。
本发明的技术方案:多通道高重频大动态范围测距测速激光雷达方法,雷达平台中激光光源产生光束经过放大后通过1×N分束器分为N个通道的输出光束,每通道输出光束进入对应的相位调制器和光学滤波器实现宽带线性调频,将射频信号源经过不同幅度的移相后驱动相位调制器,使每通道输出光束具有不同的时间延迟;将每通道输出光束再通过移频器进行移频后,经1×2分束器分为本振光束和发射光束;将N个通道的发射光束通过各自通道的光学望远镜和光束定向器并行发射至目标,并实现N个通道发射/接收视场的完全匹配;在各个通道的光学望远镜的接收端接收对应目标的回波光束,并与对应的本振光束通过光学桥接器进行相干光混频,采用平衡接收获得包含目标距离和速度信息的中频信号,对中频信号进行滤波处理和采样处理得到采样数据,使用现场可编程门阵列对获得的N个通道的采样数据进行实时的并行快速傅里叶变换和互谱处理,从而实现目标距离和速度的并行同步测量,最终由主控计算机实现N个通道点云图像的合并输出。
上述的多通道高重频大动态范围测距测速激光雷达方法,所述的实现目标距离和速度的同步测量,具体是中频信号中同相信号和正交信号经过滤波和采样处理后,分别进行傅里叶变换,再进行互谱处理取其虚部,再利用重心法提取频谱中峰值的位置和正负,得到雷达平台和目标相对运动产生的多普勒频移,然后由多普勒频移得到雷达平台与目标相对运动径向速度的大小和方向以及目标距离。
前述的多通道高重频大动态范围测距测速激光雷达方法,所述雷达平台中激光光源为窄线宽连续激光光源,经起偏器起偏后由激光放大器放大,再通过1×N分束器分为N个通道的输出光束,每通道输出光束分别经过相位调制器和光学滤波器实现宽带线性调频,产生的输出光束为频率线性调制的连续相干激光,采用对称三角波线性调制,调制信号的频率随时间成对称三角形变换,在一个周期内,前半部分为正向调频,后半部分为负向调频;其中,将射频信号源经过移相后驱动相位调制器,实现每个光束的不同时间延迟,然后每通道输出光束经过移频器移频,经过延时和移频后的第n路输出光束的光场表示为:
Figure GDA0003311751660000051
其中,t是时间,E0是振幅,T为调频周期,f0为调频初始频率,fshift_n是第n路输出光束的移频量,
Figure GDA0003311751660000052
为调频速率,B为调频带宽,
Figure GDA0003311751660000053
是第n路输出光束的移相,
Figure GDA0003311751660000054
N为正整数,Tshift_n∈[0,T),
Figure GDA0003311751660000055
φup(n)为第n个输出光束脉冲上升段的初始位相,φdown(n)为第n个输出光束脉冲下降段的初始位相,exp是以自然常数e为底的指数函数,
Figure GDA0003311751660000056
前述的多通道高重频大动态范围测距测速激光雷达方法,所述第n路输出光束经过延时和移频后,再经过1×2分束器分束,小部分能量作为本振光束,本振光束为时间延迟τL的线性调频信号,光场表示为:
Figure GDA0003311751660000057
其中,EL是本振光束振幅,φLO是本振光束的噪声位相;
大部分能量作为发射光束,经过空间光学环形器,再通过光学望远镜和光束定向器发射至目标,并由光学望远镜接收目标的回波光束,回波光束为时间延迟τS的线性调频信号,表示为:
Figure GDA0003311751660000061
其中,ES是回波光束振幅,φS是回波光束的噪声位相;
经过相同移频的目标回波光束和本振光束合束后光场表示为:
Figure GDA0003311751660000062
回波光束的时间延迟τS与本振光束的时间延迟τL的关系表示为:
Figure GDA0003311751660000063
其中,c是光速,R是目标距离,V是雷达平台与目标相对运动径向速度,fDoppler是雷达平台与目标相对运动径向速度引起的多普勒频移,
所述的回波光束和本振光束经过2×490°光学桥接器混频后的四路输出分别为:
Figure GDA0003311751660000071
Figure GDA0003311751660000072
其中,φN-n是第n路输出光束混频噪声位相,Is是和回波光束有关的直流量;Io是和本振光束有关的直流量;
光学桥接器输出的具有正交特性的同相信号和正交信号分别由光电平衡探测器进行接收,获得包含目标距离和速度信息的中频信号;所述的正向调频过程的中频信号中的同相信号和正交信号分别为:
Figure GDA0003311751660000073
所述的负向调频过程的中频信号中的同相信号和正交信号分别为:
Figure GDA0003311751660000081
其中kin是接收同相信号的光电平衡探测器响应率,kqu是接收正交信号的光电平衡探测器响应率,φi-n和φq-n分别是同相信号和正交信号的噪声位相;
将同相和正交通道的振幅用下式来代替:
Figure GDA0003311751660000082
正向调频过程的中频信号中的同相信号和正交信号简化为:
Figure GDA0003311751660000083
负向调频过程的中频信号中的同相信号和正交信号简化为:
Figure GDA0003311751660000084
同相信号和正交信号分别经过低通滤波器滤波,由模数转换器完成模数转换,再由现场可编程门阵列采集进行并行快速傅里叶变换,同相信号傅里叶变换表示为:
Figure GDA0003311751660000085
正交信号傅里叶变换表示为:
Figure GDA0003311751660000091
两通道进行互谱处理:
Figure GDA0003311751660000092
最后仅取虚部得到
Img=δ2(f-fn)-δ2(f+fn),
通过重心法提取频谱峰值位置和正负,就可以分别得到正向调频和负向调频过程中的中频频率值:
Figure GDA0003311751660000093
由上式可以得到:
Figure GDA0003311751660000094
上式中,fn-up是正向调频过程中的中频频率值,fn-down是负向调频过程中的中频频率值;
由于多普勒频率大小与雷达平台与目标相对运动的速度成正比,多普勒频移正负与相对运动径向速度的方向有关,正频移代表雷达平台与目标相向运动,负频移代表雷达平台与目标相背运动;因此由多普勒频移可以得到雷达平台与目标相对运动径向速度的大小和方向,因此由多普勒频移可以得到雷达平台与目标相对运动径向速度的大小和方向,表示为
Figure GDA0003311751660000101
式中,λ是输出光束波长,fDoppler是雷达平台与目标相对运动径向速度引起的多普勒频移;
由上式得到目标点的距离:
Figure GDA0003311751660000102
式中,
Figure GDA0003311751660000103
为调频速率,B为调频带宽,T为调频周期;
主控计算机将采集到的N个通道点云图像合并,获取一个调频周期T时间内目标坐标的空间点集:
Figure GDA0003311751660000104
P(xP_n,yP_n,zP_n)为目标P在第n个通道探测到的空间坐标系中的坐标位置;
该空间点集∑P即为最终显示的3D点云图像。
前述的多通道高重频大动态范围测距测速激光雷达方法,所述同相信号和正交信号分别经过低通滤波器滤波去除串扰信号具体是,所述串扰信号表示为:
Figure GDA0003311751660000111
,ES_m表示为串扰的回波光束振幅,ELO_n表示为第n个回波光束对应的本振光束振幅;
串扰信号数据中包含频率为|fshift_n-fshift_m|项,通过低通滤波器对其进行滤波处理,消除高频串扰信号,提高单路输出光束的探测精度。
实现如前述的多通道高重频大动态范围测距测速激光雷达方法的装置,包括激光光源,所述激光光源的输出端依次连接有起偏器和激光放大器,所述激光放大器的输出端经1×N分束器连接有N个相位调制器;每个所述的相位调制器均连接有移相器,N个移相器一同连接有射频信号源,所述相位调制器经光学滤波器连接有移频器,移频器的输出端连接有1×2分束器,所述1×2分束器连接有光学环形器;
所述光学环形器的输出端依次连接有光学望远镜和光束定向器;所述光学环形器和1×2分束器一同连有光学桥接器;所述光学桥接器的输出端依次连接有光电平衡探测器和低通滤波器;所述低通滤波器经模数转换器连接有现场可编程门阵列,所述现场可编程门阵列的输出端连接有主控计算机;所述主控计算机还与光束定向器连接。
与现有技术相比,本发明通过将激光光源产生光束经过放大后通过1×N分束器分为N个通道的输出光束,每通道输出光束经过延时和移频后,再经1×2分束器分为本振光束和发射光束;然后将N个通道的发射光束通过光学望远镜和光束定向器并行发射至目标,并实现N个通道发射/接收视场的完全匹配;在接收端接收目标的回波光束,并与对应的本振光束通过光学桥接器进行相干光混频,采用平衡接收获得包含目标距离和速度信息的中频信号,对中频信号进行滤波处理和采样处理得到采样数据,使用现场可编程门阵列对获得的N个通道采样数据进行实时的并行快速傅里叶变换和互谱处理,从而实现目标距离和速度的并行同步测量。由此本发明可以获得同一目标N个包含目标距离和速度信息的中频信号且并行输出,再使用现场可编程门阵列对获得的N个中频信号进行实时的并行快速傅里叶变换(FFT)和互谱处理,就可以实现对目标距离和速度的并行同步测量。目前的激光雷达存在高重频的问题,即脉冲式雷达、连续波式雷达通过发出脉冲序列,在目标距离一定的情况下,如果序列之间的时间间隔过短或频率太高,在接收端就难以分辨脉冲序列对应的回波,从而造成测距模糊,目前原则上只有伪随机码可以解决高重频的问题,但是采用伪随机码测距对硬件要求很高,而且需要补偿多普勒频移引入的频谱展宽问题,可实现难度大。因此本发明采用多个通道进行输出光束的发射,即按顺序将第1路输出光束、第2路输出光束直至最后1路输出光束分别延时和移频后发射,然后重复此过程,因此可以在一定的重复频率下,测距动态范围可以提高N倍,有效克服了测距动态范围模糊;或者在一定的测距动态范围下,重复频率可以提高N倍,因此具有高重频、大动态测距范围优点;本发明由主控计算机实现N个通道点云图像的合并输出,实现了高分辨3D成像;本发明采用N个通道输出光束的并行发射/接收,每通道输出光束发射/接收同轴,有利于集成小型化,降低了***的复杂性,实现每通道输出光束发射/接收视场的完全匹配,而且本发明采用分布式结构,减小了整体装置的体积,即使个别通道发生损坏也不会造成***瘫痪,从而具有体积小、质量轻、可靠性高、高分辨率、高精度等优点。本发明在车载、机载和星载等领域具有良好的发展前景。
附图说明
图1是对称三角线性调频连续波相干激光雷达回波光束-本振光束的波形关系和测距测速中频信号示意图。
图2是本发明的结构示意图。
图3是对称三角线性频率调制波形及频率调制延时的示意图;
附图中的标记为:1、激光光源;2、起偏器;3、激光放大器;4、1×N分束器;5、相位调制器;6、移频器;7、1×2分束器;8、光学环形器;9、光学望远镜;10、光束定向器;11、光学桥接器;12、光电平衡探测器;13、低通滤波器;14、模数转换器;15、现场可编程门阵列;16、主控计算机;17、移相器;18、射频信号源;19、光学滤波器。
具体实施方式
下面结合附图和实施例对本发明作进一步的说明,但并不作为对本发明限制的依据。
实施例1:雷达平台中窄线宽连续激光光源产生光束经过放大后通过1×N分束器分为N个通道的输出光束,每通道输出光束分别经过相位调制器和光学滤波器实现宽带线性调频,用于后续的测速和测距,其中,将射频信号源经过移相后驱动相位调制器,实现每通道输出光束的不同时间延迟,既在一个周期时间T内,第一路输出光束不延时发射,第二路输出光束经过
Figure GDA0003311751660000141
的延时后进行发射,第三路输出光束经过
Figure GDA0003311751660000142
的延时后进行发射,……,第n路输出光束经过
Figure GDA0003311751660000143
的延时后进行发射,直至最后一路输出光束经过
Figure GDA0003311751660000144
的延时后进行发射;将N个经过延时的输出光束分别再经过移频器进行移频,然后经1×2分束器分为本振光束和发射光束;将N个通道的发射光束通过光学望远镜和光束定向器并行发射至目标,并实现N个通道发射/接收视场的完全匹配;在各个通道的光学望远镜的接收端接收对应目标的回波光束,回波光束经过光学环形器后并与对应的本振光束通过光学桥接器进行相干光混频,采用平衡接收获得包含目标距离和速度信息的中频信号,对中频信号进行滤波处理和采样处理得到采样数据,使用现场可编程门阵列对获得的N个通道采样数据进行实时的并行快速傅里叶变换和互谱处理,从而实现目标距离和速度的并行同步测量;最终,由主控计算机实现N个通道点云图像的合并输出。实现上述方法的装置,如图2所所示,包括激光光源1,所述激光光源1的输出端依次连接有起偏器2和激光放大器3,所述激光放大器3的输出端经1×N分束器4连接有N个相位调制器5;每个所述的相位调制器5均连接有移相器17,N个移相器17一同连接有射频信号源18,所述相位调制器5经光学滤波器19连接有移频器6,移频器6的输出端连接有1×2分束器7,所述1×2分束器7连接有光学环形器8;
所述光学环形器8的输出端依次连接有光学望远镜9和光束定向器10;所述光学环形器8和1×2分束器7一同连有光学桥接器11;所述光学桥接器11的输出端依次连接有光电平衡探测器12和低通滤波器13;所述低通滤波器13经模数转换器14连接有现场可编程门阵列15,所述现场可编程门阵列15的输出端连接有主控计算机16;所述主控计算机16还与光束定向器10连接,主控计算机用于控制光学定向器10的指向和实现N个通道点云图像的合并输出,以及其他传感器数据采集和决策任务。
实施例2:多通道高重频大动态范围测距测速激光雷达方法,如图1所示,包括激光光源1,激光光源1采用人眼安全的1550nm单模窄线宽连续光纤激光器,激光器线宽10kHz,输出功率为20mW,光纤输出有隔离保护,经过起偏器2起偏,确保偏振消光比大于25dB,经过掺铒光纤放大器3放大至600mW,然后偏振光通过1×2分束器4分成2路输出光束,每个光束先后经过光纤相位调制器和光学滤波器,采用调频信号发生器产生的调频信号和基频信号发生器产生的基频信号的混频信号作为光纤相位调制器的驱动信号,驱动光纤相位调制器产生调频激光信号,并通过光学滤波器抑制谐波保留所需阶次的调频激光信号,输出光束采用对称三角波线性调制,调制信号的频率随时间成对称三角形变换,在一个周期内,前半部分为正向调频,后半部分为负向调频。调频带宽2.5GHz,调频速率0.5PHz/s,调频周期10μs,重复频率100kHz;本实施例中以两路输出光束为例,第1路输出光束不经过延时,第2路输出光束的射频信号源经过π移相后驱动相位调制器,实现光束时间延迟,第1路输出光束移频量为0MHz,第2路输出光束经过移频器6移频,移频量为150MHz,如图2所示;
第2路输出光束的光场表示为:
Figure GDA0003311751660000161
其中,n=2,t是时间,E0是振幅,T为调频周期,f0为调频初始频率,fshift_n是第n路输出光束的移频量,
Figure GDA0003311751660000162
为调频速率,B为调频带宽,
Figure GDA0003311751660000163
是第n路输出光束的移相,
Figure GDA0003311751660000164
N为正整数,Tshift_n∈[0,T),
Figure GDA0003311751660000165
φup(n)为第n个发射激光脉冲上升段的初始位相,φdown(n)为第n个发射激光脉冲下降段的初始位相,exp是以自然常数e为底的指数函数,
Figure GDA0003311751660000166
然后经过1×2分束器7出射出本振光束和发射光束,考虑到插损,每道光束的功率约250mW,输出两路光强度为1:99,小部分能量作为本振光速;
如图3所示,本振光束为时间延迟τL的线性调频信号,光场表示为:
Figure GDA0003311751660000171
其中,EL是本振光束振幅,φLO是本振光束的噪声位相;大部分能量作为发射光束,经过光学环形器8(光纤三端口环形器),再通过发射/接收光学望远镜9和光束定向器8发射至目标,实并由光学望远镜接收目标的回波光束,实现N个通道发射/接收视场的完全匹配;
回波光束为时间延迟τS的线性调频信号;表示为:
Figure GDA0003311751660000172
其中,ES是回波光束振幅,φS是回波光束的噪声位相;
经过相同移频的目标回波光束和本振光束合束后光场表示为:
Figure GDA0003311751660000181
回波光束的时间延迟τS与本振光束的时间延迟τL的关系表示为:
Figure GDA0003311751660000182
其中,c是光速,R是目标距离,V是雷达平台与目标相对运动径向速度,fDoppler是雷达平台与目标相对运动径向速度引起的多普勒频移,
所述的回波光束和本振光束经过2×490°光学桥接器11混频后的四路输出分别为:
Figure GDA0003311751660000183
Figure GDA0003311751660000184
其中,φN-n是第n路输出光束混频噪声位相,Is是和回波光束有关的直流量;Io是和本振光束有关的直流量;
光学桥接器11输出的具有正交特性的同相信号和正交信号分别由光电平衡探测器12进行接收,获得包含目标距离和速度信息的中频信号;所述的正向调频过程的中频信号中的同相信号和正交信号分别为:
Figure GDA0003311751660000192
所述的负向调频过程的中频信号中的同相信号和正交信号分别为::
Figure GDA0003311751660000193
其中kin是接收同相信号的光电平衡探测器响应率,kqu是接收正交信号的光电平衡探测器响应率,φi-n和φq-n分别是同相信号和正交信号的噪声位相;
将同相和正交通道的振幅用下式来代替:
Figure GDA0003311751660000194
正向调频过程的中频信号中的同相信号和正交信号简化为:
Figure GDA0003311751660000201
负向调频过程的中频信号中的同相信号和正交信号简化为:
Figure GDA0003311751660000202
对于其他输出光束的串扰信号表示为:
Figure GDA0003311751660000203
ES_m表示为串扰的回波光束振幅,ELO_n表示为第n个回波光束对应的本振光束振幅;
串扰信号数据中包含频率为|fshift_n-fshift_m|项,通过低通滤波器13对其进行滤波处理,消除高频串扰信号,提高单路输出光束的探测精度;
所述同相信号和正交信号分别经过低通滤波器滤波13后,由模数转换器14完成模数转换,再由现场可编程门阵列15采集进行并行快速傅里叶变换,同相信号傅里叶变换表示为:
Figure GDA0003311751660000211
正交信号傅里叶变换表示为:
Figure GDA0003311751660000212
两通道进行互谱处理:
Figure GDA0003311751660000213
最后仅取虚部得到
Img=δ2(f-fn)-δ2(f+fn),
在主控计算机16中,通过重心法提取频谱峰值位置和正负,就可以分别得到正向调频和负向调频过程中的频率值:
Figure GDA0003311751660000214
由上式可以得到:
Figure GDA0003311751660000215
上式中,fn-up是正向调频过程中的中频频率值,fn-down是负向调频过程中的中频频率值;
由于多普勒频率大小与雷达平台与目标相对运动的速度成正比,多普勒频移正负与相对运动径向速度的方向有关,正频移代表雷达平台与目标相向运动,负频移代表雷达平台与目标相背运动。因此由多普勒频移可以得到雷达平台与目标相对运动径向速度的大小和方向,因此由多普勒频移可以得到雷达平台与目标相对运动径向速度的大小和方向,表示为
Figure GDA0003311751660000221
式中,λ是输出光束波长,fDoppler是雷达平台与目标相对运动径向速度引起的多普勒频移;
由上式得到目标点的距离:
Figure GDA0003311751660000222
式中,
Figure GDA0003311751660000223
为调频速率,B为调频带宽,T为调频周期。
主控计算机将采集到的N个通道点云图像合并,获取一个调频周期T时间内目标坐标的空间点集:
Figure GDA0003311751660000224
P(xP_n,yP_n,zP_n)为目标P在第n个通道探测到的空间坐标系中的坐标位置;
该空间点集∑P即为最终显示的3D点云图像。
本实施例中的测距动态范围为1500m,重复频率200kHz,比常规方案提高1倍。
综上所述,本发明可以获得同一目标N个包含目标距离和速度信息的中频信号且并行输出,每通道输出光束具有不同的时间延迟,再使用现场可编程门阵列对获得的N个中频信号进行实时的并行快速傅里叶变换(FFT)和互谱处理,就可以实现对目标距离和速度的并行同步测量,在一定的重复频率下,测距动态范围可以提高N倍;或者在一定的测距动态范围下,重复频率可以提高N倍,因此具有高重频、大动态测距范围的优点;主控计算机实现N个通道点云图像的合并输出,实现了高分辨3D成像;本发明采用N个输出光束的并行发射/接收,每通道输出光束发射/接收同轴,采用分布式结构,实现N通道输出光束发射/接收视场的完全匹配,即使个别通道发生损坏也不会造成***瘫痪,从而具有体积小、质量轻、可靠性高、高分辨率、高精度等优点。

Claims (5)

1.多通道高重频大动态范围测距测速激光雷达方法,其特征在于:雷达平台中激光光源产生光束经过放大后通过1×N分束器分为N个通道的输出光束,每通道输出光束进入对应的相位调制器和光学滤波器实现宽带线性调频,将射频信号源经过不同幅度的移相后驱动相位调制器,使每通道输出光束具有不同的时间延迟;将每通道输出光束再通过移频器进行移频后,经1×2分束器分为本振光束和发射光束;将N个通道的发射光束通过各自通道的光学望远镜和光束定向器并行发射至目标,并实现N个通道发射/接收视场的完全匹配;在各个通道的光学望远镜的接收端接收对应目标的回波光束,并与对应的本振光束通过光学桥接器进行相干光混频,采用平衡接收获得包含目标距离和速度信息的中频信号,对中频信号进行滤波处理和采样处理得到采样数据,使用现场可编程门阵列对获得的N个通道的采样数据进行实时的并行快速傅里叶变换和互谱处理,从而实现目标距离和速度的并行同步测量,最终由主控计算机实现N个通道点云图像的合并输出;
所述的实现目标距离和速度的同步测量,具体是中频信号中同相信号和正交信号经过滤波和采样处理后,分别进行傅里叶变换,再进行互谱处理取其虚部,再利用重心法提取频谱中峰值的位置和正负,得到雷达平台和目标相对运动产生的多普勒频移,然后由多普勒频移得到雷达平台与目标相对运动径向速度的大小和方向以及目标距离。
2.根据权利要求1所述的多通道高重频大动态范围测距测速激光雷达方法,其特征在于:所述雷达平台中激光光源为窄线宽连续激光光源,经起偏器起偏后由激光放大器放大,再通过1×N分束器分为N个通道的输出光束,每通道输出光束分别经过相位调制器和光学滤波器实现宽带线性调频,产生的输出光束为频率线性调制的连续相干激光,采用对称三角波线性调制,调制信号的频率随时间成对称三角形变换,在一个周期内,前半部分为正向调频,后半部分为负向调频;其中,将射频信号源经过移相后驱动相位调制器,实现每个光束的不同时间延迟,然后每个通道的输出光束经过移频器移频,经过延时和移频后的第n路输出光束的光场表示为:
Figure FDA0003311751650000021
其中,t是时间,E0是振幅,T为调频周期,f0为调频初始频率,fshift_n是第n路输出光束的移频量,
Figure FDA0003311751650000022
为调频速率,B为调频带宽,
Figure FDA0003311751650000023
是第n路输出光束的移相,
Figure FDA0003311751650000024
K为正整数,Tshift_n∈[0,T),
Figure FDA0003311751650000025
φup(n)为第n个输出光束脉冲上升段的初始位相,φdown(n)为第n个输出光束脉冲下降段的初始位相,exp是以自然常数e为底的指数函数,
Figure FDA0003311751650000026
3.根据权利要求2所述的多通道高重频大动态范围测距测速激光雷达方法,其特征在于:所述第n路输出光束经过延时和移频后,再经过1×2分束器分束,小部分能量作为本振光束,本振光束为时间延迟τL的线性调频信号,光场表示为:
Figure FDA0003311751650000031
其中,EL是本振光束振幅,φLO是本振光束的噪声位相;
大部分能量作为发射光束,经过空间光学环形器,再通过光学望远镜和光束定向器发射至目标,并由光学望远镜接收目标的回波光束,回波光束为时间延迟τS的线性调频信号,表示为:
Figure FDA0003311751650000032
其中,ES是回波光束振幅,φS是回波光束的噪声位相;
经过相同移频的目标回波光束和本振光束合束后光场表示为:
Figure FDA0003311751650000041
回波光束的时间延迟τS与本振光束的时间延迟τL的关系表示为:
Figure FDA0003311751650000042
其中,c是光速,R是目标距离,V是雷达平台与目标相对运动径向速度,fDoppler是雷达平台与目标相对运动径向速度引起的多普勒频移,
所述的回波光束和本振光束经过2×4 90°光学桥接器混频后的四路输出分别为:
Figure FDA0003311751650000043
Figure FDA0003311751650000051
其中,φN-n是第n路输出光束混频噪声位相,Is是和回波光束有关的直流量;Io是和本振光束有关的直流量;
光学桥接器输出的具有正交特性的同相信号和正交信号分别由光电平衡探测器进行接收,获得包含目标距离和速度信息的中频信号;所述的正向调频过程的中频信号中的同相信号和正交信号分别为:
Figure FDA0003311751650000052
所述的负向调频过程的中频信号中的同相信号和正交信号分别为:
Figure FDA0003311751650000053
其中kin是接收同相信号的光电平衡探测器响应率,kqu是接收正交信号的光电平衡探测器响应率,φi-n和φq-n分别是同相信号和正交信号的噪声位相;
将同相和正交通道的振幅用下式来代替:
Figure FDA0003311751650000061
正向调频过程的中频信号中的同相信号和正交信号简化为:
Figure FDA0003311751650000062
负向调频过程的中频信号中的同相信号和正交信号简化为:
Figure FDA0003311751650000063
同相信号和正交信号分别经过低通滤波器滤波,由模数转换器完成模数转换,再由现场可编程门阵列采集进行并行快速傅里叶变换,同相信号傅里叶变换表示为:
Figure FDA0003311751650000064
正交信号傅里叶变换表示为:
Figure FDA0003311751650000065
两通道进行互谱处理:
Figure FDA0003311751650000071
最后仅取虚部得到
Img=δ2(f-fn)-δ2(f+fn),
通过重心法提取频谱峰值位置和正负,就可以分别得到正向调频和负向调频过程中的中频频率值:
Figure FDA0003311751650000072
由上式可以得到:
Figure FDA0003311751650000073
上式中,fn-up是正向调频过程中的中频频率值,fn-down是负向调频过程中的中频频率值;
由于多普勒频率大小与雷达平台与目标相对运动的速度成正比,多普勒频移正负与相对运动径向速度的方向有关,正频移代表雷达平台与目标相向运动,负频移代表雷达平台与目标相背运动;因此由多普勒频移可以得到雷达平台与目标相对运动径向速度的大小和方向,因此由多普勒频移可以得到雷达平台与目标相对运动径向速度的大小和方向,表示为
Figure FDA0003311751650000081
式中,λ是输出光束波长,fDoppler是雷达平台与目标相对运动径向速度引起的多普勒频移;
由上式得到目标点的距离:
Figure FDA0003311751650000082
式中,
Figure FDA0003311751650000083
为调频速率,B为调频带宽,T为调频周期;
主控计算机将采集到的N个通道点云图像合并,获取一个调频周期T时间内目标坐标的空间点集:
Figure FDA0003311751650000084
P(xP_n,yP_n,zP_n)为目标P在第n个通道探测到的空间坐标系中的坐标位置;
该空间点集∑P即为最终显示的3D点云图像。
4.根据权利要求3所述的多通道高重频大动态范围测距测速激光雷达方法,其特征在于:所述同相信号和正交信号分别经过低通滤波器滤波去除串扰信号,具体是,所述串扰信号表示为:
Figure FDA0003311751650000091
ES_m表示为串扰的回波光束振幅,ELO_n表示为第n个回波光束对应的本振光束振幅;
串扰信号数据中包含频率为|fshift_n-fshift_m|项,通过低通滤波器对其进行滤波处理,消除高频串扰信号,提高单路输出光束的探测精度。
5.实现如权利要求1-4任一项所述的多通道高重频大动态范围测距测速激光雷达方法的装置,其特征在于:包括激光光源(1),所述激光光源(1)的输出端依次连接有起偏器(2)和激光放大器(3),所述激光放大器(3)的输出端经1×N分束器(4)连接有N个相位调制器(5);每个所述的相位调制器(5)均连接有移相器(17),N个移相器(17)一同连接有射频信号源(18),所述相位调制器(5)经光学滤波器(19)连接有移频器(6),移频器(6)的输出端连接有1×2分束器(7),所述1×2分束器(7)连接有光学环形器(8);
所述光学环形器(8)的输出端依次连接有光学望远镜(9)和光束定向器(10);所述光学环形器(8)和1×2分束器(7)一同连有光学桥接器(11);所述光学桥接器(11)的输出端依次连接有光电平衡探测器(12)和低通滤波器(13);所述低通滤波器(13)经模数转换器(14)连接有现场可编程门阵列(15),所述现场可编程门阵列(15)的输出端连接有主控计算机(16);所述主控计算机(16)还与光束定向器(10)连接。
CN202010354911.XA 2020-04-29 2020-04-29 多通道高重频大动态范围测距测速激光雷达方法及装置 Active CN111337902B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010354911.XA CN111337902B (zh) 2020-04-29 2020-04-29 多通道高重频大动态范围测距测速激光雷达方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010354911.XA CN111337902B (zh) 2020-04-29 2020-04-29 多通道高重频大动态范围测距测速激光雷达方法及装置

Publications (2)

Publication Number Publication Date
CN111337902A CN111337902A (zh) 2020-06-26
CN111337902B true CN111337902B (zh) 2021-11-23

Family

ID=71186531

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010354911.XA Active CN111337902B (zh) 2020-04-29 2020-04-29 多通道高重频大动态范围测距测速激光雷达方法及装置

Country Status (1)

Country Link
CN (1) CN111337902B (zh)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111751834B (zh) * 2020-06-30 2024-02-20 重庆大学 基于光学调频干涉与单频干涉的高速高精度动态测距方法
CN113589304B (zh) * 2020-09-27 2022-08-16 北京一径科技有限公司 激光雷达用抗串扰处理方法、装置和存储介质
CN112363146B (zh) * 2020-11-13 2022-07-29 杭州爱莱达科技有限公司 双调制共光路线性调频激光雷达测距测速方法及装置
CN112649812B (zh) * 2020-11-24 2023-08-15 中国电子科技集团公司第十一研究所 一种非匹配滤波的单脉冲多通道数据处理方法及装置
CN113192518B (zh) * 2021-04-14 2023-03-17 上海交通大学 毫米波侦听方法及***
CN115481361A (zh) * 2021-05-31 2022-12-16 华为技术有限公司 一种光计算装置和方法
WO2023065149A1 (zh) * 2021-10-20 2023-04-27 深圳市速腾聚创科技有限公司 激光雷达及激光雷达控制方法
CN114035174A (zh) * 2021-12-15 2022-02-11 浙江科技学院 双通道双啁啾线性调频连续波激光雷达方法及装置
CN114466494A (zh) * 2022-01-20 2022-05-10 南京矽典微***有限公司 一种智能感知灯控制方法、装置及设备
CN114879208B (zh) * 2022-06-30 2022-11-01 北京摩尔芯光半导体技术有限公司 一种调频连续波激光雷达***
CN118068351A (zh) * 2023-11-22 2024-05-24 深圳市中图仪器股份有限公司 用于测距装置的延时网络
CN118033626B (zh) * 2024-04-11 2024-06-14 济南卓伦智能交通技术有限公司 一种基于双雷达的目标跟踪测速方法及***

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103412313A (zh) * 2013-07-30 2013-11-27 桂林理工大学 低空轻小型面阵激光雷达测量***
CN110988864A (zh) * 2019-12-11 2020-04-10 中国电子科技集团公司第二十研究所 一种频率捷变mti雷达测速方法

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN202281835U (zh) * 2011-10-18 2012-06-20 中国科学院上海技术物理研究所 亚毫米级线性调谐激光测距***
CN105137446A (zh) * 2015-09-22 2015-12-09 中国科学院上海技术物理研究所 基于光学正交解调的相干零差多普勒测速激光雷达***
CN205374743U (zh) * 2015-09-22 2016-07-06 中国科学院上海技术物理研究所 一种光学正交解调的相干零差多普勒测速激光雷达***
CN105510927B (zh) * 2015-11-25 2017-09-26 中国科学院上海光学精密机械研究所 机载分频激光三维成像装置及其成像方法
CN207380238U (zh) * 2017-11-22 2018-05-18 杭州爱莱达科技有限公司 一种基于衍射光学的激光雷达
US11650296B2 (en) * 2018-02-16 2023-05-16 Xiaotian Steve Yao Optical sensing based on wavelength division multiplexed (WDM) light at different wavelengths in light detection and ranging LiDAR systems
CN110531341B (zh) * 2018-05-26 2023-04-11 苏州镭智传感科技有限公司 一种使用宽带调频激光的fmcw激光雷达及其扫描方法
CN110221292B (zh) * 2019-04-01 2020-12-11 清华大学 一种微波光子多波段雷达成像***及方法
CN111045028B (zh) * 2019-12-11 2020-06-23 杭州爱莱达科技有限公司 运动多普勒频移补偿的语音相干激光雷达侦听方法及装置
CN110780310B (zh) * 2019-12-31 2020-04-14 杭州爱莱达科技有限公司 偏振分集双通道测速及测距相干激光雷达测量方法及装置
CN110806586B (zh) * 2020-01-08 2020-04-14 杭州爱莱达科技有限公司 无扫描线性调频连续波测速测距激光三维成像方法及装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103412313A (zh) * 2013-07-30 2013-11-27 桂林理工大学 低空轻小型面阵激光雷达测量***
CN110988864A (zh) * 2019-12-11 2020-04-10 中国电子科技集团公司第二十研究所 一种频率捷变mti雷达测速方法

Also Published As

Publication number Publication date
CN111337902A (zh) 2020-06-26

Similar Documents

Publication Publication Date Title
CN111337902B (zh) 多通道高重频大动态范围测距测速激光雷达方法及装置
CN105629258B (zh) 基于伪随机码相位调制和外差探测的测速测距***及方法
US11366228B2 (en) Method and system for time separated quadrature detection of doppler effects in optical range measurements
US5889490A (en) Method and apparatus for improved ranging
CN110806586B (zh) 无扫描线性调频连续波测速测距激光三维成像方法及装置
CN110780310B (zh) 偏振分集双通道测速及测距相干激光雷达测量方法及装置
US7342651B1 (en) Time modulated doublet coherent laser radar
EP3637134B1 (en) Laser radar device
CN111693988A (zh) 一种激光毫米波一体化测距测速雷达方法及装置
US5745437A (en) Method and apparatus for coherent burst ranging
CN111650601B (zh) 车载相干激光雷达高分辨3d成像方法及装置
EP3465269A1 (en) Coherent lidar system using tunable carrier-suppressed single-sideband modulation
CN110794421B (zh) 伪随机码时延自差分干涉三维成像激光雷达方法及装置
CN115015886B (zh) 双通道伪随机码调相测距测速激光雷达方法及装置
CN108303704B (zh) 一种基于偏振调制的激光测量方法及激光雷达
US11500062B2 (en) Acceleration-based fast SOI processing
CN111999739A (zh) 相位调制测距测速相干激光雷达方法及装置
WO2021036315A1 (zh) 微波与激光雷达一体化方法及装置
CN109541636B (zh) 一种无盲区高距离分辨率激光雷达测风***及方法
CA3048330A1 (en) Method for processing a signal from a coherent lidar in order to reduce noise and related lidar system
CN101788671B (zh) 应用于外差探测啁啾调幅激光测距装置的多周期调制方法
CN114035174A (zh) 双通道双啁啾线性调频连续波激光雷达方法及装置
CN112327319A (zh) 基于循环移频环的固态激光雷达探测方法及***
CN113589309B (zh) 一种基于微波光子技术的非合作目标空间位置探测方法
CN104111450A (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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20220921

Address after: Room 1803, No. 12, Lane 333, Hongshi Road, Juyuan New District, Jiading District, Shanghai, 201800

Patentee after: Sun Jianfeng

Address before: 310053 room 701, block B, 581 torch Avenue, Puyan street, Binjiang District, Hangzhou City, Zhejiang Province

Patentee before: HANGZHOU ILIDAR TECHNOLOGY CO.,LTD.