Uo=imread（‘guang.bmp‘）;??????????%調入作為物的圖像
Uo=double（Uo（::1））;????????????%取第一層，并轉為雙精度
figureimshow（Uo[]）title（‘原始像‘）
[rc]=size（Uo）;
lamda=6328*10^（-10）;k=2*pi/lamda;%賦值波長、波數
zo=0.3086;???????????????????????%物到全息記錄面的距離單位:米
Lo=5*10^（-3）?????????????????????%賦值衍射面（物）的尺寸單位:米
%=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?
%用T-FFT算法完成全息圖記錄過程的計算
xo=linspace（-Lo/2Lo/2c）;yo=linspace（-Lo/2Lo/2r）;
[xoyo]=meshgrid（xoyo）;
F0=exp（j*k*zo）/（j*lamda*zo）;
F1=exp（j*k/2/zo.*（xo.^2+yo.^2））;
fa=fft2（Uo）;?fF1=fft2（F1）;
Fuf=fa.*fF1;?
O=F0.*fftshift（ifft2（Fuf））;?????%在全息記錄面上的光場分布
I=O.*conj（O）;???????????????????%全息記錄面上的光強分布
figureimshow（I[]）colormap（pink）title（‘衍射圖‘）
%下面加入參考光
alpha=pi/2.00;??????????????????%參考光與x軸間的夾角
beita=pi/2.02;??????????????????%參考光與y軸間的夾角
R=exp（j*k*（xo*cos（alpha）+yo*cos（beita）））;?%參考光
%下面計算參、物光在全息記錄面上的干涉得到全息圖
inter=O./max（max（sqrt（I）））+R;???%調節光束比，并使參、物光干涉
II=?inter.*conj（inter）;?????????%干涉得到全息圖
figureimshow（II[]）title（‘全息圖‘）
%=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?=?
%下面計算再現
%先用S-FFT算法再現全息像
zi=0.3086;??????????????????????%全息圖到觀察面的距離單位:米
Li=r*lamda*zi/Lo????????????????%給出像面的尺寸單位:米
x=linspace（-Li/2Li/2c）;y=linspace（-Li/2Li/2r）;
[xy]=meshgrid（xy）;
F0=exp（j*k*zi）/（j*lamda*zi）*exp（j*k/2/zi*（x.^2+y.^2））;
F=exp（j*k/2/zi*（xo.^2+yo.^2））;??%用T-FFT算法得到的全息圖尺寸與物面一致
%?取再現照明光垂直入射C=1
holo=Lo/r*Lo/c*fftshift（fft2（II.*F*1））;?holo=holo.*F0;
Ii=holo.*conj（holo）;
figureimshow（Ii[0max（max（Ii））./1]）colormap（pink）title（‘S-FFT再現像‘）

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

?????文件??????21345??2016-05-17?15:58??holographic?recording?and?playback\1.jpg

?????文件??????36965??2016-05-17?15:58??holographic?recording?and?playback\2.jpg

?????文件??????44976??2016-05-17?15:58??holographic?recording?and?playback\3.jpg

?????文件??????20242??2016-05-17?15:58??holographic?recording?and?playback\4.jpg

?????文件???????1795??2014-02-12?11:41??holographic?recording?and?playback\CX11_1.m

?????文件???????1290??2014-02-12?11:42??holographic?recording?and?playback\CX11_2.m

?????文件???????1323??2014-02-12?12:09??holographic?recording?and?playback\CX11_2_1.m

?????文件???????1378??2014-02-12?12:17??holographic?recording?and?playback\CX11_2_2.m

?????文件?????786486??2008-10-26?11:59??holographic?recording?and?playback\guang.bmp

?????文件????9437238??2011-06-03?22:00??holographic?recording?and?playback\ho.bmp

?????目錄??????????0??2016-05-17?16:03??holographic?recording?and?playback

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

?????????????10353038????????????????????11