A new method of designing a robust nonlinear observer is presented for a class of neural mass models by using the Lur’e system theory and the projection lemma. The observer is robust towards input uncertainty and measurement noise. It is applied to estimate the unmeasured membrane potential of neural populations from the electroencephalogram (EEG) produced by the neural mass models. An illustrative example shows the effectiveness of the proposed method.
1. Introduction
Mathematical modelling provides a powerful tool for studying mechanisms involved in the generation of different electroencephalogram (EEG) rhythms and neuronal processes of neurological disorders. There are two types of approaches to model neural signals. One is on the basis of networks built with a large number of elementary cells to describe the activity of a given system. The other is a lumped-parameter approach in which neural populations are modeled as nonlinear oscillators. Neural mass models are based on the latter approach. These models comprise macrocolumns or cortical areas and represent the mean activity of the whole population by using one or two state variables. It is seldom tractable to model EEG signals at the neuronal level due to the complexity of real neural networks. The use of neural mass models has been the preferred approach since 1970s. Neural mass models originated from the seminal work of Lopes da Silva et al. for alpha rhythm generation [1] and redesigned by Jansen and Rit to represent the generation of evoked potentials in the visual cortex [2]. The dynamical analysis [3–6] and control [7–9] of the neural mass models have been widely studied over the years. Despite the existence of these neural mass models for simulating distinct rhythms in EEG signals, neural activity is always measured through observing just a single variable such as voltage. A combination of noise in neurons and amplifiers as well as uncertainties in recording equipment leads to uncertainty of the measurement. The observation of states therefore plays significant roles in neuroscientific studies for better understanding of the human brain [10].
In general, neural mass models can be expressed as nonlinear systems of Lur'e type [11]. Observer design for nonlinear systems of Lur'e type [12–14] and the neural mass models [15] has been widely investigated over the years. We here introduce a new method of designing a robust nonlinear observer for the neural mass models. The Lur'e system theory and new tools in linear matrix inequality (LMI) method [16] are used to obtain the new reformulation. We should mention that this new reformulation takes input uncertainty and measurement noise into account. The superiority of the proposed method is demonstrated in the last section which is devoted to numerical comparisons.
Notation. The identity matrix is denoted by I. The symmetric block component of a symmetric matrix is denoted by ⋆. The vector norm is denoted by |·|. The L2 norm is denoted by ∥·∥2. The set of positive real numbers is denoted by ℝ+.
2. Problem Formulation
Let us consider a class of neural mass models that can be formulated as the following mathematical structure:
(1)x˙=Ax+Bf(Hx)+B1u,z=Cx+Dw,
where x∈ℝn is the state vector, u∈ℝp is the input, z∈ℝq is the measurement output, w∈ℝs is the measurement noise, A∈ℝn×n, B∈ℝn×m, B1∈ℝn×p, C∈ℝq×n, H∈ℝm×n, and D∈ℝq×s are constant matrices, and f(·): ℝn→ℝm is a memoryless nonlinear vector valued function which is continuously differentiable on ℝm. Each entry of the state-dependent nonlinearity f(Hx) is a function of a linear combination of the states
(2)fl=fl(∑j=1nhljxj)=fl(hlx),l=1,2,…,m,
where hl=[hl1⋯hln]. It satisfies certain slope-restricted condition
(3)0≤fl′(σ)≤δl,∀σ∈ℝ,l=1,2,…,m,
where δl≥0. The models in David and Friston [4], Goodfellow et al. [6], Jansen and Rit [2], and Wendling et al. [3] can all be expressed as the form of (1). Let us construct the following observer for plant (1):
(4)x^˙=Ax^+Bf[Hx^+K(Cx^-z)]+L(Cx^-z)+B1(u+d),
where x^ is the estimation of state, d∈ℝp is the disturbance of input, and K∈ℝm×q, L∈ℝn×q are the observer matrices to be designed. Defining the observer error as e=x^-x, the dynamics of it are governed by
(5)e˙=(A+LC)e+B1d-LDw+Bη(e,t),
where η(e,t)=f(V)-f(U), V=Hx^+K(Cx^-z), and U=Hx. Note from (3) that each entry of the nonlinearity η(e,t) satisfies
(6)0≤ηl(e,t)vl-ul=fl(vl)-fl(ul)vl-ul≤δl,∀e∈ℝ,∀t∈ℝ+,l=1,2,…,m.
The observer design for (1) consists in finding observer matrices K and L such that the observer error e satisfies the following property for all t≥0:
(7)∥e∥2≤κ|e0|+ρd∥d∥2+ρw∥w∥2,
where scalars κ>0, ρd≥0, and ρw≥0. The disturbance gains from d and w to e are ρd and ρw.
3. Main ResultsTheorem 1.
Consider plant (1) and observer (4). Under the slope restrictions (3), if there exist a matrix P=PT>0, a diagonal matrix M=diag(m1,…,mm)≥0, matrices K and L, nonsingular matrices G and F with appropriate dimensions, and scalar constants ε>0, μw≥0, and μd≥0 such that
(8)[G+GTΓ1GBGB1-GLD⋆Γ2Γ3FB1-FLD⋆⋆-M0-12ΔMKD⋆⋆⋆-μdI0⋆⋆⋆⋆-μwI]≤0,Γ1=G(A+LC)+P-FT,Γ2=F(A+LC)+(A+LC)TFT+εI,Γ3=FB+12(H+KC)TMΔ,
then the observer error e satisfies (7) for all t≥0, whereΔ=diag(δ1,…,δm), κ=λmax(P)/ε,ρd=μd/ε, and ρw=μw/ε.
Proof.
The inequality (8) can be written as
(9)Ω+V1[GF]V2+V2T[GF]TV1T≤0,
where
(10)Ω=[0P000⋆εI12(H+KC)TMΔ00⋆⋆-M0-12ΔMKD⋆⋆⋆-μdI0⋆⋆⋆⋆-μwI],V1=[I00000I000]T,V2=[-IA+LCBB1-LD].
By using the well-known projection lemma in LMI method [16], (9) can be transformed into
(11)[Π1Π2PB1-PLD⋆-M0-12ΔMKD⋆⋆-μdI0⋆⋆⋆-μwI]≤0,
where
(12)Π1=(A+LC)TP+P(A+LC)+εI,Π2=PB+12(H+KC)TMΔ.
The derivative of V(e)=eTPe is given by
(13)V˙=e˙TPe+eTPe˙≤e˙TPe+eTPe˙-∑l=1mmlηl(e,t)[ηl(e,t)-δl(vl-ul)].
Applying (11), we have
(14)V˙≤-εeTe+μddTd+μwwTw
from which it follows that
(15)∥e∥2≤λmax(P)ε|e0|+μdε∥d∥2+μwε∥w∥2.
Hence, (7) results from κ=λmax(P)/ε,ρd=μd/ε, and ρw=μw/ε.
Theorem 1 shows that the observer design for (1) consists in finding observer matrices K and L to satisfy (8) with a symmetric matrix P>0, a diagonal matrix M≥0, nonsingular matrices G, F, and scalar constants ε>0, μw≥0, and μd≥0. The feasible solution of (8) can be obtained by solving the following optimization problem:
(16)minmax{μw,μd}s.t.(8),P>0,M≥0,ε>0,μw≥0,μd≥0.
Efficient numerical tools such as YALMIP in MATLAB are available for this task. Once the values of μw and μd are computed, the disturbance gains ρw and ρd can also be derived. When no input uncertainty and measurement noise are taken into account, Theorem 1 is simplified as follows.
Theorem 2.
Consider plant (1) and observer (4) with d=0 and w=0. Under the slope restrictions (3), if there exist a matrix P=PT>0, a diagonal matrix M=diag(m1,…,mm)≥0, matrices K, L, nonsingular matrices G, F, and scalar constants ε>0, μw≥0 such that
(17)[G+GTΓ1GB⋆Γ2Γ3⋆⋆-M]≤0,
where Γ1, Γ2, and Γ3 are defined as Theorem 1, then the origin of the observer error system (5) is globally exponentially stable.
4. Simulations
Let us consider a neural mass model developed by Jansen and Rit [2]. This type of single cortical column model with altered parameters is able to generate realistic patterns such as alpha rhythms and epileptiform spikes in EEG. It can be formulated as the form of (1) with the state vector x=[x1x2x3x4x5x6x7x8]T, where xi(i=1,3,5,7) are the mean membrane postsynaptic potentials and xj(j=2,4,6,8) are their time derivatives. The input u is the afferent influence from neighbouring or more distant columns and is modeled by a Gaussian white noise with mean value 90 and standard deviation 30. The output z is the EEG measurement available to the observer. The system matrices are as follows:
(18)A=diag(A1,…,A4),Ai=[01-κi2-2κi],κ1=κ2=a,κ3=b,κ4=ad,B1=[000θaa0000]T,C=[0010-1000],D=1,B=[0θaa00000θaad000θaaC2000000000θbbC400]T,H=[0010-1000C10000000C30000000],f(Hx)=[S(x3-x5)S(C1x1)S(C3x1)]T,S(v)=2e01+er(v0-v).
The function S(·) satisfies (3) with δl=(1/2)e0r(l=1,2,3). All values of the constants in the model are set on a physiological interpretation basis which can be found in [2]. The standard values of these constants are given anatomically as
(19)θa=3.25mV,θb=22mV,a=100s-1,b=50s-1,v0=6mV,e0=2.5s-1,r=0.56mV-1,ad=33s-1,C1=135,C2=108,C3=33.75,C4=33.75.
We design the robust nonlinear observer (4) for the neural mass model. The performance of the observer obtained from Theorem 1 is presented in what follows. Input disturbance d~N(0,0.12) and measurement noise w~N(0,0.92) are introduced in the design of robust nonlinear observer. For the robust nonlinear observer, we solve the optimization problem (16) to obtain K and L. The computed disturbance gains ρw=596 and ρd=3.53 are derived by using the YALMIP toolbox in MATLAB. They are much less than the values given in [15]. In the following simulations, the initial states of the neural mass model and the observer are chosen as x(0)=[10.510.510.510.5]T and x^(0)=[00000000]T, respectively. Figure 1 presents the time evolutions of the states x1–x8 (black lines) and their estimations, that is, the states of observer (4) proposed in this study (red lines) and in [15] (blue line). Insets are given to show the zoom-in on data. Figure 1 shows that the states of observer (4) obtained from Theorem 1 do converge to a neighbourhood of the states of the neural mass model. It also shows that the observer proposed in this study performs better than that proposed in [15].
The time evolutions of the states x1–x8 and their estimations.
5. Conclusions
We have designed a robust nonlinear observer for a class of neural mass models by using the Lur’e system theory and the projection lemma. The resulting observer inhibits input uncertainty and measurement noise. We apply this observer to the neural mass model that generates alpha rhythms to estimate the mean membrane potential of neural populations from the EEG measurement. We show that the proposed observer performs better than some existing ones. The proposed method can also be applied to other types of neural models that have the typical structure of Lur’e systems.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This research is supported by the National Natural Science Foundation of China (61004050, 61172095).
Lopes da SilvaF. H.HoeksA.SmitsH.ZetterbergL. H.Model of brain rhythmic activity197415127372-s2.0-0016138266JansenB. H.RitV. G.Electroencephalogram and visual evoked potential generation in a mathematical model of coupled cortical columns19957343573662-s2.0-002937494610.1007/BF00199471WendlingF.BellangerJ. J.BartolomeiF.ChauvelP.Relevance of nonlinear lumped-parameter models in the analysis of depth-EEG epileptic signals20008343673782-s2.0-0034305571DavidO.FristonK. J.A neural mass model for MEG/EEG: coupling and neuronal dynamics2003203174317552-s2.0-034541500910.1016/j.neuroimage.2003.07.015BabajaniA.Soltanian-ZadehH.Integrated MEG/EEG and fMRI model based on neural masses2006539179418012-s2.0-3374781641510.1109/TBME.2006.873748GoodfellowM.SchindlerK.BaierG.Self-organised transients in a neural mass model of epileptogenic tissue dynamics2012593264426602-s2.0-8485547182810.1016/j.neuroimage.2011.08.060LiuX.LiuH. J.TangY. G.GaoQ.Fuzzy PID control of epileptiform spikes in a neural mass model2013711-21323LiuX.GaoQ.MaB. W.DuJ. J.RenW. J.Analysis and control of epileptiform spikes in a class of neural mass models201320131179250710.1155/2013/792507LiuX.GaoQ.Parameter estimation and control for a neural mass model based on the unscented Kalman filter201388404290510.1103/PhysRevE.88.042905SchiffS. J.2011London, UKThe MIT PressLeonovG. A.PonomarenkoD. V.SmirnovaV. B.1996SingaporeWorld ScientificArcakM.KokotovićP.Nonlinear observers: a circle criterion design and robustness analysis20013712192319302-s2.0-003554640210.1016/S0005-1098(01)00160-1FanX.ArcakM.Observer design for systems with multivariable monotone nonlinearities20035043193302-s2.0-014188711910.1016/S0167-6911(03)00170-1ZemoucheA.BoutayebM.A unified H∞ adaptive observer synthesis method for a class of systems with both Lipschitz and monotone nonlinearities20095842822882-s2.0-6044908788610.1016/j.sysconle.2008.11.007ChongM.PostoyanR.NešićD.KuhlmannL.VarsavskyA.A robust circle criterion observer with application to neural mass models2012481129862989BoydS.GhaouiL. E.FeronE.BalakrishnanV.1994Philadelphia, USASIAM