以下程序存在以Laser2为名的.m文件中。

function dy=Laser2(t,y)
Delta=2.8e-23;%激光增益介质受激发射截面面积
Deltags=1.1e-22;%饱和吸收体的基态吸收截面面积
Deltaes=1.2e-23;%饱和吸收体的激发态吸收截面面积
tg=2.3e-4;%饱和吸收体的恢复时间
ts=3.4e-6;
C=3e8;%真空中光速
P=0.2;%泵浦光的泵浦功率
L=0.047;%谐振腔的光学长度
lg=1.2e-3;%激光增益介质长度
ls=8e-4;
WL=808e-9;%泵浦光波长
V=C/WL;%泵浦光频率
h=6.63e-34;%普朗克常数
T0=0.9;%可饱和吸收体的初始透过率
tr=2*L/C;%光在谐振腔中往返一周所需时间
asc=0.2;%激光器的散射损耗，典型值为0.002cm-1
R=0.95;%输出镜反射率
Wp=5e-5;%泵浦光的光斑半径
Rp=P/(h*V*pi*(Wp^2)*(lg+ls));
dy=zeros(3,1);
dy(1)=y(1)/tr*(2*Delta*y(2)*lg-2*Deltaes*((log(1/T0))/(Deltags*(lg+ls))-y(3))*ls-(log(1/R)+2*asc*(lg+ls)));
dy(2)=Rp-Delta*C*y(1)*y(2)-y(2)/tg;
dy(3)=-Deltags*C*y(1)*y(3)+((log(1/T0))/(Deltags*(lg+ls))-y(3))/ts;
end

以下依次程序在command命令窗口运行。

[t,y]=ode45('Laser2',[0 5],[0,0,0]);
figure,plot(t,y(:,1),'r');
figure,plot(t,y(:,2),'b');
figure,plot(t,y(:,3),'g');

我的电脑性能不是太好，貌似运行好久也出不来，不知道哪位大侠能帮忙跑一下，把结果给我啦！！感激不尽啊，急等！谢谢了！