Given a list of state-space models belonging to the StateSpace class, creates a joined system whose output is the sum of the state-space outputs. If wv is not None, this is a list of weights, such that the output is:
y = sum( wv[ii] y_ii )
Ref: equation (4.22) of Benner, P., Gugercin, S. & Willcox, K., 2015. A Survey of Projection-Based Model Reduction Methods for Parametric Dynamical Systems. SIAM Review, 57(4), pp.483–531.
system matrices must be numpy arrays
the function does not perform any check!