1. A method of model reduction and system identification of a dynamic system with multiple inputs, comprising: generating a plurality of statistically independent random numbers for use as input signals; and performing a singular-value-decomposition directly on a system response of the dynamic system due to a simultaneous excitation of the plurality of input signals.