Next Article in Journal
UAV Patrolling for Wildfire Monitoring by a Dynamic Voronoi Tessellation on Satellite Data
Next Article in Special Issue
A 3D Vision Cone Based Method for Collision Free Navigation of a Quadcopter UAV among Moving Obstacles
Previous Article in Journal
Machine Learning-Assisted Adaptive Modulation for Optimized Drone-User Communication in B5G
Previous Article in Special Issue
Computationally-Efficient Distributed Algorithms of Navigation of Teams of Autonomous UAVs for 3D Coverage and Flocking
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multiloop Multirate Continuous-Discrete Drone Stabilization System: An Equivalent Single-Rate Model

Department of Informatics and Control in Technical Systems, Sevastopol State University, 299053 Sevastopol, Russia
*
Author to whom correspondence should be addressed.
Drones 2021, 5(4), 129; https://doi.org/10.3390/drones5040129
Submission received: 6 September 2021 / Revised: 27 October 2021 / Accepted: 28 October 2021 / Published: 1 November 2021
(This article belongs to the Special Issue Conceptual Design, Modeling, and Control Strategies of Drones)

Abstract

:
The article discusses the UAV lateral motion stabilization system, as a MIMO multiloop multirate continuous-discrete system, specified in the form of an input–output model in the domain of discrete Laplace transform or in the form of a structural diagram. Approaches to the construction of equivalent T and NT single-rate models for MIMO multiloop multirate continuous-discrete systems are considered. Here, T is the largest common divisor of the sampling periods of the system, N is a natural number that is the smallest common multiple of the numbers characterizing the sampling periods of the system. The resulting impulse representations of the outputs of equivalent models are in the form of rational functions. The basis for the construction of these models is a matrix of sampling densities—a structural invariant of sampling chains. An example of the construction of the indicated matrix and an equivalent single-rate model are given. Obtaining equivalent single-rate models for MIMO multiloop multirate systems allows us to extend the methods of research and synthesis of MIMO continuous and continuous-discrete systems to a common theoretical base—the theory of polynomials and rational functions, which are typical elements of the description of these classes of systems.

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, δ N , and δ E —angles of deflection of the rudders and ailerons, respectively; ε N and ε E —stabilization errors. The corresponding transfer functions are: Φ S T ( s ) = W P T ( s ) 1 + W P T ( s ) —steering track, W K L ( s ) = a 41 b 16 s 2 μ 2 + b 16 b 33 s μ 2 + 1 s ( b 36 s μ 2 + 1 ) —kinematic link, W e ( s ) —extrapolator, W G T ( s ) = s ( T G T 2 s 2 + 2 ζ G T T G T s + 1 ) ( τ G T s + 1 ) —gyrotachometer, b i j —aerodynamic coefficients, W P T ( s ) = k P T W P M ( s ) W μ y ( s ) , W P M ( s ) = e τ P M s s ( T P M s + 1 ) , W μ y ( s ) = 1 T μ y s + 1 , T and τ with the corresponding indices–time constants. The transfer functions of the control object W δ N ψ ( s ) = ψ ( s ) δ N ( s ) , W δ N γ ( s ) = γ ( s ) δ N ( s ) , W δ E γ ( s ) = γ ( s ) δ E ( s ) , W δ E ψ ( s ) = ψ ( s ) δ E ( s ) 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 W ψ δ H y ( s ) , 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 W ψ δ H y ( s ) is given in [18]. T 1 , T 2 , T 3 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
y ( s ) = W 0 ( s ) u ( s ) , x ( s ) = E ( s ) y ( s ) + B ( s ) u ( s ) ,
where W 0 ( s ) , E ( s ) , B ( s ) are the matrices of transfer functions. Usually, B ( s ) = 0 , and in the future, we will take this into account.
x ( s ) = [ x 1 ( s ) , x 2 ( s ) , , x r ( s ) ] T is a vector of the sampling variables with periods T 1 = n 1 T , T 2 = n 2 T , , T r = n r T , respectively (among which there may be equal ones), where n 1 , n 2 , , n r are natural numbers. The control vector is a vector u ( s ) . y ( s ) —the vector of outputs.
For the system shown in Figure 2, x ( s ) = [ x 1 ( s ) , x 2 ( s ) , x 3 ( s ) , x 4 ( s ) , x 5 ( s ) , x 6 ( s ) ] T . For the variables of the vector x ( s ) , the sequence of sampling periods is T 1 , T 3 , T 1 , T 3 , T 3 , T 2 . The vector of outputs y ( s ) = [ ψ γ ] . The “k-th” component of the vector of control, u k ( s ) = i = 1 r d k i ( s ) x i T i * ( s ) + u Z k ( s ) , where d k i ( s ) —the transfer function of parallel digital circuits, x 1 T 1 * ( s ) , x 2 T 2 * ( s ) , , x r T r * ( s ) —discrete Laplace transforms sampled concerning the sampling periods T 1 , T 2 , , T r of the variables x 1 ( t ) , x 2 ( t ) , , x r ( t ) , respectively.
As considered in Figure 2, schemes u Z 1 ( s ) = ψ Z , u Z 2 ( s ) = γ Z , and the transfer matrixes, E ( s ) = [ 0 0 W G T 1 s ( W k Z + b 73 ) W k L ( K Z ( 1 + K Z s ) W G T 1 0 0 0 0 ] T , W 0 = [ Φ S T W δ N ψ Φ S T W δ E ψ Φ S T W δ N γ Φ S T W δ E γ ] .
Thus
x ( s ) = U ( s ) u ( s ) = E ( s ) W 0 ( s ) u ( s ) .
Writing the equations for the control in a vector–matrix form:
u ( s ) = W e ( s ) W * ( s ) x * ( s ) + W e ( s ) v * ( s ) ,
where W * ( s ) is a matrix with elements W k i T i * ( s ) and vectors x * ( s ) , v * ( s ) with elements x i T i * ( s ) and u k z T * ( s ) . The symbol “*” at the bottom marks the property of periodicity of the matrix W * ( s ) and the vector x * ( s ) for the corresponding periods. All elements of W * ( s ) and x * ( s ) satisfy the relation x ( s + 2 π j / T ) = x ( s ) . For W * ( s ) and x * ( s ) is fulfilled [24,25]
( c b * ) T * = c T * b * ,   ( or   ( b * c ) T * b * c T * ) .
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
u ( s ) = G ( s ) W * ( s ) x * ( s ) + G ( s ) v * ( s ) ,
x ( s ) = C ( s ) W * ( s ) x * ( s ) + C ( s ) v * ( s ) ,
y ( s ) = L ( s ) W * ( s ) x * ( s ) + L ( s ) v * ( s ) ,
where G ( s ) = W e ( s ) I ,   C ( s ) = U ( s ) G ( s ) ,   L ( s ) = W 0 ( s ) G ( s ) .
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 T , which is the greatest common divisor of the sampling periods T 1 , T 2 , T r .
We will proceed from Equation (11)
y n T * ( s ) = 1 n k = 1 n y T * ( s + 2 π j n T ( k 1 ) ) ,   ( x y n T * ) T * = x T * y n T * .
Further, until the final results, the symbol T in the designation of the corresponding impulse transformation will be omitted, i.e., for any function z ( s ) , instead of z T * ( s ) we will write z * ( s ) .
Consider Equation (4). By property (3), we have
x * ( s ) = C * ( s ) W * ( s ) x * ( s ) + C * ( s ) v * ( s ) .  
Replacing s in this relation by s + 2 π j N T ( ν 1 ) , where N , ν are natural numbers, we have
x * ν ( s ) = C * ν ( s ) W * ν ( s ) x * ν ( s ) + C * ν ( s ) v * ν ( s ) ,
where for any function φ ( s ) , we further take φ ν ( s ) = φ ( s + 2 π j N T ( ν 1 ) ) .
Setting in Equation (9), ν = 1 , , N , we write
x ^ * ( s ) = C ^ * ( s ) W ^ * ( s ) x ^ * ( s ) + C ^ * ( s ) v ^ * ( s ) ,  
where
x ^ ( s ) = [ x 1 ( s ) x 2 ( s ) x N ( s ) ] T , x ^ * ( s ) = [ x * 1 ( s ) x * 2 ( s ) x * N ( s ) ] T , v ^ ( s ) = [ v 1 ( s ) v 2 ( s ) v N ( s ) ] T ,   C ^ ( s ) = d i a g [ C 1 ( s ) , C 2 ( s ) , , C N ( s ) ] ,   W ^ ( s ) = d i a g [ W 1 ( s ) , W 2 ( s ) , , W N ( s ) ] .
Consider the components of the vector x * ν ( s )
x i T i ν ( s ) = 1 n i k = 1 n i x i ( s + 2 π j N T ( ν 1 ) + 2 π j n i T ( k 1 ) ) .
N —the least common multiple of numbers n 1 , n 2 , , n r such that N = ν 1 n 1 , N = ν 2 n 2 , , N = ν r n r , where ν 1 , ν 2 , , ν r are natural numbers. Then Equation (12) can be written in the form
x i T i ν ( s ) = 1 n i k = 1 n i x i ( s + 2 π j N T [ ν + ν i ( k 1 ) 1 ] ) , i = 1 , , r .
Note that, for any p > N we have, setting p = N + ν ,
x i * ( s + 2 π j N T p ) = x i * ( s + 2 π j N T ν )
due to the periodicity of the function x i * ( s ) in s with a period 2 π j / T .
Due to (14), the right-hand side of (13) contains only values x i * ν ( s ) , ν 1 , N ¯ , which means that there is a linear transformation of the vector x ^ * ( s ) to x ^ * ( s ) . We can write
x ^ * ( s ) = Π x ^ * ( s )
where Π r N × r N a numerical matrix, the rules for finding, which we indicate below.
Substituting Equation (15) into (10), we find
x ^ * ( s ) = ( I r N + C ^ * ( s ) W ^ * ( s ) ) 1 C ^ * ( s ) v ^ * ( s ) .
Writing the equation
y ^ ( s ) = L ^ ( s ) W ^ * ( s ) x ^ * ( s ) + L ^ ( s ) v ^ * ( s ) ,
where y ^ ( s ) = [ y 1 ( s ) y 2 ( s ) y N ( s ) ] T , L ^ ( s ) = d i a g [ L 1 ( s ) , L 2 ( s ) , , L N ( s ) ] , and using (15), (16) we find
y ^ ( s ) = L ^ ( s ) { W ^ * ( s ) ( I r N + C ^ * ( s ) W ^ * ( s ) ) 1 C ^ * ( s ) + I m N } v ^ * ( s ) .
Applying the identity [26,27,28] I A ( I + BA ) 1 B = ( I + AB ) 1 , we also find
y ^ ( s ) = L ^ ( s ) ( I m N + W ^ * ( s ) C ^ * ( s ) ) 1 v ^ * ( s ) ,   y ^ * ( s ) = L ^ * ( s ) ( I m N + W ^ * ( s ) C ^ * ( s ) ) 1 v ^ * ( s ) .
The defining matrix Π = { Π ν p } , ν , p 1 , N ¯ as block type, where Π ν p r × r matrix, and introducing the notation
W * ( s ) = W ^ * ( s ) C ^ * ( s )
we define the elements of the block representation of this matrix
W * ( s ) = { W ν p * ( s ) } , ν , p 1 , N ¯ ,   W ν p * ( s ) = W * ν ( s ) ν p C * p ( s ) , ν , p 1 , N ¯ .
where ν p = d i a g [ π ν p 1 , π ν p 2 , , π ν p r ] .
From (20). we obtain the equation for the elements of the matrix W ν p * ( s )
W ν p σ μ ( s ) = i = 1 r π ν p i W σ i ( s + ( ν 1 ) 2 π j N T ) C i μ ( s + ( p 1 ) 2 π j N T ) .
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 k T , k = 0 , 1 , .
Defining, for this purpose, the block type matrix
H * ( s ) = ( I mN + W * ( s ) ) 1 ,   H * = { H ν p * } ,   ν , p 1 ,   N ¯ ,
we find
y * ( s ) = p = 1 N L * ( s ) H 1 p * ( s ) u * p ( s ) ,
where y * ( s ) = y * 1 ( s ) ,   L * ( s ) = L * 1 ( s ) .

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 v ( s ) is a vector with elements 1 W e ( s ) u k z ( s ) , k = 1 , , r .
Instead of (4)–(6), we now consider the equations
u ( s ) = G ( s ) W * ( s ) x * ( s ) + G ( s ) v ( s ) ,
x ( s ) = C ( s ) W * ( s ) x * ( s ) + C ( s ) v ( s ) ,
y ( s ) = L ( s ) W * ( s ) x * ( s ) + L ( s ) v ( s ) .
Determining, as above, the extended vectors and matrices (14), we find
x ^ * ( s ) = C ^ * ( s ) W ^ * ( s ) x ^ * ( s ) + [ C ^ ( s ) v ^ ( s ) ] * ,
x ^ * ( s ) = ( I r N + C ^ * ( s ) W ^ * ( s ) Π ) 1 [ C ^ ( s ) v ^ ( s ) ] * ,
y ^ * ( s ) = L ^ * ( s ) W ^ * ( s ) Π ( I r N + C ^ * ( s ) W ^ * ( s ) Π ) 1 [ C ^ ( s ) v ^ ( s ) ] * + [ L ^ ( s ) v ^ ( s ) ] * ,
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 y ( k T ) , k = 0 , 1 , allows, of course, to isolate the sequence y ( k N T ) , k = 0 , 1 , . At the same time, if the last sequence is determined from the discrete Laplace transformation y N T * ( s ) , 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 N T —a single-rate equivalent model—an impulse stabilization system with a sampling period N T , where, as before, N is the least common multiple of numbers n 1 , , n r ( N = ν 1 n 1 , N = ν 2 n 2 , , N = ν r n r ) .
The starting point for building an N T —single-rate equivalent model is the [29]
y T * ( s ) = k = 1 n e ( k 1 ) T s [ y ( s ) e ( k 1 ) T s ] n T * .
Taking into account T i = N T / ν i , i = 1 , , r ,
x i T i * ( s ) = ν = 1 ν i e ( ν 1 ) T i s [ x i ( s ) e ( ν 1 ) T i s ] N T * = ν = 1 ν i e ν i [ x i ( s ) e ν i + ] N T * ,
where the numbers e ν i ± = e ± ( ν 1 ) T i s , ν = 1 , , ν i ; i = 1 , , r .
Equation (25) can be written as
x i T i * ( s ) = ν = 1 N i e ν i x ν i N T * ( s ) = E i ( s ) x ^ i N T * ( s ) ,
where E i ( s ) = [ e 1 i , e 2 i , , e ν i i ] , x ^ i ( s ) = [ x 1 i ( s ) x 2 i ( s ) x ν i i ( s ) ] T , i = 1 , , r , x ν i ( s ) = e ν i + x i ( s ) .
Using Equation (25) and the block matrix E ( s ) = d i a g [ E 1 ( s ) E 2 ( s ) E r ( s ) ] with dimension r × ( ν 1 + ν 2 + + ν r ) for a vector x * ( s ) = [ x 1 T 1 * ( s ) x 2 T 2 * ( s ) x r T r * ( s ) ] T , we can write x * ( s ) = E ( s ) x ^ N T * ( s ) and
x ( s ) = C ( s ) W * ( s ) E ( s ) x ^ N T * ( s ) + C ( s ) v * ( s ) ,
We now introduce the matrix E i + ( s ) = [ e 1 i + ( s ) e 2 i + ( s ) e ν i i + ( s ) ] T and block matrix E ( s ) = d i a g [ E 1 + ( s ) , E 2 + ( s ) , , E r + ( s ) ] with a dimension ( ν 1 + ν 2 + + ν r ) × r that allows writing E + ( s ) x ( s ) = x ^ ( s ) .
Therefore, multiplying (27) on the left by the matrix E + ( s ) , we obtain
x ^ ( s ) = E + ( s ) C ( s ) W * ( s ) E ( s ) x ^ N T * ( s ) + E + ( s ) C ( s ) v * ( s ) .
Introducing, for convenience, C + ( s ) = E + ( s ) C ( s ) , W * ( s ) = W * ( s ) E ( s ) and W N T * ( s ) = [ C + ( s ) W * ( s ) ] N T * , by Equation (28), we find
x ^ N T * ( s ) = [ I σ + W N T * ( s ) ] 1 [ C + ( s ) v * ( s ) ] N T * ,
where σ = ν 1 + ν 2 + + ν r .
We also find
y N T * ( s ) = [ L ( s ) W * ( s ) ] N T * [ I σ + W N T * ( s ) ] 1 [ C + ( s ) v * ( s ) ] N T * + [ L ( s ) v * ( s ) ] N T * .
The relations (29) and (30) define the equivalent single-rate model of the closed-loop drone stabilization system with a sampling period N T .
An alternative approach to constructing an equivalent single-rate model of the closed-loop drone stabilization system with a sampling period N T 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 n 1 = n 2 = = n r = 1 and N = 1 . The block vectors and matrices introduced above degenerate (consist of one block element). Equation (18) takes the form
y * ( s ) = L * ( s ) ( I m + W * ( s ) C * ( s ) ) 1 v * ( s ) .
Here, the matrix Π = I r 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 n 1 , , n r that determine the relationship between the sampling periods and N —is the least common multiple of the indicated numbers.
As can be seen from Equations (13) and (15), only numbers 1 n i , i = 1 , , r 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 ν p = d i a g [ π ν p 1 , π ν p 2 , , π ν p r ] . To do this, consider the equation
p = 1 N π ν p i x i * p = 1 n i k = 1 n i x i * ν + ( k 1 ) ν i ,
which takes place by Equations (13) and (15) with the notation
x i * p ( s ) = x i * ( s + ( p 1 ) 2 π j N T ) , i = 1 , , r .
Considering Equation (32), one should take into account the periodicity relation
x i * N + p = x i * p , p = 1 , 2 , ,
as property (14) of functions x i * ( s ) , i = 1 , , r .
For this purpose, consider the right-hand side of relation (32) and the set of numbers
p ( k ) = ν + ( k 1 ) ν i , k = 1 , , n i ,
with fixed indices ν , i .
The numbers of this set have the property
p ( k ) > p ( k ) , k > k .
p max be the largest of the numbers p ( k , ν , i ) = ν + ( k 1 ) ν i , k 1 , n i ¯ ; ν 1 , N ¯ ; i 1 , r ¯ . This number has the form
p max = max i 1 , r ¯ ( N + ( n i 1 ) ν i ) = N + N ν ¯ = N + γ ¯ .
where ν ¯ is the smallest of the numbers ν i , i 1 , r ¯ .
Due to the N = n i ν i , at least one of the numbers ν i , i 1 , r ¯ is less then N (at least one of the numbers n i , i 1 , r ¯ is greater than (1), which means γ ¯ < N .
As a result, all numbers of form (33) exceeding N are numbers of the form N + γ , 0 γ N 1 . Now, let, for arbitrarily taken k , k , k , such that k > k > k , we have
p ( k ) = ν + ( k 1 ) ν i = γ N ,
p ( k ) = ν + ( k 1 ) ν i = N + γ
p ( k ) = ν + ( k 1 ) ν i = N + γ
Subtracting, we find that
γ γ = ( k k ) ν i , γ γ = ( k k ) ν i + N = ( n i + k k ) ν i .
Hence,
γ < γ < γ .
Because of the above, from (35), we obtain the following rule for calculating the elements π ν l i , l = 1 , , N of the matrix Π of quantization densities for fixed indices ν , i
π ν l i = { 1 / n i , l = l ( k ) = p ( k ) , p ( k ) N 1 / n i , l = l ( k ) = p ( k ) N , p ( k ) > N 0 l l ( k ) ,
where the values p ( k ) , k = 1 , , n i 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 N , some of them are marked with an asterisk —*. The second line contains values δ ( p ) 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
p ( k ) = ν + ( k 1 ) ν i , k = 1 , , n i .
If p ( k ) N , then the marked number p * = p ( k ) , if
p ( k ) = N + γ ( k ) ,
then p * = γ k .
In this way, the bottom line defines the values δ ( p ) , l = 1 , , N , and
π ν p i = 1 n i δ ( p ) , p = 1 , , N .
Thus, this method determines the entire row of the matrix Π , i.e., diagonal elements with fixed indices ν , i of all block matrices equal to either 1 n i or 0.
Adopting ν = 1 and applying the algorithm n times for i = 1 , , r , we finally find the first row of the block representation of the density matrix Π
Π 11 , Π 12 , , Π 1 N 1 , Π 1 N .
This row defines the entire matrix Π using a sequential circular permutation. For example, the second line will be
Π 1 N , Π 11 , , Π 1 N 2 , Π 1 N 1 ,   etc .
For example, for N = 3 , the entire matrix has the form
Π = Π 11 Π 12 Π 13 Π 13 Π 11 Π 12 Π 12 Π 13 Π 11 .
Note the possibility of forming the given blocks of the matrix Π . For example, we are talking about a block Π ν p . Diagonal elements π ν p i , i = 1 , , r should be defined with, given ν , p . Then:
  • If k 1 , n i ¯ exists, satisfying one of the equations
    ν + ( k 1 ) ν i = p , ν + ( k 1 ) ν i N = p ,
    then π ν p i = 1 n i , and
  • If k 1 , n i ¯ does not exist, then π ν p i = 0 .
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: W 0 ( s ) = Φ S T ( s ) W δ E γ ( s ) , x ( s ) = [ x 1 ( s ) x 2 ( s ) ] T , y ( s ) = γ , E ( s ) = [ W 1 ( s ) W 2 ( s ) ] T , B ( s ) = 0 , W * ( s ) = [ W 1 T 1 * ( s ) W 2 T 3 * ( s ) ] , G ( s ) = W e ( s ) , L ( s ) = W e ( s ) W 0 ( s ) , C ( s ) = [ W e ( s ) W 3 ( s ) W 0 ( s ) W e ( s ) W 0 ( s ) ] .
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
y ( s ) = W e ( s ) W 0 ( s ) [ W 1 * ( s ) W 2 * ( s ) ] [ x 1 * ( s ) x 2 * ( s ) ] + W e ( s ) W 0 ( s ) u * ( s ) ,
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 r = 2 , n 1 = 1 , n 2 = 3 (defined above). Using the above equations, we obtain N = 3 , ν 1 = 3 , ν 2 = 1 .
For block Π 11 , we have the equations
i = 1 1 + ( k 1 ) 3 = 1 , 1 + ( k 1 ) 3 3 = 1 , k 1 ¯ ; i = 2 1 + ( k 1 ) 1 = 1 , 1 + ( k 1 ) 1 3 = 1 , k 1 , 3 ¯ .
Therefore, π 11 1 = 1 , π 11 2 = 1 3 .
For block Π 12 :
i = 1 1 + ( k 1 ) 3 = 2 , 1 + ( k 1 ) 3 3 = 2 , k 1 ¯ ; i = 2 1 + ( k 1 ) 1 = 2 , 1 + ( k 1 ) 1 3 = 2 , k 1 , 3 ¯ .
Therefore, π 12 1 = 0 , π 12 2 = 1 3 .
For block Π 13
i = 1 1 + ( k 1 ) 3 = 3 , 1 + ( k 1 ) 3 3 = 3 , k 1 ¯ ; i = 2 1 + ( k 1 ) 1 = 3 , 1 + ( k 1 ) 1 3 = 3 , k 1 , 3 ¯ .
Therefore, π 13 1 = 0 , π 12 2 = 1 3 .
The entire density matrix has the form
Π = [ 1 0 0 0 0 0 0 1 3 0 1 3 0 1 3 0 0 1 0 0 0 0 1 3 0 1 3 0 1 3 0 0 0 0 0 0 0 1 3 0 1 3 0 1 3 ] .
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., n 1 = n 2 = = n r = n , N = n and ν 1 = ν 2 = = ν r = 1 , then the matrix Π turns out to be maximally filled—all of its block elements Π ν p , ν , p = 1 , , n turn out to be equal
Π ν p = 1 n I r .
where I r —identity r x r matrix. Hence, π ν p i = 1 n i , i 1 , r ¯ .
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:
x ^ * ( s ) = [ x * 1 ( s ) x * 2 ( s ) x * 3 ( s ) ]   with   dimension   3   ×   1 ;   v ^ ( s ) = [ v 1 ( s ) v 2 ( s ) v 3 ( s ) ]   with   dimension   3   ×   1 ;
C ^ ( s ) = [ C 1 ( s ) 0 C 2 ( s ) 0 C 3 ( s ) ]   with   dimension   6   ×   3 ,  
where the dimension of the matrix block C i ( s ) —2 × 1;
W ^ ( s ) = [ W 1 ( s ) 0 W 2 ( s ) 0 W 3 ( s ) ]   with   dimension   3   ×   6 ;
where the dimension of the matrix block W i ( s ) —1 × 2;
L ^ ( s ) = [ L 1 ( s ) 0 L 2 ( s ) 0 L 3 ( s ) ]   with   dimension   3   ×   1 ;  
y ^ ( s ) = [ y 1 ( s ) y 2 ( s ) y 3 ( s ) ]   with   dimension   3   ×   1
.
Knowing the indicated matrices, we get
y ^ * ( s ) = L ^ * ( s ) ( I m N + W ^ * ( s ) C ^ * ( s ) ) 1 v ^ * ( s ) ,
or
[ y * 1 ( s ) y * 2 ( s ) y * 3 ( s ) ] = [ L 1 ( s ) 0 0 0 L 2 ( s ) 0 0 0 L 3 ( s ) ] × × ( [ 1 0 0 0 1 0 0 0 1 ] + [ W 1 T 1 * 1 W 2 T 2 * 1 0 0 0 0 0 0 W 1 T 1 * 2 W 2 T 2 * 2 0 0 0 0 0 0 W 1 T 1 * 3 W 2 T 2 * 3 ] [ 1 0 0 0 0 0 0 1 3 0 1 3 0 1 3 0 0 1 0 0 0 0 1 3 0 1 3 0 1 3 0 0 0 0 0 0 0 1 3 0 1 3 0 1 3 ] × × [ C 1 1 * 0 0 C 2 1 * 0 0 0 C 1 2 * 0 0 C 2 2 * 0 0 0 C 1 3 * 0 0 C 2 3 * ] ) 1 [ u 1 * u 2 * u 3 * ] .
By separating y * 1 ( s ) , 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 T 1 and T 3 . In the system T 1 = 0.1 , T 3 = 0.3 , W 0 ( s ) = 48.6 s 3 144 s 2 228 s 216 s 5 + 4.28 s 4 25.3 s 3 57.7 s 2 + 6.3 s , W 1 ( s ) = 0.708 , W 2 ( s ) = 0.13 , W 3 ( s ) = s .
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
y ( z ) = 0.09645 z 10 + 0.01635 z 9 + 0.52773 z 8 + 0.1572 z 7 0.78658 z 6 + 0.09531 z 5 + + 0 / 901137 z 4 0.016177 z 3 0.43913 z 2 + 0.005319 z + 0.087335 0.23952 z 10 + 1.107978 z 8 2.284746 z 6 + 2.20171 z 4 1.0785 z 2 + 0.210444 ,
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.
The simulation results are shown in Figure 5, Figure 6, Figure 7 and Figure 8.
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.

5. Conclusions

The input–output model of the system in the domain of discrete Laplace transform is adopted as the initial description in the article. The construction of the matrix of transfer functions of a MIMO system is a difficult task, the solution of which is based on the mathematical apparatus of the theory of polynomial matrices. Considerable attention has been paid to this problem, primarily within the framework of the theory of continuous systems [31,32]. In principle, the methods of this theory can be easily transferred to multidimensional impulse (continuous-discrete) single-rate systems. Therefore, in consideration of multirate continuous-discrete, significant attention is paid to the problem of constructing an equivalent impulse single-rate system, the complete solution of which, as shown, is possible in the case of rationally commensurate quantization periods in the digital circuits of the system.
In all cases, the impulse representations of the outputs of the equivalent T- and NT-models have the form of rational functions of the indicated variables
y T * ( s ) = f ( ψ ( e s T ) , ψ 1 ( e s T 1 ) , , ψ r ( e s T r ) ) , y N T * ( s ) = F ( Ψ ( e s T ) , Ψ 1 ( e s T 1 ) , , Ψ r ( e s T r ) ) .
The solution to this problem fundamentally allows developing methods for analyzing and synthesizing MIMO continuous and discrete systems on a common theoretical basis—the theory of polynomials and rational functions, which are typical elements of the descriptions of these classes of systems.

Author Contributions

Conceptualization, V.K. and A.K.; methodology, V.K.; software, V.A.; validation, V.K., A.K. and V.A.; formal analysis, A.K.; investigation, V.K.; resources, A.K.; writing—original draft preparation, V.K.; writing—review and editing, V.K. and V.A.; visualization, V.K.; supervision, V.K.; project administration, A.K.; funding acquisition, A.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

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

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Patel, Y. The Design and Analysis of Multirate Control Systems. Ph.D. Thesis, University of York, New York, NY, USA, 1991. [Google Scholar]
  2. Tomizuka, M. Multi-rate control for motion control applications. In Proceedings of the the 8th IEEE International Workshop on Advanced Motion Control, Kawasaki, Japan, 25–28 March 2004; pp. 21–29. [Google Scholar]
  3. Sushchenko, O.; Goncharenko, A. Design of Robust Systems for Stabilization of Unmanned Aerial Vehicle Equipment. Int. J. Aerosp. Eng. 2016, 2016, 1–10. [Google Scholar] [CrossRef] [Green Version]
  4. Fahlstrom, P.G.; Gleason, T.J. Introduction to UAV Systems, 4th ed.; Wiley: Hoboken, NJ, USA, 2012; p. 287. [Google Scholar]
  5. Li, Z.; Xue, S.-R.; Yu, X.-H.; Gao, H.-J. Controller Optimization for Multirate Systems Based on Reinforcement Learning. Int. J. Autom. Comput. 2020, 17, 417–427. [Google Scholar] [CrossRef]
  6. Carlos, M.; Vélez, S.; Andrés, A. Multirate control of an unmanned aerial vehicle. WSEAS Trans. Circuits Syst. 2005, 11, 1628–1634. [Google Scholar]
  7. Ju, C.; Son, H.I. A Hybrid Systems-based Hierarchical Control Architecture for Heterogeneous Field Robot Teams (arxiv.org). Available online: https://arxiv.org/abs/2002.08602 (accessed on 2 September 2021).
  8. Kiyashchenko, L.A. The problem of estimation in continuous/discrete dynamic systems. J. Math. Sci. 1996, 82, 3311–3315. [Google Scholar] [CrossRef]
  9. Agranovich, G.A. Models of the State of Multiple-Timing Continuous-Discrete Systems. J. Math. Sci. 2001, 103, 27–30. [Google Scholar] [CrossRef]
  10. Barabanov, A.T.; Agranovich, G.A. Linear models and optimization of continuous-discrete dynamical systems. Dyn. Syst. 1983, 2, 17–24. [Google Scholar]
  11. Berg, M.; Amit, N.; Powell, J. Multirate digital control system design. IEEE Trans. Autom. Control 1988, 33, 1139–1150. [Google Scholar] [CrossRef]
  12. Araki, M.; Yamamoto, K. Multivariable multirate sampled-data systems: State-space description, transfer characteristics, and Nyquist criterion. IEEE Trans. Autom. Control 1986, 31, 145–154. [Google Scholar] [CrossRef]
  13. Coffey, T.C.; Williams, I.J. Stability analysis of multiloop, multirate sampled systems. AIAA J. 1966, 4, 2178–2190. [Google Scholar] [CrossRef]
  14. Carlos, M.; Vélez, S.; Andrés, A. Simulation and multirate control of a mini-helicopter robot. In Proceedings of the 5th WSEAS International Conference on Wavelet Analysis and Multirate Systems, Sofia, Bulgaria, 27–29 October 2005; pp. 17–23. [Google Scholar]
  15. Feuer, A.; Goodwin, G.C. Multirate Control. Sampl. Digit. Signal Process. Control 1996, 477–500. [Google Scholar] [CrossRef]
  16. Lopez-Briones, Y.F.; Sanchez-Rivera, L.M.; Arias-Montano, A. Aerodynamic Analysis for the mathematical model of a dual-system UAV. In Proceedings of the 2020 17th International Conference on Electrical Engineering, Computing Science and Automatic Control (CCE), Mexico City, Mexico, 11–13 November 2020; pp. 1–6. [Google Scholar]
  17. Kramar, V.; Alchakov, V.; Kabanov, A.; Dudnikov, S.; Dmitriev, A. The Design of Optimal Lateral Motion Control of an UAV Using the Linear-Quadratic Optimization Method in the Complex Domain. Int. Rev. Aerosp. Eng. 2020, 13, 217. [Google Scholar] [CrossRef]
  18. Kramar, V.; Kabanov, A.; Dudnikov, S. A Mathematical Model for a Conceptual Design and Analyses of UAV Stabilization Systems. Fluids 2021, 6, 172. [Google Scholar] [CrossRef]
  19. Kranc, G.M. Multirate Sampled Systems; TR T-14/B, Electronics Research Lab, Department of Electrical Engineering; Columbia University: New York, NY, USA, 1956. [Google Scholar]
  20. Raggazzini, J.R.; Franklin, G.F. Sampled–Data Control Systems; McGraw-Hill Book: New York, NY, USA, 2015; p. 348. [Google Scholar]
  21. Kuo, B.C. Analysis and Synthesis of Sampled–Data Control Systems; Prentice–Hall Inc.: Englewood Cliffs, NJ, USA, 1963; p. 327. [Google Scholar]
  22. Jury, E.I. A note on multirate sampled-data systems. IEEE Trans. Automat. Control 1967, 12, 319–320. [Google Scholar] [CrossRef]
  23. Kramar, V. Construction of equivalent matrix model of multivaruable multi-stage control system. In Proceedings of the 2016 2nd International Conference on Industrial Engineering, Applications and Manufacturing (ICIEAM), Chelyabinsk, Russia, 19–20 May 2016; pp. 1–4. [Google Scholar]
  24. Gustafson, C.; Desoer, C.A. CAD methodology for linear multivariable feedback systems based on algebraic theory. In Proceedings of the 22nd IEEE Conference on Decision and Control, San Antonio, TX, USA, 14–16 December 1983. [Google Scholar] [CrossRef]
  25. Sandoval, J.; Salt, J. Applications of a novel tool for analysis and design in the frequency domain of Multirate Sampled-Data Control Systems. In Proceedings of the 2008 5th International Conference on Electrical Engineering, Computing Science and Automatic Control, Mexico City, Mexico, 12–14 November 2008; pp. 68–74. [Google Scholar]
  26. Gantmacher, F.R. The Theory of Matrices; AMS Chelsea Publishing: New York, NY, USA; American Mathematical Society: New York, NY, USA, 2000; p. 660. [Google Scholar]
  27. Graham, A. Matrix Theory and Applications for Scientists and Engineers; Dover Publications: Mineola, NY, USA, 2018; p. 304. [Google Scholar]
  28. Jeffrey, A. Matrix Operations for Engineers and Scientists. An Essential Guidein Linear Algebra; Springer Science Business Media: Berlin, Germany, 2010; p. 314. [Google Scholar]
  29. Jury, E.I. Theory and Application of the Z-Transform Method; John Wiley and Sons Inc.: New York, NY, USA, 1964; p. 295. [Google Scholar]
  30. Gasparyan, O.N.; Vardanyan, N.H. On behaviour of root loci of uniform multivariable feedback control systems. Model. Optim. Control 2007, 1, 3–11. [Google Scholar]
  31. Gasparyan, O.N. Linear and Nonlinear Multivariable Feedback Control: A Classical Approach; John Wiley & Sons: Chichester, UK, 2008; p. 354. [Google Scholar]
  32. Skogestad, S.; Postlethwaite, I. Multivariable Feedback Control—Analysis and Design; Wiley: Hoboken, NJ, USA, 2008; p. 590. [Google Scholar]
Figure 1. The aircraft stabilization system.
Figure 1. The aircraft stabilization system.
Drones 05 00129 g001
Figure 2. The block diagram of the complete model of the UAV lateral stabilization system.
Figure 2. The block diagram of the complete model of the UAV lateral stabilization system.
Drones 05 00129 g002
Figure 3. A modified roll stabilization contour.
Figure 3. A modified roll stabilization contour.
Drones 05 00129 g003
Figure 4. A Simulink calculation scheme.
Figure 4. A Simulink calculation scheme.
Drones 05 00129 g004
Figure 5. Comparison of continuous (red) and multirate (blue) systems.
Figure 5. Comparison of continuous (red) and multirate (blue) systems.
Drones 05 00129 g005
Figure 6. Comparison of single-rate (red) and multirate (blue) systems.
Figure 6. Comparison of single-rate (red) and multirate (blue) systems.
Drones 05 00129 g006
Figure 7. Unstable single-rate system.
Figure 7. Unstable single-rate system.
Drones 05 00129 g007
Figure 8. Comparison of equivalent and multirate systems.
Figure 8. Comparison of equivalent and multirate systems.
Drones 05 00129 g008
Table 1. The table to construct the matrix Π .
Table 1. The table to construct the matrix Π .
p12 * * N
δ ( p ) …010…010…
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kramar, V.; Kabanov, A.; Alchakov, V. Multiloop Multirate Continuous-Discrete Drone Stabilization System: An Equivalent Single-Rate Model. Drones 2021, 5, 129. https://doi.org/10.3390/drones5040129

AMA Style

Kramar V, Kabanov A, Alchakov V. Multiloop Multirate Continuous-Discrete Drone Stabilization System: An Equivalent Single-Rate Model. Drones. 2021; 5(4):129. https://doi.org/10.3390/drones5040129

Chicago/Turabian Style

Kramar, Vadim, Aleksey Kabanov, and Vasiliy Alchakov. 2021. "Multiloop Multirate Continuous-Discrete Drone Stabilization System: An Equivalent Single-Rate Model" Drones 5, no. 4: 129. https://doi.org/10.3390/drones5040129

APA Style

Kramar, V., Kabanov, A., & Alchakov, V. (2021). Multiloop Multirate Continuous-Discrete Drone Stabilization System: An Equivalent Single-Rate Model. Drones, 5(4), 129. https://doi.org/10.3390/drones5040129

Article Metrics

Back to TopTop