matlab傅里叶光学仿真,关于微透镜阵列的傅里叶光学分析

本帖最后由 欧攀 于 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.42最后