CN106574965A - 激光相位估计及校正 - Google Patents
激光相位估计及校正 Download PDFInfo
- Publication number
- CN106574965A CN106574965A CN201580045019.0A CN201580045019A CN106574965A CN 106574965 A CN106574965 A CN 106574965A CN 201580045019 A CN201580045019 A CN 201580045019A CN 106574965 A CN106574965 A CN 106574965A
- Authority
- CN
- China
- Prior art keywords
- totsegs
- phase deviation
- persistent period
- signal
- laser signal
- 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
Links
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
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/02—Systems using the reflection of electromagnetic waves other than radio waves
- G01S17/06—Systems determining position data of a target
- G01S17/08—Systems determining position data of a target for measuring distance only
- G01S17/32—Systems determining position data of a target for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated
- G01S17/34—Systems determining position data of a target for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated using transmission of continuous, frequency-modulated waves while heterodyning the received signal, or a signal derived therefrom, with a locally-generated signal related to the contemporaneously transmitted signal
-
- 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/491—Details of non-pulse systems
- G01S7/4912—Receivers
- G01S7/4915—Time delay measurement, e.g. operational details for pixel components; Phase measurement
-
- 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
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/02—Systems using the reflection of electromagnetic waves other than radio waves
- G01S17/06—Systems determining position data of a target
- G01S17/08—Systems determining position data of a target for measuring distance only
- G01S17/32—Systems determining position data of a target for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated
- G01S17/36—Systems determining position data of a target for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated with phase comparison between the received signal and the contemporaneously transmitted signal
-
- 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
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/02—Systems using the reflection of electromagnetic waves other than radio waves
- G01S17/50—Systems of measurement based on relative movement of target
- G01S17/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
-
- 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/481—Constructional features, e.g. arrangements of optical elements
- G01S7/4814—Constructional features, e.g. arrangements of optical elements of transmitters alone
-
- 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/491—Details of non-pulse systems
- G01S7/4911—Transmitters
-
- 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/497—Means for monitoring or calibrating
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)
- Measurement Of Optical Distance (AREA)
Abstract
在一个一般方面中,一种非暂时性计算机可读存储媒体可经配置以存储在被执行时造成处理器执行过程的指令。所述过程可包含:产生激光信号的区段,其中所述激光信号的所述区段具有持续时间;及基于所述激光信号而产生第一参考信号。所述过程可包含:基于所述第一参考信号而计算与所述持续时间的第一部分对应的第一相位偏差;及基于所述激光信号而产生第二参考信号。所述过程可包含:基于所述第二参考信号而计算与所述持续时间的第二部分对应的第二相位偏差;及基于所述第一相位偏差及所述第二相位偏差的组合而计算所述激光信号的所述区段的相位偏差。
Description
相关申请案
本申请案主张2015年7月23日申请的且标题为“激光相位估计及校正(LASERPHASE ESTIMATION AND CORRECTION)”的第14/807,508号美国申请案的优先权并是其接续申请案,其主张2014年7月25日申请的且标题为“激光相位估计及校正(LASER PHASEESTIMATION AND CORRECTION)”的第62/028,911号美国临时申请案的优先权及权益,所述两个申请案的全文以引用的方式并入本文中。
本申请案还主张2014年7月25日申请且标题为“激光相位估计及校正(LASERPHASE ESTIMATION AND CORRECTION)”的第62/028,911号美国临时申请案的优先权及权益,所述申请案的全文以引用的方式并入本文中。
技术领域
本发明涉及一种多参考激光检测及测距(LIDAR)系统。
背景技术
在一些已知LIDAR系统中,激光可用于估计移动物体的距离及速度。然而,距离及速度估计会因(例如)相位偏差、频率偏差及/或其它干涉而失真。因此,需要系统、方法及设备来解决现有技术的缺点且提供其它创新特征。
发明内容
在一个一般方面中,一种非暂时性计算机可读存储媒体可经配置以存储在被执行时致使处理器执行过程的指令。所述过程可包含:产生激光信号的区段,其中所述激光信号的所述区段具有持续时间;及基于所述激光信号而产生第一参考信号。所述过程可包含:基于所述第一参考信号而计算与所述持续时间的第一部分对应的第一相位偏差;及基于所述激光信号而产生第二参考信号。所述过程可包含:基于所述第二参考信号而计算与所述持续时间的第二部分对应的第二相位偏差;及基于所述第一相位偏差及所述第二相位偏差的组合而计算所述激光信号的所述区段的相位偏差。
附图说明
图1是说明激光系统的图式。
图2说明以校正偏差为目标的目标激光信号。
图3说明以校正偏差为目标的图2中所展示的目标激光信号的另一实例。
图4说明图3中所展示的目标激光信号的时间区段的位移的另一实例。
图5是说明根据实施方案的目标激光信号的另一实例的图式。
具体实施方式
图1是说明激光系统100(还可称为光检测及测距(LIDAR)系统)的图式,激光系统100经配置以使用激光源110来产生或测量可相对于激光系统100而固定或移动的物体5的距离及/或速度。在一些实施方案中,物体5可称为目标或目标物体5。激光系统100可用于调频连续波(FMCW)应用中。
激光系统100的激光源110经配置以发射(例如,产生、传播)一或多个频率下的电磁辐射,所述电磁辐射可为(例如)相干光发射(例如单色光发射)或光束。为简明起见,来自激光源110的发射将称为电磁辐射发射(例如电磁辐射发射)、发射激光信号10或发射光。
如图1中所展示,可由分离器125将激光信号10分离成多个激光信号,例如至少激光信号11-1、11-2、12A、12B、13A及13B。可由分离器125产生用于由参考系统190A、190B处理的激光信号12A、12B、13A及13B,参考系统190A、190B各自包含干涉计(其可包含经配置以将光学信号转换成电信号的一或多个光检测器或检测器(例如检测器150C))。在一些实施方案中,激光信号11可源自于分离激光信号且可称为组合激光信号。如图1中所展示,干涉计可用于产生激光信号11,所述激光信号11可由分析器170(其还可称为解调制器)分析以进行一或多个校正。在此类实施方案中,可(例如,由分离器125)将激光信号10进一步分离成激光信号11-1及激光信号11-2。激光信号11-1可从物体5反射为激光信号11-4。激光信号11-2可在延迟延时142C(其可与长度关联)之后成为激光信号11-3,且激光信号11-3可经由组合器140C而与激光信号11-4组合。来自干涉计的激光信号11(还可称为干涉计信号)可用于使用检测器150C来收集与激光信号11有关的信息。可将下文关于激光信号11的论述应用于可用于界定激光信号11的分量激光信号11-1到11-4中的任何者,激光信号11可为目标激光信号或以由分析器170分析为目标的激光信号。为简明起见,将分离器125说明为单个组件。在一些实施方案中,分离器125可包含一个以上分离器。类似地,图1中所展示的组合器中的一或多者可经组合或可包含额外组合器。
如图1中所展示,激光系统100包含频率扫描模块120。频率扫描模块120经配置以触发激光源110(例如)通过调制激光源110的驱动电流而产生各种光频(一般还可称为频率)。具体来说,频率扫描模块120经配置以触发激光源110产生光频的模式(还可称为频率模式)。举例来说,频率扫描模块120可经配置以触发激光源110以产生光频的正弦波型、光频的锯齿波型等等。在一些实施方案中,锯齿波型可具有光频连续增大(例如单调增大、线性增大、非线性增大)(还可称为渐增式线性调频(up-chirp))的部分且可具有光频连续减小(例如单调减小、线性减小、非线性减小)(还可称为渐减式线性调频(down-chirp))的部分。因此,频率模式可具有包含渐增式线性调频及渐减式线性调频的循环。
激光系统100包含组合器140C,其经配置以响应于从激光源110朝向物体5的发射激光信号11-1(从激光信号10分离)而接收从物体5反射的激光信号11-4(还可称为反射激光信号或散射激光信号)(未展示)。在一些实施方案中,来自物体5的反射激光信号(还可称为返回信号或返回光)可与发射激光信号10的部分(例如,延迟达延时142C的激光信号11-3)混合且接着由分析器170分析(在通过检测器150C而转换成电信号之后)。
激光系统100的分析器170经配置以分析来自激光源110的发射激光信号11-1及由组合器140C接收的反射激光信号11-4的组合。可根据包含渐增式线性调频及接着渐减式线性调频(或渐减式线性调频及接着渐增式线性调频)的模式而发射所发射激光信号11-1。可由分析器170分析来自激光源110的发射激光信号11-1的频率及由组合器140C接收的反射激光信号11-4的频率的组合以获得或界定拍频或差拍信号。换句话说,拍频可为在到物体5(发射激光信号)及返回(反射激光信号)的往返行程上的信号频率变化的总和,且可包含起因于激光系统100与物体5之间的相对距离运动的反射激光信号的多普勒频移。在一些实施方案中,差拍信号可具有相对较恒定频率或变频。在一些实施方案中,发射激光信号11-1的频率及反射激光信号11-4的频率的组合可称为差频、拍频或往返频率。
分析器170可经配置以计算往返时段,所述往返时段是从激光信号10的发射到反射激光信号的返回的接收的时段。发射激光信号11-1及反射激光信号11-4的组合可统称为往返激光信号。分析器170还可经配置以基于发射激光信号11-1及反射激光信号11-4的所述组合而计算距离及/或速度。
由于(例如)激光源110的驱动电流调制,激光输出的光学功率可在频率模式(例如频率扫描或渐增式线性调频/渐减式线性调频)期间显著改变。所述频率模式可由于可造成变化(例如频率、相位等等)的缺陷驱动电流信号、激光源110中的不可避免热激发等等而非理想地远离(例如,可偏离)指定频率模式。
激光系统100包含经配置以产生参考信号的参考系统190A、190B,所述参考信号可用于校正(例如)由激光源110产生的一或多个激光信号的频率偏差、相位偏差等等。换句话说,包含于激光系统100中的参考系统190A、190B可经配置以促进补偿(例如)来自激光系统100的发射激光信号11-1、反射激光信号11-4、往返激光信号及/或其它等等的频率模式的偏差(例如非线性、非理想性、误差)。参考系统190A、190B可用于实现近乎理想或理想的FMCW LIDAR实施方案。具体来说,参考系统190A、190B可用于校正偏差以获得相对较恒定拍频。以由参考系统190A、190B校正(例如,调整)为目标的激光信号可称为目标激光信号,且可至少包含发射激光信号10(或源自于其的信号)、反射激光信号(或源自于其的信号)及往返激光信号(或源自于其的信号)。
参考系统190A、190B中的每一者经配置以分别界定可用于确定(例如,识别、计算)一或多个目标激光信号(例如,激光信号11)的偏差的参考信号14A、14B。可由分离器125将激光信号10分离成用于由参考系统190A、190B处理的激光信号12A、12B。可基于激光信号12A及延迟激光信号13A'(其基于激光信号13A而产生)的组合(使用组合器140A)而产生参考信号14A。类似地,可基于激光信号12B及延迟激光信号13B'(其基于激光信号13B而产生)的组合(使用组合器140B)而产生参考信号14B。换句话说,参考信号14A、14B可为分别通过组合激光信号12A及延迟信号13A'及组合激光信号12B及延迟信号13B'而产生的差拍信号。分别通过延时142A、142B而产生延迟信号13A'、13B'。延时142A、142B可各自称为固定延迟或参考臂长度,且参考信号14A、14B可到称参考臂信号。延时142A、142B中的每一者可经配置以界定延迟时段,且可为干涉计的部分(例如,包含于干涉计中)。
偏差检测器150A、150B可经配置以分别确定与参考信号14A、14B相关联的偏差。在一些实施方案中,可将偏差检测器150A、150B组合成单个模块或分成多个模块。在一些实施方案中,偏差检测器150A、150B中的一或多者可经配置以检测包含相移等等的各种偏差。偏差检测器150A、150B中的一或多者可包含或可为光检测器。
因为分析器170无法经配置以直接测量一或多个目标激光信号的偏差,所以参考系统190A、190B可经配置以测量可与一或多个目标激光信号(例如激光信号11)的时间区段对应的参考信号14A、14B的时间区段的偏差(例如,使用偏差检测器150A、150B)。在一些实施方案中,与目标激光信号相关联的时间信号可称为目标时间区段。举例来说,参考系统190A、190B可经配置以测量与一或多个目标激光信号的时间区段对应的参考信号14A、14B的时间区段的激光相位历时。一般来说,参考信号14A、14B的希伯特(Hilbert)变换的相位历时可用于校正一或多个目标激光信号的相位历时的偏差,使得经校正的目标激光信号可为可从其确定所要频率的所要频调。
在此实施方案中,多个参考系统(参考系统190A、190B)用于测量不同时间区段的偏差。具体来说,延时142A可不同于延时142B,使得参考信号14A及14B将与不同时间区段(或延迟时段)相关联。因此,与不同时间区段中的每一者相关联的偏差可用于各种数学组合中来以相对较精确方式确定(例如,计算)一或多个目标信号的又一时间区段的偏差。图2及3中说明此概念的实例。包含于激光系统100中的多参考系统具有优于单参考信号处理系统或方法的许多优点,这是因为单参考信号处理系统可(例如)假设相位历程以估计激光相位历时。
图2说明以校正偏差为目标的目标激光信号。举例来说,所述目标激光信号在时间T0到时间T10之间的持续时间期间具有以分析为目标的时间区段20。对于参考信号A,可使用发生于时间T0到T7之间的时间区段20A来确定(例如,测量、计算)第一偏差,且对于参考信号B,可依据发生于时间T7到T10之间的时间区段20B而确定第二偏差。因此,可通过组合(例如,加总、联结)与时间区段20A及时间区段20B相关联的偏差而确定与所述目标激光信号相关联的时间区段20的总偏差。
图3说明以校正偏差为目标的图2中所展示的目标激光信号的另一实例。此实例是图2中所说明的实例的变型。在此实施方案中,目标激光信号在时间T0到时间T4之间的持续时间期间具有以分析为目标的时间区段21。对于参考信号A,可使用发生于时间T0到T7之间的时间区段21A来确定(例如,测量、计算)第一偏差,且对于参考信号B,可依据发生于时间T4到T7之间的时间区段21B而确定第二偏差。因此,可通过组合(其在此情况中是相减)与时间区段21A及时间区段21B相关联的偏差而确定与目标激光信号相关联的时间区段20的总偏差。
返回参考图1,与参考系统190A、190B相关联的时间区段可与分析器170的取样率对应。在一些实施方案中,时间区段可与取样时间间隔或取样时段的整数倍(还可称为整数)对应。举例来说,分析器170的所述取样率可包含数纳秒的取样时段。延时142A可经界定以与所述取样时段的整数倍对应(例如,5整数倍×5纳秒取样时段=25纳秒延时),使得与参考信号14A相关联的偏差可与与所述取样时段的所述整数倍对应的时间区段相关联。借此,与参考信号14A、14B的时间区段相关联的偏差可与与一或多个目标信号相关联的取样时段的时间区段匹配。因此,可精确确定目标激光信号的时间区段的偏差。
至少在(例如)图2中说明取样时段的实例。如图2中所展示,时间T0到T1之间的时段可与取样时段对应。因此,目标激光信号的区段20可包含10个取样时段。也如图2中所展示,区段20A包含7个取样时段RA1到RA7,且区段20B包含7个取样时段RB1到RB3。因此,参考信号A的区段20A可与目标激光信号的区段20的从时间T0到时间T7的取样时段匹配,且参考信号B的区段20B可与目标激光信号的区段20的从时间T7到时间T10的剩余取样时段匹配。因此,由于区段持续时间对应或匹配,所以可基于区段20A、20B的偏差而精确确定目标激光信号的时间区段20的偏差。
返回参考图1,与参考系统190A、190B相关联的延时可具有基于一或多个取样时段的差。具体来说,延时142A及延时142B可间隔达取样时段(例如,取样时间间隔)的整数倍。举例来说,如图2中所展示,区段20A可与区段20B相差4个取样时段。
在一些实施方案中,延时142A、142B中的一或多者可具有延迟时段,其为取样时段的质数倍。在一些实施方案中,可使用不具有公因子的延迟时段来界定延时142A、142B。
在一些实施方案中,就非零偏差(例如相位测量偏差)来说,两个参考系统190A、190B的延迟时段的较长者可具有基于比率(距离(距离)-LO(长度))/(长参考臂长度)的延时。距离是到物体5的往返延时,且LO长度是与延时142C关联的长度。(距离-LO)项可表示与源自于激光信号10的干涉计信号相关联的长度差。换句话说,距离项可为与激光信号10相关联的长度,其可包含到目标(例如物体5)的距离,且可为往返距离,且LO项可为与激光信号10的延迟版本相关联的长度。因此,(距离-LO)可表示源自于激光信号10的差拍及激光信号10的延迟版本的长度。在一些实施方案中,比率可小于10或应小于10。两个参考可用于有效联结(例如)相对较短时段的偏差(例如相位差)测量以获得相对较长时段的相位差估计。在一些实施方案中,可期望联结在一起的较少短时间测量用于精确长距离测量。
此外,在一些实施方案中,与两个参考系统190A、190B相关联的延迟时段的较短者可足够短,使得可精确估计相对较短延迟时段偏差,且可使用相对较少测量来估计给定(距离-LO)长度的偏差。在具有较短延时(或长度)的参考系统190A、190B中,较短延时(例如延时142A、延时142B)可足够大,使得归因于激光源110的缺陷,(例如)激光信号14B的均方根(RMS)相位测量误差(例如归因于噪声的误差)小于测量RMS相位偏差。
在一些实施方案中,偏差检测器150A、150B可经配置以确定(例如,计算、测量)与取样时段对应的时间区段的偏差。换句话说,可测量开始于一个以上取样时段的时间区段的偏差。因此,可测量具有各种持续时间且开始于不同时间(例如取样时间)的目标激光信号的时间区段的偏差。
在一些实施方案中,可基于与使用单个参考系统来测量的多个区段(例如联结区段、重叠区段)相关联的偏差而确定目标激光信号的区段的偏差。举例来说,偏差检测器可经配置以使用具有延时的参考系统来测量与开始于第一时间的第一时间区段相关联的第一偏差。所述偏差检测器可经配置以使用具有相同延时的相同参考系统来测量与开始于第二时间的第二时间区段(其不同于所述第一时段)相关联的第二偏差。在一些实施方案中,所述第一时间区段可与所述第二时间区段互斥。在一些实施方案中,所述第一时间区段可与所述第二时间区段重叠。特定来说,所述区段的重叠可发生于测量残余偏差时。结合(例如)至少图5来更详细地论述残余偏差。
图4说明图3中所展示的目标激光信号的时间区段的位移的另一实例。如图4中所展示,使时间区段22A从图3中所展示的时间区段21A位移一个取样时段。类似地,使时间区段22B从图3B中所展示的时间区段21B位移一个取样时段。因此,与时间区段22A及22B相关联的偏差可用于确定与时间区段22相关联的偏差,时间区段22从图3中所展示的时间区段21位移一个取样时段。图4中未展示与参考信号A或相关联的参考信号B相关联的额外时间区段。
返回参考图1,在一些实施方案中,可基于使用参考系统190A、190B的仅一者来检测的偏差而确定与目标激光信号的时间区段相关联的总偏差。举例来说,可组合使用与两个或两个以上时间区段(其与参考系统190A的参考信号14A相关联)相关联的两个或两个以上偏差来确定与所述目标激光信号的时间区段相关联的总偏差。
在一些实施方案中,与参考系统190A、190B相关联的时间区段无法以所要方式与目标激光系统的时间区段匹配。在此类例子中,包含于图1中所展示的激光系统100中的残余偏差计算器180可经配置以基于与参考系统190A、190B相关联的时间区段中的一或多者而使用各种方法(例如平均数、截断部分)来计算残余偏差。图5中说明此实例。
图5是说明根据实施方案的目标激光信号的另一实例的图式。如图5中所展示,所述目标激光信号具有时间S1与S12之间的目标时间区段52。与参考信号A的时间区段52A相关联的偏差可用于确定与时间S1与S8之间的目标时间区段52的第一部分相关联的偏差,且与参考信号B的时间区段52B1相关联的偏差可用于确定与时间S8与S11之间的目标时间区段52的第二部分。时间区段52A及时间区段52B1的组合使S11与S12之间的目标时间区段52中的剩余部分(或残余部分)留下。在此实施方案中,与参考信号B的时间区段52B2相关联的偏差可用于确定与时间S11与S12之间的目标时间区段52的残余部分相关联的偏差。与目标时间区段52的残余部分相关联的偏差可称为残余偏差。在一些实施方案中,所述偏差的分率(或部分)或与时间区段52B2相关联的偏差的每取样时间平均数可用于估计与时间S11到S12之间的目标时间区段52的残余部分相关联的偏差(或残余偏差)。在一些实施方案中,用于计算(或估计)与目标时间区段52的残余部分相关联的残余偏差的区段可以时间S11及S12为中心。举例来说,与参考信号B的时间S10与S13之间的时间区段相关联的偏差(或其部分)可用于估计与可以时间S11及S12为中心的目标时间区段52的残余部分相关联的残余偏差。在一些实施方案中,可使用残余偏差计算器(例如图1中所展示的残余偏差计算器180)来计算偏差。
使用适当构造的参考臂(例如图1中所展示的参考系统190A、190B),在参考臂相位测量偏差内,可将(距离-LO)延迟相位历时精确(例如,完美)估计到最接近取样长度。一般来说,子取样相位估计的偏差将相对较小,这是因为一个取样时段小于1/(激光线宽)且所述时段内的激光变化应较小。
下文是(距离-LO)相位历时估计的实例。
假设
Ref15具有15取样延时
Ref7具有7取样本延时
那么,就正整数延时来说,可将Ref15相位时间序列添加到Ref7相位时间序列以合成Ref22时间序列。如果测量在(距离-LO)延迟时间接近22个取样时段的距离处的目标,那么与目标返回信号共轭及差拍的合成Ref22将产生用于频谱处理的近乎完美频调。
下文阐述实施方案的细节:
激光相位信号:phi(t)
phi(t)=phi0+w0*t+alpha*t^2+V(t) (1),
其中t是时间,phi0是时间0处的相位,w0是时间0处的光频,alpha是界定调频斜率的常数,且V(t)表示激光频率扫描的偏差或非线性。
参考臂相位信号给出如下:
R(t,tau)=phi(t)-phi(t-tau) (2)
=w0*tau+(2*alpha*tau)*t-alpha*tau^2+v(t)-v(t-tau)
=[w0*tau-alpha*tau^2]+[2*alpha]*tau*t+v(t)-v(t-tau),
其中tau是干涉计中的两个长度之间的时间差或延时。
下文界定:
K(tau)=[w0*tau-alpha*tau^2]=恒定相位项
2*alpha=2*pi*(调频斜率)=2*pi*hzpm*c=H,
其中hzpm是路径失衡的每米赫兹,且c是光在传播介质中的速度。
调频斜率(hz/sec)=hzps=hzpm*c
R(t,tau)=K(tau)+H*tau*t+v(t)-v(t-tau) (3)
现考虑离散时间序列:
t(j)=tSamp*j,其中tSamp是取样时间。
还考虑参考臂切割以获得离散行进时间长度tSamp*k。
R(j,k)=R(j*tSamp,k*tsamp) (4)
=K(tau*tSamp)+H*k*tSamp^2*j+v(j)-v(j-k)
如果我们将外差复合到基频带以忽略常数,那么我们具有:
R0(j,k)=v(j)-v(j-k) (5)
如果我们具有对应于延迟指数k及m的两个实际参数,那么:
R0(j,k)+R0(j-k,m)
=v(j)-v(j-k)+v(j-k)-v(j-k-m)
=v(j)-v(j-(k+m))=R0(j,k+m) (6)
我们已使用R0(j,k)及R0(j,m)来正确合成R0(j,k+m)。
类似地,
R0(j,k)-R0(j-k+m,m)
=v(j)-v(j-k)-v(j-(k-m))+v(j-(k-m)-m)
=v(j)-v(j-(k-m))=R(j,k-m) (7)
定理:
给定第三整数n,对于互质(不含公因子)的任何两个非零整数k及m,存在整数因子A及B,使得:
n=A*k+B*m (8)
因此,总的来说,我们可从两个整数取样长度参考臂(其整数是互质的)合成任意整数取样长度的参考臂信号。唯一偏差可为测量噪声。
近似参考区段估计:因为无法忽略测量噪声,所以粗略估计一些参考长度区段分量(非使用涉及大整数因子A及/或B的准确实施方案)有时是有意义的。举例来说,可使用15样本参考臂来粗略估计13个样本的参考长度区段,且选择样本指数,使得13个样本位于15个样本的中心。将13样本相位差估计为(13/15)乘以15样本测量相位差。类似地,可将17样本参考区段相位差估计为(17/15)乘以具有相同中心样本区段的测量15样本参考的相位差。
可展示:对于长度k及m的参考对,准确合成小于k*m的整数样本参考长度所需的区段的最大数(方程式(8)中的abs(A)及abs(b)的总和)是:
maxSegments=[(k-1)+(m-1)], (9)
其为大于期望值的数目。通过允许选择区段长度[(k-2),k,(k+2),(m-2),m,(m+2)]而大幅减小参考区段相位测量的最大数。举例来说,为合成1个到100个样本的任何整数样本参考长度,使用长度[5,7,9,13,15,17]的区段,最大数是7。对于1个到50个样本,最大数是4。实际上,我们将具有模型化及经验性参考相位测量偏差及近似长度测量偏差估计,且相应地选择准确且近似的区段来最小化偏差。在所有情形中,最小偏差将比通过已知HS处理而获得的偏差更好。
如果系统用于定期测量长距离,那么较长参考臂及LO是恰当的。
这看似复杂的合成参考臂估计过程实际上很容易实施且产生明显优于现有HS处理实施方案的性能优势。除一致SNR改进之外,还应减小速度及距离估计的渐增式线性调频与渐减式线性调频的差。
下文是实例性1到100样本长度参考的结果。在此编纂中,准确区段并不比近似区段优先。一旦良好估计可用于准确区段及近似区段的偏差,此就将被添加。
****两个Refs模型****
nSegRef1=15 nSegRef2=7
sigmaP=0.025 sigmaM=0.050
nSegRange=1
1=-6*15+13*7 totSegs=19 ErrMeas=0.218
1=1*15+-2*7 totSegs=3 ErrMeas=0.087
nSegRange=2
2=-5*15+11*7 totSegs=16 ErrMeas=0.200
2=2*15+-4*7 totSegs=6 ErrMeas=0.122
2=-1*15+0*7+0*13+1*17+0*9+0*5GSegs=1 ASegs=1 totSegs=2
nSegRange=3
3=-4*15+9*7 totSegs=13 ErrMeas=0.180
3=3*15+-6*7 totSegs=9 ErrMeas=0.150
3=-1*15+0*7+0*13+0*17+2*9+0*5 GSegs=1 ASegs=2 totSegs=3
nSegRange=4
4=-3*15+7*7 totSegs=10 ErrMeas=0.158
4=4*15+-8*7 totSegs=12 ErrMeas=0.173
4=0*15+0*7+-1*13+1*17+0*9+0*5 GSegs=0 ASegs=2 totSegs=2
nSegRange=5
5=-2*15+5*7 totSegs=7 ErrMeas=0.132
5=5*15+-10*7 totSegs=15 ErrMeas=0.194
5=0*15+0*7+0*13+0*17+0*9+1*5 GSegs=0 ASegs=1 totSegs=1
nSegRange=6
6=-1*15+3*7 totSegs=4 ErrMeas=0.100
6=6*15+-12*7 totSegs=18 ErrMeas=0.212
6=0*15+-1*7+1*13+0*17+0*9+0*5 GSegs=1 ASegs=1 totSegs=2
nSegRange=7
7=0*15+1*7 totSegs=1 ErrMeas=0.050
nSegRange=8
8=-6*15+14*7 totSegs=20 ErrMeas=0.224
8=1*15+-1*7 totSegs=2 ErrMeas=0.071
nSegRange=9
9=-5*15+12*7 totSegs=17 ErrMeas=0.206
9=2*15+-3*7 totSegs=5 ErrMeas=0.112
9=0*15+0*7+0*13+0*17+1*9+0*5 GSegs=0 ASegs=1 totSegs=1
nSegRange=10
10=-4*15+10*7 totSegs=14 ErrMeas=0.187
10=3*15+-5*7 totSegs=8 ErrMeas=0.141
10=0*15+-1*7+0*13+1*17+0*9+0*5 GSegs=1 ASegs=1 totSegs=2
nSegRange=11
11=-3*15+8*7 totSegs=11 ErrMeas=0.166
11=4*15+-7*7 totSegs=11 ErrMeas=0.166
11=-1*15+0*7+0*13+1*17+1*9+0*5 GSegs=1 ASegs=2 totSegs=3
nSegRange=12
12=-2*15+6*7 totSegs=8 ErrMeas=0.141
12=5*15+-9*7 totSegs=14 ErrMeas=0.187
12=0*15+0*7+0*13+1*17+0*9+-1*5 GSegs=0 ASegs=2 totSegs=2
nSegRange=13
13=-1*15+4*7 totSegs=5 ErrMeas=0.112
13=6*15+-11*7 totSegs=17 ErrMeas=0.206
13=0*15+0*7+1*13+0*17+0*9+0*5 GSegs=0 ASegs=1 totSegs=1
nSegRange=14
14=0*15+2*7 totSegs=2 ErrMeas=0.071
nSegRange=15
15=1*15+0*7 totSegs=1 ErrMeas=0.050
nSegRange=16
16=-5*15+13*7 totSegs=18 ErrMeas=0.212
16=2*15+-2*7 totSegs=4 ErrMeas=0.100
16=0*15+1*7+0*13+0*17+1*9+0*5 GSegs=1 ASegs=1 totSegs=2
nSegRange=17
17=-4*15+11*7 totSegs=15 ErrMeas=0.194
17=3*15+-4*7 totSegs=7 ErrMeas=0.132
17=0*15+0*7+0*13+1*17+0*9+0*5 GSegs=0 ASegs=1 totSegs=1
nSegRange=18
18=-3*15+9*7totSegs=12 ErrMeas=0.173
18=4*15+-6*7 totSegs=10 ErrMeas=0.158
18=0*15+0*7+0*13+0*17+2*9+0*5 GSegs=0 ASegs=2 totSegs=2
nSegRange=19
19=-2*15+7*7 totSegs=9 ErrMeas=0.150
19=5*15+-8*7 totSegs=13 ErrMeas=0.180
19=-1*15+0*7+0*13+2*17+0*9+0*5 GSegs=1 ASegs=2 totSegs=3
nSegRange=20
20=-1*15+5*7 totSegs=6 ErrMeas=0.122
20=6*15+-10*7 totSegs=16 ErrMeas=0.200
20=0*15+1*7+1*13+0*17+0*9+0*5 GSegs=1 ASegs=1 totSegs=2
nSegRange=21
21=0*15+3*7 totSegs=3 ErrMeas=0.087
nSegRange=22
22=1*15+1*7 totSegs=2 ErrMeas=0.071
nSegRange=23
23=-5*15+14*7 totSegs=19 ErrMeas=0.218
23=2*15+-1*7 totSegs=3 ErrMeas=0.087
nSegRange=24
24=-4*15+12*7 totSegs=16 ErrMeas=0.200
24=3*15+-3*7 totSegs=6 ErrMeas=0.122
24=0*15+1*7+0*13+1*17+0*9+0*5 GSegs=1 ASegs=1 totSegs=2
nSegRange=25
25=-3*15+10*7 totSegs=13 ErrMeas=0.180
25=4*15+-5*7 totSegs=9 ErrMeas=0.150
25=0*15+0*7+0*13+2*17+-1*9+0*5 GSegs=0 ASegs=3 totSegs=3
nSegRange=26
26=-2*15+8*7 totSegs=10 ErrMeas=0.158
26=5*15+-7*7 totSegs=12 ErrMeas=0.173
26=0*15+0*7+0*13+1*17+1*9+0*5 GSegs=0 ASegs=2 totSegs=2
nSegRange=27
27=-1*15+6*7 totSegs=7 ErrMeas=0.132
27=6*15+-9*7 totSegs=15 ErrMeas=0.194
27=0*15+-1*7+0*13+2*17+0*9+0*5 GSegs=1 ASegs=2 totSegs=3
nSegRange=28
28=0*15+4*7 totSegs=4 ErrMeas=0.100
28=1*15+0*7+1*13+0*17+0*9+0*5 GSegs=1 ASegs=1 totSegs=2
nSegRange=29
29=1*15+2*7 totSegs=3 ErrMeas=0.087
nSegRange=30
30=2*15+0*7 totSegs=2 ErrMeas=0.071
nSegRange=31
31=-4*15+13*7 totSegs=17 ErrMeas=0.206
31=3*15+-2*7 totSegs=5 ErrMeas=0.112
31=0*15+0*7+0*13+1*17+1*9+1*5 GSegs=0 ASegs=3 totSegs=3
nSegRange=32
32=-3*15+11*7 totSegs=14 ErrMeas=0.187
32=4*15+-4*7 totSegs=8 ErrMeas=0.141
32=1*15+0*7+0*13+1*17+0*9+0*5 GSegs=1 ASegs=1 totSegs=2
nSegRange=33
33=-2*15+9*7 totSegs=11 ErrMeas=0.166
33=5*15+-6*7 totSegs=11 ErrMeas=0.166
33=0*15+1*7+0*13+1*17+1*9+0*5 GSegs=1 ASegs=2 totSegs=3
nSegRange=34
34=-1*15+7*7 totSegs=8 ErrMeas=0.141
34=6*15+-8*7 totSegs=14 ErrMeas=0.187
34=0*15+0*7+0*13+2*17+0*9+0*5 GSegs=0 ASegs=2 totSegs=2
nSegRange=35
35=0*15+5*7 totSegs=5 ErrMeas=0.112
35=0*15+0*7+0*13+1*17+2*9+0*5 GSegs=0 ASegs=3 totSegs=3
nSegRange=36
36=1*15+3*7 totSegs=4 ErrMeas=0.100
nSegRange=37
37=2*15+1*7 totSegs=3 ErrMeas=0.087
nSegRange=38
38=-4*15+14*7 totSegs=18 ErrMeas=0.212
38=3*15+-1*7 totSegs=4 ErrMeas=0.100
nSegRange=39
39=-3*15+12*7 totSegs=15 ErrMeas=0.194
39=4*15+-3*7 totSegs=7 ErrMeas=0.132
39=0*15+0*7+0*13+2*17+0*9+1*5 GSegs=0 ASegs=3 totSegs=3
nSegRange=40
40=-2*15+10*7 totSegs=12 ErrMeas=0.173
40=5*15+-5*7 totSegs=10 ErrMeas=0.158
40=0*15+-1*7+1*13+2*17+0*9+0*5 GSegs=1 ASegs=3 totSegs=4
nSegRange=41
41=-1*15+8*7 totSegs=9 ErrMeas=0.150
41=6*15+-7*7 totSegs=13 ErrMeas=0.180
41=0*15+1*7+0*13+2*17+0*9+0*5 GSegs=1 ASegs=2 totSegs=3
nSegRange=42
42=0*15+6*7 totSegs=6 ErrMeas=0.122
42=0*15+0*7+0*13+3*17+-1*9+0*5 GSegs=0 ASegs=4 totSegs=4
nSegRange=43
43=1*15+4*7 totSegs=5 ErrMeas=0.112
43=0*15+0*7+0*13+2*17+1*9+0*5 GSegs=0 ASegs=3 totSegs=3
nSegRange=44
44=2*15+2*7 totSegs=4 ErrMeas=0.100
nSegRange=45
45=3*15+0*7 totSegs=3 ErrMeas=0.087
nSegRange=46
46=-3*15+13*7 totSegs=16 ErrMeas=0.200
46=4*15+-2*7 totSegs=6 ErrMeas=0.122
46=0*15+0*7+0*13+3*17+0*9+-1*5 GSegs=0 ASegs=4 totSegs=4
nSegRange=47
47=-2*15+11*7 totSegs=13 ErrMeas=0.180
47=5*15+-4*7 totSegs=9 ErrMeas=0.150
47=0*15+0*7+1*13+2*17+0*9+0*5 GSegs=0 ASegs=3 totSegs=3
nSegRange=48
48=-1*15+9*7 totSegs=10 ErrMeas=0.158
48=6*15+-6*7 totSegs=12 ErrMeas=0.173
48=0*15+0*7+0*13+2*17+1*9+1*5 GSegs=0 ASegs=4 totSegs=4
nSegRange=49
49=0*15+7*7 totSegs=7 ErrMeas=0.132
49=1*15+0*7+0*13+2*17+0*9+0*5 GSegs=1 ASegs=2 totSegs=3
nSegRange=50
50=1*15+5*7 totSegs=6 ErrMeas=0.122
50=0*15+1*7+0*13+2*17+1*9+0*5 GSegs=1 ASegs=3 totSegs=4
nSegRange=51
51=2*15+3*7 totSegs=5 ErrMeas=0.112
51=0*15+0*7+0*13+3*17+0*9+0*5 GSegs=0 ASegs=3 totSegs=3
nSegRange=52
52=3*15+1*7 totSegs=4 ErrMeas=0.100
nSegRange=53
53=-3*15+14*7 totSegs=17 ErrMeas=0.206
53=4*15+-1*7 totSegs=5 ErrMeas=0.112
nSegRange=54
54=-2*15+12*7 totSegs=14 ErrMeas=0.187
54=5*15+-3*7 totSegs=8 ErrMeas=0.141
54=0*15+1*7+1*13+2*17+0*9+0*5 GSegs=1 ASegs=3 totSegs=4
nSegRange=55
55=-1*15+10*7 totSegs=11 ErrMeas=0.166
55=6*15+-5*7 totSegs=11 ErrMeas=0.166
55=0*15+0*7+-1*13+4*17+0*9+0*5 GSegs=0 ASegs=5 totSegs=5
nSegRange=56
56=0*15+8*7 totSegs=8 ErrMeas=0.141
56=0*15+0*7+0*13+3*17+0*9+1*5 GSegs=0 ASegs=4 totSegs=4
nSegRange=57
57=1*15+6*7 totSegs=7 ErrMeas=0.132
57=0*15+-1*7+1*13+3*17+0*9+0*5 GSegs=1 ASegs=4 totSegs=5
nSegRange=58
58=2*15+4*7 totSegs=6 ErrMeas=0.122
58=0*15+1*7+0*13+3*17+0*9+0*5 GSegs=1 ASegs=3 totSegs=4
nSegRange=59
59=3*15+2*7 totSegs=5 ErrMeas=0.112
nSegRange=60
60=4*15+0*7 totSegs=4 ErrMeas=0.100
nSegRange=61
61=-2*15+13*7 totSegs=15 ErrMeas=0.194
61=5*15+-2*7 totSegs=7 ErrMeas=0.132
61=0*15+-1*7+0*13+4*17+0*9+0*5 GSegs=1 ASegs=4 totSegs=5
nSegRange=62
62=-1*15+11*7 totSegs=12 ErrMeas=0.173
62=6*15+-4*7 totSegs=10 ErrMeas=0.158
62=1*15+0*7+1*13+2*17+0*9+0*5 GSegs=1 ASegs=3 totSegs=4
nSegRange=63
63=0*15+9*7 totSegs=9 ErrMeas=0.150
63=0*15+0*7+0*13+4*17+0*9+-1*5 GSegs=0 ASegs=5 totSegs=5
nSegRange=64
64=1*15+7*7 totSegs=8 ErrMeas=0.141
64=0*15+0*7+1*13+3*17+0*9+0*5 GSegs=0 ASegs=4 totSegs=4
nSegRange=65
65=2*15+5*7 totSegs=7 ErrMeas=0.132
65=0*15+0*7+0*13+3*17+1*9+1*5 GSegs=0 ASegs=5 totSegs=5
nSegRange=66
66=3*15+3*7 totSegs=6 ErrMeas=0.122
66=1*15+0*7+0*13+3*17+0*9+0*5 GSegs=1 ASegs=3 totSegs=4
nSegRange=67
67=4*15+1*7 totSegs=5 ErrMeas=0.112
nSegRange=68
68=-2*15+14*7 totSegs=16 ErrMeas=0.200
68=5*15+-1*7 totSegs=6 ErrMeas=0.122
68=0*15+0*7+0*13+4*17+0*9+0*5 GSegs=0 ASegs=4 totSegs=4
nSegRange=69
69=-1*15+12*7 totSegs=13 ErrMeas=0.180
69=6*15+-3*7 totSegs=9 ErrMeas=0.150
69=0*15+0*7+0*13+3*17+2*9+0*5 GSegs=0 ASegs=5 totSegs=5
nSegRange=70
70=0*15+10*7 totSegs=10 ErrMeas=0.158
70=-1*15+0*7+0*13+5*17+0*9+0*5 GSegs=1 ASegs=5 totSegs=6
nSegRange=71
71=1*15+8*7 totSegs=9 ErrMeas=0.150
71=0*15+1*7+1*13+3*17+0*9+0*5 GSegs=1 ASegs=4 totSegs=5
nSegRange=72
72=2*15+6*7 totSegs=8 ErrMeas=0.141
72=0*15+0*7+-1*13+5*17+0*9+0*5 GSegs=0 ASegs=6 totSegs=6
nSegRange=73
73=3*15+4*7 totSegs=7 ErrMeas=0.132
73=0*15+0*7+0*13+4*17+0*9+1*5 GSegs=0 ASegs=5 totSegs=5
nSegRange=74
74=4*15+2*7 totSegs=6 ErrMeas=0.122
nSegRange=75
75=5*15+0*7 totSegs=5 ErrMeas=0.112
nSegRange=76
76=-1*15+13*7 totSegs=14 ErrMeas=0.187
76=6*15+-2*7 totSegs=8 ErrMeas=0.141
76=0*15+0*7+0*13+5*17+-1*9+0*5 GSegs=0 ASegs=6 totSegs=6
nSegRange=77
77=0*15+11*7 totSegs=11 ErrMeas=0.166
77=0*15+0*7+0*13+4*17+1*9+0*5 GSegs=0 ASegs=5 totSegs=5
nSegRange=78
78=1*15+9*7 totSegs=10 ErrMeas=0.158
78=0*15+-1*7+0*13+5*17+0*9+0*5 GSegs=1 ASegs=5 totSegs=6
nSegRange=79
79=2*15+7*7 totSegs=9 ErrMeas=0.150
79=1*15+0*7+1*13+3*17+0*9+0*5 GSegs=1 ASegs=4 totSegs=5
nSegRange=80
80=3*15+5*7 totSegs=8 ErrMeas=0.141
80=0*15+0*7+0*13+5*17+0*9+-1*5 GSegs=0 ASegs=6 totSegs=6
nSegRange=81
81=4*15+3*7 totSegs=7 ErrMeas=0.132
81=0*15+0*7+1*13+4*17+0*9+0*5 GSegs=0 ASegs=5 totSegs=5
nSegRange=82
82=5*15+1*7 totSegs=6 ErrMeas=0.122
nSegRange=83
83=-1*15+14*7 totSegs=15 ErrMeas=0.194
83=6*15+-1*7 totSegs=7 ErrMeas=0.132
83=1*15+0*7+0*13+4*17+0*9+0*5 GSegs=1 ASegs=4 totSegs=5
nSegRange=84
84=0*15+12*7 totSegs=12 ErrMeas=0.173
84=0*15+1*7+0*13+4*17+1*9+0*5 GSegs=1 ASegs=5 totSegs=6
nSegRange=85
85=1*15+10*7 totSegs=11 ErrMeas=0.166
85=0*15+0*7+0*13+5*17+0*9+0*5 GSegs=0 ASegs=5 totSegs=5
nSegRange=86
86=2*15+8*7 totSegs=10 ErrMeas=0.158
86=0*15+0*7+0*13+4*17+2*9+0*5 GSegs=0 ASegs=6 totSegs=6
nSegRange=87
87=3*15+6*7 totSegs=9 ErrMeas=0.150
87=-1*15+0*7+0*13+6*17+0*9+0*5 GSegs=1 ASegs=6 totSegs=7
nSegRange=88
88=4*15+4*7 totSegs=8 ErrMeas=0.141
88=0*15+1*7+1*13+4*17+0*9+0*5 GSegs=1 ASegs=5 totSegs=6
nSegRange=89
89=5*15+2*7 totSegs=7 ErrMeas=0.132
nSegRange=90
90=6*15+0*7 totSegs=6 ErrMeas=0.122
nSegRange=91
91=0*15+13*7 totSegs=13 ErrMeas=0.180
91=0*15+-1*7+1*13+5*17+0*9+0*5 GSegs=1 ASegs=6 totSegs=7
nSegRange=92
92=1*15+11*7 totSegs=12 ErrMeas=0.173
92=0*15+1*7+0*13+5*17+0*9+0*5 GSegs=1 ASegs=5 totSegs=6
nSegRange=93
93=2*15+9*7 totSegs=11 ErrMeas=0.166
93=0*15+0*7+0*13+6*17+-1*9+0*5 GSegs=0 ASegs=7 totSegs=7
nSegRange=94
94=3*15+7*7 totSegs=10 ErrMeas=0.158
94=0*15+0*7+0*13+5*17+1*9+0*5GSegs=0 ASegs=6 totSegs=6
nSegRange=95
95=4*15+5*7 totSegs=9 ErrMeas=0.150
95=0*15+-1*7+0*13+6*17+0*9+0*5GSegs=1 ASegs=6 totSegs=7
nSegRange=96
96=5*15+3*7 totSegs=8 ErrMeas=0.141
96=1*15+0*7+1*13+4*17+0*9+0*5GSegs=1 ASegs=5 totSegs=6
nSegRange=97
97=6*15+1*7 totSegs=7 ErrMeas=0.132
nSegRange=98
98=0*15+14*7 totSegs=14 ErrMeas=0.187
98=0*15+0*7+1*13+5*17+0*9+0*5GSegs=0 ASegs=6 totSegs=6
nSegRange=99
99=1*15+12*7 totSegs=13 ErrMeas=0.180
99=0*15+0*7+0*13+5*17+1*9+1*5GSegs=0 ASegs=7 totSegs=7
nSegRange=100
100=2*15+10*7 totSegs=12 ErrMeas=0.173
100=1*15+0*7+0*13+5*17+0*9+0*5GSegs=1 ASegs=5 totSegs=6
在一些实施方案中,图1的激光系统100中所展示的组件的一或多个部分可为或可包含基于硬件的模块(例如数字信号处理器(DSP)、现场可编程门阵列(FPGA)、存储器)、固件模块及/或基于软件的模块(例如计算机代码的模块、可在计算机处执行的一组计算机可读指令)。举例来说,在一些实施方案中,激光系统100的一或多个部分可为或可包含经配置以由至少一个处理器(未展示)执行的软件模块。在一些实施方案中,组件的功能性可包含于不同于图1中所展示的模块及/或组件的模块及/或组件中。
在一些实施例中,激光系统100的组件中的一或多者可为或可包含经配置以处理存储于存储器中的指令的处理器。举例来说,分析器170(及/或其部分)可为处理器及存储器的组合,其经配置以执行与过程相关的指令以实施一或多个功能。
尽管未展示,但在一些实施方案中,激光系统100的组件(或激光系统100的部分)可经配置以在(例如)数据中心(例如云端运算环境)、计算机系统、一或多个服务器/主机装置等等内操作。在一些实施方案中,激光系统100的组件(或激光系统100的部分)可经配置以在网络内操作。因此,激光系统100(或其部分)可经配置以在可包含一或多个装置及/或一或多个服务器装置的各种类型的网络环境内运行。举例来说,网络可为或可包含局域网(LAN)、广域网(WAN)等等。网络可为或可包含无线网络及/或使用(例如)网关件、网桥、交换机等等来实施的无线网络。网络可包含一或多个区段及/或可具有基于各种协议(例如因特网协议(IP)及/或专属协议)的部分。网络可包含因特网的至少一部分。
在一些实施方案中,存储器可为任何类型的存储器,例如随机存取存储器、磁盘驱动存储器、快闪存储器及/或等等。在一些实施方案中,可将存储器实施为与激光系统100的组件相关联的一个以上存储器组件(例如一个以上RAM组件或磁盘驱动存储器)。
可在数字电子电路中或在计算机硬件、固件、软件或其组合中实施本文所描述的各种技术的实施方案。可将实施方案实施为计算机程序产品,即,在信息载体中(例如,在机器可读存储装置(计算机可读媒体、非暂时性计算机可读存储媒体、有形计算机可读存储媒体)中)或在传播信号中有形地体现的计算机程序,以用于由数据处理设备(例如可编程处理器、计算机或多个计算机)处理或用于控制数据处理设备的操作。计算机程序(例如上文所描述的计算机程序)可以任何形式的编程语言(其包含编译或解译语言)写入,且可以任何形式部署,其包含作为独立程序或作为适合用于运算环境中的模块、组件、子例程或其它单元。计算机程序可经部署以在一个计算机上或在一个地点处或分布于多个地点且由通信网络互连的多个计算机上被处理。
可由执行计算机程序的一或多个可编程处理器执行方法步骤以通过操作输入数据且产生输出而执行功能。还可由专用逻辑电路(例如FPGA(现场可编程门阵列)或ASIC(专用集成电路))执行方法步骤,且可将设备实施为所述专用逻辑电路。
适合用于处理计算机程序的处理器包含(例如)通用微处理器及专用微处理器两者、及任何类型的数字计算机的任何一或多个处理器。一般来说,处理器将从只读存储器或随机存取存储器或两者接收指令及数据。计算机的元件可包含用于执行指令的至少一个处理器及用于存储指令及数据的一或多个存储器装置。一般来说,计算机还可包含用于存储数据的一或多个大容量存储装置(例如磁盘、磁光盘或光盘),或经操作性耦合以从所述一或多个大容量存储装置接收数据或将数据传送到所述一或多个大容量存储装置,或以上两者。适合用于体现计算机程序指令及数据的信息载体包含所有形式的非易失性存储器,其包含(例如):半导体存储器装置,例如EPROM、EEPROM及快闪存储器装置;磁盘,例如内部硬盘或可抽换式磁盘;磁光盘;及CD-ROM磁盘及DVD-ROM磁盘。处理器及存储器可补充专用逻辑电路或并入专用逻辑电路中。
为提供与用户的交互,可在具有显示装置(例如液晶显示(LCD)监视器)(其用于将信息显示给用户)及键盘及指向装置(例如鼠标或轨迹球)(用户可通过所述键盘及所述指向装置而将输入提供到计算机)的计算机上实施实施方案。其它类型的装置还可用于提供与用户的交互;举例来说,提供到用户的反馈可为任何形式的感测反馈,例如视觉反馈、听觉反馈或触觉反馈;且可从用户接收任何形式的输入,其包含声音输入、语音输入或触觉输入。
可在运算系统中实施实施方案,所述运算系统包含后端组件(例如)作为数据服务器,或包含中间件组件(例如应用程序服务器),或包含前端组件(例如客户端计算机,其具有图形用户接口或网页浏览器,用户可通过所述图形用户接口或所述网页浏览器而与实施方案交互),或包含此类后端组件、中间件组件或前端组件的任何组合。组件可由任何形式或介质的数字数据通信(例如通信网络)互连。通信网络的实例包含局域网(LAN)及广域网(WAN),例如因特网。
尽管已如本文所描述那样说明所描述的实施方案的某些特征,但所属领域的技术人员现将设想许多修改、替换、改变及等效物。因此,应了解,所附权利要求书希望涵盖落入实施方案的范围内的所有此类修改及改变。应了解,所述修改及改变仅供例示,而非限制,且可对形式及细节作出各种改变。可以除互斥组合之外的任何组合方式组合本文所描述的设备及/或方法的任何部分。本文所描述的实施方案可包含所描述的不同实施方案的功能、组件及/或特征的各种组合及/或子组合。
Claims (15)
1.一种非暂时性计算机可读存储媒体,其存储在被执行时致使处理器执行过程的指令,所述过程包括:
产生激光信号的区段,所述激光信号的所述区段具有持续时间;
基于所述激光信号而产生第一参考信号;
基于所述第一参考信号而计算与所述持续时间的第一部分对应的第一相位偏差;
基于所述激光信号而产生第二参考信号;
基于所述第二参考信号而计算与所述持续时间的第二部分对应的第二相位偏差;
及
基于所述第一相位偏差及所述第二相位偏差的组合而计算所述激光信号的所述区段的相位偏差。
2.根据权利要求1所述的非暂时性计算机可读存储媒体,其中所述持续时间与所述激光信号到目标物体的往返时段对应。
3.根据权利要求1所述的非暂时性计算机可读存储媒体,其中所述持续时间的所述第一部分及所述持续时间的所述第二部分发生于互斥时段期间。
4.根据权利要求1所述的非暂时性计算机可读存储媒体,其中所述组合是和。
5.根据权利要求1所述的非暂时性计算机可读存储媒体,其中所述组合是差。
6.一种方法,其包括:
产生激光信号;
基于取样时段而分析所述激光信号的反射信号;
基于所述激光信号及第一延迟持续时间而产生第一参考信号,所述第一延迟持续时间是所述取样时段的整数倍;
基于所述第一参考信号而计算与所述第一延迟持续时间相关联的相位偏差;及
基于所述激光信号及不同于所述第一延迟持续时间的第二延迟持续时间而产生第二参考信号。
7.根据权利要求6所述的方法,其中所述相位偏差是第一相位偏差,
所述方法进一步包括:
基于所述第二参考信号而计算与所述第二延迟持续时间相关联的第二相位偏差;及
使用所述第一相位偏差及所述第二相位偏差来计算所述激光信号的往返相位偏差。
8.根据权利要求6所述的方法,其中所述相位偏差是第一相位偏差,
所述方法进一步包括:
基于所述第二参考信号而计算与所述第二延迟持续时间相关联的第二相位偏差;及
基于所述第一相位偏差及所述第二相位偏差的总和而计算所述激光信号的相位偏差,所述相位偏差与等于或大于所述第一延迟持续时间及所述第二延迟持续时间的和的往返持续时间相关联。
9.根据权利要求6所述的方法,其中所述相位偏差是第一相位偏差,
所述方法进一步包括:
基于所述第二参考信号而计算与所述第二延迟持续时间相关联的第二相位偏差;及
基于所述第一相位偏差与所述第二相位偏差之间的差而计算所述激光信号的相位偏差,所述相位偏差与小于所述第一延迟持续时间及所述第二延迟持续时间的和的往返持续时间相关联。
10.根据权利要求6所述的方法,其进一步包括:
基于到目标物体的距离而确定与所述激光信号相关联的往返持续时间;及
估计与所述往返持续时间与所述第一延迟持续时间及所述第二延迟持续时间的和之间的差对应的所述激光信号的部分的残余相位偏差。
11.根据权利要求6所述的方法,其中所述取样时段的所述整数倍是所述取样时段的第一整数倍,所述第二延迟持续时间是所述取样时段的第二整数倍。
12.根据权利要求11所述的方法,其中所述第一整数倍及所述第二整数倍不具有公因子。
13.根据权利要求11所述的方法,其中所述第一整数倍是质数,所述第二整数倍是质数。
14.一种设备,其包括:
激光源,其经配置以产生激光信号;
接收器,其经配置以基于取样时段而分析来自所述激光信号的反射信号;
第一参考模块,其包含第一延时且经配置以基于所述激光信号而产生第一参考信号;
相位检测器,其经配置以检测与由所述第一延时产生的持续时间相关联的相位偏差;
第二参考模块,其包含第二延时且经配置以基于所述激光信号而产生第二参考信号;及
相位检测器,其经配置以检测与由所述第二延时产生的持续时间相关联的相位偏差,由所述第一延时产生的所述持续时间与由所述第二延时产生的所述持续时间相差所述取样时段的整数倍。
15.根据权利要求14所述的设备,其中所述第一延时的所述持续时间是所述取样时段的整数倍。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911253119.9A CN111025263A (zh) | 2014-07-25 | 2015-07-24 | 激光相位估计及校正 |
Applications Claiming Priority (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201462028911P | 2014-07-25 | 2014-07-25 | |
US62/028,911 | 2014-07-25 | ||
US14/807,508 | 2015-07-23 | ||
US14/807,508 US10151827B2 (en) | 2014-07-25 | 2015-07-23 | Laser phase estimation and correction |
PCT/US2015/041898 WO2016014894A1 (en) | 2014-07-25 | 2015-07-24 | Laser phase estimation and correction |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911253119.9A Division CN111025263A (zh) | 2014-07-25 | 2015-07-24 | 激光相位估计及校正 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106574965A true CN106574965A (zh) | 2017-04-19 |
CN106574965B CN106574965B (zh) | 2020-01-03 |
Family
ID=53785744
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911253119.9A Pending CN111025263A (zh) | 2014-07-25 | 2015-07-24 | 激光相位估计及校正 |
CN201580045019.0A Expired - Fee Related CN106574965B (zh) | 2014-07-25 | 2015-07-24 | 激光相位估计及校正 |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911253119.9A Pending CN111025263A (zh) | 2014-07-25 | 2015-07-24 | 激光相位估计及校正 |
Country Status (6)
Country | Link |
---|---|
US (2) | US10151827B2 (zh) |
EP (2) | EP3444630B1 (zh) |
JP (1) | JP6653695B2 (zh) |
CN (2) | CN111025263A (zh) |
TW (2) | TWI648550B (zh) |
WO (1) | WO2016014894A1 (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111025320A (zh) * | 2019-12-28 | 2020-04-17 | 深圳奥锐达科技有限公司 | 一种相位式激光测距系统及测距方法 |
CN112099036A (zh) * | 2020-11-10 | 2020-12-18 | 深圳市汇顶科技股份有限公司 | 距离测量方法以及电子设备 |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10151827B2 (en) * | 2014-07-25 | 2018-12-11 | DSCG Solutions, Inc. | Laser phase estimation and correction |
US9877114B2 (en) * | 2015-04-13 | 2018-01-23 | DSCG Solutions, Inc. | Audio detection system and methods |
US10690756B2 (en) | 2016-05-10 | 2020-06-23 | Texas Instruments Incorporated | Methods and apparatus for LIDAR operation with pulse position modulation |
US10613204B2 (en) | 2016-05-10 | 2020-04-07 | Texas Instruments Incorporated | Methods and apparatus for lidar operation with sequencing of pulses |
US10908287B2 (en) * | 2016-05-10 | 2021-02-02 | Texas Instruments Incorporated | Methods and apparatus for LIDAR operation with narrowband intensity modulation |
US11221411B2 (en) | 2016-07-18 | 2022-01-11 | Texas Instruments Incorporated | Power efficient LIDAR |
TWI647470B (zh) * | 2017-09-21 | 2019-01-11 | 光寶科技股份有限公司 | 光達裝置 |
WO2019145843A1 (en) | 2018-01-23 | 2019-08-01 | Amo Wavefront Sciences, Llc | Methods and systems of optical coherence tomography with fiducial signal for correcting scanning laser nonlinearity |
US11150349B2 (en) * | 2018-08-16 | 2021-10-19 | Wei Chen | Multi-line, high-definition LiDAR device and method with integrated direct spatial reference |
DE102019135648A1 (de) * | 2019-12-21 | 2021-06-24 | Carl Zeiss Ag | Vorrichtung und Verfahren zur Abstandsermittlung eines Objekts |
US11754711B2 (en) * | 2019-12-31 | 2023-09-12 | Luminar Technologies, Inc. | Frequency chirp for lidar for high-velocity targets |
US11169246B1 (en) | 2020-10-05 | 2021-11-09 | Aeva, Inc. | Techniques to compensate for phase impairments in LIDAR systems |
CN115006826B (zh) * | 2022-08-09 | 2022-10-25 | 澳瑞特体育产业股份有限公司 | 一种运动行程采集装置、方法及气阻健身训练设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030021309A1 (en) * | 2001-06-19 | 2003-01-30 | Mattox Barry G. | Method for producing highly accurate frequency and FM of a laser |
US20080088503A1 (en) * | 2004-09-28 | 2008-04-17 | Beasley Patrick D L | Frequency Modulated Continuous Wave (Fmcw) Radar Having Improved Frequency Linearity |
CN101266293A (zh) * | 2008-04-30 | 2008-09-17 | 西安电子科技大学 | 激光合成孔径雷达成像距离向相位补偿方法 |
CN102419441A (zh) * | 2011-09-01 | 2012-04-18 | 哈尔滨工业大学 | 一种基于四通道探测技术的弱光锁相星间位移测量方法及实现该方法的装置 |
Family Cites Families (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US1219583A (en) | 1914-07-13 | 1917-03-20 | Ward Sumner Perry | Headlight-dimmer. |
US4637723A (en) * | 1982-09-30 | 1987-01-20 | Honeywell Inc. | Discriminant apparatus for ring laser angular rate sensors |
US4830486A (en) * | 1984-03-16 | 1989-05-16 | Goodwin Frank E | Frequency modulated lasar radar |
JPH01219583A (ja) * | 1988-02-26 | 1989-09-01 | Daikin Ind Ltd | 距離測定装置 |
JP2954871B2 (ja) * | 1996-03-25 | 1999-09-27 | 株式会社先進材料利用ガスジェネレータ研究所 | 光ファイバセンサ |
TW466335B (en) | 2001-02-21 | 2001-12-01 | Asia Optical Co Inc | High accuracy laser range finder device |
US6486961B1 (en) * | 2001-05-08 | 2002-11-26 | Agilent Technologies, Inc. | System and method for measuring group delay based on zero-crossings |
JP2006266797A (ja) * | 2005-03-23 | 2006-10-05 | Anritsu Corp | 光ヘテロダイン干渉装置 |
US7345744B2 (en) * | 2005-06-06 | 2008-03-18 | Raytheon Company | Optical delay line to correct phase errors in coherent ladar |
EP1994424B1 (en) * | 2006-03-15 | 2012-01-11 | QUALCOMM Incorporated | Global navigation satellite system |
JP4769668B2 (ja) * | 2006-09-05 | 2011-09-07 | 日本電信電話株式会社 | 光リフレクトメトリ測定方法および装置 |
EP2315045B1 (de) * | 2009-10-22 | 2012-08-01 | Sick Ag | Messung von Entfernungen oder Entfernungsänderungen |
CN102495411B (zh) * | 2011-10-18 | 2013-05-29 | 中国科学院上海技术物理研究所 | 亚毫米级线性调谐激光测距系统及信号处理方法 |
US8665921B2 (en) | 2012-01-11 | 2014-03-04 | Realtek Semiconductor Corp. | Method and apparatus of automatic power control for burst mode laser transmitter |
US8772731B2 (en) * | 2012-04-15 | 2014-07-08 | Kla-Tencor Corporation | Apparatus and method for synchronizing sample stage motion with a time delay integration charge-couple device in a semiconductor inspection tool |
TWI516746B (zh) * | 2012-04-20 | 2016-01-11 | 賽格股份有限公司 | 在干涉編碼系統中執行非諧循環錯誤補償的方法、裝置及計算機程式產品,以及微影系統 |
US9064977B2 (en) | 2012-08-22 | 2015-06-23 | Freescale Semiconductor Inc. | Stacked microelectronic packages having sidewall conductors and methods for the fabrication thereof |
US10151827B2 (en) | 2014-07-25 | 2018-12-11 | DSCG Solutions, Inc. | Laser phase estimation and correction |
-
2015
- 2015-07-23 US US14/807,508 patent/US10151827B2/en active Active
- 2015-07-24 CN CN201911253119.9A patent/CN111025263A/zh active Pending
- 2015-07-24 TW TW104124152A patent/TWI648550B/zh active
- 2015-07-24 EP EP18182653.8A patent/EP3444630B1/en active Active
- 2015-07-24 JP JP2017504025A patent/JP6653695B2/ja active Active
- 2015-07-24 EP EP15747906.4A patent/EP3172586B1/en active Active
- 2015-07-24 CN CN201580045019.0A patent/CN106574965B/zh not_active Expired - Fee Related
- 2015-07-24 TW TW107136951A patent/TWI676044B/zh active
- 2015-07-24 WO PCT/US2015/041898 patent/WO2016014894A1/en active Application Filing
-
2018
- 2018-11-14 US US16/191,153 patent/US10739448B2/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030021309A1 (en) * | 2001-06-19 | 2003-01-30 | Mattox Barry G. | Method for producing highly accurate frequency and FM of a laser |
US20080088503A1 (en) * | 2004-09-28 | 2008-04-17 | Beasley Patrick D L | Frequency Modulated Continuous Wave (Fmcw) Radar Having Improved Frequency Linearity |
CN101266293A (zh) * | 2008-04-30 | 2008-09-17 | 西安电子科技大学 | 激光合成孔径雷达成像距离向相位补偿方法 |
CN102419441A (zh) * | 2011-09-01 | 2012-04-18 | 哈尔滨工业大学 | 一种基于四通道探测技术的弱光锁相星间位移测量方法及实现该方法的装置 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111025320A (zh) * | 2019-12-28 | 2020-04-17 | 深圳奥锐达科技有限公司 | 一种相位式激光测距系统及测距方法 |
CN112099036A (zh) * | 2020-11-10 | 2020-12-18 | 深圳市汇顶科技股份有限公司 | 距离测量方法以及电子设备 |
CN112099036B (zh) * | 2020-11-10 | 2021-03-23 | 深圳市汇顶科技股份有限公司 | 距离测量方法以及电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN111025263A (zh) | 2020-04-17 |
US10739448B2 (en) | 2020-08-11 |
TWI648550B (zh) | 2019-01-21 |
EP3444630B1 (en) | 2020-05-27 |
JP2017526911A (ja) | 2017-09-14 |
JP6653695B2 (ja) | 2020-02-26 |
EP3172586B1 (en) | 2018-09-19 |
TWI676044B (zh) | 2019-11-01 |
TW201903430A (zh) | 2019-01-16 |
CN106574965B (zh) | 2020-01-03 |
TW201604572A (zh) | 2016-02-01 |
US10151827B2 (en) | 2018-12-11 |
EP3444630A1 (en) | 2019-02-20 |
US20160025843A1 (en) | 2016-01-28 |
US20190101630A1 (en) | 2019-04-04 |
WO2016014894A1 (en) | 2016-01-28 |
EP3172586A1 (en) | 2017-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106574965A (zh) | 激光相位估计及校正 | |
US20230161038A1 (en) | Lidar system to adjust doppler effects | |
US20190339388A1 (en) | Method and System for Time Separated Quadrature Detection of Doppler Effects in Optical Range Measurements | |
JP6903743B2 (ja) | コヒーレントライダーから生じる信号を処理する方法および関連するライダーシステム | |
JP6236369B2 (ja) | 光ファイバの温度・歪み分布測定方法および装置 | |
CN107250723A (zh) | 多束范围测量过程 | |
US10365354B2 (en) | Multi-target laser distance meter | |
JP2019052938A (ja) | 光反射測定装置及びその方法 | |
EP2799838A1 (en) | Information signal generating method | |
JP5613627B2 (ja) | レーザ光コヒーレンス関数測定方法及び測定装置 | |
JP5827140B2 (ja) | レーザ光特性測定方法及び測定装置 | |
CN112654894A (zh) | 一种雷达探测方法及相关装置 | |
JP5371933B2 (ja) | レーザ光測定方法及びその測定装置 | |
JP5470320B2 (ja) | レーザ光コヒーレンス長測定方法及び測定装置 | |
JP2020027034A (ja) | 信号処理装置 | |
JP2003185737A (ja) | 信号分析装置、および方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200103 |