非线性光纤光学

%%观察群速度色散和kell效应同时扰动作用下多孤子的传输特性

clear
clc

%?注意时间域与频域的选择问题
T = 40;
tn= 2^10;  % 时间区域个数（2^M个）                            
dt = T/tn;   %选择的时间步长  
t = ((1:tn)'-(tn+1)/2)*dt;  
w = wn(T,tn);%傅里叶域频率,调用计算频率函数 

% m=ran(0,1)%n为0-1之间随机变化数
% r=ran(-1,1)%l为-1和1之间随机变化数
% %上面两个式子调用了随机数产生函数
m0=[0,0.2,0.3,1];
r=0.2;
%由于m和r都是随机变化的量，m*r也是随机变化的量，我们假设r为定值，研究随m变化孤子的传输特性

b=0.03;
%c=-i*(1+m*r)*w'.^2/2+i*b*w'.^3;
%步长的选择影响计算效率
dz=0.1;  %迭代时选择的步长
z=0:dz:100;  % 位置坐标  
%z=0:dz:50;  % 位置坐标  

% n=ran(0,1)%n为0-1之间随机变化数
% l=ran(-1,1)%l为-1和1之间随机变化数
% %上面两个式子调用了随机数产生函数
n=0.1;
l=0.1;
%由于kt中我们已经研究了扰动系数的变化对孤子传输的影响，此时我们给kell系数一个固定的量来研究随着群速度色散的随机扰动孤子的传输特性有什么变化


% 输入初始脉冲
q=3.5;
r1=1.3;r2=1.1;r3=1.2;
o1=0;o2=0;o3=0;%q为初始间距，r为相对振幅,o为相对相位
u(1,:)=sech(t-q)+r1*sech(r1*(t+q))*exp(i*o1)+r2*sech(r2*(t+3*q))*exp(i*o2)+r3*