Dynamics of a predator-prey model with fear effects and gestation delays

Yaping Wang; Yuanfu Shao; Chuanfu Chai; Yaping Wang; Yuanfu Shao; Chuanfu Chai

doi:10.3934/math.2023378

AIMS Mathematics

2023, Volume 8, Issue 3: 7535-7559. doi: 10.3934/math.2023378

Previous Article Next Article

Research article

Dynamics of a predator-prey model with fear effects and gestation delays

College of Science, Guilin University of Technology, Guilin, Guangxi 541004, China

Received: 12 November 2022 Revised: 06 January 2023 Accepted: 10 January 2023 Published: 17 January 2023
MSC : 34C23, 92D25

Recent studies have shown that, in addition to direct predation, fear of predators alters the physiological behavior of prey. Based on this fact, this paper investigates a three-species food chain based on ratio-dependent and Beddington-DeAngelis type functional responses, which incorporates fear effects and two gestation delays. The positivity, boundedness and existence of equilibrium points of the system are investigated, and the local stability behavior of the equilibrium points and the occurrence of Hopf-bifurcation when the time lag parameters exceed the critical values are studied by analyzing the corresponding characteristic equations. The main results show that Hopf-bifurcation occurs when the time delay parameters attain the thresholds. Finally, numerical simulations are performed to verify our main results.

Keywords:

Citation: Yaping Wang, Yuanfu Shao, Chuanfu Chai. Dynamics of a predator-prey model with fear effects and gestation delays[J]. AIMS Mathematics, 2023, 8(3): 7535-7559. doi: 10.3934/math.2023378

Related Papers:

[1]	Taekyun Kim, Dae San Kim, Jin-Woo Park . Degenerate $r$ -truncated Stirling numbers. AIMS Mathematics, 2023, 8(11): 25957-25965. doi: 10.3934/math.20231322
[2]	Taekyun Kim, Dae San Kim, Hyunseok Lee, Lee-Chae Jang . A note on degenerate derangement polynomials and numbers. AIMS Mathematics, 2021, 6(6): 6469-6481. doi: 10.3934/math.2021380
[3]	Dojin Kim, Patcharee Wongsason, Jongkyum Kwon . Type 2 degenerate modified poly-Bernoulli polynomials arising from the degenerate poly-exponential functions. AIMS Mathematics, 2022, 7(6): 9716-9730. doi: 10.3934/math.2022541
[4]	Dae San Kim, Hye Kyung Kim, Taekyun Kim . New approach to $\lambda$ -Stirling numbers. AIMS Mathematics, 2023, 8(12): 28322-28333. doi: 10.3934/math.20231449
[5]	Sang Jo Yun, Jin-Woo Park . On a generation of degenerate Daehee polynomials. AIMS Mathematics, 2025, 10(5): 12286-12298. doi: 10.3934/math.2025556
[6]	Aimin Xu . Some identities connecting Stirling numbers, central factorial numbers and higher-order Bernoulli polynomials. AIMS Mathematics, 2025, 10(2): 3197-3206. doi: 10.3934/math.2025148
[7]	Feng Qi, Da-Wei Niu, Bai-Ni Guo . Simplifying coefficients in differential equations associated with higher order Bernoulli numbers of the second kind. AIMS Mathematics, 2019, 4(2): 170-175. doi: 10.3934/math.2019.2.170
[8]	Taekyun Kim, Dae San Kim, Dmitry V. Dolgy, Hye Kyung Kim, Hyunseok Lee . A new approach to Bell and poly-Bell numbers and polynomials. AIMS Mathematics, 2022, 7(3): 4004-4016. doi: 10.3934/math.2022221
[9]	Hye Kyung Kim . Note on $r$ -central Lah numbers and $r$ -central Lah-Bell numbers. AIMS Mathematics, 2022, 7(2): 2929-2939. doi: 10.3934/math.2022161
[10]	Letelier Castilla, William Ramírez, Clemente Cesarano, Shahid Ahmad Wani, Maria-Fernanda Heredia-Moyano . A new class of generalized Apostol–type Frobenius–Euler polynomials. AIMS Mathematics, 2025, 10(2): 3623-3641. doi: 10.3934/math.2025167

Abstract

1. Introduction

Many people with physical mobility impairments, such as spinal and multiple sclerosis, need to use intelligent wheelchairs for various activities in their surroundings ^[1,2,3]. Conventional wheelchairs may meet the specific needs of some patients, but most patients, however, need to use innovative chairs with advanced technologies that use smart devices and consider various safety measures ^[4]. Recent decades have seen significant advances in developing control technologies for autonomous and semi-autonomous wheelchairs. These developments have been characterized by improving comfort, safety and driving in an ideal way while avoiding obstacles in different conditions, depending on the supervisory control based on the sensory systems and the user's intention. These developments also included using smartphones, developing navigation algorithms, hands-free driving and sharing users' locations to improve system performance conditions ^[3,4].

In general, driving-level control technologies of mobile wheelchairs have been divided into two categories. The supervisory category operates the mobile wheelchairs to perform intelligent tasks, such as avoiding obstacles and following a given path, and are essential to optimum performance ^[5]. The low-level category needs to be developed to raise the level of fine motion control. This is because this level of control deals with many variables that need to design numerical algorithms for the designed systems. Therefore, this kind of research needs to develop efficient algorithms that deal with system uncertainties and external perturbations ^[6]. Nevertheless, there has been a lack of emphasis among researchers on the specific category of low-level control.

Various low-level control techniques have been proposed for driving smart wheelchairs, including proportional-derivative (PD) controller, fuzzy logic control (FLC), sliding mode control (SMC) ^[7], adaptive neural network of nonsingular terminal sliding mode control ^[8], quasi linear parameter varying control ^[9], H-infinity ^[10], Lagrangien method with FLC ^[11], and backstepping control with a proportional-integral-derivative (PID) controller ^[12]. However, the majority of these haven't considered the motorized wheelchair to be a multivariate system. The challenge of setting the control design into action for a multivariate system can be difficult because of the interactions between the various inputs and outputs of the system. The problems inherently present in a multivariate system can frequently be resolved using separation techniques, which are very effective solutions. These methods attempt to alleviate the complexity of the problem of multivariate control by breaking it down into several separate numerical control issues. Nevertheless, there has yet to be a lot of research done looking at how durable systems are when subjected to the influence of uncertainty and external perturbations.

Recently, a model reference control method for a semi-autonomous powered wheelchair has been proposed and evaluated in different condition states. This method adopted quasi-linear closed-loop behavior as a control loop to achieve robustness. To mitigate modeling uncertainties, the technique also compensated for non-linear dynamics associated with a veer rate ^[13]. A digital control and encoding method was also introduced to obtain real-time robustness tracking in smart wheelchairs. The method utilized PID with an open-loop controller to vary the speed in different proposed surfaces. Experimental results demonstrated that this method was effective in proposed surfaces using different speeds ^[14]. Another method conducted by Seki and Tanohata ^[15] implemented an FLC algorithm based on a torque scheme to drive a wheelchair on a disturbance road. In this approach, the algorithm was developed by incorporating human input torque and driving velocity. The effectiveness of this method was evaluated utilizing a real-time experimental test, and the results depicted that this control method was robust on driving on disturbance roads. Despite the evaluations made on these methods, they face several challenges in appropriately addressing the changing parameters of the system and effectively mitigating the impact of noise.

The neural network approach is one of the control tools that have been developed recently to control linear and nonlinear systems. This is done by providing it with many effective characteristics, such as the ability to generalize through experience for many planned and nonlinear inputs. This approach provided the ability to devise different ways to overcome noise and multivariate interactions. Neural networks have demonstrated efficacy in addressing complicated control operations by designing structures for those networks used in solving arithmetic sequences ^[16]. For instance, Aly et al. ^[17] developed a tracking control algorithm for the upper-limb exoskeleton using an adaptive neural backstepping in the presence of external disturbance and model uncertainty. The Lyapunov stability was obtained to achieve the system parameter convergence while evaluating the designed algorithm. The findings from the simulation study indicated that the proposed methodology effectively monitored and mitigated disturbances and uncertainties in the model. The control method for the upper-limb exoskeleton, as presented by the authors in reference [8], involved the utilization of an adaptive neural network based on fixed-time tracking. This approach incorporated an integral nonsingular terminal sliding mode. The method was evaluated using MATLAB simulink tools (The MathWorks Inc., Natick, USA), and the results show that the technique was effective in the presence of external disturbance and model uncertainty. However, the upper proposed methods were only in silico evaluated and updated, so this kind of control law is still complicated.

This paper addresses a number of challenging issues associated with the trajectory tracking control of mobile wheelchairs. This is related to stability and tracking performance, where precise trajectory monitoring is essential for advancing autonomy guidance. Therefore, this work aims to design an efficient method to track the reference trajectory in the presence of disturbances or system parameter variation. The proposed technique employed the neural pole placement sliding mode control (PP-SMC) method to achieve the system robustness. The stability analysis was studied, and the numerical evaluation approach was conducted using MATLAB simulink tools to demonstrate the system's performance. To our knowledge, this is the first study in which the placement method combined with a neural sliding mode that adopted the Gao reaching law has been employed to control the low-level dynamic model of mobile wheelchairs.

2. Materials and methods

2.1. Dynamical model

As shown in , the nonlinear dynamic system model of the mobile wheelchair can be described in the form of cartesian coordinates ( ${\chi _w}＆ {y_w}$ ) and the mobile orientation ( $\Phi$ ) ^[18]. In this model, the center of mass is denoted by ( $H$ ), and the distance between $H$ and the central point of the virtual axis ( $G$ ) is denoted by ( $l$ ). As a result, the updated kinematic model can be illustrated as:

$\left[ {\begin{array}{*{20}{c}} {\dot \chi \left( t \right)} \\ {\dot y\left( t \right)} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {\begin{array}{*{20}{c}} {\cos \Phi \left( t \right)}&{ - l\sin \Phi \left( t \right)} \end{array}} \\ {\begin{array}{*{20}{c}} {\sin \Phi \left( t \right)}&{l\cos \Phi \left( t \right)} \end{array}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{u_1}\left( t \right)} \\ {{u_2}\left( t \right)} \end{array}} \right] + \left[ {\begin{array}{*{20}{c}} {{\gamma _x}\left( t \right)} \\ {{\gamma _y}\left( t \right)} \end{array}} \right],$

(1)

Figure 1. Kinematic controlled model.

DownLoad: Full-Size Img PowerPoint

where $\chi \in {R^n}$ is the state vector of the dynamic system,

$u = \left[ {\begin{array}{*{20}{c}} {{u_1}\left( t \right)} \\ {{u_2}\left( t \right)} \end{array}} \right] \in {R^m}$

is the input vector, $\gamma \left(t \right)$ is the system disturbance,

$y = \left[ {\begin{array}{*{20}{c}} {\upsilon \left( t \right)} \\ {\omega \left( t \right)} \end{array}} \right]$

is the model output of linear velocity (v(t)) and angular velocity (ω(t)).

The dynamic output model is calculated based on the following equations:

$\left\{ {\begin{array}{*{20}{c}} {\upsilon \left( t \right) = {{R\left( {{\omega _r} + {\omega _l}} \right)} \mathord{\left/ {\vphantom {{R\left( {{\omega _r} + {\omega _l}} \right)} 2}} \right. } 2}, } \\ {\omega \left( t \right) = R{{\left( {{\omega _r} - {\omega _l}} \right)} \mathord{\left/ {\vphantom {{\left( {{\omega _r} - {\omega _l}} \right)} {D, }}} \right. } {D, }}} \end{array}} \right.$

(2)

here ${\omega _r}＆ {\omega _l}$ are right and left of angular velocities, D is the distance between the wheels, and R is the right and left wheel radius.

By using the Euler-Lagrangian approach ^[19], the dynamic model in (1) can be expressed as follows:

$\begin{gathered} \dot \chi (t) = \left( {a + \delta a} \right)\chi (t) + bu(t) + \gamma \left( t \right), \hfill \\ \dot y(t) = c\chi (t), \hfill \\ \end{gathered}$

(3)

where, $\delta a$ is the system parameter variation, $\gamma \left(t \right)$ is the system disturbance and a, b and c are matrices with compatible dimensioning.

2.2. Controller design

Consider the dynamic model in (3) is without uncertainty and disturbance is given as:

$\begin{gathered} \dot \chi (t) = a\chi (t) + bu(t), \hfill \\ \dot y(t) = c\chi (t), \hfill \\ \end{gathered}$

(4)

and we assume that the correspondence reference model is given by:

$\dot q(t) = {a_k}q(t) + {b_k}r(t),$

(5)

where $q \in {R^n}$ is the vector that representing the states of the reference model, $r \in {R^i}$ is an input vector for reference model and ( ${a_k}＆ {b_k}$ ) are compatibly dimensioned matrices.

2.2.1. Pole placement method

For the sake of implementing this strategy, we are assuming that the pairs (a & b) and ( ${a_k}＆ {b_k}$ ) are able to be controlled, thus demonstrating stability of the reference model. This implies that the eigenvalues of ${a_k}$ are located within the open unit disc. In order to develop the algorithm, we take into account the tractability of the model referenced in (5) by the model presented in (4). Therefore, this requires to determine the ( ${a_k}＆ {b_k}$ ) based on the following non-proofing theorem ^[20]:

Theorem 1. For any system $a\chi = c$ , where is $a \in {R^{n*m}}$ matrix and $c \in {\mathbb{R}^m}$ is a vector, there exists a solution $rank\left[ a \right] = rank\left[ {{a_c}} \right]$ if and only if ${a_c} \in m\left({n + 1} \right)$ is an augmented matrix.

Based on that we can write: $rank\left[ {\begin{array}{*{20}{c}} b & {{a_k} - a} \end{array}} \right] = rank\left[ b \right]$ ; $rank\left[ {\begin{array}{*{20}{c}} b & {{b_k}} \end{array}} \right] = rank\left[ b \right]$ , then the previous theorem has the significant conclusion that $K$ and $G$ have matrices with compatible dimensions in such a way as:

$\begin{gathered} {a_k} = a - bk, \hfill \\ {b_k} = bg, \hfill \\ \end{gathered}$

(6)

where $k = \left[ {{k_1}, {k_2}, \cdots, {k_{\text{n}}}} \right]$ , w represents the state feedback matrix generated with the pole placement method or Ackerman's formula to achieve asymptotically stable behavior in closed-loop dynamic systems.

The term $g$ represents the feedforward steady-state gain, which may be determined by finding the inverse of the triple components ( $c, \left({a + bk} \right), b$ ) as the following:

Consider that the proposed state feedback law is formulated as follows:

$u\left( t \right) = k\chi \left( t \right) + gr\left( t \right) ,$

(7)

the obtained state equation for the closed-loop system is as follows:

$\begin{gathered} \dot \chi (t) = \left( {a + bk} \right)\chi (t) + bgr(t), \hfill \\ \dot y(t) = c\chi (t). \hfill \\ \end{gathered}$

(8)

To obtain g, we assume that the reference input $r\left(t \right) = \Psi$ and the steady state output can be written as:

${y_{ss}} \equiv \mathop {\lim }\limits_{t \to \infty } y\left( t \right) = \Psi .$

(9)

For the constant $\Psi$ at $t \geqslant 0$ , the closed-loop state equation that involves an equilibrium state can satisfy the following condition:

$\dot \chi (t) = \left( {a + bk} \right){\chi _{ss}}\left( t \right) + bg\Psi = 0,$

(10)

${\chi _{ss}}(t) = c{\left( {a + bk} \right)^{ - 1}}bg\Psi .$

(11)

The steady-state output can be obtained using the following:

${y_{ss}} = c{\chi _{ss}}(t) = c{\left( {a + bk} \right)^{ - 1}}bg\Psi .$

(12)

From (9) we can write:

$\Psi = c{\left( {a + bk} \right)^{ - 1}}bg\Psi .$

(13)

As a result, the feedforward gain can be obtained as:

$g = {\left( {c{{\left( {a + bk} \right)}^{ - 1}}b} \right)^{ - 1}}.$

(14)

In this study, we define the trajectory error as:

$\begin{gathered} {e_\chi }\left( t \right) = \chi \left( t \right) - q\left( t \right), \hfill \\ {{\dot e}_\chi }\left( t \right) = \dot \chi \left( t \right) - \dot q\left( t \right). \hfill \\ \end{gathered}$

(15)

By substituting (4) and (5) in (15), we can write:

${\dot e_\chi }\left( t \right) = a\chi (t) + bu(t) - \left( {{a_k}q(t) + {b_k}r(t)} \right).$

(16)

By adding and subtracting the term ${a_k}\chi \left(t \right)$ we can get:

${\dot e_\chi }\left( t \right) = a\chi (t) + {a_k}\chi \left( t \right) - {a_k}\chi \left( t \right) + bu(t) - {a_k}q(t) - {b_k}r(t),$

(17)

${\dot e_\chi }\left( t \right) = {a_k}{e_\chi }\left( t \right) + \left( {a - {a_k}} \right)\chi \left( t \right) + bu(t) - {b_k}r(t).$

(18)

2.2.2. SMC method

In this study we employed the SMC method and chose the following switch function ^[21,22]:

$\begin{gathered} s\left( t \right) = \xi {e_\chi }\left( t \right), \hfill \\ \dot s\left( t \right) = \xi {{\dot e}_\chi }\left( t \right). \hfill \\ \end{gathered}$

(19)

By substituting (18) in (19) we can get:

$\dot s\left( t \right) = \xi \left( {{a_k}{e_\chi }\left( t \right) + \left( {a - {a_k}} \right)\chi \left( t \right) + bu(t) - {b_k}r(t)} \right) .$

(20)

To achieve a robust stability, we proposed to adopt a Gao reaching law ^[21] as:

$\dot s\left( t \right) = - \left( {\tau s\left( t \right) + \beta sign\left( {s\left( t \right)} \right)} \right),$

(21)

where, $\tau > 0$ , $\beta > 0$ and $sign\left({s\left(t \right)} \right)$ is the signum function given by:

$sign\left( {s\left( t \right)} \right) = \left\{ {\begin{array}{*{20}{c}} {\begin{array}{*{20}{c}} { + 1, }&{s\left( t \right) \geqslant 0, } \end{array}} \\ {\begin{array}{*{20}{c}} { - 1, }&{s\left( t \right) < 0.} \end{array}} \end{array}} \right.$

By equivalating (20) to (21) we can get:

$- \left( {\tau s\left( t \right) + \beta sign\left( {s\left( t \right)} \right)} \right) = \xi \left( {{a_k}{e_\chi }\left( t \right) + \left( {a - {a_k}} \right)\chi \left( t \right) + bu(t) - {b_k}r(t)} \right) .$

(22)

Solving (22) in terms of $u(t)$ we can write:

$u\left( t \right) = - {\left( {\xi b} \right)^{ - 1}}\left( \begin{gathered} \xi {a_k}{e_\chi }\left( t \right) + \xi \left( {a - {a_k}} \right)\chi \left( t \right) - \xi {b_k}r(t) \hfill \\ - \tau \xi {e_\chi }\left( t \right) - \beta sign\left( {\xi {e_\chi }\left( t \right)} \right) \hfill \\ \end{gathered} \right) .$

(23)

The model presented in (4) can be rewritten as follows when there is uncertainty in either the disturbances or the system parameters:

$\begin{gathered} \dot \chi (t) = \left( {a + \delta a} \right)\chi (t) + bu(t) + \gamma \left( t \right), \hfill \\ \dot y(t) = c\chi (t). \hfill \\ \end{gathered}$

(24)

The control algorithm in (24) can be written as:

(25)

The control algorithm ( ${u_t}(t)$ ) in (25) can be constructed as equivalent $\left({{u_{eqv}}(t)} \right)$ and corrective $\left({{u_{corr}}(t)} \right)$ parts as:

$u(t) = {u_{eqv}}(t) + {u_{corr}}(t) ,$

(26)

where,

${u_{eqv}}\left( t \right) = - {\left( {\xi b} \right)^{ - 1}}\left( \begin{gathered} \xi {a_k}{e_\chi }\left( t \right) + \xi \left( {a - {a_k}} \right)\chi \left( t \right) - \xi {b_k}r(t) \hfill \\ - \left( {\xi \delta a\chi \left( t \right) + \xi \gamma \left( t \right)} \right) \hfill \\ \end{gathered} \right),$

(27)

${u_{corr}}\left( t \right) = {\left( {\xi b} \right)^{ - 1}}\left( {\tau \xi {e_\chi }\left( t \right) + \beta sign\left( {\xi {e_\chi }\left( t \right)} \right)} \right) .$

(28)

2.2.3. Neural networks control method

Neural networks method has proposed to adjust the values varying of the system parameters and disturbances by estimating both equivalent and corrective laws ^[23]. In this method, the equivalent neural network keeps the system states on the $s\left(t \right) = 0$ , while the corrective neural network is trained to retune all states back to the proposed sliding surface. To eliminate any chattering, we proposed to replace the $\left({sign\left({s\left(t \right)} \right)} \right)$ function with

$\mu = {{\left( {1 - {e^{ - s\left( t \right)}}} \right)} \mathord{\left/ {\vphantom {{\left( {1 - {e^{ - s\left( t \right)}}} \right)} {\left( {1 + {e^{ - s\left( t \right)}}} \right)}}} \right. } {\left( {1 + {e^{ - s\left( t \right)}}} \right)}} .$

The estimated the equivalent control law is given by:

${\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{u} _{eqv}}(t) = {k_u}\mu \left( {s\left( t \right)} \right)\left( {\sum\limits_{i = 1}^m {{{\rm H}_i}.\mu \left( {s\left( t \right)} \right)\left( {\sum\limits_{j = 1}^n {{{\bar {\rm H}}_{j, i}}.{{\rm O}_i}} } \right)} } \right) .$

(29)

The estimated corrective control law as is given by:

${\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{u} _{corr}}(t) = {k_c}.\mu \left( {s\left( t \right)} \right)\left( {\sum\limits_{j = 1}^n {{\varphi _j}{e_j}\left( t \right)} } \right) .$

(30)

To obtain a robust estimation, the following cost functions were proposed:

$J = 0.5{u^2}_{eqv}(t) + {u_{eqv}}(t)*{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{u} _{eqv}}(t) + 0.5{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\frown}$}}{u} ^2}_{eqv}(t) ,$

(31)

${\rm Z} = 0.5{s^2}(t) .$

(32)

Based on (31), the weight updating law to reduce $J$ can be written as:

$\Delta {{\rm H}_i} = - \beta {{\partial J} \mathord{\left/ {\vphantom {{\partial J} {\partial {{\rm H}_i} = 0.5\left( {\beta {k_u}s\left( t \right) - \beta {k_u}s\left( t \right)\mu \left( {s\left( t \right)} \right){u^2}y\left( t \right)} \right)}}} \right. } {\partial {{\rm H}_i} = 0.5\left( {\beta {k_u}s\left( t \right) - \beta {k_u}s\left( t \right)\mu \left( {s\left( t \right)} \right){u^2}y\left( t \right)} \right)}} ,$

(33)

$\begin{array}{l}\Delta {\stackrel{⌢}{{\rm H}}}_{j, i} = {}^{-\beta \partial J}\!\!\diagup\!\!{}_{\partial {\stackrel{⌢}{{\rm H}}}_{j, i}}\;\\ = 0.25\left(\beta s(t){{\rm H}}_{i}-\beta s\left(t\right)\mu \left(s\left(t\right)\right){u}^{2}\left(t\right){{\rm H}}_{i}\right)\left(1-\mu \left(s\left(t\right)\right){y}^{2}\left(t\right)\right){F}_{i}\text{,}\end{array}$

(34)

where $\beta > 0$ , also based on (32), the weight updating law to reduce Z can be written as:

$\Delta {\xi _j} = - \alpha {\raise0.7ex\hbox{${\partial J}$} \!\mathord{\left/ {\vphantom {{\partial J} {\partial {\xi _j}}}}\right.} \!\lower0.7ex\hbox{${\partial {\xi _j}}$}} = - \alpha s(t){e_j}\left( t \right) ,$

(35)

where $\alpha > 0$ , k in (14) can be calculated based on the updated vector $\xi$ as:

$k = \lambda {\left( {{\xi ^T}b} \right)^{ - 1}} ,$

(36)

where, $\lambda > 0$ . In both equivalent and corrective control laws, a three-layer feedforward neural network is employed. depicts the block diagram of the proposed control system where the overall control algorithm ( ${u_t}\left(t \right)$ ) represents the summation of the neural output structure with the function $g$ in (14).

Figure 2. Block diagram of the proposed control system.

DownLoad: Full-Size Img PowerPoint

3. Simulation protocols and setting of the controller parameters

To evaluate the effectiveness of the proposed control method, we proposed to utilize various scenarios through the application of the MATLAB simulink software tools. In the first scenario, the control algorithm was examined under nominal conditions to track the desired trajectory. Here, the wheelchair motion started from the initial point. For the second scenario, we suggested to manipulate the system parameters while monitoring the dynamic system in order to follow the intended trajectory. Subsequently, uncertainties and external disturbances were introduced to the dynamic system and monitoring the performance of the controller to trace the desired trajectory. In all scenarios, we proposed the following reference signal to be tracked:

$r\left( t \right) = 2 + 2\sin \left( {2\pi + \Phi } \right).$

(37)

In addition, mean absolute error (MAE) was also used to measure the error difference between actual and reference trajectory as:

$MAE = {{\sum\nolimits_{i = 1}^n {\left| {{y_i} - {x_i}} \right|} } \mathord{\left/ {\vphantom {{\sum\nolimits_{i = 1}^n {\left| {{y_i} - {x_i}} \right|} } n}} \right. } n},$

(38)

where, ${y_i}$ is the actual trajectory and ${x_i}$ is reference trajectory. Table 1 shows the gain values and online training gains that were assigned to the parameters for the proposed control algorithm.

Table 1. The gain values of control algorithm.

Parameter	Value
$g$	-34.35
$k$	[0.067, -0.135]
$\xi$	[0.152, -2.5]
${k}_{u}$	3.4
${k}_{c}$	0.67
$\tau \beta$	0.14
$\epsilon \beta$	0.77
Training gains
${H}_{\mathrm{1, 2}}$	0.22
${H}_{\mathrm{1, 1}}$	0.46

| Show Table

DownLoad: CSV

4. Simulation results

The responses of system stability for the dynamic model was obtained using the unit step input signal with the initial condition response $\left[ {\begin{array}{*{20}{c}} 1 & 1 \end{array}} \right]$ . Figure 3 is representing the unforced states of the dynamic system model.

Figure 3. States responses of dynamic model.

DownLoad: Full-Size Img PowerPoint

4.1. The performance observed under nominal conditions

The results of the simulation of the movable wheelchair on the x-axis and y-axis under normal conditions are depicted in and . According to the figures, the proposed system successfully achieved high tracking performance while maintaining a low error rate between the reference trajectory and the actual trajectory. The findings also showed that the control algorithm was able to track the reference trajectory with the value of $MAE = 0.08$ , as shown in Figure 5. The proposed technique also demonstrated that the system stability and robustness were achieved quickly. This is evidence of the simplicity of the proposed control algorithm. This simplicity makes the proposed control algorithm suitable for tiny systems with minimal processing demands.

Figure 4. Tracking performance under nominal conditions.

DownLoad: Full-Size Img PowerPoint

Figure 5. MAE between actual and reference under nominal conditions.

DownLoad: Full-Size Img PowerPoint

4.2. The performance observed under system parameters variations

In this section, system performance is evaluated by varying the parameters of the dynamic system during the wheelchair motion. depicts the position tracking along the x-axis and y-axis. The simulation results demonstrate that the controller is capable of managing the parameters varying within an appropriate period of time. In addition, the distance error that occurs for the controllers when comparing the reference trajectory to the actual trajectory is depicted in . By applying the proposed strategy, we were able to achieve MAE values that were lower, resulting to $MAE = 0.10$ .

Figure 6. Tracking performance observed under system parameters variations.

DownLoad: Full-Size Img PowerPoint

Figure 7. MAE between actual and reference trajectory under system parameters variations.

DownLoad: Full-Size Img PowerPoint

4.3. The performance observed under uncertainties and external disturbances

In part, the performance of the robustness control system is assessed by introducing uncertainties and external disturbances into the dynamic system. The position tracking in the x-axis and y-axis under the influence of this condition is depicted in . This figure shows that the proposed control system was able to track the reference trajectory in presence of uncertainties and external disturbances. Similarly, the MAE criterion was utilized to compare between the actual and the reference trajectory. depicts the disparity in distance between the reference trajectory and the actual trajectory for the controllers, with the recorded value $MAE = 0.12.$

Figure 8. The performance observed under external disturbances.

DownLoad: Full-Size Img PowerPoint

Figure 9. MAE between actual and reference trajectory under external disturbances.

DownLoad: Full-Size Img PowerPoint

5. Discussion

Numerous techniques for optimal mobile-powered wheelchair use have been developed over the past 20 years ^[4,24,25]. As a result, there is a growing demand for developing dynamic models and control systems that meet the needs for enhancing self-efficacy and safety. Enhancing the performance of these models can be achieved by improving robust control algorithms ^[26,27]. Recently, nonlinear control approaches and computational Artificial Intelligence (AI) tools such as neural networks constitute a potentially fruitful strategy for managing electromechanical actuators. However, the properties of dynamical systems can change, and there may also be external perturbations, which makes it difficult to create efficient algorithms ^[28]. In this work, a novel control method is developed and evaluated to track the nonlinear motion of wheelchairs precisely. A sliding neural network approach is integrated with the method of PP-SMC to eliminate the uncertainties of dynamic system parameters and disturbances. The findings conducted through simulink demonstrated that the suggested system could precisely track the desired trajectory, as shown in Figures 4–9.

The primary challenge encountered in controlling these dynamic systems is the mitigation of disturbances and addressing the parameter variation. Consequently, numerous control systems are engineered to mitigate these issues. Our study employed a new technique that utilized the neural network approach to reduce the effects of disturbance and parameter variation. The design methodology introduced a model reference to construct a feedforward control part and then adopted the neural PP-SMC. The results concluded that the performance of the proposed approach was efficient and robust in tracking reference trajectory with minimal error ranging between 0.08 and 0.12. In addition, the proposed method is characterized by the results illustrated by the ability to reduce the differences between the referenced and actual trajectory. Also, the performance of the neural PP-SMC controller can be enhanced by optimizing the parameter determination process compared to other controllers.

In the context of multi-input multi-output (MIMO) control systems, certain researchers propose using decoupling techniques in servo systems to attain robustness in control systems. Nguyen et al. ^[29] conducted an empirical investigation to validate the efficacy of a real-time multivariate control strategy employing the decoupling technique. This strategy employs the sliding neural network approach as the control method to mitigate uncertainties and external disturbances. Furthermore, the examination of system stability and attainment of global convergence was also accomplished. Nevertheless, there is still chattering appearing on the system tracks for their method. Therefore, the authors are required to improve the system's performance by replacing the shifted sigmoid function with the saturation function. Compared to this work, the utilization of a shifted sigmoid function in our study has effectively eliminated the occurrence of chattering, as depicted in Figure 2.

In recent work, Zhang et al. ^[30] proposed a robust H-infinity FLC method for powered wheelchairs. In this method, an iterative linear matrix inequality algorithm was utilized via an acceleration feedback signal to eliminate disturbances and approximate errors. The results demonstrated that this method was effective to control the dynamic system. A different study conducted by Al-Mahturi et al. ^[31] developed a model-free method that utilized a type-2 FLC system for mobile robots under large uncertainties. The method adopted the structure and parameters learning to achieve the desired set point. The Lyapunov stability theory is employed to carry out the proposed control strategy. To prove the effectiveness of the proposed method, the method needed to manage the motion of a mobile robot while under the influence of various disturbances. The results illustrated that this method was efficient to eliminate the external disturbance. Nevertheless, the implementation of FLC in many applications is accompanied by several drawbacks, including the need for meticulous tuning settings and its susceptibility to noise interference.

As shown in Table 2, this paper provides a concise overview of the methodologies employed, control methods and evaluation type of tracking control algorithms in various studies. Moreover, the findings of this study are susceptible to various limitations, including an actual evaluation. The assessment of this study necessitates the incorporation of empirical experimentation alongside a comparative analysis with simulation-based research. Furthermore, it should be noted that the software model lacked the incorporation of reflex control and automatic regulatory systems.

Table 2. Technique, control method and evaluation type for different studies.

Authors	Technique	Control method	Evaluation type
Aly et al. ^[8]	Dynamic model based	Adaptive neural network	Simulation
Wang et al. ^[14]	Dynamic model based	Proportional-integral-derivative (PID)	Simulation
Martins et al. ^[18]	Dynamic model based	Adaptive controller	Simulation and experimental
Nguyen et al. ^[29]	Multivariable decoupling	Neuro-Sliding Mode	Real-time experiment
De La Cruz et al. ^[32]	Dynamic model based	Adaptive controller	Real-time experiment
Al-Mahturi et al. ^[31]	Dynamic model based	Type-2 evolving fuzzy control	Simulation

| Show Table

DownLoad: CSV

The potential influence of the reflexive and automated systems' ability to adapt to rapid changes in wheelchair movement will likely affect the current control strategy significantly. In addition, work has not been designed for underactuated mobile wheelchair systems. However, the developed technique of underactuated wheeled mobile vehicles has offered a wide range of dynamic positioning to track a path trajectory ^[33,34,35]. Therefore, this work should be designed and implemented in future work.

6. Conclusions

In this study, a tracking control method for operating the dynamic model of the powered wheelchair was developed and implemented. The technique used a reference model, pole placement method, and sliding mode control approach. A neural network-based technology was utilized to eliminate uncertainty and system disturbances. In this method, we proposed to use a sinewave signal as a modeled reference signal to be tracked by the designed control system. The proposed tracking control method was evaluated in MATLAB simulink software tools. The gain values and online training gains that were assigned to the parameters for the proposed control algorithm were obtained to achieve robust tracking. The method was examined under three scenarios: operating under normal conditions, studying the dynamic system's performance in changing parameter settings and evaluating the tracking system's performance in the presence of external and unknown disturbances. The results showed that the developed control algorithm efficiently tracks the reference trajectory. In all scenarios, the MAE was between 0.08 and 0.12 in the presence of uncertainties or perturbations.

Use of AI tools declaration

The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

Acknowledgments

The authors extend their appreciation to the King Salman Center for Disability Research for funding this work through Research Group no KSRG-2023-166.

Conflict of interest

The authors declare no conflicts of interest.

Appendix

The nomenclature for the dynamic model of a mobile wheelchair and the developed control algorithm is displayed in the following table.

Table A1. The parameters' meanings of the dynamic model of the mobile wheelchair and for the designed control algorithm.

Notations	Definitions
Dynamic model of mobile wheelchair
${x_w}＆ {y_w}$	cartesian coordinates
Φ	mobile orientation
H	center of mass
G	central point of the virtual axis
l	distance between $H$ and $G$
x $\in$ Rⁿ	state vectors of the dynamic system
u $\in$ R^m	input vector
γ(t)	system disturbance
y(t)	model output
v(t)	linear velocity
ω(t)	angular velocity
R	right and left wheel radius
D	distance between the wheels
a, b and c	compatibly dimensioned matrices
δa	system parameter variation
Controller
q $\in$ Rⁿ	reference model state vector
r $\in$ R^j	reference model input vector
a_k, b_k	compatibly dimensioned matrices for reference model
k=[k₁, k₂, $\cdots$ , k_n]	state feedback matrix
g	feed-forward steady-state gain
e_x(t)	trajectory error
s(t)	switch function
u_eqv(t)	equivalent control part
u_corr(t)	corrective control part
û_eqv(t)	estimated equivalent control part
û_corr(t)	estimated corrective control part
J & Z	weighting functions

| Show Table

DownLoad: CSV

References

[1]	A. Das, G. P. Samanta, A prey-predator model with refuge for prey and additional food for predator in a fluctuating environment, Physica A, 538 (2019), 427–450. https://doi.org/10.1016/j.physa.2019.122844 doi: 10.1016/j.physa.2019.122844
[2]	D. Pal, G. S. Mahaptra, G. P. Samanta, Optimal harvesting of prey-predator system with interval biological parameters: A bioeconomic model, Math. Biosci., 241 (2013), 181–187. https://doi.org/10.1016/j.mbs.2012.11.007 doi: 10.1016/j.mbs.2012.11.007
[3]	S. Sharma, G. P. Samanta, A Leslie-Gower predator-prey model with disease in prey incorporating a prey refuge, Chaos Soliton. Fract., 70 (2015), 69–84. https://doi.org/10.1016/j.chaos.2014.11.010 doi: 10.1016/j.chaos.2014.11.010
[4]	G. Tang, S. Tang, R. A. Cheke, Global analysis of a holling type Ⅱ predator-prey model with a constant prey refuge, Nonlinear Dynam., 76 (2014), 635–647. https://doi.org/10.1007/s11071-013-1157-4 doi: 10.1007/s11071-013-1157-4
[5]	J. P. Tripathi, S. Abbas, M. Thakur, A density dependent delayed predator-prey model with Beddington-DeAngelis type function response incorporating a prey refuge, Commun. Nonlinear Sci. Numer. Simul., 22 (2015), 427–450. https://doi.org/10.1016/j.cnsns.2014.08.018 doi: 10.1016/j.cnsns.2014.08.018
[6]	A. Q. Khan, S. A. H. Bukhari, M. B. Almatrafi, Global dynamics, Neimark-Sacker bifurcation and hybrid control in a Leslie's prey-predator model, Alex. Eng. J., 61 (2022), 11391–11404. https://doi.org/10.1016/j.aej.2022.04.042 doi: 10.1016/j.aej.2022.04.042
[7]	A. Q. Khan, F. Nazir, M. B. Almatrafi, Bifurcation analysis of a discrete Phytoplankton–Zooplankton model with linear predational response function and toxic substance distribution, Int. J. Biomath., 16 (2022). https://doi.org/10.1142/S1793524522500954 doi: 10.1142/S1793524522500954
[8]	A. Q. Khan, M. Tasneem, M. B. Almatrafi, Discrete-time COVID-19 epidemic model with bifurcation and control, Alex. Eng. J., 19 (2022), 1944–1969. https://doi.org/10.3934/mbe.2022092 doi: 10.3934/mbe.2022092
[9]	N. Mondal, H. Alrabaiah, D. Barman, J. Roy, S. Alam, Influence of predator incited fear and interference competition in the dynamics of prey-predator system where the prey species are protected in a reserved area, Ecol. Environ. Conserv., 28 (2022), 831–852. http://doi.org/10.53550/EEC.2022.v28i02.040 doi: 10.53550/EEC.2022.v28i02.040
[10]	D. Mukherjee, Role of fear in predator-prey system with intraspecific competition, Math. Comput. Simul., 177 (2020), 263–275. https://doi.org/10.1016/j.matcom.2020.04.025 doi: 10.1016/j.matcom.2020.04.025
[11]	Y. Shao, Global stability of a delayed predator-prey system with fear and Holling-type Ⅱ functional response in deterministic and stochastic environments, Math. Comput. Simul., 200 (2022), 65–77. https://doi.org/10.1016/j.matcom.2022.04.013 doi: 10.1016/j.matcom.2022.04.013
[12]	L. Y. Zanette1, A. F. White, M. C. Allen, M. Clinchy, Perceived predation risk reduces the number of offspring songbirds produce per year, Science, 334 (2011), 1398–1401. https://doi.org/10.1126/science.1210908 doi: 10.1126/science.1210908
[13]	S. Creel, D. Christianson, S. Liley, J. A. Winnie, Predation risk affects reproductive physiology and demography of elk, Science, 315 (2007), 960. https://doi.org/10.1126/science.1135918 doi: 10.1126/science.1135918
[14]	J. P. Suraci, M. Clinchy, L. M. Dill, D. Roberts, L. Y. Zanette, Fear of large carnivores causes a trophic 585 cascade, Nat. Commun., 7 (2016), 10698. https://doi.org/10.1038/ncomms10698 doi: 10.1038/ncomms10698
[15]	X. Wang, L. Zanette, X. Zou, Modelling the fear effect in predator-prey interactions, J. Math. Biol., 73 (2016), 1179–1204. https://doi.org/10.1007/s00285-016-0989-1 doi: 10.1007/s00285-016-0989-1
[16]	C. S. Holling, The components of predation as revealed by a study of small-mammal predation of the european pine sawfly, Canadian Entomol., 9 (1959), 293–320. https://doi.org/10.4039/Ent91293-5 doi: 10.4039/Ent91293-5
[17]	C. S. Holling, The functional response of predators to prey density and its role in mimicry and population regulation, Mem. Entomol. Soc. Can., 97 (1965), 5–60. https://doi.org/10.4039/entm9745fv doi: 10.4039/entm9745fv
[18]	C. S. Holling, Some characteristics of simple types of predation and parasitism, Canadian Entomol., 91 (1959), 385. https://doi.org/10.4039/Ent91385-7 doi: 10.4039/Ent91385-7
[19]	Z. Liu, R. Yuan, Stability and bifurcation in a delayed predator-prey system with Beddington-DeAngelis functional response, J. Math. Anal. Appl., 296 (2004), 521–537. https://doi.org/10.1016/j.jmaa.2004.04.051 doi: 10.1016/j.jmaa.2004.04.051
[20]	X. Shi, X. Zhou, X. Song, Analysis of a stage-structured predator-prey model with Crowley-Martin function, J. Appl. Math. Comput., 36 (2011), 459–472. https://doi.org/10.1007/s12190-010-0413-8 doi: 10.1007/s12190-010-0413-8
[21]	J. P. Tripathi, S. Tyagi, S. Abbas, Global analysis of a delayed density dependent predator-prey model with Crowley-Martin functional response, Commun. Nonlinear Sci. Numer. Simul., 30 (2016), 45–69. https://doi.org/10.1016/j.cnsns.2015.06.008 doi: 10.1016/j.cnsns.2015.06.008
[22]	R. Chinnathambi, F. A. Rihan, Stability of fractional-order prey-predator system with time-delay and monod-haldane functional response, Nonlinear Dynam., 92 (2018), 1637–1648. https://doi.org/10.1007/s11071-018-4151-z doi: 10.1007/s11071-018-4151-z
[23]	J. Roy, D. Barman, S. Alam, Role of fear in a predator-prey system with ratio-dependent functional response in deterministic and stochastic environment, Biosystems, 197 (2020), 104176. https://doi.org/10.1016/j.biosystems.2020.104176 doi: 10.1016/j.biosystems.2020.104176
[24]	M. Banerjee, Y. Takeuchi, Maturation delay for the predators can enhance stable coexistence for a class of prey-predator models, J. Theor. Biol., 412 (2017), 154–171. https://doi.org/10.1016/j.jtbl.2016.10.016 doi: 10.1016/j.jtbl.2016.10.016
[25]	B. Dubey, A. Kumar, A. P. Maiti, Global stability and Hopf-bifurcation of prey-predator system with two discrete delays including habitat complexity and prey refuge, Commun. Nonlinear Sci. Numer. Simul., 67 (2019), 528–554. https://doi.org/10.1016/j.cnsns.2018.07.019 doi: 10.1016/j.cnsns.2018.07.019
[26]	J. Zhang, Bifurcation analysis of a modified Holling-Tanner predator-prey model with time delay, Appl. Math. Model., 36 (2012), 1219–1231. https://doi.org/10.1016/j.apm.2011.07.071 doi: 10.1016/j.apm.2011.07.071
[27]	D. Hu, Y. Li, M. Liu, Y. Bai, Stability and Hopf bifurcation for a delayed predator-prey model with stage structure for prey and ivlev-type functional response, Nonlinear Dynam., 99 (2020), 3323–3350. https://doi.org/10.1007/s11071-020-05467-z doi: 10.1007/s11071-020-05467-z
[28]	Y. Bai, Y. Li, Stability and Hopf bifurcation for a stage-structured predator-prey model incorporating refuge for prey and additional food for predator, Adv. Differ. Equ., 1 (2019), 1–20. https://doi.org/10.1186/s13662-019-1979-6 doi: 10.1186/s13662-019-1979-6
[29]	F. Chen, M. You, Permanence for an integrodifferential model of mutualism, Appl. Math. Comput., 186 (2007), 30–34. https://doi.org/10.1016/j.amc.2006.07.085 doi: 10.1016/j.amc.2006.07.085
[30]	E. Beretta, Y. Kuang, Geometric stability switch criteria in delay differential systems with delay-dependent parameters, SIAM J. Math. Anal., 33 (2002), 1144–1165. https://doi.org/10.1137/S0036141000376086 doi: 10.1137/S0036141000376086
[31]	L. Deng, X. Wang, M. Peng, Hopf bifurcation analysis for a ratio-dependent predator-prey system with two delays and stage structure for the predator, Appl. Math. Comput., 231 (2014), 214–230. https://doi.org/10.1016/j.amc.2014.01.025 doi: 10.1016/j.amc.2014.01.025
[32]	X. Yang, L. Chen, J. Chen, Permanence and positive periodic solution for the single-species nonautonomous delay diffusive models, Comput. Math. Appl., 32 (1996), 109–116. https://doi.org/10.1016/0898-1221(96)00129-0 doi: 10.1016/0898-1221(96)00129-0
[33]	S. L. Ross, Differential equations, New York: John Wiley and Sons, 1984.

This article has been cited by:

Saim Ahmed, Ahmad Taher Azar, Ibraheem Kasim Ibraheem, Model-free scheme using time delay estimation with fixed-time FSMC for the nonlinear robot dynamics, 2024, 9, 2473-6988, 9989, 10.3934/math.2024489

Reader Comments

Your name:*

Email:*
© 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)