这篇文章运用简单易懂的例子给大家介绍Python中qutip的用法有哪些,内容非常详细,感兴趣的小伙伴们可以参考借鉴,希望对大家能有所帮助。
让客户满意是我们工作的目标,不断超越客户的期望值来自于我们对这个行业的热爱。我们立志把好的技术通过有效、简单的方式提供给客户,将通过不懈努力成为客户在信息化领域值得信任、有价值的长期合作伙伴,公司提供的服务项目有:域名注册、网站空间、营销软件、网站建设、白沙黎族网站维护、网站推广。一、N原子系综自旋概率分布
from qutip import * import numpy as np import matplotlib.pyplot as plt n=2#原子数 j = n//2 psi0 = spin_coherent(j, np.pi/3, 0)#设置系统的初态为自旋相干态 Jp=destroy(2*j+1).dag()#升算符 J_=destroy(2*j+1)#降算符 Jz=(Jp*J_-J_*Jp)/2#Jz H=Jz**2#系统的哈密顿量 tlist=np.linspace(0,3,100)#时间列表 result=mesolve(H,psi0,tlist)#态随时间的演化 theta=np.linspace(0, np.pi, 50) phi=np.linspace(0, 2*np.pi, 50) #分别计算四个状态下的 husimi q函数 Q1, THETA1, PHI1 = spin_q_function(result.states[0], theta, phi) Q2, THETA2, PHI2 = spin_q_function(result.states[30], theta, phi) Q3, THETA3, PHI3 = spin_q_function(result.states[60], theta, phi) Q4, THETA4, PHI4 = spin_q_function(result.states[90], theta, phi) #在四个子图中分别画出四个状态下的husimi q函数 fig = plt.figure(dpi=150,constrained_layout=1) ax1 = fig.add_subplot(221,projection='3d') ax2 = fig.add_subplot(222,projection='3d') ax3 = fig.add_subplot(223,projection='3d') ax4 = fig.add_subplot(224,projection='3d') plot_spin_distribution_3d(Q1, THETA1, PHI1,fig=fig,ax=ax1) plot_spin_distribution_3d(Q2, THETA2, PHI2,fig=fig,ax=ax2) plot_spin_distribution_3d(Q3, THETA3, PHI3,fig=fig,ax=ax3) plot_spin_distribution_3d(Q4, THETA4, PHI4,fig=fig,ax=ax4) for ax in [ax1,ax2,ax3,ax4]: ax.view_init(0.5*np.pi, 0) ax.axis('off')#不显示坐标轴 fig.show()