高斯光束复振幅分布MATLAB

%%高斯光束复振幅%%
​clc;clear;
cm=0.01;um=1e-6;mm=0.001;
lambda = 0.632*um;              %波长
x=-100*cm:0.1*cm:100*cm;
y=x;
z=0*cm;                         %截面在高斯光束中z轴位置
[X,Y]=meshgrid(x,-y);
r=sqrt(X.^2+Y.^2);         
k=2*pi/lambda;             
c=1;
%w0=sqrt(lambda*f/pi); 
w0=20*cm;                       %基模腰斑半径
f=pi*w0^2/lambda; 
R=z+f^2/z;
w=w0*sqrt(1+(z/f)^2); 
C_Amp=(c./w).*exp(-r.^2./w.^2).*exp(-1i.*(k*(z+r.^2/(2*R))-atan(z/f)))   ;
figure;imagesc(x,y,mat2gray(abs(C_Amp)));set(gca,'YDir','normal');