Я нашел в Руководстве по подсказкам Простую программу для визуализации периодического поведения JCM.
Мой код
from qutip import *
from pylab import *
import matplotlib.pyplot as plt
import numpy as np
N=10
a=tensor(destroy(N),qeye(2))
sz=tensor(destroy(N),sigmaz())
s=tensor(qeye(N),destroy(2))
wc=wq=2*np.pi
g=np.pi
psi0=tensor(basis(N,0),basis(2,0))
H = wc*a.dag()*a-0.5*wq*sz+0.5*g*(a*s.dag()+a.dag()*s)
tlist=np.linspace(0,3,100)
out=mesolve(H,psi0,tlist,[],[a.dag()*a])
plot(tlist,out.expect[0])
plt.show()
И на картинке я должен увидетьесть синусная линия, но в этой программе я получаю только одну узкую линию.Что не так?