Joint diagonalization via subspace fitting techniques