基于迭代计算的PGC相位解调误差实时补偿方法与流程

本发明属于激光干涉测量技术领域，尤其是一种基于迭代计算的pgc相位解调误差实时补偿方法。

背景技术：

相位生成载波(pgc)调制解调技术具有灵敏度高、动态范围大、抗干扰能力强等优点，被广泛应用于干涉型光纤传感器、光学反馈干涉仪等激光干涉测量系统。在pgc相位解调技术中，通常将干涉信号与本地参考信号相乘并经低通滤波处理后得到一对正交信号，进一步采用微分交叉相乘(pgc-dcm)或反正切(pgc-arctan)运算获得待测相位。其中pgc-dcm算法容易受干涉信号幅值波动、干涉条纹对比度变化、载波相位延迟和调制深度的影响，存在线性误差；pgc-arctan算法对此做了一定改进，消除了干涉信号幅值波动和干涉条纹对比度变化的影响，但是仍然会受载波相位延迟、调制深度的影响，导致正交信号幅值不相等，从而引入非线性误差。因此，正交信号的幅值的实时计算与修正是补偿pgc相位解调非线性误差与提高测量精度的关键技术问题。

技术实现要素：

为了解决背景技术中存在的问题，本发明公开了一种基于迭代计算的pgc相位解调误差的实时补偿方法，本发明解决了载波相位延迟、调制深度引起的非线性误差难以同时实时补偿的问题。

本发明解决其技术问题所采用的技术方案包括以下步骤：

一、一种基于迭代计算的pgc相位解调误差实时补偿系统：

包括频率综合器、第一乘法器、第一低通滤波器、基于迭代计算的修正系数计算模块、第一除法器、第二乘法器、第二低通滤波器、第二除法器和反正切运算器；干涉信号输入到至第一乘法器和第二乘法器的输入端，频率综合器的两个输出端分别连接至第一乘法器和第二乘法器的输入端，第一乘法器的输出端连接至第一低通滤波器的输入端，第二乘法器的输出端连接至第二低通滤波器的输入端，第一低通滤波器的输出端分别连接至基于迭代计算的修正系数计算模块的输入端和第一除法器的输入端，第二低通滤波器的输出端分别连接至基于迭代计算的修正系数计算模块的输入端和第一除法器的输入端，基于迭代计算的修正系数计算模块的两个输出端分别连接至第一除法器和第二除法器的输入端，第一除法器和第二除法器的输出端分别连接至反正切运算器的输入端，反正切运算器输出待测相位。

所述的基于迭代计算的修正系数计算模块包括平方向量计算模块、预测误差计算模块、增益矩阵计算模块、参数向量预测值计算模块、幅值修正系数计算模块；第一低通滤波器和第二低通滤波器的输出端连接到平方向量计算模块的输入端，平方向量计算模块的输出端分别连接到预测误差计算模块、增益矩阵计算模块和参数向量预测值计算模块的输入端，预测误差计算模块和增益矩阵计算模块的输出端连接到参数向量预测值计算模块的输入端，参数向量预测值计算模块的输出端反馈连接到参预测误差计算模块的输入端并连接到幅值修正系数计算模块的输入端，幅值修正系数计算模块连接到第一除法器和第一除法器。

二、一种基于迭代计算的pgc相位解调误差实时补偿方法：

步骤1)干涉信号经模数转换后进入现场可编程逻辑门阵列信号处理器(fpga)，与频率综合器产生的本地参考信号经乘法器相乘并经低通滤波器低通滤波处理，计算获得一对包含待测相位的正交信号p和q：

a＝rj1(m)cosθ

b＝rj2(m)cos2θ

其中，a和b分别表示正交信号p和q的幅值，r为干涉信号交流分量幅值，j1(m)和j2(m)分别为第一类的一阶和二阶贝塞尔函数，m表示调制深度，θ表示载波相位延迟，φ(t)表示待测相位，t表示时间；

步骤2)平方向量计算模块401中对正交信号p和q进行平方运算，进一步获得平方向量h＝(p²,q²,1)；结合正交信号p和q的幅值定义参数向量x＝(r²/a²,r²/b²,-r²)，r为干涉信号交流分量幅值，参数向量x为待求解的未知参数，平方向量和参数向量满足下列关系式：

其中，上标t表示转置；

结合平方向量h对参数向量x进行预测，k表示预测的次数，k＝0,1,2,3……，参数向量x的预测值表示为xk＝(ak,bk,ck)，其中ak,bk,ck分别表示对r²/a²、r²/b²、-r²的预测值，参数向量x的预测值初始为x0＝(1,1,-1)，第k次预测时先通过平方向量计算模块401计算平方向量hk，预测误差计算模块402中结合平方向量hk与参数向量的预测值xk计算第k次预测时的误差，得到预测误差ek：

步骤3)结合平方向量hk在增益矩阵计算模块403中处理计算求得增益矩阵gk：

gk＝2nkhk^t(hknkhk^t+(2akpk)²+(2bkqk)²)^-1

其中，i为单位矩阵，nk表示中间矩阵，其初始值等于单位矩阵i；

步骤4)参数向量预测值计算模块404运用第k次的预测误差ek、增益矩阵gk以及参数向量预测值xk计算获得第k+1次参数向量预测值xk+1，并且每次参数向量预测值的计算结果都返回预测误差计算模块402：

xk+1＝xk-gkek

步骤5)重复上述步骤进行多次迭代计算，参数向量的预测值xk将逐步逼近参数向量的真实值；当相邻两次参数向量预测值xk与xk-1的差值小于预测阈值时，幅值修正系数计算模块405中取最后的参数向量预测值xk采用以下公式计算正交信号p和q对应的幅值修正系数a*和b*：

步骤6)在除法器中采用以下公式根据幅值修正系数a*和b*对正交信号p和q进行修正，获得修正的正交信号p*和q*：

修正的正交信号p*、q*的幅值均为1，进一步根据下式在反正切运算器中获得补偿后的待测相位实现pgc相位解调误差实时补偿：

本发明具有的有益效果是：

(1)本发明方法采用迭代计算来估计pgc相位解调中正交信号的幅值修正系数，并进一步将幅值修正为1，补偿了载波相位延迟和调制深度引起非线性误差，提高了pgc相位解调精度。

(2)本发明的方法完全在现场可编程逻辑门阵列信号处理器(fpga)内实现，运算速度快，实时性高，当载波相位延迟和调制深度波动时能快速计算出新的幅值修正系数，实现了对载波相位延迟和调制深度引起非线性误差的实时补偿。

附图说明

图1是基于迭代计算的pgc相位解调误差实时补偿方法的原理框图。

图2是基于迭代计算的修正系数计算模块的原理框图。

图中：1、频率综合器，2、第一乘法器，3、第一低通滤波器，4、基于迭代计算的修正系数计算模块，5、第一除法器，6、第二乘法器，7、第二低通滤波器，8、第二除法器，9、反正切运算器。

401、平方向量计算模块，402、预测误差计算模块，403、增益矩阵计算模块，404、参数向量预测值计算模块，405、幅值修正系数计算模块。

具体实施方式

下面结合附图对本发明加以详细说明。

如图1所示，本发明系统具体实施采用现场可编程逻辑门阵列信号处理器(fpga)，现场可编程逻辑门阵列信号处理器(fpga)包括频率综合器1、第一乘法器2、第一低通滤波器3、基于迭代计算的修正系数计算模块4、第一除法器5、第二乘法器6、第二低通滤波器7、第二除法器8和反正切运算器9；干涉信号s(t)输入到至第一乘法器2和第二乘法器6的输入端，频率综合器1的两个输出端分别连接至第一乘法器2和第二乘法器6的输入端，频率综合器1用于产生参考信号，第一乘法器2的输出端连接至第一低通滤波器3的输入端，第二乘法器6的输出端连接至第二低通滤波器7的输入端，第一低通滤波器3的输出端分别连接至基于迭代计算的修正系数计算模块4的输入端和第一除法器5的输入端，第二低通滤波器7的输出端分别连接至基于迭代计算的修正系数计算模块4的输入端和第一除法器8的输入端，基于迭代计算的修正系数计算模块4的两个输出端分别连接至第一除法器5和第二除法器8的输入端，第一除法器5和第二除法器8的输出端分别连接至反正切运算器9的输入端，反正切运算器9输出待测相位。

如图2所示，基于迭代计算的修正系数计算模块4包括平方向量计算模块401、预测误差计算模块402、增益矩阵计算模块403、参数向量预测值计算模块404、幅值修正系数计算模块405；第一低通滤波器3和第二低通滤波器7的输出端分别连接到预测平方向量计算模块401的输入端，平方向量计算模块401的输出端分别连接至预测误差计算模块402、增益矩阵计算模块403和参数向量预测值计算模块404的输入端，预测误差计算模块402和增益矩阵计算模块403的输出端分别连接到参数向量预测值计算模块404的输入端，参数向量预测值计算模块404的输出端反馈连接到预测误差计算模块402的输入端并连接到幅值修正系数计算模块405，幅值修正系数计算模块405连接到第一除法器5和第一除法器8。

本发明实施过程如下：

如图1所示为在现场可编程逻辑门阵列信号处理器(fpga)内实现的基于迭代计算的pgc相位解调误差实时补偿方法的原理框图，正弦相位调制干涉仪的干涉信号s(t)经模数转换后进入现场可编程逻辑门阵列信号处理器(fpga)，干涉信号表示如下：

其中r、c分别表示干涉信号的交流分量和直流分量幅值，m表示调制深度，ω表示相位调制的角频率，θ表示载波相位延迟，表示待测相位，t表示时间。

通过第一乘法器2、第二乘法器6将干涉信号s(t)分别与频率综合器1产生的本地参考信号cos(ωt)和cos(2ωt)相乘，分别通过第一低通滤波器3和第二低通滤波器7处理后，获得一对含有待测相位的正交信号p和q，分别表示如下：

其中，a＝rj1(m)cosθ、b＝rj2(m)cos2θ分别表示正交信号p和q的幅值。理想情况下，调制深度m等于2.63rad，载波相位延迟θ等于零，此时a＝b，但在通常情况下，调制深度m会偏离2.63rad，载波相位延迟θ会偏离零，即此时a≠b，最终在求解待测相位时会引入非线性误差。

将正交信号p和q输入到基于迭代计算的修正系数计算模块4，进行多次迭代计算得到幅值修正系数a*和b*。通过第一除法器5与第二除法器8分别将正交信号p和q与幅值修正系数a*和b*相除，得到修正后的正交信号p*和q*，公式如下：

修正后的正交信号p*和q*的幅值均为1，使用反正切运算器9对修正后的正交信号p*和q*进行反正切运算，求得待测相位

如图2所示，为基于迭代计算的修正系数计算模块4的原理框图，即对图1中基于迭代计算的修正系数计算模块4内信号处理方法进一步说明：

1)在平方向量计算模块401中对正交信号p和q进行平方运算，进一步获得平方向量h＝(p²,q²,1)；结合正交信号p和q定义参数向量x＝(r²/a²,r²/b²,-r²)，参数向量x为待求解的未知参数。

结合正交信号p和q的数学关系可知其满足如下关系式：

其中，上标t表示转置。

基于迭代计算的修正系数计算模块4在每次输入新的正交信号p和q时都会对参数向量x进行预测。k表示预测的次数(k＝0,1,2,3……)，则第k次预测的参数向量预测值表示为xk＝(ak,bk,ck)，其中ak,bk,ck分别表示r²/a²,r²/b²,-r²的预测值，预测初始值为x0＝(1,1,-1)。第k次预测时输入的正交信号表示为pk和qk，平方向量计算模块401运用正交信号pk和qk计算得到第k次的平方向量hk：

hk＝(pk²,qk²,1)(8)

2)预测误差计算模块402运用平方信号hk以及参数向量预测值xk计算得到预测误差ek：

3)增益矩阵计算模块403运用第k次的平方向量hk和正交信号pk和qk计算得到增益矩阵gk：

gk＝2nkhk^t(hknkhk^t+(2akpk)²+(2bkqk)²)^-1(10)

其中，i为单位矩阵，nk表示中间矩阵，其初始值等于单位矩阵i；

4)参数向量预测值计算模块404运用第k次的预测误差ek、增益矩阵gk以及参数向量预测值xk递推计算得到新的参数向量预测值xk+1：

xk+1＝xk-gkek(12)

其中，参数向量预测值的初始值为x0＝(1,1,-1)。

5)通过对步骤1)、2)、3)、4)进行迭代计算，迭代计算时，参数向量预测值将逐步逼近参数向量的真实值。在幅值修正系数计算模块405中对相邻两次参数向量预测值xk与xk-1的差值进行检测，当差值小于一定阈值时取参数向量预测值xk来计算得到正交信号p和q对应的幅值修正系数a*和b*，表示如下：

至此，基于迭代计算的修正系数计算模块4完成幅值修正系数的快速计算，并输出幅值修正系数用于pgc相位解调误差的实时补偿。

综上，本发明方法采用迭代计算来获得pgc相位解调中正交信号的幅值，并进一步将幅值修正为1，补偿了载波相位延迟和调制深度引起非线性误差，提高了pgc相位解调精度。并且本方法完全在现场可编程逻辑门阵列信号处理器(fpga)内实现，运算速度快，实时性高，当载波相位延迟和调制深度波动时能快速计算出幅值修正系数，实现了对载波相位延迟和调制深度引起非线性误差的实时补偿。

上述具体实施方式用来解释说明本发明，而不是对本发明进行限制，在本发明的精神和权利要求的保护范围内，对本发明做出的任何修改和改变，都落入本发明的保护范围。