【光学】涡旋光产生【含Matlab源码 1927期】.zip

%vortex beam clear clc row=1080; g1=zeros(row); g2=zeros(row); w0=100; l=3; for m1=1:row for n1=1:row if (m1-row/2).^2+(n1-row/2).^2<row/2*row/2 g1(n1,m1)=sqrt((m1-row/2)^2+(n1-row/2)^2)^abs(l)*exp(-((m1-row/2)^2+(n1-row/2)^2)/w0^2)*exp(1i*l*atan2((n1-row/2),(m1-row/2))); end end end A0=max(max(abs(g1))); for m1=1:row for n1=1:row if (m1-row/2).^2+(n1-row/2).^2<row/2*row/2 g2(n1,m1)=A0.*exp(-((m1-row/2)^2+(n1-row/2)^2)/(row*0.2)^2).*exp(1i*pi*-2*((m1-row/2)^2+(n1-row/2)^2)/(row/6)^2); end end end g=g1+g2; intensity=abs(g1); phase=angle(g1); figure();imshow(intensity,[]);colormap hot; figure();imshow(phase,[]);colormap jet; intensity1=abs(g); figure();imshow(intensity1,[]);colormap hot;