A frequency domain technique for characterizing nonlinearities in biological systems*

https://doi.org/10.1016/S0022-5193(88)80323-0Get rights and content

Response asymmetry (e.g. to light ON vs. OFF), a frequently encountered property of neurons at both peripheral and central levels in sensory pathways, can be modelled as a rectifier. When fed with the sum of two sinusoids of frequencies F1 and F2 Hz, the output of a rectifier contains multiple discrete frequency components of frequency (nF1 ± mF2) where n and m are zero or integers. We describe a double Fourier series method for obtaining the amplitudes and phases of these components for physiologically relevant neural models including (1) compressive, linear, accelerating and mixed compressive/accelerating single model neurons, (2) a cascaded series of model neurons and (3) the parallel/cascaded case corresponding to dichoptic or dichotic stimulation. If one of the inputs is held constant while the other's amplitude is varied, we obtain a family of curves—one for each (nF1 ± mF2) component. The family of curves seems to be characteristic of the particular non-linear system.

Neural models may be tested by first calculating the family of curves and then comparing these theoretical data with a physiologically measured family of curves. Illustrating this approach, we compare the predictions of case (1) and case (2) models with the family of curves obtained by stimulating with F1 Hz flicker of fixed modulation depth superimposed on F2 Hz flicker of variable modulation depth. Evoked brain responses were analysed by non-destructive zoom-FFT giving resolution up to the Heisenberg-Gabor limit of ΔF = ΔT−1, where ΔT is the recording duration. Thus, for example, the spectrum could contain 50 000 lines over a D.C.-100 Hz bandwidth for a 500 second recording duration. Empirically, the bandwidth of discrete physiological components is no more than 0·004 Hz.

Like the Wiener kernel time-domain approach, this frequency-domain approach is restricted to time-invariant systems with a limited settling time and without appreciable hysteresis. On the other hand, the Weiner kernel approach seems to handle dynamic systems more conveniently. But this frequency-domain method has the important advantage that different signal components are well separated from each other and from noise, whereas signal and noise overlap in the time domain. Consequently, high-order terms are more easily recorded using the frequency domain method, thus giving sharper distinctions between candidate non-linear models.

References (53)

  • De ValoisR.L. et al.

    Vision Res.

    (1982)
  • ReganD.

    Vision Res.

    (1983)
  • ReganD. et al.

    Vision Res.

    (1973)
  • ReganD. et al.

    Vision. Res.

    (1987)
  • ReganD. et al.

    Vision Res.

    (1988)
  • ReganD. et al.

    Brain Res.

    (1988)
  • SovijarviA.R.A. et al.

    Brain Res.

    (1974)
  • VictorJ.D. et al.

    Biophys. J.

    (1980)
  • WilliamsonS. et al.

    J. Mag. Mag. Materials

    (1981)
  • BarlowH.B. et al.

    J. Physiol.

    (1965)
  • BarlowH.B. et al.

    J. Physiol.

    (1967)
  • BennettW.R.

    Bell System Tech. J.

    (1933)
  • BeverleyK.I. et al.

    J. Physiol.

    (1973)
  • BruggeJ.F. et al.

    Auditory cortex

  • ClynesM. et al.

    Ann. N. Y. Acad. Sci.

    (1964)
  • CynaderM. et al.

    J. Physiol.

    (1978)
  • EvansE.F. et al.

    J. Physiol., Lond.

    (1964)
  • FiorentiA. et al.

    Vision Res.

    (1983)
  • FosterK.H. et al.

    J. Physiol.

    (1985)
  • GaborD.

    J. Inst. elect. Eng., Lond.

    (1946)
  • GardnerR.B. et al.

    J. acoust. Soc. Am.

    (1979)
  • HobsonE.W.

    On the Theory of Functions of a Real Variable and the Theory of Fourier's Series

    (1926)
  • HubelD.H. et al.

    J. Physiol.

    (1962)
  • KingL.V.

    On the Direct Numerical Calculation of Elliptic Functions and Integrals

    (1924)
  • LeggeG.E. et al.

    J. opt. Soc. Am.

    (1980)
  • MarmarelisP.Z. et al.

    Analysis of Physiological Systems. The White Noise Approach

    (1978)
  • Cited by (0)

    *

    Parts of this research were supported by the Canadian NSERG and MRC and the NEI, and sponsored by the U.S. Air Force Office of Scientific Research.

    View full text