chenweiqiang11.m

来自「数字全息再现与优化算法研究」· M 代码 · 共 33 行

close all;
clear;
w=0.6328,z=84000,L=5120,M=512,g=0.6*pi/180;
N=L/(w*z); 
k=2*pi/w;
[x,y]=meshgrid(linspace( -1,1,M-1).*L/2);
m=length(x);%U = input ('U = ');
U=imread('true.bmp');
Un= size (U);
U=im2bw(U);
Unn=size(U);
figure(1);
imagesc (U) ,colormap (gray) %图3a
% part 2 衍射物波的复振幅分布
[u,v]=meshgrid(linspace(-1,1,M-1).*N/2);
Ui=fftshift(fft2(U));
%平行光垂直照射衍射摸板
%part 3 记录数值化的全息图
R =exp(i*k.*(x.*sin(g))); % R 参考光
Ud=Ui.*R; 
Ud =abs(Ud).^2;
figure (2)
%Ud=imread('hol2.bmp');
imagesc(Ud,[min(Ud(:)),0.001*max(Ud(:))]) ,
colormap (gray) %图3b
%part 4 再现像
Uc = R ; % 用与参考光波方向相同的光波作为再现光
U =Ud.*Uc;
U =fftshift (ifft2 (U) );
U =abs (U) ;
figure (3)
imagesc(U,[min(U(:)),0.001*max(U(:))]),
colormap (gray) %图3c ,图3d