%%初始化參數
wid=500;%圖片寬度
height=1000;%圖片高度
T=40;%條紋周期
a=100;%正弦幅值
b=120;%正弦y軸平移值
%四幅相位差的圖片
I1=zeros（widheight）;
I2=zeros（widheight）;
I3=zeros（widheight）;
I4=zeros（widheight）;
%%?產生四步相移光柵條紋
for?i=1:wid
????for?k=1:height
????????I1（ik）=a*cos（2*pi*k/T+0）+b;
????????I2（ik）=a*cos（2*pi*k/T+pi/2）+b;
????????I3（ik）=a*cos（2*pi*k/T+pi）+b;
????????I4（ik）=a*cos（2*pi*k/T+1.5*pi）+b;
????end
end
figure（1）;%顯示四步相移光柵
subplot（221）;imshow（I1[]）;title（‘初始相位0‘）;
subplot（222）;imshow（I2[]）;title（‘初始相位0.5pi‘）;
subplot（223）;imshow（I3[]）;title（‘初始相位pi‘）;
subplot（224）;imshow（I4[]）;title（‘初始相位1.5pi‘）;
%%?四步相移解相
I=atan（（I4-I2）./（I1-I3））+pi;%待修正的四步相移提取相位
figure（2）;
imshow（I[]）;title（‘提取相位圖‘）;
?%????????imwrite（I‘C:\Users\nec\Desktop\標柵提取相位圖.bmp‘‘bmp‘）;%imwrite有問題，保存的圖像是二值圖
%?顯示解出來的一列
figure（3）;
plot（1:heightI（3:））;grid?on;title（‘相位圖結果分析（第三列）‘）;
hold?on;
figure（4）;mesh（I）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1045??2013-05-31?10:36??zhou3.m