Uo=imread（‘guang.bmp‘）;?????????????%調(diào)入作為物的圖像
Uo=double（Uo?（::1））;??????????????%取第一層，并轉(zhuǎn)為雙精度
[rc]=size（Uo）;
Uo=ones（rc）*0.98-Uo/255*0.5;???????%將物轉(zhuǎn)換為高透射率射系數(shù)體
figureimshow（Uo[01]）title（‘物‘）
lamda=6328*10^（-10）;k=2*pi/lamda;???%賦值波長(zhǎng)和波數(shù)
Lo=5*10^（-3）????????????????????????%賦值衍射面（物）的尺寸
xo=linspace（-Lo/2Lo/2r）;yo=linspace（-Lo/2Lo/2c）;
[xoyo]=meshgrid（xoyo）;????????????%生成衍射面（物）的坐標(biāo)網(wǎng)格
zo=0.20;????????????????????????????%全息記錄面到衍射面的距離單位:米
%下面用T-FFT算法完成物面到全息記錄面的衍射計(jì)算
F0=exp（j*k*zo）/（j*lamda*zo）;
F1=exp（j*k/2/zo.*（xo.^2+yo.^2））;
fF1=fft2（F1）;
fa1=fft2（Uo）;
Fuf1=fa1.*fF1;?
Uh=F0.*fftshift（ifft2（Fuf1））;?
Ih=Uh.*conj（Uh）;
figureimshow（Ih[0max（max（Ih））/1]）title（‘全息圖‘）
%下面用T-FFT算法完成全息面到觀察面的衍射計(jì)算（重構(gòu)再現(xiàn)像）
zi=0.20????????????????????????????%賦值再現(xiàn)距離（可以調(diào)整）
F0i=exp（j*k*zi）/（j*lamda*zi）;
F1i=exp（j*k/2/zi.*（xo.^2+yo.^2））;??%T-FFT算法，物面、全息圖和再現(xiàn)像尺寸相同
fF1i=fft2（F1i）;
fIh=fft2（Ih）;?
FufIh=fIh.*fF1i;?
Ui=F0i.*fftshift（ifft2（FufIh））;?
Ii=Ui.*conj（Ui）;
figureimshow（Ii[0max（max（Ii））/1]）title（‘再現(xiàn)像‘）

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件??????14601??2016-05-07?13:11??Bo?coaxial?holographic\1.jpg

?????文件??????34556??2016-05-07?13:11??Bo?coaxial?holographic\2.jpg

?????文件??????35657??2016-05-07?13:11??Bo?coaxial?holographic\3.jpg

?????文件???????1180??2014-02-12?11:37??Bo?coaxial?holographic\CX10_1.m

?????文件???????1517??2014-02-12?11:51??Bo?coaxial?holographic\CX10_1_1.m

?????文件?????786486??2008-10-26?11:59??Bo?coaxial?holographic\guang.bmp

?????目錄??????????0??2016-05-07?13:14??Bo?coaxial?holographic

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

???????????????873997????????????????????7