Читать книгу EEG Signal Processing and Machine Learning - Saeid Sanei - Страница 42

Neuronal synchronization processes, measured with brain imaging data, can be described using weakly coupled oscillator (WCO) models. Dynamic causal modelling (DCM) is used to fit the WCOs to brain imaging data and so make inferences about the structure of neuronal interactions [2]. The complex behaviours are mediated by the interaction of particular brain regions. Recent studies agreed that such interactions may be instantiated by the transient synchronization of oscillatory neuronal ensembles [3]. For example, contour detection is accompanied by gamma band synchronization in distant parts of visual cortex, multimodal object processing by parieto‐temporal synchronization in the beta band [4] and spatial memory processes by hippocampal–prefrontal synchronization in the theta band [2, 5]. DCM allows for different model structures to be compared using Bayesian model selection [6]. In [2] DCM has been extended to the study of phase coupling. One direction is based on the WCO models in which the rate of change of phase of one oscillator is related to the phase differences between itself and other oscillators [7].

The WCO theory applies to system dynamics close to limit cycles. By assuming that weak coupling leads to only small perturbations away from these cycles, one can reduce a high‐dimensional system of differential equations to one based solely on the phases of the oscillators, and pairwise interactions between them [2].

Dynamics on the limit cycle are given by [2]:

Figure 3.2 A pair of oscillators weakly coupled via the perturbation function p(φ ₁, φ ₂).

(3.3)

Then, for any perturbation of phase p(ϕ) this changes to [2]:

(3.4)

(3.5)

where:

(3.6)

For a pair of oscillators of phases φ ₁ and φ ₂, using the same analysis (see Figure 3.2, assuming p ₁₂ = p ₂₁):

(3.7)

(3.8)

If it is further assumed that the phase difference φ ₂ − φ ₁ = φ changes slowly, then [2]:

(3.9)

(3.10)

(3.11)

Γ_ij (ϕ) is called phase interaction function (PIF). Similarly, for N_R regions the rate of change of phase of the ith oscillator is given by:

(3.12)

where f_i is the intrinsic frequency of the ith oscillator. In these formulations there are two key assumptions: the first one is that the perturbations are sufficiently small that the differentiations can equivalently be evaluated at X ₀ rather than X. The second assumption is that the relative changes in the oscillator phase are sufficiently slow with respect to the oscillation frequency, that the phase offset term can be replaced by a time average.

In [2] an extension of the above WCO model is used to describe the dynamic phase changes in a network of oscillators. The use of Bayesian model comparison allows one to infer the mechanisms underlying synchronization processes in the brain. This has been applied to synthetic bimanual finger movement data from physiological models and on magnetoencephalogram (MEG) data from a study of visual working memory. The WCO approach accommodates signal nonstationarity by using differential equations which describe how the changes in phase are driven by pairwise differences in instantaneous phase.