单模光纤三维模场分布的MATLAB实现

在上一篇中我们已经计算出了模场分布的数据
三维图形和二维图形的区别在于绘图函数，我们直接换用surf函数即可：

clear
close all
clc

V = 2.4000;
U = 1.6453;
W = 1.7473;

Npoint = 501;
Rx = 5;
Ry = 5;
x = linspace(-Rx,Rx,Npoint);
y = linspace(-Ry,Ry,Npoint);
[X Y] = meshgrid(x,y);
R = sqrt(X.^2+Y.^2);

E1 = besselj(0,U*R);
E2 = besselj(0,U).*besselk(0,W.*R)./besselk(0,W);
I1 = E1.^2;
I2 = E2.^2;
I = I1;
E = E1;
pos = find(R >= 1);
E(pos) = E2(pos);
I(pos) = I2(pos);

surf(x,y,E)
shading interp
colorbar
xlabel('x')
xlabel('y')
zlabel('E')
title('LP_{01} mode, V = ',num2str(V))

figure
surf(x,y,I)
shading interp
colorbar
xlabel('x')
xlabel('y')
zlabel('I')
title('LP_{01} mode, V = ',num2str(V))

分别绘制出电场强度和光场强度的分布图，结果如下：

参考：《高等光学仿真——光波导、激光》