airy_airy_Airy光束_

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 = -3; Xmax = 3; Ymin = -3; Ymax = 3; x = linspace(Xmin,Xmax,N); y = linspace(Ymin,Ymax,N); L=1.1*10^(-2); kXmin = 0; kXmax = (2*pi)/L; kYmin = 0; kYmax = (2*pi)/L; kx = linspace(kXmin,kXmax,N); ky = linspace(kYmin,kYmax,N); figure; [x,y]=meshgrid(x,y); z=airy(x/x0).*exp(a*x/x0).*airy(y/y0).*exp(a*y/y0); mesh(x,y,abs(z).^2) colormap pink shading interp xlabel('Sx'); ylabel('Sy'); view(2);