untitled3.m

Study on LD Pumped Yb3+ -doped Double-clad Fiber Lasers数值分析对线形腔LD泵浦掺镱的双包层光纤激光器进行研究

clc;
clear;
p0=0.5533;%赋初值
R1=0.98;%赋初值
R2=0.04;%赋初值

for var1=0:0.001:15
var2=var1/R1;
[x,y]=ode23('fun0088',[0,60],[var1;var2]);
pp=(-p0.*(1+R1)+sqrt(p0^2.*(1+R1)^2+4.*R1.*y(:,1).*(p0+R2.*y(:,1)+p0.*R2)))/(2.*R1);
ppp=(-p0.*(1+R2)+sqrt(p0^2.*(1+R2)^2+4.*R2.*pp.*(p0+R1.*pp+p0.*R1)))/(2.*R2);
pppp=ppp*R2;
if abs(mean((y(:,2)-pppp)))<1;
    break
end
end
g0=8  ;%赋初值
a= 0.299;      ;%赋初值
p0=0.5553 ;%赋初值
psat=1.908e-2;    ;%赋初值
as=0.0025 
p013=300;
p014=0;
R3=0.96;
R4=0.98;
pout11=p013*(exp((-a).*x));
pout12=p014*(exp((-a).*(60-x)));
pout13=pout11+pout12;
go=g0*pout13;
g=go/(1+((y(1)+y(2))/psat));
plot(x,g);
hold on
clc;
clear;
p0=0.5533;%赋初值
R1=0.98;%赋初值
R2=0.04;%赋初值

for var1=0:0.001:15
var2=var1/R1;
[x,y]=ode23('fun006',[0,60],[var1;var2]);
pp=(-p0.*(1+R1)+sqrt(p0^2.*(1+R1)^2+4.*R1.*y(:,1).*(p0+R2.*y(:,1)+p0.*R2)))/(2.*R1);
ppp=(-p0.*(1+R2)+sqrt(p0^2.*(1+R2)^2+4.*R2.*pp.*(p0+R1.*pp+p0.*R1)))/(2.*R2);
pppp=ppp*R2;
if abs(mean((y(:,2)-pppp)))<1;
    break
end
end
g0=8  ;%赋初值
a= 0.299;      ;%赋初值
p0=0.5553 ;%赋初值
psat=1.908e-2;    ;%赋初值
as=0.0025 
p02=150;
p03=150;
pout3=p02*(exp((-a)*x));
pout4=p03*(exp((-a)*(60-x)));
pout5=pout3+pout4;
go=g0*pout5;
plot(x,g);
hold on
clc;
clear;
p0=0.5533;%赋初值
R1=0.98;%赋初值
R2=0.04;%赋初值

for var1=0:0.001:15
var2=var1/R1;
[x,y]=ode23('fun8888',[0,20],[var1;var2]);
pp=(-p0.*(1+R1)+sqrt(p0^2.*(1+R1)^2+4.*R1.*y(:,1).*(p0+R2.*y(:,1)+p0.*R2)))/(2.*R1);
ppp=(-p0.*(1+R2)+sqrt(p0^2.*(1+R2)^2+4.*R2.*pp.*(p0+R1.*pp+p0.*R1)))/(2.*R2);
pppp=ppp*R2;
if abs(mean((y(:,2)-pppp)))<1;
    break
end
end

g0=8  ;%赋初值
a= 0.299;      ;%赋初值
p0=0.5553 ;%赋初值
psat=1.908e-2;    ;%赋初值
as=0.0025 
p021=60;
p031=240;
pout31=p021*(exp((-a).*x));
pout41=p031*(exp((-a).*(60-x)));
pout51=pout31+pout41;
go=g0.*pout51;
g=go/(1+((y(1)+y(2))/psat));
plot(x,g)
hold on
clc;
clear;
p0=0.5533;%赋初值
R1=0.98;%赋初值
R2=0.04;%赋初值

for var1=0:0.001:15
var2=var1/R1;
[x,y]=ode23('fun998',[0,60],[var1;var2]);
pp=(-p0.*(1+R1)+sqrt(p0^2.*(1+R1)^2+4.*R1.*y(:,1).*(p0+R2.*y(:,1)+p0.*R2)))/(2.*R1);
ppp=(-p0.*(1+R2)+sqrt(p0^2.*(1+R2)^2+4.*R2.*pp.*(p0+R1.*pp+p0.*R1)))/(2.*R2);
pppp=ppp*R2;
if abs(mean((y(:,2)-pppp)))<1;
    break
end
end

g0=8  ;%赋初值
a= 0.299;      ;%赋初值
p0=0.5553 ;%赋初值
psat=1.908e-2;    ;%赋初值
as=0.0025 
p021=0;
p031=300;
pout31=p021.*(exp((-a).*x));
pout41=p031*(exp((-a).*(60-x)));
pout51=pout31+pout41;
go=g0.*pout51;
g=go/(1+((y(1)+y(2))/psat));
plot(x,g)