Wigner functions and density matrices are computer simulated for various quantum mechanical states of light. Wigner function and density matrices are evaluated by filtered back projection which includes inverse Radon transform from the distribution function of the photocurrents, which are calculated in the balanced homodyne detection scheme. The density matrix is also directly obtained by using the pattern function from the simulated phase independent photocurrent distribution function. ction.