%四步相移輪廓術的仿真實驗
%仿真的步驟分為以下五步

%--------------第一步，三維圖形的模擬-------------------
%函數表達式z（xy）=0.4*sqrt[100*100-（x-150）^2-（y-150）^2]
%模擬的時候以像素作為單位
clear
clc
x?=?1:300;
y?=?1:300;?
[XY]=meshgrid（xy）;
Z=0.4*sqrt（10000-（X-150）.^2-（Y-150）.^2）;
Z（find（angle（Z）））=0;
figure;mesh（xyZ）;
title（‘實際物體的三維輪廓‘）
clear?indexs;


%--------------第二步，正弦光柵的制作-------------------
f?=?1/5;%一個光柵空間頻率
X?=?ones（3001）*x;
fai?=2*pi*f.*X;
I00?=?1?+?cos（fai）;?
imwrite（I00‘shatter00.bmp‘‘bmp‘）;%圖片的寫入，沒有發生移位

%系統參數的設置
len?=?100;%相機的出瞳中心到參考面的距離
d?=20;?%投影中心到照相機中心的距離
p?=1/f;%投影光柵的空間周期，6個像素
detafai?=?2*pi*f*d.*Z/len;%物體表面增加的相位
clear?X?Y?Z?I00;

%變形光柵的寫入
I11?=?1?+?cos（fai+detafai）;%1+cos（）
imwrite（I11‘shatter11.bmp‘‘bmp‘）;%圖片的寫入，第一個光柵
I12?=?1?+?cos（fai+detafai+pi/2）;%?1-sin（）
imwrite（I12‘shatter12.bmp‘‘bmp‘）;%圖片的寫入，第二個光柵
I13?=?1?+?cos（fai+detafai+pi）;%?1-cos（）
imwrite（I13‘shatter13.bmp‘‘bmp‘）;%圖片的寫入，第三個光柵
I14?=?1?+?cos（fai+detafai+3*pi/2）;%1+sin（）
imwrite（I14‘shatter14.bmp‘‘bmp‘）;%圖片的寫入，第四個光柵

%原始光柵的寫入
I01?=?1?+?cos（fai）;%1+cos（）
imwrite（I01‘shatter01.bmp‘‘bmp‘）;%圖片的寫入，第一個光柵
I02?=?1?+?cos（fai+pi/2）;%?1-sin（）
imwrite（I02‘shatter02.bmp‘‘bmp‘）;%圖片的寫入，第二個光柵
I03?=?1?+?cos（fai+pi）;%?1-cos（）
imwrite（I03‘shatter03.bmp‘‘bmp‘）;%圖片的寫入，第三個光柵
I04?=?1?+?cos（fai+3*pi/2）;%1+sin（）
imwrite（I04‘shatter04.bmp‘‘bmp‘）;%圖片的寫入，第四個光柵
clear?I01?I02?I03?I04?I11?I12?I13?I14;

%--------------第三步，圖片的讀取以及初始化濾波-------------------

h=fspecial（‘gaussian‘91）;%產生高斯低通濾波器

%參考圖片與調制圖片的讀取
%參考圖片讀取，為了獲取載波相位
I01=imread（‘shatter01.bmp‘‘bmp‘）;?I01=double（filter2（hdouble（I01）））;
I02=imread（‘shatter02.bmp‘‘bmp‘）;?I02=double（filter2（hdouble（I02）））;
I03=imread（‘shatter03.bmp‘‘bmp‘）;?I03=double（filter2（hdouble（I03）））;
I04=imread（‘shatter04.bmp‘‘bmp‘）;?I04=double（filter2（hdouble（I04）））;
%調制光柵圖片的讀取，為了獲取變形后的相位
I11=imread（‘shatter11.bmp‘‘bmp‘）;?I11=double（filter2（hdouble（I11））