%基于HNLF的全光邏輯與門主程序
clear;clf;
global?tmax?tmin?delta_t?delta_r?????????????%全局變量
zmin=0;zmax=1450;delta_z=1;
tmin=-25;tmax=25;delta_t=1;
rmin=0;rmax=18;delta_r=1;
Leff=135;
z=zmin:delta_z:zmax;?????????????????????%長度單位:mm
t=tmin:delta_t:tmax;??
r=rmin:delta_r:rmax;
%時(shí)間單位:ps
c=2.998e8;??
a=0;%光速單位:m/s

Pp0=0.8;????????????????????????????????%泵浦峰值功率?單位:W
Ps0=0.08;

Ap=zeros（（zmax-zmin）/delta_z+1（tmax-tmin）/delta_t+1）;??%分配空間，脈沖振幅
As=zeros（（zmax-zmin）/delta_z+1（tmax-tmin）/delta_t+1）;
Ac=zeros（（zmax-zmin）/delta_z+1（tmax-tmin）/delta_t+1）;
Pp=zeros（（zmax-zmin）/delta_z+1（tmax-tmin）/delta_t+1）;??%脈沖功率
Ps=zeros（（zmax-zmin）/delta_z+1（tmax-tmin）/delta_t+1）;
Pc=zeros（（zmax-zmin）/delta_z+1（tmax-tmin）/delta_t+1）;
%Gpeak=zeros（（zmax-zmin）/delta_z+1（tmax-tmin）/delta_t+1）;
Jp=zeros（1（zmax-zmin）/delta_z+1）;??????????????????????%脈沖能量
Js=zeros（1（zmax-zmin）/delta_z+1）;
Jc=zeros（1（zmax-zmin）/delta_z+1）;

tau_0=5;????????????????????????????????????????????????%初始脈寬
Ap（1:）=sqrt（Pp0）*sech（1.76*t/tau_0）;
As（1:）=sqrt（Ps0）*sech（1.76*t/tau_0）;
Ac（1:）=0;
Pp（1:）=（abs（Ap（1:）））.^2;
Ps（1:）=（abs（As（1:）））.^2;
bp=1.5528;??????????????????????????????%波長單位:um
bs=1.548;
bc=bp*bs/（2*bs-bp）;
bc

u0=4*pi*1e-7;???????????????????????????%磁導(dǎo)率?單位:N*s^2/C^2
Aeff=11;????????????????????????????????%有效相互作用面積?單位:um^2
%r=20;
rp=r*1.55/bp;?????