matlab-利用MATLAB程序实现单缝衍射，完成拉盖尔高斯光束的光场、单缝孔以及LG光场通过单缝孔衍射后的衍射图-源码

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%LG光单缝衍射%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% clc; clear; close all; warning off; addpath(genpath(pwd)); l=1; na=0.6328*10^-6; z=57; w=7.05*10^-3; x0=-0.01995*10:0.00010*10:0.01995*10; size(x0) y0=-0.01995*10:0.00010*10:0.01995*10; x1=x0/(na*z); y1=y0/(na*z); [x0,y0]=meshgrid(x0,y0); A=zeros(400,400); A(:,198:203)=1; % A(:,193:194)=1; % A(:,203:204)=1; fai=1/i*log((x0+i*y0)./(x0.^2+y0.^2).^0.5); U0=(2/(pi*factorial(abs(l))))^0.5*(w^-1)*(2*(x0.^2+y0.^2)/w^2).^(abs(l)/2).*exp(-(x0.^2+y0.^2)/w^2).*exp(i*l*fai-i*abs(l)*pi/2); size(A) size(U0) U0=U0.*A; I0=abs(U0.*conj(U0)); figure(1) imshow(I0,[]); axis on; [x1,y1]=meshgrid(x1,y1); fai=1/i*log((x1+i*y1)./(x1.^2+y1.^2).^0.5); U00=(2/(pi*factorial(abs(l))))^0.5*(w^-1)*(2*(na*z)^2*(x1.^2+y1.^2)/w^2).^(abs(l)/2).*exp(-(na*z)^2*(x1.^2+y1.^2)/w^2).*exp(i*l*fai-i*abs(l)*pi/2).*A; U1=U00.*exp(i*na*z*pi*(x1.^2+y1.^2)); figure(2) imshow(U1,[]);axis on; f=fft2(U1); f=fftshift(f); x=-0.01995*10:0.00010*10:0.01995*10; y=-0.01995*10:0.00010*10:0.01995*10; [x,y]=meshgrid(x,y); U=exp(i*2*pi*z/na)*(-i*na*z)*exp(i*pi/(na*z)*(x.^2+y.^2)).*f; I=abs(U.*conj(U)); %figure(2) %imshow(U,[]);axis on; figure(3) %axis([150 250 100 300]); imshow(I,[]);axis on; figure(4) mesh(x,y,I); %colormap(gray); %变灰度图命令 shading interp; axis off; xlabel('x/m','Fontsize',14); ylabel('y/m','Fontsize',14); view([0,0,1]) [az,el]=view;