Recursive Computation for Error Covariance Matrix in Subspace Identification Methods