clear?all
close?all
format?long


%?parameters

zmax=0.1;???????????????????%?fiber?segment?[m]??
Tw=24e-12;????????????????%?temporal?window?[s]
n=2048*2;???????????????????%?fft?points
lambda=1550*1e-9;???????????%?pump?wavelength?[m]
sm=1;???????????????????????%?SPM?[1/W/m]
cm=0;???????????????????????%?XPM?fraction
iter=5000;??????????????????%?split-step?iterations?for?0-zmax
usa=200;???????????????????%?Store?“usa“?intermediate?results
Pp=300;????????????????????%?peak?power?[W]
pwidth=0.4e-12;?????????????%?FWHM?pulsewidth?[s]
lin_pol_angle=0;????????????%?polarization?angle?[rad]
seedamp=0;	????????????????%?seed?spectral?power?[W]?（0=no?noise）
fr=0.18;	????????????????????%?Raman?fraction
tau1=12.2e-15;	????????????%?Raman?response?[s]
tau2=32e-15;	????????????%?Raman?response?[s]
loss=0;?????????????????????%?Fiber?loss?[dB/m]
cohere_yes=0;???????????????%?include?coherent?terms?yes=1?no=0;
slfstp_yes=1;???????????????%?include?self?steepening?yes=1?no=0;
ramand_yes=1;???????????????%?include?Raman?nl?yes=1;?no=0;


%?loading?dispersion?curve
[pl_fib_dompl_fib_dx]=dispersion_tellurite_pcf_1;
[pl_fib_dompl_fib_dy]=dispersion_tellurite_pcf_1;


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%?do?not?change?%%%%%%%%%%%%%%%%%%%%%%%%%%
%?definition?for?phase?birefringence
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%?definition?of?domains?%%%%%%%%%%%%%%%%%%
unds=1;??????????????????%?undersampling?factor?for?contourplot??
poly_ord=14;?????????????%?polynomial?expansion?of?GVD
c=3e8;???????????????????%?speed?of?light?[m/s]
ssfram_yes=1;????????????%?Self?steepening?molecular?yes=1;?no=0;
dnf_ut=0.0;
%?definition?for?group?birefringence
dnfg_ut=0.0;
Deltan=dnfg_ut;
Deltabeta=Deltan/（2*c）;
csi=zmax*sm;
tmax=Tw/2;
Tnorm=tmax;
deltat=Tw/n;
De=Deltabeta/（Tnorm*sm）;
bfase=dnf_ut*2*pi/lambda/sm;
t=-tmax/Tnorm:deltat/Tnorm:（tmax-deltat）/Tnorm;
deltaz=csi/iter;
indfreq=-n/2:1:n/2-1;
kt=（Tnorm*pi/（tmax））.*indfreq;
ft=kt/（2*pi）;
ftot=c/lambda+ft/Tnorm;
dolam=c./ftot*1e9;
w0=2*pi*c/lambda*Tnorm*n;
tau1_norm=tau1/Tnorm;
tau2_norm=tau2/Tnorm;
loss_norm=1/2/sm*log（10^（loss/10））;
z_axis=0:zmax/usa:zmax;
Qx=[];
zv=[];
Qxf=[];
mat_x=[];
mat_y=[];
mat_dom=[];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%?To?change?dispersion?coefficient?into?the?group?velocity?dispersion
%?and?to?fit?the?available?data?with?the?BPM?frequency?domain
%?please?pay?attention?that?the?BPM?frequency?range?is?included?in?your
%?available?dispersion?dat

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????7706??2011-11-04?12:52??supercontinuum.m