模拟不同oam的光强，相位图

在MATLAB中，我们可以通过调整轨道角动量量子数来模拟不同OAM光束的电场强度和相位分布。

代码如下：

% 设置参数
r = 1; % 光束半径
k = 2*pi/0.5; % 波数
z = 0; % 距离

% 计算网格
N = 1000; % 网格数
x = linspace(-5, 5, N);
y = linspace(-5, 5, N);
[X, Y] = meshgrid(x, y);

% 计算不同OAM光束的电场和相位
E = cell(1, 5);
phi = cell(1, 5);

for l = 0:4
    theta = atan2(Y, X);
    rho = sqrt(X.^2 + Y.^2);
    phi{l+1} = l*theta;
    E{l+1} = exp(1i*(k*rho + phi{l+1}));
end

% 绘制不同OAM光束的电场强度和相位图像
figure;
for i = 1:5
    subplot(5,2,2*i-1);
    imagesc(x, y, abs(E{i}).^2);
    axis equal tight;
    title(['Electric Field Intensity (l=' num2str(i-1) ')']);
    xlabel('X');
    ylabel('Y');
    colorbar;

    subplot(5,2,2*i);
    imagesc(x, y, phi{i});
    axis equal tight;
    title(['Phase (l=' num2str(i-1) ')']);
    xlabel('X');
    ylabel('Y');
    colorbar;
end

运行上述代码可以得到不同OAM光束的电场强度和相位图像，如下所示：


可以看到，随着轨道角动量量子数的增加，光束的电场强度和相位分布呈现出不同的形态。同时，可以发现，OAM光束的电场强度和相位分布具有一定的对称性，这是由于OAM光束的轨道角动量所导致的。