Open Access is an initiative that aims to make scientific research freely available to all. To date our community has made over 100 million downloads. It’s based on principles of collaboration, unobstructed discovery, and, most importantly, scientific progression. As PhD students, we found it difficult to access the research we needed, so we decided to create a new Open Access publisher that levels the playing field for scientists across the world. How? By making research easy to access, and puts the academic needs of the researchers before the business interests of publishers.
We are a community of more than 103,000 authors and editors from 3,291 institutions spanning 160 countries, including Nobel Prize winners and some of the world’s most-cited researchers. Publishing on IntechOpen allows authors to earn citations and find new collaborators, meaning more people see your work not only from your own field of study, but from other related fields too.
When the manipulated variable (MV) has significantly large time delay in changing the control variable (CV), use of the currently measured CV in the feedback error can result in very deficient feedback control (FBC). However, control strategies that use forecast modeling to estimate future CV values and use them in the feedback error have the potential to control as well as a feedback controller with no MV deadtime using the measured value of CV. This work evaluates and compares FBC algorithms using discrete-time forecast modeling when MV has a large deadtime. When a feedforward control (FFC) law results in a physically unrealizable (PU) controller, the common approach is to use approximations to obtain a physically realizable feedforward controller. Using a discrete-time forecast modeling method, this work demonstrates an effective approach for PU FFC. The Smith Predictor is a popular control strategy when CV has measurement deadtime but not MV deadtime. The work demonstrates equivalency of this discrete-time forecast modeling approach to the Smith Predictor FBC approach. Thus, this work demonstrates effectiveness of the discrete-time forecast modeling approach for FBC with MV or DV deadtime and PU FFC.
Model Predictive Control
Nonlinear Dynamic Modeling
chapter and author info
Derrick K. Rollins*
Department of Chemical, Biological Engineering and Statistics, Iowa State University, Ames, Iowa, USA
*Address all correspondence to: email@example.com
From the Edited Volume
Model Predictive Control - Recent Design and Implementations for Varied Applications [Working Title]
Dr. Umar Zakir Abdul Hamid and Dr. Ahmad 'Athif Mohd Faudzi
Modeling data is critical to the advancement of information and data science on many levels and in many areas. Accurately modeling data is often important to system monitoring, understanding, and control; and thus, ultimately to the advancement of technology.
A characteristic of data that that is not well understood, even by those in the physical sciences, is dynamic behavior. However, the behavior of Covid-19, which is inherently dynamic, has forced wide-spread conversations from even the non-science community about such terms as lag and deadtime. Just as understanding the attributes of Covid-19 dynamically can lead to intelligent and thus, safe behavior, good decision making and preparation, and even save lives, the lack of understanding can do the opposite. Thus, the better our understanding of dynamic behavior is, physically and biologically, the better our understanding of data will be, leading only to better solutions to many problems facing society.
Forecast modeling is a type of predictive modeling that uses current and anticipated future input values to predict values for outputs in the future. For example, forecast modeling is used to predict the wind velocity in a certain region five days into the future. Another example is forecast modeling to predict the number of deaths caused by a virus a week into the future. This chapter focuses on the application of forecast modeling to enhance feedback control (FBC) and feedforward control (FFC) when the problem is physically unrealizable (PU).
A PU system is a mathematical phenomenon of a dynamic system. It occurs in two ways. The first one is when the order of the differential equation for the output is less than the order of the differential equation for the input. An example is the development of a FFC law determined from a load transfer function divided by the process transfer of a lower order. The second one is when an output depends on deadtime that has a negative value, in effect causing a dependence on future values of a time dependent variable(s). This also occurs in FFC when the load transfer function has a smaller deadtime than the process transfer function. The most common approach for addressing a PU system is to use approximation(s) to make it physically realizable. However, such approximations can lead to large modeling errors, thus leading to unacceptable control.
The dynamic modeling literature [1, 2] defines causalitysomewhat differently than the statistics literature. More specifically, “if a system[‘s] output depends on the future input values … the system is noncausal .” This definition is synonymous with PU, it seems. In forecast modeling, all values of inputs are before the forecast distance in the future. They can be in the future, but not a distance beyond the forecast time. Another description for PU in the dynamic modeling literature is improper transfer function.
In FFC, MV is the output and depends on input changes. When MV has a deadtime of θMV, for example, it takes this time before a change in MV affects CV. Within this period, other inputs may change that impact CV before the change in MV does. An example is the control of blood glucose concentration (BGC) in type 1 diabetes. The deadtime for insulin infusion is much larger than the deadtime for carbohydrate consumption. The best approach for control of BGC is to consider the timing and the amount of carbohydrates consumed and to bolus this with a determined amount of insulin, a calculated amount of the time before the meal. This procedure is just a manual type of PU FFC practiced by people with type 1 diabetes.
The use of causality in the statistics literature seeks to distinguish it from correlation. Thus, in the statistics literature, causality is not focused exclusively on dynamic systems (e.g., only those with lag or deadtime) but a cause-and-effectrelationship between input and output, that can be nondynamic . For forecast modeling, cause-and-effectis not essential if the model is accurate. However, in control, cause-and-effectmodeling is essential. A PU system does not have an exact solution which would be a continuous-time solution. However, a discrete-time solution can be determined directly from the PU continuous-time structure. Hence, this work uses highly structured discrete-time forecasting and FFC models.
1.1 Objective and contribution of the work
Moreover, the primary objective of this chapter is to apply a novel discrete-time forecasting modeling methodology to systems with large MV deadtime in FBC and PU FFC without physical realizable approximations. For this scope, FBC is examined and evaluated under three prediction horizons: 1. None -- Classical FBC that uses the currently measured value of CV [4, 5]; 2. θMV -- Feedback Predictive Control (FBC)  and; greater than θMV – Model Predictive Control (MPC) [7, 8]. The Smith Predictor (SP)  is a novel FBC approach when θMV = 0 and there is deadtime in the measurement of CV. This work shows that FBPC gives equivalent control of the SP, but also has the advantage that it is applicable when θMV > 0, which the SP is not. In addition, this work reveals the detrimental use of the bias correction as given in the block diagram of the SP and used widely in process control . Thus, this work proposes a better bias correction method. Finally, this works presents a novel discrete-time PU FFC algorithm that is multiple-input and single-output and hence, is able to treat complex multiple-input feedforward model structures. Although MPC is a FBC approach, comparison is made to illustrate the potential improvement of FFPC over model-based predictive FBC.
Physical unrealizability (PU) is an anomaly that is strictly an artifact of a dynamic system. A dynamic system has at least one process state (i.e., output or response) that does not change to its new value immediately when input changes occur that cause its value (i.e., level) to change. This behavior contrasts with a nondynamic system that changes to its new state immediately when inputs change (also, called “disturbances”).
There are two basic dynamic phenomena – time lag and time delay (also called “dead time”) which are shown in Figure 1. This figure illustrates nondynamic and dynamic relationships for the response, y, to a step change in the input, x, occurring at time (t) = 0. As shown, for the nondynamic response (a), ychanges immediately to its new steady state value. Lag (b) is shown by the change in ystarting to occur at t = 0, but monotonically increasing over time to its new steady state value. Time delay (c) is shown by the change in yto its new value occurring θtime later. Lag and time delay (d) are shown by ystarting to change θtime later and then monotonically increasing over time to its new steady state value. “Everyday” examples of nondynamic changes are eyes opening (x) and immediate sight (y) and turning on the radio in a car (x) and hearing it (y) immediately. A dynamic time delay example is lightening occurring very far away. It occurs when one sees the lightening, but the thunder is delayed and occurs at a significant time after seeing the lightening. It does not build up to its final value, there is just a big boom that occurs, essentially, at once. A dynamic change with lag occurs when a person has been out in the cold for a while and their skin temperature is quite cold and when they move to a warmer environment, their temperature starts to rise but it takes time for the temperature to reach its new level in this warmer environment.
When a system is dynamic, its mass and/or internal energy changes over time, being driven to a new state due to input changes, arriving there at a time different than when the input was changed. Mathematically (and theoretically) this is seen as a time-order differential equation. Such an equation is given in terms of input xand output yin Eq. (1).
The dynamic form represented by Eq. (1) is PU if n < m, if θis negative, or if both are true. More specifically, the output, y, which depends on the input, x, cannot have a time dependent derivative structure that is of lower order than the variable that causes it to change. The response of a system to a disturbance also cannot have negative time delay. A system cannot respond to a disturbance before it occurs. There is no true solution for these conditions.
However, there are ways to address these PU cases in practice. For the first case, n < m, discrete-time backwards different derivatives can be used to approximate the continuous-time derivatives. This approach should provide adequate accuracy when the sampling time is constant and sufficiently small, and sensor noise is not too great. Eq. (2) illustrates this approximation when m = 2, n = 1, a constant sampling time, Δt, and with θ = 0 (for simplicity). Note, since y(t) cannot be immediately affected by x(t), xt-Δtis used to approximate x(t).
Thus, digital and sensor technologies, among other advancements, have significantly contributed to an ability to approximate Eq. (1) when n < m.
For the other case, i.e., when the time delay is a negative value such as a -5Δt(e.g., in FFC when deadtime for the disturbance variable (θDV) is smaller than the deadtime for the manipulated variable (θMV)), Eq. (2) becomes
Thus, as shown by Eq. (3), the output depends on future values of the input. However, discrete-time modeling provides a means to express an approximate solution in a PU form where knowledge of future changes allows approximation of the output. An example where this type of approximation is applied,is the control of blood glucose concentration (BGC) for people with type 1 diabetes . The manipulated variable (MV) for the automatic regulation of the exogenous insulin infusion from a servo-mechanical pump can have a deadtime of 60 minutes and carbohydrates from meals can have a deadtime of 30 minutes , resulting in x(t) becoming x(t + 30) in the numerator of the FFC law, with MVas the output variable. Moreover, for these values, a change in insulin flow rate will take one hour to begin to lower BGC. During this period, eating can increase BGC. People with type 1 diabetes understand this phenomenon and will bolus their insulin infusion, based on when they will eat and how many carbohydrates they will eat. This is called “a meal announcement” [12, 13]. But just as this idea is applied to carbohydrates, it can be applied to other variables with dead times less than MVsuch as stress, exercise, etc. As one can imagine, the relationships of such a set of variables on BGC is quite complex and accurately modeling their relationships and automatic feedback/feedforward control (FBFFC), with accurate announcements, appears to be the most viable one for success. In the content to follow, we focus on dynamic modeling with application to feedback control (FBC) and FFC when time delay in MV is significantly larger than time delay in disturbances. For this situation, when CV is the output, forecast modeling is necessary and when MV is the output, future announcement (i.e., knowledge) is needed for any variable with a dead time less than that of MV.
3. Discrete-time forecast dynamic modeling
Accurate forecast dynamic modeling in the context of process control has two critical applications. One is accurately forecastingCV at least θMV distance into the future, depending on the type of model-based control algorithm being used in FBC. The other one is an accurate cause-and-effectmodel for CV that is inverted for determining MV as a function of disturbances in FFC. Empirical modeling methods (EMM) (i.e., the so-called “data-driven” methods such as linear regression and artificial neural networks) are fit to a correlation structure and should not be used for cause-and-effectmodeling unless the modeling data are generated from a statistical experimental design covering the full range of the operating (i.e., input) space. This input space will be orthogonal and prevent extrapolation, which is risky for EM. For “freely existing data” or any data not generated from a statistical experimental design, accurate EM for CV forecasting is possible. However, when modeling data are not generated by a statistical experimental design, it would not be wise to use EM when cause-and-effectmodeling is needed since EMM are data driven and not knowledge driven, rely on high levels of parametrization, do not have structures or parameters that are physically interpretable based on first principles modeling, and are typically very risky for, even slight, extrapolation. In contrast, first principles model structures are: 1. nonlinear and thus, naturally break down correlation structures in the input data; 2. have physically interpretable parameters; and 3. often physical constraints with a theoretical basis. Nonetheless, theoretically based modeling of real data outside a controlled environment such as a lab, is often some combination of empiricism and first principles knowledge, which is essentially a “hybrid model” that is often called “gray box” or “semi-empirical” models. Models that are fully theoretical in derivation and structure but use data to obtain unknown physically interpretable model parameters, are classified in this document as semi-theoretical models.
Theoretically structured dynamic systems can be linearized (i.e., approximated) in time dependent variables (i.e., x = x(t)) while maintaining their time derivative structures (i.e., the order of derivatives will remain intact) and physical parametrization. For example, Eq. (4) represents the result of a dynamic overall mass balance on a process tank with one inlet stream with flow rate, q1(t), and one outlet stream though a hand valve with flow rate, q(t) = h2(t)/Rv, where his the tank level, Rvis the resistance to flow through the valve, and Ais the cross-sectional area of the tank. The density and temperature of the fluid in the tank is constant in this example. Using a 1st order Taylor Series approximation to linearize all time dependent variables in Eq. (4), gives the solution in Eq. (5), where the “′” represents a variable as a deviation from its initial steady state at t = 0.
where Eq. (6) is a first-order dynamic relationship with time constant, τ, and steady-state gain, K,and is represented in “standard form [4, 5].” Many dynamic processes can be approximated accurately by either a first-order-plus-deadtime (FOPDT) or second-order-plus-deadtime (SOPDT) structure [4, 5].
Eq. (7) gives a second-order version of Eq. (1) for inputs and unity gain, with yreplaced by Thus, Eq. (7) represents the dynamic response due to in the units of Eq. (8) represent the dynamic response of yas a function of where f(V) is an unrestricted mathematical function that maps each to the units of the output variable in standard form. Thus, it is f(V) that transforms the linear dynamic inputs into the nonlinear dynamic and static response for the output y.
where is a vector of the and the estimate of y(t), denoted as is equal to the estimate of This hybrid dynamic modeling structure is called a Wiener network  and is in a class of structures that are called “block-oriented models.” The block diagram for this network is shown in Figure 2.
Rollins, et al.  developed a multiple-input, single-output, discrete-time, non-linear Wiener dynamic approach using backwards difference derivatives based on Eqs. (7) and (8). Using a backward difference approximation applied to a sampling interval of Δt, an approximate discrete-time form of Eq. (7) is obtained (for pinputs):
where to satisfy the unity gain constraint with
After obtaining for each input i(i = pis MV), the modeled output value is determined by substituting these results into f(V), such as
Modification of Eq. (13) for forecasting a distance θMVinto the future with p = 3, for example, gives
where ai, bi, and ci,j, denote the linear, quadratic and interaction parameters for i = 1, 2, 3 and j = 2 and 3, and
where the θ’s are integer multiples of Δt.Depending on the rate of change of CV, forecast accuracy (and hence, control) can suffer significantly by setting θMV- θDVto zero. Developers of BGC devices that use current sensor glucose measurements in the feedback error restrict these devices for use only during long sleeping periods when BGC changes very slowly. For an application such as automatic BGC control, with a very large deadtime for MV and many disturbances with smaller deadtime than MV, that are nonlinear and interactive, the required accuracy for forecasting BGC is quite challenging. However, as health monitoring sensor technology continues to advance, forecast modeling accuracy continues to improve. The strengths of the method presented in this chapter are the use of dynamic structures that are embedded in first principles modeling; that is, they have physically interpretable parameters embedded in highly nonlinear structures (Eqs. (10)–(12)) with physical constraints such as , and for all i. While the method has these strengths for forecasting [16, 17], these strengths are quite critical in FFC applications where cause-and-effectmodeling is essential. Next PU is examined from a control perspective – FBC first and then feedback feedforward (FBFFC).
4. FBC when θMVis large
As discussed above, a change in MV will not affect CV until θMVtime into the future. When θMVis 0, the feedback error for FBC is, rightly, et = Yset– yt, where ytis the measured value of CV at the current time, t.When θMVis not 0, the equivalent feedback error is et = Yset– yt+θMVwhich is unknown because yt+θMVis not obtained until time t + θMV. This section describes and compares three FBC approaches when θMVis not 0. The first one is classical FBC [4, 5] which uses et = Yset– yt. The second one is feedback predictive control (FBPC)  which uses et = Yset–The third one is model predictive control (MPC) [18, 19, 20]. The MPC control law is for CV to be equal to Yset,time steps after t + θMVwhile holding the current value of MV fixed . Thus, is the only controller tuning parameter for MPC. More specifically, its feedback error is: et = Yset–Hence, the MPC prediction horizon is longer than FBPC by the amount For MPC, the optimal value of will tend to increase as the time lag of ytincreases for changes in MV. Since forecasting accuracy typically decreases as the distance into the future increases, MPC control can significantly deteriorate as increases.
These three control algorithms were compared by  in their ability to automatically control CV for a true FOPDT process with K = 1, τ = 10 and 50 min, θMV = 3 min, and a sampling time, Δt, equal to 1 min. FBPC and FBC were PI-Controllers with tuning parameters to give the best response with little, to no, overshoot for a unit step change in the set point at time t = 0.
Figure 3 presents the results of this study found in . As shown, CV (y) is on the left and MV (M) is on the right. The top row represents τ = 10 (J = 3, 8, and 20) and the bottom row represents τ = 50 (J = 8, 20, and 30). As shown, as Jdecreases, yreaches the set point faster and overshoots it for the lowest values of J. FBPC reaches the set point much faster than MPC, even when MPC overshoots the set point. As τincreases, MPC takes longer to reach the set point, but this is not the case for FBPC and FBC. FBC reaches the set point faster without overshooting than MPC for the case with the larger τ. Moreover, depending on Jand τ, FBC and MPC can reach the set point about the same time without overshooting the set point. However, FBPC has a faster response and reaches the set point much earlier than FBC and MPC in all cases. MV for FBPC has an initial “kick” much greater than FBC or MPC. However, its MV quickly drops below that of FBC and MPC and has significantly less movement in both cases as shown in Figure 3. Thus, as expected, because of the longer control horizon, which increases as τincreases, MPC responded slower than FBPC in reaching and staying at the new set point. Similar conclusions were seen in a comparison of FBPC and MPC in this article  for a simulated CSTR. Nonetheless, the main conclusion is that there are model-based forecasting FBC algorithms that are viable alternative to classical FBC when θMVis appreciably large.
5. FFC when θMVis large
All classical process control textbooks (e.g., [4, 5]) derive the FFC algorithm from a block diagram and give the FFC transfer function for each DV, Gf, as -GDV/GMV. The outputs from each Gfare added together to form the multiple-disturbance feedforward control law. An example when both GDVand GMVare FOPDT is:
Thus, Gfwill be PU when θMV > θDV, i.e., when Δθ > 0. Typically, this limitation is addressed by just setting Δθto 0 or increasing to and setting to 1 . This approximation is usually acceptable when Δθis small, as commonly found in chemical processes. However, modern applications of process control have gone beyond chemical processes to biological processes where transport is cellular, slow, and complex (i.e., not well understood). A common example is exogenous insulin taken by people with diabetes as mentioned above . Insulin deadtime is significantly greater than the deadtime for carbohydrate intake and other disturbances. Recent advancements in activity trackers measure multiple variables that likely affect BGC, and most, if not, all have a smaller deadtime than insulin [15, 16, 17].
Process Control textbooks commonly describe additive and linear dynamic FFC and present the algorithms in the continuous-time Laplace (s-) domain. This section presents a FFC approach that is: 1. given in the time domain; 2. discrete-time; 3. able to treat all types of non-additive behavior as well as nonlinear dynamic and static behavior and; 4. combines all disturbances functionally into one FFC law (i.e., all the DV’s enter one FFC equation). A block diagram of this FFC approach based on the Wiener network is given in Figure 4. As shown, the modeled disturbances are x1to xp-1and xpis MV. Inputs x1to xp-1pass through their dynamic blocks to produce v1to vp-1. The FFC law associated with Figure 4 is:
where is defined in Eq, (8) with the superscript eassociating it with the FFC law, i.e., value of MV that makes effc = 0, and thus, satisfying Eq. (19).
For p = 2, i.e., one disturbance and MV, and FOPDT structures for both inputs, and application of linear forms for Eq. (13) (for simplicity) into Eq. (19) gives:
Eq. (20) gives an explicit solution for the FFC signal, xMV, in this example. When f(V) has terms higher than first order, numerical root solving methods may be required to find xMV, as illustrated in .
Eq. (20) is evaluated now to determine if it can meet the standard of perfect control for x1load changes. For a frame of reference, MPC is also included in this study although it is a FBC method. For this example, Yset = 100 and remains constant. Input changes are made in x1(t) and its dynamic response to these input changes, v1(t), are given in Figure 5. The tuning parameter for MPC, J, has values of 1, 2 and 10. The model parameters are: a1 = 1, τ1 = 5 min, θ1 = 5 min; a2 = −1, τ2 = 10 min, θ2 = 10 min; and the sampling time, Δt = 1 min. The results for CV and MV for both FFPC and MPC are given in Figure 5. FFPC gives perfect control, as anticipated, and MPC does not, as anticipated. The response of MV that gives perfect control is the heavy black line in Figure 5. MPC with J = 1 appears to match the FFPC MV profile the best in terms of shape and time of changes, but it is also the most extreme. Thus, this example illustrates the ability of FFPC to meet the requirement of theoretically perfect control. Figure 6 gives a general multiple-input, block-oriented model FBFF block diagram similar to the one in Figure 4. For more information see .
FFPC is now evaluated on the in silicocontinuous stirred tank reactor (CSTR) in Figure 7 and described in  and taken from  with some minor modifications. This study has two DVs -- feed composition, CAi(x1), and temperature of the coolant entering the jacket, TCi(x2)). MV is the flowrate of the coolant entering the jacket, FC(x3). The output (CV) is the measured tank temperature, Tm(y). The model for each input is SOPDT, as shown in Eq. (7). The output, y, follows Eq. (8) and measurement noise was added to the true tank temperature (T) to produce Tm. Modeling this process was an application of Eqs. (9)–(17) with θ1 = θ2 = 5 seconds (s) and θ3 = θMV = 10 s. Thus, Δθ = 5 sfor each input and both are PU in the FFC law. Therefore, the objective of this study is to compare announcement of input changes 5 sahead versus no announcement.
A proportional-integral (PI) feedback controller was implemented in this study. Thus, FBPC was not used for FBC in any case to evaluate FFPC exclusive of FBPC. For this controller KC = 1.40, τI = 11.0 and MFBis the FBC signal to the control valve. MFF= is the FFPC signal. Thus, the signal to the valve, M, in Figure 7 is M = MFB + MFF. The input sequence used for training the model is given in Figure 8. The excellent fit of the model to Tmfor these input changes is also shown in Figure 8. The testing sequence (not shown) fit the response as well as the training sequence.
The results of FBC with FFPC for the two disturbances is shown in Figure 9. The left plot is for FBC only. The right plot is FBC with announcements for TCiand CAi. As shown, the variation of Tmaround its set point decreased greatly with FFC and announcements for both disturbances. More specifically, the standard deviation about the set point temperature dropped from 0.4352°C to 0.1131°C, a 74% reduction. Thus, modeling disturbances effectively and implementing them into FFC algorithms that can take advantage of announcements of future changes for critical disturbances can have a significant impact in reducing variation around the set point of CV.
As demonstrated above, FBPC is an effective FBC strategy when θMVis large. The Smith Predictor (SP) [4, 9] is a widely accepted FBC strategy when there is no deadtime in MV, and CV is measured, not at t, when MV changes, but at t + θCV(see Figure 10). The SP idea is to obtain a predictive model for CV with deadtime, remove the deadtime, and use this estimator in the feedback error term for CV at t. In Figure 11 an example of a SP process is illustrated with a block diagram representing the process given in Figure 10. As shown in Figure 11, MV = M(t) (Mis the signal of MV) immediately affects CV = B(t) (Bis the signal of CV), but the sensor is at B1(t) = B(t-θCV). Thus, SP uses a prediction of CV at twith an MV that changes CV immediately. In contrast, FBPC uses a forecast prediction of CV at t + θMVwhen it takes a time of θMVfor a change in MV to affect CV. Moreover, the SP does not use a forecast estimator and is not applicable to cases with deadtime in MV. However, this section will show that FBPC, using the forecasting estimator of CV at t + θCV, gives the same result as the SP. Consequently, FBPC can be used in place of the SP. However, the opposite is not true. More specifically, the SP is not applicable when a change in MV has time delay in changing CV, whereas FBPC is applicable for this case.
For a Figure 10 type process, the SP should compensate for the deadtime (i.e., reduce its effect) and respond quicker using an accurate estimate of Btthan using an accurate measurement of B1,t. The block diagram for the SP [4, 9] shows feedback control using with bias correction (BC) to address measurement bias. BC is the current measurement of where is the estimated value of Thus, in the SP block diagram,
Similarly, for FBPC,
The same simulated CSTR used above was used in this study to compare classical FBC, FBPC and the SP control algorithm with and without feedback correction. A step test in Mwas done and obtained B1over time from the initial steady state to a final steady state. These values are: M0 = 0.2569 and corresponding to B10 = 0.4000 and The input change was large enough to cover the change in Tmfor the test data (i.e., a 4°C change in the set point temperature). A FOPDT model was fit to the data and the fitted response is given in Figure 12. As shown, the fit is excellent with the following estimates: and with and Δt = 0.1. A plot of the response over time is given in Figure 12. The fitted forecast equation for is derived as follows:
The SP estimate without BC is obtained from Eq. (23) as
With BC, the SP estimate is
For a step change in the set point temperature of 4°C, the responses for tank temperature (Tm) for all four cases are given in Figure 13. The proportional-integral (PI) controller is the slowest to get to the new set point of 92°C. This is no surprise since the deadtime is quite large as shown in Figure 12. The SP with BC gives a modest improvement over PI which is surprising. This is also reflected in the modest increase in Kcfrom 0.20 to 0.41 as shown in Table 1. However, when the BC was removed, the SP response improved considerably as did Kcto 1.50. It gives the same response as FBPC, which also has no BC. These two cases give the same results, supporting the conclusion that FBPC and SP are equivalent for the SP application. However, when MV has deadtime with respect to CV, SP is not applicable, but FBPC is applicable. In  this BC method also did quite poorly with MPC being 132% worse (This BC method was not applied to FBPC in this study). Moreover, for BC, the following time series approach is recommended where the ϕ’s are estimated with all other parameters for fitted model :
SP w BC
SP w/o BC
Tuning values in the SP study.
This chapter has focused on the use of discrete-time dynamic forecast modeling to enhance FBC and all types of FFC. Discrete-time modeling has the advantage of obtaining solutions to PU systems without having to make assumptions to make the system an approximation of a physically realizable system. Models do not have to be cause-and-effectfor forecasting but need to be as FFC models. Cause-and-effect models result from statistical design of experiments because input changes are orthogonal (i.e., uncorrelated) and for theoretical structured models because they will be nonlinear in one or more physically based parameters, have physical constraints that must be met, and physically interpretable unknown model parameters .
When MV has deadtime with-respects-to CV (e.g., θMV), a change in MV will not begin to change CV until a time distance of θMVin the future. Three FBC approaches where evaluated in this scenario: classical FBC, FBPC and MPC, using the current measured value of CV, forecast estimate of CV, θMVin the future, and forecast estimate of CV, θMV + JΔtin the future, respectively. In the simulation study, FBC responded quicker than MPC when the process lag was large and MPC responded quicker when the process lag was small. FBPC responded much faster than both under small and large lag. FBPC control has a prediction horizon of θMVbut for MPC it is JΔtlonger. Since the optimal value of Jincreases as the lag increases, MPC can be significantly more sluggish than FBPC when Jis large. A definite advantage of MPC is that Jis its only tuning parameter.
A discrete-time FFC approach (FFPC) was presented in this chapter that can be effective when θMVis large and the multiple-input FFPC model is PU for any reason (i.e., the order of the differential equation or negative deadtime). FFPC was shown to satisfy perfect theoretical control in a simulated data study. A critical strength of the approach presented in this work is that the FFC variables enter one mathematical function that simultaneously solves for one FFC control signal. This contrasts with classical FFC that has a FFC algorithm for each input and combines their values to determine the value of the FFC control signal for MV. The classical approach cannot treat complex interactive and nonlinear behavior of the disturbances in determining the optimal value of the FFC signal for MV. Block diagrams of this novel FFC approach were shown for the Wiener Network and a general block-oriented modeling approach. When FFC inputs have a PU impact, knowing how their values will change over the control horizon (i.e., announcements), can significantly improve FFC as demonstrated in the CSTR simulation study.
The SP is a model-based feedback control algorithm that can be quite effective when there is no deadtime between a change in MV and its impact on CV, and the measured value of CV has deadtime. For this situation, FBPC, that uses a forecast value for CV based on a model developed from the measured value of CV with deadtime, gives the equivalent result of the SP. However, the SP is limited to this case, but FBPC is not. More specifically, FBPC is applicable when there is deadtime for changes in MV and its effect on CV but the SP is not. Finally, one should exercise care when using the bias correction (BC) method in the block diagram for the SP. It can lead to a significantly suboptimal SP as shown in this work. A better alternative is to use one that is obtained from modeling the serially correlated structure in the process as given in this chapter. A flowchart illustrating the complete process of the proposed framework of this chapter is shown in Figure 14.
I thank Kendra Kreienbrink for helping with the references and for proofreading this document.
cross-sectional area of the tank
the nth constant
the signal of CV in the SP algorithm
bias correction in the SP algorithm
the signal of CV with deadtime θCV
the mth constant
inlet concentration to the CSTR (x1)
the feedforward control law criterion for perfect control
the feedback error = Yset – yt,
inlet flow rate of the coolant to the jacket of the CSTR
the true function in the FFC law that satisfy effc=0
is an unrestricted mathematical function that maps each vi′ to the units of the output variable in standard form. Thus, it is f(V) that transforms the linear dynamic inputs into the nonlinear dynamic response for the output y.
feedforward transfer function
DV transfer function
transfer function with output vi
transfer function for CV in signal
MV transfer function
process transfer function
transfer function for MV
the tank level
the only controller tuning parameter for MPC.
controller gain for FBC
mth order derivative
the input signal to MV transfer function
FBC input signal to MV transfer function
FFC input signal to MV transfer function
nth order derivative
inlet flow rate of Stream 1
flowrate in the Laplace domain entering the system
flow rate of the outlet stream
current time or just time
the physical value of CV in SP case
inlet temperature of the coolant to the jacket of the CSTR
measured tank temperature
dynamic output for ith input, xi, in the same units as xi
is a vector of the vi″s
the value of x at time t -Δt
deviation of xi at time t from xi at time t = 0
xMVe equals the value of MV in the FFC law that satisfy effc=0
Derrick K. Rollins (August 19th 2021). Use of Discrete-Time Forecast Modeling to Enhance Feedback Control and Physically Unrealizable Feedforward Control with Applications [Online First], IntechOpen, DOI: 10.5772/intechopen.99340. Available from:
Over 21,000 IntechOpen readers like this topic
Help us write another book on this subject and reach those readers
We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.