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