半波片对偏振态的影响： import numpy as npimport qutip# 这是半波片对应的SU(2)操作def hwp(state, theta):rn_mat = -1*qutip.Qobj([[np.cos(2*theta),np.sin(2*theta)],[np.sin(2*theta),-np.cos(2*theta)]])r_state = (rn_mat*sta

上一篇文章说到python的qutip包可以实现布洛赫球的可视化，并附上了产生动画的代码。 这一篇文章可以对SO(3)旋转操作进行可视化：(需要用到上一篇的new_bloch3d) import numpy as npimport qutip# SO(3)旋转矩阵总共包含三个参数theta, phi, deltadef rn(state, theta, phi, delta):rn_mat