这个程序只合适计算k=1 (也就是Lcd=1）,对于其他k值，一些参数要做适当变化。这也是soliton format 不宜用于非线性耦合开关的一个原因。它使问题复杂化。这个程序是计算孤子格式化下的非线性薛定格方程组，程序中的方程与我们之前的博文提供的程序的方程不同。我们当初做这个仿真的目的，一方面要重复 S. Trillo, S. Wabnitz, E. M. Wright and G. I. Stegman, 在他们论文“Soliton switching in fiber nonlinea