Design of PID Controllers for Decoupled MIMO Systems

DOI : 10.17577/IJERTCONV2IS04012

Download Full-Text PDF Cite this Publication

Text Only Version

Design of PID Controllers for Decoupled MIMO Systems

Vaishali Killarikar

Assistant Professor, Electronics Department KCCEMSR , Thane (E)

Jayshree Katkar,

Assistant Professor, Electronics & Tele. Department KCCEMSR , Thane (E)

Abstract- The design of PID controllers for systems with interacting loops is discussed. It is important to deal with the interaction at the lower-level loops, since supervisory control based on for instance MPC seldom has sufficient bandwidth. A new scheme based on modified scalar PID design and static decoupling is developed, where the frequency characteristics of the coupling between the lower-level loops is taken into account. This leads to a design method emphasizing the trade-off between the individual loop performances and the so called interaction indices. The controller is easily implemented, due to its simple configuration based on standard components.

Keywords: PID controllers, multivariable control systems, decoupling


Model predictive control (MPC) is becoming the standard technique to solve multivariable control problems in the process industry [12, 5,8 ]. Practically all MPC systems are however operating in a supervisory mode with PID controllers at the lower level. A substantial portion of the performance improvement credited toMPC is actually due to improvements in the lower-level PID loops. Interaction among the loops causes difficulties when the lower-level loops are closed. There are some difficulties in dealing with the interaction at the MPC level because the bandwidths of the MPC loops are limited: they operate in supervisory mode with sampling intervals that are longer than the PID loops.It is consequently of interest to investigate ways of dealing with interaction at the loop level. The approach we take is to investigate standard PID tuning

  1. and see what can be achieved by adding simple interactions between the feedback loops. In many cases the performance of the system can improved, particularly if the coupling in the process is not severe. The proposed scheme is based on a simple decoupling, which implies that it can be easily implemented at the loop level. The advantage by doing this is that it gives

    performance enhancement in a frequency range that is normally not dealt with by MPC.

    1. THE DESIGN PROCEDURE The PID controller is described by

      where K is the proportional gain, Ti is the integral time and Td is the derivative time.

      The design of a decentralized control system with a decoupling matrix can be done combining a diagonal controller Kd (s) with a block compensator D(s), sothat the controller manipulates the variable ui insteadof the ui, as can be appreciated in figure 1, for the 2 x2 case.With this configuration the controller see theprocess as a set of n completely independent processes or with the interaction minimized. The essence of decoupling is the imposition of a calculation net that cancels the existent process interaction, allowing the independent control of the loops. In decentralized design, the question is not to eliminate interaction, but to take it into account. A multivariable system may still experience interactions and responds poorly. The objective in decoupling is to compensate for the effect of interactions brought about by cross coupling of the process variables.

      In literature there are different decoupling methods: Lineal decoupling (Desphande, 1989) is most extended method. In this, the decoupling matrixes try to eliminate interactions from all loops, obtaining following elements for a 2 x 2 system

      Implementation of this decoupling matrix has some problems: What happens if numerator has bigger order

      than denominator?, what if delays exist?, and what if these delays only appear in the denominator? The different works propose some solutions like Pade approximations for delays or steady state decoupling as a first option to prevent this problems. Other solutions consist of applying partial decoupling, setting d12 or d21 null. This fact avoid problems originated in one loop reach the other.

      There are some other methods: one of them looks for diagonal dominance in the system as ALIGN algorithm designed by McFarlane and Kouvaritakis in 1974 and described in (Maciejowski, 1989). Other design by means of singular value decomposition (Desphande, 1989) or by means of inverse decoupling (Wade, 1997) that produces the process input signals by combining one controller output with the other process input signals.Even if decouplers are incorporated, the interaction effects cannot be completely eliminated because of model mismatch. Then, single-loop controllers cannot be tuned independently, and a sequential tuning algorithm that takes interactions into account should be used. Next section describes one of these algorithms.

    2. THE METHOD:

      Analysis of decentralized methods described in section 1 show several interesting points:

      • Obtaining different design specifications for each loop implies solution methods based on some kind of iteration. If designs are carried out loop-by-loop, tuning one of them can detune the others.

      • In spite of abundance of methods, there are no simple and general solutions for the multivariable tuning problem. Or methods are used by optimization (Wang et al, 2000) that do not guarantee an appropriate solution and they require several iteration processes, or the solutions are too particular, for concrete models(Ho et al, 1996) or too simple (Shinskey, 1996),(Luyben, 1992) that only get a first approach to the problem.In

        (Vázquez et al, 1999), a method of tuning PID controllers for systems with decentralized control is presented. Its fundamental characteristics can be summarized in following points:

      • It is a generic method for n x n systems with decentralized control.

      • It is a method based on successive SISO tuning that does not suppose any additional constraintto obtained controllers (except for their decentralized structure) neither to the transfer functions matrix. This matrix can include the transfer functions of only the process or also including the decoupling net. Present work exploits this characteristic.

      • The SISO methodology integrated in the algorithm has certain imposed limitations: it should allow the controller be designed from a frequential description of the process (their frequency response). And it should quantify the specifications achievement by means of a quality index, J. This index controls the iteration evolution of the tuning algorithm.

      • The method allows a decoupling matrix to be included between the plant and the controllers.

      • The design methodology is divided into two phases that are carried out in a sequential and combined way: the structural decomposition and the controller design. These two phases are described next.

      1. System structural decomposition

        This phase consists of the decomposition of a n x n multivariable system into n SISO systems. To solve this problem the structural decomposition, introduced by Zhu (1996), is used. Let be a n input system controlled by means of a decentralized control strategy. Now, n-1 loops have been closed by means of n-1 controllers. The process seen from the free input to the free output is needed. The scheme corresponds to the figure 2, where loop i has not yet been closed (the pairing problem has been already solved and input-output pairs are in the diagonal).

        In the structural decomposition schme, K1 is the controller that closes the loop between input i and output i. The set of n-1 controllers, whose loops are already been closed, is K2. This way, the controller matrix is:

        The process element gii has its feedback loop open. The rest of elements of its same row are called G12. The rest of elements of its same column are called G21. And the rest of elements of G are called G22. Then, the process matrix G can be written as:

        The process seen from input i to output i when the other loops has been closed is given by

        The structural decomposition implies some important considerations for the decentralized control: some global system properties, as interaction, stability, integrity, etc, could be deduced from the properties of the obtained n SISO subsystems. Also, no supposition about G has been made, and then, structural decomposition could be applied to every transfer functions matrix. Applying Nyquist theorem to each individual loop obtained with structural decomposition, stability criterion is more interesting from the point of view of applicability to decentralized multivariable systems. Then, the n SISO systems stability implies MIMO stability

        Theorem 1: Supposing that individual elements of G(s) and their SISO independent subsystems do not have poles in the right hand plane, the system with

        decentralized control is stable if and only if the Nyquist contour of the equivalent open loop transfer function i ig

        ~ k does not encircle point (-1,0), i.e.

        The structural direct Nyquist arrays (SDNA) are the representation of the n direct Nyquist diagrams of

      2. Controllers design

      The proposed is a frequency-based method and stability measures as phase and gain margins could be used. But most of existent methods requires a model for obtaining the design parameters. Carrying out the structural decomposition, the resulting transfer function i g ~ could be not a rational function, possibly due to the inclusion of delays in model systems. In order to solve this problem a frequential representation of this transfer function i g ~ it is used, i.e. a frequency response array. But this take implicit another problem: what SISO design method to use? The particular solution has been the mono variable methodology described in (Morilla and Dormido, 2000). This methodology supposes, and in fact requires, a model with the description of the system dynamics. However it leaves open the possibility to be adapted to solve the problems outlined previously it can be modified allowing the controllers to be designed from an array with the frequential representation of the frequency response and also, the obtained results can be compared with the design specifications to present an indicative of quality, indispensable for the proposed solution search method. These two phases cannot be considered isolated, since the design process consists of an iteration of both. This is because, as equation (4) shows, the gi transfer functions, obtained with structural decomposition, incorporate the controller transfer functions in their definitions. Thus, every time a design is done for one of these transfer functions, it will be needed to recalculate the others.From the new i gi, a new controller could be designed, and then, definitely, the algorithm needs an alternation between the two phases


      One of advantages that present this methodology of multivariable controller design is to be independent of the process model. This has been shown in numerous examples in (Vazquez et al, 1999), which can be repeated with TITO tool. Also, the model does not have reason to be a rational expression in the Laplace operator, but rather it can be a frequency response array. These characteristics allow the immediate extension of this method, using it accompanied by some decoupling

      strategies (one of the described in section 2 or another), so that in a first step the decoupling matrix is obtained and, in a second step, the controllers for the decoupler+process block designed. Some of the possibilities are shown in following examples.

      Example 1: Following shows a system model with a RGA next to 1.7. Some decoupling strategy could be tested. It is a distillation column (Vinante and Luyben, 1972) that describes the existent dynamics between reflux and vapor flow and the temperatures of plates 4 and 17:

      In order to analyze the effects of designing with decouplers, expressions (1) are used. In this case,dynamic and steady state decoupling coincide:

      Design specifications are phase margin (PM) of 45º and gain margin (GM) of 4 for both loops, and looking for a PID with a limit of the relationship between the derivative and integral constant equal to 0.01. The design using iterative algorithm without decouplers gets following controllers:

      Design specifications are met with combined (PM and GM) tuning design in three iterations. Interaction impedes using other methods based on Gershgorin bands, as Hos method Next, a complete decoupling is chosen. Design specifications are the same as before (PM=45º and GM=4 for both loops). The algorithm converges insame number of iterations as before, obtaining following controllers

      With this design, specifications are also met, but now, controlled system was decoupler + process block. Time responses are shown in figure 3. Response without decoupler is also superimposed.

      A slight modification can be appreciated, mainly in second loop, and an interaction decrease on the first one.Although modifications in time response are minimum, design is good to show methodology effectiveness when decouplers are incorporated.

      Example 2: Table 1 shows the results of PID controller design for the system

      for different decoupling nets (dynamic, in steady state, partial and total) and for different controllers (PI, PID). All the designs have been carried out with specifications of PM=60º and GM=4 for both loop. It is a water- methanol distillation column (Wood and Berry, 1973) analyzed in numerous later works (Ho et al, 1996), (Toh and Devanathan, 1993), (Vázquez et al, 1999) (it appears in the practical whole of multivariable control references).


      Figure 4 shows superimposed time responses of the system without decouplers and with total dynamic decouplers, with a PI and a PID (the two last designs of table 1). The PID design has been carried out limiting

      ,to 0.25

      Note that interaction effects decrease not very significantly. This same reduction could be obtained tuning the second loop in a less aggressive way (for example PM=80º), which would reduce its effect on first loop.

      Example 3: This third example studies a process proposed by Niederlinski (1971) and analyzed with a set of decouplers proposed by Shiu and Hwang (1998). The process function transfer matrix is:

      Shiu proposes the following decoupler, obtained after identifying the i g ~ functions applying a relay method:

      In order to compare time responses with and without decoupler, the tuning algorithm has been used with the same design specifications (PM = 60º for both loops). Without decoupler, controllers are

      obtained in six iterations of the algorithm. With decouplers, the controller equations are

      obtained in only three iterations, because, in general,the algorithm converges quicker the more dominant the system is. Figure 5 shows time response, with and without decouplers. Note that interaction decreases. If this fact is quantified by means of some measure as the IAE test, it can be proven that in the response without decouplers, IAE is 3.42, while when decouplers are applied it is 0.98.


The suggested control scheme is a quite simple multivariable controller. However, it handles a large number of practical control problems, although it can be implemented with regular PID controllers together with a proportional controller D. It is desirable to have efficient dsign methods for decoupled PID control. Such a method has been developed in the paper, under the assumption that the interaction is not too severe. It takes advantage of methodology proposed in (Vázquez et al, 1999), and it does not need a transfer function matrix with the system model but a representation with frequency response. Independently of how decouplers

had been obtained, they can be included between process and controllers. Then, a decentralized PID controller is designed, using some of the different possibilities of proposed algorithm: design with only phase margin specifications for both loops, only gain margin, and a combination of PM and GM specifications.


  1. K.J.Astrom.Limitations on control system performance. European Journal on Control, 6:2-20, 2000.

  2. K. J. Wstrom and T. Hagglund. PID Controllers:Theory, Design, and Tuning. Instrument Society of America, Research Triangle Park, NC, 1995.

  3. K. J. bstrom, H. Panagopoulos, and T. Hagglund. Design of PI controllers based on non-convex optimization. Automatica, 34(5):585-601, 1998.

  4. E. Bristol. On a new measure of interaction for multivariable process control. IEEE lhns. Autom.Control, 11:133, 1966.

  5. E. F. Camacho and C. Bordons. Model Predictive Control. Springer, London, 1999.

  6. W. L. Luyben. Process modeling, simulation, and control for chemical engineers. McGraw-Hill, New York, 1990.

  7. J. M. Maciejowski. Predictive Control with Constraints. Pearson Education, 2000.

  8. A. Niederlinski. A heuristic approach to the design of h e a r multivariable interacting control systems. Automatica, 7:691- 701, 1971.

  9. H. Panagopoulos. PID-Control. Design, Eztension, Application. PhD thesis, Dept. of Automatic Control, Lund University, Sweden, 2000.

  10. P. Persson and K. J. Astrom. Dominant pole design – a unified view of PID controller tuning. In Preprints 4th IFAC Symposium on Adaptive Systems in Control and Signal Processing, pages 127-132, Grenoble, France 1992.

  11. Deshpande, B. (1989): Multivariable process control. Instrument Society of America.

  12. Halevy, Y, Palmor, Z., Efrati,T. (1997): Automatic tuning of decentralized PID controllers for MIMO processes, J. Proc. Cont. Vol. 7 Nº 2, pp 119-128 Elsevier Science, Ltd.

  13. Ho, W.K., T.H. Lee, O.P. Gan (1996): Tuning of multiloop PID controllers based on gain and phase margins specifications. 13th IFAC World Congress, pp 211-216.

  14. Morilla, F.; Dormido, S. (2000): Methodologies for the tuning of PID controllers in the frequency domain. Proceedings of PID00. Tarrasa, pp.155-160

  15. Vázquez, F.; Morilla, F.; Dormido, S. (1999): An iterative method for tuning decentralized PID controllers. Proceeding of the 14th IFAC World Congress, pp. 491-496.

Leave a Reply