# Exponential Stabilization of Coupled Hybrid Stochastic Delayed BAM Neural Networks: A Periodically Intermittent Control Method

Authors: Yunjian Peng; Birong Zhao; Weijie Sun; Feiqi Deng
Journal: Complexity
 (2018)

---

## Abstract

This paper considers exponential stabilization for a class of coupled hybrid stochastic delayed bidirectional associative memory neural networks (HSD-BAM-NN) with reaction-diffusion terms. A periodically intermittent controller is proposed to exponentially stabilize such an unstable HSD-BAM-NN, and sufficient conditions of the closed-loop BAM-NN system with exponential stabilization are derived by using Lyapunov-Krasovskii functional method, stochastic analysis techniques, and integral inequality property, which decide the basic parameters of the proposed controller. Furthermore, a framework to establish simulation algorithm with sampled states is presented to implement the stabilization controller. With a HSD-BAM-NN model of power synchronization in a photovoltaic (PV) array field, we illustrate numerical simulation results to verify the correctness and effectiveness of the proposed controller.


---

## Body

## 1. Introduction

Spurred by pioneering works on the neural networks models with BAM in Kosko [1–3], much attention has been paid to BAM-NN owning to its various applications in image processing, pattern recognition, automatic control [4, 5], associate memory, parallel computing [6], and optimization [7]. It is worth noting that the results of global stability for BAM models obtained in [1–3] require severe constraint conditions of symmetric connection weight matrix. In the neural networks with very large scale circuits, it is difficult for a practical NN system to satisfy the absolutely symmetric conditions in BAM models. Recently, the stability analysis of BAM neural networks has been paid considerable attention, and many stability conditions of such NN models have been reported in the published literature [8–22]. Since diffusion effects cannot be avoided in the neural networks, in a physical sense, when electrons are moving in asymmetric electromagnetic fields, it is more valuable to consider the alternative activation of neurons in an available space as well as in a given time interval. Therefore, the model of BAM neural networks could be formulated as partial differential equations (PDE) instead of only ordinary differential equations (ODE) [9, 11]. Based on PDE models, many contributions have been published focusing on the stability of BAM neural networks (BAM-NN) with reaction-diffusion terms by using Lyapunov functional method and LMI techniques [7, 9, 11–15, 17, 18, 23, 24].Haykin [25] proposed that in real neural networks systems, synaptic transmission is a stochastic process with noise released by random electronic fluctuations of the neurotransmitters and other disturbance. Following the experimental conclusions, Hossain and Anagnostou [26] further stated that noise posed a basic problem for information processing which affected all aspects of nervous system function in a nervous system. Hence, it is also important to study the effects of noise perturbations existing in neutral networks dynamics. Literature [8] investigated the stability of stochastic NN for the first time, which guided followers exploring novel results of such problems, especially for stochastic NN with reaction-diffusion terms, and some striking published contributions can be found in [7, 11–15, 17, 18, 27].It is worth mentioning that the hybrid BAM-NN driven by continuous-time Markov chains have been used to model many practical systems where they may experience abrupt changes in their structure and parameters. To consider the issues of system structures’ abrupt changes, hybrid disturbance, and unreliable subsystem interconnections, the evolution dynamics of BAM-NN could be modeled as jump systems [28, 29]. As Markovian jump system was first introduced in [10], two fundamental components have formed original ideas of considering a jumping system as continuous state described by differential equations, and discrete-time state described by a continuous-time finite-state Markovian process [28–32]. For theory and techniques recently developed to analyze BAM-NN with Markovian jump factors, here we mention [11, 22, 24, 28, 29].In recent literature, a variety of approaches have been published for the stabilization control of BAM neural networks with or without delays and reaction-diffusion terms which include feedback control [9, 21, 23, 32, 33], impulsive control [22, 24, 34], and intermittent control [35–38]. The type of control considered in this paper is intermittent control, which was first introduced to control nonlinear dynamical systems in [39] and has aroused much interest of researchers due to its merits in engineering applications. Different from continuous control approaches, intermittent control is more effective because the system output is measured intermittently rather than continuously [35]. Some novel contributions of intermittent control [35, 36, 38] based on the stability analysis of BAM-NN have been achieved in recent years. Nevertheless, the complex effects on BAM-NN with stimulation time-varying delays and stochastic reaction-diffusion terms have not been considered to use intermittent control in published results.Motivated by the above discussion, the main purpose of this paper is to investigate BAM-NN with delays and reaction-diffusion terms and focus on its exponential stabilization by designing a periodically intermittent controller [35, 36, 38]. The main novelties in this paper can be highlighted as follows: Firstly, the model under consideration covers the frequently investigated models which often are characterized by special cases in structures or systematic functions. Secondly, most available results have been concerned with the stability problems, and in this technical note, stochastic stabilization is taken into account and a periodically intermittent controller with more flexible conditions is first proposed for such BAM-NN with stochastic reaction-diffusion terms. Finally, a numerical simulation method is designed to simulate the behavior of the time-varying delays stochastic hybrid partial differential equations, which enriches the theory of delayed stochastic partial differential equation.Throughout the paper, we take(Ω,F,{Ft}t≥0,P) as a complete probability space with a filtration {F}t≥0 satisfying the usual conditions; i.e., {F}t≥0 is right-continuous and {F}0 contains all P-null sets. Let W(t), t≥t0, be one-dimensional Brownian motion defined on the probability space. Let r(t), t≥0, be right-continuous Markov chain on the probability space taking values in a finite-state space S={1,2,…,N} with generator Γ=(rij)N×N. And the transition probability from model i at time t to model j at time t+Δ is (1)Prt+Δ=j∣rt=i=γijΔ+oΔ,ifi≠j1+γiiΔ+oΔ,ifi=jwhere Δ>0 and limΔ→0o(Δ)/Δ=0. Here, γij≥0, i≠j, is the transition probability from model i to model j and γii=-∑j=1,j≠iN‍γij. We assume that the Markovian chain r(·) is independent of the Brownian motion W(·). It is well known that almost every sample path of r(·) is right-continuous step function with a finite number of simple jumps in any finite subinterval R+.


## 2. Preliminaries and Problem Formulation

In this paper, we consider the hybrid stochastic BAM neural network model with reaction-diffusion terms [11–13, 15], which are formulated as a couple of neurons’ dynamics indexed by a pair superscripts of (n,m)(2)dxint,λn=∑k=1ln∂∂λknDiknrt∂xint,λn∂λkn-ainrtxint,λn+∑j=1Nmbjinrtfjnxjmt,λm+∑j=1Nmcjinrtfjnxjmt-τmt,λm+uint,λn,λmdt+∑j=1Nmhjinxint,λn,xjmt,λm,xjmt-τmt,λmdWjtwhere n,m denote the index of a couple of neuron networks (NN) with N(n) and N(m) neurons, respectively; i.e., (n,m)=(1,2) or (2,1). The notations i=1,2,…,N(n) and j=1,2,…,N(m) are, respectively, the indices of the neurons in the nth NN and the mth NN. l(n) is the dimension of space variable vector λ(n), and λ(n)=(λ1(n), λ2(n),…,λl(n)(n))∈Ω0⊂Rl(n) where Ω0 is a compact set and measurable with smooth boundary ∂Ω0 in space Rl(n), 0<mesΩ0<+∞. The definition of λ(m) is the same as λ(n) and we let λ=(λ(n),λ(m)). xi(n)=xi(n)(t,λ(n)) and xj(m)=xj(m)(t,λ(m)) are the states of the ith neuron in nth NN and the jth neuron in mth NN at time t(t≥0) and in space of λ(n), λ(m), respectively. For simplicity, we define x(n)=(x1(n),…,xN(n)(n))T∈RN(n), x(m)=(x1(m),…,xN(m)(m))T∈RN(m), and x=((x(n))T,(x(m))T)T to denote the state vectors of the nth, mth, and whole coupled neural networks, respectively. τi(n)(t) denotes the time delay satisfying 0≤τ(n)(t)≤τ(n),τ˙(n)(t)≤τ0(n)<1 with constant τ(n),τ0(n). fj(n)(·) denotes the activation function of the jth neuron in mth NN stimulating the ith neurons in the nth NN. ai(n)(r(t))>0 denotes changing rate of the ith neuron under the condition that neural network is disconnected and no external additional activation exists. bji(n)(r(t)) is the connection weight of neurons in coupled NN and cji(n)(r(t)) denotes the corresponding delayed connection weight. The smooth function Dik(n)((r(t)))≥0 is the transmission diffusion operator of neurons, hji(n) is the stochastic disturbance function of neurons, and Wj(t) is the Brownian motions as noise acting on the transmission from the jth neuron to the ith neuron.According to [9, 10], we present the following initial value condition for (2):(3)xins,λn=ϕis,λn,s,λn∈-τn,0×Ω0which is with Dirichlet boundary value.(4)xint0,λn=0,t,λn∈-τn,+∞×∂Ω0For convenience, we denote all the possible models of BAM-NN as A(n)(r(t))=diag(a1(n)(r(t)), ⋯,aN(n)(n)(r(t))), B(n)(r(t))=(bji(n)(r(t)))N(n)×N(m), and C(n))(r(t))=(cji(n)(r(t)))N(n)×N(m).Definition 1.
A stochastic vectorx is the solution of system (2)-(4) if it satisfies the following conditions:
(i)x is adapted to {Ft}{t≥0},
(ii) forT0∈R+,x∈CF0b([0,T0]×Ω0;R(N(n)+N(m))), (5)Emaxλn∈Ω0∫0T0∑n=12xn2+∇xn2dt<+∞,(iii) for t∈R+, (6)∫Ω0xindλn=∫Ω0ϕi0,λndλn+∫Ω0∫0t-ainrsxins,λn+∫Ω0∫0t∑k=1ln∂∂λknDiknrs∂xins,λn∂λkndsdλn+∑j=1Nmbjinrsfjnxjms,λm+∑j=1Nmcjinrsfjnxjms-τms,λmdsdλ+∫Ω0∫0t∑j=1Nmhjinxins,λn,xjms,λm,xjms-τms,λmdWjsdλ.Definition 2.
The HSD-BAM-NN (2)-(4) is exponentially stable in pth moment if, for arbitrary model combinations of (A(r~), B(r~), C(r~))(r~∈S), the states of the system satisfy (7)limsupt→∞log⁡E∑n=12xnpt<0where x(n)p=∫Ω0xnpdλ(n)1/p, p≥2. To ensure the existence and uniqueness of the solution to system (2)-(4), we suppose the following assumptions.Assumption 3.
Fori=1,…,N(n), arbitrary s1,s2∈R, the activation function fi(n)(·) is bounded with fi(n)(0)=0 and (8)Li-≤fins1-fins2s1-s2≤Li+ where s1≠s2,Li-,Li+ are constants.Remark 4.
In the above assumption,Li-,Li+ could be positive, negative, or zero. Such an assumption is weaker than the one in [14] where it demands Li-≡-Li+.Assumption 5.
Fors1,s2,s3∈R, there exits nonnegative constant σji(n) satisfying the following. (9)hjins1,s2,s32≤σjins12+s22+s32
As noted in [6, 10], even if the parameters or time-varying delay in neural networks are appropriately chosen, neural networks may lead to some phenomena such as instability, divergence, oscillation, and chaos [10, 11, 40]. In order to stabilize the BAM-NN system (2)-(4), we introduce the following periodically intermittent controller(10)uint,λn=∑l=1NnKilnxln,μT≤t<μT+δ0,μT+δ≤t<μ+1Twhere T denotes the control period, μ=0,1,2,⋯ is the control periods number, δ(0<δ<T) is called the control time width, and Kil(n) are the control gains.
Under the periodically intermittent controller (10), the closed-loop systems of (2)-(4) can be described as follows.(11)dxin=∑k=1ln∂∂λknDiknrt∂xin∂λkn-ainrtxin+∑j=1Nmbjinrtfjnxjm+∑j=1Nmcjinrtfjnxjmt-τmt,λm+∑l=1NnKilnxlndt+∑j=1Nmhjinxin,xjm,xjmt-τmt,λmdWjt
For thepth moment of x(t,λ), we cite the following lemma to derive the stability conditions of the system.Lemma 6 (see [36]).
LetΩ0 be a super cuboid set and λ=(λk)1×l, λk≤θ(θ>0), k=1,2,…,l. If x(t,λ) is a real-valued continuous derivable function as x∈C((R+,Ω0);R), xt,λ∂Ω0=0, then(12)∫Ω0xt,λpdλ≤p2θ4∫Ω0xt,λp-2∂xt,λ∂λk2dλ.


## 3. Main Results

In this section, we present the sufficient conditions for stability of the controlled HDS-BAM-NN (11). For simplicity, the following notations are used to state and prove the main results.(13)ςin=minr~∈Sμr~n∑k=1ln4p-1D_iknpθk2+pa~in-∑j=1Nm∑l¯=1p-1b~jinpαl¯jinL~jpβl¯jin+c~jinpξl¯jinL~jpζl¯jin-p-12∑j=1Nm∑l¯=1p-2σjinpϵl¯jin+σjinpϵp-1jin+σjinpϵpjin-minr~∈Sμr~m∑j=1Nmb~ijmpα¯pjinN~ipβ¯pjin+p-12σijmpϵ¯p-1ijn+σijmpϵ¯pijn(14)kin=minr~∈Sμr~n∑k=1ln4p-1D_iknpθk2+pa~in-∑j=1Nm∑l¯=1p-1b~jinpαl¯jinL~jpβl¯jin+c~jinpξl¯jinL~jpζl¯jin+p-12∑l¯=1p-2σjinpϵl¯jin+σjinpϵp-1jin+σjinpϵpjin-minr~∈Sμr~m∑j=1Nmb~ijmpα¯pijnN~ipβ¯pijn+p-12σijnpϵ¯p-1ijn+σijnpϵ¯pijn(15)νin=maxr~∈Sμr~npKiin+∑l=1l≠iNn∑l¯=1p-1Kilnpη¯l¯il+∑l=1l≠iNnKlinpη¯pli(16)ηin=maxr~∈Sμr~m∑j=1Nmp-12σijnpϵ¯p-1ijn+σijnpϵ¯pijn+c~ijnpξ¯pijnN~ipζ¯pijnIn the above equations,a~i(n)=minr~∈Sai(n)(r~), b~ji(n)=maxr~∈Sbjinr~, c~ij(n)=maxr~∈Scjinr~, D_ik(n)=minr~∈SDik(n)(r~), L~j=maxLj-,Lj+, N~i=maxNi-,Ni+,μr~(n)>0,μr~(m)>0, and other variable parameters are given as αl¯ji(n),βl¯ji(n),ξl¯ji(n),ζl¯ji(n),ϵl¯ji(n),α¯l¯ji(n),β¯l¯ji(n),ξ¯l¯ji(n),ζ¯l¯ji(n),ϵ¯l¯ji(n)∈(0,1) and satisfy that ∑l¯=1p‍αl¯ji(n)=∑l¯=1p‍βl¯ji(n)=∑l¯=1p‍ξl¯ji(n)=∑l¯=1p‍ζl¯ji(n)=∑l¯=1p‍ϵl¯ji(n)=∑l¯=1p‍ϵ¯l¯ij(n)=∑l¯=1p‍α¯l¯ji(n)=∑l¯=1p‍β¯l¯ji(n)=∑l¯=1p‍ξ¯l¯ji(n)=∑l¯=1p‍ζ¯l¯ij(n)=1.In order to establish the sufficient conditions, the following two assumptions are further introduced.Assumption 7.
The following inequalities holds:(17)ςin-νin-maxr~∈S∑q=1Nγr~qnμqn-ηin1-τ0n>0kin+ϱin-maxr~∈S∑q=1Nγr~qnμqn-ηin1-τ0n>0where ϱi(n)>0.
Here we define a function(18)Hiε^i=ςin-νin-maxr~∈S∑q=1Nγr~qnμqn-ε^imaxr~∈Sμr~n-ηineε^iτn1-τ0nwhere ε^i(n)∈R+. For Hi(ε^i) being continuous derivable, by intermediate value theorem, we could conclude that there exists ε^i0>0 satisfying Hi(ε^i)>0 with ε^i∈(0,ε^i0). Similarly, choose ε^i1>0 and ε^i∈(0,ε^i1) to have the following function. (19)Fiε^i=kin+ϱin-maxr~∈S∑q=1Nγr~qnμqn-ε^imaxr~∈Sμr~n-ηieε^iτn1-τ0n Let ε=min{ε^i0,ε^i1}; we have the uniform inequality as Hi(ε)>0, Fi(ε)>0.Assumption 8.
The following inequality holds:ε-ϱ(T-δ)/μ¯T>0, where ϱ=maxmax1≤i≤N(n)ϱi(n),max1≤j≤N(m)ϱj(m), μ¯=minminr~∈Sμr~(n),minr~∈Sμr~(m).Theorem 9.
Under Assumptions3–5 and 7-8, the closed-loop system of (2)-(4) with the periodically intermittent controller (10) is exponentially stable in pth moment.Proof.
Choose a candidate average Lyapunov-Krasovskii function [41] (20)V^:C0,∞×Ω0;RNn+Nm×S×R+→R+V^x,rt,t=∫Ω0Vx,rt,tdλ,with(21)Vx,rt,t=∑n=12μrtn∑i=1Nneεtxinp+eετnt1-τ0n∑i=1Nnηin∫t-τntteεsxins,λnpdswhere μr(t)(n)>0.
By the generalized Itô formula [42], we have the following.(22)EV^x,rt,t=EV^ϕ,ψ,r0,0+E∫0t∫Ω0LVx,rs,sdλdsFor short expression, denote r(t)=r~, xj,τ(m)=xj(m)(t-τ(m)(t),λ(m)), xi,τ(n)=xi(n)(t-τ(n)(t),λ(n)). By Lemma 2.7 in [10], we can get the following.(23)LVx,r~,t=∑n,m=1,22,1εμr~n∑i=1Nneεtxinp+pμr~neεt∑i=1Nnxinp-1∑k=1ln∂∂λknDiknr~∂xin∂λkn-ainr~xin+∑j=1Nmbjinr~fjnxjm+∑l=1NnKilnxln+∑j=1Nmcjir~fjnxj,τm+μr~neεtpp-12∑i=1Nnxinp-2∑j=1Nmhjinxin,xjm,xj,τm2+∑q=1Nγr~qnμqneεt∑i=1nxinp+eετn1-τ0n∑i=1Nnηineεtxinp-eετn1-τ0n∑i=1Nnηieεt-τntxi,τnp1-τ˙nt+∑q=1Nγr~qneετn1-τ0n∑i=1Nnηin∫t-τntteεsxins,λnpds
By employing the absolute value inequality and noticing that∑q=1N‍γr~q(n)=0, we can obtain the following.(24)LVx,r~,t≤∑n,m=1,22,1∑i=1Nneεtεμr~nxinp+pμr~nxinp-1∑k=1ln∂∂λknDiknr~∂xin∂λkn+Kii-ainr~xin+pμr~nxinp-1∑j=1Nmbjinr~fjnxjm+pμr~nxinp-1∑j=1Nmcjinr~fjnxj,τm+μr~npp-12xinp-2∑j=1Nmσjinxin2+xim2+xj,τm2+∑l=1l≠iNnpμr~nxinp-1Kilxln+∑q=1Nγr~qnμqnxinp+eετn1-τ0nηinxinp-ηinxj,τmp
Applying the fundamental inequality, i.e.,a1p+a2p+⋯+app≥pa1a2⋯ap, (ah≥0,h=1,2,…,p), yields the following.(25)pμr~nxinp-1∑j=1Nmbjinr~fjnxjm≤μr~n∑j=1Nmpxinp-1b~jinL~jxjm=μr~n∑j=1Nmp∏l¯=1p-1b~jinαl¯jinL~jβl¯jinxinb~jinαpjinL~jβpjinxjm≤μr~n∑j=1Nm∑l¯=1p-1b~jinpαl¯jinL~jpβl¯jinxinp+μr~n∑j=1Nmb~jinpαpjinL~jpβpjinxjmp(26)μr~n∑l=1l≠iNnpxinp-1Kilxln=μr~n∑l=1l≠iNnp∏l¯=1p-1Kilη¯l¯ilxinKilη¯pilxln=μr~n∑l=1l≠iNn∑l¯=1p-1Kilpη¯l¯ilxinp+∑l=1l≠iNnKilpη¯l¯ilxlnp(27)pμr~nxinp-1∑j=1Nmcjinr~fjnxj,τm≤μr~n∑j=1Nmpxinp-1c~jinL~jxj,τm=μr~n∑j=1Nmp∏l¯=1p-1c~jinξl¯jinL~jζl¯jinxinc~jinξpjinL~jζpjinxj,τm≤μr~n∑j=1Nm∑l¯=1p-1c~jinpξl¯jinL~jpζl¯jinxinp+μr~n∑j=1Nmc~jinpξpjinL~jpζpjinxj,τmp(28)μr~npp-12xinp-2∑j=1Nmσ1ijxin2+xim2+xj,τm2=μr~npp-12∑j=1Nmp∏l¯=1p-2σjinϵl¯jinxinσjinϵp-1jinxinσjinϵpjinxin+σjinϵp-1jinxjmσjinϵpjinxjm+σjinϵp-1jinxj,τmσjinϵpjinxj,τm=μr~np-12∑j=1Nm∑l¯=1p-2σjinpϵl¯jinxinp+μr~np-12∑j=1Nmσjinpϵp-1jin+σjinpϵpjinxinp+xjmp+xj,τmp
Substituting (25)-(28) into (24) leads to the following.(29)LVx,r~,t≤∑n,m=1,22,1∑i=1Nneεtεμr~n+pμr~nkii-pμr~nain~+∑q=1Nγr~qnμqn+μr~np-12∑j=1Nm∑l¯=1p-2σjinpϵl¯jin+μr~n∑j=1Nm∑l¯=1p-1b~jinpαl¯jinL~jpβl¯jin+c~jinpξl¯jinL~jpζl¯jin+μr~n∑l=1l≠iNn∑l¯=1p-1Kilpη¯l¯il+∑j=1Nmμr~mb~ijmpα¯pijnN~ipβ¯pijn+μr~nσjinpϵp-1jin+σjinpϵpjin+μr~mσijnpϵ¯p-1ijn+σijnpϵ¯pijnp-12xinp+μr~n∑l=1l≠iNnKilpη¯pilxlnp+eεtηin1-τ0n+μr~m∑j=1Nmc~ijmpξ¯pijnN~ipζ¯pijnxj,τmp+μr~mp-12∑j=1Nmσijnpϵ¯p-1ijn+σijnpϵ¯pijnxj,τmp-ηinxj,τmp+∑i=1Nneεtpμr~nxinp-1∑k=1ln∂∂λknDiknr~∂xin∂λkn
Let(30)
ς
¯
i
n
=
m
i
n
r
~
∈
S
μ
r
~
n
p
a
~
i
n
-
∑
j
=
1
N
m
∑
l
¯
=
1
p
-
1
b
~
j
i
n
p
α
l
¯
j
i
n
L
~
j
p
β
l
¯
j
i
n
+
c
~
j
i
n
p
ξ
l
¯
j
i
n
L
~
j
p
ξ
l
¯
j
i
n
-
p
-
1
2
∑
j
=
1
N
m
∑
l
¯
=
1
p
-
2
σ
j
i
n
p
ϵ
l
¯
j
i
n
-
p
-
1
2
∑
j
=
1
N
m
σ
j
i
n
p
ϵ
p
-
1
j
i
n
+
σ
j
i
n
p
ϵ
p
j
i
n
-
m
i
n
r
~
∈
S
μ
r
~
m
∑
j
=
1
N
m
b
~
i
j
m
p
α
¯
p
i
j
n
N
~
i
p
β
¯
p
i
j
n
+
p
-
1
2
σ
i
j
n
p
ϵ
¯
p
-
1
i
j
n
+
σ
i
j
n
p
ϵ
¯
p
i
j
n
and by (22), we get the following inequality.(31)EV^x,rt,t≤EV^ϕ,ψ,r0,0-∑n,m=1,22,1E∫0t∫Ω0∑i=1Nneεsς¯in-εmaxr~∈Sμr~n-νin-maxr~∈S∑q=1Nγr~qnμqn-eετnηin1-τ0nxins,λnpdλnds-E∫0t∑i=1neεs∫Ω0pμrsxins,λnp-1∑k=1ln∂∂λknDiknr~∂xin∂λkndλndsAccording to Dirichlet boundary conditions, by Lemma 6, we can derive the integral part of evolution dynamics as follows.(32)
∫
Ω
0
p
x
i
n
s
,
λ
n
p
-
1
∑
k
=
1
l
n
∂
∂
λ
k
n
D
i
k
n
r
~
∂
x
i
n
s
,
λ
n
∂
λ
k
n
d
λ
n
≤
-
∑
k
=
1
l
n
4
p
-
1
D
i
k
n
r
~
p
θ
k
2
∫
Ω
0
x
i
n
s
,
λ
n
p
d
s
≤
-
∑
k
=
1
l
n
4
p
-
1
D
_
i
k
n
p
θ
k
2
∫
Ω
0
x
i
n
s
,
λ
n
p
d
sUsing Ito^’s formula EVt′=E(LV(t)), we obtain(33)EV^x,rt,t≤EV^ϕ,ψ,r0,0-∑n,m=1,22,1E∫0t∫Ω0∑i=1Nneεsς¯in+minr~∈Sμr~n∑k=1ln4p-1D_iknpθk2-εmaxr~∈Sμr~n-νin-maxr~∈S∑q=1Nγr~qnμqn-eετnηin1-τ0nxins,λnpdλnds≤EV^ϕ,ψ,r0,0-∑n,m=1,22,1E∫0t∫Ω0∑i=1Nneεsς¯in-εmaxr~∈Sμr~n-νin-maxr~∈S∑q=1Nγr~qnμqn-eετnηin1-τ0nxins,λnpdλnds≤EV^ϕ,ψ,r0,0where (t,λ(n))∈[μT,μT+δ]×Ω0.
Similarly, when(t,λ(n))∈[μT+δ,(μ+1)T]×Ω0, we have(34)EV^x,rt,t≤EV^xμT+δ,λ,rμT+δ,μT+δ-E∫0t∫Ω0∑i=1Nneεskin+ϱin-εmaxr~∈Sμr~n-νin-maxr~∈S∑q=1Nγr~qnμqn-eετnηin1-τ0nxins,λnpdλds≤EV^xμT+δ,λ,rμT+δ,μT+δ+ϱμ¯E∫0tV^x,rs,sdswhere ϱ=max{max1≤i≤N(n)ϱi,max1≤j≤N(m)ϱj(m)}, μ¯=min{minl∈Sμr~(n),minl∈Sμr~(m)}.
By using the Gronwall inequality to (34) it yields the following.(35)EV^x,rt,t≤EV^xμT+δ,λ,rμT+δ,μT+δeϱ/μ¯t-μT-δ
From (33)-(35), we can conclude that
(I) for(t,λ)∈([0,δ),Ω0) and (t,λ)∈([T,T+δ),Ω0), by (33), we can, respectively, have (36)EV^x,rt,t≤EV^x0,λ,r0,0, and(37)EV^x,rt,t≤EV^xT,λ,rT,T≤EV^x0,λ,r0,0eϱ/μ¯T-δ(II) for (t,λ)∈([δ,T),Ω0) and (t,λ)∈([T+δ,2T),Ω0), by (35), we can, respectively, have (38)EV^x,rt,t≤EV^x0,λ,r0,0eϱ/μ¯t-δ, and(39)EV^x,rt,t≤EV^xT+δ,λ,rT+δ,T+δeϱ/μ¯t-T-δ≤EV^x0,λ,r0,0eρ/μ¯t-2δ.Repeating above procedure (I)-(II), for (t,λ)∈[μT,μT+δ)×Ω0, then μ≤t/T and(40)EV^x,rt,t≤EV^xμT,λ,rμT,μT≤EV^x0,λ,r0,0eϱ/μ¯TT-δt and for (t,λ)∈[μT+δ,(μ+1)T)×Ω0, then t/T<μ+1 and(41)EV^x,rt,t≤EV^xμT+δ,λ,rμT+δ,μT+δeϱ/μ¯t-μT-δ≤EV^x0,λ,r0,0eϱ/μ¯TT-δt. Hence, for arbitrary (t,λ)∈[0,∞)×Ω0, (42)EV^x,rt,t≤EV^x0,λ,r0,0eϱ/μ¯TT-δt. By calculation, we obtain(43)EV^x,rt,t≥eεtμ¯E∫Ω0∑i=1Nnxinpdλn+∫Ω0∑j=1Nmxjmpdλn(44)EV^x0,λ,r0,0≤maxr~∈Sμr~nE∫Ω0∑i=1Nnϕi0,λnpdλn+sup-τn≤s≤0max1≤i≤Nnηinτneετn1-τ0n×E∫Ω0∑i=1Nnϕis,λnpdλn+maxr~∈Sμr~mE∫Ω0∑j=1Nmψj0,λnpdλn+sup-τn≤s≤0max1≤j≤Nmηmjτneετn1-τ0nE∫Ω0∑j=1Nmψjs,λnpdλn=μ0so(45)E∫Ω0∑i=1Nnxinpdλn+∫Ω0∑j=1Nmxjmpdλn≤μ0μ¯e-ε-ϱ/μ¯TT-δt.Under Assumption 8, Theorem 9 holds.
In Theorem9, if we let ςi(n)=ki(n), then αl¯ji(n)=αl¯ji(n), βl¯ji(n)=βl¯ji(n), ξl¯ji(n)=ξl¯ji(n), ζl¯ji(n)=ζl¯ji(n), ϵl¯ji(n)=ϵl¯ji(n), ϵ¯l¯ji(n)=ϵ¯l¯ji(n), α¯l¯ji(n)=α¯l¯ji(n), β¯l¯ji(n)=β¯l¯ji(n), ξ¯l¯ji(n)=ξ¯l¯ji(n), ζ¯l¯ji(n)=ζ¯l¯ji(n), ϱi=-νi(n). Thus, we can obtain the following corollary.Corollary 10.
Under Assumptions3–5, system (2)-(4) with periodically intermittent controllers (10) is exponentially stable in pth moment if the following conditions hold:
(I)νi(n)<0,ςi(n)-νi(n)-maxr~∈S∑q=1N‍γr~q(n)μq(n)-ηi(n)/1-τ0(n)>0.
(II)ε-ν¯(T-δ)/μ¯T>0, where ν¯=maxmax1≤i≤Nnνin,max1≤j≤Nmνjm,μ¯=min{minr~∈Sμr~(n),minr~∈Sμr~(m)}.
In Theorem9, letting ςi(n)=1/p, ςj(m)=1/p which means that αl¯ji(n)=1/p, βl¯ji(n)=1/p, ξl¯ji(n)=1/p, ζl¯ji(n)=1/p and ϵl¯ji(n)=1/p, we can get the following.(46)ς~in=k~in=minr~∈Sμr~n∑k=1ln4p-1D_iknpθk2+pa~in-p-1∑j=1Nmb~jinL~j+c~jinL~j-pp-12∑j=1Nmσjin-minr~∈Sμr~ml∑j=1Nmb~ijmN~i+p-1σijn
We can also have a further extended corollary.Corollary 11.
Under Assumptions3 to 5, the BAM-NN system (2)-(4) with periodically intermittent controllers (10) is exponentially stable in pth moment if the following conditions hold:
(I)ς~i(n)-ν~i(n)-maxr~∈S∑q=1N‍γr~q(n)μq(n)-η~i(n)/1-τ0(n)>0.
(II)ς~i(n)+ϱ~i(n)-maxr~∈S∑q=1N‍γr~q(n)μq(n)-η~i(n)/1-τ0(n)>0, ϱ~i(n)>0.
(III)ε~-ϱ~(T-δ)/μ¯T>0, where ϱ~=maxmax1≤i≤Nnϱ~in,max1≤j≤Nmϱ~jm,μ¯=minminr~∈Sμr~n,minr~∈Sμr~m.Remark 12.
Unlike the result in [22], which only considered the mean square exponential stability of stochastic BAM-NN with time-varying delays and reaction-diffusion terms, we design a periodically intermittent controller to exponentially stabilize the unstable neural network in pth moment. Moreover, the controller (10) is linear which can be easily implemented in practice.Algorithm 13.
For the periodically intermittent controller (10) of the BAM-NN system (2)-(4), we summarize the following algorithm to implement the controller.(1) By instigating a practical system’s structure and dynamics’ characterization with its disturbance, we can develop a model of the BAM-NN system, i.e., S , γ i j . Supposing and examining the transition probability P ( · ) as well as the main parameters of each mode, we can have A ( n ) ( r ( t ) ) , B ( n ) ( r ( t ) ) , C ( n ) ) ( r ( t ) ) and other relations of their states.(2) By using (13)-(16) and Corollary 11, we can calculate the values of stabilization indices of the system and the controller gains K i j with the parameters of BAM-NN.(3) Choose a numerical solution to the stochastic partial differential equation (SPDE) to simulate the sample states of BAM-NN. Here we use a so-called estimation-correction method, which is based on the following main steps [43].Firstly , with given steps H , L of simulation time interval t ∈ 0 , t f and space λ k ∈ [ λ min ⁡ k , λ max ⁡ k ] ∈ Ω 0 , we can get the grids of time and space, i.e., t h = h Δ t , h = 0,1 , 2 , … , H , λ l = l Δ λ , l = 0,1 , 2 , … , L , by which we define a numerical solution to BAM-NN as x l h = x ( t h , λ l ) = x 1 t h , λ l , x 2 t h , λ l T . Then we denote ΣF(xlh)=A(r(th))xlh+B(r(th))f(xlh)+C(r(th))f(xlh-τh)+Kxlh+h(xlh,xlh-τh)W(th)-W(th-1)/Δt, where τh≈τ(th)/Δt is an integer times of τ(th) to Δt.Secondly , according to the basic formulation of SPDE (2) with its discretized time and space, we have an estimation formula for the numerical solution x l h (47) x l h + Δ t 2 Σ F x l h = - ρ 2 D x l + 1 h + 1 / 2 + I + ρ D x l h + 1 / 2 - ρ 2 D x l - 1 h + 1 / 2 where x l h + 1 / 2 = x ( t h + 1 / 2 Δ t , λ l ) , ρ = Δ t / Δ λ 2 , I is an identity matrix with suitable dimensions. By the boundary function (4), we can calculate the values of x 0 h + 1 / 2 and x L h + 1 / 2 . Then we can solve three diagonal linear equations from (47) to get an estimated value of x l h + 1 / 2 .Finally , based on the estimated value of x l h + 1 / 2 , we have a further correction formula as follows. (48) - ρ 2 D x l + 1 h + 1 + I + ρ D x l h + 1 - ρ 2 D x l - 1 h + 1 = ρ 2 D x l + 1 h + I - ρ D x l h + ρ 2 D x l - 1 h + Δ t Σ F x l h + 1 / 2 Similarly, we also transfer formulation (48) to three diagonal linear equations and calculate the numerical solution x l h with permitted calculation errors.(4) Comparing the sampled data profiles of the simulation states from (1)-(3) in Algorithm13 with the measured data profiles from the aimed practical system, we can identify the validness of the model and performance of the controlled BAM-NN system. Furthermore, we can improve the parameters and control gains to obtain a better model by repeating steps (1)-(3).


## 4. Numerical Simulation of an Illustrative Application

For a grid-connected photovoltaic (PV) power generation system with series-connected inverters, which is illustrated by it is equivalent topology structure in Figure1, every branch-circuit with a group of PV panels penetrates power into utility grid via a inverter. Since the current/power (Ii,Pi) generated by the ith group of PV panels is strongly dependent on operating conditions and field factors, such as sun geometric locations, their irradiation levels (Ri) of the sun and the ambient temperature (Ti) stochastically fluctuate with the environmental factors of PV power fields. Thus it is technically necessary to maintain power synchronization of series-connected inverters to improve their output power and standard voltage and current [44]. For this purpose, we need to develop a model with power and current difference between every couple-connected PV inverter based on basic photovoltaic model of PV panels [44, 45], which can be modeled as a HSD-BAM-NN ((2)-(4)), and the power/current difference (ΔPij,ΔIij) is taken as states x(t) with irradiation levels (Ri) as the space variable λ. It is assumed that the temperature difference can be ignored. By [44, 45], we formulate every part of BAM-NN ((2)-(4)) and calculate the parameters in the model; i.e., the modeled BAM-NN (2) has the same coefficient matrix, which is denoted as follows.(49)Dr=4.966.34008.056.76Ar=0.28000.18000.42000.35Br=0.30.4120.3000.20.210-0.5Cr=0.10.30.20.300-0.20.2-0.4-0.1Figure 1
Equivalent topology structure of a grid-connected PV power system with series multiple inverters.Firstly, we get the initial value conditions(50)x1s,λ1=escos⁡2πλ1-1x2s,λ2=e2ssin⁡4πλ2 where (s,λ(2))∈[-2,0]×Ω0, Ω0=[0.5,8.6]. And the boundary value functions are taken as follows.(51)x1t,0.5=2.3,x1t,8.6=6.95x2t,0.5=1.8,x2t,8.6=0.34,t≥-2 Denote τ1(t)=0.8sin⁡(π/4t+0.2), τ(2)(t)=1/1+0.3et; the generator of the Markovian chain is as follows.(52)Γ=-454512-12The activation functions aref(x)=3/4sin⁡x+1/4x, g(x)=1/2(x+1-x-1), and the stochastic disturbed functions are h12(s1,s2,s3)≤0.01(s12+s22+s32), h22(s1,s2,s3)≤0.02(s12+s22+s32).According to the given structure and parameters of the BAM-NN (2)-(4), we can perform the numerical simulations in instability and stabilization cases.(1) Instability. Using Algorithm 13 in Section 3, the sample states x(t,λ) of BAM-NN (2) are calculated, and the surfaces of (t,λ) versus x(t,λ) are shown in Figures 2 and 3, while the t versus x(t,λ) profile curves are shown in Figures 4 and 5. These figures show the instability behavior.Figure 2
Unstable sample state surface ofx(1)(t,λ) in HSD-BAM-NN’s simulation with no control.Figure 3
Unstable sample state surface ofx(2)(t,λ) in HSD-BAM-NN’s simulation with no control.Figure 4
Unstable sample state profiles ofx(1)(t,λ) in HSD-BAM-NN’s simulation with no control.Figure 5
Unstable sample state profiles ofx(2)(t,λ) in HSD-BAM-NN’s simulation with no control.(2) Stabilization. Let p=2,μ1(1)=μ1(2)=1,μ2(1)=μ2(2)=2. By (13)-(16), we can calculate that L-=-1/2,L+=1,N-=0,N+=1,τ(n)=1,τ0(n)=1/4, ς~i(n)=k~i=-1.43, η~i=2.84, ς~j(2)=k~j(2)=-1,33, η~j(2)=2.82,j=1,2. The gain coefficients of periodically intermittent controllers (10) are obtained as follows by sufficient conditions (53)K=-2.952.863.421.6Furthermore, we can get the indices’ values of the stable conditions of HSD-BAM-NN with ν~1(1)=ν~2(1)=-18.56, ν~1(2)=ν~2(2)=-25.02, ς~i(1)-ν~i(1)-maxr~∈S∑q=1N‍γr~q(1)μq(1)-ηi(1)/1-τ0(1)>0, ς~i(2)-ν~i(2)-maxr~∈S∑q=1N‍γr~q(2)μq(2)-ηi(2)/1-τ0(2)>0, i=1,2. Also, we have ε~1(1)=ε~2(1)=3.15,ε^1(2)=ε^2(2)=8.399. From (III) in Corollary 11, we know that ε~-ν^(T-δ)/μ¯T>0. Letting δ=3.6,T=5, we can calculate the sample stabilized states of the controlled system. The trend surfaces of (t,λ) versus x(t,λ) are shown in Figures 6 and 7 with the profile curves of t versus x(t,λ) in Figures 8 and 9. From the states’ trends, we can see that an unstable system (2)-(4) can be stabilized by the controller (10) with appropriate parameters, which theoretically shows that the output power of connected PV panels can be synchronized to a given level with its BAM-NN model, even though the power stochastically fluctuate.Figure 6
Sample state surface ofx(1)(t,λ) in HSD-BAM-NN’s simulation with controller (10).Figure 7
Sample state surface ofx(2)(t,λ) in HSD-BAM-NN’s simulation with controller (10).Figure 8
Sample state profiles ofx(1)(t,λ) in HSD-BAM-NN’s simulation with controller (10).Figure 9
Sample state profiles ofx(2)(t,λ) in HSD-BAM-NN’s simulation with controller (10).


## 5. Conclusion

In this paper, a hybrid stochastic delayed BAM neural network is considered for its stabilization problem and a periodically intermittent controller is designed to stabilize an unstable HSD-BAM-NN with an exponential convergence property. The sufficient conditions of exponential stabilization of HSD-BAM-NN are derived by Lyapunov-Krasovskii functional method, stochastic analysis techniques, and integral inequality. And the framework is established to give a solution algorithm to the sufficient conditions. The simulation results of the grid-connected photovoltaic (PV) power generation system verify the effectiveness of the proposed controller.



---

*Source: 1708935-2018-12-02.xml*