1. Introduction
The present paper introduces a new design requirement in the family of noninteracting control problems. The design requirement is that of a common noninteracting control with simultaneous common partial output zeroing for multi-model linear time-invariant systems. The special case of the present design requirement, where the number of system inputs is greater than or equal to the number of system outputs, was presented in [
1].
The present design requirement belongs to the class of common design problems. The common control design is an indispensable step that should be checked before proceeding to switching. The common control design has a long history (see [
2,
3,
4,
5,
6,
7,
8,
9,
10,
11]) and several applications; indicatively, see [
9,
10,
11], where the case of robot manipulators was studied. In the framework of the present common design requirement, studied for the class of multi-model normal linear time-invariant systems, the set of the outputs of each model of the multi-model system is divided into two subsets. Each output belonging to the first subset is required to be controlled independently by a respective external command. The rest outputs, namely the outputs belonging to the second subset, are required to be equal to zero, independently from the choice of the external commands. It is important to point out that the partition of the set of the system outputs into two subsets depends upon the respective model of the multi-model description. So, the two subsets may be different for different models of the multi-model system description. After an appropriate row rearrangement, the transfer matrices of the closed-loop systems, resulting after the application of the common controller, are expressed in four block forms. One block form is a rational diagonal matrix, and the remaining three are zero blocks. The main motivation for the derivation of the results of the present paper was the observation that, in several MIMO system applications, the deviations of some output performance variables from their nominal values are required to be controlled independently, while the other performance variables are required to remain at their nominal value.
Here, the design requirement of common noninteracting control with simultaneous common partial output zeroing is solved using a regular static measurement output controller. The solvability conditions of this design problem are established, and the general solution of the common controllers, solving the problem, is determined. Additionally, the design requirement of the approximate asymptotic command following is studied by the appropriate choice of the free parameters of the precompensator matrix. The remaining free parameters of the common controller are derived using a metaheuristic algorithm (see [
12,
13]) that is executed towards the minimization of a 2-norm cost, evaluating the performance of the transient responses of the closed-loop systems. Finally, the results of the proposed common design scheme are illustrated through a numerical example, where a nonlinear process in two different operating points is controlled. It is mentioned that the pure (noncommon) output zeroing problem, presented in [
14], requires the deviations of the system outputs from their nominal values to be equal to zero, i.e., the original performance variables are required to remain at their desired operating values.
Before closing the presentation of the design results of the present paper, it is important to mention that the design requirement of a common noninteracting control with simultaneous common partial output zeroing was studied in [
1] for the special case where the number of inputs is greater than or equal to the number of outputs. Finally, it is mentioned that in [
1], a rather abstractive presentation of the design procedure and the proofs of the results was presented. Here, the complete design procedure, using all the required formalism, is presented.
Another motivation for the derivation of the present design results is the open control problem of the two-model robot-tracked vehicle description presented in [
15,
16]. Tracked Unmanned Ground Vehicles (UGVs) are important for open land applications, where conventional robotic vehicles or two-legged robots are not appropriate. Indicative applications (see [
15,
16,
17,
18,
19,
20,
21,
22]) are met in agriculture, surveillance, patrols, the transfer of goods, emergency situations, the military, underwater operations, and constructions. It is important to mention that the control problem of the two-model robot-tracked vehicle description in [
15,
16] has not yet been studied.
The precision and accuracy of the operations of Tracked UGVs depend greatly upon control efficiency. Several approaches and models were proposed (indicatively, see [
15,
16,
22,
23,
24,
25,
26,
27,
28,
29,
30,
31,
32,
33]) to handle the issue. In [
15,
16], a two-model variable switching description of a robot-tracked vehicle was proposed. The models include the simplified dynamics of the rigid body and the motors driving the vehicle. The switching condition between the models is in the form of a sharp constraint upon the forces that are generated to be the left-side and right-side motors. In [
22], the mathematical description of an unmanned tracked excavator vehicle that is equipped with an excavation arm was developed. In [
23], a switching dynamic description of a tracked vehicle was presented, and a speed-regulating control strategy and a torque-regulating control strategy were developed. In [
24], the mathematical description of a tracked vehicle was presented in modular form, and a Stanley-type controller was derived. In [
25], the mathematical description model of a 3DOF-tracked mobile robot, including the dynamics and kinematics of the skid-steering mechanism, was presented, and a PID controller was applied to control the speed of the vehicle. In [
26], based on the kinematics of a robot-tracked vehicle, a control law for trajectory tracking under slip conditions, system constraints, and varying dynamics were presented. In [
27], the problem of localization and trajectory tracking control under slip conditions was investigated based on the kinematics of a robot-tracked vehicle. The approach was based on feedback linearization and indirect Kalman filtering. In [
28], the dynamic system description of a tracked vehicle undergoing skid-steering on horizontal hard terrain was presented. The vehicle’s path was controlled via a modified PID computed torque control scheme. In [
29,
30], the problem of trajectory tracking control of an autonomous tracked vehicle was investigated through a backstepping kinematic controller and an integral sliding mode controller based on vehicle dynamics. In [
31], the problem of position control for a tracked vehicle was studied using driving force control and virtual-turning velocity control. In [
32], the problem of the steering control of tracked vehicles was studied through differential steering control rules and a Fuzzy PID control scheme. In [
33], three control strategies were proposed for the investigation of the closed-loop dynamical performance of tracked vehicles using different desired trajectories. The three control schemes were a robust nonlinear controller, a speed compensation-based fuzzy logic controller, and a speed compensation-based proportional–integral controller.
In the present paper, based on the two-model description of the Tracked UGV presented in [
15,
16] and the herein-derived theoretical results, the independent forward and angular motion of the robot vehicle is achieved for the case where the performance outputs are the linear velocity and the heading angle. Here, the measurement variables are the two motor currents and the heading angle of the vehicle. These measurements are clearly offered for acquisition by a small-scale embedded controller. The problem of common noninteracting control with partial output zeroing is proven to be solvable, and the controller matrices’ general solution is derived. This general solution achieves approximate command following and I/O stability using the previously mentioned metaheuristic algorithm. The present application was studied in [
1]. Furthermore, in the case where the four mechanical resistances of the vehicle, as well as the lateral friction of the tracks, are uncertain, it is shown through extensive simulation experiments that the efficiency of the derived common controller remains satisfactory. This result provides promising perspectives for physical experimentation using the derived controller. Finally, the performance of the derived common controller for the Tracked UGV is compared to a switching controller designed to satisfy noninteracting control with simultaneous partial output zeroing for each model of the Tracked UGV. The superiority of the performance of the common controller is illustrated through the simulation results for a composite maneuver.
3. Common Noninteracting Control with Simultaneous Common Partial Output Zeroing for Multi-Model Descriptions
3.1. Definition of the Common Design Requirement
Consider the multi-model linear time-invariant system description (see also [
1,
9,
10])
where
is the state vector,
is the control input vector,
is the performance output vector, and
is the measurement output vector. The parameters
,
, and
were already presented in the previous section. Moreover, similarly to the previous section and without a loss of generality, the following assumptions are considered:
In the present multi-model case, the outputs are also grouped into two sets, depending upon the index of the respective model. Let
be the predefined set of outputs. The set of the rest of the outputs is defined to be
. Let
and
. Clearly, it holds that
Let
be a row rearrangement matrix, setting first the elements of
belonging to
. Thus, the elements of
are rearranged, and so the rearranged vector has the form
and the rows of the output matrix are rearranged accordingly, i.e.,
Here, the controller in (2), being independent of the number of the system model, is also considered. Following the second case, presented in the previous section, it is assumed that the following equalities hold true:
where (11) was used. Regarding
, it is also mentioned that it may have any relative value with respect to
. Additionally, the controller in (2) is considered to be regular, i.e., the precompensator matrix, being square (see (15)), is invertible, i.e.,
Using the definitions (13) and (14) as well as the assumptions in (15) and (16), the design requirement of common noninteracting control with simultaneous common partial output zeroing is expressed as follows: find the appropriate
,
, and
such that
where
Regarding (15), it is observed that it can be derived from (18) and (19) as a necessity. Moreover, it is observed that (15) can be covered by (17). Hence, the only assumption considered here is (17), in the sense that all the system assumptions are covered by the definition of the problem in (18) and (19) and the regularity of the controller (17). However, the inequalities in (16) are quite useful to exclude system cases before testing the solvability conditions that will be presented in
Section 3.3.
3.2. Preliminary Results
In this subsection, a set of necessary conditions, as well as a transformation of the controller matrices, are presented.
Lemma 1. For common noninteracting control with simultaneous common partial output zeroing to be satisfied, it is necessary for the following conditions to be satisfied:where denotes the rank of the argument matrix over the field rational functions of with real coefficients. Remark 1. The first condition in Lemma 1 is the right invertibility of the subsystem of system (1), including only the independently controlled outputs. The second condition in Lemma 1 is a structural dependence relation of the rest of the outputs of (1). □
From (17), the following transformation of the controller matrices and the diagonal elements of the closed-loop system is proposed:
If (21) is satisfied, then the problem is reduced to that of solving only (53) (see
Appendix A) under the constraints (17) and (19). Thus, using the definitions in (22), the problem is reduced to that of solving the following set of linear equations with respect to the unknowns
,
, and
,
where
is the
-th row of
, and where the unknowns are constrained to satisfy the inequalities
Following the approach in [
10,
11], it can be proven that (24)–(26) is satisfied if, and only if, the following set of equations is solvable with respect to
,
,
and
, under the constraints (25) and (26),
where
and
denote the left inverse and the left orthogonal of the argument matrix, while
and
denote the right inverse and the right orthogonal of the argument matrix. As is known, there are several methods to compute the orthogonal (right or left) and the inverse (right or left) of the argument matrix. Moreover, the determination of these matrices is not unique. However, their participation in the solution of a linear nonhomogeneous vector equation does not modify the general class of the solutions of the equation.
3.3. Necessary and Sufficient Conditions
In order to present the solvability conditions of the problem at hand, a set of recursive definitions related to Markov parameter vectors of the system (10a) and (10b) after the application of the row rearrangement matrix
will be presented. These definitions provide a left transformation of the rows of the output matrix. Particularly for
and
, the following recursive definitions, being modifications of the respective definitions in [
10], are presented:
Clearly, the existence of an appropriate
, such that
, is a consequence of the first condition in Lemma 1. The nonnegative integers
are the essential orders corresponding to the decoupled parts of the closed-loop transfer matrices. Moreover, from the definition of
, the following equalities are derived:
Finally, the following set is introduced:
where
The above positive integer denotes the maximum number of models requiring independent control of the
-
th output. From (30), (32), and (16), it is observed that
The set
can be expressed, in terms of its elements, as follows:
In the above expression, the elements of are considered to be distinct and ordered by their magnitude, i.e., , where and .
Consider the following block matrix structure:
where
,
, and
is a nonnegative integer. The quantities
and
are positive integers. Using (34), the following matrix, which will be used in the proof of Lemma 2 and the presentation of the solvability conditions, is defined:
The following lemma reveals the necessary dependence relations among the Markov parameter row vectors of the multi-model description.
Lemma 2. If common noninteracting control with simultaneous common partial output zeroing is satisfied, then the following conditions must be satisfied: A set of analytic formulas of the respective dependence coefficients is determined in the following corollary using condition (37).
Corollary 1. If (37) is satisfied, then for everyand everythere exist real nonzero real coefficients, denoted by, satisfying the following dependence relationThe coefficients are uniquely determined by the relations Before presenting the necessary and sufficient conditions for the solvability of the problem at hand, the definition of a vector and a matrix, using data from all models in (10a) and (10b), will be presented. First, the following two vectors are defined:
The vector and the matrix, using data from all models, are defined as follows:
We are now in a position to present the main result of the paper.
Theorem 1. The common noninteracting control problem with simultaneous common partial output zeroing is solvable via a regular measurement output controller of the form (2), if, and only if, the conditions (20), (21), (36), (37), and the following conditions, are satisfied: 3.4. General Solution of the Controller
The following definitions are introduced to present the general solution of the controller matrices solving the problem at hand:
Theorem 2. If the conditions in Theorem 1 are satisfied, then the general solution of the controller matrices, solving the problem of common noninteracting control with simultaneous common partial output zeroing, is The parameters, where, are arbitrary. The vectors arearbitrary row vectors where. Moreover, the following matrices are arbitrary: The second arbitrary matrix is constrained to satisfy the inequality: Remark 2. The results of Theorems 1 and 2 cover two important special cases. The first is the case of the common design of the problem at hand via regular static state feedback, i.e., . The second is the case of single model descriptions, i.e., the case where . □
3.5. Numerical Example
To illustrate the proposed controller design procedure, consider the following nonlinear system
where
,
,
, and
are the state, input, performance output, and measurable output vectors, respectively. Furthermore, assume that the matrix
is equal to the 4-by-4 identity matrix as well as that of
and
where
Consider the following set of operating points for the actuatable inputs and state variables corresponding to two different modes of operation:
Mode 1: , , , , ,
Mode 2: , , , , ,
The linear approximants of the nonlinear model are of the form (10), and the respective system matrices are computed to e:
Note that the system matrices
,
,
, and
are identical to the respective matrices of the original nonlinear model. Let
From the form of the linear approximant system matrices as well as the above d above definitions, it can be observed that
,
,
,
,
,
,
,
, and
. Using the above data, the following matrices are computed:
Using the system matrices of the linear approximants, the following matrices are computed to be
The following quantities are computed to be
,
, and
. Furthermore, the following vectors are computed to be
From the above matrices, it is observed that the rank conditions in (36) and (37) are satisfied since computed to be
,
,
, and
. Moreover, the following sets are derived:
,
. Finally, the following matrices are computed:
Clearly, the condition in (40) is satisfied. Thus, the necessary and sufficient conditions for the solvability of the common noninteracting control problem with simultaneous common partial output zeroing are satisfied.
It can readily be verified that
and that the general form of the controller matrices solving the problem are
where
,
,
,
, and
are free parameters, with
and
. It is important to mention that the method for the computation of the general solution of the common controllers in Theorem 2 is purely algebraic, requiring only basic linear algebra multi-member and single-member operations.
The closed-loop transfer matrices for the first and second mode of operation are computed to be of the following forms:
where
Furthermore, it can be observed that, in Mode 1, there is a canceled-out polynomial of the form
The free controller parameters
,
,
,
, and
can be used to satisfy linear approximant closed-loop stability and the appropriate adjustment of the gains of the closed-loop transfer matrices. Indicatively, the stability of the linear approximant closed-loop system in both modes of operation is guaranteed if, and only if, the following inequalities are satisfied:
If the above inequalities are satisfied, by setting the and asymptotic command following is achieved for and .
3.6. An Optimization Criterion for Approximate Command following and I/O Closed-Loop Stability
For the specifications of the approximate command following and I/O asymptotic stability of the closed-loop system, the goal is to determine the degrees of freedom in (42a) and (42b), such that (a) all transfer functions in diagonal elements of the closed-loop transfer matrices are stable and (b) the steady state value of the step responses of these transfer functions are as close as possible to one. In addition to these two requirements, two more requirements are imposed. The first is to minimize the imaginary part of the poles of the diagonal transfer function, and the second is to achieve stable system poles of the closed-loop system. The satisfaction of all the above requirements, through appropriately bounded input variables, is translated here to the minimization of an appropriate cost function. To introduce this cost function, a set of definitions will be presented.
First, consider the following “decoupled” unitary step external command:
where
is the unitary step signal. The respective response of the
-th output of the closed-loop system is denoted by
The respective response of the
-th input of the closed-loop system is denoted by
Note that, in the closed-loop system, the input vector is produced by the controller. So, it depends upon the index of the system description model. Moreover, from the form of the closed-loop transfer matrix, it is observed that
The following cost criterion is defined using the above definitions:
where
is the denominator polynomial of the argument rational function, and where
where
is the absolute value of the imaginary part of the argument complex number, and
is the ring of polynomials of
. The 2-norm type cost function is defined to be of the following finite time horizon form
where the parameter
is the simulation time of the computational experiment. The weighting factors
,
,
satisfy the equality
Note that, for different
, the cost functions in (44) are decoupled among themselves, with respect to the free parameters of the controller matrices. In other words, different cost functions have different free parameters. It is mentioned that the rest of the free parameters in (42a) and (42b), namely the elements of the matrices
and
, do not participate in the diagonal elements of the closed-loop transfer matrices. To verify this property, post multiply (27) by
and use (A19) and (A21) (see
Appendix A) to compute
. Thus, the following optimization problems are defined:
under the algebraic constraints
where
is a small enough, nonnegative real number set by the designer. The constraint (46) formulates the design requirement of the approximate command following.
To solve the present optimization problem, the metaheuristic algorithm in [
34], being an extension of the algorithm in [
12,
13], will be used to compute the respective free parameters of the controller. Clearly, the solution derived using the metaheuristic algorithm is suboptimal.
4. Robot Vehicle Application
This section presents the robot dynamics and controllers achieving common noninteracting control with simultaneous common partial output zeroing. Additionally, approximate command following and closed-loop stability are derived.
4.1. Dynamics of Robot-Tracked Vehicles
According to [
15,
16], the dynamics of a tracked vehicle in horizontal planar motion have two distinct modes. In the first mode, the vehicle follows a straight path, while in the second, the vehicle follows a curved one. Switching from the first mode (not turning) to the second mode (turning) and vice versa is governed by an appropriate switching function. In the description of the vehicle, the dynamics of the driving electric motors are also included. According to [
15,
16], the dynamic mode of the tracked vehicle for each mode of operation can be described by the following sets of equations:
Straight Path Mode (
)
Curved Path Mode (
)
Note that
is the linear velocity of the vehicle,
is the angular velocity of the vehicle,
is the heading angle of the vehicle,
and
are the left and the right track motor currents,
and
are the left and the right track motor angular velocities,
and
are the left and right track generated forces, and
and
are the left and the right motor voltage supplies (actuatable inputs). The vector of the angular velocity of the tracked vehicle is perpendicular to the motion plane, and the vector of the linear velocity is along the heading axis of the vehicle. It is important to mention that switching between modes of operation occurs through an appropriate switching condition variable
, where
is accomplished through the following rule:
Regarding the notation of the parameters of the switching model, is the mass of the vehicle, is the angular mass of the vehicle, is the width of the vehicle from track to track, is the mechanical resistance of the tracks to rolling forward, is the mechanical resistance of the tracks to sliding forward, is the mechanical resistance of the tracks to turning, is the inductance of the motors, is the resistance of the motors, is the angular mass of the gears, is the mechanical resistance of the gears to rotation, is the gear ratio of the gearbox, is the current–torque ratio of the electric motors, is the radius of the sprocket wheels, is the compliance of the tracks, and is the lateral friction of the tracks.
4.2. Two-Model State Space Representation
The switching model of the robot-tracked vehicle is in the form (10a) with
and
In the present UGV case, the matrices in (10a) have the following forms:
where
where
and where the nonzero elements of
,
, and
are:
4.3. Performance Variables and Design Goals
A favorite control objective in vehicle control systems is the independent control of the translation and the rotation of the vehicle. Here, the problem of the independent control of the forward velocity and the heading angle of the tracked vehicle will be examined. In this case, the performance output matrix
is a
matrix of the form
The performance output vector is of the form .
4.4. The Controller
In general, the most accurate real-time measurements of the variables of such a vehicle are the left and the right track motor currents. Additionally, the orientation of the vehicle, i.e., the heading angle with respect to the earth frame, can accurately be measured in real-time using optical sensors (indicatively, see [
35]). In what follows, it will be assumed that the motor currents and orientation angle of the vehicle are measurable, i.e., it holds that
. The measurement output matrix
is a three-by-nine real matrix of the form
As already presented in
Section 3, the controller is of the static measurement output regular form of relation (2). In the present case, it holds that
and
are constrained to be invertible.
The controller design procedure described in
Section 3 will be used for the independent regulation of the translation and the rotation of the vehicle. In addition to the design requirement in
Section 3, it is required for the controller to satisfy the following: (i) approximate asymptotic command following, (ii) stability of appropriate variables, and (iii) suboptimal transient response.
Recall that in the present UGV case, the multi-model state space representation in (10a) and (10b) is used with
,
, and
; since the tracked vehicle has two distinct modes of operation, it also holds that
. Clearly, in the first mode, the linear speed must be independently controlled while the heading angle must remain at its initial value. In the second mode, both performance outputs must be independently controlled. This design requirement translates to
,
,
, and
, while
After appropriate algebraic manipulations, it can be verified that the necessary and sufficient conditions in Theorem 1 are satisfied. Hence, the common noninteracting control problem with simultaneous common partial output zeroing is solvable. According to Theorem 2, the general form of the controller matrices solving the problem are
where
,
,
,
, and
are free controller parameters to be selected by the designer, and
,
.
The transfer matrices mapping the external commands to the performance outputs, when using the above solutions for the controller matrices for the two modes of motion, are derived to be
where
and where
From the above expressions, it is observed that in the second model, the transition poles are the system poles, while in the first model, there exists a fifth-order polynomial being canceled out of the closed-loop transfer matrix. The canceled-out polynomial is computed to be
, where
The double pole at zero of
is an inherent characteristic of the second model resulting from the conditions of constant heading angle and zero angular velocity. Thus, with regard to stability, the roots of
,
,
, and
must lie on the left half complex plane. To this end, using the Hurwitz criterion (see [
36]), the appropriate inequality constraints of
are derived.
Here, the following nonsymmetric version of the approximate command following is selected:
The above selection is translated to
,
, and
. Choosing
the equalities
and
are satisfied and
Clearly, is constrained to satisfy the inequality , which is not in contradiction with the Hurwitz inequalities mentioned above. Thus, an appropriate can be derived.
Overall, the controller parameters guarantee asymptotic command following and stability for the performance outputs in the second mode, while in the first mode, the steady state error of the forward velocity is stable and appropriately bounded. The remaining controller parameters will be used to improve the closed-loop performance through a 2-norm performance criterion and a metaheuristic algorithm.
4.5. Optimized Selection of the Free Controller Parameters
The goal is to determine the free controller parameters
,
, and
such that the transient performance of the closed-loop system is improved. The degree of improvement will be evaluated through the cost criterion presented in (45). The controller parameters will be determined using the metaheuristic algorithm in [
34], being an extension of the algorithm in [
12,
13]. This way, a suboptimal solution of the parameters will be derived. According to
Section 3.6, the step responses of the closed-loop transfer functions
,
, and
are denoted by
,
, and
, respectively. Furthermore, let
,
,
,
,
, and
be the respective actuatable input signals generated by the controller. Considering that the parameters
and
have been predefined towards the asymptotic command following in Mode 2, the cost functions in (44) are specified as follows:
Note that
and
are decoupled with respect to the controller parameters, i.e.,
depends only upon
while
depends only upon
and
. The goal of the metaheuristic algorithm is to minimize
and
under the stability constraints presented in the previous subsection and
. The metaheuristic algorithm will be executed two times, one for determining
and one for determining
and
. Let (see [
16])
The parameters of the cost criteria and the algorithm are
where
,
, and
are the initial central value, initial half width of the search area, and the convergence bound of the argument parameter, respectively;
is the number of simulation experiments per repetition;
is the total number of repetitions per parameter update; and
is the total of simulation experiment after which the algorithm terminates without parameters converge (see [
12,
13]). The metaheuristic algorithm provides the following values of the controller parameters using the above parameter:
yielding
The resulting feedback matrix is
Using the derived general solution of the precompensator for the UGV case, the precompensator is specified to be
4.6. Simulation Results
The above controller will be used to execute the simulation providing the closed-loop response of the tracked vehicle for a complex maneuver to demonstrate the performance of the proposed control scheme. All computational experiments were produced using the MATLAB
® R2021a/Simulink simulation software. The variable step ODE45 (Dormand–Prince) solver was selected to produce the simulation results, setting the relative and absolute tolerance and the relative tolerance of the solver to
. The maximum, minimum, and initial step sizes were set to “automatic”. The elements of the external command vector
will be chosen to be of the smooth sigmoid form, namely
where
denotes the inverse Laplace transform of the argument signal,
(
,
) are distinct positive real numbers,
(
) are appropriate time delays being nonnegative real numbers,
(
) are the amplitudes of the respective external commands, and
is the unit step signal.
In order to produce a unified dynamic model, the mathematical description of the motion of the tracked vehicle will be modified to include a logistic function (see [
37]) and an additional friction coefficient, incorporating both modes of operation to be used for the simulation and to avoid numerical errors around the switching area. This description covers both models of the vehicle, and the transition between models is approximated by sharp but continuous responses, interpreted as friction terms. The switching nature of the model is guaranteed by a rapid decrease in the angular velocity of the vehicle when it enters the first mode of operation. To this end, consider the function:
where
defines the flex points of each term in (52),
is a positive real number adjusting the steepness of the curve, and
is a positive real number defining the maximum value of the curve (see
Figure 1).
Define
. Using this signal as an argument variable of the above-presented modified logistic function, the two modes of operation of the tracked vehicle are unified, reducing to the respective form of Mode 2, while the linear and angular velocity equations are modified in the form
where
is the additional friction term.
Consider the model parameters and the controller parameters presented in
Section 4.5. The logistic function parameters, the additional friction term, and the external commands amplitudes are selected to be
Figure 2,
Figure 3,
Figure 4,
Figure 5,
Figure 6 and
Figure 7 present the simulation results for the robot-tracked vehicle’s closed-loop response. In particular, in
Figure 2 and
Figure 3, the forward and angular velocities are presented, respectively; in
Figure 4, the heading angle is presented; in
Figure 5 and
Figure 6, the motor currents and angular velocities are presented; in
Figure 7, the track forces are presented; and, in
Figure 8, the voltage supply to the motors of the vehicle are presented. With respect to the performance variables
and
, it can be verified that the external commands are followed accurately (see
Figure 2 and
Figure 4). With respect to the voltage supply to the electric motors (see
Figure 8), it is observed that they are smooth and appropriately bounded, being easily implementable. Note that the voltage supplies, motor currents, and motor velocities are within the respective technical specifications of the motor. For more details, see the Data Sheet of the FA-130 Mabuchi Motor and [
16]. Besides, the performance outputs, the remaining state variables (see
Figure 3,
Figure 5,
Figure 6 and
Figure 7) are appropriately bounded.
To demonstrate the effectiveness of the proposed control scheme, with respect to uncertainties, what follows will be considered the mechanical resistances of the tracks to rolling (
and
) and the mechanical resistance of sliding forward (
). Moreover, the mechanical resistance of the gears to rotation (
) and the lateral friction of the tracks (
) are not perfectly known to the designer. So, these parameters are considered to be uncertain, and they are expressed as follows:
,
,
,
, and
, where
,
,
,
, and
are the respective nominal values of the parameters upon which the controller was evaluated, and where
,
,
,
, and
are the respective uncertainties. Let
and
be the nonlinear closed-loop responses of the performance variables. As already mentioned, the controller is designed using the parameters
,
,
,
, and
instead of their real values. Furthermore, it is mentioned that the external commands and the rest model parameters are chosen as in the previous paragraphs. Finally, let
and
be the closed-loop response of the performance variables produced after applying the controller to the original uncertain system. The following metrics will be used to compare the respective responses:
where
and
denote the finite horizon 2-norm and infinity norm (indicatively, see [
38,
39]) of the argument signal
, i.e., it holds that
where
denotes computational experiment simulation time. For demonstration purposes, it will be assumed that
.
Performing a series of computational experiments for
,
,
,
, and
, where
, it was observed that the second performance variable, namely the heading angle of the vehicle, is minimally affected by the presence of the uncertainties. The maximum values of
and
were found to be
and
. Regarding the first performance variable, namely the speed of the UGV, it is observed that the uncertainties affect the closed-loop response at acceptable levels. The maximum values of
and
were computed to be
and
, respectively. In
Table 1, the maximum values of
,
,
, and
are presented for different values of
. The present simulation results demonstrate some of the perspectives of the proposed controller in a real environment.
4.7. Comparison to Switching Controllers
In this subsection, the common control design scheme derived in the previous subsections for the robot-tracked UGV is compared to a switching control design scheme satisfying the same design requirements with two different controllers, one for each mode of the UGV. These two controllers are switched by the external commands and measurement of the heading angle of the vehicle, in the sense that the first is triggered when the external vector command corresponds to a straight path and constant heading angle, while the second is triggered when the external vector command corresponds to a curved path. The curved path period is completed, and consequently, the straight path controller is switched on when the heading angle has reached 97% of the final value of its variation. Clearly, this time period covers the settling time of the heading angle.
Regarding the design requirements of the switching scheme, the first requirement is to satisfy (51), i.e., the first controller satisfies the first equality in (51), while the second controller satisfies the second equality in (51). Following Remark 2 for
, the following general forms of the controller matrices are derived:
where
and
are the
and
elements of the feedback and precompensator matrices, respectively, for
. The free elements of the feedback matrices
and
, i.e., the elements
,
,
,
,
,
,
, and
, will be chosen such that the eigenvalues of the matrices
(
) are stable. The free elements of the precompensator matrices
and
, i.e., the elements
,
,
, and
, must be chosen such that
,
. In addition to the above constraints, the asymptotic command following is also achieved for the nonzero elements for the closed-transfer matrices of the two modes. In what follows, the determination of the free controller matrix elements will be carried out through the metaheuristic algorithm mentioned in
Section 4.5, where the model parameters are those presented in
Section 4.5. After applying a series of computations, the two sets of controller matrices are determined to be
Here, the simulations are performed with zero initial conditions. Furthermore, the two external commands are chosen to be
The simulations begin with the second controller of the switching scheme. At
the heading angle has reached the 97% of its final value, i.e.,
, and so the first controller is switched on. From that point on, since no change in the orientation command is imposed, the straight path controller of the switching scheme is preserved. The closed-loop responses of the performance outputs for both the common controller case and the switching controller case are presented in
Figure 9 and
Figure 10, where it is observed that, for the present metaheuristic evaluation of the controller-free parameters, the closed-loop response of the common controller is preferable as compared to the switching controller. Indicatively, consider the performance metric
,
. According to the simulation results, in the common controller case, the metric is evaluated to be
and
, while in the switching controller case, the metric is evaluated to be
and
.
5. Discussion
In its general form, the problem of common noninteracting control with simultaneous common partial output zeroing is applicable in the sense that some of the outputs are independently controlled while some of them are retained at their nominal values. Such a requirement can be met in redundant manipulators, multi-mode processes, and robotic vehicles.
The approach followed to solve the problem at hand is purely algebraic, facilitating the computation of the quantities in the solvability criteria as well as the general solution of the controller.
The application of the present results to the case of a UGV is quite useful and provides useful insight into the properties of the problem at hand. From the point of view of the availability of measurements, the measurement variables used by the controller appear to be the preferable choice.
Significant issues, such as the switching stability of the controlled UGV and the robustness of the derived results with respect to measurement noise, are currently under investigation.
A future perspective of the results of the present paper is the extension to other problems in the wider class of noninteracting control; indicatively, see the common noninteracting control problems in [
40,
41], as well as the seminal works in [
42,
43,
44], where pure noninteracting control problems are studied for nonsquare systems with static controllers, including nonsquare precompensators.