Nonlinear analysis of biological systems using short M-sequences and sparse-stimulation techniques

The m-sequence pseudorandom signal has been shown to be a more effective probing signal than traditional Gaussian white noise for studying nonlinear biological systems using cross-correlation techniques. The effectiveness is evidenced by the high signal-to-noise (S/N) ratio and speed of data acquisition. However, the “anomalies” that occur in the estimations of the cross-correlations represent an obstacle that prevents m-sequences from being more widely used for studying nonlinear systems. The sparse-stimulation method for measuring system kernels can help alleviate estimation errors caused by anomalies. In this paper, a “padded sparse-stimulation” method is evaluated, a modification of the “inserted sparse-stimulation” technique introduced by Sutter, along with a short m-sequence as a probing signal. Computer simulations show that both the “padded” and “inserted” methods can effectively eliminate the anomalies in the calculation of the second-order kernel, even when short m-sequences were used (length of 1023 for a binary m-sequence, and 728 for a ternary m-sequence). Preliminary experimental data from neuromagnetic studies of the human visual system are also presented, demonstrating that the system kernels can be measured with high signal-to-noise (S/N) ratios using short m-sequences.

[1]  I. Ohzawa,et al.  Stereoscopic depth discrimination in the visual cortex: neurons ideally suited as disparity detectors. , 1990, Science.

[2]  J. P. Jones,et al.  The two-dimensional spatial structure of simple receptive fields in cat striate cortex. , 1987, Journal of neurophysiology.

[3]  C C Wood,et al.  Retinotopic organization of human visual cortex: departures from the classical model. , 1996, Cerebral cortex.

[4]  B. Knight,et al.  Nonlinear analysis with an arbitrary stimulus ensemble , 1979 .

[5]  Vasilis Z. Marmarelis,et al.  Applications of minimum-order Wiener modeling to retinal ganglion cell spatiotemporal dynamics , 1989 .

[6]  M. Schetzen.  The Volterra and Wiener Theories of Nonlinear Systems , 1980 .

[7]  E. Adelson,et al.  Directionally selective complex cells and the computation of motion energy in cat visual cortex , 1992, Vision Research.

[8]  V. Marmarelis.  Random versus pseudorandom test signals in nonlinear-system identification , 1978 .

[9]  D. Regan.  Human brain electrophysiology: Evoked potentials and evoked magnetic fields in science and medicine , 1989 .

[10]  M. Korenberg,et al.  Exact orthogonal kernel estimation from finite data records: Extending Wiener's identification of nonlinear systems , 1988, Annals of Biomedical Engineering.

[11]  R. Srebro.  WOULD AN EVALUATION OF BINOCULARITY USING ALGEBRAIC INTERACTION BETWEEN THE TWO EYES BE MODIFIED USING THE PSEUDORANDOM BINARY SEQUENCE METHOD? , 1980, Annals of the New York Academy of Sciences.

[12]  D. Pollen,et al.  Cross-correlation analyses of nonlinear systems with spatiotemporal inputs (visual neurons) , 1993, IEEE Transactions on Biomedical Engineering.

[13]  H. Barker,et al.  Performance of antisymmetric pseudorandom signals in the measurement of 2nd-order Volterra kernals by crosscorrelation , 1972 .

[14]  R. Srebro.  An analysis of the VEP to luminance modulation and of its nonlinearity , 1992, Vision Research.

[15]  L. Chua,et al.  Measuring Volterra kernels (II) , 1989 .

[16]  Y. W. Lee,et al.  Measurement of the Wiener Kernels of a Non-linear System by Cross-correlation† , 1965 .

[17]  Leon O. Chua,et al.  Measuring Volterra kernels , 1983 .

[18]  Lowell D. Jacobson,et al.  Structural testing of multi-input linear—nonlinear cascade models for cells in macaque striate cortex , 1993, Vision Research.

[19]  C J Aine,et al.  Temporal dynamics of visual-evoked neuromagnetic sources: effects of stimulus parameters and selective attention. , 1995, The International journal of neuroscience.

[20]  N. Wiener,et al.  Nonlinear Problems in Random Theory , 1964 .

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

[22]  Vito Volterra,et al.  Theory of Functionals and of Integral and Integro-Differential Equations , 2005 .

[23]  V. Marmarelis.  Identification of nonlinear biological systems using laguerre expansions of kernels , 1993, Annals of Biomedical Engineering.

[24]  N. Ream.  Nonlinear identification using inverse-repeatm sequences , 1970 .

[25]  R. Shapley,et al.  Spatial structure of cone inputs to receptive fields in primate lateral geniculate nucleus , 1992, Nature.

[26]  J.L. Wyatt,et al.  Nonlinear system theory. The Volterra/Weiner approach , 1982, Proceedings of the IEEE.

[27]  W. Wright,et al.  PSEUDORANDOM SEQUENCES IN THE STUDY OF EVOKED POTENTIALS * , 1980, Annals of the New York Academy of Sciences.

[28]  Vasilis Z. Marmarelis,et al.  Nonlinear analysis of renal autoregulation under broadband forcing conditions , 1993, Annals of Biomedical Engineering.

[29]  Hai-Wen Chen,et al.  Modeling and identification of parallel nonlinear systems: structural classification and parameter estimation methods , 1995, Proc. IEEE.

[30]  S. Klein,et al.  The topography of visual evoked response properties across the visual field. , 1994, Electroencephalography and clinical neurophysiology.

[31]  I. Sandberg.  The mathematical foundations of associated expansions for mildly nonlinear systems , 1983 .

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

[33]  William H. Press,et al.  Numerical recipes in C , 2002 .

[34]  Jonathan D. Victor,et al.  AN EXTENSION OF THE M-SEQUENCE TECHNIQUE FOR THE ANALYSIS OF MULTI-INPUT NONLINEAR SYSTEMS , 1994 .

[35]  Katsuhiko Ogata,et al.  Modern Control Engineering , 1970 .

[36]  W. D. T. Davies,et al.  System identification for self-adaptive control , 1970 .

[37]  Vasilis Z. Marmarelis,et al.  Practicable identification of nonstationary nonlinear systems , 1981 .

[38]  Ee Sutter,et al.  A deterministic approach to nonlinear systems analysis , 1992 .

[39]  H. A. Barker,et al.  High-order autocorelation functions of pseudorandom signals based on m sequences , 1970 .

[40]  Hiroko M. Sakai,et al.  Generation and transformation of second-order nonlinearity in catfish retina , 2006, Annals of Biomedical Engineering.

[41]  Leon O. Chua,et al.  Measuring volterra kernels III: How to estimate the highest significant order , 1991, Int. J. Circuit Theory Appl..

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

[43]  D. Pollen,et al.  Space-time spectra of complex cell filters in the macaque monkey: A comparison of results obtained with pseudowhite noise and grating stimuli , 1994, Visual Neuroscience.