-
本文提出的二次强度调制测距系统可以实现绝对距离的测量, 其利用马赫曾德尔强度调制器代替二次偏振调制测距中的电光相位调制器, 通过对信号的光强进行二次调制来进行测距. 相比于二次偏振调制测距, 二次强度调制测距无需考虑测距系统中的偏振态问题, 简化了系统结构, 提高了系统的稳定性. 经过相关理论推导以及实验验证: 二次强度调制测距系统的输出光强与调制频率成余弦关系, 并且可以直接测量调制器到目标物体之间的绝对距离, 系统的频率稳定度、相对测距精度皆达到10–7 量级. 本文提出的测距系统量程达到100 m, 相对测距精度稳定在10–7 量级. 采用摇摆法快速测距, 避免了直接扫频寻找光强极小值点对应的频率, 数据刷新率达到2 kHz. 二次强度调制测距系统测距速度快, 同时兼顾了较大的量程与较好的测距精度, 系统结构简单, 易于搭建, 具有广阔的应用前景.Long-range, high-precision, and high-refresh rate absolute distance measurement based on double intensity modulation is proposed and experimentally demonstrated. In this scheme, a Mach-Zehnder modulator is utilized to perform bidirectional modulation by a reversible optical path. In the Mach-Zehnder modulator, interference demodulation is completed by bidirectional modulated light with time difference. By adjusting the driving frequency of the modulator, the curve of light intensity versus driving frequency is achieved. Consequently, the distance to be measured can be obtained by extracted the frequency interval between two adjacent light intensity minimum points. In the traditional double polarization modulation ranging, the optical carrier is polarized by a polarizing beam splitter (PBS) before phase modulator. Moreover, a quarter wave or Faraday rotating mirror need to be placed to adjust the polarization in front of the target object. Therefore, the polarization state is an indispensable factor in the traditional double polarization modulation ranging. Due to the advantage of the intensity modulation, absolute distance measurement is achieved without additional polarization control, greatly simplifying the system. Theoretical analysis of the system is developed, which is then demonstrated by experiments. In the experiments, we achieved the following results. Firstly, the relationship between the intensity of the output light of the system and the modulation frequency is theoretically analyzed, which proved to be a cosine form. Secondly, swing method is introduced to realize high-speed absolute distance measurement during the analytical distance algorithm, and we achieved a refresh rate of 2 kHz in the experiments. Thirdly, the relationship between measurement ranging precision and frequency stability is analyzed. When the modulation frequency is 8.9 GHz, the experimentally measured frequency stability is on the order of 10–7. And when the distance to be measured is 2.73 m, the standard deviation of ranging reaches 1 μm, which represents the precision of the system. That is, the relative measurement precision is also on the order of 10–7, which is consistent with theoretical analysis. Finally, when the distance to be measured increases from 1.57 m to 100.83 m, the measurement precision increases from 1 μm to 30 μm. It is worth mentioning that the relative measurement precision of the system is always stable in the order of 10–7. Our scheme has some significant advantages, such as long-range, high-precision, high-refresh rate, and a simple and compact configuration. Moreover, our method can be used in important applications such as precision instruments, metrology, and aerospace.
1. 引 言
近年来, 高精度、大量程的快速测距在大型设备制造、精密仪器装配、航空航天、军事国防等科学工程领域的应用愈加广泛. 在我国的火星探测工程中, 祝融号火星车成功软着陆火星地表后, 需要其相对定位精度在距离10 m左右时优于3%[1]; 在大型客机的生产制造过程中, 各个部位工件加工精度与定位精度都需要达到毫米量级[2]; 在载人空间站的交会对接中, 目标相对测距机高速移动, 不仅需要高精度的绝对距离测量(毫米量级), 还需要足够的测距速率(3 Hz)以保证实时有效的跟踪控制[3].
目前, 绝对距离激光测距技术分为相干测量和非相干测量两大类. 相干测量以基于飞秒光频率梳的干涉测距为研究热点[4-10], 它在大量程的高精度测量中具有很大的应用潜力, 但由于飞秒光频率梳的高昂成本, 尚不能广泛应用于工业测量. 此外, 激光干涉也被应用在引力波探测中, 并取得了巨大的成功. 当引力波经过干涉仪的测量臂时, 通过干涉激光相位变化推导出测试质量间的距离变化, 从而反演得到引力波信息. 欧洲航天局(ESA)和美国航空航天局(NASA)合作的引力波探测计划LISA, 通过利用空间卫星构建空间激光干涉仪, 可以测量百万公里量级的皮米级距离变化[11]. 国内的太极计划和天琴计划也将激光干涉技术应用于引力波探测, 并分别于2019年8月和2019年12月发射了实验验证卫星. 经过实验测试, 他们构建的空间激光干涉仪也都可以实现皮米量级的超高精度距离测量[12,13]. 非相干测量中常用的技术手段包括脉冲测距法[14-16]、微波相位测距法[17-19]. 脉冲测距法通过测量激光往返待测路径的飞行时间来测得往返距离, 现阶段该方法主要用于超长距离测距, 测距精度一般在毫米甚至厘米量级[16]. 微波相位法测距是对测量光进行正弦调制, 通过比较激光往返待测距离的相位变化获取距离信息. 微波相位法测距是一种比较成熟的测距方案, 在几十厘米至几百米的测距范围内测距精度可以达到毫米甚至亚毫米量级[17,18], 并且可以实现0.03 s/次[19]的测距速度. 但要想进一步提升测距性能, 还需要进一步优化系统的鉴相模块.
二次偏振调制测距结合了相位法和干涉法的测距优点, 利用相位调制器对光信号进行两次偏振调制, 再经由干涉解调, 最终通过输出信号的相邻光强极值点对应的调制频率变化获取距离信息. 黑克非[20]、 肖洋[21]、高书苑等[22,23]都曾先后进行过相关研究, 并取得了一定的研究成果; 2022年, 高超等[24]尝试的波导式相位调制器测距系统, 利用光纤延时线测量得到的测距误差为0.39 mm. 二次偏振调制测距方法由传统的对时间、相位的测量转变为对光强的测量, 同时兼顾了较大的量程与较好的测距精度, 在工业领域中有很大的应用潜力. 但为了实现光的偏振调制, 需要用到偏振分光棱镜(PBS)进行起偏, 测量过程中偏振态的扰动也会在调制时引入额外的噪声, 影响测试精度. 此外, 测距过程中还需要在目标物体之前放置1/4波片或者法拉第旋光镜进行偏振态的调整, 这使得系统的复杂度增加, 限制了应用场景.
二次强度调制测距系统利用马赫曾德尔强度调制器代替二次偏振调制测距系统中的相位调制器, 与二次偏振调制测距一样, 二次强度调制测距也是通过解调后的光强反演得出距离信息, 理论上可以达到和二次偏振调制相同的测距效果. 而相比二次偏振调制测距, 二次强度调制测距通过对光的强度调制代替偏振调制, 激光器输出的光信号不用PBS进行起偏, 可以直接输入到强度调制器中进行调制, 两次调制间的光信号也无需进行偏振态的调整. 整个系统的结构可得到充分简化, 能适应更多不同的测量场景, 更有利于后续的工业化设计.
本文对二次强度调制测距原理进行了理论分析, 推导证明了系统输出光强度和调制频率满足余弦函数关系, 并进行了实验验证. 实验结果表明: 二次强度调制测距系统可以实现调制器到目标物体之间的绝对距离测量, 系统的频率稳定度以及相对测距精度都在10–7量级, 测距量程达到100 m. 此外, 本文提出通过摇摆法间接找寻光强极小值点对应的频率, 不仅提高了测距数据刷新率(2 kHz), 还降低了噪声的影响.
2. 二次强度调制测距系统
2.1 测距原理
图 1是本文提出的二次强度调制测距原理图, 激光器输出信号由环形器1端口输入, 2端口输出进入马赫-曾德尔强度调制器, 假设激光的光场强度为
$ E_{0} $ , 频率为$ \omega_{0} $ , 初始相位为$ \phi _{0} $ , 则进入调制器的光信号可表示为$$ \begin{array}{l} E_{\rm in} = E_{0}{\rm exp}[{\rm j}(\omega_{0}t + \phi _{0})]. \end{array} $$ (1) 马赫-曾德尔强度调制器的调制原理是将输入光分为振幅、相位相等的两束光. 通过对其中一束光施加外部电场进行相位调制, 使上下两束光信号产生相位差. 两束光在输出端干涉叠加时的光强会随着两束光信号的相位差改变, 从而实现外部电场对光信号的强度调制[25]. 定义调制器的半波电压参数为
$ V_{\pi} $ , 表示两路信号相位差为$ \pi $ 时引入的外电场. 则当外部电场的直流电压为$ V_{\rm dc} $ , 射频部分为$ V_{\rm ac}{\rm cos}(\omega_{\rm ac}t) $ 时, 经过马赫-曾德尔调制器后的输出光信号可表示为[26]$$ E_{1} = \frac{E_{0}}{2}{\rm exp}[{\rm j}(\omega_{0}t + \phi _{0})][{\rm exp}({\rm j}\varphi_{1}) + 1], \ $$ (2) $$ \begin{array}{l} \varphi_{1} = \varphi_{\rm bias} + m{\rm cos}(\omega_{\rm ac}t), \end{array} $$ (3) 式中,
$ \varphi_{\rm bias} = {\pi V_{\rm dc}}/{V_{\pi}} $ 表示调制器的外部直流偏置,$ m = {\pi V_{\rm ac}}/{V_{\pi}} $ 表示调制器的调制深度.调制器第一次调制后的光信号经过一段待测距离后到达反射镜反射, 光信号沿原路返回, 返回的光信号再次经过调制器进行第2次调制:
$$ \begin{split} E_{2} =\;& \frac{E_{0}}{4}{\rm exp[j}(\omega_{0}t + \phi _{0})][{\rm exp(j}\varphi_{1}) + 1]\\ &\times[{\rm exp(j}\varphi_{2}) + 1], \end{split} $$ (4) $$ \begin{array}{l} \varphi_{2} = \varphi_{\rm bias} + m'{\mathrm{cos}}[\omega_{\rm ac}(t + \Delta t)]. \end{array} $$ (5) 第2次调制时, 由于光信号从调制器的输出端进入, 调制器的调制深度与第1次有所不同[27]:
$ m'=\dfrac{\pi V\mathrm{_{ac}}}{V_{\pi}'} $ ,$ V_{\pi}'=V_{\pi}\dfrac{\omega\mathrm{_{ac}}\tau}{\mathrm{sin}(\omega_{\mathrm{ac}}\tau)} $ , 其中$ \tau $ 为马赫-曾德尔强度调制器中信号的渡越时间.可得信号光强为
$$ \begin{split} \;&P_{{\mathrm{out}}} = E_{2}\cdot E_{2}^{*} \\ = \;&\frac{E_{0}^{2}}{16}{\rm exp[j}(\omega_{0}t + \phi_{0})] [{\rm exp(j}\varphi_{1}) + 1][{\rm exp(j}\varphi_{2}) + 1] \\ &\times {\rm exp[j}( - \omega_{0}t - \phi _{0})][{\rm exp(-j}\varphi_{1}) + 1 ]\\ &\times[{\rm exp(-j}\varphi_{2}) + 1] \\ = \;&\frac{E_{0}^{2}}{16}[2{\rm cos}(\varphi_{1} + \varphi_{2}) + 2{\rm cos}(\varphi_{1} - \varphi_{2}) \\ &+ 4{\rm cos}(\varphi_{1}) + 4{\rm cos}(\varphi_{2}) + 4]\\ = \;&\frac{E_{0}^{2}}{4}[{\rm cos}(\varphi_{1}){\rm cos}(\varphi_{2}) + {\rm cos}(\varphi_{1}) + {\rm cos}(\varphi_{2}) + 1].\\[-1pt] \end{split} $$ (6) 马赫-曾德尔强度调制器根据不同的直流偏置, 会在不同状态工作, 图 2为调制器在不同直流偏置下的传输曲线, 其中有3种典型的工作状态: 最大输出点、最小输出点和正交传输点. 在本实验中, 调节直流偏置电压
$ V\mathrm{_{dc}}= {3V_{\pi}}/{2} $ , 使其工作在斜率为正的正交传输点, 此时理论上调制器的输出光信号功率$ P_{1} $ 的变化频率等于调制的射频信号$ U_{\mathrm{ac}} $ 的频率, 推导如下:$$ \begin{split} &{\rm cos}(\varphi_{1}) = {\rm cos}\Big[\frac{3\pi}{2} + m{\rm cos}(\omega_{\rm ac}t)\Big] = {\rm sin} \Big[m{\rm cos}(\omega_{\rm ac}t)\Big]\\ =\;& -2\sum\limits_{n = 1}^{\infty}( - 1)^{n}\mathrm{J}_{2n - 1}(m) {\rm cos}[(2n - 1)\omega_{\rm ac}t]. \\[-1pt] \end{split} $$ (7) 图 3为
$ \rm{J}_{2 n - 1}(m) $ 表示的第一类贝塞尔函数, 实验中选取的强度调制器的半波电压$ V_{\pi} $ 为6 V, 射频信号幅度$ V_{\rm ac} $ 为600 mV, 调制深度m约为0.31, 此时贝塞尔函数一阶项系数为0.15, 三阶项系数为0.0006. 贝塞尔函数三阶及以上的高阶项远小于一阶项, 因此可以忽略:$$ \begin{array}{l} {\rm cos}(\varphi_{1}) = 2\mathrm{J}_{1}(m){\rm cos}(\omega_{\rm ac}t). \end{array} $$ (8) 同样, 可以得到:
$$ \begin{array}{l} {\rm cos}(\varphi_{2}) = 2\mathrm{J}_{1}(m'){\rm cos}[\omega_{\rm ac}(t + \Delta t)]. \end{array} $$ (9) $$ \begin{split} P_{\rm out} =\;& {E_{0}^{2}} \{4\mathrm{J}_{1}(m){\rm cos}(\omega_{\rm ac}t)\mathrm{J}_{1}(m'){\rm cos}[\omega_{\rm ac}(t + \Delta t)] \\ &+ 2\mathrm{J}_{1}(m){\rm cos}(\omega_{\rm ac}t) \\ &+ 2\mathrm{J}_{1}(m'){\rm cos}[\omega_{\rm ac}(t + \Delta t)] + 1\}/4\\ =\;& {E_{0}^{2}} \{2\mathrm{J}_{1}(m){\mathrm{J}}_{1}(m'){\rm cos}(\omega_{\rm ac}\Delta t)\\ &+ 2\mathrm{J}_{1}(m)\mathrm{J}_{1}(m'){\rm cos}[\omega_{\rm ac}(2t + \Delta t)]\\ &+ 2\mathrm{J}_{1}(m){\rm cos}(\omega_{\rm ac}t) \\ &+ 2\mathrm{J}_{1}(m'){\rm cos}[\omega_{\rm ac}(t + \Delta t)] + 1\}/4.\\[-1pt] \end{split} $$ (10) 两次调制后的光信号再次经过环形器, 由2端口输入, 3端口输出进入光电探测器(PD)进行光电转换. 当待测距离固定时, 探测光的往返时间
$ \Delta t $ 是固定不变的, 包含$ \mathrm{cos}(\omega_{\rm ac}\Delta t) $ 的部分为直流分量, 而包含$ \mathrm{cos}(\omega\mathrm{_{ac}}t) $ ,$ \mathrm{cos}(2\omega_{\mathrm{ac}}t) $ 的部分为交流分量. 由于PD的响应带宽远小于调制信号的频率, 因此信号中除了直流分量以外, 其他交流分量都被滤除, 最终探测的信号光强$ P\mathrm{_{out}} $ 为$$ P_{\rm out} = \frac{E_{0}^{2}}{4}[1 + 2{\mathrm{J}}_{1}(m){\mathrm{J}}_{1}(m'){\rm cos}(\omega_{\rm ac}\Delta t)]. $$ (11) 2.2 连续扫频测距
(11)式表示PD探测光强与调制信号角频率
$ \omega\mathrm{_{ac}} $ 以及信号飞行时间$ \Delta t $ 的乘积满足余弦函数关系.$ \Delta t $ 可用待测距离D、真空中光速c以及空气的等效折射率n表示为$$ \begin{array}{l} \Delta t = \dfrac{2Dn}{c}. \end{array} $$ (12) 将其代入(11)式, 并将角频率
$ \omega_{\mathrm{ac}} $ 替换为$ 2\pi f\mathrm{_{ac}} $ :$$ P_{\rm out} = \frac{E_{0}^{2}}{4}\left[1 + 2{\mathrm{J}}_{1}(m){\mathrm{J}}_{1}(m') {\rm cos}\left(\frac{4\pi f_{\rm ac}Dn}{c}\right)\right]. $$ (13) (13)式表明, 当待测距离D确定, 在
$ {\mathrm{J}}_{1}(m) $ 和$ {\mathrm{J}}_{1}(m') $ 不变的前提下, PD探测到的光功率随着$ f\mathrm{_{ac}} $ 的变化呈余弦变化. 通过连续扫频的方法, 找到相邻的两个光强极小值点对应的频率$ f_{1} $ ,$ f_{2} $ , 分别对应余弦曲线的$ (2 N - 1)\pi $ ,$ (2 N + 1)\pi $ , 则有:$$ \frac{4\pi f_{1}Dn }{c} = (2N - 1)\pi, $$ (14) $$ \frac{4\pi f_{2}Dn }{c} = (2N + 1)\pi. $$ (15) 联立两个方程可以得到待测距离D的表达式:
$$ D=(2N-1)\frac{c}{4nf_1}=\left[\frac{2f_1}{f_2-f_1}\right]\frac{c}{4nf_1}, $$ (16) 式中[·]表示取整运算. 通过连续扫频测距方法, 只需找到相邻两个光强极小值点对应的频率
$ f_{1} $ 和$ f_{2} $ 即可求出最终的待测距离, 无需考虑传统相位法测距中模糊距离的限制.3. 实验系统搭建
3.1 测距系统可行性分析
根据上述二次强度调制的实验原理, 搭建实验测距装置如图 4所示. 光学结构部分由激光器, 环形器(Cir, 中心波长1550 nm, 插入损耗0.7 dB), 可调光衰减器(VOA, 衰减范围0.8—60 dB), 准直镜(中心波长1550 nm, 出射损耗0.5 dB)和反射镜(插入损耗0.6 dB)构成. 激光器使用半导体激光器, 中心波长为1550 nm, 为应对远距离测量要求, 其输出功率为14 dBm. 因为最后需要根据光强信号得出距离信息, 因此激光器的输出功率要尽可能的稳定, 激光器的短期功率波动小于0.01 dB. 激光器产生的光信号由环形器1 端口输入, 2 端口输出进入到马赫-曾德尔调制器进行第1次调制, 直流电压源提供调制器的直流偏置, 射频源提供调制器的交流驱动. 调制器输出的光信号经过准直镜结构进入测量空间并经过一段距离后被反射镜反射沿原路返回, 再次进入调制器进行第2次调制. 两次调制后的信号由环形器2 端口输入, 3 端口输出, 最后经由可调光衰减器进行功率调整. 光接收系统由光电探测器(PD, 可承受光功率18 μW, 响应度0.9 A/W)和模数转换(A/D, 转换位数16 bit)组成. PD 先将光信号转换为模拟电信号, A/D再将模拟电信号转换为数字信号. 控制与处理系统包括电脑接收处理A/D转换后的数据以及利用STM32微处理器控制射频源的射频输出两部分. 实验开始时需要保证调制器工作在正交传输点, 设置方式是改变直流偏置电压, 确定调制器输出的最大功率. 在最大功率的基础上, 再次改变直流偏置电压, 最大功率减3 dB就是调制器工作在正交传输点时的输出. 根据(16)式可知, 待测距离D还与空气的等效折射率n有关, 空气扰动会直接影响测量精度, 因此要求外界环境的温度、湿度以及大气压强等要尽可能的稳定. 为了减少由于空气扰动造成的误差, 实验测试时使用玻璃罩对系统的光学结构部分进行封闭.
系统搭建完成后, 对系统的测距可行性进行实验验证, 控制射频源以100 kHz为步进, 进行100—400 MHz范围内的扫频, 再由电脑对AD转换的数字信号进行分析处理. 根据AD转换器测量得到的光强信息绘制扫频曲线如图 5所示, PD的输出电压与调制频率满足三角函数关系, 这与(13)式的推导一致.
3.2 摇摆法快速测距
根据(16)式可知, 二次强度调制测距的关键是找到两个相邻的光强极小值点对应的频率. 如果利用直接测量的方法寻找, 需要从某个频率点开始进行多次频率扫描, 极大地延长了测距时间. 为了尽可能地缩短测距时间, 提高数据刷新率, 本文提出利用摇摆法进行快速测距.
根据(13)式以及图 5可知, 探测光强和调制频率成余弦关系. 测量时选取中心频率
$ f $ 和偏移频率$ \Delta f $ , 由于余弦函数的对称性, 可以根据调制频率分别为$ f + \Delta f $ 和$ f - \Delta f $ 的系统输出来判断中心频率$ f $ 的位置. 具体用频率为$ f + \Delta f $ 时的系统输出减去频率为$ f - \Delta f $ 时的系统输出, 这个差值记作摇摆差值$ \Delta U $ , 则当中心频率$ f $ 位于不同位置时$ \Delta U $ 的变化情况如图 6所示.$ \Delta U$ 的表达式为图 6 不同位置下f与 的关系 (a) f在光强极小值点对应频率的左侧; (b) f为光强极小值点对应的频率; (c) f在光强极小值点对应频率的右侧$ \Delta U $ Fig. 6. Relationship between f and at different positions: (a) f is on the left side of the frequency corresponding to the light intensity minimum point; (b) f is the frequency corresponding to the minimum point of light intensity; (c) f is on the right side of the frequency corresponding to the light intensity minimum point.$ \Delta U $ $$ \begin{split} & \Delta U \\=\;& \frac{E_{0}^{2}}{4} \left\{ 1 + 2{\mathrm{J}}_{1}(m){\mathrm{J}}_{1}(m'){\rm cos}\left[\frac{4\pi (f + \Delta f)Dn}{c}\right] \right\}\\ &- \frac{E_{0}^{2}}{4}\left\{1 + 2{\mathrm{J}}_{1}(m){\mathrm{J}}_{1}(m'){\rm cos}\left[\frac{4\pi (f - \Delta f)Dn}{c}\right]\right\}\\ = \;&- E_{0}^{2}{\mathrm{J}}_{1}(m){\mathrm{J}}_{1}(m'){\rm sin}\left(\frac{4\pi\Delta fDn}{c}\right){\rm sin}\left(\frac{4\pi fDn}{c}\right). \end{split} $$ (17) 根据上述摇摆差值的计算方法, 当偏移频率
$ \Delta f $ 固定, 摇摆差值$ \Delta U $ 与中心频率$ f $ 满足正弦函数关系. 随着中心频率$ f $ 从扫频曲线极小值点的左侧移动到右侧, 摇摆差值$ \Delta U $ 的值从小于0变化到大于0, 具体关系曲线如图 7(a)所示, 扫频曲线的极小值点对应的调制频率$ f_{0} $ 就是摇摆差值曲线的零点对应的调制频率. 图 7(a)的摇摆差值曲线是一个正弦曲线, 在其零点附近局部放大, 可以近似为一条直线(如图 7(b)所示), 选取中心频率$ f_{1} $ 与$ f_{2} $ 分别位于摇摆差值曲线零点对应调制频率$ f_{0} $ 的左右两侧, 相应的摇摆差值$ \Delta U_{1} $ 小于0,$ \Delta U_{2} $ 大于0. 通过$ f_{1} $ ,$ f_{2} $ ,$ \Delta U_{1} $ 以及$ \Delta U_{2} $ 的大小, 就可以拟合得到摇摆差值曲线零点对应的调制频率$ f_{0} $ 的具体数值:$$ \frac{ f_{0} - f_{1} }{ f_{2} - f_{0}} = \frac{0 - \Delta U_{1} }{\Delta U_{2} - 0}, $$ (18) $$ f_{0} = f_{1} + \frac{ f_{2} - f_{1}}{1 - \dfrac{\Delta U_{2}}{\Delta U_{1}}}. $$ (19) 根据(19)式可知, 找到摇摆差值曲线零点对应频率两侧的两个频率值, 再测量其相应的摇摆差值, 即可得到扫频曲线极小值点对应的调制频率. 具体方法: 设置合适的扫频间隔, 根据
$ \Delta U $ 的正负判断当前频率与摇摆差值曲线零点对应频率的相对位置, 最终找到4组频率, 确保它们包括两个相邻的摇摆差值曲线零点对应频率. 通过这4组频率及其摇摆差值计算得到扫频曲线连续两个光强极小值点对应的频率, 将这两个频率值代入(16)式即可求出待测距离. 扫频间隔的设置有两点原则, 一是要尽可能的大, 以减少扫频次数; 二是待测距离大致确定时, 摇摆差值曲线的周期也就大致确定, 要使扫频间隔小于摇摆差值曲线周期的一半, 确保每次扫频时不会漏掉摇摆差值曲线的零点.相比于通过直接测量的方式找寻光强极小值点对应的频率, 摇摆法所需要的跳频次数明显减少, 这无疑极大地减少了测距时间. 在当前的实验条件下, 通过摇摆法快速测距, 最高可以实现2 kHz的数据刷新率. 除了快速测距, 摇摆法的另一个优势是把光强极小值点的系统输出测量转变为对光强极小值点附近的输出测量. 而相比于系统在光强极小值点时的输出, 附近点的输出响应更大, 受噪声的影响相对更小.
4. 分析及讨论
4.1 绝对距离测量分析
由上述分析可知, 根据系统的扫频曲线, 结合摇摆法, 可以实现调制器到待测目标之间的绝对距离测量. 为了验证绝对距离测量的测距精度及其影响因素, 实验选取待测距离2.73 m, 对这段距离重复测量100 次, 得到的100 组数据整合为图 8所示的折线.
根据图 8的测距结果可知: 100次重复测距的结果均稳定在2.73 m, 且波动不超过5 μm. 对这些数据进一步进行分析计算, 得到距离结果的标准差为1 μm, 系统的相对测距精度为
$ 3.66 \times 10^{ - 7} $ .根据(16)式, 测距结果与相邻光强极小值点对应的频率值有关, 因此频率的稳定与否决定了系统的测距能力. 为了进一步探究频率稳定度与测距精度的关系, (16)式对频率求导得
$$ \begin{split} D &= \left\{ \left[ \dfrac{2 f_{1}}{ f_{2}- f_{1}}\right]\frac{c}{4n f_{1}}\right\}'\\ & = -\left[\dfrac{2 f_{1}}{ f_{2}- f_{1}}\right]\frac{c}{4n f^{2}_{1}}\Delta f_{1}. \end{split} $$ (20) $$ \Delta D = - \frac{D}{ f_{1}}\Delta f_{1}, $$ (21) $$ \frac{\Delta D}{D} = - \frac{\Delta f_{1}}{ f_{1}}. $$ (22) 可以发现, 系统相对测距精度的大小等于第一个频率极值点f1的相对稳定度大小. 为了进一步验证它们的关系, 实验进行3组重复距离测试, 每组进行30次扫频测距, 根据扫频曲线算出两个相邻光强极小值点对应的频率值. 分析第一个频率极值点f1的平均值与标准差, 得到的实验结果见表 1.
表 1 f1的相对稳定度Table 1. Relative stability of f1.Group Average of $ f_{1}$/Hz Standard deviation of $ f_{1}$/Hz Relative accuracy/10–7 Ⅰ 8965550359 2737 3.05 Ⅱ 8965546127 2760 3.07 Ⅲ 8965538910 3736 4.17 根据表1的结果可以发现,
$ f_{1} $ 的相对稳定度在10–7量级, 与测试得到的相对测距精度大致相等, 这与理论分析相一致. 说明二次强度调制测距系统可以实现绝对距离的测量, 相对测距精度与频率的相对稳定度有关, 可以达到10–7量级.4.2 测距量程分析
(16)式也可写成如下形式:
$$ \begin{split} D =\;& \left [\dfrac{2 f_{1}}{ f_{2} - f_{1}}\right] \frac{c}{4n f_{1}} \approx \dfrac{2 f_{1}}{ f_{2} - f_{1}} \times \frac{c}{4n f_{1}} \\ =\;& \frac{c}{2n( f_{2} - f_{1})}. \end{split} $$ (23) 在光速c以及空气的等效折射率n确定的情况下, 随着待测距离的增大, 两个光强极小值点对应的频率间隔将减小. 根据(23)式, 当待测距离从1 m变化到100 m, 两个光强极小值点对应的频率间隔从150 MHz变化到1.5 MHz. 因此为了满足1—100 m范围内的测距, 所需要的扫频范围需涵盖300—3 MHz的变化区间, 这个扫频范围在目前的实验条件下是可以实现的, 这意味着本系统的测距量程可以达到100 m.
为了验证系统在不同待测距离下测距精度的变化, 选取几段固定长度的光纤进行重复性测距. 我们分别选取等效空间距离为1.57, 9.14, 23.31, 38.79, 54.11, 64.20, 79.67和100.83 m长的8段光纤, 对每段光纤进行30次重复距离测试, 分析测距结果的平均值以及标准差, 计算相对测距精度, 实验结果见表 2.
表 2 相对测距精度Table 2. Relative distance measurement accuracy.Group Average of D/m Standard deviation of D/μm Relative accuracy/10–7 Ⅰ 1.575537 1 6.35 Ⅱ 9.139261 5 5.47 Ⅲ 23.313040 8 3.43 Ⅳ 38.786468 10 2.57 Ⅴ 54.107275 15 2.77 Ⅵ 64.199287 16 2.55 Ⅶ 79.673707 26 3.49 Ⅷ 100.830172 30 2.98 由表2可见, 随着待测距离从1.57 m增至100.83 m, 测距结果的标准差从1 μm增至30 μm, 但其相对测距精度始终在10–7量级. 说明在100 m量程范围内, 该测距系统的相对测距精度能够达到10–7量级.
5. 结 论
本文提出了二次强度调制系统, 利用马赫-曾德尔强度调制器对光信号的二次调制来获取距离与光强的信息. 我们搭建了实验系统, 获得了探测光强与调制频率的余弦函数关系曲线, 与理论推导相一致, 并通过摇摆法间接找寻光强极小值点对应的频率, 不仅降低了噪声的影响, 还提高了测距速度, 实现2 kHz的数据刷新率. 二次强度调制系统可以直接测量调制器到目标物体之间的绝对距离, 测得的频率稳定度以及相对测距精度都在10–7量级, 与理论分析结果一致. 实验可实现1—100 m量程内的测距, 且相对测量精度始终在10–7 量级. 综上, 二次强度调制测距系统可以实现绝对距离的测量, 不仅在测距精度以及测距量程上保留了二次偏振调制测距的优势, 还省略掉了很多偏振元器件, 系统结构更加简化, 应用场景适应性更强.
[1] 王镓, 李达飞, 何锡明, 成子青, 许倩, 钱雪茹, 万文辉 2022 深空探测学报 9 62
Wang J, Li D F, He X M, Cheng Z Q, Xu Q, Qian X R, Wan W H 2022 J. Deep Space Explor. 9 62
[2] 于勇, 陶剑, 范玉青 2009 航空制造技术 11 56
Google Scholar
Yu Y, Tao J, Fan Y Q 2009 Aeronaut. Manuf. Technol. 11 56
Google Scholar
[3] 邵珠法 2005 硕士学位论文 (成都: 电子科技大学)
Shao Z F 2005 M. S. Thesis (Chengdu: University of Electronic Science and Technology of China
[4] Minoshima K, Matsumoto H 2000 Appl. Opt. 39 5512
Google Scholar
[5] Ye J 2004 Opt. Lett. 29 1153
Google Scholar
[6] Cui M, Zeitouny M G, Bhattacharya N, van den Berg S A, Urbach H P, Braat J J M 2009 Opt. Lett. 34 1982
Google Scholar
[7] Lee J, Kim Y J, Lee K, Lee S, Kim S W 2010 Nat. Photonics 4 716
Google Scholar
[8] Minoshima K, Arai K, Inaba H 2011 Opt. Express 19 26095
Google Scholar
[9] Lee J, Han S, Lee K, Bae E, Kim S, Lee S, Kim S W, Kim Y J 2013 Meas. Sci. Technol. 24 045201
Google Scholar
[10] Wu G H, Liao L, Xiong S L, Li G Y, Cai Z J, Zhu Z 2018 Sci. Rep. 8 4362
Google Scholar
[11] Armano M, Audley H, Baird J, et al. 2018 Phys. Rev. Lett. 120 061101
Google Scholar
[12] Gao R H, Liu H S, Luo Z R, Jin G 2019 Chin. Opt. 12 425
Google Scholar
[13] Gong Y G, Luo J, Wang B 2021 Nat. Astron. 5 881
Google Scholar
[14] Nissinen J, Nissinen I, Kostamovaara J 2009 IEEE J. Solid-State Circuits 44 1486
Google Scholar
[15] Wang T S, Huo J, Wang P C, Dong P, Yang R, Li M F 2017 Proceedings of the 2017 Symposium on Quantum Information Technology and Applications Beijing, China, August 6–15, 2017 p50
[16] Kilpela A, Pennala R, Kostamovaara J 2001 Rev. Sci. Instrum. 72 2197
Google Scholar
[17] 刘玉周 2015 硕士学位论文 (武汉: 华中科技大学)
Liu Y Z 2015 M. S. Thesis (Wuhan: Huazhong University of Science and Technology
[18] 许贤泽, 翁名杰, 徐逢秋, 白翔 2017 光学精密工程 8 1979
Xu X Z, Weng M J, Xu F Q, Bai X 2017 Opt. Precis. Eng. 8 1979
[19] 郑大青, 陈伟民, 陈丽, 李存龙 2015 光电子·激光 26 303
Zheng D Q, Chen W M, Chen L, Li C L 2015 J. Optoelectron. Laser 26 303
[20] 黑克非, 于晋龙, 王菊, 王文睿, 贾石, 吴穹, 薛纪强 2014 物理学报 63 100602
Google Scholar
Hei K F, Yu J L, Wang J, Wang W R, Jia S, Wu Q, Xue J Q 2014 Acta Phys. Sin. 63 100602
Google Scholar
[21] 肖洋, 于晋龙, 王菊, 王文睿, 王子雄, 谢田元 2016 物理学报 65 100601
Google Scholar
Xiao Y, Yu J L, Wang J, Wang W R, Wang Z X, Xie T Y 2016 Acta Phys. Sin. 65 100601
Google Scholar
[22] 高书苑, 石俊凯, 纪荣祎, 黎尧, 周维虎 2018 中国激光 45 4005
Gao S Y, Shi J K, Ji R W, Li Y, Zhou W H 2018 Chin. J. Lasers 45 4005
[23] 高书苑, 黎尧, 纪荣祎, 石俊凯, 胡哲文, 周维虎 2019 光学精密工程 27 279
Google Scholar
Gao S Y, Li Y, Ji R W, Shi J K, Hu Z W, Zhou W H 2019 Opt. Precis. Eng. 27 279
Google Scholar
[24] 高超, 纪荣祎, 高书苑, 董登峰, 周维虎 2022 光学精密工程 30 246
Gao C, Ji R Y, Gao S Y, Dong D F, Zhou W H 2022 Opt. Precis. Eng. 30 246
[25] 亢洋 2021 硕士学位论文 (北京: 北京邮电大学)
Kang Y 2021 M. S. Thesis (Beijing: Beijing University of Posts and Telecommunications
[26] 张伟婷 2020 硕士学位论文 (成都: 电子科技大学)
Zhang W T 2020 M. S. Thesis (Chengdu: University of Electronic Science and Technology of China
[27] Frankel M Y, Esman R D 1998 J. Lightwave Technol. 16 859
Google Scholar
-
图 6 不同位置下f与
$ \Delta U $ 的关系 (a) f在光强极小值点对应频率的左侧; (b) f为光强极小值点对应的频率; (c) f在光强极小值点对应频率的右侧Fig. 6. Relationship between f and
$ \Delta U $ at different positions: (a) f is on the left side of the frequency corresponding to the light intensity minimum point; (b) f is the frequency corresponding to the minimum point of light intensity; (c) f is on the right side of the frequency corresponding to the light intensity minimum point.表 1 f1的相对稳定度
Table 1. Relative stability of f1.
Group Average of $ f_{1}$/Hz Standard deviation of $ f_{1}$/Hz Relative accuracy/10–7 Ⅰ 8965550359 2737 3.05 Ⅱ 8965546127 2760 3.07 Ⅲ 8965538910 3736 4.17 表 2 相对测距精度
Table 2. Relative distance measurement accuracy.
Group Average of D/m Standard deviation of D/μm Relative accuracy/10–7 Ⅰ 1.575537 1 6.35 Ⅱ 9.139261 5 5.47 Ⅲ 23.313040 8 3.43 Ⅳ 38.786468 10 2.57 Ⅴ 54.107275 15 2.77 Ⅵ 64.199287 16 2.55 Ⅶ 79.673707 26 3.49 Ⅷ 100.830172 30 2.98 -
[1] 王镓, 李达飞, 何锡明, 成子青, 许倩, 钱雪茹, 万文辉 2022 深空探测学报 9 62
Wang J, Li D F, He X M, Cheng Z Q, Xu Q, Qian X R, Wan W H 2022 J. Deep Space Explor. 9 62
[2] 于勇, 陶剑, 范玉青 2009 航空制造技术 11 56
Google Scholar
Yu Y, Tao J, Fan Y Q 2009 Aeronaut. Manuf. Technol. 11 56
Google Scholar
[3] 邵珠法 2005 硕士学位论文 (成都: 电子科技大学)
Shao Z F 2005 M. S. Thesis (Chengdu: University of Electronic Science and Technology of China
[4] Minoshima K, Matsumoto H 2000 Appl. Opt. 39 5512
Google Scholar
[5] Ye J 2004 Opt. Lett. 29 1153
Google Scholar
[6] Cui M, Zeitouny M G, Bhattacharya N, van den Berg S A, Urbach H P, Braat J J M 2009 Opt. Lett. 34 1982
Google Scholar
[7] Lee J, Kim Y J, Lee K, Lee S, Kim S W 2010 Nat. Photonics 4 716
Google Scholar
[8] Minoshima K, Arai K, Inaba H 2011 Opt. Express 19 26095
Google Scholar
[9] Lee J, Han S, Lee K, Bae E, Kim S, Lee S, Kim S W, Kim Y J 2013 Meas. Sci. Technol. 24 045201
Google Scholar
[10] Wu G H, Liao L, Xiong S L, Li G Y, Cai Z J, Zhu Z 2018 Sci. Rep. 8 4362
Google Scholar
[11] Armano M, Audley H, Baird J, et al. 2018 Phys. Rev. Lett. 120 061101
Google Scholar
[12] Gao R H, Liu H S, Luo Z R, Jin G 2019 Chin. Opt. 12 425
Google Scholar
[13] Gong Y G, Luo J, Wang B 2021 Nat. Astron. 5 881
Google Scholar
[14] Nissinen J, Nissinen I, Kostamovaara J 2009 IEEE J. Solid-State Circuits 44 1486
Google Scholar
[15] Wang T S, Huo J, Wang P C, Dong P, Yang R, Li M F 2017 Proceedings of the 2017 Symposium on Quantum Information Technology and Applications Beijing, China, August 6–15, 2017 p50
[16] Kilpela A, Pennala R, Kostamovaara J 2001 Rev. Sci. Instrum. 72 2197
Google Scholar
[17] 刘玉周 2015 硕士学位论文 (武汉: 华中科技大学)
Liu Y Z 2015 M. S. Thesis (Wuhan: Huazhong University of Science and Technology
[18] 许贤泽, 翁名杰, 徐逢秋, 白翔 2017 光学精密工程 8 1979
Xu X Z, Weng M J, Xu F Q, Bai X 2017 Opt. Precis. Eng. 8 1979
[19] 郑大青, 陈伟民, 陈丽, 李存龙 2015 光电子·激光 26 303
Zheng D Q, Chen W M, Chen L, Li C L 2015 J. Optoelectron. Laser 26 303
[20] 黑克非, 于晋龙, 王菊, 王文睿, 贾石, 吴穹, 薛纪强 2014 物理学报 63 100602
Google Scholar
Hei K F, Yu J L, Wang J, Wang W R, Jia S, Wu Q, Xue J Q 2014 Acta Phys. Sin. 63 100602
Google Scholar
[21] 肖洋, 于晋龙, 王菊, 王文睿, 王子雄, 谢田元 2016 物理学报 65 100601
Google Scholar
Xiao Y, Yu J L, Wang J, Wang W R, Wang Z X, Xie T Y 2016 Acta Phys. Sin. 65 100601
Google Scholar
[22] 高书苑, 石俊凯, 纪荣祎, 黎尧, 周维虎 2018 中国激光 45 4005
Gao S Y, Shi J K, Ji R W, Li Y, Zhou W H 2018 Chin. J. Lasers 45 4005
[23] 高书苑, 黎尧, 纪荣祎, 石俊凯, 胡哲文, 周维虎 2019 光学精密工程 27 279
Google Scholar
Gao S Y, Li Y, Ji R W, Shi J K, Hu Z W, Zhou W H 2019 Opt. Precis. Eng. 27 279
Google Scholar
[24] 高超, 纪荣祎, 高书苑, 董登峰, 周维虎 2022 光学精密工程 30 246
Gao C, Ji R Y, Gao S Y, Dong D F, Zhou W H 2022 Opt. Precis. Eng. 30 246
[25] 亢洋 2021 硕士学位论文 (北京: 北京邮电大学)
Kang Y 2021 M. S. Thesis (Beijing: Beijing University of Posts and Telecommunications
[26] 张伟婷 2020 硕士学位论文 (成都: 电子科技大学)
Zhang W T 2020 M. S. Thesis (Chengdu: University of Electronic Science and Technology of China
[27] Frankel M Y, Esman R D 1998 J. Lightwave Technol. 16 859
Google Scholar
计量
- 文章访问数: 3140
- PDF下载量: 193
- 被引次数: 0