Search

This paper considers the use of cubic splines, instead of polynomials, to represent the static nonlinearities in block structured models. It introduces a system identification algorithm for the Hammerstein structure, a static nonlinearity followed by a linear filter, where cubic splines represent the static nonlinearity and the linear dynamics are modeled using a finite impulse response filter. The algorithm uses a separable least squares Levenberg-Marquardt optimization to identify Hammerstein cascades whose nonlinearities are modeled by either cubic splines or polynomials. These algorithms are compared in simulation, where the effects of variations in the input spectrum and distribution, and those of the measurement noise are examined. The two algorithms are used to fit Hammerstein models to stretch reflex electromyogram (EMG) data recorded from a spinal cord injured patient. The model with the cubic spline nonlinearity provides more accurate predictions of the reflex EMG than the polynomial based model, even in novel data. [ABSTRACT FROM AUTHOR]

The Hammerstein cascade, consisting of a zero-memory nonlinearity followed by a linear filter, is often used to model nonlinear biological systems. This structure can represent some high-order nonlinear systems accurately with relatively few parameters. However, it is not possible, in general, to estimate the parameters of a Hammerstein cascade in closed form. The most effective method available to date uses an iterative approach, which alternates between estimating the linear element from a crosscorrelation, and then fitting a polynomial to the nonlinearity via linear regression. This paper proposes the use of separable least squares optimization methods to estimate the linear and nonlinear elements simultaneously in a least squares framework. A separable least squares algorithm for the identification of Hammerstein cascades is developed and used to analyze stretch reflex electromyogram data from two experimental subjects. The results show that in each case the proposed algorithm produced a better model, in that it predicted the system’s response to novel inputs more accurately, than did models estimated using the traditional iterative algorithm. Monte-Carlo simulations demonstrated that when the input is a non-Gaussian, nonwhite signal, as is often the case experimentally, the traditional iterative identification approach produces biased models, whereas the separable least squares approach proposed in this paper does not. © 2001 Biomedical Engineering Society. PAC01: 8719Ff, 8719Nn, 8710+e [ABSTRACT FROM AUTHOR]