这个程序只合适计算k=1 (也就是Lcd=1）,对于其他k值，一些参数要做适当变化。这也是soliton format 不宜用于非线性耦合开关的一个原因。它使问题复杂化。这个程序是计算孤子格式化下的非线性薛定格方程组，程序中的方程与我们之前的博文提供的程序的方程不同。我们当初做这个仿真的目的，一方面要重复 S. Trillo, S. Wabnitz, E. M. Wright and G. I. Stegman, 在他们论文“Soliton switching in fiber nonlinear directional couplers.” Opt. Lett. Vol.13, No.8, pp672-674(1988) 的仿真结果，以确保我们程序是正确的。另一方面要以己之矛攻己之盾。当时，孤子开关用孤子归一化的非线性薛定格方程描述已经广泛被接受和使用。很多人因为他们孤子开关的工作，也获得了重要的学术地位，要想发一篇跟他们之前的工作不同的观点的文章是很困难的。那些审稿的人，不是他们本人，就是学生或粉丝。图1中的曲线跟我们其它文中的曲线有很大不同，这是因为在这里输入脉冲是bsech(bt),因而输入脉冲的脉宽跟脉冲幅度一起变化，而其他地方我们只考虑脉冲幅度的变化，因为脉宽的变化已经包含在参量Lcd 中了。-This Matlab file solves the nonlinear Schrodinger equation in soliton format