相位跟踪算法计算条纹图相位

function [thitan0mm] = regulation（xymthitan0）
%UNtitleD3 此处显示有关此函数的摘要
%   此处显示详细说明
%% 定义初值
Nx=7;  %%领域的大小
n=（Nx-1）/2;
lamda=2; %% 控制调制项的参数 1或者2
tao=10;%%τ的范围是在（0.110）之间

%% 设置Ur的初值，即在（xy）处的梯度

k=1;%%k表示循环次数，k=1是初值
errf=1;
while errf>1e-8&&k<5
for kesi=（x-n+1）:（x+n-1） %%x表示行，y表示列
    for yita =（y-n+1）:（y+n-1）
        
        p（kesiyita）=thitan0（1）+thitan0（2）*（x-kesi）+thitan0（3）*（y-yita）;
        U_phi0（kesiyita）=2*（In（kesiyita）-cos（p（kesiyita）））*sin（p（kesiyita））-2*lamda*（thitan0（1）-p（kesiyita）