計算脈沖在非線性耦合器中演化的Matlab 程序 ~K96y$ DTE E?mp6R]}% % This Matlab script file solves the coupled nonlinear Schrodinger equations of
B|=maz:_ % soliton in 2 cores coupler. The output pulse evolution plot is shown in Fig.1 of
5r<(Z0 % Youfa Wang and Wenfeng Wang, “A simple and effective numerical method for nonlinear
+I n"OR% % pulse propagation in N-core optical couplers”, IEEE Photonics Technology lett. Vol.16, No.4, pp1077-1079, 2004
zrO|L|F&P ug,|'<G+ %fid=fopen('e21.dat','w');
V^z;^mdd N = 128; % Number of Fourier modes (Time domain sampling points)
D,j5k3< # M1 =3000; % Total number of space steps
wKlCx J =100; % Steps between output of space
{Yk20Zn T =10; % length of time windows:T*T0
} XU:DE T0=0.1; % input pulse width
--YUiNhh MN1=0; % initial value for the space output location
S1`0d9ds# dt = T/N; % time step
Iq]6] n = [-N/2:1:N/2-1]'; % Index
~ p.W*skD t = n.*dt;
"T%'Rp`j| u10=1.*sech(1*t); % input to waveguide1 amplitude: power=u10*u10
-!>ZATL<B u20=u10.*0.0; % input to waveguide 2
%QgAilj, u1=u10; u2=u20;
|a$w;s>\ U1 = u1;
eKz~viM' U2 = u2; % Compute initial condition; save it in U
SWsv, ww = 4*n.*n*pi*pi/T/T; % Square of frequency. Note i^2=-1.
*.c9$`s w=2*pi*n./T;
Q|KD/s?? g=-i*ww./2; % w=2*pi*f*n./N, f=1/dt=N/T,so w=2*pi*n./T
XHM"agrhSQ L=4; % length of evoluation to compare with S. Trillo's paper
&[Zap6] dz=L/M1; % space step, make sure nonlinear<0.05
zZc@;S# for m1 = 1:1:M1 % Start space evolution
T;r];Y(b* u1 = exp(dz*i*(abs(u1).*abs(u1))).*u1; % 1st sSolve nonlinear part of NLS
64'
涿州市|
富民县|
武汉市|
龙江县|
宁蒗|
吴桥县|
成武县|
平谷区|
从江县|
北宁市|
大埔县|
广南县|
七台河市|
孝昌县|
潞西市|
贡嘎县|
宁化县|
喀喇|
玉龙|
德惠市|
丹东市|
武穴市|
乐昌市|
农安县|
民乐县|
江华|
兴城市|
荔浦县|
吴川市|
荆门市|
新巴尔虎左旗|
高尔夫|
湘潭县|
托里县|
郁南县|
玉龙|
麻栗坡县|
分宜县|
南京市|
西贡区|
温宿县|