本帖最后由 欧攀 于 2013-1-29 10:01 编辑
53页到56页,主要阐述物理背景。
我主要的编程工作是根据公式3.42,以及公式3.45。
clear
clc
f=25; %焦距
a=0.3;
b=0.3; %微透镜尺寸
lamda=0.6238e-3; %波长
k=2*pi/lamda;
x=-1.2:0.01:1.2; %对微透镜表面或者说入射波前取离散点
y=-1.2:0.01:1.2;
n0=1.75; %透镜折射率
Delta0=2; %透镜厚度
n=length(x);
u0=ones(n);
[x1 y1]=meshgrid(x,y);
fx=x1./lamda./f; %根据公式3.46
fy=y1./lamda./f;
x_org=-1.05:0.3:1.05; %微透镜中心坐标
y_org=-1.05:0.3:1.05;
nx=length(x_org);
ny=length(y_org);
Ft2=zeros(n,n);
% uf=1./j./lamda./f.*exp(j.*k./2./f.*(x1.*x1+y1.*y1));
Ft1=fft2(u0.*exp(j.*k.*n0.*Delta0)); %公式3.4