; 参数：
ls：= 1060 nm {中心波长}
T_range：= 25 ps {时间范围的宽度}
N_t：= 2 ^ 10 {网格点数}
dt：= T_range / N_t {时间分辨率}
E0：= 10 nJ {初始脉冲能量}
tau0：= 1 ps {初始脉冲持续时间}
chirp0：= 0 GHz / ps {初始chi}
 
calc {计算以下复合表达式}
  开始
  set_pulse_grid（T_range，N_t，l_s）; {定义脉冲网格}
  startpulse_G（E0，tau0，chirp0）; {产生起始脉冲}
  结尾

calc {定义纤维}
  开始
    {主动光纤：}
    set_device（1）;
    set_fiber（L_active，20，'Yb'）;
    add_ring（r_co，N_Yb）;
    泵：= addinputchannel（P_pump_in，l_p，'I_p'，0，转发）;
    signal_active_fw：= addinputchannel（0，l_s，'I_s'，0，forward）;
    signal_active_bw：= addinputchannel（0，l_s，'I_s'，0，向后）;
    set_GVD（GVD_af）;
    set_n2（n2_f）; {非线性指标}
    finish_fiber（）;
    {无源光纤：}
    set_device（2）;
    set_fiber（L_passive，5，'-'）;
    signal_passive_fw：= addinputchannel（0，l_s，'I_s'，0，forward）;
    signal_passive_bw：= addinputchannel（0，l_s，'I_s'，0，向后）;
    set_GVD（GVD_pf）;
    set_n2（n2_f）; {非线性指标}
    finish_fiber（）;
  结尾;

DoResonatorRoundTrip（）：= 
  {模拟脉冲的一个谐振器往返行程。}
  开始
    全球允许所有人；
    pp_multiply_expr_f（'r_FBG％（f）'）;
      {光纤布拉格光栅作为输出耦合器}
    pp_fiber（1，signal_active_fw）; {主动光纤}
    pp_fiber（2，signal_passive_fw）; {无源光纤}
    pp_sat_abs（dR_S，tau_S，E_sat_S）; {SESAM}
    pp_loss（loss_S）; {SESAM损失}
    pp_fiber（2，signal_passive_bw）; {无源光纤}
    pp_fiber（1，signal_active_bw）; {主动光纤}
    pp_center（1）; {使脉动居中}
    calc_dyn（0，T_rt，T_rt）;
      {增益恢复的动态模拟}
  结尾

计算
  对于j：= 1到1000
  开始
    DoResonatorRoundTrip（）;
    store_pulse（j）;
  结尾
显示“能量：”，E_p（）：d3：“ J”
显示“持续时间：”，tau_p（）：d3：“ s”
显示“带宽：”，dl_p（）：d3：“ m”

FindSteadyState（tol，N_max）：= 
  {找到稳定状态，
    由相对变化的时间定义
    E_out和tau的连续三个步骤的总和小于tol。
    如果在N_max个步骤内找到稳定状态，则返回1，
    否则为0。}
  开始
    var n_sc，N，E1，E2，tau1，tau2，df1，df2;
    N：= 0; {往返次数}
    n_sc：= 0;
      {脉冲参数变化很小的往返次数}
    E1：= E_p（）; tau1：= tau_p（）; df1：= df_p（）;
    重复
      inc（N）;
      DoResonatorRoundTrip（）;
      E2：= E_p（）; tau2：= tau_p（）; df2：= df_p（）;
      如果abs（E2 / E1-1）<= tol
         和abs（tau2 / tau1-1）<= tol
         和abs（df2 / df1-1）<= tol
      然后inc（n_sc）
      否则n_sc：= 0;
      E1：= E2; tau1：= tau2; df1：= df2;
    直到n_sc = 3或N> = N_max;
    （n_sc = 3）；
  结尾

“泵功率的变化”
x：400、600
“泵浦功率（mW）”，@ x
y：0、20
框架
 
！对于x：= CS_x1至CS_x2步骤20
  开始
    set_device（1）; {主动光纤}
    set_P_in（泵，x * mW）;
    如果FindSteadyState（0.01，1000）然后
    开始
      点（x + i * E_p（）/ nJ，“ R”）;
        {实心矩形，以nJ为单位的脉冲能量}
      点（x + i * tau_p（）/ ps，“ t”）;
        {脉冲持续时间以ps为单位的空心三角形}
    结尾;
  结尾