Input Shaping Predictive Functional Control for Different Types of Challenging Dynamics Processes

Predictive functional control (PFC) is a fast and effective controller that is widely used for processes with simple dynamics. This paper proposes some techniques for improving its reliability when applied to systems with more challenging dynamics, such as those with open-loop unstable poles, oscillatory modes, or integrating modes. One historical proposal considered is to eliminate or cancel the undesirable poles via input shaping of the predictions, but this approach is shown to sometimes result in relatively poor performance. Consequently, this paper proposes to shape these poles, rather than cancelling them, to further enhance the tuning, feasibility, and stability properties of PFC. The proposed modification is analysed and evaluated on several numerical examples and also a hardware application.


Introduction
Predictive functional control (PFC) is a simple controller that is effective for small-scale single-input-single-output (SISO) applications, especially for low-order and stable processes [1][2][3]. The main advantages of PFC compared to its prime competitor-that is, proportional integral derivative (PID)-are its ability to handle constraints and its transparent tuning parameters. Indeed, it must be emphasised that the performance of PFC should not be benchmarked against more advanced model predictive control (MPC) strategies [4], because the implementation is relatively much cheaper and requires only low computation with very straightforward coding [5,6].
Nevertheless, despite its apparently attractive benefits, PFC has received relatively little attention in the academic literature because of the lack of a priori stability guarantees [7,8], which are possible with more advanced MPC approaches. Consequently, several recent works have developed the basic concept of PFC to improve its overall tuning properties while providing a confident assurance in the resulting closed-loop performance and constraint handling [9][10][11][12]. However, most of these modifications perform well only with low-order and simple dynamical processes. For a system with open-loop unstable poles, significant underdamping, or integrating dynamics, PFC is quite challenging to tune effectively [13], and the resulting divergent or oscillating predictions may give rise to infeasibility and/or robustness issues.
PFC practitioners often deploy a type of cascade structure to handle a challenging dynamics process, where an inner loop is used to improve the dynamics for an outer loop to control [14,15]. This modification enables a user to retain an independent model (IM) structure that handles uncertainties, while retaining a similar tuning concept as the conventional approach. Nevertheless, the inner PFC can only accept a proportional-type controller to avoid any overcomplication when handling constraints [6]. This restriction makes it difficult to have a systematic selection of gain, rather than an ad hoc approach. Moreover, a back-calculation or anti-wind-up technique is required to avoid any saturation while satisfying the constraints [6], which also can be conservative and thus affect performance.
A fundamental conceptual weakness of a simplistic PFC approach is that one is basing decisions on an open-loop prediction, which may have undesirable, possibly divergent, dynamics; matching this to a desirable closed-loop dynamic will lead to an ill-posed approach. Hence, an alternative approach is to pre-stabilise/pre-shape the output predictions by shaping the future input dynamics so that the effect of unwanted poles on the predictions are alleviated or cancelled [16]. This modification can retain the systematic tuning concept of PFC, in addition to facilitating a recursive feasibility guarantee feature for constraint handling. However, the performance of this control law is not always desirable, since the cancellation of specified modes (poles) within the predictions using a minimal number of control input changes often requires an aggressive input trajectory [16][17][18], which is not implementable or ideal for some cases.
In practice, less aggressive shaping to remove the undesirable modes from the prediction [19] is preferable, and this forms the main thrust of the proposal in this paper. More specifically, the proposal here is that, rather than using a small, finite number of control moves (effectively, finite impulse response (FIR) parameterisation), such as in Generalized Predictive Control (GPC) Added the defination and conventional PFC, the future input moves will be parameterised using an infinite impulse response (IIR) instead. The preference for IIR over FIR is due to IIR being more convenient to manipulate and define than a high-order FIR, in general. In turn, by allowing the output modes to evolve over many more samples, the required input will be less aggressive. Nevertheless, due to the desire for simplicity and transparency that is a central tenet of PFC, in this paper, we choose not to generalise the parameterisation for different dynamics. Instead, this work seeks to provide some rigour on how to effectively and systematically shape the future input for a given dynamic and, moreover, how to ensure some recursive feasibility properties during constraint handling. Section 2 of this paper presents a brief formulation of conventional PFC. Section 3 introduces the concept of input shaping PFC, together with the constraint handling approach. Section 4 demonstrate the proposed algorithms on several numerical examples and also on some laboratory hardware. Section 5 provides the conclusions.

Conventional PFC
This section presents, in brief, the main concepts, notation, and formulation of PFC, together with a systematic constraint handling approach. For more detailed derivations, theory, and concepts of PFC, an interested reader may refer to the references, e.g., [5,6,13,20]. Without loss of generality, a controlled autoregressive and integrated moving average (CARIMA)-based model is used instead of the standard independent model (IM) structure to derive the required unbiased predictions, as this form is more amenable to the algebra required to implement the shaping. Hence, the model will take the form: where b(z) = b 1 z −1 + ... , a(z) = 1 + a 1 z −1 + ..., y k , u k are the outputs and inputs, respectively, at sample k, and ζ k is an unknown zero mean random variable used to capture uncertainty.

Control Law
PFC design is based on the assumption that a closed-loop response should follow a first-order dynamic from the current state y k to the desired target r [20]. In practice, one aims to achieve this by ensuring a matching using the open-loop prediction, but only at a single point. In other words, the predicted output trajectory is chosen to satisfy the following equality: where y k+n y |k is the n-step ahead system prediction at sample time k, and λ is the desired closed-loop pole that will determine the speed of convergence. PFC practitioners typically select the desired closed-loop time response (CLTR) which has an explicit link with the corresponding closed-loop pole, that is, λ = e −3T CLTR , where T is the sampling period [6]. The predictions for the CARIMA model (1) are standard in the literature (e.g., [4,21]), so only the final form is given here. For input increments ∆u k+i , the n-step ahead future output prediction is formed as: where the left and right arrow underlying vectors represent past and future variables, respectively. The parameters H, P, Q depend on the model parameters, and for a model of order m: In conventional PFC [6,20], within the predictions we select ∆u k+i = 0, i > 0. Uusing this and substituting the n-step ahead prediction from (3) into equality (2) gives: where e i is the ith standard basic vector and H n = e T n H, P n = e T n P, Q n = e T n Q. Hence, the PFC control law is given as: Remark 1. Figure 1 shows a comparison of simplified flow diagrams, where both PFC and MPC share the same structure, yet have a different optimisation process, where the constraint handling is embedded inside the main block. As for PID, the control input is obtained simply by tuning the gains, while the constraints are handled via a rule base [22].

Remark 2.
It is noted that PFC performs well for a system with close to a monotonic step response, such as a first-order system and overdamped second-order dynamics, assuming, of course, a sensible choice of the tuning parameters λ and n [11][12][13]. However, the same tuning procedure may not work for systems with less simple open-loop dynamics, leading to ill-posed decision making and unreliable control.

Constraint Handling
The constraint handling approach presented here is adapted from the standard MPC literature [21,23] but with simpler code, and it is more systematic and less conservative than the back-calculation typically used in conventional PFC algorithms [24]. Similarly, it will be more efficient than the ad hoc approaches used with PID. Assume constraints, at every sample, on input rate, input, and states, as follows: The user needs to: 1. Find a suitable horizon m [23] over which to compute the entire set of future output predictions in a single vector: y →k+1|k = He 1 ∆u k + P ∆u k ← + Q y k ← . The horizon for output predictions y →k+1 , and thus the row dimension of H, should be long enough to capture all core dynamics! 2. Combine the input constraints, rate constraints, and output predictions into a single set of linear inequalities: where f k depends on past data in ∆u k ← , y k ← and on the limits. 3. The input/output predictions will satisfy constraints if inequalities (8) are satisfied, and thus the PFC algorithm should consider these explicitly.
Next, a single simple loop is utilised to find the ∆u k closest to the unconstrained solution of (6), while satisfying (8).

Algorithm 1. Simple PFC algorithm with systematic constraint handling
At each sample: 1. Define the unconstrained value for ∆u k from (6). 2. Define the vector f k of (8) (it is noted that C does not change). 3. Use a simple loop covering all the rows of C as follows: Remark 3. For a simple and stable open-loop process and for suitably large m, Algorithm 1 is guaranteed to be recursively feasible and, moreover, to converge to a possible value for ∆u k that is closest to the unconstrained choice [12]. However, this benefit does not apply to systems with more challenging dynamics, such as when the open-loop prediction is divergent.

Input Shaping PFC
This section presents the concept of input shaping to pre-stabilise (or pre-condition) the open-loop predictions of processes with undesirable dynamics. The shaping of input predictions can be done either via explicit pole cancellation or pole shaping, and both methods can be used to formulate a PFC control law which retains a recursively feasible constrained solution. A key issue, however, is whether one deploys FIR or IIR parameterisations of the predictions for the future input increments ∆u k+i|k , i ≥ 0.

Pre-Stabilisation of Predictions via Pole Cancellation
For systems with poor open-loop dynamics, the constant input assumption of typical PFC does not provide enough flexibility within the predicted input to both cancel the effect of an undesirable pole and to get nice convergent behaviour [13,17,19]. Thus, it is crucial to first stabilise the prediction before implementing it into a control law [25]. The first step is to factorise the poles in the denominator: where a + (z) contains the undesirable poles. Utilising the Toeplitz/Hankel form [21], the future output predictions can be computed as: where for general polynomial f (z) Rearranging prediction (10) in more compact form, we get: from which the presence of the undesirable modes are transparent through the factor [C a + ] −1 .

Lemma 1.
Selection of the future input sequence ∆ u →k , at each sample, such that the following equality is satisfied: where γ is a convergent sequence or a polynomial, will ensure that the corresponding output predictions in (11) do not contain the undesirable modes in a + .
Proof. This is self-evident by substitution of (12) into (11), which gives: so that only the acceptable modes in a − (z) remain in the predictions, along with any components in the convergent sequence γ. It is noted that this choice automatically includes the initial conditions within p and thus updates each sample as required. (12) can be solved by a small number of simultaneous equations [21], where the minimal-order solution can be represented as:

Remark 4. Requirement
for suitable P 1 , P 2 . The required dimension of non-zero elements in vector ∆ u →k corresponds to at least one more than the number of undesirable modes (n a + ), while the order of γ is usually taken as n γ = n p − n a + , where n p is the effective dimension of p (which depends upon the column dimensions of H b , H A ).
To ensure the future manipulated control moves are convergent, while adding some flexibility for modifying the output predictions, the input requirement in (14) can be enhanced to: where the vector parameter φ denotes the degrees of freedom (d.o.f.) within the predictions.

Theorem 1.
Using the new shaped input (15) ensures that the undesirable modes do not appear in the output predictions, irrespective of the choice of φ. The output predictions are convergent if φ is finite-dimensional or a convergent infinite dimensional sequence.
Proof. Substitute input prediction (15) into output prediction (12), and the predictions become: The prediction can be represented with an equivalent z-transform: It is known from Lemma 1 that the contribution from γ gives a convergent prediction, and thus overall convergence is obvious as long as φ(z) is convergent (or an FIR).

Remark 5.
Noting the definition of p in (11), the n-step ahead output prediction with prediction class (15) and (17) can be put in a more common form as: where H s,n , P s,n , and Q s,n are suitable matrices, and the additional subscript 's' is used to denote shaping and φ is taken to be FIR (equivalently a finite dimensional vector). Note, however, that typically for PFC, φ is a scalar. Also, it is easy to show [18] that choosing φ = 0 will automatically give the same input predictions as those deployed at the previous sample, which enables consistency of predictions from one sample to the next.

Pre-Stabilisation via Pole Shaping
It is known that dead-beat pole cancellation can require aggressive inputs, and the minimal-order solutions to (12) are in effect dead-beat input predictions [16,19]. Although dead-beat solutions are easy to define and thus have some advantages in terms of computation and transparency, in practice, a user may desire a less aggressive shaping that is more implementable in a real system. Alongside this, the popularity of dual mode approaches in the literature [26] is partially because they allow the implied input predictions to converge to the steady state asymptotically, rather than in a small, finite number of steps. Thus, a logical question to ask is whether a smoother solution to (12)-that is, one where the implied solutions for γ(z), φ(z) used in (17) have some poles, say α(z)-would work better for PFC.
The mainstream MPC community has focussed on optimal control solutions, but, given that PFC is intended to be simple and low-dimensional, the proposal here is that it is more reasonable to investigate the potential of simple default choices for the asymptotic dynamics α(z) within the input and output predictions. Clearly, this choice can be strongly linked to the target closed-loop behaviour and/or system knowledge. Proposal 1. By definition, the integrator has a pole on the unit circle-that is, factor (1 − z −1 )-and, conversely, cancelling the pole as in (12) is equivalent to enforcing a pole on the origin-that is, factor (1 − 0 z −1 ). Hence, the choice of pole factor α = (1 − 0.5 z −1 ) represents a simple half-way house trade-off between these two choices.

Proposal 2.
For a process with significant underdamping, the implied α(z) will have only real poles which are chosen to be close to the real parts of the oscillatory poles. This will reduce the undesirable oscillation in the output predictions, but not change the convergence speed, albeit the input may then be somewhat oscillatory.

Proposal 3.
For open-loop unstable systems, a simple default solution simply inverts the unstable poles, that is, defining α(z) such that a + (z i ) = 0 ⇒ α(1/z i ) = 0. Lemma 2. The dynamics α(z) will be present in the predictions if the following Diophantine equation is used to solve the input/output prediction pairing. (19) is equivalent to solving: and, moreover, Equation (20) follows directly from enforcing (12) while assuming ∆ u →k = C −1 α w. Hence, substituting this ∆ u →k into (11) gives: It is evident, therefore, that the desired poles are in the predictions for both the input and output.

Remark 6.
The new requirement (20) can be solved similarly to (12), where the minimal-order solution for w andγ are: Theorem 2. A convergent prediction class which embeds both the desired asymptotic poles and some degrees of freedom (d.o.f.) can be defined from: where convergent IIR or FIR φ constitutes the d.o.f.
Proof. This is analogous to Theorem 1 and is based on superposition. The additional component in w-that is, C a + φ-necessarily cancels the undesirable poles and gives overall convergent output predictions. So, using (21), then: Added hat on top of gamma in Equation (24).
Remark 7. By extracting the n th row and noting the definition of p in (11), the n-step ahead prediction from (24) can be rearranged in a more general form as: for suitable h n,α , P n,α , Q n,α , and it is noted that as is conventional for PFC, φ has just a single non-zero parameter in order to retain computational simplicity and to have just a single d.o.f. for satisfying the control law (2).

Proposed Shaping PFC Control Laws
Since the shaped predictions of (18) and (25) are derived in a general form, two new control laws-Pole Cancellation PFC (PC-PFC) and Pole Shaping PFC (PS-PFC)-can be formulated in a conventional manner after selecting a suitable coincidence horizon n and desired closed-loop pole λ.
[PC-PFC] The d.o.f φ is computed by substituting prediction (18) of PC-PFC into equality (2), and thus: then, the current input increment ∆u k is determined simply by inserting φ into the predicted input of (15).
[PS-PFC] The d.o.f φ is computed by substituting prediction (25) of PS-PFC into equality (2), and thus: then, the current input increment ∆u k is determined simply by inserting φ into the predicted input of (23).

Constraint Handling Approaches with Recursive Feasibility
A core advantage of MPC, in general, is that the optimised predictions can be restricted to ones which satisfy constraints; the d.o.f. within the predictions are used to ensure constraint satisfaction. For PC-PFC and PS-PFC, the d.o.f. in the predictions is the variable φ. This section gives a brief overview of how the constraint inequalities ensuring (7) depend upon φ. We will use PS-PFC and assume that the reader can easily find the equivalent matrices for PC-PFC (for which, in effect, α = 1).
Noting the definition of future input increments in (23) and output predictions in (25), the constraints inequalities for (7) can be defined as: where C I/∆ is a lower triangular matrix one ones, and L is a vector of ones with an appropriate dimension (typically a horizon long enough to capture the core dynamics in the predictions). The reader should note that the horizon for the predictions used in (28) will, in general, be much longer than the coincidence horizon used in (27), as one needs to ensure that the implied long-range predictions satisfy constraints. The inequalities can be combined for convenience as follows, although this is not necessary for online coding where efficient alternatives may exist:

Algorithm 2. [PS-PFC with constraint handling]
At each sample: 1. Define the unconstrained value for φ from (27). 2. Update the vector f k of (29) (it is noted that C does not change). 3. Use a simple loop covering all the rows of C as follows: (a) Check satisfaction of the ith constraint using: e T i Cφ ≤ f k,i . (b) If e T i Cφ > f k,i , then set φ = ( f k,i )/[e T i C], else leave φ unchanged.

Theorem 3.
In the presence of constraints, Algorithm 2 is recursively feasible where the computed φ will not only enforce the input/output predictions to satisfy constraints at the current sample, but will also guarantee that one can make the same statement at the next sample.
Proof. By definition, the choice of φ = 0 ensures feasibility in the nominal case because the input componentP 1 p is the unused part of the input prediction from the previous sample, and this is known to satisfy constraints by assumption. One can ensure feasibility at start-up by beginning from a sensible state.
Readers should note that using the pre-stabilised/shaped predictions is essential for this recursive feasibility result, which is not available for more conventional PFC approaches, for which the implied long-range predictions may be divergent. Thus, this Theorem is an important contribution of this paper.

Remark 8.
Although Algorithm 2 allows recursive feasibility, which is a strong result, ironically, the use of PC-PFC or PS-PFC does not not give any a priori stability and/or performance guarantees in general, which is a well-understood weakness of PFC approaches [7] and a consequence of wanting a very simple and cheap control approach.

Numerical Examples
This

Description of Case Studies
The case studies presented here are inspired from the real process applications. However, for clarity of presentation, the model parameters are not specific to a given piece of apparatus, but rather are generic to attain suitable dynamics which enable an explicit comparison between the control laws. In this work, it is assumed that there is no plant model mismatch. The robustness properties of these controllers will remain as future work, although it is noted that, as with most predictive controllers, disturbance rejection and offset free tracking is implicit.

Case Study 1: Boiler Level Control
In the process industry, the use of a boiler is frequent, and the level of water needs to be controlled within the manufacturer's specified limits. Exceeding the allowable limits may lead to water overflow, overheating, and/or damage to many components. Conversely, if the level is low, the water wall tubes may overheat and cause tube ruptures, resulting in expensive repairs and other potential hazards. Hence, the prime control objective is to raise the water level at the boiler start-up point while retaining it at a constant steam load. Since the conversion process from water to steam is very slow, a typical model for this process is usually a first-order system with an integrator and stable zero [27]. In a discrete form, one of the poles should reside in a unit circle. The relationship between the output water level (m), y(z), and the input water flow rate (m 3 s −1 ), u(z), can be represented by a representative model, such as G 1 :

Case Study 2: Depth Control of Unmanned Free-Swimming Submersible (UFSS)
In the marine application, the depth of an unmanned submarine can be controlled by deflecting its elevator surface, whereby the vehicle will rotate about its pitch axis; the associated vertical forces due to the water flow beside the vehicle enable the vehicle to sink or rise. Since a step input deflection may create an oscillatory angle of dive due to the water current, typical dynamics to represent this system often consist of at least one stable pole and two complex poles with stable zeros [22]. Thus a representative third-order underdamped process G 2 can be assumed to represent this pitch control system: with the output as the pitch angle (rad) and input as the input elevator deflection (m).

Case Study 3: Temperature Control of Fluidised Bed Reactor
A fluidised bed reactor is used to produce a variety of multiphase chemical reactions that are highly exothermic and can be considered as unstable. The reactor bed temperature needs to be controlled by manipulating the coolant flow rate to avoid overheating and other potential hazards. In this case, a drastic change in flow rate will trigger a reaction between the chemicals that releases extra energy and increases the bed temperature. In fact, the change in flow rate needs to follow specific dynamics to avoid this reaction while stabilising the temperature. A reduced control model includes at least one stable pole, typically, and one unstable pole [28] to relate the dynamics between the output temperature ( • C) and input coolant flow rate (m 3 s −1 ). Inspired by this system, a representative second-order unstable process G 3 is considered as a good case study:

The Impact of Input Shaping on Predictions and Feasibility
The prime purpose of shaping the future input dynamics is to eliminate the effect of unwanted poles in the future predictions. Nevertheless, it is also undesirable to have an overaggressive input activity, which may not be implementable in a real plant. To analyse this issue, the prediction behaviour of PFC, PC-PFC, and PS-PFC for processes G 1 , G 2 , and G 3 are plotted in Figure 2. From these results, it can be observed that: • For an integrating process, such as G 1 , the constant input prediction of PFC leads to a divergent output prediction, and thus, output constraints can only be satisfied if the input is selected to be zero! Hence, the PFC plots do not appear in this example, as the constraint handling forces a choice of u k = 0, ∀k. • For G 1 , the default input prediction (Equation (15)) for PC-PFC (blue-dotted line) is of dead-beat form and aggressive, whereas the input prediction (Equation (23)) for PS-PFC moves smoothly to the steady state and is less aggressive. There is no significant difference in the corresponding output predictions. • For the underdamped system G 2 , the output prediction from PFC includes a significant oscillation, which is undesirable and will also cause constraint handling to be conservative. The differences between PC-PFC and PS-PFC predictions are similar to those noted for G 1 , that is, PS-PFC has a much smoother and less aggressive input prediction, albeit slow, and output prediction, due to the choice of α. Of course, this difference means that the constraint handling for PS-PFC will be far more preferable and less conservative, in general. Conversely, since PC-PFC cancelled out two of their oscillatory open-loop poles, a sudden spike or aggressive damping is expected in the output response. • For the unstable process G 3 , a conventional PFC cannot be used because the divergent predictions will automatically violate constraints so that no feasible choice for u k will exist. Once again, it is seen that the predictions for PS-PFC are preferable to those from PC-PFC. In summary, PS-PFC produces the best prediction behaviour because it ensures convergent predictions which will satisfy constraints with less aggressive input predictions, and thus less conservative constraint handling, than given by PC-PFC/PFC.

Tuning Efficacy and Closed-Loop Performance: The Unconstrained Case
First we give a brief discussion on PFC tuning for completeness. In general terms, a good practice guidance is to select the coincidence horizon in between 40% and 80% rise of the step input response to the steady-state value [13]: here, G 1 (4 ≤ n ≤ 9), G 2 (8 ≤ n ≤ 15), and G 3 (11 ≤ n ≤ 21). Selecting a smaller horizon will lead to a more aggressive input, while larger horizons reduce the efficacy of λ as a tuning parameter.
In general therefore, the main designer choice is the desired closed-loop pole; for simplicity of illustration, we take the desired pole to be λ = 0.7. Figure 3 demonstrates the closed-loop performance of PFC, PC-PFC and PS-PFC on the three case study processes with these tunings. It is noted that: • For all cases, PS-PFC (using a default choice of α) gives the best trade-off between the rate of convergence and the aggressiveness of input activity compared to PFC and PC-PFC. Changes to α could offer a further tuning parameter for varying this trade-off.
• It is notable that PS-PFC gives similar or better output behaviour to PFC/PC-PFC while using a far smoother and less aggressive input trajectory. • For processes G 2 , G 3 , the input and output behaviour of PC-PFC is extremely aggressive and would not be implementable in a real application. • For process G 3 , the conventional PFC cannot be stabilised with the given choice of n.

Constraint Handling
As noted in Section 4.2, for many dynamics a conventional PFC approach is infeasible or highly conservative because the output predictions inevitably violate constraints beyond a given horizon. Thus, conventional PFC can only be implemented by using short constraint horizons and thus with the loss of any recursive feasibility assurance; if it does work this is luck not good design and thus should be avoided. For the case studies here, conventional PFC could only be used safely with output constraints for G 2 , although in that case we would expect some conservatism due to the oscillations in the output predictions.
PC-PFC and PS-PFC pre-stabilise the output predictions and thus can be used safely and with a recursive feasibility assurance. Figure 4 compares the performance of PFC, PC-PFC, and PS-PFC when constraints are added to the process. Several observations can be noted: • As expected, PS-PFC and PC-PFC satisfy constraints, retain recursive feasibility throughout and converge safely. • For process G 1 , the constrained performance of the controllers are almost the same, but PS-PFC provides a smoother input transition. • For process G 2 it is clear that handling the under-damping will cause some challenges to any control law, but clearly PS-PFC provides the best responses. • For process G 3 , the inherent dead-beat input predictions deployed by PC-PFC mean the performance is poor and slow to converge whereas PS-PFC performs well.
Closed-loop constrained behaviour for system G 1 with n = 4 and −0.

Application of PS-PFC on Laboratory Hardware
This section demonstrates the implementation of PS-PFC on laboratory hardware, that is, a Quanser SRV02 servo-based unit (Quanser, Markham, ON, Canada) (see Figure 5). This device is powered by a Quanser VoltPAQ-X1 amplifier and operates by National Instrument ELVIS II+ (National Instruments, Austin, TX, USA) multifunctional data acquisition via USB connection. The control objective is to track the servo position θ(t), measured in radians, by manipulating the supplied voltage V(t). This servo will rotate counter-clockwise with positive supplied voltage and vice versa. A second-order model of (33) with an integrator is used to represent the system dynamics (refer to [29] for a more formal derivation) as: θ(t) = 1.53 s(0.254s + 1) V(s) (33) To implement the proposed control law, the continuous model (33) is discretised with sampling time 0.02s to obtain the discrete model of: The plant is set to have a CLTR of 0.2s (which is equivalent to λ = 0.89). Using a similar procedure to that in the previous section, the coincidence horizon is selected at n = 4. Figure 6 demonstrates the unconstrained and constrained performances of PS-PFC in tracking the servo position. In the unconstrained case, the controller manages to track the target position with the desired convergence speed. As for the constrained case, all the implied input limits (−8 ≤ u k ≤ 8), rate limits (−3 ≤ ∆u k ≤ 3), and output limits (−0.8 ≤ y k ≤ 0.8) are satisfied systematically without any conflict.

Conclusions
This paper proposes two potential simple modifications to a conventional PFC algorithm to improve the constraint handling properties for processes with challenging dynamics, such as integrating modes, underdamping, or unstable modes. Both proposals use relatively simple algebra-in effect, the solution of a small number of linear simultaneous equations-to parameterise the future input trajectories which lead to convergent and desirable output behaviours; this is done in terms of a component to deal with the current state (or initial condition) and a free component for control purposes. A core contribution is to show that using the proposed parameterisations allows a simple proof of recursive feasibility so that the constraint handling can be performed more safely and reliably.
A specific novelty of this paper is the proposed PS-PFC algorithm, which gives a pragmatic and simple proposal for deriving input and output prediction pairs which do not require aggressive inputs during transients; the more classical alternative approach of PC-PFC, in general, deploys very aggressive inputs in transients and thus cannot be used in practice. Simulation evidence on a variety of simulation case studies and hardware demonstrates that the proposed PS-PFC algorithm significantly outperforms both a conventional PFC approach and PC-PFC.
Although, as with all predictive control laws, both PC-PFC and PS-PFC are robust to some parameter uncertainty and disturbances, a detailed sensitivity analysis is an important next step. Also, it is of particular interest to compare these approaches with the alternative feedback formulations for PFC [6,15] by way of systematic design, nominal performance, constraint handling, and sensitivity.
Author Contributions: This paper is a collaborative work between both authors. J.A.R. provided initial proposals and accurate communication of the concepts employed in previous MPC and PFC control laws while reviewing the whole project. M.A. developed the code and analysed the concept in various challenging dynamics process within PFC framework.
Funding: This research received no external funding.