空间光-单模光纤耦合效率计算

风扇考察

已于 2024-07-14 22:20:07 修改

阅读量1.4k

点赞数 39

分类专栏： zemax 光纤文章标签：改行学it 笔记

于 2024-07-14 22:17:07 首次发布

本文链接：https://blog.csdn.net/sinat_35238851/article/details/140404680

版权

zemax 同时被 2 个专栏收录

2 篇文章 0 订阅

订阅专栏

光纤

1 篇文章 0 订阅

订阅专栏

一、耦合效率的定义

光纤耦合效率被定义为耦合进光纤的光功率与耦合系统接收平面上接收的光功率之比；

可以用下面的积分公式来描述：

$\eta=\frac{\left|\iint E_{fiber, b}^* E_{lens, b} \mathrm{~d} s\right|^2}{\iint\left|E_{fiber, b}\right|^2 \mathrm{~d} s \cdot \iint\left|E_{lens, b}\right|^2 \mathrm{~d} s}$

测试一下，贴图片和csdn自带Latex公式编辑器的清晰程度

文章引用：

【Massimo Lazzaroni, Fabio E. Zocchi,Optical coupling from plane wave to step-index single-mode fiber,Optics Communications】

分子：两种不同电场的重叠积分；(代表耦合进光纤的能量,模场匹配)

分母：两种电场强度的模的乘积，自己积自己；(代表空间光经过光学系统耦合到光纤端面的全部能量)

$E_{fiber,b}$ 指的是单模光纤端面的电场分布，一般来说单模光纤的强度分布为零阶贝塞尔强度轮廓，或者强度接近于高斯分布。

$E_{fiber,b}=\exp \left[-\left(r / r_0\right)^2\right]$

r为光纤端面任意一点到中心的距离，r0为单模光纤模场半径，康宁SMF-28e+光纤（1550nm），r0一般为5.25um，模场直径10.5um。

$E_{lens,b}$ 指的是光学系统在光纤端面的光场分布，无穷远的空间光传播到光学镜头第一个面，或者说入瞳的时候，可以看做是平行光，平面波。依据夫琅禾费-圆孔-衍射成像的理论，在像面上形成一个艾里斑。

圆孔远场夫琅禾费公式：

整个公式和二元的傅里叶变换表达式接近，那么可以用贝塞尔公式来求解积分，于是科学家们又得到了如下的表达式：

二、直白化的理解

公式确实比较难以理解……，其实可以换一个角度来考虑这个问题，如果要达到耦合效率最大，或者说在什么情况下两种不同的光场是最匹配的，那最直接的形式就是两个光场完全相同，如果做不到两个光场相同，那就让他们尽可能的相似。

简而言之，光纤模场-光学系统艾里斑，尽可能相似。继续从工程应用的角度出发，我们一般用束腰直径（1/e2），模场直径、发散角等参数来描述高斯光束。从光纤的手册中可以得到，1550nm波段使用的单模光纤模场直径MFD=10.5um。

2.1 光纤MFD-1550nm-10.5um

2.2 艾里斑的计算公式

从衍射的公式来看，艾里斑直径为 $D_{Air}=2R=2.44*\tfrac{\lambda }{(\tfrac{D}{F})}$ ，其中D/F是我们比较关心的光学系统孔径，即F数。下图中的艾里斑表达是半径的表达，所以常数系数是1.22。

注意，这里的 $D_{Air}$ 是第一暗环的中心直径，是艾里斑光强第一次等于0的中心直径，和我们高斯光束1/e2，13.5%的强度轮廓是不同的，需要进行光斑直径的转换。

2.3 MATLAB代码的艾里斑计算

很快我们可以得到 $D_{Air}=3.8317,D_{Air{0.135}}=1.5291,D_{Air}/D_{Air{0.135}}=2.5059$ ，这里表示光强的13.%，用电场积分的面积的13.5%来计算光斑轮廓。

计算得到Dair=1.47*MFD时，两种光斑匹配最佳（这里需要补充文献再看看实际的设计）

%%
Z = -10:0.1:10;
I1 = 4*((besselj(1,Z))./Z).^2; %艾里斑表达式

%这一部分用来计算艾里斑中心的值，求极限得到
syms X
eqn2 = ((besselj(1,X))/X)^2;
I1(101) = 4*limit(eqn2,X,0);

plot(Z,I1)
hold on

I2 = exp(-(Z/2.588).^2).^2;
plot(Z,I2)
% 2.588 = S2/1.48 经验公式

%% 这里用来数值计算，13.5%的艾里斑光强轮廓对应的光束直径Dair（13.5%）
syms Z
eqn1 = ((besselj(1,Z))/Z)^2 == 0.135
S1 = vpasolve(eqn1,Z,1)

% 这里用来数值计算，0%的艾里斑光强轮廓对应的光束直径Dair（光强为零）
eqn2 = ((besselj(1,Z))/Z)^2 == 0
S2 = vpasolve(eqn2,Z,1)

三、关于耦合的讨论

将艾里斑和模场直径的式子做一个变换，可以得到光学系统孔径D，和系统焦距F，波长λ，光纤MFD之间的关系，满足如下的经验公式时，耦合效率理论上可以达到最大值。例如，考虑在车载激光雷达系统中，受限于系统高度要求，D<35mm，波长1550nm或者905nm是固定的，想要或者一个焦距较小的F，就需要减小MFD的数值，例如选用6.5um模场MFD的光纤，能够减小系统总长，短焦距也可以设计较大的视场的镜头，相对而言比较简单。

当然，总体来说需要按照系统参数要求，合理分配其中的光学参数。

$D_{Air}=2.44*\tfrac{\lambda }{(\tfrac{D}{F})}=1.47*MFD$