被动调Q_高等光学MATLAB仿真_激光速率方程_

function dy=rate_eq(t,y,flag,Rp,T0,R) sigma=5.4e-23; %增益介质的受激发射截面 sigma_gs=8.7e-23; %可饱和吸收体基态的吸收截面 sigma_es=2.2e-23; %可饱和吸收体激发态的吸收截面 N_T=1.68e26; %增益介质的总粒子数密度 tao_a=750e-6; %增益介质的上能级寿命 tao_gs=3e-6; %可饱和吸收体的恢复时间 n1=1.82; %Nd3+:YAG的折射率 n2=1.80; %Cr4+:YAG的折射率 delta=0.02; %激光器腔体的耗散性损耗 l=0.001; %增益介质的长度 ls=0.001; %可饱和吸收体沿光腔轴线的长度 gamma=1; %反转因子，对于四能级系统为1 c=2.997963e8; %真空中的光速 lc=n1*l+n2*ls; %谐振腔等效光程长度 tr=lc/c; %光在腔中往返一周的时间 n0s=-log(T0)/(sigma_gs*ls); %求可饱和吸收体粒子数密度 y(1)=max(y(1),1); %光子数密度的最小值 %被动调Q耦合方程组 dy=[y(1)*(2*sigma*y(2)*l-2*sigma_gs*y(3)*ls-2*sigma_es*... (n0s-y(3))*ls-(log(1/R)+delta))/tr; Rp*(1-y(2)/N_T)-gamma*sigma*c*y(1)*y(2)-y(2)/tao_a; (n0s-y(3))/tao_gs-sigma_gs*c*y(1)*y(3)];