Dynamic analysis of a mosquito population model with a stage structure and periodic releases of sterile males

1.   Introduction

Owing to the distinctive characteristics of high power density and high efficiency, permanent magnet synchronous motor (PMSM) has drawn extensive popularity over the last several decades ^[1,2,3]. For this reason, PMSM has been widely applied for industry, especially for electric vehicles, vacuum cleaners, electric ship propulsion and more. Nevertheless, the traditional proportional integral (PI) control algorithm is typically difficult to offer the high-precision control property ^[4,5,6,7]. To improve the unsatisfactory performance of the PI controller, lots of advanced algorithms have been applied to further optimize the control precision, such as predictive control ^[8,9], fuzzy control ^{[10,11,12,13]}, adaptive control ^{[14,15,16,17,18]} and sliding mode control (SMC) ^{[19,20,21,22,23]}.

Among the previously mentioned modern control approaches, the SMC is an effective strategy to enhance the anti-disturbance capability of the PMSM system ^[24]. A traditional linear SMC method was presented in ^[25] to obtain a fast dynamic response and strong robustness. However, the parameter variation has not been completely considered in theory, which limits the stable operation of motors. In ^[26], a new piecewise but differentiable switching controller was introduced with a proper design idea to avert the singularity problem. Moreover, in order to make the speed error converge faster, the nonsingular terminal sliding mode (NTSM) control algorithm was proposed ^[27] to make sure that the speed error will finite-time converge to the origin.

The NTSM strategy can handle the singular issue and provide satisfactory tracking accuracy. However, the NTSM controller still brings severe chattering^[28,29], and a fast NTSM controller with no switching function was used to weaken the chattering ^[30]. However, this kind of controller has poor anti-disturbance performance and only converges the speed error to a region instead of to zero. At the same time, a brand-new adaptive terminal sliding mode (TSM) reaching law combined with a fast TSM control method was utilized to construct the speed controller for the PMSM system in ^[31]. Even if an adaptive law was adopted to reduce chattering, this method still cannot adjust the control gain automatically with the change of total disturbance.

Notably, utilizing the disturbance observer (DOB) to estimate the unknown time-varying disturbance is a valid strategy to avoid the overestimated switching gain in SMC ^[32,33,34]. Compensating the precise estimates to the baseline controller, the PMSM drive system will achieve significant steady-state performance and strong disturbance rejection property synchronously ^[35]. Without the accurate information of the PMSM mathematical model, a conventional extended state observer (ESO) was applied ^[36] to restrain the property deterioration of PMSM in the existence of external disturbance and parameter variation. Unfortunately, the baseline controller only adopted the linear PI method. It is extremely hard for this category of composite controllers to achieve strong disturbance rejection property. On this basis, the traditional ESO combined with a simple SMC algorithm was constructed ^[37] to address the above tough issue. Although the traditional ESO can estimate the total disturbance exactly, it just guarantees that the estimated error converges asymptotically to zero, which may result in poor estimation precision. In order to obtain higher estimation accuracy and robustness, various improved ESOs have been proposed in recent years^[38,39]. In ^[40], the natural evolution theory has been applied to graph structure learning, where an evolutionary method was constructed to evolve a population of graph neural network (GNN) models to adapt to dynamical environments. In ^[41], a heterogeneous network representation learning method was reported to characterize implicitly inside Ethereum transactions. In ^[42], a center-based transfer feature learning with classifier adaptation for the surface defect recognition was proposed. In ^[43], the authors proposed an integrated triboelectric nanogenerator and tribovoltaic nanogenerator in the air cylinder as difunctional pneumatic sensors for simultaneous position and velocity monitoring. In ^[44], the recent developments in the area of arc fault detection were studied.

In this paper, a novel adaptive nonsingular terminal sliding mode (ANTSM) controller combined with a modified ESO (MESO) is proposed to enhance the disturbance rejection property of the PMSM system. First of all, considering the parameter variation and time-varying load torque, an ANTSM controller is constructed to provide the desired steady-state and dynamic performance. Furthermore, in order to solve the trouble of unsatisfactory control effect due to large switching gain, one novel MESO is utilized to observe the unknown time-varying disturbance and compensate the estimates to the ANTSM controller simultaneously. Lastly, the validity of the ANTSM + ESO composite control algorithm is proved by comprehensive experiments. The main contribution of this article could be summarized in the following points:

$1)$ A novel NTSM controller is designed. The control gain is automatically tuned by the proposed adaptive law. It avoids unsatisfactory control performance caused by excessive control gain.

$2)$ A MESO is proposed to improve the anti-disturbance capability of the PMSM system. By using the finite-time technique, the estimation error can converge to zero in finite time. The proposed MESO has a higher estimation accuracy and faster estimation speed.

$3)$ The proposed composite controller combines the ANTSM algorithm with MESO to improve the system robustness. Since the high disturbance estimation accuracy of MESO, the rejection ability to disturbances of the PMSM speed regulation system can be improved.

The rest of the paper is organized as follows. In Section 2, the PMSM mathematical model with parameter perturbation and the traditional NTSM control method are described. Section 3 shows the design of the proposed MESO-based ANTSM controller in detail. Comprehensive experiments are illustrated in Section 4. The summary of this paper is presented in Section 5.

Notations: Throughout the paper, the symbol ${\left\lfloor {x}\right\rceil}^{{m}}$ is used to present the $|x|^{m} \cdot \text{sign}(x)$ for a real number $m$, and the symbol ${\left\lfloor {x}\right\rceil}^{{*}}$ is defined as

2.   Problem formulation

2.1.   Mathematical model of PMSM

In the ideal case, the motion equation of PMSM can be given by

where $\Omega$ is the rotor angular speed, $n_p$ is the number of pole pairs, ${\psi _f }$ is the flux linkage, $J$ is the moment of inertia, $B$ is the viscous friction coefficient, $i_{q}$ is the stator current in the $q$-axis and $T_{L}$ is the load torque.

Because of load variation in practical applications, the value of inertia may be mismatched. Therefore, we introduce $\Delta J = J - J_0$, where $J_0$ is the nominal value and $\Delta J$ is parameter variation. Then, replacing current $i_q$ by the reference current $i_{q}^{*}$, system (1) can be rewritten as

where $d_0\left(t \right) = -\frac{T_L}{J_0} + b\left(i_q-i_{q}^{*} \right) -\frac{\Delta J}{J_0}\dot{\Omega}$ and $b = \frac{1.5n_p\psi _{f}}{J_0}$.

2.2.   Traditional NTSM controller design

The traditional NTSM algorithm is used to design the speed controller for the PMSM system in this subsection. First of all, setting $\Omega_r$ as the reference speed, the speed error $\Omega_e$ is expressed as

With the help of (2.2), one has

where $d(t) = -d_0\left(t \right)+\dot{\Omega}_r+\frac{B}{J_0}\Omega _r$ is the total disturbance.

Assumption 1: ^[45,46,47] The total disturbance $d(t)$ is bounded and differentiable, and there exist known positive constants $l_1$ and $l_2$ such that $|d(t)| \le {l_1}$ and $|\dot d(t)| \le {l_2}$.

Remark 1: From (2.2) and (2.4), we can get that the disturbance $d(t)$ consists of the load torque $T_L$, the stator current $i_q$, the stator reference current $i_q^*$, the speed $\Omega$, the setting speed $\Omega_r$ and other components. It can be concluded that the above variables are bounded and differentiable. Therefore, it is reasonable that the total disturbance satisfies $|d(t)| \le {l_1}$ and $|\dot d(t)| \le {l_2}$, where $l_1$ and $l_2$ are known positive constants.

According to ^[48], the nonsingular terminal sliding manifold is selected as

where $\beta$ is positive constant, $p$ and $q$ are positive odd integers and $1 < p/q < 2$.

Based on this, the NTSM controller will be designed as

with a positive constant $k$.

Choose the widely-used Lyapunov function as

Differentiating $V(s)$ gets

Substituting NTSM controller (2.6) into (2.8), one obtains

Similar to ^[27], if the control gain satisfies $k > l_1$, then the speed error will converge to origin in a finite time.

Remark 2: To guarantee the stability of the PMSM, the value of $k$ in the traditional NTSM controller should be chosen to be larger than the upper bound of the total disturbance. However, the total disturbance in practical engineering is time-varying, and its upper bound may be much smaller than the fixed control gain. Since the control gain directly affects the chattering amplitude, conservative control gain may lead to unsatisfactory control performance in PMSM. Thus, it is urgent to design a novel controller with variable and a small enough control gain to suppress the time-varying disturbance.

3.   Proposed controller design

To address the forgoing disadvantages of the conventional NTSM controller, a novel ANTSM is developed in this section to automatically tune the value of control gains. Based upon this, compensating the disturbance estimation value derived by MESO to the ANTSM controller, the control gain can get a further reduction.

3.1.   Construction of ANTSM controller

With the aid of controller (2.6), the ANTSM controller is constructed as

The adaptive law $k(t)$ is given as

where $k_M > l_1$ and $k_m > 0$ are the maximums and minimums of the control gain, constant $\varepsilon \in (0, 1)$, $\eta > \frac{l_2}{\varepsilon k_m}$, $N > \eta k_M$, and $[{\rm sign}(s)]_{av}$ is produced by the signal $z(t)$ of the low-pass filter

with $\lambda$ being the tunable parameter. The ANTSM controller structure is depicted by Figure 1.

Remark 3: Since the discontinuous switching function ${\rm sign}(s)$ varies at negative one and one, the function $[{\rm sign}(s)]_{av}$ is continuous and its value belongs to $(-1, 1)$. Therefore, it can be assumed that the second derivative of function $[{\rm sign}(s)]_{av}$ satisfies $|\frac{d^2}{dt^2}[{\rm sign}(s)]_{av}| \leq C$ and $C > 0$.

Next, we will prove that $k(t)$ can converge in a finite time to the minimum absolute value of the total uncertainty. In addition, the equivalent control theory is crucial for the proof in the paper.

The equivalent controller is constructed as

Since the exact information about the total disturbance $d(t)$ is not available, the equivalent controller (3.4) cannot be directly applied to the PMSM system. Therefore, an average control of ANTSM controller (3.1) is employed to track controller (3.4). Meanwhile, to enhance the accuracy of the average control, a continuous function $[{\rm sign}(s)]_{av}$ is used to replace the discontinuous function ${\rm sign}(s)$ in equivalent controller. Then, the average control of ANTSM controller (3.1) is constructed as

According to the above analysis, if the average control (3.5) can track the equivalent controller (3.4) well, one has

If $[{\rm sign}(s)]_{av}$ approaches one, the $k(t)$ is sufficiently large to counteract the $d(t)$. The basic idea of the ANTSM control approach is to ensure that $k(t)$ can converge to $\frac{|d(t)|}{|[{\rm sign}(s)]_{av}|}$ with $\varepsilon$ being close to one in a finite time. Meanwhile, the continuous function $[{\rm sign}(s)]_{av}$ should converge to parameter $\varepsilon$.

Consider the following Lyapunov function:

Evaluating the derivative of $V_1(\delta)$ yields

Since the range of control gain $k\left(t \right) \in \left[k_m, k_M \right]$, one has $\frac{\left| d\left(t \right) \right|}{\varepsilon} > k_m$. With the help of Assumption 1, (3.8) can be rewritten as

Combining $\eta > \frac{l_2}{\varepsilon k_m}$ and (3.9), one has

where $\kappa = \frac{\sqrt2 (\varepsilon k_m\eta-l_2)}{k_M} > 0$.

Therefore, inequity (3.10) satisfies the finite-time stability theorem ^[49], and it indicates that $k(t)$ will converge to $\frac{\left| d\left(t \right) \right|}{\varepsilon}$.

3.2.   Construction of MESO

In this subsection, the conventional ESO and the presented MESO are exploited to estimate the disturbance $d_0(t)$, respectively. First, provided that $|{\dot d_0}| \le {l_0}, l_0 > 0$, the traditional ESO designed for system (2.2) can be expressed as

where $\tilde{\Omega} = \hat{\Omega}-\Omega$, $\hat{\Omega}$ and $\hat{d}_0$ are the estimated values of $\Omega$ and $d_0$, $h_1$ and $h_2$ are the gains of traditional ESO.

Nevertheless, the traditional ESO just ensures that the estimation of error converges to zero asymptotically, which causes poor estimation precision. Therefore, it is urgent to accelerate the estimation speed of traditional ESO. For (2.2), the MESO is constructed as

where functions $\phi _1\left(\tilde{\Omega} \right)$ and $\phi _2\left(\tilde{\Omega} \right)$ are given by

Subtracting (3.12) from (2.2) obtains

where $\tilde{d}_0 = \hat{d}_0-d_0$.

By defining $\xi ^T = \left[\phi _1\left(\tilde{\Omega} \right), \tilde{d}_0 \right]$, the time derivative of $\xi$ can be written as

where $A = \left[\begin{matrix} -h_1 & 1\\ -h_2 & 0\\ \end{matrix} \right]$, $B = \left[\begin{array}{c} 0\\ 1\\ \end{array} \right]$, $\varPhi \left(\tilde{\Omega} \right) = \frac{3}{2}\left| \tilde{\Omega} \right|^{\frac{1}{2}}+1$ and $\psi = \dot{d}_{0}/ \varPhi \left(\tilde{\Omega} \right)$.

With the help of $|{\dot d_0}| \le {l_0}$, one yields $|\psi | \le {l_0}$. Then, we define

Theorem 1: For a positive constant $\gamma$ and the symmetric and positive definite matrix $Q$, if the following inequality holds

then the estimation error of the proposed MESO (3.12) will converge to zero in a finite time.

Proof. The Lyapunov function is chosen as

From (3.15) and (3.16), the differential coefficient of $V_2$ can be given as

From (3.17), it is derived that

where $\lambda \left\{ \cdot \right\}$ is the eigenvalue of matrix $\left\{ \cdot \right\}$ and $\left\| \xi \right\|_2^2$ is the Euclidean norm of $\xi$.

Considering $\left\| \xi \right\|_2^2 = \left| \tilde{\Omega} \right|+2\left| \tilde{\Omega} \right|^{\frac{3}{2}}+\tilde{\Omega}^2+\tilde{d}_{0}^{2}$, it follows that

In accordance with (3.18) and (3.20), one obtains

Clearly, (3.21) satisfies the finite-time stability theory ^[50]. To this end, we have proved that under the proposed MESO (3.12) the estimation error will be guaranteed to converge to zero in a finite time.

The proposed MESO can estimate the disturbance precisely. Then, compensating the estimate $\hat{d}_0(t)$ to the baseline controller, the output signal of the MESO-based ANTSM controller can be designed as

4.   Experimental validation

Aiming to validate the performance of the ANTSM + MESO control algorithm, comprehensive experimental results are given in this section. The complete schematic of the proposed ANTSM + MESO control system is illustrated in Figure 4. The experimental platform mainly consists of a three-phase motor, a controller, a three-phase inverter, a host computer and more. The controlled motor is a permanent magnet synchronous motor with a power of 1.5 KW, and the motor parameters are illustrated in Table 1. The controller is based on the RTU-BOX204 real-time digital control platform. The magnetic powder brake is used to generate the load torque. The parameters of the ANTSM controller are chosen as $\beta = 600$, $q = 11$, $p = 17$, $\eta = 1.5$, $\varepsilon = 0.99$, $N = 80$, $k_m = 1$ and $k_M = 30$, and the parameters of MESO are selected as $h_1 = 30$ and $h_2 = 225$.

4.1.   Comparison of experimental results among PI, traditional NTSM and ANTSM

The start-up results of $\omega$, $i_q$ and $i_a$ under the PI, traditional NTSM and the proposed ANTSM controller are respectively depicted in Figure 5. The given speed is 500 $rpm$, and no additional load torque was added. One can notice that the start-up transient response under the PI has the longest convergence time and largest speed overshoot. In Figure 5(c), the developed ANTSM controller has the smaller start speed overshoot. Figure 6 shows a sudden load experiment. From the experimental results, one can see that the proposed ANTSM controller has the smallest speed fluctuation and current ripple.

4.2.   Comparison of experimental results among ANTSM+ESO and ANTSM+MESO

In this subsection, the ESO-based ANTSM controller is utilized to compare with the presented MESO-based ANTSM controller. Figure 7 exhibits the sudden load change responses at the speed of 500 $rpm$. After loading, the speed recovery time under the proposed ANTSM+MESO controller is shorter than that under the ANTSM+ESO controller, and the speed fluctuation under the proposed ANTSM+MESO controller is also smaller. The starting responses are depicted by Figure 8. One can undoubtedly observe that the speed overshoot of the proposed ANTSM+MESO controller is smaller than that under the ANTSM+ESO. Additionally, the time for the q-axis and a-phase currents to achieve the system stability was shorter than that under the ANTSM+ESO controller.

4.3.   Experimental results of dynamic performance tests

To further validate the availability of the proposed controller, numerous additional experiments were carried out. The moment of inertia has a straightforward impact on the starting and braking performance of the motor. As a result, the controller performance was tested by changing the value of inertia in the presented ANTSM+MESO controller. The inertia $J{\rm{ = }}\left({1/3} \right){J_0}$, $J{\rm{ = }}\left({1/2} \right){J_0}$ and $J{\rm{ = }}{J_0}$ were chosen for comparative experiments. Figure 9 shows the step speed responses at the three different inertia. The sudden load change responses at the speed of 500 $rpm$ are exhibited by Figure 10. As can be seen in Figures 9 and 10, the inertia mismatch imposes an impact on the stable operation of the motor. Nevertheless, the system can still operate stably under the proposed ANTSM+MESO controller.

The experimental results of speed, current $i_q$ and current $i_a$ during the speed change are shown in Figure 11(a). The anti-disturbance property of the ANTSM + MESO control strategy at 500 $rpm$ is exhibited in Figure 11(b). From Figure 11, one can summarize that the ANTSM + MESO algorithm can better control the motor operation at a wide range of speed.

5.   Conclusions

In this paper, a novel MESO-based ANTSM controller was constructed to improve the anti-disturbance performance of the PMSM drive system with parameter variation and unknown disturbance. To prevent the unsatisfactory control effect due to overestimating switching gain, an adaptive law was combined with the traditional NTSM strategy to achieve the expected performance. Furthermore, compensating the disturbance estimation value obtained by MESO to the ANTSM controller, the switching gain can be further reduced. Comprehensive experimental results demonstrate that the property of the MESO-based ANTSM algorithm outperforms both traditional PI and NTSM control methods. In future work, we will work on solving the stability problem caused by parameter variations.

Use of AI tools declaration

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

Acknowledgments

This paper was funded by the Natural Science Foundation of the Jiangsu Higher Education Institutions of China under grant number 21KJB510019, the National Natural Science Foundation of China under Grant 62203188, the Natural Science Foundation of Jiangsu Province under Grant BK20220517 and the China Postdoctoral Science Foundation under Grant 2022M721386.

Conflict of interest

The authors declare there is no conflict of interest.