Next Article in Journal
Differential Transform Method as an Effective Tool for Investigating Fractional Dynamical Systems
Next Article in Special Issue
Determination of the Function of the Course of the Static Property of PAMs as Actuators in Industrial Robotics
Previous Article in Journal
A Combined Non-Destructive and Micro-Destructive Approach to Solving the Forensic Problems in the Field of Cultural Heritage: Two Case Studies
Previous Article in Special Issue
Chimney Sweeping Robot Based on a Pneumatic Actuator
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Methods of Pre-Identification of TITO Systems

1
Faculty of Mechanical Engineering, University of Zilina, Univerzitna 8215/1, 01026 Zilina, Slovakia
2
Faculty of Applied Informatics, Tomas Bata University in Zlin, Nad Stráněmi 4511, 760 05 Zlin, Czech Republic
3
VIPO, a. s., Gen. Svobodu 1069/4, 95801 Partizanske, Slovakia
*
Author to whom correspondence should be addressed.
Appl. Sci. 2021, 11(15), 6954; https://doi.org/10.3390/app11156954
Submission received: 21 April 2021 / Revised: 9 July 2021 / Accepted: 23 July 2021 / Published: 28 July 2021

Abstract

:
The content of this article is the presentation of methods used to identify systems before actual control, namely decentralized control of systems with Two Inputs, Two Outputs (TITO) and with two interactions. First, theoretical assumptions and reasons for using these methods are given. Subsequently, two methods for systems identification are described. At the end of this article, these specific methods are presented as the pre-identification of the chosen example. The Introduction part of the paper deals with the description of decentralized control, adaptive control, decentralized control in robotics and problem formulation (fixing the identification time at the existing decentralized self-tuning controller at the beginning of control and at the beginning of any set-point change) with the goal of a new method of identification. The Materials and methods section describes the used decentralized control method, recursive identification using approximation polynomials and least-squares with directional forgetting, recursive instrumental variable, self-tuning controller and suboptimal quadratic tracking controller, so all methods described in the section are those ones that already exist. Another section, named Assumptions, newly formulates the necessary background information, such as decentralized controllability and the system model, for the new identification method formulated in Pre-identification section. This section is followed by a section showing the results obtained by simulations and in real-time on a Coupled Drives model in the laboratory.

1. Introduction

Most processes in practice are processes that have multiple inputs and multiple outputs, and these are influenced by interactions. These systems can be controlled by a centralized or decentralized controller. The main advantages of decentralization include simplifying the overall task by dividing it into a set of sub-tasks. Decentralized control is very often used in practice. Decentralized charge control of electric vehicles is a nice application. There was introduced a fully decentralized and participatory learning mechanism for privacy-preserving coordinated charging control of electric vehicles that regulates three Smart Grid socio-technical aspects: (i) reliability, (ii) discomfort and (iii) fairness [1]. Another good application of decentralized control is a quadrocopter, namely outdoor flocking of quadcopter drones with decentralized model predictive control [2]. From the theoretical point of view, there was proposed a decentralized explicit (closed-form) iterative formula that solves convex programming problems with linear equality constraints and interval bounds on the decision variables [3], or decentralized control problem for non-affine large-scale systems with nonaffine functions possibly being discontinuous [4], or decentralized adaptive tracking control strategy consisting of a steady-state controller and modified optimal feedback controller. Design parameters-dependent feasibility conditions were formulated by using Lyapunov theory to guarantee the existence of the proposed decentralized control scheme [5]. Decentralized voltage control is another example of a decentralized strategy. It includes network partitioning strategy for the optimal voltage control of Active Distribution Networks actuated by means of a limited number of Distributed Energy Storage Systems [6]. Another paper is concerned with the problem of decentralized event-triggered dynamic output feedback control for large-scale systems with unknown interconnections. By using a modified cyclic-small gain condition and introducing a free-matrix-based integral inequality, a sufficient condition was derived to ensure that the overall closed-loop system is asymptotically stable with a prescribed H∞ performance [7]. There was also implemented fuzzy decentralized control, for example an adaptive fuzzy decentralized control approach for a class of nonlinear systems with unknown control direction and different performance constraints. In the control design, the different performance constraints, that were the prescribed performance error constraints for some subsystems and the asymmetric time-varying output constraints for the others, could be unified as one form by selecting proper performance functions [8]. Adaptive control is another area that has expanded the use of decentralized control, and here are at least a few such examples. A minimal-neural-networks-based design approach was presented for the decentralized output feedback tracking of uncertain interconnected strict-feedback nonlinear systems with unknown time varying delayed interactions unmatched in control inputs [9]. The decentralized output-feedback adaptive backstepping control scheme was also proposed for a class of interconnected nonlinear systems with unknown actuator failures by introducing a kind of high gain K-filters [10], or decentralized output-feedback adaptive control scheme for a class of interconnected nonlinear systems with input quantization. Both logarithmic quantizers and improved hysteretic quantizers were studied, and a linear time-varying model was introduced to handle the difficulty caused by quantization [11]. A decentralized adaptive backstepping control scheme was also proposed for a class of interconnected systems with nonlinear multisource disturbances and actuator faults. The nonlinear multisource disturbances comprised two parts: one was the time-varying parameterized uncertainty; the other was the dynamic unexpected signal formulated by a nonlinear exogenous system [12]. Additionally, the problem of decentralized adaptive backstepping control for a class of large-scale stochastic nonlinear time-delay systems with asymmetric saturation actuators and output constraints was also solved [13]. It can be mentioned that a backstepping-based robust decentralized adaptive neural H ∞ tracking control method was addressed for a class of large-scale strict feedback nonlinear systems with uncertain disturbances [14]. Decentralized control was implemented in robotics. The example where a discrete-time decentralized neural identification and control for large-scale uncertain nonlinear systems at a two degree of freedom planar robot was implemented can be mentioned [15], or the work that investigated the use of a decentralized control system for suppressing the vibration of a multi-link flexible robotic manipulator using embedded smart piezoelectric transducers [16]. Decentralized motion coordination algorithms were proposed for the robots so that they collectively moved in a rectangular lattice pattern from any initial position [17]. Mobile robot formations differ in accordance to the mission, environment and robot abilities. In the case of decentralized control, the ability to achieve the shapes of these formations has to be built in the controllers of each autonomous robot [18]. A decentralized control algorithm for the robots to accomplish the sweep coverage was also proposed. The sweep coverage was achieved by coordinating the robots to move along a given path that was unknown to the vehicles a priori [19].
During the simulation experiments in the real-time laboratory in recent years, we reveal the fact that some time is needed to get the appropriate behaviour of control when the self-tuning controller is used. This time depends on the type of the controlled system. This is the known problem of self-tuning control because the controller needs the adequate model of the system. This is one problem that we solve by a new approach described in this paper, by the method named as pre-identification. Another problem comes from the fact that we used the decentralized controllers for the control of multivariable systems. If one set-point changes its value, it influences all other main subsystems by interactions and therefore the model of subsystems changes. This could be also fixed by a self-tuning controller but some time is needed to obtain the stable model. Therefore, by the new method described in this paper, named as pre-identification, we also solved this problem.

2. Materials and Methods

2.1. Decentralized Control

Using the decentralized approach, the control is divided into a set of sub-tasks that are matched by simple controllers. These partial tasks will then give us the overall course of control. The main advantages of decentralized control are primarily that a more complex system is divided into a set of simple tasks and the resulting controller is more flexible [20].
A special example of multidimensional systems is a system with two inputs and two outputs. This can be realized by the so-called P structure, see Figure 1. In this case, the inputs to the systems describing the interactions are the values of the action signals of the SISO controllers and their outputs are added to the opposite outputs of the main diagonal systems. From this figure, we get the transfer function equations of the model in the form
G S 1 = Y 1 U 1 = G S 11 G S 21 G S 12 G R 22 1 + G S 22 G R 22
G S 2 = Y 2 U 2 = G S 22 G S 21 G S 12 G R 11 1 + G S 11 G R 11

2.2. Recursive Identification Using Approximation Polynomials

A prerequisite for good control is the most accurate description of the regulated system. Identification is the procedure by which the mathematical model of a system is obtained. The beginnings of identification based on continuous models date back to the middle of the 20th century. For continuous-time identification, the identified model is in the form of the differential equations. Differential equations contain expressions with derivatives over time that are not measurable. It is possible to replace the segment by an approximation polynomial whose derivatives can be calculated analytically in advance and then calculated numerically, see Figure 2. This approach was for example used by Perutka [20].

2.3. Least Squares Method with Exponential Forgetting

The estimation of model parameters is computed as
θ ^ ( k ) = θ ^ ( k 1 ) + K ( k ) e ^ ( k )
The gain vector is calculated as
K ( k ) = C ( k 1 ) ϕ ( k ) 1 + ϕ T ( k ) C ( k 1 ) ϕ ( k )
and covariance matrix
C ( k ) = C ( k 1 ) C ( k 1 ) ϕ ( k ) ϕ T ( k ) C ( k 1 ) 1 + ϕ T ( k ) C ( k 1 ) ϕ ( k )
The following applies to the calculation of the prediction error
e ^ ( k ) = y ( k ) ϕ T ( k ) θ ^ ( k 1 )
In the case of exponential forgetting, the criterion of identification is
J = i = k 0 k ( φ k i e ( i ) ) 2
where the exponential forgetting factor is chosen in the range of 0 to 1, the most common near 1.
If
ϕ T ( k ) C ( k 1 ) ϕ ( k ) > 0
Then
C ( k ) = C ( k 1 ) C ( k 1 ) ϕ ( k ) ϕ T ( k ) C ( k 1 ) η 1 + ϕ T ( k ) C ( k 1 ) ϕ ( k )  
where
η ( k ) = φ ( k ) 1 φ ( k ) ξ ( k )
If
ϕ T ( k ) C ( k 1 ) ϕ ( k ) = 0
Then
C ( k ) = C ( k 1 )
Furthermore
φ ( k ) = { 1 + ( 1 + ρ ) [ ln ( 1 + ξ ( k 1 ) ) ] + [ ( v ( k 1 ) + 1 ) η ( k 1 ) 1 + ξ ( k 1 ) + η ( k 1 ) 1 ] ξ ( k 1 ) 1 + ξ ( k 1 ) } 1
η ( k ) = e 2 ( k ) λ ( k )
v ( k ) = φ ( k ) [ v ( k 1 ) + 1 ]
λ ( k ) = φ ( k ) [ λ ( k 1 ) + e 2 ( k ) 1 + ξ ( k 1 ) ]
ξ ( k ) = ϕ T ( k ) C ( k 1 ) ϕ ( k )
The parameters estimation vector is in the form
θ ^ T ( k ) = ( a ^ 0 , a ^ 1 ,     ,   a ^ deg ( a ) ,   b ^ 0 , b ^ 1 ,     ,   b ^ deg ( b ) ,   d )
and regressor
ϕ T ( k ) = ( y ( t k ) ,   , y L ( n 1 ) ( t k ) ,   u ( t k ) ,   ,   u L ( m ) ( t k ) ,   1 )

2.4. Self-Tuning Controller

The main reason for using adaptive control is that the systems change over time or the characteristics of the controlled system are unknown. The basic principle of adaptive systems is to change the characteristics of the controller based on the characteristics of the controlled process [21]. The general scheme of the self-tuning controller is shown in Figure 3.

2.5. Suboptimal Linear Quadratic Tracking Controller

The method was introduced by Dostál [22]. If the system of Figure 4 is considered
Let us minimize a quadratic functional with two penalty constants
J = 0 { μ e 2 ( t ) + φ u ˜ 2 ( t ) } d t
The Laplace image of the set point holds
w ( s ) = h w ( s ) s f w ( s )
It holds for degrees of polynomials
deg ( h w ) deg ( f w ) ,   f w ( 0 ) 0
We calculate stable polynomials g and n as results of spectral factorizations
( a s ) φ a s + b μ b = g g , n n = a a
We solve the following diophantine equation
a s p + b q = g n
Considering the transfer function of the system
G ( s ) = b 0 s 2 + a 1 s + a 0
then the controller is
F ( s ) Q ( s ) = q 2 s 2 + q 1 s + q 0 s ( p 2 s 2 + p 1 s + p 0 )
In this case, the polynomials have the form
g ( s ) = g 3 s 3 + g 2 s 2 + g 1 s + g 0
n ( s ) = s 2 + n 1 s + n 0
and to calculate their coefficients obtained by spectral factorization
g 0 = μ b 0 2
g 1 = 2 g 2 g 0 + φ a 0 2
g 2 = 2 g 3 g 1 + φ ( a 1 2 2 a 0 )
g 3 = φ
n 0 = a 2 0
n 1 = 2 n 0 a 1 2 2 a 0

2.6. Calculation of Derivatives Using Approximation Functions

To calculate the derivatives, we approximate the closest neighborhood for a given time by the approximation function. For example, we will use the Lagrange polynomial in the form
P 2 ( x ) = ( x b ) ( x c ) ( a b ) ( a c ) f ( a ) + ( x a ) ( x c ) ( b a ) ( b c ) f ( b ) + ( x b ) ( x b ) ( c a ) ( c b ) f ( a c )
whose first derivative is
f ( x ) P 2 ( x ) = 2 x ( b + c ) ( a b ) ( a c ) f ( a ) + 2 x ( a + c ) ( b a ) ( b c ) f ( b ) + 2 x ( a + b ) ( c a ) ( c b ) f ( c )
and second derivative is
f / / ( x ) P 2 / / ( x ) = 2 f ( a ) ( a b ) ( a c ) + 2 f ( b ) ( b a ) ( b c ) + 2 f ( c ) ( c a ) ( c b )

2.7. Recursive Instrumental Variable Method

The instrumental variable method is a modification of the least squares method. The least squares method uses the quadratic criterion and the existence of one global minimum. However, a prerequisite for successful least-squares modelling is that the fault is represented by white noise with zero mean value. The instrumental variable method does not make it possible to determine the noise properties, but is based on weaker assumptions than the least squares method. The identification proceeds according to number 5. As with the least squares method, the method of the instrumental variable can also be formulated recursively [22,23,24,25]. The parameter estimation vector has the form
θ ^ T ( k ) = ( a ^ 0 ,   a ^ 1 ,   ,   a ^ deg ( a ) ,   b ^ 0 ,   b ^ 1 ,     ,   b ^ deg ( b ) ,   d )
and data vector
θ ^ T ( k ) = ( y ( t k ) ,   y L ( n 1 ) ( t k ) , u ( t k ) ,   ,   u L ( m ) ( t k ) , 1
The gain vector is calculated by relation
L ( k ) = C ( k 1 ) z ( k ) 1 + ϕ T ( k ) C ( k 1 ) z ( k 1 )
In addition to the data vector, it is necessary to know the covariance matrix
C ( k ) = C ( k 1 ) C ( k 1 ) z ( k ) ϕ T ( k ) C ( k 1 ) 1 + ϕ T ( k ) C ( k 1 ) z ( k )
and instrument vector
z ( k ) = ( u ( t k ) ,   u ( t k 1 ) ,     ,   u ( t k n m ) ,    
which we choose as a set of delayed inputs. The prediction error is calculated by
e ^ ( k ) = y ( k ) ϕ T ( k ) θ ^ ( k 1 )
and estimating the parameters according to
θ ^ ( k ) = θ ^ ( k 1 ) + L ( k ) e ^ ( k )

3. Assumptions

3.1. Decentralized Controllability

Assume the existence of a stable minimum phase Linear Time Invariant (LTI) in time of a continuous multidimensional system of the dimension N × N. Its Laplace transformation S ( s ) : S ( t ) > S ( s ) , which we call the transformed system is in the form
S ( s ) = ( S 11 ( s ) S 12 ( s ) S 1 N ( s ) S N 1 ( s ) S N 2 ( s ) S N N ( s ) )
where S i j ( s ) ,           i = 1 , 2 ,   , N ,           j = 1 , 2 ,    is Laplace transformation of the ij-th subset of Sij(t) of the transformed system S(s). The transformed system S(s) je is decentrally controllable only when its main diagonal is dominant.

3.2. System Model and Shape of Reference Signal

Suppose there exists a system S(t) and a transformed system S(s) as described above. Then we formulate a model created for the purpose of decentralized control, which we call M(t), and its Laplace transformation M(s). Consider M(s) as a stable minimum phase linear time invariant multivariate diagonal matrix in the form
M ( s ) = ( M ( s ) 0 0 0 M ( s ) 0 0 0 M N ( s ) )
where M i ( s ) ,         i = 1 , 2 ,   ,   N  is the Laplace transformation of the i-th submodule Mi(t) of the model M(t)of the transformed system S(s). This assumes minimal impact of extra-diagonal transmissions, which is important because of the deployment of a decentralized controller. Simplification of the N-dimensional system to N-dimensional systems is simplified.
Furthermore, suppose that the Laplace transformation of the reference signal vector r(s) is always in the form
r ( s ) = ( R 1 ( s ) ,   R 2 ( s ) ,   ,   R N ( s ) )
where Ri(s), i = 1, 2..., N is the i-th Laplace reference signal of Laplace transformation of the reference signal vector r(s) and has the form
R i ( s ) = h i s
where h i     R ,   i = 1 ,   2 ,   ,   N , is the i-th part by constant function, i.e., reference signal, which is only a combination of p-l step changes of the signal of its different constant values
h i = { j i 1   f o r   t     0 ,   t i 1 ) j i 2   f o r   t     t i 1 ,   t i 2 ) j i p   f o r   t     t i p 1 ,   t i p ) ,   0 < t i 1 < t i 2 < < t i p 1 < t i p
where j i l   R ,         i = 1 ,   2 ,   ,   N ,         l = 1 , 2 ,   ,   p ,   the l-th constant function i-th in parts by the constant function hi, t is the time t i l R + ,         i = 1 , 2 ,   ,   N ,         l = 1 , 2 ,   , p , the l-th moment of the i-th in portions of the constant function hi. This means that each non-zero element of the matrix M(s) has exactly one non-zero element of the vector r(s), i.e., that each partial transmission of the overall system model has a reference signal defined for it. As for the form of the reference signal, it is a constant function in parts. This function is approximated from an arbitrary but predetermined number of p segments of a different but concise value, i.e., it varies over time.

4. Pre-Identification

Consider the validity of the assumption of decentralized controllability, system description and system model. Then, the transformed system S(s) can be viewed as a black box model, and let it be identified by direct and/or indirect time-continuous algorithms. In time, continuous model identification can be realized by following steps: The controller is not connected in the closed circuit. The values of the vector of difference of output quantities and reference signals E(t) are sent to the input of the system S(t). The values of the reference signals are the same and at the same time as those that will be used during regulation.
  • The controller is not connected in the closed circuit. The values of the vector of difference of output quantities and reference signals E(t) are sent to the input of the system S(t). The values of the reference signals are the same and at the same time as those that will be used during regulation.
  • If switching control is considered, each time interval of the control of the system S(s) at which all reference signals have a constant value is identified separately, in so-called Identification Elements (IE).
  • Each identification element is identified several times, each time by a different identification algorithm, and the obtained model can be verified by comparison with the measured data. The obtained model, which is most consistent with the measured data, is then used for control. Let us call this method of Identification More Than One Method (IMTOM).
The system S(s) is completely identified by the above procedure before the actual regulation begins, therefore identification during the regulation is not necessary. This procedure is suitable for processes where the same procedure is repeated many times.

5. Results

5.1. Simulation Results

The verified system is described as
G S ( s ) = ( 3.7 s 2 + 5.2 s + 4.6 0.4 s 2 + 4.4 s + 3.8 0.6 s 2 + 10.6 s + 10.2 8.7 s 2 + 7.4 s + 8.2 )
Since it is the system with two inputs and two outputs, we obtained two suboptimal linear quadratic controllers in the form that was described in the previous section. These controllers had the following penalty constants
μ 1 = 1 ,   φ 1 = 30 , μ 2 = 1 ,   φ 2 = 30
We used our algorithm, pre-identification, at this system and we obtained the following results, see Figure 5, Figure 6, Figure 7 and Figure 8. First, we performed system response on the given set-points depicted in Figure 5. Together with this system response, we obtained the system pre-identified parameters shown in Figure 6., for the first subsystem in the left one and for the second subsystem in the right one. According to these pre-identified parameters we performed the simulation of control which is shown in Figure 7. During the control, the controller parameters were changing, and they are recorded in Figure 8., for the first controller in left one and and for the second controller in the right one.

5.2. Results in Real-Time at Laboratory Model

We verified the presented method in real-time control using MATLAB at CE108 Coupled Drives Apparatus Model [26] which is shown in the Figure 9. The laboratory model CE108 allows solving practical problems of tensioning and speed of material movement in continuous processes. An example is the speed and tension of the thread when rewinding from one spool to another, which must be controlled. This situation is adapted for laboratory measurements so that the elastic band is mounted on three wheels. The lower two wheels are fixed, their speed can be measured and regulated. The third wheel can move (located on a movable arm suspended on a spring) and simulates a workstation along with tension and speed measurements. Two servomotors control the speed of movement and tensioning of the belt. The speed is 0–3000 rpm, which corresponds to a voltage of 0–10 V. Tension measurement is indirect through the angle of the movable arm, from −10 ° to 10 °, which corresponds to a voltage from −10 V to 10 V. Inputs and outputs are located on the front panel of the device; it is the control voltage to the servomotor amplifiers, which are bidirectional, and which are inputs. There are four outputs, the voltage corresponding to the speed of the two servomotors and the voltage corresponding to the tension and the speed of the belt.
Using the pre-identification method and fully implementing interactions in the main plants, we obtained the following mathematical model to be used at control of speed
G S ( s ) = ( 0.78 s 2 + 2.66 s + 1.33 0 0 4.16 s 2 + 11.66 s + 16.66 )
Since it is the system with two inputs and two outputs, we obtained two suboptimal linear quadratic controllers in the form that was described in the previous section. These controllers had the following penalty constants
μ 1 = 1 , φ 1 = 1 , μ 2 = 1 , φ 2 = 0.85
We used our algorithm, pre-identification, at this system and we obtained the following results, see Figure 10, Figure 11, Figure 12 and Figure 13. First, we performed system response on the given set-points depicted in Figure 10. Together with this system response, we obtained the system pre-identified parameters shown in Figure 11, for the first subsystem in the left one and for the second subsystem in the right one. According to these pre-identified parameters we performed the simulation of control which is shown in Figure 12. During the control, the controller parameters were changing, and they are recorded in Figure 13, for the first controller in left one and and for the second controller in the right one.

6. Conclusions

This paper presents the new method of identification named as pre-identification on the theoretical level and subsequently verified it by simulations and in the real-time experiments at Coupled Drives model in the laboratory. The results confirm that the method can be successfully used and future work will focus on the verification of this method on more examples both in simulation and in laboratory conditions. This new method enhances the usage of a decentralized self-tuning controller in the way that it fixes the time the adaptive controller needs to adapt its model.

Author Contributions

Conceptualization, K.P. and I.K.; methodology, K.P. and I.K.; software, I.Z.; validation, V.B., K.P. and V.T.; formal analysis, M.B.; investigation, K.P. and V.T.; resources, K.P.; data curation, K.P. and M.B.; writing—original draft preparation, K.P.; writing—review and editing, I.Z.; visualization, K.P. and I.Z.; supervision, K.P. and M.S.; project administration, K.P. and I.K.; funding acquisition, M.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by by the European Regional Development Fund under the project CEBIA-Tech No. CZ.1.05/2.1.00/03.0089. This research was funded by the Ministry of Education, Science, research and Sport of the Slovak Republic under the project STIMULY MATADOR 1247/2018. Project title: Research and development of modular reconfigurable production systems using Smart Industry principles for automotive with pilot application in MoBearing Line industry.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This work was supported by the European Regional Development Fund under the project CEBIA-Tech No. CZ.1.05/2.1.00/03.0089. This article was made under support of project: STIMULY MATADOR 1247/2018. Project title: Research and development of modular reconfigurable production systems using Smart Industry principles for automotive with pilot application in MoBearing Line industry.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pournas, E.; Jung, S.; Yadhunathan, S.; Zhang, H.; Fang, X. Socio-technical smart grid optimization via decentralized charge control of electric vehicles. Appl. Soft Comput. J. 2019, 82, 10573. [Google Scholar]
  2. Yuan, Q.; Zhan, J.; Li, X. Outdoor flocking of quadcopter drones with decentralized model predictive control. ISA Trans. 2017, 71, 84–92. [Google Scholar] [CrossRef] [PubMed]
  3. Blanchini, F.; Casagrande, D.; Fabiani, F.; Giordano, G.; Pesenti, R. Network-decentralized optimization and control: An explicit saturated solution. Automatica 2019, 103, 379–389. [Google Scholar] [CrossRef]
  4. Wang, Y.; Hu, J.; Zheng, Y. Improved decentralized prescribed performance control for non-affine large-scale systems with uncertain actuator nonlinearity. J. Frankl. Inst. 2019, 356, 7091–7111. [Google Scholar] [CrossRef]
  5. Qu, Q.; Zhang, H.; Feng, T.; Jiang, H. Decentralized adaptive tracking control scheme for nonlinear large-scale interconnected systems via adaptive dynamic programming. Neurocomputing 2017, 225, 1–10. [Google Scholar] [CrossRef]
  6. Bahramipanah, M.; Cherkaoui, R.; Paolone, M. Decentralized voltage control of clustered active distribution network by means of energy storage systems. Electr. Power Syst. Res. 2016, 136, 370–382. [Google Scholar] [CrossRef]
  7. Liu, D.; Yang, G.H. Decentralized event-triggered output feedback control for a class of interconnected large-scale systems. ISA Trans. 2019, 93, 156–164. [Google Scholar] [CrossRef] [PubMed]
  8. Shi, W.; Yan, F.; Li, B. Adaptive fuzzy decentralized control for a class of nonlinear systems with different performance constraints. Fuzzy Sets Syst. 2019, 374, 1–22. [Google Scholar] [CrossRef]
  9. Choi, Y.H.; Yoo, S.J. Decentralized adaptive output-feedback control of interconnected nonlinear time-delay systems using minimal neural networks. J. Frankl. Inst. 2018, 355, 81–105. [Google Scholar] [CrossRef]
  10. Wang, C.; Wen, C.; Guo, L. Decentralized output-feedback adaptive control for a class of interconnected nonlinear systems with unknown actuator failures. Automatica 2016, 71, 187–196. [Google Scholar] [CrossRef]
  11. Wang, C.; Wen, C.; Lin, Y.; Wang, W. Decentralized adaptive tracking control for a class of interconnected nonlinear systems with input quantization. Automatica 2017, 81, 359–368. [Google Scholar] [CrossRef]
  12. Wang, Z.; Zhang, B.; Yuan, J. Decentralized adaptive fault tolerant control for a class of interconnected systems with nonlinear multisource disturbances. J. Frankl. Inst. 2018, 355, 4493–4514. [Google Scholar] [CrossRef]
  13. Si, W.; Dong, X.; Yang, F. Decentralized adaptive neural control for interconnected stochastic nonlinear delay-time systems with asymmetric saturation actuators and output constraints. J. Frankl. Inst. 2018, 355, 54–80. [Google Scholar] [CrossRef]
  14. Li, X.; Liu, X. Backstepping-based decentralized adaptive neural H tracking control for a class of large-scale nonlinear interconnected systems. J. Frankl. Inst. 2018, 355, 4533–4552. [Google Scholar] [CrossRef]
  15. Tellez, F.O.; Loukianov, A.G.; Sanchez, E.N.; Corrochano, J.B. Decentralized neural identification and control or uncertain nonlinear systems: Application to planar robot. J. Frankl. Inst. 2010, 347, 1015–1034. [Google Scholar] [CrossRef]
  16. Halim, D.; Luo, X.; Trivailo, P.M. Decentralized vibration control of a multi-link flexible robotic manipulator using smart piezoelectric transducers. Acta Astronaut. 2014, 104, 186–196. [Google Scholar] [CrossRef]
  17. Cheng, T.M.; Savkin, A.V. Decentralized Control of Multi-robot Systems for Rectangular Aggregation. IFAC Proc. Vol. 2011, 44, 11574–11579. [Google Scholar] [CrossRef] [Green Version]
  18. Necsulescu, D.; Pruner, E.; Kim, B.; Sasiadek, J. Decentralized Control of Autonomous Mobile Robots Formations using Velocity Potentials. IFAC Proc. Vol. 2012, 45, 31–36. [Google Scholar] [CrossRef]
  19. Cheng, T.M.; Savkin, A.; Javed, F. Decentralized Control of a group of mobile robots for deployment in sweep coverage. Robot. Auton. Syst. 2011, 59, 497–507. [Google Scholar] [CrossRef]
  20. Perutka, K. Adaptive LQ control with pre-identification of two tanks laboratory model. In Annals of DAAAM and Proceedings of the International DAAAM Symposium; Danube Adria Association for Automation and Manufacturing: Wienna, Austria, 2009; pp. 439–440. [Google Scholar]
  21. Bobal, V.; Böhm, J.; Fessl, J.; Machacek, J. Digital Self-Tuning Controllers; Springer: London, UK, 2009. [Google Scholar]
  22. Dostal, P.; Bobal, V. The suboptimal tracking problem in linear systems. In Proceedings of the 7th Conference on Control and Automation, Haifa, Israel, 28–30 June 1999; pp. 667–673. [Google Scholar]
  23. Perutka, K. Decentralized Adaptive Control. Ph.D. Thesis, Tomas Bata University iz Zlin, Zlín, Czech Republic, 2007. [Google Scholar]
  24. Velíšek, K.; Holubek, R.; Sobrino, D.R.; Ružarovský, R.; Vetríková, N. Design of a robotized workstation making use of the integration of CAD models and Robotic Simulation software as way of pairing and comparing real and virtual environments. In MATEC, 4th International Conference on Computing and Solutions in Manufacturing Engineering 2016—CoSME´16; EDP Sciences: Youlis, France, 2017. [Google Scholar]
  25. Holubek, R.; Ružarovský, R.; Delgado Sobrino, D.R.; Košťál, P.; Švorc, A.; Velíšek, K. Novel trends in the assembly process as the results of human—The industrial robot collaboration. In MATEC Web of Conferences; EDP Sciences: Youlis, France, 2017; Volume 137. [Google Scholar]
  26. TQ Tecquioment Academia-CE108 Coupled Drives Apparatus. November 2018. Available online: https://www.tecquipment.com/assets/documents/datasheets/CE108-Coupled-Drives-Datasheet.pdf (accessed on 28 July 2021).
Figure 1. Decentralized system control with two inputs and two outputs, the so-called P structure.
Figure 1. Decentralized system control with two inputs and two outputs, the so-called P structure.
Applsci 11 06954 g001
Figure 2. Identification scheme for continuous-time systems.
Figure 2. Identification scheme for continuous-time systems.
Applsci 11 06954 g002
Figure 3. The general scheme of the self-tuning controller.
Figure 3. The general scheme of the self-tuning controller.
Applsci 11 06954 g003
Figure 4. System with feedback controller.
Figure 4. System with feedback controller.
Applsci 11 06954 g004
Figure 5. System response (green and cyan) on the given set-points (blue and red).
Figure 5. System response (green and cyan) on the given set-points (blue and red).
Applsci 11 06954 g005
Figure 6. System pre-identification using recursive intrumental variable—1st subsystem (left), 2nd (right).
Figure 6. System pre-identification using recursive intrumental variable—1st subsystem (left), 2nd (right).
Applsci 11 06954 g006
Figure 7. Output of control.
Figure 7. Output of control.
Applsci 11 06954 g007
Figure 8. Controllers parameters during simulation—1st (left), 2nd (right).
Figure 8. Controllers parameters during simulation—1st (left), 2nd (right).
Applsci 11 06954 g008
Figure 9. Photo of CE108 Coupled Drives Apparatus model connected with PC with MATLAB.
Figure 9. Photo of CE108 Coupled Drives Apparatus model connected with PC with MATLAB.
Applsci 11 06954 g009
Figure 10. System response (green and cyan) on the given set-points (blue and red).
Figure 10. System response (green and cyan) on the given set-points (blue and red).
Applsci 11 06954 g010
Figure 11. System pre-identification using recursive instrumental variable—1st subsystem (left), 2nd (right).
Figure 11. System pre-identification using recursive instrumental variable—1st subsystem (left), 2nd (right).
Applsci 11 06954 g011
Figure 12. Output of control.
Figure 12. Output of control.
Applsci 11 06954 g012
Figure 13. Controllers parameters during simulation—1st (left), 2nd (right).
Figure 13. Controllers parameters during simulation—1st (left), 2nd (right).
Applsci 11 06954 g013
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Saga, M.; Perutka, K.; Kuric, I.; Zajačko, I.; Bulej, V.; Tlach, V.; Bezák, M. Methods of Pre-Identification of TITO Systems. Appl. Sci. 2021, 11, 6954. https://doi.org/10.3390/app11156954

AMA Style

Saga M, Perutka K, Kuric I, Zajačko I, Bulej V, Tlach V, Bezák M. Methods of Pre-Identification of TITO Systems. Applied Sciences. 2021; 11(15):6954. https://doi.org/10.3390/app11156954

Chicago/Turabian Style

Saga, Milan, Karel Perutka, Ivan Kuric, Ivan Zajačko, Vladimír Bulej, Vladimír Tlach, and Martin Bezák. 2021. "Methods of Pre-Identification of TITO Systems" Applied Sciences 11, no. 15: 6954. https://doi.org/10.3390/app11156954

APA Style

Saga, M., Perutka, K., Kuric, I., Zajačko, I., Bulej, V., Tlach, V., & Bezák, M. (2021). Methods of Pre-Identification of TITO Systems. Applied Sciences, 11(15), 6954. https://doi.org/10.3390/app11156954

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

Article Metrics

Back to TopTop