拉曼光纤激光器超连续谱（supercontinuum）仿真程序

拉曼光纤激光器超连续谱（supercontinuum）仿真程序

clear all
close all
format long


% parameters

zmax=0.1;                   % fiber segment [m]  
Tw=24e-12;                % temporal window [s]
n=2048*2;                   % fft points
lambda=1550*1e-9;           % pump wavelength [m]
sm=1;                       % SPM [1/W/m]
cm=0;                       % XPM fraction
iter=5000;                  % split-step iterations for 0-zmax
usa=200;                   % Store “usa“ intermediate results
Pp=300;                    % peak power [W]
pwidth=0.4e-12;             % FWHM pulsewidth [s]
lin_pol_angle=0;            % polarization angle [rad]
seedamp=0;	                % seed spectral power [W] （0=no noise）
fr=0.18;	                    % Raman fraction
tau1=12.2e-15;	            % Raman response [s]
tau2=32e-15;	            % Raman response [s]
loss=0;                     % Fiber loss [dB/m]
cohere_yes=0;               % include coherent terms yes=1 no=0;
slfstp_yes=1;               % include self steepening yes=1 no=0;
ramand_yes=1;               % include Raman nl yes=1; no=0;


% loading dispersion curve
[pl_fib_dompl_fib_dx]=dispersion_tellurite_pcf_1;
[pl_fib_dompl_fib_dy]=dispersion_tellurite_pcf_1;


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% do not change %%%%%%%%%%%%%%%%%%%%%%%%%%
% definition for phase birefringence
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% definition of domains %%%%%%%%%%%%%%%%%%
unds=1;                  % undersampling factor for contourplot  
poly_ord=14;             % polynomial expansion of GVD
c=3e8;                   % speed of light [m/s]
ssfram_yes=1;            % Self steepening molecular yes=1; no=0;
dnf_ut=0.0;
% definition for group birefringence
dnfg_ut=0.0;
Deltan=dnfg_ut;
Deltabeta=Deltan/（2*c）;
csi=zmax*sm;
tmax=Tw/2;
Tnorm=tmax;
deltat=Tw/n;
De=Deltabeta/（Tnorm*sm）;
bfase=dnf_ut*2*pi/lambda/sm;
t=-tmax/Tnorm:deltat/Tnorm:（tmax-deltat）/Tnorm;
deltaz=csi/iter;
indfreq=-n/2:1:n/2-1;
kt=（Tnorm*pi/（tmax））.*indfreq;
ft=kt/（2*pi）;
ftot=c/lambda+ft/Tnorm;
dolam=c./ftot*1e9;
w0=2*pi*c/lambda*Tnorm*n;
tau1_norm=tau1/Tnorm;
tau2_norm=tau2/Tnorm;
loss_norm=1/2/sm*log（10^（loss/10））;
z_axis=0:zmax/usa:zmax;
Qx=[];
zv=[];
Qxf=[];
mat_x=[];
mat_y=[];
mat_dom=[];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% To change dispersion coefficient into the group velocity dispersion
% and to fit the available data with the BPM frequency domain
% please pay attention that the BPM frequency range is included in your
% available dispersion dat

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        7706  2011-11-04 12:52  supercontinuum.m