MATLAB代码：锥形光纤纤芯周围功率分布轮廓

% 定义常量 n_core = 1.45; % 纤芯折射率 n_clad = 1; % 包层折射率 d_core = 400e-9; % 纤芯直径 lambda = 633e-9; % 入射光波长

% 定义空间范围和离散步长 r_max = d_core/2; % 最大半径 dr = r_max/100; % 离散步长 z_max = 10*d_core; % 最大长度 dz = z_max/100; % 离散步长

% 定义网格 r = 0:dr:r_max; % 径向坐标 z = 0:dz:z_max; % 纵向坐标 [R, Z] = meshgrid(r, z); % 网格

% 计算模场分布 V = (2pin_core*d_core/lambda) * sqrt((n_core^2-n_clad^2)R.^2 + n_clad^2Z.^2); % 归一化模场

% 绘制功率分布轮廓 P = V.^2.R; % 功率分布 figure; contourf(r1e9, z*1e9, P, 20); xlabel('Radius (nm)'); ylabel('Z (nm)'); title('Power Distribution of Conical Fiber');