simulation_激光光纤_光纤激光器_光纤锁模_激光器_锁模光纤

% Run this scrip to start the simulation % The round trip loop is defined as: Output - > SESAM -> Fiber -> Output parameter; disp(sprintf('Now counting...')) if Clear_flag delete('datas\f\*.mat'); delete('datas\rt\*.mat'); disp(sprintf('Warning: All former datas have been cleared!')) disp(sprintf('\n')) end tic for rt = Start_loop : Cycle_time disp(sprintf('Round trip: %d',rt)) % Output coupler A_oc = sqrt(1-Attenu_sys)*A_in; %A_in初始脉冲 U=sqrt(rho)*A_oc; V=i*sqrt(1-rho)*A_oc; %fiber [V_1]=SMF_ccw(V,T_span,Beta_smf,Gama,T_2,L,H,Save_step); [V_2]=PCF_ccw(V_1,T_span,Beta_2,Gama,G_0,E_l,T_2,L,H,Save_step); [U_1]=PCF_cw(U,T_span,Beta_2,Gama,G_0,E_l,T_2,L,H,Save_step); [U_2]=SMF_cw(U_1,T_span,Beta_smf,Gama,T_2,L,H,Save_step) %coupler A_c=sqrt(rho).*U_2+sqrt(1-rho).*V_2 ; [A_final]=SMF_left(A_c,T_span,Beta_2smf,Gama,T_2,L,H,Save_step); E_next = trapz(T,abs(A').^2); disp(sprintf('Pulse energy after DDL is: %g nJ',E_next/1000)) Differ = abs(sum(abs(A').^2-abs(A_in).^2)/sum(abs(A_in).^2)); disp(sprintf('Relative change: %g',Differ)) disp(sprintf('\n')) A_in=A_final; clear U_1 U_2 V_1 V_2 A_final U V A_oc end toc