[0055] 下面结合附图和实施例对本发明加以详细说明。
[0056] 如图1所示，本发明具体实施包括光源系统I、正弦相位调制干涉系统II和信号处理系统III；光源系统I输出端连接到正弦相位调制干涉系统II的输入端，正弦相位调制干涉系统II输出干涉光束到信号处理系统III；光源系统I输出包含两束波长不同的激光光束射向正弦相位调制干涉系统II，正弦相位调制干涉系统II中形成干涉光束并被信号处理系统III处理，进行PGC相位解调，根据频率扫描干涉FSI和多波长干涉MWI测量原理得到待测距离的测量结果。
[0057] 如图1所示，光源系统I包括飞秒光频梳1、第一激光器2、第二激光器3、反射镜4和第一分色镜5；输出波长为λa的第一激光器2与输出波长为λb的第二激光器3均锁定至飞秒光频梳1上，第二激光器3输出端发出波长为λb的激光经反射镜4反射后与第一激光器2输出波长为λa的激光一起在第一分色镜5处合光，第一分色镜5合束光输出作为所述的光源系统I的输出光。
[0058] 如图1所示，正弦相位调制干涉系统II包括分光镜6、参考角锥棱镜7、相位调制器8、测量角锥棱镜11、信号发生器9和高压放大器10；光源系统I输出光射向分光镜6发生透射和反射，分光镜6的透射光经参考角锥棱镜7反射到电光相位调制器8，经电光相位调制器8调制后输出返回入射到分光镜6；分光镜6的反射光射向测量角锥棱镜11，经测量角锥棱镜
11反射后与电光相位调制器8输出光回到分光镜6处合光输出，电光相位调制器8输出光到分光镜6发生反射，测量角锥棱镜11反射光到分光镜6发生透射，形成干涉信号；电光相位调制器8经高压放大器10与信号发生器9的输出端连接。
[0059] 如图1所示，信号处理系统III包括第二分色镜13、第一滤光片14、第一光电探测器15、第二滤光片16、第二光电探测器17、数据采集模块18和计算机19；正弦相位调制干涉系统II的干涉信号入射到第二分色镜13发生透射和反射，波长为λa的光被第二分色镜13透射，第二分色镜13的透射光经第一滤光片14后被第一光电探测器15探测接收，波长为λb的光被第二分色镜13反射，第二分色镜13的反射光经第二滤光片16后被第二光电探测器17探测接收。
[0060] 具体实施还包括数据采集模块18和计算机19，所述的信号发生器9具有两个输出端，信号发生器9一个输出端经高压放大器10电光相位调制器8连接，信号发生器9另一个输出端连接到数据采集模块18的输入端，第一光电探测器15和第二光电探测器17连接到数据采集模块18的输入端，数据采集模块18的输出端连接到计算机19；分光镜6和参考角锥棱镜7之间以及分光镜6和测量角锥棱镜11之间设有用于检测温度、湿度、压强、二氧化碳等环境参数值的环境参数测量模块12，温度、湿度、压强、二氧化碳等环境参数值由环境参数测量模块12实时监测，环境参数测量模块12连接到计算机19。
[0061] 信号发生器9产生高频正弦信号一分为二从两个输出端输出，一路经高压放大器10放大后施加于电光相位调制器8，另一路输入至数据采集模块18用于信号解调。
[0062] 信号处理系统III中通过数据采集模块18同时采集信号发生器9输出信号以及第一光电探测器15和第二光电探测器17的信号，然后传输至计算机19进行信号处理，根据PGC解调原理对测量干涉仪和补偿干涉仪对应的干涉信号进行实时相位解调，得到待测绝对距离的值。
[0063] 具体实施中，飞秒光频梳1采用德国MenloSystem公司的FC1500-250-WG，包括HMP633模块和SHG780模块，第一激光器2采用美国NewFocus公司的TLB6700型可调谐半导体激光器，输出固定波长λa为633nm，第二激光器3采用另一台美国NewFocus公司的TLB6700型可调谐半导体激光器，输出波长λb的范围为765-781nm，环境参数测量模块12为美国Fluke公司的1586A高精度多通道环境参数测量仪，数据采集模块18为凌华科技有限公司的PCI-9812型四通道数据采集卡，计算机19采用惠普公司的Pro4500台式机。
[0064] 本发明基于飞秒光频梳的正弦相位调制干涉绝对测距的具体实施过程如下：
[0065] 1)t0时刻，将第一激光器2和第二激光器3分别锁定至飞秒光频梳1上，通过飞秒光频梳1锁定第一激光器2和第二激光器3输出激光的频率分别为νa和νb0；
[0066] 通过计算机19记录t0时刻第一激光器2和第二激光器3输出频率为νa、νb0的两路激光各自对应的干涉条纹小数εa0、εb0，0≤εa0,εb0<1；并且由环境参数测量模块12采集到的环境参数值(温度、湿度、压强、CO2含量)，由环境参数值转换计算获得t0时刻第一激光器2和第二激光器3输出频率为νa和νb0的两路激光各自对应的空气折射率na0、nb0；
[0067] 同时采用以下公式求解计算可获得两路激光对应的干涉条纹整数ma0，mb0：
[0068]
[0069] 其中，c为真空光速，L表示t0时刻的待测绝对距离。
[0070] 2)t0时刻后，解除飞秒光频梳1对第二激光器3的频率锁定，对第二激光器3的PZT调制输入端施加一个扫描电压，使第二激光器3的输出频率进行正向扫描，正向扫描具体是频率逐渐连续增加；
[0071] t1时刻，将第二激光器3重新锁定至飞秒光频梳1，锁定后第二激光器3输出激光的频率为νb1，νb1>νb0，由于第一激光器2一直被飞秒光频梳1锁定，第一激光器2输出激光的频率为νa，通过计算机19记录t0至t1时刻第二激光器3输出激光对应的干涉条纹整数变化Δm10、t1时刻第二激光器3输出激光对应的干涉条纹小数εb10≤εb1<1、t0至t1时刻第一激光器2输出激光对应的干涉条纹变化量ΔNa1以及t1时刻由环境参数测量模块12采集到的环境参数值，由环境参数值转换计算获得t1时刻第一激光器2和第二激光器3输出频率为νa和νb1的两路激光各自对应的空气折射率na1和nb1；
[0072] 然后采用以下两个公式联立进行求解获得FSI测距法对应的绝对距离粗测值L1和t0至t1时刻待测绝对距离的漂移量为ΔL1：
[0073]
[0074]
[0075] 其中，λS1为第二激光器3频率正向扫描所构建的首级合成波长，λS1＝c/(nb1νb1-nb0νb0)；
[0076] 3)t1时刻后，再次解除飞秒光频梳1对第二激光器3的频率锁定，通过控制第二激光器3内部的直流电机使第二激光器3的频率发生增加跳变；
[0077] t2时刻，将第二激光器3重新锁定至飞秒光频梳1，锁定后第二激光器3输出激光的频率为νb2，νb2>νb1，通过计算机19记录t2时刻第二激光器3输出激光对应的干涉条纹小数εb2、t0至t2时刻第一激光器2输出激光对应的干涉条纹变化量ΔNa2以及t2时刻由环境参数测量模块12采集到的环境参数值，由环境参数值转换计算获得t2时刻第一激光器2和第二激光器3输出频率为νa和νb2的两路激光各自对应的空气折射率na2和nb2，0≤εb2<1；
[0078] 然后采用以下公式求解获得t0至t2时刻待测绝对距离的漂移量为ΔL2：
[0079]
[0080] 具体实施中，na2-na0为一极小量。
[0081] 实施例的上述三个步骤过程中，绝对距离测量中第二激光器3频率变化情况，以及第一激光器2和第二激光器3对应的干涉条纹变化情况如图2所示。
[0082] 4)采用以下公式计算获得t0至t2时刻待测绝对距离漂移引起的第二激光器3输出激光的干涉条纹变化量ΔNb20：
[0083]
[0084] 再采用以下公式根据干涉条纹变化量ΔNb20对干涉条纹小数εb2进行修正：
[0085]
[0086] 其中，frac()表示求取小数部分，εb′2表示修正后的干涉条纹小数εb2；
[0087] 5)t0至t2时刻第二激光器3输出激光的频率从νb0跳变至νb2对应的二级合成波长为λS2，第二激光器3输出激光的频率在预先选择时满足多波长干涉的级间过渡条件u(L1)<λS2/4，u(L1)为频率扫描干涉对应的绝对距离测量不确定度，采用以下公式计算二级合成波长λS2对应的干涉条纹小数εS2：
[0088]
[0089] 其中，εS2表示二级合成波长λS2对应的干涉条纹小数；
[0090] 再采用以下公式计算二级合成波长λS2对应的干涉条纹整数mS2：
[0091]
[0092] 其中，mS2表示二级合成波长λS2对应的干涉条纹整数，floor()表示向下取整；
[0093] 最后采用以下公式的MWI测距法计算获得待测绝对距离的精测值L2：
[0094]
[0095] 其中，λS2表示t0至t2时刻第二激光器3输出激光的频率从νb0跳变至νb2对应的二级合成波长，λS2＝c/(nb2νb2-nb0νb0)。
[0096] 无论是FSI粗测，还是MWI精测，测量距离都可表示为如下的形式
[0097]
[0098] 其中， 为真空中的合成波长，n为空气折射率，m和ε分别为合成波长对应的干涉条纹整数和小数，测距不确定度为：
[0099]
[0100] 其中，λ为测量干涉仪对应的激光波长，波长值在780nm附近。
[0101] 本发明实施例中，测量距离为10m，空气折射率的测量不确定度为u(n)/n＝5×10-8，激光锁定至飞秒光频梳波长不确定度可达u(λ)/λ＝5×10-12，小数条纹测量不确定度u(ε)＝1×10-2。FSI粗测对应的频率扫描范围为85GHz，对应的首级合成波长为λS1＝3.53mm，代入至公式(11)得FSI粗测对应的测量不确定度为u(L1)＝5.3μm。选择二级合成波长大小为λS2＝60μm(由波长差为10nm两激光波长构建)，满足u(L1)<λS2/4的级间过渡条件。由公式(11)可知，激光波长不确定度对总的测距不确定度的影响可以忽略不计，绝对距离的综合测量不确定度可表示为：
[0102]
[0103] 由实施例可见，本发明将FSI与MWI相结合进行大范围的绝对距离测量，实现了10m测量范围内，亚微米级的测量不确定度。通过飞秒光频梳标定FSI中频率扫描的起点和终点以及MWI中每个波长的值，不仅将距离测量直接溯源至时间频率基准，还将绝对距离测量中由激光波长不确定度引起的测量不确定度降低至可以忽略的程度。
[0104] 构建基于EOM正弦相位调制的测量干涉仪和补偿干涉仪，利用PGC解调原理实现FSI和MWI过程中的实时相位解调，并通过补偿干涉仪实时监测和补偿由于环境参数改变和目标漂移引起的相位改变，光路结构简单、测量范围大、精度高，具有其突出显著的技术效果。
