# A New Hybrid BFOA-PSO Optimization Technique for Decoupling and Robust Control of Two-Coupled Distillation Column Process

Authors: Noha Abdelkarim; Amr E. Mohamed; Ahmed M. El-Garhy; Hassen T. Dorrah
Journal: Computational Intelligence and Neuroscience
 (2016)

---

## Abstract

The two-coupled distillation column process is a physically complicated system in many aspects. Specifically, the nested interrelationship between system inputs and outputs constitutes one of the significant challenges in system control design. Mostly, such a process is to be decoupled into several input/output pairings (loops), so that a single controller can be assigned for each loop. In the frame of this research, the Brain Emotional Learning Based Intelligent Controller (BELBIC) forms the control structure for each decoupled loop. The paper’s main objective is to develop a parameterization technique for decoupling and control schemes, which ensures robust control behavior. In this regard, the novel optimization technique Bacterial Swarm Optimization (BSO) is utilized for the minimization of summation of the integral time-weighted squared errors (ITSEs) for all control loops. This optimization technique constitutes a hybrid between two techniques, which are the Particle Swarm and Bacterial Foraging algorithms. According to the simulation results, this hybridized technique ensures low mathematical burdens and high decoupling and control accuracy. Moreover, the behavior analysis of the proposed BELBIC shows a remarkable improvement in the time domain behavior and robustness over the conventional PID controller.


---

## Body

## 1. Introduction

Controller design for multi-input multiple-output (MIMO) systems is of significant importance, as they constitute the majority of physical systems. However, the design of such controllers is confronted with the challenge to overcome the influence of the nested interrelationship between system inputs and outputs. Based on the level of interrelationship between system inputs and outputs, the control structure is to be selected. The systems, which are characterized with moderate interaction between control loops, can be controlled with the decentralized control structure so that a single controller is assigned for each control loop [1–5]. On the other hand, a centralized controller has to be designed for systems with higher level of interrelationship between control loops [6–9], so that a central matrix of controllers is allocated for the whole MIMO system. In spite of this, the centralized controllers are not widely utilized due to their high complexity. Alternatively, the issue of high interrelationship can be alleviated by decoupling the MIMO system into several relatively independent single-input single-output (SISO) control loops [10–15]. Accordingly, the multivariable process can be controlled based on independent loop structure. In this paper, the decoupled control structure is utilized for controlling the presented MIMO system. For the decoupler to be designed, the adequate input-output pairings are primarily determined through the evaluation of the system relative gains [16–19]. Accordingly, the decoupling network has to be designed, so that the interaction between control loops is minimized [20, 21]. For control purposes, conventional control schemes as PID can be utilized [22–24]. It is to mention that the nonlinearity and the model imprecision that characterize the majority of physical systems reduce the robustness and accuracy of such controllers.Control design approaches based on intelligent algorithms as fuzzy logic, neural network, and genetic algorithms are of increasing spread due to their proven ability to overcome system model uncertainty [25–27]. The computational model of emotional learning in mammalian brain, which is introduced in [28, 29], inspired a new learning algorithm. Afterwards, this learning technique is deployed in the system control design presenting the novel Brain Emotional Learning Based Intelligent Controller (BELBIC) [30]. Thus, the Brain Emotional Learning Based Intelligent Controller (BELBIC) has shown robustness in the control of the nonlinear and uncertain systems such as Van Der Pol oscillator, Duffing forced oscillator and automatic self-balancing scale [31], washing machine [32], microheat exchanger [33], switched reluctance motor [34], unmanned aerial vehicle [35], path tracking of a vehicle [36], two-coupled distillation column process [37], multiple-area power systems [38], and continuous stirred tank reactor [39, 40].The parameterization of the decoupling compensation network and the BELBIC are still open research topics. Currently, the detailed analytical methods used in purpose of the parameterization of the decoupling compensation network cost excessive mathematical burdens. Regarding the parameterization of BELBIC, the trial and error is widely used for controller parameters estimation. Thus, the utilization of various optimization techniques based on artificial intelligence introduces an efficient alternative in both cases [11, 37].The utilization of the biologically inspired algorithms as Ant Colony Optimization (ACO) [41], Genetic Algorithm (GA) [42], and Particle Swarm Optimization (PSO) [43, 44] in optimization problems constitutes currently a promising solution. Recently, the new evolutionary computation technique depending on the behavior of foraging ofE. coli bacteria, which is named as Bacterial Foraging Optimization Algorithm (BFOA), is proposed by [45]. This technique has been successfully deployed in many applications as power systems [46–49], stock market prediction [50, 51], and design of PI/PID controllers [52, 53]. The key drawback of the BFOA is the delay in reaching the global solution because it is based on random searching directions. For this delay to be overcome, the BFOA is integrated with the Particle Swarm Optimization (PSO) technique, introducing the Bacterial Swarm Optimization (BSO) algorithm [54]. One of the major characteristics of PSO technique, which is inherited to the BSO algorithm, is the idea of velocity updating. Accordingly, the process of searching for the global solution in the BSO technique depends on the individual and global best positions concurrently. As discussed in [55–58], The Bacterial Swarm Optimization (BSO) algorithm grants better performance in determining the optimum solution compared with PSO and BFOA algorithms.This research paper aims to investigate the feasibility of applying the Bacterial Swarm Optimization (BSO) algorithm in the field of MIMO control system. As an example for MIMO system the two-coupled distillation column is to be studied. Primarily, the mathematical model of the system is decoupled into several independent loops. The parameters of the decoupling compensation network are determined based on the summation minimization of the integral time-weighted squared outputs (ITSOs) of unpaired outputs regarding a particular input. In this regard, the BSO technique is utilized. On the other hand, an optimal BELBIC is designed for each decoupled control loop through the summation minimization of the integral time-weighted squared errors (ITSEs) of all loops using BSO as well. For the comparison purpose the PID controller is also implemented for the same application, so that the strength of each of considered control structures can be studied.


## 2. Two-Coupled Distillation Column Process

As formerly stated, the two-coupled distillation column is handled in this paper as an example of MIMO system. In this section, the physical system is discussed. Primarily, the function of the system and its components are presented. Afterwards, the system mathematical modeling and decoupling are handled.


### 2.1. System Description

Distillation units are mainly utilized for the separation of fluid mixture components. The distillation column major components can be listed as follows:(i)
A vertical shell in which the separation of fluid substances is accomplished(ii)
A cascade of trays for improving component separation(iii)
A reboiler for maintaining the required heat energy for the distillation process(iv)
A condenser for liquefying the vapor leaving the column(v)
A reflux drum in which portion of the condensed liquid is recycled back to the vertical shell.The physical process considered in this research comprises two-coupled distillation columns as shown in Figure1. Thereby, ternary petrochemical mixtures can be separated. While the main glass column is composed of 40 bubble cap trays excluding those of the boiler and condenser, the side glass column contains 10 bubble cap trays. At the 22nd stage, the system intake port is located, at which the fluid mixture is to be fed. The three separated fluid components are to be extracted from the process, so that the heaviest and the lightest components can be provided from the bottom and the top of the main column, respectively, and the intermediate component is derived from the top of the side column.Figure 1
The two-coupled distillation columns process.


### 2.2. Mathematical Modeling of the Two-Coupled Distillation Column

Primarily, the manipulated variables are to be determined. Thus, the selected manipulated variables are the heat input to the reboiler (QE), the vapor flow rate in the vapor transfer line (SAB), the reflux ratio in the main column (RL1), and the reflux ratio in the second column (RL2) [59]. Regarding the system technical constraints, the input heat energy rate QE and the vapor flow rate SAB have to not exceed the values 8.2 KW and 3.95 m3/h, respectively, in order to avoid flooding of the smaller side column [59]. Secondly, the controlled values are selected to be the fluid temperature at four different stages. The selection of these stages is based on the following requirements:(i)
The temperatures should be sensitive to the major disturbances as changes in feed rate and feed concentrations.(ii)
The temperatures should exhibit a fairly linear behavior.(iii)
They should be a good indication of product quality.Accordingly, the temperature at the 11th, 30th, 34th, and 48th trays are chosen to form the controlled variables [59]. The mathematical model of realistic two-coupled distillation column, mentioned in [59], is to be used in this research. The transfer function matrix of the considered MIMO system is stated in (1).(1)Gs=2.61.69s+1-6.0983.5s+1-4.990.2s+14.5s+10.06s+10.0713.5s+17.321.05s+110.4s+10.14s+1-1.450.4s+1-1.570.23s+11.34s+10.2s+1-0.141.92s+14.60.53s+12.78s+10.09s+1-2.370.23s+12s+10.3s+1-2.71.75s+1-0.360.02s+12.47s+10.04s+12.110.92s+1-2.110.06s+12.38s+10.05s+1-1.752.16s+1-0.31.89s+14.35s+10.16s+1.As shown in the transfer function matrix, the physical system is characterized by the high interrelationship between all system inputs and outputs. As discussed in the introduction, the control structure adopted in this paper is based on the mathematical model decoupling, so that each output is independently controlled by a single input, forming four independent SISO control loops.Primarily, the noninteracting design is preceded by a relative gain analysis to determine the most suitable input-output pairings. In this regard, the relative gain array (RGA) method is used. Afterwards, a decoupling compensation network is designed for the reduction of residual interactions [11].


#### 2.2.1. The Relative Gain Array

The RGA matrix forN×N system represented in (2) is used for the valuation of the input influences on each system output [16–19], so that output “i” is to be paired with input “j” for which γij is a positive value and as close to unity as possible.(2)RGAN×N=Gs=0·∗Gs=0-T=γ11γ12⋯γ1Nγ21γ22⋯γ2N⋮⋮⋮γN1γN2⋯γNN,where G(s) is the transfer function matrix of the MIMO process and “·∗” operator implies element by element multiplication.


#### 2.2.2. Decoupling Compensation Network Design Procedure

The steady state decoupling compensation matrix illustrated in (3) is to be integrated into the given MIMO system in the way demonstrated in (4), so that interrelationships between each input and unpaired outputs is minimized [20, 21]. In this paper, the simplest form of the decoupling compensation matrix introduced in the previous study [11], which has unity diagonal elements, is replaced with the general form presented in (3). This modification introduces a remarkable improvement in the system response as discussed in Section 5. Figure 2 illustrates the decoupling compensation network scheme integrated with the required controllers.(3)Λss=λ11λ12λ13⋯λ1N-1λ1Nλ21λ22λ23⋯λ2N-1λ2N⋮⋮⋮⋮⋮λN1λN2λN3⋯λNN-1λNN,(4)Ys=GsΛssWs,where Ys is the output vector, Ws is the input vector, Gs is the transfer function matrix of the MIMO process, and Λss is the steady state decoupling compensation matrix of MIMO process.Figure 2
The decoupled MIMO system including controllers.


#### 2.2.3. Fitness Function Design

As discussed in the former two sections, the determination of the adequate system input/output pairing as well as the elements’ values of the decoupling compensation matrix define the effectiveness of interrelationship minimization between system inputs and outputs. The selection of the proper input/output pairing is already tackled by the relative gain array concept. For the elements’ values of the decoupling compensation matrix to be estimated using optimization techniques, the fitness functions have to be designed that minimize the interrelationship between system inputs and their unpaired outputs.The four commonly used performance criteria for fitness function design are the integral absolute error (IAE), integral squared error (ISE), integral time-weighted squared error (ITSE), and Integral time-weighted absolute error (ITAE). It’s to mention that in the context of decoupling system design, the performance criteria are concerned with the minimization of the system outputs towards the unpaired inputs rather than control system error. Accordingly, the corresponding criteria are integral absolute output (IAO), integral squared output (ISO), integral time-weighted squared output (ITSO), and Integral time-weighted absolute output (ITAO). Although the criteria ISO and IAO grant less overshoot in the system dynamics, none of these criteria is adopted in this research due to the long settling time [56]. This drawback could be overcome by utilizing the ITAO. However, this criterion is not used in this research for the difficulty of its analytical tracking [60]. Thus, the ITSO performance criterion is employed in this work for ensuring the minimization of settling time without being confronted with unnecessary analytical complications.In the former research [11], the decoupling compensation matrix elements are estimated by applying the PSO technique, so that the integral squared outputs (ISOs) of unpaired outputs with respect to a specific input are minimized. On the other hand, the proposed technique in this paper is based on the minimization of the integral time-weighted squared outputs (ITSOs) by utilizing the BSO technique.The integral time-weighted squared outputs (ITSOs) for each input are calculated such thatifW1(s)=1/s and W2,W3,…,WN are zeroes, then(5)ITSO11=∫0∞tY12·dt,ITSO21=∫0∞tY22·dt,⋮ITSON1=∫0∞tYN2·dtand ifW2(s)=1/s and W1,W3,…,WN are zeroes, then(6)ITSO12=∫0∞tY12·dt,ITSO22=∫0∞tY22·dt,⋮ITSON2=∫0∞tYN2·dt.Consecutively, ifWN(s)=1/s and W1,W2,…,WN-1 are zeroes, then(7)ITSO1N=∫0∞tY12·dt,ITSO2N=∫0∞tY22·dt,⋮ITSONN=∫0∞tYN2·dt.Thus the fitness function for specific inputWj can be described as follows:(8)Fitnessj=∑i=1NITSOij,j=1,2,3,…,N,i≠q,where i is the specific output subscript, j is the specific input subscript, and q is the subscript of the output that has been paired with input j.


## 2.1. System Description

Distillation units are mainly utilized for the separation of fluid mixture components. The distillation column major components can be listed as follows:(i)
A vertical shell in which the separation of fluid substances is accomplished(ii)
A cascade of trays for improving component separation(iii)
A reboiler for maintaining the required heat energy for the distillation process(iv)
A condenser for liquefying the vapor leaving the column(v)
A reflux drum in which portion of the condensed liquid is recycled back to the vertical shell.The physical process considered in this research comprises two-coupled distillation columns as shown in Figure1. Thereby, ternary petrochemical mixtures can be separated. While the main glass column is composed of 40 bubble cap trays excluding those of the boiler and condenser, the side glass column contains 10 bubble cap trays. At the 22nd stage, the system intake port is located, at which the fluid mixture is to be fed. The three separated fluid components are to be extracted from the process, so that the heaviest and the lightest components can be provided from the bottom and the top of the main column, respectively, and the intermediate component is derived from the top of the side column.Figure 1
The two-coupled distillation columns process.


## 2.2. Mathematical Modeling of the Two-Coupled Distillation Column

Primarily, the manipulated variables are to be determined. Thus, the selected manipulated variables are the heat input to the reboiler (QE), the vapor flow rate in the vapor transfer line (SAB), the reflux ratio in the main column (RL1), and the reflux ratio in the second column (RL2) [59]. Regarding the system technical constraints, the input heat energy rate QE and the vapor flow rate SAB have to not exceed the values 8.2 KW and 3.95 m3/h, respectively, in order to avoid flooding of the smaller side column [59]. Secondly, the controlled values are selected to be the fluid temperature at four different stages. The selection of these stages is based on the following requirements:(i)
The temperatures should be sensitive to the major disturbances as changes in feed rate and feed concentrations.(ii)
The temperatures should exhibit a fairly linear behavior.(iii)
They should be a good indication of product quality.Accordingly, the temperature at the 11th, 30th, 34th, and 48th trays are chosen to form the controlled variables [59]. The mathematical model of realistic two-coupled distillation column, mentioned in [59], is to be used in this research. The transfer function matrix of the considered MIMO system is stated in (1).(1)Gs=2.61.69s+1-6.0983.5s+1-4.990.2s+14.5s+10.06s+10.0713.5s+17.321.05s+110.4s+10.14s+1-1.450.4s+1-1.570.23s+11.34s+10.2s+1-0.141.92s+14.60.53s+12.78s+10.09s+1-2.370.23s+12s+10.3s+1-2.71.75s+1-0.360.02s+12.47s+10.04s+12.110.92s+1-2.110.06s+12.38s+10.05s+1-1.752.16s+1-0.31.89s+14.35s+10.16s+1.As shown in the transfer function matrix, the physical system is characterized by the high interrelationship between all system inputs and outputs. As discussed in the introduction, the control structure adopted in this paper is based on the mathematical model decoupling, so that each output is independently controlled by a single input, forming four independent SISO control loops.Primarily, the noninteracting design is preceded by a relative gain analysis to determine the most suitable input-output pairings. In this regard, the relative gain array (RGA) method is used. Afterwards, a decoupling compensation network is designed for the reduction of residual interactions [11].


### 2.2.1. The Relative Gain Array

The RGA matrix forN×N system represented in (2) is used for the valuation of the input influences on each system output [16–19], so that output “i” is to be paired with input “j” for which γij is a positive value and as close to unity as possible.(2)RGAN×N=Gs=0·∗Gs=0-T=γ11γ12⋯γ1Nγ21γ22⋯γ2N⋮⋮⋮γN1γN2⋯γNN,where G(s) is the transfer function matrix of the MIMO process and “·∗” operator implies element by element multiplication.


### 2.2.2. Decoupling Compensation Network Design Procedure

The steady state decoupling compensation matrix illustrated in (3) is to be integrated into the given MIMO system in the way demonstrated in (4), so that interrelationships between each input and unpaired outputs is minimized [20, 21]. In this paper, the simplest form of the decoupling compensation matrix introduced in the previous study [11], which has unity diagonal elements, is replaced with the general form presented in (3). This modification introduces a remarkable improvement in the system response as discussed in Section 5. Figure 2 illustrates the decoupling compensation network scheme integrated with the required controllers.(3)Λss=λ11λ12λ13⋯λ1N-1λ1Nλ21λ22λ23⋯λ2N-1λ2N⋮⋮⋮⋮⋮λN1λN2λN3⋯λNN-1λNN,(4)Ys=GsΛssWs,where Ys is the output vector, Ws is the input vector, Gs is the transfer function matrix of the MIMO process, and Λss is the steady state decoupling compensation matrix of MIMO process.Figure 2
The decoupled MIMO system including controllers.


### 2.2.3. Fitness Function Design

As discussed in the former two sections, the determination of the adequate system input/output pairing as well as the elements’ values of the decoupling compensation matrix define the effectiveness of interrelationship minimization between system inputs and outputs. The selection of the proper input/output pairing is already tackled by the relative gain array concept. For the elements’ values of the decoupling compensation matrix to be estimated using optimization techniques, the fitness functions have to be designed that minimize the interrelationship between system inputs and their unpaired outputs.The four commonly used performance criteria for fitness function design are the integral absolute error (IAE), integral squared error (ISE), integral time-weighted squared error (ITSE), and Integral time-weighted absolute error (ITAE). It’s to mention that in the context of decoupling system design, the performance criteria are concerned with the minimization of the system outputs towards the unpaired inputs rather than control system error. Accordingly, the corresponding criteria are integral absolute output (IAO), integral squared output (ISO), integral time-weighted squared output (ITSO), and Integral time-weighted absolute output (ITAO). Although the criteria ISO and IAO grant less overshoot in the system dynamics, none of these criteria is adopted in this research due to the long settling time [56]. This drawback could be overcome by utilizing the ITAO. However, this criterion is not used in this research for the difficulty of its analytical tracking [60]. Thus, the ITSO performance criterion is employed in this work for ensuring the minimization of settling time without being confronted with unnecessary analytical complications.In the former research [11], the decoupling compensation matrix elements are estimated by applying the PSO technique, so that the integral squared outputs (ISOs) of unpaired outputs with respect to a specific input are minimized. On the other hand, the proposed technique in this paper is based on the minimization of the integral time-weighted squared outputs (ITSOs) by utilizing the BSO technique.The integral time-weighted squared outputs (ITSOs) for each input are calculated such thatifW1(s)=1/s and W2,W3,…,WN are zeroes, then(5)ITSO11=∫0∞tY12·dt,ITSO21=∫0∞tY22·dt,⋮ITSON1=∫0∞tYN2·dtand ifW2(s)=1/s and W1,W3,…,WN are zeroes, then(6)ITSO12=∫0∞tY12·dt,ITSO22=∫0∞tY22·dt,⋮ITSON2=∫0∞tYN2·dt.Consecutively, ifWN(s)=1/s and W1,W2,…,WN-1 are zeroes, then(7)ITSO1N=∫0∞tY12·dt,ITSO2N=∫0∞tY22·dt,⋮ITSONN=∫0∞tYN2·dt.Thus the fitness function for specific inputWj can be described as follows:(8)Fitnessj=∑i=1NITSOij,j=1,2,3,…,N,i≠q,where i is the specific output subscript, j is the specific input subscript, and q is the subscript of the output that has been paired with input j.


## 2.2.1. The Relative Gain Array

The RGA matrix forN×N system represented in (2) is used for the valuation of the input influences on each system output [16–19], so that output “i” is to be paired with input “j” for which γij is a positive value and as close to unity as possible.(2)RGAN×N=Gs=0·∗Gs=0-T=γ11γ12⋯γ1Nγ21γ22⋯γ2N⋮⋮⋮γN1γN2⋯γNN,where G(s) is the transfer function matrix of the MIMO process and “·∗” operator implies element by element multiplication.


## 2.2.2. Decoupling Compensation Network Design Procedure

The steady state decoupling compensation matrix illustrated in (3) is to be integrated into the given MIMO system in the way demonstrated in (4), so that interrelationships between each input and unpaired outputs is minimized [20, 21]. In this paper, the simplest form of the decoupling compensation matrix introduced in the previous study [11], which has unity diagonal elements, is replaced with the general form presented in (3). This modification introduces a remarkable improvement in the system response as discussed in Section 5. Figure 2 illustrates the decoupling compensation network scheme integrated with the required controllers.(3)Λss=λ11λ12λ13⋯λ1N-1λ1Nλ21λ22λ23⋯λ2N-1λ2N⋮⋮⋮⋮⋮λN1λN2λN3⋯λNN-1λNN,(4)Ys=GsΛssWs,where Ys is the output vector, Ws is the input vector, Gs is the transfer function matrix of the MIMO process, and Λss is the steady state decoupling compensation matrix of MIMO process.Figure 2
The decoupled MIMO system including controllers.


## 2.2.3. Fitness Function Design

As discussed in the former two sections, the determination of the adequate system input/output pairing as well as the elements’ values of the decoupling compensation matrix define the effectiveness of interrelationship minimization between system inputs and outputs. The selection of the proper input/output pairing is already tackled by the relative gain array concept. For the elements’ values of the decoupling compensation matrix to be estimated using optimization techniques, the fitness functions have to be designed that minimize the interrelationship between system inputs and their unpaired outputs.The four commonly used performance criteria for fitness function design are the integral absolute error (IAE), integral squared error (ISE), integral time-weighted squared error (ITSE), and Integral time-weighted absolute error (ITAE). It’s to mention that in the context of decoupling system design, the performance criteria are concerned with the minimization of the system outputs towards the unpaired inputs rather than control system error. Accordingly, the corresponding criteria are integral absolute output (IAO), integral squared output (ISO), integral time-weighted squared output (ITSO), and Integral time-weighted absolute output (ITAO). Although the criteria ISO and IAO grant less overshoot in the system dynamics, none of these criteria is adopted in this research due to the long settling time [56]. This drawback could be overcome by utilizing the ITAO. However, this criterion is not used in this research for the difficulty of its analytical tracking [60]. Thus, the ITSO performance criterion is employed in this work for ensuring the minimization of settling time without being confronted with unnecessary analytical complications.In the former research [11], the decoupling compensation matrix elements are estimated by applying the PSO technique, so that the integral squared outputs (ISOs) of unpaired outputs with respect to a specific input are minimized. On the other hand, the proposed technique in this paper is based on the minimization of the integral time-weighted squared outputs (ITSOs) by utilizing the BSO technique.The integral time-weighted squared outputs (ITSOs) for each input are calculated such thatifW1(s)=1/s and W2,W3,…,WN are zeroes, then(5)ITSO11=∫0∞tY12·dt,ITSO21=∫0∞tY22·dt,⋮ITSON1=∫0∞tYN2·dtand ifW2(s)=1/s and W1,W3,…,WN are zeroes, then(6)ITSO12=∫0∞tY12·dt,ITSO22=∫0∞tY22·dt,⋮ITSON2=∫0∞tYN2·dt.Consecutively, ifWN(s)=1/s and W1,W2,…,WN-1 are zeroes, then(7)ITSO1N=∫0∞tY12·dt,ITSO2N=∫0∞tY22·dt,⋮ITSONN=∫0∞tYN2·dt.Thus the fitness function for specific inputWj can be described as follows:(8)Fitnessj=∑i=1NITSOij,j=1,2,3,…,N,i≠q,where i is the specific output subscript, j is the specific input subscript, and q is the subscript of the output that has been paired with input j.


## 3. Controlling of Two-Coupled Distillation Column Process

In this section, the proposed BELBIC and the conventional PID controller are handled regarding the control structure.


### 3.1. Brain Emotional Learning Based Intelligent Controller (BELBIC) Model

As formerly stated, the considered MIMO system is to be split into several decoupled SISO systems with the minimum possible interrelationship between them. Thus, each of these systems can be independently controlled. In this research, the BELBIC adaptive control structure is utilized for each single control loop. This control structure is based on the functional model of brain emotional learning introduced by [28, 29]. Over the past decade, this control scheme has proven its robustness in many complex control applications [31–40]. Apart from the application in control systems, the computational model of brain emotional learning in its discrete and continuous form is discussed here, respectively. Afterwards, the methodology of the proposed BSO-BELBIC scheme, which assigns one BELBIC for each decoupled loop, is presented.The emotional learning computational model designed by [28, 29] is graphically illustrated in Figure 3. In mammalian brains, the emotional learning process occurs in a part of the brain called the limbic system, which consists of four main components corresponding to the amygdala, orbitofrontal cortex, thalamus, and the sensory cortex. As shown in the figure, sensory input signals are primarily passed to the thalamus model, which relay the sensory information from the peripheral sensory systems to the sensory cortices. Moreover, the sensory input maximum value is passed directly to the amygdala for the fast response to be insured. Then, the available sensory data are to be processed by sensory cortex model. Hence, highly analyzed data are to be sent to the amygdala and orbitofrontal cortex models. The emotional evaluation of stimuli and the formulation of long-term memories are carried out by the amygdala. Finally, the orbitofrontal cortex is supposed to inhibit inappropriate responses from the amygdala.Figure 3
Graphical depiction of the brain emotional learning (BEL) process [28, 29].The outputs of the model two major components amygdala and orbitofrontal cortex are described in (9) and (10), respectively. The feedback element MO′ is defined in (11) as the subtraction of the orbitofrontal cortex inhibitory outputs (Oi) from the summation of amygdala nodes (Ai) excluding Ath node. As illustrated in relation (12), the output (MO) of the brain emotional learning model (BEL) constitutes the subtraction of the orbitofrontal cortex inhibitory outputs (Oi) from the summation of amygdala nodes (Ai) including the Ath node.(9)Ai=SiGAi,(10)Oi=SiGOi,(11)MO′=∑iAi-∑iOiexcludingAthnode.,(12)MO=∑iAi-∑iOiincludingAthnode.,where Si forms the ith sensory input, and GA and GO are the plastic connection weights of the amygdala and orbitofrontal cortex, respectively. These plastic connection weights are responsible for the emotional change towards specific object characteristics. Thus, they constitute the adaptive component in the model structure, as their values are to be updated continuously. While the formulas (13) and (14) represent the discrete form for the change in plastic connection weights, (15) and (16) constitute the continuous one. It is to mention that the continuous form is utilized in this research.(13)ΔGAi=αSimax⁡0,Rew-∑jAj,(14)ΔGOi=βSiMO′-Rew,(15)dGAidt=α·Si·Rew-Ai,(16)dGOidt=β·Si·Ai-Oi-Rew,where α and β are learning rate constants, and the symbol Rew forms the reward signal. The operator “max” in the formula (13) is the responsible for maintaining the monotonic learning change of amygdala. This characteristic models the incapability of the amygdala to unlearn the formerly learned emotions [28, 29]. In the continuous form, the operator “max” is eliminated for analytical simplicity [61].The BELBIC internal structure as well as its interface with the controlled physical plant is illustrated in Figure4. As shown, sensory input block as well as the reward signal builder manipulates orbitofrontal cortex and amygdala based on the control error signal according to (17) and (18), respectively.(17)S=K·e,(18)Rew=Kp·e+Ki·∫e·dt+Kd·dedt,where K, Kp, Ki, and Kd besides the learning rate constants (α and β) constitute the controller parameters, which characterize the controlled dynamic system behavior.Figure 4
Control system configuration using BELBIC.The BELBIC parameterization is one of the current challenges confronting such a novel control structure. In this regard, the utilization of optimization techniques like Particle Swarm Optimization (PSO) has shown a robust behavior [37]. In the frame of this research, the feasibility of applying the Bacterial Swarm Optimization (BSO) algorithm regarding the tuning process of BELBIC parameters is to be investigated. As mentioned in the previous section, the integral time-weighted squared errors (ITSEs) of all control loops are selected to be minimized rather than the integral-squared-errors (ISEs) for its faster settling response. Accordingly, the fitness function is represented in (19).(19)Fitness  function=∫0∞tT30desired-T302·dt+∫0∞tT11desired-T112·dt+∫0∞tT34desired-T342·dt+∫0∞tT48desired-T482·dt.Twenty-four parameters (six for each loop) should be tuned simultaneously with the aim of minimizing the fitness function.


### 3.2. PID Control

The main terms that constitute the conventional PID controller are the proportional, the integral, and the derivative terms. The three terms are added to each other as shown in Figure5. The transfer function of the conventional PID controller is stated in (20).(20)Gcs=Kp+Kis+sKd,where Kp, Ki, and Kd are proportional gain, integral gain, and a derivative gain, respectively.Figure 5
Block diagram of the conventional PID controller.The Bacterial Swarm Optimization (BSO) algorithm is utilized regarding the parameterization of PID controllers with the same policies as BELBICs.


## 3.1. Brain Emotional Learning Based Intelligent Controller (BELBIC) Model

As formerly stated, the considered MIMO system is to be split into several decoupled SISO systems with the minimum possible interrelationship between them. Thus, each of these systems can be independently controlled. In this research, the BELBIC adaptive control structure is utilized for each single control loop. This control structure is based on the functional model of brain emotional learning introduced by [28, 29]. Over the past decade, this control scheme has proven its robustness in many complex control applications [31–40]. Apart from the application in control systems, the computational model of brain emotional learning in its discrete and continuous form is discussed here, respectively. Afterwards, the methodology of the proposed BSO-BELBIC scheme, which assigns one BELBIC for each decoupled loop, is presented.The emotional learning computational model designed by [28, 29] is graphically illustrated in Figure 3. In mammalian brains, the emotional learning process occurs in a part of the brain called the limbic system, which consists of four main components corresponding to the amygdala, orbitofrontal cortex, thalamus, and the sensory cortex. As shown in the figure, sensory input signals are primarily passed to the thalamus model, which relay the sensory information from the peripheral sensory systems to the sensory cortices. Moreover, the sensory input maximum value is passed directly to the amygdala for the fast response to be insured. Then, the available sensory data are to be processed by sensory cortex model. Hence, highly analyzed data are to be sent to the amygdala and orbitofrontal cortex models. The emotional evaluation of stimuli and the formulation of long-term memories are carried out by the amygdala. Finally, the orbitofrontal cortex is supposed to inhibit inappropriate responses from the amygdala.Figure 3
Graphical depiction of the brain emotional learning (BEL) process [28, 29].The outputs of the model two major components amygdala and orbitofrontal cortex are described in (9) and (10), respectively. The feedback element MO′ is defined in (11) as the subtraction of the orbitofrontal cortex inhibitory outputs (Oi) from the summation of amygdala nodes (Ai) excluding Ath node. As illustrated in relation (12), the output (MO) of the brain emotional learning model (BEL) constitutes the subtraction of the orbitofrontal cortex inhibitory outputs (Oi) from the summation of amygdala nodes (Ai) including the Ath node.(9)Ai=SiGAi,(10)Oi=SiGOi,(11)MO′=∑iAi-∑iOiexcludingAthnode.,(12)MO=∑iAi-∑iOiincludingAthnode.,where Si forms the ith sensory input, and GA and GO are the plastic connection weights of the amygdala and orbitofrontal cortex, respectively. These plastic connection weights are responsible for the emotional change towards specific object characteristics. Thus, they constitute the adaptive component in the model structure, as their values are to be updated continuously. While the formulas (13) and (14) represent the discrete form for the change in plastic connection weights, (15) and (16) constitute the continuous one. It is to mention that the continuous form is utilized in this research.(13)ΔGAi=αSimax⁡0,Rew-∑jAj,(14)ΔGOi=βSiMO′-Rew,(15)dGAidt=α·Si·Rew-Ai,(16)dGOidt=β·Si·Ai-Oi-Rew,where α and β are learning rate constants, and the symbol Rew forms the reward signal. The operator “max” in the formula (13) is the responsible for maintaining the monotonic learning change of amygdala. This characteristic models the incapability of the amygdala to unlearn the formerly learned emotions [28, 29]. In the continuous form, the operator “max” is eliminated for analytical simplicity [61].The BELBIC internal structure as well as its interface with the controlled physical plant is illustrated in Figure4. As shown, sensory input block as well as the reward signal builder manipulates orbitofrontal cortex and amygdala based on the control error signal according to (17) and (18), respectively.(17)S=K·e,(18)Rew=Kp·e+Ki·∫e·dt+Kd·dedt,where K, Kp, Ki, and Kd besides the learning rate constants (α and β) constitute the controller parameters, which characterize the controlled dynamic system behavior.Figure 4
Control system configuration using BELBIC.The BELBIC parameterization is one of the current challenges confronting such a novel control structure. In this regard, the utilization of optimization techniques like Particle Swarm Optimization (PSO) has shown a robust behavior [37]. In the frame of this research, the feasibility of applying the Bacterial Swarm Optimization (BSO) algorithm regarding the tuning process of BELBIC parameters is to be investigated. As mentioned in the previous section, the integral time-weighted squared errors (ITSEs) of all control loops are selected to be minimized rather than the integral-squared-errors (ISEs) for its faster settling response. Accordingly, the fitness function is represented in (19).(19)Fitness  function=∫0∞tT30desired-T302·dt+∫0∞tT11desired-T112·dt+∫0∞tT34desired-T342·dt+∫0∞tT48desired-T482·dt.Twenty-four parameters (six for each loop) should be tuned simultaneously with the aim of minimizing the fitness function.


## 3.2. PID Control

The main terms that constitute the conventional PID controller are the proportional, the integral, and the derivative terms. The three terms are added to each other as shown in Figure5. The transfer function of the conventional PID controller is stated in (20).(20)Gcs=Kp+Kis+sKd,where Kp, Ki, and Kd are proportional gain, integral gain, and a derivative gain, respectively.Figure 5
Block diagram of the conventional PID controller.The Bacterial Swarm Optimization (BSO) algorithm is utilized regarding the parameterization of PID controllers with the same policies as BELBICs.


## 4. Bacterial Swarm Optimization Algorithm

Bacterial Swarm Optimization Algorithm (BSO) forms a hybrid between two efficient optimization techniques. These algorithms are the Particle Swarm Optimization (PSO) and the Bacterial Foraging Optimization Algorithm (BFOA).PSO is a stochastic optimization approach inspired from the behavior of the flock of birds, insects, and fish. Every single particle in the search space adjusts its own direction based on its own experience as well as the experience of the most successful particle in the swarm [43, 44]. Nevertheless, the optimization technique based on the PSO algorithm may lead to entrapment in local optimum solution rather than catch the global one due to the rapidity and simplicity of the algorithm [62].On the other hand, BFOA is a new bio-inspired algorithm depending on foraging behavior ofEscherichia coli (E. coli) bacteria [45]. Bacteria have the tendency to group around the nutrient-rich regions by the activity named chemotaxis. The bacteria which fail to reach nutrient-rich areas may die due to the nutrient lack. However, the ones that survived reproduce the next generation in nutrient-rich areas. Once the current living environment becomes inconvenient to the bacteria, it tends to disperse randomly to search for an alternative environment. Consequently, the optimization technique that simulates the foraging behavior of these bacteria requires a long time for achieving the global optimum solution due to the dependence on random search directions [57].The time consumed by BFOA finding the global optimum solution can be reduced by granting theE. coli bacteria the ability of exchanging social information. This ability is inherited from the PSO technique forming the BSO algorithm. Therefore, BSO algorithm requires less time for the optimum solution determination, while maintaining the BFOA ability in finding a new solution with elimination and dispersal. Thus, the BSO algorithm solves the insufficient scattering problem that confronted the PSO algorithm. Furthermore, the chemotaxis step of BSO technique safeguards against the PSO shortcoming regarding the weak search ability.As mentioned, the BSO technique is utilized for the parameterization of the decoupling compensation network as well as the BELBIC. As an overview, the procedure of applying the proposed technique on the case study is demonstrated in the process chart in Figure6. Afterwards, the basic flowchart of the BSO algorithm is presented in Figure 7. Moreover, the pseudocode of the BSO algorithm is stated delivering more details.Figure 6
Process chart for the main steps of applying the proposed technique on the case study.Figure 7
Flowchart of BSO algorithm.The pseudocode of the BSO algorithm is as follows.Step 1.
The initialization of all stated parameters:p, S, Sr, Nc, Ns, Nre, Ned, Ped, C(i)(i=1,2,…,S), Delta, w, C1, C2, R1, and R2, where(i)
p is dimension of search space,(ii)
S is the number of bacteria,(iii)
Sr is the number of bacteria splits per generation,(iv)
Nc is the number of chemotactic steps,(v)
Ns is the limits of the length of a swim,(vi)
Nre is the reproduction steps number,(vii)
Ned is the amount of elimination-dispersal events,(viii)
Ped is the elimination-dispersal probability,(ix)
C(i) is the bacteria step size length,(x)
Delta is the direction of each bacteria,(xi)
w is the weight of inertia,(xii)
C1 is the weight of local information,(xiii)
C2 is the weight of global information,(xiv)
R1, R2 are two Random numbers.Step 2.
Loop of elimination and dispersal:l=l+1.Step 3.
Loop of reproduction:k=k+1.Step 4.
Loop of chemotaxis:j=j+1.
Substep  4.1. For i=1,2,…,S each bacterium i moves a chemotactic step as follows.(a)
Calculate cost function,J(i,j,k,l).(b)
LetJlast=J(i,j,k,l) to save the current value of the cost function in order to be able to compare it with the one determined in the next swim.(c)
LetJlocal(i,j)=Jlast; the better cost per each bacteria is going to be chosen to be the local best Jlocal.(d)
Update positionP(i,j+1,k,l)=P(i,j,k,l)+C(i)∗Delta(i).(e)
Calculate cost function,J(i,j+1,k,l).(f)
Swim
(i)
Letm=0 (counter for swim length).(ii)
Whilem<Ns (if have not climbed down too long).
(1)
Letm=m+1.(2)
IfJ(i,j+1,k,l)<Jlast (if doing better), let Jlast=J(i,j+1,k,l) and let(21)Pi,j+1,k,l=Pi,j+1,k,l+Ci∗Deltaiand use thisP(i,j+1,k,l) to calculate the new J(i,j+1,k,l).(3)
For each bacteria estimate the current position and local cost(22)Pcurrenti,j+1=Pi,j+1,k,l,Jlocali,j+1=Ji,j+1,k,l.(4)
Else, let(23)Pcurrenti,j+1=Pi,j+1,k,l,Jlocali,j+1=Ji,j+1,k,l,m=Ns.(5)
While statement end.(g)
Go to next bacterium(i+1) if i≠S (i.e., go to (b) to execute the next bacterium).
Substep  4.2. For each bacteria estimate the local best position (PLbest) and global best position (PGbest).
Substep  4.3. For each bacteria estimate the new direction as(24)V=w∗V+C1∗R1PLbest-Pcurrent+C2∗R2PGbest-Pcurrent,Delta=V.Step 5 (ifj<Nc, go to Step 4).
In this situation keep on chemotaxis since the life of the bacteria is not terminated.Step 6 (reproduction).
Substep  6.1. For k and l and for each i=1,2,…,S, the health of the bacterium i is given by(25)Jhealthi=∑j=1Nc+1Ji,j,k,l.Sort bacteria cost Jhealth in ascending order (greater cost indicates lower health).
Substep  6.2. The Sr bacteria with the lowest Jhealth values split and the rest of the Sr bacteria die.Step 7 (ifk<Nre, return to Step 3).
In this instance, the specified maximum number of reproduction steps is not over; therefore the bacteria begin a new generation of a chemotactic loop.Step 8 (elimination dispersal).
Fori=1,2,…,S with probability Ped, eliminate and then disperse one to a random place. If l<Ned, then go to Step 2, otherwise end.


## 5. Simulation and Results

The decoupled physical system and its controllers are designed and simulated with MATLAB™ Simulnik®.The RGA method is applied on the mathematical model of the two-coupled distillation column process, and the resulted matrix is given by(26)RGA4×4=-0.04290.56290.40830.07171.45590.5558-0.9138-0.0980-0.4647-3.11104.48320.09260.05172.9923-2.97770.9337.According to the RGA matrix the suitable input-output pairing is(27)T11-SAB,T30-QE,T34-RL1,T48-RL2.Depending on the suitable pairing, the fitness functions of the 1st, 2nd, 3rd, and 4th input are described in (28), (29), (30), and (31), respectively.(28)
F
i
t
n
e
s
s
1
=
I
T
S
O
11
+
I
T
S
O
31
+
I
T
S
O
41
=
16185764.6415
∗
λ
11
2
+
-
31208202.2860
∗
λ
11
∗
λ
21
+
-
29085704.6263
∗
λ
11
∗
λ
31
+
-
2104353.3818
∗
λ
11
∗
λ
41
+
23626420.1763
∗
λ
21
2
+
40518661.8661
∗
λ
21
∗
λ
31
+
1053225.6938
∗
λ
21
∗
λ
41
+
17625377.5018
∗
λ
31
2
+
1142698.7804
∗
λ
31
∗
λ
41
+
112317.1448
∗
λ
41
2
(29)
F
i
t
n
e
s
s
2
=
I
T
S
O
22
+
I
T
S
O
32
+
I
T
S
O
42
=
39587638.7663
∗
λ
12
2
+
-
25965747.1866
∗
λ
12
∗
λ
22
+
-
27602372.6861
∗
λ
12
∗
λ
32
+
-
3313543.6120
∗
λ
12
∗
λ
42
+
6085664.7540
∗
λ
22
2
+
12367844.1767
∗
λ
22
∗
λ
32
+
1689163.6136
∗
λ
22
∗
λ
42
+
6408628.3541
∗
λ
32
2
+
1716766.7412
∗
λ
32
∗
λ
42
+
119596.6264
∗
λ
42
2
(30)
F
i
t
n
e
s
s
3
=
I
T
S
O
13
+
I
T
S
O
23
+
I
T
S
O
43
=
32387849.9550
∗
λ
13
2
+
-
30918315.4051
∗
λ
13
∗
λ
23
+
-
28156040.3138
∗
λ
13
∗
λ
33
+
-
1472985.7141
∗
λ
13
∗
λ
43
+
21869260.2750
∗
λ
23
2
+
36396234.4688
∗
λ
23
∗
λ
33
+
403039.5905
∗
λ
23
∗
λ
43
+
15212859.0003
∗
λ
33
2
+
390512.5127
∗
λ
33
∗
λ
43
+
57318.4117
∗
λ
43
2
(31)
F
i
t
n
e
s
s
4
=
I
T
S
O
14
+
I
T
S
O
24
+
I
T
S
O
34
=
40741561.5731
∗
λ
14
2
+
-
37368063.2195
∗
λ
14
∗
λ
24
+
-
36883365.2610
∗
λ
14
∗
λ
34
+
-
2495963.4210
∗
λ
14
∗
λ
44
+
22451663.5575
∗
λ
24
2
+
39102720.2746
∗
λ
24
∗
λ
34
+
623243.6320
∗
λ
24
∗
λ
44
+
17326594.9351
∗
λ
34
2
+
837511.9391
∗
λ
34
∗
λ
44
+
77118.8721
∗
λ
44
2
.The BSO algorithm parameters’ values which are utilized to implement an optimized decoupling network as well as an optimized BELBIC are summarized in Table1. The resulting values of the steady state decoupling compensation elements that minimize the above fitness functions are summarized in Table 2.Table 1
Gains of BSO.
Parameter Symbol BSO for decoupling network implementation BSO for BELBIC implementation Number of bacteria in the population S 50 50 Dimension of search space p 4 24 Maximum number of swim length N s 4 4 Maximum number of chemotactic steps N c 100 20 Number of reproductive steps N r e 4 2 Number of elimination dispersal events N e d 2 2 Elimination dispersal probability P e d 0.25 0.25 Step size C ( i ) 0.05 0.05 Cognitive factor C 1 1.2 1.2 Social acceleration factors C 2 0.5 0.5 Momentum/inertia w 0.9 0.9Table 2
Resulting values of steady state decoupling compensation elements based on BSO.
Fitness function Final value of fitness function Values of steady state decoupling elements (λs) Fitness1 9.8506 λ 11 = 0.2453 λ 21 = - 0.4755 λ 31 = 0.7210 λ 41 = 0.8564 Fitness2 0.5941 λ 12 = - 0.0195 λ 22 = - 0.1090 λ 32 = - 0.0967 λ 42 = 1.1934 Fitness3 7.6138 λ 13 = - 0.0481 λ 23 = 0.6249 λ 33 = - 0.7904 λ 43 = - 0.1218 Fitness4 0.3040 λ 14 = - 0.0026 λ 24 = 0.1492 λ 34 = - 0.1790 λ 44 = 0.3273The resulting best gains’ values for different BELBICs that minimize the summation of the integral time-weighted squared errors (ITSEs) for different decoupled loops are presented in Table3. The best gains’ values of the designed conventional PID controllers are shown in Table 4. The PID controllers’ parameters are determined utilizing the same algorithm utilized for the design of BELBICs. In this regard, proposed BSO algorithm with the same gains given in Table 1 is used to minimize the same fitness function.Table 3
The proper gains of the BELBICs optimized by the BSO for different loops.
Loop Gains α β K K p K i K d ( Q E , T 30 ) 0.2321 0.8914 −0.4830 3.4765 0.0028 −0.6471 ( S A B , T 11 ) 0.3418 −1.3651 0.3936 −21.5658 −1.0949 −0.7771 ( R L 1 , T 34 ) 87.3587 −8.9268 −0.1020 29.4672 17.3183 14.6593 ( R L 2 , T 48 ) 17.4282 22.6201 0.2650 −62.2460 −0.6042 10.8162Table 4
The proper gains of the PID controllers optimized by the BSO for different loops.
Loop Gains K p K i K d ( Q E , T 30 ) 1.9662 0.3524 0.7445 ( S A B , T 11 ) 1.0643 1.9718 0.2621 ( R L 1 , T 34 ) 1.6276 1.4202 −0.0366 ( R L 2 , T 48 ) 0.0784 −3.0889 0.8859The dynamic behavior of the system is analyzed based on its step response based on the sequential step input shown in Figure8. In Figure 9, the decoupled system response on the given step sequence is illustrated. As shown, the experimental results demonstrate the efficiency of the designed decoupling technique. Thus, the interrelationship between system inputs and their unpaired outputs is noticeably minimized compared to the formerly deployed techniques [11]. Accordingly, the spikes in the response of the decoupled system outputs on the unpaired inputs are significantly reduced. The step response of the designed closed loop control system is shown in Figure 10.Figure 8
Step changes in system inputs [37].Figure 9
The outputs of different decoupled loops in case of no controllers.Figure 10
The outputs of different decoupled loops in the presence of the controller.For verification purposes, the response of the designed control scheme, which is based on the minimization of ITSEs of all control loops utilizing the BSO algorithm, is compared to that of the latest research that considers the same application. It is to mention that the control scheme of the previous research is mainly based on the minimization of the ISE using the PSO algorithm [37]. In Figure 11, the proposed technique is compared to the former study in terms of step response. Moreover, the steady state errors for both controllers are stated in Table 5.Table 5
The steady state errors for all loops after control.
Loop Steady state errors Former technique Proposed technique ( Q E , T 30 ) 0.087510 0.1189 ( S A B , T 11 ) 0.019915 0.0117 ( R L 1 , T 34 ) 0.036152 0.0041 ( R L 2 , T 48 ) 0.189710 0.0207Figure 11
Comparison between the step response of the proposed technique and that of the formerly studied technique.As shown, the control scheme presented in this research offers a valuable improvement in the last three control loops (T11, T34, and T48) in terms of minimizing steady state errors.The step response of the proposed BELBICs and conventional PID controllers is compared in the presence of a disturbance step with the final value of “1” at the 500th second in the first and the third decoupled loop. The simulation results are presented in Figure12. In the figure it is clearly recognizable that the robustness of the proposed BELBICs is higher than that of the conventional PID controllers regarding the handling of the unexpected disturbance. The controlled system using BELBICs is better damped as well. On the other side, the PID controllers in all control loops achieve remarkably less steady state error.Figure 12
Comparison between the step response of the designed BELBIC and that of the PID controller in the presence of disturbance att=500 seconds.For comparison purpose, the PID controllers are designed by utilizing the PSO technique for minimizing the summation of the integral time-weighted squared errors (ITSEs) of system control loops. This controller is to be compared with the one designed using the BSO algorithm regarding the integral time-weighted squared errors (ITSEs) for each control loop, which are listed in Table6. The remarkable difference between both algorithms regarding ITSEs indicates that the BSO technique is more efficient in determining the global best solution in the field of MIMO control system.Table 6
Comparison between PID controllers designed by PSO and BSO algorithms regarding integral time-weighted squared error ITSE.
Loop ITSE PSO technique BSO technique ( Q E , T 30 ) 11.5703 3.2552 ( S A B , T 11 ) 0.4473 1.8333 ( R L 1 , T 34 ) 2.7937 1.8592 ( R L 2 , T 48 ) 11.8501 5.9165  Summation 26.6614 12.8642


## 6. Conclusions

The challenge of decoupling and controlling higher order multi-input multioutput (MIMO) process is tackled in this research. The Bacterial Swarm Optimization (BSO) technique is used to develop an optimized scheme for decoupling highly interactive 4 input/4 output two-coupled distillation column processes. The scheme consists of two stages. In the first stage, the optimum group of fitness functions is determined through the analysis of precalculated proper pairing, which is based on the derived relative gain array (RGA). The derivation of the RGA is based on the transfer function matrix of the physical process. In the second stage, the values of decoupling compensation elements (λs) that minimize the interactions are estimated based on the formerly driven fitness functions. Designing the decoupled system based on the general form of the decoupling compensation matrix showed a remarkable improvement in the system dynamics compared to the utilization of the simple form of the compensation matrix. The simulation results showed the efficiency of the proposed BSO technique in estimating steady state decoupling compensation elements values. For control purpose a scheme of Brain Emotional Learning Based Intelligent Controller (BELBIC) is designed and optimized using BSO algorithm to obtain the optimal values of controllers’ parameters. Furthermore, PID controllers are developed using the same optimization technique, in order to validate the robustness of the BELBIC. The robust control behavior of the designed BELBIC-based control scheme is validated in the simulation results. The BELBIC designed using the BSO algorithm showed a remarkable improvement in the transient and steady state errors of the last three control loops compared with the controller designed utilizing the PSO technique.



---

*Source: 8985425-2016-10-11.xml*