CN105487067B - 粗测和精测距离信号处理方法、处理模块及基于该模块的啁啾调制光子计数激光雷达系统 - Google Patents
粗测和精测距离信号处理方法、处理模块及基于该模块的啁啾调制光子计数激光雷达系统 Download PDFInfo
- Publication number
- CN105487067B CN105487067B CN201510990679.8A CN201510990679A CN105487067B CN 105487067 B CN105487067 B CN 105487067B CN 201510990679 A CN201510990679 A CN 201510990679A CN 105487067 B CN105487067 B CN 105487067B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- signal
- mfrac
- tau
- 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
- 238000005259 measurement Methods 0.000 title claims abstract description 120
- 238000012545 processing Methods 0.000 title claims abstract description 21
- 238000003672 processing method Methods 0.000 title abstract description 5
- 238000001228 spectrum Methods 0.000 claims abstract description 18
- 238000000034 method Methods 0.000 claims abstract description 8
- 230000003287 optical effect Effects 0.000 claims description 24
- 238000004364 calculation method Methods 0.000 claims description 13
- 230000010354 integration Effects 0.000 claims description 10
- 238000001914 filtration Methods 0.000 claims description 7
- 230000004927 fusion Effects 0.000 claims description 6
- 238000012805 post-processing Methods 0.000 abstract description 5
- 238000001514 detection method Methods 0.000 description 8
- 229920006395 saturated elastomer Polymers 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000035945 sensitivity Effects 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
- G01S7/4802—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Optical Radar Systems And Details Thereof (AREA)
- Measurement Of Optical Distance (AREA)
Abstract
粗测和精测距离信号处理方法、处理模块及基于该模块的啁啾调制光子计数激光雷达系统，涉及激光雷达技术领域，具体涉及基于相位后处理方法的啁啾调制光子计数激光雷达测距技术领域。为了解决现有啁啾调制光子计数雷达测距误差大的问题。先通过中频频谱质心算法获得一个粗测的距离值，再通过对中频波形的相位后处理获得一个精细的距离测量值，将粗测值和精细值的互补在一起，从而有效的提高测距精度。本发明适用于啁啾调制光子计数激光雷达测距系统。
Description
技术领域
本发明涉及激光雷达技术领域，具体涉及基于相位后处理方法的啁啾调制光子计数激光雷达测距技术。
背景技术
啁啾调制光子计数雷达是一种新型的雷达体制，它结合了光子计数和啁啾调制两大技术，这使得它既拥有Gm-APD的单光子响应的极高探测灵敏度，能极大的增加探测距离，还拥有啁啾调制外差探测高精度的特点。但是啁啾调制光子雷达探测器采用的是Gm-APD，它由于工作在盖革模式下，信号的到来会引起雪崩效应造成饱和输出电流，如果不及时抑制，饱和电流将会击穿探测器，这就需要一定的时间来抑制饱和电流、并将探测器重置到盖革模式以准备下一次的探测，这个时间就是死时间。由于死时间的存在造成了Gm-APD啁啾调制光子雷达的探测是离散采样的，因此经过傅里叶变换的处理，中频频谱也是离散的，对应的距离间隔为δR＝c/2B，其中，δR为中频频率的固有间隔，c为光速，B为啁啾调制信号的带宽。由于无法准确的确定目标位于间隔内的位置，从而造成大的测距误差。
发明内容
本发明是为了解决现有啁啾调制光子计数雷达测距误差大的问题，从而提供粗测和精测距离信号处理方法、处理模块及基于该模块的啁啾调制光子计数激光雷达系统。
粗测和精测距离信号处理方法，该方法包括以下步骤：
参数设定步骤：
设定f0＝B，f0是啁啾调制信号的载波频率，即基频，B是啁啾调制信号的带宽；
粗测距离值计算步骤：
获得中频频域信号，从中频频谱上获得中频峰的数据(wm，Pm)，m是整数，表示第m个测量点，wm表示第m个测量点的频率位置，Pm表示第m个测量点的中频频谱强度；
采用质心算法估计粗测距离值，中频峰值的频率fIF|WCLA表示为：
其中，l是中频信号的半峰宽度；
粗测距离值Rraw为：
Rraw＝(fIF|WCLA/k)·c/2 (2)
其中，k＝B/T，k是啁啾调制信号的斜率，T是啁啾调制信号的时间长度，c为光速；
精测距离值计算步骤：
获得中频时域信号SIF(t)，
其中，M是激光信号衰减的系数，I0是发射激光的强度，ε(t)表示噪声；
根据粗测距离值计算步骤获得的Rraw，得到粗测的回波延迟时间τ，τ＝2Rraw/c；根据回波延迟时间τ，产生I/Q信号，该信号的Q分量信号、I分量信号分别为
SQ(t)＝cos(kτt) (4a)
SI(t)＝cos(kτt+π/2) (4b)
中频时域信号SIF(t)分别与I分量信号和Q分量信号相乘，并经滤波，得到的Q分量信号和I分量信号的积分结果为：
Q分量信号和I分量信号的积分结果相除，得到
其中，n是非负整数，表示重复循环的周期数；则相位差从而相位精测的距离值Rfine为：
其中，ΔR＝c/(2f0)，ΔR是相位测距的周期；
粗测和精测的距离值融合步骤：
利用粗测距离值将表示真正回波峰值的精测的距离值即目标的距离值R挑选出来，目标的距离值R是一系列精测的距离值Rfine中距离粗测距离值Rraw最近的一个精测的距离值。
粗测和精测距离信号处理模块，该模块包括以下模块：
参数设定模块：
设定f0＝B，f0是啁啾调制信号的载波频率，即基频，B是啁啾调制信号的带宽；
粗测距离值计算模块：
获得中频频域信号，从中频频谱上获得中频峰的数据(wm，Pm)，m是整数，表示第m个测量点，wm表示第m个测量点的频率位置，Pm表示第m个测量点的中频频谱强度；
采用质心算法估计粗测距离值，中频峰值的频率fIF|WCLA表示为：
其中，l是中频信号的半峰宽度；
粗测距离值Rraw为：
Rraw＝(fIF|WCLA/k)·c/2 (2)
其中，k＝B/T，k是啁啾调制信号的斜率，T是啁啾调制信号的时间长度，c为光速；
精测距离值计算模块：
获得中频时域信号SIF(t)，
其中，M是激光信号衰减的系数，I0是发射激光的强度，ε(t)表示噪声；
根据粗测距离值计算模块获得的Rraw，得到粗测的回波延迟时间τ，τ＝2Rraw/c；根据回波延迟时间τ，产生I/Q信号，该信号的Q分量信号、I分量信号分别为
SQ(t)＝cos(kτt) (4a)
SI(t)＝cos(kτt+π/2) (4b)
中频时域信号SIF(t)分别与I分量信号和Q分量信号相乘，并经滤波，得到的Q分量信号和I分量信号的积分结果为：
经过I/Q除法器得到
其中，n是非负整数，表示重复循环的周期数；
则相位差从而相位精测的距离值Rfine为：
其中，ΔR＝c/(2f0)，ΔR是相位测距的周期；
粗测和精测的距离值融合模块：
利用粗测距离值将表示真正回波峰值的精测的距离值即目标的距离值R挑选出来，目标的距离值R是一系列精测的距离值Rfine中距离粗测距离值Rraw最近的一个精测的距离值。
基于上述模块的啁啾调制光子计数激光雷达系统，它包括啁啾信号发生器、激光器、发射光学系统、接收光学系统、Gm-APD探测器、混频模块、低通滤波器、傅里叶变换器及信号处理器；
啁啾信号发生器的控制信号输出端一连接激光器的控制信号输入端，激光器出射的激光经发射光学系统的准直和扩束后发射，接收光学系统接收目标反射回的激光，接收光学系统的输出端连接Gm-APD探测器的输入端，Gm-APD探测器的输出端连接混频模块的光信号输入端，啁啾信号发生器的控制信号输出端二连接混频模块的电信号输入端，混频模块的输出端连接低通滤波器的输入端，低通滤波器的输出端同时连接傅里叶变换器的输入端及信号处理器的输入端一，傅里叶变换器的输出端连接信号处理器的输入端二；
信号处理器内嵌入有软件实现的粗测和精测距离信号处理模块。
由于Gm-APD的死时间，Gm-APD的采样探测是离散的，因此外差和傅里叶变化后的中频频谱也是离散的。当中频频谱的峰值在两个离散测量点之间，由于无法准确的给出峰值的位置从而造成了测距误差，限制了测距精度。本发明所述的粗测和精测距离信号处理方法，先通过中频频谱质心算法获得一个粗测的距离值，再通过对中频波形的相位后处理获得一个精细的距离测量值，将粗测值和精细值的互补在一起，从而有效的提高测距精度。
本发明所述的粗测和精测距离信号处理模块，先通过粗测距离值模块获得一个粗测的距离值，再通过精测距离值模块获得一个精细的距离测量值，将粗测值和精细值的互补在一起，得到的目标的距离值即为所求距离值，本发明能准确的确定目标位于频谱间隔内的位置，有效的提高测距精度。
本发明所述的基于粗测和精测距离信号处理模块的啁啾调制光子计数激光雷达系统，首先由啁啾信号发生器产生啁啾调制的电信号控制激光器产生幅度受啁啾调制的激光信号，该被调制的激光信号经过发射光学系统的准直和扩束后发射出去，经过往返大气的衰减，由目标反射回来的激光信号到达接收光学系统，然后通过接收光学系统将回波信号收集到Gm-APD探测器上，回波信号光子的到达率受啁啾信号发生器调制，Gm-APD响应回波光子到达率输出一系列疏密相间的脉冲串，脉冲密集的地方信号光子到达率高，反之脉冲稀疏的地方信号光子到达率低，这样Gm-APD探测器探测的结果携带着调制信息，与啁啾信号发生器产生的另一路啁啾调制电信号在混频器进行混频，再经过低通滤波器滤除高频的信号得到中频信号的时域信号，再经过傅里叶变换器的转换又可以得到中频信号的频域信号，最后将中频信号时域和频域信号同时输入到信号处理器进行数据后处理。本发明所述的系统将粗测值和精细值互补在一起，从而有效的提高啁啾调制光子计数激光雷达测距精度。
附图说明
图1是具体实施方式一所述的粗测和精测距离信号处理方法的原理图；
图2是具体实施方式一中的中频频谱和相位差精测值曲线图；
图3是具体实施方式三所述的基于粗测和精测距离信号处理模块的啁啾调制光子计数激光雷达系统的结构示意图。
具体实施方式
具体实施方式一：参照图1和图2具体说明本实施方式，本实施方式所述的粗测和精测距离信号处理方法，该方法包括以下步骤：
参数设定步骤：
设定f0＝B，f0是啁啾调制信号的载波频率，即基频，B是啁啾调制信号的带宽；
粗测距离值计算步骤：
获得中频频域信号，从中频频谱上获得中频峰的数据(wm，Pm)，m是整数，表示第m个测量点，wm表示第m个测量点的频率位置，Pm表示第m个测量点的中频频谱强度；
采用质心算法估计粗测距离值，中频峰值的频率fIF|WCLA表示为：
其中，l是中频信号的半峰宽度；
粗测距离值Rraw为：
Rraw＝(fIF|WCLA/k)·c/2 (2)
其中，k＝B/T，k是啁啾调制信号的斜率，T是啁啾调制信号的时间长度，c为光速；
精测距离值计算步骤：
获得中频时域信号SIF(t)，
其中，M是激光信号衰减的系数，I0是发射激光的强度，ε(t)表示噪声；
根据粗测距离值计算步骤获得的Rraw，得到粗测的回波延迟时间τ，τ＝2Rraw/c；根据回波延迟时间τ，产生I/Q信号，该信号的Q分量信号、I分量信号分别为
SQ(t)＝cos(kτt) (4a)
SI(t)＝cos(kτt+π/2) (4b)
中频时域信号SIF(t)分别与I分量信号和Q分量信号相乘，并经滤波，得到的Q分量信号和I分量信号的积分结果为：
Q分量信号和I分量信号的积分结果相除，得到
其中，n是非负整数，表示重复循环的周期数；则相位差从而相位精测的距离值Rfine为：
其中，ΔR＝c/(2f0)，ΔR是相位测距的周期；
粗测和精测的距离值融合步骤：
利用粗测距离值将表示真正回波峰值的精测的距离值即目标的距离值R挑选出来，目标的距离值R是一系列精测的距离值Rfine中距离粗测距离值Rraw最近的一个精测的距离值。
计算机产生一个I/Q信号，它的频率为kτ，它的初始相位0对应于啁啾信号发射的0时刻，I/Q分量经过正交检相探测得到Q分量信号和I分量信号的积分结果，再经反正切的解算得到相位差。由于三角函数是周期性函数，所以相位测距量程周期性重复，它的周期为ΔR＝c/(2f0)，为了将相位和频率两种测距方法结合在一起，令相位方法的重复周期ΔR等于中频频率的固有间隔δR，即ΔR＝δR，由于δR＝c/(2B)以及ΔR＝c/(2f0)，那么需要满足f0＝B，如图2所示，图2中a为b的局部放大图，A、B、C、D和E所在曲线为中频频谱图，A、B、C、D和E为五个测量点，相邻测量点频率的间隔为δR，a为相位差精测值曲线图，相位差精测值能给出中频峰值在两个测量点间隔内的具体位置，但是无法给出中频峰值的相位差精测值具体在哪个间隔(周期)内，从而得到一系列相位差精测值，如图a中的F、G和H，最后利用粗测距离值将真正中频峰值的相位差精测值挑选出来,，该相位差精测值对应的距离值即为目标距离值R。
具体实施方式二：粗测和精测距离信号处理模块，该模块包括以下模块：
参数设定模块：
设定f0＝B，f0是啁啾调制信号的载波频率，即基频，B是啁啾调制信号的带宽；
粗测距离值计算模块：
获得中频频域信号，从中频频谱上获得中频峰的数据(wm，Pm)，m是整数，表示第m个测量点，wm表示第m个测量点的频率位置，Pm表示第m个测量点的中频频谱强度；
采用质心算法估计粗测距离值，中频峰值的频率fIF|WCLA表示为：
其中，l是中频信号的半峰宽度；
粗测距离值Rraw为：
Rraw＝(fIF|WCLA/k)·c/2 (2)
其中，k＝B/T，k是啁啾调制信号的斜率，T是啁啾调制信号的时间长度，c为光速；
精测距离值计算模块：
获得中频时域信号SIF(t)，
其中，M是激光信号衰减的系数，I0是发射激光的强度，ε(t)表示噪声；
根据粗测距离值计算模块获得的Rraw，得到粗测的回波延迟时间τ，τ＝2Rraw/c；根据回波延迟时间τ，产生I/Q信号，该信号的Q分量信号、I分量信号分别为
SQ(t)＝cos(kτt) (4a)
SI(t)＝cos(kτt+π/2) (4b)
中频时域信号SIF(t)分别与I分量信号和Q分量信号相乘，并经滤波，得到的Q分量信号和I分量信号的积分结果为：
经过I/Q除法器得到
其中，n是非负整数，表示重复循环的周期数；则相位差从而相位精测的距离值Rfine为：
其中，ΔR＝c/(2f0)，ΔR是相位测距的周期；
粗测和精测的距离值融合模块：
利用粗测距离值将表示真正回波峰值的精测的距离值即目标的距离值R挑选出来，目标的距离值R是一系列精测的距离值Rfine中距离粗测距离值Rraw最近的一个精测的距离值。
在软件设计中，从一系列精测的距离值Rfine中选取距离粗测距离值Rraw最近的一个精测的距离值，是从一系列精测的距离值Rfine|n＝0,1,2,…中在由Rraw-ΔR/2至Rraw+ΔR/2的范围内取交集，即：R＝{Rfine|n＝0,1,2,…}∩(Rraw-ΔR/2,Rraw+ΔR/2)，从而实现从一系列精测的距离值Rfine中选取距离粗测距离值Rraw最近的一个精测的距离值。
具体实施方式三：结合图3具体说明本实施方式，基于粗测和精测距离信号处理模块的啁啾调制光子计数激光雷达系统，它包括啁啾信号发生器1、激光器2、发射光学系统3、接收光学系统4、Gm-APD探测器5、混频模块6、低通滤波器7、傅里叶变换器8及信号处理器9；
啁啾信号发生器1的控制信号输出端一连接激光器2的控制信号输入端，激光器2出射的激光经发射光学系统3的准直和扩束后发射，接收光学系统4接收目标反射回的激光，接收光学系统4的输出端连接Gm-APD探测器5的输入端，Gm-APD探测器5的输出端连接混频模块6的光信号输入端，啁啾信号发生器1的控制信号输出端二连接混频模块6的电信号输入端，混频模块6的输出端连接低通滤波器7的输入端，低通滤波器7的输出端同时连接傅里叶变换器8的输入端及信号处理器9的输入端一，傅里叶变换器8的输出端连接信号处理器9的输入端二；
信号处理器9内嵌入有软件实现的粗测和精测距离信号处理模块。
具体实施方式四：本实施方式是对具体实施方式一所述的基于粗测和精测距离信号处理模块的啁啾调制光子计数激光雷达系统作进一步说明，本实施方式中，接收光学系统4还设置有窄带滤光片。窄带滤光片能够滤除非工作波长的背景噪声。
Claims (4)
1.粗测和精测距离信号处理方法，其特征在于，该方法包括以下步骤：
参数设定步骤：
设定f0＝B，f0是啁啾调制信号的载波频率，即基频，B是啁啾调制信号的带宽；
粗测距离值计算步骤：
获得中频频域信号，从中频频谱上获得中频峰的数据(wm，Pm)，m是整数，表示第m个测量点，wm表示第m个测量点的频率位置，Pm表示第m个测量点的中频频谱强度；
采用质心算法估计粗测距离值，中频峰值的频率fIF|WCLA可以表示为：
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mi>I</mi>
<mi>F</mi>
</mrow>
</msub>
<msub>
<mo>|</mo>
<mrow>
<mi>W</mi>
<mi>C</mi>
<mi>L</mi>
<mi>A</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mo>-</mo>
<mi>l</mi>
</mrow>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mi>l</mi>
</mrow>
</munderover>
<msub>
<mi>P</mi>
<mi>m</mi>
</msub>
<msub>
<mi>w</mi>
<mi>m</mi>
</msub>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mo>-</mo>
<mi>l</mi>
</mrow>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mi>l</mi>
</mrow>
</munderover>
<msub>
<mi>P</mi>
<mi>m</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，l是中频频域信号的半峰宽度；
粗测距离值Rraw为：
Rraw＝(fIF|WCLA/k)·c/2 (2)
其中，k＝B/T，k是啁啾调制信号的斜率，T是啁啾调制信号的时间长度，c为光速；
精测距离值计算步骤：
获得中频时域信号SIF(t)，
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>I</mi>
<mi>F</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msub>
<mi>MI</mi>
<mn>0</mn>
</msub>
<mi>r</mi>
<mi>e</mi>
<mi>c</mi>
<mi>t</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mrow>
<mi>T</mi>
<mo>-</mo>
<mi>&tau;</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mi>&tau;</mi>
<mo>+</mo>
<mi>k</mi>
<mi>t</mi>
<mi>&tau;</mi>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mi>k&tau;</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>&epsiv;</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，M是激光信号衰减的系数，I0是发射激光信号的强度，ε(t)表示噪声；
根据粗测距离值计算步骤获得的Rraw，得到粗测的回波延迟时间τ，τ＝2Rraw/c；根据回波延迟时间τ，产生I/Q信号，该信号的Q分量信号、I分量信号分别为
SQ(t)＝cos(kτt) (4a)
SI(t)＝cos(kτt+π/2) (4b)
中频时域信号SIF(t)分别与I分量信号和Q分量信号相乘，并经滤波，得到的Q分量信号和I分量信号的积分结果为：
<mrow>
<mi>Q</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msub>
<mi>MI</mi>
<mn>0</mn>
</msub>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mi>&tau;</mi>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mi>k&tau;</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mi>a</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>I</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msub>
<mi>MI</mi>
<mn>0</mn>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mi>&tau;</mi>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mi>k&tau;</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mi>b</mi>
<mo>)</mo>
</mrow>
</mrow>
Q分量信号和I分量信号的积分结果相除，可以得到
<mrow>
<mi>t</mi>
<mi>a</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mi>&tau;</mi>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mi>k&tau;</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>I</mi>
<mo>/</mo>
<mi>Q</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，n是非负整数，表示重复循环的周期数；则相位差从而相位精测的距离值Rfine为：
其中，ΔR＝c/(2f0)，ΔR是相位测距的周期；
粗测和精测的距离值融合步骤：利用粗测距离值将目标的距离值R挑选出来，目标的距离值R是一系列精测的距离值Rfine中距离粗测距离值Rraw最近的一个精测的距离值。
2.粗测和精测距离信号处理模块，其特征在于，该模块包括以下模块：
参数设定模块：
设定f0＝B，f0是啁啾调制信号的载波频率，即基频，B是啁啾调制信号的带宽；
粗测距离值计算模块：
获得中频频域信号，从中频频谱上获得中频峰的数据(wm，Pm)，m是整数，表示第m个测量点，wm表示第m个测量点的频率位置，Pm表示第m个测量点的中频频谱强度；
采用质心算法估计粗测距离值，中频峰值的频率fIF|WCLA表示为：
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mi>I</mi>
<mi>F</mi>
</mrow>
</msub>
<msub>
<mo>|</mo>
<mrow>
<mi>W</mi>
<mi>C</mi>
<mi>L</mi>
<mi>A</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mo>-</mo>
<mi>l</mi>
</mrow>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mi>l</mi>
</mrow>
</munderover>
<msub>
<mi>P</mi>
<mi>m</mi>
</msub>
<msub>
<mi>w</mi>
<mi>m</mi>
</msub>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mo>-</mo>
<mi>l</mi>
</mrow>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mi>l</mi>
</mrow>
</munderover>
<msub>
<mi>P</mi>
<mi>m</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，l是中频频域信号的半峰宽度；
粗测距离值Rraw为：
Rraw＝(fIF|WCLA/k)·c/2 (2)
其中，k＝B/T，k是啁啾调制信号的斜率，T是啁啾调制信号的时间长度，c为光速；
精测距离值计算模块：
获得中频时域信号SIF(t)，
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>I</mi>
<mi>F</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msub>
<mi>MI</mi>
<mn>0</mn>
</msub>
<mi>r</mi>
<mi>e</mi>
<mi>c</mi>
<mi>t</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mrow>
<mi>T</mi>
<mo>-</mo>
<mi>&tau;</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mi>&tau;</mi>
<mo>+</mo>
<mi>k</mi>
<mi>t</mi>
<mi>&tau;</mi>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mi>k&tau;</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>&epsiv;</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，M是激光信号衰减的系数，I0是发射激光的强度，ε(t)表示噪声；
根据粗测距离值计算模块获得的Rraw，得到粗测的回波延迟时间τ，τ＝2Rraw/c；根据回波延迟时间τ，产生I/Q信号，该信号的Q分量信号、I分量信号分别为
SQ(t)＝cos(kτt) (4a)
SI(t)＝cos(kτt+π/2) (4b)
中频时域信号SIF(t)分别与I分量信号和Q分量信号相乘，并经滤波，得到的Q分量信号和I分量信号的积分结果为：
<mrow>
<mi>Q</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msub>
<mi>MI</mi>
<mn>0</mn>
</msub>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mi>&tau;</mi>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mi>k&tau;</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mi>a</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>I</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msub>
<mi>MI</mi>
<mn>0</mn>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mi>&tau;</mi>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mi>k&tau;</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mi>b</mi>
<mo>)</mo>
</mrow>
</mrow>
经过I/Q除法器得到
<mrow>
<mi>t</mi>
<mi>a</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mi>&tau;</mi>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mi>k&tau;</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>I</mi>
<mo>/</mo>
<mi>Q</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，n是非负整数，表示重复循环的周期数；则相位差从而相位精测的距离值Rfine为：
其中，ΔR＝c/(2f0)，ΔR是相位测距的周期；
粗测和精测的距离值融合模块：
利用粗测距离值将目标的距离值R挑选出来，目标的距离值R是一系列精测的距离值Rfine中距离粗测距离值Rraw最近的一个精测的距离值。
3.基于权利要求2所述的粗测和精测距离信号处理模块的啁啾调制光子计数激光雷达系统，其特征在于，所述啁啾调制光子计数激光雷达系统包括啁啾信号发生器(1)、激光器(2)、发射光学系统(3)、接收光学系统(4)、Gm-APD探测器(5)、混频模块(6)、低通滤波器(7)、傅里叶变换器(8)及信号处理器(9)；
啁啾信号发生器(1)的控制信号输出端一连接激光器(2)的控制信号输入端，激光器(2)出射的激光经发射光学系统(3)的准直和扩束后发射，接收光学系统(4)接收目标反射回的激光，接收光学系统(4)的输出端连接Gm-APD探测器(5)的输入端，Gm-APD探测器(5)的输出端连接混频模块(6)的光信号输入端，啁啾信号发生器(1)的控制信号输出端二连接混频模块(6)的电信号输入端，混频模块(6)的输出端连接低通滤波器(7)的输入端，低通滤波器(7)的输出端同时连接傅里叶变换器(8)的输入端及信号处理器(9)的输入端一，傅里叶变换器(8)的输出端连接信号处理器(9)的输入端二；
信号处理器(9)内嵌入有软件实现的粗测和精测距离信号处理模块。
4.根据权利要求3所述的啁啾调制光子计数激光雷达系统，其特征在于，接收光学系统(4)还设置有窄带滤光片。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510990679.8A CN105487067B (zh) | 2015-12-25 | 2015-12-25 | 粗测和精测距离信号处理方法、处理模块及基于该模块的啁啾调制光子计数激光雷达系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510990679.8A CN105487067B (zh) | 2015-12-25 | 2015-12-25 | 粗测和精测距离信号处理方法、处理模块及基于该模块的啁啾调制光子计数激光雷达系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105487067A CN105487067A (zh) | 2016-04-13 |
| CN105487067B true CN105487067B (zh) | 2017-12-05 |
Family
ID=55674176
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510990679.8A Active CN105487067B (zh) | 2015-12-25 | 2015-12-25 | 粗测和精测距离信号处理方法、处理模块及基于该模块的啁啾调制光子计数激光雷达系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105487067B (zh) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105891840A (zh) * | 2016-06-29 | 2016-08-24 | 长春理工大学 | 激光测距系统 |
| CN106813784B (zh) * | 2017-03-03 | 2023-10-20 | 浙江工业大学 | 一种实时微波脉冲啁啾检测装置及其检测方法 |
| CN108445471B (zh) * | 2018-03-26 | 2019-08-02 | 武汉大学 | 一种单光子激光雷达多探测器条件下的测距精度评估方法 |
| CN108444948A (zh) * | 2018-04-10 | 2018-08-24 | 中国科学院上海技术物理研究所 | 测量大气二氧化碳浓度的差分吸收激光雷达系统及方法 |
| IL259190A (en) * | 2018-05-07 | 2018-06-28 | Arbe Robotics Ltd | System and method for frequency hopping MIMO FMCW imaging radar |
| WO2020107234A1 (zh) * | 2018-11-27 | 2020-06-04 | 华为技术有限公司 | 雷达系统及波形生成方法 |
| CN110542907B (zh) * | 2019-08-19 | 2021-04-20 | 杭州电子科技大学 | 双测距模块复合快速高精度激光三维测量方法 |
| CN111162839A (zh) * | 2019-12-31 | 2020-05-15 | 华东师范大学重庆研究院 | 一种远距离高精度光纤断点位置检测方法及系统 |
| CN113447946B (zh) * | 2021-06-28 | 2022-08-05 | 哈尔滨工业大学 | 微弱激光回波信号的微多普勒信息测量系统 |
| CN113534103B (zh) * | 2021-06-29 | 2024-04-09 | 昂纳科技(深圳)集团股份有限公司 | 激光频率自适应调制方法、系统、探测设备和存储介质 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN201876545U (zh) * | 2010-09-02 | 2011-06-22 | 淄博职业学院 | 脉冲相位式激光测距仪 |
-
2015
- 2015-12-25 CN CN201510990679.8A patent/CN105487067B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105487067A (zh) | 2016-04-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105487067B (zh) | 粗测和精测距离信号处理方法、处理模块及基于该模块的啁啾调制光子计数激光雷达系统 | |
| EP2198323B9 (en) | Time delay estimation | |
| US8848171B2 (en) | Highly accurate distance measurement device | |
| CN107907878B (zh) | 高精度获取fmcw雷达距离测量值的方法 | |
| CN109964143B (zh) | 用于处理由相干激光雷达引起的信号的方法及相关激光雷达系统 | |
| US10151827B2 (en) | Laser phase estimation and correction | |
| CA2583337C (en) | Electro-optical method for measuring distance and detecting a non-ideal chirp profile | |
| JP2016029369A (ja) | タイミング同期を備えた周波数変調連続波（ｆｍｃｗ）レーダー | |
| US11500062B2 (en) | Acceleration-based fast SOI processing | |
| CN101788671B (zh) | 应用于外差探测啁啾调幅激光测距装置的多周期调制方法 | |
| JP5235737B2 (ja) | パルスドップラレーダ装置 | |
| CN110187353A (zh) | 一种远目标激光测距方法 | |
| RU2688921C2 (ru) | Способ измерения дальности и радиальной скорости в РЛС с зондирующим составным псевдослучайным ЛЧМ импульсом | |
| CN110726995B (zh) | 激光雷达高精度测距方法及系统 | |
| KR100780525B1 (ko) | 레이저 광파기 | |
| CN116338708A (zh) | 飞行时间的测量方法、装置和系统 | |
| KR101359344B1 (ko) | Ｆｍｃｗ 기반 거리 측정 장치 | |
| Müller et al. | A method for the analysis of ramp-inherent linearity distortions in automotive radar applications | |
| JP7445107B2 (ja) | レーザ測距装置 | |
| Nejad et al. | A new design of laser phase-shift range finder independent of environmental conditions and thermal drift | |
| Sorochan et al. | Radio altimeter with J-correlation signal processing | |
| RU164241U1 (ru) | Устройство управления коэффициентом усиления передающего устройства рлс на основе восстановления параметров усиливаемого сигнала | |
| RU2439592C1 (ru) | Чм-дальномер непрерывного слежения с дробно-дифференцирующим фильтром | |
| CN115825975A (zh) | 基于脉冲调制型激光调频连续波的超分辨测距方法及系统 | |
| CN116299410A (zh) | 一种针对mfsk波形雷达的近距离目标模拟方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |