Joint Sensitivity Estimation and Image Reconstruction in Parallel Imaging Using Pre-learned Subspaces of Coil Sensitivity Functions