or m2=1: fix(z/h)2 a=(abs(u)).2+(abs(v)).2 w=sum((abs(u)).2)+sum((abs(v)).2)*Dt %gw=sanl(a, T, NUM) gk=(1+gw/Es) gp=gO /(1+(w(1: NUM)/(omega)).6)/gk: according to equation(10) Dl=exp(- deltas*(i*W1)+i*0.5*dsp2*(*1).2)*h/2); ul=fft(u).**DI 2=ifft(ul) D2=exp(+ deltas(i*W1)+1*0.5*disp2*(i*w1).2)*h/2); v1=fft(v).*D2 v2=ifft(v1) N1=i*gama*(abs(u)).2+2*(abs(v).2/3+conj(u).*V.2./(uteps+i*eps)/3)+i*beta+gp 12=i*gama*((abs(v)). 2+2*(abs(u)).2/3+conj(v).*u. 2. /(v+eps+isps)/3)-i*betatgp u3=u12.*eXp(N1.*h) v3=v2. *exp (N2 *h) D4=exp((- deltas*(i*W1)+1*0.5*disp2*(i*w1).2)*h/2); u4=fft(u3).*D4; u5=ifft(u4 D5=exp((deltaj* (i*w1)+i*0. 5*disp2*(i=w1).2)*h/2) v5=ifft(v4) u=u5 %% propogation in standard single-mode fiber, 3m long Lb=2*1e-3;%12/5*1e-3Lb=12/5*1e-3;%10m1ong beta=pi/ Lb delta j=beta*1. 554*le-9/(2*pi*3e-7)for m2=1:fix(z/h)2 a=(abs(u)).^2+(abs(v)).^2; gw=sum((abs(u)).^2)+sum((abs(v)).^2)*Dt; %gw=san1(a,T,NUM); gk=(1+gw/Es); gp=g0./(1+(w(1:NUM)/(omegag)).^6)/gk; % according to equation (10) D1=exp((-deltaj*(i*w1)+i*0.5*disp2*(i*w1).^2)*h/2); u1=fft(u).*D1; u2=ifft(u1); D2=exp((+deltaj*(i*w1)+i*0.5*disp2*(i*w1).^2)*h/2); v1=fft(v).*D2; v2=ifft(v1); N1=i*gama*((abs(u)).^2+2*(abs(v)).^2/3+conj(u).*v.^2./(u+eps+i*eps)/3)+i*beta+gp /2; N2=i*gama*((abs(v)).^2+2*(abs(u)).^2/3+conj(v).*u.^2./(v+eps+i*eps)/3)-i*beta+gp /2; u3=u2.*exp(N1.*h); v3=v2.*exp(N2.*h); D4=exp((-deltaj*(i*w1)+i*0.5*disp2*(i*w1).^2)*h/2); u4=fft(u3).*D4; u5=ifft(u4); D5=exp((deltaj*(i*w1)+i*0.5*disp2*(i*w1).^2)*h/2); v4=fft(v3).*D5; v5=ifft(v4); u=u5; v=v5; end %%% propogation in standard single-mode fiber, 3m long z=3*1e-3; Lb=2*1e-3;%12/5*1e-3Lb=12/5*1e-3; % 10m long beta=pi/Lb; deltaj=beta*1.554*1e-9/(2*pi*3e-7); disp2=23;