EDSQ Operator on 2DS and Limit Behavior

Abstract

This paper evaluates the limit behavior for symmetry interactions networks of set points for nonlinear mathematical models. Nonlinear mathematical models are being increasingly applied to most software and engineering machines. That is because the nonlinear mathematical models have proven to be more efficient in processing and producing results. The greatest challenge facing researchers is to build a new nonlinear model that can be applied to different applications. Quadratic stochastic operators (QSO) constitute such a model that has become the focus of interest and is expected to be applicable in many biological and technical applications. In fact, several QSO classes have been investigated based on certain conditions that can also be applied in other applications such as the Extreme Doubly Stochastic Quadratic Operator (EDSQO). This paper studies the behavior limitations of the existing 222 EDSQ operators on two-dimensional simplex (2DS). The created simulation graph shows the limit behavior for each operator. This limit behavior on 2DS can be classified into convergent, periodic, and fixed.

1. Introduction

Discrete dynamical systems have been a subject of interest in many scientific and engineering research applications. The theoretical framework of such systems has been presented in [1,2,3,4,5] using stochastic analysis for population genetics. These dynamical systems consist of set point networks that require symmetrical interactions [6,7]. The work has then been extended in [8] by deploying the majorization technique to derive new notions and properties of the related systems.

A specific case based on linear dynamic systems has been used in [9] to address the agreement issue in the technology science and engineering areas, and has subsequently been reformulated to a nonlinear dynamic system for the same problem [10,11]. In effect, the focus on applications-related studies has shifted to nonlinear dynamical systems in order to address such potential problems. The main reason for this shift is merited by their advantage of being more efficient than their linear counterparts. The quadratic stochastic operator (QSO) constitutes one particular class of nonlinear dynamic systems of particular interest, as introduced in [12], along with its related research history. As a matter of fact, many classes of QSO have been established for the application of consensus [13,14,15,16,17,18,19,20]. This class has been extensively examined and investigated in various works yielding many sub-classes in the process. These QSO sub-classes have been thoroughly investigated as volterra [21,22], dissipative [23,24], Doubly Stochastic Quadratic Operator (DSQOs), and Extreme Doubly Stochastic Quadratic Operator (EDSQOs) [5]. The limit conduct study constitutes the major issue for most nonlinear dynamic models [22]. The mapping of such finite dimensional simplex in this case is what defines the QSO [1]. Lyubich has studied the limit conduct of QSO on one-dimensional simplex (1DS) in [1], albeit with the exception of the high dimensional simplex. Other studies have introduced additional properties for certain QSO sub-classes that are deemed more significant than others [22]. One of these classes is the DSQO.

The points of DSQO are presented in the form of a polyhedron in dimensional space. This polyhedron of DSQO will consist of extreme points which are called EDSQO. Since DSQO constitutes a sub-class of QSO, the QSO is DSQO if the operator of $s$ is majorized by $Vs\prec s$, where $Vs$ is the operator (output), $"\prec "$ signifies majorization, and $s$ is the value of inputs [25]. Consensually, the values of $s$ is arranged in a non-increasing sequence [25]. Mirsky [26] identified the necessary conditions for special stochastic and doubly stochastic matrices. In addition, an overview of the concept of the doubly stochastic matrix along with earlier open problems were addressed in [27]. Ryff [28] extended the related work by introducing the notion of the majorization theory in regard to the convex hull for some functions of doubly stochastic transformations. Moreover, a clear illustration on the significance of these operators and the underlying conceptual perception of majorization were presented in [29], while the theory and necessary conditions for DSQO were described in [23]. Since the dimensional space can be found in finite sets, studying the limit is expected to be rather complicated [23]. Due to this reason, the study on the EDSQOs’ limit behavior has been examined in specific dimensions, in particular from the two-dimensional simplex (2DS) point of view, where such operators are defined using majorization techniques. Such studies have proved the existence of 37 extreme points based on the permutation of the DSQO sets on 2DS. Therefore, it can be used to form 222 EDSQOs on 2DS, as shown in [22]. However, the existing studies have failed to simulate the limit of all 222 EDSQOs on 2DS. Hence, our work is designed to study the limit of these 222 EDSQOs on 2DS and generalize the outcomes for finite-dimensional simplex.

Many previous studies have considered complex non-linear models for synchronization to a complex network. The issues of these models have to do with synchronizing and convergence, as well as in their complex calculations. For example, the Duffiny oscillator [30] involves complex, nonlinear, dynamic equations that convert to convergence and synchronization. However, the downside of this model that it requires highly complicated computations in the form of second-order differential equations as compared to the EDSQ model, which is a more flexible nonlinear model that can deal with a large number of nodes for the networks.

This work presents the limit of all 222 EDSQ operators on 2DS. The paper consists of an introductory part (Section 1), methodology (Section 2), the theoretical result (Section 3), and the discussion and numerical solution (Section 4), followed by the research conclusion (Section 5).

2. Methods

It is appropriate to start this discussion by introducing the concept of majorization and its rules for DSQO operators.

Definition 1:

QSO is a nonlinear dynamical model defined on finite dimensional simplex$\left(m-1\right)$ [14,22,31]:

$$S^{m-1}=\{s=\left(s_1,s_2, \dots ,s_m\right)\in R^m : s_i\ge 0, \mathrm{for}\mathrm{all}\left[1,m\right]$$

(1)

Definition 2:

The set of$s_i$falls within simplex if$s_i>0$. The set consists of$\left\{0,\dots ,1,\dots ,0\right\}$and is called a set of extreme points, while a simplex center set consists of$\left\{\frac{1}{m},\frac{1}{m},\frac{1}{m}\right\}$.

Definition 3:

The nonlinear model of QSO is defined as [32]:

$${\left(Vs\right)}_k={\displaystyle \sum }_{i,j,k=1}^m{P}_{ij,k} s_i{s}_j, s_i>0,$$

(2)

where $p_{ij,k}$ is a transaction matrix under conditions [22,32]:

$$p_{ij,k}=p_{ij,k} \ge 0, {\displaystyle \sum }_{k=1}^m{p}_{ij,k}=1,$$

(3)

The QSO has applications in population genetics, where the population is defined as $m$ species. Here, $s^0=\left(s_1^0,s_2^0, \dots , s_m^0\right)$ is defined as independent members of the population in their initial states. The communication from member $i_{th}$ to member $j_{th}$ called ’parents’ allows the transition to produce a new member $k$. This communication among the members configures a transition matrix $P_{ij,k}$ and produces an inherited factor. Each individual member has a separate transition matrix eventually consisting of the total of these separated matrices with the sum of all its elements being equal to 1. This means that the transition matrix is denoted by $P_{ij,k}$ which has distributed to $m$ matrices $P_{ij,k}=\left(p_{ij,1}\left|p_{ij,2}\right|\dots |p_{ij,m}\right)$, where $p_{ij,m}$ is symmetric and nonnegative.

In this case, a steady state $s=\left(s_1,s_2, \dots ,s_m\right)$ of the potential parents xx appears $s_i{s}_j$ [32].

Definition 4:

The matrix$P_{ij,k}$is doubly stochastic if:

$$p_{ij}\ge 0, and {{\displaystyle \sum }}_{i=1}^m{p}_{ij}={{\displaystyle \sum }}_{j=1}^m{p}_{ij}=1, \mathrm{for}\mathrm{all}i,j=\left[1,m\right]$$

(4)

Definition 5:

The operator$V:S^{m-1}$is called a doubly stochastic operator (DSQO) if:

$$Vs\prec s\mathrm{for}\mathrm{all}s\in S^{m-1}$$

(5)

The set of elements $(s_{\left[1\right]}, \dots , s_{\left[m\right]}$) in DSQO are re-arranged in a non-increasing order $\left(s_{\left[1\right]}\ge \dots \ge s_{\left[m\right]}\right)$, then:

$${\displaystyle \sum }_{i=1}^k{y}_{\left[i\right]} \ge {\displaystyle \sum }_{i=1}^k{s}_{\left[i\right]}, k=\left[1,m\right].$$

(6)

From this point, the majorization concept stipulates that $y$ majorizes $s$, which is denoted by $y\succ s$ or $s$ is majorized by $y$. The doubly stochastic matrix $P$ has been defined according to [31] stating that if $y=Ps$ then $y\succ s$. It follows that in DSQO all $s$ elements belong to the simplex $Ps\succ s$ if $P$ is a doubly stochastic matrix.

Therefore, the particular operator of $V\left(s\right)$ is called DSQO if the transition matrix $P$ is doubly stochastic in accordance with the majorization theory [31].

In other words, the operator of $V\left(s\right)$ is defined a DSQO if the transition matrix $P_{ij,k}$ in Equation (2) satisfies the condition $V\left(s\right)\prec s$, which means:

$$P=\left\{P=\left(a_{ij}\right):a_{ij}=a_{ji}\ge 0, {{\displaystyle \sum }}_{i,j\in \alpha }{a}_{ij}\le \left|\alpha \right|, {{\displaystyle \sum }}_{i,j,k}{a}_{ij}=m\right\}.$$

(7)

where $\alpha$ is a subblock in the matrix that contains elements $\left(m-1\right)$ and its sum should be $\le \left(m-1\right)$ [22].

Definition 6:

As mentioned above, DSQO has vertices points which are called extremes, which is why we call it EDSQO. The additional conditions for DSQO to become EDSQO [15,32] are that the transition matrix should be:

a. 

$p_{ij,k}=0 or \frac{1}{2} or 1;$

b. 

$p_{ii,k}=0 or 1;$

Definition 7:

To define the limit behavior of EDSQO we should first insert the initial values of$s^0\in S^{m-1}$to render it$V\left(s^0\right)$, while$V\left(s^0\right)$is inserted into$V^2\left(s^0\right)$, and so on. The series of$\left\{s^0,V\left(s^0\right),V^2\left(s^0\right), \dots , V^n\left(s^0\right)\right\}$is termed the limit of the trajectory of the operator$V\left(s^0\right)$. Therefore, the QSO, DSQO, and EDSQO can be called the evolutionary process, where$s^0$is the initial generation and$V^{n+1}\left(s^0\right)$the next generation. To study the limit for the operator we find the limit as$n\to \infty$until the trajectory reaches the fixed points. Since the set of points are fixed on$\left(\frac{1}{m},\frac{1}{m},\frac{1}{m}\right)$, the set has converged on the center of the simplex. These sequences of generations are also called iterations.

Definition 8:

If the transition matrix of an individual member$s_d$has interactions only with itself$s_{dd}$and with itself and others$s_{dk}$, then the interaction here is called a selfish interaction.

Definition 9:

If the transition matrix of an individual member$s_d$has at least an interaction between other members$s_{kl}$, then the interaction is called an unselfish interaction.

Definition 10:

If the trajectory of the set points does not change through the limit of$n\to \infty$, then the set points are fixed. Even if the set points start from different values and are fixed after several iterations, the set of points are called fixed where $V^i\left(x^0\right)=x^0$,for all$i=n\to \infty$.

Definition 11:

If the set points of the trajectory have exchanged among themselves through the limit of$n\to \infty$, without changing the values, then the set points are called periodic where$V^i\left(s^0\right)=s^0$,for all$i=n\to \infty$.

3. Theoretical Result

The limit behavior of EDSQO has been studied theoretically in this section.

Theorem 1:

Considering DSQO as an evaluation operator$V$that has been defined on 2DS, there exists a limit of each operator of DSQO on 2DS.

Proof: 

Suppose that $V^{\left(n\right)}\left(s_m\right)$ is EDSQO on 2DS where $n$ is the number of iterations and $m$ is the number of set points, then series of $V^{\left(n\right)}\left(s_m\right)$ is = $V^{\left(1\right)}\left(s_m\right),\dots ,V^{\left(n\right)}\left(s_m\right))$. Based on the majorization concept for DSQO that $s\prec V\left(s\right)$, it means that:

$$V^{\left(0\right)}\left(s_m\right)\succ V^{\left(1\right)}\left(s_m\right)\succ V^{\left(2\right)}\left(s_m\right)\succ \dots \succ V^{\left(n\right)}\left(s_m\right)$$

□

In a non-increasing re-arrangement for $s_m$, we obtain:

$$V^{\left(n\right)}\left(s_1\right)\le \dots \le V^{\left(2\right)}\left(s_1\right)\le V^{\left(1\right)}\left(s_1\right)<V^{\left(0\right)}\left(s_1\right)$$

$$\begin{array}{c}{V}^{\left(n\right)}\left(s_1\right)+V^{\left(n\right)}\left(s_2\right)\le \dots \le V^{\left(2\right)}\left(s_1\right)+V^{\left(2\right)}\left(s_2\right)\le V^{\left(1\right)}\left(s_1\right)+V^{\left(1\right)}\left(s_2\right)<V^{\left(0\right)}\left(s_1\right)+V^{\left(0\right)}\left(s_2\right)\\ ⋮ ⋮ ⋮\\ V^{\left(n\right)}\left(s_1\right)+V^{\left(n\right)}\left(s_2\right)+\dots +V^{\left(n\right)}\left(s_m\right)\le \dots \le V^{\left(2\right)}\left(s_1\right)+V^{\left(2\right)}\left(s_2\right)+\dots +V^{\left(n\right)}\left(s_m\right)\\ \le V^{\left(1\right)}\left(s_1\right)+V^{\left(1\right)}\left(s_2\right)+\dots +V^{\left(n\right)}\left(s_m\right)\le V^{\left(0\right)}\left(s_1\right)+V^{\left(0\right)}\left(s_2\right)+\dots +V^{\left(n\right)}\left(s_m\right)\\ {\displaystyle {\displaystyle \sum }_{k=1}^m}{V}^{\left(n\right)}\left(s\right)\le \dots \le {\displaystyle {\displaystyle \sum }_{k=1}^m}{V}^{\left(2\right)}\left(s\right)\le {\displaystyle {\displaystyle \sum }_{k=1}^m}{V}^{\left(1\right)}\left(s\right)\le {\displaystyle {\displaystyle \sum }_{k=1}^m}s\end{array}$$

Therefore, the maximum of the set points increases, while the minimum of the set points decreases monotonically through the series of $V^{\left(n\right)}\left(s_m\right)$ and becomes fixed.

This means that $li{m}_{n\to \infty }{\displaystyle \sum }_{i=1}^m{V}^{\left(n\right)}\left(s_m\right)$ exists.

We can thus conclude that the limit of EDSQO on 2DS exists.

4. Discussion and Numerical Solution

In this section, we present the simulation results for all operators of EDSQO on 2DS using the MATLAB software.

First, we outline all 37 EDSQOs on 2DS as follows:

$$\begin{array}{cccc}{V}_1\{\begin{array}{c}V\left(s_1\right)=s_1{s}_3+2s_2{s}_3 \\ V\left(s_2\right)=s_1{s}_3+2s_1{s}_2 \\ V\left(s_3\right)=s_1{}^2+s_2{}^2+s_3{}^2\end{array}& V_{11}\{\begin{array}{c}V\left(s_1\right)=s_1{s}_2+2s_2{s}_3 \\ V\left(s_2\right)=s_2{}^2+2s_1{s}_3 \\ V\left(s_3\right)=s_1{}^2+s_3{}^2+s_1{s}_2 \end{array}& V_{21}\{\begin{array}{c}V\left(s_1\right)=s_1{s}_2+2s_1{s}_3 \\ V\left(s_2\right)=s_2{}^2+s_3{}^2+s_1{s}_2 \\ V\left(s_3\right)=s_1{}^2+2s_2{s}_3 \end{array}& V_{31}\{\begin{array}{c}V\left(s_1\right)=s_3{}^2+s_1{s}_2+s_2{s}_3 \\ V\left(s_2\right)=s_2{}^2+2s_1{s}_3 \\ V\left(s_3\right)=s_1{}^2+s_3{s}_2+s_1{s}_2 \end{array}\\ V_2\{\begin{array}{c}V\left(s_1\right)=s_3{}^2+2s_1{s}_2 \\ V\left(s_2\right)=s_1{s}_3+2s_2{s}_3 \\ V\left(s_3\right)=s_1{}^2+s_2{}^2+s_1{s}_3\end{array}& V_{12}\{\begin{array}{c}V\left(s_1\right)=s_1{s}_2+2s_2{s}_3 \\ V\left(s_2\right)=s_2{}^2+s_3{}^2+s_1{s}_3 \\ V\left(s_3\right)=s_1{}^2+s_1{s}_2+s_1{s}_3 \end{array}& V_{22}\{\begin{array}{c}V\left(s_1\right)=2s_1{s}_2+s_2{s}_3 \\ V\left(s_2\right)=s_3{}^{2 }+s_1{s}_3+s_2{s}_3 \\ V\left(s_3\right)=s_1{}^2+s_2{}^2+s_1{s}_3 \end{array}& V_{32}\{\begin{array}{c}V\left(s_1\right)=s_3{}^2+s_1{s}_2+s_2{s}_3 \\ V\left(s_2\right)=s_2{}^2+s_1{s}_3+s_1{s}_2 \\ V\left(s_3\right)=s_1{}^2+s_3{s}_2+s_1{s}_3 \end{array}\\ V_3\{\begin{array}{c}V\left(s_1\right)=s_2{}^2+s_1{s}_2+s_1{s}_3 \\ V\left(s_2\right)=s_1{s}_3+2s_2{s}_3 \\ V\left(s_3\right)=s_1{}^2+s_3{}^2+s_1{s}_2 \end{array}& V_{13}\{\begin{array}{c}V\left(s_1\right)=s_1{s}_2+s_1{s}_3+s_2{s}_3 \\ V\left(s_2\right)=s_3{}^2+s_1{s}_2+s_2{s}_3 \\ V\left(s_3\right)=s_1{}^2+s_2{}^2+s_1{s}_3 \end{array} & V_{23}\{\begin{array}{c}V\left(s_1\right)=2s_1{s}_2+s_2{s}_3 \\ V\left(s_2\right)=s_2{}^{2 }+2s_1{s}_3 \\ V\left(s_3\right)=s_1{}^2+s_3{}^2+s_2{s}_3 \end{array}& V_{33}\{\begin{array}{c}V\left(s_1\right)=s_3{}^2+s_1{s}_2+s_1{s}_3 \\ V\left(s_2\right)=s_2{}^2+s_1{s}_3+s_3{s}_2 \\ V\left(s_3\right)= s_1{}^2+s_3{s}_2+s_1{s}_2 \end{array}\\ V_4\{\begin{array}{c}V\left(s_1\right)=s_2{}^2+s_3{}^2+s_1{s}_2 \\ V\left(s_2\right)=s_1{s}_3+2s_2{s}_3 \\ V\left(s_3\right)=s_1{}^2+s_1{s}_3+s_1{s}_2 \end{array}& V_{14}\{\begin{array}{c}V\left(s_1\right)=s_1{s}_2+s_1{s}_3+s_2{s}_3 \\ V\left(s_2\right)=s_3{}^2+s_1{s}_2+s_1{s}_3 \\ V\left(s_3\right)=s_1{}^2+s_2{}^2+s_2{s}_3 \end{array} & V_{24}\{\begin{array}{c}V\left(s_1\right)=2s_1{s}_2+s_2{s}_3 \\ V\left(s_2\right)=s_2{}^{2 }+s_3{}^2+s_1{s}_3 \\ V\left(s_3\right)=s_1{}^2+s_1{s}_3+s_2{s}_3 \end{array}& V_{34}\{\begin{array}{c}V\left(s_1\right)=s_3{}^2+s_1{s}_2+s_1{s}_3 \\ V\left(s_2\right)=s_2{}^2+s_1{s}_2+s_3{s}_2 \\ V\left(s_3\right)=s_1{}^2+s_1{s}_3+s_2{s}_3 \end{array}\\ V_5\{\begin{array}{c}V\left(s_1\right)=s_2{s}_3+2s_1{s}_3 \\ V\left(s_2\right)=s_2{s}_3+2s_1{s}_2 \\ V\left(s_3\right)=s_1{}^2+s_2{}^2+s_3{}^2 \end{array}& V_{15}\{\begin{array}{c}V\left(s_1\right)=s_1{s}_2+s_1{s}_3+s_2{s}_3 \\ V\left(s_2\right)=s_2{}^2+s_1{s}_3+s_2{s}_3 \\ V\left(s_3\right)=s_1{}^2+s_3{}^2+s_1{s}_2 \end{array}& V_{25}\{\begin{array}{c}V\left(s_1\right)=2s_1{s}_2+s_1{s}_3 \\ V\left(s_2\right)=s_3{}^2+s_1{s}_3+s_3{s}_2\\ V\left(s_3\right)=s_1{}^{2 }+s_2{}^2+s_2{s}_3 \end{array}& V_{35}\{\begin{array}{c}V\left(s_1\right)=s_3{}^2+s_1{s}_2+s_1{s}_3 \\ V\left(s_2\right)=s_2{}^2+s_1{s}_3+s_1{s}_2\\ V\left(s_3\right)=s_1{}^2+2s_2{s}_3 \end{array} \\ V_6\{\begin{array}{c}V\left(s_1\right)=s_2{s}_3+2s_1{s}_3 \\ V\left(s_2\right)=s_3{}^2+2s_1{s}_2 \\ V\left(s_3\right)=s_1{}^2+s_2{}^2+s_2{s}_3 \end{array}& V_{16}\{\begin{array}{c}V\left(s_1\right)=s_1{s}_2+s_1{s}_3+s_2{s}_3 \\ V\left(s_2\right)=s_2{}^2+s_1{s}_2+s_1{s}_3 \\ V\left(s_3\right)=s_1{}^2+s_3{}^2+s_2{s}_3 \end{array}& V_{26}\{\begin{array}{c}V\left(s_1\right)=2s_1{s}_2+s_1{s}_3 \\ V\left(s_2\right)=s_2{}^2+s_1{s}_3+s_2{s}_3 \\ V\left(s_3\right)=s_1{}^{2 }+s_3{}^{2 }+s_2{s}_3 \end{array}& V_{36}\{\begin{array}{c}V\left(s_1\right)=s_3{}^2+2s_1{s}_2 \\ V\left(s_2\right)=s_2{}^2+s_1{s}_3+s_2{s}_3 \\ V\left(s_3\right)=s_1{}^2+s_2{s}_3+s_1{s}_3 \end{array}\\ V_7\{\begin{array}{c}V\left(s_1\right)=s_2{s}_3+2s_1{s}_3 \\ V\left(s_2\right)=s_2{}^2+s_1{s}_2+s_2{s}_3 \\ V\left(s_3\right)=s_1{}^2+s_3{}^2+s_1{s}_2 \end{array}& V_{17}\{\begin{array}{c}V\left(s_1\right)=s_1{s}_2+s_1{s}_3+s_2{s}_3 \\ V\left(s_2\right)=s_2{}^2+s_3{}^2+s_1{s}_3 \\ V\left(s_3\right)=s_1{}^2+s_1{s}_2+s_2{s}_3 \end{array}& V_{27}\{\begin{array}{c}V\left(s_1\right)=2s_1{s}_2+s_1{s}_3 \\ V\left(s_2\right)=s_2{}^2+s_3{}^2+s_1{s}_3 \\ V\left(s_3\right)=s_1{}^2+2s_2{s}_3 \end{array}& V_{37}\{\begin{array}{c}V\left(s_1\right)=s_3{}^2+2s_1{s}_{2 } \\ V\left(s_2\right)=s_2{}^2+2s_1{s}_3 \\ V\left(s_3\right)=s_1{}^2+2s_2{s}_3 \end{array}\\ V_8\{\begin{array}{c}V\left(s_1\right)=s_2{}^2+s_3{}^2+s_1{s}_2 \\ V\left(s_2\right)=s_2{s}_3+2s_1{s}_3 \\ V\left(s_3\right)=s_1{}^2+s_2{s}_3+s_1{s}_2 \end{array}& V_{18}\{\begin{array}{c}V\left(s_1\right)=s_1{s}_2+s_1{s}_3+s_2{s}_3 \\ V\left(s_2\right)=s_2{}^2+s_3{}^2+s_1{s}_2 \\ V\left(s_3\right)=s_1{}^2+s_1{s}_3+s_2{s}_{3 } \end{array}& V_{28}\{\begin{array}{c}V\left(s_1\right)=s_3{}^2+s_1{s}_3+s_2{s}_3 \\ V\left(s_2\right)=s_2{}^2+s_1{s}_2+s_2{s}_{3 } \\ V\left(s_3\right)=s_1{}^2+s_1{s}_2+s_1{s}_3, \end{array}\\ V_9\{\begin{array}{c}V\left(s_1\right)=s_1{s}_2+2s_2{s}_3 \\ V\left(s_2\right)=s_1{s}_2+2s_1{s}_3 \\ V\left(s_3\right)=s_1{}^2+s_2{}^2+s_3{}^2 \end{array}& V_{19}\{\begin{array}{c}V\left(s_1\right)=s_1{s}_2+2s_1{s}_3 \\ V\left(s_2\right)=s_3{}^2+s_1{s}_2+s_2{s}_3 \\ V\left(s_3\right)=s_1{}^2+s_2{}^2+s_2{s}_3 \end{array}& V_{29}\{\begin{array}{c}V\left(s_1\right)=s_3{}^2+s_1{s}_3+s_2{s}_3 \\ V\left(s_2\right)=s_2{}^2+s_1{s}_2+s_1{s}_3 \\ V\left(s_3\right)=s_1{}^2+s_1{s}_2+s_2{s}_3 \end{array}\\ V_{10}\{\begin{array}{c}V\left(s_1\right)=s_1{s}_2+2s_2{s}_3 \\ V\left(s_2\right)=s_3{}^2+s_1{s}_2+s_1{s}_3 \\ V\left(s_3\right)=s_1{}^2+s_2{}^2+s_1{s}_3 \end{array}& V_{20}\{\begin{array}{c}V\left(s_1\right)=s_1{s}_2+2s_1{s}_3 \\ V\left(s_2\right)=s_2{}^2+s_1{s}_2+s_2{s}_3 \\ V\left(s_3\right)=s_1{}^2+s_3{}^2+s_2{s}_3 \end{array}& V_{30}\{\begin{array}{c}V\left(s_1\right)=s_3{}^2+s_1{s}_2+s_2{s}_3 \\ V\left(s_2\right)=s_2{}^2+s_1{s}_3+s_2{s}_3 \\ V\left(s_3\right)=s_1{}^2+s_1{s}_2+s_1{s}_3 \end{array}\end{array}$$

where $V_{1}to{V}_{37}$ are EDSQOs names, and each of the 37 EDSQOs will have $3!$ permuted operators. For instance, for operator $V_1$,

$$\begin{array}{c}{V}_1\{\begin{array}{c}{V}_{1a}\{\begin{array}{c}V\left(s_1\right)=s_1{s}_3+2s_2{s}_3 \\ V\left(s_2\right)=s_1{s}_3+2s_1{s}_2 \\ V\left(s_3\right)=s_1{}^2+s_2{}^2+s_3{}^2 \end{array}\\ V_{1b}\{\begin{array}{c}V\left(s_1\right)=s_1{s}_3+ 2s_2{s}_3 \\ V\left(s_2\right)=s_1{}^2+s_2{}^2+s_3{}^2 \\ V\left(s_3\right)=s_1{s}_3+2s_1{s}_2 \end{array}\\ V_{1c}\{\begin{array}{c}V\left(s_1\right)=s_1{s}_3+2s_1{s}_2 \\ V\left(s_2\right)=s_1{s}_3+ 2s_2{s}_3 \\ V\left(s_3\right)=s_1{}^2+s_2{}^2+s_3{}^2 \end{array}\\ V_{1d}\{\begin{array}{c}V\left(s_1\right)=s_1{s}_3+2s_1{s}_2 \\ V\left(s_2\right)=s_1{}^2+s_2{}^2+s_3{}^2 \\ V\left(s_3\right)=s_1{s}_3+ 2s_2{s}_3 \end{array}\\ V_{1e}\{\begin{array}{c}V\left(s_1\right)=s_1{}^2+s_2{}^2+s_3{}^2 \\ V\left(s_2\right)=s_1{s}_3+ 2s_2{s}_3 \\ V\left(s_3\right)=s_1{s}_3+2s_1{s}_2 \end{array}\\ V_{1f}\{\begin{array}{c}V\left(s_1\right)=s_1{}^2+s_2{}^2+s_3{}^2 \\ V\left(s_2\right)=s_1{s}_3+2s_1{s}_2 \\ V\left(s_3\right)=s_1{s}_3+ 2s_2{s}_3 \end{array}\end{array}\end{array}$$

where $a,b,c,d,e,andf$ are the names of $3!$ permuted EDSQOs. Consequently, we generate $3!\ast 37=222$ EDSQOs. These 222 distributed permuted matrices and operators of 37 EDSQOs have been figures on Appendix A.

The limit behavior of each EDSQO’s permutation is as follows (see Figure 1, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12 and Figure 13):

In this section, we explain the Figure 1, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12 and Figure 13. The initial statuses $\left(s_1,s_2,s_3\right)$ are random, $I$ is the number of iterations for convergence, and $T$ is the computation time to converge. Each figure is assigned a name under EDSQO, followed by the number of operators from $V_{1}to{V}_{37}$, where each operator from $V_{1}to{V}_{37}$ is given one row of Figures $\left(a,b,c,d,e,f\right)$ which are the permutations of 3!.

From the simulation result it can be observed that are altogether 222 EDSQOs, whereby 198 EDSQOs always converge to the average $\left(\frac{1}{3}\right)$ or the center (see Figure 1, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12 and Figure 13). In the case of 18 EDSQOs $\left(V_{4c},V_{4d},V_{7a},V_{7c},V_{12e},V_{12f}, V_{20a},V_{20f},V_{22b},V_{22e},V_{25b},V_{25e},V_{29d},V_{29f},V_{30e},V_{30f},V_{34a},V_{34f}\right)$, it was found that each has one fixed point and two converged points $\left(\frac{s_i^0+s_i^0}{2}\right)$ (see Figure 2, Figure 3 and Figure 4, Figure 7, Figure 8, Figure 9 and Figure 10, Figure 12). Moreover, five EDSQOs $\left(V_{28a},V_{28b},V_{28c},V_{28d},V_{28e}\right)$ had periodic points (see Figure 10), and one EDSQO had only fixed points (see Figure 10).

This simulation analysis proves that the limit behavior of DSQOs in the form of converged or periodic or fixed points. That is due to the fact that if the EDSQO has positive points and none of them has a selfish interaction, these points converge (the 198 EDSQOs). However, if one of these points has a selfish interaction (the 18 EDSQOs and $\left(V_{28f}\right)$), they do not change and are fixed. Consequently, if there are two or more points that have selfish communications, they are periodized (the five EDSQOs).

In conclusion, the study of the limit behavior of all DSQOs on finite dimensional simplex is considerably complex. This is due to the fact that the set points of DSQOs represent a relatively large class. Thus, we have taken a subclass of DSQO that can be studied more conveniently and be used to create a general theory for the entire DSQO class. Therefore, the study has been confined to the EDSQOs class on 2DS due to its low-complexity computation. The limit behavior of EDSQOs on 2DS has been determined and the analysis result has been generalized to the DSQO class.

5. Conclusion and Future Work

We have attempted to analyze a considerably challenging issue related to the nonlinear mathematical model. We have deliberated the trajectories limit behavior of all extreme points for DSQO on 2DS. Subsequently, we have generalized the implications of the result of those EDSQO trajectories’ limit behavior by applying them to the whole class of DSQO. In the process, our work has demonstrated the existence of the EDSQOs’ limit behavior. Through simulation, we have been able to observe the set points’ convergence to the center if they have positive statuses and unselfish communications or their fixation if they have positive statuses and selfish communications or their periodization if they have two or more points with selfish communications. For any future work on this subject, we recommend the study of EDSQ operators on infinite dimensional simplex. Furthermore, it is recommended to evaluate the model for certain specific applications and compare it with existing similar models for the consensus problem.