四步相移法 matlab

I1=double（l5）;
I2=double（l6）;
I3=double（l7）;
I4=double（l8）;
 [M N]=size（I1）;


%参考平面的
Ii1=double（m5）;
Ii2=double（m6）;
Ii3=double（m7）;
Ii4=double（m8）;
 [Mi Ni]=size（Ii1）;
for i=1:M
     I1（:i）=smooth（I1（:i）30）;
      I2（:i）=smooth（I2（:i）30）;
       I3（:i）=smooth（I3（:i）30）;
        I4（:i）=smooth（I4（:i）30）;
        Ii1（:i）=smooth（Ii1（:i）30）;
        Ii2（:i）=smooth（Ii2（:i）30）;
        Ii3（:i）=smooth（Ii3（:i）30）;
        Ii4（:i）=smooth（Ii4（:i）30）;
 end
 
for ji=1:512
    for ii=2:Mi
phasei（ii