随着耦合强度的变化，给定的代码可以找到耦合振荡器系统的最大李雅普诺夫指数 (LLE)。 全局变量“be”和“gm”是表示两个耦合振荡器的耦合微分方程中的系统特定参数。 'gm' 是变化的耦合强度。 每个轨迹的初始条件包括两个振荡器的初始坐标和速度。 变量“ t0”和“ tf”代表每个时间序列的起点和终点，而“ N”是每个时间序列的步数。 在第 32 行和第 33 行中，调用函数 ode_RK4_Yang() 以通过 Runge-Kutta (RK4) 方法对两组初始条件中的每组的耦合微分方程进行数值求解。 其余的描述以注释的形式给出了代码本身。

matlab色散补偿程序，光谱仪采集到的数据进行一系列变换，画出耦合曲线的图