clc
clear;
L=0:5:30;%增益光纤长度
Ces=2*10^(-25);%激光的发射截面
As=0.8;%激光光场重叠积分
R1=0.99;%左腔镜反射率
R2=0.04;%右腔镜反射率
Cap=2.5*10^(-24);%泵浦光吸收截面
Cep=2.5*10^(-24);%泵浦光发射截面
h=6.62606896*10^(-34);%普朗克常数
Cas=10^(-27);%激光的吸收截面
tao=770*10^(-9);%上能级粒子寿命
A=6.25*pi;%纤芯有效截面积
lamda1=975*10^(-9);%当泵浦波长为975nm。
Vp=(3*10^8)/lamda1;%当泵浦波长为975nm时对应的泵浦频率.
Ap=1/17^2;%泵浦光场重叠积分，Ap=Sd/Sic Sd和Sic分别为增益光纤纤芯和内包层的截面积,纤芯直径5um，内包层直径85um。
Rou=3/100;%镱离子掺杂浓度（未知）
bita=Ap*Rou*(Cap+Cep)*((log(1/R1/R2))/(2*As*Rou)-(Ces*Cap-Cep*Cas).*L/(Cap+Cep))/(Cas+Ces);
P=(h*Vp*A)/Ap/(Cap+Cep)/tao;
Pth=P*(bita+Ap*Rou*Cap.*L)/(1-exp(bita));%Pth为泵浦阈值功率
plot(L,Pth);
