CN111595468B - 一种补偿载波相位延迟非线性误差的pgc相位解调方法 - Google Patents

一种补偿载波相位延迟非线性误差的pgc相位解调方法 Download PDF

Info

Publication number
CN111595468B
CN111595468B CN202010396826.XA CN202010396826A CN111595468B CN 111595468 B CN111595468 B CN 111595468B CN 202010396826 A CN202010396826 A CN 202010396826A CN 111595468 B CN111595468 B CN 111595468B
Authority
CN
China
Prior art keywords
phase
input end
output end
multiplier
order
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
Application number
CN202010396826.XA
Other languages
English (en)
Other versions
CN111595468A (zh
Inventor
严利平
陈本永
张倚得
谢建东
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University of Technology ZJUT
Original Assignee
Zhejiang University of Technology ZJUT
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang University of Technology ZJUT filed Critical Zhejiang University of Technology ZJUT
Priority to CN202010396826.XA priority Critical patent/CN111595468B/zh
Publication of CN111595468A publication Critical patent/CN111595468A/zh
Application granted granted Critical
Publication of CN111595468B publication Critical patent/CN111595468B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J9/00—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength
    • G01J9/02—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength by interferometric methods
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D5/00—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable
    • G01D5/26—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light
    • G01D5/32—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light
    • G01D5/34—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells
    • G01D5/353—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre
    • G01D5/35306—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre using an interferometer arrangement
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J9/00—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength
    • G01J9/02—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength by interferometric methods
    • G01J2009/0249—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength by interferometric methods with modulation

Landscapes

  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

本发明公开了一种补偿载波相位延迟非线性误差的PGC相位解调方法。采用高频正弦和低频正弦扫描复合的调制方法，对采样得到的复合正弦调制数字干涉信号进行一阶、二阶正交下混频，获得两对正交幅值信号，运用一阶正交幅值信号及一阶正交微分信号准确提取出相位延迟并对相位延迟进行补偿，根据相位延迟大小对一阶、二阶正交幅值信号进行优选获得两路不受相位延迟影响的新幅值信号，再通过反正切与相位解包裹运算求得复合相位，最后通过滑动平均计算获得待测相位。本发明解决了PGC相位解调技术中相位延迟难以准确测量及相位延迟引起的非线性误差难以实时补偿的问题，提高了相位测量精度，能广泛应用于正弦相位调制干涉技术领域。

Description

一种补偿载波相位延迟非线性误差的PGC相位解调方法
技术领域
本发明属于激光干涉测量技术领域，特别是一种补偿载波相位延迟非线性误差的PGC相位解调方法。
背景技术
相位生成载波(PGC)解调技术因其抗低频干扰、灵敏度高、动态范围大等优点，被广泛应用于干涉型光纤传感器和正弦相位调制干涉仪中。PGC解调技术主要包括微分交叉相乘算法(PGC-DCM)和反正切算法(PGC-Arctan)。PGC-DCM法通过对正交信号进行微分交叉相乘以及积分等运算获得待测相位，容易受激光光强波动、载波相位延迟、相位调制深度的影响。PGC-Arctan法将正交分量进行除法以及反正切运算直接得到待测相位，消除了光强波动带来的影响，但仍会受载波相位延迟和调制深度的影响。其中载波相位延迟的补偿是进行调制深度补偿的前提，为了补偿载波相位延迟带来的影响，通常先使用正交解调方法求得载波相位延迟值，并在参考载波信号中增加相位补偿器使得干涉信号的载波项与参考载波信号保持同相位。正交解调法在待测相位为特定值时，无法求得载波相位延迟。此外，在实际中，相位延迟会随环境变化出现漂移，现有方法难以对相位延迟实现准确的快速补偿，导致非线性误差的产生，这限制了测相精度的提高。所以，准确地提取PGC相位解调算法中的相位延迟并补偿相位延迟的漂移是提高正弦调制干涉测量精度需要解决的关键技术问题。
发明内容
为了克服现有技术中的不足，本发明公开了一种补偿载波相位延迟非线性误差的PGC相位解调方法，解决了PGC解调中相位延迟带来的影响，解决了PGC相位解调技术中相位延迟难以准确测量及相位延迟引起的非线性误差难以实时补偿的问题，提高了相位测量精度，能广泛应用于正弦相位调制干涉技术领域。
本发明采用的技术方案包括以下步骤：
对正弦相位调制干涉仪中的电光相位调制器施加包含高频正弦调制信号和低频正弦扫描信号的复合调制信号，实现对干涉信号相位的复合调制；
通过带通滤波器滤除正弦相位调制干涉仪输出的复合相位调制干涉信号中的直流成分和高频噪声，然后对复合相位调制干涉信号进行模数采样得到正弦相位调制数字干涉信号S(t)，表达式如下：
Figure BDA0002487888390000021
其中，A为正弦相位调制数字干涉信号的幅值，m为相位调制深度，J0(m)为零阶第一类贝塞尔函数，J2n(m)和J2n-1(m)分别为偶数阶和奇数阶第一类贝塞尔函数，n表示阶数，ωc为高频正弦调制信号的角频率，θ为载波相位延迟，t表示时间，
Figure BDA0002487888390000022
表示正弦相位调制数字干涉信号的复合相位；
相位
Figure BDA0002487888390000023
包含t时刻正弦相位调制干涉仪中的待测相位和扫描相位，公式如下：
Figure BDA0002487888390000024
其中，
Figure BDA0002487888390000025
为低频正弦扫描相位，
Figure BDA0002487888390000026
为待测相位，B为扫描相位的幅值，ωs为扫描相位的角频率；
所述的待测相位为正弦相位调制干涉仪中待测对象进行位移引起的偏移相位。扫描相位为正弦相位调制干涉仪中利用低频正弦扫描信号进行扫描时施加的相位。
数字频率合成器(1、18)产生的一阶参考信号(cosωct，sinωct)及二阶参考信号(cos2ωct，sin2ωct)分别与正弦相位调制数字干涉信号S(t)相乘，并进行低通滤波，得到一阶正交幅值信号(P1，Q1)和二阶正交幅值信号(P2，Q2)，公式如下：
Figure BDA0002487888390000027
Figure BDA0002487888390000028
Figure BDA0002487888390000029
Figure BDA00024878883900000210
其中，LPF[]表示低通滤波操作；sin(ωct)、cos(ωct)分别表示一阶参考信号的正弦分量和余弦分量，sin(2ωct)、cos(2ωct)分别表示二阶参考信号的正弦分量和余弦分量，P1,Q1分别表示一阶正交幅值信号的余弦幅值分量和正弦幅值分量，P2,Q2分别表示二阶正交幅值信号的余弦幅值分量和正弦幅值分量；
然后一阶正交幅值信号(P1，Q1)经过微分运算后得到一阶正交微分信号(DP，DQ)：
Figure BDA00024878883900000211
Figure BDA0002487888390000031
其中，
Figure BDA0002487888390000032
为复合相位
Figure BDA0002487888390000033
对时间t的偏微分，DP和DQ分别表示一阶正交微分信号的余弦微分分量和正弦微分分量；
利用一阶正交幅值信号(P1，Q1)、一阶正交微分信号(DP，DQ)进行运算，得到载波相位延迟θc，计算公式如下：
Figure BDA0002487888390000034
式中，sign()表示符号函数，当括号中的数值大于等于零时，其值为1，当括号中的数值小于零时，其值为-1；
计算得到的载波相位延迟θc的取值范围为-π/2到π/2。
运用步骤中计算得到的载波相位延迟θc、一阶正交幅值信号(P1，Q1)和二阶正交幅值信号(P2，Q2)，重构出一对幅值不受载波相位延迟影响的新幅值信号(R1，R2)计算公式如下：
Figure BDA0002487888390000035
Figure BDA0002487888390000036
其中，R1和R2分别表示新幅值信号的正弦幅值分量和余弦幅值分量；
对新幅值信号(R1，R2)进行四象限反正切运算，得到包裹相位
Figure BDA0002487888390000037
计算公式如下：
Figure BDA0002487888390000041
式中计算得到的包裹相位
Figure BDA0002487888390000042
包裹在-π到+π之间。
具体实施中，假定调制深度m＝2.63，J1(m)＝J2(m)。
对包裹相位
Figure BDA0002487888390000043
进行相位解包裹得到连续变化的复合相位
Figure BDA0002487888390000044
根据一个低频正弦扫描信号的周期内的扫描相位
Figure BDA0002487888390000045
均值为零的性质，在一个低频正弦扫描信号的周期内共有M个数据，采用长度为M的队列存储复合相位
Figure BDA0002487888390000046
对存储的M个数据进行求和运算，并将求和运算的结果除以M完成滑动平均运算，消除复合相位中的扫描相位，最终得到待测相位
Figure BDA0002487888390000047
公式如下：
Figure BDA0002487888390000048
其中，U[]表示相位解包裹运算，∑[]表示M个数据的求和运算。
方法采用以下PGC相位解调系统，第一乘法器、第二乘法器、第三乘法器和第四乘法器的输入端均连接数字干涉信号S(t)；第一数字频率合成器的输出端分别连接至第一乘法器和第二乘法器的输入端，第一乘法器的输出端连接至第一低通滤波器的输入端，第二乘法器的输出端连接至第二低通滤波器的输入端；第一低通滤波器的输出端分别连接至第二平方运算器的输入端、第一符号提取器的输入端、第一微分运算器的输入端和相位延迟补偿模块的输入端，第一微分运算器的输出端连接至第一平方运算器的输入端；第二低通滤波器的输出端连接至第三平方运算器的输入端、第二符号提取器的输入端、第二微分运算器的输入端和相位延迟补偿模块的输入端，第二微分运算器的输出端连接至第四平方运算器的输入端；第一平方运算器的输出端和第二平方运算器的输出端均连接至第一加法器的输入端，第三平方运算器的输出端和第四平方运算器的输出端均连接至第二加法器的输入端，第一加法器的输出端经第一开方运算器后和第一符号提取器的输出端一起连接至第五乘法器的输入端，第二加法器的输出端经第二开方运算器后和第二符号提取器的输出端一起连接至第五乘法器的输入端，第五乘法器的输出端和第六乘法器的输出端均连接至第一反正切运算器的输入端；第二数字频率合成器的输出端连接至第三乘法器和第四乘法器的输入端，第三乘法器和第四乘法器的输出端分别经第三低通滤波器和第四低通滤波器连接至相位延迟补偿模块的输入端，第一反正切运算器的输出端连接至相位延迟补偿模块的输入端，相位延迟补偿模块的两个输出端均连接至第二反正切运算器的输入端，第二反正切运算器输出端连接至相位解包裹处理器的输入端，相位解包裹处理器的输出端连接至滑动平均处理器的输入端，滑动平均处理器的输出端输出待测相位。
所述的相位延迟补偿模块具体为：第一反正切运算器的输出端分别连接至一阶新幅值信号选择器的输入端、二阶新幅值信号选择器的输入端、第一正/余弦运算器的输入端，第一正/余弦运算器的两个输出端分别连接至第一除法器的输入端和第二除法器的输入端，同时第一反正切运算器的输出端经倍乘器连接至第二正/余弦运算器的输入端，第二正/余弦运算器的两个输出端分别连接至第三除法器的输入端和第四除法器的输入端，第一低通滤波器的输出端和第二低通滤波器的输出端分别连接至第一除法器的输入端和第二除法器的输入端，第三低通滤波器的输出端和第四低通滤波器的输出端分别连接至第三除法器的输入端和第四除法器的输入端，第一除法器的输出端和第二除法器的输出端均连接至一阶新幅值信号选择器的输入端，第三除法器的输出端和第四除法器的输出端均连接至二阶新幅值信号选择器的输入端。
与背景技术相比，本发明具有的有益效果是：
(1)本发明方法运用一阶正交幅值信号及一阶正交微分信号来准确提取出载波相位延迟，提取出的相位延迟不受待测相位的影响，可以实现相位延迟的实时提取和补偿；
(2)本发明方法通过在高频正弦电压中附加低频正弦扫描电压实现复合相位调制，使得待测对象静止时的复合相位仍然连续变化，解决了待测相位为特定值时，相位延迟难以求解和补偿的问题；
(3)本发明通过对相位延迟的计算和补偿消除了相位延迟带来的非线性误差，提高了PGC相位解调的精度，可以广泛应用于正弦相位调制干涉技术领域。
附图说明
图1是本发明方法和装置原理框图。
图2是相位延迟补偿模块的原理框图。
图3是本发明仿真实验数据结果图。
图中：1、第一数字频率合成器，2、第一乘法器，3、第二乘法器，4、第一低通滤波器，5、第二低通滤波器，6、第一微分运算器，7、第二微分运算器，8、第一平方运算器，9、第二平方运算器，10、第三平方运算器，11、第四平方运算器，12、第一加法器，13、第二加法器，14、第一开方运算器，15、第二开方运算器、16、第一反正切运算器，17、相位延迟补偿模块，18、第二数字频率合成器，19、第三乘法器，20、第四乘法器，21、第三低通滤波器，22、第四低通滤波器，23、第二反正切运算器，24、相位解包裹处理器，25、滑动平均处理器，26、第一符号提取器，27、第二符号提取器，28、第五乘法器，29、第六乘法器，1701、倍乘器，1702、第一正/余弦运算器，1703、第二正/余弦运算器，1704、第一除法器，1705、第二除法器，1706、第三除法器，1707、第四除法器，1708、一阶新幅值信号选择器，1709、二阶新幅值信号选择器。
具体实施方式
下面结合附图和实施例对本发明加以详细说明。
如图1所示，方法具体实施采用以下PGC相位解调系统，第一乘法器2、第二乘法器3、第三乘法器19和第四乘法器20的输入端均连接数字干涉信号S(t)；第一数字频率合成器1的输出端分别连接至第一乘法器2和第二乘法器3的输入端，第一乘法器2的输出端连接至第一低通滤波器4的输入端，第二乘法器3的输出端连接至第二低通滤波器5的输入端；第一低通滤波器4的输出端分别连接至第二平方运算器9的输入端、第一符号提取器26的输入端、第一微分运算器6的输入端和相位延迟补偿模块17的输入端，第一微分运算器6的输出端连接至第一平方运算器8的输入端；第二低通滤波器5的输出端连接至第三平方运算器10的输入端、第二符号提取器27的输入端、第二微分运算器7的输入端和相位延迟补偿模块17的输入端，第二微分运算器7的输出端连接至第四平方运算器11的输入端；第一平方运算器8的输出端和第二平方运算器9的输出端均连接至第一加法器12的输入端，第三平方运算器10的输出端和第四平方运算器11的输出端均连接至第二加法器13的输入端，第一加法器12的输出端经第一开方运算器14后和第一符号提取器26的输出端一起连接至第五乘法器28的输入端，第二加法器13的输出端经第二开方运算器15后和第二符号提取器27的输出端一起连接至第五乘法器28的输入端，第五乘法器28的输出端和第六乘法器29的输出端均连接至第一反正切运算器16的输入端。
第二数字频率合成器18的输出端连接至第三乘法器19和第四乘法器20的输入端，第三乘法器19和第四乘法器20的输出端分别经第三低通滤波器21和第四低通滤波器22连接至相位延迟补偿模块17的输入端，第一反正切运算器16的输出端连接至相位延迟补偿模块17的输入端，相位延迟补偿模块17的两个输出端均连接至第二反正切运算器23的输入端，第二反正切运算器23输出端连接至相位解包裹处理器24的输入端，相位解包裹处理器24的输出端连接至滑动平均处理器25的输入端，滑动平均处理器25的输出端输出待测相位。
如图2所示，相位延迟补偿模块具体为：第一反正切运算器16的输出端输出θ分别连接至一阶新幅值信号选择器1708的输入端、二阶新幅值信号选择器1709的输入端、第一正/余弦运算器1702的输入端，第一正/余弦运算器1702的两个输出端分别连接至第一除法器1704的输入端和第二除法器1705的输入端，同时第一反正切运算器16的输出端经倍乘器1701连接至第二正/余弦运算器1703的输入端，第二正/余弦运算器1703的两个输出端分别连接至第三除法器1706的输入端和第四除法器1707的输入端，第一低通滤波器4的输出端输出P1和第二低通滤波器5的输出端输出Q1分别连接至第一除法器1704的输入端和第二除法器1705的输入端，第三低通滤波器21的输出端输出P2和第四低通滤波器22的输出端输出Q2分别连接至第三除法器1706的输入端和第四除法器1707的输入端，第一除法器1704的输出端和第二除法器1705的输出端均连接至一阶新幅值信号选择器1708的输入端，第三除法器1706的输出端和第四除法器1707的输出端均连接至二阶新幅值信号选择器1709的输入端。
本发明实施例的实施情况和过程如下：
步骤1：对正弦相位调制干涉仪中的电光相位调制器施加包含高频正弦调制信号和低频正弦扫描信号的复合调制信号，实现对干涉信号相位的复合调制，最终干涉仪输出复合相位调制干涉信号。
步骤2：通过带通滤波器滤除复合相位调制干涉信号中的直流成分和高频噪声，然后对复合相位调制干涉信号进行模数采样得到干涉信号S(t)，其中采样频率大于等于参考载波信号频率的十倍，数字干涉信号S(t)的表达式为：
Figure BDA0002487888390000071
相位
Figure BDA0002487888390000072
包含t时刻的待测相位和扫描相位，公式如下：
Figure BDA0002487888390000073
步骤3：通过第一乘法器2、第二乘法器3、第三乘法器19和第四乘法器20将数字干涉信号S(t)与第一数字频率合成器1产生的一阶正交参考信号和第二数字频率合成器18产生的二阶正交参考信号相乘，经第一低通滤波器4、第二低通滤波器5、第三低通滤波器21和第四低通滤波器22滤波后，获得一阶正交幅值信号(P1，Q1)和二阶正交幅值信号(P2，Q2)，公式如下：
Figure BDA0002487888390000074
Figure BDA0002487888390000081
Figure BDA0002487888390000082
Figure BDA0002487888390000083
然后一阶正交幅值信号(P1，Q1)经过第一微分运算器6和第二微分运算器7进行微分运算后得到正交微分信号(DP，DQ)，公式如下：
Figure BDA0002487888390000084
Figure BDA0002487888390000085
步骤4：幅值信号P1和微分信号DP分别经第一平方运算器8、第二平方运算器9进行平方运算，并通过第一加法器12相加，相加后的结果经第一开方运算器14进行开方运算后通过第五乘法器28与第一符号提取器26的输出相乘，幅值信号Q1和微分信号DQ分别经第三平方运算器10、第四平方运算器11进行平方运算，并通过第二加法器13相加，相加后的结果经第二开方运算器15进行开方运算后通过第六乘法器29与第二符号提取器27的输出相乘，第五乘法器28与第六乘法器29运算得到的结果通过第一反正切运算器16进行反正切运算，构建出不受待测相位影响的相位延迟计算公式，公式如下：
Figure BDA0002487888390000086
计算得到的载波相位延迟θc的取值范围为-π/2到π/2。
其中，附加低频正弦相位调制使得待测物体在静止时相位仍连续变化，使得在待测物体静止时也可以不受初相位影响而准确地计算出相位延迟值；当待测物体开始移动时，构建的计算公式也不会像传统计算方法那样出现因待测相位等于特定值而分子分母同为0的情况。
步骤5：第一反正切运算器16运算得到的载波相位延迟通过第一正/余弦运算器1702和第二正/余弦运算器1703进行正/余弦运算，运算得到的结果分别通过第一除法器1704、第二除法器1705、第三除法器1706、第四除法器1707与第一低通滤波器4、第二低通滤波器5、第三低通滤波器21、第四低通滤波器22输出的结果相除，第一除法器1704和第二除法器1705输出的结果通过一阶新幅值信号选择器1708进行选择，选择后输出不受载波相位延迟影响的一阶新幅值信号R1，第三除法器1706和第四除法器1707输出的结果通过二阶新幅值信号选择器1709进行选择，选择后输出不受载波相位延迟影响的二阶新幅值信号R2，计算公式如下：
Figure BDA0002487888390000091
Figure BDA0002487888390000092
其中，为避免相位延迟修正系数趋近于0从而导致后续计算出现较大误差，根据正/余弦三角函数的特性，限定不同相位延迟下选取正弦项还是余弦项作为相位延迟修正系数；同理，为了避免选取的幅值信号趋近于0，也根据正/余弦三角函数的特性，限定不同相位延迟下选取包含正弦项还是余弦项的幅值信号作为进行后续计算的幅值信号。
步骤6：假定调制深度m＝2.63rad，此时J1(m)＝J2(m)；对一阶新幅值信号选择器1708输出的一阶新幅值信号R1和对二阶新幅值信号选择器1709输出的二阶新幅值信号R2通过第二反正切运算器23进行四象限反正切运算，得到包裹相位
Figure BDA0002487888390000093
计算公式如下：
Figure BDA0002487888390000094
式中，计算得到的包裹相位
Figure BDA0002487888390000095
包裹在-π到+π之间。
步骤7：通过相位解包裹处理器24对包裹相位
Figure BDA0002487888390000096
进行相位解包裹，得到连续变化的复合相位
Figure BDA0002487888390000097
并根据一个低频正弦扫描信号的周期内
Figure BDA0002487888390000098
均值为零的性质，通过采用先入先出的存储结构进行存储复合相位
Figure BDA0002487888390000099
的滑动平均处理器25，对相位解包裹处理器24输出的连续变化的复合相位
Figure BDA00024878883900000910
进行滑动平均运算来消除连续变化的复合相位
Figure BDA00024878883900000911
中的扫描相位，最终得到待测相位
Figure BDA00024878883900000912
公式如下：
Figure BDA0002487888390000101
实际仿真中，按照数字干涉信号S(t)的公式，在MATLAB中产生相同的仿真的复合正弦相位调制干涉信号，其中调制深度设置为2.63，载波相位延迟设置为15°，采用本发明提出的补偿载波相位延迟非线性误差的PGC相位解调方法和传统的PGC-Arctan方法对待测相位进行解调，最终得到如图3所示的实验数据。图3所示的数据中，红线表示未补偿载波相位延迟的PGC-Arctan相位解调算法测量的相位与待测相位的差值(包含非线性误差)。显然该非线性误差随着待测相位呈现正弦规律变化，峰峰值约为6°。蓝线所示为本发明提出的方法测量的相位与待测相位的差值，显然该结果不存在非线性误差并且几乎等于零(误差小于0.02°)。该实验数据表明本发明提出的补偿载波相位延迟非线性误差的PGC相位解调方法可以有效消除由载波相位延迟带来的非线性误差，实现高精度的相位解调。

Claims (3)

1.一种补偿载波相位延迟非线性误差的PGC相位解调方法，其特征在于：
(1)对正弦相位调制干涉仪中的电光相位调制器施加包含高频正弦调制信号和低频正弦扫描信号的复合调制信号，实现对干涉信号相位的复合调制；
(2)通过带通滤波器滤除正弦相位调制干涉仪输出的复合相位调制干涉信号中的直流成分和高频噪声，然后对复合相位调制干涉信号进行模数采样得到正弦相位调制数字干涉信号S(t)，表达式如下：
Figure FDA0002909214140000011
其中，A为正弦相位调制数字干涉信号的幅值，m为相位调制深度，J0(m)为零阶第一类贝塞尔函数，J2n(m)和J2n-1(m)分别为偶数阶和奇数阶第一类贝塞尔函数，n表示阶数，ωc为高频正弦调制信号的角频率，θ为载波相位延迟，t表示时间，
Figure FDA0002909214140000012
表示正弦相位调制数字干涉信号的复合相位；
相位
Figure FDA0002909214140000013
包含t时刻正弦相位调制干涉仪中的待测相位和扫描相位，公式如下：
Figure FDA0002909214140000014
其中，
Figure FDA0002909214140000015
为低频正弦扫描相位，
Figure FDA0002909214140000016
为待测相位，B为扫描相位的幅值，ωs为扫描相位的角频率；
(3)数字频率合成器(1、18)产生的一阶参考信号(cosωct，sinωct)及二阶参考信号(cos2ωct，sin2ωct)分别与正弦相位调制数字干涉信号S(t)相乘，并进行低通滤波，得到一阶正交幅值信号(P1，Q1)和二阶正交幅值信号(P2，Q2)，公式如下：
Figure FDA0002909214140000017
Figure FDA0002909214140000018
Figure FDA0002909214140000019
Figure FDA00029092141400000110
其中，LPF[]表示低通滤波操作；sin(ωct)、cos(ωct)分别表示一阶参考信号的正弦分量和余弦分量，sin(2ωct)、cos(2ωct)分别表示二阶参考信号的正弦分量和余弦分量，P1,Q1分别表示一阶正交幅值信号的余弦幅值分量和正弦幅值分量，P2,Q2分别表示二阶正交幅值信号的余弦幅值分量和正弦幅值分量；
然后一阶正交幅值信号(P1，Q1)经过微分运算后得到一阶正交微分信号(DP，DQ)：
Figure FDA0002909214140000021
Figure FDA0002909214140000022
其中，
Figure FDA0002909214140000023
为复合相位
Figure FDA0002909214140000024
对时间t的偏微分，DP和DQ分别表示一阶正交微分信号的余弦微分分量和正弦微分分量；
(4)利用一阶正交幅值信号(P1，Q1)、一阶正交微分信号(DP，DQ)进行运算，得到载波相位延迟θc，计算公式如下：
Figure FDA0002909214140000025
式中，sign()表示符号函数，当括号中的数值大于等于零时，其值为1，当括号中的数值小于零时，其值为-1；
(5)运用步骤(4)中计算得到的载波相位延迟θc、一阶正交幅值信号(P1，Q1)和二阶正交幅值信号(P2，Q2)，重构出一对幅值不受载波相位延迟影响的新幅值信号(R1，R2)计算公式如下：
Figure FDA0002909214140000026
Figure FDA0002909214140000027
其中，R1和R2分别表示新幅值信号的正弦幅值分量和余弦幅值分量；
(6)对新幅值信号(R1，R2)进行四象限反正切运算，得到包裹相位
Figure FDA0002909214140000028
计算公式如下：
Figure FDA0002909214140000031
(7)对包裹相位
Figure FDA0002909214140000032
进行相位解包裹得到连续变化的复合相位
Figure FDA0002909214140000033
采用长度为M的队列存储复合相位
Figure FDA0002909214140000034
对存储的M个数据进行求和运算，并将求和运算的结果除以M完成滑动平均运算，消除复合相位中的扫描相位，最终得到待测相位
Figure FDA0002909214140000035
公式如下：
Figure FDA0002909214140000036
其中，U[]表示相位解包裹运算，∑[]表示M个数据的求和运算。
2.根据权利要求1所述的一种补偿载波相位延迟非线性误差的PGC相位解调方法，其特征在于：方法采用以下PGC相位解调系统，第一乘法器(2)、第二乘法器(3)、第三乘法器(19)和第四乘法器(20)的输入端均连接数字干涉信号S(t)；第一数字频率合成器(1)的输出端分别连接至第一乘法器(2)和第二乘法器(3)的输入端，第一乘法器(2)的输出端连接至第一低通滤波器(4)的输入端，第二乘法器(3)的输出端连接至第二低通滤波器(5)的输入端；第一低通滤波器(4)的输出端分别连接至第二平方运算器(9)的输入端、第一符号提取器(26)的输入端、第一微分运算器(6)的输入端和相位延迟补偿模块(17)的输入端，第一微分运算器(6)的输出端连接至第一平方运算器(8)的输入端；第二低通滤波器(5)的输出端连接至第三平方运算器(10)的输入端、第二符号提取器(27)的输入端、第二微分运算器(7)的输入端和相位延迟补偿模块(17)的输入端，第二微分运算器(7)的输出端连接至第四平方运算器(11)的输入端；第一平方运算器(8)的输出端和第二平方运算器(9)的输出端均连接至第一加法器(12)的输入端，第三平方运算器(10)的输出端和第四平方运算器(11)的输出端均连接至第二加法器(13)的输入端，第一加法器(12)的输出端经第一开方运算器(14)后和第一符号提取器(26)的输出端一起连接至第五乘法器(28)的输入端，第二加法器(13)的输出端经第二开方运算器(15)后和第二符号提取器(27)的输出端一起连接至第六乘法器(29)的输入端，第五乘法器(28)的输出端和第六乘法器(29)的输出端均连接至第一反正切运算器(16)的输入端；第二数字频率合成器(18)的输出端连接至第三乘法器(19)和第四乘法器(20)的输入端，第三乘法器(19)和第四乘法器(20)的输出端分别经第三低通滤波器(21)和第四低通滤波器(22)连接至相位延迟补偿模块(17)的输入端，第一反正切运算器(16)的输出端连接至相位延迟补偿模块(17)的输入端，相位延迟补偿模块(17)的两个输出端均连接至第二反正切运算器(23)的输入端，第二反正切运算器(23)输出端连接至相位解包裹处理器(24)的输入端，相位解包裹处理器(24)的输出端连接至滑动平均处理器(25)的输入端，滑动平均处理器(25)的输出端输出待测相位。
3.根据权利要求2所述的一种补偿载波相位延迟非线性误差的PGC相位解调方法，其特征在于：所述的相位延迟补偿模块具体为：第一反正切运算器(16)的输出端分别连接至一阶新幅值信号选择器(1708)的输入端、二阶新幅值信号选择器(1709)的输入端、第一正/余弦运算器(1702)的输入端，第一正/余弦运算器(1702)的两个输出端分别连接至第一除法器(1704)的输入端和第二除法器(1705)的输入端，同时第一反正切运算器(16)的输出端经倍乘器(1701)连接至第二正/余弦运算器(1703)的输入端，第二正/余弦运算器(1703)的两个输出端分别连接至第三除法器(1706)的输入端和第四除法器(1707)的输入端，第一低通滤波器(4)的输出端和第二低通滤波器(5)的输出端分别连接至第一除法器(1704)的输入端和第二除法器(1705)的输入端，第三低通滤波器(21)的输出端和第四低通滤波器(22)的输出端分别连接至第三除法器(1706)的输入端和第四除法器(1707)的输入端，第一除法器(1704)的输出端和第二除法器(1705)的输出端均连接至一阶新幅值信号选择器(1708)的输入端，第三除法器(1706)的输出端和第四除法器(1707)的输出端均连接至二阶新幅值信号选择器(1709)的输入端。
CN202010396826.XA 2020-05-12 2020-05-12 一种补偿载波相位延迟非线性误差的pgc相位解调方法 Active CN111595468B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010396826.XA CN111595468B (zh) 2020-05-12 2020-05-12 一种补偿载波相位延迟非线性误差的pgc相位解调方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010396826.XA CN111595468B (zh) 2020-05-12 2020-05-12 一种补偿载波相位延迟非线性误差的pgc相位解调方法

Publications (2)

Publication Number Publication Date
CN111595468A CN111595468A (zh) 2020-08-28
CN111595468B true CN111595468B (zh) 2021-03-30

Family

ID=72185293

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010396826.XA Active CN111595468B (zh) 2020-05-12 2020-05-12 一种补偿载波相位延迟非线性误差的pgc相位解调方法

Country Status (1)

Country Link
CN (1) CN111595468B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112929008A (zh) * 2021-01-25 2021-06-08 思特尔智能检测系统(苏州)有限公司 一种指定幅值波形信号生成方法和系统
CN113098812B (zh) * 2021-03-22 2023-08-18 青岛农业大学 一种4psk相关解调与时钟纠偏方法
CN113375785B (zh) * 2021-06-08 2022-04-15 中国人民解放军国防科技大学 一种光纤水听器全海深高稳定光电信号检测方法
CN113834561B (zh) * 2021-07-30 2023-05-05 西安石油大学 一种pgc相位解调中提取及补偿相位调制深度的系统及方法
CN115834312B (zh) * 2022-11-18 2024-04-16 西安石油大学 基于PGC-Arctan消除载波相位延迟且补偿调制深度的解调方法及装置
CN116679988B (zh) * 2023-08-02 2023-10-27 武汉芯必达微电子有限公司 硬件加速单元、硬件加速方法、芯片及存储介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09203665A (ja) * 1996-01-25 1997-08-05 Oki Electric Ind Co Ltd 光ファイバセンサシステム
CN106850076A (zh) * 2017-03-27 2017-06-13 吉林大学 一种采用模拟前端的数字相位生成载波解调方法及装置
CN109459070A (zh) * 2018-11-15 2019-03-12 浙江理工大学 一种pgc相位解调法中相位延迟提取与补偿方法
CN110307780A (zh) * 2019-06-25 2019-10-08 浙江理工大学 基于迭代计算的pgc相位解调误差实时补偿方法
CN110411486A (zh) * 2019-07-26 2019-11-05 浙江理工大学 对相位延迟和调制深度不敏感的pgc-dcdm解调方法
CN110429987A (zh) * 2019-07-31 2019-11-08 中船海洋探测技术研究院有限公司 一种基于fpga的相位生成载波解调系统

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09203665A (ja) * 1996-01-25 1997-08-05 Oki Electric Ind Co Ltd 光ファイバセンサシステム
CN106850076A (zh) * 2017-03-27 2017-06-13 吉林大学 一种采用模拟前端的数字相位生成载波解调方法及装置
CN109459070A (zh) * 2018-11-15 2019-03-12 浙江理工大学 一种pgc相位解调法中相位延迟提取与补偿方法
CN110307780A (zh) * 2019-06-25 2019-10-08 浙江理工大学 基于迭代计算的pgc相位解调误差实时补偿方法
CN110411486A (zh) * 2019-07-26 2019-11-05 浙江理工大学 对相位延迟和调制深度不敏感的pgc-dcdm解调方法
CN110429987A (zh) * 2019-07-31 2019-11-08 中船海洋探测技术研究院有限公司 一种基于fpga的相位生成载波解调系统

Also Published As

Publication number Publication date
CN111595468A (zh) 2020-08-28

Similar Documents

Publication Publication Date Title
CN111595468B (zh) 一种补偿载波相位延迟非线性误差的pgc相位解调方法
CN109459070B (zh) 一种pgc相位解调法中相位延迟提取与补偿方法
CN110411486B (zh) 对相位延迟和调制深度不敏感的pgc-dcdm解调方法
CN111609792B (zh) 一种pgc相位解调法中相位延迟补偿方法
Xie et al. Extraction of carrier phase delay for nonlinear errors compensation of PGC demodulation in an SPM interferometer
CN110487313B (zh) 光频域反射技术中光源扫频非线性自校正方法
CN110307780A (zh) 基于迭代计算的pgc相位解调误差实时补偿方法
CN108458654B (zh) 基于双通道正交锁相解调的光学非线性误差测量方法及装置
CN105486225A (zh) 一种抑制光强波动噪声的相位解调装置及解调方法
CN102435186B (zh) 一种光纤陀螺的数字信号处理方法、装置及光纤陀螺仪
CN101686087B (zh) 基于3×3耦合器两路输出的波长解调系统及方法
CN111609791B (zh) 一种pgc相位解调法中调制深度的提取与补偿方法
Li et al. A method to remove odd harmonic interferences in square wave reference digital lock-in amplifier
US20120320381A1 (en) Measurement apparatus and measurement method
CN108007340B (zh) 相位生成载波反正切解调中非线性误差的实时计算方法
Link et al. Amplitude and phase measurement of sinusoidal vibration in the nanometer range using laser interferometry
CN105865500B (zh) 一种光纤激光干涉型传感器的移相解调法
CN112462380A (zh) 一种基于激光调频连续波远距离测距的色散补偿方法
CN114964001A (zh) 一种基于混合相位调制的线性化pgc解调信号处理方法
JP3407413B2 (ja) エンコーダのオフセット補正方法
JP7048870B2 (ja) 変位計測装置
CN118551137B (zh) 一种改进的相位解调算法及装置
CN117201250B (zh) 一种相位生成载波解调方法、装置、电子设备及存储介质
CN115993087A (zh) 单频激光干涉测量中一种实时相位延迟补偿方法
CN113834561A (zh) 一种pgc相位解调中提取及补偿相位调制深度的系统及方法

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