高斯光束的聚焦仿真

clear all
w0=200e-9; %高斯光束的光腰半径
wavelength=532e-9; %高斯光束的光波长
ZR=piw0^2/wavelength; %高斯光束的瑞利长度
L=3ZR; %高斯光束的传输范围
N=200;
z=linspace(-L,L,N);
r=linspace(-4w0,4w0,N);
[pr,pz]=meshgrid(r,z);
wz=w0sqrt(1+(wavelengthpz/pi/w0^2).2);
I=w0^2./wz.2.exp(-2pr.^2./wz.2);

surf(pz,pr,I);
shading interp;
xlabel(‘Z(λ) /m’);
ylabel(‘r(λ) /m’);
zlabel(‘intensity /a.u.’);
axis([-3ZR 3ZR -Inf Inf -Inf Inf])
title(‘焦点附近三维光场分布’)
colorbar;
colormap hot; %colormap gray线性灰度图 colormap hot热色颜色图
view(2) %显示二维视图