
The escalating penetration of renewable energy, notably wind and solar, introduces unique complexities into power systems, particularly in frequency and voltage regulation. Current grid codes are inadequate for these emerging dynamics, necessitating significant enhancements to cope with the evolving energy landscape. The paper highlights several pivotal areas of improvement, including increased flexibility within power systems, integration of energy storage systems, expansion of ancillary services, the inclusion of grid-forming inverters, and the international harmonization of grid codes. The research underscores that conventional power system methodologies, primarily reliant on traditional power plants, fail to manage the fluctuating supply-demand dynamics of renewable energy effectively. By proposing improvements in grid codes, the research contributes towards resolving this issue. Furthermore, the paper underlines the paramount importance of international harmonization of grid codes for system interoperability, efficient operation, and exchange of best practices across diverse regions. Through its exhaustive exploration and recommendations, the study empowers policymakers, grid operators, and energy producers to advance grid code frameworks. Consequently, this facilitates renewable energy integration, ensures grid stability, and paves the way for a more sustainable energy future.
Citation: Anshu Murdan, Iqbal Jahmeerbacus, S Z Sayed Hassen. Challenges of existing grid codes and the call for enhanced standards[J]. Clean Technologies and Recycling, 2023, 3(4): 241-256. doi: 10.3934/ctr.2023015
[1] | Huaying Liao, Zhengqiu Zhang . Finite-time stability of equilibrium solutions for the inertial neural networks via the figure analysis method. Mathematical Biosciences and Engineering, 2023, 20(2): 3379-3395. doi: 10.3934/mbe.2023159 |
[2] | Pan Wang, Xuechen Li, Qianqian Zheng . Synchronization of inertial complex-valued memristor-based neural networks with time-varying delays. Mathematical Biosciences and Engineering, 2024, 21(2): 3319-3334. doi: 10.3934/mbe.2024147 |
[3] | Biwen Li, Xuan Cheng . Synchronization analysis of coupled fractional-order neural networks with time-varying delays. Mathematical Biosciences and Engineering, 2023, 20(8): 14846-14865. doi: 10.3934/mbe.2023665 |
[4] | Zhen Yang, Zhengqiu Zhang, Xiaoli Wang . New finite-time synchronization conditions of delayed multinonidentical coupled complex dynamical networks. Mathematical Biosciences and Engineering, 2023, 20(2): 3047-3069. doi: 10.3934/mbe.2023144 |
[5] | Na Zhang, Jianwei Xia, Tianjiao Liu, Chengyuan Yan, Xiao Wang . Dynamic event-triggered adaptive finite-time consensus control for multi-agent systems with time-varying actuator faults. Mathematical Biosciences and Engineering, 2023, 20(5): 7761-7783. doi: 10.3934/mbe.2023335 |
[6] | Ruoyu Wei, Jinde Cao . Prespecified-time bipartite synchronization of coupled reaction-diffusion memristive neural networks with competitive interactions. Mathematical Biosciences and Engineering, 2022, 19(12): 12814-12832. doi: 10.3934/mbe.2022598 |
[7] | Yanjie Hong, Wanbiao Ma . Sufficient and necessary conditions for global attractivity and stability of a class of discrete Hopfield-type neural networks with time delays. Mathematical Biosciences and Engineering, 2019, 16(5): 4936-4946. doi: 10.3934/mbe.2019249 |
[8] | Rui Ma, Jinjin Han, Li Ding . Finite-time trajectory tracking control of quadrotor UAV via adaptive RBF neural network with lumped uncertainties. Mathematical Biosciences and Engineering, 2023, 20(2): 1841-1855. doi: 10.3934/mbe.2023084 |
[9] | Hai Lin, Jingcheng Wang . Pinning control of complex networks with time-varying inner and outer coupling. Mathematical Biosciences and Engineering, 2021, 18(4): 3435-3447. doi: 10.3934/mbe.2021172 |
[10] | Zixiao Xiong, Xining Li, Ming Ye, Qimin Zhang . Finite-time stability and optimal control of an impulsive stochastic reaction-diffusion vegetation-water system driven by Lˊevy process with time-varying delay. Mathematical Biosciences and Engineering, 2021, 18(6): 8462-8498. doi: 10.3934/mbe.2021419 |
The escalating penetration of renewable energy, notably wind and solar, introduces unique complexities into power systems, particularly in frequency and voltage regulation. Current grid codes are inadequate for these emerging dynamics, necessitating significant enhancements to cope with the evolving energy landscape. The paper highlights several pivotal areas of improvement, including increased flexibility within power systems, integration of energy storage systems, expansion of ancillary services, the inclusion of grid-forming inverters, and the international harmonization of grid codes. The research underscores that conventional power system methodologies, primarily reliant on traditional power plants, fail to manage the fluctuating supply-demand dynamics of renewable energy effectively. By proposing improvements in grid codes, the research contributes towards resolving this issue. Furthermore, the paper underlines the paramount importance of international harmonization of grid codes for system interoperability, efficient operation, and exchange of best practices across diverse regions. Through its exhaustive exploration and recommendations, the study empowers policymakers, grid operators, and energy producers to advance grid code frameworks. Consequently, this facilitates renewable energy integration, ensures grid stability, and paves the way for a more sustainable energy future.
Hopfield neural network (HNN) is in the powerful class of artificial neural networks (ANNs), which has been introduced in the literature by John Hopfield in 1982. Since then, its investigation has become a worldwide focus [1,2,3,4,5]. All of those models are real-valued neural networks (NNs). However, complex-valued neural networks (CVNNs) can handle problems that cannot be handled with real-valued parallel networks [6,7]. Recently, some authors started the dynamical analysis of CVNNs. For example, in [8], Ali et al. studied the finite time stability analysis of delayed fractional-order memristive CVNNs by using the Gronwall inequality, H¨older inequality and inequality scaling skills. In [9], Zhang and Cao investigated the existence and global exponential stability of periodic solutions of neutral type CVNNs by collecting the Lyapunov functional method with the coincidence degree theory as well as graph theory. In article [21], the authors dealt with the problem stability of impulsive CVNNs with time delay.
In the functioning of ANNs, time-delays always exist in the signal transmission between neurons on account to the limited pace of signal exchange and transmission. Thus, it is one of the main reasons poor performance and instability in a system is the presence of delays. Therefore, in recent years, the analysis of the stability of delayed ANNs has been the great attention of researchers, and various results have been discussed in the literature [16,17,18,19,20,21,22].
ANNs with neutral-type delays were also very seldom discussed in the existing literature; such models are called neutral-type ANNs. They have been the subject of in-depth studies by many researchers, for instance, Guo and Du [17], concerned with the global exponential stability of periodic solution for neutral-type CVNNs.
The theory of passivity is a significant notion of automatic for the analysis and control of models whose certain input/output characteristics are established in terms of energy criteria. The notions of passivity are adapted to several scientific fields and are effective for the regulation of electrical, mechanical, and electromechanical systems present in several fields of engineering, such as robotics, power electronics, aeronautics, etc. Many results are concentrated in the research of the passivity of NN systems. In [11], Li and Zheng proved the global exponential passivity for delayed quaternion-valued memristor-based NNs. Ge et al. [12] studied the robust passivity analysis for a class of uncertain NNs subject to mixed delays. The problem of passivity analysis for uncertain bidirectional associative memory NNs in the presence of mixed delays are discussed in [13]. In [35], Khonchaiyaphum et al. studied FTP analysis of neutral-type NNs.
Many researchers have been interested in finite-time control problems of dynamical systems based on the Lyapunov theory of stability to give big attention to the asymptotic model of systems over an infinite time interval. Recently, many interested results have been published. In [14], Thuan et al. discussed the robust FTP for a fractional order NNs with uncertainties. In paper [19], Wei et al. are focused on a class of coupled quaternion-valued NNs with several delayed couplings to study fixed-time passivity.
The main objective of this manuscript is to deal with the FTP for the following neutral-type CVNNs:
{˙z(t)=−Az(t)+Bf(z(t))+Cf(z(t−τ(t)))+D˙z(t−h(t))+ω(t),u(t)=K1z(t)+K2f(z(t)). | (1.1) |
Here, we have
∘z(⋅)=[z1(⋅),z2(⋅),⋯,zn(⋅)]T∈Cn is the complex valued state vector,
∘ω(⋅)∈Cn is the disturbance input which belongs to L2[0,+∞),p
∘u(⋅)∈Cn is the control output,
∘f(z(⋅))=[f(z1(⋅)),f(z2(⋅)),⋯,f(zn(⋅))]T∈Cn is the complex-valued activation function,
∘A=diag{a1,a2,⋯,an}>0∈Rn×n is a diagonal matrix,
∘B=(bjk)∈Cn×n,C=(cjk)∈Cn×n,D=(djk)∈Cn×n are respectively the connection weight matrix and the connection weight matrix with delays,
∘K1 is a known real constant matrice with appropriate dimension,
∘K2 is a known complex matrice with appropriate dimension.
The initial condition of (1.1) is given as
z(s)=ψ(s),s∈[−ρ,0], | (1.2) |
where ρ=max{supt∈R(τ(t),supt∈R(h(t))}.
The main contributions of our article are :
∘ We use a more adequate hypothesis for the complex-valued activation functions (CVAF) considered in our model. Based on this assumption, a controllable model is formed by divised the real and imaginary parts, and we give a sufficient condition to realize FTP. This result is more general than the existing passivity results on real-valued NNs [35,36,37].
∘ A Lyapunov–Krasovskii function that support triple, four, and five integral terms is introduced, and the Wirtinger-type inequality technique and convex combination approach are espoused.
∘ A new set of sufficient conditions in terms of LMIs is derived to guarantee FTB and FTP results. These conditions can be simply found by the Matlab LMI toolbox.
The main contents of the manuscript are outlined as follows: in Section 2, we established new assumptions, definitions, and lemmas for the dynamic systems (1.1), which will be used later. New sufficient conditions on the FTB and FTP are discussed in Section 3. Two examples are presented in Section 4 to verify our results. At last, in Section 5, we end with the conclusion and perspective.
In this section, we present new necessary assumptions and some definitions and lemmas, which are utilized in the next section.
Assumption 1: The delays τ(⋅) and h(⋅) are differentiable functions satisfies the inequalities below:
0≤τ(t)≤ˉτ, ˙τ(t)≤μ,0≤h(t)≤ˉh, ˙h(t)≤h1. |
Assumption 2: The neuron activation function fk(⋅) satisfy the below Lipschitz condition:
|fk(z1)−fk(z2)|≤lk|z1−z2|,lk>0(k=1,2,⋯,n). | (2.1) |
Furthermore, by means of Assumption 2, it is clear to get
(f(z1)−f(z2))∗(f(z1)−f(z2))≤(z1−z2)∗LTL(z1−z2), | (2.2) |
where L=diag{l1,l2,⋯,ln}.
Assumption 3: The neuron activation functions f(⋅) can be divided into two parts as follows:
f(z)=fR(x,y)+ifI(x,y), |
where z=x+iy,i shows the imaginary unit, and the real imaginary parts check the below conditions:
1) The partial derivatives ∂f∂x,∂f∂y exist and are continuous.
2) There exist positive constant numbers γRRk,γRIk,γIRk,γIIk such that
∂fRk∂x≤γRRk,∂fRk∂y≤γRIk,∂fIk∂x≤γIRk,∂fIk∂y≤γIIk. |
Thus, one can obtain that for any x1,x2,y1,y2∈R,
|fRk(x1,y1)−fRk(x2,y2)|≤γRRk|x1−x2|+γRIk|y1−y2|,|fIk(x1,y1)−fIk(x2,y2)|≤γIRk|x1−x2|+γIIk|y1−y2|, |
for all k=1,2,⋯,n.
Remark 1: As we know, in many literatures [29,30,31], after the separation of the activation functions into real-imaginary type, ∂f∂x and ∂f∂y are still supposed to exist and be continuous and limited. Via these conditions, we can point that they checked the same inequalities given in hypothesis 3 by considering the mean value theorem for multivariable functions. Here, we delete these constraints of the partial derivatives and make the hypothesis into an appropriate one that imaginary and real parts of the activation functions just check the inequalities in hypothesis 3. Consequently, our study can help a broader class of CVAF.
Remark 2: By applying the modulus of a complex number to some simple inequalities, it is worth mentioning that hypothesis 2 is equivalent to hypothesis 3. Therefore, for writing convenience, we will assume that activation functions satisfy hypothesis 2 to study our main results.
Assumption 4: The neuron activation function fk(⋅) may be divided into two parts according to the complex number z as follows
fk(z)=fRk(Re(z))+ifIk(Im(z)) |
where fRk(⋅),fIk(⋅):R→R, then for any k=1,2,⋯,n, there exist constants ˇlRk,ˆlRk,ˇlIkˆlIk such that
ˇlRk≤fRk(κ1)−fRk(κ2)κ1−κ2≤ˆlRk,ˇlIk≤fIk(κ1)−fIk(κ2)κ1−κ2≤ˆlIk, | (2.3) |
where fRk(0)=0,fIk(0)=0,κ1,κ2∈R,κ1≠κ2.
For presentation convenience, we denote
ˆLR=diag{ˆlR1,⋯,ˆlRn},ˇLR=diag{ˇlR1,⋯,ˇlRn},ˆLI=diag{ˆlI1,⋯,ˆlIn},ˇLI=diag{ˇlI1,⋯,ˇlIn},L1=ˆLRˇLR,L2=ˆLIˇLI,L3=ˆLR+ˇLR2,L4=ˆLI+ˇLI2,ˉL1=diag{L1,L2},ˉL2=diag{L3,L4}. |
Remark 3: In [23,24,25], hypotheses 3 and 4 are given on the activation functions which separate between the real and imaginary parts. Moreover, we can even simply say that hypothesis 3 is quite strict and that it is a special case of hypothesis 2. This fact has already been mentioned in [26] and [27]. In addition, hypothesis 4 is also a strong constraint. For instance, if the activation functions fk(z)(k=1,2,⋯,n) checked hypothesis 4, we have:
|fRk(x1)−fRk(x2)|≤lRk|x1−x2|,|fIk(y1)−fIk(y2)|≤lIk|y1−y2| |
where z1=x1+iy1,z2=x2+iy2, lRk=max{|ˇlRk|,|ˆlRk|}, and lIk=max{|ˇlIk|,|ˆlIk|}, then one can have that
|fk(z1)−fk(z2)|=|(fRk(x1)+ifIk(y1))−(fRk(x2)+ifIk(y2))|≤|fRk(x1)−(fRk(x2)|+|fIk(y1)−fIk(y2)|≤lRk|x1−x2|+lIk|y1−y2|≤√(lRk)2+(lIk)2|z1−z2|. |
Thus, the activation function fk(z)(k=1,2,⋯,n) satisfies the Lipschitz condition of hypothesis 2. Therefore, hypotheses 3 and 4 imply hypothesis 2, that is, the condition of hypothesis 2 is a special case of hypotheses 3 and 4.
Remark 4: Indeed, when the activation function is given as follows: fk(z)=fRk(x,y)+ifIk(x,y) includes fk(z)=fRk(Re(z))+ifIk(Im(z)) as a special case. Therefore, in the following section, we will indicate that fk(z) is of the type of distinct real-imaginary activation functions. Let z(t)=x(t)+iy(t),z(t−τ(t))=zτ(t),x(t−τ(t))=xτ(t), y(t−τ(t))=yτ(t),x(t−h(t))=xh(t), y(t−h(t))=yh(t), B=BR+iBI,C=CR+iCI,f(z(t))=fR(x(t),y(t))+ifI(x(t),y(t)), f(zτ(t))=fR(xτ(t),yτ(t))+ifI(xτ(t),yτ(t)), ω(t)=ωR(t)+iωI(t), and u(t)=uR(t)+iuI(t), then model (1.1) can be divided into imaginary and real parts as the following form
{˙x(t)=−Ax(t)+BRfR(x(t),y(t))−BIfI(x(t),y(t))+CRfR(xτ(t),yτ(t))−CIfI(xτ(t),yτ(t))+DR˙xh(t)−DI˙yh(t)+ωR(t),˙y(t)=−Ay(t)+BIfR(x,y)(t))+BRfI(x(t),y(t))+CIfR(xτ(t),yτ(t))+CRfI(xτ(t),yτ(t))+DR˙yh(t)+DI˙xh(t)+ωI(t),uR(t)=K1x(t)+KR2fR(x(t),y(t))−KI2fI(x(t),y(t)),uI(t)=K1y(t)+KI2fR(x(t),y(t))+KR2fI(x(t),y(t)), | (2.4) |
The initial condition of (2.4) is given by:
{x(s)=φR(s),s∈[−ρ,0],y(s)=φI(s),s∈[−ρ,0], |
where φR(s),φI(s)∈C([−ρ,0],Rn).
Let θ=(xy), ˉf(θ)=(fR(x,y)fI(x,y)), ˉf(θτ)=(fR(xτ,yτ)fI(xτ,yτ)), ˉω=(ωRωI), ˉA=(A00A), ˉB=(BR−BIBIBR), ˉC=(CR−CICICR), ˉD=(DR−DIDIDR), ˉu=(uRuI),ˉK1=(K100K1),ˉK2=(KR1−KI1KI1KR1).
Model (2.4) can be rewritten as:
{˙θ(t)=−ˉAθ(t)+ˉBˉf(θ(t))+ˉCˉf(θτ(t))+ˉD˙θh(t)+ˉω(t),ˉu(t)=ˉK1θ(t)+ˉK2ˉf(θ(t)), | (2.5) |
with
θ(s)=ϕ(s),s∈[−ρ,0], |
where ϕ(s)=(φR(s)φI(s))∈C([−ρ,0],R2n).
The norm is defined as follows: For every ϕ(⋅)∈C([−ρ,0],R2n), ‖ϕ‖ˆh=max{‖ϕ‖,‖ϕ′‖}=max{supς∈[−ρ,0]|ϕ(θ)|,supς∈[−ρ,0]|ϕ′(ς)|}.
It is clear from (2.3) that:
ˉL−k≤ˉfk(κ1)−ˉfk(κ2)κ1−κ2≤ˉL+k, | (2.6) |
where, ˉL−=(ˇLR00ˇLI),ˉL+=(ˆLR00ˆLI).
Assumption 5: Given a positive value bR,bI, then ωR(t),ωI(t) checks
∫T10(ωR)T(t)ωR(t)dt≤bR,∫T10(ωI)T(t)ωI(t)dt≤bI,T1≥0,bR≥0,bI≥0,⇒∫T10ˉωT(t)ˉω(t)dt≤ˉb, | (2.7) |
where ˉb=(bRbI).
Definition 1 (FTB): For a known constant T1>0, system (1.1) is FTB with regard to (ˉc1,ˉc2,T1,ˉL,ˉb), and ˉω(⋅) checked (2.7) if: supt0∈[−ˉτ,0]{θT(t0)ˉLθ(t0),˙θT(t0)ˉL˙θ(t0)}≤ˉc1⇒θT(t)ˉLθ(t)≤ˉc2, for t∈[0,T1], where ˉc1=(c11c12),ˉc2=(c21c22), 0<c11<c21,0<c12<c22, ˉL>0 is a matrix.
Definition 2(FTP): For a known constant T1>0, system (1.1) is said to be FTP, with regard to (ˉc1,ˉc2,T1,ˉL,ˉb), if the beelow conditions are checked:
(ⅰ) Model (1.1) is FTB with regard to (ˉc1,ˉc2,T1,ˉL,ˉb).
(ⅱ) For a known constant γ>0, and since the zero initial condition, the below relation is true:
2∫T10ˉuT(t)ˉω(t)dt≥−γ∫T10ˉωT(t)ˉω(t)dt, |
when ˉω(⋅) checked (2.7).
Lemma 1 [28]: For a symmetric matrix ϑ=ϑT≥0, scalars e1>e2>0, such that the integrations below are well defined. We have
−(e1−e2)∫t−e2t−e1θT(s)ϑθ(s)ds≤−(∫t−e2t−e1θ(s)ds)Tϑ(∫t−e2t−e1θ(s)ds),−(e21−e22)2∫−e2−e1∫tt+sθT(u)ϑθ(u)duds≤−(∫−e2−e1∫tt+sθ(u)duds)Tϑ(∫−e2−e1∫tt+sθ(u)duds),−(e31−e32)6∫−e2−e1∫0s∫tt+uθT(v)ϑθ(v)dvduds≤−(∫−e2−e1∫0s∫tt+uθ(v)dvdudλds)Tθ×(∫−e2−e1∫0s∫tt+uθ(v)dvdudλds). | (2.8) |
Lemma 2: For a symmetric matrix ϑ=ϑT≥0, scalars e1>e2>0, such that the integrations below are well defined. We have
−(e41−e42)24∫−e2−e1∫0s∫0λ∫tt+uθT(v)ϑθ(v)dvdudλds≤−(∫−e2−e1∫0s∫0λ∫tt+uθ(v)dvdudλds)Tϑ×(∫−e2−e1∫0s∫0λ∫tt+uθ(v)dvdudλds). | (2.9) |
Proof: The proof of Lemma 2 is inspired by the proof of Lemma 2 in [28].
For any symmetric matrix ϑ=ϑT≥0, we have
[θT(v)ϑθ(v)θT(v)θ(v)ϑ−1]≥0, |
then after integration of it from t+u to t, from λ to 0, from s to 0, and from −e1 to −e2 in turn, we can obtain
(Π11Π12ΠT12Π22)≥0, |
where
Π11=∫−e2−e1∫0s∫0λ∫tt+uθT(v)ϑθ(v)dvdudλds,Π12=∫−e2−e1∫0s∫0λ∫tt+uθT(v)dvdudλds,Π22=∫−e2−e1∫0s∫0λ∫tt+uϑ−1dvdudλds=∫−e2−e1∫0s∫0λ−uϑ−1dudλds=∫−e2−e1∫0s[−u22ϑ−1]0λdλds=∫−e2−e1[λ36ϑ−1]0sds=[−s424M−1]−e2−e1=e41−e4224ϑ−1. |
According to Schur complement, [16] is equivalent to the below condition: Π22>0,Π11−Π12Π−122ΠT12≥0. ⇒Π11≥Π12Π−122ΠT12, which is equivalent to inequality (2.9).
This completes the proof.
Lemma 3 [38]: For a given symmetric definite matrix ξ>0 and for any differentiable function ζ(⋅):[e1,e2]→Rn, the below inequality holds:
∫e2e1˙ζT(s)ξ˙ζ(s)ds≥1e2−e1(ζ(e2)ζ(e1)χ)T×Ξ2(ξ)(ζ(e2)ζ(e1)χ), |
where χ=1e2−e1∫e2e1ζ(s)ds, Ξ2(ξ)=(ξ−ξ0⋆ξ0000)+π24(ξξ−2ξ⋆ξ−2ξ⋆⋆4ξ).
Remark 5: In this manuscript, our goal is to know how to determine the sufficient condition to study the FTB and FTP of the proposed model. So, to get the desired objectives, a processable whole model (2.5) is first formed by dividing the initial model into imaginary and real parts and founding an equivalent real-valued model. Second, using the Lyapunov function approach, the design procedure can be easily performed by checked the LMIs, so in the next section the expected conditions will be obtained.
In this section, we will concentrate on the problem of FTB and FTP.
In this first part, our goal is to study the FTB in the following system:
˙θ(t)=−ˉAθ(t)+ˉBˉf(θ(t))+ˉCˉf(θτ(t))+ˉD˙θh(t)+ˉω(t). | (3.1) |
Theorem 1: Suppose that Assumptions 1–5 hold. Let ˉτ,μ,ˉh,h1, and δ be positive scalars, then system (3.1) is FTB, with respect to (ˉc1, ˉc2, T1, ˉL, ˉb), if there exists symmetric positive definite matrices P, Q1, Q2, R01, R1, R02, R2, R3, R4 R5 R6, T01, T2, T3 and diagonal matrices M1>0, M2>0, M3>0, such that the following LMIs hold:
Ω=[η1,1ˉL1M3η1,3η1,4η1,5G1Dη1,7η1,80η1,11Π22ˉhR6ˉτ22T2ˉτ36T3G1⋆η2,20000−MT3FT2η2,8000000⋆⋆η3,3000000Π22ˉτR40000⋆⋆⋆η4,4000000Π22ˉhR6000⋆⋆⋆⋆η5,5G2ˉDη5,7G2ˉC00000G2⋆⋆⋆⋆⋆η6,600000000⋆⋆⋆⋆⋆⋆η7,7M3000000⋆⋆⋆⋆⋆⋆⋆η8,8000000⋆⋆⋆⋆⋆⋆⋆⋆η9,900000⋆⋆⋆⋆⋆⋆⋆⋆⋆η10,100000⋆⋆⋆⋆⋆⋆⋆⋆⋆⋆η11,11000⋆⋆⋆⋆⋆⋆⋆⋆⋆⋆⋆−T200⋆⋆⋆⋆⋆⋆⋆⋆⋆⋆⋆⋆−T30⋆⋆⋆⋆⋆⋆⋆⋆⋆⋆⋆⋆⋆−δI]<0, | (3.2) |
ˉc1Γ+ˉb(1−e−δT1)<ˉc2λ1e−δT1, | (3.3) |
where
η1,1=R01+R1−R4−Π24R4+ˉτ2R5−R6−Π24R6−ˉτ2T01−ˉτ44T2−ˉτ66T3−ˉL1M1−ˉL1M3−G1ˉA−ˉATGT1−δP,η1,3=R4−π24R4,η1,4=R6−Π24R6,η1,5=P−ˉL−Q1+ˉL+Q2−G1−ˉATGT2,η1,7=ˉL2M1+ˉL2M3+G1ˉB,η1,11=Π22ˉτR4+ˉτT01,η2,2=−(1−μ)R01−ˉL1M2−ˉL1M3,η2,8=ˉL2M2+MT3LT2,η3,3=−R4−π24R4−R1,η4,4=−R6−Π24R6,η5,5=R3+ˉτ2R4+ˉh2R6+ˉτ44T01+ˉτ636T2+ˉτ8576T3−G2−GT2,η5,7=QT1−QT2+G2ˉB,η6,6=−(1−h1)R3,η7,7=R02+R2−M1−M3,η8,8=−(1−μ)R02−M2−M3,η9,9=−R2,η10,10=−Π2ˉτ2R4−R5−T01,η11,11=−Π2ˉh2R6. |
Proof: We consider the Lyapunov functional below:
V(θ(t))=5∑i=1Vi(θ(t)), | (3.4) |
where
V1(θ(t))=θT(t)Pθ(t),V2(θ(t))=2∫θ(t)0Q1(ˉf(s)−ˉL−s)+Q2(ˉL+s−ˉf(s))ds,V3(θ(t))=∫tt−τ(t)θT(s)R01θ(s)ds+∫tt−τθT(s)R1θ(s)ds+∫tt−τ(t)ˉfT(θ(s))R02ˉf(θ(s))ds+∫tt−ˉτˉfT(θ(s))R2ˉf(θ(s))ds,V4(θ(t))=∫tt−h(t)˙θT(s)R3˙θ(s)ds+ˉτ∫0−ˉˉτ∫tt+β˙θT(s)R4˙θ(s)dsdβ+ˉτ∫0−ˉτ∫tt+βθT(s)R5θ(s)dsdβ+ˉh∫0−ˉh∫tt+β˙θT(s)R6˙θ(s)dsdβ,V5(θ(t))=ˉτ22∫0−ˉτ∫0γ∫tt+β˙θT(s)T01˙θ(s)dsdβdγ+ˉτ36∫0−ˉτ∫0λ∫0γ∫tt+β˙θT(s)T2˙θ(s)dsdβdγdλ+ˉτ424∫0−ˉτ∫0λ∫0α∫0γ∫tt+β˙θT(s)T3˙θ(s)dsdβdγdαdλ. |
Calculating the time-derivative of V(θ(⋅)) along any trajectory of model (3.1), we obtain
˙V(θ(t))=5∑i=1˙Vi(θ(t)), | (3.5) |
where
˙V1(θ(t))=2θT(t)P˙θ(t),˙V2(θ(t))=2(ˉf(θ(t))−ˉL−θ(t))TQ1˙θ(t)+2(ˉL+θ(t)−ˉf(θ(t)))TQ2˙θ(t),=2ˉf(θ(t))TQ1˙θ(t)−2θ(t)TˉL−Q1˙θ(t)+2θT(t)ˉL+Q2˙θ(t)−2ˉf(θ(t))TQ2˙θ(t),˙V3(θ(t))≤θT(t)[R01+R1]θ(t)+ˉfT(θ(t))[R02+R2]ˉf(θ(t))−(1−μ)ˉfT(θ(t−τ(t)))R02ˉf(θ(t−τ(t)))−θT(t−ˉτ)R1θ(t−ˉτ)−(1−μ)θT(t−τ(t))R01θ(t−τ(t))−ˉfT(θ(t−ˉτ))R2ˉf(θ(t−ˉτ)), |
˙V4(θ(t))=˙θT(t)R3˙θ(t)−(1−˙h(t))˙θT(t−h(t))R3˙θ(t−h(t))+ˉτ∫0−ˉτ˙θT(t)R4˙θ(t)−˙θT(t+β)R4˙θ(t+β)dβ+ˉτ∫0−ˉτθT(t)R5θ(t)−θT(t+β)R5θ(t+β)dβ+ˉh∫0−ˉh˙θT(t)R6˙θ(t)−˙θT(t+β)R6˙θ(t+β)dβ≤˙θT(t)R3˙θ(t)−(1−h1)˙θT(t−h(t))R3˙θ(t−h(t))+ˉτ2˙θT(t)R4˙θ(t)−ˉτ∫tt−ˉτ˙θT(s)R4˙θ(s)ds+ˉτ2θT(t)R5θ(t)−ˉτ∫tt−ˉτθT(s)R5θ(s)ds+ˉh2˙θT(t)R6˙θ(t)−ˉh∫tt−ˉh˙θT(s)R6˙θ(s)ds. |
Applying Lemma 3 to the above inequality, we get
−ˉτ∫tt−ˉτ˙θT(s)R4˙θ(s)ds≤−(θ(t)θ(t−ˉτ)1ˉτ∫tt−ˉτθ(s)ds)T×Ξ2(R4)(θ(t)θ(t−ˉτ)1ˉτ∫tt−ˉτθ(s)ds), |
where
Ξ2(R4)=(R4−R40⋆R40000)+Π24(R4R4−2R4⋆R4−2R4⋆⋆4R4), |
−ˉτ∫tt−ˉτ˙θT(s)R4˙θ(s)ds≤−[θT(t)(R4+Π24R4)θ(t)+2θT(t)(−R4+Π24R4)θ(t−ˉτ)+2θT(t)(−Π22ˉτR4)∫tt−ˉτθ(s)ds+θT(t−ˉτ)(R4+Π24R4)θ(t−ˉτ)+2θT(t−ˉτ)(−Π22ˉτR4)∫tt−ˉτθ(s)ds+(∫tt−ˉτθT(s)ds)Π2ˉτ2R4∫tt−ˉτθ(s)ds]. |
Similary,
−ˉh∫tt−ˉh˙θT(s)R6˙θ(s)ds≤−[θT(t)(R6+Π24R6)θ(t)+2θT(t)(−R6+Π24R6)θ(t−ˉh)+2θT(t)(−Π22ˉhR6)∫tt−ˉhθ(s)ds+θT(t−ˉh)(R6+Π24R6)θ(t−ˉh)+2θT(t−ˉh)(−Π22ˉhR6)∫tt−ˉhθ(s)ds+∫tt−ˉhθT(s)ds(Π2ˉh2R6)∫tt−ˉhθ(s)ds], |
and by applying Lemma 1,
−ˉτ∫tt−ˉτθT(s)R5θ(s)ds≤−(∫tt−ˉτθT(s)ds)TR5(∫tt−ˉτθ(s)ds). |
˙V5(θ(t))=ˉτ44˙θT(t)T01˙θ(t)−ˉτ22∫0−ˉτ∫tt+γ˙θT(s)T01˙θ(s)dsdγ+ˉτ636˙θT(t)T2˙θ(t)−ˉτ36∫0−ˉτ∫0λ∫tt+γ˙θT(s)T2˙θ(s)dsdγdλ+ˉτ8576˙θT(t)T3˙θ(t)−ˉτ424∫0−ˉτ∫0λ∫0α∫tt+γ˙θT(s)T3˙θ(s)dsdγdαdλ. |
Applying Lemma 1 and 3 to the above inequality, we get
A1=−ˉτ22∫0−ˉτ∫tt+γ˙θT(s)T01˙θ(s)dsdγA1≤−(∫0−ˉτ∫tt+γ˙θT(s)ds)T01(∫0−ˉτ∫tt+γ˙θ(s)dsdγ)=−[ˉτθT(t)−∫tt−ˉτθT(s)ds]T1[ˉτθ(t)−∫tt−ˉτθ(s)ds]=−ˉτ2θT(t)T01θ(t)+2ˉτθT(t)T01∫tt−ˉτθ(s)ds−(∫tt−ˉτθT(s)ds)T1(∫tt−ˉτθ(s)ds).A2=−ˉτ36∫0−ˉτ∫0λ∫tt+γ˙θT(s)T2˙θ(s)dsdγdλ≤(∫0−ˉτ∫0λ∫tt+γ˙θT(s)dsdγdλ)×T2(∫0−ˉτ∫0λ∫tt+γ˙θ(s)dsdγdλ)=−[ˉτ22θT(t)−∫0−ˉτ∫tt+λθT(s)dsdλ]T2[ˉτ22θ(t)−∫0−ˉτ∫tt+λθ(s)dsdλ]=−ˉτ44θT(t)T2θ(t)+2θT(t)(ˉτ22T2)∫0−ˉτ∫tt+λθ(s)dsdλ−(∫0−ˉτ∫tt+λθT(s)dsdλ)T2(∫0−ˉτ∫tt+λθ(s)dsdλ). |
A3=−ˉτ424∫0−ˉτ∫0λ∫0α∫tt+γ˙θT(s)T3˙θ(s)dsdγdαdλA3≤−(∫0−ˉτ∫0λ∫0α∫tt+γ˙θT(s)dsdγdαdλ)T3(∫0−ˉτ∫0λ∫0α∫tt+γ˙θ(s)dsdγdαdλ)=−[ˉτ36θT(t)−∫0−ˉτ∫0λ∫tt+αθT(s)dsdαdλ]T3[ˉτ36θ(t)−∫0−ˉτ∫0λ∫tt+αθ(s)dsdαdλ]=−ˉτ636θT(t)T3θ(t)+2θT(t)(ˉτ36T3)(∫0−ˉτ∫0λ∫tt+αθ(s)dsdαdλ)−(∫0−ˉτ∫0λ∫tt+αθT(s)dsdαdλ)T3(∫0−ˉτ∫0λ∫tt+αθ(s)dsdαdλ). | (3.6) |
For any ρ1k>0,ρ2k>0,ρ3k>0,k=1,2,⋯,n, it follows from (2.3) that
[ˉfk(θk(t))−ˉL−kθk(t)]Tρ1k[ˉL+kθk(t)−ˉfk(θk(t))]≥0,[ˉfk(θk(t−τ(t)))−ˉL−kθk(t−τ(t))]Tρ2k[ˉL+kθk(t−τ(t))−ˉfk(θk(t−τ(t)))]≥0,[ˉfk(θk(t))−ˉfk(θk(t−τ(t)))−ˉL−k(θk(t)−θk(t−τ(t)))]Tρ3k[ˉL+k(θk(t)−θk(t−τ(t)))−(ˉfk(θk(t))+ˉfk(θk(t−τ(t))))]≥0, |
which implies
(θT(t)ˉfT(θ(t)))(−ˉL1M1ˉL2M1⋆−M1)×(θ(t)ˉf(θ(t)))≥0, | (3.7) |
(θT(t−τ(t))ˉfT(θ(t−τ(t))))(−ˉL1M2ˉL2M2⋆−M2)×(θ(t−τ(t))ˉf(θ(t−τ(t))))≥0, | (3.8) |
and
(θT(t)ˉfT(θ(t))θT(t−τ(t))ˉfT(θ(t−τ(t))))×(−ˉL1M3ˉL2M3ˉL1M3−ˉL1M3⋆−M3−ˉL2M3ˉL3⋆⋆−ˉL1M3ˉL1M3⋆⋆⋆−M3)(θ(t)ˉf(θ(t))θ(t−τ(t))ˉf(θ(t−τ(t))))≥0, | (3.9) |
where M1=diag{ρ11,ρ12,⋯ρ1n},M2=diag{ρ21,ρ22,⋯ρ2n}, M3=diag{ρ31,ρ32,⋯ρ3n}.
Moreover, for any matrices G1 and G2 with appropriate dimensions, it is true that,
2[θT(t)G1+˙θT(t)G2][−˙θ(t)−ˉAθ(t)+ˉBˉf(θ(t))+ˉCˉf(θ(t−τ(t)))+ˉD˙θ(t−h(t))+ˉω(t)]=0=−2θT(t)G1˙θ(t)−2θT(t)G1ˉAθ(t)+2θT(t)G1ˉBˉf(θ(t))+2θT(t)G1ˉCˉf(θ(t−τ(t)))+2θT(t)G1ˉD˙θ(t−h(t))+2θT(t)G1ˉω(t)−2˙θT(t)G2˙θ(t)−2˙θT(t)G2ˉAθ(t)+2˙θT(t)G2ˉBˉf(θ(t))+2˙θT(t)G2ˉCf(θ(t−τ(t)))+2˙θT(t)G2ˉD˙θ(t−h(t))+2˙θT(t)G2ˉω(t). | (3.10) |
We can say that
˙V(θ(t))−δV1(θ(t))−δˉωT(t)ˉω(t)≤ξT(t)Ωξ(t)<0, | (3.11) |
where
ξ(t)=[θT(t)θT(t−τ(t))θT(t−ˉτ)θT(t−ˉh)˙θT(t)˙θT(t−h(t))ˉfT(θ(t))ˉfT(θ(t−τ(t)))ˉfT(θ(t−ˉτ))∫tt−ˉτθT(s)ds∫tt−ˉhθT(s)ds∫0−ˉτ∫tt+λθT(s)dsdλ∫0−ˉτ∫0λ∫tt+αθT(s)dsdαdλˉω(t)]T, |
and Ω is given in (3.2).
˙V(θ(t))≤δV1(θ(t))+δˉωT(t)ˉω(t)<δV(θ(t))+δˉωT(t)ˉω(t) | (3.12) |
Multiplying by e−δt, we can obtain
e−δt˙V(θ(t))−δe−δtV(θ(t))<δe−δtˉωT(t)ˉω(t),ddt(e−δtV(θ(t)))<δe−δtˉωT(t)ˉω(t). | (3.13) |
By integrating (3.13] between 0 to t, such as t∈[0,T1], we can write
e−δtV(θ(t))−V(x(0))<δ∫t0e−δsˉωT(s)ˉω(s)ds,V(θ(t))<eδt[V(θ(0))+∫t0e−δsˉωT(s)ω(s)ds]. | (3.14) |
So,
V(θ(0))=θT(0)Pθ(0)+2∫θ(0)0Q1(ˉf(s)−ˉL−s)+Q2(ˉL+s−ˉf(s))ds+∫0−τ(t)θT(s)R01θ(s)ds+∫0−ˉτθT(s)R1θ(s)ds+∫0−ˉτˉfT(θ(s))R2ˉf(θ(s))ds+∫0−τ(t)ˉfT(θ(s))R02ˉf(θ(s))ds+∫0−h(t)˙θT(s)R3˙θ(s)ds+ˉτ∫0−ˉτ∫tβ˙θT(s)R4˙θ(s)dsdβ+ˉτ∫0−ˉτ∫0βθT(s)R5θ(s)dsdβ+ˉh∫0−ˉh∫0β˙θT(s)R6˙θ(s)dsdβ+ˉτ22∫0−ˉτ∫0γ∫0β˙θT(s)T01˙θ(s)dsdβdγ+ˉτ36∫0−ˉτ∫0λ∫0γ∫0β˙θT(s)T2˙θ(s)dsdβdγdλ+ˉτ424∫0−ˉτ∫0λ∫0α∫0γ∫0β˙θT(s)T3˙θ(s)dsdβdγdαdλ. |
Letting
ˆP=ˉL−12PˉL−12,^Q1=ˉL−12Q1ˉL−12,^Q2=ˉL−12Q2ˉL−12,^R01=ˉL−12R01ˉL−12,^R1=ˉL−12R1ˉL−12,^R02=ˉL−12R02ˉL−12,^R2=ˉL−12R2ˉL−12,^R3=ˉL−12R3ˉL−12,^R4=ˉL−12R4ˉL−12,^R5=ˉL−12R5ˉL−12,^R6=ˉL−12R6ˉL−12,ˆT01=ˉL−12T01ˉL−12,ˆT2=ˉL−12T2ˉL−12,ˆT3=ˉL−12T3ˉL−12,λ1=λmin(ˆP),λ2=λmax(ˆP),λ3=λmax(ˆQ1),λ4=λmax(ˆQ2),λ5=λmax(ˆR01),λ6=λmax(ˆR1),λ7=λmax(ˆR02),λ8=λmax(ˆR2),λ9=λmax(ˆR3),λ10=λmax(ˆR4),λ11=λmax(ˆR5),λ12=λmax(ˆR6),λ13=λmax(ˆT01),λ14=λmax(ˆT2),λ15=λmax(ˆT3), |
we obtain
V(θ(0))=θT(0)ˉL12ˆPˉL12θ(0)+2∫θ(0)0ˉL12ˆQ1ˉL12(ˉf(s)−ˉL−s)+ˉL12ˆQ2ˉL12(ˉL+s−ˉf(s))ds+∫0−τ(t)θT(s)ˉL12ˆR01ˉL12θ(s)ds+∫0−ˉτθT(s)ˉL12ˆR1ˉL12θ(s)ds+∫0−ˉτ((ˉL+)2θT(s)ˉL12ˆR2ˉL12θ(s)ds+∫0−τ(t)((ˉL+)2θT(s)ˉL12ˆR02ˉL12θ(s)ds+∫0−h(t)˙θT(s)ˉL12ˆR3ˉL12˙θ(s)ds+ˉτ∫0−ˉτ∫0β˙θT(s)ˉL12ˆR4ˉL12˙θ(s)dsdβ+ˉτ∫0−ˉτ∫0βθT(s)ˉL12ˆR5ˉL12θ(s)dsdβ+ˉh∫0−ˉh∫0β˙θT(s)ˉL12ˆR6ˉL12˙θ(s)dsdβ+ˉτ22∫0−ˉτ∫0γ∫0β˙θT(s)ˉL12ˆT1ˉL12˙θ(s)dsdβdγ+ˉτ36∫0−ˉτ∫0λ∫0γ∫0β˙θT(s)ˉL12ˆT2ˉL12˙θ(s)dsdβdγdλ+ˉτ424∫0−ˉτ∫0λ∫0α∫0γ∫0β˙θT(s)ˉL12ˆT3ˉL12˙θ(s)dsdβdγdαdλ≤{λmax(ˆP)θT(0)ˉLθ(0)+2λmax(^Q1)[max{|ˉL+,ˉL−|2−ˉL−}]+2λmax(^Q2)[max{ˉL+−|ˉL+,ˉL−|2}]+ˉτ[λmax(^R01)+λmax(^R1)]+ˉτ[max{|ˉL+,ˉL−|2}][λmax(^R02)+λmax(^R2)]+ˉhλmax(^R3)+ˉτ32λmax(^R4)+ˉτ32λmax(^R5)+ˉh32λmax(^R6)+ˉτ512λmax(^T1)+ˉτ7144λmax(^T2)+ˉτ92880λmax(^T3)}×supt0∈[−ˉτ,0]{θT(t0)ˉLθ(t0),˙θT(t0)ˉL˙θ(t0)},V(θ(0))≤ˉc1Γ, |
where
Γ=[λ2+λ3[˜L2−ˉL−]+λ4[ˉL+−˜L2]+ˉτ[λ5+λ6]+ˉτ˜L2[λ7+λ8]+ˉhλ9+ˉτ32[λ10+λ11]+ˉh32λ12+ˉτ512λ13+ˉτ7144λ14+ˉτ92880λ15, |
where ˜L=max{|ˉL+,ˉL−|2.
Furthermore, it follows from (3.14) that
V(θ(t))≥θT(t)Pθ(t)≥λmin(ˆP)θT(t)ˉLθ(t)=λ1θT(t)ˉLθ(t). | (3.15) |
Because of inequalities (3.14) and (3.15), we obtain
λ1θT(t)Pθ(t)≤eδt[V(θ(0))+∫t0e−δtˉωT(s)ˉω(s)ds],⇒θT(t)ˉLθ(t)≤eδT1[ˉc1Γ+b(1−e(−δT1))]λ1. | (3.16) |
From condition (3.3), we arrive at θT(t)Lθ(t)<ˉc2. From Definition 1, the model (3.1) is FTB with regard to (ˉc1,ˉc2,T1,ˉL,ˉb).
This allowed the proof to be obtained.
Remark 6: In Theorem 1, sufficient conditions are met to verify that the model (3.1) is FTB, then Theorem 2 will present FTP conditions.
In this second part, we study the FTP analysis for the below model
{˙θ(t)=−ˉAθ(t)+ˉBˉf(θ(t))+ˉCˉf(θ(t−τ(t)))+ˉD˙θ(t−h(t))+ˉω(t),ˉu(t)=ˉK1θ(t)+ˉK2ˉf(θ(t)). | (3.17) |
Theorem 2: Suppose that Assumptions 1–5 hold. Let ˉτ,μ,ˉh,h1, and δ be scalars, then system (3.17) is FTP with keeping the parameter (ˉc1, ˉc2, T1, ˉL, ˉb), if there exists symmetric positive definite matrices P, Q1, Q2, R01, R1, R02, R2, R3, R4, R5, T01, T2, T3 and diagonal matrices M1>0, M2>0, M3>0, and scalar β>0 such that the following LMIs (3.18) hold:
˜Ω=[η1,1ˉL1M3η1,3η1,4η1,5G1Dη1,7η1,80η1,10Π22hR6ˉτ22T2ˉτ36T3G1−βI−ˉK1⋆η2,20000−MT3ˉLT2η2,8000000⋆⋆η3,3000000Π22ˉτR40000⋆⋆⋆η4,4000000Π22hR6000⋆⋆⋆⋆η5,5G2ˉDη5,7G2ˉC00000G2⋆⋆⋆⋆⋆η6,600000000⋆⋆⋆⋆⋆⋆η7,7M300000−ˉK2⋆⋆⋆⋆⋆⋆⋆η8,8000000⋆⋆⋆⋆⋆⋆⋆⋆η9,900000⋆⋆⋆⋆⋆⋆⋆⋆⋆η10,100000⋆⋆⋆⋆⋆⋆⋆⋆⋆⋆η11,11000⋆⋆⋆⋆⋆⋆⋆⋆⋆⋆⋆−T200⋆⋆⋆⋆⋆⋆⋆⋆⋆⋆⋆⋆−T30⋆⋆⋆⋆⋆⋆⋆⋆⋆⋆⋆⋆⋆−βI]<0, | (3.18) |
ˉc1Γ+ˉb(1−e−δT1)<ˉc2λ1e−δT1, |
where the parameters are kept the same as Theorem 3.1.
Proof: Select the Lyapunov function used in the Theorem 3.1. We get
˙V(θ(t))−δV1(θ(t))−2ˉuT(t)ˉω(t)−βˉωT(t)ˉω(t)≤ξT(t)˜Ωξ(t)<0,⇒˙V(θ(t))−δV(θ(t))−2ˉuT(t)ˉω(t)−βˉωT(t)ˉω(t)<0, |
where ˜Ω is given in (3.18), then
˙V(θ(t))−δV(θ(t))≤2ˉuT(t)ˉω(t)+βˉωT(t)ˉω(t). | (3.19) |
Multiplying (3.19) by e−δt, we can get
e−δt˙V(θ(t))−δe−δtV(θ(t))≤2e−δtˉuT(t)ˉω(t)+βe−δtˉωT(t)ˉω(t), |
ddt(e−δtV(θ(t)))≤2e−δtˉuT(t)ˉω(t)+βe−δtˉωT(t)ˉω(t). |
Integrating the above inequality between 0 to T1, we can write
e−δT1V(θ(t))−V(θ(0))≤2∫T10e−δtˉuT(t)ˉω(t)dt+β∫T10e−δtˉωT(t)ˉω(t)dt. |
Consider the zero initial condition for θ0=0 and we have V(θ(0))=0, then
e−δT1V(θ(t))≤2∫T10e−δtˉuT(t)ˉω(t)dt+β∫T10e−δtˉωT(t)ˉω(t)dt, |
V(θ(t))≤eδT1∫T10[2ˉuT(t)ˉω(t)+βˉωT(t)ˉω(t)]dt. | (3.20) |
Since V(θ(t))≥0, we can say from (3.20) that
2∫T10ˉuT(t)ˉω(t)dt≥−β∫T10ˉωT(t)ˉω(t)dt. | (3.21) |
Finally, we can say that the system (3.17) is FTP. This allowed the proof to be obtained.
Remark 7: It should be emphasized that the LMI approach proposed in this paper is more useful for reducing the conservatism of the delay system, which may lead to obtain less conservative results. This proves the advantage of our proposed method.
Remark 8: We note that the optimal value of ˉc2 depends on the parameter δ, then the optimal minimum value of ˉc2 is determined from the minimum value of δ, such that the LMIs matrix solution remains feasible.
Remark 9: In the majority of published work, the Lyapunov function theory is the most effective approach to investigating the problem of stability and passivity for divers dynamical models. Moreover, it can be seen that the existing literature [11,12,13,15,19,32,33,34,35] contains the Lyapunov function with single, double, tripe, and four integral terms. However, in this article, we give the Lyapunov–Krasovskii function with five integral terms such as ˉτ424∫0−ˉτ∫0λ∫0α∫0γ∫tt+β˙θT(s)T3˙θ(s)dsdβdγdαdλ. Different from the published work, this is the first time studying the problem of FTP of neutral-type CVNNs. Via a Lyapunov–Krasovskii function with triple, four and five integral terms, by utilizing Jensons inequality and the Wirtinger-type inequality technique, new sufficient conditions for FTB and FTP are taken in terms of LMIs, which is effective on reducing conservatism.
In this section, two examples are given to show the feasibility of our results.
Example 1: Consider the following model
{˙z(t)=−Az(t)+Bf(z(t))+Cf(z(t−τ(t)))+D˙z(t−h(t))+ω(t),z(s)=ψ(s),s∈[−ρ,0], | (4.1) |
where
A=(1.9001.2),B=(0.5+0.5i0.03+0.03i−0.4−0.4i−0.1−0.1i),C=(1.1+1.1i0.03+0.03i0.07+0.07i0.1+0.1i),D=(0.1+0.1i000.1+0.1i),ˉA=(1.900001.200001.900001.2),ˉB=(0.50.03−0.5−0.03−0.4−0.10.40.10.50.030.50.03−0.4−0.1−0.4−0.1),ˉC=(1.10.03−1.1−0.030.070.1−0.07−0.11.10.031.10.030.070.10.070.1),ˉD=(0.10−0.1000.10−0.10.100.1000.100.1), |
fj(z)=0.3(|xj+1|−|xj−1|)+i0.3(|yj+1|−|yj−1|). Thus, L1=L2=02×2 and L3=L4=(0.3000.3),⇒ ˉL1=04×4 and ˉL2=0.3I4 τ(t)=0.15(1−sin(2t)),h(t)=0.15(1−cos(2t)),ω(t)=0.9sin(πt)e−0.5t+sin(πt)e−0.3ti, δ=1,ˉc1=(0.50.5)0.5,ˉb=(0.810.60),T1=15, matrix ˉL=I. Using the Matlab LMI toolbox to solve the LMIs (3.2) and (3.3) in Theorem 1, we obtained the feasible solutions for an optimal minimum value of ˉc2=(9.82089.4072). The trajectories of the solution of system (4.1) in Example 1 are shown in Figures 1 and 2, and the time history of (Re(z))TL(Re(z)) and (Im(z))TL(Im(z)) are given in Figure 3. Hence, it can be concluded that the system (4.1) is FTB.
Example 2: Consider the following system :
{˙z(t)=−Az(t)+Bf(z(t))+Cf(z(t−τ(t)))+D˙z(t−h(t))+ω(t),u(t)=K1z(t)+K2f(z(t)),z(s)=ψ(s),s∈[−ρ,0], | (4.2) |
where
A=(3.8002.4),B=(1+i0.06+0.06i−0.8−0.8i−0.2−0.2i),C=(2.2+2.2i0.06+0.06i0.14+0.14i0.2+0.2i),D=(0.2+0.2i000.2+0.2i),K1=(1.21.61.251),K2=(1+i0.5+0.6i0.5−0.4i1−i),ˉA=(3.800002.400003.800002.4),ˉB=(10.06−1−0.06−0.8−0.20.80.210.0610.06−0.8−0.2−0.8−0.2),ˉC=(2.20.06−2.2−0.060.140.2−0.14−0.22.20.062.20.06−0.140.20.140.2),ˉD=(0.20−0.2000.20−0.20.200.2000.200.2),ˉK1=(1.2−1.600−1.25100001.2−1.600−1.251),ˉK2=(10.51−0.60.510.41−10.610.5−0.4−10.51), |
fj(z)=0.6tanh(xj)+i0.7tanh(yj). Thus, L1=L2=02×2 and L3=(0.3000.3),L4=(0.350.35)⇒ ˉL1=04×4 and ˉL2=(0.300000.300000.3500000.35), τ(t)=0.15(1−cos(2t)),h(t)=0.15(1−sin(2t)),ω(t)=0.2e−0.5t+0.3e−0.3ti, δ=1,ˉc1=(0.10.1)0.5, ˉb=(0.0190.053),T1=12, matrix ˉL=I. Using the Matlab LMI toolbox to solve the LMIs (3.2) and (3.3) in Theorem 2, we obtained the feasible solutions for an optimal minimum value of ˉc2=(12.443115.0562). From Table 1, The results presented in this manuscript are significantly better than those of [32,33,34,35], which confirms the validity of our work.
μ | 0.8 | 0.9 |
[32] | 3.9212 | 2.3901 |
[33] | 5.2403 | 3.5211 |
[34] | 5.6384 | 3.7718 |
[35] | 6.5411 | 4.5074 |
ˉτ of our result | 6.9307 | 5.2113 |
The trajectories of the solution of model (4.2) in Example 4 are given in Figure 4 and Figure 5. The time history of (Re(z))TL(Re(z)) and (Im(z))TL(Im(z)) are given in Figure 6. Hence, it can be concluded that the system (4.2) is FTP.
This paper is focused on the FTP problem of neutral-type complex-valued NNs in the presence of time-varying delays. By using Lyapunov functionals, the Wirtinger inequality, and LMIs, new sufficient conditions are gotten to guarantee the finite-time boundedness and finite-time passivity of our model. Finally, two examples are presented to prove the effectiveness of our main results. In the future work, we will take the challenge to study the finite-time dissipativity of stochastic complex NNs with mixed delays via a non-separation approach and the fixed-time passivity of coupled clifford-valued NNs subject to multiple delayed couplings.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
The authors declare there is no conflict of interest.
[1] | IRENA, Grid codes for renewable powered systems. IRENA, 2022. Available from: https://www.irena.org/publications/2022/Apr/Grid-codes-for-renewable-powered-systems. |
[2] | Altın M, Göksu Ö, Teodorescu R, et al. (2010) Overview of recent grid codes for wind power integration. 2010 12th international conference on optimization of electrical and electronic equipment, IEEE, 1152–1160. https://doi.org/10.1109/OPTIM.2010.5510521 |
[3] | Central Electricity Regulatory Commission, Draft Indian Electricity Grid Code 2022. Central Electricity Regulatory Commission, 2022. Available from: https://cercind.gov.in/2022/draft_reg/Draft-IEGC-07062022.pdf. |
[4] | IRENA, Renewable capacity highlight. IRENA, 2022. Available from: https://www.irena.org/-/media/Files/IRENA/Agency/Publication/2022/Apr/IRENA_-RE_Capacity_Highlights_2022.pdf?la = en & hash = 6122BF5666A36BECD5AAA2050B011ECE255B3BC7. |
[5] | ENTSO-E, Frequency Stability in Long-Term Scenarios and Relevant Requirements. ENTSO-E, 2021. Available from: https://eepublicdownloads.azureedge.net/clean-documents/Publications/ENTSO-E general publications/211203_Long_term_frequency_stability_scenarios_for_publication.pdf. |
[6] |
Zhao F, Bai F, Liu X, et al. (2022) A review on renewable energy transition under China's carbon neutrality target. Sustainability 14: 15006. https://doi.org/10.3390/su142215006 doi: 10.3390/su142215006
![]() |
[7] | FERC, Department of Energy Federal Energy Regulatory Commission. FERC, 2018. Available from: https://ferc.gov/sites/default/files/2020-06/Order-845.pdf. |
[8] |
Robles E, Haro-Larrode M, Santos-Mugica M, et al. (2019) Comparative analysis of European grid codes relevant to offshore renewable energy installations. Renewable Sustainable Energy Rev 102: 171–185. https://doi.org/10.1016/j.rser.2018.12.002 doi: 10.1016/j.rser.2018.12.002
![]() |
[9] |
Foley AM, McIlwaine N, Morrow DJ, et al. (2020) A critical evaluation of grid stability and codes, energy storage and smart loads in power systems with wind generation. Energy 205: 117671. https://doi.org/10.1016/j.energy.2020.117671 doi: 10.1016/j.energy.2020.117671
![]() |
[10] |
Zhao X, Flynn D (2022) Stability enhancement strategies for a 100% grid-forming and grid-following converter-based Irish power system. IET Renew Power Gener 16: 125–138. https://doi.org/10.1049/rpg2.12346 doi: 10.1049/rpg2.12346
![]() |
[11] | IEA, Status of Power System Transformation 2018: Advanced Power Plant Flexibility. International Energy Agency, 2018. Available from: https://iea.blob.core.windows.net/assets/ede9f1f7-282e-4a9b-bc97-a8f07948b63c/Status_of_Power_System_Transformation_2018.pdf. |
[12] | de Jong J, Hassel A, Jansen J, et al. (2017) Improving the market for flexibility in the electricity sector. CEPS 2017: 13093. |
[13] | Enerdata, World Energy & Climate Statistics—Yearbook 2023. Enerdata, 2023. Available from: https://yearbook.enerdata.net/renewables/wind-solar-share-electricity-production.html. |
[14] | European Commission, Summary of the responses to the targeted stakeholder consultation on the priority list for the development of network codes and guidelines on electricity for the period 2020–2023 and on gas for 2020 and beyond. European Commission, 2020. Available from: https://ec.europa.eu/energy/sites/default/files/summary_for_publication_ares.pdf. |
[15] | SmartEn, Open letter: A Network Code for Distributed Flexibility. Smart Energy Europe, 2020. Available from: https://smarten.eu/wp-content/uploads/2020/04/Open-Letter-on-a-Network-Code-for-Distributed-Flexibility.pdf. |
[16] | European Commission, Recommendations on energy storage. European Commission, 2023. Available from: https://energy.ec.europa.eu/topics/research-and-technology/energy-storage/recommendations-energy-storage_en. |
[17] | ENTSOE-E, Storage Expert Group: Phase Ⅱ final report, Belgium. ENTSOE-E, 2020. Available from: https://www.entsoe.eu/network_codes/cnc/expert-groups/. |
[18] | Service Public Federal Economie P.M.E., Classes Moyennes et Energie, Arrêté royal établissant un règlement technique pour la gestion du réseau de transport de l'électricité et l'accès à celui-ci. Moniteur belgisch, 2019. Available from: https://www.creg.be/sites/default/files/assets/Publications/Advices/A1661Annex.pdf. |
[19] | Fingrid, Grid code specifications for grid energy storage systems SJV2019. Fingrid, 2019. Available from: https://www.fingrid.fi/globalassets/dokumentit/en/customers/grid-connection/grid-energy-storage-systems-sjv2019.pdf. |
[20] | OFGEM, Modification proposal: Grid Code GC0096: Energy Storage. OFGEM, 2020. Available from: https://www.ofgem.gov.uk/. |
[21] | VDE FNN, VDE-AR-N 4110: Technical requirements for the connection and operation of customer installations to the high voltage network (TAR high voltage). VDE FNN, 2018. Available from: https://www.vde.com/en/fnn/topics/technical-connection-rules/tcr-for-medium-voltage. |
[22] |
Gulotta F, Daccò E, Bosisio A, et al. (2023) Opening of ancillary service markets to distributed energy resources: A review. Energies 16: 2814. https://doi.org/10.3390/en16062814 doi: 10.3390/en16062814
![]() |
[23] | SmartNet, Ancillary service provision by RES and DSM connected at distribution level in the future power system. SmartNet, 2016. Available from: https://smartnet-project.eu/wp-content/uploads/2016/12/D1-1_20161220_V1.0.pdf. |
[24] | Arnold GW (2012) Realizing an interoperable and secure smart grid on a national scale, In: Sorokin A, Rebennack S, Pardalos PM, et al., Handbook of Networks in Power Systems I, Berlin: Springer Science & Business Media, 487–504. https://doi.org/10.1007/978-3-642-23193-3_19 |
[25] |
Singh B, Singh SN (2009) Wind power interconnection into the power system: A review of grid code requirements. Electr J 22: 54–63. https://doi.org/10.1016/j.tej.2009.04.008 doi: 10.1016/j.tej.2009.04.008
![]() |