
Citation: Jalil Manafian, Onur Alp Ilhan, Sizar Abid Mohammed. Forming localized waves of the nonlinearity of the DNA dynamics arising in oscillator-chain of Peyrard-Bishop model[J]. AIMS Mathematics, 2020, 5(3): 2461-2483. doi: 10.3934/math.2020163
[1] | Saulo Orizaga, Maurice Fabien, Michael Millard . Efficient numerical approaches with accelerated graphics processing unit (GPU) computations for Poisson problems and Cahn-Hilliard equations. AIMS Mathematics, 2024, 9(10): 27471-27496. doi: 10.3934/math.20241334 |
[2] | Mingliang Liao, Danxia Wang, Chenhui Zhang, Hongen Jia . The error analysis for the Cahn-Hilliard phase field model of two-phase incompressible flows with variable density. AIMS Mathematics, 2023, 8(12): 31158-31185. doi: 10.3934/math.20231595 |
[3] | Hyun Geun Lee . A mass conservative and energy stable scheme for the conservative Allen–Cahn type Ohta–Kawasaki model for diblock copolymers. AIMS Mathematics, 2025, 10(3): 6719-6731. doi: 10.3934/math.2025307 |
[4] | Martin Stoll, Hamdullah Yücel . Symmetric interior penalty Galerkin method for fractional-in-space phase-field equations. AIMS Mathematics, 2018, 3(1): 66-95. doi: 10.3934/Math.2018.1.66 |
[5] | Alain Miranville . The Cahn–Hilliard equation and some of its variants. AIMS Mathematics, 2017, 2(3): 479-544. doi: 10.3934/Math.2017.2.479 |
[6] | Naveed Iqbal, Mohammad Alshammari, Wajaree Weera . Numerical analysis of fractional-order nonlinear Gardner and Cahn-Hilliard equations. AIMS Mathematics, 2023, 8(3): 5574-5587. doi: 10.3934/math.2023281 |
[7] | Aymard Christbert Nimi, Daniel Moukoko . Global attractor and exponential attractor for a Parabolic system of Cahn-Hilliard with a proliferation term. AIMS Mathematics, 2020, 5(2): 1383-1399. doi: 10.3934/math.2020095 |
[8] | Joseph L. Shomberg . Well-posedness and global attractors for a non-isothermal viscous relaxationof nonlocal Cahn-Hilliard equations. AIMS Mathematics, 2016, 1(2): 102-136. doi: 10.3934/Math.2016.2.102 |
[9] | Saulo Orizaga, Ogochukwu Ifeacho, Sampson Owusu . On an efficient numerical procedure for the Functionalized Cahn-Hilliard equation. AIMS Mathematics, 2024, 9(8): 20773-20792. doi: 10.3934/math.20241010 |
[10] | Taishi Motoda . Time periodic solutions of Cahn-Hilliard systems with dynamic boundary conditions. AIMS Mathematics, 2018, 3(2): 263-287. doi: 10.3934/Math.2018.2.263 |
In the past few years, investigations of complex dynamical networks (CDNs) have drawn much attention from scholars due to their applications in real-world schemes such as the World Wide Web, electric power grids, food webs, and so on [1,2,3]. Complex networks comprise a vast multitude of interconnected nodes and edges that depict the dynamic behaviors exhibited by various individuals and capture the relationships between them. Synchronization is a common and significant collective phenomenon observed in complex networks. It refers to the gradual convergence of node trajectories toward a shared state over time, despite varying initial conditions and interactions between nodes. This convergence of node states forms the basis of synchronization analysis in complex networks [4]. There are two types of synchronization: Internal synchronization and external synchronization. External synchronization refers to synchronization between nodes in different networks while ignoring synchronization between nodes within the same CDNs. On the other hand, internal synchronization involves adjusting the movements of two or more nodes to achieve a shared behavior. Various studies have focused on internal synchronization, such as finite-time synchronization [8], cluster synchronization [9], projective synchronization [10], lag synchronization [11], and more. Therefore, investigating internal synchronization is an important issue. Dynamic systems in the real world are commonly acknowledged for their intricate, uncertain, and nonlinear characteristics, making them challenging to control or manipulate. The fuzzy logic theory is an intriguing and efficient approach to addressing the complexities of analyzing and synthesizing such complex nonlinear systems. Among the various fuzzy methods, the Takagi-Sugeno (T-S) fuzzy model [5] is an effective tool that merges principles from linear and fuzzy theories to resolve control problems in nonlinear systems. By utilizing fuzzy membership functions, these linear subsystems are integrated. Consequently, the author merged the T-S fuzzy model with complex networks to study the synchronization challenge described in the references [6,7].
It is widely recognized that in complex networks, the nodes must share information with their neighboring nodes through an information channel. Nevertheless, time delay is an unavoidable consequence caused by long-distance communication, transmission rate, and transmission blockage. This time delay has significant implications for the instability of complex networks and the decline in system performance, as highlighted in previous studies [12,13,14]. Nevertheless, these sources only considered scenarios where neighboring nodes were either fully connected or completely disconnected. However, in many realistic networks, incomplete information is commonly encountered. For instance, research [15] discovered that only 5 of all excitatory synapses in the human brain network, connecting cortical regions, can transmit information. Additionally, references [16,17] established that the transmitter cannot receive information from all users' information channels. Hence, it becomes essential to investigate and analyze the problem of partial coupling in CDNs, which has been explored in reference [18]. Scholars in [19,20,21] have examined a different form of connection called distributed coupling. The spatial nature of complex networks introduces variations in axon sizes and the length of parallel paths. It is crucial to account for the continuous distributed coupling delay when representing such networks. This consideration brings the model closer to real-world scenarios and aids in addressing practical issues more effectively. Previous research on synchronizing CDNs with time-varying delays primarily assumes that differentiable functions can represent these delays. However, this article focuses on CDNs that incorporate continuous time-varying non-differentiable delays.
When all nodes within a network cannot synchronize independently, it presents a challenge to devise a practical controller capable of achieving network synchronization. CDNs can utilize various control schemes, such as pinning control [22,23], impulsive control [17], intermittent control [24,25], and adaptive control [26,27]. Nevertheless, due to frequent perturbations in the controller gain and potential actuator deterioration, inaccuracies inevitably arise in the controller's operation. Therefore, researchers have explored nonfragile controllers, which are precise controllers capable of tolerating acceptable faults in the controller's functioning, and this area of study has been investigated by several researchers [28,29,30,31,32,33,34]. For example, in [33] a non-fragile controller was utilized to tackle the synchronization problem in discrete-time CDNs with interval time-varying delays. The investigation in reference [31] focused on the problem of achieving exponential synchronization in CDNs with time-varying coupling delays by integrating non-fragile control with sampled-data control. In a subsequent study, the authors in reference [29] presented a method for achieving non-fragile synchronization in CDNs with hybrid delays and stochastic disturbance using sampled-data control. Following that, concern [34] addressed the issue of attaining H∞ performance synchronization in a specific class of T-S fuzzy complex dynamical networks with hybrid coupling delays. The hybrid delay in this context encompassed both discrete-time coupling and distributed coupling delays. In contrast to these references, this paper adopts a combination of interval time-varying delays, which include discrete-time coupling delay, partial coupling delay, and distributed coupling delay. Consequently, the problem of researching synchronization in CDNs with hybrid time-varying delays poses a significant challenge.
However, during the actual control process, certain unavoidable uncertainties, external disturbances, or changes in control parameters may occur, leading to a decline in the performance of complex networks and even jeopardizing their stability. In such cases, the dissipativity theory, which is based on the energy properties of dynamic systems, offers valuable tools for system analysis. Dissipative systems especially posit that the energy dissipated within a dynamic system is always less than the energy supplied from an external source. Recently, researchers have applied dissipative analysis to address the impact of external disturbances in dynamical systems [35,36]. Additionally, in [4,37,38,39,40], the notion of extended dissipativity performance was introduced as a comprehensive framework that unifies various robust performance measures such as H∞, L2−L∞, passivity, and dissipativity performances. These robust performances are achieved by appropriately tuning the weighting parameters in the extended dissipativity criteria. Despite these developments, to our knowledge, there have yet to be any reports on the extended dissipativity analysis of T-S fuzzy CDNs with mixed interval time-varying delay. Inspired by earlier research, we aim to conduct an extended dissipativity analysis of the conceded CDNs framework with hybrid interval time-varying delays.
The main findings of this paper, drawn from the discussions mentioned earlier, can be summarized in the following aspects:
● For the first time, we conducted a novel study exploring the extended dissipative performance and non-fragile control concerning the synchronization of T-S fuzzy complex dynamical networks. These networks incorporate interval discrete and distributed coupling time-varying delay, which may involve non-differentiable functions and the lower bound of the delays is not limited to zero.
● Compared with the model of the complex dynamical networks with coupling delays presented in references [20,41], we opted for a T-S fuzzy complex dynamical networks model with a combination of hybrid coupling delays. This includes partial coupling delay, distributed coupling delay, and discrete-time coupling delay. Furthermore, we employed a regrouping approach to handle these coupling delays efficiently to replace the Kronecker product, thus avoiding complex calculations.
● We have developed a synchronization criterion for hybrid delayed CDNs, utilizing a non-fragile state-feedback control design to assess the extended dissipativity performance. To make a comparison with the controller designed in [34], we utilized a single free-matrix parameter F in Eq (3.17) instead of the two parameters used in the zero equation. Importantly, we were able to achieve this without introducing any additional variables.
● We can achieve a less conservative condition by creating multiple integral Lyapunov functionals and employing more stringent scaling methods, an integral inequality based on Jensen's inequality, Wirtinger's inequality, and various mathematical techniques. Furthermore, our findings can be used to validate Chua's circuit system.
The remaining sections are structured as follows: Section 2 presents the model description and some preliminary concepts. Section 3 discusses the results of the extended dissipative synchronization under the non-fragile control. Section 4 provides two simulation examples to demonstrate the feasibility and validity of our obtained results. Finally, in Section 5, we present the overall conclusion of our main work.
Notation: Rm×n and Rn correspond to the collection of all real matrices with dimensions m×n and the n-dimensional Euclidean space, respectively. The symbol XT denotes the transpose of the matrix X. The notation diag{…} indicates the formation of a block diagonal matrix. The symbol ∗ is used to denote symmetry. col{…} denotes a column vector. In a matrix X, the expression X>0 signifies that X is positive definite. sym{X} denotes X+XT. The space Ln2[0,∞) is the Hilbert space containing square-integrable functions defined over the interval [0,∞).
This paper studies T-S fuzzy complex dynamical networks involving hybrid coupling delays. Within the network, every T-S fuzzy node is equipped with dth plant rules:
Plant Rule d: IF χ1(t) is ψd1,χ2(t) is ψd2,…,χp(t) is ψdp.
Then,
˙xi(t)=Adxi(t)+Bdf(xi(t))+c1N∑j=1,j≠iFijΥdij(xj(t)−xi(t))+c2N∑j=1,j≠iFijΥdij(xj(t−γ(t))−xi(t−γ(t)))+c3N∑j=1,j≠iFijΥdij∫t−ς1(t)t−ς2(t)(xj(s)−xi(s))ds+udi(t)+ωi(t)zi(t)=ℑxi(t),i=1,2,…,N, | (2.1) |
where χ1(t),…,χp(t) are the premise variables;ψd1,…,ψdp are the fuzzy sets identified by their membership values, d∈{1,2,…,r},r denotes the number of fuzzy rules;xi(t)=col{xi1(t),xi2(t),…,xin(t)}∈Rn denotes the state vector of the ith node;zi(t)∈Rn represents output of each node;Ad=diag{ad1,ad2,…,adn},Bd=[bdij]n×n∈Rn×n and ℑ=diag{ℑ1,ℑ2,…,ℑn} are given constant matrices. The function f(xi(t))=col{f1(xi1(t)),f2(xi2(t)),…,fn(xin(t))}∈Rn is a continuous nonlinear function with a global Lipschitz condition. ci>0,(i=1,2,3) is the constant coupling strength. F=[Fij]N×N is the outer coupling matrix, if there is a connection from node j to node i, then the coupling Fij≠0(j≠i); otherwise, Fij=0 and the diagonal elements Fii=−N∑j=1,j≠iFij;Υdij=diag{Υ1dij,Υ2dij,…,Υndij} is an inner coupling matrix with Υldij=0 or 1.Υldij denotes the information transmission of channel, if the lth channel of state node j can be connected to node i, then Υldij=1; otherwise, Υldij=0. The real function γ(t) and ςi(t) (i=1,2) represent the interval discrete and distributed time-varying delays, respectively. These delays adhere to the conditions 0≤γ1≤γ(t)≤γ2 and 0≤ς1≤ς1(t)≤ς2(t)≤ς2, where γ1, γ2, ς1 and ς2 are known real constants. Using the weighted average defuzzifier, the complex dynamical networks system can be described as follows:
˙xi(t)=r∑d=1ρd(χ(t)){Adxi(t)+Bdf(xi(t))+c1N∑j=1,j≠iFijΥdij(xj(t)−xi(t))+c2N∑j=1,j≠iFijΥdij(xj(t−γ(t))−xi(t−γ(t)))+c3N∑j=1,j≠iFijΥdij∫t−ς1(t)t−ς2(t)(xj(s)−xi(s))ds+udi(t)+ωi(t)}zi(t)=r∑d=1ρd(χ(t))ℑxi(t)i=1,2,…,N, | (2.2) |
where χ(t)=col{χ1(t),χ2(t),…,χp(t)} defines the normalized membership function that satisfies the condition:
ρd(χ(t))=ηd(χ(t))r∑d=1ηd(χ(t)),ηd(χ(t))=p∏j=1ψdj(χ(t)), | (2.3) |
where ηd(χ(t)) denotes the membership grade of χi(t) in ρd. As mentioned in the fuzzy sets definition, it is presumed that ηd(χ(t))≥0,r∑d=1ρd(χ(t))=1.
Assume that s(t)=col{s1(t),s2(t),…,sn(t)} represents the state of isolated nodes and satisfies the following:
˙s(t)=r∑d=1ρd(χ(t)){Ads(t)+Bdf(s(t))}. | (2.4) |
Consider νi(t) as the output of an isolated node, defined as νi(t)=ℑsi(t). The output error ˜zi(t) can be represented as ˜zi(t)=zi(t)−νi(t) for i=1,2,…,N. Moreover, to address the issue of partial couplings, researchers have adopted the regrouping method to investigate the impact of channels on T-S fuzzy CDNs. Let Gldij=FijΥdij=diag{g1dij,g2dij,…,gndij}(i≠j) and Gldii=−N∑j=1,j≠iGldij=diag{g1dii,g2dii,…,gndii}.
The error system is defined as ei(t)=xi(t)−s(t), where ei(t)=col{ei1,ei2(t),…,ein(t)}. Subsequently, according to (2.2), it can be inferred that the error system can be represented as:
˙ei(t)=r∑d=1ρd(χ(t)){Adei(t)+Bdf(ei(t))+c1N∑j=1FijΥdijej(t)+c2N∑j=1FijΥdijej(t−γ(t))+c3N∑j=1FijΥdij∫t−ς1(t)t−ς2(t)ej(s)ds+udi(t)+ωi(t)}˜zi(t)=ℑei(t),i=1,2,…,N, | (2.5) |
where f(ei(t))=col{f1(ei1(t)),f2(ei2(t)),…,fn(ein(t))}∈Rn with f(ei(t))=f(xi(t))−f(s(t)),˜zi(t) as the output for ith node.
Based on the regrouping method, the matrix Gldij∈RN×N for the lth channel can be described as follows:
Gld=[gld11gld12⋯gld1Ngld21gld22⋯gld2N⋮⋮⋱⋮gldN1gldN2⋯gldNN],l=1,2,…,n. | (2.6) |
Let el(t)=col{e1l(t),e2l(t),…,eNl(t)} be the error state vector of lth channel, and the synchronization error system can be established as below:
˙el(t)=r∑d=1ρd(χ(t)){adlel(t)+n∑j=1bdljˆFj+c1Gldel(t)+c2Gldel(t−γ(t))+c3Gld∫t−ς1(t)t−ς2(t)el(s)ds+udl(t)+ωl(t)}˜zl(t)=ˆℑel(t),l=1,2,…,n, | (2.7) |
where udl(t)∈RN represents the control input of the lth channel, and
ˆFj=[fj(e1j(t))fj(e2j(t))⋮fj(eNj(t))]=[fj(x1j(t))−fj(s(t))fj(x2j(t))−fj(s(t))⋮fj(xNj(t))−fj(s(t))]. |
The main objective of this paper is to attain extended dissipative synchronization to the target node s(t) in T-S fuzzy CDNs with hybrid coupling delays by the non-fragile controller. The controller is defined as follows:
udl(t)=(Kdl+ΔKl)el(t), | (2.8) |
where Kdl=diag{kd1l,kd2l,…,kdNl} is the control gain matrix. Moreover, the perturbed matrix ΔK(t) can be expressed as the norm-bounded additive in the following form
ΔKl(t)=ClDl(t)El, | (2.9) |
where Cl and El are fixed matrices provided and Dl(t)∈Rn×n, where l=1,2,…,n, is a time-varying unknown matrix that satisfies the condition DTl(t)Dl(t)≤I.
Remark 1. The controller application introduces perturbations in the controller gains, which could cause system errors or adjustments made to modify the controller gains. Hence, creating a controller system that remains unaffected by these fluctuations becomes essential. The non-fragile controller embodies precision and insensitivity to permissible operational errors of the controller, a topic that has garnered significant focus in previous studies [28,29,30,31,32,33,34]. Given these aspects, the rationale for exploring the adoption of a non-fragile state feedback controller becomes apparent.
Next, the paper will utilize various definitions, assumptions and lemmas to demonstrate the proof of the theorem.
Assumption 1. [42] Each component of the nonlinear function f(xi(t)) satisfies a global Lipschitz condition, acquired through the following way:
(fj(x1)−fj(x2))2≤Lj(x1−x2)2, | (2.10) |
where j=1,2,…,n, and Lj is a positive constant.
Remark 2. This study deals with several communication restrictions within T-S fuzzy complex dynamical networks. A decoupling approach is employed to explore the communication channel challenge in the analyzed system rather than relying on the Kronecker product. This shift significantly diminishes the computational load, especially when dealing with large-scale CDNs.
Definition 1. The T-S fuzzy complex dynamical networks (2.1) with ω(t)=0 is said to be asymptotic synchronization if
limt→∞ei(t)=limt→∞xi(t)−s(t)=0,i=1,2,…,N. | (2.11) |
Definition 2. [39] If matrices Θ1,Θ2,Θ3 and Θ4 are given, where Θ1,Θ3 and Θ4 are real symmetric, system (2.7) is considered to be extended dissipative when, for any nonnegative time tf and all ω(t)∈L2[0,∞), the following inequality holds:
∫tf0J(t)dt≥sup0≤t≤tf˜zT(t)Θ4˜z(t), | (2.12) |
where J(t)=˜zT(t)Θ1˜z(t)+2˜zT(t)Θ2ω(t)+ωT(t)Θ3ω(t).
In this paper, the general assumptions regarding Θ1,Θ2,Θ3 and Θ4 are applied consistently.
Assumption 2. [39] The real symmetric matrices Θ1,Θ2,Θ3 and Θ4 satisfy the following conditions:
(1) Θ1≤0,Θ3>0andΘ4≥0;
(2) (‖
Remark 3. The concept of extended dissipativity performance introduced in Definition 2 encompasses various performance measures as specific instances, achieved by adjusting the weighting matrices \Theta_1, \Theta_2, \Theta_3 and \Theta_4, in the following ways:
● When \Theta_1 = -I, \Theta_2 = 0, \Theta_3 = \delta^{2}I and \Theta_4 = 0, Definition 2 simplifies to the H_\infty performance.
● If \Theta_1 = 0, \Theta_2 = 0, \Theta_3 = \delta^{2}I and \Theta_4 = I, Definition 2 transforms into the L_2-L_\infty performance.
● In the case of \Theta_1 = 0, \Theta_2 = I, \Theta_3 = \delta I and \Theta_{4} = 0, Definition 2 represents passivity performance.
● When \Theta_{1} = Q, \Theta_{2} = S, \Theta_{3} = R-\delta I and \Theta_4 = 0, Definition 2 corresponds to (Q, S, R)-dissipativity performance.
Lemma 1. [34] For a matrix \mathcal{P} > 0 and a vector function \psi(t), t\in[a, b], the inequality is formulated as follows:
\begin{align} \int_{a}^{b}\dot{\psi}^{T}(t)\mathcal{P}\dot{\psi}(t)dt\geq\frac{1}{b-a}\Gamma^{T}_{1}\mathcal{P}\Gamma_{1}+\frac{3}{b-a}\Gamma^{T}_{2}\mathcal{P}\Gamma_{2}+\frac{5}{b-a}\Gamma^{T}_{3}\mathcal{P}\Gamma_{3}, \end{align} | (2.13) |
where
\begin{align*} \Gamma_{1} & = \psi(b)-\psi(a),\\ \Gamma_{2} & = \psi(b)+\psi(a)-\frac{2}{b-a}\int_{a}^{b}\psi(t)dt,\\ \Gamma_{3} & = \psi(b)-\psi(a)+\frac{6}{b-a}\int_{a}^{b}\psi(t)dt-\frac{12}{(b-a)^2}\int_{a}^{b}\int_{\beta}^{b}\psi(t)dtd\beta. \end{align*} |
Lemma 2. [43] For a matrix \mathcal{M} > 0, the following inequality holds:
\begin{align} -\frac{(b-a)^2}{2}\int_{a}^{b}\int_{u}^{b}\psi^T(s)\mathcal{M}\psi(s)dsdu\leq-\left(\int_{a}^{b}\int_{u}^{b}\psi(s)dsdu\right)^{T}\mathcal{M}\left(\int_{a}^{b}\int_{u}^{b}\psi(s)dsdu\right). \end{align} | (2.14) |
Lemma 3. [43] (Jensens Inequality) For a matrix \mathcal{Q} > 0 and a continuously differentiable function \psi(t) : [a, b] \rightarrow \; \mathbb{R}^{n} , the following statement is true:
\begin{align} (b-a)\int_{a}^{b}\psi^{T}(t)\mathcal{Q}\psi(t)dt\geq\int_{a}^{b}\psi^T(t)dt\mathcal{Q}\int_{a}^{b}\psi(t)dt. \end{align} | (2.15) |
Lemma 4. [44] For a matrix \mathcal{T} > 0, constant \gamma_{1}\leq\gamma(t)\leq\gamma_{2} and a function \dot{x}(t): [-\gamma_{2}, \gamma_{1}]\rightarrow \mathbb{R}^{n} such that the following integration is well defined, then the following inequality holds:
\begin{align} -(\gamma_{2}-\gamma_{1})\int_{t-\gamma_{2}}^{t-\gamma_{1}}\dot{x}^{T}(s)\mathcal{T}\dot{x}(s)ds\leq\begin{bmatrix} x(t-\gamma_{1})\\ x(t-\gamma(t))\\ x(t-\gamma_{2})\end{bmatrix}^{T}\begin{bmatrix} -\mathcal{T} & \mathcal{T} & 0 \\ 0 & -2\mathcal{T} & \mathcal{T} \\ 0 & 0 & -\mathcal{T} \end{bmatrix}\begin{bmatrix} x(t-\gamma_{1})\\ x(t-\gamma(t))\\ x(t-\gamma_{2})\end{bmatrix}. \end{align} | (2.16) |
Lemma 5. [29] Given matrices \mathcal{H} and \mathcal{N} , and a matrix \mathcal{F}(t) that satisfies \mathcal{F}(t)\mathcal{F}^{T}(t)\leq I , for any positive value of \epsilon , the following statement is true:
\begin{align} \mathcal{H}\mathcal{F}(t)\mathcal{N} + \mathcal{N}^{T}\mathcal{F}^{T}(t)\mathcal{H}^T\leq\frac{1}{\epsilon}\mathcal{HH}^{T} + \epsilon \mathcal{NN}^{T}. \end{align} | (2.17) |
In this section, we present certain notations to aid in expressing matrices and vectors more conveniently. The notations are described below:
\begin{align} \zeta_{1}(t) = &\left[ e^T_{l}(t)\; \dot{e}^T_{l}(t)\; e^T_{l}(t-\gamma_{1})\;e^T_{l}(t-\gamma_{2})\;\frac{1}{\gamma_{1}}\int_{t-\gamma_{1}}^{t}e^{T}_{l}(s)ds\;\frac{1}{\gamma^2_{1}}\int_{t-\gamma_{1}}^{t}\int_{u}^{t}e^{T}_{l}(s)dsdu\;\frac{1}{\gamma_{2}}\int_{t-\gamma_{2}}^{t}e^{T}_{l}(s)ds\right.\\ &\left. \frac{1}{\gamma^2_{2}}\int_{t-\gamma_{2}}^{t}\int_{u}^{t}e^{T}_{l}(s)dsdu\;\int_{t-\varsigma_{2}(t)}^{t-\varsigma_{1}(t)}e^{T}_{l}(s)ds\;e^{T}_{l}(t-\gamma(t))\right]^{T},\\ \zeta_{2}(t) = & \left[\zeta^{T}_{1}(t)\;\omega^{T}_{l}(t)\right]^{T},\\ e_{i} = & \left[0_{N\times(i-1)N},\;\;I_{N},\;\;0_{N\times(10-i)N}\right]\quad \text{for}\; i = 1,2,\dots ,10. \end{align} | (3.1) |
We will introduce a novel sufficient synchronization criterion for the analyzed complex dynamical networks (2.1) when \omega_{l}(t) = 0.
Theorem 1. Given positive scalar \epsilon and nonnegative scalars \gamma_1, \gamma_2, c_{i}, (i = 1, 2, 3), if there exists positive matrices P_{il}\in\mathbb{R}^{N\times N}, (i = 1, 2, 3), \; Q_{il}\in\mathbb{R}^{N\times N}, (i = 1, 2, 3, 4), \; R_{l}\in\mathbb{R}^{N\times N} and any matrices \Lambda_{1l}, \Lambda_{2l}\in\mathbb{R}^{N\times N}, F_{l} = diag\bigg\{F_{1}, F_{2}\dots, F_{N}\bigg\}, \Gamma = diag\bigg\{\Gamma_{1}, \Gamma_{2}\dots, \Gamma_{N}\bigg\}\; ,l = 1, 2, \dots, n, such that the LMIs hold as follows:
\begin{align} \begin{bmatrix} \hat{\Pi}^{(1)} & \sqrt{n}e^T_{1}F^T_{l} & \sqrt{n}e^T_{2}F^T_{l} & N_{1} & N_{2}\\ * & -\Lambda_{1l} & 0 & 0 & 0\\ * & * & -\Lambda_{2l} & 0 & 0\\ * & * & * & -\epsilon I & 0\\ * & * & * & * & -\epsilon I \end{bmatrix} < 0, \end{align} | (3.2) |
where
\begin{align} \hat{\Pi}^{(1)} = & sym\big\{e^T_{1}P_{1l}e_{2}\big\}+e^T_{1}P_{2l}e_{1}-e^T_{3}P_{2l}e_{3}+e^T_{1}P_{3l}e_{1}-e^T_{4}P_{3l}e_{4}+\gamma^2_{1}e^T_{2}Q_{1l}e_{2}+\gamma^2_{2}e^T_{2}Q_{2l}e_{2}\\ &-(e_{1}-e_{3})^{T}Q_{1l}(e_{1}-e_{3})-3(e_{1}+e_{3}-2e_{5})^{T}Q_{1l}(e_{1}+e_{3}-2e_{5})\\ &-5(e_{1}-e_{3}+6e_{5}-12e_{6})^{T}Q_{1l}(e_{1}-e_{3}+6e_{5}-12e_{6})\\ &-(e_{1}-e_{4})^{T}Q_{2l}(e_{1}-e_{4})-3(e_{1}+e_{4}-2e_{7})^{T}Q_{2l}(e_{1}+e_{4}-2e_{7})\\ &-5(e_{1}-e_{4}+6e_{7}-12e_{8})^{T}Q_{2l}(e_{1}-e_{4}+6e_{7}-12e_{8})\\ &+\frac{\gamma^4_{1}}{4}e^T_{2}Q_{3l}e_{2}-(e_{1}-e_{5})^{T}\gamma^2_{1}Q_{3l}(e_{1}-e_{5})+(\varsigma_{2}-\varsigma_{1})^{2}e^T_{1}Q_{4l}e_{1}-e^T_{9}Q_{4l}e_{9}\\ &+(\gamma_{2}-\gamma_{1})^{2}e^T_{2}R_{l}e_{2}+\begin{bmatrix} e_{3}\\ e_{10}\\ e_{4} \end{bmatrix}^{T}\begin{bmatrix} -R_{l} & R_{l} & 0\\ * & -2R_{l} & R_{l}\\ * & * & -R_{l} \end{bmatrix}\begin{bmatrix} e_{3}\\ e_{10}\\ e_{4} \end{bmatrix}\\ &+sym\bigg\{(e_{1}+e_{2})^{T}F^T_{l}\left(a_{dl}e_{1}+c_{1}\mathcal{G}^{l}_{d}e_{1}+c_{2}\mathcal{G}^{l}_{d}e_{10}+c_{3}\mathcal{G}^{l}_{d}e_{9}+K^{d}e_{1}-e_{2}\right)\bigg\}\\ &+e^T_{1} \sum\limits_{j = 1}^{n}b^2_{dlj}L_{j}(\Lambda_{1l}+\Lambda_{2l})e_{1}\\ N_{1}& = e^T_{1}E^T_{l}\\ N_{2}& = \epsilon(e_{1}+e_{2})^{T}F^T_{l}C_{l}, \end{align} | (3.3) |
then, the error system (2.7) achieves asymptotic synchronization. Additionally, the control gain can be determined as K_{l}^{d} = (F_{l}^{T})^{-1}\Gamma .
Proof. The candidate Lyapunov functional for system (2.7) is formulated as follows:
\begin{align} V(t) = \sum\limits_{k = 1}^{6}V_{k}(t), \end{align} | (3.4) |
where
\begin{align} V_{1}(t) & = \sum\limits_{l = 1}^{n}e^{T}_{l}(t)P_{1l}e_{l}(t),\\ V_{2}(t) & = \sum\limits_{l = 1}^{n}\left[\int_{t-\gamma_{1}}^{t}e^T_{l}(s)P_{2l}e_{l}(s)ds+\int_{t-\gamma_{2}}^{t}e^T_{l}(s)P_{3l}e_{l}(s)ds\right],\\ V_{3}(t)& = \sum\limits_{l = 1}^{n}\left[\gamma_{1}\int_{-\gamma_{1}}^{0}\int_{t+u}^{t}\dot{e}^T_{l}(s)Q_{1l}\dot{e}_{l}(s)dsdu+\gamma_{2}\int_{-\gamma_{2}}^{0}\int_{t+u}^{t}\dot{e}^T_{l}(s)Q_{2l}\dot{e}_{l}(s)dsdu\right],\\ V_{4}(t)& = \sum\limits_{l = 1}^{n}\left[\frac{\gamma^2_{1}}{2}\int_{t-\gamma_{1}}^{t}\int_{u}^{t}\int_{\theta}^{t}\dot{e}^T_{l}(s)Q_{3l}\dot{e}_{l}(s)dsd\theta du\right],\\ V_{5}(t)& = \sum\limits_{l = 1}^{n}\left[(\varsigma_{2}-\varsigma_{1})\int_{-\varsigma_{2}}^{-\varsigma_{1}}\int_{t+u}^{t}e^T_{l}(s)Q_{4l}e_{l}(s)dsdu\right],\\ V_{6}(t)& = \sum\limits_{l = 1}^{n}\left[(\gamma_{2}-\gamma_{1})\int_{-\gamma_{2}}^{-\gamma_{1}}\int_{t+u}^{t}\dot{e}^T_{l}(s)R_{l}\dot{e}_{l}(s)dsdu\right]. \end{align} | (3.5) |
By computing the time derivative of V(t), we have
\begin{align} \dot{V}_{1}(t)& = \sum\limits_{l = 1}^{n}2e^T_{l}(t)P_{1l}\dot{e}_{l}(t) = \sum\limits_{l = 1}^{n} 2\zeta_{1}^T(t)e^T_{1}P_{1l}e_{2}\zeta_{1}(t), \end{align} | (3.6) |
\begin{align} \dot{V}_{2}(t)& = \sum\limits_{l = 1}^{n}\left[e^T_{l}(t)P_{2l}e_{l}(t)-e^T_{l}(t-\gamma_{1})P_{2l}e_{l}(t-\gamma_{1})+e^T_{l}(t)P_{3l}e_{l}(t)-e^T_{l}(t-\gamma_{2})P_{3l}e_{l}(t-\gamma_{2})\right]\\ & = \sum\limits_{l = 1}^{n}\left[\zeta^T_{1}(t)\left(e^T_{1}P_{2l}e_{1}-e^T_{3}P_{2l}e_{3}+e^T_{1}P_{3l}e_{1}-e^T_{4}P_{3l}e_{4}\right)\zeta_{1}(t)\right], \end{align} | (3.7) |
\begin{align} \dot{V}_{3}(t) = & \sum\limits_{l = 1}^{n}\left[\gamma^2_{1}\dot{e}^T_{l}(t)Q_{1l}\dot{e}_{l}(t)-\gamma_{1}\int_{t-\gamma_{1}}^{t}\dot{e}^T_{l}(s)Q_{1l}\dot{e}_{l}(s)ds+\gamma^2_{2}\dot{e}^T_{l}(t)Q_{2l}\dot{e}_{l}(t)\right.\\ &\left.-\gamma_{2}\int_{t-\gamma_{2}}^{t}\dot{e}^T_{l}(s)Q_{2l}\dot{e}_{l}(s)ds\right]. \end{align} | (3.8) |
By using Lemma 1, we gained:
\begin{align} -\gamma_{1}\int_{t-\gamma_{1}}^{t}\dot{e}^T_{l}(s)Q_{1l}\dot{e}_{l}(s)ds\leq& -(e_{1}-e_{3})^{T}Q_{1l}(e_{1}-e_{3})-3(e_{1}+e_{3}-2e_{5})^{T}Q_{1l}(e_{1}+e_{3}-2e_{5})\\ &-5(e_{1}-e_{3}+6e_{5}-12e_{6})^{T}Q_{1l}(e_{1}-e_{3}+6e_{5}-12e_{6}) \end{align} | (3.9) |
\begin{align} -\gamma_{2}\int_{t-\gamma_{2}}^{t}\dot{e}^T_{l}(s)Q_{2l}\dot{e}_{l}(s)ds\leq& -(e_{1}-e_{4})^{T}Q_{2l}(e_{1}-e_{4})-3(e_{1}+e_{4}-2e_{7})^{T}Q_{2l}(e_{1}+e_{4}-2e_{7})\\ &-5(e_{1}-e_{4}+6e_{7}-12e_{8})^{T}Q_{2l}(e_{1}-e_{4}+6e_{7}-12e_{8}). \end{align} | (3.10) |
\dot{V}_{4} can be calculated as follows:
\begin{align} \dot{V}_{4}(t)& = \sum\limits_{l = 1}^{n}\left[\frac{\gamma^4_{1}}{4}\dot{e}^T_{l}(t)Q_{3l}\dot{e}_{l}(t)-\frac{\gamma^2_{1}}{2}\int_{t-\gamma_{1}}^{t}\int_{u}^{t}\dot{e}^T_{l}(s)Q_{3l}\dot{e}_{l}(s)dsdu\right]. \end{align} | (3.11) |
By utilizing Lemma 2, we have
\begin{align} -\frac{\gamma^2_{1}}{2}\int_{t-\gamma_{1}}^{t}\int_{u}^{t}\dot{e}^T_{l}(s)Q_{3l}\dot{e}_{l}(s)dsdu\leq-\left(e_{1}-e_{5}\right)^{T}\gamma^2_{1}Q_{3l}\left(e_{1}-e_{5}\right). \end{align} | (3.12) |
\dot{V}_{5} can be calculated as follows:
\begin{align} \dot{V}_{5}(t)& = \sum\limits_{l = 1}^{n}\left[(\varsigma_{2}-\varsigma_{1})^{2}e^T_{l}(t)Q_{4l}e_{l}(t)-(\varsigma_{2}-\varsigma_{1})\int_{t-\varsigma_{2}}^{t-\varsigma_{1}}e^T_{l}(s)Q_{4l}e_{l}(s)ds\right]. \end{align} | (3.13) |
After utilizing Lemma 3, we gain
\begin{align} -(\varsigma_{2}-\varsigma_{1})\int_{t-\varsigma_{2}}^{t-\varsigma_{1}}e^T_{l}(s)Q_{4l}e_{l}(s)ds&\leq-(\varsigma_{2}(t)-\varsigma_{1}(t))\int_{t-\varsigma_{2}(t)}^{t-\varsigma_{1}(t)}e^T_{l}(s)Q_{4l}e_{l}(s)ds\\ &\leq-e^T_{9}Q_{4l}e_{9}. \end{align} | (3.14) |
\dot{V}_{6} can be calculated as follows:
\begin{align} \dot{V}_{6}(t) & = \sum\limits_{l = 1}^{n}\left[(\gamma_{2}-\gamma_{1})^{2}\dot{e}^T_{l}(t)R_{l}\dot{e}_{l}(s)ds-(\gamma_{2}-\gamma_{1})\int_{t-\gamma_{2}}^{t-\gamma_{1}}\dot{e}^T_{l}(s)R_{l}\dot{e}_{l}(s)ds\right]. \end{align} | (3.15) |
By applying Lemma 4, we have
\begin{align} -(\gamma_{2}-\gamma_{1})\int_{t-\gamma_{2}}^{t-\gamma_{1}}\dot{e}^T_{l}(s)R_{l}\dot{e}_{l}(s)ds\leq\begin{bmatrix} e_{3}\\ e_{10}\\ e_{4} \end{bmatrix}^{T}\begin{bmatrix} -R_{l} & R_{l} & 0\\ * & -2R_{l} & R_{l}\\ * & * & -R_{l} \end{bmatrix}\begin{bmatrix} e_{3}\\ e_{10}\\ e_{4} \end{bmatrix}. \end{align} | (3.16) |
Based on the error system (2.7), we can deduce that for any suitable matrix F_{l}, the following equation is true:
\begin{align} 0 = & 2 \sum\limits_{l = 1}^{n}(e^T_{l}(t)+\dot{e}^T_{l}(t))F^T_{l} \sum\limits_{d = 1}^{r}\rho_{d}(\chi(t))\bigg\{a_{dl}e_{l}(t)+ \sum\limits_{j = 1}^{n}b_{dlj}f_{j}(e_{lj}(t))+c_{1}\mathcal{G}^{l}_{d}e_{l}(t)+c_{2}\mathcal{G}^{l}_{d}e_{l}(t-\gamma(t))\\ &+ c_{3}\mathcal{G}^{l}_{d}\int_{t-\varsigma_{2}(t)}^{t-\varsigma_{1}(t)}e_{l}(s)ds+(K^{d}_{l}+C_{l}D_{l}(t)E_{l})e_{l}(t)-\dot{e}_{l}(t)\bigg\}, \end{align} | (3.17) |
\begin{align} &2(e^T_{l}(t)+\dot{e}^T_{l}(t))F^T_{l} \sum\limits_{j = 1}^{n}b_{dlj}f_{j}(e_{lj}(t))\leq ne^T_{l}(t)F^T_{l}\Lambda^{-1}_{1l}F_{l}e_{l}(t) + n\dot{e}^T_{l}(t)F^T_{l}\Lambda^{-1}_{2l}F_{l}\dot{e}_{l}(t)\\ &\hskip6cm + \sum\limits_{j = 1}^{n}b^2_{dlj}L_{j}e^T_{j}(t)(\Lambda_{1l}+\Lambda_{2l})e_{j}(t). \end{align} | (3.18) |
Then, by combining (3.6) – (3.18), we get
\begin{align} \dot{V}(t)\leq \sum\limits_{d = 1}^{r}\rho_{d}(\chi(t)) \sum\limits_{l = 1}^{n}\zeta_{1}^T(t)\Omega\zeta_{1}(t), \end{align} | (3.19) |
where
\begin{align} \Omega = & \hat{\Pi}^{(1)} + ne^T_{1}F^T_{l}\Lambda^{-1}_{1l}F_{l}e_{1}+ne^T_{2}F^T_{l}\Lambda^{-1}_{2l}F_{l}e_{2} + (e_{1}+e_{2})^T\left(F^T_{l}C_{l}D_{l}(t)E_{l}\right)e_{1} +e^T_{1}\left(E^T_{l}D^T_{l}(t)C^T_{l}F_{l}\right)(e_{1}+e_{2}). \end{align} | (3.20) |
By using Lemma 5,
\begin{align} \Omega \leq \hat{\Pi}^{(1)} + \epsilon(e_{1}+e_{2})^T\left(F^T_{l}C_{l}\right)\left(C^T_{l}F_{l}\right)(e_{1}+e_{2})+\epsilon^{-1}\left(e^T_{1}E^T_{l}\right)\left(E_{l}e_{1}\right). \end{align} | (3.21) |
By utilizing the Schur complement, we can derive \Omega < 0 from (3.2), which in turn ensures \dot{V}(t) < 0. According to the Lyapunov stability theory, the T-S fuzzy complex dynamical network (2.1) can achieve asymptotic synchronization. This concludes the proof.
We conduct an extended dissipative performance analysis of the considered CDNs (2.1) under various initial conditions, where \omega_{l}(t)\not = 0 and \omega_{l}(t)\in\mathcal{L}_{2}[0, \infty) .
Theorem 2. Given positive scalar \epsilon and nonnegative scalars \gamma_1, \gamma_2, c_{i}, i = 1, 2, 3, \; \Theta_{i}\; (i = 1, 2, \dots, 4) satisfying Assumption 2, if there exists positive matrices P_{il}\in\mathbb{R}^{N\times N}, (i = 1, 2, 3), \; Q_{il}\in\mathbb{R}^{N\times N}, (i = 1, 2, 3, 4), \; R_{l}\in\mathbb{R}^{N\times N} and any matrices \Lambda_{1l}, \Lambda_{2l}\in\mathbb{R}^{N\times N}, F_{l} = diag\bigg\{F_{1}, F_{2}\dots, F_{N}\bigg\}, \Gamma = diag\bigg\{\Gamma_{1}, \Gamma_{2}\dots, \Gamma_{N}\bigg\}\; ,l = 1, 2, \dots, n, such that the LMIs hold as follows:
\begin{align} &\begin{bmatrix} \hat{\Pi}^{(2)} & \sqrt{n}\tilde{e}^T_{1}F^T_{l} & \sqrt{n}\tilde{e}^T_{2}F^T_{l} & N_{1} & N_{2}\\ * & -\Lambda_{1l} & 0 & 0 & 0\\ * & * & -\Lambda_{2l} & 0 & 0\\ * & * & * & -\epsilon I & 0\\ * & * & * & * & -\epsilon I \end{bmatrix} < 0, \end{align} | (3.22) |
\begin{align} &P_{1l} - \hat{\Im}^{T}\Theta_{4}\hat{\Im} > 0, \end{align} | (3.23) |
where
\begin{align} \hat{\Pi}^{(2)} = & sym\big\{\tilde{e}^T_{1}P_{1l}\tilde{e}_{2}\big\}+\tilde{e}^T_{1}P_{2l}\tilde{e}_{1}-\tilde{e}^T_{3}P_{2l}\tilde{e}_{3}+\tilde{e}^T_{1}P_{3l}\tilde{e}_{1}-\tilde{e}^T_{4}P_{3l}\tilde{e}_{4}+\gamma^2_{1}\tilde{e}^T_{2}Q_{1l}\tilde{e}_{2}+\gamma^2_{2}\tilde{e}^T_{2}Q_{2l}\tilde{e}_{2}\\ &-(\tilde{e}_{1}-\tilde{e}_{3})^{T}Q_{1l}(\tilde{e}_{1}-\tilde{e}_{3})-3(\tilde{e}_{1}+\tilde{e}_{3}-2\tilde{e}_{5})^{T}Q_{1l}(\tilde{e}_{1}+\tilde{e}_{3}-2\tilde{e}_{5})\\ &-5(\tilde{e}_{1}-\tilde{e}_{3}+6\tilde{e}_{5}-12\tilde{e}_{6})^{T}Q_{1l}(\tilde{e}_{1}-\tilde{e}_{3}+6\tilde{e}_{5}-12\tilde{e}_{6})\\ &-(\tilde{e}_{1}-\tilde{e}_{4})^{T}Q_{2l}(\tilde{e}_{1}-\tilde{e}_{4})-3(\tilde{e}_{1}+\tilde{e}_{4}-2\tilde{e}_{7})^{T}Q_{2l}(\tilde{e}_{1}+\tilde{e}_{4}-2\tilde{e}_{7})\\ &-5(\tilde{e}_{1}-\tilde{e}_{4}+6\tilde{e}_{7}-12\tilde{e}_{8})^{T}Q_{2l}(\tilde{e}_{1}-\tilde{e}_{4}+6\tilde{e}_{7}-12\tilde{e}_{8})\\ &+\frac{\gamma^4_{1}}{4}\tilde{e}^T_{2}Q_{3l}\tilde{e}_{2}-(\tilde{e}_{1}-\tilde{e}_{5})^{T}\gamma^2_{1}Q_{3l}(\tilde{e}_{1}-\tilde{e}_{5})+(\varsigma_{2}-\varsigma_{1})^{2}\tilde{e}^T_{1}Q_{4l}\tilde{e}_{1}-\tilde{e}^T_{9}Q_{4l}\tilde{e}_{9}\\ &+(\gamma_{2}-\gamma_{1})^{2}\tilde{e}^T_{2}R_{l}\tilde{e}_{2}+\begin{bmatrix} \tilde{e}_{3}\\ \tilde{e}_{10}\\ \tilde{e}_{4} \end{bmatrix}^{T}\begin{bmatrix} -R_{l} & R_{l} & 0\\ * & -2R_{l} & R_{l}\\ * & * & -R_{l} \end{bmatrix}\begin{bmatrix} \tilde{e}_{3}\\ \tilde{e}_{10}\\ \tilde{e}_{4} \end{bmatrix}\\ &+sym\bigg\{(\tilde{e}_{1}+\tilde{e}_{2})^{T}F^T_{l}\left(a_{dl}\tilde{e}_{1}+c_{1}\mathcal{G}^{l}_{d}\tilde{e}_{1}+c_{2}\mathcal{G}^{l}_{d}\tilde{e}_{10}+c_{3}\mathcal{G}^{l}_{d}\tilde{e}_{9}+K^{d}\tilde{e}_{1}-\tilde{e}_{2}\right)\bigg\}\\ &+\tilde{e}^T_{1} \sum\limits_{j = 1}^{n}b^2_{dlj}L_{j}(\Lambda_{1l}+\Lambda_{2l})\tilde{e}_{1}-\tilde{e}^T_{1}\Im^T\Theta_{1}\Im \tilde{e}_{1}-sym\bigg\{\tilde{e}^T_{1}\Im^{T}\Theta_{2}\tilde{e}_{11}\bigg\}-\tilde{e}^T_{11}\Theta_{3}\tilde{e}_{11}\\ N_{1}& = \tilde{e}^T_{1}E^T_{l}\\ N_{2}& = \epsilon(\tilde{e}_{1}+\tilde{e}_{2})^{T}F^T_{l}C_{l},\\ \tilde{e}_{i} & = \left[0_{N\times(i-1)N},\;\;I_{N},\;\;0_{N\times(11-i)N}\right]\quad \mathit{\text{for}}\; i = 1,2,\dots ,11. \end{align} | (3.24) |
Then, the error system (2.7) is extended dissipative, with all other parameters unchanged as described in Theorem 1.
Proof. By applying the Schur complement and Lemma 5, it is clear that (3.22) is equivalent to
\begin{align} \sum\limits_{d = 1}^{r}\rho_{d}(\chi(t)) \sum\limits_{l = 1}^{n}\zeta^T_{2}(t)\tilde{\Omega}\zeta_{2}(t) < 0, \end{align} | (3.25) |
where
\tilde{\Omega} = \hat{\Pi}^{(2)} + n\tilde{e}^T_{1}F^T_{l}\Lambda^{-1}_{1l}F_{l}\tilde{e}_{1}+n\tilde{e}^T_{2}F^T_{l}\Lambda^{-1}_{2l}F_{l}\tilde{e}_{2} + (\tilde{e}_{1}+\tilde{e}_{2})^T\left(F^T_{l}C_{l}D_{l}(t)E_{l}\right)\tilde{e}_{1} +\tilde{e}^T_{1}\left(E^T_{l}D^T_{l}(t)C^T_{l}F_{l}\right)(\tilde{e}_{1}+\tilde{e}_{2}).
From (3.19), we have
\begin{align} \dot{V}(t)\leq \sum\limits_{d = 1}^{r}\rho_{d}(\chi(t)) \sum\limits_{l = 1}^{n}\zeta_{1}^T(t)\Omega\zeta_{1}(t), \end{align} | (3.26) |
and it is clear that
\begin{align} \sum\limits_{d = 1}^{r}\rho_{d}(\chi(t)) \sum\limits_{l = 1}^{n}\zeta_{2}^T(t)\tilde{\Omega}\zeta_{2}(t)& = \sum\limits_{d = 1}^{r}\rho_{d}(\chi(t)) \sum\limits_{l = 1}^{n}\zeta_{1}^T(t)\Omega\zeta_{1}(t)-\mathcal{J}(t), \end{align} | (3.27) |
where \mathcal{J}(t) is provided in Definition 2.
Then, we have
\begin{align} \dot{V}(t) \leq \sum\limits_{d = 1}^{r}\rho_{d}(\chi(t)) \sum\limits_{l = 1}^{n}\zeta_{2}^T(t)\tilde{\Omega}\zeta_{2}(t) + \mathcal{J}(t) \leq \mathcal{J}(t). \end{align} | (3.28) |
By integrating both sides of inequality (3.28) from 0 to t \; (t\geq 0) , we obtain the following result:
\begin{align} \int_{0}^{t}\mathcal{J}(s)ds\geq V(t)-V(0)\geq e^T_{l}(t)P_{1l}e_{l}(t). \end{align} | (3.29) |
Our attention is specifically directed toward two cases: \Theta_{4} = 0 and \Theta_{4}\geq 0 . The rationale behind this focus is that the extended dissipative performance can encompass the conditions of strictly (Q, S, R) -dissipativity, H_{\infty} and passive performance when \Theta_{4} = 0 or the L_{2} - L_{\infty} performance when \Theta_{4} > 0 .
Taking \Theta_{4} = 0 into account, the inequality (3.29) implies the following:
\begin{align} \int_{0}^{t_{f}}\mathcal{J}(s)ds\geq 0. \end{align} | (3.30) |
Simultaneously, in the case where \Theta_{4} > 0 as indicated in Assumption 2, this implies that the matrices are \Theta_{1} = 0 , \Theta_{2} = 0 and \Theta_{3} > 0 . Consequently, for t\in[0, t_{f}] , the inequality (3.29) leads to the following chain of relations: \int_{0}^{t_{f}}\mathcal{J}(s)ds\geq\int_{0}^{t}\mathcal{J}(s)ds\geq e^T_{l}(t)P_{1l}e_{l}(t) . Subsequently, utilizing (3.23), we can deduce the following expression:
\begin{align} \hat{z}^T(t)\Theta_{4}\hat{z}(t)& = e^T_{l}(t)\hat{\Im}^{T}\Theta_{4}\hat{\Im}e_{l}(t)\\ &\leq e^T_{l}(t)P_{1l}e_{l}(t)\leq\int_{0}^{t_{f}}\mathcal{J}(s)ds. \end{align} | (3.31) |
By combining (3.30) and (3.31), system (2.1) is extended dissipative synchronization. This concludes the proof.
Remark 4. To make a comparison with the controller designed in [34], we utilized a single free-matrix parameter F in Eq (3.17) instead of the two parameters used in the zero equation in [34]. Importantly, we were able to achieve this without introducing any additional variables.
Remark 5. Employing multiple integral components within the Lyapunov function addresses the challenges of mixed interval time-varying delays. This paper introduces a less conservative stability criterion using auxiliary function-dependent integral inequalities and mathematical approaches. In real-world scenarios, if the system possesses initial energy, this energy will progressively diminish as time elapses, ultimately guiding the system toward a state of equilibrium. However, determining the energy function of the system can be inconvenient. Consequently, Lyapunov introduced the function V(x, t), serving as an energy representation known as the Lyapunov equation.
Remark 6. The external disturbance \omega(t) can generally impact the system's stability. Various robust performance measures, including H_\infty, L_{2}-L_{\infty}, passivity, and dissipativity, are utilized to assess system stability. These measures offer improved ways of gauging stability. By adjusting weighting matrices and scalars, extended dissipative performance can be transformed into four standard individual performance criteria, enhancing the comprehensiveness and generality of system analysis. A recent study [34] established H_\infty synchronization conditions for uncertain CDNs with hybrid delays using inequalities such as Jensens inequality and Wirtinger inequality. However, there has been no prior investigation into the extended dissipativity criteria through a synchronization approach in T-S fuzzy CDNs.
This section will present two simulation examples to showcase reduced conservatism and validate the proposed method.
Example 1. Let us examine the T-S fuzzy complex dynamical networks with four nodes. Each node has three information vectors, and the membership functions are \rho_{1}(\chi(t)) = \sin^{2}\frac{t}{4} and \rho_{2}(\chi(t)) = \cos^{2}\frac{t}{4}. The following network parameters are provided:
\begin{align*} A_{1} & = \begin{bmatrix} -0.1 & 0 & 0\\ 0 & -0.1 & 0\\ 0 & 0 & -0.1 \end{bmatrix},\;A_{2} = \begin{bmatrix} -0.1 & 0 & 0\\ 0 & -0.1 & 0\\ 0 & 0 & -0.1 \end{bmatrix},\; B_{1} = \begin{bmatrix} 1.1 & 0 & -1.2\\ 1.2 & 1.4 & 0.7\\ 1.1 & 0 & 0.5 \end{bmatrix}\\ B_{2} & = \begin{bmatrix} -1.3 & 0 & 0.9\\ 0.8 & 0.5 & 1.2\\ 1 & 1 & 0 \end{bmatrix},\; \mathcal{F} = \begin{bmatrix} -0.8 & 0.4 & 0.2 & 0.2\\ 0.1 & -2.4 & 0.3 & 2\\ 3 & 0.1 & -4.1 & 1\\ 0.5 & 0.5 & 2.2 & -3.2 \end{bmatrix}. \end{align*} |
The diagonal matrix \Upsilon_{1ij} represents the channel matrix associated with Rule 1 and is defined as follows:
\begin{align*} \Upsilon_{112} & = diag\{1,1,0\},\;\Upsilon_{113} = diag\{1,0,0\},\; \Upsilon_{114} = diag\{0,0,0\}, \\ \Upsilon_{121} & = diag\{0,1,1\},\; \Upsilon_{123} = diag\{1,0,0\},\;\Upsilon_{124} = diag\{0,0,0\},\\ \Upsilon_{131}& = diag\{0,0,0\},\;\Upsilon_{132} = diag\{0,0,1\},\;\Upsilon_{134} = diag\{1,0,1\},\\ \Upsilon_{141} & = diag\{1,0,0\},\; \Upsilon_{142} = diag\{0,0,0\},\;\Upsilon_{143} = diag\{0,0,1\}. \end{align*} |
Additionally, in the context of Rule 2, the diagonal matrix \Upsilon_{2ij} is defined as follows:
\begin{align*} \Upsilon_{212} & = diag\{1,0,0\},\;\Upsilon_{213} = diag\{0,0,0\},\; \Upsilon_{214} = diag\{0,1,0\},\\ \Upsilon_{221} & = diag\{0,0,1\},\; \Upsilon_{223} = diag\{0,0,0\},\;\Upsilon_{224} = diag\{0,0,1\},\\ \Upsilon_{231}& = diag\{1,1,0\},\;\Upsilon_{232} = diag\{0,0,0\},\;\Upsilon_{234} = diag\{0,0,1\},\\ \Upsilon_{241} & = diag\{0,0,0\},\; \Upsilon_{242} = diag\{1,0,0\},\;\Upsilon_{243} = diag\{0,0,1\}. \end{align*} |
Figure 1 illustrates the information pathways of the network for Rule 1. It uses two-way arrows (i \leftrightarrow j) to indicate the activation of information channels between node i and j and one-way arrows (i \rightarrow j) to represent the activation of information channels from node i to j.
Next, a decoupling method addresses the partial coupling between linked nodes. The N \times N matrix \mathcal{G}^{l}_{1}, (l = 1, 2, 3) represents the l^{th} channel within a network of N nodes, following the guidelines of Rule 1. It can be readily observed that:
\begin{align*} \mathcal{G}^{1}_{1} & = \begin{bmatrix} -0.6 & 0.4 & 0.2 & 0\\ 0 & -0.3 & 0.3 & 0\\ 0 & 0 & -1 & 1\\ 0.5 & 0 & 0 & -0.5 \end{bmatrix},\;\mathcal{G}^{2}_{1} = \begin{bmatrix} -0.4 & 0.4 & 0 & 0\\ 0.1 & -0.1 & 0 & 0\\ 0 & 0 & 0 & 0\\ 0 & 0 & 0 & 0 \end{bmatrix},\;\mathcal{G}^{3}_{1} = \begin{bmatrix} 0 & 0 & 0 & 0\\ 0.1 & -0.1 & 0 & 0\\ 0 & 0.1 & -1.1 & 1\\ 0 & 0 & 2.2 & -2.2 \end{bmatrix}. \end{align*} |
Figure 2 depicts the network's channel configuration according to Rule 2. Using a similar analysis as mentioned earlier, the computation of matrix \mathcal{G}^{l}_{2} (l = 1, 2, 3) is determined as follows:
\begin{align*} \mathcal{G}^{1}_{2} & = \begin{bmatrix} -0.4 & 0.4 & 0 & 0\\ 0 & 0 & 0 & 0\\ 3 & 0 & -3 & 0\\ 0 & 0.5 & 0 & -0.5 \end{bmatrix},\;\mathcal{G}^{2}_{2} = \begin{bmatrix} -0.2 & 0 & 0 & 0.2\\ 0 & 0 & 0 & 0\\ 3 & 0 & -3 & 0\\ 0 & 0 & 0 & 0 \end{bmatrix},\;\mathcal{G}^{3}_{2} = \begin{bmatrix} 0 & 0 & 0 & 0\\ 0.1 & -2.1 & 0 & 2\\ 0 & 0 & -1 & 1\\ 0 & 0 & 2.2 & -2.2 \end{bmatrix}. \end{align*} |
The nonlinear function f(x_{i}(t)) is chosen as:
\begin{align*} f(x_{i}(t)) = \begin{bmatrix} f_{1}(x_{i1}(t))\\ f_{2}(x_{i2}(t))\\ f_{3}(x_{i3}(t)) \end{bmatrix} = \begin{bmatrix} \frac{\vert x_{i1}(t) + 1\vert-\vert x_{i1}(t) - 1\vert}{2}\\ \frac{\vert x_{i2}(t) + 1\vert-\vert x_{i2}(t) - 1\vert}{2}\\ \frac{\vert x_{i3}(t) + 1\vert-\vert x_{i3}(t) - 1\vert}{2} \end{bmatrix}. \end{align*} |
It is possible to determine that any element f_l(x_{il}) (l = 1, 2, 3; i = 1, 2, 3, 4) meets the Lipschitz condition with L_l = 1.
In this case, we take \omega(t) = 0, and the uncertain parameter matrices are taken as: C_{l} = diag\{0.1, 0.1, 0.1, 0.1\}, \; E_{l} = diag\{0.2, 0.2, 0.2, 0.2\} and D_{l}(t) = diag\{\sin t, \sin t, \cos t, \sin t\}. Moreover, c_{1} = c_{2} = c_{3} = 1, \gamma(t) = 0.2 \vert\tanh{t}\vert, \; \varsigma_{1}(t) = 0.1+0.1\sin{t} and \varsigma_{2}(t) = 0.2\sin{t}, which means \gamma_{1} = 0.1 and \gamma_{2} = 0.2, \sigma_{1} = 0.1, \sigma_{2} = 0.2 .Then, the non-fragile control matrices can be calculated as:
\begin{align*} K^{1}_{1} & = diag\{-16.0001,-14.0399,-56.4629,-49.6643\},\\ K^{1}_{2} & = diag\{-20.3844,-18.0041,-68.7015,-61.4408\},\\ K^{1}_{3} & = diag\{-21.5908,-19.0994,-72.0947,-64.6059\},\\ K^{2}_{1} & = diag\{-27.5577,-24.5250,-88.5943,-79.5904\},\\ K^{2}_{2} & = diag\{-9.5513,-8.5177,-36.4208,-31.3356\},\\ K^{2}_{3} & = diag\{-21.9717,-19.4468,-73.1039,-65.5740\}. \end{align*} |
In the given example, we set the initial conditions as follows: x_{1}(0) = \left[-1, 1.5, 3\right]^T, \; x_{2}(0) = \left[2.5, -1.2, 1\right]^T\; x_{3}(0) = \left[3.5, -1, 0.5\right]^T, \; x_{4}(0) = \left[-2, -3, 4\right]^T, and s(0) = \left[3.7, 2.4, -1.3\right]^{T}. Based on the observation from Figure 3, it is evident that the error system (2.7) cannot achieve synchronization without using a non-fragile controller. However, Figure 4 illustrates that the trajectory of the error system converges toward zero. The application of the non-fragile controller enables the error to converge to zero, facilitating synchronization in the T-S fuzzy complex dynamical networks. The control input u(t) for the synchronization system is depicted in Figures 5 and 6, where the trajectories converge to zero. To provide more evidence of the effectiveness of the proposed approach, we define e(t) = \sqrt{ \sum\limits_{i = 1}^{4}\Vert x_{i}(t)-s(t)\Vert^2} as the overall synchronization error between the network and the target node. Figure 7 demonstrates that the trajectories of the synchronization error can reach a convergence to 0.
Remark 7. We choose continuous time-delay functions, denoted as \gamma(t), which fulfill the criterion \gamma_{1} \leq \gamma(t) \leq \gamma_{2}. In the context of Example 1, we adopt the specific function \gamma(t) = 0.2\vert \tanh{t}\vert. As illustrated in Figure 8, the time-delay functions \gamma(t) that we employ are continuous and adhere to the condition \gamma_{1} \leq \gamma(t) \leq \gamma_{2}. Notably, our chosen delay function is not required to be differentiable. This contrasts with prior research [4,6,8,34,39], where the time-delay function was consistently differentiable.
Example 2. To implement the networks in practical applications, Chua's circuit system is carefully chosen as the isolated node within the networks to demonstrate the effectiveness of the outcome. This can be explained as follows:
\begin{align*} \dot{s}_{1}(t) & = \sigma_{1}(-s_{1}(t)+s_{2}(t)-\varphi(s_{1}(t)))\\ \dot{s}_{2}(t) & = s_{1}(t) - s_{2}(t) + s_{3}(t)\\ \dot{s}_{3}(t) & = -\sigma_{2}s_{2}(t), \end{align*} |
where \sigma_{1} = 10, \sigma_{2} = 14.87, \varphi(s_{1}(t)) = \eta_{1}s_{1}(t) + 0.5(\eta_{1}-\eta_{2})(\vert s_{1}(t) + 1\vert-\vert s_{1}(t) - 1\vert), \; \eta_{1} = -1.27 and \eta_{2} = -0.68. Let s(t) = \left[s_{1}(t), s_{2}(t), s_{3}(t)\right]^{T} . As a result, the function f(s(t)) can be expressed as the following matrix equation:
\begin{align*} f(s(t)) = \begin{bmatrix} -\eta_{1}-\eta_{1}\sigma_{1} & \eta_{1} & 0\\ 1 & -1 & 1\\ 0 & -14.87 & 0 \end{bmatrix} + \begin{bmatrix} -0.5(\eta_{1}-\eta_{2})(\vert s_{1}(t) + 1\vert -\vert s_{1}(t) - 1\vert)\\ 0\\ 0 \end{bmatrix}. \end{align*} |
The example focuses on five nodes, where the outer coupling matrix and channel matrices are given as follows:
\begin{align} \mathcal{F} = \begin{bmatrix} -3 & 1 & 1 & 0 & 1\\ 1 & -4 & 1 & 1 & 1\\ 1 & 1 & -3 & 1 & 0\\ 0 & 1 & 1 & -3 & 1\\ 1 & 1 & 0 & 1 & -3 \end{bmatrix}, \end{align} | (4.1) |
\begin{align*} \Upsilon_{12} & = diag\{1,1,1\},\;\Upsilon_{13} = diag\{1,1,1\},\;\Upsilon_{14} = diag\{0,0,0\},\;\Upsilon_{15} = diag\{1,1,1\},\\ \Upsilon_{21} & = diag\{1,1,1\},\;\Upsilon_{23} = diag\{1,1,1\},\;\Upsilon_{24} = diag\{1,1,1\},\;\Upsilon_{25} = diag\{0,0,0\},\\ \Upsilon_{31} & = diag\{1,1,1\},\;\Upsilon_{32} = diag\{1,1,1\},\;\Upsilon_{34} = diag\{1,1,1\},\;\Upsilon_{35} = diag\{0,0,0\},\\ \Upsilon_{41} & = diag\{0,0,0\},\;\Upsilon_{42} = diag\{1,1,1\},\;\Upsilon_{43} = diag\{1,1,1\},\;\Upsilon_{45} = diag\{1,1,1\},\\ \Upsilon_{51} & = diag\{1,1,1\},\;\Upsilon_{52} = diag\{0,0,0\},\;\Upsilon_{53} = diag\{0,0,0\},\;\Upsilon_{54} = diag\{1,1,1\}, \end{align*} |
from which one has \mathcal{G}^{k} = \mathcal{F} (k = 1, 2, 3).
Additionally, we consider the perturbation \omega(t) = 1/(10t^2 + 1), and the uncertain parameter matrices are provided as follows: C_{l} = \text{diag}\{1, 1, 1, 1, 1\}, E_{l} = \text{diag}\{0.3, 0.3, 0.3, 0.3, 0.3\}, D_{l}(t) = \text{diag}\{\sin t, \sin t, \sin t, \cos t, \sin t\} and \hat{\Im} = \text{diag}\{0.1, 0.1, 0.1, 0.1, 0.1\}. Furthermore, we set c_{1} = 0, c_{2} = 1, c_{3} = 0, \gamma(t) = 0.04, \varsigma_{1} = \varsigma_{2} = 0 and the performance index for H_{\infty} performance is \delta = 0.3. By utilizing Theorem 2, we can derive the non-fragile control matrices as shown below:
\begin{align*} K^{1}_{1} & = diag\{-269.2818,-321.2147,-269.2818,-269.2818,-269.2818\},\\ K^{1}_{2} & = diag\{-80.4129,-117.1156,-80.4129,-80.4129,-80.4129\},\\ K^{1}_{3} & = diag\{-289.7346,-342.6869,-289.7346,-289.7346,-289.7346\}. \end{align*} |
We select matrices corresponding to each scenario outlined in Table 1 to investigate the extended dissipative behavior. We employed this approach to examine and analyze the extended dissipativity performance in this particular case. Using the Matlab control toolbox, we solved the LMIs stated in Theorem 2. The solutions obtained, ensuring extended dissipative synchronization, are presented in Table 2, encompassing the optimal minimum value \gamma and the maximum value \alpha.
Performance | \Theta_{1} | \Theta_{2} | \Theta_{3} | \Theta_{4} |
H_\infty | -I | 0 | \delta^{2}I | 0 |
L_{2}-L_\infty | 0 | 0 | \delta^{2}I | I |
Passivity | 0 | I | \delta I | 0 |
Dissipativity | -I | I | (2-\alpha)I | 0 |
\gamma_{2}=0.04 | L_{2}-L_{\infty} | H_{\infty} | Passivity | (Q, S, R)-dissipativity |
Performance | 3.1788\times 10^{-10} | 6.7701\times 10^{-11} | 2.3737\times 10^{-11} | 1.9999 |
In the given illustration, we consider the initial conditions: x_{1}(0) = [4, 4, 2]^T, x_{2}(0) = [3, -0.7, 1.5]^T, x_{3}(0) = [2.8, 2.8, -1]^T, x_{4}(0) = [2.9, -0.7, 5.7]^T and x_{5}(0) = [1, 2, 3]^T along with s(0) = [1, -1, 0]^T . In the context of the non-fragile controller, Figure 11 visually demonstrates the convergence of the synchronization error system (2.7) toward zero. Conversely, Figure 10 highlights that the error system (2.7) cannot accomplish synchronization unless a non-fragile controller is utilized. The control input u_{i}(t) , where (i = 1, 2, 3, 4, 5) , is depicted in Figure 12 and consistently exhibits zero trajectories. Furthermore, Figure 9 graphically represents the chaotic behavior of Chua's circuit.
Remark 8. This study introduces a novel criterion for achieving non-fragile control in T-S fuzzy complex networks. Innovative and practical scaling techniques are also applied to address derivatives within the multi-integral Lyapunov functional. Compared with [29,34], we observe that the convergence speed of the state trajectory of the error system under a non-fragile controller in Figure 11 is faster than the findings in [29,34]. In summary, the method proposed in this study offers a combination of reduced computational complexity and less conservatism.
Remark 9. In the past few years, there has been a significant surge in fascination with the management and synchronization of chaotic behaviors in dynamic systems. Chaotic systems exhibit complex behaviors characterized by distinct traits, including their high sensitivity to minor variations in initial conditions and the confinement of their trajectories within phase space. Researchers have delved into the control of chaos within nonlinear systems like Chua's, Lur'e, and Chen's systems. Consequently, Chua's circuit is utilized as unforced, isolated nodes in (2.7) to illustrate its effectiveness through a practical example.
This article has focused on achieving extended dissipative synchronization in T-S fuzzy complex dynamical networks, considering interval hybrid coupling delays and additional disturbances. To begin with, a non-fragile controller was formulated to address the challenge of extended dissipative synchronization. Subsequently, developing a multi-integral Lyapunov functional and applying a mathematical method established novel synchronization criteria for the studied system, ensuring the desired extended dissipative behavior. Finally, the proposed methodology's simulation results validated the synchronization criterion's rationality. Future research might extend these findings to singular systems featuring hybrid time-varying delays.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
This research was supported by the Fundamental Fund of Khon Kaen University. The research on Novel non-fragile extended dissipative synchronization of T-S fuzzy complex dynamical networks with interval hybrid coupling delays by Khon Kaen University has received funding support from the National Science, Research and Innovation Fund or NSRF.
The authors declare that there are no conflicts of interest regarding the publication of this paper.
[1] | M. Dehghan, J. Manafian, A. Saadatmandi, Solving nonlinear fractional partial differential equations using the homotopy analysis method, Num. Meth. Partial Diff. Eq. J., 26 (2010), 448-479. |
[2] | M. Dehghan, J. Manafian, A. Saadatmandi, Application of semi-analytic methods for the Fitzhugh-Nagumo equation, which models the transmission of nerve impulses, Math. Meth. Appl. Sci., 33 (2010), 1384-1398. |
[3] |
M. Dehghan, J. Manafian, A. Saadatmandi, Application of the Exp-function method for solving a partial differential equation arising in biology and population genetics, Int. J. Num. Meth. Heat Fluid Flow, 21 (2011), 736-753. doi: 10.1108/09615531111148482
![]() |
[4] |
X. G. Geng, Y. L. Ma, N-soliton solution and its wronskian form of a (3+1)-dimensional nonlinear evolution equation, Phys. Lett. A, 369 (2007), 285-289. doi: 10.1016/j.physleta.2007.04.099
![]() |
[5] |
J. Manafian, M. Lakestani, Optical solitons with Biswas-Milovic equation for Kerr law nonlinearity, Eur. Phys. J. Plus, 130 (2015), 1-12. doi: 10.1140/epjp/i2015-15001-1
![]() |
[6] |
J. Manafian, On the complex structures of the Biswas-Milovic equation for power, parabolic and dual parabolic law nonlinearities, Eur. Phys. J. Plus, 130 (2015), 1-20. doi: 10.1140/epjp/i2015-15001-1
![]() |
[7] | J. Manafian, M. Lakestani, Solitary wave and periodic wave solutions for Burgers, Fisher, Huxley and combined forms of these equations by the (G'/G)-expansion method, Pramana, 130 (2015), 31-52. |
[8] | J. Manafian, M. Lakestani, New improvement of the expansion methods for solving the generalized Fitzhugh-Nagumo equation with time-dependent coefficients, Int. J. Eng. Math., 2015 (2015), Article ID 107978. |
[9] |
J. Manafian, M. Lakestani, Application of tan(φ/2)-expansion method for solving the Biswas-Milovic equation for Kerr law nonlinearity, Optik, 127 (2016), 2040-2054. doi: 10.1016/j.ijleo.2015.11.078
![]() |
[10] | J. Manafian, M. Lakestani, Dispersive dark optical soliton with Tzitzéica type nonlinear evolution equations arising in nonlinear optics, Opt. Quan. Elec., 48 (2016), 16. |
[11] |
J. Manafian, M. Lakestani, Optical soliton solutions for Schrödinger type nonlinear evolution equations by the tan(φ/2)-expansion method, Optik, 127 (2016), 4222-4245. doi: 10.1016/j.ijleo.2016.01.078
![]() |
[12] | H. M. Baskonus, H. Bulut, Exponential prototype structures for (2+1)-dimensional Boiti-Leon-Pempinelli systems in mathematical physics, Waves in Random and Complex Media, 26 (2016), 201-208. |
[13] | H. M. Baskonus, D. A. Koç, H. Bulut, New travelling wave prototypes to the nonlinear Zakharov-Kuznetsov equation with power law nonlinearity, Nonlinear Sci. Lett. A, 7 (2016), 67-76. |
[14] |
M. Peyrard, A. R. Bishop, Statistical mechanics of a nonlinear model for DNA denaturation, Phys. Rev. Lett., 62 (1989), 2755-2758. doi: 10.1103/PhysRevLett.62.2755
![]() |
[15] |
R. Abazari, S. Jamshidzadeh, G. Wang, Mathematical modeling of DNA vibrational dynamics and its solitary wave solutions, Revista Mexicana de Fisica, 64 (2018), 590-597. doi: 10.31349/RevMexFis.64.590
![]() |
[16] |
S. Dusuel, P. Michaux, M. Remoissenet, From Kinks to compactonlike Kinks, Phys. Rev. E, 57 (1998), 2320-2326. doi: 10.1103/PhysRevE.57.2320
![]() |
[17] |
A. Alvarez, S. R. Romero, J. F. R. Archilla, et al. Breather trapping and breather transmission in a DNA model with an interface, Eur. Phys. J. B, 51 (2006), 119-130. doi: 10.1140/epjb/e2006-00191-0
![]() |
[18] | L. Yakushevich, Nonlinear Physics of DNA, Wiley and Sons, 1998. |
[19] | Mika Gustafsson, Coherent waves in DNA within the Peyrard Bishop model, Master thesis, Linopings Universitet, 2003. |
[20] | E. Villagran, Estructuras solitonicas y su influencia en la dinamica vibraional del ADN, PhD Thesis, Universidad Autonoma del Estado de Mexico, Mexico, 2007. |
[21] |
Z. Wang, B. Zineddin, J. Liang, et al. A novel neural network approach to cDNA microarray image segmentation, Comput Methods Programs Biomed., 111 (2013), 189-198. doi: 10.1016/j.cmpb.2013.03.013
![]() |
[22] |
M. Aguero, M. Najera, M. Carrillo, Nonclassic solitonic structures in DNA's vibrational dynamics, Int. J. Mod. Phys. B, 22 (2008), 2571-2582. doi: 10.1142/S021797920803968X
![]() |
[23] |
G. Gaeta, Results and limitations of the soliton theory of DNA transcription, J. Biol. Phys., 24 (1999), 81-96. doi: 10.1023/A:1005158503806
![]() |
[24] |
J. B. Okaly, A. Mvogo, R. L. Woulaché, et al. Nonlinear dynamics of damped DNA systems with long-rangeinteractions, Commun. Nonlinear Sci. Numer. Simulat., 55 (2018), 183-193. doi: 10.1016/j.cnsns.2017.06.017
![]() |
[25] | G. Miloshevich, J. P. Nguenang, T. Dauxois, et al. Traveling solitons in long-range oscillator chains, J. Phys. A: Math. Theor., 50 (2017), 12LT02. |
[26] |
J. B. Okaly, A. Mvogo, R. L. Woulaché, et al. Semi-discrete breather in a helicoidal DNA double chain-model, Wave Motion, 82 (2018), 1-15. doi: 10.1016/j.wavemoti.2018.06.005
![]() |
[27] | J. B. Okaly, F. L. Ndzana, R. L. Woulaché, et al. Base pairs opening and bubble transport in damped DNA dynamics with transport memory effects, Chaos, 29 (2019), 093103. |
[28] | M. Peyrard, Nonlinear dynamics and statistical physics of DNA, Nonlinearity, 17 (2004), R1-R40. |
[29] |
J. B. Okaly, A. Mvogo, R. L. Woulaché, et al. Nonlinear dynamics of DNA systems with inhomogeneity effects, Chin. J. of Phys., 56 (2018), 2613-2626. doi: 10.1016/j.cjph.2018.07.006
![]() |
[30] | A. Mvogo, G. H. Ben-Bolie, T. C. Kofané, Solitary waves in an inhomogeneous chain of α-helical proteins, Int. J. Mod. Phys B., 28 (2014), 1-14. |
[31] | S. Zdravković, D. Chevizovich, A. N. Bugay, et al. Stationary solitary and kink solutions in the helicoidal Peyrard-Bishop model of DNA molecule, Chaos, 29 (2019), 053118. |
[32] |
S. Zdravković, S. Zeković, Nonlinear dynamics of microtubules and series expansion unknown function method, Chin. J. Phys., 55 (2017), 2400-2406. doi: 10.1016/j.cjph.2017.10.009
![]() |
[33] |
E. Tala-Tebue, Z. I. Djoufack, D. C. Tsobgni-Fozap, et al. Traveling wave solutions along microtubules and in theZhiber-Shabat equation, Chin. J. Phys., 55 (2017), 939-946. doi: 10.1016/j.cjph.2017.03.004
![]() |
[34] | J. Manafian, M. Lakestani, Lump-type solutions and interaction phenomenon to the bidirectional Sawada-Kotera equation, Pramana, 92 (2019), 41. |
[35] |
J. Manafian, Novel solitary wave solutions for the (3+1)-dimensional extended Jimbo-Miwa equations, Comput. Math. Appl., 76 (2018), 1246-1260. doi: 10.1016/j.camwa.2018.06.018
![]() |
[36] | J. Manafian, B. Mohammadi-Ivatlo, M. Abapour, Lump-type solutions and interaction phenomenon to the (2+1)-dimensional Breaking Soliton equation, Appl. Math. Comput., 13 (2019), 13-41. |
[37] |
O. A. Ilhan, J. Manafian, M. Shahriari, Lump wave solutions and the interaction phenomenon for a variable-coefficient Kadomtsev-Petviashvili equation, Comput. Math. Appl., 78 (2019), 2429-2448. doi: 10.1016/j.camwa.2019.03.048
![]() |
[38] |
S. T. R. Rizvi, I. Afzal, K. Ali, et al. Stationary Solutions for Nonlinear Schrödinger Equations by Lie Group Analysis, Acta Physica Polonica A., 136 (2019), 187-189. doi: 10.12693/APhysPolA.136.187
![]() |
[39] | K. Ali, S. T. R. Rizvi, B. Nawaz, et al. Optical solitons for paraxial wave equation in Kerr media, Modern Phys. Let. B, 33 (2019), 1950020. |
[40] | S. Ali, M. Younis, M. O. Ahmad, et al. Rogue wave solutions in nonlinear optics with coupled Schrodinger equations, Opt. Quan. Elec., 50 (2018), 266. |
[41] | A. Arif, M. Younis, M. Imran, et al. Solitons and lump wave solutions to the graphene thermophoretic motion system with a variable heat transmission, Eur. Phys. J. Plus, 134 (2019), 303. |
[42] |
S. Zdravković, Helicoidal PeyrardBishop Model of DNA Dynamics, J. Nonlinear Math. Phys., 18 (2011), 463-484. doi: 10.1142/S1402925111001635
![]() |
[43] |
M. Peyrard, A. R. Bishop, Statistical mechanics of a nonlinear model for DNA denaturation, Phys. Rev. Lett., 62 (1989), 2755-2758. doi: 10.1103/PhysRevLett.62.2755
![]() |
[44] |
T. Dauxois, Dynamics of breather modes in a nonlinear helicoidal model of DNA, Phys. Lett. A, 159 (1991), 390-395. doi: 10.1016/0375-9601(91)90367-H
![]() |
[45] |
M. A. Aguero, M. D. L. Najera, M. Carrillo, Nonclassic solitonic structures in DNA's vibrational dynamics, Int. J. Modern Physics B, 22 (2008), 2571-2582. doi: 10.1142/S021797920803968X
![]() |
[46] | L. Najera, M. Carrillo, M. A. Agüero, Non-classical solitons and the broken hydrogen bonds in DNA vibrational dynamics, Adv. Studies Theor. Phys., 4 (2010), 495-510. |
[47] | S. Zdravković, J. A. Tuszyński, M. V. Satarić, Peyrard-Bishop-Dauxois model of DNA dynamics and impact of viscosity, J. Comput. Theor. Nanosci., 2 (2005), 1-9. |
[48] |
S. Zdravković, M. V. Satarić, Parameter selection in a PeyrardBishopDauxois model for DNA dynamics, Phys. Let. A, 373 (2009), 2739-2745. doi: 10.1016/j.physleta.2009.05.032
![]() |
1. | Yongsheng Rao, Jianwei Su, Behrouz Kheirfam, A Full-Newton Step Interior-Point Method for Weighted Quadratic Programming Based on the Algebraic Equivalent Transformation, 2024, 12, 2227-7390, 1104, 10.3390/math12071104 | |
2. | Tiantian Fan, Jingyong Tang, New smooth weighted complementarity functions and a cubically convergent method for wLCP, 2024, 1862-4472, 10.1007/s11590-024-02139-4 |
Performance | \Theta_{1} | \Theta_{2} | \Theta_{3} | \Theta_{4} |
H_\infty | -I | 0 | \delta^{2}I | 0 |
L_{2}-L_\infty | 0 | 0 | \delta^{2}I | I |
Passivity | 0 | I | \delta I | 0 |
Dissipativity | -I | I | (2-\alpha)I | 0 |
\gamma_{2}=0.04 | L_{2}-L_{\infty} | H_{\infty} | Passivity | (Q, S, R)-dissipativity |
Performance | 3.1788\times 10^{-10} | 6.7701\times 10^{-11} | 2.3737\times 10^{-11} | 1.9999 |
Performance | \Theta_{1} | \Theta_{2} | \Theta_{3} | \Theta_{4} |
H_\infty | -I | 0 | \delta^{2}I | 0 |
L_{2}-L_\infty | 0 | 0 | \delta^{2}I | I |
Passivity | 0 | I | \delta I | 0 |
Dissipativity | -I | I | (2-\alpha)I | 0 |
\gamma_{2}=0.04 | L_{2}-L_{\infty} | H_{\infty} | Passivity | (Q, S, R)-dissipativity |
Performance | 3.1788\times 10^{-10} | 6.7701\times 10^{-11} | 2.3737\times 10^{-11} | 1.9999 |