【光学】涡旋光与球面波的干涉【含Matlab源码 597期】.zip

%vortex beam clear clc row=1080; %像素点设置成1080pix g1=zeros(row); %生成1080阶零矩阵，初始化g1 g2=zeros(row); %生成1080阶零矩阵，初始化g1 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;