Summary: We obtain expansions of weighted pseudoinverse matrices with singular weights into matrix power products with negative exponents and arbitrary positive parameters. We show that the rate of convergence of these expansions depends on a parameter. On the basis of the proposed expansions, we construct and investigate iteration methods with quadratic rate of convergence for the calculation of weighted pseudoinverse matrices and weighted normal pseudosolutions. Iteration methods for the calculation of weighted normal pseudosolutions are adapted to the solution of least-squares problems with constraints.