The aim is to build a model for the S-parameter matrix of a linear, time invariant multiple-port system operating at RF frequencies. The proposed method delivers stable, accurate models starting from measurements as well as from simulation data. The use of orthogonal polynomials allows high model orders.