%菲涅耳算法計算同軸菲涅耳全息圖的再現象
%lambda=632.8nm再現距離z，dxh=dyh=8.33um512*512像素
close?all;
clear;
lamda=632.8*10^-9;?????????%?首先要計算機模擬再現光，使之與全息圖相乘，得到重構波前在全息圖平面的復振幅分布，再利用菲涅爾變換法重建
z=-110*10^-3;???????????????%再現距離
k=2*pi/lamda;??????????????%波矢量
A=exp（1i*k*z）/（1i*lamda*z）;%相位因子
dxh=8.33*10^-6;????????????%像素大小
dyh=8.33*10^-6;
xh=ones（5121）*（-256:255）*dxh;
yh=（ones（5121）*（-256:255）*dyh）‘;
dxi=lamda*z/（512*dxh）;????%再現像素大小
dyi=lamda*z/（512*dyh）;
xi=ones（5121）*（-256:255）*dxi;
yi=（ones（5121）*（-256:255）*dyi）‘;
cosx=0.0180;??????????????%夾角的余弦值
cosy=0.0180;
I=double（rgb2gray（imread（‘512p1.bmp‘）））;
ref=exp（-1i*k*（xh*cosx+yh*cosy））;??????%再現光
FI=fftshift（fft2（I））;
mesh（log（abs（FI）））;
%imwrite（mat2gray（log（abs（FI）））‘23.bmp‘）;
%figureimagesc（log（abs（FI）））colormap（gray）;
a=zeros（5

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1589??2012-05-05?16:18??h\fresnell.m

?????文件????????965??2012-05-09?11:22??h\jiaopu.m

?????目錄??????????0??2012-05-09?11:59??h

-----------?---------??----------?-----??----

?????????????????2554????????????????????3