CN104990495B - 高分辨率频率扫描干涉仪中基于峰值演化消畸变的色散相位补偿方法 - Google Patents
高分辨率频率扫描干涉仪中基于峰值演化消畸变的色散相位补偿方法 Download PDFInfo
- Publication number
- CN104990495B CN104990495B CN201510443526.1A CN201510443526A CN104990495B CN 104990495 B CN104990495 B CN 104990495B CN 201510443526 A CN201510443526 A CN 201510443526A CN 104990495 B CN104990495 B CN 104990495B
- Authority
- CN
- China
- Prior art keywords
- compensation
- frequency
- interferometer
- phase
- 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.)
- Active
Links
Abstract
高分辨率频率扫描干涉仪中基于峰值演化消畸变的色散相位补偿方法,涉及扫描干涉仪色散补偿技术领域。本发明是为了解决辅助干涉仪光纤色散效应引起的校正后的测量干涉仪信号拍频随调频带宽和被测距离增大而产生线性变化导致的测量分辨率较低和测距误差较大的问题。本发明将光纤马赫泽德干涉仪频率采样法校正非线性后的信号乘以复相位补偿项后,得到根据补偿相位为将Ib表示为选择相位补偿系数αcomp,使‑πσdispn2+παcompn2最小,得到经过色散相位补偿的测量信号Ib,完成对频率扫描干涉仪色散影响的补偿。本发明适用于扫描干涉仪的色散补偿。
Description
技术领域
本发明涉及扫描干涉仪色散补偿技术领域。
背景技术
随着航空航天、船舶工业、重型机械、发电设备等行业的发展,大尺寸高精度测量问题突出涌现,如飞机大部件安装位置的测量、水汽轮机主轴长度的测量、大型精密机床床身的测量及水电站发电机组转子定子直径的测量等。三坐标机、经纬仪以及视觉测量等传统大型工件检测方法存在着接触被测工件、精度低、成本高及效率低等缺点,难以满足加工制造过程与现代科学研究中大尺寸高精度绝对距离测量的需求。频率扫描干涉仪具有非接触测量、测量精度高、低信噪比探测等优点,因此,在表面形貌测量、计量、及三维激光雷达领域中得到广泛应用。现代工业生产及计量要求测量具有高精度高分辨率的特点,高分辨率频率扫描干涉仪的分辨率能够达到亚毫米甚至几十微米。
频率扫描干涉仪采用频率采样法校正测量信号拍频非线性,但实际发现,提供外部采样时钟信号的辅助干涉仪光纤色散效应将引起校正非线性后的测量干涉仪信号拍频随调频带宽和被测距离增大产生线性变化,导致测量分辨率下降和测距误差增大。该现象主要由于辅助干涉仪由单模光纤构成,而测量干涉仪在空气中,两者色散系数不匹配,辅助干涉仪的光纤色散无法与测量干涉仪抵消,导致测量信号拍频频谱展宽。
发明内容
本发明的目的是为了解决辅助干涉仪光纤色散效应引起的校正后的测量干涉仪信号拍频随调频带宽和被测距离增大而产生线性变化导致的测量分辨率较低和测距误差较大的问题。
高分辨率频率扫描干涉仪中基于峰值演化消畸变的色散相位补偿方法,包括以下步骤:
经光纤马赫泽德干涉仪频率采样法校正非线性后的信号可表示为式(14)
其中,Ib为测量信号,AT为经目标反射回来的信号能量,AR是本振光信号能量;ηH为干涉效率,n为采样点序号;为采样后信号的畸变相位,主要由辅助干涉仪光纤色散引起;σdisp为相位畸变系数;τ和τaux分别表示高分辨率频率扫描干涉仪中测量干涉仪和辅助干涉仪对应的时间延迟;
将式(14)表示为复数形式,则有
其中,j为复数;
将采样后信号乘以复相位补偿项补偿采样后信号的相位畸变,如式(16):
其中,补偿相位为αcomp表示相位补偿系数;
选择相位补偿系数αcomp,使-πσdispn2+παcompn2最小,得到经过色散相位补偿的测量信号Ib′,完成对频率扫描干涉仪色散影响的补偿。
由于并不能直接获得σdisp的具体值,所以不能直接选取αcomp对采样后信号的相位畸变进行补偿;从式(16)可以看出当αcomp值逐渐趋近于σdisp时,采样后信号的啁啾分量将会减小,则其对应的频谱谱峰FWHM也将减小,当αcomp=σdisp时,则采样后信号的啁啾分量为0,表明相位畸变被完全消除;提高测量分辨率的关键在于选择合理的αcomp值,使其满足αcomp=σdisp的条件。
通过测量谱峰FWHM是否接近测量理论分辨率,可以判断信号相位畸变是否得到消除;当谱峰FWHM接近理论分辨率时,相位畸变逐渐得到补偿,反之,则相位畸变增大;但在实际中,由于目标表面特性影响,多个峰值发生严重频谱干涉时,将增加测量频谱FWHM的困难,且多个峰值无法区分时,只能测量整体峰值的FWHM,无法实现对每个峰值的最佳补偿。本发明通过分析峰值幅度变化与峰值FWHM之间的关系,发现在峰值FWHM减小的同时,峰值幅度呈现逐渐增大的趋势,从而为通过测量峰值幅度来判断峰值FWHM是否达到最小提供了依据。
所述选择相位补偿系数αcomp,使-πσdispn2+παcompn2最小,得到经过色散相位补偿的测量信号Ib′,完成对频率扫描干涉仪色散影响的补偿;包括以下步骤:
分析采样后信号函数的特点,将采样后信号的频谱幅度平方表示为式(17):
其中,ξ1表示频谱成分,Fresnelc(·)表示菲涅尔余弦积分函数、Fresnels(·)表示菲涅尔正弦积分函数,T表示采样后信号的时间长度;
令由式(17)可知,的基本形状只与有关;由于菲涅尔积分较为复杂,为方便分析,采用数值仿真的方法研究M与测量信号频谱之间的关系;假设T=1s,分别设置σdisp=6.5,σdisp=4.9,σdisp=3.3,σdisp=1.7,σdisp=0.1,则对应的M=3.6056,M=3.1305,M=2.5690,M=1.8439,M=0.4472。谱峰随M值的变化如图5所示。
图5中,σdisp越大,则M越大,谱峰展宽越严重,说明该信号含有的色散越大,此时,对应的谱峰峰值越小;而当M逐渐变小,谱峰展宽同时变小,此时,对应的谱峰峰值变大;因而可以通过谱峰变化判断色散是否消除。
令σd′isp=|αcomp-σdisp|,则经过每次色散相位补偿后的采样后信号的频谱幅度平方表示为:
调整αcomp,将每次色散补偿后的采样后信号以公式(18)的频谱形式表现出来,即根据公式(18)绘制随变化波形图,在绘制的波形图中选择峰值最高的波形图作为补偿效果最好的频谱,选取此时对应的αcomp对相位畸变进行消除,即选取此时对应的αcomp使得-πσdispn2+παcompn2最小;得到经过色散相位补偿的测量信号Ib′,完成对频率扫描干涉仪色散影响的补偿。
本发明提出基于峰值演化消畸变的色散相位补偿方法,将补偿后信号频谱峰值幅度变化作为判断标准并结合相位补偿方法,可实现对色散影响补偿,提高测量分辨率,降低误差。该方法具有对所有频谱峰值进行最佳补偿的优点,且算法复杂度低,判断方便。
对自由空间中2.4352m距离处的200μm高度差台阶进行测量,对采样后信号进行色散影响补偿前,峰值发生展宽,由于频谱干涉效应,谱峰形状产生畸变,无法区分真实目标峰值,导致测量错误。采用本发明提出的方法对采样后信号进行色散影响补偿后,可清晰分辨台阶两个面形成的距离峰值,且谱峰FWHM变窄,提高了测量分辨率,测量台阶高度为199.5μm,实现了大尺寸高精度高分辨率测量。
附图说明
图1为频率扫描干涉仪测量系统示意图
图2为光纤色散效应对辅助干涉仪本振光及发射光的影响图;
图3为采样后信号随调频带宽增大的距离峰变化图;
图4(a)-图4(f)分别为采样后信号在20nm扫频带宽下谱峰随距离增加的变化图;
图5为谱峰随不同M值的变化图;
图6(a)为由四个点目标形成的拍信号时域图;
图6(b)为由四个点目标形成的拍信号频谱图;
图7(a)为选择补偿系数48800对测量信号进行相位补偿后的频谱图;
图7(b)为选择补偿系数50000对测量信号进行相位补偿后的频谱图;
图7(c)为选择补偿系数50400对测量信号进行相位补偿后的频谱图;
图7(d)为选择补偿系数50800对测量信号进行相位补偿后的频谱图;
图8(a)为采样后信号的未进行色散影响补偿的距离峰效果图;
图8(b)为采样后信号的进行色散影响补偿后距离峰效果图;
图9(a)为对自由空间中2.4352m距离处的200μm高度差台阶测量结果图(未进行色散影响补偿的距离峰效果图);
图9(b)为对自由空间中2.4352m距离处的200μm高度差台阶测量结果图(进行色散影响补偿后距离峰效果图)。
具体实施方式
分析频率扫描干涉仪测量原理:
采用锯齿波对激光器进行线性调频,激光器出射角频率可表示为式(1):
ωL(t)=2π(f0+μt) (1)
其中,f0为激光调频起始频率,μ为调频斜率,t为激光器调频时间;
本振光信号可表示为式(2):
测量光到达被测目标并沿原路返回后对应的时延表示为τ,并与本振光形成拍频干涉,则测量信号可表示为式(3):
Sb∝exp(j2πμτt) (3)
啁啾线性调频激光测量方法的理论分辨率为式(4):
其中c为光速,Ω为激光器调频带宽;
由式(4)可以看出,增加调频带宽可以提高测量分辨率;
实际中,激光器调频特性无法满足理想线性调频,将导致测量信号的拍频频谱发生展宽,无法提取目标谱峰,使得测量结果产生很大误差;为避免该影响,采用频率采样法对测量信号拍频进行非线性校正。频率扫描干涉仪测量系统如图1所示,频率扫描干涉仪包括光学发射接收系统(a)、光学环形器(b)、第二3dB耦合器(c)、第二平衡探测器(d)、采集卡(e)、第一平衡探测器(f)、第一3dB耦合器(g)、第一普通耦合器(h)、第二普通耦合器(i)、第三普通耦合器(j)、光学隔离器(k)和激光器(l)。激光器(l)为外腔调频激光器。
采用锯齿波对激光器进行线性调频,激光器(l)发出的光束经第一普通耦合器(h)后分为两路,一路为马赫泽德干涉光路(辅助干涉仪,图1中下方的虚线框所示),经两臂光程差220m后,在第一平衡探测器(f)上形成触发信号,将触发信号的相位每经过2π的时刻作为外部采样时钟接入采集卡;另一路为为测量光路(测量干涉仪,图1中上方的虚线框所示),光路经第二普通耦合器(i)后分为两路,第一路经光学环形器(b)、光学发射接收系统(a)到达目标,并返回与经第二普通耦合器(i)后的第二路的本振光在第二平衡探测器(d)上形成测量信号,并接入到采集卡;
触发信号对测量信号同步采样后,测量信号的拍频非线性得到校正,测量信号被触发信号同步采样后的信号称作采样后信号,通过对采样后信号进行Fourier变换,可以计算出目标距离。其中光纤端面为测量零点,光学系统中除光学发射接收系统外,均由单模光纤构成。
频率采样法的原理是通过构造辅助干涉仪产生触发信号,并将触发信号相位每经过2π的时刻作为外部时钟对测量信号进行同步采样,则采样后信号的幅度过零点之间的间隔变为相等,经Fourier变换后,在频谱图上形成单一峰值,完成测量信号频率非线性校正。采用频率采样法可以校正测量信号的非线性主要是由于辅助干涉仪和测量干涉仪的光来源于同一个激光器,它们的调频非线性函数形式是相同的,因此,它们的比值仅和两光路的光程差有关。
经频率采样法校正(消除非线性)后的测量信号为采样后信号,可表示为:
其中τ和τaux分别表示测量干涉仪和辅助干涉仪对应的时间延迟,根据奈奎斯特采样定理,需要满足τaux≥2τ。由式(5)可以看出,通过预先标定出辅助干涉干涉仪的光程,可以得到τaux,进而通过Fourier变换得到被测目标时延τ,实现测距。
分析式(5),若辅助干涉仪的光路部分在空气中,则τaux为固定值,τ和τaux之间成正比例关系,通过Fourier变换计算该比值,并乘以τaux,可以得到τ,进而完成测量,此时采样后信号的频率为单一频率,其频率不随时间变化。
但在实际中,为实现大尺寸测量,需要对辅助干涉仪设置较长的光程差,通常采用单模光纤构成辅助干涉仪,为实现高分辨率测量,要求激光器调频带宽很大,在采用宽带外腔调频激光器和220m光程差的光纤辅助干涉仪后,实验发现,由单模光纤色散效应引起辅助干涉仪的时延τaux随着调频带宽的增加而变化,此时,τ和τaux之间将不成比例关系,采样后信号的频率随时间发生变化,表现为频谱展宽,测量分辨率下降。对该过程进行理论分析,如图2所示,其中ω0为调频起始角频率,ΔΩ为调频带宽,(ω0+Δω)表示调频任意时刻角频率,Tm为调频周期。通过辅助干涉仪光纤短臂的本振光频率表示为ωL(t),通过辅助干涉仪光纤长臂的传输光频率表示为ωT(t),经过两臂时延τaux后,在第一平衡探测器(f)上形成触发信号。若不考虑光纤色散,如辅助干涉仪光路在空气中,则本振光和传输光形成的拍频不随调频带宽发生变化,为恒定频率。当辅助干涉仪由光纤构成时,光纤色散导致在第一平衡探测器(f)上形成的触发信号频率随调频带宽发生变化,等效于校正测量信号非线性的外部时钟时间间隔随调频带宽发生变化,进而导致测量信号拍频随调频带宽发生变化,频谱展宽,测量分辨率下降。
首先分析在宽带调频情况下,由单模光纤构成的辅助干涉仪形成的拍频频率。
由激光器(l)发出的线性调频光经辅助干涉仪的短臂后,其频率及光信号可分别表示为式(1)和式(2)。激光器(l)发出的光经辅助干涉仪的长臂后,单模光纤的色散效应将引起辅助干涉仪长臂的时延发生变化,可表示为τGVD;
τGVD=β2ΔωLaux (6)
式中β2=-20ps2/km表示单模光纤的群速度色散,Δω=2πμt为某一时刻的调频带宽,t为调频时间,Laux为辅助干涉仪光纤长度;经过辅助干涉仪长臂后的光频可表示为式(7):
选用G.652.D标准的单模光纤,1550nm的有效群折射率为ng=1.4682,vg=2.0437×108m/s为群速度,τaux为不考虑光纤色散情况下的辅助干涉仪长臂对应的时延,t是调频时间;则经过光纤长臂后的光相位可表示为:
激光器(l)发出的光经过光纤短臂后的光频可表示为式(1),其与长臂发生干涉形成触发信号的相位可表示为式(9):
对式(9)相位求导,辅助干涉仪的拍频频率表示为式(10):
由式(10)可以看出光纤色散效应导致辅助干涉仪的拍频发生啁啾变化,即2πμβ2vgt项,且随着调频时间t增加,即调频带宽相应增大,啁啾量变大。
采用频率采样法对测量信号进行采样,参考公式(5),测量干涉仪的时延与辅助干涉仪的时延比值变为:
调谐波长选择1542nm-1562nm,换算为调频带宽为μt=2.4911THz,由于2πμβ2vgt<<1,忽略二次以上的项后,式(11)可表示为如下:
对应于220m光程差的时延为τaux=733ns,则|2πμβ2vgτaux|<10-9,其影响可以忽略,对式(12)整理后,如式(13)所示。
当调频时间t增加,调频带宽相应增大。由式(13)可以看出在考虑辅助干涉仪光纤色散情况下,测量信号经频率采样法校正非线性后,测量干涉仪时延τ和辅助干涉仪时延τaux(t)之间将不具备比例关系,而是在采样后信号中引入了啁啾量2πτμβ2vgt,这表明采样后信号的频率不是单一频率,而是随调频带宽增大而线性变化,从而导致对采样后信号进行Fourier变换后,频谱峰值发生展宽,测量分辨率降低。对采样后信号距离谱图的仿真如图3所示,图中真实距离为4.55231257m,当选择扫频波长分别为1542-1547nm,1542-1552nm,1542-1557nm,1542-1562nm时,测距峰值展宽越来越严重,测量其与真实值之间的偏差分别为36.48μm,73.00μm,109.24μm,147.34μm。表明在测量大距离目标时,随着扫频带宽增大,测量分辨率下降的同时,测距值与真实值的误差在增大。由式(13)可以看出,当被测目标距离增大时,采样后信号拍频含有的啁啾分量变大,导致频谱展宽增大。
采样后信号的距离谱图随被测距离增加变化的仿真如图4(a)-图4(f)所示,由图中可以看出在扫频范围1542-1562nm下,随着被测距离增大,测距峰值展宽越来越严重,且与真实距离的误差也变大。采用宽带扫频光源进行大尺寸测量时,光纤色散效应导致辅助干涉仪的光程发生变化,引起拍频随扫频带宽变化,等效于外部采样时钟的时间间隔产生线性变化,测量信号经外部时钟采样后,采样后信号拍频也将发生线性变化,从而引起距离峰值发生展宽,降低测量分辨率,同时测量误差增大。
为实现大尺寸高分辨率测量,需要对频率扫描干涉仪测量系统中辅助干涉仪光纤色散影响进行补偿。
具体实施方式一:高分辨率频率扫描干涉仪中基于峰值演化消畸变的色散相位补偿方法,包括以下步骤:
经光纤马赫泽德干涉仪频率采样法校正非线性后的信号可表示为式(14)
其中,Ib为测量信号,AT为经目标反射回来的信号能量,AR是本振光信号能量;ηH为干涉效率,n为采样点序号;为采样后信号的畸变相位,主要由辅助干涉仪光纤色散引起;σdisp为相位畸变系数;τ和τaux分别表示高分辨率频率扫描干涉仪中测量干涉仪和辅助干涉仪对应的时间延迟;
将式(14)表示为复数形式,则有
其中,j为复数;
将采样后信号乘以复相位补偿项补偿采样后信号的相位畸变,如式(16):
其中,补偿相位为αcomp表示相位补偿系数;
选择相位补偿系数αcomp,使-πσdispn2+παcompn2最小,得到经过色散相位补偿的测量信号Ib′,完成对频率扫描干涉仪色散影响的补偿。
由于并不能直接获得σdisp的具体值,所以不能直接选取αcomp对采样后信号的相位畸变进行补偿;从式(16)可以看出当αcomp值逐渐趋近于σdisp时,采样后信号的啁啾分量将会减小,则其对应的频谱谱峰FWHM也将减小,当αcomp=σdisp时,则采样后信号的啁啾分量为0,表明相位畸变被完全消除;提高测量分辨率的关键在于选择合理的αcomp值,使其满足αcomp=σdisp的条件。
通过测量谱峰FWHM是否接近测量理论分辨率,可以判断信号相位畸变是否得到消除;当谱峰FWHM接近理论分辨率时,相位畸变逐渐得到补偿,反之,则相位畸变增大;但在实际中,由于目标表面特性影响,多个峰值发生严重频谱干涉时,将增加测量频谱FWHM的困难,且多个峰值无法区分时,只能测量整体峰值的FWHM,无法实现对每个峰值的最佳补偿。本发明通过分析峰值幅度变化与峰值FWHM之间的关系,发现在峰值FWHM减小的同时,峰值幅度呈现逐渐增大的趋势,从而为通过测量峰值幅度来判断峰值FWHM是否达到最小提供了依据。
具体实施方式二:所述选择相位补偿系数αcomp,使-πσdispn2+παcompn2最小,得到经过色散相位补偿的测量信号Ib′,完成对频率扫描干涉仪色散影响的补偿;包括以下步骤:
分析采样后信号函数的特点,将采样后信号的频谱幅度平方表示为式(17):
其中,ξ1表示频谱成分,Fresnelc(·)表示菲涅尔余弦积分函数、Fresnels(·)表示菲涅尔正弦积分函数,T表示采样后信号的时间长度;
令由式(17)可知,的基本形状只与有关;由于菲涅尔积分较为复杂,为方便分析,采用数值仿真的方法研究M与测量信号频谱之间的关系;假设T=1s,分别设置σdisp=6.5,σdisp=4.9,σdisp=3.3,σdisp=1.7,σdisp=0.1,则对应的M=3.6056,M=3.1305,M=2.5690,M=1.8439,M=0.4472。谱峰随M值的变化如图5所示。
图5中,σdisp越大,则M越大,谱峰展宽越严重,说明该信号含有的色散越大,此时,对应的谱峰峰值越小;而当M逐渐变小,谱峰展宽同时变小,此时,对应的谱峰峰值变大;因而可以通过谱峰变化判断色散是否消除。
令σd′isp=|αcomp-σdisp|,则经过每次色散相位补偿后的采样后信号的频谱幅度平方表示为:
调整αcomp,将每次色散补偿后的采样后信号以公式(18)的频谱形式表现出来,即根据公式(18)绘制随变化波形图,在绘制的波形图中选择峰值最高的波形图作为补偿效果最好的频谱,选取此时对应的αcomp对相位畸变进行消除,即选取此时对应的αcomp使得-πσdispn2+παcompn2最小;得到经过色散相位补偿的测量信号Ib′,完成对频率扫描干涉仪色散影响的补偿。
本发明提出基于峰值演化消畸变的色散相位补偿方法,将补偿后信号频谱峰值幅度变化作为判断标准并结合相位补偿方法,可实现对色散影响补偿,提高测量分辨率。该方法具有对所有频谱峰值进行最佳补偿的优点,且算法复杂度低,判断方便。
实施例
仿真过程如下:假设由四个点目标形成的测量信号拍频依次为5500Hz、5800Hz、5900Hz、6000Hz,为便于说明原理,增加对比效果,其相应的最大畸变系数分别设为4.88×104,5×104,5.04×104,5.08×104,该测量信号的时域图和频域图分别如图6(a)和图6(b)所示,由图6(b)可以看出信号发生严重频谱干涉,且产生伪峰,4个点目标无法区分和测量。为消除该影响,将采样后的信号与相乘,并分别设置补偿系数αcomp=48800,αcomp=50000,αcomp=50400,αcomp=50800,然后对补偿后的信号进行Fourier变换,则其频谱图如图7所示。
对比图7(a)-图7(d)和图6(b),可以看出经补偿后的采样后信号的频谱图四个点目标形成的谱峰可清晰分辨。进一步由图7(a)可以看出,当选择αcomp=48800时,频率为5500Hz的拍信号相位畸变得到完全补偿,其对应的频谱峰值1得到最佳补偿,其高度相对于自身的变化达到最大值,而谱峰2、3、4的补偿效果依次变差。改变αcomp的值,使αcomp=50000时,频率为5800Hz的拍信号相位畸变得到完全补偿,其对应的频谱峰值2得到最佳补偿,其高度相对于自身的变化达到最大值,而谱峰1、3、4的补偿效果变差,如图7(b)所示。以此类推,进一步改变αcomp的值,使其分别为αcomp=50400,αcomp=50800,则频率分别为5.04×104,5.08×104的拍信号相位畸变得到完全补偿,分别使频率3和频率4得到最佳补偿,如图7(c)和图7(d)所示。由图7(a)-7(d)可知,通过连续改变补偿系数αcomp的值,可以实现对所有目标形成的拍信号相位畸变的最佳补偿,提高测量分辨率。
频率扫描干涉仪采用下调频的方式,输出功率1.5mW,波长调谐范围选择1542-1562nm,调频速度为100nm/s,采用频率采样法校正测量信号拍频非线性,其中辅助干涉仪的光程差为220m,设置采集卡为外部触发和外部采样时钟的工作模式。
辅助干涉仪光纤色散对频率扫描干涉仪测量系统影响补偿实验,将采样后信号按照式(16)进行相位补偿,并根据频谱峰值幅度变化调节补偿系数,当峰值幅度相对于自身变化达到最大时,即认为完成色散影响补偿。以光纤端面4.5242m为测量零点,对自由空间中2.53m距离处的量块进行测量,色散影响补偿前后的距离峰值对比如图8所示,图8(a)为补偿前峰值高度为965.3,峰值FWHM为166.8μm,图8(b)为补偿后峰值高度为1396,峰值FWHM为64.5μm,接近测量理论分辨率60.2μm。
对自由空间中2.4352m距离处的200μm高度差台阶进行测量,其峰值进行相位补偿前后结果如图9所示,图9(a)显示对采样后信号进行色散影响补偿前,峰值发生展宽,由于频谱干涉效应,谱峰形状产生畸变,无法区分真实目标峰值,导致测量错误。图9(b)显示采用我们提出的方法对采样后信号进行色散影响补偿后,可清晰分辨台阶两个面形成的距离峰值,且谱峰FWHM变窄,提高了测量分辨率,测量台阶高度为199.5μm,实现了大尺寸高精度高分辨率测量。
Claims (2)
1.高分辨率频率扫描干涉仪中基于峰值演化消畸变的色散相位补偿方法,其特征在于,包括以下步骤:
经光纤马赫泽德干涉仪频率采样法校正非线性后的信号可表示为式(14)
其中,Ib为测量信号,AT为经目标反射回来的信号能量,AR是本振光信号能量;ηH为干涉效率,n为采样点序号;为采样后信号的畸变相位;σdisp为相位畸变系数;τ和τaux分别表示高分辨率频率扫描干涉仪中测量干涉仪和辅助干涉仪对应的时间延迟;
将式(14)表示为复数形式,则有
其中,j为复数;
将采样后信号乘以复相位补偿项补偿采样后信号的相位畸变,如式(16):
其中,补偿相位为αcomp表示相位补偿系数;
选择相位补偿系数αcomp,使-πσdispn2+παcompn2最小,得到经过色散相位补偿的测量信号Ib′,完成对频率扫描干涉仪色散影响的补偿:
将采样后信号的频谱幅度平方表示为
其中,σ′disp=|αcomp-σdisp|;ξ1表示频谱成分,Fresnelc(·)表示菲涅尔余弦积分函数、Fresnels(·)表示菲涅尔正弦积分函数,T表示采样后信号的时间长度;
调整αcomp,将每次色散补偿后的采样后信号以公式(18)的频谱形式表现出来,即根据公式(18)绘制随变化波形图,在绘制的波形图中选择峰值最高的波形图作为补偿效果最好的频谱,选取此时对应的αcomp对相位畸变进行消除,即选取此时对应的αcomp使得-πσdispn2+παcompn2最小;得到经过色散相位补偿的测量信号Ib′,完成对频率扫描干涉仪色散影响的补偿。
2.根据权利要求1所述的高分辨率频率扫描干涉仪中基于峰值演化消畸变的色散相位补偿方法,其特征在于:所述选择相位补偿系数αcomp,使-πσdispn2+παcompn2最小,得到经过色散相位补偿的测量信号Ib′,完成对频率扫描干涉仪色散影响的补偿中的得到过程包括以下步骤:
将采样后信号的频谱幅度平方表示为式(17):
其中,ξ1表示频谱成分,Fresnelc(·)表示菲涅尔余弦积分函数、Fresnels(·)表示菲涅尔正弦积分函数,T表示采样后信号的时间长度;
令σ′disp=|αcomp-σdisp|,则经过每次色散相位补偿后的采样后信号的频谱幅度平方表示为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510443526.1A CN104990495B (zh) | 2015-07-24 | 2015-07-24 | 高分辨率频率扫描干涉仪中基于峰值演化消畸变的色散相位补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510443526.1A CN104990495B (zh) | 2015-07-24 | 2015-07-24 | 高分辨率频率扫描干涉仪中基于峰值演化消畸变的色散相位补偿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104990495A CN104990495A (zh) | 2015-10-21 |
CN104990495B true CN104990495B (zh) | 2017-07-28 |
Family
ID=54302338
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510443526.1A Active CN104990495B (zh) | 2015-07-24 | 2015-07-24 | 高分辨率频率扫描干涉仪中基于峰值演化消畸变的色散相位补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104990495B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105259548B (zh) * | 2015-10-29 | 2017-10-27 | 哈尔滨工业大学 | 一种用于fmcw绝对距离测量技术中色散失配校正方法 |
CN105954735B (zh) * | 2016-07-15 | 2018-06-22 | 哈尔滨工业大学 | 一种用于fmcw绝对距离测量技术中改进的高速色散失配校正方法 |
CN108803006B (zh) * | 2017-09-18 | 2021-01-05 | 成都理想境界科技有限公司 | 光纤扫描成像系统、设备及其畸变检测与矫正系统 |
CN110068828B (zh) * | 2018-08-07 | 2023-10-17 | 天津大学 | 基于激光调频连续波远距离测距的装置及色散补偿方法 |
CN110730055B (zh) * | 2019-10-22 | 2022-07-05 | 上海创远仪器技术股份有限公司 | 基于信号分析仪实现5g信号发射调制质量测量的方法 |
CN111239657B (zh) * | 2020-01-20 | 2022-05-06 | 上海东软医疗科技有限公司 | 谱图的相位校正方法、装置及设备 |
CN112051583B (zh) * | 2020-08-25 | 2022-06-14 | 哈尔滨工业大学 | Fmcw距离测量系统中拍频信号非线性校正方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6160627A (en) * | 1998-05-21 | 2000-12-12 | Electronics And Telecommunications Research Institute | Optical fiber Mach-Zehnder interferometer filter |
CN101181153A (zh) * | 2007-12-12 | 2008-05-21 | 中国科学院上海光学精密机械研究所 | 高分辨率光学相干层析成像方法 |
CN101839686A (zh) * | 2010-03-26 | 2010-09-22 | 中国计量科学研究院 | 激光干涉仪非线性误差修正方法、装置及应用其的干涉仪 |
CN202617124U (zh) * | 2012-05-31 | 2012-12-19 | 陕西理工学院 | 具有光信号色散补偿和相位补偿功能的干涉装置 |
CN104771144A (zh) * | 2015-04-13 | 2015-07-15 | 中国科学院上海光学精密机械研究所 | 频域光学相干层析成像系统的色散补偿方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE10130159C2 (de) * | 2001-06-22 | 2003-06-18 | Litef Gmbh | Verfahren zur Vermeidung von Bias-Fehlern aufgrund synchroner Einstreuung bei faseroptischen Gyroskopen |
-
2015
- 2015-07-24 CN CN201510443526.1A patent/CN104990495B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6160627A (en) * | 1998-05-21 | 2000-12-12 | Electronics And Telecommunications Research Institute | Optical fiber Mach-Zehnder interferometer filter |
CN101181153A (zh) * | 2007-12-12 | 2008-05-21 | 中国科学院上海光学精密机械研究所 | 高分辨率光学相干层析成像方法 |
CN101839686A (zh) * | 2010-03-26 | 2010-09-22 | 中国计量科学研究院 | 激光干涉仪非线性误差修正方法、装置及应用其的干涉仪 |
CN202617124U (zh) * | 2012-05-31 | 2012-12-19 | 陕西理工学院 | 具有光信号色散补偿和相位补偿功能的干涉装置 |
CN104771144A (zh) * | 2015-04-13 | 2015-07-15 | 中国科学院上海光学精密机械研究所 | 频域光学相干层析成像系统的色散补偿方法 |
Non-Patent Citations (1)
Title |
---|
Ultrahigh-resolution,high-speed,Fourier domain optical coherence tomography and methods for dispersion compensation;M.Wojtkowski et al;《Optics Express》;20041231;第12卷(第11期);第2404-2422页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104990495A (zh) | 2015-10-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104990495B (zh) | 高分辨率频率扫描干涉仪中基于峰值演化消畸变的色散相位补偿方法 | |
CN105259548B (zh) | 一种用于fmcw绝对距离测量技术中色散失配校正方法 | |
CN110068828B (zh) | 基于激光调频连续波远距离测距的装置及色散补偿方法 | |
CN106895790B (zh) | 一种光频域反射中提升分布式光纤传感分辨率方法 | |
Zhang et al. | Simultaneous measurements of velocity and distance via a dual-path FMCW lidar system | |
CN105136021B (zh) | 基于调焦清晰度评价函数的激光频率扫描干涉仪色散相位补偿方法 | |
CN102332956B (zh) | 一种宽带光源的色散补偿方法 | |
JP2019045200A (ja) | 光学的距離測定装置および測定方法 | |
CN103176173B (zh) | 基于光纤采样技术的lfmcw激光雷达调频的非线性校正方法 | |
CN105091740B (zh) | 一种基于标定高分辨率频率扫描干涉仪色散啁啾斜率补偿色散的方法 | |
CN112923960B (zh) | 用于校正非线性调谐效应的光纤参数测量装置 | |
CN105137444B (zh) | 双干涉光路调频连续波激光测距信号处理方法 | |
CN113029034A (zh) | 激光扫频干涉测量的非线性校正与量程扩展装置及方法 | |
CN102819022A (zh) | 一种基于实时色散傅立叶变换的飞秒测距激光雷达和方法 | |
CN107144537B (zh) | 一种可见光傅里叶变换吸收光谱测量方法及系统 | |
CN108663684A (zh) | 一种基于等光频间隔重采样的相位差测距方法 | |
CN103163513A (zh) | 基于相位解调方法的fmcw激光雷达高精度信号测量方法 | |
CN105954735A (zh) | 一种用于fmcw绝对距离测量技术中改进的高速色散失配校正方法 | |
CN104914444A (zh) | 一种远距离激光外差干涉测距结构 | |
Zehao et al. | FMCW LiDAR with an FM nonlinear kernel function for dynamic-distance measurement | |
CN109541621B (zh) | 一种频率扫描干涉绝对测距系统的振动补偿方法 | |
CN116930995B (zh) | 调频连续波激光高速目标的速度和距离测量系统及方法 | |
CN101738167A (zh) | 基于谐振腔稳频的绝对距离测量系统及实现方法 | |
Niu et al. | Arbitrary distance measurement without dead zone by chirped pulse spectrally interferometry using a femtosecond optical frequency comb | |
CN104792269B (zh) | 一种对线性相移误差不敏感的光纤端面高度值的解算方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |