A new algorithm for the identification of multiple input Wiener systems

Multiple-input Wiener systems consist of two or more linear dynamic elements, whose outputs are transformed by a multiple-input static non-linearity. Korenberg (1985) demonstrated that the linear elements of these systems can be estimated using either a first order input-ouput cross-covariance or a slice of the second, or higher, order input-output cross-covariance function. Korenberg's work used a multiple input LNL structure, in which the output of the static nonlinearity was then filtered by a linear dynamic system. In this paper we show that by restricting our study to the slightly simpler Wiener structure, it is possible to improve the linear subsystem estimates obtained from the measured cross-covariance functions. Three algorithms, which taken together can identify any multiple-input Wiener system, have been developed. We present the theory underlying these algorithms and detail their implementation. Simulation results are then presented which demonstrate that the algorithms are robust in the presence of output noise, and provide good estimates of the system dynamics under a wide set of conditions.

[1]  M. Korenberg Identifying Noisy Cascades of Linear and Static Nonlinear Systems , 1985 .

[2]  I. Hunter,et al.  Dynamics of human ankle stiffness: variation with mean ankle torque. , 1982, Journal of biomechanics.

[3]  M. J. Korenberg,et al.  The identification of nonlinear biological systems: LNL cascade models , 1986, Biological Cybernetics.

[4]  K. Naka,et al.  Nonlinear analysis and synthesis of receptive-field responses in the catfish retina. 3. Two-input white-noise analysis. , 1973, Journal of neurophysiology.

[5]  Michael J. Korenberg,et al.  Parallel cascade identification and kernel estimation for nonlinear systems , 2006, Annals of Biomedical Engineering.

[6]  W. Rugh Nonlinear System Theory: The Volterra / Wiener Approach , 1981 .

[7]  P Z Marmarelis,et al.  Nonlinear analysis and synthesis of receptive-field responses in the catfish retina. II. One-input white-noise analysis. , 1973, Journal of neurophysiology.

[8]  J. Bendat,et al.  Random Data: Analysis and Measurement Procedures , 1987 .

[9]  J. Gaska,et al.  Structural classification of multi-input nonlinear systems , 1990, Biological Cybernetics.

[10]  I. Hunter,et al.  Nonlinear identification of stretch reflex dynamics , 2006, Annals of Biomedical Engineering.

[11]  Gene H. Golub,et al.  Matrix computations , 1983 .

[12]  I. Hunter,et al.  Two-sided linear filter identification , 1983, Medical and Biological Engineering and Computing.

[13]  Vasilis Z. Marmarelis,et al.  Analysis of Physiological Systems , 1978, Computers in Biology and Medicine.

[14]  K. Naka,et al.  Identification of multi-input biological systems. , 1974, IEEE transactions on bio-medical engineering.

[15]  P Z Marmarelis,et al.  Nonlinear analysis and synthesis of receptive-field responses in the catfish retina. I. Horizontal cell leads to ganglion cell chain. , 1973, Journal of neurophysiology.

[16]  M. J. Korenberg,et al.  The identification of nonlinear biological systems: Wiener and Hammerstein cascade models , 1986, Biological Cybernetics.

[17]  Syozo Yasui,et al.  Spatio-Temporal Receptive Field Measurement of Retinal Neurons by Random Pattern Stimulation and Cross Correlation , 1979, IEEE Transactions on Biomedical Engineering.

[18]  I. Hunter,et al.  Dynamics of human ankle stiffness: variation with displacement amplitude. , 1982, Journal of biomechanics.

[19]  J. Bendat,et al.  Random Data: Analysis and Measurement Procedures , 1971 .