Analysis of stability and convergence for L-type formulas combined with a spatial finite element method for solving subdiffusion problems