% 参数设置
wavelength = 1.55e-6; % 波长(单位:米)
fiber_length = 10e-3; % 光纤长度(单位:米)
mode_order = 1; % OAM的轴模式阶数
beam_radius = 1e-3; % 光束半径(单位:米)
% 光纤参数
fiber_attenuation = 0.2; % 光纤衰减系数(单位:1/米)
% 计算光纤中的OAM传输
z = linspace(0, fiber_length, 1000); % 在光纤长度上进行均匀分布的采样点
omega = 2*pi*wavelength; % 光的角频率
% 计算光束的复振幅分布
u = sqrt(2*fiber_attenuation*z); % 光纤中的衰减系数
% 计算相位分布
phase = mode_order * atan2(z, beam_radius^2);
% 计算复振幅
E = u .* exp(1i * phase);
% 绘制光束复振幅分布
figure;
plot(z, abs(E));
xlabel('光纤长度(米)');
ylabel('光束复振幅');
title('光纤中OAM的仿真结果');