In this paper, a parametric model order reduction approach based on scaled Laguerre expansion series is investigated in frequency domain for large‐scale parametric systems. The proposed approach is implemented by two algorithms: the one‐sided algorithm and the two‐sided algorithm. Based on modified multi‐Arnoldi process for construction of projection matrices, reduced parametric systems with high fidelity are achieved in the full range of parametric domain. The reduced parametric systems obtained by the proposed algorithms preserve not only the parametric dependence of the original parametric system but also the first several Laguerre coefficients. The proposed algorithms are verified by two benchmarks (a synthetic parametric system and a microthruster unit). Comparisons of the results of original and reduced systems show that the obtained reduced systems have high fidelity.