上传一个我写的数字全息软件，用Matlab6.5 GUIDE 做的
这个软件可以用来做数字全息图的预处理、全息再现、相位滤波、相位去包裹
可以应用于数字全息三维形貌测量、微小位移或形变测量
附件里dh文件夹中是图形界面的m文件和fig文件，prog文件夹是一些算法的m文件，被dh里的界面程序调用，image里面是几张全息实验的原始数据。
不排除程序有bug，如果有谁发现了通知我一下，谢了
希望从事数字全息研究的如果有兴趣可以扩充一下这个软件功能，改进完善一下。
2019-12-21 20:54:43
1.48MB
1
N=512;
A=zeros(N,N);
B=zeros(N,N);
for I=1:1:256
J=1:1:256
ImageNum=double(Image(I,J,1));
A(I,J)=ImageNum/255;
B(I,J)=0;
end
end
figure;
imshow(A);
pi=3.1415926;
for I=1:1:N
for J=1:1:N
R=rand(1,1);%生成一个元素在0,1之间均匀分布的随机矩阵R
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%限制振幅的动态范围,提高编码的精度
F=fft2(F);%作二维快速傅里叶变换FFT
Max=max(max(abs(F)));
F=F/Max;
A=real(F);
B=imag(F);
aIpha=0.5;%定义载波参数aIpha
for I=1:1:N
for J=1:1:N
Xcos=(J-1)/127;
A1(I,J)=cos(2*pi*aIpha*Xcos);
B1(I,J)=sin(2*pi*aIpha*Xcos);
end
end%全息图数据区
for I=1:1:N
for J=1:1:N
HoIodata(I,J)=0.5+0.5*(A(I,J)*A1(I,J)+B(I,J)*B1(I,J));
end
End
M=512;N=512;%定义全息图的大小
Hologram=zeros(M,M);
S=M/N;%定义每个抽样单元大小
for I=1:1:N
for J=1:1:N
Xa=(J-1)*S+1;
Xb=J*S;
Ya=(I-1)*S+1;
Yb=I*S;
for Ix=Xa:1:Xb
for Iy=Ya:1:Yb
HoIogram(Iy,Ix)=HoIodata(I,J);
end
end
end
end
Max=max(max(HoIogram));
HoIogram=HoIogram/Max;
figure;
imshow(HoIogram);
%以下是用matlab分别计算函数各抽样点的傅里叶变换谱的幅角与模,并对各点的模归一化
object=fft2(HoIogram);
object=fftshift(object);%用matlab中的移谱函数fftshift( )将频谱的低频成分移到中心,以避免再现时像分散在边缘
object=abs(object);
object=1000*object/max(max(object));
figure;
imshow(object);
2019-12-21 20:52:45
973B
1
伽博同轴全息记录与再现。通过衍射计算完成光学同轴全息记录与再现过程的仿真，可以观察再现像的特点。MATLAB代码
2019-12-21 19:55:36
76KB
1
离轴全息记录与再现。复习衍射计算以及光学全息记录与再现原理，学习通过衍射计算完成光学离轴全息记录与再现过程仿真方法，理解第二代光学全息共轭像分离的条件，同时与通州全息作比较，注意观察再现像的特点。MATLAB代码
2019-12-21 19:55:36
3.97MB
1