%建立airy光束ticclearlamda=0.532;%波长k=2*pi/lamda;A=1[x,y]=meshgrid(-20:40/200:20,20:-40/200:-20);theta=atan2(y,x)+pi*(abs(sign(y))-sign(y));%极坐标旋向r=sqrt(x.^2+y.^2);%极坐标径向z=0;U=A*airy((x+y)-(z/2).^2).*airy((x-y)-(z/2).^2).*exp(1i*z/2*((x+y)+(x-y)-z^3/3));I=U.*conj(U);imshow(I,[0,max(max(I))]),colormap(parula)toc
得到结果为
免责声明：本公众号旨在传递更多科研相关内容及分享，所有文章图片来源均注明出处，如涉及版权问题，请作者第一时间后台联系，我们将协调进行处理，对于投稿内容文责自负。
