Research Article  Open Access
Quantitative Feedback Control of Multiple Input Single Output Systems
Abstract
This paper presents a robust feedback control solution for systems with multiple manipulated inputs and a single measurable output. A structure of parallel controllers achieves robust stability and robust disturbance rejection. Each controller uses the least possible amount of feedback at each frequency. The controller design is carried out in the Quantitative Feedback Theory framework. The method pursues a smart load sharing along the frequency spectrum, where each branch must either collaborate in the control task or be inhibited at each frequency. This reduces useless fatigue and saturation risk of actuators. Different examples illustrate the ability to deal with complex control problems that current MISO methodologies cannot solve. Main control challenges arise due to the uncertainty of plant and disturbance models and when a fastslow hierarchy of plants cannot be uniquely established.
1. Introduction
This paper deals with systems where multiple inputs are used to govern a single output. Although scientific literature refers to them with diverse names, here they will simply be called MISO systems. In some cases, each individual output is accessible, as in distributed energy generation systems [1], decentralised production systems [2], or unmanned fleets [3]. Usually, this leads to complex multiloop control structures [4]. However, a lot of engineering MISO systems lack physical individual outputs or sensors to measure them. Such systems are common in process industry [5], where the MISO control sometimes pursues the management of the global production system [6], whereas other times it governs lowlevel process variables (typical examples are two pumps or a pump and a valve, used as actuators in pressure or flow control [7–9]). In other cases the MISO control attends certain subsystems in a process, such as chemical reactors [10–12] or biological reactors [13–15]. More specific usages can be found in drying sections of paper machines [7, 16] or in aerobic digesters of waste water treatment plants [17]. Heat exchangers [18–20], chemical reactors in polymerization processes [9, 18], or distillation columns [9, 20–23] are repeated references in the scientific literature as MISO control applications. The automotive industry has also adopted these principles, firstly for the government of internal combustion engines [24–27] and recently for HCCI (Homogeneous Charge Compression Ignition) engines [28–30]. Another area devoted to MISO control is the consumer electronics, and particularly the massive data storage devices [31–33]. And finally, biological engineering applications can be found in [18, 34].
Within those MISO systems with nonindividual measurable outputs, the control strategies can be divided into noncollaborative and collaborative ones. Noncollaborative control selects a plant inside a battery of them, which covers a wide range of operating points for the output. The selection criterion is based on the stationary capacity of each plant. Thus, the control law is designed for an equivalent SISO system. A selector splits online the control action to the plant or plants with capacity to regulate the output in the actual operating point. The splitrange control [35, 36] is the most representative of this methodology. A simpler method reduces to a pure SISO control system, which closes a single feedback loop around a plant. The inputs to the other plants are manipulated manually or are left constant [37].
On the other hand, collaborative strategies benefit from the dynamic strengths of each plant to improve the controlled output performance considering the restrictions of manipulated variables and individual outputs. In the scientific literature collaborative MISO systems appear under diverse designations: VPC (Valve Position Control) [8, 9, 21–23], habituating control [10–12, 17, 18, 34, 38], mainvernier control [39–42], cooperativefeedback control [20, 43, 44], load sharing control [2, 45–47], or PQ design method [27, 48, 49]. Midranging control [7, 13, 15, 16, 35, 50–54] is frequently found in process control literature.
The controlled collaboration of parallel plants can be performed by serial (Figure 1(a)) or parallel (Figure 1(b)) structures of controllers. Serial collaboration is based on a qualitative organization of plants: the main loopcorresponds to the fastest plant, the first one to react to output deviations. The vernier loops are arranged to produce a chained intervention of gradually slower plants. The last loop, which incorporates the slowest plant, is in charge of the steady state. In this way, each system collaborates to the extent it provides benefits, avoiding useless actuator fatigue and risk of saturation. This is also the goal of parallel structures. However, the involvement of each plant is determined by its preceding controller, instead of by its position in the overall scheme. Thus, the design of each openloop transfer functionbecomes more complex, since it has to define both the actuation domain and the action itself. However, unlike serial structures, parallel designs can cope with systems where the fastslow hierarchy cannot be uniquely established.
(a) Serial structure or mainvernier
(b) Parallel structure
In addition to all this, a proper design must pay some attention to the pervasive presence of uncertainty. In this sense, some solutions for the MISO control problem come from the Quantitative Feedback Theory (QFT). Horowitz [4, 55] did not pay much attention to the MISO problem with no individual measurable outputs. His work suggests the reduction of the MISO system to a SISO equivalent by plant summation. Then, a master controller (in Figure 1(b)) is designed for the equivalent plant, and no individual controller is used (). The master control strategy can extract the maximum dynamic potential of each plant under certain restrictions for their phases. Nevertheless, its main drawback is that the same control action reaches all plant actuators. This entails using more feedback than necessary at certain frequencies on each branch, which leads to different negative effects in reallife actuators. In particular, an excess of noise amplification arises in the control actions to the slowest plants and steadystate offsets reduce the available actuation range for the fastest ones. Another drawback of plant summation strategy is that feedback is allocated by the plant instead of by the designer, losing some flexibility. This might have some importance in certain systems [7] whose high operation expenses suggest to arrange plant interventions according to monetary criteria rather than dynamic ones. Further QFT developments include the approaches in [25, 56], which detail a collaborative control of parallel plants with uncertainties and restrictions for particular applications, and the approach in [3], which adds feedback loops of actuations inside a parallel structure as described in [57].
However, a general robust methodology to deal with MISO collaborative systems through QFT tools has not been presented yet. Such is the goal of this paper. In particular, it focuses on the disturbance rejection or regulation problem. The new technique is fitted to the parallel structure. This is the most versatile one since other arrangements can always be transformed into it while the inverse transformation is not always feasible. Besides, some examples will show that mainvernier design philosophy cannot cope with certain types of MISO plants. The new methodology seeks a quantitative division of feedback amongst parallel loops in the frequency domain. In absence of restrictions, this distribution leads to the accomplishment of the specifications with minimum use of feedback, that is, each loop only contributes at those frequencies where its plant favours the output performance.
2. MISO QFT Design Methodology
2.1. Robust MISO Regulation Problem
Consider Figure 1(b) with. Each parallel plant,, defines the effect of each manipulated variableover the single measurable output. The plantdefines the way in which the nonmeasurable disturbancedeviates the single outputfrom its desired constant setpoint. A set of parallel controllers,are designed to minimize such deviation. Eachfixes its correspondingbased on a common output measurement, which is corrupted by the noise signal. For this regulation problem, the closed loop functions are whereare the individual openloop transfer functions and is the global openloop transfer function.
All plants are assumed to present parametric uncertainty. Let us defineas the vector of all parameters appearing in the transfer functions. Eachvaries within certain lower and upper limits. Therefore, the uncertainty vectorbelongs to a hyperrectangle incalled the uncertainty space ; that is,
The purpose of QFT control is to enforce the specifications for all elements in the uncertainty space. In this case, the specifications include robust stability, and robust disturbance rejection, whereandare upper tolerances for each closedloop frequency response.
As long as feedback control guarantees expected specifications at a certain price, QFT advocates for that control solution with the minimum amount of feedback, which will be denoted as the QFT optimum. In SISO control, the bounds express the minimal gain requirements on the nominal openloop transfer function for the accomplishment of the specifications. Consequently, an QFT optimal loopshaping tries first a narrow fulfilment of the bounds, and then an abrupt rolloff to minimize the presence of noise in the control signals [58].
In the MISO case, the search for the QFT optimum is not so obvious. There is a minimal gain of the nominal openloop global functionthat meets (4) and (5). However, it can be achieved by several combinations ofin (2). That combination of minimum gains of eachat each frequency is the solution that demands the minimum amount of feedback. Then, one option is splitting the global problem in several designsin such a way that each bound arrangement shows the minimum amount of feedback solution. If so, each loopshaping reduces to a narrow fulfilment of the bounds and a maximal reduction of gain in the frequencies where no restrictions are imposed.
In practice, the new methodology looks for a wise distribution of feedback, which exploits the power of each plant frequency response to achieve controllers of minimum gain. This minimizes the presence of noise in the control signals to the slowest plants, which usually exhibit the most powerful actuators and provide the bias action. In this way, faster plant actuators preserve their full range of operation for output regulation in response to disturbances.
2.2. Parallel Plants and the Amount of Feedback
The core of the method is to quantify the amount of feedback, which is evaluated with the controller gain. For a better understanding, let us take a single frequency, exclude the uncertainty, and consider solely the specification of robust disturbance rejection. Then, certainwill be a solution to (5). Two parallel structures from Figure 1(b) will be studied: a master control, which sintonizes, and does each, and a strictly parallel control, which sintonizes eachand does.
To achieve, we are firstly comparing the controller gain necessities for plants of similar characteristics. For simplification purposes, let us take two plants . As expected, if the system power to regulate the output was increased fromto, feedback necessities would reduce in the same proportion fromto, independently if a master controller or two parallel controllers were used.
Secondly, control necessities to achieveare being evaluated for plants with very different gains at the same phase. Let us take two plants as . Note that, ifplant worked alone, the feedback demand would be considerably inferior to the case ofworking alone:. Therefore, there is no sense in using the actuation unless you consider other criteria different from meeting the performance with the minimum feedback. In the case of a master controller, its control demand is, which approaches to. Thus, is the only plant contributing to the performance andactuator is being excited unnecessarily. A parallel controller structure can overcome this drawback as follows. Basing on the formulation,can be achieved by different combinations of. However, as long aspowerfulness is negligible, any attempt of collaboration withto the performance will require a huge control demand in the second branch. Then, the best option is to achieve the performance with thebranch, whose control demand will be. Simultaneously, tuningavoids that useless signals reach the actuators.
One of the above two situations (plants that should collaborate or not) arises in the loopshaping of controllers at each frequency. In a first stage of the proposed methodology, special QFT bounds guide the designer to define when several plants are similar enough to collaborate at certain frequency or when they are different enough such that some must take the whole responsibility and others must be inhibited. Model uncertainty and the whole set of specifications are included in the method. In a second stage, the feedback is shared accordingly through the sequential loopshaping of controllers in the parallel structure.
2.3. First Stage: Load Sharing Planning at Each Frequency of Design
The proposed method firstly evaluates the amount of controller gainthat each plantwould demand to fulfil the specifications (4), and (5) by itself, that is, assuming. The required gain is expressed graphically by a QFT bound at each frequency, denoted by. Note that, unlike common QFT bounds, these are bounds on the controller gain, not on the nominal openloop gain, in the vertical axis. This is so, because their purpose is not to serve as design guidelines, but to allow a graphical comparison among the controller gains that are being demanded at subsequent design phaseby each plant. Once all the bounds for a single frequencyare computed and a design phaseis selected for comparisons, a conclusion is reached according to the following criteria.(a)If one of the bounds, say, is sensibly lower than the others at, then regulation at that frequency must be accomplished through the openloop function, and the other branches must be inhibited by decreasing the gain ofas much as possible. The quantity of is a practical guideline to establish the bound height difference. It is funded on the ideal collaboration ofequal plants, which would reducetimes the amount of feedback, as Section 2.2 details. Figure 2 shows the analysis for a 2branch structure, which suggests that the controllershould carry the burden of control atwhatever the chosen is.(b)If several bounds are located at the bottom and their heights do not differ too much from each other at, then their associated controllers must share the regulating task. The branches whose bounds locate above this bound group must be inhibited. In this case, the practical rule of thumb is that any set ofbounds within a range of are suitable for cooperation. According to this, Figure 2 suggests that controllersandshould collaborate atwhatever the chosen is to fulfil the specification.
(a)
(b)
Remark 1. The first stage of the new methodology arranges the best load sharing amongst parallel plants at each frequency for minimum amount of feedback in each branch to meet the control specifications. This approach is completely new inside the MISO control literature.
2.4. Second Stage: Feedback Sharing in Parallel Structure of Controllers
Once it is decided which plants or plants must deal with the specifications at each frequency, theparallel controllers are designed sequentially. The procedure presented here aims to minimize the number of iterations required to reach a feasible solution. At each step, one of the controllers is tuned by QFT loopshaping, while the remaining ones are considered constant components of the overall uncertain system. This rotating nature of the tunable controller is enabled by the genbnds command of the Terasoft QFT Toolbox [59], which admits general specifications in the form
Thus, if at some point the controlleris to be adjusted, the robust stability (4) bounds can be computed by choosing, , , , and , while the robust disturbance rejection (5) bounds require the choice, , , , , and . Note that these are common QFT bounds, since they are imposed on the openloop transfer function, and consequently denoted, which represent all the specifications to meet at.
The design procedure begins with all controllers set to zero. Subsequently, in the first step of design, the loopshaping of the branch/branches that should work at the lowest frequency is performed fully. In an increasing order of frequencies, the last step corresponds with the design of branch/branches working at the highest frequency. During any loopshaping design, the course of action at certain frequency depends on the result of the load sharing planning stage. As common in QFT for each loopshaping, the design begins satisfying the bounds in an increasing order of frequencies, and finally, a fast rolloff of the openloop gain is shaped in the high frequency range. A novelty of MISO control is that the openloop gain must also be cut or reduced in the frequencies where the branch should be noncollaborating. Following, there are some hints of guide for both goals of design: when the branch must achieve the performance atand when the branch must be inhibited atbecause other branches do the work or because is the rolloff frequency for the branch.
At the frequencies where the branch/branches must perform the regulation task, two different cases are detailed next. Figures 3 and 4 show examples of both cases. Superscripts in bound designation indicate the step in the sequential design.(a)A single plantwill be responsible for the regulating task at(see the twobranch example in Figure 3). Thus,is placed above the bound(seeand). After this, thebound arrangement does not change (see). However, since the accomplishment of the specifications has just been secured by, the boundsdelimit now closed forbidden regions, which will allow a reduction of thegains for feedback saving atwhen these designs are faced (comparewith). Note also that, an excess of gain to meetwould imply a major contraction ofif this was necessary in the design of the other loops.(b)A group ofplants will share the control task at(see the fourbranch example in Figure 4, where). Their loop shaping is carried out jointly as follows. For convenience, let us denote the group. Let us assume as well thatis tuned first. The purpose of the designer is to placearound decibels below the boundat the design phase of interest(see and ). The consequences are thatbound does not change (see) and that a dip appears atin bounds(seeand). Then, tuning oftakes place, and the goal is to placearound,, decibels below the bound; (see and ). Once achieved, thebound does not change (see), but it extends the depth of the dip in(seeand). The process goes on in this manner, temporarily violating their bounds, until the tuning ofis reached. In this case, the purpose is simply to meet the bound, whose dip has been repeatedly enlarged by each one of the controllers(see and ). Once this step is completed, all previously violated bounds are now met (see, and). If theplants were identical atand, the gain of each controllerwould be lower. Nevertheless, differences between plant gains up to are acceptable for successful cooperation and controller gain reduction, as Section 2.3(b) reports. In consequence, some extra adjustments are usually necessary to meet each. After the design of thecontrollers with control responsibilities at, the bounds of the remaining loops will delimit closed forbidden regions (see).
(a)
(b)
At the frequencies where the branch/branches has/have no regulation responsibilities, a great deal of different situations may arise, specially when there are more than two inputs in the MISO system. Some of the most frequent ones are being described. In any case, the bound recalculation after any design step contributes to value the tradeoffs, which are inherent to any design with multiple degrees of freedom.(i)Design hints for the rolloff frequenciesof theloop. Thegain reduction is usually a tradeoff between a maximum noise amplification and a minimum controller order. One case is that the boundsdemand certain feedback at. This happens if there are other loops with control responsibilities at that have not been designed yet. In this circumstance, thedesign is performed temporarily violating its bounds at the rolloff frequencies. After the design of the loops with control responsibilities at , the boundsare automatically reshaped and now demarcate forbidden closed regions. In the new picture, will usually meet its bounds. Nevertheless, further iterations may be required to optimize the gain reduction at these rolloff frequencies. Other case is if theloop is the last one to be designed. In this circumstance,bounds, which already delimit forbidden closed regions, are the definitive restrictions to be met in thegain reduction at its rolloff frequencies.(ii)Design hints for the frequencieswhere theplant participation must be inhibited. The boundscurrently demarcate fobidden closed regions at, since other loops have already achieved the prescribed specifications. The restrictions for thegain reduction at frequencies are the bounds at frequencies higher than, where theplant will have to achieve or collaborate in the regulation task. The controller complexity is another reason to take into account in the gain saving procedure.
Remark 2. The second stage of quantitative MISO control details the QFT framework and procedure for the sequential loopshaping of parallel controllers in order to meet the predefined loadsharing arrangement. A set of branches collaborate in the control tasks or are inhibited for feedback saving. A standard method is given to deal withinput singleoutput systems from a quantitative feedback control point of view, which remained unsolved till now.
3. Design Examples
3.1. Example 1
A simple MISO system illustrates the new methodology and its potential benefits. The uncertaintyfree plant models are
Inside a parallel disposition (Figure 1(b) with), the regulatory problem aims to designto meet certain stability (4) and performance (5) specifications with minimum controller gains at each frequency. For a minimum stability phase margin of45°, the tolerance choice is And for disturbance rejection, the performance upper model is
For bound computation and loopshaping guide, a suitable vector of discrete frequencies is chosen:
3.1.1. First Stage: Load Sharing Planning at Each Frequency
According to explanations in Section 2.3, the required amount of feedback at is computed in case that a single system regulated the output:or, in Figure 5(a) and Figure 5(b), respectively (only three frequencies inare drawn to illustrate the three potential outcomes). Compare controller bounds of both systems as follow.(i)At low frequencies (), the boundsare over lower than the bounds; see. Thus,branch should assumeregulation andbranch should be switched off. Subsequent loopshaping goals are ,at.(ii)At medium frequencies (), the bound heights differ less than20 log 2; see. Thus, the collaboration of both plants is suggested. Subsequent loopshaping goals are at.(iii)At high frequencies (), the boundsare over lower than the bounds; see. Subsequent loopshaping goals are ,at.
(a)
(b)
Accordingly, Table 1 summarises the load sharing planning at design frequencies.

3.1.2. Second Stage: Design of Parallel Controllers
Sequential loopshaping ofandis performed with Table 1 aims. Figure 6 illustrates the procedure in detailed steps. Steps 0, I, and II belong to the same first iteration, which explains the design ofand (marked with superscripts). Further iterations (step III) are required to optimize both designs looking for the strictly minimum amount of feedback at each frequency, which usually trades off with a reasonable order of the controllers. General guidelines were described in Section 2.4.(0)Initial bound computation is drawn withfor the loop 1 (Figure 6(a)) and withfor the loop 2 (Figure 6(b)).(I)Loopshaping of(see Figure 6(c)). To achieveat,is located onto bound. For plant collaboration at,is shaped belowand. At, the goal becomes. Then,are the rolloff frequencies for. Its gain is reduced as much as possible, even momentary violating bounds, and, since control requirements will be achieved by. As this reduction is not bound guided, it could require later areshaping at these high frequencies. After step I, the controllers are As a consequence, suitable boundsarise for the subsequent loopshaping of. Let us compare new bounds in Figure 6(d) with initial bounds in Figure 6(b): closed exclusion regions appear inforfeedback saving at, sincehas already achieved the feedback requirements; dips arise inand, wherehas to collaborate with previouscontribution; the bounds at are not modified since must assume the whole control tasks.(II)Loopshaping of(see Figure 6(f)). Firstly, reducegain atto meet sincewas achieved. This example does not need adding a lag network, which favours a lower order of controller. Secondly, locateonto its boundsat; note that the purpose of an excess of feedback at is the illustration of subsequent optimization iterations. Finally, reduce the high frequency gain ofat(rolloff frequencies) as much as possible as in any classical QFT design. After step II, the controllers are As a consequence, let us compare new boundsin Figure 6(e) with bounds of the previous step in Figure 6(c). Note thatnow meets all its bounds. Pay special attention to collaboration frequencies in Figures 6(e) and 6(f): the feedback sharing is close to the optimum at(andlie on their respective bounds), but it is not the case at. Thus, small adjustments are required to economize feedback.(III)Both loops are iteratively redesigned in Figures 6(g) and 6(h). Thanks to the help of software tools, this is usually a simple procedure. The final controllers are Note thatachieves by itself the control requirements at and shares this task withat. Thus,lies onat. Alongside, achieves by itself the performance load at. Thus,lies onat. It’s also necessary to remark how the bound disposition eases the reduction of the openloop gains, in particular at for and at , for . Nevertheless, the feedback saving at these frequencies trades off with a minimum order of controllers.
3.1.3. Analysis and Comparisons
Parallel control,(13) and (14), is being compared with master control,, which is designed for the equivalent plantin order to satisfy the same performance and stability specifications. QFT design yields the single controller
For a fairer comparison,lies exactly on its bounds and reduces its high frequency gain as much as possible. Thus, it is the solution of minimum amount of feedback at each frequency if the system was single input; Figure 7(a) depicts the minimum. This QFT optimal loopshaping is easily achievable as there is no uncertainty and there is a single controller. For the MISO solution, the exact meeting of bounds with controllers (13) and (14) required a second iteration to reshape the former controllers (11) and (12). Nevertheless, the main feedback saving is not so much to lie exactly on the bounds but to share conveniently the feedback between branches. Note thatat low frequencies, whereregulates the output, and at high frequencies, whereworks. In this way,filters low frequencies, which avoids unnecessary offsets ofactuation. Andfilters high frequencies, which reduces the amplification of high frequency signals (sensor noise) inactuator (usually labelled as cost of feedback in SISO QFT [58]).
(a)
(b)
Figure 8 clearly demonstrates these benefits in the time domain. External inputs are a unity step disturbance at s and the sensor noisebuilt with a whitenoise source (noise power = 0.001, sample time = 0.01). All response signals are depicted with and without noise intervention in blue and green lines, respectively. As long as(a better agreement could be achieved with higher order of controllers), the closedloop transfer functionsandare similar for both control strategies. Accordingly, theperformance coincides for master and parallel control (only one is depicted in Figure 8), in presence (blue line) or absence (green line) of noise. However, the control effort strikes the differences as setting forward by closedloop frequency responses in Figure 7(b). Let us compare the pairwith the signalin Figure 8. Actuationis hardly affected by sensor noise amplifications, which avoids unnecessary fatigue ofactuator; the exact quantification of this benefit is on the distance betweenand, which mainly depends on the affordable bandwidth saving during theshaping. On the other hand, the noise affects similarly to both and . However,recovers the initial stationary state anddoes not: see,in timedomain, and,in the frequency domain. Thus, for accumulative disturbances, the fast time actuation ofstarts out always from zero, and consequently, the full range ofsignal is available to speed up the transient response with a less saturation risk. Both improvements (the noise release ofand the zero steady state of) are more valuable in practice, where reallife plants, actuators, and sensors have physical limitations.
As some final comparisons, onceandare achieved, there are equivalences in other MISO control topologies (Eitelberg’s masterslaves [46] and Lurie’s serial structure [39]) for this particular example. Designs according to those other structures and methodologies can always be transformed into the strictly parallel structure, while the inverse transformation is not always feasible. In any case, the same QFT optimal solutions (13) and (14) would not probably be attainable since those other design methodologies do not quantify neither the amount of feedback nor the collaboration degree. Besides, some MISO plants or designcondition arrangements are only affordable with the new parallel design methodology; as following examples are discussing.
Remark 3. Firstly, this example details the two stages of the proposed MISO control methodology. Then analysis and comparisons show the expected benefits of a parallel structure of controllers in comparison with a single master controller: the reduction of useless fatigue and saturation risk of actuators. In comparison with other MISO control strategies (masterslaves or serial structures), the new one offers the minimum amount of feedback to govern each plant for prescribed specifications.
3.2. Example 2
Heat exchangers are common engineering applications of MISO control [18–20]. Figure 9 outlines a heatexchanger where the output temperature ,can be regulated through two manipulated variables: the hot flowand the heating power control(which conditions the input temperature of the hot flow). The cold flowacts as a disturbance since the production rate changes on demand. A laboratoryscale plant commercialised by Feedback Systems Inc. agrees with this arrangement [60].
Around certain operating conditions, linear plant models of the form are identified [61]. They express the relationships from the two manipulated variables to the output (and) and from the disturbance to the output (). A total amount of 112 cases are defined for each output/input relation, which generate the uncertainty of numerator and denominator coefficients in (16). Figure 10 depicts the envelopes of the plant frequency responses. As long as all themodels are minimum phase, only magnitude Bode diagrams are included.
The robust stability control specification is (4) with tolerance (8) for a minimum phase margin of45°. The specification for robust disturbance rejection is (5) with upper tolerance model
After computing the demand of feedback iforworked alone, Table 2 resumes the loadsharing planning at discrete frequencies to meet the robust control specifications with the minimum amount of feedback.

Consequently, the loopshaping of and are performed (see Figure 11) and yield
(a)
(b)
In Figure 12, closedloop frequency responses illustrate the achievement of control specifications and the control effort that is used for it. Only the envelopes of a set of 112 frequency responses are depicted for each plant. In summary, the originality of this example relapses in the engineering character of the MISO plant, the robustness of the new methodology (model uncertainty consideration), and the challenge of disturbance incorporation through a dynamical model with uncertainty.
(a)
(b)
(c)
(d)
For comparison purposes, a master controller is firstly designed for the summation of plantsin order to meet the same robust control specifications
The parallel arrangement and the master controller behave equally attending to the prescribed stability and performance. However, the parallel one is superior in terms of the actuation signals. The same reasoning and benefits of example 1 can be extended to this example.
Secondly, comparisons of the parallel arrangement and the masterslaves structures are shown. The masterslaves structure (Figure 1(b) with,,) can be considered a particular case of the parallel structure (), since this last one will be demonstrated more versatile. Attending their functionality as filters, let us rename slaveas filterand slaveas filter. Let us suppose that(18) and(19) have already been designed for a parallel control structure. Then, there would not be a single equivalence for them in the masterslaves structure. Equivalents should meet:and. A practical solution would be firstly to assume the masterof (20), as long as its output performance matches that of the set(18) and(19). And secondly, to overcome the drawbacks of master control for parallel plants, filters are computed asandwith(18) and (19). This equivalence procedure is also reliable in example 1. The inverse procedure, that is, the direct design of,andcan be successfully done in example 1; however, it would become extremely difficult in example 2 for the following reasons. The masterslaves method firstly addresses the design of the filters and then the master. As a first difficulty, the filters have to perform the load sharing in frequency terms, but their manual openloop design would not meet quantitative closedloop criteria. Besides, the uncertainty in parallel plant models (,) hinders the choice of the cutoff frequency and the rolloff of the filters. Another difficulty is that the disturbance participation through dynamical models () makes impossible to guess the best load sharing along the frequencies between the parallel plants (,), and then it hinders a suitable parametrization of filters. At a second step of the masterslave method of design, the master would be designed, taking into account the initial MISO plant and the branch filters. In this stage, filter rolloffs may add serious problems: a high rolloff improves the filtering power but the branch phase deviates far from the original plant phase. Thus, branches can reach counterphase conditions at certain frequencies. In such a case, a master controller spoils plant collaborations, since the vectorial sumyields a smaller modulus thanand .
Remark 4. A common engineering problem is solved: the robust regulation of outputflow temperature in a heat exchanger manipulating the inputflow heating power and the input heating flow itself, for different production rates (outputflows). Model uncertainty is fully considered to obtain the minimum amount of feedback controllers, whose benefits are detailed in previous examples. Additionally, the difficulties exhibited by the masterslaves strategy are overcome. The challenge is due to plant uncertainty and to the dynamical incorporation of disturbances.
3.3. Example 3
Despite the large number of engineering systems that have more than two manipulated inputs and single regulated output, few applications of this kind have been solved [40] so far. In fact, only a few scientific works incorporate a control theory forinputs [25, 26, 45, 48, 56] because of its complexity. In this sense, the following example tries to demonstrate the potential of the new methodology. Consider a system with three manipulated inputs and a single regulated output, whose transfer function models are which are particularised for eachinput by,,,,,,,,. Disturbances come into the single output straightforward,. As the difficulty of this example relapses on dealing with three actuators, parameter uncertainty has been intentionally omitted in the plant models. Stability and disturbance rejection have to be attended meeting (4), (8) and (5), (9), respectively.
Upper section of Table 3 shows the amount of feedbackthat a single branchwould demand at the design phaseto meet the specifications by itself. Note that the selection of design phase is part of the designer tasks. Previous examples depicted control necessities on the Nichols chart, that is, for all possible design phases. Consequently with the upper part, the bottom of Table 3 shows the planning on load sharing along the frequencies; revise the criteria in Section 2.3 to decide plant collaborations or not as a function of feedback demand.

Accordingly to results and choices on Table 3, the loopshaping of each,is sequentially performed (see the final arrangements in Figure 13). The achieved controllers are
(a)
(b)
(c)
Figure 14 depicts the expected feedback sharing along the frequency band (Table 3) in terms of branch openloop functions and their contribution to the global openloop functionin dB units.
The originality of this theoretical example relapses mainly in two challenges: facing regulatory problems with more than two manipulated inputs and dealing with parallel plants when some of them must work in disjoint bands of frequencies. Note thatjoins first atand later on at. This causes that there is not an equivalent for(22),(23), and (24) in a serial structure (Figure 1(a)). Let us rename serial controllers as,, and. The equivalents can be computed as , and ,, which would yield a nonproper controller. On the other hand, it is not feasible either a straight away design of the serial controllers using a mainvernier strategy, since the load sharing planning (Table 3) does not allow sorting the plants from faster to slower (joins at discontinuous intervals of frequencies). That is, the serial strategy is conceived for that each plant works in one and only one frequency interval.
Remark 5. This example tries to illustrate the new methodology ability (i) to solve regulatory problems of more than two manipulated inputs, which are frequent in engineering applications but are not so common in the scientific literature and (ii) to deal with parallel plants when some of them must work in disjoint bands of frequencies. For this last case, a serial structure of controllers (mainvernier) fails.
4. Conclusions
In attention to engineering application demands, a parallel structure of controllers has been proposed to solve MISO regulatory problems, which concern robust stability and robust disturbance rejection. MISO plant models included parametric uncertainty. Inside a QFT framework, a new methodology of design was proposed with the aim of tuning a set of controllers that met the robust specifications, and simultaneously, each controller used the least possible gain at each frequency. The keypoint was to quantify the feedback demanded by each plant at each frequency to meet the specifications by itself. This allowed a planning of the best load sharing along the frequency band. At each frequency, either one or several plants in collaboration must assume feedback responsibilities to meet the specifications or the plant/s participation must be inhibited for feedback saving. To accomplish it, a sequential design of the parallel controllers was fully detailed, using current QFT software tools for bound computations and for loopshapings.
The abilities and engineering interest of the new methodology were demonstrated through several examples. The first one illustrated the methodology and detailed the benefits of a parallel structure of robust controllers in comparison with a single robust master controller. The main advantages are (i) to reduce the sensor noise amplification in the actuators of the slowest plants, which are conceived to work in the band of lower frequencies and provide the stationary and (ii) to avoid unnecessary offsets in the actuators of the faster plants, and consequently, that their full range of actuation is available to speed up the response with less risk of saturation. Nevertheless, these are also the aims of current MISO methodologies, whose weak points have been also overcome by the new technique. In particular, a second engineering example exhibited the difficulties of tuning a parallel control using the masterslave philosophy. The challenge is mainly the filter (slaves) design due to plant uncertainties and the dynamic incorporation of disturbances, which are common in reallife plants. A third example illustrated the weak points of a chained serial structure of controllers (mainvernier), which fails when plants must work in disjoint bands of frequencies. In any case, the proposed parallel structure can always be converted to both masterslaves or mainvernier structures. The inverse procedure was demonstrated as not always being possible. In addition, the third example showed the difficulty of facing a MISO problem with three manipulabled variables, which are not common in the scientific literature but are frequent in engineering applications.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
The authors are grateful for the assistance provided by La Rioja Government through Project IMPULSA 2010/01 and through ADER Grant 2012IIDD00093.
References
 M. N. Marwali, J.W. Jung, and A. Keyhani, “Control of distributed generation systems—part II: load sharing control,” IEEE Transactions on Power Electronics, vol. 19, no. 6, pp. 1551–1561, 2004. View at: Publisher Site  Google Scholar
 E. Eitelberg, “Macrodynamic feedback interaction between trade and production,” International Journal of Robust and Nonlinear Control, vol. 17, no. 23, pp. 203–224, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. GarciaSanz and F. Y. Hadaegh, “Loadsharing robust control of spacecraft formations: deep space and low Earth elliptic orbits,” IET Control Theory and Applications, vol. 1, no. 2, pp. 475–484, 2007. View at: Publisher Site  Google Scholar
 I. Horowitz, Synthesis of Feedback Systems, 1963.
 M. Kano and M. Ogawa, “The state of the art in chemical process control in Japan: good practice and questionnaire survey,” Journal of Process Control, vol. 20, no. 9, pp. 969–982, 2010. View at: Publisher Site  Google Scholar
 S. K. Jha and N. Kaistha, “Valve positioning control for process throughput maximization,” Chemical Engineering Research and Design, vol. 85, no. 11, pp. 1465–1475, 2007. View at: Publisher Site  Google Scholar
 B. J. Allison and A. J. Isaksson, “Design and performance of midranging controllers,” Journal of Process Control, vol. 8, no. 56, pp. 469–474, 1998. View at: Google Scholar
 B. J. Allison and S. Ogawa, “Design and tuning of valve position controllers with industrial applications,” Transactions of the Institute of Measurement and Control, vol. 25, no. 1, pp. 3–16, 2003. View at: Publisher Site  Google Scholar
 C.C. Yu and W. L. Luyben, “Analysis of valveposition control for dualinput processes,” Industrial and Engineering Chemistry Fundamentals, vol. 25, no. 3, pp. 344–350, 1986. View at: Google Scholar
 R. B. McLain, M. J. Kurtz, M. A. Henson, and F. J. Doyle III, “Habituating control for nonsquare nonlinear processes,” Industrial and Engineering Chemistry Research, vol. 35, no. 11, pp. 4067–4077, 1996. View at: Google Scholar
 R. MonroyLoperena, R. Solar, and J. AlvarezRamirez, “Balanced control scheme for reactor/separator processes with material recycle,” Industrial and Engineering Chemistry Research, vol. 43, no. 8, pp. 1853–1862, 2004. View at: Google Scholar
 A. VelascoPerez, J. AlvarezRamirez, A. Velasco, and R. SolarGonzalez, “Control múltiple entrada una salida (miso) de un cstr,” Revista Mexicana De Ingeniería Química, vol. 10, no. 2, pp. 321–331, 2011. View at: Google Scholar
 O. A. PradoRubio, S. B. Jorgensen, and G. Jonsson, “Ph control structure design for a periodically operated membrane separation process,” Computers and Chemical Engineering, vol. 43, pp. 120–129, 2012. View at: Google Scholar
 S. Haugwitz, P. Hagander, and T. Norén, “Modeling and control of a novel heat exchange reactor, the open plate reactor,” Control Engineering Practice, vol. 15, no. 7, pp. 779–792, 2007. View at: Publisher Site  Google Scholar
 S. Velut, L. de Maré, and P. Hagander, “Bioreactor control using a probing feeding strategy and midranging control,” Control Engineering Practice, vol. 15, no. 2, pp. 135–147, 2007. View at: Publisher Site  Google Scholar
 M. Karlsson, O. Slätteke, B. Wittenmark, and S. Stenström, “Reducing moisture transients in the papermachine drying section with the midranging control technique,” Nordic Pulp and Paper Research Journal, vol. 20, no. 2, pp. 150–155, 2005. View at: Google Scholar
 A. VelascoPérez and J. AlvarezRamírez, “Algorithm for parallel control for an aerobic reactor with a recirculation current,” Revista Mexicana De Ingeniería Química, vol. 6, no. 2, pp. 229–236, 2007. View at: Google Scholar
 M. A. Henson, B. A. Ogunnaike, and J. S. Schwaber, “Habituating control strategies for process control,” AIChE Journal, vol. 41, no. 3, pp. 604–618, 1995. View at: Google Scholar
 D. F. Marselle, M. Morari, and D. F. Rudd, “Design of resilient processing plantsII design and control of energy management systems,” Chemical Engineering Science, vol. 37, no. 2, pp. 259–270, 1982. View at: Google Scholar
 L. Giovanini, “Cooperativefeedback control,” ISA Transactions, vol. 46, no. 3, pp. 289–302, 2007. View at: Publisher Site  Google Scholar
 G. F. Fauth and F. G. Shinskey, “Advanced control of distillation columns,” Chemical Engineering Progress, vol. 71, no. 6, pp. 49–54, 1975. View at: Google Scholar
 F. G. Shinskey, “Energyconserving control systems for distillation units,” Chemical Engineering Progress, vol. 72, no. 5, pp. 73–78, 1976. View at: Google Scholar
 F. G. Shinskey, “Control systems can save energy,” Chemical Engineering Progress, vol. 74, no. 5, pp. 43–46, 1978. View at: Google Scholar
 S. J. Williams, D. Hrovat, C. Davey, D. Maclay, J. W. V. Crevel, and L. F. Chen, “Idle speed control design using an Hinfinity approach,” in Proceedings of the American Control Conference, pp. 1950–1956, June 1989. View at: Google Scholar
 S. Jayasuriya and M. A. Franchek, “A qfttype design methodology for a parallel plant structure and its application in idle speed control,” International Journal of Control, vol. 60, no. 5, pp. 653–670, 1994. View at: Google Scholar
 V. Besson and A. T. Shenton, “Interactive parameter space design for robust performance of MISO control systems,” IEEE Transactions on Automatic Control, vol. 45, no. 10, pp. 1917–1924, 2000. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 P. Gorzelic, E. Hellstrom, A. Stefanopoulou, L. Jiang, and S. Gopinath, “A coordinated approach for throttle and wastegate control in turbocharged spark ignition engines,” in Proceedings of the 24th Control and Decision Conference (CCDC '12), pp. 1524–1529, 2012. View at: Google Scholar
 N. Ravi, H.H. Liao, A. F. Jungkunz, and J. C. Gerdes, “Midranging control of a multicylinder HCCI engine using split fuel injection and valve timings,” in Proceedings of the 6th IFAC Symposium Advances in Automotive Control (AAC '10), pp. 797–802, July 2010. View at: Publisher Site  Google Scholar
 S. Jade, E. Hellstrom, L. Jiang, and A. G. Stefanopoulou, “Fuel governor augmented control of recompression hcci combustion during large load transients,” in Proceedings of the American Control Conference, pp. 2084–2089, 2012. View at: Google Scholar
 S. Jade, J. Larimore, E. Hellstrom, L. Jiang, and A. G. Stefanopoulou, “Enabling large load transitions on multicylinder recompression hcci engines using fuel governors,” in Proceedings of the American Control Conference, pp. 4423–4428, 2013. View at: Google Scholar
 J. Zhai, Y. Huang, S. Schroeck, W. Messner, D. D. Stancil, and T. E. Schlesinger, “High bandwidth electrooptic scanner for optical data storage,” Japanese Journal of Applied Physics, vol. 39, no. 2, pp. 883–887, 2000. View at: Google Scholar
 A. Al Mamun, I. Mareels, T. H. Lee, and A. Tay, “Dual stage actuator control in hard disk drive—a review,” in Proceedings of the 29th Annual Conference of the IEEE Industrial Electronics Society, pp. 2132–2137, November 2003. View at: Publisher Site  Google Scholar
 T. Gao, C. Du, C. P. Tan, Z. He, J. Yang, and L. Xie, “High bandwidth control design and implementation for a dualstage actuation system with a microthermal actuator,” IEEE Transactions on Magnetics, vol. 49, no. 3, pp. 1082–1087, 2013. View at: Google Scholar
 J. Sun, F. Cameron, and B. W. Bequette, “A habituating blood glucose control strategy for the critically ill,” Journal of Process Control, vol. 22, no. 8, pp. 1411–1421. View at: Google Scholar
 K. Aström and T. Hägglund, Control Pid Avanzado, Pearson Prentice Hall, 2009.
 V. Lersbamrungsuk, T. Srinophakun, S. Narasimhan, and S. Skogestad, “Control structure design for optimal operation of heat exchanger networks,” AIChE Journal, vol. 54, no. 1, pp. 150–162, 2008. View at: Publisher Site  Google Scholar
 J. AlvarezRamirez and H. Puebla, “On classical PI control of chemical reactors,” Chemical Engineering Science, vol. 56, no. 6, pp. 2111–2121, 2001. View at: Publisher Site  Google Scholar
 J. AlvarezRamirez, A. Velasco, and G. FernandezAnaya, “A note on the stability of habituating process control,” Journal of Process Control, vol. 14, no. 8, pp. 939–945, 2004. View at: Publisher Site  Google Scholar
 B. J. Lurie and P. Enright, Classical Feedback Control: With MATLAB, Control Engineering, Taylor & Francis, 2000.
 R. L. Grogan, G. H. Blackwood, and R. J. Calvet, “Optical delay line nanometer level pathlength control law design for spacebased interferometry,” in Proceedings of the International Society for Optical Engineering (SPIE '98), pp. 14–25, March 1998. View at: Publisher Site  Google Scholar
 G. W. Neat and A. Abramovici, “Control technology lessons learned: case study using the microprecision interferometer testbed,” in Proceedings of the American Control Conference, vol. 2, 1998. View at: Google Scholar
 J. Daeges, B. Lurie, and A. Bhanji, “An improved highvoltage dc regulator for a radar and communication transmitter,” in Proceedings of the 18th IEEE Conference Record of the Power Modulator Symposium, pp. 106–108, 1988. View at: Google Scholar
 Q.G. Wang, Y. Zhang, W.J. Cai, Q. Bi, and C.C. Hang, “Cooperative control of multiinput singleoutput processes: online strategy for releasing input saturation,” Control Engineering Practice, vol. 9, no. 5, pp. 491–500, 2001. View at: Publisher Site  Google Scholar
 L. L. Giovanini, “Flexiblestructure control: a strategy for releasing input constraints,” ISA Transactions, vol. 43, no. 3, pp. 361–376, 2004. View at: Google Scholar
 E. Eitelberg, Load Sharing Control, NOYB Press, 1999.
 E. Eitelberg, “Load sharing in a multivariable temperature control system,” Control Engineering Practice, vol. 7, no. 11, pp. 1369–1377, 1999. View at: Publisher Site  Google Scholar
 E. Eitelberg, “Some peculiarities of load sharing control,” International Journal of Robust and Nonlinear Control, vol. 13, no. 7, pp. 607–618, 2003. View at: Publisher Site  Google Scholar
 S. J. Schroeck, W. C. Messner, and R. J. McNab, “On compensator design for linear timeinvariant dualinput singleoutput systems,” IEEE/ASME Transactions on Mechatronics, vol. 6, no. 1, pp. 50–57, 2001. View at: Publisher Site  Google Scholar
 S. Brennan and A. Alleyne, “Integrated vehicle control via coordinated steering and wheel torque inputs,” in Proceedings of the American Control Conference, pp. 7–12, June 2001. View at: Google Scholar
 B. J. Allison and J. B. Ball, “Constrained model predictive control of blow tank consistency,” Control Engineering Practice, vol. 12, no. 7, pp. 837–845, 2004. View at: Publisher Site  Google Scholar
 S. Haugwitz, M. Karlsson, S. Velut, and P. Hagander, “Antiwindup in midranging control,” in Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference (CDCECC '05), pp. 7570–7575, December 2005. View at: Publisher Site  Google Scholar
 M. Karlsson, O. Slatteke, T. Hagglund, and S. Stenstrom, “Feedforward control in the paper machine drying section,” in Proceedings of the American Control Conference, 2006. View at: Google Scholar
 S. Gayadeen and W. Heath, “An internal model control approach to midranging control,” in Proceedings of the 7th IFAC International Symposium on Advanced Control of Chemical Processes (ADCHEM '09), pp. 542–547, July 2009. View at: Google Scholar
 W. P. Heath and S. Gayadeen, “Simple robustness measures for control of miso and simo plants,” in IFAC Proceedings, vol. 18, pp. 11356–11361, 2011. View at: Google Scholar
 I. Horowitz, Quantitative Feedback Design Theory (QFT), QFT Publications, Boulder, Colorado, USA, 1993.
 P.O. Gutman, E. Horesh, R. Guetta, and M. Borshchevsky, “Control of the aeroelectric power station—an exciting qft application for the 21st century,” International Journal of Robust and Nonlinear Control, vol. 13, no. 7, pp. 619–636, 2003. View at: Publisher Site  Google Scholar
 F. G. Shinskey and J. H. P. Castellanos, “Sistemas de control de procesos: aplicación, diseño y sintonización,” in Sistemas De Control De Procesos: Aplicación, Diseño Y Sintonización, McGrawHill, 1996. View at: Google Scholar
 A. Gera and I. Horowitz, “Optimization of the loop transfer function,” International Journal of Control, vol. 31, no. 2, pp. 389–398, 1980. View at: Google Scholar
 C. Borghesani, Y. Chait, and O. Yaniv, Quantitative Feedback Theory Toolbox. For Use With Matlab, Terasoft, 2nd edition, 2002.
 F. Instruments, Limited Proces control, http://www.fbk.com.
 J. Rico and M. GilMartinez, “Multivariable QFT robust control of a heat exchanger,” in Proceedings of the 19th Mediterranean Conference on Control and Automation (MED '11), pp. 588–593, June 2011. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Javier RicoAzagra et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.