This is something like an ‘appendix’ to two blogposts I wrote, going into some technical details. It won’t make a lot of sense without that context - read them first!
I want to show that the procedure I outlined in those posts for calculating the Wigner function of a qubit matches the normal calculation method described in e.g. Wootters 1987.
In the blog posts, I write the Wigner function out as follows:
where the are the answers to the three ‘questions’ - the expectation values of the measurements
where the are the Pauli matrices. The square-bracketed matrices are supposed to represent the 2x2 phase space graphs I used in the blog posts. There’s an annoying notational clash where matrices normally have their indices labelled starting in the top left, whereas graphs have their origin in the bottom left. I’ve kept with the graph convention (as did Wootters), so that, for example would be the bottom left component, and used square brackets rather than round brackets to visually distinguish them from normal matrices.
Simplifying, we get
Now to make contact with the standard version. Wootters calculates the four components of the Wigner function from the density matrix of the qubit as follows:
where the are the phase point operators
Written out in components these become
Any density matrix can be decomposed in the following way (see e.g. here):
where the s are components of a vector with . Plugging this decomposition into the standard formula  for the Wigner function gives
Finally we need to relate the s back to the s. We have , which decomposes as
Putting these back into  gives the formula  that I used for the Wigner function in my blog posts.