非线性薛定谔方程(NLSE)模拟光在非线性光纤和平面波导中的传播。NLSE具有多孤子解,只能通过离散谱来定义,这些谱由特征值和范数常数组成,不包含连续谱。经典的Darboux变换是生成多孤子的常用方法[1]。根据[2]中的算法2进行了一些修改,以改善数值条件,并根据[3]中的排序结果优化特征值的添加顺序,从而显著减少了舍入误差。