单模光纤二维强度分布的MATLAB仿真

类似于上一篇文章中的二维模场分布，这里的强度分布同样基于电场的表达式

我们知道，光强等于电场强度的平方。

对应代码如下：

clear
close all

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);

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

imagesc(x,y,I)
colorbar
xlabel('x')
xlabel('y')
title('I')

本次仿真中使用的V参数也和上次仿真相同。

结果：

如果需要，也可以同样绘制出等高线图与梯度箭头，这里不再重复讲解，只绘制二维彩色图形。

可以看到能量更多的集中在中央区域，由于我们设置的绘图范围较大，看上去只有中间的亮点。

我们可以观察不同的V参数对光纤中光场分布的影响，比如更改上面的程序中V值为0.8：

结果：

可以看到，将V参数变小之后，中央的亮点变大

这不是一件好事，因为这意味着更多的能量进入了包层，对于光的传输是不利的

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