matlab 全息图,全息图matlab程序.doc

全息图matlab程序

%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;

A=real(F);

B=imag(F);

%对物函数做二维FFT变换%

%定义载波参数

alpha=0.5;