Next Article in Journal
Six-Dimensional Manifold with Symmetric Signature in a Unified Theory of Gravity and Electromagnetism
Next Article in Special Issue
Algebraic Theory of Patterns as Generalized Symmetries
Previous Article in Journal
Generic Modeling Method of Quasi-One-Dimensional Flow for Aeropropulsion System Test Facility
Previous Article in Special Issue
A New 4D Hyperchaotic System with Dynamics Analysis, Synchronization, and Application to Image Encryption
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamics of Shunting Inhibitory Cellular Neural Networks with Variable Two-Component Passive Decay Rates and Poisson Stable Inputs

by
Marat Akhmet
1,*,
Madina Tleubergenova
2,3 and
Akylbek Zhamanshin
1,2
1
Department of Mathematics, Middle East Technical University, Ankara 06531, Turkey
2
Department of Mathematics, Aktobe Regional University, Aktobe 030000, Kazakhstan
3
Institute of Information and Computational Technologies CS MES RK, Almaty 050010, Kazakhstan
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(6), 1162; https://doi.org/10.3390/sym14061162
Submission received: 10 May 2022 / Revised: 29 May 2022 / Accepted: 1 June 2022 / Published: 5 June 2022
(This article belongs to the Special Issue Symmetry in Nonlinear Dynamics and Chaos)

Abstract

:
Shunting inhibitory cellular neural networks with continuous time-varying rates and inputs are the focus of this research. A new model is considered with compartmental passive decay rates which consist of periodic and Poisson stable components. The first component guarantees the Poisson stability of the dynamics, and the second one causes irregular oscillations. The inputs are Poisson stable to take into account the more sophisticated environment of the networks. The rates and inputs are synchronized to obtain Poisson stable outputs. A new efficient technique for checking the recurrence, the method of included intervals, is applied. Sufficient conditions for the existence of a Poisson stable solution and its asymptotic stability were obtained. Numerical simulations of Poisson stable outputs as well as inputs are provided. Examples of the model with Poisson stable rates, inputs and outputs confirm the feasibility of theoretical results. Discussions were undertaken to provide additional light on the relation of the obtained results with practical and theoretical potentials of neuroscience. Quantitative characteristics are suggested, which can be useful for the future applications of the results. In particular, the center of antisymmetry for the degree of periodicity is determined.

1. Introduction

Cellular neural networks (CNNs), introduced in [1], immediately gained wide recognition to focus on parallel asynchronous computing processes. In recent years, many algorithms for image processing, pattern recognition, evaluation of the dynamics of mechanical systems, etc., were successfully implemented with CNNs.
Shunting inhibitory cellular neural networks (SICNNs) were proposed in [2] and have been extensively applied in psychophysics, optimization, speech, robotics, adaptive pattern recognition, vision and image processing, among other tasks.
In SICNNs, neighboring cells exert mutual inhibitory interactions of the shunting type. In its original formulation, they describe the dynamics of the cell C i j , i = 1 , , m ,   j = 1 , , n , by the following nonlinear ordinary differential equation,
x i j ( t ) = a i j x i j ( t ) C k l N r ( i , j ) C i j k l f ( x k l ( t ) ) x i j ( t ) + v i j ( t ) .
Here, x i j indicates the i j th neuron state, coefficient a i j represents the passive decay rate, C i j k l 0 denotes the strength of the coupling of the cell C k l conveyed to the cell C i j , f ( x k l ) is the activation function, and v i j ( t ) is the external input to cell C i j .
SICNNs are a grid of cells, such that each cell C i j is connected with cells of its r-neighborhood. Here, we follow the definition in [1], where the r-neighborhood N r ( i , j ) is defined as
N r ( i , j ) = { C k l : max ( | k i | , | l j | ) r , 1 k m ,   1 l n } ,
with fixed natural numbers m and n .
It is of strong interest to replace the constant rates in SICNNs by time-varying ones, and many researchers have investigated the dynamic behavior of SICNNs, wherein passive decay rates are periodic [3], anti-periodic [4,5], almost periodic [6,7,8] and pseudo-periodic [9,10] functions.
We investigate the following SICNNs with continuous time-varying coefficients:
x i j ( t ) = ( a i j ( t ) + b i j ( t ) ) x i j ( t ) C k l N r ( i , j ) C i j k l f ( x k l ( t ) ) x i j ( t ) + v i j ( t ) ,
where the rates are compartmental and consist of two components such that a i j ( t ) are periodic whilst b i j ( t ) and inputs v i j ( t ) are Poisson stable functions.
The following concept is fundamental for the present study.
Definition 1
([11]). A function ψ ( t ) : R R , bounded and continuous, is said to be Poisson stable if there is a sequence of moments t p , t p as p , such that the uniform convergence ψ ( t + t p ) ψ ( t ) is valid on each bounded interval of the real axis.
Consider the sequence τ p , p = 1 , 2 , , such that t p = τ p ( m o d ω ) , and ω > 0 is a fixed number. It is clear that 0 τ p < ω . Hence, there exists a subsequence τ p l , l = 1 , 2 , , which converges to a real number τ ω . Consequently, one can find a subsequence t p l , l = 1 , 2 , , of the convergence sequence t p , such that t p l τ ω ( m o d ω ) as l . The number τ ω is said to be Poisson shift for the convergence sequence. The set of all Poisson shifts, T ω , is not empty, and is of several or even infinite number elements. In what follows, κ ω = i n f T ω , 0 κ ω < ω , is called the Poisson number for t p , p = 1 , 2 , , and ω .

2. Methods

It is historical fact that Poisson stability emerged in works [12,13] and then spread to many areas of applied mathematics, similarly to how periodic, quasi-periodic and almost periodic motions [14,15,16,17,18,19] spread. In recent decades, the dynamics were intensively considered in neuroscience [3,6,7,8,20,21,22,23,24]. Currently, Poisson stable motions are the most sophisticated type of recurrence [25,26,27,28].
It is natural to assume that Poisson stable motions should be considered in neuroscience due to the arguments common with those for other types of recurrence, such as periodicity and almost periodicity. The dynamics admit an additional significant reason to be involved in research since it is strongly adjoint with chaos [29].
A new technique, the method of included intervals, was introduced and developed in papers [30,31,32,33,34,35] and book [36] to approve Poisson stability for solutions in differential equations and neural networks. The approach strongly differs from the comparability method in [37,38,39,40].
In papers [33,34,35], it is shown that the Poisson stable motions became chaotic if one equipped them with the unpredictability property. It is a very remarkable result since it was proven that a chaotic trajectory can be Poisson stable and vice versa. Moreover, the dynamics are enriched with infinitely many unstable recurrent motions [36,41,42] which strongly relate to the cognitive activity [43,44,45,46]. This is why one can apply to the dynamics of neural networks effective methods of chaos control such as OGY method and Piragas method, which are very useful in neural models study and their applications [47].
In this paper, SICNNs with Poisson stable passive decay rates and inputs are studied. Moreover, the rates of cell activity have a special structure consisting of two components. One of them is periodic, responsible for the overall dynamics and stability whilst the second is Poisson stable, causing irregular behavior in SICNNs. The passive decay rates and inputs are synchronized by the common moments of convergence on bounded intervals. The relationship between periodic and Poisson stable functions has been studied in [32].

3. Main Results

Using the result for differential equations [48], one can convince that the following lemma is valid.
Lemma 1.
A function y ( t ) = ( y i j ( t ) ) , i = 1 , , m ,   j = 1 , , n , bounded on the real axis, is a solution of (2) if and only if the following equations are true,
y i j ( t ) = t e s t a i j ( u ) d u b i j ( s ) y i j ( s ) + C k l N r ( i , j ) C i j k l f ( y k l ( s ) ) y i j ( s ) v i j ( s ) d s .
Throughout the paper, the norm for bounded functions u : R R , u = sup t R | u ( t ) | , will be used.
Denote by B the set of rectangle matrix-functions ϕ ( t ) = ( ϕ i j ( t ) ) ,   i = 1 , 2 , , m ,   j = 1 , 2 , , n , with Poisson stable entries such that | ϕ i j ( t ) | < H , i = 1 , 2 , , m ,   j = 1 , 2 , , n , where H is a positive number. All the functions are with the common convergence sequence t p ,   p = 1 , 2 , . Moreover, determine the norm ϕ ( t ) 0 = max i , j ϕ i j ( t ) in the set.
Let us define the operator Π ϕ ( t ) = ( Π i j ϕ ( t ) ) , ϕ B , i = 1 , , m ,   j = 1 , , n , where
Π i j ϕ ( t ) t e s t a i j ( u ) d u b i j ( s ) ϕ i j ( s ) + C k l N r ( i , j ) C i j k l f ( ϕ k l ( s ) ) ϕ i j ( s ) v i j ( s ) d s .
We shall need the following assumptions:
(C1)
The functions a i j ( t ) are ω -periodic, such that 0 ω a i j ( u ) d u > 0 , i = 1 , , m , j = 1 , , n ;
(C2)
Functions b i j ( t ) and v i j ( t ) , i = 1 , 2 , , m ,   j = 1 , 2 , , n , are Poisson stable with the common convergence sequence t p , p = 1 , 2 , ;
(C3)
The Poisson number κ ω is equal to zero;
(C4)
m f > 0 such that sup | s | < H | f ( s ) | = m f ;
(C5)
L > 0 such that | f ( s 1 ) f ( s 2 ) | L | s 1 s 2 | if | s 1 | < H , |   s 2 | < H .
Condition ( C 1 ) implies that there exist constants K i j 1 , λ i j > 0 which satisfy
e t s a i j ( u ) d u K i j e λ i j ( t s ) ,
for all i = 1 , , m ,   j = 1 , , n .
We will use notations:
m i j b = sup t R | b i j ( t ) | , m i j v = sup t R | v i j ( t ) | ,   i = 1 , 2 , , m ,   j = 1 , 2 , , n .
The additional assumptions are required.
(C6)
K i j ( m i j b + ( m f + L H ) C k l N r ( i , j ) C i j k l ) < λ i j , i = 1 , , m ,   j = 1 , , n ;
(C7)
K i j m i j v λ i j K i j m f C k l N r ( i , j ) C i j k l K i j m i j b < H , i = 1 , , m ,   j = 1 , , n .
Lemmas A1 and A2 for the following proof are from Appendix A.
Theorem 1.
Suppose that assumptions (C1)–(C7) hold; then, the system (1) admits a unique asymptotically stable Poisson stable solution.
Proof. 
Firstly, we prove the completeness of B . Let ϕ k ( t ) = ( ϕ i j k ( t ) ) , i = 1 , 2 , , m , j = 1 , 2 , , n , be a sequence converging on R to a limit function ϕ ( t ) . Consider a section [ a , b ] , < a < b < , and i = 1 , 2 , , m ,   j = 1 , 2 , , n . We have that
| ϕ i j ( t + t p ) ϕ i j ( t ) | | ϕ i j ( t + t p ) ϕ i j k ( t + t p ) | + | ϕ i j k ( t + t p ) ϕ i j k ( t ) | + | ϕ i j k ( t ) ϕ i j ( t ) | .
For sufficiently large p and k, it can be attained that each term on the right hand-side of (5) is less than ϵ 3 for an arbitrary small ϵ > 0 and t [ a , b ] . The inequalities | ϕ i j ( t + t p ) ϕ i j ( t ) | < ϵ , i = 1 , 2 , , m ,   j = 1 , 2 , , n , imply that ϕ ( t + t p ) ϕ ( t ) uniformly on [ a , b ] . The completeness of B is shown.
Under conditions (C1)–(C7), the operator Π is invariant in B (see Lemma A1) and contractive (see Lemma A2). According to the contraction mapping theorem, there exists a unique Poisson stable solution, z ( t ) = ( z i j ( t ) ) , i = 1 , 2 , , m ,   j = 1 , 2 , , n , of the neural network (2).
Now, consider the stability of solution z ( t ) . It is correct that for all i = 1 , , m , j = 1 , , n ,
z i j ( t ) = e t 0 t a i j ( u ) d u z i j ( t 0 ) t 0 t e s t a i j ( u ) d u b i j ( s ) z i j ( s ) + C k l N r ( i , j ) C i j k l f ( z k l ( s ) ) z i j ( s ) v i j ( s ) d s .
Let x ( t ) = ( x i j ( t ) ) , i = 1 , , m ,   j = 1 , , n , be another solution of model (2). One can write that
x i j ( t ) = e t 0 t a i j ( u ) d u x i j ( t 0 ) t 0 t e s t a i j ( u ) d u b i j ( s ) x i j ( s ) + C k l N r ( i , j ) C i j k l f ( x k l ( s ) ) x i j ( s ) v i j ( s ) d s .
Using the formula
x i j ( t ) z i j ( t ) = e t 0 t a i j ( u ) d u x i j ( t 0 ) z i j ( t 0 ) t 0 t e t 0 t a i j ( u ) d u ( b i j ( s ) x i j ( s ) b i j ( s ) z i j ( s ) + C k l N r ( i , j ) C i j k l f ( x k l ( s ) ) x i j ( s ) C k l N r ( i , j ) C i j k l f ( z k l ( s ) ) z i j ( s ) ) d s = e t 0 t a i j ( u ) d u x i j ( t 0 ) z i j ( t 0 ) t 0 t e t 0 t a i j ( u ) d u ( b i j ( s ) x i j ( s ) b i j ( s ) z i j ( s ) + C k l N r ( i , j ) C i j k l f ( x k l ( s ) ) x i j ( s ) C k l N r ( i , j ) C i j k l f ( z k l ( s ) ) x i j ( s ) + C k l N r ( i , j ) C i j k l f ( z k l ( s ) ) x i j ( s ) C k l N r ( i , j ) C i j k l f ( z k l ( s ) ) z i j ( s ) ) d s ,
we have that
| x i j ( t ) z i j ( t ) | e t 0 t a i j ( u ) d u | x i j ( t 0 ) z i j ( t 0 ) | + t 0 t e t 0 t a i j ( u ) d u ( | b i j ( s ) | | x i j ( s ) z i j ( s ) | + C k l N r ( i , j ) C i j k l | f ( x k l ( s ) ) f ( z k l ( s ) ) | | x i j ( s ) | d s + C k l N r ( i , j ) C i j k l | f ( z k l ( s ) ) | | x i j ( s ) z i j ( s ) | ) d s K i j e λ i j ( t t 0 ) | x i j ( t 0 ) z i j ( t 0 ) | + t 0 t K i j e λ i j ( t t 0 ) m i j b + ( L H + m f ) C k l N r ( i , j ) C i j k l | x i j ( s ) z i j ( s ) | d s ,
for all i = 1 , 2 , , m ,   j = 1 , 2 , , n .
Apply the Gronwall–Bellman Lemma to get
| x i j ( t ) z i j ( t ) | K i j | x i j ( t 0 ) z i j ( t 0 ) | e ( K i j ( m i j b + ( m f + L H ) C k l N r ( i , j ) C i j k l ) λ i j ) ( t t 0 ) ,
for each i = 1 , 2 , , m ,   j = 1 , 2 , , n .
The condition (C6) confirms that the solution z ( t ) is asymptotically stable. □
We shall then illustrate the theoretical results by numerical simulations.
It is known [49] that there exists a Poisson stable sequence ρ i , i Z , which is the solution for the logistic discrete equation
λ i + 1 = μ λ i ( 1 λ i ) ,
where i Z , and μ [ 3 + ( 2 / 3 ) 1 / 2 , 4 ] is a fixed parameter. The sequence is a discrete version of the Poisson stable function. The section [ 0 , 1 ] is invariant with respect to (6) for the considered values of μ .
In the examples below, we will use the Poisson stable function Θ ( t ) [49] which is the unique bounded on the real axis solution of differential equation
Θ = 3 Θ + W ( t ) ,
where W ( t ) = ρ i for t [ h i , h ( i + 1 ) ) , i Z and h is a positive real number. The number h is called the length of step for the piecewise constant function W ( t ) , and for the function Θ ( t ) .
Example 1.
Let us consider the system,
d x i j d t = ( a i j ( t ) + b i j ( t ) ) x i j C k l N 1 ( i , j ) C i j k l f ( x k l ) x i j + v i j ( t ) ,
where i = 1 , 2 , j = 1 , 2 , 3 , f ( s ) = 0.5 a r c t a n ( s ) , C i j 11 = 0.04 , C i j 12 = 0.03 , C i j 13 = 0.05 , C i j 21 = 0.01 , C i j 22 = 0.06 , C i j 23 = 0.02 . For each fixed k = 1 , 2 , and l = 1 , 2 , 3 , the coefficients C i j k l are the same for all i , j , such that the cell C k l belongs to N 1 ( i , j ) . The functions a i j ( t ) are periodic with period ω = π , such that a 11 ( t ) = 3 + cos ( 4 t ) , a 12 ( t ) = 2 + sin ( 2 t ) , a 13 ( t ) = 4 + sin ( 4 t ) , a 21 ( t ) = 5 + cos ( 2 t ) , a 22 ( t ) = 3 + sin ( 6 t ) , a 23 ( t ) = 6 + sin ( 4 t ) . The functions b i j ( t ) and v i j ( t ) are Poisson stable, b 11 ( t ) = 0.5 Θ ( t ) , b 12 ( t ) = 0.3 Θ ( t ) , b 13 ( t ) = 0.2 Θ ( t ) , b 21 ( t ) = 0.7 Θ ( t ) , b 22 ( t ) = 0.8 Θ ( t ) , b 23 ( t ) = 0.3 Θ ( t ) , v 11 ( t ) = 5 Θ ( t ) , v 12 ( t ) = 4 Θ 2 ( t ) , v 13 ( t ) = Θ ( t ) , v 21 ( t ) = 2 Θ ( t ) , v 22 ( t ) = 7 Θ 3 ( t ) , v 23 ( t ) = 2 Θ 2 ( t ) , where Θ ( t ) is the solution of Equation (7) with the step of length h = 3 π . Below, we shall need the degree of periodicity, κ = ω / h , which in the present case is equal to 1 / 3 .
Since the elements of the convergence sequence t i , i Z are multiples of h and the period ω is equal to the Poisson number π , κ ω is equal to zero. Furthermore, sup t R | b i j ( t ) | = m i j b , sup t R | v i j ( t ) | = m i j v , where m 11 b = 0.17 ,   m 12 b = 0.1 ,   m 13 b = 0.07 ,   m 21 b = 0.24 ,   m 22 b = 0.27 , m 23 b = 0.1 ,   m 11 v = 1.67 ,   m 12 v = 0.45 ,   m 13 v = 0.34 ,   m 21 v = 0.67 ,   m 22 v = 0.26 ,   m 23 v = 0.34 . The conditions ( C 1 ) ( C 7 ) of Theorem 1 hold for the network (8) with K i j = 1 , for all i = 1 , 2 , j = 1 , 2 , 3 , L = 0.5 , λ 11 = 3 π ,   λ 12 = 2 π ,   λ 13 = 4 π ,   λ 21 = 5 π ,   λ 22 = 3 π ,   λ 23 = 2 π . As is the case for Theorem 1, Equation (8) admits a unique asymptotically Poisson stable solution z ( t ) = ( z i j ( t ) ) , i = 1 , 2 , j = 1 , 2 , 3 . In Figure 1, the solution x ( t ) = ( x i j ( t ) ) , i = 1 , 2 , j = 1 , 2 , 3 , of SICNNs (8) is shown, which asymptotically converges to z ( t ) .
Figure 2 depicts the coordinates x 11 x 12 x 13 of the solution of (8).
To emphasize the role of the variable passive rates for the behavior of the outputs, in the next simulations, we consider the constant and Poisson stable components of the passive rates. In Figure 3, the outputs of the neural network (8) with a 11 ( t ) + b 11 ( t ) = 3 , a 12 ( t ) + b 12 ( t ) = 2 , a 13 ( t ) + b 13 ( t ) = 4 , a 21 ( t ) + b 21 ( t ) = 5 , a 22 ( t ) + b 22 ( t ) = 3 , a 23 ( t ) + b 23 ( t ) = 6 are demonstrated. One can observe that the constancy is combined with the irregular contribution of the inputs.
Example 2.
Now, we consider SICNNs (8) where i = 1 , 2 , j = 1 , 2 , 3 , activation function f ( s ) = 0.5 a r c t a n ( s ) , the coupling strength C i j 11 = 0.04 , C i j 12 = 0.03 , C i j 13 = 0.05 , C i j 21 = 0.01 , C i j 22 = 0.06 , C i j 23 = 0.02 . For each fixed k = 1 , 2 , and l = 1 , 2 , 3 , the coefficients C i j k l are the same for all i , j , such that the cell C k l belongs to N 1 ( i , j ) . Functions a i j ( t ) are π-periodic such that a 11 ( t ) = 3 + cos ( 4 t ) , a 12 ( t ) = 2 + sin ( 2 t ) , a 13 ( t ) = 4 + sin ( 4 t ) , a 21 ( t ) = 5 + cos ( 2 t ) , a 22 ( t ) = 3 + sin ( 6 t ) , a 23 ( t ) = 6 + sin ( 4 t ) . The functions b i j ( t ) and v i j ( t ) are Poisson stable such that b 11 ( t ) = 0.5 Θ ( t ) , b 12 ( t ) = 0.3 Θ ( t ) , b 13 ( t ) = 0.2 Θ ( t ) , b 21 ( t ) = 0.7 Θ ( t ) , b 22 ( t ) = 0.8 Θ ( t ) , b 23 ( t ) = 0.3 Θ ( t ) , v 11 ( t ) = sin ( 0.2 t ) + 0.5 Θ ( t ) , v 12 ( t ) = sin ( 0.4 t ) + 0.4 Θ 2 ( t ) , v 13 ( t ) = cos ( 0.1 t ) 0.1 Θ ( t ) , v 21 ( t ) = sin ( 0.2 t ) + 0.2 Θ ( t ) , v 22 ( t ) = cos ( 0.2 t ) + 0.7 Θ 3 ( t ) , v 23 ( t ) = cos ( 0.1 t ) + 0.2 Θ 2 ( t ) , where Θ ( t ) is the solution of Equation (7) with the step of length h = π / 10 . The common period of the functions a i j ( t ) and v i j ( t ) equals to 20 π , and the degree of periodicity is equal to κ = 200 .  Figure 4 depicts simulation results of SICNNs (8) with initial values x 11 ( 0 ) = 0.8 , x 12 ( 0 ) = 0.5 , x 13 ( 0 ) = 0.1 , x 21 ( 0 ) = 0.2 , x 22 ( 0 ) = 0.1 , x 23 ( 0 ) = 0.1 . The coordinates x 11 x 12 x 13 of the trajectory with the same initial data are shown in Figure 5.
Comparing Figure 1 and Figure 4, one can understand the importance of the periodicity degree for the models. That is, the periodicity characteristic is seen to be better if the degree is larger than one, and the Poisson stability is better visualized for κ smaller than one. This is why κ = 1 is the center of antisymmetry for degree of periodicity.
Example 3.
Let us consider the following neural network
d y i j d t = ( a i j ( t ) + b i j ( t ) ) y i j C k l N 1 ( i , j ) C i j k l f ( y k l ) y i j + v i j ( t ) ,
where i = 1 , 2 , j = 1 , 2 , 3 , activation function f ( s ) = 0.25 t a n h ( s ) , the coupling strength C i j 11 = 0.003 , C i j 12 = 0.002 , C i j 13 = 0 , C i j 21 = 0.005 , C i j 22 = 0.004 , C i j 23 = 0.001 . For each fixed k = 1 , 2 and l = 1 , 2 , 3 , the coefficients C i j k l are the same for all i , j such that the cell C k l belongs to N 1 ( i , j ) . The functions a i j ( t ) are 2-periodic such that a 11 ( t ) = 4 + cos ( 2 π t ) , a 12 ( t ) = 4 + sin ( 5 π t ) , a 13 ( t ) = 3 + sin ( 4 π t ) , a 21 ( t ) = 5 + cos ( 4 π t ) , a 22 ( t ) = 5 + sin ( 2 π t ) , a 23 ( t ) = 3 + sin ( π t ) . Functions b i j ( t ) and v i j ( t ) are Poisson stable, b 11 ( t ) = 0.1 x 11 ( t ) , b 12 ( t ) = 0.4 x 21 ( t ) , b 13 ( t ) = 0.3 x 13 ( t ) , b 21 ( t ) = 0.6 x 22 ( t ) , b 22 ( t ) = 0.4 x 12 ( t ) , b 23 ( t ) = 0.5 x 23 ( t ) , v 11 ( t ) = 0.2 x 12 ( t ) , v 12 ( t ) = 0.4 x 22 ( t ) , v 13 ( t ) = 0.6 x 21 ( t ) , v 21 ( t ) = 0.2 x 13 ( t ) , v 22 ( t ) = 0.5 x 11 ( t ) , v 23 ( t ) = 0.3 x 22 ( t ) . The function x ( t ) = ( x i j ( t ) ) , i = 1 , 2 , j = 1 , 2 , 3 is the solution of SICNNs (8) obtained for Θ ( t ) with h = 2 . In Figure 6, the solution y ( t ) = ( y i j ( t ) ) , i = 1 , 2 , j = 1 , 2 , 3 , of (9) with y 11 ( 0 ) = 0.05 , y 12 ( 0 ) = 0.01 , y 13 ( 0 ) = 0.04 , y 21 ( 0 ) = 0 , y 22 ( 0 ) = 0.07 , y 23 ( 0 ) = 0.03 , is depicted. In the model, the periodicity degree is equal to 1 and the periodicity is clearly not observable.

4. Discussion

In this part of the paper, we provide the ideas behind the theoretical research of the dynamics of neural networks which are useful for applications in engineering. We found that the simulations of our study are similar to those modeled by neural networks in industrial problems. For example, consider tracking of periodical or square waves dynamics. We are confident that one should better learn what can happen in the dynamics of neural networks with a sophisticated environment and structure. Furthermore, the information of the paper is useful in that sense. One can assume that in addition to the tracking periods, amplitudes of the periodical and square wave dynamics, the neural networks’ motions can also be considered for the Poisson sequences. Thus, the instruments for synchronization and control can be significantly improved, if one consider our proposals.
To shed more light for the suggestions and provide several issues for the next observations, let us look not only at the simulations results for the networks in the examples above but simulations for Poisson stable functions, especially driven to compare with the industrial simulations in papers [50,51,52,53,54]. For this reason, let us consider the two-component function F ( t , κ ) = a cos ( 2 π t ω ) + b Θ ( t ) , where a and b are real coefficients, Θ ( t ) is the Poisson stable function with step h described in the last section, and the first component is an ω -periodic function. The parameter κ is the periodicity degree, which was introduced above in Example 1. It can be very useful for the analysis of time series obtained in different areas—both theoretical and engineering. In Figure 7, Figure 8 and Figure 9, we have the graphs of functions F ( t , 80 ) ,   F ( t , 1 / 5 ) and F ( t , 1 ) , respectively. The graph of function F ( t , 80 ) with the periodicity degree κ > 1 is shown in Figure 7. We obtain the function that admits a clear periodic shape which is enveloped by irregular Poisson stability—which can be said about the dominance of periodicity in this case.
For the function F ( t , 1 / 5 ) , where κ < 1 , and observing its graph in Figure 8, we have the opposite effect, when the irregular stability dominates and the periodicity, which appears locally on separated intervals, envelopes the Poisson stability.
If κ = 1 , then the graph of function F ( t , 1 ) in Figure 9 demonstrates that both phenomena—periodicity and Poisson stability—are equally present in the dynamics.
To estimate the future prospects of our suggestions, the simulations of the paper were compared with some industrial applications of neural networks. First of all, it was found that examples from an historical hourly weather dataset, Figure 9 in [50], had the same shape as the solution x ( t ) = ( x i j ( t ) ) , i = 1 , 2 , j = 1 , 2 , 3 , of neural network (8) with the periodicity degree κ = 200 and the initial values x 11 ( 0 ) = 0.8 , x 12 ( 0 ) = 0.5 , x 13 ( 0 ) = 0.1 , x 21 ( 0 ) = 0.2 , x 22 ( 0 ) = 0.1 , x 23 ( 0 ) = 0.1 in Figure 4. Moreover, it is very similar to the graph of function F ( t , 80 ) in Figure 7. It is remarkable that both mathematical objects are with κ > 1 . It is important that the dominance is numerically indicated for the value of the periodicity degree larger than one.
Our observations then relate to the values of the degree less than one. In Figure 8 of the paper [51] (track control of 1 Hz square wave (no load)), one can see oscillations similar to Poisson stable motions with κ < 1 in Figure 1 and Figure 8, which are obtained on the basis of our theoretical issues. Analyzing simulations with a periodicity degree smaller than one, one can say that there is a loss of periodicity and dominance of irregularity, which is specified as Poisson stability, while the result in [51] might be accepted as coexistence of the square wave periodicity and the Poisson stability in our comprehension. Nevertheless, we believe that one can use the present theoretical results for qualitative characteristics of the square wave tracking with SICNNs [52,53].
Figure 36(a) in [54] is very similar to our simulations in Figure 6 and Figure 9 with the unit degree of periodicity. This confirms that specific Poisson stable motions can be useful to represent the dynamics of linear and nonlinear single-degree-of-freedom systems, as well as full-scale three-story multi-degree of freedom steel frame [54].

5. Conclusions

In this paper, we considered SICNNs with variable two-component rates and Poisson stable inputs. Sufficient conditions were established to provide the existence of Poisson stable solutions for the SICNNs. The line of periodic, almost periodic and recurrent oscillations in neural networks is continued. Introducing the concept of the degree of periodicity, we showed by examples how it affects the dynamics of neural networks. The suggestions are constructive and suitable for simulations.

Author Contributions

M.A.: conceptualization; methodology; investigation. M.T.: investigation; supervision; writing—review and editing. A.Z.: software; investigation; writing—original draft. All authors have read and agreed to the published version of the manuscript.

Funding

M. Tleubergenova and A. Zhamanshin have been supported by the Science Committee of the Ministry of Education and Science of the Republic of Kazakhstan (grant No. AP08856170 and No. AP09258737). M. Akhmet and A. Zhamanshin have been supported by 2247-A National Leading Researchers Program of TUBITAK, Turkey, N 120C138.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors are very thankful to the editor and referees for their valuable comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

In what follows, the basic constructive properties of the operator Π are provided, which are needed to prove Theorem 1.
Lemma A1.
Suppose that assumptions (C1)–(C7) are true. Then, operator Π is invariant in the set B .
Proof. 
Let a function ϕ ( t ) belong to B , then we obtain
| Π i j ϕ ( t ) | | t e s t α i j ( u ) d u b i j ( s ) ϕ i j ( s ) + C k l N r ( i , j ) C i j k l f ( ϕ k l ) ( s ) ϕ i j ( s ) v i j ( s ) d s | t K i j e λ i j ( t s ) | b i j ( s ) | | ϕ i j ( s ) | + C k l N r ( i , j ) C i j k l | f ( ϕ k l ) ( s ) | | ϕ i j ( s ) | + | v i j ( s ) | d s K i j λ i j ( m i j b H + m f H C k l N r ( i , j ) C i j k l + m i j v ) ,
for all i = 1 , 2 , , m ,   j = 1 , 2 , , n . By conditions (C6) and (C7), we obtain that | Π i j ϕ ( t ) | < H , i = 1 , 2 , , m ,   j = 1 , 2 , , n .
We then show that Π i j ϕ ( t ) , i = 1 , 2 , , m , j = 1 , 2 , , n , are Poisson stable functions. Consider a fixed number ϵ > 0 and bounded interval [ α , β ] R . One can choose numbers γ , ξ such that α > γ and ξ > 0 which obey the following relations,
2 K i j λ i j e λ i j ( α γ ) ( m f H C k l N r ( i , j ) C i j k l + m i j b H + m i j v ) < ϵ 4 ,
K i j λ i j ( ( 2 L H 2 + 2 H m f ) C k l N r ( i , j ) C i j k l + 2 m i j b H + 2 m i j v ) e λ i j ( α γ ) < ϵ 4 ,
K i j λ i j ( e ξ ( β γ ) 1 ) ( m f H C k l N r ( i , j ) C i j k l + m i j b H + m i j v ) < ϵ 4 ,
K i j ξ λ i j ( ( L H + m f ) C k l N r ( i , j ) C i j k l + H + m i j b + 1 ) < ϵ 4 ,
for all i = 1 , 2 , , m , j = 1 , 2 , , n . Since the interval [ α , β ] [ γ , β ] , the present technique of discussion is called the method of included intervals.
Because the functions b i j ( t ) and v i j ( t ) , i = 1 , 2 , , m , j = 1 , 2 , , n , are Poisson stable, the function ϕ ( t ) = ( ϕ i j ( t ) ) , i = 1 , 2 , , m , j = 1 , 2 , , n , belongs to B , whilst the convergence sequence, t p , is common for all of them, and the Poisson number κ ω is equal to zero, by Lemma A2 in [32], it is true that | b i j ( t + t p ) b i j ( t ) | < ξ , | v i j ( t + t p ) v i j ( t ) | < ξ , | ϕ i j ( t + t p ) ϕ i j ( t ) | < ξ for t [ γ , β ] , and | a i j ( t + t p ) a i j ( t ) | < ξ for t R , sufficiently large p and i = 1 , 2 , , m ,   j = 1 , 2 , , n . We have that
| Π i j ϕ ( t + t p ) Π i j ϕ ( t ) | | t e s t a i j ( u + t p ) d u ( b i j ( s + t p ) ϕ i j ( s + t p ) + C k l N r ( i , j ) C i j k l f ( ϕ k l ( s + t p ) ) ϕ i j ( s + t p ) v i j ( s + t p ) ) d s + t e s t a i j ( u ) d u b i j ( s ) ϕ i j ( s ) + C k l N r ( i , j ) C i j k l f ( ϕ k l ( s ) ) ϕ i j ( s ) v i j ( s ) d s | t | e s t a i j ( u + t p ) d u e s t a i j ( u ) d u | ( | b i j ( s + t p ) ϕ i j ( s + t p ) + C k l N r ( i , j ) C i j k l f ( ϕ k l ( s + t p ) ) ϕ i j ( s + t p ) v i j ( s + t p ) | ) d s +
t e s t a i j ( u ) d u ( | b i j ( s + t p ) ϕ i j ( s + t p ) b i j ( s ) ϕ i j ( s ) + C k l N r ( i , j ) C i j k l ( f ( ϕ k l ( s + t p ) ) ϕ i j ( s + t p ) f ( ϕ k l ( s ) ) ϕ i j ( s ) ) v i j ( s + t p ) + v i j ( s ) | ) d s .
Let us present the sum of integrals in the last inequality as I 1 + I 2 , such that
I 1 = γ | e s t a i j ( u + t p ) d u e s t a i j ( u ) d u | ( | b i j ( s + t p ) ϕ i j ( s + t p ) + C k l N r ( i , j ) C i j k l f ( ϕ k l ( s + t p ) ) ϕ i j ( s + t p ) v i j ( s + t p ) | ) d s + γ e s t a i j ( u ) d u ( | b i j ( s + t p ) ϕ i j ( s + t p ) b i j ( s ) ϕ i j ( s ) + C k l N r ( i , j ) C i j k l ( f ( ϕ k l ( s + t p ) ) ϕ i j ( s + t p ) f ( ϕ k l ( s ) ) ϕ i j ( s ) ) v i j ( s + t p ) + v i j ( s ) | ) d s ,
and
I 2 = γ t | e s t a i j ( u + t p ) d u e s t a i j ( u ) d u | ( | b i j ( s + t p ) ϕ i j ( s + t p ) + C k l N r ( i , j ) C i j k l f ( ϕ k l ( s + t p ) ) ϕ i j ( s + t p ) v i j ( s + t p ) | ) d s + γ t e s t a i j ( u ) d u ( | b i j ( s + t p ) ϕ i j ( s + t p ) b i j ( s ) ϕ i j ( s ) + C k l N r ( i , j ) C i j k l ( f ( ϕ k l ( s + t p ) ) ϕ i j ( s + t p ) f ( ϕ k l ( s ) ) ϕ i j ( s ) ) v i j ( s + t p ) + v i j ( s ) | ) d s .
One can find that
I 1 γ 2 K i j e λ i j ( t s ) m i j b H + C k l N r ( i , j ) C i j k l m f H + m i j v d s + γ K i j e λ i j ( t s ) 2 m i j b H + C k l N r ( i , j ) C i j k l ( 2 L H 2 + 2 H m f ) + 2 m i j v d s ,
and
I 2 γ t K i j e λ i j ( t s ) ( e ξ ( β γ ) 1 ) m i j b H + C k l N r ( i , j ) C i j k l m f H + m i j v d s + γ t K i j e λ ( t s ) m i j b ξ + H ξ + C k l N r ( i , j ) C i j k l ( L H ξ + m f ξ ) + ξ d s .
Now, applying inequalities (A2)–(A5), we obtain that the following estimates are correct for each i = 1 , 2 , , m ,   j = 1 , 2 , , n :
I 1 2 K i j λ i j e λ i j ( α γ ) ( m i j b H + m f H C k l N r ( i , j ) C i j k l + m i j v ) + K i j λ i j ( ( 2 L H 2 + 2 H m f ) C k l N r ( i , j ) C i j k l + 2 m i j b H + 2 m i j v ) e λ i j ( α γ ) < ϵ 4 + ϵ 4 = ϵ 2 ,
and
I 2 K i j λ i j ( e ξ ( β γ ) 1 ) ( m i j b H + m f H C k l N r ( i , j ) C i j k l + m i j v ) + K i j λ i j ( ( L H ξ + m f ξ ) C k l N r ( i , j ) C i j k l + m i j b ξ + H ξ + ξ ) < ϵ 4 + ϵ 4 = ϵ 2 .
From inequalities (A6) and (A7) it follows that | Π i j ϕ ( t + t p ) Π i j ϕ ( t ) | I 1 + I 2 < ϵ for t [ α , β ] , and all i = 1 , 2 , , m ,   j = 1 , 2 , , n . Therefore, Π B B .
Lemma A2.
Assumptions (C1) and (C4)–(C6) imply that Π is a contraction operator in B .
Proof. 
Let us take two functions φ and ψ that belong to B . It is true that for fixed i and j that the inequality
| Π i j φ ( t ) Π i j ψ ( t ) | t e s t a i j ( u ) d u ( | b i j ( s ) | | φ i j ( s ) ψ i j ( s ) | + C k l N r ( i , j ) C i j k l ( | f ( φ k l ( s ) ) φ i j ( s ) f ( φ k l ( s ) ) ψ i j ( s ) | + | f ( φ k l ( s ) ) ψ i j ( s ) f ( ψ k l ( s ) ) ψ i j ( s ) ) | ) d s K i j λ i j m i j b + ( m f + L H ) C k l N r ( i , j ) C i j k l φ ψ 0
is valid. This is why Π φ Π ψ 0 K i j ( m i j b + ( m f + L H ) C k l N r ( i , j ) C i j k l ) λ i j φ ψ 0 . By condition ( C 6 ) , the operator Π is a contraction in B . □

References

  1. Chua, L.; Yang, L. Cellular neural networks: Applications. IEEE Trans. Circuits Syst. 1988, 35, 1273–1290. [Google Scholar] [CrossRef]
  2. Bouzerdoum, A.; Pinter, R. Shunting inhibitory cellular neural networks: Derivation and stability analysis. IEEE Trans. Circuits Syst. I Fundam. Theory Appl. 1993, 40, 215–221. [Google Scholar] [CrossRef]
  3. Li, Y.; Wang, L.; Fei, Y. Periodic Solutions for Shunting Inhibitory Cellular Neural Networks of Neutral Type with Time-Varying Delays in the Leakage Term on Time Scales. J. Appl. Math. 2014, 2014, 496396. [Google Scholar] [CrossRef]
  4. Huang, C.; Wen, S.; Huang, L. Dynamics of anti–periodic solutions on shunting inhibitory cellular neural networks with multi-proportional delays. Neurocomputing 2019, 357, 47–52. [Google Scholar] [CrossRef]
  5. Peng, G.; Li, L. Anti–periodic solutions for shunting inhibitory cellular neural networks with continuously distributed delays. Nonlinear Anal. Real World Appl. 2009, 10, 2434–2440. [Google Scholar] [CrossRef]
  6. Ou, C. Almost periodic solutions for shunting inhibitory cellular neural networks. Nonlinear Anal. Real World Appl. 2019, 10, 2652–2658. [Google Scholar] [CrossRef]
  7. Li, Y.; Meng, X. Almost periodic solutions for quaternion-valued shunting inhibitory cellular neural networks of neutral type with time delays in the leakage term. Int. J. Syst. Sci. 2018, 49, 2490–2505. [Google Scholar] [CrossRef]
  8. Li, Y.; Wang, C. Almost periodic solutions of shunting inhibitory cellular neural networks on time scales. Commun. Nonlinear Sci. Numer. Simul. 2018, 17, 3258–3266. [Google Scholar] [CrossRef]
  9. Lu, Y.; Ji, D. Pseudo almost periodic solutions for shunting inhibitory cellular neural networks with continuously distributed delays. J. Inequalities Appl. 2017, 2017, 242. [Google Scholar] [CrossRef] [Green Version]
  10. Zhang, A. Pseudo Almost Periodic Solutions for SICNNs with Oscillating Leakage Coefficients and Complex Deviating Arguments. Neural Process. Lett. 2017, 45, 183–196. [Google Scholar] [CrossRef]
  11. Sell, G. Topological Dynamics and Ordinary Differential Equations; Van Nostrand Reinhold Company: London, UK, 1971. [Google Scholar]
  12. Poincare, H. New Methods of Celestial Mechanics; Dover Publications: New York, NY, USA, 1957. [Google Scholar]
  13. Birkhoff, G. Dynamical Systems; American Mathematical Society: Providence, RI, USA, 1927. [Google Scholar]
  14. Nayfeh, A.; Mook, D. Nonlinear Oscillations; Wiley: New York, NY, USA, 1979. [Google Scholar]
  15. Hoppensteadt, F.; Izhikevich, E. Weakly Connected Neural Networks; Springer: New York, NY, USA, 1997. [Google Scholar]
  16. Fink, A. Almost Periodic Differential Equations; Springer: New York, NY, USA, 1974. [Google Scholar]
  17. Levitan, B.; Zhikov, V. Almost Periodic Functions and Differential Equations; Cambridge University Press: Cambridge, UK, 1983. [Google Scholar]
  18. Corduneanu, C. Almost Periodic Oscillations and Waves; Springer: New York, NY, USA, 2009. [Google Scholar]
  19. Besicovitch, A. Almost Periodic Functions; Dover: Cambridge, UK, 1954. [Google Scholar]
  20. Liu, B. Almost periodic solutions for Hopfield neural networks with continuously distributed delays. Math. Comput. Simul. 2007, 73, 327–335. [Google Scholar] [CrossRef]
  21. Cao, J. New results concerning exponential stability and periodic solutions of delayed cellular neural networks. Phys. Lett. A 2003, 307, 136–147. [Google Scholar] [CrossRef]
  22. Zhang, L.; Du, B. Periodic solution for inertial neural networks with variable parameters. AIMS Math. 2021, 6, 13580–13591. [Google Scholar] [CrossRef]
  23. Xu, C.; Pang, Y.; Li, P. Anti–periodic solutions of Cohen-Grossberg shunting inhibitory cellular neural networks on time scales. J. Nonlinear Sci. Appl. 2016, 9, 2376–2388. [Google Scholar] [CrossRef] [Green Version]
  24. Huang, Z. Almost periodic solutions for fuzzy cellular neural networks with multi-proportional delays. Int. J. Mach. Learn. Cybern. 2017, 8, 1323–1331. [Google Scholar] [CrossRef]
  25. Bender, P. Recurrent solutions to systems of ordinary differential equations. J. Differ. Equ. 1969, 5, 271–282. [Google Scholar] [CrossRef] [Green Version]
  26. Kumar, A.; Bhagat, R. Poisson stability in product of dynamical systems. Int. J. Math. Math. Sci. 1987, 10, 613–614. [Google Scholar] [CrossRef] [Green Version]
  27. Hino, Y. Recurrent solutions for linear almost periodic systems. Funkc. Ekvacioj 1985, 28, 117–119. [Google Scholar]
  28. Knight, R. Recurrent and Poisson stable flows. Proc. Am. Math. Soc. 1981, 83, 49–53. [Google Scholar] [CrossRef]
  29. Holmes, P. Poincare, celestial mechanics, dynamical-systems theory and chaos. Phys. Rep. 1990, 193, 137–163. [Google Scholar] [CrossRef]
  30. Akhmet, M.; Tleubergenova, M.; Zhamanshin, A. Poincare chaos for a hyperbolic quasilinear system. Miskolc Math. Notes 2019, 20, 33–44. [Google Scholar] [CrossRef]
  31. Akhmet, M.; Tleubergenova, M.; Zhamanshin, A. Quasilinear differential equations with strongly unpredictable solutions. Carpathian J. Math. 2020, 36, 341–349. [Google Scholar] [CrossRef]
  32. Akhmet, M.; Tleubergenova, M.; Zhamanshin, A. Modulo periodic Poisson stable solutions of quasilinear differential equations. Entropy 2021, 23, 1535. [Google Scholar] [CrossRef] [PubMed]
  33. Akhmet, M.; Seilova, R.; Tleubergenova, M.; Zhamanshin, A. Shunting inhibitory cellular neural networks with strongly unpredictable oscillations. Commun. Nonlinear Sci. Numer. Simul. 2020, 89, 105287. [Google Scholar] [CrossRef]
  34. Akhmet, M.; Tleubergenova, M.; Zhamanshin, A. Inertial neural networks with unpredictable oscillations. Mathematics 2020, 8, 1797. [Google Scholar] [CrossRef]
  35. Akhmet, M.; ÇinÇin, D.A.; Tleubergenova, M.; Nugayeva, Z. Unpredictable oscillations for Hopfield–type neural networks with delayed and advanced arguments. Mathematics 2020, 9, 571. [Google Scholar] [CrossRef]
  36. Akhmet, M. Domain Structured Dynamics: Unpredictability, Chaos, Randomness, Fractals, Differential Equations and Neural Networks; IOP Publishing: Bristol, UK, 2021. [Google Scholar]
  37. Shcherbakov, B.A. Classification of Poisson—stable motions. Pseudo—recurrent motions. Dokl Akad Nauk. SSSR 1962, 146, 322–324. [Google Scholar]
  38. Shcherbakov, B.A. Poisson stable solutions of differential equations, and topological dynamics. Differ. Uravn 1969, 5, 2144–2155. [Google Scholar]
  39. Cheban, D.; Liu, Z. Poisson stable motions of monotone nonautonomous dynamical systems. Sci. China Math. 2019, 62, 1391–1418. [Google Scholar] [CrossRef] [Green Version]
  40. Cheban, D.; Liu, Z. Periodic, quasi–periodic, almost periodic, almost automorphic, Birkhoff recurrent and Poisson stable solutions for stochastic differential equations. J. Differ. Equ. 2020, 269, 3652–3685. [Google Scholar] [CrossRef] [Green Version]
  41. Alligood, K.; Sauer, T.; Yorke, J. CHAOS: An Introduction to Dynamical Systems; Springer: New York, NY, USA, 1996. [Google Scholar]
  42. Devaney, R. An Introduction to Chaotic Dynamical Systems; Addison–Wesley: Menlo Park, CA, USA, 1990. [Google Scholar]
  43. Skarda, C.; Freeman, W.J. How brains make chaos in order to make sense of the world? Behav. Brain Sci. 1987, 10, 161–173. [Google Scholar] [CrossRef] [Green Version]
  44. Boccaletti, S.; Grebogi, C.; Lai, Y.; Mancini, H.; Maza, D. The control of chaos: Theory and applications. Phys. Rep. 2000, 329, 103–197. [Google Scholar] [CrossRef] [Green Version]
  45. Wright, J. Simulation of EEG: Dynamic changes in synaptic efficacy, cerebral rhythms, and dissipative and generative activity in cortex. Biol. Cybern. 1999, 81, 131–147. [Google Scholar] [CrossRef]
  46. Iasemidis, L.; Sackellares, J. REVIEW: Chaos Theory and Epilepsy. Neuroscientist 1996, 2, 118–126. [Google Scholar] [CrossRef]
  47. Hoff, A. Chaos Control and Neural Classification. Z. Naturforschung A 1994, 49, 589–593. [Google Scholar] [CrossRef]
  48. Hartman, P. Ordinary Differential Equations; Birkhauser: Boston, MA, USA, 2002. [Google Scholar]
  49. Akhmet, M.; Fen, M. Poincare chaos and unpredictable functions. Commun. Nonlinear Sci. Numer. Simul. 2017, 48, 85–94. [Google Scholar] [CrossRef] [Green Version]
  50. Zhang, H.; Lu, H.; Nayak, A. Periodic time series data analysis by deep learning methodology. IEEE Access 2020, 8, 078–088. [Google Scholar] [CrossRef]
  51. Li, X.; Xu, F.; Zhang, J.; Wang, S. A multilayer feed forward small-world neural network controller and its application on electrohydraulic actuation system. J. Appl. Math. 2013, 1, 211–244. [Google Scholar] [CrossRef]
  52. Wei, Y.; Zhang, Q. Square wave analysis. In Common Waveform Analysis; Cai, K.Y., Ed.; Springer: Boston, MA, USA, 2000; pp. 13–40. [Google Scholar]
  53. Mohammad, U.; Yasin, M.; Yousuf, R.; Anwar, I. A novel square wave generator based on the translinear circuit scheme of second generation current controlled current conveyor–CCCII. SN Appl. Sci. 2019, 1, 587. [Google Scholar] [CrossRef] [Green Version]
  54. Wu, U.; Yasin, M.; Yousuf, R.; Anwar, I. Deep convolutional neural network for structural dynamic response estimation and system identification. J. Eng. Mech. Appl. Sci. 2019, 145, 04018125. [Google Scholar] [CrossRef]
Figure 1. The solution x ( t ) = ( x i j ( t ) ) , i = 1 , 2 , j = 1 , 2 , 3 , of the neural network (8) with the degree of periodicity κ = 1 / 3 and the initial values x 11 ( 0 ) = 0.8 , x 12 ( 0 ) = 0.5 , x 13 ( 0 ) = 0.1 , x 21 ( 0 ) = 0.2 , x 22 ( 0 ) = 0.1 , x 23 ( 0 ) = 0.1 .
Figure 1. The solution x ( t ) = ( x i j ( t ) ) , i = 1 , 2 , j = 1 , 2 , 3 , of the neural network (8) with the degree of periodicity κ = 1 / 3 and the initial values x 11 ( 0 ) = 0.8 , x 12 ( 0 ) = 0.5 , x 13 ( 0 ) = 0.1 , x 21 ( 0 ) = 0.2 , x 22 ( 0 ) = 0.1 , x 23 ( 0 ) = 0.1 .
Symmetry 14 01162 g001
Figure 2. The coordinates x 11 x 12 x 13 of the solution x ( t ) = ( x i j ( t ) ) , i = 1 , 2 , j = 1 , 2 , 3 , of the neural network (8) with the degree of periodicity κ = 1 / 3 and the initial values x 11 ( 0 ) = 0.8 , x 12 ( 0 ) = 0.5 , x 13 ( 0 ) = 0.1 , x 21 ( 0 ) = 0.2 , x 22 ( 0 ) = 0.1 , x 23 ( 0 ) = 0.1 .
Figure 2. The coordinates x 11 x 12 x 13 of the solution x ( t ) = ( x i j ( t ) ) , i = 1 , 2 , j = 1 , 2 , 3 , of the neural network (8) with the degree of periodicity κ = 1 / 3 and the initial values x 11 ( 0 ) = 0.8 , x 12 ( 0 ) = 0.5 , x 13 ( 0 ) = 0.1 , x 21 ( 0 ) = 0.2 , x 22 ( 0 ) = 0.1 , x 23 ( 0 ) = 0.1 .
Symmetry 14 01162 g002
Figure 3. The solution x ( t ) = ( x i j ( t ) ) , i = 1 , 2 , j = 1 , 2 , 3 , of the neural network (8) with constant passive decay rates and initial values x 11 ( 0 ) = 0.8 , x 12 ( 0 ) = 0.5 , x 13 ( 0 ) = 0.1 , x 21 ( 0 ) = 0.2 , x 22 ( 0 ) = 0.1 , x 23 ( 0 ) = 0.1 . The degree of periodicity κ = 1 / 3 .
Figure 3. The solution x ( t ) = ( x i j ( t ) ) , i = 1 , 2 , j = 1 , 2 , 3 , of the neural network (8) with constant passive decay rates and initial values x 11 ( 0 ) = 0.8 , x 12 ( 0 ) = 0.5 , x 13 ( 0 ) = 0.1 , x 21 ( 0 ) = 0.2 , x 22 ( 0 ) = 0.1 , x 23 ( 0 ) = 0.1 . The degree of periodicity κ = 1 / 3 .
Symmetry 14 01162 g003
Figure 4. The solution x ( t ) = ( x i j ( t ) ) , i = 1 , 2 , j = 1 , 2 , 3 , of the neural network (8) with the degree of periodicity κ = 200 and the initial values x 11 ( 0 ) = 0.8 , x 12 ( 0 ) = 0.5 , x 13 ( 0 ) = 0.1 , x 21 ( 0 ) = 0.2 , x 22 ( 0 ) = 0.1 , x 23 ( 0 ) = 0.1 .
Figure 4. The solution x ( t ) = ( x i j ( t ) ) , i = 1 , 2 , j = 1 , 2 , 3 , of the neural network (8) with the degree of periodicity κ = 200 and the initial values x 11 ( 0 ) = 0.8 , x 12 ( 0 ) = 0.5 , x 13 ( 0 ) = 0.1 , x 21 ( 0 ) = 0.2 , x 22 ( 0 ) = 0.1 , x 23 ( 0 ) = 0.1 .
Symmetry 14 01162 g004
Figure 5. The coordinates x 11 x 12 x 13 of the trajectory x ( t ) = ( x i j ( t ) ) , i = 1 , 2 , j = 1 , 2 , 3 , of the neural network (8). The degree of periodicity is equal to 200.
Figure 5. The coordinates x 11 x 12 x 13 of the trajectory x ( t ) = ( x i j ( t ) ) , i = 1 , 2 , j = 1 , 2 , 3 , of the neural network (8). The degree of periodicity is equal to 200.
Symmetry 14 01162 g005
Figure 6. The solution y ( t ) = ( y i j ( t ) ) , i = 1 , 2 , j = 1 , 2 , 3 , of the SICNNs (9), with a degree of periodicity κ = 1 and initial data y 11 ( 0 ) = 0.05 , y 12 ( 0 ) = 0.01 , y 13 ( 0 ) = 0.04 , y 21 ( 0 ) = 0 , y 22 ( 0 ) = 0.07 , y 23 ( 0 ) = 0.03 .
Figure 6. The solution y ( t ) = ( y i j ( t ) ) , i = 1 , 2 , j = 1 , 2 , 3 , of the SICNNs (9), with a degree of periodicity κ = 1 and initial data y 11 ( 0 ) = 0.05 , y 12 ( 0 ) = 0.01 , y 13 ( 0 ) = 0.04 , y 21 ( 0 ) = 0 , y 22 ( 0 ) = 0.07 , y 23 ( 0 ) = 0.03 .
Symmetry 14 01162 g006
Figure 7. The graph of function F ( t , 80 ) with a = 0.6 , b = 1 , ω = 40 , κ = 80 .
Figure 7. The graph of function F ( t , 80 ) with a = 0.6 , b = 1 , ω = 40 , κ = 80 .
Symmetry 14 01162 g007
Figure 8. The graph of function F ( t , 1 / 5 ) with a = 0.4 , b = 5 , ω = 2 , κ = 1 / 5 .
Figure 8. The graph of function F ( t , 1 / 5 ) with a = 0.4 , b = 5 , ω = 2 , κ = 1 / 5 .
Symmetry 14 01162 g008
Figure 9. The graph of function F ( t , 1 ) with a = 0.5 , b = 2 , ω = π , κ = 1 .
Figure 9. The graph of function F ( t , 1 ) with a = 0.5 , b = 2 , ω = π , κ = 1 .
Symmetry 14 01162 g009
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Akhmet, M.; Tleubergenova, M.; Zhamanshin, A. Dynamics of Shunting Inhibitory Cellular Neural Networks with Variable Two-Component Passive Decay Rates and Poisson Stable Inputs. Symmetry 2022, 14, 1162. https://doi.org/10.3390/sym14061162

AMA Style

Akhmet M, Tleubergenova M, Zhamanshin A. Dynamics of Shunting Inhibitory Cellular Neural Networks with Variable Two-Component Passive Decay Rates and Poisson Stable Inputs. Symmetry. 2022; 14(6):1162. https://doi.org/10.3390/sym14061162

Chicago/Turabian Style

Akhmet, Marat, Madina Tleubergenova, and Akylbek Zhamanshin. 2022. "Dynamics of Shunting Inhibitory Cellular Neural Networks with Variable Two-Component Passive Decay Rates and Poisson Stable Inputs" Symmetry 14, no. 6: 1162. https://doi.org/10.3390/sym14061162

APA Style

Akhmet, M., Tleubergenova, M., & Zhamanshin, A. (2022). Dynamics of Shunting Inhibitory Cellular Neural Networks with Variable Two-Component Passive Decay Rates and Poisson Stable Inputs. Symmetry, 14(6), 1162. https://doi.org/10.3390/sym14061162

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop