h 数字全息的菲涅耳再现算法和角谱算法

数字全息的菲涅耳再现算法和角谱算法，包括频谱滤波。

%菲涅耳算法计算同轴菲涅耳全息图的再现象
%lambda=632.8nm再现距离z，dxh=dyh=8.33um512*512像素
close all;
clear;
lamda=632.8*10^-9;         % 首先要计算机模拟再现光，使之与全息图相乘，得到重构波前在全息图平面的复振幅分布，再利用菲涅尔变换法重建
z=-110*10^-3;               %再现距离
k=2*pi/lamda;              %波矢量
A=exp（1i*k*z）/（1i*lamda*z）;%相位因子
dxh=8.33*10^-6;            %像素大小
dyh=8.33*10^-6;
xh=ones（5121）*（-256:255）*dxh;
yh=（ones（5121）*（-256:255）*dyh）‘;
dxi=lamda*z/（512*dxh）;    %再现像素大小
dyi=lamda*z/（512*dyh）;
xi=ones（5121）*（-256:255）*dxi;
yi=（ones（5121）*（-256:255）*dyi）‘;
cosx=0.0180;              %夹角的余弦值
cosy=0.0180;
I=double（rgb2gray（imread（‘512p1.bmp‘）））;
ref=exp（-1i*k*（xh*cosx+yh*cosy））;      %再现光
FI=fftshift（fft2（I））;
mesh（log（abs（FI）））;
%imwrite（mat2gray（log（abs（FI）））‘23.bmp‘）;
%figureimagesc（log（abs（FI）））colormap（gray）;
a=zeros（5

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1589  2012-05-05 16:18  h\fresnell.m

     文件        965  2012-05-09 11:22  h\jiaopu.m

     目录          0  2012-05-09 11:59  h

----------- ---------  ---------- -----  ----

                 2554                    3