贝塞尔光束光强模拟程序

贝塞尔光束的横向光强分布表现为一个中心光斑和一系列同心圆环。在物理上可以实现的贝塞尔光束，其无衍射传播范围是有限的。在贝塞尔光束的无衍射范围内，贝塞尔光束保持横向光强分布，即使在遇到不透明障碍物后也可以恢复到原来的横向光强分布。贝塞尔光束独特的光强分布和传播性质使其得到了广泛应用，例如光学成像，微细加工，光学互联和校直，粒子操控，微缩平板印刷，非线性光学等。

% This program calculates the Azimuthal Bessel-Gauss beam based on the
% formula of J. Opt. Soc. Am. A Vol.15 No.12 1998

tic;
% parameters for Azimuthal Bessel-Gauss （ABG） beam.
% all length are in the unit of mm.
mm=1.0e-3;
lambda=632.8e-6*mm;              % wavelength of incident BG beam
k=2*pi/lambda;                   % wave number of incident BG beam
omega0=1.0*mm;                   % waist of incident BG beam
beta=0.0005*k;                   % transverse wave number of ABG beam
L=k*omega0^2/2;                  % Rayleigh range

x=linspace（-2*mm2*mm401）;      
y=x;
[XY]=meshgrid（xy）;
X=X‘;Y=Y‘;
r=sqrt（X.^2+Y.^2）;
theta=atan2（YX）;
z=lambda;

omega=omega0*sqrt（1+