The Wigner Function of $x(t)$ is
$$ W(t,f) = \frac{1}{2\pi}\int x\Big(t+\frac{\tau}{2}\Big)x^*\Big(t-\frac{\tau}{2}\Big) e^{-j2\pi f\tau}\;d\tau $$
I know how to get the $W(t,f)$ of $$x(t)=\cos(2\pi f_i t)$$ which is
$$[\delta(f+f_i)+\delta(f-f_i) + 2\delta(f) \cos(4 \pi f_i t)].$$
However, I found a difficulty in deriving the cosine when the phase $\theta$ is included, $x(t)=\cos(2\pi f_i t + \theta)$.