韩玉祥,丛至诚,高丙坤,崔翔宇
(1.东北石油大学计算机与信息技术学院,黑龙江 大庆 163318;2.东北石油大学电气信息工程学院,黑龙江 大庆 163318)
激光自混合干涉系统结构简单、装置紧凑、鲁棒性好,能够实现物体位移、距离、振动、速度、表面轮廓、透明物体厚度等高精度测量[1-12]。在干扰较小的环境中,大都采用激光器发射的调幅(AM)信号对目标进行测量,在满足高精度测量要求的同时,结构更加轻量化,产品也更加实用化。但在大噪声环境中,待测物体反射到激光自混合干涉系统中产生的相干调幅信号难以检测,为此通过调频(FM)信号测量目标物体状态的方法陆续提出[13-16]。微位移测量已大量应用于大型土木结构、航空航天、健康监测等,故本文在激光自混合干涉平台上,采用马赫-曾德尔干涉仪,进行了调幅/调频信号转换与检测;针对大噪声干扰环境下的调频信号,通过多次希尔伯特变换进行相位解卷,重构了待测信号,实现了微位移的高精度测量。
激光自混合干涉调幅/调频转换原理如图1所示。其中,DFB为激光器,PD为光电探测器,PZT为微纳米运动平台,马赫-曾德尔干涉仪用于实现AM/FM信号转换[16]。
图1 调频信号产生原理图Fig.1 Frequency modulation signal generation schematic diagram
激光自混合干涉系统调幅信号的功率方程为[10]:
P=P0[1+βcos(ωτ)]
(1)
式中,P是有光反馈时的光功率;P0是无光反馈时候的初始光功率;β表示激光自混合干涉中的条纹可见度;ω为角频率;τ为光束在外腔往返一次所需要的时间。
自混合干涉系统的频率方程为:
ωτ=ω0τ-Csin(ωτ+arctanα)
(2)
式中,ω0为无光反馈时候的初始角频率;C为光反馈强度系数;α为线宽增强因子,系激光器的属性参数。用光频率表示的频率方程为:
2π(f-f0)τ=-Csin(2πfτ+arctanα)
(3)
令Δf=f-f0,有:
(4)
当自混合干涉信号经过马赫-曾德尔干涉仪后,通过线性拟合,干涉仪输出功率PFM为:
PFM=P0SΔf=Asin(2πfτ+arctanα)
(5)
式中,A为光束电矢量的振幅;S为干涉仪的强度传输系数T对光频率f的导数,即:
(6)
式中,n为光纤纤芯的折射率;c为光在真空中的速度;ΔL为振动位移。
由式(6)可知,调节马赫-曾德尔干涉仪的S可对光束滤波,其强度传输系数T与光频率f的对应变化曲线如图2所示。
图2 滤波器强度传输系数与光频率的对应变化曲线以及FM信号转换原理图Fig.2 The corresponding change curve of filter intensity transmission coefficient and light frequency and the schematic diagram of FM signal conversion
由图2可知,当T=0.5时,S取得最大值,即:
(7)
可见马赫-曾德尔干涉仪滤波后得到的调频信号为正弦信号,因为相位信息难以解调,故通过多重希尔伯特变换[17],得到与原信号PFM正交的信号PN:
PN=Acos(2πvτ+arctanα)
(8)
综合式(5)和(8),可计算出在进行位移重构时所需要的相位:
(9)
根据自混合干涉调幅信号与经过干涉仪滤波后的调频信号之间对应的相位关系,结合式(1)和(5),两信号之间的相位关系可以表示为:
(10)
分别对调幅和调频信号进行探测,然后对其进行相位提取,即可对线宽增强因子进行求解。
激光自混合干涉相干项的相位与外腔长度变化的关系为:
(11)
式(11)反映了相位与外腔长度变化的对应关系,实现了相位解卷。结合上面有关相位公式,即可进行位移重构,恢复出待测物体微位移轨迹。
基于上述原理,首先通过数值模拟验证方法的有效性,分析微位移测量精度。选择外部物体运动频率为100 Hz,振幅为3.1 μm,初始相位为0的简谐运动。激光器模拟波长为1550 nm,线宽增强因子α为4.6,光反馈因子C取为0.1,采样点数取为4000,可得到无噪声时调频信号仿真结果如图3所示。图3(a)中,“a”线为调频信号,“b”线调频信号经过希尔伯特变换后的信号。在图3(b)中,“a”线为模拟外部物体的简谐运动,“b”线为重构曲线,二者几乎重叠。图3(c)给出了重构位移与模拟位移的误差,可见该方法能够实现微位移的高精度测量。
图3 无噪声下调频信号仿真Fig.3 Simulation of frequency modulation signal without noise
对前述系统加入10 dB的高斯白噪声,数值仿真结果如图4所示。图4(a)中,“a”线为加噪调频信号,“b”线为希尔伯特变换后的信号。在图4(b)中,“a”线为模拟外部物体进行的简谐运动,“b”线为去噪后位移重构曲线。图4(c)给出了重构位移与模拟位移之间的误差。可见,在大噪声环境下,基于调频信号的重构仍具有很高的精度。
图4 加入10 dB噪声的数值仿真Fig.4 Numerical simulation with 10 dB noise
数值模拟和分析证明了本文提出的高精度微位移测量方法的有效性,为此开展实验测试,装置如图5所示。为模拟大噪声干扰环境,该实验选用的光源为无光隔离器无温控的光纤耦合的DFB激光器(THORLABS,S1FC1550PM),波长为1550 nm,使用的三个耦合器均为3 dB标准耦合度(分光比为50∶50),耦合器的四个传感臂的插入损耗分别为3.32 dB、3.37 dB、3.46 dB、3.60 dB。第一个耦合器一条输出臂连接准直镜头(THORLABS、F220FC-1550)照射在振动物体表面,振动物体选用镜面,放置在PZT(P753.1CD,PI)上,反射光再通过耦合器反射回激光腔形成自混合干涉现象,通过PD1(THORLABS,PDA20CS2)进行探测。耦合器另一条输出臂连接马赫-曾德尔干涉仪,干涉仪由两个3 dB耦合器和一段1 m的光纤(THORLABS,P1-1550A-FC-1)组成,再连接到PD2(THORLABS,PDA20CS2)进行调频信号的探测。PD1和PD2都通过数据采集模块(USB-4431,NI)获取电压信号,并在PC机上进行处理。在本次实验中,所有实验器材均放置在光学隔振平台(T1225QK,THORLABS)上,以避免外界振动的干扰。
图5 振动物体调幅/调频信号检测实验装置图Fig.5 An experimental device for detecting AM/FM signals of vibrating objects
设置目标物体振动频率为5 Hz,振动峰峰值为3.1 μm,采样频率设为50 kHz,两个光电探测器同一时刻探测到的调幅和调频信号如图6所示。可以看出,调频信号的振幅远大于调幅信号且比较光滑,有利于相位解卷。
图6 同时采集的调幅和调频信号Fig.6 Amplitude and frequency modulation signals are collected simultaneously
根据前述激光自混合干涉调频测量原理,基于采集到的调频信号进行微位移重构,结果如图7所示。图7(a)中,“a”线为采集到的调频信号,“b”线为调频信号经过希尔伯特变换后的信号。在图7(b)中,“a”线为外部物体振动曲线,“b”线为重构曲线。图7(c)给出了调频重构信号与原信号之间的误差,可见,基于调频信号可直接实现位移重构,误差在100 nm以内。
图7 基于调频信号的微位移测量实验结果Fig.7 Experimental results of micro displacement measurement based on frequency modulation signal
因为使用调幅信号难以直接进行位移重构,故首先对调幅信号做低通滤波,再进行位移重构,结果如图8所示。图8(a)中,“a”线为滤波后的调幅信号,“b”线为希尔伯特变换变换后的信号。在图8(b)中,“a”线为外部物体振动曲线,“b”线为重构曲线。图8(c)给出了调幅重构信号与原信号之间的误差,可见,滤波后的调幅信号进行位移重构的误差达到了520 nm。
图8 基于调幅信号的微位移测量实验结果Fig.8 Experimental results of micro-displacement measurement based on amplitude modulated signal
调整物体的位移,同时获取两调频和调幅信号,按照前述方法分析其信噪比,得到三组数据示于表1。可见,调频信号的信噪比远高于调幅信号的信噪比。
表1 调频和调幅信号的信噪比Tab.1 The signal-to-noise ratio of FM and AM signals
(1)基于激光自混合干涉系统,利用调频信号进行微位移测量,误差在100 nm以内,明显低于利用调幅信号测量的误差,适用于超精密测量。
(2)基于激光自混合干涉系统测量微位移,调频信号的信噪比远高于调幅信号的信噪比,适用于大噪声环境下高精度测量。