A method to recover the density matrix of radiation states expressed b
y a finite superpositions of number states is presented. It starts fro
m realistic, not fully efficient, heterodyne or double homodyne detect
ion. Examples are given by means of numerical simulations.