Airy 光束的MATLAB仿真

目录

一维Airy光束

仿真代码

输出结果

二维Airy光束

仿真代码

输出结果

基于旋转矩阵的Airy光束阵列

​​​​​​​仿真代码

输出结果

写在最后

---

一维Airy光束

经截肢的有限能量一维Airy光束光强表达式如下：

式中各参数及相关表达式的含义：φ  表示电场包络，s=x/x0  为一维无量纲横向坐标，x0  为任意的横向尺度，ξ=z/k0x02  为归一化传输距离，k0=2πn/λ0  为波数。

仿真代码

clc,clear
%% 有限能量艾里光束传播轨迹
% 参数，单位是 mm
a = 0.05;
x = linspace(-0.5, 0.5, 1000);x0 = 8E-3; % 8um
z = linspace(0, 4, 1000); lambda = 850E-6;k = 2 * pi / lambda; 
e = z / (k * x0^2);sx = x / x0;[SX, E] = meshgrid(sx, e);
% 计算光场分布
U = airy(SX-((E/2).^2)+1i*a.*E)...
        .*exp(a.*SX-a.*E.^2/2 - ...
        1i*E.^3/12+1i*a^2.*E./2+1i*SX.*E./2);
I = abs(U).^2;
% 绘图
imagesc(x, z, I);
xlabel('x/mm'); ylabel('z/mm');
grid on;axis xy;
colorbar;colormap jet;
unit = 'um';
title(sprintf('衰减因子 a = %s  横向尺度 x0 = %s %s', ...
         num2str(a),num2str(x0*1E3), unit));

输出结果

二维Airy光束

二维Airy光束光强表达式如下：

​​​​​​​仿真代码

clc,clear
%% 二维艾里光束
% 参数
a = 0.05;z = 0;
x = linspace(-0.2, 0.2, 1000);x0 = 8E-3; 
y = linspace(-0.2, 0.2, 1000);y0 = 8E-3; 
lambda = 850E-6;k = 2*pi/lambda;e = z/(k*x0^2);
sx = x/x0; sy = y/y0; [SX,SY] = meshgrid(sx,sy);E = meshgrid(e);
% 计算光场分布
U = airy(SX-(E.^2/4)+1i*a*E) ...
        .*exp(a*SX-a*E.^2/2-1i*E.^3/12+1i*a^2*E/2+1i*SX.*E/2)...
        .*airy(SY-(E.^2/4)+1i*a*E)...
        .*exp(a*SY-a*E.^2/2-1i*E.^3/12+1i*a^2*E/2+1i*SY.*E/2);
I = abs(U).^2;
% 绘制图像
imagesc(x, y, I); 
colorbar; colormap jet;
xlabel('x/mm'); ylabel('y/mm');
grid on; axis xy;
unit = 'mm';
title(sprintf('传输距离 z = %s %s ', num2str(z),unit));

输出结果

基于旋转矩阵的Airy光束阵列

包含n个二维Airy光束的Airy光束阵列的光场表达式如下：

​​​​​​​仿真代码

clc,clear
%% 艾里光束阵列，基于旋转矩阵
% 参数,单位mm
a = 0.05; z = 0; c = 0.1; n = 4;
x = linspace(-0.3, 0.3, 1000); x0 = 8E-3;  
y = linspace(-0.3, 0.3, 1000); y0 = 8E-3; 
lambda = 850E-6; k = 2 * pi / lambda; 
e = z / (k * x0^2);E = meshgrid(e);
I_sum = zeros(size(x, 2), size(y, 2)); % 分配空间
for i = 1:n
    sx = x/x0; sy = y/y0;
    theta = 2*pi*(i-1)/n + pi;
    M = [cos(theta), -sin(theta); sin(theta), cos(theta)];
    N = [sx;sy];
    C = [c/x0;c/y0];
    S = M * N + C; 
    [SX,SY] = meshgrid(S(1, :),S(2, :));
    U = airy(SX-((E/2).^2)+1i*a*E)...
        .*exp(a*SX-a*E.^2/2-1i*E.^3/12+1i*a^2*E/2+1i*SX.*E/2)...
        .*airy(SY-((E/2).^2)+1i*a*E)...
        .*exp(a*SY-a*E.^2/2-1i*E.^3/12+1i*a^2*E/2+1i*SY.*E/2);
    I_sum = I_sum + abs(U).^2;    
end   
% 绘制图像
imagesc(x, y, I_sum); 
colorbar; colormap jet;
xlabel('x/mm'); ylabel('y/mm');
axis xy; grid on;
unit = 'mm';
title(sprintf('传输距离 z = %s %s 偏移量 c = %s %s', ...
    num2str(z),unit,num2str(c), unit));

输出结果

写在最后

本文为毕业论文的干货整理，欢迎大家交流学习，批评指正。

参考文献

1. 赵超越.艾里光束的新型数字化调控与测量[D].哈尔滨理工大学,2019.
2. 鲁强.基于LCOS的Airy光束阵列生成及其光强闪烁特性的研究[D].中国科学院长春光学精密机械与物理研究所,2017.