Search

The study of the dynamics of the size of a population via mathematical modelling is a problem of interest and widely studied. Traditionally, continuous deterministic methods based on differential equations have been used to deal with this problem. However, discrete versions of some models are also available and sometimes more adequate. In this paper, we randomize the Pielou logistic equation in order to include the inherent uncertainty in modelling. Taking advantage of the method of transformation of random variables, we provide a full probabilistic description to the randomized Pielou logistic model via the computation of the probability density functions of the solution stochastic process, the steady state, and the time until a certain level of population is reached. The theoretical results are illustrated by means of two examples: The first one consists of a numerical experiment and the second one shows an application to study the diffusion of a technology using real data. [ABSTRACT FROM AUTHOR]

In this work, a system of non‐linear difference equations is employed to model the opinion dynamics between a small group of agents (the target group) and a very persuasive agent (the indoctrinator). Two scenarios are investigated: the indoctrination of a homogeneous target group, in which each agent grants the same weight to his (or her) partner's opinion and the indoctrination of a heterogenous target group, in which each agent may grant a different weight to his or her partner's opinion. Simulations are performed to study the required times by the indoctrinator to convince a group. Initially, these groups are in a consensus about a doctrine different to that of the ideologist. The interactions between the agents are pairwise. [ABSTRACT FROM AUTHOR]

Some errors exist in the above paper. [ABSTRACT FROM AUTHOR]

The main work of this paper is to investigate two kinds of generalized focal surfaces and two kinds of evolutes generated by spacelike curve γ lying in lightlike surfaces in Minkowski three‐space. Applying the method of unfolding theory in singularity theory to our study, it is shown that there exist the cuspidal edge and the swallowtail types of singularities in each of two classes of generalized focal surfaces under certain conditions; the only cusps will appear in each of evolutes. Two new geometric invariants are presented to classify the singularities of generalized focal surfaces and evolutes. Much more importantly, we reveal the correspondence among the geometric invariants, the types of singularities on generalized focal surfaces and evolutes, the singularities of two kinds of evolutes, and the contact of γ with the osculating spheres. Finally, several examples are presented to demonstrate the correctness of the theoretical results. [ABSTRACT FROM AUTHOR]

The present note presents some errors in the aforementioned paper published in Mathematical Methods in the Applied Sciences. Two errors are found in the definition of the non-dimensional parameters and correct results are presented for temperature profiles included in figure 10 of the previous paper. [ABSTRACT FROM AUTHOR]

The purpose of the present note is to point out mathematical mistakes in some published papers on boundary value problems to prevent usage of mistaken results and methods by mathematicians. Copyright © 2013 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]

There is no correct form of the term HT ht in the literature. In Physics it is not allowed to add quantities with different units and the term HT ht is wrong. [Extracted from the article]

No abstract is available for this article. [ABSTRACT FROM AUTHOR]

Dutch graphic artist M.C. Escher created many famous drawings with a deep mathematical background based on wallpaper symmetry, hyperbolic geometry, spirals, and regular polyhedra. However, he did not attempt any spiral drawings in hyperbolic space. In this paper, we consider a modified hyperbolic geometry by removing the condition that a geodesic is orthogonal to the unit circle in the Poincaré model. We show that spiral symmetry and the similarity property exist in this modified geometry so that the creation of uncommon hyperbolic spiral drawings is possible. To this end, we first establish the theoretical foundation for the proposed method by deriving a contraction mapping and a rotation for constructing modified hyperbolic spiral tilings (MHSTs) and introduce symmetry groups to analyze the structure of MHSTs. Then, to embed a pre‐designed wallpaper template into the tiles, we derive a one‐to‐one mapping between a tile of MHST and a rectangle. Finally, we specify some technical implementation details and give a gallery of the resulting MHST drawings. Using existing wallpaper templates, the proposed method is able to generate a great variety of exotic Escher‐like drawings. [ABSTRACT FROM AUTHOR]

This paper deals with the uniformly stability (US), globally uniformly asymptotically stability (GUAS) and instability of zero solution as well as integrability and boundedness of nonzero solutions of certain nonlinear integro‐differential equations (IDEs). We prove five new theorems, which include sufficient conditions related to these fundamental qualitative properties of solutions to the IDEs considered. The main tools used in the proof are two new and suitable Lyapunov–Krasovskiiˇ functionals (LKFs). In particular cases, two numerical examples are given and solved via the fourth order Runge–Kutta method (RKM) in MATLAB to illustrate the theoretical results of this paper. Compared with the existing results on the fundamental qualitative behaviors of scalar IDEs, our results are new, original, more effective and convenient for tests and applications. [ABSTRACT FROM AUTHOR]

The uniform stabilization of a coupled system arising in the active control of noise in a cavity with a flexible boundary (strings under thermal effects) is considered. Unlike most articles on this subject, which employ the scalar wave equation when analyzing the asymptotic behavior of structural acoustic models, in this paper, we consider classical equations in terms of flow velocity and pressure to describe the acoustic vibrations of the fluid which fills the cavity. This allows to consider, for example, more realistic boundary conditions to model the coupling on the interface between the acoustic chamber and the wall. The main result of this paper, concerning the exponential stability of the model, is established by means of the frequency domain method and the semigroup theory. This method can be adapted to other first‐order hyperbolic dissipative systems as well. [ABSTRACT FROM AUTHOR]

In this paper, we provide an extension result for the existence and uniqueness of solutions for a coupled system of fractional Langevin differential equations with antiperiodic boundary conditions. The system involves two different Caputo fractional derivatives defined on different intervals and associated with boundary conditions described by sequential fractional derivatives. As a conclusion for our main result, we deduce the results of H. Fazli, J.J. Nieto, Fractional Langevin equation with anti‐periodic boundary conditions, Chaos Soliton Fract, 114 (2018),332–337 under less restrictive conditions. The consistency of the main results is demonstrated by two numerical examples. For the sake of completeness, we end the paper by a concluding discussion. [ABSTRACT FROM AUTHOR]

The paper is devoted to the study of 2D version special functions generalizing the famous Slepian functions. As a consequence, many desirable spectral properties of the corresponding weighted finite Fourier transform are deduced from the rich literature. In particular, a similar aspect related to Slepian's seminal papers is the investigation of differential operators that commute with appropriate integral operators. Finally, we provide some analytic expressions for the finite Fourier transforms of the disk polynomials and the two variables Gegenbauer polynomials. [ABSTRACT FROM AUTHOR]

In the present paper, we construct the existence of nontrivial solutions to a new kind of Hadamard fractional boundary value problem on an unbounded domain. With the contribution of some fixed point theorems in cone and the corresponding Green function, we ensure sufficient conditions for the Hadamard fractional boundary value problem. Also, the paper concludes with two examples to demonstrate our results. [ABSTRACT FROM AUTHOR]

The dynamics of the propagation and outspread of infectious diseases are eminently intricate, mainly due to the heterogeneity of the host individuals. In this paper, an age‐stratified SEIR (susceptible‐exposed‐infected‐recovered) epidemiological model incorporating saturated treatment function and heterogeneous contact rates is developed to study infectious disease transmission dynamics among various age groups. The expression for the basic reproduction number R0$$ {R}_0 $$ and conditions for the global stability of the system have been derived by a recently developed graph‐theoretic (GT) approach. Digraph reduction creates a GT version of the Gauss elimination method for computing the R0$$ {R}_0 $$. The global dynamics results have been formed by constructing the Lyapunov function using a GT approach. The endemic equilibrium exists uniquely if R0>1$$ {R}_0>1 $$, whereas the disease‐free equilibrium is observed to be globally stable if R0≤1$$ {R}_0\le 1 $$. The numerical simulations are demonstrated by extracting the daily reported COVID‐19 cases for the second wave in Italy. The age‐dependent contact matrix for the Republic of Italy (data sourced from the POLYMOD study) is computed via paper–diary methodology (PDM) grounded on a population‐prospective survey in European countries. Our numerical findings imply that (i) for the age group (20–49) years and (50–69) years, the number of infected persons is quite double as compared with the exposed cases; (ii) approximately 50% of positive cases lies in (20–69) years age group; (iii) for the (00–19) years age group, only half of the exposed individuals got infected; and (iv) a consistent graph is detected for the age group of (70–99) years in both cases; it shows that almost all the exposed got infected. [ABSTRACT FROM AUTHOR]

This paper explores an inverse problem pertaining to the determination of a source function in non‐divergence parabolic equations, where the solution is known at a discrete set of points. Being different from other ordinary inverse source problems, which are often dependent on only one variable, the unknown coefficient in this paper not only depends on the space variable x$$ x $$ but also depends on the time t$$ t $$. On the basis of the optimal control framework, the existence of the optimal solution of the control function is proved. The necessary conditions to be satisfied by the optimal solution are given. The convergence of the optimal solution when the mesh parameters tend to zero is obtained. The conjugate gradient method is applied to the inverse problem and some numerical results are presented for various typical test examples. [ABSTRACT FROM AUTHOR]

In this paper steady‐state vibrations of the two‐degrees‐of‐freedom oscillatory systems with van der Pol coupling are investigated. The model is a system of two differential equations with weak nonlinearity. A new solving procedure based on D′Alembert's method and the method of time‐variable amplitude and phase is developed. The main advantage of the method in comparison to others is that it gives the solution of the system of two coupled weak nonlinear equations in the form that is simple to analyze, as it has the same form as the solution of the corresponding system of linear equations. In the paper two types of systems are considered: one, a two‐mass system with two degrees of freedom, and second, the one‐mass system with two degrees of freedom. The torsional vibrations of a two‐mass system and vibrations of a Jeffcott rotor with two‐degrees‐of‐freedom are analyzed. Analytically obtained results are numerically tested. It is obtained that the difference between analytic and numeric results is small and almost negligible. As the accuracy of the analytic solution is high, it is suitable for application in technics and engineering. Conclusions about steady‐state self‐sustainable oscillators, orbital, and limit cycle motions are given. [ABSTRACT FROM AUTHOR]

This paper focuses on the three‐dimensional (3D) incompressible anisotropic Boussinesq system while the velocity of fluid only involves horizontal dissipation and the temperature has a damping term. By utilizing the structure of the system, the energy methods and the means of bootstrapping argument, we prove the global stability property in the Sobolev space Hk(ℝ3)(k≥3)$$ {H}^k\left({\mathrm{\mathbb{R}}}^3\right)\left(k\ge 3\right) $$ of perturbations near the hydrostatic equilibrium. Moreover, we take an effective approach to obtain the optimal decay rates for the global solution itself as well as its derivatives. In this paper, we aim to reveal the mechanism of how the temperature helps stabilize the fluid. Additionally, exploring the stability of perturbations near hydrostatic equilibrium may provide valuable insights into specific severe weather phenomena. [ABSTRACT FROM AUTHOR]

In this paper, we study the jerk vector that is the rate of change of the acceleration vector over time. In three‐dimensional space, the decomposition of the jerk vector is a new concept in the field. This decomposition expresses the jerk vector as the sum of three unique components in specific directions: the tangential direction, the radial direction in the osculating plane, and the radial direction in the rectifying plane. The snap vector is the rate of change of the jerk vector over time. In this paper, the authors examine non‐relativistic particles moving along non‐lightlike Frenet curves at low speeds compared to the speed of light in Minkowski 3‐space. They resolve the jerk and snap vectors using Frenet–Serret frames. Additionally, the cases for motion along non‐lightlike Frenet planar curves in the Minkowski 3‐space are given as corollaries. To help understand these results, the paper provides some illustrative examples [ABSTRACT FROM AUTHOR]

In this paper, a Hopfield‐type neural network system with Caputo–Hadamard fractional derivative is discussed. The importance of the existence of the equilibrium point in the analysis of artificial neural networks is well known. Another important investigation is the stability properties. So, the stability of a neural network system is dealt with in the present paper. First, a theorem that asserts the existence and uniqueness of the equilibrium point of the system is proven. Later, the conditions that ensure the Mittag‐Leffler stability of the equilibrium point is obtained by using the Lyapunov's direct method. In addition, an example is given with numerical simulations to show the effectiveness of our theoretical results. [ABSTRACT FROM AUTHOR]

In this paper, we develop and analyze a mathematical model for spreading malaria, including treatment with transmission‐blocking drugs (TBDs). The paper's main aim is to demonstrate the impact the chosen model for demographic growth has on the disease's transmission and the effect of its treatment with TBDs. We calculate the model's control reproduction number and equilibria and perform a global stability analysis of the disease‐free equilibrium point. The mathematical analysis reveals that, depending on the model's demography, the model can exhibit forward, backward, and even some unconventional types of bifurcation, where disease elimination can occur for both small and large values of the reproduction number. We also conduct a numerical analysis to explore the short‐time behavior of the model. A key finding is that for one type of demographic growth, the population experienced a significantly higher disease burden than the others, and when exposed to high levels of treatment with TBDs, only this population succeeded in effectively eliminating the disease within a reasonable timeframe. [ABSTRACT FROM AUTHOR]

The primary goal of this paper is to suggest a fully discrete numerical solution approach for the time‐fractional Burgers' equation. This paper will consider the fractional derivative in the Caputo sense. The time derivative of this equation will be discretized using the L2‐type discretization formula. The spatial variable is approximated by using the nonsymmetric interior penalty discontinuous Galerkin method. The proposed method is globally and unconditionally stable. The accuracy of the solution is evaluated using a convergence analysis. Computational experiments further confirm the accuracy and stability of the suggested strategy. [ABSTRACT FROM AUTHOR]

This paper deals with the probabilistic analysis of a class of compartmental models formulated via a system of linear differential equations with time‐dependent non‐homogeneous terms. For the sake of generality, we assume that initial conditions and rates between compartments are random variables with arbitrary distributions while the source terms defining the flows entering the compartments are stochastic processes. We then take extensive advantage of the so‐called random variable transformation (RVT) technique to determine the first probability distribution of the solution of such a randomized model under very general hypotheses. In the simplest but relevant case of time‐independent source terms, we also obtain the probability distribution of the equilibrium, which is a random vector. Furthermore, we first particularize the aforementioned results for different types of integrable source terms that permit obtaining an explicit solution of the compartmental model: random constants, a train of Dirac delta impulses with random intensities, and a Wiener process. Secondly, we show an alternative approach based on the Liouville–Gibbs equation, which is useful when dealing with source terms that do not admit a closed‐form primitive. All the previous theoretical results are first illustrated through several numerical examples and simulations where a wide range of different probability distributions are assumed for the model parameters. The paper concludes by applying a compartmental model that describes the dynamics of oral drug administration through multiple chronologically spaced doses using synthetic data generated according to pharmacological references. [ABSTRACT FROM AUTHOR]

This paper deals with compatibility of space‐time kernels with (either) full, spatially dynamical, or space‐time compact support. We deal with the dilemma of statistical accuracy versus computational scalability, which are in a notorious trade‐off. Apparently, models with full support ensure maximal information but are computationally expensive, while compactly supported models achieve computational scalability at the expense of loss of information. Hence, an inspection of whether these models might be compatible is necessary. The criterion we use for such an inspection is based on equivalence of Gaussian measures. We provide sufficient conditions for space‐time compatibility. As a corollary, we deduce implications in terms of maximum likelihood estimation and misspecified kriging prediction under fixed domain asymptotics. Some results of independent interest relate about the space‐time spectrum associated with the classes of kernels proposed in the paper. [ABSTRACT FROM AUTHOR]

In this article, we introduce a modified class of Bernstein–Kantorovich operators depending on an integrable function ψα$$ {\psi}_{\alpha } $$ and investigate their approximation properties. By choosing an appropriate function ψα$$ {\psi}_{\alpha } $$, the order of approximation of our operators to a function f$$ f $$ is at least as good as the classical Bernstein–Kantorovich operators on the interval [0,1]$$ \left[0,1\right] $$. We compared the operators defined in this study not only with Bernstein–Kantorovich operators but also with some other Bernstein–Kantorovich type operators. In this paper, we also study the results on the uniform convergence and rate of convergence of these operators in terms of the first‐ and second‐order moduli of continuity, and we prove that our operators have shape‐preserving properties. Finally, some numerical examples which support the results obtained in this paper are provided. [ABSTRACT FROM AUTHOR]

It is well known that investigation on exact solutions of nonlinear fractional partial differential equations (PDEs) is a very difficult work compared with integer‐order nonlinear PDEs. In this paper, based on the separation method of semifixed variables and integral bifurcation method, a combinational method is proposed. By using this new method, a class of generalized time‐fractional thin‐film equations are studied. Under two kinds of definitions of fractional derivatives, exact solutions of two generalized time‐fractional thin‐film equations are investigated respectively. Different kinds of exact solutions are obtained and their dynamic properties are discussed. Compared to the results in the existing references, the types of solutions obtained in this paper are abundant and very different from those in the existing references. Investigation shows that the solutions of the model defined by Riemann–Liouville differential operator converge faster than those defined by Caputo differential operator. It is also found that the profiles of some solutions are very similar to solitons, but they are not true soliton solutions. In order to visually show the dynamic properties of these solutions, the profiles of some representative exact solutions are illustrated by 3D graphs. [ABSTRACT FROM AUTHOR]

In a lot of engineering applications, one has to recover the temperature of a heat body having a surface that cannot be measured directly. This raises an inverse problem of determining the temperature on the surface of a body using interior measurements, which is called the sideways problem. Many papers investigated the problem (with or without fractional derivatives) by using the Cauchy data u(x0,t),ux(x0,t)$$ u\left({x}_0,t\right),{u}_x\left({x}_0,t\right) $$ measured at one interior point x=x0∈(0,L)$$ x={x}_0\in \left(0,L\right) $$ or using an interior data u(x0,t)$$ u\left({x}_0,t\right) $$ and the assumption limx→∞u(x,t)=0$$ {\lim}_{x\to \infty }u\left(x,t\right)=0 $$. However, the flux ux(x0,t)$$ {u}_x\left({x}_0,t\right) $$ is not easy to measure, and the temperature assumption at infinity is inappropriate for a bounded body. Hence, in the present paper, we consider a fractional sideways problem, in which the interior measurements at two interior point, namely, x=1$$ x=1 $$ and x=2$$ x=2 $$, are given by continuous data with deterministic noises and by discrete data contaminated with random noises. We show that the problem is severely ill‐posed and further constructs an a posteriori optimal truncation regularization method for the deterministic data, and we also construct a nonparametric regularization for the discrete random data. Numerical examples show that the proposed method works well. [ABSTRACT FROM AUTHOR]

In this paper, we study a higher order Kirchhoff problem with variable exponent of type M∫Ω|Dru|m(x)m(x)dxΔm(x)ru=f(x,u)inΩ,Dαu=0,on∂Ω,for eachα∈ℝNwith|α|≤r−1,$$ \left\{\begin{array}{ll}M\left({\int}_{\Omega}\frac{{\left|{\mathcal{D}}_ru\right|}^{m(x)}}{m(x)} dx\right){\Delta}_{m(x)}^ru=f\left(x,u\right)& \mathrm{in}\kern0.30em \Omega, \\ {}{D}^{\alpha }u=0,\kern0.30em & \mathrm{on}\kern0.30em \mathrm{\partial \Omega },\kern0.30em \mathrm{for}\ \mathrm{each}\kern0.4em \alpha \in {\mathrm{\mathbb{R}}}^N\kern0.4em \mathrm{with}\kern0.4em \mid \alpha \mid \le r-1,\end{array}\right. $$where Ω⊂ℝN$$ \Omega \subset {\mathrm{\mathbb{R}}}^N $$ is a smooth bounded domain, r∈ℕ∗,m∈C(Ω‾),1

Stochastic differential equations (SDEs) are very useful in modeling many problems in biology, economic data, turbulence, and medicine. Fractional Brownian motion (fBm) and variable‐order fractional Brownian motion (vofBm) are suitable alternatives to standard Brownian motion (sBm) for describing and modeling many phenomena, since the increments of these processes are dependent of the past and for H>12$$ \mathcal{H}>\frac{1}{2} $$ these increments have the property of long‐term dependence. Classical mathematical techniques such as Ito's calculus do not work for stochastic computations on fBm and vofBm due to they are not semi‐Martingale for H(ξ)≠12$$ \mathcal{H}\left(\xi \right)\ne \frac{1}{2} $$. Therefore, solving these equations is much more difficult than solving SDEs with sBm. On the other hand, these equations do not have an analytical solution, so we have to use numerical methods to find their solution. In this paper, a computational approach based on hybrid of block‐pulse and parabolic functions (HBPFs) has been introduced for simulating vofBm and solving a modern class of SDEs. The mechanism of this approach is based on stochastic and fractional integration operational matrices, which transform the intended problem to a nonlinear system of algebraic equations. Thus, the complexity of solving the mentioned problem is reduced significantly. Also, convergence analysis of the expressed method has been theoretically examined. Finally, the accuracy and efficiency of the proposed algorithm have been experimentally investigated through some test problems and comparison of obtained results with results of previous papers. High accurate numerical results are obtained by using a small number of basic functions. Therefore, this method deals with smaller matrices and vectors, which is one of the most important advantage of our suggested method. Also, presenting an applicable procedure to construct vofBm is another innovation of this work. To gain this aim, at first, discretized sBm is generated via fundamental features of this process, and afterward, block‐pulse functions (BPFs) and HBPFs are utilized for simulating discretized vofBm. Finally, spline interpolation method has been employed to provide a continuous path of vofBm. [ABSTRACT FROM AUTHOR]

Geometrically, a complex dynamical network (CDN) can be regarded as the interconnected system composed of the node subsystem (NS) and the link subsystem (LS) coupled with each other. Guided by this idea, in order to achieve the goal of each node following asymptotically its own reference target in a CDN, this paper investigates the model following adaptive control (MFAC) problem of NS via the dynamics of links, which implies that the LS plays the important dynamic auxiliary role in the MFAC realization of nodes. Meanwhile, we focus on the condition that the links state information is unavailable, due to sensor practical application and measurement cost constraints. To obviate this restriction, we construct the asymptotical state observer for the LS. Next, to achieve the control goal of this paper, an appropriate Lyapunov candidate function is selected, by which the MFAC scheme for NS is synthesized based on the state observer of LS. Finally, the simulation example is performed to demonstrate the theoretical results in this paper. [ABSTRACT FROM AUTHOR]

This paper discusses the issue of passivity analysis of the neutral‐type Cohen–Grossberg neural networks (NTCGNNs) with distributed, transmission, and neutral time‐varying delays. A new method based on the system solutions was proposed to analyze passivity. It is important to note that this paper focuses solely on introducing the proposed method and does not provide a comprehensive analysis of NTCGNNs. This method does not require for the construction of Lyapunov–Krasovskii functional (LKF) or model transformation, thus reducing computational effort and complexity. By using this new way and the inequality technique, a passivity criterion for NTCGNNs with multiple time‐varying delays is derived. The applicability of the derived criterion was demonstrated through a numerical example. [ABSTRACT FROM AUTHOR]

In this paper, we find a new regularity criterion of the 3D tropical climate model in terms of one component of the barotropic velocity. Our criteria of the velocity field improve the previous result by Berti et al. Moreover, we extend the regularity criteria of the Lebesgue spaces to that of larger Lorentz spaces. The key of this paper is by using interaction between partial and usual H1$$ {H}^1 $$ estimate to get upper bound of velocities. [ABSTRACT FROM AUTHOR]

This paper introduces and establishes new Hermite‐Hadamard type inequalities specifically tailored for m$$ m $$‐preinvex functions within the (p,q)$$ \left(p,q\right) $$‐calculus framework. These newly developed inequalities come with accompanying left‐right estimates, which enhance their practical utility. The primary objective of this research is to investigate the properties of (p,q)$$ \left(p,q\right) $$‐differentiable m$$ m $$‐preinvex functions and derive inequalities that extend and generalize existing results in the domain of integral inequalities. The techniques employed in this study hold broader implications, finding relevance in various fields where symmetry is paramount. The findings presented in this paper make a significant contribution to the field of analytic inequalities, offering valuable insights into the behavior and characteristics of m$$ m $$‐preinvex functions. Moreover, the established results demonstrate the wider applicability and generalization of analogous findings from prior literature. The techniques and inequalities introduced herein pave the way for further exploration and research in the realm of integral inequalities. [ABSTRACT FROM AUTHOR]

The stability analysis of dynamic continuous structural system (DCSS) has often been investigated by discretizing it into several low‐dimensional elements. The integrated results of all elements are employed to describe the whole dynamic behavior of DCSS. In this paper, DCSS is regarded as the complex dynamic network with the discretized elements as the dynamic nodes and the time‐varying stiffness as the dynamic link relations between them, by which the DCSS can be regarded to be the large‐scale system composed of the node subsystem (NS) and link subsystem (LS). Therefore, the dynamic model of DCSS is proposed as the combination of dynamic equations of NS and LS, in which their state variables are coupled mutually. By using the model, this paper investigates the stability of DCSS. The research results show that the state variables of NS and LS are uniformly ultimately bounded (UUB) associated with the synthesized coupling terms in LS. Finally, the simulation example is utilized to demonstrate the validity of method in this paper. [ABSTRACT FROM AUTHOR]