Estimation of multiple large covariance matrices by a partially common diagonal and low-rank matrix decomposition