CN113687329B - 结合扫频和外差干涉仪的非合作目标测距系统及方法 - Google Patents
结合扫频和外差干涉仪的非合作目标测距系统及方法 Download PDFInfo
- Publication number
- CN113687329B CN113687329B CN202111063094.3A CN202111063094A CN113687329B CN 113687329 B CN113687329 B CN 113687329B CN 202111063094 A CN202111063094 A CN 202111063094A CN 113687329 B CN113687329 B CN 113687329B
- Authority
- CN
- China
- Prior art keywords
- beam splitter
- interferometer
- laser
- frequency
- target
- 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
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000005259 measurement Methods 0.000 claims abstract description 34
- 239000013307 optical fiber Substances 0.000 claims description 25
- 238000006243 chemical reaction Methods 0.000 claims description 16
- 238000001228 spectrum Methods 0.000 claims description 14
- 238000005070 sampling Methods 0.000 claims description 12
- 230000008878 coupling Effects 0.000 claims description 11
- 238000010168 coupling process Methods 0.000 claims description 11
- 238000005859 coupling reaction Methods 0.000 claims description 11
- 230000003287 optical effect Effects 0.000 claims description 11
- 238000012937 correction Methods 0.000 claims description 9
- 238000001514 detection method Methods 0.000 claims description 9
- CLOMYZFHNHFSIQ-UHFFFAOYSA-N clonixin Chemical compound CC1=C(Cl)C=CC=C1NC1=NC=CC=C1C(O)=O CLOMYZFHNHFSIQ-UHFFFAOYSA-N 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- PNEYBMLMFCGWSK-UHFFFAOYSA-N Alumina Chemical compound [O-2].[O-2].[O-2].[Al+3].[Al+3] PNEYBMLMFCGWSK-UHFFFAOYSA-N 0.000 claims description 5
- 238000001914 filtration Methods 0.000 claims description 5
- 238000000691 measurement method Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 9
- 238000010586 diagram Methods 0.000 description 6
- 239000000835 fiber Substances 0.000 description 4
- 238000004088 simulation Methods 0.000 description 3
- 238000005305 interferometry Methods 0.000 description 2
- 230000000712 assembly Effects 0.000 description 1
- 238000000429 assembly Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- 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
-
- 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
-
- 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/4818—Constructional features, e.g. arrangements of optical elements using optical fibres
-
- 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/495—Counter-measures or counter-counter-measures using electronic or electro-optical means
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
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)
- Instruments For Measurement Of Length By Optical Means (AREA)
Abstract
结合扫频和外差干涉仪的非合作目标测距系统及方法,本发明涉及非合作目标测距领域。解决了现有技术中对于非合作目标的测距系统存在测量精度低的问题。本发明中外差干涉仪与扫频干涉仪采用共路结构,可以同步测量目标的振动和距离,振动测量结果用于校正扫频干涉信号I1(k)中的由于多普勒效应引起的振动,校正振动对测距结果的影响。本发明主要用于对非合作目标距离的测量。
Description
技术领域
本发明涉及非合作目标测距领域。
背景技术
在大型组件正在组装或机器正在运行的工业环境(航空航天工业或其他重型制造)中执行高精度绝对距离测量的能力变得越来越重要。由于频率扫描激光器的发展,频率扫描干涉测量(FSI)作为一种具有高精度和无测距盲区的技术,在绝对距离测量中越来越受到关注。
对于基本FSI,其中测量干涉仪的光程差(OPD)变化,测量误差比实际OPD变化大数百倍。这是由振动引入多普勒效应引起的,也即由于振动导致测距准确度低。
为了克服这个问题,2001年,Schneider等人。使用两个同时在频率上向上和向下调谐的激光二极管,并使用辅助干涉仪来校正频率扫描非线性。
2014年,Matthew·Warden等人。使用两个外腔激光器开发了一种类似于RichardSchneider的设置。外腔激光器具有较大的扫频范围,可以提高FSI的精度。双频扫描方法可以获得较好的校正结果,但成本太高。
2018年,Zhang等人提出结合FSI和多波长干涉测量距离,其中HeNe激光器是用于监测测量距离的波动。
2020年,Bin等人提出使用定频干涉仪来消除多普勒效应误差。
2016年,我们使用外差干涉仪来校正FSI中的多普勒效应系统。
然而,上面显示的所有方法都旨在测量合作目标,例如:球面反射目标或光滑镜面,并且忽略了多径干扰的影响。在这种情况下,当测量粗糙的物体(非合作目标)时,由于不可忽略多径干扰的影响,受多径干扰的影响导致对非合作目标测量精度低,特别是当距离大于10m时精度不超过10μm。因此,针对现有技术中对于非合作目标的测距系统存在测量精度低的缺陷亟需解决。
发明内容
本发明目的是为了解决现有技术中对于非合作目标的测距系统存在测量精度低的问题,本发明提供了一种结合扫频和外差干涉仪的非合作目标测距系统及方法。
结合扫频和外差干涉仪的非合作目标测距系统,包括外腔式激光器、单频激光器、第一辅助干涉仪、第二辅助干涉仪、测量干涉仪、1号分束器、2号分束器、3号分束器、4号分束器、5号分束器、声光调制器AOM、数据采集器和处理器;
第一辅助干涉仪和测量干涉仪构成扫频干涉仪;
第二辅助干涉仪和测量干涉仪构成外差干涉仪;
外腔式激光器输出的激光经1号分束器分束后,生成的两路激光,其中,一路送至第一辅助干涉仪进行干涉及光电转化后形成的信号送至数据采集器,作为数据采集器的采样时钟,另一路送入2号分束器进行分束后形成的两路激光均送至测量干涉仪,且2号分束器输出的一路激光作为探测光,其输出的另一路激光作为参考光;
单频激光器输出的激光经3号分束器分束后,生成两路激光,其中,
3号分束器输出的一路激光经过4号分束器分束成两路激光,且4号分束器输出的一路激光送至测量干涉仪,4号分束器输出的另一路激光送至第二辅助干涉仪;
3号分束器输出的另一路激光经声光调制器AOM调频后送至5号分束器,5号分束器对接收的激光进行分束后形成两路激光分别送至测量干涉仪和第二辅助干涉仪;
测量干涉仪,用于根据接收的4路激光对被测目标进行探测,获得携带有振动信息的距离信号,并将该信号通过数据采集器送至处理器;
第二辅助干涉仪,用于对接收的两路激光进行干涉及光电转化后,生成差频信号Ia(k),并将该差频信号Ia(k)通过数据采集器送至处理器;其中,k为整数,表示采样点;
处理器,利用差频信号Ia(k)对携带有振动信息的距离信号进行振动补偿,从而获得目标真实距离Rm0。
优选的是,第一辅助干涉仪包括6号分束器、1号耦合器和1号光电探测器PD1;
6号分束器,用于对1号分束器输出的一路激光进行分束后,形成的两束激光通过光纤送至1号耦合器进行耦合后,送至1号光电探测器PD1进行光电转化,并将形成的信号送至数据采集器。
优选的是,测量干涉仪包括1号波分复用器、2号波分复用器、2号耦合器、光纤循环器、光纤尾纤FP、汇聚透镜组和2号光电探测器PD2;
1号波分复用器,用于对2号分束器输出的探测光和4号分束器输出的一路激光进行耦合后,送至光纤循环器,光纤循环器的第一输出端输出激光通过光纤尾纤FP送至汇聚透镜组汇聚后入射至被测目标,经被测目标反射的光按原路依次经过汇聚透镜组、光纤尾纤FP返回至光纤循环器,并通过光纤循环器的第二输出端送至2号耦合器;
2号波分复用器,用于对2号分束器输出的参考光和5号分束器输出的一路激光进行耦合后送至2号耦合器;
2号耦合器对接收的两路激光进行耦合后,送至2号光电探测器PD2进行光电转化后,生成携带有振动信息的距离信号。
优选的是,第二辅助干涉仪包括3号耦合器和3号光电探测器PD3;
3号耦合器,用于对接收的4号分束器输出的另一路激光和5号分束器输出的另一路激光耦合后,送至3号光电探测器PD3进行光电转化生成差频信号Ia(k)。
优选的是,1号分束器、2号分束器、3号分束器和4号分束器的分光比均为99:1,5号分束器的分光比为50:50的分束器;
6号分束器和1号耦合器的分光比均为50:50;
2号耦合器的分光比为50:50;
3号耦合器的分光比为50:50。
采用所述的结合扫频和外差干涉仪的非合作目标测距系统实现的测距方法,该方法包括如下步骤:
S1、处理器接收到差频信号Ia(k)和携带有振动信息的距离信号之后,先对携带有振动信息的距离信号进行滤波,获得扫频干涉信号I1(k)和外差干涉信号I2(k);
S2、利用差频信号Ia(k)对外差干涉信号I2(k)进行解调,从而获得II(k)和IQ(k),再利用最小二乘法对II(k)和IQ(k)进行圆拟合,从而获得α和β的值;
其中,
II(k)为外差干涉信号I2(k)的同相分量;
IQ(k)为外差干涉信号I2(k)的正交分量;
α为同相分量II(k)的直流偏置;
β为正交分量IQ(k)的直流偏置;
S3、利用S2中获得的α和β的值,对II(k)和IQ(k)进行修正及处理,从而获得目标振动信息ΔRm(k);
S4、利用目标振动信息ΔRm(k)对扫频干涉信号I1(k)进行振动补偿,从而获得目标绝对距离Icor(k),再对目标绝对距离Icor(k)进行线性调频Z变换,从而获得距离谱,并将距离谱中的最大值作为目标真实距离Rm0。
优选的是,S1中,差频信号Ia(k)的表达式为:Ia(k)=Aacos[2πΔfAOMt(k)] (公式一);
其中,
Aa为第二辅助干涉仪的幅值;
ΔfAOM为声光调制器AOM的差频;
t(k)为第k个采样点所对应的时间;
S1中,对携带有振动信息的距离信号进行滤波,获得扫频干涉信号I1(k)和外差干涉信号I2(k)的实现方式为:
其中,
A1为扫频干涉信号I1(k)的幅值;
R0为第一辅助干涉仪的光程差;
c为光速;
f0为外腔式激光器的起始频率;
A2为外差干扰信号I2(k)的幅值;
fAOM为声光调制器AOM的频移;
t(k)为第k个采样点所对应的时间;
f1为单频激光器的光频。
优选的是,S2中,
其中,
A2为外差干扰信号I2(k)的幅值;
f为外腔式激光器(1)的瞬时频率;
A3为杂散光干扰信号的幅值;
为杂散光干扰信号的初相位;
α为同相分量II(k)的直流偏置,
β为正交分量IQ(k)的直流偏置,
优选的是,S3中,利用S2中获得的α和β的值,对II(k)和IQ(k)进行修正及处理,从而获得目标振动信息ΔRm(k)的实现方式为:
其中,f1为单频激光器的光频。
c为光速。
优选的是,S4中,利用目标振动信息ΔRm(k)对扫频干涉信号I1(k)进行振动补偿,从而获得目标绝对距离Icor(k)的实现方式为:
其中,HT[·]为希尔伯特变换;
所述的目标真实距离Rm0的表达式为:Rm0=arcmax(czt(Icor(k))),其中,
czt(·)为线性调频Z变换。
本发明提供了一种结合扫频和外差干涉仪的非合作目标测距系统,该系统与外差干涉仪相结合。外差干涉仪与扫频干涉仪采用共路结构,可以同步测量目标的振动和距离,振动测量结果用于校正扫频干涉信号I1(k)中的由于多普勒效应引起的振动。
为了校正扫频干涉信号I1(k),本发明采用附加外差干涉仪来监测目标的振动,其中,单频激光器的输出分为两个部分,其中一部分进入第一辅助干涉仪和测量干涉仪,另一部分经过声光调制器AOM后再进入第二辅助干涉仪和测量干涉仪。
本发明所述采用结合扫频和外差干涉仪的非合作目标测距系统实现的测距方法中,使用与扫频干涉仪平行的附加外差干涉仪来测量目标的振动,分析多径干扰的影响,然后采用后处理算法校正多普勒效应并提取距离,具体为,在使用圆拟合来校正多径干扰的影响,多径干扰会损害振动测量结果,尤其是在非合作目标测量中。然后,通过振动补偿消除了含多径效应的振动干扰,最终提高了测距精度。
附图说明
图1是本发明所述结合扫频和外差干涉仪的非合作目标测距系统的结构示意图;
图2是利用差频信号Ia(k)对外差干涉信号I2(k)进行解调获得II(k)和IQ(k)的原理示意图;其中,LPF为低通滤波器,HT为希尔伯特变换;
图3是获得目标真实距离的原理示意图;
图4为应用本发明测距方法对由于多普勒效应引起的振动效应的校正仿真图;
图4(a)为外差干扰信号I2(k)的波形图;
图4(b)为同相分量II(k)和正交分量IQ(k)的波形图;
图4(c)为同相分量II(k)和正交分量IQ(k)的李萨如曲线;
图4(d)为未校正多径干涉得到的目标振动曲线;
图4(e)为采用圆拟合校正多径干涉后的目标振动曲线;
图4(f)为三种状态下的距离谱,其中,三种状态下的距离谱分别为理想距离谱,采用圆拟合校正多径干涉影响后的距离谱,未校正多径干涉影响的距离谱。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
需要说明的是,在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互组合。
具体实施方式一:下面结合图1说明本实施方式,本实施方式所述的结合扫频和外差干涉仪的非合作目标测距系统,包括外腔式激光器1、单频激光器2、第一辅助干涉仪3、第二辅助干涉仪4、测量干涉仪5、1号分束器6、2号分束器7、3号分束器8、4号分束器9、5号分束器10、声光调制器AOM11、数据采集器和处理器;
第一辅助干涉仪3和测量干涉仪5构成扫频干涉仪;
第二辅助干涉仪4和测量干涉仪5构成外差干涉仪;
外腔式激光器1输出的激光经1号分束器6分束后,生成的两路激光,其中,一路送至第一辅助干涉仪3进行干涉及光电转化后形成的信号送至数据采集器,作为数据采集器的采样时钟,另一路送入2号分束器7进行分束后形成的两路激光均送至测量干涉仪5,且2号分束器7输出的一路激光作为探测光,其输出的另一路激光作为参考光;
单频激光器2输出的激光经3号分束器8分束后,生成两路激光,其中,
3号分束器8输出的一路激光经过4号分束器9分束成两路激光,且4号分束器9输出的一路激光送至测量干涉仪5,4号分束器9输出的另一路激光送至第二辅助干涉仪4;
3号分束器8输出的另一路激光经声光调制器AOM11调频后送至5号分束器10,5号分束器10对接收的激光进行分束后形成两路激光分别送至测量干涉仪5和第二辅助干涉仪4;
测量干涉仪5,用于根据接收的4路激光对被测目标进行探测,获得携带有振动信息的距离信号,并将该信号通过数据采集器送至处理器;
第二辅助干涉仪4,用于对接收的两路激光进行干涉及光电转化后,生成差频信号Ia(k),并将该差频信号Ia(k)通过数据采集器送至处理器;其中,k为整数,表示采样点;
处理器,利用差频信号Ia(k)对携带有振动信息的距离信号进行振动补偿,从而获得目标真实距离Rm0。
本实施方式中,本发明提供了一种结合扫频和外差干涉仪的非合作目标测距系统,该系统与外差干涉仪相结合。外差干涉仪与扫频干涉仪采用共路结构,可以同步测量目标的振动和距离,振动测量结果用于校正扫频干涉信号I1(k)中的由于多普勒效应引起的振动。
进一步的,第一辅助干涉仪3包括6号分束器3-1、1号耦合器3-2和1号光电探测器PD1;
6号分束器3-1,用于对1号分束器6输出的一路激光进行分束后,形成的两束激光通过光纤送至1号耦合器3-2进行耦合后,送至1号光电探测器PD1进行光电转化,并将形成的信号送至数据采集器。
本优选实施方式中1号光电探测器PD1为平衡光电探测器。
更进一步的,测量干涉仪5包括1号波分复用器5-1、2号波分复用器5-2、2号耦合器5-3、光纤循环器5-4、光纤尾纤FP、汇聚透镜组5-5和2号光电探测器PD2;
1号波分复用器5-1,用于对2号分束器7输出的探测光和4号分束器9输出的一路激光进行耦合后,送至光纤循环器5-4,光纤循环器5-4的第一输出端输出激光通过光纤尾纤FP送至汇聚透镜组5-5汇聚后入射至被测目标,经被测目标反射的光按原路依次经过汇聚透镜组5-5、光纤尾纤FP返回至光纤循环器5-4,并通过光纤循环器5-4的第二输出端送至2号耦合器5-3;
2号波分复用器5-2,用于对2号分束器7输出的参考光和5号分束器10输出的一路激光进行耦合后送至2号耦合器5-3;
2号耦合器5-3对接收的两路激光进行耦合后,送至2号光电探测器PD2进行光电转化后,生成携带有振动信息的距离信号。
本优选实施方式中,给出了测量干涉仪5的一种具体结构,结构简单,便于实现,且2号光电探测器PD2优选为平衡光电探测器。
更进一步的,第二辅助干涉仪4包括3号耦合器4-1和3号光电探测器PD3;
3号耦合器4-1,用于对接收的4号分束器9输出的另一路激光和5号分束器10输出的另一路激光耦合后,送至3号光电探测器PD3进行光电转化生成差频信号Ia(k)。
本优选实施方式中,给出了第二辅助干涉仪4的结构,结构简单,便于实现,且3号光电探测器PD3优选为平衡光电探测器。
更进一步的,1号分束器6、2号分束器7、3号分束器8和4号分束器9的分光比均为99:1,5号分束器10的分光比为50:50的分束器;
6号分束器3-1和1号耦合器3-2的分光比均为50:50;
2号耦合器5-3的分光比为50:50;
3号耦合器4-1的分光比为50:50。
本优选实施方式中,给出了各分束器和耦合器的最优分光比。
具体实施方式二:下面结合图1说明本实施方式,采用具体实施方式一所述的结合扫频和外差干涉仪的非合作目标测距系统实现的测距方法,该方法包括如下步骤:
S1、处理器接收到差频信号Ia(k)和携带有振动信息的距离信号之后,先对携带有振动信息的距离信号进行滤波,获得扫频干涉信号I1(k)和外差干涉信号I2(k);
S2、利用差频信号Ia(k)对外差干涉信号I2(k)进行解调,从而获得II(k)和IQ(k),再利用最小二乘法对II(k)和IQ(k)进行圆拟合,从而获得α和β的值;
其中,
II(k)为外差干涉信号I2(k)的同相分量;
IQ(k)为外差干涉信号I2(k)的正交分量;
α为同相分量II(k)的直流偏置;
β为正交分量IQ(k)的直流偏置;
S3、利用S2中获得的α和β的值,对II(k)和IQ(k)进行修正及处理,从而获得目标振动信息ΔRm(k);
S4、利用目标振动信息ΔRm(k)对扫频干涉信号I1(k)进行振动补偿,从而获得目标绝对距离Icor(k),再对目标绝对距离Icor(k)进行线性调频Z变换,从而获得距离谱,并将距离谱中的最大值作为目标真实距离。
进一步的,S1中,差频信号Ia(k)的表达式为:Ia(k)=Aacos[2πΔfAOMt(k)] (公式一);
其中,
Aa为第二辅助干涉仪(4)的幅值;
ΔfAOM为声光调制器AOM(11)的差频;
t(k)为第k个采样点所对应的时间;
进一步的,S1中,对携带有振动信息的距离信号进行滤波,获得扫频干涉信号I1(k)和外差干涉信号I2(k)的实现方式为:
其中,
A1为扫频干涉信号I1(k)的幅值;
R0为第一辅助干涉仪3的光程差;
c为光速;
f0为外腔式激光器1的起始频率;
A2为外差干扰信号I2(k)的幅值;
fAOM为声光调制器AOM 11的频移;
t(k)为第k个采样点所对应的时间;
f1为单频激光器2的光频。
更进一步的,具体参见图2,S2中,
其中,
Aa为第二辅助干涉仪4的幅值;
A2为外差干扰信号I2(k)的幅值;
c为光速;
f为外腔式激光器1的瞬时频率;
A3为杂散光干扰信号的幅值;
为杂散光干扰信号的初相位;
α为同相分量II(k)的直流偏置,
f1为单频激光器2的光频;
β为正交分量IQ(k)的直流偏置,
更进一步的,具体参见图2和图3,S3中,利用S2中获得的α和β的值,对II(k)和IQ(k)进行修正及处理,从而获得目标振动信息ΔRm(k)的实现方式为:
其中,f1为单频激光器2的光频。
c为光速。
更进一步的,参见图3,S4中,利用目标振动信息ΔRm(k)对扫频干涉信号I1(k)进行振动补偿,从而获得目标绝对距离Icor(k)的实现方式为:
其中,HT[·]为希尔伯特变换;
更进一步的,所述的目标真实距离Rm0的表达式为:Rm0=arcmax(czt(Icor(k))),
其中,czt(·)为线性调频Z变换。
具体应用时,应用本发明所述的采用结合扫频和外差干涉仪的非合作目标测距系统的测量方法进行应用时仿真图,具体参见图4,在仿真中,目标位于10m处,多径干扰信号的振幅是有用信号的0.8倍。外差干扰信号如图4(a)所示。为了提取变化,应首先获得II(k)和IQ(k)信号,如图4(b)所示。受多径干扰的影响,II(k)和IQ(k)信号有明显的偏差。在这种影响下,李萨如曲线如图4(c)所示,圆心偏离零点。采用最小二乘法得到拟合圆,拟合中心为(-0.21,0.34)。如果忽略多径干扰的影响,计算的振动曲线如图4(d)所示。利用曲线的拟合中心,可通过公式六获得精确计算的振动曲线,也即:目标振动信息ΔRm(k),结果如图4(e)所示,其与预设振动信息相似。计算出的振动可用于校正扫频干涉信号I1(k)中的多普勒效应,如公式七所示,校正多径干扰影响后的距离谱如图4(f)中的采用圆拟合校正多径干涉影响后的距离谱,图4(f)中没有校正多径干扰影响的距离谱,多径干扰具有明显的失真。
虽然在本文中参照了特定的实施方式来描述本发明,但是应该理解的是,这些实施例仅仅是本发明的原理和应用的示例。因此应该理解的是,可以对示例性的实施例进行许多修改,并且可以设计出其他的布置,只要不偏离所附权利要求所限定的本发明的精神和范围。应该理解的是,可以通过不同于原始权利要求所描述的方式来结合不同的从属权利要求和本文中所述的特征。还可以理解的是,结合单独实施例所描述的特征可以使用在其他所述实施例中。
Claims (9)
1.采用结合扫频和外差干涉仪的非合作目标测距系统实现的测距方法,结合扫频和外差干涉仪的非合作目标测距系统包括外腔式激光器(1)、单频激光器(2)、第一辅助干涉仪(3)、第二辅助干涉仪(4)、测量干涉仪(5)、1号分束器(6)、2号分束器(7)、3号分束器(8)、4号分束器(9)、5号分束器(10)、声光调制器AOM(11)、数据采集器和处理器;
第一辅助干涉仪(3)和测量干涉仪(5)构成扫频干涉仪;
第二辅助干涉仪(4)和测量干涉仪(5)构成外差干涉仪;
外腔式激光器(1)输出的激光经1号分束器(6)分束后,生成的两路激光,其中,一路送至第一辅助干涉仪(3)进行干涉及光电转化后形成的信号送至数据采集器,作为数据采集器的采样时钟,另一路送入2号分束器(7)进行分束后形成的两路激光均送至测量干涉仪(5),且2号分束器(7)输出的一路激光作为探测光,其输出的另一路激光作为参考光;
单频激光器(2)输出的激光经3号分束器(8)分束后,生成两路激光,其中,
3号分束器(8)输出的一路激光经过4号分束器(9)分束成两路激光,且4号分束器(9)输出的一路激光送至测量干涉仪(5),4号分束器(9)输出的另一路激光送至第二辅助干涉仪(4);
3号分束器(8)输出的另一路激光经声光调制器AOM(11)调频后送至5号分束器(10),5号分束器(10)对接收的激光进行分束后形成两路激光分别送至测量干涉仪(5)和第二辅助干涉仪(4);
测量干涉仪(5),用于根据接收的4路激光对被测目标进行探测,获得携带有振动信息的距离信号,并将该信号通过数据采集器送至处理器;
第二辅助干涉仪(4),用于对接收的两路激光进行干涉及光电转化后,生成差频信号Ia(k),并将该差频信号Ia(k)通过数据采集器送至处理器;其中,k为整数,表示采样点;
处理器,利用差频信号Ia(k)对携带有振动信息的距离信号进行振动补偿,从而获得目标真实距离Rm0;
其特征在于,测距方法包括如下步骤:
S1、处理器接收到差频信号Ia(k)和携带有振动信息的距离信号之后,先对携带有振动信息的距离信号进行滤波,获得扫频干涉信号I1(k)和外差干涉信号I2(k);
S2、利用差频信号Ia(k)对外差干涉信号I2(k)进行解调,从而获得II(k)和IQ(k),再利用最小二乘法对II(k)和IQ(k)进行圆拟合,从而获得α和β的值;
其中,
II(k)为外差干涉信号I2(k)的同相分量;
IQ(k)为外差干涉信号I2(k)的正交分量;
α为同相分量II(k)的直流偏置;
β为正交分量IQ(k)的直流偏置;
S3、利用S2中获得的α和β的值,对II(k)和IQ(k)进行修正及处理,从而获得目标振动信息ΔRm(k);
S4、利用目标振动信息ΔRm(k)对扫频干涉信号I1(k)进行振动补偿,从而获得目标绝对距离Icor(k),再对目标绝对距离Icor(k)进行线性调频Z变换,从而获得距离谱,并将距离谱中的最大值作为目标真实距离Rm0。
2.根据权利要求1所述的采用结合扫频和外差干涉仪的非合作目标测距系统实现的测距方法,其特征在于,第一辅助干涉仪(3)包括6号分束器(3-1)、1号耦合器(3-2)和1号光电探测器PD1;
6号分束器(3-1),用于对1号分束器(6)输出的一路激光进行分束后,形成的两束激光通过光纤送至1号耦合器(3-2)进行耦合后,送至1号光电探测器PD1进行光电转化,并将形成的信号送至数据采集器。
3.根据权利要求2所述的采用结合扫频和外差干涉仪的非合作目标测距系统实现的测距方法,其特征在于,测量干涉仪(5)包括1号波分复用器(5-1)、2号波分复用器(5-2)、2号耦合器(5-3)、光纤循环器(5-4)、光纤尾纤FP、汇聚透镜组(5-5)和2号光电探测器PD2;
1号波分复用器(5-1),用于对2号分束器(7)输出的探测光和4号分束器(9)输出的一路激光进行耦合后,送至光纤循环器(5-4),光纤循环器(5-4)的第一输出端输出激光通过光纤尾纤FP送至汇聚透镜组(5-5)汇聚后入射至被测目标,经被测目标反射的光按原路依次经过汇聚透镜组(5-5)、光纤尾纤FP返回至光纤循环器(5-4),并通过光纤循环器(5-4)的第二输出端送至2号耦合器(5-3);
2号波分复用器(5-2),用于对2号分束器(7)输出的参考光和5号分束器(10)输出的一路激光进行耦合后送至2号耦合器(5-3);
2号耦合器(5-3)对接收的两路激光进行耦合后,送至2号光电探测器PD2进行光电转化后,生成携带有振动信息的距离信号。
4.根据权利要求3所述的采用结合扫频和外差干涉仪的非合作目标测距系统实现的测距方法,其特征在于,第二辅助干涉仪(4)包括3号耦合器(4-1)和3号光电探测器PD3;
3号耦合器(4-1),用于对接收的4号分束器(9)输出的另一路激光和5号分束器(10)输出的另一路激光耦合后,送至3号光电探测器PD3进行光电转化生成差频信号Ia(k)。
5.根据权利要求4所述的采用结合扫频和外差干涉仪的非合作目标测距系统实现的测距方法,其特征在于,1号分束器(6)、2号分束器(7)、3号分束器(8)和4号分束器(9)的分光比均为99:1,5号分束器(10)的分光比为50:50的分束器;
6号分束器(3-1)和1号耦合器(3-2)的分光比均为50:50;
2号耦合器(5-3)的分光比为50:50;
3号耦合器(4-1)的分光比为50:50。
6.根据权利要求1所述的采用结合扫频和外差干涉仪的非合作目标测距系统实现的测距方法,其特征在于,
S1中,差频信号Ia(k)的表达式为:Ia(k)=Aacos[2πΔfAOMt(k)] (公式一);
其中,
Aa为第二辅助干涉仪(4)的幅值;
ΔfAOM为声光调制器AOM(11)的差频;
t(k)为第k个采样点所对应的时间;
S1中,对携带有振动信息的距离信号进行滤波,获得扫频干涉信号I1(k)和外差干涉信号I2(k)的实现方式为:
其中,
A1为扫频干涉信号I1(k)的幅值;
R0为第一辅助干涉仪(3)的光程差;
c为光速;
f0为外腔式激光器(1)的起始频率;
A2为外差干扰信号I2(k)的幅值;
fAOM为声光调制器AOM(11)的频移;
t(k)为第k个采样点所对应的时间;
f1为单频激光器(2)的光频。
7.根据权利要求6所述的采用结合扫频和外差干涉仪的非合作目标测距系统实现的测距方法,其特征在于,S2中,
其中,
A2为外差干扰信号I2(k)的幅值;
f为外腔式激光器(1)的瞬时频率;
A3为杂散光干扰信号的幅值;
为杂散光干扰信号的初相位;
α为同相分量II(k)的直流偏置,
β为正交分量IQ(k)的直流偏置,
8.根据权利要求1所述的采用结合扫频和外差干涉仪的非合作目标测距系统实现的测距方法,其特征在于,S3中,利用S2中获得的α和β的值,对II(k)和IQ(k)进行修正及处理,从而获得目标振动信息ΔRm(k)的实现方式为:
其中,f1为单频激光器(2)的光频;
c为光速。
9.根据权利要求6所述的采用结合扫频和外差干涉仪的非合作目标测距系统实现的测距方法,其特征在于,S4中,利用目标振动信息ΔRm(k)对扫频干涉信号I1(k)进行振动补偿,从而获得目标绝对距离Icor(k)的实现方式为:
其中,HT[·]为希尔伯特变换;
所述的目标真实距离Rm0的表达式为:Rm0=arcmax(czt(Icor(k))),其中,
czt(·)为线性调频Z变换。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111063094.3A CN113687329B (zh) | 2021-09-10 | 2021-09-10 | 结合扫频和外差干涉仪的非合作目标测距系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111063094.3A CN113687329B (zh) | 2021-09-10 | 2021-09-10 | 结合扫频和外差干涉仪的非合作目标测距系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113687329A CN113687329A (zh) | 2021-11-23 |
CN113687329B true CN113687329B (zh) | 2023-12-05 |
Family
ID=78586010
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111063094.3A Active CN113687329B (zh) | 2021-09-10 | 2021-09-10 | 结合扫频和外差干涉仪的非合作目标测距系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113687329B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115327515B (zh) * | 2022-08-10 | 2023-04-18 | 哈尔滨工业大学 | 一种基于相位传递的双扫频干涉动态测量系统及测量方法 |
CN115327514B (zh) * | 2022-08-10 | 2023-05-05 | 哈尔滨工业大学 | 一种基于相位传递的扫频干涉动态测量系统及测量方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011118255A1 (ja) * | 2010-03-26 | 2011-09-29 | 株式会社日立製作所 | 距離測定装置および距離測定方法 |
CN109029246A (zh) * | 2018-09-11 | 2018-12-18 | 哈尔滨工业大学 | 基于光学分频锁相非线性校正的动态扫频干涉测距系统及测距方法 |
CN109188454A (zh) * | 2018-09-11 | 2019-01-11 | 哈尔滨工业大学 | 基于数字锁相非线性校正的动态扫频干涉测距系统及方法 |
CN109188453A (zh) * | 2018-09-11 | 2019-01-11 | 哈尔滨工业大学 | 基于锁相非线性校正的动态扫频干涉测距系统及测距方法 |
CN110132138A (zh) * | 2019-05-14 | 2019-08-16 | 宁波核芯光电科技有限公司 | 基于级联干涉仪的双扫频光源测距系统及方法 |
CN110716208A (zh) * | 2019-10-18 | 2020-01-21 | 中国科学院光电技术研究所 | 一种基于声光调制器的扫频干涉动态绝对距离测量系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060061769A1 (en) * | 2004-09-03 | 2006-03-23 | California Institute Of Technology | Homodyne based optical coherence microscope |
-
2021
- 2021-09-10 CN CN202111063094.3A patent/CN113687329B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011118255A1 (ja) * | 2010-03-26 | 2011-09-29 | 株式会社日立製作所 | 距離測定装置および距離測定方法 |
CN109029246A (zh) * | 2018-09-11 | 2018-12-18 | 哈尔滨工业大学 | 基于光学分频锁相非线性校正的动态扫频干涉测距系统及测距方法 |
CN109188454A (zh) * | 2018-09-11 | 2019-01-11 | 哈尔滨工业大学 | 基于数字锁相非线性校正的动态扫频干涉测距系统及方法 |
CN109188453A (zh) * | 2018-09-11 | 2019-01-11 | 哈尔滨工业大学 | 基于锁相非线性校正的动态扫频干涉测距系统及测距方法 |
CN110132138A (zh) * | 2019-05-14 | 2019-08-16 | 宁波核芯光电科技有限公司 | 基于级联干涉仪的双扫频光源测距系统及方法 |
CN110716208A (zh) * | 2019-10-18 | 2020-01-21 | 中国科学院光电技术研究所 | 一种基于声光调制器的扫频干涉动态绝对距离测量系统 |
Non-Patent Citations (1)
Title |
---|
"基于宽带扫描干涉的高精度绝对距离测量方法研究";路程;《中国博士学位论文全文数据库工程科技II辑》(第1期);正文38-48页,图3-9,图3-11 * |
Also Published As
Publication number | Publication date |
---|---|
CN113687329A (zh) | 2021-11-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113687329B (zh) | 结合扫频和外差干涉仪的非合作目标测距系统及方法 | |
CN110068828B (zh) | 基于激光调频连续波远距离测距的装置及色散补偿方法 | |
US8144334B2 (en) | Fiber-optic, digital system for laser Doppler vibrometers (LDVs) | |
CN108873007B (zh) | 一种抑制振动效应的调频连续波激光测距装置 | |
CN108120378B (zh) | 基于飞秒光频梳的正弦相位调制干涉绝对测距装置与方法 | |
GB2529752A (en) | Dual laser frequency sweep interferometry system and method | |
CN110132138B (zh) | 基于级联干涉仪的双扫频光源测距系统及方法 | |
WO2009017291A1 (en) | Scanning micrometer using heterodyne interferometer | |
CN113029034A (zh) | 激光扫频干涉测量的非线性校正与量程扩展装置及方法 | |
Shao et al. | Dynamic clearance measurement using fiber-optic frequency-swept and frequency-fixed interferometry | |
CN110995341A (zh) | 基于光载微波干涉的光纤时延测量方法及装置 | |
CN116428966A (zh) | 一种三自由度激光干涉测量装置 | |
CN109541621B (zh) | 一种频率扫描干涉绝对测距系统的振动补偿方法 | |
CN109084884B (zh) | 一种零差激光测振装置及其振动检测方法 | |
CN114754689A (zh) | 一种基于双电光外差调制的相位式距离测量装置及方法 | |
Tomic et al. | Low-coherence interferometric method for measurement of displacement based on a 3× 3 fibre-optic directional coupler | |
CN108845333B (zh) | 一种抑制振动效应的调频连续波激光测距方法 | |
CN116819551A (zh) | 一种基于差分波前测量物体横向位移的装置 | |
CN115452015A (zh) | 一种双标度参考干涉的相噪精确校正光频域反射计 | |
CN113607277B (zh) | 一种窄线宽激光器线宽测量系统的解调方法 | |
CN115711633A (zh) | 一种环路结构参考干涉仪的相位噪声精确校正光频域反射计 | |
CN114720998A (zh) | 基于模态分解的非接触式动态位移测量系统 | |
US20140111804A1 (en) | Heterodyne Optical Spectrum Analyzer | |
CN112816180A (zh) | 光纤色散测量方法及测量装置 | |
KR20100077861A (ko) | 헤테로다인 레이저 간섭계 오차 보정 장치 및 그 방법 |
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 |