[发明专利]一种基于误差补偿的移相干涉信息处理方法

权利要求书

1.一种基于误差补偿的移相干涉信息处理方法，其特征是，包括以下步骤，

第一步，确定移相干涉检测的误差源的敏感方程；

所述误差源为移相不准误差、探测器的响应误差、光源的光强不稳定和频率不稳定；

各种误差源的敏感方程如下：

移相不准误差，实际移相量和移相不准可用理想移相量的多项式表示为

δn′=(1+ϵ1)δn+ϵ2(δn)2+ϵ3(δn)3+...·Δδn=ϵ1δn+ϵ2(δn)2+ϵ3(δn)3+...·=ϵ1(n-1)δ+ϵ2[(n-1)δ]2+ϵ3[(n-1)δ]3+...·---(8)]]>

对待测相位等式求偏微分，并带入待测相位等式，可得到移相不准的敏感方程

Δφ=cos2(φ)Σn=1Nwnsin2(δn)Δδn+sin2(φ)Σn=1Nwncos2(δn)Δδn+2sin(φ)cos(φ)Σn=1Nwnsin(δn)cos(δn)Δδn---(9)]]>

探测器响应误差，探测器存在二阶非线性响应误差时，探测器输出的干涉图信号与干涉 图光强信号之间的关系以及响应ΔIn偏差可表示为，

I′_n(x,y)＝I_n(x,y)+ηI_n²(x,y)

ΔI_n＝I′_n-I_n＝η[A²+2ABcos(φ+δ_n)+B²cos²(φ+δ_n)](10)

对待测相位等式求偏微分，并带入等式(10)，可得到探测器响应非线性的敏感方程：

Δφ=-1Bcos(φ)ηA2Σn=1Nwnsin(δn)+sin(φ)ηA2Σn=1Nwncos(δn)+2cos2(φ)ηABΣn=1Nwnsin(δn)cos(δn)-sin2(φ)ηABΣn=1Nwncos(δn)sin(δn)+2sin(φ)cos(φ)ηABΣn=1Nwncos2(δn)-2sin(φ)ηABΣn=1Nwnsin2(δn)+cos(φ)ηB22Σn=1Nwnsin(δn)+cos(φ)ηB2cos(2φ)2Σn=1Nwnsin(δn)cos(2δn)cos(φ)ηB2sin(2φ)2Σn=1Nwnsin(δn)sin(2δn)+sin(φ)ηB22Σn=1Nwncos(δn)+sin(φ)ηB2cos(2φ)2Σn=1Nwncos(δn)cos(2δn)sin(φ)ηB2sin(2φ)2Σn=1Nwncos(δn)sin(2δn)---(11)]]>

光源的光强不稳定和频率不稳定，存在光强波动时，干涉图可以用帧数的多项式表示， 第n帧干涉图以及探测器输出的干涉图信号与理想干涉图光强信号之间的偏差可表示为可分 别表示为，

In′(x,y)=[1+ϵ1(n-1)+ϵ2(n-1)2+......+ϵk(n-1)k]In(x,y)ΔIn′=In′-In=[ϵ1(n-1)+ϵ2(n-1)2+......+ϵk(n-1)k]In---(12)]]>

对待测相位等式求偏微分，并带入等式(12)，可得到光强不稳定的敏感方程：

当仅存在激光器的频率波动时，引入的额外相位差可表示为

Δδ=[(OPD+ΔOPD)×(υ+Δυ)C-OPD×υC)×2π≈[OPD×ΔυC]×2π]]>

这种移相不准与激光器的频率波动大小成正比，当存在频率波动时，频率波动引起的移相不 准可用帧数的多项式表示，

Δδ_n＝ε₁(n-1)δ+ε₂(n-1)²δ+ε₃(n-1)³δ+…(14)

对待测相位等式求偏微分，并带入等式(14)，可得到光源频率不稳定引入的测量误差ΔΦ：

Δφ=cos2(φ)Σn=1Nwnsin2(δn)Δδn+sin2(φ)Σn=1Nwncos2(δn)Δδn+2sin(φ)cos(φ)Σn=1Nwnsin(δn)cos(δn)Δδn---(15);]]>

第二步：基于权重待定最小二乘法和敏感方程，获取约束方程组；

综合敏感方程(9)、(11)、(13)、和(15)，要实现移相算法对上述误差源不敏感，需要 满足约束条件为，

Σn=1Nwn=1Σn=1Nwnsin(δn)=Σn=1Nwncos(δn)=0Σn=1Nwnsin(2δn)=Σn=1Nwncos(2δn)=0Σn=1N(n-1)wnsin(2δn)=Σn=1N(n-1)2wnsin(2δn)=...=0Σn=1N(n-1)wncos(2δn)=Σn=1N(n-1)2wncos(2δn)=...=0Σn=1N(n-1)wnsin(δn)=Σn=1N(n-1)2wnsin(δn)=...=0Σn=1N(n-1)wncos(δn)=Σn=1N(n-1)2wncos(δn)=...=0Σn=1Nwnsin(δn)cos(2δn)=Σn=1Nwncos(δn)cos(2δn)=0Σn=1Nwncos(δn)sin(2δn)=Σn=1Nwnsin(δn)sin(2δn)=0---(16);]]>

第三步：解约束方程组，确定权重，获取对误差不敏感的相位提取算法；

当误差源伟PZT的移相不准、CCD的二阶响应非线性和光源的不稳定性不敏感，权重为 正实数且对称分布，移相量为π/2的等间隔多步移相，移相步数为13，解方程组获得权重如 下

w1=w13=1256;w2=w12=4256]]>

w3=w11=10256;w4=w10=20256]]>

w5=w9=31256;w6=w8=40256]]>

w7=44256---(17)]]>

确定加权最小二乘相位提取算法为

φ=arctan[-4(I~2-I~12)-20(I~4-I~10)+40(I~6-I~8)(I~1+I~13)-10(I~3+I~11)+31(I~5+I~9)-44I~7]---(18).]]>