1. Introduction
The MIMO multiloop multirate continuous-discrete system is a mathematical model that describes the stabilization system of a drone with a multiprocessor onboard computer. Such systems are considered due to the need to optimize the performances of systems with different dynamic modes and efficiently use available hardware resources [
1,
2]. Usually, MIMO stabilization systems are considered, which have vector characters of control, disturbances, parameters of the control system, and outputs [
3,
4,
5,
6,
7]. There are two approaches to the description of the multirate continuous-discrete system: models in the time domain and models obtained using the discrete Laplace transform or z-transform. The advantages of describing the considered class of systems in the time domain are detailed in [
1]. A large number of publications have solved the problem of describing multiloop multirate continuous-discrete systems in the time domain [
1,
8,
9,
10,
11,
12].
The description of the considered system in the domain of discrete Laplace transform or z-transform does not lose its relevance.
The drone stabilization system is considered for the first time [
13] as a multiloop multirate control system (
Figure 1). Similar systems are considered in [
14,
15,
16].
We consider the stabilization system of the lateral movement of the aircraft-type drone when solving the problem of a UAV landing on a non-aircraft-carrying ship [
17]. The stabilization system of the lateral movement for the autonomous drone, with an onboard computer in the control loop, is considered in [
18], in the form of a structural diagram (
Figure 2). The system in
Figure 2 is a multiloop multirate continuous-discrete system.
In this figure,
—yaw angle,
—roll angle,
, and
—angles of deflection of the rudders and ailerons, respectively;
and
—stabilization errors. The corresponding transfer functions are:
—steering track,
—kinematic link,
—extrapolator,
—gyrotachometer,
—aerodynamic coefficients,
,
,
,
and
with the corresponding indices–time constants. The transfer functions of the control object
,
,
,
can be determined by transforming, using the Laplace transform for a system of differential equations, describing the rotational motion of the aircraft around the center of mass in the horizontal plane [
18]. The transfer function
, which characterizes the effect on the dynamics of stabilization of elastic vibrations of a moving UAV, is not taken into account in further consideration. The model of the
is given in [
18].
are the sampling periods of digital feedback circuits of control signals.
To solve the problem of analysis or synthesis of the system shown in
Figure 2, we need to obtain an adequate mathematical model.
The mathematical models obtained using the discrete Laplace transform or z-transform are considered with sufficient simplification. The studies were limited to SISO systems and generally resulted in very cumbersome formulas [
1].
Although the problems of analyzing and synthesizing multirate systems appeared a long time ago, there is no explicit and complete methodology for constructing their equivalent mathematical model [
1].
Development in the issue of constructing mathematical models of multi-loop control systems with different sampling periods specified using transfer functions was achieved, thanks to the works in references [
19,
20,
21,
22]. In these works, it was proposed to replace the pulse elements (keys), which sample signals with different sampling periods, with an equivalent system containing one key and some delay elements. Such an equivalent system operates with a uniform time sampling period. Further description of the considered system is carried out using modified z-transformations, which makes it possible to introduce the time delays of the equivalent system into the description [
21]. A description of multiloop multirate systems, in the form of transfer functions and two different algorithms for calculating the characteristic polynomial, was proposed in [
19,
21]. These two methods, for constructing the characteristic polynomial of a multiloop multirate control system, are equivalent [
13].
The main disadvantage of the approach proposed in [
19] is the complexity of the implementation associated with the inconvenience of numerical calculations of topological operations [
13].
As in most other works, attention in [
13] was paid to the construction of characteristic equations of multiloop discrete control systems, containing impulse elements with different sampling periods, and not to the construction of the entire mathematical model of the system suitable for analysis and synthesis.
Thus, existing approaches to the construction of the models of multiloop multirate continuous-discrete systems, in the form of transfer functions, allow one to obtain equations suitable only for studying stability. The existing approaches do not allow a wide range of analysis and synthesis of the specified class of control systems and reduce the efficiency of their research.
Note that the problem of synthesis and analysis of linear continuous-discrete multiloop UAV stabilization systems contains two essential aspects of complexity: the system’s multidimensionality and the multirate of sampling chains. In the case of commensurate periods (their multiplicity to a certain “effective” period), the second aspect is equivalent to the first, significantly strengthening it at large numbers of the multiplicity of sampling periods to the “effective” sampling period. In this case (rationally commensurate sampling periods), it is possible to transform the mathematical model of a multirate system to an equivalent mathematical model of a single-rate system of increased dimension. In the theory of impulse systems, such transformations are usually performed based on the “parallelization” of impulse circuits [
1].
In this article, we propose a complete description of a general approach to the construction of an equivalent single-rate model, which makes it possible to transform the vector–matrix model of the original MIMO multiloop multirate system into a vector–matrix model of a single-rate system, in the form typical for vector–matrix models of continuous systems. The construction includes complete descriptions of the equivalent T-single-rate and NT-single-rate, where T is the largest common divisor of the sampling periods of the system, and N is a natural number that is the least common multiple of the numbers characterizing the sampling periods of the system. Detailed descriptions of these characteristics are given below, describing the central element of the construction—a structural invariant of sampling chains—a matrix, which, in the general case, describes a complex situation of the mutual influence of sampling processes in digital circuits of a closed-loop MIMO multirate stabilization system with different sampling periods. The dimension and complexity of this matrix are due to the mutual numerical characteristics of digital circuits.
2. The Mathematical Model of the Stabilization System
To build an equivalent mathematical model of the multirate system, we write the fundamental equation in the input–output form of the Laplace transform. The complete construction of the basic mathematical model in the input–output form of the Laplace transform for the MIMO multiloop multirate continuous-discrete system is given in [
18,
23].
The continuous object equation and analog circuit stabilization system equation from the object to the sampling keys are
where
are the matrices of transfer functions. Usually,
, and in the future, we will take this into account.
is a vector of the sampling variables with periods , respectively (among which there may be equal ones), where are natural numbers. The control vector is a vector . —the vector of outputs.
For the system shown in
Figure 2,
. For the variables of the vector
, the sequence of sampling periods is
. The vector of outputs
. The “k-th” component of the vector of control,
, where
—the transfer function of parallel digital circuits,
—discrete Laplace transforms sampled concerning the sampling periods
of the variables
, respectively.
As considered in
Figure 2, schemes
,
, and the transfer matrixes,
,
.
Writing the equations for the control in a vector–matrix form:
where
is a matrix with elements
and vectors
with elements
and
. The symbol “*” at the bottom marks the property of periodicity of the matrix
and the vector
for the corresponding periods. All elements of
and
satisfy the relation
. For
and
is fulfilled [
24,
25]
Finally, we obtain a mathematical model of the MIMO multiloop multirate continuous-discrete drone stabilization system, which takes into account both the influence of all sampling periods of the system and its multiloop, including the cross-connections of control channels
where
.
Due to the presence in the description of factors that take into account the different sampling periods, the presented mathematical model is not suitable for solving problems of analysis and synthesis. It is a basic model for constructing an equivalent mathematical model, which, in turn, can be used to extend analysis and synthesis methods for multidimensional single-rate systems to the class of multirate systems.
3. Results
As noted earlier, for a fundamental possible transfer of the synthesis and analysis methods of MIMO continuous and discrete systems, to the class of continuous-discrete multirate drone stabilization systems, it is necessary to construct equivalent mathematical models of a single-rate system in a form characteristic for vector–matrix models of continuous systems. Approaches to the construction of equivalent mathematical models are briefly presented in [
23]. Now, we give their full description.
3.1. The Equivalent T-Single-Rate Model of the Stabilization System
The given model of a multirate drone stabilization system can be transformed into a model of a single-rate system with a sampling period , which is the greatest common divisor of the sampling periods .
We will proceed from Equation (11)
Further, until the final results, the symbol in the designation of the corresponding impulse transformation will be omitted, i.e., for any function , instead of we will write .
Consider Equation (4). By property (3), we have
Replacing
in this relation by
, where
,
are natural numbers, we have
where for any function
, we further take
.
Setting in Equation (9),
, we write
where
Consider the components of the vector
—the least common multiple of numbers
such that
, where
are natural numbers. Then Equation (12) can be written in the form
Note that, for any
we have, setting
,
due to the periodicity of the function
in
with a period
.
Due to (14), the right-hand side of (13) contains only values
which means that there is a linear transformation of the vector
to
. We can write
where
a numerical matrix, the rules for finding, which we indicate below.
Substituting Equation (15) into (10), we find
Writing the equation
where
,
, and using (15), (16) we find
Applying the identity [
26,
27,
28]
, we also find
The defining matrix
as block type, where
matrix, and introducing the notation
we define the elements of the block representation of this matrix
where
.
From (20). we obtain the equation for the elements of the matrix
The construction of the matrix (20) of an equivalent T-single-rate system allows one to find the outputs of the closed original system, sampled at the moments .
Defining, for this purpose, the block type matrix
we find
where
. 3.2. The Equivalent T-Single-Rate Model of the Stabilization System: The Analog Control
In order to not change the notation of the matrices, we assume that is a vector with elements .
Instead of (4)–(6), we now consider the equations
Determining, as above, the extended vectors and matrices (14), we find
Hence, what follows is the corresponding representation for the impulse transformation of the vector of outputs, similar to Equation (22).
3.3. The Equivalent NT-Single-Rate Model of the Stabilization System
The construction of the sequence allows, of course, to isolate the sequence . At the same time, if the last sequence is determined from the discrete Laplace transformation , then this one can see some additional possibilities of analyzing the dynamics of the initial multirate stabilization system. Therefore, it is of interest to consider —a single-rate equivalent model—an impulse stabilization system with a sampling period , where, as before, is the least common multiple of numbers .
The starting point for building an
—single-rate equivalent model is the [
29]
Taking into account
,
where the numbers
.
Equation (25) can be written as
where
,
,
.
Using Equation (25) and the block matrix
with dimension
for a vector
, we can write
and
We now introduce the matrix and block matrix with a dimension that allows writing .
Therefore, multiplying (27) on the left by the matrix
, we obtain
Introducing, for convenience,
,
and
, by Equation (28), we find
where
.
The relations (29) and (30) define the equivalent single-rate model of the closed-loop drone stabilization system with a sampling period .
An alternative approach to constructing an equivalent single-rate model of the closed-loop drone stabilization system with a sampling period
is given in [
28].
3.4. The Structural Invariant of Sampling Chains of the Stabilization System
3.4.1. The Construction Method of the Structural Invariant of Sampling Chains of the Stabilization System
When, in stabilization system (4)–(6), the same sampling period is
T, then
and
. The block vectors and matrices introduced above degenerate (consist of one block element). Equation (18) takes the form
Here, the matrix degenerates into one block with diagonal elements equal to one.
In the general case of modeling multidimensional multiloop multirate stabilization systems, there is a complex situation of the mutual influence of sampling processes in digital circuits of a closed system with different sampling periods. A characteristic element of the description of this complex situation is the matrix , i.e., transformation (15).
The dimension and complexity of this matrix determine the dimension and complexity of the processes in the stabilization system, and is due to the mutual numerical characteristics of the digital circuits—the numbers that determine the relationship between the sampling periods and —is the least common multiple of the indicated numbers.
As can be seen from Equations (13) and (15), only numbers are nonzero entries of this matrix. Since these numbers have the meaning of relative sampling period densities (ratios of the number of “rare” samples to the number of “frequent” ones), the matrix will be called the matrix of tilt densities. Since the structure of the system, in particular, the arrangement of digital chains and sampling keys, does not affect the matrix in any way, it is a structural invariant of the sampling chains.
To complete the formation of an equivalent single-rate model of the stabilization system, we now indicate a method for calculating the elements of the matrix
. We will look for them in the form of the matrix
. To do this, consider the equation
which takes place by Equations (13) and (15) with the notation
Considering Equation (32), one should take into account the periodicity relation
as property (14) of functions
.
For this purpose, consider the right-hand side of relation (32) and the set of numbers
with fixed indices
.
The numbers of this set have the property
be the largest of the numbers
,
. This number has the form
where
is the smallest of the numbers
.
Due to the , at least one of the numbers is less then (at least one of the numbers is greater than (1), which means .
As a result, all numbers of form (33) exceeding
are numbers of the form
. Now, let, for arbitrarily taken
, such that
, we have
Subtracting, we find that
Because of the above, from (35), we obtain the following rule for calculating the elements
of the matrix
of quantization densities for fixed indices
where the values
are calculated by (33).
3.4.2. The Algorithm of Construction the Structural Invariant of Sampling Chains of the Stabilization System
In this part, we present the corresponding algorithm. We build the table to construct the matrix
(
Table 1).
The top line of the table is a series of natural numbers from 1 to , some of them are marked with an asterisk —*. The second line contains values equal to 1 under the marked upper numbers and 0 otherwise.
The rule for marking the numbers on the top line is as follows. The values are calculated
If
, then the marked number
, if
then
.
In this way, the bottom line defines the values
, and
Thus, this method determines the entire row of the matrix , i.e., diagonal elements with fixed indices of all block matrices equal to either or 0.
Adopting
and applying the algorithm n times for
, we finally find the first row of the block representation of the density matrix
This row defines the entire matrix
using a sequential circular permutation. For example, the second line will be
For example, for
, the entire matrix has the form
Note the possibility of forming the given blocks of the matrix . For example, we are talking about a block . Diagonal elements should be defined with, given . Then:
Thus, the construction algorithm is fully formulated.
3.4.3. Example
Consider a part of the lateral motion stabilization system—a modified roll stabilization contour (
Figure 3).
For the system shown in
Figure 3, according to the notation introduced in
Section 2, we have:
,
,
,
,
,
,
,
,
.
From the above diagram, it is clear that the methods based on the analysis of the location of the poles are not available with such a representation of the system.
Using Equations (1)–(6), we obtain an expression for the output of the system
To construct an equivalent model, the usual matrix operations are used. A separate issue is the construction of a matrix of sampling densities.
For the considered in a
Figure 3 system
(defined above). Using the above equations, we obtain
.
For block
, we have the equations
Therefore, .
Therefore, .
Therefore, .
The entire density matrix has the form
The difference in the sampling densities is reflected in the matrix
by the inhomogeneous filling of its blocks. When the sampling density in all circuits is the same, i.e.,
and
, then the matrix
turns out to be maximally filled—all of its block elements
turn out to be equal
where
—identity
matrix. Hence,
.
For the system under consideration, to construct a T—single-rate equivalent model, we will have the following vectors and matrices used to construct the model:
where the dimension of the matrix block
—2 × 1;
where the dimension of the matrix block
—1 × 2;
.
Knowing the indicated matrices, we get
or
By separating , we get an equivalent single-rate model as a result. The rest of the equivalent models given in the article are obtained in the same way.
A comparative analysis of the step response of stabilization systems for the different types of models was carried out in [
18]: a complete continuous-discrete system, taking into account the presence of different sampling periods, a model of a system in which all quantization periods are considered small, and are excluded from the model, and the model in which the periods of quantization are taken the same. In this case, the transfer functions were not transformed and their parameters remained the same. It was shown in [
18] that a simple change of sampling periods without performing equivalent transformations leads to various transient processes of systems.
Now we simulate the continuous-discrete system considered in
Figure 3 with different sampling periods
and
. In the system
,
,
,
,
,
.
The expression for the obtained equivalent single-rate representation of the system in the z domain with a sample period
T, which is the greatest common divisor of the sampling periods, has the form
A Simulink calculation scheme is shown in
Figure 4.
As equivalent models, we will take a continuous system with the same transfer functions, but with continuous circuits in the feedback, a discrete single-rate system with equal sampling periods in feedback loops, and an equivalent system with a sampling period T, which is the greatest common divisor of the sampling periods.
Analyzing the results, we conclude that a simple replacement of a multirate system with a continuous and single-rate system leads to unsatisfactory results. In the case of a single-rate system, the system becomes unstable at certain values of the sampling period. The transient process of the constructed equivalent system practically coincides with the transient process of the original model.
4. Discussion
If the system contains several digital information-processing circuits with different periods, then its study is greatly complicated. In the general case of incommensurate quantization periods, the mathematical apparatus of the study is practically not developed. In the case of rationally commensurate quantization periods (the periods are integer multiples of a certain “effective” period), it is possible to transform a multirate system to an equivalent single-rate of the increased dimension. In the theory of impulse systems, such transformations are usually performed based on parallelizing impulse circuits.
As a result, vector–matrix equivalent single-rate models are obtained using the discrete Laplace transform of MIMO multirate multiloop stabilization systems given by relations (4)–(6), or by structural diagrams, an example of which is shown in
Figure 2.
Because of the construction of equivalent mathematical models, it becomes fundamentally possible to transfer the synthesis and analysis methods of continuous MIMO systems to the class of continuous-discrete MIMO multiloop multirate systems.
On the one hand, we obtained an effective tool for constructing an equivalent single-rate form for describing MIMO multiloop multirate stabilization systems, to which the methods of synthesis and analysis of continuous MIMO systems can be extended. On the other hand, implementing the proposed technique leads to cumbersome calculations and work with large transfer function matrices. The dimension of the equivalent system is influenced by the mutual numerical characteristics of digital circuits. Thus, the general construction of the system model complication is indicated due to the increase in the dimension of the equivalent model at significant values of the clock multiplicity numbers and the number of quantization chains. This is illustrated by an example, including constructing a matrix of quantization densities—a structural invariant of quantization chains.
Linear MIMO control systems—systems with multiple inputs and multiple outputs —are part of the linear system class. Therefore, their theory is based on the same foundation as the theory of SISO systems [
30,
31,
32]. If we talk about models of linear systems of a complex domain (domains of Laplace transformations, z-transformations), then this dynamics is determined by the input–output relations in the form of the corresponding rational functions-scalar for a one-dimensional matrix for multidimensional systems. At the same time, the class of multidimensional systems (more complex by the property of multidimensionality) in comparison with the class of one-dimensional systems requires further development of the theory of linear systems, in connection with special mathematical problems of research, which are absent in the class of one-dimensional systems or appear in a much simpler form. Thus, for example, the solution of several problems requires the use of equivalent transformations of polynomial matrices that determine the transfer matrix of the system to a relatively simple form.
Implementing the corresponding computational procedures is a difficult task, the solution of which becomes more complex with an increase in the dimension of the system. Thus, for example, the solution of several problems requires using equivalent transformations of polynomial matrices that determine the transfer matrix of the system to a relatively simple form. Methods for constructing the most significant common divisors of polynomial matrices also play an essential role in implementing these transformations.