A hybrid mode matching-finite element method and spectral decomposition approach for the analysis of large finite phased arrays of waveguides