Iterative alternating least-aquares approach to jointly estimate the RETFs and the diffuse PSD