cIc tic 1=3*pi/8 %input( theta=') detphi=0*pi/50: %input( detphi') Gi P theta2=thetal+pi/2-0 NUM=1024 t=linspace(-T, T, NUM+1) %qu=0. 5kexp(t(1: NUM). 2/100000/2): %quasi-CW, pulse width Ins=1000ps qu=0. 1*rand (1, NUM) %qu=0.1* randn(1,1024) %qu=0. 1*imnoise(1: 1024, gaussian') %qu=wgn(1, 1024, 0.01):% Gaussian noise Dt=2*T/NUM w=linspace(-pi /Dt, pi. /Dt, NUM+1) wl=fftshift(w(1: NUM)) tt=linspace(T, T, NUM) the units are W, m, ps h=0. 1*le-3: calculation step unit ki for m=l: nn nn number of round trips etal) v=qu*sin(theta1)*exp (-i*detphi) %% propogation in standard fiber, 5m long Lb=2*le-3: %12/5*1e-3: according to the paper it's 10m long beta=pi/Lbclc clear all tic nn=input('nn=') theta1=3*pi/8; %input('theta1=') detphi=0*pi/50; %input('detphi') % Given Parameters gama=3; % 1/W/km theta2=theta1+pi/2-0.05; NUM=1024; T=30; %ps t=linspace(-T,T,NUM+1); %qu=0.5*exp(-t(1:NUM).^2/100000/2); %quasi-CW, pulse width 1ns=1000ps qu=0.1*rand(1,NUM); %qu=0.1*randn(1,1024); %qu=0.1*imnoise(1:1024,'gaussian'); %qu=wgn(1,1024,0.01); % Gaussian noise Dt=2*T/NUM; w=linspace(-pi./Dt,pi./Dt,NUM+1); w1=fftshift(w(1:NUM)); tt=linspace(-T,T,NUM); % the units are W,m,ps h=0.1*1e-3; % calculation step unit km for m=1:nn % nn number of round trips u=qu*cos(theta1); v=qu*sin(theta1)*exp(-i*detphi); %%% propogation in standard fiber, 5m long z=5*1e-3; Lb=2*1e-3;%12/5*1e-3; % according to the paper it's ~10m long beta=pi/Lb;