Er Yb codoped fier laser.zip_Er激光器_doped fiber_laser_matlab_激光速率

function laser_New %铒镱共掺光纤激光器的传输方程 global sigma_12 sigma_21 sigma_13 sigma_56 sigma_65 Ner.... Nyb alpha_p alpha_s gamma_p gamma_s h c nu_p nu_s Ccr Acore tau_er tau_yb lambda_s Es Ep Pp0 Ps0 %参量设置 lambda_p=980*1e-9; %泵浦光波长1 lambda_s=1550*1e-9; %激光波长 tau_er=10*1e-3; %Er上能级粒子寿命 tau_yb=1.35*1e-3; %Yb上能级粒子寿命 sigma_12=1.75*1e-25; %铒离子的吸收截面面积 sigma_21=2.4*1e-25; %铒离子的发射截面面积 sigma_13=1.97*1e-27; sigma_56=2.19*1e-25; %镱离子的吸收截面面积 sigma_65=3.14*1e-26; %镱离子的发射截面面积 Ner=4.33*1e25; %铒离子的浓度 Nyb=3.61*1e26; %镱离子的浓度 alpha_p=0.15; %双包层光纤对泵浦光的损耗 alpha_s=0.1; %双包层光纤对激光的损耗 gamma_p=0.0058; %泵浦光重叠因子 gamma_s=1; %激光重叠因子 Ccr=1e-22; %交叉弛豫系数 Acore=25*pi*1e-12; %纤芯的截面积 Pp0=0:0.01:20; Ps0=0.15; %物理常量及中间过程参数计算 c=3*1e8; h=6.626*1e-34; nu_s=c/lambda_s; nu_p=c/lambda_p; Es=h*nu_s; Ep=h*nu_p; G=zeros(2001,1); %options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]); %OPTION = odeset('Stats','ON'); for i=1:2001 [t,y]=ode45(@singlepump,[0,10],[Pp0(i),Ps0]); %指令ode45解方程 G(i)=10*log10(y(2,end)/Ps0); end plot(Pp0,G)