厄米-高斯光束.rar_matlab高斯光束_光束_厄米-高斯光束_厄米光束MATLAB_高斯光束

%% main.m lambda = 500; k = 2*pi/lambda; z0 = 1; A0 = 1; W0 = sqrt(lambda*z0/pi); W = @(z) W0*sqrt(1+(z/z0)^2); R = @(z) z*(1+(z/z0)^2); Zeta = @(z) atan(z/z0); A = [ 1 0 0 0; 1 1 0 0; 0 0 0 0; 0 0 .2 0]; res = 800; z = 1e-9; x = linspace(-2*W(z),2*W(z),res); y = linspace(-2*W(z),2*W(z),res); [X Y] = meshgrid(x,y); X = X(:); Y = Y(:); U = zeros(length(X),1); Utemp = zeros(length(X),1); Utemp2 = zeros(length(X),1); lpf = exp(-1i*k*z - 1i*k*(X.^2 + Y.^2)/(2*R(z))); %lateral phase factor u = sqrt(2)*X/W(z); v = sqrt(2)*Y/W(z); for l = 1:size(A,1) if sum(A(l,:) ~= 0) ~= 0 Utemp2 = (W0/W(z))*polyval(HermitePoly(l-1),u).*exp(-u.^2/2); else continue; end for m = 1:size(A,2) if A(l,m) ~= 0 Utemp = Utemp2.*(polyval(HermitePoly(m-1),v)).*exp(-v.^2/2); Utemp = A(l,m)*Utemp.*lpf*exp(1i*(l+m+1)*Zeta(z)); U = U + Utemp; end end end %% DRAW PLOTS figure; U = reshape(U,res,res); imagesc(x,y,abs(U).^2); axis square; set(1,'color','w'); title('厄米-高斯光束'); xlabel('x (nm)'); ylabel('y (nm)');