FWM 耦合方程 matlab代码

%基于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;
%时间单位: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;                           %磁导率 单位:N*s^2/C^2
Aeff=11;                                %有效相互作用面积 单位:um^2
%r=20;
rp=r*1.55/bp;