Next Article in Journal
Indoor Visible Light Positioning System Based on Point Classification Using Artificial Intelligence Algorithms
Previous Article in Journal
Design, Development and Testing of a Monitoring System for the Study of Proton Exchange Fuel Cells and Stacks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonlinear Robust Control by a Modulating-Function-Based Backstepping Super-Twisting Controller for a Quadruple Tank System

by
Italo Aranda-Cetraro
1,*,
Gustavo Pérez-Zúñiga
1,*,
Raul Rivas-Pérez
1,2 and
Javier Sotomayor-Moriano
1
1
Departamento de Ingeniería, Pontificia Universidad Católica del Perú (PUCP), Avenida Universitaria 1801, San Miguel, Lima 15088, Peru
2
Departamento de Automática y Computación, Universidad Tecnológica de la Habana José Antonio Echeverría (CUJAE), La Habana 19390, Cuba
*
Authors to whom correspondence should be addressed.
Sensors 2023, 23(11), 5222; https://doi.org/10.3390/s23115222
Submission received: 21 April 2023 / Revised: 15 May 2023 / Accepted: 27 May 2023 / Published: 31 May 2023
(This article belongs to the Topic Industrial Control Systems)

Abstract

:
In this paper, a robust nonlinear approach for control of liquid levels in a quadruple tank system (QTS) is developed based on the design of an integrator backstepping super-twisting controller, which implements a multivariable sliding surface, where the error trajectories converge to the origin at any operating point of the system. Since the backstepping algorithm is dependent on the derivatives of the state variables, and it is sensitive to measurement noise, integral transformations of the backstepping virtual controls are performed via the modulating functions technique, rendering the algorithm derivative-free and immune to noise. The simulations based on the dynamics of the QTS located at the Advanced Control Systems Laboratory of the Pontificia Universidad Católica del Perú (PUCP) showed a good performance of the designed controller and therefore the robustness of the proposed approach.

1. Introduction

Currently, there are a considerable number of industrial multivariable processes with complex nonlinear dynamic behavior [1,2,3]. Therefore, multivariable control strategies are highly important and have received significant attention from the international scientific community; see, for example, [4,5,6,7,8,9].
The quadruple tank system (QTS) has been widely used in academia to design multi-input, multi-output (MIMO) control schemes for liquid level regulation in the presence of complex nonlinear dynamics [10]. Classical control approaches to this problem involve linearizing the model about an operating point (OP) and using decouplers to eliminate or reduce the interaction sensitivity between inputs and outputs, in order to apply decentralized proportional-integral (PI) control [11]. In [12], an actuator fault-tolerant decentralized PI-controller based on the design of simplified decouplers was proposed, allowing feed-forward control, as if it were two independent single input single output (SISO) systems [13]. The fault-tolerant part of the controller had the ability to additively compensate for the magnitude of the fault. Yet, one of the downsides of this approach was that the designed decouplers were not always physically feasible and were prone to modeling errors [12].
In [14], the performance of a decentralized PI controller was compared to robust control strategies such as a multivariable internal model control (IMC) and an H-infinity control ( H ), concluding that these control strategies yielded better performance than the decentralized PI controller. In [15], a reconfigurable model predictive control (MPC) approach was followed, where the feasibility to compensate for the control signal was shown when an actuator (pump) failure was detected by switching its control signal to another actuator and switching to the output of an observer when a sensor fault was detected. In general, all these methods limited their overall performance to an operating point, since they worked with a linearized model of the system.
However, other authors proposed working directly with the variable structure control (VSC) of nonlinear multivariable systems. Variable structure control systems (VSCS) theory, based on the phase-plane method of the oscillatory theory, was proposed by Alexander Andronov in the 1940s [16]. This control theory was formally presented by Stanislav Emelyanov in the 1960s and developed by Vadim Utkin [17] and other authors afterwards. A VSCS is composed of continuous subsystems and a commutation law between these subsystems, originating discontinuous or bang–bang control efforts to stabilize or regulate the response of the processes.
Sliding mode control (SMC) is a special type of VSCS, since it introduces an error variable named a “sliding variable” to steer the trajectory of the system to a sliding manifold and maintain motion on the manifold by means of discontinuous control, regardless of disturbances or uncertainties to the process [18]. Nonetheless, due to its inherent robustness, the SMC introduces chattering or high-frequency oscillations to the control, which is undesirable for actuators. Various chattering-suppression methods have been proposed [19], with the equivalent control method [17] and higher-order sliding modes (HOSM) [20] as the most representative. Moreover, in [21] it was proved that the first-degree SMC and HOSM were sensitive to non-Gaussian measurement noise, making it necessary to test and implement special HOSM differentiator algorithms and filters to mitigate the influence of the nonlinear noise distributions in the sliding variable for real applications.
Concerning the use the use of the SMC and HOSM algorithms for liquid level regulation of the QTS, in [22], an SMC technique through feedback linearization was proposed, yielding better performance than the conventional PI controller. In [23], a second-order sliding mode (SOSM) controller based on the twisting algorithm (TA) was designed to regulate the liquid levels while considerably reducing the chattering level of the control effort. Even though the SOSM and HOSM methods proved to be effective in reducing the chattering level, they did not completely suppress it in some cases.
In [24], the active disturbance rejection control (ADRC) was introduced as an alternate approach to PID control, based on the design of a tracking differentiator, which provided the error signal, the derivative of the error signal, and a sliding mode controller (SMC). Recently, improvements were made in [25] with regard to the ADRC based on the work provided in [24] with promising results, although the performance of the tracking differentiator was not addressed in the presence of measurement noise.
The main contribution of this paper consists in the proposal of a robust nonlinear approach for control of liquid levels in a quadruple tank system (QTS) based on the combination of a backstepping controller and a super-twisting controller (BSSTC), implementing a multivariable sliding surface, where the error trajectories converge to the origin at any operating point of the system. Because the backstepping algorithm is dependent on the derivatives of the state variables and is sensitive to measurement noise, integral transformations of the backstepping virtual controls are carried out using modulating functions technique, making the algorithm derivative-free and immune to noise. All the modeling and control methodologies shown in this paper have been developed for a real laboratory QTS.
This paper proceeds as follows. In Section 2, a mathematical model of the process studied is obtained using modeling techniques. The theoretical background is explained in Section 3. The design of the MF-BSSTC controller is developed in Section 4. The discussions of the attained results are presented in Section 5. Lastly, a few conclusions are given in Section 6.

2. Quadruple Tank System Modeling

The study presented in this paper is based on the QTS located at PUCP’s Advanced Control Systems Laboratory. This QTS consists of four coupled tanks, a reservoir tank, four ball valves, two proportional valves, and two pumps. Figure 1 shows a view of this QTS, and its schematic representation is exhibited in Figure 2.

QTS Modeling under Disturbances

The control objective is to regulate the liquid levels of tanks n°3 and n°4. The inputs to the system are the flow rates u 1 and u 2 from pumps n°1 and n°2, respectively, measured by the flow transmitters F T 1 and F T 2 . The measured outputs are the liquid levels h 1 , h 2 , h 3 , and h 4 of tanks n°1, n°2, n°3, and n°4, respectively, which are proportional to the voltages generated by the level transmitters L T 1 , L T 2 , L T 3 , and L T 4 . The system has two proportional valves V p 1 and V p 2 with aperture percentages k 1 and k 2 [ 0 , 1 ] , respectively, two ball valves V b 2 and V b 4 , which generate flow towards tanks n°2 and n°4 with a split constant γ 1 [ 0 , 1 ] , and two ball valves V b 1 and V b 3 which generate flow towards tanks n°1 and n°3 with a split constant γ 2 [ 0 , 1 ] . Since the output flow rates of pump n°1 split between tank n°2 and tank n°4 and the output flow of pump n°2 between tank n°1 and tank n°3, the position of valves γ 1 and γ 2 controls the split ratio. For instance, if 0 < γ 1 + γ 2 < 1 holds, one transmission zero out of the two transmission zeros of the system locates at the right half plane (RHP) of the root locus, generating nonminimum phase dynamics. On the other hand, if 1 < γ 1 + γ 2 < 2 holds, both transmission zeros locate at the left half plane (LHP) of the locus, rendering the QTS with minimum phase dynamics.
In [14], mass balances and Bernoulli’s law yield the following multivariable QTS model with disturbances:
h ˙ 1 = a 1 2 g h 1 A 1 + ( 1 γ 2 ) k 2 u 2 A 1 + ξ 1 , h ˙ 2 = a 2 2 g h 2 A 2 + ( 1 γ 1 ) k 1 u 1 A 2 + ξ 2 , h ˙ 3 = a 3 2 g h 3 A 3 + a 2 2 g h 2 A 3 + γ 2 k 2 u 2 A 3 + ξ 3 , h ˙ 4 = a 4 2 g h 4 A 4 + a 1 2 g h 1 A 4 + γ 1 k 1 u 1 A 4 + ξ 4 , y = h 1 h 2 h 3 h 4 + v ,
where x = h 1 h 2 h 3 h 4 represents the state variables vector, h i i = 1 , 2 , 3 , 4 is the tank i liquid level in cm, y is the measurement vector with v as additive noise, u = u 1 u 2 is the control vector, A i is the cross section of tank i in cm 2 , a i is the cross section of the outlet hole of the tank i in cm 2 , d i is the diameter of the tank i in cm, h i m a x is the maximum liquid level of tank i, and g is the acceleration of gravity in cm/s 2 . Moreover, ξ 1 = ( 1 γ 2 ) f p 2 A 1 , ξ 2 = ( 1 γ 1 ) f p 1 A 2 , ξ 3 = γ 2 f p 2 A 3 , and ξ 4 = γ 1 f p 1 A 4 are disturbances to the system generated by the flow losses f p 1 and f p 2 of pump n°1 and pump n°2, respectively, which are modeled as step functions, such that ξ i ˙ = 0 . In Table 1, the parameter values of the four-tank coupled system are presented.

3. Theoretical Background

Sliding mode control (SMC) has the ability to reject bounded matched uncertainties at the cost of introducing chattering to the control input. This could be detrimental to the performance if the mechanical systems are controlled, as is the case with QTS. In fact, one of the many methodologies available to suppress chattering is the design of super-twisting algorithm (STA) based controllers [20], which has been extensively used in recent years for closed-loop control of the QTS.
A super-twisting controller (STC) is a type of SOSM applicable to a system, where the control appears in the first derivative of the sliding variable [26], which has the ability to compensate for disturbances or uncertainties with only the knowledge of the measured output or sliding variable σ , while suppressing or attenuating chattering. For instance, the STC could be applied to the QTS if a multivariable STC scheme was designed, considering that there should exist as many sliding surfaces as independent controls [27], and some “hierarchy of controls” [17] should be established.
Definition 1 
(Multivariable Sliding Surface [27]). Let the general state-space representation of a system with multiple controls be
x ˙ = f ( x ) + g ( x ) u , y = h ( x ) ,
where x R n is the state variable vector of dimension n, u R m is the control vector of dimension m, and y R n is the output vector of dimension n, where n and m represent the number of controlled systems and the number of sliding surface coordinate functions defined as system outputs, respectively. Then, the multivariable sliding surface is represented by the simultaneous satisfaction of m smooth algebraic state restrictions, summarized in the equation σ ( x ) = 0 , which represents the intersection manifold,
S = { x R n | σ i ( x ) = 0 , i = 1 , 2 , . . , m } = i = 1 m S i .
For instance, the following vector of sliding surfaces was chosen for the QTS,
σ ( x ) = σ 1 ( x ) σ 2 ( x ) T = c 1 e 4 c 2 e 3 T ,
where e i = h i h i are the errors between the actual liquid level h i and the desired level h i , c i > 0 are constants, and i = 1 , 2 , 3 , 4 . h i must be a continuously differentiable reference trajectory. If the reference trajectory is a step, then h i = h ˙ i = = 0 .
Therefore, the simultaneous satisfaction of algebraic constraints σ 1 ( x ) and σ 2 ( x ) , which geometrically represents the existence of a smooth intersection manifold S, ideally produces a desired closed-loop behavior of the system, where xS holds in finite time [26].
In order to satisfy the algebraic constraints shown above, a first-degree SMC and a BSSTC are designed.

Synthesis of a First-Degree Sliding Mode Controller (SMC)

In this subsection a first-degree SMC is designed using the equivalent control method proposed in [17,27].
Definition 2 
(Equivalent Control and Ideal Sliding Dynamics [27]). Let the lie derivative L g σ ( x ) be locally invertible, where σ ( x ) is a vector of the sliding surfaces that satisfies σ ( x ) = 0 . Then,
σ ˙ ( x ) = σ ( x ) x T ( f ( x ) + g ( x ) u e q ( x ) ) = 0 ,
or
σ ˙ ( x ) = L f σ ( x ) + [ L g σ ( x ) ] u e q ( x ) | σ = 0 = 0 .
Therefore, the equivalent control is expressed as
u e q ( x ) = [ L g σ ( x ) ] 1 L f σ ( x ) | σ = 0 ,
with ideal sliding dynamics given by:
x ˙ = f ( x ) G ( x ) [ L G σ ( x ) ] 1 L f σ ( x ) .
The equivalent control is the smooth feedback control law, denoted by u e q ( x ) , which ideally locally holds the state evolution in the smooth manifold S for any initial state of the system located locally in S [26]. However, the closed-loop controller obtained with the equivalent control method generates a phenomena called “chattering”, as mentioned above, which can potentially wear out actuators when applied. Thus, in this paper, two second-order super-twisting controllers (2-STC) are designed through the recursive nonlinear backstepping technique to generate smooth control actions to regulate the liquid level of the selected subsystems, in order to reach the smooth intersection manifold S at any point of operation without the need to linearize the system.
Taking the model (1) into the state-space representation introduced in (2), the matrices f ( x ) and g ( x ) are shown:
f ( x ) = a 1 2 g x 1 A 1 a 2 2 g x 2 A 2 a 3 2 g x 3 A 3 + a 2 2 g x 2 A 3 a 4 2 g x 4 A 4 + a 1 2 g x 1 A 4 ,
g ( x ) = 0 ( 1 γ 2 ) k 2 A 1 ( 1 γ 1 ) k 1 A 2 0 0 γ 2 k 2 A 3 γ 1 k 1 A 4 0 ,
where x i = h i i = 1 , 2 , 3 , 4 .
Let the sliding surfaces of the multivariable closed-loop control be σ 1 = x 4 x 4 and σ 2 = x 3 x 3 . So, the lie derivatives of σ 1 and σ 2 along the direction of vector fields f ( x ) and g ( x ) are computed as follows:
L f σ ( x ) = L f σ 1 ( x ) L f σ 2 ( x ) = σ 1 ( x ) x T f ( x ) σ 2 ( x ) x T f ( x ) = 0 0 0 1 0 0 1 0 a 1 2 g x 1 A 1 a 2 2 g x 2 A 2 a 3 2 g x 3 A 3 + a 2 2 g x 2 A 3 a 4 2 g x 4 A 4 + a 1 2 g x 1 A 4 , L f σ ( x ) = a 4 2 g x 4 A 4 + a 1 2 g x 1 A 4 a 3 2 g x 3 A 3 + a 2 2 g x 2 A 3 ,
L g σ ( x ) = L g σ 1 ( x ) L g σ 2 ( x ) = σ 1 ( x T ) x g ( x ) σ 2 ( x T ) x g ( x ) = 0 0 0 1 0 0 1 0 0 ( 1 γ 2 ) k 2 A 1 ( 1 γ 1 ) k 1 A 2 0 0 γ 2 k 2 A 3 γ 1 k 1 A 4 0 , L g σ ( x ) = γ 1 k 1 A 4 0 0 γ 2 k 2 A 3 .
Then, replacing (11) and (12) in Equation (7), the equivalent control u e q to ensure σ ( x ) 0 in finite time is found as follows:
u e q ( x ) = a 4 2 g x 4 a 1 2 g x 1 γ 1 k 1 a 3 2 g x 3 a 2 2 g x 2 γ 2 k 2 ,
with the following control law:
u ( x ) = | η | 1 2 ( 1 1 s i g n ( σ ( x ) T L g σ ( x ) ) T ,
where | η | = 1718   cm 3 / s is the upper saturation threshold value for the flow rate, which equates to 60% of the maximum capacity of the pump. Furthermore, this control law u ( x ) = 0 takes the lower saturation threshold value of 10% of the pump capacity when it is equal to 0.
Proof. 
Robust closed-loop stability: Replacing Equation (13) in model (1), the following closed-loop dynamics are obtained:
x ˙ 1 = γ 2 a 1 2 g x 1 + ( a 3 γ 2 a 3 ) 2 g x 3 ( a 2 γ 2 a 2 ) 2 g x 2 A 1 γ 2 , x ˙ 2 = γ 1 a 2 2 g x 2 + ( a 4 γ 1 a 4 ) 2 g x 4 ( a 1 γ 1 a 1 ) 2 g x 1 A 2 γ 1 , x ˙ 3 = 0 , x ˙ 4 = 0 .
The closed-loop system (15) shows that x 3 and x 4 will reach the intersection manifold S in finite time. □

4. Design of the Modulating-Function-Based Backstepping Super-Twisting Controller

Multivariable control of the QTS based on the synthesis of a controller by the equivalent control method is feasible. Nonetheless, the implementation of the first-order sliding modes in the QTS could be detrimental to the plant actuators due to the chattering and the sensitivity to the measurement noise. On the other hand, using the backstepping and modulating functions technique can in fact render the controller chattering-free and robust to nonlinear noise. In this section, a modulating-function-based backstepping super-twisting controller is designed for the first time.
The backstepping technique [28] is a recursive back-deduced Lyapunov-based approach for systems transformable in their parametric-strict-feedback form or their pure parametric feedback form. This technique uses some of the system state variables as “virtual controls” at each step of the algorithm, implementing intermediate control laws to stabilize the system energy [28]. The advantage of applying backstepping control (BSC) is that it avoids the cancellation of the nonlinearities that are useful for controlling the system, easing the control effort. Moreover, the BS technique provides a framework to develop adaptive laws to unknown process parameters [29,30] and calibrate the gains of the online sliding mode controller, to efficiently compensate for parametric uncertainties and disturbances. However, its dependence on the derivatives of the state variables limits its application in plants or processes where these derivatives are not measurable. In general, recent nonlinear SMC and BSC approaches for liquid level control of the QTS [22,23,24,25,31,32] have been proved to be efficient only in conditions where the derivatives of the state variables are available and in the absence of measurement noise.
Design of the backstepping super-twisting controller (BSSTC) for tank n°4. The dynamics of the QTS are reduced only by considering the dynamics of the tank n°4. In order to find a suitable STC through the backstepping technique, it is necessary to extend the dynamics of tank n°4 by adding an auxiliary input w 1 as follows:
h ˙ 4 = A h 4 + B h 1 + C u 1 + ξ 4 , u 1 ˙ = w 1 ,
where A = a 4 2 g A 4 , B = a 1 2 g A 4 , and C = γ 1 k 1 A 4 . For the extended subsystem above, the following coordinate change is introduced:
x ˙ 1 = A 2 + B 2 x 1 x 2 + C 2 x 1 u 1 + 1 2 x 1 ξ 4 , u 1 ˙ = w 1 ,
where x 1 = h 4 , x 2 = h 1 , and u ˙ 1 = w 1 is an auxiliary input. By adding u 1 ˙ = w 1 , system (19) is already in pure parametric feedback form [28]. The control objective is to design a continuous controller that regulates the liquid level of tank n°4 at any operating point of the plant, achieving a SOSM ( σ i = σ ˙ i = 0 i = 1 , 2 ) in finite time over the sliding surface σ i ( x ) .
Step 1. Starting with the dynamics of x 1 in system (17), let us define the new error coordinate,
z 1 = x 1 x 1 ,
where x 1 is a twice-continuously differentiable reference trajectory. However, if the reference trajectory x 1 is a step function, x ˙ 1 = x ¨ 1 = x 1 = 0 holds. So, the derived dynamics of the new coordinate are:
z ˙ 1 = A 2 + B 2 x 1 x 2 + C 2 x 1 u 1 + 1 2 x 1 ξ 4 x ˙ 1 .
Let x 2 be an internal control variable, α 1 a virtual control law, and z 2 = x 2 α 1 the error between the actual control variable and the virtual control. Then, the control objective is to design the virtual control law α 1 , such that z 1 0 in finite time. In order to stabilize σ 1 0 , the following candidate Lyapunov function is introduced:
V 1 = 1 2 z 1 2 ,
and its time derivative:
V ˙ 1 = z 1 z ˙ 1 .
Replacing Equation (21) in (23) yields:
V ˙ 1 = z 1 ( A 2 + B 2 x 1 α 1 + C 2 x 1 u 1 + 1 2 x 1 ξ 4 x ˙ 1 ) + B 2 x 1 z 2 z 1 .
If z ˙ 1 = c 1 z 1 , where c 1 > 0 must hold, then
V ˙ 1 = c 1 z 1 2 + B 2 x 1 z 2 z 1 < 0 ,
if and only if z 2 = 0, and z 1 is locally asymptotically stable. To this end, virtual control law α 1 is found:
α 1 = A B x 1 C B u 1 1 B ξ 4 + 2 B Φ 1 2 c 1 B Φ 2 ,
with dynamics:
α ˙ 1 = A B x ˙ 1 C B u ˙ 1 1 B ξ ˙ 4 + 2 B Φ ˙ 1 2 c 1 B Φ ˙ 2 ,
where
Φ 1 = x 1 x ˙ 1 , Φ 2 = x 1 z 1 , Φ ˙ 1 = x ˙ 1 x ˙ 1 + x 1 x ¨ 1 , Φ ˙ 2 = x ˙ 1 z 1 + x 1 z ˙ 1 .
Step 2. The dynamics of z 2 are derived such as:
z ˙ 2 = x ˙ 2 α ˙ 1 , z ˙ 2 = x ˙ 2 + A B x ˙ 1 + C B w 1 + 1 B ξ ˙ 4 2 B Φ ˙ 1 + 2 c 1 B Φ ˙ 2 ,
where w 1 = u ˙ 1 . Then, to make z 2 0 to hold in finite time, the following candidate Lyapunov function is designed:
V 2 = V 1 + 1 2 z 2 2 ,
and its derivative:
V ˙ 2 = V ˙ 1 + z 2 z ˙ 2 , V ˙ 2 = c 1 z 1 2 + B 2 x 1 z 2 z 1 + z 2 z ˙ 2 , V ˙ 2 = c 1 z 1 2 + z 2 ( B 2 x 1 z 1 + z ˙ 2 ) .
To stabilize the energy in the system, B 2 x 1 z 1 + z ˙ 2 = σ 1 c 2 z 2 must be fulfilled, where σ 1 = c 2 z 2 + z ˙ 2 is the sliding surface, since z 1 0 holds in finite time. So,
V ˙ 2 = c 1 z 1 2 + z 2 σ 1 c 2 z 2 2 < 0 .
Step 3. The sliding surface,
σ 1 = c 2 z 2 + z ˙ 2 + B 2 Φ 3 ,
with dynamics
σ ˙ 1 = c 2 z ˙ 2 + z ¨ 2 + B 2 Φ ˙ 3 ,
where
z ¨ 2 = x ¨ 2 + A B x ¨ 1 + C B w ˙ 1 + 1 B ξ ¨ 4 2 B Φ ¨ 1 + 2 c 1 B Φ ¨ 2 , Φ 3 = z 1 x 1 , Φ ˙ 3 = z ˙ 1 x 1 z 1 x ˙ 1 x 1 2 , Φ ¨ 1 = x ¨ 1 x ˙ 1 + x ˙ 1 x ¨ 1 + x ˙ 1 x ¨ 1 + x 1 x 1 , Φ ¨ 2 = x ¨ 1 z 1 + x ˙ 1 z ˙ 1 + x ˙ 1 z ˙ 1 + x 1 z ¨ 1 ,
must be designed, such that V ˙ 2 decreases its energy in a finite time. Then, the following candidate Lyapunov function is introduced:
V 3 = V 2 + 1 2 σ 1 2 ,
with time derivative:
V ˙ 3 = V ˙ 2 + σ 1 σ ˙ 1 , V ˙ 3 = c 1 z 1 2 + z 2 σ 1 c 2 z 2 2 + σ 1 σ ˙ 1 , V ˙ 3 = c 1 z 1 2 c 2 z 2 2 + σ 1 ( z 2 + σ ˙ 1 ) .
Thus, to obtain
V ˙ 3 = c 1 z 1 2 c 2 z 2 2 c 3 σ 1 2 < 0 ,
or
V ˙ 3 = i = 1 2 c i z i 2 c 3 σ 1 2 < 0 ,
the following relationship must hold:
z 2 + σ ˙ 1 = c 3 σ 1 ,
where c 3 > 0 . Replacing Equations (31)–(33) in (38) yields:
z 2 + c 2 z ˙ 2 + x ¨ 2 + A B x ¨ 1 + C B w ˙ 1 + 1 B ξ ¨ 4 2 B Φ ¨ 1 + 2 c 1 B Φ ¨ 2 + B 2 Φ ˙ 3 = c 3 ( c 2 z 2 + z ˙ 2 ) c 3 B 2 Φ 3 .
Later, term w ˙ 1 is isolated:
w ˙ 1 = B C ( z 2 c 2 z ˙ 2 x ¨ 2 A B x ¨ 1 + 2 B Φ ¨ 1 2 c 1 B Φ ¨ 2 B 2 Φ ˙ 3 c 3 c 2 z 2 c 3 z ˙ 2 c 3 B 2 Φ 3 ) 1 C ξ ¨ 4 .
Then, by performing the double integration of term w ˙ 1 , the actual control u 1 is obtained:
u 1 = 0 t 0 t w ˙ 1 d τ d τ 1 C ξ 4 λ 1 | σ 1 | 1 / 2 s i g n ( σ 1 ) 0 t λ 2 s i g n ( σ 1 ) d τ .
The super-twisting terms λ 1 | σ 1 | 1 / 2 s i g n ( σ 1 ) and 0 t λ 2 s i g n ( σ 1 ) d τ , where λ 1 = 1.5 | Δ | 1 / 2 , λ 2 = 1.1 | Δ | > 0 , and | Δ | is an upper-bound of the expected disturbance to the system, are added to provide robust compensation for the disturbances to and uncertainties in the system. Replacing (40) in (41), the variable dependencies of the BSSMC control law are taken into account:
u 1 = 0 t 0 t ( B C ( z 2 c 2 z ˙ 2 x ¨ 2 A B x ¨ 1 + 2 B Φ 1 ¨ 2 c 1 B Φ ¨ 2 B 2 Φ ˙ 3 c 3 c 2 z 2 c 3 z ˙ 2 c 3 B 2 Φ 3 ) ) d τ d τ λ 1 | σ 1 | 1 / 2 s i g n ( σ 1 ) 0 t λ 2 s i g n ( σ 1 ) d τ .
It should be noted that the same procedure is carried out to obtain the control law u 2 . It is important to state, after obtaining the control law u 1 , that for the classical backstepping technique only two steps would be needed to obtain an adequate control law. However, by taking an additional step in the algorithm, the control law is integrated twice, giving additional integral action for the removal of the rate of change of the error variables at the steady state, as well as giving increased robustness against the disturbances and parameter uncertainties. Since the control law u 1 depends on unmeasurable derivatives x ˙ 1 , x ˙ 2 , x ¨ 1 , and x ¨ 2 , the state-of-the-art approach would force the implementation of differentiators, such as in [21,23], that would amplify the noise coming from the measured liquid levels y 1 and y 4 . For this, the technique of modulating functions is proposed to obtain the sliding surfaces σ i ( x ) that do not depend on the derivatives of the state variables and are noise-free.
The modulating functions (MF) technique performs integral transformations to compute the derivative-free algebraic relations between the inputs and outputs of a system, which allows the estimation of the internal states, unknown parameters, and fault detection with the advantage of not relying on the derivatives of the state variables and filtering noise while the modulation operation is carried out. To understand the mathematical background of the MF technique, the following definitions are presented.
Definition 3 
(Total Modulating Function [33]). Consider a sufficiently smooth function R × R R , with partial derivatives as
φ ( i ) ( t , t 1 ) : = i φ τ i ( τ , t 1 ) τ = t .
Then, function φ is called a modulating function of order k , if there exists t 0 < t 1 , such that
φ ( i ) ( t 0 , t 1 ) · φ ( i ) ( t 1 , t 1 ) = 0 , i = 0 , 1 , , n 1 .
A modulating function whose boundaries satisfy φ ( i ) ( t 0 , t 1 ) = φ ( i ) ( t 1 , t 1 ) = 0 is called a total modulating function.
Definition 4 
(Modulation Functional [34]). The modulation functional is defined as:
M [ h ] = h , φ Ω , I = Ω 0 t h ( x , τ + t T ) φ ( x , τ ) d τ d x ,
where h : Ω × R 0 + R and φ : Ω × I R represent the signal to be modulated and the modulating function, respectively, spatially defined on the n-dimensional rectangular region Ω : = x R n : 0 < x i < L i , i = 1 , 2 , , n and temporally on the moving time horizon I = [ t T , t ] of length T > 0 . Moreover, if the integration concerns only the spatial or temporal variable, the inner product notations h , φ Ω or h , φ I are used, respectively.
When implementing the MFs for filtering, an FIR filter with a modulating receding horizon can be realized, integrating only in the time dimension as follows:
M i [ h ] = ( 1 ) i t T t ϕ i ( t τ + T ) h ( τ ) d τ , M i [ h ] ( 1 ) i T s k = 0 N W k ϕ i ( k T s ) h ( ( l N + k ) T s ) ,
and we define it in its matrix form
M i [ h ] = K M F i h ( l N ) . . . h ( l 1 ) h ( l ) ,
where
K M F i = ( 1 ) i T s W 0 φ ( i ) ( 0 ) . . . W N 1 φ ( i ) ( ( N 1 ) T s ) W N φ ( i ) ( N T s ) T
is an i-dimension vector of modulating gains sampled with t = l T s , l N .
Lemma 1 
(Shift of derivatives property). Using integration by parts, a shift of derivatives can be attained as follows:
M 0 [ h ( n ) ] = t T t φ ( t τ + T ) h ( n ) ( τ ) d τ , M 0 [ h ( n ) ] = M n [ h ] ( 1 ) n t T t φ ( n ) ( t τ + T ) h ( τ ) d τ .
After presenting these definitions, the task is to modulate the error coordinates z 1 and z 2 , the virtual control law α 1 , and the sliding surface σ 1 , such that control u 1 is resilient to noise and does not depend on the derivatives of the system’s state variables. For this, the error coordinates z 1 ,
M 0 [ z 1 ] = M 0 [ x 1 x 1 ] M 0 [ z 1 ] = M 0 [ h 4 ] M 0 [ h 4 ] ,
and z 2 ,
M 0 [ z 2 ] = M 0 [ x 2 α 1 ] M 0 [ z 2 ] = M 0 [ y 1 ] + A B M 0 [ x 1 ] + C B M 0 [ u 1 ] 2 B M 0 [ Φ 1 ] + 2 c 1 B M 0 [ Φ 2 ] ,
are modulated over a prescribed modulation receding horizon. The nonlinearities h 1 , h 4 , Φ 1 , and Φ 2 in Equations (52) and (53) cannot be directly modulated. However, if h 1 and h 4 are replaced by measurements y 1 and y 4 , respectively, and h 4 and h ˙ 4 are known trajectories, then these nonlinear terms can be computed numerically and modulated on each iteration. Moreover, after modulating z 2 , the following modulated virtual control law α 1 is obtained:
M 0 [ α 1 ] = A B M 0 [ y 4 ] C B M 0 [ u 1 ] + 1 B M 0 [ Φ 1 ] 2 c 1 B M 0 [ Φ 2 ] ,
with modulated dynamics:
M 0 [ α ˙ 1 ] = M 0 [ A B x ˙ 1 C B u ˙ 1 + 2 B Φ ˙ 1 2 c 1 B Φ ˙ 2 ] , M 1 [ α 1 ] = A B M 1 [ y 4 ] C B M 1 [ u 1 ] + 2 B M 1 [ Φ 1 ] 2 c 1 B M 1 [ Φ 2 ] .
Finally, the sliding surface is modulated as follows:
M 0 [ σ 1 ] = M 0 [ c 2 z 2 + z ˙ 2 + B 2 Φ 3 ] , M 0 [ σ 1 ] = c 2 M 0 [ z 2 ] + M 1 [ z 2 ] + B 2 M 0 [ Φ 3 ] ,
with dynamics:
M 0 [ σ ˙ 1 ] = M 0 [ c 2 z ˙ 2 + z ¨ 2 + B 2 Φ ˙ 3 ] , M 1 [ σ 1 ] = c 2 M 1 [ z 2 ] + M 2 [ z 2 ] + B 2 M 1 [ Φ 3 ] .
With the modulating functions technique, all the virtual controls of u 1 have derivative-free input–output algebraic relations. Even though u 1 depends on the noise-measured state variables y 1 and y 4 , the MF technique filters noise when performing integration of these signals. Figure 3 shows the block diagram of the modulating-function-based backstepping super-twisting control system for a QTS.

5. Results and Discussion

5.1. Multivariable Sliding Mode Control of the Liquid Levels in Tanks n°3 and n°4

The control objective is to regulate the liquid levels of tank n°3 and n°4 at different operating points with a nonlinear controller. To this end, the performance of a first-degree sliding mode controller synthesized with the equivalent control method and a BSSTC were compared in the absence of measurement noise and disturbances, for a total simulation time T s i m = 1000 s and sample time T s = 0.01 s. The model parameters and liquid level operating points, the latter modeled as step inputs to the system, are shown in Table 1 and Table 2, respectively. The lower and upper saturation thresholds for the pumps control efforts were set at 10 % and 60 % , respectively. The design specifications of the BSSTC are shown in Table 3.
Figure 4 shows the control system time responses of all tanks at the operating points described in Table 2, where the liquid levels of tanks n°3 and n°4 were satisfactorily controlled by both the first-degree SMC and by the BSSTC. In addition, since the system valves were calibrated for the minimum phase response, the liquid levels in tanks n°1 and n°2 remained at safe operating points without overflow.
Figure 5 and Figure 6 show the control system time responses of tanks n°3 and n°4 in the different time intervals, where the SMC had a faster settling time than the BSSTC. The drawback of using the SMC for this specific application is seen in Figure 7, where the control efforts u 1 and u 2 of the SMC exhibited a considerable energy effort with undesirable chattering that oscillated from 10 % to 60 % of the pumping effort. This condition would be detrimental to the pumps actuators. On the other hand, the BSSTC did not present chattering in the absence of noise, and the energy effort was much lower than the SMC.
Figure 8 shows the errors between the actual and the desired liquid level over time, where it is noted that the error with the SMC converged faster than with the BSSTC. Moreover, Table 4 shows the results related to the following performance indices: the integral time absolute error (ITAE), integral absolute error (IAE), and the integral square error (ISE), where the SMC obtained better scores since it had a faster convergence, without overshoot, compared to the BSSTC, which showed a slower response with overshoot in some intervals of the trajectory. This result is especially important as it shows that although the SMC obtains slightly better dynamic error scores, it does so at the expense of more energy used in the control effort and the introduction of chattering.
Although, the BSSTC would work in conditions where the derivatives of the system state variables were accessible, and there was no measurement noise, this is not the case in most industrial applications, where sensors are not available due to budget constraints or they are not physically realizable. For example, for this plant, the velocity and acceleration of the liquid level could not be directly measured.

5.2. Modulating-Function-Based Backstepping Super-Twisting Control of Liquid Levels in Tanks n°3 and n°4

For this end, the goal was to test the robustness of the modulating-function-based BSSTC against nonlinear measurement noise. Noise with a mixture of Gaussian distributions was generated for the level transmitters L T 1 , L T 2 , L T 3 , and L T 4 , according to the mixture probability p ( x ) = ( 1 ϵ ) p g 1 + ϵ p g 2 , where ϵ is the probability for distribution p g 2 N ( 0 , ( σ 2 ) 2 ) , and σ 1 and σ 2 are the standard deviations. Table 5 shows the specifications of the Gaussian mixture model (GMM). In Figure 9 and Figure 10, the quantile–quantile plots (QQ plots) and sample vs. amplitude plots are shown, respectively, where it is evident that the GMM exhibited nonlinear characteristics.
In Figure 3, the modulating function filter was synthesized with a modulation function kernel order of one and two kernel derivatives, with a polynomial waveform type, a time horizon interval of 3.00 s, and a sample time of 0.01 s. The polynommial waveforms of the kernel and its two derivatives are shown in Figure 11. It received the measurement vector y and the control vector u as inputs, and sent the modulated vectors M i [ u ] and M i [ y ] to the backstepping sliding surface block. The backstepping sliding surface block generated the modulated sliding surfaces σ 1 and σ 2 . Finally, the backstepping super-twisting control yielded the modulated virtual control w i and output a modulated backstepping control signal along the super-twisting terms for the plant.
For this result, the performance of the MF-BSSTC in the presence of measurement noise was tested, setting up the gains according to Table 6.
In Figure 12 and Figure 13, the closed-loop time response of the QTS control system with the MF-BSSTC (green line) is plotted against the closed-loop time response of the system with measurement noise (red line). It is remarkable that the MF technique adequately filtered out the noise by modulating all the virtual controls of the backstepping algorithm, rendering it free from the derivatives of the system state variables. Nonetheless, the MF technique generated a delay of 1.5 s. The method to compensate for this delay, either a Smith predictor or any other delay-compensation technique, is outside the scope of this paper.
The controller gains were set to minimize the amplitude of the oscillations in the closed-loop response of the system. Since the MF-BSSTC was tested with non-Gaussian measurement noise, which is a mixture of two different Gaussian distributions, the MF-BSSTC effectively filtered both distributions, but a minimal amplitude variation over time was still visible when the Gaussian mixture changed from one distribution to another. Figure 14 shows the control efforts u 1 and u 2 , which exhibited minimum oscillations at steady state that were caused by the varying amplitude of the filtered noise. If the gains were increased to improve the settling time and robustness of the response, then the oscillations would increase, and with this, the control effort itself would show greater oscillations. Then, the MF-BSSTC must be synthesized based on the required closed-loop behavior and in terms of a relationship between the noise filtering capability and the desired robustness of the controller.
To test the robustness of the controller in the presence of measurement noise, a disturbance ξ 4 = 300 m 3 / s was applied to the system in the time interval t = [700 900] s, corresponding to a loss of flow in pump n°2. Then, the MF-BSSTC gains were increased to λ 1 = 29.24 and λ 2 = 418 , for | Δ | = 380 , so that the controller offered more robustness. Figure 15 shows that the controller was almost insensitive to the disturbance. Moreover, this robustness came with the chattering of the control signal u 2 at time interval t = [310 461] s, as shown in Figure 16. At the mentioned time interval, there was an overestimation of the controller gains, set as fixed for the worst-case disturbance. Perhaps, this high-gain condition of the controller could be solved by adapting the gains online.

5.3. Comparison with Other Controllers

The MF-BSSTC was compared to two of the most used controllers with the QTS, the decoupled PI controller [12] and the model predictive controller (MPC) [15].
A disturbance ξ 4 = 300 m 3 / s at time interval t = [600 1000] s was introduced to the system, corresponding to the flow loss in pump n°2. The MF-BSSTC gains for subcontrollers n°1 and n°2 were c 1 = 1.65 , c 2 = 1.05 , c 3 = 1.0 , λ 1 = 29.24 , and λ 2 = 418 and c 1 = 1.70 , c 2 = 1.05 , c 3 = 1.0 , λ 1 = 29.24 , and λ 2 = 418 , respectively. The decoupled PI controller was synthesized using the MATLAB PID tuner app, with gains K p 1 = 19.01 and k i 1 = 0.74 , and K p 1 = 26.91 and k i 1 = 1.18 , respectively, which maximized the robust transient behavior. On the other hand, the MPC was set with a prediction horizon H p = 200 and a control horizon H c = 200 . Moreover, both linear controllers were linearized with respect to the operating point ( h o 3 = 12.4 , h o 4 = 12.7 ).
Figure 17 shows that the MF-BSSTC had more robustness to disturbances and a faster response than the decoupled PI and MPC controllers. When disturbed, the PI controller and MPC showed strong set-point deviations from the set point, whereas the MF-BSSTC showed more insensitivity to the disturbance. As shown in Figure 18, the control efforts of the PI and MPC were lower than 10% for all the time responses, indicating that their robustness and speed of convergence were related to this low energy consumption.
Table 7 shows that the best results related to the ITAE, IAE, and ISE criteria were obtained with the MF-BSSTC. The PI and MPC controllers obtained high ITAE and ISE values for subcontrollers u 1 due to their underdamped responses to the disturbance. Even though they converged to the desired liquid level, their transient responses were very slow.
Regarding the mean computational cost, it was calculated based on the time it took to execute the loop controls. For instance, the mean computational efforts of the MF-BSSTC, PI, and MPC, for a total of 100 , 000 iterations, were 1.1685 × 10 4 s, 3.9600 × 10 4 s, and 2.8959 × 10 5 s, respectively. Even though the MF-BSSTC had memory requirements related to storing measured samples during the modulating receding horizon, its computational cost was lower than that of the PI controller. In addition, the MPC achieved a lower computational cost.

6. Conclusions

In this paper, it has been proved that it is feasible to adopt a nonlinear approach to the QTS liquid level control problem, using a multivariable sliding mode controller. For this, a first-degree SMC and a BSSTC were proposed based on the design of a multivariable sliding surface through the equivalent control method and backstepping virtual controls, respectively.
The results in the absence of noise showed that the liquid level control can be attained at any operating point without the need to linearize the system over a specific operating point. Although, the first-degree SMC exhibited a faster and more robust response, it introduced chattering into the control signal that was detrimental to the actuators of the plant.
On the other hand, the BSSTC showed almost no chattering at the expense of a slower settling time. In addition, a method to make the BSSTC derivative-free and noise-immune based on the modulation of the backstepping virtual controls was proposed, demonstrating that it is possible to use the backstepping algorithm and sliding modes in the presence of non-Gaussian noise through the modulating functions technique at the expense of introducing a time delay equal to half the receding horizon used.
Finally, with respect to future work, we will try to extend the current results considering the time delay in the dynamic behavior of the QTS and the conditions in which this system operates in non-minimum phase mode, as well as leakages and actuator failures.

Author Contributions

All the authors contributed to the development of the experiments, the result analysis, and the writing and review of the paper. Specifically, I.A.-C. and R.R.-P. were in charge of preparing the introduction, G.P.-Z. and J.S.-M. of the modeling of the QTS system, I.A.-C. of the design and the simulation of the QTS control system, and I.A.-C., R.R.-P. and J.S.-M. of the overall ideas of the exposed research and the general conception of the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Dirección de Fomento de la Investigación at the Pontificia Universidad Católica del Perú (PUCP) through grant DFI-PI0764.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Albertos, P.; Salas, A. Multivariable Control Systems. An Engineering Approach; Springer: Berlin, Germany, 2004. [Google Scholar]
  2. Che, Y.; Zhao, Z.; Wang, Z.; Liu, F. Iterative learning model predictive control for multivariable nonlinear batch processes based on dynamic fuzzy PLS model. J. Process Control 2022, 119, 1–12. [Google Scholar] [CrossRef]
  3. Sotomayor-Moriano, J.; Pérez-Zúñiga, G.; Soto, M. A virtual laboratory environment for control design of a multivariable process. IFAC-PapersOnLine 2019, 52, 15–20. [Google Scholar] [CrossRef]
  4. Altabey, W.A. Model optimal control of the four tank system. Int. J. Syst. Sci. Appl. Math. 2016, 8, 30–41. [Google Scholar]
  5. Osman, A.; Kara, T.; Arıcı, M. Robust adaptive control of a quadruple tank process with sliding mode and pole placement control strategies. IETE J. Res. 2021, 1–14. [Google Scholar] [CrossRef]
  6. Rivas-Perez, R.; Sotomayor-Moriano, J.; Pérez-Zuñiga, G.; Soto-Angles, M.E. Real-time implementation of an expert model predictive controller in a pilot-scale reverse osmosis plant for brackish and seawater desalination. Appl. Sci. 2019, 9, 2932. [Google Scholar] [CrossRef]
  7. Shah, D.H.; Patel, D.M. Design of sliding mode control for quadruple-tank MIMO process with time delay compensation. J. Process Control 2019, 76, 46–61. [Google Scholar] [CrossRef]
  8. Rivas-Perez, R.; Sotomayor-Moriano, J.; Pérez-Zúñiga, C.G.; Calderon-Mendoza, E.M. Design of a multivariable GPC based on an industrial PC for control of a reverse osmosis unit of a pharmaceutical industry. Rev. Mex. Ing. Quim. 2016, 15, 259–273. [Google Scholar]
  9. Aravindan, M.; Chilambuchelvan, A.; Tamilselvi, S. A supervised model of multivariable control in quadruple tank system. Appl. Artif. Intell. 2023, 37, 606–622. [Google Scholar]
  10. Johansson, K.H. The quadruple-tank process: A multivariable laboratory process with an adjustable zero. IEEE Trans. Control Syst. Technol. 2000, 8, 456–465. [Google Scholar] [CrossRef]
  11. Garrido, J.; Vazquez, F.; Morilla, F. Centralized multivariable control by simplified decoupling. J. Process Control 2012, 22, 1044–1062. [Google Scholar] [CrossRef]
  12. Cordova, J.; Sotomayor, J.; Perez-Zuniga, G. Actuator fault tolerant PI-control design on TITO process. In Proceedings of the XVIII Latin-American Congress of Automatic Control, Quito, Ecuador, 24–26 October 2018; pp. 1–6. [Google Scholar]
  13. Bequette, B. Process Control: Modeling, Design, and Simulation; Prentice Hall: Upper Saddle River, NJ, USA, 2010. [Google Scholar]
  14. Vadigepalli, R.; Gatzke, E.; Doyle, F. Robust control of a multivariable experimental four-tank system. Ind. Eng. Chem. Res. 2001, 40, 1916–1927. [Google Scholar] [CrossRef]
  15. Sanchez-Zurita, V.; Perez-Zuñiga, G.; Sotomayor-Moriano, J. Reconfigurable model predictive control applied to the quadruple tank process. In Proceedings of the 15th European Workshop on Advanced Control and Diagnosis, Bologna, Italy, 21–22 November 2019; pp. 937–957. [Google Scholar]
  16. Andronov, A.; Chaikin, C. Theory of Oscillations, 1st ed.; Princeton University Press: Princeton, NJ, USA, 1949. [Google Scholar]
  17. Utkin, V. Sliding Modes and Their Applications in Variable Structure Systems; MIR: Moscow, Russia, 1978. [Google Scholar]
  18. Shtessel, Y.; Levant, A.; Fridman, L.; Edwards, C. Sliding Mode Control and Observation; Springer: New York, NY, USA, 2014. [Google Scholar]
  19. Lee, H.; Utkin, V.I. Chattering suppression methods in sliding mode control systems. Annu. Rev. Control 2007, 31, 179–188. [Google Scholar] [CrossRef]
  20. Levant, A. Robust exact differentiation via sliding mode technique. Automatica 1998, 34, 379–384. [Google Scholar] [CrossRef]
  21. Aranda, I.A.; Pérez-Zúñiga, G. Highly maneuverable target tracking under glint noise via uniform robust exact filtering differentiator with intrapulse median filter. IEEE Trans. Aerosp. Electron. Syst. 2022, 58, 2541–2559. [Google Scholar] [CrossRef]
  22. Biswas, P.; Srivastava, R.; Ray, S.; Samanta, A. Sliding mode control of quadruple tank process. Mechatronics 2009, 19, 548–561. [Google Scholar] [CrossRef]
  23. Chaudhari, V.; Tamhane, B.; Kurode, S. Robust liquid level control of quadruple tank system—Second order sliding mode approach. IFAC-PapersOnLine 2020, 53, 7–12. [Google Scholar] [CrossRef]
  24. Han, J. From PID to active disturbance rejection control. IEEE Trans. Ind. Electron. 2009, 56, 900–906. [Google Scholar] [CrossRef]
  25. Hashim, Z.S.; Khani, H.I.; Azar, A.T.; Khan, Z.I.; Smait, D.A. Robust liquid level control of quadruple tank system: A nonlinear model-free approach. Actuators 2023, 12, 119. [Google Scholar] [CrossRef]
  26. Chalanga, A.; Kamal, S.; Fridman, L.; Bandyopadhyay, B.; Moreno, J. Implementation of super-twisting control: Super-twisting and higher order sliding-mode observer-based approaches. IEEE Trans. Ind. Electron. 2016, 63, 3677–3685. [Google Scholar] [CrossRef]
  27. Sira-Ramirez, H. Sliding Mode Control: The Delta-Sigma Modulation Approach; Birkhauser Verlag AG: Basel, Switzerland, 2016. [Google Scholar]
  28. Kanellakopoulos, I.; Kokotovic, P.; Morse, A. Systematic design of adaptive controllers for feedback linearizable systems. IEEE Trans. Automat. Control 1991, 36, 1241–1253. [Google Scholar] [CrossRef]
  29. Zhou, J.; Wen, C. Adaptive Backstepping Control of Uncertain Systems; Springer: Berlin, Germany, 2008. [Google Scholar]
  30. Sira-Ramirez, H.; Llanes-Santiago, O. Adaptive dynamical sliding mode control via backstepping. In Proceedings of the 32nd IEEE Conference on Decision and Control, San Antonio, TX, USA, 15–17 December 1993; pp. 1422–1427. [Google Scholar]
  31. Turki, A.; Said, H.; M’Sahli, F. Backstepping control for a quadruple tank process based on adaptive observer. In Proceedings of the 12th IEEE International Multi-Conference on Systems, Signals and Devices, Mahdia, Tunisia, 16–19 March 2015; pp. 1–5. [Google Scholar]
  32. Gouta, H.; Said, H.; M’Sahli, F. Observer-based backstepping liquid level controller for a quadruple tank process. In Proceedings of the 16th IEEE Conference on Sciences and Techniques of Automatic Control and Computer Engineering, Monastir, Tunisia, 21–23 December 2015; pp. 351–356. [Google Scholar]
  33. Jouffroy, J.; Reger, J. Finite-time simultaneous parameter and state estimation using modulating functions. In Proceedings of the IEEE Conference on Control Applications, Sydney, NSW, Australia, 21–23 September 2015; pp. 394–399. [Google Scholar]
  34. Ghaffour, L.; Noack, M.; Reger, J.; Laleg-Kirati, T.M. Non-asymptotic state estimation of linear reaction diffusion equation using modulating functions. IFAC-PapersOnLine 2020, 53, 4196–4201. [Google Scholar] [CrossRef]
Figure 1. Quadruple tank system of the PUCP’s Advanced Control Systems Laboratory.
Figure 1. Quadruple tank system of the PUCP’s Advanced Control Systems Laboratory.
Sensors 23 05222 g001
Figure 2. Quadruple tank system diagram.
Figure 2. Quadruple tank system diagram.
Sensors 23 05222 g002
Figure 3. Block diagram of the modulating-function-based backstepping super-twisting control system for a QTS.
Figure 3. Block diagram of the modulating-function-based backstepping super-twisting control system for a QTS.
Sensors 23 05222 g003
Figure 4. Control system time responses of all tanks: (a) tank n°1; (b) tank n°2; (c) tank n°3; (d) tank n°4; QTS with SMC (red line); QTS with BSSTC (green line); set point (blue line).
Figure 4. Control system time responses of all tanks: (a) tank n°1; (b) tank n°2; (c) tank n°3; (d) tank n°4; QTS with SMC (red line); QTS with BSSTC (green line); set point (blue line).
Sensors 23 05222 g004
Figure 5. Control system time responses of tank n°3: (a) t = [0 1000] s; (b) t = [0 300] s; (c) t = [300 600] s; (d) t = [600 1000] s; QTS with SMC (red line); QTS with BSSTC (green line); set point (blue line).
Figure 5. Control system time responses of tank n°3: (a) t = [0 1000] s; (b) t = [0 300] s; (c) t = [300 600] s; (d) t = [600 1000] s; QTS with SMC (red line); QTS with BSSTC (green line); set point (blue line).
Sensors 23 05222 g005
Figure 6. Control system time responses of tank n°4: (a) t = [0 1000] s; (b) t = [0 300] s; (c) t = [300 600] s; (d) t = [600 1000] s; QTS with SMC (red line); QTS with BSSTC (green line); set point (blue line).
Figure 6. Control system time responses of tank n°4: (a) t = [0 1000] s; (b) t = [0 300] s; (c) t = [300 600] s; (d) t = [600 1000] s; QTS with SMC (red line); QTS with BSSTC (green line); set point (blue line).
Sensors 23 05222 g006
Figure 7. Control efforts u 1 and u 2 : (a) control effort u 1 ; (b) control effort u 2 ; QTS with SMC (red line); QTS with BSSTC (green line).
Figure 7. Control efforts u 1 and u 2 : (a) control effort u 1 ; (b) control effort u 2 ; QTS with SMC (red line); QTS with BSSTC (green line).
Sensors 23 05222 g007
Figure 8. Control system errors: (a) QTS with BSSTC; (b) QTS with SMC; error h 3 - h 3 (green dotted line); error h 4 - h 4 (blue line).
Figure 8. Control system errors: (a) QTS with BSSTC; (b) QTS with SMC; error h 3 - h 3 (green dotted line); error h 4 - h 4 (blue line).
Sensors 23 05222 g008
Figure 9. Gaussian mixture quantile–quantile plot: (a) level transmitter L T 1 ; (b) level transmitter L T 2 ; (c) level transmitter L T 3 ; (d) level transmitter L T 4 .
Figure 9. Gaussian mixture quantile–quantile plot: (a) level transmitter L T 1 ; (b) level transmitter L T 2 ; (c) level transmitter L T 3 ; (d) level transmitter L T 4 .
Sensors 23 05222 g009
Figure 10. Noise amplitude: (a) level transmitter L T 1 ; (b) level transmitter L T 2 ; (c) level transmitter L T 3 ; (d) level transmitter L T 4 .
Figure 10. Noise amplitude: (a) level transmitter L T 1 ; (b) level transmitter L T 2 ; (c) level transmitter L T 3 ; (d) level transmitter L T 4 .
Sensors 23 05222 g010
Figure 11. Modulating functions: modulating function kernel (green line); modulating function kernel first derivative (blue line); modulating function kernel second derivative (red line).
Figure 11. Modulating functions: modulating function kernel (green line); modulating function kernel first derivative (blue line); modulating function kernel second derivative (red line).
Sensors 23 05222 g011
Figure 12. Control system time responses of tank n°3: (a) t = [0 1000] s; (b) t = [0 300] s; (c) t = [300 600] s; (d) t = [600 1000] s; control system time responses with measurement noise (red line); control system time responses with MF-BSSTC (green line); set point (blue line).
Figure 12. Control system time responses of tank n°3: (a) t = [0 1000] s; (b) t = [0 300] s; (c) t = [300 600] s; (d) t = [600 1000] s; control system time responses with measurement noise (red line); control system time responses with MF-BSSTC (green line); set point (blue line).
Sensors 23 05222 g012
Figure 13. Control system time responses of tank n°4: (a) t = [0 1000] s; (b) t = [0 300] s; (c) t = [300 600] s; (d) t = [600 1000] s; control system time responses with measurement noise (red line); control system time responses with MF-BSSTC (green line); set point (blue line).
Figure 13. Control system time responses of tank n°4: (a) t = [0 1000] s; (b) t = [0 300] s; (c) t = [300 600] s; (d) t = [600 1000] s; control system time responses with measurement noise (red line); control system time responses with MF-BSSTC (green line); set point (blue line).
Sensors 23 05222 g013
Figure 14. Control efforts u 1 and u 2 : (a) control effort u 1 ; (b) control effort u 2 ; MF-BSSTC (green line).
Figure 14. Control efforts u 1 and u 2 : (a) control effort u 1 ; (b) control effort u 2 ; MF-BSSTC (green line).
Sensors 23 05222 g014
Figure 15. Control system time responses of tank n°3 in the presence of a disturbance to the system: (a) t = [0 1000] s; (b) t = [0 300] s; (c) t = [300 600] s; (d) t = [600 1000] s; control system time responses with measurement noise (red line); control system time responses with MF-BSSTC (green line); set point (blue line).
Figure 15. Control system time responses of tank n°3 in the presence of a disturbance to the system: (a) t = [0 1000] s; (b) t = [0 300] s; (c) t = [300 600] s; (d) t = [600 1000] s; control system time responses with measurement noise (red line); control system time responses with MF-BSSTC (green line); set point (blue line).
Sensors 23 05222 g015
Figure 16. Control efforts u 1 and u 2 in the presence of a disturbance to the system: (a) control effort u 1 ; (b) control effort u 2 ; MF-BSSTC (green line).
Figure 16. Control efforts u 1 and u 2 in the presence of a disturbance to the system: (a) control effort u 1 ; (b) control effort u 2 ; MF-BSSTC (green line).
Sensors 23 05222 g016
Figure 17. Control system time responses of tank n°3 in the presence of a disturbance to the system: (a) t = [0 1000] s; (b) t = [0 300] s; (c) t = [300 600] s; (d) t = [600 1000] s; control system time responses with measurement noise (red line); control system time responses with MF-BSSTC (green line); control system time responses with decoupled PI (cyan line); control system time responses with MPC (magenta line); set point (blue line).
Figure 17. Control system time responses of tank n°3 in the presence of a disturbance to the system: (a) t = [0 1000] s; (b) t = [0 300] s; (c) t = [300 600] s; (d) t = [600 1000] s; control system time responses with measurement noise (red line); control system time responses with MF-BSSTC (green line); control system time responses with decoupled PI (cyan line); control system time responses with MPC (magenta line); set point (blue line).
Sensors 23 05222 g017
Figure 18. Control efforts u 1 and u 2 in the presence of a disturbance to the system: (a) control effort u 1 ; (b) control effort u 2 ; MF-BSSTC (green line); decoupled PI (cyan line); MPC (magenta line).
Figure 18. Control efforts u 1 and u 2 in the presence of a disturbance to the system: (a) control effort u 1 ; (b) control effort u 2 ; MF-BSSTC (green line); decoupled PI (cyan line); MPC (magenta line).
Sensors 23 05222 g018
Table 1. QTS parameters.
Table 1. QTS parameters.
ParameterUnitsValue
A 1 , A 2 , A 3 , A 4 cm 2 706.85
a 1 , a 2 , a 3 , a 4 cm 2 1.89, 1.89, 5.39, 5.39
gcm/s 2 981
γ 1 , γ 2 unitless0.80, 0.90
k 1 , k 2 unitless1.00
u 1 m a x , u 1 m a x cm 3 /s2863
d 1 , d 2 , d 3 , d 4 cm30
h 1 m a x , h 2 m a x , h 3 m a x , h 4 m a x cm45
Table 2. Liquid level operating points.
Table 2. Liquid level operating points.
Operating PointTime SpanUnitsValue
h 3 , h 4 T s i m = 0 300 scm40, 7
h 3 , h 4 T s i m = 300 600 scm10, 34
h 3 , h 4 T s i m = 600 1000 scm24, 15
Table 3. The design specifications of the BSSTC for simulation in the absence of measurement noise.
Table 3. The design specifications of the BSSTC for simulation in the absence of measurement noise.
SubcontrollerParameterValue (Dimensionless)
h 3 c 1 , c 2 , c 3 , | Δ | , λ 1 , λ 2 1.50, 1.20, 1.00, 1.00, 1.50, 1.10
h 4 c 1 , c 2 , c 3 , | Δ | , λ 1 , λ 2 2.50, 1.50, 1.00, 1.00, 1.50, 1.10
Table 4. Performance indices.
Table 4. Performance indices.
ControllerSubcontrollerITAE (cm)IAE (cm)ISE (cm 2 )
SMC h 3 1813.81267.92525.7
h 4 2864.3782.91028.9
BSSTC h 3 1919.51346.02549.2
h 4 3203.0861.11050.3
Table 5. Gaussian mixture model.
Table 5. Gaussian mixture model.
Distr.ParametersUnitsValues
p g 1 Mean μ h 1 , μ h 2 , μ h 3 , μ h 4 cm0.01, 0.02, 0.03, 0.02
Standard deviation σ h 1 , σ h 2 , σ h 3 , σ h 4 cm0.15
Probability (1 −  ϵ )unitless0.95
p g 2 Mean μ h 1 , μ h 2 , μ h 3 , μ h 4 cm0.05, 0.04, 0.07, 0.04
Standard deviation σ h 1 , σ h 2 , σ h 3 , σ h 4 cm0.25
Probability ( ϵ )unitless0.05
Table 6. MF-BSSTC design specifications for the simulation under measurement noise.
Table 6. MF-BSSTC design specifications for the simulation under measurement noise.
SubcontrollerParameterValue (Dimensionless)
h 3 c 1 , c 2 , c 3 , | Δ | , λ 1 , λ 2 1.20, 1.02, 1.00, 1.00, 1.50, 1.10
h 4 c 1 , c 2 , c 3 , | Δ | , λ 1 , λ 2 1.20, 1.05, 1.00, 1.00, 1.50, 1.10
Table 7. Performance indices.
Table 7. Performance indices.
ControllerControl InputITAEIAEISE
MF-BSSTC u 1 4711.877.1659.0
u 2 4176.1101.4786.7
PI u 1 33,243.0911.94564.3
u 2 5549.4387.21935.1
MPC u 1 9997.1250.91414.0
u 2 2804.9243.61176.9
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Aranda-Cetraro, I.; Pérez-Zúñiga, G.; Rivas-Pérez, R.; Sotomayor-Moriano, J. Nonlinear Robust Control by a Modulating-Function-Based Backstepping Super-Twisting Controller for a Quadruple Tank System. Sensors 2023, 23, 5222. https://doi.org/10.3390/s23115222

AMA Style

Aranda-Cetraro I, Pérez-Zúñiga G, Rivas-Pérez R, Sotomayor-Moriano J. Nonlinear Robust Control by a Modulating-Function-Based Backstepping Super-Twisting Controller for a Quadruple Tank System. Sensors. 2023; 23(11):5222. https://doi.org/10.3390/s23115222

Chicago/Turabian Style

Aranda-Cetraro, Italo, Gustavo Pérez-Zúñiga, Raul Rivas-Pérez, and Javier Sotomayor-Moriano. 2023. "Nonlinear Robust Control by a Modulating-Function-Based Backstepping Super-Twisting Controller for a Quadruple Tank System" Sensors 23, no. 11: 5222. https://doi.org/10.3390/s23115222

APA Style

Aranda-Cetraro, I., Pérez-Zúñiga, G., Rivas-Pérez, R., & Sotomayor-Moriano, J. (2023). Nonlinear Robust Control by a Modulating-Function-Based Backstepping Super-Twisting Controller for a Quadruple Tank System. Sensors, 23(11), 5222. https://doi.org/10.3390/s23115222

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop