四步相移中，相移误差对相位求解的影响

1.四步相移公式：

$I_0(x,y)=A(x,y)+B(x,y)cos[\phi (x,y)]$
$I_1(x,y)=A(x,y)+B(x,y)cos[\phi (x,y)+\pi /2]$
$I_2(x,y)=A(x,y)+B(x,y)cos[\phi (x,y)+2\pi /2]$
$I_3(x,y)=A(x,y)+B(x,y)cos[\phi (x,y)+3\pi /2]$

2.变形：

$I_0(x,y)=A(x,y)+B(x,y)cos[\phi (x,y)]$
$I_1(x,y)=A(x,y)+B(x,y)sin[\phi (x,y)]$
$I_2(x,y)=A(x,y)-B(x,y)cos[\phi (x,y)]$
$I_3(x,y)=A(x,y)-B(x,y)sin[\phi (x,y)]$

3.求解相位：

$\phi (x,y)=arctan(\frac{I_1-I_3}{I_0-I_2})$

4.引入误差：

一般地，没有一个系统可以保证无限精确的相移，实际的相移量会和理想的有一个微小的差异，实际的相移为：
${\delta }_i^{{}^{\prime }}={\delta }_i+{\epsilon }_i$

5.光强分布为：

$I_0(x,y)=A(x,y)+B(x,y)cos[\phi (x,y)]$
$I_1^{{}^{\prime }}(x,y)=A(x,y)+B(x,y)sin[\phi (x,y)+{\epsilon }_1]$
$I_2^{{}^{\prime }}(x,y)=A(x,y)-B(x,y)cos[\phi (x,y)+{\epsilon }_2]$
$I_3^{{}^{\prime }}(x,y)=A(x,y)-B(x,y)sin[\phi (x,y)+{\epsilon }_3]$

则：

$\frac{I_1^{{}^{\prime }}-I_3^{{}^{\prime }}}{I_0^{{}^{\prime }}-I_2^{{}^{\prime }}}$
$=\frac{A+Bsin[\phi (x,y)+{\epsilon }_1]-A+Bsin[\phi (x,y)+{\epsilon }_3]}{A+Bsin[\phi (x,y)]-A+Bsin[\phi (x,y)+{\epsilon }_2]}$
$=\frac{sin\phi cos{\epsilon }_1+cos\phi sin{\epsilon }_1+sin\phi cos{\epsilon }_3+cos\phi sin{\epsilon }_3}{cos\phi +cos\phi cos{\epsilon }_2-sin\phi sin{\epsilon }_2}$
$=\frac{sin\phi (cos{\epsilon }_1+cos{\epsilon }_3)+cos\phi (sin{\epsilon }_1+sin{\epsilon }_3)}{cos(\phi (1+cos{\epsilon }_2)-sin\phi sin{\epsilon }_2}$
$=\frac{tan\phi (cos{\epsilon }_1+cos{\epsilon }_3)+(sin{\epsilon }_1+sin{\epsilon }_3)}{1+cos{\epsilon }_2-tan\phi sin{\epsilon }_2}$
其中，${\epsilon }_i$一般较小，有$sin{\epsilon }_i\approx \epsilon ,cos{\epsilon }_i\approx 1$

则：

$tan\phi =\frac{2(I_1^{{}^{\prime }}-I_3^{{}^{\prime }})-({\epsilon }_1+{\epsilon }_3)(I_0^{{}^{\prime }}-I_2^{{}^{\prime }})}{2(I_0^{{}^{\prime }}-I_2^{{}^{\prime }})+{\epsilon }_2(I_1^{{}^{\prime }}-I_3^{{}^{\prime }})}$

即：

$\phi =arctan\frac{2(I_1^{{}^{\prime }}-I_3^{{}^{\prime }})-({\epsilon }_1+{\epsilon }_3)(I_0^{{}^{\prime }}-I_2^{{}^{\prime }})}{2(I_0^{{}^{\prime }}-I_2^{{}^{\prime }})+{\epsilon }_2(I_1^{{}^{\prime }}-I_3^{{}^{\prime }})}$

总结：

四步相移的核心方法是将干涉图间的相移量调整为π/2，从而在计算中消去直流项，约去常数，方便求解出相位值。有额外的相移量（误差）时，求解公式应作出调整（如最后一个公式）。如果使用原来的公式，求解出的相位会有偏差。

参考文献

[1]王吉元,田爱玲,王红军.四步相移干涉测量的相位补偿及仿真[J].西安工业学院学报, 2005, 25(1):5.DOI:10.3969/j.issn.1673-9965.2005.01.007.
[1]蔡长青,贺玲凤.基于四步相移的相位差提取方法[J].华南理工大学学报：自然科学版, 2011, 39(9):4.DOI:10.3969/j.issn.1000-565X.2011.09.016.