# Active Disturbance Rejection Synchronization of Morris-Lecar Neurons

Authors: Wei Wei; Yanjie Shao; Min Zuo
Journal: Complexity
 (2018)

---

## Abstract

Synchronization of biological neurons is not only a hot topic, but also a difficult issue in the field of bioelectrical physiology. Numerous reported synchronization algorithms are designed on the basis of neural model, but they have deficiencies like relatively complex and poor robustness and are difficult to be realized. Morris-Lecar neuron is considered, and linear active disturbance rejection control (LADRC) is designed. Only one control input signal is utilized to synchronize membrane potentials of biological neurons. Meanwhile, in order to verify the robustness of synchronization, sinusoidal signal and parameter perturbations are introduced in numerical simulations. LADRC can still achieve satisfactory synchronization. Both theoretical and numerical simulation results show that LADRC is capable of estimating and cancelling disturbances and uncertainties. Neither accurate neural models nor concrete disturbance signal models are indispensable. A more practical and effective thought is provided to address the synchronization between neurons.


---

## Body

## 1. Introduction

Complex and powerful biological nervous systems are composed of a great number of connected neurons. By organizing and analysing experimental data of membrane potentials of biological neurons, it is found that synchronization of neuronal membrane potentials plays a crucial role in biological information processing and transmission. The synchronization and synchronized degree between biological neurons are vital for physiological function [1]. Actualizing the synchronization of biological neurons has become a hot and difficult topic in bioelectricity physiology, nonlinear system theory, control, and many other areas.For quantitative analysis of discharging characteristics, Hodgkin-Huxley (HH) model [2] was proposed by Hodgkin and Huxley in 1950s. It is the first mathematic model that describes dynamics of biological neuronal membrane potentials, which has turned into a milestone of bioelectric physiology. Subsequently, various neuronal models describing different biological neuronal discharging modes have been proposed, such as FitzHugh-Nagumo (FHN) model [3], Ghostburster model [4], Hindmarsh-Rose (HR) model [5], and Chay model [6].Based on the above models, in order to assist the abnormal nervous systems to operate normally, scholars from nonlinear control area began to study how to exert an external control signal, i.e., an external electrical stimulation, to drive different discharging modes between neurons to be synchronized. A variety of quantitative analyses, control, and synchronization studies have been carried out. For instance, according to HH neuron model, accurate feedback [7] has been applied to realize the discharging synchronization between two HH neurons; under assumptions that all states are available, nonlinear control is proposed [8], additionally, when only membrane potential is available, linear adaptive control has also been designed to realize the discharging synchronization of HR neurons [8]; when system states and model are known, stable feedback control [9] based on Lyapunov stability theory has been designed and the synchronization of FitzHugh–Nagumo (FHN) neurons is achieved; adaptive neural networkH∞ approach proposed in [10] is utilized to synchronize two Ghostburster neurons.Morris-Lecar neuron is a kind of simplified biological neuron model, its asymptotic behaviour has been analysed in [11, 12]. González has explained the pacemaker dynamics in special parameter space [13]. For the purpose of describing memory features of Morris-Lecar neuron, fractal dimension Morris-Lecar neuron model has also been proposed in [14]. In order to achieve the synchronization of Morris-Lecar neurons, linear feedback control [12], adaptive H∞ control [15], and other control algorithms have been reported, respectively.Actually, two points are of great importance in the synchronization of neurons by external stimulations. Firstly, the synchronization approach should be robust enough to disturbance. Owing to the weakness of electrical signals created by neurons, the neuronal membrane potentials are sensitive to disturbance. Model based control, whose performance largely depends on the controlled system model, cannot keep the satisfactory performance in presence of disturbance. Therefore, control algorithms which are robust enough to disturbance are able to guarantee the synchronization. Secondly, low power electrical stimulations play a critical role in the synchronization [16]. As we know, a high power electrical stimulation is harmful to biological tissues. Therefore, reducing the power of a stimulation will guarantee the safe of tissues and lengthen the lifetime of neurological implant.Active disturbance rejection control (ADRC) is the very kind of control, which has strong ability to reject disturbance. By rethinking the essence of control and making a discussion on control theory (i.e., a model analysis approach or a direct control approach?) Han clarified the nature of control [17] in the late 1980s and proposed ADRC [18]. It does not need faithful model information. Before disturbance corrupting system output, ADRC estimates and cancels disturbance to guarantee the system output to be still desired. Therefore, in presence of disturbance, less oscillations of a control signal are needed to regulate system output. In other words, less power is consumed in regulation. So far, theoretical and applied researches on ADRC have been becoming a hot topic and numerous results have been reported in such area [19–21].However, for ADRC, many parameters have to be determined. Performance of a closed-loop system, to a great extent, depends on the experience of tuning parameters. It does limit the applications of ADRC. In order to reduce difficulty of parameter tuning, Gao proposed linear active disturbance rejection control (LADRC) and provided bandwidth-parameterization based tuning method for LADRC [22].In this paper, the synchronization of Morris-Lecar neuronal membrane potentials is considered and LADRC is utilized to reject disturbance and achieve desired synchronization. The rest of the paper is organized as follows: problem description is presented in Section2, basic LADRC and the design of linear active disturbance rejection synchronization are provided in Section 3, simulation results are given in Section 4, and, finally, the concluding remarks are given in Section 5.


## 2. Problem Description

### 2.1. Morris-Lecar Neuron Model

Morris-Lecar neuron model can be expressed as follows [12]:(1)CdVdt=Iext-gLV-VL-gCaβVV-VCa-gKnV-VK,dndt=τVαV-n,where(2)αV=0.51+tanh⁡V-v3v4,βV=0.51+tanh⁡V-v1v2,τV=ϕcosh⁡V-v3v4V is the membrane potential, Iext is the externally applied DC current, t denotes the time measured in milliseconds, n is the activation variable for potassium, VK,VCa,VL are equilibrium potentials of potassium, calcium, and leakage currents, respectively. gK,gCa,gL represent corresponding maximum conductance of corresponding ionic currents and C stands for membrane capacitance. v1 is the midpoint potential at which the calcium current is half-activated, v2 is a constant that corresponds to the steepness of the activation voltage dependence, v3 is the activation midpoint potential at which the potassium current is half-activated, and v4 denotes the slope factor of the potassium activation. ϕ represents a constant. Here, values of each variable are listed in Table 1.Table 1
Variables of a Morris-Lecar neuron.
Variables Values C 5 g L 2 V L −60 g C a 4 V C a 120 g K 8 v K −80 ϕ 1/15 v 1 −1.2 v 2 18 v 3 2 v 4 17.4


### 2.2. Discharge Characteristic of a Morris-Lecar Neuron

Neuron owns different discharge characteristics according to different values of external electric fields. Discharge characteristics of a Morris-Lecar neuron, whenIext=50 and Iext=200, are shown in Figure 1, respectively.Figure 1
Membrane potential of a Morris-Lecar neuron whenIext are different.
(a)
I e x t = 50 (b)
I e x t = 200


### 2.3. Synchronization of Morris-Lecar Neurons

Master and slave Morris-Lecar neurons are defined as follows:(3)dVmdt=Iextm-gLVm-VL-gCaβVmVm-VCa-gKnmVm-VK1C,dnmdt=τVmαVm-nm,dVsdt=1CIexts-gLVs-VL-gCaβVsVs-VCa-gKnsVs-VK+u,dnsdt=τVsαVs-ns,where footnotes m and s represent the master neuron and the slave neuron, respectively.u is the external control input, i.e., the external electrical stimulation, which drives the slave neuron to act the same as the master neuron.In this paper, the synchronization of Morris-Lecar neuronal membrane potentials is to realizelimt→∞⁡(Vm-Vs)=0 by LADRC.


## 2.1. Morris-Lecar Neuron Model

Morris-Lecar neuron model can be expressed as follows [12]:(1)CdVdt=Iext-gLV-VL-gCaβVV-VCa-gKnV-VK,dndt=τVαV-n,where(2)αV=0.51+tanh⁡V-v3v4,βV=0.51+tanh⁡V-v1v2,τV=ϕcosh⁡V-v3v4V is the membrane potential, Iext is the externally applied DC current, t denotes the time measured in milliseconds, n is the activation variable for potassium, VK,VCa,VL are equilibrium potentials of potassium, calcium, and leakage currents, respectively. gK,gCa,gL represent corresponding maximum conductance of corresponding ionic currents and C stands for membrane capacitance. v1 is the midpoint potential at which the calcium current is half-activated, v2 is a constant that corresponds to the steepness of the activation voltage dependence, v3 is the activation midpoint potential at which the potassium current is half-activated, and v4 denotes the slope factor of the potassium activation. ϕ represents a constant. Here, values of each variable are listed in Table 1.Table 1
Variables of a Morris-Lecar neuron.
Variables Values C 5 g L 2 V L −60 g C a 4 V C a 120 g K 8 v K −80 ϕ 1/15 v 1 −1.2 v 2 18 v 3 2 v 4 17.4


## 2.2. Discharge Characteristic of a Morris-Lecar Neuron

Neuron owns different discharge characteristics according to different values of external electric fields. Discharge characteristics of a Morris-Lecar neuron, whenIext=50 and Iext=200, are shown in Figure 1, respectively.Figure 1
Membrane potential of a Morris-Lecar neuron whenIext are different.
(a)
I e x t = 50 (b)
I e x t = 200


## 2.3. Synchronization of Morris-Lecar Neurons

Master and slave Morris-Lecar neurons are defined as follows:(3)dVmdt=Iextm-gLVm-VL-gCaβVmVm-VCa-gKnmVm-VK1C,dnmdt=τVmαVm-nm,dVsdt=1CIexts-gLVs-VL-gCaβVsVs-VCa-gKnsVs-VK+u,dnsdt=τVsαVs-ns,where footnotes m and s represent the master neuron and the slave neuron, respectively.u is the external control input, i.e., the external electrical stimulation, which drives the slave neuron to act the same as the master neuron.In this paper, the synchronization of Morris-Lecar neuronal membrane potentials is to realizelimt→∞⁡(Vm-Vs)=0 by LADRC.


## 3. Linear Active Disturbance Rejection Synchronization

### 3.1. Brief Introduction to Linear Active Disturbance Rejection Control

In this paper, second-order linear active disturbance rejection control is designed to realize the synchronization. Assuming that a first-order unknown nonlinear controlled system can be described as(4)y˙=fy,y˙,w+bu,where y is the system output, w is external disturbance, f denotes unknown system dynamics, b is the control coefficient, and u is the control input signal. Here, LADRC is utilized, and the control input can be designed as(5)u=kpyr-z1-z2b0,where kp is the proportional coefficient, yr is the desired output, b0 is a tunable parameter, z1,z2 are outputs of the extended state observer (ESO), and ESO is designed as(6)z˙1=z2+l1y-z1+b0u,z˙2=l2y-z1,where z1 estimates the system output y, z2 estimates the integrated disturbance including internal uncertainties and unmodeled dynamics, and external disturbances, y is the system output, l1,l2 are gains of ESO, and u is the control input defined in (5).If we letz=[z1,z2]T,eo=y-z1, then extended state observer (6) can be rewritten as(7)z˙=Aoz+Bou+Leo,where(8)Ao=0100,Bo=0b0,L=l1l2.For the first-order unknown nonlinear controlled system, it can be rewritten as(9)y˙=fy,y˙,w+bu=fy,y˙,w+b-b0u+b0u≜f-+b0u.Here,f- is defined as f-≜f(y,y˙,w)+(b-b0)u, which represents the integrated disturbance.According to control input (5), when ESO performs well, i.e., z1≈y,z2≈f-, we have the closed-loop system as(10)y˙=f-+b0u=f-+kpyr-z1-z2≈kpyr-z1.Here, a nonlinear controlled system, by LADRC, is approximately linearized. By setting proper parameters, one can get the desired performance.


### 3.2. Closed-Loop Stability of Linear Active Disturbance Rejection Synchronization

LADRC is taken to realize the synchronization of Morris-Lecar neurons; synchronization structure of Morris-Lecar neurons is given in Figure2.Figure 2
Synchronization structure of Morris-Lecar neurons by LADRC.In Figure2, Vm is the membrane potential of the master neuron, and Vs is the membrane potential of the slave neuron, respectively. Let x1=[Vm,nm]T,x2=[Vs,ns]T, system model of the master and slave neurons can be rewritten as(11)x˙1=Ax1+Fx1,x˙2=Ax2+Fx2+Bu,where(12)A=-gLC000,B=10,Fx1=f1x1f2x1,Fx2=f1x2f2x2,f1x1=Iextm+gLVL-gCaβVmVm-VCa-gKnmVm-VKC,f1x2=Iexts+gLVL-gCaβVsVs-VCa-gKnsVs-VKC,f2x1=τVmαVm-nm,f2x2=τVsαVs-ns.Consider the disturbance, system (11) can be rewritten as(13)x˙1=Ax1+Fx1,x˙2=Ax2+Fx2+Bd+Bu,where d is bounded disturbance signal.Letec=x2-x1=even=Vs-Vmns-nm be synchronization error vector, and synchronization error system of Morris-Lecar neuronal membrane potentials can be obtained from (13):(14)e˙c=Aec+Fx2-Fx1+Bd+Bu;i.e.,(15)e˙v=-gLCev+f1x2-f1x1+d+u,e˙n=f2x2-f2x1.The first synchronization error equation in (15) can be written as(16)e˙v=-gLCev+f1x2-f1x1+d+u=-gLCev+f1x2-f1x1+d+1-b0u+b0u=Gev,x1,x2,u,d+b0u,where G(ev,x1,x2,u,d) can be regarded as the integrated disturbance that includes internal uncertainties, un-modeled dynamics, and external uncertainties.From Figure2, we can see that y=Vs-Vm. The output z1 of ESO estimates y=Vs-Vm, i.e., the synchronization error ev, and z2 estimates integrated disturbance G(∙).Leth=G˙(∙), and then(17)e˙v=Gev,x1,x2,u,d+b0u,G˙·=h.Letε=ε1ε2=ev-z1G-z2 be observation error of ESO; then observation error equation of ESO can be obtained by subtracting system (6) from system (17); i.e.,(18)ε˙=Aεε+Eh,where(19)Aε=-l11-l20,E=01.Lemma 1 (see [23]).
Ifh is globally Lipschitz on x1,x2, then extended state observer gain L=[l1,l2]T can be selected to make observation errors converge to zero asymptotically; i.e., limt→∞⁡ε1=limt→∞⁡ε2=0.Lemma 2 (see [23]).
If observation errors of the extended state observer are asymptotically convergent to zero, and then there exist appropriate parameters, for the closed-loop system, which can force the synchronization error to converge to zero asymptotically.According to Lemmas1 and 2, we have the following theorem.Theorem 3.
Ifh is globally Lipschitz on x1,x2, appropriate parameters of controller and extended state observer, i.e., kp and L=[l1,l2]T, can be chosen so as to make synchronization error system (15) be asymptotically stable; i.e., asymptotic synchronization of the membrane potentials can be realized.Proof.
Substitute linear active disturbance rejection control law (5) into system (16); considering yr=0 in synchronization, we have(20)e˙v=Gev,x1,x2,u,d+b0u=-kpz1+ε2=-kpev+kpε1+ε2;i.e.(21)e˙v=-kpev+kpε1+ε2.
Choosing suitable ESO gain vectorL=[l1,l2]T, one can make observation errors converge to zero. At the same time, suitable control parameter kp will drive the tracking error (21) to converge to zero asymptotically; i.e., synchronization error between Vs and Vm is asymptotic to zero.
Whenev=Vs-Vm=0, i.e., Vm=Vs, we have(22)e˙n=f2x2-f2x1=τVsαVs-ns-τVmαVm-nm=τVmαVm-ns-αVm+nm=τVmnm-ns=-τVmen.
Sinceτ(V)=ϕcosh⁡((V-v3)/v4), then(23)τV=ϕcosh⁡·=ϕexp⁡·+exp-1⁡·2≥ϕ>0.
Hence, the zero dynamic of synchronization error system (15) is also asymptotically stable.
Therefore, selecting suitable ESO gain vectorL=[l1,l2]T and control parameter kp, one can synchronize the membrane potentials of Morris-Lecar neurons asymptotically.


## 3.1. Brief Introduction to Linear Active Disturbance Rejection Control

In this paper, second-order linear active disturbance rejection control is designed to realize the synchronization. Assuming that a first-order unknown nonlinear controlled system can be described as(4)y˙=fy,y˙,w+bu,where y is the system output, w is external disturbance, f denotes unknown system dynamics, b is the control coefficient, and u is the control input signal. Here, LADRC is utilized, and the control input can be designed as(5)u=kpyr-z1-z2b0,where kp is the proportional coefficient, yr is the desired output, b0 is a tunable parameter, z1,z2 are outputs of the extended state observer (ESO), and ESO is designed as(6)z˙1=z2+l1y-z1+b0u,z˙2=l2y-z1,where z1 estimates the system output y, z2 estimates the integrated disturbance including internal uncertainties and unmodeled dynamics, and external disturbances, y is the system output, l1,l2 are gains of ESO, and u is the control input defined in (5).If we letz=[z1,z2]T,eo=y-z1, then extended state observer (6) can be rewritten as(7)z˙=Aoz+Bou+Leo,where(8)Ao=0100,Bo=0b0,L=l1l2.For the first-order unknown nonlinear controlled system, it can be rewritten as(9)y˙=fy,y˙,w+bu=fy,y˙,w+b-b0u+b0u≜f-+b0u.Here,f- is defined as f-≜f(y,y˙,w)+(b-b0)u, which represents the integrated disturbance.According to control input (5), when ESO performs well, i.e., z1≈y,z2≈f-, we have the closed-loop system as(10)y˙=f-+b0u=f-+kpyr-z1-z2≈kpyr-z1.Here, a nonlinear controlled system, by LADRC, is approximately linearized. By setting proper parameters, one can get the desired performance.


## 3.2. Closed-Loop Stability of Linear Active Disturbance Rejection Synchronization

LADRC is taken to realize the synchronization of Morris-Lecar neurons; synchronization structure of Morris-Lecar neurons is given in Figure2.Figure 2
Synchronization structure of Morris-Lecar neurons by LADRC.In Figure2, Vm is the membrane potential of the master neuron, and Vs is the membrane potential of the slave neuron, respectively. Let x1=[Vm,nm]T,x2=[Vs,ns]T, system model of the master and slave neurons can be rewritten as(11)x˙1=Ax1+Fx1,x˙2=Ax2+Fx2+Bu,where(12)A=-gLC000,B=10,Fx1=f1x1f2x1,Fx2=f1x2f2x2,f1x1=Iextm+gLVL-gCaβVmVm-VCa-gKnmVm-VKC,f1x2=Iexts+gLVL-gCaβVsVs-VCa-gKnsVs-VKC,f2x1=τVmαVm-nm,f2x2=τVsαVs-ns.Consider the disturbance, system (11) can be rewritten as(13)x˙1=Ax1+Fx1,x˙2=Ax2+Fx2+Bd+Bu,where d is bounded disturbance signal.Letec=x2-x1=even=Vs-Vmns-nm be synchronization error vector, and synchronization error system of Morris-Lecar neuronal membrane potentials can be obtained from (13):(14)e˙c=Aec+Fx2-Fx1+Bd+Bu;i.e.,(15)e˙v=-gLCev+f1x2-f1x1+d+u,e˙n=f2x2-f2x1.The first synchronization error equation in (15) can be written as(16)e˙v=-gLCev+f1x2-f1x1+d+u=-gLCev+f1x2-f1x1+d+1-b0u+b0u=Gev,x1,x2,u,d+b0u,where G(ev,x1,x2,u,d) can be regarded as the integrated disturbance that includes internal uncertainties, un-modeled dynamics, and external uncertainties.From Figure2, we can see that y=Vs-Vm. The output z1 of ESO estimates y=Vs-Vm, i.e., the synchronization error ev, and z2 estimates integrated disturbance G(∙).Leth=G˙(∙), and then(17)e˙v=Gev,x1,x2,u,d+b0u,G˙·=h.Letε=ε1ε2=ev-z1G-z2 be observation error of ESO; then observation error equation of ESO can be obtained by subtracting system (6) from system (17); i.e.,(18)ε˙=Aεε+Eh,where(19)Aε=-l11-l20,E=01.Lemma 1 (see [23]).
Ifh is globally Lipschitz on x1,x2, then extended state observer gain L=[l1,l2]T can be selected to make observation errors converge to zero asymptotically; i.e., limt→∞⁡ε1=limt→∞⁡ε2=0.Lemma 2 (see [23]).
If observation errors of the extended state observer are asymptotically convergent to zero, and then there exist appropriate parameters, for the closed-loop system, which can force the synchronization error to converge to zero asymptotically.According to Lemmas1 and 2, we have the following theorem.Theorem 3.
Ifh is globally Lipschitz on x1,x2, appropriate parameters of controller and extended state observer, i.e., kp and L=[l1,l2]T, can be chosen so as to make synchronization error system (15) be asymptotically stable; i.e., asymptotic synchronization of the membrane potentials can be realized.Proof.
Substitute linear active disturbance rejection control law (5) into system (16); considering yr=0 in synchronization, we have(20)e˙v=Gev,x1,x2,u,d+b0u=-kpz1+ε2=-kpev+kpε1+ε2;i.e.(21)e˙v=-kpev+kpε1+ε2.
Choosing suitable ESO gain vectorL=[l1,l2]T, one can make observation errors converge to zero. At the same time, suitable control parameter kp will drive the tracking error (21) to converge to zero asymptotically; i.e., synchronization error between Vs and Vm is asymptotic to zero.
Whenev=Vs-Vm=0, i.e., Vm=Vs, we have(22)e˙n=f2x2-f2x1=τVsαVs-ns-τVmαVm-nm=τVmαVm-ns-αVm+nm=τVmnm-ns=-τVmen.
Sinceτ(V)=ϕcosh⁡((V-v3)/v4), then(23)τV=ϕcosh⁡·=ϕexp⁡·+exp-1⁡·2≥ϕ>0.
Hence, the zero dynamic of synchronization error system (15) is also asymptotically stable.
Therefore, selecting suitable ESO gain vectorL=[l1,l2]T and control parameter kp, one can synchronize the membrane potentials of Morris-Lecar neurons asymptotically.


## 4. Simulation Studies

In order to verify LADRC in the synchronization of Morris-Lecar neurons, initial conditions are chosen the same as the ones taken in [12], and three groups of simulations are performed. In the first group, no disturbance is considered, and synchronization performance has been compared. In the second group, numerical results have been compared, and the synchronization of LADRC has been confirmed in presence of sinusoidal disturbance. In the last group, parameters of slave neurons are changed to simulate parameter perturbations of a controlled plant. Results also verify LADRC.Parameters of LADRC are chosen according to the approach proposed in [22]. Let control bandwidth ωc, bandwidth of ESO ωo, and control parameter kp satisfy kp=ωc,ωo=10ωc.Group  I. Synchronization is achieved in absence of disturbance. Control parameters and integral of absolute error (IAE, which can be calculated as ∫evdt) are present in Table 2.Table 2
Control parameters and IAE values (in absence of disturbance).
Control parameters Values IAE indexes Values k 1 20 IAE1 54.8163 k 2 20 ω c 26 IAE0 8.6921 b 0 −50In Table2, k1,k2 are the same values chosen as those given in [12]. IAE1 is the IAE value generated by controllers taken in [12], IAE0 is the IAE value produced by LADRC. ωc,b0 are parameters of LADRC. Synchronous response and comparisons of control inputs are shown in Figure 3. Figure 3(a) presents the synchronous response of LADRC. Figure 3(b) exhibits the comparison of synchronization errors between controllers presented in [12] and LADRC. In Figure 3(b), ev is the synchronization error generated by LADRC, while ev1 represents synchronization error generated by controllers proposed in [12]. Figure 3(c) is the comparison of en, which is produced by LADRC, and en1, which is generated by controllers utilized in [12]. Figure 3(d) is the comparisons among control input u of LADRC and two control inputs u1 and u2 designed in [12].Figure 3
Synchronous performance of Morris-Lecar neurons and the comparison of control inputs (in absence of disturbance).
(a) (b) (c) (d)From Table2 and Figure 3, it can be seen that LADRC, compared with two controllers used in [12], is able to get smaller synchronization error and IAE values with less control energy.Group  II. For the purpose of verifying synchronization in presence of disturbance, sinusoidal disturbance (d(t)=10sin⁡t) is taken. Synchronizations have been compared between controllers designed in [12] and LADRC. Controller parameters and IAE values are shown in Table 3.Table 3
Control parameters and IAE values (in presence of sinusoidal disturbance).
Control parameters Values IAE indexes Values k 1 20 IAE1 275.2547 k 2 20 ω c 26 IAE0 21.0479 b 0 −50In Table3, k1, k2, ωc, b0 are the same values taken as those given in Table 2. IAE1 and IAE0 are IAE values generated by controller taken in [12] and LADRC, respectively. Control is activated at 200 ms, and the sinusoidal disturbance signal is introduced from 400 ms. Simulation lasts 600 ms. Synchronization response and the comparisons are given in Figure 4.Figure 4
Synchronous response of Morris-Lecar neurons and the comparisons of control inputs (in presence of sinusoidal disturbance).
(a) (b) (c) (d) (e) (f)Figure4(a) presents trajectories of neurons as time passes; Figure 4(b) is the phase trajectories betweenVm andVs. It will be a line when synchronization is achieved. Obviously, Figure 4(b) shows that synchronization is obtained. Figures 4(c), 4(d), and 4(e) and IAE values given in Table 3 depict the fact that, with less control energy, smaller synchronization errors and IAE values can be obtained by LADRC. Simulation results manifest that ESO can estimate the disturbance signal effectively and then disturbance will be compensated to guarantee satisfactory synchronization.Group  III. For the purpose of verifying synchronization, parameters of slave neurons have also been changed. According to [13], gCa,gK can be chosen from [0,12] and [0,24], respectively. In simulations, gCa,gK of slave neurons have been changed by 100%; i.e., gCa=8 and gK=16, from 400 ms to the end of simulation. Controller parameters and IAE values are listed in Table 4.Table 4
Control parameters and IAE values (in presence of parameter perturbations).
Control parameters Values IAE indexes Values k 1 20 IAE1 259.0081 k 2 20 ω c 26 IAE0 22.7972 b 0 −50In Table4, k1, k2, ωc, b0 are the same values taken as those listed in Table 3. Control is also activated from 200 ms, and the changed parameters (see Figure 5(f)) are introduced from 400 ms. Simulation lasts 600 ms. Synchronization response and the comparisons are given in Figure 5.Figure 5
Synchronous response of Morris-Lecar neurons and the comparisons of control inputs (in presence of parameter perturbations).
(a) (b) (c) (d) (e) (f)Figure5 shows that when gCa and gK are changed, smaller synchronization errors can be achieved by LADRC. Both Figure 5 and IAE values given in Table 4 show the same fact that, with the help of ESO, LADRC is able to get better synchronization.Numerical results show that LADRC is able to achieve desired synchronization with stronger robustness and less power. It does satisfy two important issues in the synchronization of neurons with external electrical stimulations.


## 5. Conclusion

In this paper, LADRC is adopted to discuss the synchronization of Morris-Lecar neuronal membrane potentials. Theoretically, it has been analysed that LADRC can get asymptotic synchronization of neuronal membrane potentials by selecting appropriate parameters of ESO and controller. Meanwhile, from numerical results, it has been verified that LADRC can perform well in the synchronization, no matter disturbance and (or) uncertainties exist or not.Precise model is not necessary and less control energy is needed, both of the facts guarantee that LADRC is able to achieve better synchronization. Broader prospects in practical synchronization of membrane potentials in nervous systems can be expected.



---

*Source: 7656252-2018-06-05.xml*