The General Realization Algorithm (GRA) is developed to identify modal parameters of linear multi-degree-of-freedom dynamic systems subjected to measured (known) arbitrary dynamic loading from known initial conditions. The GRA extends the well known Eigensystem Realization Algorithm (ERA) based on Hankel matrix decomposition by allowing an arbitrary input signal in the realization algorithm. This generalization is obtained by performing a weighted Hankel matrix decomposition, where the weighting is determined by the loading. The state-space matrices are identified in a two-step procedure that includes a state reconstruction followed by a least squares optimization to get the minimum prediction error for the response. The statistical properties (i.e., bias, variance, and robustness to added output noise introduced to model measurement noise and modeling errors) of the modal parameter estimators provided by the GRA are investigated through numerical simulation based on a benchmark problem with non-classical damping.