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

这个程序用于仿真孤子的演变波形，这里Lcd 可以取不同的值-This Matlab scr ipt file solves the nonlinear Schrodinger equation

这两个程序用于仿真 图4和图5.这里的方程用新的归一化方程，不是孤子归一化方程。也就是论文中的方程3-study input peak power vursus output total power for different Lcd, P~P CURVES Simulate Fig.4 Using wangs equations