代码：
clcclear allclose all%% 振幅光学相位元件产生完美涡旋光场N = 300;row = linspace(-1.5e-6,1.5e-6,N); col = linspace(-1.5e-6,1.5e-6,N);[x,y] = meshgrid(row,col);[phi,rho] = cart2pol(x,y);lambda = 632e-9; %波长632nmk = 2*pi/lambda; %波数f = 0.5; %透镜焦距0.5mR = 0.5; %圆孔光阑半径0.3mX = k*rho/f;alpha = 20;%径向波矢 - 常数figure;for m = 0 : 8subplot(3,3,m+1)E1 = (-1i)^(m+1)*(k*R/f)*exp(1i*m*phi).*...((alpha*besselj(m+1,alpha*R)*besselj(m,X*R)-X.*besselj(m,alpha*R).*besselj(m+1,X*R))./(alpha^2-X.^2)); I1 = E1.*conj(E1); I1 = I1/max(max(I1));%二维h1 = pcolor(x,y,I1);colorbar;set(h1,'edgecolor','none','facecolor','interp');title(['m = ',num2str(m)]);%colormap(gray); %输出灰度图像axis square;%三维%mesh(x,y,I1) %三维%set(gca,'fontname','times new roman','fontsize',16);%title(['m = ',num2str(m)],'fontname','华文中宋','fontsize',16);%xlabel('x/m','fontname','times new roman','fontsize',16);%ylabel('y/m','fontname','times new roman','fontsize',16);%zlabel('归一化强度','fontname','华文中宋','fontsize',16);endsuptitle('振幅光学相位元件产生完美涡旋光场') %为图一添加总标题%% 利用锥透镜产生完美涡旋光场row = linspace(-0.8,0.8,N); col = linspace(-0.8,0.8,N);[x,y] = meshgrid(row,col);[phi,rho] = cart2pol(x,y);lambda = 632e-9; %波长632nmk = 2*pi/lambda; %波数f = 0.5; %透镜焦距0.5mR = 0.5; %圆孔光阑半径0.3malpha = 20;%径向波矢 - 常数w0 = 0.3; %透镜焦平面上的束腰w_g = 0.5;%高斯项的束腰figure;for m = 0 : 8subplot(3,3,m+1)E2 = (1i)^(m+1)*(w_g/w0)*exp(1i*m*phi).*exp(-(rho.^2+R^2)/w0^2).*besseli(m,2*R*rho/w0^2);I2 = E2.*conj(E2); I2 = I2/max(max(I2));%二维h2 = pcolor(x,y,I2);colorbar;set(h2,'edgecolor','none','facecolor','interp');title(['m = ',num2str(m)]);%colormap(gray); %输出灰度图像axis square;%三维%mesh(x,y,I2) %三维%set(gca,'fontname','times new roman','fontsize',16);%title(['m = ',num2str(m)],'fontname','华文中宋','fontsize',16);%xlabel('x/m','fontname','times new roman','fontsize',16);%ylabel('y/m','fontname','times new roman','fontsize',16);%zlabel('归一化强度','fontname','华文中宋','fontsize',16);endsuptitle('利用锥透镜产生完美涡旋光场') %为图二添加总标题
结果：
参考文献：
[1]王亚军, 李新忠, 李贺贺, et al. 完美涡旋光场的研究进展[J]. 激光与光电子学进展, , 054(009):67-74.
如果觉得《Matlab：完美涡旋光束》对你有帮助，请点赞、收藏，并留下你的观点哦！
