The paper describes a method that uses a multivariate model in the form of a nonhomogeneous finite difference equation to identify modal parameters of a mechanical structure. The modal parameters of this equation are estimable using a model that involves multiple input, multiple output vibration data. Thus, improved global estimates of modal parameters can be obtained, including the identification of highly coupled and pseudo-repeated modes of vibration. When the data are in the form of impulse or free decay responses, then the parameters of the homogeneous part of the equation can be estimated separately, and the method is then related to the Least Squares Complex Exponential method, the Polyreference Time Domain method and the Ibrahim Time Domain method.