引言

光纤耦合器是一种在多根光纤之间、或在光源与光纤之间实现光功率定向传输的无源器件。该器件可广泛应用于光纤通信、光纤检测及光纤传感等领域。一般情况下，光纤耦合器中输入的光能量大小可以决定耦合器的工作状态。当输入光能量较低时，光纤耦合器处于线性工作状态；当输入光能量高到一定程度后，光纤中将产生非线性效应，全光开关是非线性光纤耦合器的一个非常重要的应用实例，同时也是目前的热门研究课题。

光脉冲在定向耦合器中的传输可由一组耦合模方程来描述，方程的个数可以由定向耦合器中波导的个数来确定。方程组在一般情况下没有解析解，但可以通过一定方法得到其数值解。虽然在1990年就有人首次使用变分法得到了方程的近似解，但是，当波导数较多时，其计算十分繁琐。事实上，耦合模方程组也可以通过傅里叶级数展开法将耦合模方程转换到频域中，然后再采用龙格库塔法进行求解。

另外，广泛应用到解非线性色散介质脉冲传输问题的另一种方法是分步傅里叶法，这种方法由于采用了快速傅里叶变换(FFT)算法而具有较快的计算速度。而采用迭代算法对该算法作进一步改进，则可得到更高的计算精度，改进后的算法称为对称分步傅里叶法。本文给出了用对称分步傅里叶法求解耦合模方程组的方法，并通过该方法采用Matlab对光脉冲在双芯和三芯非线性耦合器中的传输进行了仿真。

1传输分析

描述光脉冲在N芯耦合器中传输的耦合模方程组如下：

为二、三阶色散和损耗，γn为非线性系数，knj是纤芯n和j之间的线性耦合系数，ηnj为模间色散。

一般来说，沿光纤长度方向的色散和非线性是同时作用的。分步傅里叶方法通过假定在传输过程中，光场每通过一小段距离dz，色散和非线性效应可分别作用，来得到近似结果。当步长dz足够小时，这种分析有足够的精度。对称分步傅里叶法的原理图如图1所示。

根据对称分步傅里叶法原理，可将方程(1)分为非线性和色散部分分别进行求解。