laser18.m

fiber-laser.m为光纤激光器中的功率分布；laser18.m为fiber-laser.m中的调用函数

% 微分方程  
function f=laser4(z,p)
global input 
wavelength1=975e-9;    
wavelength2=1086e-9;   
t=1e-3;                
aap=2.5e-20;         
aep=2.5e-20;        
aes=2.5e-21;         
aas=1.6e-23;      
A=pi*(29.7/2)*1e-8;            
n=8.93e19           %掺杂浓度
     ap=5e-5; 
tp=0.0043       %重叠因子
ts=0.8;          %重叠因子
pp=input*exp(-(tp*aap*n+ap)*(z));
n2=n*(pp*aap*tp/(6.62e-34*3e8/wavelength1*A)+ts*aas*(p(1)+p(2))/(6.62e-34*3e8/wavelength2*A))/(pp*(aap+aep)*tp/(6.62e-34*3e8/wavelength1*A)+1/t+ts*(aes+aas)*wavelength2*(p(1)+p(2))/(6.62e-34*3e8*A));
f=[ts*((aes+aas)*n2-aas*n)*p(1)-5e-5*p(1) -ts*((aes+aas)*n2-aas*n)*p(2)+5e-5*p(2)]';