怎么用python和qutip写wigner函数（相空间）？

例如稳态的密度矩阵为ρ=0.5*(|α><α|+|-α><-α|)，在n=20的情况下，写出他的魏格纳函数。我需要考虑哈密顿量演化吗？还算是可以直接写？α=1.3+1.2j