光波导仿真：传输特性仿真_（1）.光波导基础理论.docxVIP

PAGE1

PAGE1

光波导基础理论

在上一节中，我们介绍了光波导的基本概念和分类。本节将深入探讨光波导的基础理论，包括光波导的传输模式、波导中的光波传播方程、以及传输特性的基本分析方法。通过本节的学习，您将能够理解光波在波导中的传播机制，并为后续的仿真工作打下坚实的理论基础。

光波导的传输模式

光波导的传输模式是指光波在波导中传播时所呈现的不同状态。传输模式可以分为导模（GuidedModes）、辐射模（RadiationModes）和表面模（SurfaceModes）。其中，导模是最常见的模式，也是光波导传输的主要模式。

导模

导模是指光波在波导中传播时，其场分布主要集中在波导核心区域内的模式。导模的存在条件可以通过求解波导中的麦克斯韦方程组来确定。对于一个简单的单模光纤，导模的场分布可以用以下公式表示：

E

其中，Ezx,y是沿z方向的电场分布，E0

辐射模

辐射模是指光波在波导中传播时，其场分布主要集中在波导外部的模式。这些模态通常不被波导所限制，而是向周围空间辐射。辐射模的存在条件可以通过求解波导中的麦克斯韦方程组并考虑边界条件来确定。

表面模

表面模是指光波在波导表面附近传播的模式。这些模态主要存在于波导与周围介质的交界处，通常与表面粗糙度和不均匀性有关。表面模的存在条件可以通过求解波导表面的边界条件来确定。

波导中的光波传播方程

光波在波导中的传播可以通过求解麦克斯韦方程组来描述。对于平面波导，麦克斯韦方程组可以简化为波动方程。波动方程的一般形式为：

?

其中，E是电场，k是波数。对于特定的波导结构，波动方程需要结合边界条件来求解。

平面波导的波动方程

平面波导是常见的波导结构之一，其核心和包层可以视为无限长的平面。平面波导的波动方程可以通过以下步骤求解：

假设电场分布：假设电场在z方向上是均匀的，即Ex

代入波动方程：将假设的电场分布代入波动方程，得到：

?

求解特征方程：求解上述方程，得到特征值k2?β2

例子：单模光纤的波动方程

考虑一个单模光纤，其核心和包层的折射率分别为n1和n2。假设光波的波长为λ，波数为

?

?

Python代码示例：求解平面波导的特征值和特征函数

以下是使用Python和SciPy库求解平面波导特征值和特征函数的代码示例：

importnumpyasnp

fromscipy.sparseimportdiags

fromscipy.sparse.linalgimporteigs

importmatplotlib.pyplotasplt

#定义参数

n1=1.5#核心折射率

n2=1.4#包层折射率

lambda_=1.55e-6#光波波长(米)

k=2*np.pi/lambda_#波数

a=5e-6#核心半径(米)

L=20e-6#计算区域长度(米)

N=1000#网格点数

#定义网格

x=np.linspace(-L,L,N)

dx=x[1]-x[0]

#定义折射率分布

n=np.where(np.abs(x)a,n1,n2)

#构建离散化的波动方程矩阵

D2=diags([1,-2,1],[-1,0,1],shape=(N,N))/dx**2

K=k**2*n**2

H=D2+diags(K,0)

#求解特征值和特征函数

eigenvalues,eigenvectors=eigs(H,k=10,which=LM)

#绘制前几个特征函数

foriinrange(5):

plt.plot(x,np.abs(eigenvectors[:,i])**2,label=fEigenmode{i+1})

plt.xlabel(Position(x)[m])

plt.ylabel(Intensity[a.u.])

plt.legend()

plt.title(IntensityDistributionofEigenmodesinaPlanarWaveguide)

plt.show()

代码解释

定义参数：设置核心和包层的折射率、光波波长、核心半径和计算区域长度。

定义网格：生成计算区域内的网格点。

定义折射率分布：根据网格点的位置，设置核心和包层的折射率。

构建离散化的波动方程矩阵：使用中心差分法构建二阶导数矩阵，并结合折射率分布构造波动方程矩阵。

求解特征值和特征函数：使用eigs函数求解矩阵的前10个特征值和特征函数。

绘制特征函数：绘制前5个特征函数的强度分布。

传输特性的基本分析方法