Special Session 66: Advances in discrete-time dynamical systems with applications

Optimal Control Approaches for Managing Infectious Diseases with Behavioral Dynamics

Mo`tassem Al-arydah
Khalifa University
United Arab Emirates
Co-Author(s):    Omar Forrest
Abstract:
This study develops a simplified SVIR (Susceptible, Vaccinated, Infected, Recovered) model to examine COVID-19 transmission, incorporating the nonlinear effects of population caution on infection rates. The model`s validity is confirmed by demonstrating the existence of positive bounded solutions. The basic reproduction number is calculated, and the local stability of both the disease-free equilibrium (DFE) and endemic equilibrium (EE) is assessed, revealing that the EE only exists when the basic reproduction number exceeds a critical threshold. Global stability for both equilibria is established using Lyapunov functions. Additionally, an optimal control strategy for vaccination is proposed, proving its existence and uniqueness, with simulations showing that the strategy effectively minimizes infection rates and associated costs. The impact of integrating public education into the model is also explored, emphasizing its critical role in enhancing vaccination coverage and reducing overall transmission.

Strong local asymptotic stability enhances global stability techniques

Ziyad AlSharawi
American University of Sharjah and Universidad Politecnica de Cartagena
United Arab Emirates
Co-Author(s):    Ziyad AlSharawi, Jose Canovas and Sadok Kallel
Abstract:
In this talk, we consider the scalar difference equation $x_{n+1}=F_0(x_n,\ldots,x_{n-k+1})$ in which $F_0$ is sufficiently differentiable and has a fixed point $\bar x.$ When $F_0(x_{n-1},\ldots,x_{n-k})$ replaces $x_n,$ we obtain a new system with higher delay, namely $$y_{n+1}=F_1(y_{n-1},\ldots,y_{n-k})=F_0(F_0(y_{n-1},\ldots,y_{n-k}),y_{n-1},\ldots,y_{n-k+1}).$$ The authors define the expansion strategy as successively repeating the above process, i.e., repeating the process $j$-times gives the system $$u_{n+1}=F_j(u_{n-j},u_{n-j-1},\ldots,u_{n-j-k+1}).$$ The fixed point $\bar x$ of $F_0$ is a fixed point of $F_j$ for all $j.$ In this talk, we discuss the relationship between the local stability of $F_0$ and $F_j$ at $\bar x.$ In particular, $\bar x$ is locally asymptotically stable (LAS) under $F_0,$ if and only if $\|\nabla F_j\|_1

Enveloping in difference equations of order greater than one

Jose S Canovas
Department of Applied Mathematics and Statistics, Technical University of Cartagena
Spain
Co-Author(s):    Ziyad AlSharawi
Abstract:
In this talk, we consider the scalar difference equation $x_{n+1}=f(x_n,x_{n-1})$ in which $f$ is sufficiently differentiable and has a fixed point $x_0.$ We say that $x_0$ is strongly locally stable (SLAS for short) whenever the sum of the absolute values of $f_x(x_0, x_0)$ and $f_y(x_0, x_0)$ is smaller than $1$. We show the advantages of the SLAS on developing the notion of dominance condition as introduced in [1] and the similar enveloping notion for one-dimensional difference equations given in [2]. We give some geometric results that make finding an enveloping one-dimensional map more practical, which makes proving global stability a manageable task. When the map $f$ is of mixed monotonicity, we establish the connection between the enveloping and the embedding techniques. In particular, we prove that embedding is enough to give the existence of an enveloping function and provide ideas for finding enveloping functions. Some interesting questions will be posed. References. [1] H. A. El-Morshedy and V. Jimenez Lopez. Global attractors for difference equations dominated by one-dimensional maps. J. Difference Equ. Appl., 14(4):391-410, 2008 [2] Paul Cull. Enveloping implies global stability. In Difference equations and discrete dynamical systems, pages 71-85. World Sci. Publ., Hackensack, NJ, 2005. [3] Ziyad AlSharawi, Jose S. Canovas. Integrating the expansion strategy with the enveloping technique to establish stability 2024.

The simplest neural network does solve the simplest classification problem

Victor Jimenez Lopez
Universidad de Murcia
Spain
Co-Author(s):    V\`{\i}ctor Jim\`{e}nez L\`{o}pez (coauthor: Jes\`{u}s Molina Rodr\`{\i}guez de Vera)
Abstract:
A single-layer neural network, usually referred to as a \emph{perceptron}, is the most elementary of all machine learning tools to deal with the \emph{classificacion problem}. Roughly speaking, if disjoint finite sets $A,B\subset \mathbb{R}^h$ are given, we look for a (computationally efficient) way to distinguish points $\mathbf{x}=(x_1,\ldots,x_h)$ from $A$ and $B$. More precisely, we look for a \emph{weight vector} $\mathbf{w} =(w_1,\ldots,w_h)$ and a \emph{bias} $b$ so that, when applying the \emph{sigmoid activation function} $\sigma(z)=1/(1+e^{-z})$ to the scalar product of $\overline{\mathbf{x}}=(1,\mathbf{x})$ and $\mathbf{u}=(b,\mathbf{w})$, the resultant function $Z(\mathbf{x})=Z(\mathbf{x},\mathbf{u})=\sigma(\overline{\mathbf{x}}\cdot\mathbf{u})$ takes values close to $1$ (respectively, to $0$) when $\mathbf{x}\in A$ (respectively, $\mathbf{x}\in B$). To do this, and after rewriting $A\cup B=\{\mathbf{x}_k\}_{k=1}^m$, the natural way is to minimize the \emph{squared error function} $E(\mathbf{u})=\sum_{k=1}^m (Z(\overline{\mathbf{x}_k},\mathbf{u})-y_k)^2$, where $y_k=1-\epsilon$ or $y_k=\epsilon$ according to, respectively, $\mathbf{x}_k\in A$ or $\mathbf{x}_k\in B$, and $\epsilon$ is a fixed small positive number (the natural choice $\epsilon=0$ must be avoided for technical reasons; typically, $\epsilon=0.1$ is used). Hopefully, the error function will have exactly one minimum $\mathbf{u}^*=(b^*, \mathbf{w}^*)$ having the prescribed properties. Moreover, it should be found by using the \emph{gradient descent algorithm} with a \emph{learning rate} $\alpha>0$ sufficiently small. In dynamical systems terms, this means that $\textbf{u}_{n+1}=\textbf{f}(\textbf{u}_n)$ given by $\textbf{f}(\textbf{u})=\textbf{u} - \alpha \nabla E(\textbf{u})$ has $\mathbf{u}^*$ as a stable global attractor. While it is very well-known that this dream scenario needs not hold (thus the need to resort to multilayer networks), there are no examples in the literature, to the best of our knowledge, where global convergence of a perceptron, in the previous sense, has been rigorously proved. In the present work we succeed in doing so in the simplest of all cases: when $A=\{1\}$, $B=\{0\}$.

Historic Behavior in Replicator Equations

Mansur Saburov
Department of Mathematics and Natural Science, College of Arts and Sciences (CAS), Center for Applied Mathematics and Bioinformatics (CAMB), Gulf University for Science and Technology (GUST)
Kuwait
Co-Author(s):    Mansur Saburov
Abstract:
Evolutionary game dynamics provide an explanation for how the collective behavior of a large population of players changes over time. This approach consists of two fundamental elements: a population game that describes the strategic interaction that is to occur repeatedly, and a replicator equation that outlines the procedure used by players to determine when and how to adopt new strategies. Drastically different dynamics of zero-sum and nonzero-sum games can be observed under replicator equations. In zero-sum games, heteroclinic cycles occur naturally when species of the population supersede each other in a cyclic fashion, like for the Rock-Paper-Scissors game. The dynamics in the vicinity of a stable heteroclinic cycle is marked by intermittency, where an orbit remains close to the heteroclinic cycle, repeatedly approaching and lingering at the saddles for increasing periods of time, and quickly transitioning from one saddle to the next. This causes the time spent near each saddle to increase at an exponential rate. This highly erratic behavior causes the time averages of the orbit to diverge, a phenomenon known as historic behavior. F.Takens made a prediction in his work that certain dynamical systems found in evolutionary game theory and population dynamics could exhibit persistent historic behavior through the presence of attractive heteroclinic cycles. Recently, the problem of describing persistent families of systems exhibiting historic behavior, known as Takens` Last Problem, has been widely studied in the literature. In this talk, our objective is to confirm Takens` prediction by proposing a persistent and broad class of replicator equations which exhibit historic behavior wherein the slow oscillation of time averages of the orbit ultimately causes the divergence of higher-order repeated time averages.