我试图使用mathematica解决两级的光学Bloch方程。
这是我写的代码。我想解决Rho12,看看Rho12的虚构部分在失谐[Capital Delta]方面的变化。
显示错误“他们是超验方程式” 有人可以帮助我......
\[CapitalOmega]21 = Conjugate[\[CapitalOmega]12];
n := \[Rho]11 + \[Rho]22 + \[Rho]12 + \[Rho]21 = 1;
Solution == DSolve[{
\[Rho]11'[t] == (\[CapitalGamma]2*\[Rho]22[t]) +
i/2 (\[CapitalOmega]21*\[Rho]12[t] - \[CapitalOmega]12*\[Rho]21[
t]),
\[Rho]22'[t] == -\[CapitalGamma]2*\[Rho]22[t] +
i/2 (\[CapitalOmega]12*\[Rho]21[t] - \[CapitalOmega]21*\[Rho]12[
t]),
\[Rho]12'[
t] == (-(\[CapitalGamma]2/2) - i*\[CapitalDelta])*\[Rho]12[t] +
i/2*\[CapitalOmega]12*(\[Rho]11[t] - \[Rho]22[t]),
\[Rho]21'[t] == Conjugate[\[Rho]12'[t]] ,
\[Rho]22[0] == 0, \[Rho]11[0] == 1, \[Rho]12[0] ==
0, \[Rho]21[0] == 0}, {\[Rho]11[t], \[Rho]22[t], \[Rho]12[
t], \[Rho]21[t]}, t];