clear
clc
close?all
N?=1000;
i=sqrt（-1）;
lamda=632.8*10^（-6）;
k=2*pi/lamda;
f=1/lamda;
x0=0.5*lamda*f/pi;
y0=0.5*lamda*f/pi;
w0=5;%光束直徑
a=1/（w0^2）;

Xmin?=?-5;
Xmax?=?5;
Ymin?=?-5;
Ymax?=?5;
x?=?linspace（XminXmaxN）;
y?=?linspace（YminYmaxN）;

L=1.1*10^（-2）;
kXmin?=?0;
kXmax?=?（2*pi）/L;
kYmin?=?0;
kYmax?=?（2*pi）/L;
kx?=?linspace（kXminkXmaxN）;
ky?=?linspace（kYminkYmaxN）;

figure;
[xy]=meshgrid（xy）;
z=airy（x/x0）.*exp（a*x/x0）.*airy（y/y0）.*exp（a*y/y0）;
mesh（xyabs（z）.^2）
xlabel（‘Sx‘）;
ylabel（‘Sy‘）;
view（2）;

[kxky]=meshgrid（kxky）;
figure;
d=2;
FI=fftshift（fft2（z））;
h=exp（i*k*d*sqrt（1-lamda^2*（（kx.^2）+（ky.^2））））;
u=FI.*h;
u1=（ifft2（u））;
u2=abs（u1）.^2;
mesh（xyu2）
view（2）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????747??2013-04-24?23:45??Airy_pro.m

?????文件?????541319??2013-04-23?14:57??衍射理論對局域空心光束及無衍射光束重建的描述_劉嵐.pdf

-----------?---------??----------?-----??----

???????????????542534????????????????????4