Nonlinear Dyn (2014) 76:1291–1299 DOI 10.1007/s11071-013-1210-3
ORIGINAL PAPER
Fuzzy adaptive unscented Kalman filter control of epileptiform spikes in a class of neural mass models Xian Liu · Hui-Jun Liu · Ying-Gan Tang · Qing Gao · Zhan-Ming Chen
Received: 9 December 2012 / Accepted: 17 December 2013 / Published online: 12 January 2014 © Springer Science+Business Media Dordrecht 2014
Abstract A new closed-loop control method based on the fuzzy adaptive unscented Kalman filter (FAUKF) is proposed to suppress epileptiform spikes in a class of neural mass models with uncertain measurement noise. The FAUKF is used to estimate the nonlinear system states of the underlying models and amend measurement noise adaptively. The control law is constructed via the estimated states. Numerical simulations illustrate the efficiency of the proposed method. Keywords Neural mass model · Epileptiform spikes · The FAUKF control
1 Introduction The neural mass models play important roles in modeling of electroencephalography (EEG), most notably in dynamic modeling of epilepsy seizures [1,2]. These models describe the activity of neurons at the macroX. Liu · H.-J. Liu · Y.-G. Tang · Q. Gao Key Lab of Industrial Computer Control Engineering of Hebei Province, Institute of Electrical Engineering, Yanshan University, Qinhuangdao 066004, China X. Liu e-mail:
[email protected] Z.-M. Chen (B) School of Economics, Renmin University of China, Beijing 100872, China e-mail:
[email protected]
scopic level, i.e., the activity of populations of neurons. In comparison with the models aiming at the neuronal levels, such as the Hodgkin–Huxley model [3], the neural mass models are empirical priors to emulate realistic EEG because the categories of neurons are manifold and the connections of neurons are complex. It is demonstrated that the class of models can simulate the dynamics of real EEG measured during the transition from interictal to fast ictal activity of epilepsy seizures in the hippocampus. An epileptic neural mass exhibits the features of excitability which is reflected in its use to produce interictal spiking activities in a single population in response to random input [1]. The spiking activity can be propagated between populations in networks of two- and three-coupled neural populations, and the ever-increasing coupling strength may increase the rhythm of synchronized activities for multiple neural populations [1]. A new method for a parameter-driven transition to epilepsy seizure based on the neural mass models is provided by Goodfellow et al. [4] to produce slow spike–wave discharges as well as sinusoidal background oscillations. The effects of local functional heterogeneities, which are considered crucial for the mechanism of epilepsy, are investigated in a spatially extended neural mass formulation with local excitatory and inhibitory circuits [5]. The dynamic evolution of EEG activity during epileptic seizures is characterized as a path through the parameter space of a neural mass model, reflecting gradual changes in underlying physiological mechanisms [6].
123
1292
These investigations provide profound insight into understanding the mechanism of epilepsy seizures. However, the model-based closed-loop control for epilepsy seizures is also an important research subject. It can provide a theoretical basis for clinical electrotherapy of seizures, which still depends heavily on the empirical tuning of parameters and protocols until now [7–9]. The closed-loop control has been shown to be able to optimize the whole course of seizures treatment including efficacy of treatment, minimization of side effects, improvement of response by providing intermittent or minimal stimulation, reduction of damage, and minimization of power consumption [10]. A fuzzy PID controller is designed for a class of neural mass models to quench epileptiform spikes and make the output track an expected one [11]. In the closed-loop control strategy, the development of feedback algorithms requires precise metric of system state [12]. Recently, different types of state observers have been designed for different nonlinear systems, such as descriptor observer [13], acceleration observer [14], and adaptive sliding mode observer [15]. A robust circle criterion observer is designed and applied to the neural mass models [16]. The unscented Kalman filtering (UKF) is another powerful tool for the state estimate of nonlinear systems [17,18]. It is applied to observe state and track dynamics from a spatiotemporal computational model of cortical dynamics [19]. It is also applied to the parameter estimation and control of epileptiform spikes for the neural mass model [20,21]. The UKF is able to achieve good performance if the complete information of measurement noise distribution is assumed to be known. This assumption is inconsistent with practical applications because the environment changes all the time and statistics of noise cannot be accurately obtained. The UKF combining with the fuzzy inference system (FIS) is developed to amend measurement noise adaptively [22–24]. The measurement of EEG in clinical practice is quite easily to be contaminated by many factors such as high-frequency interference and power interference. The measurement noise is of uncertainty and difficult to be measured. In this study, a fuzzy adaptive UKF (FAUKF) control based on the UKF and the FIS is proposed to suppress epileptiform spikes in a class of neural mass models with uncertain measurement noise. The FAUKF estimates the system states through the UKF and adjusts measurement noise adaptively through the FIS. Thus, it enhances the filtering accu-
123
X. Liu et al.
racy in comparison with the UKF. The efficiency of the FAUKF control is demonstrated through plentiful simulation tests. 2 Model description As mentioned in the introduction, our results apply to a class of neural mass models put forward by Jansen and Rit [25]. In this model, each neuron population is composed of the pyramidal cells that receive excitatory and inhibitory feedbacks from inter-neurons (other pyramidal cells, stellate, or basket cells) that only receive excitatory input. The delay from a connection population is modeled by a similar postsynaptic potential block. This type of model can generate the epileptiform spikes in EEGs [1]. It can be described by the following differential equations: ⎧ x˙1l (t) ⎪ ⎪ ⎪ ⎪ x˙2l (t) ⎪ ⎪ ⎪ ⎪ x˙3l (t) ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎨ x˙4l (t) ⎪ ⎪ ⎪ ⎪ ⎪ x˙5l (t) ⎪ ⎪ ⎪ ⎪ x˙6l (t) ⎪ ⎪ ⎪ ⎪ x˙7l (t) ⎪ ⎪ ⎩ x˙8l (t)
= x2l (t) = AaS(x3l (t) − x5l (t)) − 2ax2l (t) − a 2 x1l (t) = x4l (t) j=l = Aa[ pl (t) + C2 S(C1 x1l (t)) + j=1,2,h¯ , N
= = = =
kl j x7l (t)] − 2ax4l (t) − a 2 x3l (t) x6l (t) BbC4 S(C3 x1l (t)) − 2bx6l (t) − b2 x5l (t) x8l (t) Aad S(x3l (t)−x5l (t)) − 2ad x8l (t) − ad2 x7l (t),
(1) where N is an integer greater than 1; l is the population under consideration; the state variables xkl (t), k = 1, 3, 5, 7 are the mean membrane potentials of the neuronal populations, while xkl (t), k = 2, 4, 6, 8 are their derivatives; the excitatory input pl (t) is the afferent influence from other populations, which is assumed to be a uniformly distributed signal; and the connectivity constants C1 , C2 , C3 , and C4 are used to model interactions between main cells and interneurons. The output of population l is y(l) (t) = x3l (t) − x5l (t), which is used to simulate EEG signals. All parameters in (1) are set on a physiological interpretation basis. The standard value of the parameters is given anatomically as [25] A = 3.25 mv, B = 22 mv, a = 100 s−1 , b = 50 s−1 , v0 = 6 mv, e0 = 2.5 s−1 , r = 0.56 mv−1 , ad = 33 s−1 , (2) C1 = 135, C2 = 108, C3 = 33.75, C4 = 33.75.
In our simulations, N is set at three, and (1) is solved by using a fourth-order Runge–Kutta differ-
Fuzzy adaptive unscented Kalman filter control
1293
Fig. 1 The model of three coupled neural populations
k13 1
3 Fuzzy adaptive unscented Kalman filter control
2 k21
3
The discrete-time filtering form of the neural mass model is described as xk = F(xk−1 ) (3) yk = Hk xk + Vk ,
k32
ential solver. The concrete connection among threecoupled neural populations is shown in Fig. 1. If the parameters are set at standard values and the coupling constants k21 , k32 , k13 are set at 0, then the simulated signals (Fig. 2a) resemble the spontaneous EEG recorded from neocortical structure electrodes during interictal periods [6]. Thus, they reflect normal activities resembling those reflected by real SEEG signals. Increasing the value of A of population 1–3.4 mv and keeping the coupling constants invariant, the first population produces sporadic spikes, and the two others exhibit normal activities, as shown in Fig. 2b. Keeping A of population 1 invariant and increasing the coupling constants k21 , k32 , k13 to 100, sustained discharges of spikes that resemble real SEEG signals during the propagation of temporal lobe seizures are observed (see Fig. 2c).
1 2n
X˜ j ,
y˜k =
j=1
15
20
4 2
5
10
15
−10 0 6
2
5
10 t(s)
15
20
y(1)(t)(mv) 5
10
15
20
4 2 0 0 6
20
4
0
5
10
15
20
4 2 0 0
Y˜ j
(5) (6)
j=1
5
(a) Fig. 2 The output of each population. a All parameters keep standard for all populations, k21 = k32 = k13 = 0. b All parameters keep standard except for A of population 1 which is 3.4 mv,
10 t(s)
(b)
15
20
10 0
−10 0 20 y(2)(t)(mv)
10
1 2n
2n
(4)
20
5
10
15
20
5
10
15
20
5
10 t(s)
15
20
10 0
−10 0 20 y(3)(t)(mv)
(1)
5
y(2)(t)(mv)
y(2)(t)(mv) y(3)(t)(mv)
y (t)(mv)
2
0 0
x˜k =
2n
10
y(3)(t)(mv)
y(1)(t)(mv)
4
0 0 6
X j = xˆk−1 ± ( n Px x ) j , j = 1, 2, . . . , 2n X˜ j = F(X j ), Y˜ j = H j ( X˜ j )
20
6
0 0 6
where xk is the system state vector at timestep k; yk is the measurement vector; F(·) is the system transition function; Hk is the measurement matrix; and Vk is the measurement noise vector and assumed to be Gaussian white noise satisfying E(Vk ) = 0, E(Vk , V jT ) = Rk δk j , where δk j is the Kronecker delta function. The UKF is summarized as the following recursive equations:
10 0
−10 0
(c)
k21 = k32 = k13 = 0. c All parameters keep standard except for A of population 1 which is 3.4 mv, k21 = k32 = k13 = 100
123
1294
X. Liu et al.
2n 1 ˜ ( X j − x˜k )( X˜ j − x˜k )T P˜x x = 2n
(7)
j=1
Px y =
2n 1 ˜ ( X j − x˜k )(Y˜ j − y˜k )T 2n
(8)
2n 1 ˜ (Y j − y˜k )(Y˜ j − y˜k )T + Rk 2n
(9)
j=1
Pyy =
j=1
−1 K = Px y Pyy , xˆk = x˜k + K (yk − y˜k )
(10)
Pˆx x = P˜x x − K PxTy ,
(11)
where n is the number of system states. When the information of measurement noise distribution is known completely, the UKF is able to achieve good performance. In practical applications, the statistics of noise can not be accurately obtained, and the filtering performance of UKF is affected seriously. So we combine the UKF with FIS to amend measurement noise adaptively. Define the residual vector as rk = yk − y˜k .
(12)
It is obvious that the value of rk reflects the degree to which the model fits the data, and the covariance of rk determines the degree of divergence. The theoretical residual variance is Pr = Pyy . The actual variance is defined as Cr [26] 1 Cr = M
k
ri ∗ riT ,
(13)
i=k−M+1
where M denotes the number of the latest residual vectors. A ratio of trace between the actual variance and the theoretic variance is qk =
T r (Cr ) , T r (Pr )
(14)
where T r (·) means the matrix trace. As seen in the previous chapter, the value of qk should be around 1 if the model is accurate. When the actual measurement noise suddenly increases, Cr increases and qk > 1. Hence, Pr should increase to reduce the influence from the measurement noise. On the contrary, when the measurement noise reduces, Cr reduces and qk < 1. Hence, Pr should decrease. The adjustment of Pr can be converted to the adjustment of Rk . From this point of view, Rk can be adjusted as follows Rk = s v Rk ,
123
(15)
where s is the regulation factor which is calculated by FIS, and v ≥ 0 indicates the extent of the amplification of s. If v = 0 is chosen, Rk keeps invariant. If v < 1, the adjustment action of s decreases and Rk is slowly changed. On the contrary, if v > 1, the adjustment action of s is amplified and Rk is greatly changed. Thus, the FAUKF algorithm is derived by replacing (9) in the recursive equations of the UKF with the following formulas: 2n 1 ˜ (Y j − y˜k )(Y˜ j − y˜k )T + Rk Pyy = 2n j=1
Rk = s v Rk 2n 1 ˜ Pyy = (Y j − y˜k )(Y˜ j − y˜k )T + Rk . 2n
(16)
j=1
It is noted that if v = 0, then the FAUKF is equivalent to the regular UKF. The regulation factor s is calculated by FIS. The specific process is as follows. FIS contains four parts including fuzzification of input qk , fuzzy databases, fuzzy inference, and defuzzification of output s. The fuzzy sets for qk and s include less, equal, and more. “less” denotes the set for qk and s are less than 1, “equal” denotes the set for qk and s are equal to 1, “more” denotes the set for qk and s are more than 1. The triangulares membership functions are used mainly in the premise which is specified by three parameters: (α, β, γ ) ⎧ 0 x ≤α ⎪ ⎪ ⎪ ⎨ x−α α ≤ x ≤ β (17) μ = f (x, α, β, γ ) = β−α γ −x ⎪ β ≤ x ≤ γ ⎪ γ −β ⎪ ⎩ 0 x ≥ γ, where α ≤ β ≤ γ , μ is the membership grade, β is the membership function’s center, and α and γ determine the endpoints. The membership functions for qk and s are shown in Fig. 3. The fuzzy rule is usually knowledge bases constituted by the if–then rule. It is given as follows: If qk ∈ less, then s ∈ less. If qk ∈ equal, then s ∈ equal. If qk ∈ more, then s ∈ more. The Mamdani minimax reasoning method is used for fuzzy logic inference, and the centroid method is used for the defuzzification. A closed-loop control strategy based on the FAUKF is given here to suppress the epileptiform spikes in the
Fuzzy adaptive unscented Kalman filter control equal
less
more
1
Degree of membership
1
Degree of membership
Fig. 3 The membership functions of FIS. a The membership functions of input of FIS. b The membership functions of output of FIS
1295
0.8 0.6 0.4 0.2
equal
less
0.8 0.6 0.4 0.2
0
0
0
0.5
1
1.5
2
3
3.5
4
(a) Fig. 4 The FAUKF control scheme
uk
more
xk
1
Unit delay
xk
H
5
yˆ k
yk
Cr Pr
F
4.5
(b)
qk
FIS
s
UKF
uk
L considered neural mass models, as shown in Fig. 4, where yk is the noisy measurement influenced by the random noise η; yˆk is the estimated output from the FAUKF; qk is the ratio of trace between the actual variance and the theoretic variance; s is the regulation factor; L is the feedback gain matrix; and u k is the control law with the form of u k = L ∗ yˆk .
(18)
4 Simulation In this section, simulations are provided to demonstrate the filtering effect of the FAUKF and the ability of the FAUKF-based closed-loop control (Fig. 4) to suppress epileptiform spikes in the neural mass models.
4.1 Filtering results The filtering effect of the FAUKF is conducted in the model of three coupled neural populations. The parameters of each population are kept at standard values, and the coupling constants k21 , k32 , k13 are set at 100. The variance of measurement noise is chosen to guarantee that qk is around 1 from the beginning to 5 s.
It is increased to four times from 5 to 10 s and eight times from 10 to 15 s, and then allowed to recover to the original value. Figure 5a presents the output without the interference of measure noise (the black trace) and the noisy measurements (the blue trace). An inset is given there to show the zoom-in on data. Figure 5b–f present the comparison of the application of FAUKF with different values of v (v = 0, 1.5, 3, 11, 12). The black trace is the output without the interference of measure noise, and the red trace is the estimated output. Observation shows that tshe filtering accuracy has a meaningful improvement as v is varied from 0 to 3. However, no meaningful improvement of the filtering accuracy is verified when v is varied from 3 to 11. The estimated output is divergent seriously from the output without the interference of measurement noise when v is increased to 12 or a larger value, as shown in Fig. 5f. Furthermore, it is known that the FAUKF is equivalent to the regular UKF when v = 0 is chosen. One can observe from Fig. 5b and d that the filtering accuracy of the UKF and the FAUKF has no obvious differences at the beginning and last 5 s. However, the filtering accuracy of FAUKF is higher than that of the UKF from the 5 to 15 s. This is due to the adaptive adjustment of the FAUKF for the measurement noise covariance.
123
1296
X. Liu et al. 6
4
5
10
15
(2)
4 2 0 0 6
10
15
4 2 0 0
5
10 t(s)
15
(1)
5
5
5
15
5
10
15
20
5
10
15
20
5
10 t(s)
15
20
4 2
4
2
0 0
20
15.5
(c)
6 4
2 0 14.5
15
15.5
2
5
10
15
y(1)(t)(mv)
(1)
10 t(s)
15
2
0 0 6
20
2
4
y (t)(mv)
15
0 14.5
(b)
6
0 0 6
10
4
0 0 6
20
4
(a)
4
15
2
0 0
20
10
4
0 0 6
20
y(3)(t)(mv)
5
15.5
2 0 0 6
20
15
y (t)(mv)
2
4 2
0 14.5
y(2)(t)(mv)
15.5
y (t)(mv)
y(2)(t)(mv)
15
4
y(3)(t)(mv)
0 14.5
0 0 6
y(3)(t)(mv)
2
y(1)(t)(mv)
y(1)(t)(mv)
4
6
4
2
4
0 14.5
15
15.5
2 0 0 6
20
10
2
5
10
15
y(1)(t)(mv)
6
5
0 0
20
5
10
15
20
5
10
15
20
5
10 t(s)
15
20
0 0
5
10
15
4 2 0 0 6
20
6
(2)
2
y (t)(mv)
y(2)(t)(mv)
(2)
y (t)(mv)
10 4
5
10
15
5
0 0
20
0 0
5
10 t(s)
(d)
15
20
4
(3)
2
y (t)(mv)
y(3)(t)(mv)
y(3)(t)(mv)
10 4
2
0 0
5
10 t(s)
(e)
15
20
5
0 0
(f)
Fig. 5 The output and its estimation from the FAUKF with different values of v. a The outputs with and without interference of measurement noise. b v = 0, c v = 1.5, d v = 3, e v = 11, f v = 12
123
Fuzzy adaptive unscented Kalman filter control
1297
It is possible to make the following conjectures based on above results: the increase in v significantly improves the filtering accuracy up to a limit. The increase of v beyond this limit has no further meaningful effect of improving the filtering accuracy or even deteriorates the filtering accuracy. The adaptive adjustment for the uncertain measurement noise results in the higher filtering accuracy of the FAUKF in comparison with the UKF.
is increased to 15 times from the 5 to 10 s, and 30 times from the 10 to 15 s , and then allowed to recover to the original value at the last 5 s. Figure 6a shows the noisy measurements. The feedback gain matrix L is chosen as L = [L 1 L 2 L 3 ]T , where L 1 0 0 0) L 1 = (0 0 0 0
L 2 0 0 0) L 2 = (0 0 0 0
L 3 0 0 0) L 3 = (0 0 0 0
the superscript T denotes the transpose of a matrix,
L 1,
L 2 , and
L 3 are constant numbers. The control action is imposed on population 1 by setting
L 2 and
L 3 at 0, and v is set at 2.5 in simulations. The control T energy is defined as u k ∗ u k . Several simulations are provided to demonstrate the effectiveness of the proposed FAUKF control. Figure 6b and c present the results with the application of the UKF control and the FAUKF control with
L 1 = −0.95 which is the least value that suppresses the epileptiform spikes for the FAUKF. The black trace is the output, and the red trace is its estimation. An inset of the output is given in Fig. 6b to show the zoom-in on data. An inset of the output and its estimation is also given in Fig. 6c to show the zoom-in on data. It is
4.2 Controlling results In this section, the FAUKF control given in Fig. 4 is used to control epileptiform spikes in the model of three coupled neural populations. Populations 2, 3 are assumed to exhibit normal activity by keeping all parameters standard, and population 1 is assumed to be hyperexcitable by keeping all parameters standard except for A which is 3.4 mv. The coupling constants k21 , k32 , k13 are set at 100. The corresponding outputs are shown in Fig. 2c, where the epileptiform spikes are observed. The variance of measurement noise is chosen to guarantee that qk is around 1 at the beginning 5 s. It
6
1000
−500 5
10
15
20
−1000 0 300 y(2)(t)(mv)
20 0
5
10
15
20
y(3)(t)(mv)
(3)
y (t)(mv)
40 20 0
15
15.5
5
10
15
20
5
10
15
20
5
10 t(s)
15
20
10 0 0
5
10
15
20
5
10
15
0
2 0 0 6
20
50
4
4 2
−50
−20 −40 0
−100 0 100
2
2 0 0 6
20
20
−10
100
15
0
−20 −40 0
200
10
(3)
y(2)(t)(mv)
40
5
4
0 14.5
y(2)(t)(mv)
−20 −40 0
0
(1)
0
4
y (t)(mv)
(1)
20
500
y (t)(mv)
y(1)(t)(mv)
y (t)(mv)
40
5
10 t(s)
(a)
15
20
−100 0
5
10 t(s)
(b)
15
20
0 0
(c)
Fig. 6 The output and its estimation with the application of different control methods. a The noisy measurement output. b The UKF control. c The FAUKF control
123
1298
X. Liu et al.
Table 1 The total control energy (mv2 ) over 20-s simulation realizations for different control methods The varied times of measurement noise
L1
The total control enenrgy
UKF
FAUKF
UKF
FAUKF
4.8
−1.058
−0.94
87729.49
73748.53
15.30
−1.08
−0.95
151026.92
74651.21
20.40
−1.30
−0.955
315081.23
75184.51
observed that the sustained spikes do not occur under the FAUKF control (comparing Figs. 2c, 6c), while they still occur under the UKF control (comparing Fig. 2c and the inset of Fig. 6b). Meanwhile, the estimated output from the UKF control has a serious deviation from the output. When the magnitude of
L 1 is increased to 1.08, the sustained spikes are suppressed just enough by using the UKF control. The FAUKF control can suppress the epileptiform spikes certainly at the same gain condition. However, the total control energy required in the latter control is less than that required in the former one, as shown in Table 1. Table 1 presents the critical value of
L 1 and the total control energy over 20 s simulation realizations for different measurement noises, for the UKF control and the FAUKF control. The variance of measurement noise is chosen to guarantee that qk is around 1 at the beginning 5 s and the last 5 s, and is varied for different times from 5 to 10 s and from the 10 to 15 s in simulations. It is observed that the magnitude of
L 1 and the total control energy in the FAUKF control are less than those in the UKF control under identical level of measurement noise. The gaps of the magnitude of
L 1 and the total control energy between two control methods are increased with the increasing level of measurement noise. As is expected, suppression of the epileptiform spikes in the underlying neural mass models with uncertain measurement noise can be reached using the FAUKF control. In comparison with the UKF control, the adaptive adjustment of the FAUKF makes the total control energy reduce greatly, especially when the varied amplitude of measurement noise is large.
5 Conclusion The problem of suppression of the epileptiform spikes in the neural mass models with uncertain measurement
123
noise has been studied by using the FAUKF closed-loop control strategy. It has been shown by simulations that the epileptiform spikes in the underlying models can be suppressed using the FAUKF control. In comparison with the UKF control, the application of the FAUKF control can reduce the total control energy. This is due to the role of the adaptive adjustment to the uncertain measurement noise that the FAUKF plays. The FAUKF that we described is generic and can be applied in any setting where there is a need to estimate the states of systems with uncertain measurement noise. Acknowledgments This research was supported by the National Natural Science Foundation of China (61004050, 61273260, 61172095, 51207144), the Specialized Research Fund for the Doctoral Program of Higher Education (20101333110006), and the Humanities and Social Sciences Fund, the Ministry of Education Foundation of China (12YJCZH021).
References 1. Wendling, F., Bellanger, J.J., Bartolomei, F., Chauvel, P.: Relevance of nonlinear lumped-parameter models in the analysis of depth-EEG epileptic signals. Biol. Cybern. 83, 367–378 (2000) 2. Gao, J.B., Hu, J., Tung, W.W.: Entropy measures for biological signal analyses. Nonlinear Dyn. 68, 431–444 (2012) 3. Ding, L., Hou, C.: Stabilizing control of Hopf bifurcation in the Hodgkin–Huxley model via washout filter with linear control term. Nonlinear Dyn. 60, 131–139 (2010) 4. Goodfellow, M., Schindler, K., Baier, G.: Intermittent spikewave dynamics in a heterogeneous, spatially extended neural mass model. NeuroImage 55, 920–932 (2011) 5. Goodfellow, M., Schindler, K., Baier, G.: Self-organised transients in a neural mass model of epileptogenic tissue dynamics. NeuroImage 59, 2644–2660 (2012) 6. Nevado-Holgado, A.J., Marten, F., Richardson, M.P., Terry, J.R.: Characterising the dynamics of EEG waveforms as the path through parameter space of a neural mass model: application to epilepsy seizure evolution. NeuroImage 59, 2374– 2392 (2012) 7. Krauss, G.L., Koubeissi, M.Z.: Cerebellar and thalamic stimulation treatment for epilepsy. Acta Neurochir. Suppl. 97, 347–356 (2007) 8. Pollo, C., Villemure, J.G.: Rationale, mechanisms of efficacy, anatomical targets and future prospects of electrical deep brain stimulation for epilepsy. Acta Neurochir. Suppl. 97, 311–320 (2007) 9. Benabid, A.L.: What the future holds for deep brain stimulation. Expert Rev. Med. Devices 4, 895–903 (2007) 10. Sunderam, S., Gluckman, B., Reato, D., Bikson, M.: Toward rational design of electrical stimulation strategies for epilepsy control. Epilepsy Behav. 17, 6–22 (2010) 11. Liu, X., Liu, H.J., Tang, Y.G., Gao, Q.: Fuzzy PID control of epileptiform spikes in a neural mass model. Nonlinear Dyn. 71, 13–23 (2013)
Fuzzy adaptive unscented Kalman filter control 12. Kalman, R.E.: A new approach to linear filtering and prediction problems. J. Basic Eng. 82, 35–45 (1960) 13. Wang, H., Han, Z.Z., Zhang, W., Xie, Q.Y.: Chaotic synchronization and secure communication based on descriptor observer. Nonlinear Dyn. 57, 69–73 (2009) 14. Choi, H.H., Jung, J.W.: Fuzzy speed control with an acceleration observer for a permanent magnet synchronous motor. Nonlinear Dyn. 67, 1717–1727 (2012) 15. Gholami, A., Markazi, A.H.D.: A new adaptive fuzzy sliding mode observer for a class of MIMO nonlinear systems. Nonlinear Dyn. 70, 2095–2105 (2012) 16. Chong, M., Postoyan, R., Nesic, D., Kuhlmann, L., Varsavsky, A.: A robust circle criterion observer with application to neural mass models. Automatica 48, 2986–2989 (2012) 17. Mariani, S., Ghisi, A.: Unscented Kalman filtering for nonlinear structural dynamics. Nonlinear Dyn. 49, 131–150 (2007) 18. Ali, J., Mirza, M.R.U.B.: Performance comparison among some nonlinear filters for a low cost SINS/GPS integrated solution. Nonlinear Dyn. 61(3), 491–502 (2010) 19. Schiff, S.J., Sauer, T.: Kalman filter control of a model of spatiotemporal cortical dynamics. J. Neural Eng. 5, 1–8 (2008) 20. Liu, X., Gao, Q.: Parameter estimation and control for a neural mass model based on the unscented Kalman filter. Phys. Rev. E 88(4), 042905 (2013)
1299 21. Liu, X., Gao, Q., Li, X.L.: Control of epileptiform spikes based on nonlinear unscented Kalman filter. Chin. Phys. B. 23, 010202 (2014) 22. Ren, H.W., Zhao, Y.D.: A kind of fuzzy adaptive Kalman filter and its application on INS/Doppler integrated navigation. Aerosp. Control 29, 3–6 (2011) 23. Xu, T.L., You, W.H., Cui, P.Y.: Research on INS/GPS integrated navigation system based on fuzzy adaptive Kalman filters. J. Astronaut. 26, 571–575 (2005) 24. Jwo, D.J., Yang, C.F., Chuang, C.H., Lee, T.Y.: Performance enhancement for ultra-tight GPS/INS integration using a fuzzy adaptive strong tracking unscented Kalman filter. Nonlinear Dyn. 73, 377–395 (2013) 25. Jansen, B.H., Rit, V.G.: Electroencephalogram and visual evoked potential generation in a mathematical model of coupled cortical columns. Biol. Cybern. 73, 357–366 (1995) 26. Ma, Z.X., Liu, Z.Y., Chen, M.: Application of adaptive Kalman filtering based on fuzzy logic to the integrated GPS/INS navigation. Inf. Control 35, 457–461 (2006)
123