全息图matlab程序.docVIP

%Modified off-axis reference beam CGH of Burch% %读入原始图象 PI=3.14159; Image=imread(lena256.jpg); N=512; A=zeros(N,N); B=zeros(N,N); for I=1:1:256 for J=1:1:256 ImageNum=double(Image(I,J,1)); A(I,J)=ImageNum/255; B(I,J)=0; end;end; figure; imshow(A); %为降低全息图的动态范围，乘一随机位相因子% for I=1:1:N for J=1:1:N R=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); end;end; %为降低全息图的动态范围，乘一随机位相因子% %对物函数做二维FFT变换% F=fft2(F); Max=max(max(abs(F))); F=F/Max;