Have a personal or library account? Click to login
A Custom Robotic Vehicle Prototype for Disinfecting Workplaces Cover
Open Access
|Sep 2025

Full Article

1.
Introduction

Control methods have been developed to stabilize practical systems. The control problems include regulation (stabilization), tracking, and path following. Linear quadratic regulation (LQR), (Proportional Integral Derivative) PID, and model predictive control (MPC) are the most popular control methods. These methods have some advantages and some disadvantages. One of the most essential parameters in control systems is defining the settling time and settling point. Some algorithms are developed to predict, determine, and calculate the settling time or settling point [14]. The predefined time algorithms can provide the upper bound of the settling (stabilization) time before applying the control signal to the system [59]. But no method provides a way to define the settling time and settling point simultaneously and determining the accurate settling time is another challenge.

Tuning the control parameters is another challenge in control systems. The control parameters should be tuned to satisfy the control goals. The PID controller is one of the most common control methods, but it has a tuning control parameters problem. There are some methods outlined in the literature that solve this problem [1014]. Another issue in PID that the user cannot define the settling time and be sure that the system’s output will reach the settling is point at the predefined settling time. The LQR technique provides an optimal solution by tuning the control parameters. However, the settling time of the controlled system remains an open research challenge.

A self-tuning control (STC) approach has been defined as a method that predicts/adjusts control parameters. These methods are known as “parameter tuners.” The parameter tuners are of two types: offline and online [1517]. The main contribution of the parameter tuners is increasing the control performance. In some cases, the parameter tuners are the main part of the control method within the control system. In these cases, the control method is called an STC method. One of the challenges in the STC methods is the system’s stability. In the parameter tuners, the system’s stability is not a problem, since the system’s stability should be guaranteed by the main control method. Fuzzy logic is one of the most common methods used to design the controllers as the parameter tuner [18]. The main controller in this study is a nonsingular PID finite-time sliding mode control (FSMC). The FSMC guarantees system stability, and the fuzzy logic system tunes the control parameters. In other studies [1820], fuzzy logic is used as the tuner, and the main controller is PID. In studies [2124], other types of tuners have been introduced to adjust the control parameters, which can be categorized as STC methods. These works have developed relay feedback systems for parameter tuning.

The settling time is an important parameter in the performance of a system. Different stabilization methods can guarantee various settling times [2530]. An asymptotic stabilization method requires an infinite amount of time to guarantee stability. However, a finite-time method provides stability with a limited upper-bound settling time [31]. Fixed-time stabilization methods are terminal methods that provide an upper bound on the settling time, independent of the system’s initial condition [32, 33]. Predefined-time stabilization methods provide fixed-timestability, where the settling time is a controller parameter [34]. Preselected-time control is defined as a method where the settling time is chosen in advance, ensuring that the system’s state reaches the preselected-state at the exact time. The time selected in the preselected-time approach can be considered the settling time when the preselected state (point) is chosen to be very close to zero.

In most cases in control systems, it is necessary to know/define the settling point and settling time together. Some studies have been published that present ways of calculating, determining, and defining the settling time and settling point separately, but defining the settling time and settling point together is an excellent feature for users. In this study, an STC method will be designed to tune the control parameter of the scalar systems. The proposed STC method provides the possibility of defining the settling point and settling time simultaneously. The Hamiltonian conditions will be used to prove the optimality of the method. The features of the proposed STC method are as follows:

  • -

    It is possible to select the settling point (preselected state) and settling time (preselected time) simultaneously.

  • -

    It uses a Lyapunov-based stabilization method.

  • -

    It is a practicable method that can be used in different applications.

  • -

    The proposed STC method provides an optimal solution for single-input single-output (SISO) linear (scalar) systems.

  • -

    It is a linear predefined-time solution for linear systems.

  • -

    It provides four possible scenarios.

  • -

    It can be updated for the tracking problem and (Multi-Input Multi-Output) MIMO systems.

  • -

    It provides an exponential solution without any unwanted overshoot and lower shoot for scalar systems where y = x.

2.
Problem Statement

A suitable controller for practical systems should have the following features.

Consider the following SISO linear scalar system: 1{ x˙=ax+buy=x,x(0)=x0 ,\left\{ {\matrix{ {\dot x = ax + bu} \cr {y = x} \cr } ,x(0) = {x_0}} \right., where x ∈ ℝ is the system’s state, y ∈ ℝ is the system’s output, a ∈ ℝ, b ∈ {ℝ-{0}} are the system’s parameters, u ∈ ℝ is the control input, and x0 ∈ ℝ is the initial condition of the system. The control input should guarantee the system’s stability and satisfy the discussed features. It is assumed that the system’s state is available for designing the controller. On the other hand, the system uses the feedback of the state.

Remark 1:

In practical and controllable systems, the system’s parameter b is not equal to 0.

The STC will guarantee that the system’s state reaches a specific point at a predefined time, which will be selected by the system’s user. This specific point and predefined time will be chosen simultaneously. Also, the solution will be optimal, which is proven by the Hamiltonian equations. The nonoptimal solution of the STC was published in [17] as an offline self-tuning controller (OSTC).

In OSTC, two performance criteria [integral of absolute error (IAE) and integral time absolute error (ITAE)] are considered for preparing some formulas to use in tuning the control parameters. In the OSTC method, some rules are predefined, and the control parameters are tuned accordingly before applying the controller to the system. However, in STC (online) methods, the control parameters are adjusted dynamically while controlling the system.

In this paper, the Hamiltonian equations will be used to design the optimal solution of the controller. Considering the system (1) and the following cost function: 2J=Qx2+Ru2dt,J = \int Q {x^2} + R{u^2}dt, where Q ∈ ℝ≥0 and R ∈ ℝ>0 are the weight parameters. The Hamiltonian equation for this system and cost function is as follows: 3H=Qx2+Ru2+λ(ax+bu),H = Q{x^2} + R{u^2} + \lambda (ax + bu), where λ is the costate variable. The Hamiltonian conditions are presented as follows [3537]: 4{ Hx=λHλ=x˙Hu=0. \left\{ {\matrix{ {{{\partial H} \over {\partial x}} = - \lambda } \cr {{{\partial H} \over {\partial \lambda }} = \dot x} \cr {{{\partial H} \over {\partial u}} = 0} \cr } .} \right.

The STC control input will be designed by using Hamiltonian conditions.

3.
Self-Tuning Controller Design

Considering the SISO scalar linear system of Equation (1) and the cost function of Equation (2), in order to have an optimal self-tuning solution for these conditions, the following equation is used: 5{ u=1b(ax+uSTC)uSTC=Px;P=ln(x0xS)TS, \left\{ {\matrix{ {u = {1 \over b}\left( { - ax + {u_{STC}}} \right)} \cr {{u_{STC}} = - Px;P = {{\ln \left( {{{{x_0}} \over {{x_S}}}} \right)} \over {{T_S}}}} \cr } ,} \right. where P is the STC parameter that is a positive constant as P=ln(x0xS)TsP = {{\ln \left( {{{{x_0}} \over {{x_S}}}} \right)} \over {{T_s}}}, where x0 is the initial condition of the system’s state, xs is the settling point (which we can choose and want to reach at the predefined settling time), and Ts is the predefined settling time. Figure 1 shows the block diagram of the proposed STC method.

Figure 1.

Block diagram of the STC method

Remark 2:

The concept of the STC method is that the user can choose the settling point (position) and settling time together. This is one of the greatest advantages of this method.

Remark 3:

A settling point is a value with the same sign as the initial condition, less than a positive initial condition, and greater than a negative initial condition. This fact causes the system not to have any unwanted overshoot or lower shoot.

Theorem 1:

The scalar linear system (1) will be stabilized as an optimal solution considering the cost function (2) by the self-tuning controller (5).

Proof:

The cost function of Equation (2) can be written as 6J=Qx2+Ru2dt=Qx2+Rb2(ax+uSTC)2dt.J = \int Q {x^2} + R{u^2}dt = \int Q {x^2} + R{b^{ - 2}}{\left( { - ax + {u_{STC}}} \right)^2}dt.

The Hamiltonian equation of system (1) and the cost function (6) is as follows: 7H=Qx2+Rb2(ax+uSTC)2+λuSTC.H = Q{x^2} + R{b^{ - 2}}{\left( { - ax + {u_{STC}}} \right)^2} + \lambda {u_{STC}}.

The Hamiltonian conditions of Equation (7) can be written as 8{ Hx=λ˙=2Qx+Rb2(2auSTC+2a2x),(A)Hλ=x˙=uSTC,(B)HuCTC=0=Rb2(2ax+2uSTC)+λ.(C) \left\{ {\matrix{ {{{\partial H} \over {\partial x}} = - \dot \lambda = 2Qx + R{b^{ - 2}}\left( { - 2a{u_{STC}} + 2{a^2}x} \right),} \hfill & {(A)} \hfill \cr {{{\partial H} \over {\partial \lambda }} = \dot x = {u_{STC}},} \hfill & {(B)} \hfill \cr {{{\partial H} \over {\partial {u_{CTC}}}} = 0 = R{b^{ - 2}}\left( { - 2ax + 2{u_{STC}}} \right) + \lambda .} \hfill & {(C)} \hfill \cr } } \right.

These conditions can be simplified as follows: 9 from (C):λ=Rb2(2ax+2uSTC)λ˙=Rb2(2ax˙+2u˙STC);(D)\matrix{ \hfill {{\rm{ from }}(C):\lambda = - R{b^{ - 2}}\left( { - 2ax + 2{u_{STC}}} \right)} \cr \hfill { \Rightarrow \dot \lambda = - R{b^{ - 2}}\left( { - 2a\dot x + 2{{\dot u}_{STC}}} \right);\;(D)} \cr } from – (D) = (A) we have 10Rb2(2auSTC+2u˙STC)=2Qx+Rb2(2auSTC+2a2x);R{b^{ - 2}}\left( { - 2a{u_{STC}} + 2{{\dot u}_{STC}}} \right)\quad = 2Qx + R{b^{ - 2}}\left( { - 2a{u_{STC}} + 2{a^2}x} \right); after a simplification can be written 11Rb2u˙STC=Qx+Rb2a2x.R{b^{ - 2}}{{\dot u}_{STC}} = Qx + R{b^{ - 2}}{a^2}x.

As u˙STC=Px˙{{\dot u}_{STC}} = - P\dot x, and x˙=uSTC\dot x = {u_{STC}}, so we can write 12Rb2Px˙=Qx+Rb2a2xRb2PuSTC=Qx+Rb2a2x. - R{b^{ - 2}}P\dot x = Qx + R{b^{ - 2}}{a^2}x \Rightarrow - R{b^{ - 2}}P{u_{STC}} = Qx + R{b^{ - 2}}{a^2}x.

After a simplification, 13Rb2P2x=Qx+Rb2a2xRb2P2=Q+Rb2a2.R{b^{ - 2}}{P^2}x = Qx + R{b^{ - 2}}{a^2}x \Rightarrow R{b^{ - 2}}{P^2} = Q + R{b^{ - 2}}{a^2}.

The STC parameter will be achieved as follows: 14P=QRb2+a2.P = \sqrt {{Q \over R}{b^2} + {a^2}} .

As defined previously, the STC parameter is equal to P=ln(x0xs)TsP = {{\ln \left( {{{{x_0}} \over {{x_s}}}} \right)} \over {{T_s}}}, and from the Hamiltonian, it is calculated as P=QRb2+a2P = \sqrt {{Q \over R}{b^2} + {a^2}} . These equations can be equal as 15P =ln(x0xs)Ts=QRb2+a2QR =((ln(x0xs)Ts)2a2)b2.P = {{\ln \left( {{{{x_0}} \over {{x_s}}}} \right)} \over {{T_s}}} = \sqrt {{Q \over R}{b^2} + {a^2}} \Rightarrow {Q \over R} = \left( {{{\left( {{{\ln \left( {{{{x_0}} \over {{x_s}}}} \right)} \over {{T_s}}}} \right)}^2} - {a^2}} \right){b^{ - 2}}.

This result is the self-tuning tool that allows us to define different scenarios. These scenarios will be discussed in the next section.

The Lyapunov theory is employed to prove the stability of the self-tuning controller. The following Lyapunov function can be defined: 16V=12x2V˙=xx˙V˙=x(ax+bu).V = {1 \over 2}{x^2} \Rightarrow \dot V = x\dot x \Rightarrow \dot V = x(ax + bu).

After applying the control input to the system, we have the following: 17V˙=xuSTCV˙=Px2V˙=2PVV˙=ρ1V0,\matrix{ {\dot V = x{u_{STC}} \Rightarrow \dot V = - P{x^2}} \cr { \Rightarrow \dot V = - 2PV \Rightarrow \dot V = - {\rho _1}V \le 0,} \cr } where ρ1 = 2P is a positive value. Therefore, the stability of the system is proven by the Lyaponuv stability theory [38].

To prove the fact that the self-tuning controller guarantees that the chosen settling point will be reached at the predefined settling time, we can solve the system’s differential equation as follows: 18x˙=uSTC=Pxx(t)=x0ePt{ xS=x(Ts)=x0ePTsorTs=ln(x0xs)PorP=ln(x0xs)Ts .\matrix{ {\dot x} \hfill & { = {u_{STC}} = - Px \Rightarrow x(t)} \hfill \cr {} \hfill & { = {x_0}{e^{ - Pt}}\left\{ {\matrix{ {{x_S} = x\left( {{T_s}} \right) = {x_0}{e^{ - P{T_s}}}} \cr {or} \cr {{T_s} = {{\ln \left( {{{{x_0}} \over {{x_s}}}} \right)} \over P}} \cr {or} \cr {P = {{\ln \left( {{{{x_0}} \over {{x_s}}}} \right)} \over {{T_s}}}} \cr } } \right..} \hfill \cr }

Therefore, theorem 1 is proven, and the settling point will be reached at the settling time by applying the STC control signal to the system.

The controlled system’s response is x(t) = x0ePt, which exhibits exponential behavior and ensures stability, as limt→∞ x = 0.

Remark 4:

Relations in Equation (18) show the relationship between the prede^ined settling time, settling point, and STC parameter. These relations can be used to calculate the proper values of the parameters by de^ining (selecting) two other parameters.

Remark 5:

As evident in Equation (18), the system’s state and control signal of the STC method will change exponentially. It is provable that for every STC solution of the SISO linear system, there is an optimal solution with a speci^ied cost function. By the terms of Equation (15), the weight parameters of the speci^ied cost function can be calculated.

Remark 6:

The presented STC method can be extended to MIMO linear systems by decoupling the first-order SISO subsystems as follows: 19A=diag(a1,a2,,an)=[ a100an ]B=diag(b1,b2,,bn)=[ b100bn ].\matrix{ {A = {\mathop{\rm diag}\nolimits} \left( {{a_1},{a_2}, \ldots ,{a_n}} \right) = \left[ {\matrix{ {{a_1}} & \cdots & 0 \cr \vdots & \ddots & \vdots \cr 0 & \cdots & {{a_n}} \cr } } \right]} \hfill \cr {B = {\mathop{\rm diag}\nolimits} \left( {{b_1},{b_2}, \ldots ,{b_n}} \right) = \left[ {\matrix{ {{b_1}} & \cdots & 0 \cr \vdots & \ddots & \vdots \cr 0 & \cdots & {{b_n}} \cr } } \right]} \hfill \cr } .

4.
Different Possible Scenarios

Considering Equation (15), four scenarios can be envisioned for the proposed STC method. In this section, these scenarios will be discussed.

4.1.
First Scenario

The main goal of the STC method is that the settling point (xs) and predefined settling time (Ts) occur simultaneously.

In the first scenario, the user should select the settling point, settling time, and weight parameter R in the cost function; then the STC will calculate the control parameter P and the weight parameter of Q.

4.2.
Second Scenario

The second scenario is similar to the first one, but in this scenario, the weight parameter Q should be selected by the user, and R and P will be calculated by the STC.

Remark 7:

In general, the weight parameters of the cost function can be positive or negative, but in most cases, a positive value for these parameters is required. In order to have positive values of the weigh parameters (Q and R) in the first and second scenarios, one of the following conditions is necessary: 20{ Ts<| ln(x0xS)a |orxs<| x0e| aTs | | .\left\{ {\matrix{ {{T_s} < \left| {{{\ln \left( {{{{x_0}} \over {{x_S}}}} \right)} \over a}} \right|} \cr {or} \cr {{x_s} < \left| {{{{x_0}} \over {{e^{\left| {a{T_s}} \right|}}}}} \right|} \cr } } \right..

4.3.
Third Scenario

In this scenario, the weight parameters of the cost function Q and R and the settling point xs should be defined by the user, and the settling time Ts and control parameter f will be calculated by the STC. In this scenario, the weight parameters should be positive.

4.4.
Fourth Scenario

This scenario is similar to the third scenario, but in this scenario, the settling time and weight parameters should be selected, and the control parameter and settling point will be calculated.

In all scenarios, the control parameter will be calculated by the defined condition (other parameters); then the control parameter will be applied to the system (1). In the previous section, the system stability of the STC method was proven.

5.
Simulations

In this section, three practices will be defined and simulated for every scenario. First, it is necessary to determine the system’s parameters and the desired parameters; then the STC will calculate the rest of the parameters and the value of the cost function.

MATLAB software has been used to simulate the scenarios with the ode4 solver, and the step time equals 0.01. The calculations have been done up to two decimal places.

For all practices, the values of the cost functions for the STC, LQR, and proportional (P) controllers are presented, and the results are compared with those of the LQR and P controllers tuned using the PID block in Simulink/MATLAB. It is expected that the comparison results for the STC and LQR will be the same.

5.1.
Simulations for the First Scenario

The determined scalar linear system is as follows: 21x˙=4.5x+0.6u,x(0)=10.\dot x = - 4.5x\; + \;0.6u,\;x(0)\; = \;10.

Three practices are defined to test the proposed STC algorithm. The simulation results are presented in Table 1. Figures 2, 3, and 4 show the simulation results of the different practices. In the figures, the selected settling point and settling time are shown. The P controller has been tuned using the PID block Simulink/MATLAB. The P controller is tuned as uPID = Kpx, where Kp = 15.3191 for this scenario.

Figure 2.

Simulation results of the first practice of the first scenario

Figure 3.

Simulation results of the second practice of the first scenario

Figure 4.

Simulation results of the third practice of the first scenario

Table 1.

Parameters in different practices of the first scenario

PracticesParameters Selected by the UserParameters Calculated by the STC MethodCalculated Values of the Cost Functions
xsTsRQPJSTCJLQRJPID
1st80.11013268.9022.315.11 × 1045.11 × 1045.70 × 104
2nd0.50.5143.475.991428.9428.91015.8
3rd0.0110.17.636.90868.0868.08113.57
5.2.
Simulations for the Second Scenario

The determined scalar linear system for the second scenario is as follows: 22x˙=0.01x+0.2u,x(0)=5.\dot x\; = \;0.01x\; + \;0.2u,\;x(0)\; = \; - 5.

The simulation results of the second scenario are presented in Table 2. Figures 5, 6, and 7 show the simulation results of the different practices. In the figures, the selected settling point and settling time are shown.

Figure 5.

Simulation results of the first practice of the second scenario

Figure 6.

Simulation results of the second practice of the second scenario

Figure 7.

Simulation results of the third practice of the second scenario

Table 2.

Parameters in different practices of the second scenario

PracticesParameters Selected by the UserParameters Calculated by the STC MethodCalculated Values of the Cost Functions
xsTsQRPJSTCJLQRJPID
1st–128.0930.50.8047263.9263.9513.41
2nd–40.02311206.210011.16720678.15720678.151366900
3rd–0.011.54.2910.014.14326.5326.5327.1698

The P controller has been tuned using the PID block Simulink/MATLAB. The P controller is tuned as uPID = –Kpx, where Kp = 15.3191 for this scenario.

5.3.
Simulations for the Third Scenario

The selected scalar linear system for this scenario is as follows: 23x˙=0.025x+2.6u,x(0)=1.5.\dot x = \;0.025x\; + \;2.6u,\;x(0)\; = \;1.5.

The simulation results are presented in Table 3. Figures 8, 9, and 10 show the simulation results of the different practices. In the figures, the selected settling point and settling time are shown. The P controller has been tuned using the PID block Simulink/MATLAB. The P controller is tuned as uPID = –Kpx, where Kp = 15.3191 for this scenario.

Figure 8.

Simulation results of the first practice of the third scenario

Figure 9.

Simulation results of the second practice of the third scenario

Figure 10.

Simulation results of the third practice of the third scenario

Table 3.

Parameters in different practices of the third scenario

PracticesParameters Selected by the UserParameters Calculated by the STC MethodCalculated Values of the Cost Functions
xsQRTsPJSTCJLQRJPID
1st0.5110.402.60.89970.89976.6644
2nd0.0140.20.4211.630.80220.80221.4403
3rd10.040.20.351.1630.081110.081111.3284
5.4.
Simulations for the Fourth Scenario

The determined scalar linear system for the fourth scenario is as follows: 24x˙=0.1x0.6u,x(0)=3.\dot x\; = \; - 0.1x\; - \;0.6u,\;x(0)\; = \; - 3.

The simulation results are presented in Table 4. Figures 11, 12, and 13 show the simulation results of the different practices. In the figures, the selected settling point and settling time are shown. The P controller has been tuned using the PID block Simulink/MATLAB. The P controller is tuned as uPID = –Kpx, where Kp = –0.3404252 for this scenario.

Figure 11.

Simulation results of the first practice of the fourth scenario

Figure 12.

Simulation results of the second practice of the fourth scenario

Figure 13.

Simulation results of the third practice of the fourth scenario

Table 4.

Parameters in different practices of the fourth scenario

PracticesParameters Selected by the UserParameters Calculated by the STC MethodCalculated Values of the Cost Functions
QRTsxsPJSTCJLQRJPID
1st423–0.230.854438.3938.3962.4464
2nd0.016010–1.10.10030.45630.4563102.9896
3rd110.01–2.9820.608313.0513.0516.4667

For the four possible scenarios, three practices are defined and tested under various conditions. In all practices across all scenarios, the selected or calculated settling time and settling point are accurately achieved. The STC is compared with the LQR and P controllers, and the cost function calculations indicate that the STC achieves the same optimal results as the LQR while performing significantly better than the P controller.

6.
Conclusion

This paper presents a new control strategy named a self-tuning controller. The proposed STC can stabilize the scalar linear system. This method can be updated for MIMO linear systems. The simulation results show the power of the proposed method. The STC is compared with the LQR method, and the results demonstrate that STC provides an optimal solution for the system. Also, it can control the system to reach the settling point at the predefined settling time. Future works can focus on developing a MATLAB toolbox for the STC method and developing the STC for MIMO systems and tracking problems.

DOI: https://doi.org/10.14313/jamris-2025-022 | Journal eISSN: 2080-2145 | Journal ISSN: 1897-8649
Language: English
Page range: 25 - 37
Submitted on: Aug 9, 2024
Accepted on: Feb 28, 2025
Published on: Sep 10, 2025
Published by: Łukasiewicz Research Network – Industrial Research Institute for Automation and Measurements PIAP
In partnership with: Paradigm Publishing Services
Publication frequency: 4 issues per year

© 2025 Ali Soltani Sharif Abadi, published by Łukasiewicz Research Network – Industrial Research Institute for Automation and Measurements PIAP
This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License.