Optimal Estimation of Shared Singular Subspaces across Multiple Noisy Matrices