分步傅立叶法解广义非线性薛定谔方程

clear all; 
%输入脉冲参数
P0=1e-3; %input pwr in watts
T0=5e-12; %initial pulse half width in second
C=0; %Input chirp parameter for first calculation
alpha=0; % Fiber loss value in dB/km  当传输长度很短时可以不考虑损耗。
% gamma=2*pi*n2/lamda0/Aeffn2=3.2*1e-20w0=2*pic/lamdac=3e8Aeff=20um2;
gamma=0.25; %fiber non linearity in /W/m 用高非线性光纤实现脉冲压缩，非线性系数能达到的量级
b0=1e-27;%1ps2/km=1e-27s2/m
b2=-1*b0; %2nd order disp. （s2/m） 采用DSF色散补偿光纤，即低的色散系数，反常色散。

%特征参数与离散步长
Ld=（T0^2）/（abs（b2））; %dispersion length in meter
Ln=1/（gamma*P0）;
dz=Ln/100;%离散步长可相对非线性长度来取，一般取非线性长度的十分之一以下，即相差一个数量级。
% L=10*Ln;%让非线性发挥作用，实现脉冲压缩。
% nz=L/dz;
%----------------------------------------------------------
A0=sqrt（P0）; %Amplitude
T=-40e-12:1e-12: 39e-12;%  dt=T/T0
dt=1e-12;
u0=A0*exp（-（（1+1i*（-C））/2）*（T/T0）.^2）;%page#47 G.P.AGrawal
U0=u0;
hmw_in=find（abs（u0）>abs（max（u0）/2））;
hmw_in=length（hmw_in）-1;

tol=1e-5;
maxiter=4;
nt = length（u0）;
w = 2*pi*[（0:nt/2-1）（-nt/2:-1）]/（dt