Stability Preservation for Parametric Model Order Reduction by Matrix Interpolation