計算傅里葉變換全息

《計算傅里葉變換全息》由會員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。

1、計算傅里葉變化全息原圖（秦時明月-衛(wèi)莊）離散圖（512*512）計算全息圖再現(xiàn)圖實(shí)現(xiàn)程序%用imread()函數(shù)讀入物波函數(shù)抽樣數(shù)據(jù)并保存到矩陣A中,這樣原始物體就以離散數(shù)據(jù)形式保存在一個二值矩陣中(矩陣中每個元素非0即1)Image=imread('zcy6202.jpg');N=512;A=zeros(N,N);B=zeros(N,N);forI=1:1:256forJ=1:1:256ImageNum=double(Image(I,J,1));A(I,J)=ImageNum/255;B(I,J)=0;endendfigure;ims

2、how(A);%為降低全息圖的動態(tài)范圍，通過乘隨機(jī)相位因子的辦法來平滑物函數(shù)的傅里葉變換譜(起到光學(xué)全息中毛玻璃漫射光的作用)PI=3.14159;forI=1:1:NforJ=1:1:NR=rand(1,1);B(I,J)=A(I,J)*sin(R*2*PI);A(I,J)=A(I,J)*cos(R*2*PI);F(I,J)=A(I,J)+j*B(I,J);endend%對物函數(shù)做二維快速傅里葉變換（FFT）變換%F=fft2(F);Max=max(max(abs(F)));F=F/Max;A=real(F);B=imag(F);%定

3、義載波參數(shù)alpha=0.5;forI=1:1:NforJ=1:1:NXcos=(J-1)/127;A1(I,J)=cos(2*PI*alpha*Xcos);B1(I,J)=sin(2*PI*alpha*Xcos);endend%全息圖數(shù)據(jù)區(qū)forI=1:1:NforJ=1:1:NHolodata(I,J)=0.5+0.5*(A(I,J)*A1(I,J)+B(I,J)*B1(I,J));endend%繪制全息圖，制作灰階全息圖像，灰度等級256M=512;%定義全息圖的大小，Hologram=zeros(M,M);S=M/N;%定義每個

4、抽樣單元大小為S,SforI=1:1:NforJ=1:1:NXa=(J-1)*S+1;Xb=J*S;Ya=(I-1)*S+1;Yb=I*S;forIx=Xa:1:XbforIy=Ya:1:YbHologram(Iy,Ix)=Holodata(I,J);endendendendMax=max(max(Hologram));Hologram=Hologram/Max;figure;imshow(Hologram);%計算傅里葉變換全息的再現(xiàn)Object=fft2(Hologram);Object=fftshift(Object);Objec

5、t=abs(Object);Object=1000*Object/max(max(Object));figure;imshow(Object);