=Paper= {{Paper |id=Vol-3790/paper8 |storemode=property |title=Generalized integro - differentiating controller for mechatronic devices of mobility nodes of humanoid robots |pdfUrl=https://ceur-ws.org/Vol-3790/paper08.pdf |volume=Vol-3790 |authors=Oleksandr Lysenko,Olena Tachinina,Oleksandr Guida,Iryna Alekseeva,Vladyslav Kutiepov |dblpUrl=https://dblp.org/rec/conf/icst2/LysenkoTGAK24 }} ==Generalized integro - differentiating controller for mechatronic devices of mobility nodes of humanoid robots== https://ceur-ws.org/Vol-3790/paper08.pdf
                                Generalized integro - differentiating controller for
                                mechatronic devices of mobility nodes of humanoid
                                robots
                                Oleksandr Lysenko1,, Olena Tachinina2 , Oleksandr Guida3, Iryna Alekseeva1 and Vladyslav
                                Kutiepov2
                                1
                                                                                                                                                   , 37, Prosp. Peremohy, Kyiv, 03056,
                                Ukraine
                                2
                                  National Aviation University, 1, Liubomyra Huzara ave., Kyiv, 03058, Ukraine
                                3



                                                 Abstract
                                                 The article considers the generalized integro-differentiating controller (GID-controller) as an alternative
                                                 to the PID-controller for use in cascaded SISO LTI systems for automatic control of mechatronic devices
                                                 of mobility nodes of humanoid robots. GID - controller is set by a generalized integro-differentiating
                                                 circuit or a connection of an ideal integrator with generalized prejudice-delay compensators. The article
                                                 shows that the main positive property of the GID controller compared to the PID controller is that, in the
                                                 presence of the SISO LTI mathematical model of the control object, the primary parametric setting of the
                                                 GID controller gives a practically acceptable rational result of controlling the robot movements. That is,
                                                 the initial parametric setting of the GID controller does not require further additional adjustment of the
                                                 controller parameters. This positive quality of the GID - controller allows you to significantly reduce the
                                                 time for adjusting the controller parameters on a real object. Therefore, the method of parametric
                                                 adjustment of the GID - controller was called the method of express adjustment of the generalized
                                                 integro-differentiating controller (MEA GID - controller). The result of a computer experiment is
                                                 presented, which showed that the MEA GID - controller provides quality, simplicity, convenience and
                                                 time saving during parametric adjustment of the controller, which justifies the expediency of using the
                                                 GID controller for controlling mechatronic devices of robot mobility nodes in general and, in particular,
                                                 humanoid robots with increased requirements to human-like movements.

                                                 Keywords
                                                 Automatic control system, PID-controller, integro-differentiating circuit, prejudice-delay
                                                 compensators 1


                                1. Introduction
                                When making humanoid robots designed to work next to human (housework, nursing robots,
                                service jobs in customer service areas), one of the main requirements is the implementation of the
                                principle of safe interaction between humanoid robots and human [1-6]. This principle is
                                implemented thanks to an approach that can be called "smoothness + sensuality" ("S + S"): the
                                robots perform movements that resemble (practically do not differ from) human movements, that
                                means, that they are smooth and sensual in terms of strength [5-7]. The structure of cascade
                                (multi-loop) SISO LTI automatic control systems with the properties of quasi-invariance (quasi-
                                adaptability) to the action of external disturbances in the best way ensures the implementation of
                                the "S + S" approach [8-11] (see Figure 1). Usually, in each of the cascades, controller with
                                parametric adjustment are used with a structure of PID-controller varieties (from a proportional
                                controller to a full structure with proportional, integral, and differential signals) [12, 13]. The
                                presence of an integrating link provides the property of quasi-invariance to external disturbances
                                such as step action [8]. Considering the fact that in humanoid robots the number of mobility nodes
                                exceeds hundreds [6, 7, 9], and in cascade automatic control systems of drives in mobility nodes, at
                                least two cascades are used, the total number of controllers that need to be adjusted can exceed
                                thousands. It is clear that an urgent engineering problem arises regarding the rapid adjustment (or

                                ICST-2024: Information Control Systems & Technologies, September , 23 25, 2024, Odesa, Ukraine
                                   lysenko.a.i.1952@gmail.com (O. Lysenko); tachinina5@gmail.com (O. Tachinina); guydasg@ukr.net (O. Guida);
                                alexir1@ukr.net (I. Alekseeva), vladcorvt@gmail.com (V. Kutiepov)
                                   0000-0002-7276-9279 (O. Lysenko); 0000-0001-7081-0576 (O. Tachinina); 0000-0002-2019-2615 (Guida O.); 0000-0002-
                                2878-6514 (I. Alekseeva), 0000-0002-1055-9698 (V. Kutiepov)
                                          Β© 2024 Copyright for this paper by its authors. Use permitted under Creative Commons License Attribution 4.0 International (CC BY 4.0).

CEUR
                  ceur-ws.org
Workshop      ISSN 1613-0073
Proceedings
re-adjustment) of controllers. Let us emphasize that, if computer numerical tuning programs are
used for parametric tuning according to algorithmically set criteria, then the actual task of quick
adjustment of the controller turns into the actual task of finding the first successful approximation,
that is, finding the initial conditions from which the computer adjustment algorithms "starts" [14
17].

2. Problem Statement
Two-cascade control systems (Figure 1) have better indicators of control quality compared to
single-cascade systems [8]. Therefore, they should be used where they do not exist yet. If two-
cascade control is abandoned in those mobility nodes, where it exists, in order to reduce the
number of controllers to be adjusted, the following positive properties of cascade control will be
lost:

   1.   External disturbances acting on the part of the control object that is covered by local
        feedback will directly affect the output coordinate, and won`t be reduced in the internal
        auxiliary loop.
   2.   Parametric disturbances that occur in the internal circuit will significantly affect the output
        signal.
   3.   The time of the transient response at the output of the system will increase significantly if
        the internal feedback is turned off, due to which the dynamic properties of the mechatronic
        devices of the mobility node are corrected.




Figure 1: Structural diagram of the cascade (multi-circuit) SISO LTI system of automatic control of
mechatronic devices in the mobility nodes of humanoid robots: π‘Šπ‘…π‘— (𝑠), π‘Šπ‘‚π‘— (𝑠) ( 𝑗 = 1,2) - scalar
continuous transfer functions that reflect the algorithm of the controller (index Rj) and
mathematical model of the control object (index Oj) corresponding to the first cascade (internal
circuit, j=1) and the second cascade (external circuit, j=2) of the automatic control system
    Currently, in mechatronic devices of mobility nodes of humanoid robots are used (almost 100%)
as controllers in both cascades PID-controllers [8-12]. As known, the PID-controller forms its
output signal as the sum of proportional, integral and differential signals from the error applied to
its input. We will use the so-called standard form to display the mathematical model of the PID-
controller operation algorithm
                                                 1       𝑇𝑑𝑗 βˆ™ 𝑠
                      π‘Šπ‘ƒπΌπ·π‘— (𝑠) = π‘˜π‘π‘— βˆ™ (1 +          +            ),                         (1)
                                               𝑇𝑖𝑗 βˆ™ 𝑠 πœπ‘‘π‘— βˆ™ 𝑠 + 1
where (according to recommendations [8]) we assume that the additional time constant can be
calculated from the following relation πœπ‘‘π‘— = 0.15 βˆ™ 𝑇𝑑𝑗 (𝑗 = 1,2) .
   Let us assume that the mathematical model of the control object is known π‘Šπ‘‚π‘— (𝑠) ( 𝑗 = 1,2). As
a controller in both cascades, it is planned to use a PID controller, that is, the structure of the
controller is known: π‘Šπ‘…π‘— (𝑠) = π‘Šπ‘ƒπΌπ·π‘— (𝑠) (𝑗 = 1,2). As a rule, a two-stage procedure is used to set
the parameters of both PID controllers, in which the smoothness (human-likeness) of movements is
implemented in the automatic control system: at the first stage, the initial adjustment is performed
using the Ziegler-Nichols or Cohen-Kun methods; at the second stage, the result of the initial
adjustment is improved using computer simulation. The duration and effectiveness of the second
stage significantly depends on the initial adjustment. The experience of adjusting the PID -
controller as a whole, shows that 90% of the time (and at the same time not always with the desired
result) is spent on the second stage of adjustment [8-12].
   A scientific-technical problem arises: to reduce the time spent on such a setting of the
cascade system of automatic control of mechatronic devices in mobility nodes, which ensures the
smoothness of the movements of the humanoid robot, that means that almost human-like
movements are achieved.
   The engineering experience of solving the problems of adjusting the structure and parameters
of the controllers indicates two effective approaches: first, to ensure a successful first
approximation to the acceptable structure and parameters of the controller (initial adjustment);
secondly, after the first approximation, adjust the minimum number of parameters.

3. Method of express adjustment of the generalized integro-
   differentiating controller (MEA GID- controller)
The initial data for solving the scientific-technical problem is: the structure of the cascaded SISO
LTI automatic control system (ACS) (Figure 1) and the mathematical model of the control object
π‘Šπ‘‚π‘— (𝑠) ( 𝑗 = 1,2).
    An ACS with PID controllers is considered as a prototype ACS, which should be improved by
usage of a GID - controller. PID - controllers are connected to the internal and external cascades
respectively (Figure 1).
    The task of synthesizing a cascade ACS (Figure 1) is considered solved if the algorithms of the
                                                                                π‘Šπ‘…π‘— (𝑠) (𝑗 = 1,2) .
    The general approach to the synthesis of transfer functions of regulators π‘Šπ‘…π‘— (𝑠) (𝑗 = 1,2) is as
follows. First, the regulator for the internal cascade (circuit) is synthesized with the mathematical
model of its operation algorithm, which is specified by the transfer function π‘Šπ‘…1 (𝑠) , where as a
mathematical model of the control object is considered π‘Šπ‘‚1 (𝑠) (Figure1). After that, the synthesis
of the regulator of the external cascade (circuit) is performed, which means that the transfer
function is found π‘Šπ‘…2 (𝑠). As a mathematical model of the control object, the serial connection of
the transfer functions of the internal cascade (circuit) and π‘Šπ‘‚2 (𝑠) .
    In order to solve the scientific-technical problem set above, it is proposed: replace the PID-
controllers in the cascade system-prototype with generalized integro-differentiating controllers
(GID - controllers), for the initial adjustment of which use a special method of initial express
adjustment with increased adjustment quality.
    Content of the main material: structure and parameters of GID - controller; method of express
adjustment of GID - controller; an example of the initial setting of PID and GID - controllers for a
cascade ACS and comparative modeling of transient response in a cascade ACS with synthesized
regulators.
    Structure and parameters of GID - controllers.
    As an alternative to the PID-controller algorithm, it is proposed to apply the algorithm, which is
given by the transfer function of the GID-controller:
                                                                                               (2)
                                   π‘˜πΊπΌπ· (𝑇1 𝑠 + 1)π‘š1 βˆ™ (𝑇3 𝑠 + 1)π‘š3
                       π‘Šπ‘ˆπΌπ· (𝑠) = 𝑣 βˆ™
                                    𝑠     (𝑇2 𝑠 + 1)π‘š2 βˆ™ (𝑇4 𝑠 + 1)π‘š4
parametric synthesis of which (search of parameters π‘˜π‘ˆπΌπ· > 0; 𝑣, π‘š1,2,3,4 ∈ {0; 1; 2; … }; 𝑇1,2,3,4 >
0) is proposed to be carried out by the method of express adjustment of the GID - controller (MEA
GID-controller). We emphasize once again that the MEA GID-controller is considered as a method
of primary parametric adjustment of the regulator.
    We will remind, that according to the terminology used in the national or English-language
scientific literature, the GID - controller is also called a generalized integro-differentiating circuit
or a connection of an ideal integrator with generalized prejudice-delay compensators.

    3.1. Stages of MEA GID-controller

   Stage 1. Create a mathematical model of an open circuit.
   Consider the open circuit in the form of a serial connection of mathematical models of the GID -
controller and control object and calculate the transfer function of the open circuit:
                                        (𝑇 𝑠+1)π‘š1 βˆ™(𝑇 𝑠+1)π‘š3               1
                                                                           βˆ— (𝑠)
 π‘Šπ‘‘ (𝑠) = π‘ŠπΊπΌπ· (𝑠) βˆ™ π‘Šπ‘œ (𝑠) = π‘˜πΊπΌπ· βˆ™ (𝑇1 𝑠+1)π‘š2 βˆ™(𝑇3 𝑠+1)π‘š4 βˆ™ 𝑠𝑣 π‘Šπ‘œ (𝑠) = π‘ŠπΊπΌπ·   βˆ™ π‘Šπ‘œβˆ— (𝑠), (3)
                                         2                 4

where
                                 βˆ— (𝑠)                         (𝑇 𝑠+1)π‘š1 βˆ™(𝑇 𝑠+1)π‘š3
                                π‘Šπ‘ˆπΌπ·   = π‘˜πΊπΌπ· βˆ™ (𝑇1 𝑠+1)π‘š2 βˆ™(𝑇3 𝑠+1)π‘š4 ;
                                                                2              4
                                                                 1
                                                 π‘Šπ‘œβˆ— (𝑠) = 𝑠𝑣 π‘Šπ‘œ (𝑠) ;
                                             1     𝑏 π‘ π‘š +𝑏              π‘ π‘šβˆ’1 +β‹―+𝑏 𝑠+𝑏
                               π‘Šπ‘œ (𝑠) = 𝑠𝑙 βˆ™ π‘Žπ‘š 𝑠𝑛 +π‘Žπ‘šβˆ’1π‘ π‘›βˆ’1 +β‹―+π‘Ž 1𝑠+π‘Ž 0,
                                                       𝑛        π‘›βˆ’1                1       0

π‘Šπ‘œ transfer function of the control object. We will remind that considered as known: 𝑙 the
number of ideal integrators in the mathematical model of control object; π‘š and 𝑛 - orders of
polynomials in the numerator and denominator π‘Šπ‘œ (𝑠) and the coefficients of these polynomials in
the corresponding powers 𝑠.
   Stage 2. Set the structural parameters of the transfer function π‘ŠπΊπΌπ· (𝑠) of GID-controller 𝑣 and
π‘š1,2,3,4 .
   Usually, 𝑣 and π‘š1,2,3,4 ∈ {0; 1; 2}.
   We choose 𝑣 (𝑣 the number of ideal integrators in the controller) to fulfill the requirement to
ensure the given order of astatism of the closed circuit (Figure 1) taking into account the number of
ideal integrators in the control object.
   We choose π‘š1,2,3,4 taking into account the properties of prejudice-delay compensators (PDC)
                                                  𝑇 𝑠+1                            𝑇 𝑠+1
                                   π‘Š12 (𝑠) = 𝑇1 𝑠+1 and π‘Š34 (𝑠) = 𝑇3 𝑠+1 .
                                                   2                                   4

    Analysis of the amplitude and phase-frequency characteristics of bias-delay compensators
allows us to draw the following conclusions:
        1. A PDC with the properties of an advance link allows you to increase the phase margin.
        2. A PDC with the properties of a delay link allows you to reduce the impact of high-
            frequency disturbances.
        3. The sequential inclusion of the prejudice-delay links with raising the binomials to the
            appropriate power π‘š1,2,3,4 allows you to obtain and strengthen both positive effects.
    Stage 3. Calculate the cutoff frequency πœ”π‘§ for the transfer function π‘Šπ‘œβˆ— (𝑠), that means, the
frequency at which |π‘Šπ‘œβˆ— (𝑗 βˆ™ πœ”π‘§ )| = 1.
    For designing the transfer function
                                                       1        𝑏 π‘ π‘š +𝑏          π‘ π‘šβˆ’1 +β‹―+𝑏 𝑠+𝑏
                                     π‘Šπ‘œβˆ— (𝑠) = 𝑠𝑣+𝑙 βˆ™ π‘Žπ‘š 𝑠𝑛 +π‘Žπ‘šβˆ’1π‘ π‘›βˆ’1 +β‹―+π‘Ž 1𝑠+π‘Ž 0                      (4)
                                                                    𝑛      π‘›βˆ’1                 1   0

and calculating the cutoff frequency, corresponding functions of the computer mathematics system
MATLAB+Simulink can be used.
   Stage 4. Calculation of the GID-controller parameters.
   Based on the known cut-off frequency πŽπ’› , calculate the parameters of the transfer function of
the GID-controller:
                                   3.3           1            1            1
                            𝑇1 =       ; 𝑇3 =        ; 𝑇2 =       ; 𝑇4 =        .
                                   πœ”π‘§         3.3βˆ™πœ”π‘§        33βˆ™πœ”π‘§        330βˆ™πœ”π‘§

   Written down ratios make it possible to obtain a successful first approximation to the
acceptable values of the GID-controller parameters at any values of its structural parameters 𝒗 and
π’ŽπŸ,𝟐,πŸ‘,πŸ’ . These ratios are obtained as a result of empirical generalization of the experience of
synthesis of GID - controllers.
   Stage 5. Selection of structural parametersπ‘š1,2,3,4 and the gain of the regulator π‘˜πΊπΌπ· .
   First step: set π‘š1,2,3,4 = 1.
   Design a computer mathematical model and perform a simulation experiment for selection π‘˜πΊπΌπ·
. Recommended: Start with a value π‘˜πΊπΌπ· ∈ [0.1 βˆ™ π‘˜π‘πΎπ‘βˆ’π‘ ; 0.5 βˆ™ π‘˜π‘πΎπ‘βˆ’π‘ ], where π‘˜π‘πΎπ‘βˆ’π‘            gain
coefficient, which is calculated by the method of Ziegler-Nichols oscillations [8], and choose such a
value of π‘˜πΊπΌπ· , at which the duration of the transient response, oscillation and overregulation will
have acceptable values for the specific task.
   If it was successful in choose π‘˜πΊπΌπ· , in which the above stated parameters of the transient
response satisfy the requirements of a specific task, then we consider that stage 5 is completed.
   If it was not possible to meet the requirements of a specific task, then we change one of the
structural parameters π‘š1,2,3,4 and repeat the selection π‘˜πΊπΌπ· .
   During the initial setup, stage 5 is performed in the interactive "manual" mode. Experience
shows that no more than a few "runs" of a computer mathematical model are enough to obtain a
result acceptable for practical use, that is, before stopping the adjustment as a whole. Let us
emphasize that the "human-likeness" of movements is assessed by experts, although formally the
smoothness of movements can be specified using well-known standard forms [12].

4. Results and discussions
   In this section, we will consider an example of the initial setting of the PID- and GID-
controllers for a cascade ACS (Figure 1) and perform a comparative simulation of transient
responses in a cascade ACS with synthesized controllers. Let us start with the initial setting of the
PID - controllers for the cascade ACS (Figure 1).
   The initial setting of the PID - controllers for the cascade ACS (Figure 1). For the initial setting
of the PID-controller parameters (that is, determination of the gain π‘˜π‘ , constant of integration 𝑇𝑖
and constant of differentiation 𝑇𝑑 ) we will use well-known methods [8]: Ziegler-Nichols
oscillations (KZ-N); Ziegler - Nichols when using the transient response (Z-NPP); Cohen Kuhn
using the transient response (K-KPP). To demonstrate the exact method by which the parameters
of the PID controller were determined, we will use the corresponding indices for each of the
methods: π‘˜π‘πΎπ‘βˆ’π‘ , π‘‡π‘–πΎπ‘βˆ’π‘ , π‘‡π‘‘πΎπ‘βˆ’π‘ ; π‘˜π‘π‘βˆ’π‘π‘ƒπ‘ƒ , π‘‡π‘–π‘βˆ’π‘π‘ƒπ‘ƒ , π‘‡π‘‘π‘βˆ’π‘π‘ƒπ‘ƒ ; π‘˜π‘πΎβˆ’πΎπ‘ƒπ‘ƒ , π‘‡π‘–πΎβˆ’πΎπ‘ƒπ‘ƒ , π‘‡π‘‘πΎβˆ’πΎπ‘ƒπ‘ƒ .
   The above-mentioned methods require conducting a computer or field experiments [8-12].
According to the results of the experiment, auxiliary parameters are determined, which we will
denote as 𝐾0 , 𝑃0 ; π‘₯0 (𝑑0 ), π‘₯0∞ = lim π‘₯0 (𝑑) , 𝑑0 , 𝑑1 , 𝑑2 . Physical and mathematical content of these
                                    π‘‘β†’βˆž
parameters is illustrated with the help of figures 2 and 3.




Figure 2: A transient response with constant oscillations with a period of 𝑃0 at the output of a
closed system consisting of a proportional controller with a gain factor of 𝐾0 and a control object:
the values 𝐾0 , 𝑃0 are used to calculate the parameters (initial setting) of the PID - controller using
the Ziegler-Nichols oscillation method

   According to the known 𝐾0 , 𝑃0 π‘Žπ‘›π‘‘ π‘₯0 (𝑑0 ), π‘₯0∞ = lim π‘₯0 (𝑑) , 𝑑0 , 𝑑1 , 𝑑2 parameters of the PID-
                                                       π‘‘β†’βˆž
controller are calculated using the above-mentioned methods [8]:
                                                                               𝑃
                          π‘˜π‘πΎπ‘βˆ’π‘ = 0.6 βˆ™ 𝐾0 , π‘‡π‘–πΎπ‘βˆ’π‘ = 0.5 βˆ™ 𝑃0 , π‘‡π‘‘πΎπ‘βˆ’π‘ = 80;
                                         𝐷
                                        2
                     π‘˜π‘π‘βˆ’π‘π‘ƒπ‘ƒ = 1.2 βˆ™ 𝐷 βˆ™π· , π‘‡π‘–π‘βˆ’π‘π‘ƒπ‘ƒ = 2 βˆ™ 𝐷1 , π‘‡π‘‘π‘βˆ’π‘π‘ƒπ‘ƒ = 0.5 βˆ™ 𝐷1;
                                         0   1

                                                   𝐷2               𝐷1
                                 π‘˜π‘πΎβˆ’πΎπ‘ƒπ‘ƒ =               βˆ™ (0.9 +         ),
                                                 𝐷0 βˆ™ 𝐷1          12 βˆ™ 𝐷2
                                                    30 βˆ™ 𝐷2 + 3 βˆ™ 𝐷1
                                   π‘‡π‘–πΎβˆ’πΎπ‘ƒπ‘ƒ = 𝐷1 βˆ™                    ,
                                                    9 βˆ™ 𝐷2 + 20 βˆ™ 𝐷1
                                                        4βˆ™π· βˆ™π·
                                                        1 2
                                         π‘‡π‘‘πΎβˆ’πΎπ‘ƒπ‘ƒ = 11βˆ™π· +0.2βˆ™π· ,
                                                          2       1
where 𝐷0 = π‘₯0∞ βˆ’ π‘₯0 (𝑑0 ), 𝐷1 = 𝑑1 βˆ’ 𝑑0 , 𝐷2 = 𝑑2 βˆ’ 𝑑1 .




Figure 3: Transient response at the output of the control object: value π‘₯0 (𝑑0 ), π‘₯0∞ =
lim π‘₯0 (𝑑) , 𝑑0 , 𝑑1 , 𝑑2 are used to calculate parameters (initial setting) of the PID - controller by
π‘‘β†’βˆž
methods of Ziegler-Nichols and Cohen Kuhn when using a transient response (a tangent is
drawn at point A (the inflection point of the transient response))
   We consider that the mathematical model of the mechatronic device in each of the mobility
nodes is given by the continuous transfer functions π‘Šπ‘œ1 (𝑠) and π‘Šπ‘œ2 (𝑠) (Figure 1). The transfer
function π‘Šπ‘œ1 (𝑠) can be calculated as the transfer function of serially connected links of directed
action. These links are a power transformator with a transfer function π‘Šπ‘“ (𝑠) and a direct current
motor (DC) with a transfer function π‘Šπ·π‘ƒπ‘† (𝑠), where the output signal is considered to be the
angular speed of the rotor, i.e.:
                                       π‘Šπ‘œ1 (𝑠) = π‘Šπ‘“ (𝑠) βˆ™ π‘Šπ·π‘ƒπ‘† (𝑠).
    The first cascade (circuit) consists of a controller with a transfer function π‘Šπ‘…1 (𝑠) and a control
object with a transfer function π‘Šπ‘œ1 (𝑠) (Figure 1). The second cascade consists of a controller with
a transfer function π‘Šπ‘…2 (𝑠) and a control object with a transfer function equal to the product of the
                                                                π‘Š (𝑠)βˆ™π‘Šπ‘œ1 (𝑠)
transfer functions of the first cascade (circuit), 𝐻1 (𝑠) = 𝑅1 (𝑠)βˆ™π‘Š
                                                                 1+π‘Šπ‘…1
                                                                               to the transfer function
                                                                         π‘œ1 (𝑠)
          1
π‘Šπ‘œ2 (𝑠) = 𝑠 , where the output signal is the rotation angle of the DC motor rotor.
  When performing computer experiments, typical transfer functions of mechatronic devices of
mobility nodes of humanoid robots were used [8-12]:
                                           10                             1
                         π‘Šπ‘“ (𝑠) = 0.0003βˆ™π‘ 2 +0.04βˆ™π‘ +1 ;π‘Šπ·π‘ƒπ‘† (𝑠) = 0.4βˆ™π‘ 2 +1.3βˆ™π‘ +1 .
    Let us perform the initial setting of the PID - controller of the first cascade using the above
stated methods. During the initial setting of the PID - controller of the first cascade (internal
circuit), the following auxiliary values of parameters were obtained:
             𝐾0 = 3.5, 𝑃0 = 0.72 𝑠; π‘₯0 (𝑑0 ) = 0, π‘₯0∞ = 10, 𝑑0 = 0 𝑠, 𝑑1 = 0.25 𝑐, 𝑑2 = 2 𝑠 .
   The result of calculating the parameters of PID - controller of the first cascade (internal circuit):
                    π‘˜π‘πΎπ‘βˆ’π‘1 = 2.100 , π‘‡π‘–πΎπ‘βˆ’π‘1 = 0.3600 𝑠 , π‘‡π‘‘πΎπ‘βˆ’π‘1 = 0.0900 𝑠;
                   π‘˜π‘π‘βˆ’π‘π‘ƒπ‘ƒ1 = 0.8400 , π‘‡π‘–π‘βˆ’π‘π‘ƒπ‘ƒ1 = 0.500 𝑠 , π‘‡π‘‘π‘βˆ’π‘π‘ƒπ‘ƒ1 = 0.1250 𝑠;
                  π‘˜π‘πΎβˆ’πΎπ‘ƒπ‘ƒ1 = 0.6383 , π‘‡π‘–πΎβˆ’πΎπ‘ƒπ‘ƒ1 = 0.6416 𝑠, π‘‡π‘‘πΎβˆ’πΎπ‘ƒπ‘ƒ1 = 0.0907 𝑠.
   Let us perform the initial setting of the GID - controller of the first cascade using the MEN UID
- regulator. Assume that 𝑣 and π‘š1,2,3,4 are equal to one. We will calculate the cut-off frequency
and parameters of the GID - regulator of the first cascade (internal circuit). We will use the
positioning of the MEA GID-controller and as a result we will get:
                                π‘Ÿπ‘Žπ‘‘
                   πœ”π‘§1 = 2.61 𝑠 ,𝑇11 = 1.2644 𝑠, 𝑇31 = 0.1161 𝑠, 𝑇21 = 0.0116 𝑠,
                    𝑇41 = 0.0012 𝑠 . We assume, that π‘˜π‘ˆπΌπ·1 = 0.1 βˆ™ π‘˜π‘πΎπ‘βˆ’π‘1 =0.21 .
    Before moving on to adjusting the parameters of the controller of the second cascade (Figure 1),
consider the transient response at the output of the two-cascade ACS under the condition that
 π‘Šπ‘…2 (𝑠) = 1. Let us clarify the issue of the necessity to complicate the general procedure for
setting up a two-stage ACS by adjusting the controller for the second cascade.
    Let us perform a visual analysis of the transient responses (see Figure 4) at the output of the
two-stage ACS (Figure 1).
    This analysis shows: when using the GID-controller, the initial setting of which is performed
using the proposed method, the smoothness of movements in the mobility node of the humanoid
robot will be better than when using the PID-controller.
    Thus, if there is no need to reduce the readjustment and the duration of the transient response,
then it can be assumed that the cascade ACS can be designed with the connection of only one GID-
controller in the internal cascade of the ACS (we consider the external cascade as having a
proportional controller with a gain factor connected to it, which is equal to one).
    We especially emphasize that it is not necessary to adjust the GID - controller: only the initial
setting is enough.
    Let us assume that there is still a necessity to reduce the duration of the transient response
while maintaining the smoothness of the movements in the mobility node. Then we will perform
                                                                                                     1
the initial setting of the PID and GID-controllers for the second cascade (circuit). Since π‘Šπ‘‚2 (𝑠) = 𝑠
, then only the Ziegler-Nichols oscillation method can be used to adjust PID-controller.




   Figure 4. Transient response at the output of the two-cascade ACS (Figure 1) under the
condition that π‘Šπ‘…2 (𝑠) = 1 and different π‘Šπ‘…1 (𝑠): 1, 2, 3, transfer function π‘Šπ‘…1 (𝑠) is equal to the
transfer function of the PID - controller, the parameters of which are adjusted by the methods of
Ziegler-Nichols oscillations, Ziegler-Nichols when using a transient response, and Cohen-Kuhn
when using a transient response respectively; 4 - transfer function π‘Šπ‘…1 (𝑠) is equal to the transfer
function of the GID - controller.
    During the initial setting of the PID - controller of the second cascade (external circuit), the
following auxiliary values of parameter were obtained
                                        𝐾0 = 3.5, 𝑃0 = 0.72 𝑠 .
   Based on these values, the parameters of the PID-controller of the second cascade were
calculated:
                    π‘˜π‘πΎπ‘βˆ’π‘2 = 1.800 , π‘‡π‘–πΎπ‘βˆ’π‘2 = 0.3300 𝑠 , π‘‡π‘‘πΎπ‘βˆ’π‘2 = 0.0825 𝑠.
    Let us move on to setting up the GID-controller. Assume, that 𝑣 and π‘š1,2,3,4 equal to one.
During the initial setting of the GID-controller of the second cascade (external circuit), the
                                                               π‘Ÿπ‘Žπ‘‘
following value of the cutoff frequency was obtained πœ”π‘§2 = 0.99 𝑠 and the values of the GID-
controller were calculated:
                𝑇12 = 3.3333 𝑠, 𝑇32 = 0.3061 𝑠, 𝑇22 = 0.0306 𝑠, 𝑇42 = 0.0031 𝑠 .
    Let us assume, that π‘˜π‘ˆπΌπ·2 = 0.2 βˆ™ π‘˜π‘πΎπ‘βˆ’π‘2 =0.36 . As we can see (see Figure 5), the use in the
second cascade of the GID-controller with the initial setting of parameters allows to obtain
smoother movements in the mobility node of the humanoid robot compared to the use of the
initially configured PID-controller. We note, that there will be almost no readjustment. In order to
improve the result of the GID-controllers application, we will change the structural parameters of
the GID-regulators in both cascades (circuits). Let us assume, that 𝑣 = 1 and π‘š1,2 = 1, π‘š3,4 = 2
provided that the parameter values 𝑇1𝑗 , 𝑇3𝑗 , 𝑇2𝑗 , 𝑇4𝑗 (𝑗 = 1,2) and π‘˜π‘ˆπΌπ·1 remained unchanged.
Let us assume, that π‘˜π‘ˆπΌπ·2 = 0.5 βˆ™ π‘˜π‘πΎπ‘βˆ’π‘2 . It was possible to reduce the duration of the transient
response by almost two times while maintaining the smoothness of the movement and the absence
of readjustment (Figure 6).




Figure 5. Transient responses at the output of the two-cascade ACS (Figure 1) under the condition
that: 1 transfer functions π‘Šπ‘…1 (𝑠)and π‘Šπ‘…2 (𝑠)are equal to the transfer functions of the PID -
controllers, parameters of which are adjusted by the Ziegler - Nichols oscillation method; 2
transfer functions π‘Šπ‘…1 (𝑠) and π‘Šπ‘…2 (𝑠) are equal to the transfer functions of the GID - controllers,
under the condition that the structural parameters of the GID - controllers 𝑣 and π‘š1,2,3,4 are equal
to one
Figure 6. Transient responses at the output of the two-cascade ACS (Figure 1) under the condition
that: 1 transfer functions π‘Šπ‘…1 (𝑠) and π‘Šπ‘…2 (𝑠) are equal to the transfer functions of the PID -
controllers, parameters of which are adjusted by the Ziegler - Nichols oscillation method; 2 the
transfer functions π‘Šπ‘…1 (𝑠) and π‘Šπ‘…2 (𝑠) are equal to the transfer functions of the GID - controllers,
under the condition that the structural parameters of the GID - controllers 𝑣 v and π‘š1,2,3,4 are
equal to one, and π‘š3,4 = 2

5. Conclusions
As an alternative to the PID-controller for usage in cascaded SISO LTI systems for automatic
control of the mobility nodes of humanoid robots, it is proposed to use a controller, which is set by
a generalized integro-differentiating circuit or a connection of an ideal integrator with generalized
prejudice-delay compensators. This controller was called a generalized integro-differentiating
controller (GID - controller).
   The main positive property of the GID-controller compared to the PID -controller is that when
SISO LTI mathematical model of the control object is available, the primary parametric setting of
the GID-controller gives an almost acceptable rational result of control that does not require
additional adjustment. Therefore, the method of parametric adjustment of the GID - controller was
named the method of express adjustment of the generalized integro-differentiating regulator (MEA
GID-controller).
   MEA GID-controller provides quality, simplicity, convenience and time saving during
parametric adjustment of the controller, which justifies the expediency of using the GID -
controller for controlling mechatronic devices of robot mobility nodes in general and, in particular,
humanoid robots with increased requirements for human-like movements.
   The computer experiment presented in the article illustrated and confirmed the positive
properties of the GID - controller, which give it an advantage compared to the PID - controller.
   Further research will be aimed at developing a method of structural-parametric synthesis of the
GID - controller and its application in automatic control systems with significant nonlinearities.

6. References
[1] S. V. Viraktamath, Rashmi N. Raikar, Humanoid Robot: A Review, International Journal for
     Research     in    Applied      Science and  Engineering   Technology      (2021) 2884
     2894. doi:10.22214/ijraset.2021.37890.
[2] A. Pon Bharathi et al., An Investigation on Humanoid Robots with Biped Locomotion and
     Walking, Design, Modelling and Fabrication of Advanced Robots (2022) 55 61.
     doi:10.46632/dmfar/1/1/9.
[3] S H. Ugiuchi, J.Yoneshima, S. Kohama, Development of control system for small humanoid
     robot, in: Proceedings of JSME annual Conference on Robotics and Mechatronics(Robomec),
     2017, pp. 2A1 C09. doi:10.1299/jsmermd.2017.2a1-c09
[4] L. M. CabΓ‘s Ormaechea, Humanoid Robots, Advanced Mechanics in Robotic Systems. London
     (2011) 1 18. doi:10.1007/978-0-85729-588-0_1.
[5] S. Aswini et al., Biomechanics-Inspired Control Strategies for Humanoid Robots, in: IEEE Uttar
     Pradesh Section International Conference on Electrical, Electronics and Computer Engineering
     (UPCON),         Gautam      Buddha       Nagar,      India,     1 3       December,        2023.
     doi:10.1109/upcon59197.2023.10434421
[6] E. Mogos, The fuzzy-based systems in the communication between a human and a humanoid
     robot, Journal of Physics: Conference Series 2251 1 (2022) 012003. doi:10.1088/1742-
     6596/2251/1/012003
[7] S. N. Bhatlawande, Prachi Prashant Kulkarni. Humanoid Robot. International Journal of
     Advanced Research in Science, Communication and Technology (2023) 32.
[8] G.C. Goodwin, S. F. Graebe, M. E., Salgado Control system design, 1 st Edishen by Prentice Hall,
     2001.
[9] Yu. Irlyk, A. Stopakevych, Analysis of the perspectives of application of artificial intelligence
     technologies for building autonomous industrial systems of automatic control, Automation of
     technological and business processes, 15 4 (2024) 8-13.
[10] doi:10.15673/atbp.v15i4.2578
[11] Y. Wang et al., Development of a humanoid robot control system based on AR-BCI and SLAM
     navigation, Cognitive Neurodynamics (2024). doi:10.1007/s11571-024-10122-z .
[12] O. Cetin, Present Applications of Humanoid Robots and Fuzzy Control, Toward Humanoid
     Robots: The Role of Fuzzy Sets. Cham (2021) 237 251. doi:10.1007/978-3-030-67163-1_9.
[13] O.I. Lysenko, Technique for tuning digital controllers of robot drives, taking into account the
     correction of their dynamic characteristics, Problems of informatization and control, 4 68
     (2021) 44 51. doi:10.18372/2073-4751.68.16526.
[14] M. O. Morneva, S. M. Golubeva, A. S. Toropov, Digital automatic control systems in
     mechatronic systems, Visnik of the Volodymyr Dahl East Ukrainian National University 5 269
     (2021) 20 23. doi:10.33216/1998-7927-2021-269-5-20-23.
[15] L. Zhao et al., Compliant Motion Control of Wheel-Legged Humanoid Robot on Rough
     Terrains,     in:    IEEE     Transactions     on    Mechatronics       (ASME)      (2023) 1 11.
     doi:10.1109/tmech.2023.3320762
[16] O.V. Solomentsev, M.Yu. Zaliskyi, O.V. Zuiev, M.M. Asanov, Data processing in exploitation
     system of unmanned aerial vehicles radioelectronic equipment, in: IEEE International
     Conference Actual Problems of Unmanned Air Vehicles Developments (APUAVD), Kyiv,
     Ukraine, 2013, pp. 77 80. doi: 10.1109/APUAVD.2013.6705288.
[17] O.V. Solomentsev, M. Yu. Zaliskyi, M. M. Asanov, V. H. Melkumyan, UAV operation system
     designing, in: IEEE 3rd International Conference on Actual Problems of Unmanned Aerial
     Vehicles Developments (APUAVD), Kyiv, Ukraine, 2015, pp. 95 98. doi:
     10.1109/APUAVD.2015.7346570.
[18]                                         Methodology of Synthesizing Digital Regulators in
     Precision Electric Drives for Orientation and Stabilization Target Tracking System of Mobile
                                       in: Proceedings of the 11-th International Conference
     "Information Control Systems & Technologies", Odesa, Ukraine, 2023, pp. 51-63. https://ceur-
     ws.org/Vol-3513/.