Disturbance rejection for multivariable processes with time delays based on partial inverted decoupling control

Abstract

Motivated by the thought of using couplings to reject disturbances, the authors propose a novel partial inverted decoupling control for improving the disturbance rejection of multivariable processes with time delays. With the relative disturbance gain (RDG) for n × n processes, whether the couplings among loops are favorable is evaluated. The inverted decoupling usually used in full decoupling control is then extended to partial decoupling control to retain the favorable couplings and remove the unfavorable ones. The elements of the partial inverted decoupler are derived and the implementation structure is provided. Based on the partially decoupled processes, the decentralized controllers for both the decoupled and non-decoupled loops are designed. An analytical design method for the decentralized controllers is proposed to improve disturbance rejection. Controller parameters are determined in terms of both performance and robustness. Two simulation examples show that the presented partial inverted control strategy can provide superior disturbance rejection to several other classical multivariable control strategies.

Keywords

Multivariable process time delay partial inverted decoupling disturbance rejection

Introduction

In process control industries, most controlled processes are multivariable processes in nature. An important characteristic of multivariable processes is the couplings between the inputs and outputs.¹ Couplings make the control problem for multivariable processes much more complicated than that for univariate processes. Therefore, how to handle couplings is a key issue in multivariable control.

Decentralized control, with no doubt is the most commonly used control strategy for multivariable processes in practice due to its easy implementation.² The primary step for decentralized control is usually to determine the main channel of each loop. Therefore, how to pair the inputs and outputs so as to reduce relative couplings is of great importance³ and many indicators have been proposed, such as relative gain array (RGA)⁴ and some relevant methods with process dynamics considered.^5,6 With the well-paired system, the decentralized controllers can be tuned by detuning method,⁷ auto-tuning method⁸ and iterative design method,⁹ etc. Recently, Mahapatro et al.¹⁰ designed decentralized PI controllers based on frequency response matching of a predefined reference model and the actual model. Mahapatro and Subudhi¹¹ proposed a robust decentralized PID controller based on complementary sensitivity function through a graphical tuning method. For minimizing the undesirable effects of the couplings, Euzébio et al.² optimized decentralized controllers by determining the cost function as the sum of the integral absolute errors (IAE) from each output due to set-point changes in coupled loops. Jeng and Lee¹² proposed a direct synthesis method for decentralized controllers on the basis of specifying the frequency response of the desired closed-loop transfer function matrix. Apparently, the decentralized controllers designed only with the main channels may not achieve sufficiently satisfactory control performances since the coupling information is ignored. As a result, the effective open-loop transfer function (EOTF) was proposed for independent design of multi-loop controllers.¹³ In the formulation of the EOTF, the coupling channels were contained. With the concept of EOTF, the multivariable control system can be decomposed into several equivalent single-loops for decentralized controller design.^14,15 However, since the calculation of EOTF is complex, especially for high-dimensional processes, model approximation is usually inevitable to simplify controller design. Approximation errors will affect the control performance. In sum, decentralized control cannot fully remove the effect of couplings in spite of its simple structure, which means it is suitable for the process with modest couplings.

When significant couplings exist, an intuitional way to improve control performance is to remove or reduce the loop couplings by inserting a decoupler. When the process is well decoupled, the control strategies for single-input-single-output (SISO) processes can be used directly. Hence, the design of the decoupler is the key issue in decoupling control. Generally, decoupling strategies are classified into two categories, that is, static decoupling and dynamic decoupling.¹⁶ A static decoupler is designed by only considering the steady state of a process. It has the advantages of simple calculation and easy implementation while it may not provide satisfactory decoupling performance because of the neglected process dynamics. A dynamic decoupler can achieve better decoupling performance but it needs more accurate process information. As the most representative multivariable processes, many dynamic decouplers are firstly studied for 2 × 2 processes. Among them, the most classical three methods are ideal decoupling, simplified decoupling and inverted decoupling.¹⁷ Ideal decoupling which uses the inverse of the process to design the decoupler is rarely used in practice because it suffers from complex decouplers. Compared with ideal decoupling, simplified decoupling is more popular as it provides a simple and implementable decoupler. However, a simplified decoupler usually results in a complex decoupled process, which is not desirable for controller design. This will be more evident when simplified decoupling is extended to high-dimensional processes.¹⁸ Unlike the above two, inverted decoupling can achieve simple decouplers and decoupled processes simultaneously, but the limitations of achieving stable and realizable decoupler elements should be considered when it is applied for high-dimensional processes.¹⁹ To overcome the effect of the mismatched model, Aguiar et al.^20,21 proposed closed-loop evaluation and redesign methodologies for simplified decoupling and inverted decoupling. In essence, these dynamic decoupling methods seek to make the decoupled process matrix a diagonal matrix at more frequencies. Motivated by this, Jin and Liu²² use the Nyquist set to describe the dynamics of the decoupler and use a complex-curve fitting technique to obtain an approximate decoupler matrix. Moreover, Xiong et al.²³ proposed the equivalent transfer function (ETF) based on the inherent relation between the inverse process and the dynamic relative gain array (DRGA) to simplify the design of the decoupler. Naik et al.²⁴ approximated the ETF using the static relative gains and the average residence time. Khandelwal and Detroja²⁵ presented a novel ETF parameterization procedure for large scale processes or systems with higher-order elements. To simplify the structure of multivariable system, the multivariable controller can have the abilities of both decoupling and control, which is called the centralized controller. Usually, a centralized controllers can be realized by combining a decoupler and a set of decentralized controllers.²⁶ Besides, a centralized controllers can also be synthesized directly based on a desired closed-loop transfer function matrix.²⁷ Garrido et al.²⁸ presented an iterative methodology for designing centralized PID controllers based on the proposed equivalent loop transfer functions. Gündeş²⁹ gave the necessary and sufficient conditions for the existence of decoupling controllers. Actually, the couplings can be regarded as the internal disturbances of the multivariable systems, which means the methods for handling couplings is to suppress the internal disturbances.

Apart from couplings, external disturbances and process uncertainties are inevitable in many industries, which makes the primary goal of the controllers to mitigate the effects of external disturbances and process uncertainties.^30,31 Feedback control is an important strategy for suppressing disturbances. As the classical feedback controller, PID controller still dominates the engineering systems. The key issue of PID controller is to tune the controller parameters. The traditional tuning formula, such as Z-N method, are experiential methods, which are easy to implement while only achieve limited control performance. When the process model is known, the synthesis of PID controller can be transformed to solving the optimization problem with given objectives and constraints.^32,33 With the development of machine learning techniques, some are incorporated to PID tuning to improve control performance such as neural network PID,³⁴ fuzzy PID³⁵ and swarm intelligent PID³⁶, etc. To reduce the dependence of accurate models, data-driven techniques, like iterative feedback tuning³⁷ and virtual reference feedback tuning,³⁸ are employed to tune fix-structure controller, which is suitable for the processes with big uncertainties. Apart from feedback control, feedforward control is another effective scheme for disturbance rejection, which cancels disturbances with the reverses of disturbance signals. In light of the main drawback of feedforward control that it can only deal with measurable disturbances, many disturbance observers are designed to estimate the disturbances by using the measurable states and the process dynamics, which results in the disturbance observer-based (DOB) control. DOB control contains several related methods, such as frequency domain DOB control,³⁹ active disturbance rejection control (ADRC),⁴⁰ equivalent input disturbance (EID),⁴¹ etc. Although under different names, these methods share the same basic idea, that is, the lumped disturbances including the external disturbances and uncertainties are estimated and corresponding compensations are generated based on the estimations. Due to realizing the compensation of uncertainties, DOB control does not require an accurate process model.³¹ Meanwhile, its inherent 2 degree-of-freedom (DOF) structure allow the DOB control system to achieve a satisfactory performance in both set-point tracking and disturbance rejection. As a result, DOB control has been applied into more and more industries.^42–44

As the typical characteristic of the multivariable system, the couplings can also be used for disturbance rejection when they have the opposite effect of the external disturbances. Motivated by this, Stanley et al.⁴⁵ defined the relative disturbance gain (RDG) to evaluate couplings for distillation processes. Similarly, Lin et al.⁴⁶ proposed the relative load gain (RLD) as the ratio of the gains of effective disturbances and open-loop disturbances. Moreover, Chang and Yu⁴⁷ systematically described the controller structure via internal model control (IMC) structure and proposed the relative disturbance gain array (RDGA) to synthesize the controller structure straightforwardly. Based on the evaluation of couplings, the partial decoupling structure can be constructed and some control strategies have been studied.^48,49 As well, the thought of using couplings for disturbance rejection has been extended to other control structures.^50,51

Motivated by the thought of using couplings to reject disturbances, this work is aimed to propose a novel partial decoupling control strategy. Unlike the existing methods which cannot remove the unfavorable couplings thoroughly due to the approximation or simplification of the decoupler, the proposed strategy extends the inverted decoupler which used in full decoupling control to partial decoupling control. The designed partial inverted decoupler can realize better elimination of the unfavorable couplings. Besides, in the existing methods, the controllers of the non-decoupled loops are usually designed based on the main channel models and disturbance channel models. The dynamic characteristics of non-decoupled loops are not analyzed sufficiently. In the proposed method, the equivalent process model and the equivalent disturbance channel model are obtained for designing the controllers of the non-decoupled loops. Finally, to simplify the design of the controllers, an analytical design method for the decentralized controllers is proposed by considering the performance of disturbance rejection and robustness. In sum, the main contributions of this work are concluded as follows:

(1) The RDG for general n × n processes is explored to evaluate whether the couplings are favorable or not.

(2) The inverted decoupling is extended to partial decoupling control to retain the favorable couplings and remove the unfavorable ones. The elements of the partial inverted decoupler are derived and the implementation structure is provided.

(3) The equivalent process models and the equivalent disturbance channel models are obtained to design the controllers of the non-decoupled loops.

(4) An analytical design method for the decentralized controllers is proposed with the consideration of both the performance of disturbance rejection and robustness.

This paper is organized as follows. In Section “Relative disturbance gain (RDG),” the RDG for general n × n processes is derived to evaluate the couplings. Based on the evaluation result, the design and implementation of the partial inverted decoupler are discussed in Section “Partial inverted decoupling control structure.” In Section “Design of the decentralized controller,” the decentralized controllers for decoupled and non-decoupled loops are designed respectively. An analytical design method is proposed. Two simulations examples are provided in Section “Simulation study.” Section “Conclusion” is the conclusion.

Relative disturbance gain (RDG)

With the purpose of using the couplings between loops to reject disturbances, the first but important step is to judge whether the couplings are favorable. Hence, an index for evaluating the couplings is necessary. With this index, the couplings beneficial to disturbance rejection should be retained while the unfavorable ones have to be decoupled. Herein, the index is chosen as RDG.

The RDG was proposed by Stanley et al.⁴⁵ for 2 × 2 processes, like dual-composition control. Equation (1) shows the definition of the RDG of the ith loop, that is, β_i, where u_i and y_j are the input of the ith loop and the output of the jth loop, respectively; d is a disturbance signal.

β_{i} = \frac{{(\frac{\partial u_{i}}{\partial d})}_{perfect control over y_{j} (j \neq i)}}{{(\frac{\partial u_{i}}{\partial d})}_{no control over y_{j} (j \neq i)}}

(1)

Obviously, β_i is the ratio of the required input changes for rejecting the disturbances in two different cases. The numerator indicates the case where the decentralized control is used for the multivariable system while y_j (j ≠ i) is under perfect control. In this case, the couplings are retained. The denominator indicates the case where the system has been well decoupled since y_j (j ≠ i) has no effect on y_i. Therefore, it can be concluded that the couplings are favorable to disturbance rejection if |β_i| < 1, while harmful to disturbance rejection if |β_i| > 1. According to equation (1), the expression of β_i for n × n processes can be derived.

Figure 1(a) shows the commonly used decentralized control structure for general n × n processes, where C(s) is the decentralized controller, expressed as C(s) = diag(c₁(s), c₂(s), …, c_n(s)); P(s) is the controlled process, whose element is described as p_ij(s), i, j = 1, 2, …, n; L(s) is the disturbance transfer function vector, expressed as L(s) = [L₁(s), L₂(s), …, L_n(s)]^T;d(s) is a specific disturbance signal; r, u and y are the reference vector, input vector and output vector, respectively. Assuming that the loops have been paired, according to Figure 1(a), the ith disturbance response, $y_{i}^{d}$ can be expressed by (the Laplace operator s is omitted)

\begin{matrix} y_{i}^{d} = (p_{ii} - P_{i *} H_{- ii} {(P_{- ii})}^{- 1} P_{* i}) u_{i} \\ + (L_{i} - P_{i *} H_{- ii} {(P_{- ii})}^{- 1} L_{- i}) d \end{matrix}

(2)

where

H_{- ii} = {(I + C_{- i} P_{- ii})}^{- 1} C_{- i} P_{- ii}

(3)

Figure 1.

Control structures of multivariable systems: (a) decentralized control structure and (b) decoupling control structure.

$L_{- i}$ is the disturbance transfer function vector with L_i removed; $P_{i *}$ and $P_{* i}$ are the ith row and ith column of transfer function matrix P with p_ii removed, respectively; C_−i and P_−ii are controller matrix C and transfer function matrix P with the ith row and ith column removed, respectively; I is the unit matrix. Since other loops are assumed to be under perfect control condition, H_−ii satisfies H_−ii≈I. Substituting H_−ii = I into equation (2), the numerator of β_i can be derived as

{\frac{\partial u_{i}}{\partial d} |}_{y_{i}^{d} = 0} = - \frac{L_{i} - P_{i *} {(P_{- ii})}^{- 1} L_{- i}}{p_{ii} - P_{i *} {(P_{- ii})}^{- 1} P_{* i}}

(4)

On the other hand, when a multivariable system is well-decoupled, $y_{i}^{d}$ can be expressed as $y_{i}^{d} = p_{ii} u_{i} + L_{i} d$ . Hence,

{\frac{\partial u_{i}}{\partial d} |}_{y_{i}^{d} = 0} = - \frac{L_{i}}{p_{ii}}

(5)

Based on equations (4) and (5), β_i at steady state can be expressed as

β_{i} = \frac{L_{i} (0) p_{ii} (0) - P_{i *} (0) {(P_{- ii} (0))}^{- 1} L_{- i} (0) p_{ii} (0)}{L_{i} (0) p_{ii} (0) - P_{i *} (0) {(P_{- ii} (0))}^{- 1} P_{* i} (0) L_{i} (0)}

(6)

Next, β_i will be used to evaluate the couplings, which is essential to determine the constitution of the partial inverted decoupler.

Partial inverted decoupling control structure

The basic decoupling control structure is shown in Figure 1(b). In Figure 1(b), the decoupler Q(s), between the controller and the process, is the main part to be designed. The partial inverted decoupling control shares the same structure with Figure 1(b). The primary step is to design the partial inverted decoupler.

Design of the partial inverted decoupler

In decoupling control, the decoupler can be synthesized directly in terms of the desired decoupled process,

Q (s) = P^{- 1} (s) T (s)

(7)

where T(s) is the desired decoupled process. For partial decoupling, T(s) here should be determined based on RDG. For n × n processes, assuming that there are m (m ≤ n) loops with unfavorable couplings, the process transfer function matrix can be reformed as

P (s) = [\begin{matrix} P_{11} (s) & P_{12} (s) \\ P_{21} (s) & P_{22} (s) \end{matrix}]

(8)

where P₁₁(s), P₁₂(s), P₂₁(s), and P₂₂(s) are an m × m, m × (n−m), (n−m) × m and (n−m) × (n−m) matrix, respectively. Since the first m loops need decoupling, the decoupled process should have the following form,

T (s) = [\begin{matrix} T_{11} (s) & 0 \\ T_{21} (s) & T_{22} (s) \end{matrix}]

(9)

where T₁₁(s) = diag{t₁(s), t₂(s), …, t_m(s)}. To facilitate analysis, the decoupler matrix Q(s) is written as a block matrix, that is,

Q (s) = [\begin{matrix} Q_{11} (s) & Q_{12} (s) \\ Q_{21} (s) & Q_{22} (s) \end{matrix}]

(10)

The design of Q(s) is transformed to determine Q₁₁(s), Q₁₂(s), Q₂₁(s), and Q₂₂(s). In view of T(s) = P(s)Q(s), the following equations can be obtained:

{\begin{matrix} P_{11} (s) Q_{11} (s) + P_{12} (s) Q_{21} (s) = T_{11} (s) \\ P_{11} (s) Q_{12} (s) + P_{12} (s) Q_{22} (s) = 0 \\ P_{21} (s) Q_{11} (s) + P_{22} (s) Q_{21} (s) = T_{21} (s) \\ P_{21} (s) Q_{12} (s) + P_{22} (s) Q_{22} (s) = T_{22} (s) \end{matrix}

(11)

Considering that there are (n−m) loops using decentralized control, let Q₂₁(s) = 0 and Q₂₂(s) = I to simplify the design of the decoupler. As a result, equation (11) is rewritten as

{\begin{matrix} P_{11} (s) Q_{11} (s) = T_{11} (s) \\ P_{11} (s) Q_{12} (s) + P_{12} (s) = 0 \\ P_{21} (s) Q_{11} (s) = T_{21} (s) \\ P_{21} (s) Q_{12} (s) + P_{22} (s) = T_{22} (s) \end{matrix}

(12)

where the remaining undetermined parts are Q₁₁(s) and Q₁₂(s). According to equation (12), it is obvious that $Q_{11} (s) = P_{11}^{- 1} (s) T_{11} (s)$ , which is a problematic solve as it requires the inverse of a transfer function matrix. To deal with this problem, Q₁₁(s) is designed with the inverted decoupling which is implemented by the V canonical structure shown in Figure 2.

Figure 2.

Structure of inverted decoupling.

where Q₁₁(s) is realized by two parts, that is, V₁(s) and V₂(s). According to Figure 2, Q₁₁(s) can be expressed as

Q_{11} (s) = V_{1} (s) {[I - V_{2} (s) V_{1} (s)]}^{- 1}

(13)

Substituting $Q_{11} (s) = P_{11}^{- 1} (s) T_{11} (s)$ into equation (13), V₁(s) and V₂(s) should satisfy

V_{1}^{- 1} (s) - V_{2} (s) = T_{11}^{- 1} (s) P_{11} (s)

(14)

Since T₁₁(s) is a diagonal matrix, it is straightforward to compute its inverse. Similarly, V₁(s) should be determined as a diagonal matrix, that is, V₁(s) =diag{δ₁(s), δ₂(s), …, δ_m(s)}. For further simplification, V₂(s) is determined as

V_{2} (s) = [\begin{matrix} 0 & ρ_{12} (s) & \dots & ρ_{1 m} (s) \\ ρ_{12} (s) & 0 & \dots & ρ_{1 m} (s) \\ \dots & \dots & \dots & \dots \\ ρ_{m 1} (s) & ρ_{m 2} (s) & \dots & 0 \end{matrix}]

(15)

To satisfy equation (14), the elements of V₁(s) and V₂(s) should satisfy

{\begin{matrix} δ_{i} (s) = \frac{t_{i} (s)}{p_{i i} (s)} & i = 1, 2, \dots, m \\ ρ_{i j} (s) = - \frac{p_{i j} (s)}{t_{i} (s)} & i, j = 1, 2, \dots, m \end{matrix}

(16)

Equation (16) shows that the elements of V₁(s) and V₂(s) rely on the selection of the elements of T₁₁(s).

After determining Q₁₁(s), the design of Q₁₂(s) will be undertaken. In light of equation (12) and Q₁₁(s), Q₁₂(s) is derived as

Q_{12} (s) = V_{1} (s) {[I - V_{2} (s) V_{1} (s)]}^{- 1} V_{3} (s)

(17)

where $V_{3} (s) = - T_{11}^{- 1} (s) P_{12} (s)$ and its element η_ij(s) satisfies

η_{ij} (s) = - \frac{p_{i (m + j)} (s)}{t_{i} (s)} \begin{matrix} i = 1, 2, \dots, m, j = 1, 2, \dots, n - m \end{matrix}

(18)

Obviously, the elements of V₃(s) are also related to the elements of T₁₁(s). Therefore, how to determine T₁₁(s) needs a discussion.

The principle of determining T₁₁(s)

The above procedure reveals that the partial decoupler can be accomplished after T₁₁(s) is determined. However, considering the realizability of V₁(s), V₂(s), and V₃(s), there are some constraints for the elements of T₁₁(s). Concretely, the elements of V₁(s), V₂(s), and V₃(s) should be causal, proper and stable, which results in the constraints for the time delay, relative degree and right-half-plane (RHP) zeros of t_i(s).

Constraint for the time delay

Equations (16) and (18) show that each element of V₁(s), V₂(s), and V₃(s) is a ratio of two transfer functions. Hence, to make sure the elements of V₁(s), V₂(s), and V₃(s) are causal, the time delay of t_i(s) should satisfy

θ (p_{ii} (s)) \leq θ (t_{i} (s)) \leq min_{j \neq i} θ (p_{ij} (s)) \begin{matrix} i = 1, 2, \dots, m, j = 1, 2, \dots, n \end{matrix}

(19)

where θ(x) is the time delay of a transfer function x(s).

Constraint for the relative degree

Similarly, the elements of V₁(s), V₂(s), and V₃(s) will be proper if they have nonnegative relative degrees. Hence, the relative degree of t_i(s) should satisfy

\begin{matrix} φ (p_{ii} (s)) & \leq φ (t_{i} (s)) \leq min_{j \neq i} φ (p_{ij} (s)) \begin{matrix} i = 1, 2, \dots, m, j = 1, 2, \dots, n \end{matrix} \end{matrix}

(20)

where φ(x) is the relative degree of a transfer function x(s).

Constraint for the RHP zeros

Moreover, the elements of V₁(s), V₂(s), and V₃(s) should be stable. Therefore, the numerator and denominator in equations (16) and in (18) should have the same RHP zeros and the multiplicity of each zero should satisfy

ψ^{k} (p_{ii} (s)) \leq ψ^{k} (t_{i} (s)) \leq min_{j \neq i} ψ^{k} (p_{ij} (s)) \begin{matrix} i = 1, 2, \dots, m, j = 1, 2, \dots, n \end{matrix}

(21)

where ψ^k(x) is the multiplicity of the kth zero of a transfer function x(s).

From Inequalities (19) to (21), it can be seen that p_ii(s) should have the smallest values in terms of time delay, relative degree and multiplicity of RHP zero in the ith row. However, the original process matrix in equation (8) may not meet these requirements. If so, the process matrix P₁₁(s) can firstly be reformed by interchanging the rows to meet these requirements. Unfortunately, in some cases, reformation may make no difference for the above problem. An effective solution is to introduce a compensator expressed as equation (22) to modify the characteristic of the process.

N_{i} (s) = e^{- α_{i} s} \frac{1}{{(μ s + 1)}^{ε_{i}}} Π_{j = 1}^{k} {(\frac{- s + z_{j}}{s + z_{j}})}^{Γ_{i}^{j}}

(22)

N _i(s) is the compensator inserted in the ith loop and it is composed of three parts, which are used to compensate the time delay, relative degree and multiplicity of RHP zero, respectively. In equation (22), μ is a positive real constant; z_j denotes the jth RHP zero; α_i, ε_i and $Γ_{i}^{j}$ are undetermined compensator parameters. With the compensators, the compensated process transfer function matrix can be expressed as P_C(s) = P(s)N(s), where N(s) is the compensator matrix, that is, N(s) = diag{N₁(s), N₂(s), …, N_n(s)}. The compensator parameters α_i, ε_i and $Γ_{i}^{j}$ should make sure the compensated process satisfies the abovementioned requirements, that is, p_ii(s)N_i(s) has the smallest values in terms of time delay, relative degree and multiplicity of RHP zero in the ith row of the compensated process transfer function matrix. Consequently, the elements of V₁(s), V₂(s), and V₃(s) are designed based on the compensated process:

{\begin{matrix} δ_{i} (s) = \frac{t_{i} (s)}{p_{ii} (s) N_{i} (s)} \begin{matrix} i = 1, 2, \dots, m \end{matrix} \\ ρ_{ij} (s) = - \frac{p_{ij} (s) N_{j} (s)}{t_{i} (s)} \begin{matrix} i, j = 1, 2, \dots, m \end{matrix} \\ η_{ij} (s) = - \frac{p_{i (m + j)} (s) N_{m + j} (s)}{t_{i} (s)} \begin{matrix} i = 1, 2, \dots, m, j = 1, 2, \dots, n - m \end{matrix} \end{matrix}

(23)

Remark 1. The compensators are solely used for meeting the realizability requirements. Therefore, there is no need to select large compensator parameters. In fact, smaller compensator parameters will be preferable as it can result in a simpler compensated process, which is friendly to the design of the controller later. In a specific design problem, the calculation of the compensator parameters can be described as the following linear programming problem:

\begin{array}{l} \min X = \sum_{i = 1}^{n} x_{i} \\ s . t . {\begin{matrix} \sum_{k = 1}^{n} a_{i k} x_{k} \leq b_{i} & i = 1, 2, \dots, m \times (n - 1) \\ x_{i} \geq 0 & i = 1, 2, \dots, n \end{matrix} \end{array}

(24)

where x_i can be parameter α_i, ε_i or $Γ_{i}^{j}$ ; a_ik and b_i are coefficients of a specific problem. The above linear programming problem can be solved with the simplex method.

Remark 2. Equation (22) shows the general expression of the compensator. In most cases, the final compensator has a much simpler expression than equation (22). The reason is that, in process control, the elements of the process transfer function matrix are mostly described as first-order-plus-time-delay (FOPTD) models. Therefore, compensations for relative degree and RHP zeros are usually not needed. Generally, the compensators are introduced to compensate the time delays.

Finally, it is time to determine t_i(s) to finish the design of V₁(s) V₂(s), and V₃(s). Given the constraints of t_i(s), an intuitive and simple selection is t_i(s) = p_ii(s)N_i(s). This selection can result in the simplest form of V₁(s), that is, V₁(s) = I. Certainly, this is not the only way to select t_i(s). Sometimes, other simpler forms of t_i(s) are also available, but they will undoubtedly complicate V₁(s). With t_i(s) = p_ii(s)N_i(s), the elements of V₂(s) and V₃(s) can be derived accordingly as

{\begin{matrix} ρ_{ij} (s) = - \frac{p_{ij} (s) N_{j} (s)}{p_{ii} (s) N_{i} (s)} \begin{matrix} i, j = 1, 2, \dots, m \end{matrix} \\ η_{ij} (s) = - \frac{p_{i (m + j)} (s) N_{m + j} (s)}{p_{ii} (s) N_{i} (s)} \begin{matrix} i = 1, 2, \dots, m, j = 1, 2, \dots, n - m \end{matrix} \end{matrix}

(25)

In the end, the implementation structure of the partial inverted decoupling control is shown in Figure 3, where C_d(s) and C_m(s) are the decentralized controller matrix of decoupled loops and non-decoupled loops, respectively; L_d(s) and L_m(s) are the disturbance transfer function vector of decoupled loops and non-decoupled loops, respectively; Y_d(s) and Y_m(s) are the output vector of decoupled loops and non-decoupled loops, respectively; R_d(s) and R_m(s) are the reference vector of Y_d(s) and Y_m(s), respectively.

Figure 3.

Implementation structure of partial inverted decoupling.

In Figure 3, the decentralized controllers C_d(s) and C_m(s) are still the unknown parts, which will be studied in the next section.

Design of the decentralized controller

Decentralized controllers C_d(s) and C_m(s) are located in the decoupled loops and the non-decoupled loops, respectively, to improve the performance of disturbance rejection. The design of C_d(s) and C_m(s) will rely on the partially decoupled process. According to Figure 3, the disturbance response [Y_d(s), Y_m(s)]^T can be expressed as

\begin{matrix} [\begin{matrix} Y_{d} (s) \\ Y_{m} (s) \end{matrix}] = \\ [\begin{matrix} {(I + T_{11} (s) C_{d} (s))}^{- 1} L_{d} (s) \\ {(I + T_{22} (s) C_{m} (s))}^{- 1} (L_{m} (s) - T_{21} (s) C_{d} (s) {(I + T_{11} (s) C_{d} (s))}^{- 1} L_{d} (s)) \end{matrix}] d (s) \end{matrix}

(26)

From equation (26), it can be observed that Y_d(s) is only related to C_d(s) while Y_m(s) is determined by both C_d(s) and C_m(s). Thus, the decentralized controllers C_d(s) and C_m(s) should be designed in sequence. Firstly, C_d(s) will be designed by optimizing Y_d(s). And then the design of C_m(s) can be undertaken for optimizing Y_m(s) with the determined C_d(s).

Design of the controllers in the decoupled loops

As the couplings in the decoupled loops have been removed, the decoupled system can be regarded as a series of independent subsystems. The ith element of C_d(s), ${[C_{d} (s)]}_{ii}$ is the sub-controller of the ith subsystem. The design of ${[C_{d} (s)]}_{ii}$ is carried out by using the internal model control (IMC) structure shown in Figure 4, where ${[C_{imc} (s)]}_{ii}$ is the equivalent internal model controller. ${[C_{d} (s)]}_{ii}$ can be obtained based on equation (27) after ${[C_{imc} (s)]}_{ii}$ is determined.

{[C_{d} (s)]}_{ii} = \frac{{[C_{imc} (s)]}_{ii}}{1 - t_{i} (s) {[C_{imc} (s)]}_{ii}}

(27)

Figure 4.

Internal model control structure of the ith subsystem.

According to the IMC structure, the disturbance response $[Y_{d} (s)]_{i}$ is

{[Y_{d} (s)]}_{i} = (1 - t_{i} (s) {[C_{imc} (s)]}_{ii}) {[L_{d} (s)]}_{i} d (s)

(28)

where t_i(s) and $[L_{d} (s)]_{i}$ are assumed as FOPTD models:

\begin{matrix} t_{i} (s) = \frac{K_{ti}}{T_{ti} s + 1} e^{- τ_{i} s}, & {[L_{d} (s)]}_{i} = \frac{K_{Li}}{T_{Li} s + 1} e^{- ϑ_{i} s} \end{matrix}

(29)

Assuming d(s) = 1/s, the necessary condition for steady-state performance of disturbance rejection is

lim_{s \to 0} (1 - t_{i} (s) {[C_{imc} (s)]}_{ii}) {[L_{d} (s)]}_{i} = 0

(30)

Based on equation (30) and IMC principle, the potential internal model controller can be designed as

{[C_{imc} (s)]}_{ii} = t_{i -}^{- 1} (s) \cdot J_{i} (s)

(31)

where $t_{i -}^{- 1} (s)$ is the inverse of the minimum phase part of $t_{i} (s)$ and $J_{i} (s)$ is the undetermined low-pass filter which is used to make the controller realizable and improve the control performance.

From equation (28), obviously, if slow modes exist in the disturbance channel, the performance for disturbance rejection will be deteriorated. To speed up disturbance rejection, $J_{i} (s)$ is designed using equation (32), which has two undetermined parameters, that is, χ_i and λ_i.

J_{i} (s) = \frac{χ_{i} s + 1}{{(λ_{i} s + 1)}^{2}}

(32)

In order to eliminate the slow mode, the undesired poles of ${[L_{d} (s)]}_{i}$ should be cancelled. As a result, χ_i should satisfy

χ_{i} = \frac{2 λ_{i} \cdot T_{Li} + T_{Li} \cdot τ_{i} - λ_{i}^{2}}{T_{Li} + τ_{i}}

(33)

The proof of equation (33) is shown in the Appendix.

From equation (33), χ_i can be calculated when a concrete λ_i is given. Hence, the selection of λ_i will be focused on next. Obviously, in terms of speeding up disturbance rejection, a small value of λ_i is preferable. Meanwhile, a small value of λ_i can also provide a small gain for the disturbance output, which can suppress the peak value of the disturbance output (The gain of disturbance output can be approximately described as $K_{Li} (2 λ_{i} + τ_{i} - χ_{i}) = K_{Li} {(λ_{i} + τ_{i})}^{2} / (T_{Li} + τ_{i})$ ). Besides performance, robustness should also be considered due to the inevitable uncertainties. A widely-used robustness index is the maximum sensitivity (M_s):

M_{si} = ‖ 1 - t_{i} (s) {[C_{imc} (s)]}_{ii} ‖_{\infty}

(34)

A large value of M_s reveals bad robustness and a typical range of M_s is [1.4, 2.0]. For ease of analyzing the relationship between λ_i and M_si, define Λ_i = λ_i/T_Li, Ξ_i = T_Li/τ_i, and Ω_i = ω·τ_i, where ω is the frequency. Then, M_si in equation (34) can be rewritten as

\begin{matrix} M_{si} = max_{Ω_{i}} \\ | 1 - \frac{1 + \frac{j Ω_{i} (2 Λ_{i} Ξ_{i}^{2} + Ξ_{i} - Λ_{i}^{2} Ξ_{i}^{2})}{1 + Ξ_{i}}}{1 - {(Λ_{i} \cdot Ξ_{i} \cdot Ω_{i})}^{2} + 2 j (Λ_{i} \cdot Ξ_{i} \cdot Ω_{i})} e^{- j Ω_{i}} | \begin{matrix} Ω_{i} \in [0, + \infty) \end{matrix} \end{matrix}

(35)

where j is the imaginary unit. Based on equation (35), the relationship between system parameters and M_s is depicted in Figure 5.

Figure 5.

Relationship between system parameters and M_s: (a) relationship between Λ, Ξ and Ms and (b) relationship between Λ and Ξ for M_s = 2.0.

According to Figure 5(a), for a specific Ξ, M_s decreases while Λ increases. This reveals that a large λ (Λ) can improve the robustness. Therefore, with respect to the selection of λ, there is a tradeoff between robustness and performance. Meanwhile, for a specific Λ, M_s decreases while Ξ increases. In light of this, it can be inferred that if a system has a relatively large time delay, controller parameter λ ought to be increased to guarantee sufficient robustness. Figure 5(b) shows the relationship between Λ and Ξ for M_s = 2.0 and an approximate expression has been obtained as equation (36). In this way, the value of M_s will lie in the typical range, when λ ≥ 1.254T_L/(1.672 + T_L/τ).

(λ / T_{L}) = \frac{1.254}{1.672 + (T_{L} / τ)}

(36)

The design of internal model controller can be finished along with the determination of λ_i. Furthermore, based on equation (27), the equivalent feedback controller which is in the form of a PID controller can be fulfilled as

{[C_{d} (s)]}_{ii} = K_{Pi} (1 + \frac{1}{T_{Ii} s} + T_{Di} s) \frac{1}{T_{Fi} s + 1}

(37)

where

{\begin{matrix} K_{Pi} = \frac{(T_{ti} + χ_{i})}{K_{ti} (2 λ_{i} + τ_{i} - χ_{i})} \\ T_{Ii} = T_{ti} + χ_{i} \\ T_{Di} = \frac{T_{ti} \cdot χ_{i}}{T_{ti} + χ_{i}} \\ T_{Fi} = \frac{λ_{i}^{2} + τ_{i} \cdot χ_{i}}{(2 λ_{i} + τ_{i} - χ_{i})} \end{matrix}

(38)

Design of the controllers in the non-decoupled loops

After C_d(s) is obtained, the decentralized controllers in the non-decoupled loops, C_m(s) will be designed then. Similarly, the disturbance output Y_m (s) can be denoted as

Y_{m} (s) = {(I + T_{22} (s) C_{m} (s))}^{- 1} L_{e} (s) d (s)

(39)

where L_e(s) and T₂₂(s) are the equivalent disturbance channel model and the equivalent controlled process model, respectively. L_e(s) and T₂₂(s) have the following expressions:

{\begin{matrix} T_{22} (s) = (P_{22} (s) - P_{21} (s) P_{11}^{- 1} (s) P_{12} (s)) N_{M} (s) \\ L_{e} (s) = L_{m} (s) - P_{21} (s) P_{11}^{- 1} (s) ϒ (s) J (s) L_{d} (s) \end{matrix}

(40)

where

{\begin{matrix} N_{M} (s) = diag {N_{m + 1} (s), N_{m + 2} (s), \dots, N_{n} (s)} \\ ϒ (s) = diag {e^{- τ_{1} s}, e^{- τ_{2} s}, \dots, e^{- τ_{m} s}} \\ J (s) = diag {J_{1} (s), J_{2} (s), \dots, J_{m} (s)} \end{matrix}

(41)

According to equation (39), the design of C_m(s) relies on L_e(s) and T₂₂(s). For the independent design of the decentralized controllers, the idea of EOTF is used.¹³ The EOTFs of the process model and disturbance model are

{\begin{matrix} g_{i}^{eff} (s) = {[T_{22} (s)]}_{ii} - {[T_{22} (s)]}_{i *} \cdot {[T_{22} (s)]}_{- ii}^{- 1} \cdot {[T_{22} (s)]}_{* i} \\ L_{i}^{eff} (s) = {[L_{e} (s)]}_{i} - {[T_{22} (s)]}_{i *} \cdot {[T_{22} (s)]}_{- ii}^{- 1} \cdot {[L_{e} (s)]}_{- i} \end{matrix}

(42)

where ${[T_{22} (s)]}_{ii}$ is the ith row and ith column element of T₂₂(s); ${[T_{22} (s)]}_{i *}$ and ${[T_{22} (s)]}_{* i}$ are the ith row and ith column of T₂₂(s) with ${[T_{22} (s)]}_{ii}$ removed, respectively; ${[T_{22} (s)]}_{- ii}$ is T₂₂(s) removing its ith row and ith column; ${[L_{e} (s)]}_{i}$ is the ith element of L_e(s) and ${[L_{e} (s)]}_{- i}$ is L_e(s) removing ${[L_{e} (s)]}_{i}$ . Apparently, both $g_{i}^{eff} (s)$ and $L_{i}^{eff} (s)$ in equation (42) are too complex to be used for controller design. Therefore, model reduction technique is usually used to get low-order models. The approximate FOPDT model can be obtained based on the following method¹⁰:

g_{i}^{eff} (s) \approx \frac{{\bar{K}}_{gi} e^{- {\bar{τ}}_{gi} s}}{{\bar{T}}_{gi} s + 1}

(43)

where

{\begin{matrix} {\bar{K}}_{gi} = g_{i}^{eff} (0) \\ {\bar{T}}_{gi} = \sqrt{\frac{{\bar{K}}_{gi}^{2} - {| g_{i}^{eff} (j ω_{ci}) |}^{2}}{ω_{ci}^{2} {| g_{i}^{eff} (j ω_{ci}) |}^{2}}} \\ {\bar{τ}}_{gi} = \frac{π - \arctan ({\bar{T}}_{gi} ω_{ci})}{ω_{ci}} \end{matrix}

(44)

and ω_ci is the phase cross-over frequency of $g_{i}^{eff} (s)$ . Similarly, $L_{i}^{eff} (s)$ can be approximated as $L_{i}^{eff} (s) \approx {\bar{K}}_{Li} e^{- {\bar{τ}}_{Li} s} / ({\bar{T}}_{Li} s + 1)$ . Then, the ith sub-controller ${[C_{m} (s)]}_{ii}$ will be designed on the basis of the approximate models of $g_{i}^{eff} (s)$ and $L_{i}^{eff} (s)$ . The design procedure can consult the design of ${[C_{d} (s)]}_{ii}$ .

Simulation study

In this section, two examples are provided to demonstrate the effectiveness of the proposed method.

Example 1. Industrial scale polymerization (ISP) reactor

The process transfer function matrix of ISP reactor is

\begin{matrix} P (s) = [\begin{matrix} \frac{22.89 e^{- 0.2 s}}{4.572 s + 1} & \frac{- 11.64 e^{- 0.4 s}}{1.807 s + 1} \\ \frac{4.689 e^{- 0.2 s}}{2.174 s + 1} & \frac{5.8 e^{- 0.4 s}}{1.801 s + 1} \end{matrix}] \end{matrix}

(45)

The disturbance transfer function vector is assumed to be

L (s) = {[\begin{matrix} \frac{22.89 e^{- 0.2 s}}{4.572 s + 1} & \frac{5.8 e^{- 0.4 s}}{1.801 s + 1} \end{matrix}]}^{T}

(46)

First, the couplings should be evaluated based on RDG. According to equation (6), the values of RDG are computed as |β₁| = 1.069 and |β₂| = 0.1357, which indicates that the couplings in loop 1 need to be removed while the couplings in loop 2 need not. Therefore, the decoupled process is determined as

T (s) = [\begin{matrix} t_{11} (s) & 0 \\ t_{21} (s) & t_{22} (s) \end{matrix}]

(47)

Then, t₁₁(s) will be determined to design the components of the partial inverted decoupler, that is, V₁(s), V₂(s), and V₃(s). Fortunately, p₁₁(s) meets the requirements given in Section “The principle of determining T₁₁(s),” which means the extra compensation is unnecessary. In this case, N(s) = I and t₁₁(s) = p₁₁(s). Substituting t₁₁(s) into equations (16) and (18), the components are obtained as V₁(s) = 1, V₂(s) = 0 and V₃(s) = $- p_{12} (s) / p_{11} (s)$ . The next step is to design the decentralized controllers based on the decoupled process. For the decoupled loop, C_d(s) is designed based on t₁₁(s) and L_d(s), where t₁₁(s) = L_d(s) = 22.89e^−0.2s/(4.572s + 1). According to our proposed method, the controller parameter λ_d₁ should satisfy λ_d₁ > 0.2337, where λ_d₁ is determined as λ_d₁ = 0.7 to achieve better robustness. With λ_d₁, parameter χ_d₁ = 1.4303 is obtained by equation (33) and C_d(s) is finally determined as

C_{d} (s) = (1.545 + \frac{0.2574}{s} + 1.683 s) \frac{1}{4.572 s + 1}

(48)

For the non-decoupled loop, C_m(s) can be designed with the approximately equivalent models $g_{1}^{eff} (s) \approx 8.184 e^{- 0.4036 s} / (1.360 s + 1)$ and $L_{1}^{e f f} (s) \approx 1.111 e^{- 0.3776 s} / (0.3469 s + 1)$ . Similarly, the final C_m(s) is determined as

C_{m} (s) = (0.3157 + \frac{0.1853}{s} + 0.0867 s) \frac{1}{0.3469 s + 1}

(49)

with λ_m₁ = 0.3 for better robustness. With the decentralized controllers, the output response is shown in Figure 6(a). The proposed method is compared with several classical control strategies for multivariable processes, such as decentralized control,⁵² centralized control⁵³ and decoupling control.²² In Figure 6(a), it can be seen that the proposed method can achieve better performances for disturbance rejection. The peak of disturbance outputs of the proposed method are smaller. Meanwhile, the proposed method can bring the outputs back to their set-points in shorter time, especially in loop 2. Moreover, the integral absolute error (IAE) defined as

IAE = \sum_{i = 1}^{n} IA E_{i}

(50)

where

IA E_{i} = \int_{0}^{\infty} | y_{i} (t) | dt

(51)

is also used to evaluate the control performance. The IAE values of different control strategies are listed in Table 1 for comparison. It is noticed that the smallest IAE value is provided by the proposed strategy. In addition to the nominal control performance, it is also important to evaluate the robustness of the control system. Therefore, +20% parameter perturbations are introduced to the controlled process to evaluate the robust performance. Figure 6(b) depicts the disturbance outputs of the system with parameter perturbations. It can be seen that the proposed method can keep the better performances in terms of the output peak, response time and IAE value (The IAE values of the perturbed system are also listed in Table 1 and denoted by IAE_p). Then, the robust stability will be evaluated. For a multivariable system, both input uncertainty and output uncertainty should be considered. According to the M-Δ structure, in order to ensure robust stability, the multiplicative input uncertainty Δ_I(s) and the multiplicative output uncertainty Δ_O(s) should satisfy

\begin{matrix} {\begin{matrix} ‖ Δ_{I} (j ω) ‖ < 1 / \bar{σ} (M_{I} (j ω)) & \forall ω \\ ‖ Δ_{O} (j ω) ‖ < 1 / \bar{σ} (M_{O} (j ω)) & \forall ω \end{matrix} \end{matrix}

(52)

where

{\begin{matrix} M_{I} (j ω) = {(I + Q (j ω) C (j ω) P (j ω))}^{- 1} Q (j ω) C (j ω) P (j ω) \\ M_{O} (j ω) = {(I + P (j ω) Q (j ω) C (j ω))}^{- 1} P (j ω) Q (j ω) C (j ω) \end{matrix}

(53)

and C(s) = diag{C_d(s), C_m(s)}. The curves of $1 / \bar{σ} (M_{I} (j ω))$ and $1 / \bar{σ} (M_{O} (j ω))$ are plotted in Figure 6(c). The minimums of $1 / \bar{σ} (M_{I} (j ω))$ and $1 / \bar{σ} (M_{O} (j ω))$ at all frequencies denote the upper bounds of $‖ Δ_{I} (j ω) ‖$ and $‖ Δ_{O} (j ω) ‖$ . Therefore, generally, the minimums of $1 / \bar{σ} (M_{I} (j ω))$ and $1 / \bar{σ} (M_{O} (j ω))$ , γ_I and γ_O are used as robustness indexes.²⁶ Larger values of γ_I or γ_O reveal better robust stability. From Table 1, the proposed method provides the biggest γ_O and an acceptable γ_I, which means a good robust stability is achieved.

Example 2. Industrial distillation column

Figure 6.

Control performance for example 1: (a) output responses for example 1 on nominal condition, (b) output responses for example 1 on perturbed condition, and (c) curves of 1/||M(jω)|| for example 1.

Table 1.

Performance indexes of example 1.

Method	IAE₁	IAE₂	IAE	IAE_p	γ _I	γ _O
Proposed	3.8895	1.0560	4.9455	5.3752	0.4443	0.6669
Decoupling	3.4244	14.8361	18.2605	16.2709	0.5825	0.5973
Decentralized	9.4711	4.1524	13.6235	16.9606	0.2740	0.3775
Centralized	39.2103	4.6024	43.8127	42.1012	0.7989	0.5850

Here, a widely-studied 3 × 3 industrial distillation column is considered as

\begin{matrix} P (s) = [\begin{matrix} \frac{1.986 e^{- 0.71 s}}{66.7 s + 1} & \frac{- 5.24 e^{- 60 s}}{400 s + 1} & \frac{- 5.984 e^{- 2.24 s}}{14.29 s + 1} \\ \frac{- 0.0204 e^{- 0.59 s}}{{(7.14 s + 1)}^{2}} & \frac{0.33 e^{- 0.68 s}}{{(2.38 s + 1)}^{2}} & \frac{- 2.38 e^{- 0.42 s}}{{(1.43 s + 1)}^{2}} \\ \frac{- 0.374 e^{- 7.75 s}}{22.22 s + 1} & \frac{11.3 e^{- 3.79 s}}{{(21.74 s + 1)}^{2}} & \frac{9.811 e^{- 1.59 s}}{11.36 s + 1} \end{matrix}] \end{matrix}

(54)

with an assumed disturbance transfer function vector

L (s) = {[\begin{matrix} \frac{1.986 e^{- 0.71 s}}{66.7 s + 1} & \frac{0.33 e^{- 0.68 s}}{{(2.38 s + 1)}^{2}} & \frac{9.811 e^{- 1.59 s}}{11.36 s + 1} \end{matrix}]}^{T}

(55)

First, RDGs are calculated as |β₁| = 3.5827 and |β₂| = 1.0120 and |β₃| = 0.0290 to evaluate couplings. Consequently, the couplings of the first two loops need to be removed while that in the last loop need not. As a result, the decoupled process is determined as

T (s) = [\begin{matrix} t_{11} (s) & 0 & 0 \\ 0 & t_{22} (s) & 0 \\ t_{31} (s) & t_{32} (s) & t_{33} (s) \end{matrix}]

(56)

In this case, a compensator designed as N(s) =diag{e^−0.09s, 1, e^−0.26s} is needed as the requirement for the time delay cannot be met if t₁₁(s) = p₁₁(s) and t₂₂(s) = p₂₂(s) are selected directly. Based on the compensated process, V₁(s), V₂(s), and V₃(s) are achieved as V₁(s) = I,

V_{2} (s) = [\begin{matrix} 0 & \frac{2.6385 (66.7 s + 1) e^{- 59.2 s}}{400 s + 1} \\ \frac{0.06182 {(2.38 s + 1)}^{2}}{{(7.14 s + 1)}^{2}} & 0 \end{matrix}]

(57)

and V₃(s) = $[\begin{matrix} 3.0131 e^{- 1.7 s} (66.7 s + 1) / (14.29 s + 1), & 7.2121 {(2.38 s + 1)}^{2} / {(1.43 s + 1)}^{2} \end{matrix}]^{T}$ , respectively. Next, the decentralized controllers for the partially decoupled process are designed. According to the proposed method, the decentralized controllers for the decoupled loops are

\begin{matrix} C_{d} (s) = diag {(110.8718 + \frac{1.4752}{s} + 832.2807 s) \frac{1}{66.7 s + 1} \\ (2.3626 + \frac{0.3493}{s} + 3.9828 s) \frac{1}{3.57 s + 1}} \end{matrix}

(58)

Then, the controller of the non-decoupled loop, C_m(s) is designed based on the approximate equivalent models, $g_{1}^{eff} (s) \approx 99.856 e^{- 2.530 s} / (123.746 s + 1)$ and $L_{1}^{eff} (s) \approx - 2.90 (- 124.855 s + 1) e^{- 5.597 s} / (410.506 s^{2} + 40.045 s + 1)$ . Here, a second-order-plus-time-delay (SOPTD) model is used to approximate $L_{1}^{eff} (s)$ , because the dynamics of $L_{1}^{eff} (s)$ are too complex. With the similar thought, C_m(s) is finally determined as

\begin{matrix} C_{m} (s) = (0.07578 + \frac{0.002637}{s} + 0.5824 s) \\ \frac{123.746 s + 1}{410.506 s^{2} + 40.045 s + 1} \end{matrix}

(59)

The output responses of the nominal process are shown in Figure 7(a) compared with Liu’s⁵⁴ analytical decoupling control strategy, ETF-based partial decoupling control strategy⁴⁹ and controller synthesis based on RDGA.⁴⁷ The corresponding IAE values are listed in Table 2. Obviously, the proposed method provide better output performances than other methods. Similarly, to evaluate robust performance, −20% parameter perturbations are introduced to the controlled process and the corresponding output responses are depicted in Figure 7(b). The proposed method can still achieve the best IAE value for the perturbed process. Besides, the curves of $1 / \bar{σ} (M_{I} (j ω))$ and $1 / \bar{σ} (M_{O} (j ω))$ and the robustness indexes, γ_I and γ_O are also shown in Figure 7(c) and Table 2, respectively. It can be seen that the proposed method has better robustness than ETF-based partial decoupling control strategy. Liu’s method achieve bigger γ_I and γ_O, but it suffers from complex calculation and controller form. The controllers based on RDGA achieve the biggest γ_I, but it results in worse output responses.

Figure 7.

Control performance for example 2: (a) output responses for example 2 on nominal condition, (b) output responses for example 2 on perturbed condition, and (c) curves of 1/||M(jω)|| for example 2.

Table 2.

Performance indexes of example 2.

Method	IAE₁	IAE₂	IAE₃	IAE	IAE_P	γ _I	γ _O
Proposed	0.737	2.889	72.066	75.692	102.555	0.05592	0.2746
Liu	12.292	3.244	77.625	93.161	119.853	0.1559	0.7212
ETF	85.242	1.174	1172.6	1259.016	1137.189	0.01044	0.1495
RDGA	496.268	3.418	544.710	1044.396	1314.852	0.6539	0.2628

Conclusion

In this work, a novel partial inverted decoupling control strategy to improve the disturbance rejection of multivariable process with time delays is presented. Two main issues: the design of the partial decoupler and the decentralized controllers, are focused on. Motivated by the thought of using favorable couplings to reject disturbances, the couplings are firstly evaluated with RDGs. Then the inverted decoupling usually used in full decoupling control is extended to partial decoupling control. The elements of the partial inverted decoupler are designed in detail and the implementation structure is presented. The advantage of the proposed partial inverted decoupler is that the decoupler is implemented without approximation, which can remove the unfavorable couplings thoroughly. Based on the partial decoupled processes, the decentralized controllers for both the decoupled and non-decoupled loops are designed. An analytical design method is proposed with due consideration of the performance of disturbance rejection and robustness. Moreover, to improve the control performance of the non-decoupled loop, the equivalent process model and the equivalent disturbance channel model are obtained for controller design. Two simulation examples are provided and the results show that the proposed method can achieve better nominal control performance and robustness compared with several existing multivariable control strategies.

Footnotes

Appendix

Substituting equations (29), (31), and (32) into equation (28), the transfer function from d(s) to $[Y_{d} (s)]_{i}$ is

(A.1)

\frac{{[Y_{d} (s)]}_{i}}{d (s)} = [1 - \frac{(χ_{i} s + 1)}{{(λ_{i} s + 1)}^{2}} e^{- τ_{i} s}] \frac{K_{Li} e^{- ϑ_{i} s}}{T_{Li} s + 1}

Using the first-order Taylor approximation, $e^{- τ_{i} s} \approx 1 - τ_{i} s$ , the above equation can be rewritten as

(A.2)

\begin{matrix} [1 - \frac{(χ_{i} s + 1)}{{(λ_{i} s + 1)}^{2}} e^{- τ_{i} s}] \frac{K_{Li} e^{- ϑ_{i} s}}{T_{Li} s + 1} \\ \approx \frac{[(2 λ_{i} + τ_{i} - χ_{i}) s (\frac{λ_{i}^{2} + χ_{i} \cdot τ_{i}}{2 λ_{i} + τ_{i} - χ_{i}} s + 1)]}{{(λ_{i} s + 1)}^{2}} \frac{K_{Li} e^{- ϑ_{i} s}}{T_{Li} s + 1} \end{matrix}

To cancel the undesired pole of $[L_{d} (s)]_{i}$ , χ_i should satisfy

(A.3)

\frac{λ_{i}^{2} + χ_{i} \cdot τ_{i}}{2 λ_{i} + τ_{i} - χ_{i}} = T_{Li}

Equation (33) can be obtained based on equation (A.3).

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work is supported by the National Natural Science Foundation of China [62003152] and Xuzhou Science and Technology Project of China [KC23008].

ORCID iDs

Xinghan Du

Wu Cai

Data availability statement

Data sharing not applicable to this article as no datasets were generated or analyzed during the current study.

References

Zhang

Quantitative process control theory. Boca Raton: CRC Press, 2012.

Euzébio

TAM

Silva

MTD

Yamashita

AS.

Decentralized PID controller tuning based on nonlinear optimization to minimize the disturbance effects in coupled loops. IEEE Access 2021; 9: 156857–156867.

Mappas

Vassiliadis

Dorneanu

, et al. Automated control loop selection via multistage optimal control formulation and nonlinear programming. Chem Eng Res Des 2023; 195: 76–95.

Bristol

EH.

On a new measure of interaction for multivariable process control. IEEE Trans Autom Control 1966; 11(1): 133–134.

Cai

, et al. RNGA based control system configuration for multivariable processes. J Process Control 2009; 19(6): 1036–1042.

Liao

Cai

, et al. Interaction analysis and loop pairing for MIMO processes described by T–S fuzzy models. Fuzzy Sets Syst 2012; 207: 64–76.

Luyben

WL.

Simple Method for Tuning SISO controllers in multivariable systems. Ind Eng Chem Process Des Dev 1986; 25(3): 654–660.

Wang

Zou

Lee

, et al. Auto-tuning of multivariable PID controllers from decentralized relay feedback. Automatica 1997; 33(3): 319–330.

Lee

Cho

Edgar

TF.

Multiloop PI controller tuning for interacting multivariable processes. Comput Chem Eng 1998; 22(11): 1711–1723.

10.

Mahapatro

Subudhi

Ghosh

Design and experimental realization of a robust decentralized PI controller for a coupled tank system. ISA Trans 2019; 89: 158–168.

11.

Mahapatro

Subudhi

A Robust Decentralized PID controller based on complementary sensitivity function for a multivariable system. IEEE Trans Circuits Syst II Express Briefs 2020; 67(10): 2024–2028.

12.

Jeng

Lee

MW.

Multi-loop PID controllers design with reduced loop interactions based on a frequency-domain direct synthesis method. J Franklin Inst 2023; 360: 2476–2506.

13.

Huang

Jeng

Chiang

, et al. A direct method for multi-loop PI/PID controller design. J Process Control 2003; 13(8): 769–786.

14.

Zhang

Xiao

Energy-efficient decentralized control method with enhanced robustness for multi-evaporator air conditioning systems. Appl Energy 2020; 279: 115732.

15.

Chu

MH.

IMC-PID tuning method based on maximum sensitivity for uncertain multivariable systems. Int J Robust Nonlinear Control 2023; 33(13): 7395–7414.

16.

Liu

Tian

Xue

, et al. A review of industrial MIMO decoupling control. Int J Control Autom Syst 2019; 17(5): 1246–1254.

17.

Gagnon

Pomerleau

Desbiens

Simplified, ideal or inverted decoupling?

ISA Trans 1998; 37(4): 265–276.

18.

Garrido

Vázquez

Morilla

Centralized multivariable control by simplified decoupling. J Process Control 2012; 22(6): 1044–1062.

19.

Garrido

Vázquez

Morilla

An extended approach of inverted decoupling. J Process Control 2011; 21(1): 55–68.

20.

Aguiar

APVDA

Júnior

Barros

. Evaluation and redesign of simplified decouplers for TITO processes using relay experiments. J Control Autom Electr Syst 2021; 32: 830–841.

21.

Aguiar

APVDA

Júnior

Barros

. Evaluation and redesign of the inverted decoupler: open and closed-loop approaches. Int J Control Autom Syst 2020; 18(6): 1435–1444.

22.

Jin

Liu

Decoupling proportional–Integral–Derivative controller design for multivariable processes with time delays. Ind Eng Chem Res 2014; 53(2): 765–777.

23.

Xiong

Cai

MJ.

Equivalent transfer function method for PI/PID controller design of MIMO processes. J Process Control 2007; 17(8): 665–673.

24.

Naik

Kumar

DVA

Sujatha

Independent controller design for MIMO processes based on extended simplified decoupler and equivalent transfer function. Ain Shams Eng J 2020; 11: 343–350.

25.

Khandelwal

Detroja

KP.

The optimal detuning approach based centralized control design for MIMO processes. J Process Control 2020; 96: 23–36.

26.

Jin

Jiang

BY.

Novel centralized IMC-PID controller design for multivariable processes with multiple time delays. Ind Eng Chem Res 2017; 56(15): 4431–4445.

27.

Ghosh

Pan

Centralized PI controller design method for MIMO processes based on frequency response approximation. ISA Trans 2021; 110: 117–128.

28.

Garrido

Ruz

Morilla

, et al. Iterative design of centralized PID controllers based on equivalent loop transfer functions and linear programming. IEEE Access 2021; 10: 1440–1450.

29.

Gündeş

AN.

Controller design for diagonal decoupling and integral action. IEEE Trans Autom Control 2021; 66(3): 1361–1366.

30.

Emre

Roberto

Kouhei

Disturbance observer-based robust control and its applications: 35th anniversary overview. IEEE Trans Ind Electron 2020; 67(3): 2042–2053.

31.

Sun

You

FQ.

Machine learning and data-driven techniques for the control of smart power generation systems: an uncertainty handling perspective. Engineering 2021; 7: 1239–1247.

32.

Sun

Lee

KY.

Optimal disturbance rejection for PI controller with constraints on relative delay margin. ISA Trans 2016; 63: 103–111.

33.

Zhang

, et al. Design of modified 2-degree-of-freedom proportional–integral–derivative controller for unstable processes. Meas Control 2020; 53: 1465–1471.

34.

Damour

Benne

Lebreton

, et al. Real-time implementation of a neural model-based self-tuning PID strategy for oxygen stoichiometry control in PEM fuel cell. Int J Hydrogen Energy 2014; 39(24): 12819–12825.

35.

Asgharnia

Shahnazi

Jamali

Performance and robustness of optimal fractional fuzzy PID controllers for pitch control of a wind turbine using chaotic optimization algorithms. ISA Trans 2018; 79: 27–44.

36.

Wang

Luo

A novel HTD-CS based PID controller tuning method for time delay continuous systems with multi objective and multi-constraint optimization. Chem Eng Res Des 2016; 115: 98–106.

37.

Jeng

Jeon

Kang

, et al. Iterative feedback tuning of cascade control of two-inertia system. IEEE Control Syst Lett 2021; 5(3): 785–790.

38.

William

Marcello

Virtual reference feedback tuning for linear discrete-time systems with robust stability guarantees based on set membership. Automatica 2023; 157: 111228.

39.

Liu

Jin

, et al. Design of DOB-based control system in the presence of uncertain delays for low-order processes. IEEE Trans Control Syst Technol 2020; 28(2): 558–565.

40.

Sun

Zhu

, et al. Robust tuning of active disturbance rejection controller for time-delay systems with application to a factual electrostatic precipitator. IEEE Trans Control Syst Technol 2022; 30(5): 2204–2211.

41.

Cao

She

, et al. Disturbance rejection and robustness of improved equivalent-input-disturbance-based system. IEEE Trans on Cybern 2022; 52(8): 8537–8546.

42.

Sun

Lee

KY.

Enhanced decentralized PI control for fluidized bed combustor via advanced disturbance observer. Control Eng Pract 2015; 42: 128–139.

43.

Sun

Xue

, et al. Quantitative tuning of active disturbance rejection controller for FOPTD model with application to power plant control. IEEE Trans Ind Electron 2022; 69(1): 805–815.

44.

Sun

Zhang

, et al. Tuning of active disturbance rejection control with application to power plant furnace regulation. Control Eng Pract 2019; 92: 104122.

45.

Stanley

Maria

Mcavoy

TJ.

Short cut operability analysis. 1. The relative disturbance gain. Ind Eng Chem Process Des Dev 1985; 24(4): 1181–1188.

46.

Lin

Jeng

Huang

HP.

Multivariable control with generalized decoupling for disturbance rejection. Ind Eng Chem Res 2009; 48(20): 9175–9185.

47.

Chang

CC.

Relative disturbance gain array. AIChE J 1992; 38(4): 521–534.

48.

Gilbert

Yousef

Natarajan

, et al. Tuning of PI controllers with one-way decoupling in 2×2 MIMO systems based on finite frequency response data. J Process Control 2003; 13: 553–567.

49.

Shen

Ning

, et al. Partial decoupling control for multivariable processes. Ind Eng Chem Res 2011; 50(12): 7380–7387.

50.

Jin

Cai

Disturbance observer-based control with partially coupled inverse model for stable TITO processes with time delays. ISA Trans 2021; 110: 148–159.

51.

Cai

Shen

Shang

, et al. Partially coupled inverse model design for multivariable disturbance observer-based control based on generalized relative input disturbance gain. IEEE Access 2022; 10: 90388–90400.

52.

Jin

Liu

Multi-loop PI/PID controllers design for disturbance rejection based on non-parametric effective model and non-convex optimization. IET Control Theory Appl 2014; 8(15): 1499–1512.

53.

Besta

Chidambaram

Tuning of multivariable PI controllers by BLT method for TITO systems. Chem Eng Commun 2016; 207: 527–538.

54.

Liu

Zhang

Gao

Analytical decoupling control strategy using a unity feedback control structure for MIMO processes with time delays. J Process Control 2007; 17: 173–186.

Disturbance rejection for multivariable processes with time delays based on partial inverted decoupling control

Abstract

Keywords

Introduction

Relative disturbance gain (RDG)

Partial inverted decoupling control structure

Design of the partial inverted decoupler

The principle of determining T11(s)

Constraint for the time delay

Constraint for the relative degree

Constraint for the RHP zeros

Design of the decentralized controller

Design of the controllers in the decoupled loops

Design of the controllers in the non-decoupled loops

Simulation study

Conclusion

Footnotes

Appendix

Declaration of conflicting interests

Funding

ORCID iDs

Data availability statement

References

The principle of determining T₁₁(s)