Matlab:产生艾里光束

据说艾里光束是一类不沿着直线传播的光束，这里用matlab简单仿真了一个光束的产生。当然，还没有开始研究其中光束的传输特性。

由于艾里光束的理论二维光强分布的傅里叶变换是一个高斯振幅项×立方相位项，那么反过来，传统产生艾里光束的方法便是利用高斯光束与立方相位作用，并用焦距为f的透镜做傅里叶变换，那么在焦平面上便可以得到理想的艾里光斑。而这种方式所构建的实验系统必须要将透镜带来的2f距离考虑在内，这就会带来实验系统冗长的问题。而在器件微型化和小型化的当下，2f的距离简直是不能忍的，利用直径不到1mm的超构表面承载立方相位与透镜相位，可以直接解决问题。这也是后面两种仿真的思路，一种是不考虑光场的振幅分布，将立方相位叠加透镜相位，经过一定距离的衍射之后得到准艾里光束。而另一种则是将高斯振幅分布与立方相位点乘编码，再叠加上透镜相位，同样衍射之后得到艾里光束。

话不多说，上代码：

%% ***************************************
% author: Wangjing  
% 2024/8/14

clear,clc,close all;

% parameters setting**************
lamda = 1550e-9;
k = 2*pi/lamda;
xlens = 0.8e-3;
ylens = 0.8e-3;
w = 1.8e-4;        % waist radius
Nx = 1024;
Ny = 1024;
z = 0;

% coordinates generation**********
xpos = linspace(-xlens/2,xlens/2,Nx);
ypos = linspace(-ylens/2,ylens/2,Ny);
[x,y] = meshgrid(xpos,ypos);

% Gaussian_beams generation*******
E_Gauss = f_Gaussian_beams(w,z,lamda,xlens,ylens,Nx,Ny);
figure;
subplot(1,2,1),imagesc(abs(E_Gauss).^2),title('gauss_intensity'),colormap(subplot(1,2,1),"hot"),axis square,axis off;
subplot(1,2,2),imagesc(angle(E_Gauss)),title('gauss_phase'),colormap(subplot(1,2,2),"gray"),axis square,axis off;

% Cubic phase & lens phase*********
s = 5e-7;
A = exp(1i*k*(x.^3+y.^3)/3/s);
figure;imagesc(angle(A)),colormap gray,title('Cubic phase'),axis square,axis off;
f = 5e-4;
Fresnel = exp(1i*k*(x.^2+y.^2)/2/f);
figure;imagesc(angle(Fresnel)),colormap gray,title('lens phase'),axis square,axis off;

% 传统透镜傅里叶变换产生艾里光束*********
E1 = E_Gauss.*A;
E2 = fftshift(fft2(E1));      % 光斑较小
figure;
subplot(1,2,1),imagesc(abs(E2).^2),title('Airy_beams'),colormap(subplot(1,2,1),"hot"),axis square,axis off;
subplot(1,2,2),imagesc(angle(E2)),title('bemas_phase'),colormap(subplot(1,2,2),"gray"),axis square,axis off;

% 平面波入射到cubic及lens复合相位产生准艾里光束********
E_plane = ones(Nx,Ny);
figure;
subplot(1,2,1),imagesc(abs(E_plane).^2),title('平面波'),colormap(subplot(1,2,1),"hot"),axis square,axis off;
subplot(1,2,2),imagesc(angle(E_plane)),title('phase'),colormap(subplot(1,2,2),"gray"),axis square,axis off;

E_qusi = E_plane.*(exp(1i*angle(A+Fresnel)));
E_qusi_Airy = f_angular_spectrum(lamda,E_qusi,xlens/Nx,2*f,'forward');
figure;
subplot(1,2,1),imagesc(abs(E_qusi_Airy).^2),title('qusi_Airy_beams'),colormap(subplot(1,2,1),"hot"),axis square,axis off;
subplot(1,2,2),imagesc(angle(E_qusi_Airy)),title('qusi_phase'),colormap(subplot(1,2,2),"gray"),axis square,axis off;

% 平面波入射到振幅编码的cubic_phase并结合透镜相位产生艾里光束********
m = 1e-7;
Amp = -(x.^2 + y.^2)./m;
figure;imagesc(Amp),title('Amplitude'),colormap gray,axis square,axis off;
phase = angle(Amp.*A + Fresnel);
E_middle = E_plane.*exp(1i*phase);
E_Airy = f_angular_spectrum(lamda,E_middle,xlens/Nx,f,'forward');
figure;
subplot(1,2,1),imagesc(abs(E_Airy).^2),title('Airy_beams'),colormap(subplot(1,2,1),"hot"),axis square,axis off;
subplot(1,2,2),imagesc(angle(E_Airy)),title('phase'),colormap(subplot(1,2,2),"gray"),axis square,axis off;

下面是相应的仿真结果

传统傅里叶变换法的仿真光斑：

这是放大的焦面光斑，不放大看其实只有一个小点，实验上也应该是这个样子。

这是对应直接将立方相位与透镜相位叠起来的仿真结果，即准艾里光束（言外之意就是没那么标准的艾里光束）：

确实这个光斑效果都差很多，没有上面fft那么光滑，而且还有很多杂斑。

最后一种是对应振幅编码的艾里光束仿真结果：

emm,这个纯粹手误，出现了一个很明显的十字衍射，有可能是衍射距离没有设置好，或者是传输函数有点疏忽啥的，意思就那个意思。意思意思就可以了。

至于为什么要研究艾里光束，以及艾里光束有什么用，其实我也不知道，纯粹是写着玩的，那什么，如果文章或者代码有错误有疏忽，欢迎补充和批评指正。
祝大家工作科研顺利，暑假愉快~~~~