两种方法,一个是使用matlab内嵌ode45函数,一个是使用数值迭代runge-butte的方法。
1.常规的三能级激光系统速率方程表示如下:
{ d N 3 d t = N 1 W 13 − N 3 ( S 32 + A 31 ) d N 2 d t = − ( N 2 − g 2 g 1 N 1 ) σ 21 v ϕ − N 2 ( A 21 + S 21 ) + N 3 S 32 N 1 + N 2 + N 3 = N T d ϕ d t = ( N 2 − g 2 g 1 N 1 ) σ 21 v ϕ − ϕ τ L \left\{ \begin{aligned} \frac{\mathrm{d} N_{3}}{\mathrm{~d} t}&=N_{1} W_{13}-N_{3}\left(S_{32}+A_{31}\right) \\ \frac{\mathrm{d} N_{2}}{\mathrm{~d} t}&=-\left(N_{2}-\frac{g_{2}}{g_{1}} N_{1}\right) \sigma_{21} v \phi- N_{2}\left(A_{21}+S_{21}\right)+N_{3} S_{32} \\ N_{1}&+N_{2}+N_{3}=N_{\mathrm{T}} \\ \frac{\mathrm{d} \phi}{\mathrm{d} t}&=\left(N_{2}-\frac{g_{2}}{g_{1}} N_{1}\right) \sigma_{21} v \phi-\frac{\phi}{\tau_{\mathrm{L}}} \end{aligned} \right. ⎩⎪⎪⎪⎪⎪⎪⎪⎪⎪⎨⎪⎪⎪⎪⎪⎪⎪⎪⎪⎧ dtdN3