Elsevier

Automatica

Volume 50, Issue 4, April 2014, Pages 1087-1099
Automatica

Observability limits for networked oscillators

https://doi.org/10.1016/j.automatica.2014.02.035Get rights and content

Abstract

Inspired by the neuro-scientific problem of predicting brain dynamics from electroencephalography (EEG) measurements of the brain’s electrical activity, this paper presents limitations on the observability of networked oscillators sensed with quantised measurements. The problem of predicting highly complex brain dynamics sensed with relatively limited amounts of measurement is abstracted to a study of observability in a network of oscillators. It is argued that a low-dimensional quantised measurement is in fact, by itself, an exceptionally poor observer for a large-scale oscillator network, even for the case of a completely connected graph. The main rational is based on (i) an information-theoretic argument based on ideas of entropy in measure preserving maps, (ii) a linear deterministic observability argument, and (iii) a linear stochastic approach using Kalman filtering. For prediction of brain network activity, the findings indicate that the classic EEG signal is just not precise enough to be able to provide reliable prediction and tracking in a clinical setting in view of the complexity of underlying neural dynamics.

Introduction

Observability issues in networked oscillators have implications for tracking and control in a wide range of distributed systems that exhibit self-organising behaviour through synchronisation, from modern telecommunications (Prehofer & Bettstetter, 2005) to future power systems (Butler, 2007, Rohden et al., 2012). Coupled oscillatory networks are particularly prevalent in biological systems where there is increasing interest in tracking and predicting dynamics for applications in medical bionics.

Our particular motivational interest is the human brain, an oscillatory system with an estimated 86G neuron cells (Azevedo et al., 2009) networked with 1P synaptic connections.2 EEG recordings of the brain’s electrical activity typically provide our output signal from which to observe the underlying activity. Observability of the activity of neurons in the brain system from EEG measurements is essential for advancement of medical treatment through systems control in a range of neural conditions including epilepsy, Parkinson’s disease and depression. Strong observability is also necessary for prediction in neurology, for example of epileptic episodes or to determine which patient are likely to respond well to classes of pharmaceutical drugs.

Current network observability analyses typically seek to exploit redundancies in the connection pathways between network nodes to determine the minimum number of sensor measurements required such that all nodes are either directly or indirectly reachable. Early work includes finding observable islands (or sub-networks) when the network as a whole is unobservable and further determining where to place additional measurements to reach network areas beyond these islands (Monticelli and Wu, 1985, Wu and Monticelli, 1985). This early work is iterative in nature and therefore computationally intractable for large-scale networks. Recent work by Liu et al. in network observability, using graphical approaches to find the minimum sensor set, does however cater for large-scale networks (Liu, Slotine, & Barabási, 2013). These methods all rely on exploiting topological clustering. By contrast, the approach in this paper considers a network with fully-connected graph, where clustering is unavailable.

Liu et al. also highlight interest in the applicability of their methods to coupled oscillator systems (Liu et al., 2013); a similar question was recently considered, in the small network case, for synchronous neuronally-inspired networks (Whalen, Brennan, Sauer, & Schiff, 2012). In this case observability was found to be quite limited and heavily influenced by topology and symmetry. In this paper we consider observability in large scale coupled oscillator systems.

Liu and Bitmead consider network observability in non-linear stochastic networks, defining observability in information theoretic terms by comparing the entropy of the state with the conditional entropy of the state given the measurement (Liu & Bitmead, 2011). Although stochastic systems are not considered in this paper, similar ideas of linking observability to entropy are applied here (although it is entropy of the state map that is considered in this work rather than a condition related to the mutual information between state and measurement signal). In both cases, however, positive entropy effectively allows for a reduction of state uncertainty with an increasing sequence of measurements. As articulated by Liu and Bitmead, the power in an entropy definition of observability is that it can apply to both linear and non-linear systems.

The Takens–Aeyels embedding theorem (Aeyels, 1981, Takens, 1981) states that observability is a generic property in non-linear (autonomous) systems. Moreover, a state may be reconstructed from the measurement vectors for sufficiently large measurement time series. This state can then be used to infer dynamics. This remarkable theorem of delay reconstruction provides an elegant and resourceful tool, but, it is limited to autonomous, stationary, noise-free systems (Kantz & Schreiber, 2004, Chapter 3). Despite such limitations, delay reconstruction is widely applied to real-world systems, for example to the brain (Iasemidis et al., 1990, Lehnertz and Elger, 1998, Le Van Quyen et al., 2001) where dynamics are certainly not noise-free or indeed autonomous and can only be considered quasi-stationary on short time scales (10s) (Niedermeyer & Lopes Da Silva, 2005). How realistic is the Takens–Aeyels embedding theorem here, particularly in light of such a large-scale system as the brain?

In this paper the question of what one can observe from an EEG record is reduced to the generalised question of how quantisation of an output, from a large-scale system, effects the observability. This effect of quantisation has implications beyond the reconstruction of brain dynamics for the observability of any practical system.

A “synthetic” brain-like situation is presented in Section  2, fully under our control, within which the limits of a brain-like EEG recording can be investigated. While Section  3 demonstrates that we have a theoretically observable system, practical considerations reveal a severe lack of observability using arguments from (i) information-theoretic ideas of entropy in measure preserving maps in Section  4.1, (ii) linear deterministic observability in Section  4.2 and (iii) linear stochastic Kalman filtering in Section  4.3. The implications for prediction and tracking of brain dynamics are discussed in Section  5 followed by concluding comments in Section  6.

Section snippets

Abstraction to networked clocks

A generic and scalable coupled oscillator model (with origins in the 1985 work of Wright, Kydd, & Lees, 1985) is proposed as a synthetic brain. It is important to highlight that this is not a model which can tell us anything about the nature of brain function, but this approach is suited, however, to the specific task of investigating what underlying information one may expect to recover from the EEG signal.

The model abstracts the problem to the study of a simple network of second order

Observable in theory

Proposition 3.1 Networked Clocks with an EEG-Like Measurement are Theoretically Observable

The pair (C,A) of form given in Eqs.   (8), (3)   is observable under the condition that the clocks oscillate at distinct frequencies (i.e. ωiωj,ij ).

Proposition 3.1 can be proven using the PBH test (see footnote 4) on the basis of Lemma 3.1 combined with Lemma 3.2.

Lemma 3.1

An uncoupled clock network, (C,A) of form given in Eqs.   (8), (3)   with αij=0i,j, can be show to be observable by the PBH test.

A detailed proof of Lemma 3.1 is contained in Appendix A showing that an uncoupled clock network

Poor observability in practice

Although Section  3 shows that our system is indeed observable in the normal sense, it turns out that the simple limitation of a finite precision measurement, that is the quantisation of the EEG measurement, reduces the ability to distinguish all states drastically. Evidence for practical non-observability are presented under three approaches, an information theoretic argument, a linear deterministic argument and a stochastic argument.

Discussion

The brain system is approximately oscillatory with information communicated across neural networks by oscillatory electric field potentials. Any brain model (and indeed all existing models) attempting to characterise the link between brain tissue and EEG needs to represent these oscillatory dynamics. The coupled oscillator model used in this paper is very suitable for the express purpose of investigating observability properties. Its suitability is based on the reasoning that if the most basic

Conclusion

Networks of oscillators, despite theoretical observability and strongly connected (indeed complete) graphs, are practically unobservable given realistic constraints, such as quantised, noisy measurements and non-stationarity.

So what can be done to generate sufficiently good prediction and tracking in large-scale oscillatory networks, particular for therapeutic medical applications in neural networks? More electrodes? It is not currently feasible to dramatically increase electrode numbers in an

Acknowledgements

This work is based on ongoing collaboration with St. Vincent’s Hospital, Melbourne. We want to acknowledge in particular Prof. Mark Cook, Dr. Karen Fuller, Dr. Andrea Varsavsky, Dr. Dean Freestone, Andre Peterson, Michelle Chong and A/Prof. David Grayden. Many thanks to Prof Jonathan Manton for sifting out useful trigonometric identities from papers past. We thank National ICT Australia (NICTA) and the Victorian Life Sciences Computation Initiative (VLSCI), for financial and computational

Elma O’Sullivan-Greene received the B.E. in Electrical & Electronic Engineering from University College Cork, Ireland in 2005 and the Ph.D. degree in Biomedical Engineering from the University of Melbourne, Australia, in 2011. Currently, she is an Early Career Academic Fellow in the Department of Electrical & Electronic Engineering, The University of Melbourne, Australia. Her research interests are in biological systems and control, particularly related to neuro-engineering problems including

References (32)

  • Paul R. Halmos

    Lectures on ergodic theory

    (1956)
  • Paul Richard Halmos

    Entropy in ergodic theory

    (1959)
  • L.D. Iasemidis et al.

    Phase space topography and the Lyapunov exponent of electrocorticograms in partial seizures

    Brain Topography

    (1990)
  • Rudolph Emil Kalman

    A new approach to linear filtering and prediction problems

    Transactions of the ASME—Journal of Basic Engineering

    (1960)
  • Cited by (10)

    • Revisiting Seizure Prediction with Path Signatures: the Limitations of System Identification

      2024, 2024 Australian and New Zealand Control Conference, ANZCC 2024
    View all citing articles on Scopus

    Elma O’Sullivan-Greene received the B.E. in Electrical & Electronic Engineering from University College Cork, Ireland in 2005 and the Ph.D. degree in Biomedical Engineering from the University of Melbourne, Australia, in 2011. Currently, she is an Early Career Academic Fellow in the Department of Electrical & Electronic Engineering, The University of Melbourne, Australia. Her research interests are in biological systems and control, particularly related to neuro-engineering problems including epilepsy and vision.

    Iven Mareels was born in Aalst, Belgium, in 1959. He received the (ir) degree in electromechanical engineering from Gent University, Ghent, Belgium, in 1982 and the Ph.D. degree in systems engineering from the Australian National University, Canberra, Australia, in 1987. He is a Professor with the Department of Electrical and Electronic Engineering, University of Melbourne, Melbourne, Australia, where he is presently the Dean of the School of Engineering. His research interest focuses on the modelling and control of large scale systems.

    Levin Kuhlmann received his B.Sc. degree with honours in Neuroscience from the University of Melbourne in 2000. He obtained his Ph.D. degree in Cognitive and Neural Systems from Boston University in 2007. He has been a research fellow in the Department of Electrical Engineering and the Centre for Neural Engineering at the University of Melbourne since 2007. His research areas include signal processing and control theory applications to neural engineering, computational neuroscience, anaesthesia, epilepsy and vision.

    Anthony Burkitt holds the Chair in Bio-Signals and Bio-Systems in the Department of Electrical and Electronic Engineering at the University of Melbourne and is the Director of Bionic Vision Australia. He completed his undergraduate studies in 1980 in Physics at the Australian National University and his Ph.D. in 1983 in Theoretical Physics at Edinburgh University. His research is in the area of medical bionics, including cochlear-implant speech processing, neuro-engineering, computational neuroscience and bio-signal processing for epilepsy.

    The material in this paper was partially presented at the 48th IEEE Control and Decision Conference (CDC) and 28th Chinese Control Conference (CCC), December 16–18, 2009 in Shanghai, China and submitted to the 19th IFAC World Congress, August 24–29, 2014 Cape Town, South Africa. This paper was recommended for publication in revised form by Associate Editor Peng Shi under the direction of Editor Toshiharu Sugie.

    1

    Tel.: +61 4 06478403; fax: +61 383447707.

    View full text