Open access peer-reviewed chapter

# Adjustment of the PID Gains Vector Due to Parametric Variations in the Plant Model in Terms of Internal Product

Written By

José Pinheiro de Moura and João Viana da Fonseca Neto

Submitted: July 15th, 2020 Reviewed: November 17th, 2020 Published: January 5th, 2021

DOI: 10.5772/intechopen.95051

From the Edited Volume

## Control Based on PID Framework

Edited by Wei Wang

Chapter metrics overview

View Full Metrics

## Abstract

The tuning of the gains of a controller with proportional-integral-derivative (PID) actions has been prevalent in the industry. The adjustment of these gains in PID controllers is often determined by classical methods, such as Ziegler-Nichols, and trial and error. However, these methods fail to deliver satisfactory performance and often do not meet specific project demands because of the inherent complexity of industrial processes, such as plant parameter variations. To solve the tuning problem in highly complex industrial processes, a controller adjustment method based on the internal product of PID terms is proposed, and a propagation matrix (PM) is generated by the numerator coefficients of the plant transfer function (TF). In the proposed method, each term of the PID controller is influenced by each of the numerator and the denominator coefficients. Mathematical models of practical plants, such as unloading and resumption of bulk solids by car dumpers and bucket wheel resumption, were employed to evaluate the proposed method. The obtained results demonstrated an assertive improvement in the adjustment gains from PID actions, thereby validating it as a promising alternative to conventional methods.

### Keywords

• parameter variations
• industrial processes
• internal product
• propagation matrix
• PID actions

## 1. Introduction

The tuning of the parameters of PID controllers is challenging and requires expertise to achieve superior performance [1]. PID controllers are extensively used in the industries. However, the controllers are often implemented without a derivative action because of the highly sensitive tuning of parameters, which affects the efficiency of the controller [2]. This study presents a methodology for tuning three terms of the PID controller simultaneously to ensure overall efficiency of the controller [3].

The advantage of the PID controller tuning methodology, which is based on the internal product of the PID terms that generates the propagation matrix (PM), is that a vector of the specified parameters of a characteristic polynomial can be projected, and an error vector is obtained on comparison with the parameters of the characteristic plant polynomial [4, 5]. The method minimizes the error from the specified parameters, thereby facilitating the design of a high-performance PID controller. The method also enables the allocation of the poles by direct replacement using specified parameters, thereby ensuring the desired operating point of the control system.

This chapter presents a formulation proposal to resolve the PID controller tuning problem. The proposal is based on the dot product of the gain vector parameters of the controller and the rows of the propagation matrix. The dot product represents the changes in the behavior of the plant that are determined by the parametric variations in the coefficients of the TF polynomial characteristic.

The following topics and proposal development are presented in the remainder of this chapter. In Section 2, a preliminary on the transfer functions of the plant and PID controller are presented. In terms of internal product, the main properties of PID controllers and the development of proposed method are presented in Section 3. Taking into account three industrial plants of the mining sector, computational evaluation experiments of the PID tuning proposal are presented in Section 4. Finally, the conclusion of the work is presented in Section 5.

## 2. Preliminaries

Adjusting the PID controller gain parameters is not a trivial task and requires in-depth knowledge from the experts. In this work, the problem of tuning PID controllers is based on original studies regarding polynomial compensators and problems strongly related to the specification of parameters to meet operational constraints in plant dynamics, presented as a particular form of compensators in the s domain [6, 7].

### 2.1 Mathematical model of the plant in terms of transfer function

The plant’s dynamic system is represented by ordinary differential equations (ODE), described by TFs in the s domain (Laplace transform). The ODE concept in terms of TFs established in this chapter is in accordance with the block diagram shown in Figure 1, where the closed loop system relates the input and output signals: Rs is reference input, Ws is disturbance signal and Vs is noise signal, Ys is plant output and Yms is plant output measured by the sensor, U(s) and E(s) are the control effort and closed loop error, respectively, that are internal control system performance variables.

Applying the Laplace transforms to the control elements of Figure 1, the generalized TF with a polynomial structure in the s domain is obtained that is given by

Ggs=bosm+b1sm1++bm1s+bmaosn+a1sn1++an1+an,E1

where Ggs is the general TF of the control system blocks diagram, related to Figure 1, n is the order of the plant model and the number of poles that are entered into the system and m is the number of zeros, which is associated with the PID actions of the controller. As an imposition of the controller gains values, the coefficients ai and bk are the adjustable parameters to compensate for the parametric variations of the plant. In the proposed formulation, the bk coefficients are kept constant and adjustment are made only to the ai parameters.

### 2.2 PID controller model

The controller model associated with the TF given in Eq. (1) is customized to perform the actions of the controller’s PID terms, where n=1, m=2, a0=1, the PID actions are represented by the transfer function that is given by

CKpids=KDs2+KPs+KIs,E2

where CKpids is the controller model associated with the TF given in Eq. (1).

Adjustments of parameters that meet the project specifications, can be found in a large number of scientific and technical publications in controle specialized books, conferences and high quality journals [8]. The importance of developing methods for adjusting parameters of PID controllers and systematizing applications in industrial processes of real-world plants, has the objective of meeting the project specifications contained in technological advances, in order to guarantee the optimal adjustment of the parameters of the PID term of the controller [9, 10]. The challenge of tuning with optimal performance of the parameters of a PID controller, started around 1920 and continues to the present days [11, 12, 13].

The parameters of the PID controllers are adjusted to adapt to the tuning needs in a combination of proportionality associated with the proportional action, lead associated with the derivative action, and delay associated with the integral action of the error signal. However, there are still many problems that can be solved with computational intelligence-based algorithms. The purpose of this work is to contribute with a method of tuning PID controllers, which can support the development of electronic devices that contribute to technological advancement and the evolution of industry 4.0 with logical planning units, for optimal, robust decision-making and adaptability [14, 15]. Such units must be based on digital control technologies and embedded systems [16] in real time [17], to be reliably deployed in real-world systems [18].

To meet the demands of design specitifications, the proposed solution contributes to the evolution in approaches of optimal and adaptive control, providing the optimization of the figures of merit [19], ensuring a solution with satisfactory performance, meeting the requirements specified in projects, in a way that minimizes efforts of computational cost and control.

TF is specified in the factored form, that is, by the roots of the numerator and denominator polynomials associated with Eq. (1). TF in the factored form is represented in terms of product, where the designer inserts the specified or desired parameters. TF in the form of a product is given by

Ggs=Kk=1msszki=1nsspi,E3

where spi=ωdi are the poles and szk=σkωdk are the zeros of the dynamic system. The poles and zeros of the system are represented by the pair (ζ, ωn), the first component is the damping factor and the second is the undamped natural frequency, σi=ζiωni, and ωdi=ωni1ζi2.

The PID controller adjustment problem is formulated based on the parametric difference between the specified coefficients and the original coefficients of the TF denominator polynomial. This formulation is based on the references [20, 21]. Where the authors present the development of models for optimized online optimization that is based on computational intelligence approaches.

The formulation of the proposed PID controller adjustment problem is presented in this section and is illustrated by the block diagram of Figure 1. Thus, in the context of the proposal, the performance matrix of the control system provides the means to determine the values of the parameters of the gain vector of the PID controller Kpid.

The models are represented in the form of internal produc <,>, which is a notation widely used in this text as the product of two vectors (internal product) and the models are called internal product models of the plant. The internal product is the appropriate form for analysis, allowing the designer to observe the impact of the earnings vector parameters Kpid (KD, KP and KI), in the output of the dynamic system associated with the polynomial coefficients of the TF denominator.

The PID controller model, in terms of the ODE equations and the Laplace transform, obtains the TF of the PID controller in terms of the dot product, which is given by

CKpids=Kpidspids,E4
Kpid=KDKPKIT,E5

and spid is vector-powered in s of the PID gains associated with transfer function numerator that is given by

spid=s2s1s0T.E6

### 3.1 PID model in the form of internal product

Inserting the characteristics of the plant, through the values of the coefficients of the polynomials of the numerator and the denominator (poles and zeros), associated with the mathematical models in terms of TFs given in Eq. (1) in internal product form and in Eq. (3) in generic form, which is given by

Gps=Kbksmsn+aisn1,E7

where Gps is the TF in the form of internal product and the coefficients ai with i=1,,n and bk with k=0,,m are a combination of the spi poles and szk zeros.

### 3.2 Open-loop transfer function

The open-loop FT or direct branch of the control system is given by

GpOLs=KKpidspidbksmsn+1+saisn1,E8

where GpOLs is the TF of the plant in the open loop and K the gain of the plant and n>m.

The structure of TF is determined by the relationship n+12+m. For n+1=2+m, the system is proper and for n+1<2+m the system is strictly proper, thereby establishing a general relationship between the order of the PID controller and the order of plant dynamics. This relationship ensures that the system structure is adequate, not allowing the system to present a nonpractical structure. In this way, it establishes that the relationship of the closed loop system is given by

ncl=mPID+n,E9

where mPID can only assume 0 (zero) or 1 (one) values. The PID is observed to impose a proper TF, if the closed-loop system is of order n+1 and the numerator is of order mPID+m.

According to the block diagram of Figure 1, the TFs Ys/Rs, Ws=0, and Vs=0 are given by

GpCLs=CKpidsGs1+CKpidsGsHs.E10

### 3.3 Propagation of PID terms x bk coefficients

The development of the polynomials of the numerator (zeros) and the denominator (poles) consists of the propagation of the gain vector Kpid of the controller by the numerator coefficients (bk) associated with the coefficients of the denominated (ai) TF of the plant. The equationing of the problem is given in the form of an internal product that weights the coefficients of the polynomial of zeros in the closed-loop and additive to the dynamics of the closed-loop transfer function.

#### 3.3.1 Polynomial of zeros

When replacing Eqs. (1) and (2) in Eq. (10), the numerator polynomial of the closed-loop TF is obtained, which is given by

NCLs=CKpidsGs.E11

Expanding and ordering Eq. (11), one obtains

NCLs= b0KD+b1KP+b2KIsm+mpid+b1KD+b0KP+b1KIsm+mpid1+b2KD+b1KP+b0KIsm+b3KD+b2KP+b1KIsm1.E12
.+bmKD+bm1KP+bm2KIs2+bm+1KD+bmKP+bm1KIs1+bm+2KD+bm+1KP+bmKIs0.

In terms of inner product, the general polynomial form of the closed-loop numerator polynomial is given by

NCLs=i=0mclKpidb¯k1smcli,E13

where mcl=m+mpid and vector b¯k of the polynomial of zeros of the closed-loop system is given by

b¯k=bkbk1bk2.E14

In similar way Eq. (13), one obtains the closed-loop denominator polynomial is given by

DCLs=sn+nDpid+i=0n+nDpid1bkKD+bk1KP+bk2KI+ai+1,E15

where nDpid=0 or 1. When nDpid=0, the PID controller structure have the terms derivative and proportional. When nDpid=1, the structure of the PID controller has an integrator term that increases the order of the system by 1, starting with the three terms: proportional, derivative and integrative [21]. In this work, when nDpid=0, the proposal is to specify an additional ais coefficient., to ensure that the PID controller has the three terms.

#### 3.3.2 Characteristic polynomial

The general form of the closed-loop denominator polynomial is given by

PCLs=sncl+i=1ncl1ai+Kpidb¯k1sncli,E16

where PCLs is the general form of the characteristic polynomial of the closed-loop plant and ncl is the order of referred polynomial.

the characteristic closed-loop polynomial for unit feedback (Hs=1) is obtained in a similar way and given by

PpCLs=1+CKpidsGs.E17

PpCLs is the characteristic polynomial of the closed loop. The representation of the problem in the form of an internal product that relates the coefficients of the zero polynomial with the coefficients of the closed plant dynamics is the basis for obtaining the characteristic closed-loop polynomial.

### 3.4 Proposed method

The problem is formulated based on the propagation matrix generated from the dot product between the terms of the earnings vector Kpid with the coefficients of the numerator bk associated with the coefficients of the denominator ai of the TF of plant. The propagation matrix product associated with the TF numerator coefficients of the plant, give rise to a new characteristic polynomial based on new specified operating points, which are imposed by new zeros and new poles.

#### 3.4.1 Propagation matrix of PID design

The design is based on the propagation matrix, allowing the designer to specify new points of operation that improve the performance of the controller acting on the plant dynamics, where changes in the order and coefficients of the characteristic polynomial can be observed through the internal product of the zero and gain coefficients of PID controllers.

The problem is formulated based on the propagation matrix B¯, which is a consequence of the interaction between the parameters of the gain vector Kpid (KD, KP, KI) with the coefficients of the plant TF numerator, this matrix is represented by

B¯=b000b1b00b2b1b0b3b2b1bmbm1bm20bmbm100bm.E18

One case notice in [20] that the diagonals are not repeated, and they vary according to the order n of the system’s characteristic polynomial. The propagation of the gains is weighted by the coefficients of the numerator polynomial. The closed-loop TF of the plant is given in terms of the product of the gains KD, KP, and KI with the coefficients bk, bk1, and bk2 of the closed-loop TF.

The law of formation of the propagation Matrix (18) is ruled by m+2 rows and 3 columns. The rows represent the order of the system, starting with the propagation in the dynamics of order sncl and ending in the dynamics of order s0 zero. The columns represent the gains of the controller in the poles and zeros of the plant dynamics.

#### 3.4.2 Proposed characteristic polynomial

The proposed characteristic polynomial based on propagation matrix of PID controller gains idea is presented. From the system of equations that represents the actions of the PID controller in the plant dynamics, the formulation of the adjustment problem is established from the perspective of the inner product of the gains and the coefficients of the polynomials of the zeros of the closed-loop TF. In the case of the characteristic polynomial, the inner product is added to its coefficients. This way, the mechanism of gain adjustment is represented for allocations of zeros or poles.

From Eq. (17), the equation system that has an unknown vector Kpid and design specifications ais, i=1,2,,n+1 is assembled. In scalar form, this system of equations is represented by

ai+Kpidb¯k=aisKpidb¯k=aisaiKpidb¯k=aie.E19

where b¯k vector is assembled with the rows of the B¯ matrix.

Expanding the scalar representation of Eq. (19), the system of equations to be solved is given by

Kpidb¯k=aieKDb0+KP0+KI0=a1eKDb1+KPb0+KI0=a2eKDb2+KPb1+KIb0=a3eKDbm+KPbm1+KIbm2=a4eKDb0+KPbm+KIbm1=a5e++=KD0+KP0+KIb0=aneE20

The formulation of the problem presented in Eq. (19) and expanded in Eq. (20) is the starting point for the development of forms of parametric variation problems of TFs, as well as, for the establishment of operational points.

To determine the numerical values of the parameters Kpid, the following rules are presented: rule-1) the B¯ matrix is assembled via Eq. (18), where bk is the coefficients of TF numerator polynomial; rule-2) the new ais parameters of the characteristic polynomial are specified; rule-3) the dot product of the parameters of the gain vector Kpid is made with the rows of the matrix B¯, associated with the original ai parameters of the characteristic polynomial and with the specified ais parameters; rule-4) the system of equations given in Eq. (20) and rule-5 the system of equations given in rule-4 (Eq. (20)) is solved to determine the numerical values of the parameters of the gains vector Kpid.

## 4. Experiments

The experimental results are evaluated in three plants with mathematical models in terms of TF obtained with real data, being: Plant I of second order, with a zero; Third order plan II, with two zeros and fourth order plant III, with three zeros.

### 4.1 Plant I

Plant I, is a car dumper, which is used to unload solids in bulk, this equipment has the capacity to move up to 4,000 tons per hour (t/h). The general mathematical model of Plant I in TF is given by

GPIGs=b0s2+a1s+a2,E21

where, GPIGs is the TF of Plant I, it is a second order plant with zero at infinity.

The product of the TF numerator of Plant I given in Eq. (21) associated with the TF numerator of the controller given in Eq. (2) is given by

CPIKpidsGPIGsN=KDs2+KPs+KIb0=KDb0s2+KPb0s+KIb0.E22

The product of the TF denominator of Plant I given in Eq. (21) associated with the TF denominator of the controller given in Eq. (2) is given by

CPIKpidsGPIsD=ss2+a1s2+a2s=s3+a1s2+a2s.E23

The characteristic polynomial of Plant I is given by

PIs=s3+a1+KDb0s2+a2+KPb0s+KIb0.E24

System equations of Plant I related to Eq. (19) in the form Ax=b is given by

Kpidb¯i=aiea1+KDb0=a1s;a2+KDb0=a2s;a3+KDb0=a3s;KDb0=a1sa1;KPb0=a2sa2;KIb0=a3sa3;1)KDb0=a1e;2)KPb0=a2e;3)KIb0=a3e.E25

Placing the systems of equations given in (25) in matrix form, we have

KpidB¯=aieb0000b0000b0×KDKPKI=a1ea2ea3e.E26

The transfer function of Plant I related to Eq. (21) is given by

GPIs=0.4381s2+0.0861s+0.0421s.E27

The ai coefficients of the transfer function of the Plant - I related to Eq. (27) are given by

aia1=0.0861;a2=0.0421;a3=0.E28

The specified coefficients ais of Plant I are given by

aisa1s=0.8604;a2s=0.421;a3s=0.0641.E29

The error calculation aie for Plant I associated with the coefficients ais given in (29) and with the coefficients ai given in (28) is given by

aiea1sa1=a1e=0.86040.0861=0.7743;a2sa2=a2e=0.4210.0421=0.3789;a3sa3=a3e=0.06410=0.0641.E30

The calculation of the Kpid gains vector of Plant I, is done by replacing the numerical values of Eq. (27) in the system of equations given in (25) and (26), is given by.

0.4380000.4380000.438×KDKPKI=0.77430.37890.0641.E31

The Plant I given in Eq. (21) related to Eq. (21) has only the coefficient b0, with that, the system of equations generated, related to the system of equations given in (25), has three equations and three unknowns.

Kpidi)0.438KD=0.7749KD=0.7749/0.431KD=1.78;ii)0.438KP=0.3789KP=0.3789/0.431KP=0.87;iii)0.438KI=0.0641KI=0.0641/0.438KI=0.15.E32

Solving the system of equations given in (32), you can start with any of the equations to find the numerical values of KD, KP and KI, since they are independent. With the numerical values of the earnings KD, KP and KI, it replaces in the simulator developed in the MATLAB/SIMULINK software to monitor the performance of the proposed method.

Figure 2 shows the performance of the PID-Specified controller, which has the transfer function parameters specified by the designer and the KSpecifiedpid gain vector determined by the internal product of the vector of gains with the propagation matrix in purchase with the controller with the gains determined by the second method of ZN.

### 4.2 Plant II

Plant II, is a solid bulk reclaimer, which is used to recover bulk for ship loading, this equipment has the capacity to move up to 8,000 tons per hour (t/h).

GPIIs=b1s+b0s3+a1s2+a2s+a3,E33

where GPIIGs is the TF of Plant II, it is a third order plant with zero at infinity.

The product of the TF numerator of Plant II given in Eq. (33) associated with the TF numerator of the controller given in Eq. (2) is given by

CPIIpidsGPIIsN=KDs2+KPs+KIb1+b0=KDb1s3+KDb0+KIb1s2+KIb0.E34

The product of the TF denominator of Plant II given in Eq. (33) associated with the TF denominator of the controller given in Eq. (2) is given by

CPIIpidsGPIIsD=ss3+a1s2+a2s+a3=s4+a1s3+a2s2+a3s.E35

The characteristic polynomial of Plant II is given by

PPIIs=s4+a1+KDb1s3+a2+KDb0+KPb1s2+(a3+KPb0+KIb1s+KIb0.E36

System equations of Plant II related to Eq. (19) in the form Ax=b is given by

Kpidb¯i=aiea1+KDb1=a1s;a2+KDb0+KPb1=a2s;a3+KPb0+KIb1=a3s;a4+KIb0=a4s;KDb1=a1sa1;KDb0+KPb1=a2sa2;KPb0+KIb1=a3sa3;KIb0=a4sa4;1)KDb1=a1e;2)KDb0+KPb1=a2e;3)KPb0+KIb1=a3e;4)KIb0=a4e.E37

Placing the systems of equations given in (37) in the matrix form, we have

KpidB¯=aieb000b1b000b1b000b1×KDKPKI=a1ea2ea3ea4e.E38

The transfer function of Plant II related to Eq. (33) is given by

GPIIs=0.1812s+0.087s3+0.3853s2+0.117s+0.01567.E39

The ai coefficients of the Plant TF - II related to Eq. (39) are given by

aia1=0.3853;a2=0.117;a3=0.01567;a4=0.E40

The specified coefficients ais of Plant II are given by

a4eisa1s=1.7133;a2s=0.8542;a3s=0.2670;a4s=0.0478.E41

The error calculation aie for Plant II associated with the coefficients ais given in (41) and with the coefficients ai given in (40) is given by

a4eiea1e=a1sa1=1.71330.3853=1.328;a2e=a2sa2=0.73720.117=0.6202;a3e=a3sa3=0.25140.01567=0.2358;a4e=a4sa4=0.04780=0.0478.E42

The calculation of the Kpid gains vector is performed by replacing the numerical values of Eq. (39) in the system of equations given in (37) and (38).

0.087000.18120.087000.18120.087000.1812×KDKPKI=1.32800.95380.09960.0478.E43

The plant II given in Eq. (39) related to Eq. (33) has the b0 coefficients and b1, with this, the system of equations generated, for the system of equations related to the system of equations given in (37), has four equations and three unknowns.

Kpidi)KDb1=a1eKD=1.328/0.1812KD=7.329;ii)KDb0+KPb1=a2eKP=0.95380.6376/0.1812KP=1.745;iii)KPb0+KIb1=a3eKI=25140.1518/0.1812KI=0.549;iv)KIb0=a4eKI=0.0478/0.087KI=0.549.E44

Solving the system of equations given in (44), first, solve Equation i) to find the numerical value of KD. Then, replace the numerical value of KD in Equation ii) and find the numerical value of KP. To find the numerical value of KI, solve equation iv) or replace the values of KP in equation iii). With the numerical values of the gains KD, KP and KI, it replaces in the simulator developed in the MATLAB/SIMULINK software to monitor the performance of the proposed method.

Figure 3 shows the performance of the PID-Specified controller, which has the transfer function parameters specified by the designer and the KSpecifiedpid gain vector determined by the internal product of the vector of gains with the propagation matrix in purchase with the controller with the gains determined by the second method of ZN.

### 4.3 Plant III

Plant III, is a car dumper with two feeders, which is used to unload solids in bulk, this equipment has the capacity to move up to 8,000 tons per hour (t/h). The general mathematical model of Plant III in TF is given by

GPIIIGs=b2s2+b1s+b0s4+a1s3+a2s2+a3s+a4.E45

where, GPIIIGs is the TF of Plant III, it is a fourth order plant with zero at infinity.

The product of the TF numerator of Plant III given in Eq. (45) associated with the TF numerator of the controller given in Eq. (2) is given by

CPIIIpidsGPIIIsN=KDs2+KPs+KIb2s2+b1s+b0=KDb2s4+KDb1+KPb2s3.E46
+KDb0+KPb1+KIs2+KPb0+KIb1s+KIb0.

The product of the TF denominator of Plant III given in Eq. (45) associated with the TF numerator of the controller given in Eq. (2) is given by

CPIIIpidsGPIIIsD=ss4+a1s3+a2s2+a3s+a4=s5+a1s4+a2s3+a3s2+a4s.E47

The characteristic polynomial of Plant III is given by

PPIIIs=s5+a1+KDb2s4+a2+KDb1+KPb2s3+a3+KDb0+KPb1+KIb2s2+a4+KPb0+KIb1s+KIb0.E48

System equations da Planta III related to Eq. (19) in the form Ax=b is given by

Kpidb¯i=aiea1+KDb2=a1s;a2+KDb1+KPb2=a2s;a3+KDb0+KPb1+KIb2=a3s;a4+KPb0+KIb1=a4s;a5+KIb0=a5s.KDb2=a1sa1;KDb1+KPb2=a2sa2;KDb0+KPb1+KIb2=a3sa3;KPb0+KIb1=a4sa4;KIb0=a5sa5.1)KDb2=a1e;2)KDb1+KPb2=a2e;3)KDb0+KPb1+KIb2=a3e;4)KPb0+KIb1=a4e;5)KIb0=a5e.E49

Placing the systems of equations given in (49) in the matrix form, we have

KpidB¯=aieb000b1b00b2b1b00b2b100b2×KDKPKI=a1ea2ea3ea4ea5e.E50

The transfer function of Plant III related to Eq. (45) is given by

GPIIIs=0.959s2+0.1698s+0.1593s4+0.1767s3+0.3463s2+0.029s+0.02331.E51

The ai coefficients of the Plant TF - III related to Eq. (51) are given by

aia1=0.1767;a2=0.3463;a3=0.029;a4=0.2331;a5=0.E52

The specified coefficients ais of Plant III are given by

aisa1s=1.8358;a2s=8.619;a3s=8.9947;a4s=1.4277;a5s=1,3254.E53

The error calculation aie for Plant III associated with the coefficients ais given in (53) and with the coefficients ai given in (52) is given by

aiea1e=a1sa1=1.83580.1767=1.659;a2e=a2sa2=8.6190.3463=4.9066;a3e=a3sa3=8.99470.029=8.9657;a4e=a4sa4=1.42770.0233=1.4044;a5e=a5sa5=1.32540=1.3254.E54

The calculation of the gain vector Kpid is done by replacing the numerical values of Eq. (51) in the system of equations given in (49) and (50).

0.1593000.16980.159300.9590.16980.159300.9590.1698000.959×KDKPKI=1.65914.90669.07122.73811.3254.E55

The Plant III given in Eq. (51) related to Eq. (45), has the coefficients b0, b1 and b2, with that, the generated system of equations, for the system of equations related to the system of equations given in (49), has five equations and three unknowns.

Kpidi)KDb2=a1eKD=1.6591/959KD=1.73;ii)KDb1+KPb2=a2eKP=4.6128/0.959KP=4.81;iii)KDb0+KPb1+KIb2=a3eKI=7.978/0.959KI=8.32;iv)KPb0+KIb1=a4eKI=1.4127/0.1698KI=8.32;v)KIb0=a5eKI=1.3254/0.1593KI=8.32.E56

Solving the system of equations given in (56), first, solve Equation i) to find the numerical value of KD. Then, replace the numerical value of KD in Equation ii) and find the numerical value of KP. To find the numerical value of KI, solve Equation v) or replace the values of KD and KP in Equation iii) or you can substitute the value of KP in Equation iv). With the numerical values of the gains KD, KP and KI, it replaces in the simulator developed in the MATLAB/SIMULINK software to monitor the performance of the proposed method.

Figure 4 shows the performance of the PID-Specified controller, which has the transfer function parameters specified by the designer and the KSpecifiedpid gain vector determined by the internal product of the vector of gains with the propagation matrix in purchase with the controller with the gains determined by the second method of ZN.

## 5. Conclusions

The study presented a methodology for adjusting the gains of a PID controller in terms of the internal product of the gains and the propagation matrix. In addition, the relevance of the matrix was shown, which enabled impact assessment of the PID actions associated with the plant parameters. The proposed methodology complied with the project specifications and ensured high controller efficiency without suppressing the PID terms caused by the adjustments. The three PID controller terms were adjusted. Therefore, this methodology can be considered as an alternative to conventional methods for the computation of KD, KP, and KI gains parameters in practical applications and highly complex control plants.

## Acknowledgments

The authors would like to thank PPGEE of the UFMA for the resources to develop this work. We are especially grateful to FAPEMA for research incentive and infrastructure. We acknowledge the Department of Computer Engineering of the UEMA for making this research possible. Finally, we also acknowledge CAPES and CNPq for promoting and supporting the advanced studies that contributed to this work.

## References

1. 1. Sutikno Juwari Purwo and Hidayah Nur and Handogo Renanto. Maximum Peak-Gain Margin (Mp-GM) Tuning Method for Two Degree of Freedom PID Controller. ID Control for Industrial Processes. BoD–Books on Demand, 2018
2. 2. Bucz Štefan and Kozáková Alena. Advanced methods of PID controller tuning for specified performance. PID Control for Industrial Processes. IntechOpen London, 2018
3. 3. Jeng Jyh-Cheng. Data-based tuning of PID controllers: a combined model-reference and VRFT method. PID Control for Industrial Processes. BoD–Books on Demand, 2018
4. 4. Moura J P, Fomseca Neto J V, Rêgo P H M. A Neuro-Fuzzy Model for Online Optimal Tuning of PID Controllers in Industrial Systems Applications to the Mining Sector. IEEE - Transactions on Fuzzy Systems, 2019
5. 5. Moura J P, Fomseca Neto, Ferreira E F M, Araujo Filho E. M. On the Design and Analysis of Structured-ANN for Online PID-Tuning to Bulk Resumption Process in Ore Mining System. Elsevier - Neurocomputing, 2020
6. 6. Chen, Chi-Tsong. Linear System Theory and Design. Inc. 3rd. Oxford University Press, New York, NY, USA, 199
7. 7. Åström K J, Wittenmark B. Computer-Controlled Systems: Theory and Design, Third Edition, isbn: 9780486284040, series: Dover Books on Electrical Engineering. Dover Publications, 2013
8. 8. Vilanova R, Visioli A. PID Control in the Third Millennium: Lessons Learned and New Approaches, isbn: 9781447124252, series: Advances in Industrial Control. Springer London, 2012
9. 9. Tavakoli, Ali Reza and Seifi, Ali Reza. Adaptive self-tuning PID fuzzy sliding mode control for mitigating power system oscillations. Neurocomputing, volume 218, pages 146–153. Elsevier, 2016
10. 10. Liu P, Gu, Haibo, Kang, Yu and Lü, Jinhu. Global synchronization under PI/PD controllers in general complex networks with time-delay. Neurocomputing. Elsevier, 2019
11. 11. Aidan O’Dwyer. Handbook Of Pi And Pid Controller Tuning Rules, isbn: 1848162421,9781848162426. Imperial College Press, 2009
12. 12. Alfaro, V.M. and Vilanova, R. Model-Reference Robust Tuning of PID Controllers, isbn: 9783319282138, lccn: 2016935592, series: Advances in Industrial Control. Springer International Publishing, 2016
13. 13. S. Bennett. A brief history of automatic control. doi 10.1109/37.506394, issn: 1066-033X. IEEE Control Systems, volume 16, number 3, pages 17–25, 1996
14. 14. Lee, Jay and Kao, Hung-An and Yang, Shanhu. Service innovation and smart analytics for industry 4.0 and big data environment. Procedia Cirp, volume 16, pages 3–8. Elsevier, 2014
15. 15. Lasi, Heiner and Fettke, Peter and Kemper, Hans-Georg and Feld, Thomas and Hoffmann, Michael. Industry 4.0. Business & Information Systems Engineering, volume 6, number 4, pages 239–242, Springer, 2014
16. 16. Marwedel, P. Embedded System Design: Embedded Systems Foundations of Cyber-Physical Systems, and the Internet of Things, isbn: 9783319560458, series: Embedded Systems. Springer International Publishing, 2017
17. 17. Kopetz, H. Real-Time Systems: Design Principles for Distributed Embedded Applications, isbn: 9781441982377, series: Real-Time Systems Series. Springer US, 2011
18. 18. Fadali, M.S. and Visioli, A. Digital Control Engineering: Analysis and Design, isbn: 9780123983244, lccn: 2012021488. lsevier Science, 2012
19. 19. Quddious Abdul and Antoniades Marco A and Vryonides Photos and Nikolaou Sym. Voltage-Doubler RF-to-DC Rectifiers for Ambient RF Energy Harvesting and Wireless Power Transfer Systems. Wireless Energy Transfer Technology. IntechOpen, 2019
20. 20. Moura, José Pinheiro, Fonseca Neto, João Viana and Rêgo, Patrícia Helena M. Models for Optimal Online Tuning Based on Computational Intelligence of PID Controllers Applied to Operational Processes of Bulk Reclaimers, issn: 2195–3899, doi:10.1007/s40313-018-00438-7. Journal of Control, Automation and Electrical Systems. Springer, 2019
21. 21. Moura, José Pinheiro and Neto, João Viana Fonseca and Ferreira, Ernesto Franklin Marçal and Araujo Filho, Evandro Martins. On the Design and Analysis of Structured-ANN for Online PID-Tuning to Bulk Resumption Process in Ore Mining System. Neurocomputing. Elsevier, 2020

Written By

José Pinheiro de Moura and João Viana da Fonseca Neto

Submitted: July 15th, 2020 Reviewed: November 17th, 2020 Published: January 5th, 2021