Model Predictive Control Tuning by Inverse Matching for a Wave Energy Converter

: This paper investigates the application of a method to ﬁnd the cost function or the weight matrices to be used in model predictive control (MPC) such that the MPC has the same performance as a predesigned linear controller in state-feedback form when constraints are not active. This is potentially useful when a successful linear controller already exists and it is necessary to incorporate the constraint-handling capabilities of MPC. This is the case for a wave energy converter (WEC), where the maximum power transfer law is well-understood. In addition to solutions based on numerical optimization, a simple analytical solution is also derived for cases with a short prediction horizon. These methods are applied for the control of an empirically-based WEC model. The results show that the MPC can be successfully tuned to follow an existing linear control law and to comply with both input and state constraints, such as actuator force and actuator stroke.


Introduction
Control design for wave energy converters (WECs) is a critical component in the broader WEC design problem. A WEC is generally designed to achieve a desired performance goal, such as maximization of electrical power transferred to the grid while minimizing levelized cost of energy and should survive in extreme ocean environments, under various design constraints. This imposes numerous engineering challenges, comprising environmental characterization, modeling, failure detection, and fatigue analysis. Active control strategies are usually developed to achieve the desired goals, but they are significantly related to the WEC design process in that they alter the overall WEC dynamics in terms of system performance and development costs. Additionally, various types of constraints resulting from WEC design or environmental conditions must be considered when developing control systems.
A number of control design approaches and frameworks have been applied to WECs in the past. In regular waves, maximum energy is captured when the natural frequency of the WEC device is close to the dominant frequency of the incoming wave [1]. At resonance, the velocity of the WEC oscillator is in phase with the excitation force induced by the wave. Complex conjugate (CC) control or impedance matching described in the frequency domain enables this resonance and leads to the theoretically optimal device velocity guaranteeing maximum energy extraction [2]. In realizing this CC control strategy in practice, the main challenge to be overcome is that its general solution is noncausal and thus requires a good prediction for the upcoming wave force. Since the requirement under real irregular sea conditions is that the device should behave as if resonant over the wide range of frequencies, Budal and Falnes [3] proposed the so-called latching control in which the motion of the device is stalled at the extremes of its movement and released when its phase is in good agreement with the wave force. In Korde [4] an efficiency improvement using latching control was demonstrated by conducting an experimental study. In order to expand the efficiency range on either side of the resonant frequency, reactive control [5] was suggested by dynamically adjusting the WEC parameters for maximum energy absorption. Korde [4] showed that the damping coefficient could be adjusted by velocity feedback to maximize permissible energy absorption. In Scruggs et al. [6] the optimal controller was found by solving a nonstandard linear quadratic Guassian optimal control problem. Various other control algorithms have been also applied, including dynamic programming [7], pseudo-spectral methods [8], fuzzy logic control [9], shape-based control [10], moment matching [11], and model predictive control [12].
Among practical control considerations, model predictive control (MPC) possesses a unique feature, in that it takes into account constraints from the beginning of the controller design process [13]. Conversely, the analysis of the stability, robustness, and sensitivity of MPC is more intricate to handle in contrast to linear feedback controllers [14]; these factors often act as the main hurdles for implementation of MPC in real-life applications. Moreover, it is often difficult to find an appropriate cost function for MPC optimization. The vast experience gained from linear controller design can provide many benefits for MPC design and tuning the MPC such that it has the same behavior as a well-designed linear controller when constraints are not active is an attractive approach. Several researchers have attempted to link MPC to an existing linear control law to improve performance due to the constraint-handling capabilities of MPC. In Maciejowski [15] a procedure was presented to construct an observer-based MPC realization of an existing output-feedback controller such that the MPC is equivalent to the given controller when all constraints are inactive. This methodology was extended in Foo and Weyer [16] using a reduced-order observer and state augmentation. Cairano and Bemporad [17] presented approximate solutions for the weight matrices of MPC that behaves as a given linear controller in state-feedback or linear quadratic regulator form in case constraints are inactive by recasting the problem as a convex optimization problem with linear matrix inequality constraints. The global stability of the closed-loop system in the presence of constraints is also guaranteed. Kong et al. [18] further generalized this kind of inverse matching by regarding MPC as a desirable upgrade of an existing control law for performance improvement. Introducing a confidence parameter λ ∈ [0, 1], Kong et al. [18] sought how to move smoothly from the existing controller (λ = 0) to the full MPC strategy (λ = 1).
In a wave energy perspective, CC control guarantees maximum power absorption with the limitation of noncausality such that a prediction for the wave elevation or the excitation force must be known or sufficiently accurately estimated in advance. In practice, CC control based on prediction is often accomplished via a velocity tracking MPC-this has the added benefit of handling constraints. On the other hand, a feedback controller such as a proportional-integral (PI) controller circumvents this issue and requires only signals/sensors which are readily available, but with the lack of constraint-handling capabilities. Hence, in this paper, we incorporate the linear feedback control designs within MPC frameworks to take into account constraints in an optimal manner by revisiting the inverse matching method in Cairano and Bemporad [17]. A new analytical solution is also proposed that immediately generates optimal weight matrices in MPC for a short prediction horizon case. This approach has a number of advantages compared to the current WEC control strategies. First, predictions of future wave elevation or excitation force are not required for power maximization. The MPC behaves exactly as the predesigned linear controller when constraints are inactive, while satisfying the constraints in an optimal manner when they are active. When constraints are active, MPC acts as a supervisory controller to prevent constraint violations. For a WEC, these constraints can be torque limits of the generator or displacement limits of a linear piston or stage. Second, this approach, differently from all the previous approaches for MPC for WECs (see, e.g., [19]), is in a standard MPC form and thus can be implemented by the software already available (e.g., the Model Predictive Control Toolbox from MathWorks). These are significant practical advantages, as the structure of the controller is simpler and implementing such a solver for real-time control is far from trivial.
To investigate the potential of this approach, we consider a single body one degree-of-freedom (1DOF) WEC. The WEC plant model was estimated from real experimental data including power take-off system dynamics and was described in a linear transfer function form with the control force as the input and the velocity of the buoy as the output [20]. Based on this model, an MPC matching problem is solved to yield proper weight matrices. The MPC designed via the proposed approach behaves exactly as a given PI controller when the control system is unconstrained. When constraints are applied, the MPC retains the general impedance matching dynamics, but also avoids constraint violation.

Controller Matching Problem
In this section, an MPC algorithm is derived that will behave exactly as the predesigned state-feedback controller when constraints are inactive. In Cairano and Bemporad [17], it was shown that a linear state-feedback controller can be realized via MPC with specially tuned weight matrices-this approach is referred to as MPC inverse matching. For example, a PI controller is a simple example of state feedback control, so we can apply the inverse matching method to obtain an MPC that exactly behaves as the PI controller when constraints are inactive.
In MPC formulation, the following finite-horizon optimal control problem is solved at each sampling time [13,17] V (x(k), u(k)) = min where the notation denotes the transpose of a vector, such that where Equation (2) is the discrete-time state space equation of the plant model, x ∈ R n is the state vector, u ∈ R m is the control input vector, A and B are the state matrix and the input matrix, respectively, N is the prediction horizon, U(k) = [u (0|k) ... u (N − 1|k)] ∈ R Nm is the vector to be optimized, and V : R n → R 0+ is the value function. The matrices P ∈ R n×n , Q ∈ R n×n , and R ∈ R m×m are the weight matrices that should be tuned. Assuming that the current state x(k) is given, the finite horizon optimal control problem Equation (1) can be recast to the following quadratic program (QP) with respect to U(k) min U(k) U (k)HU(k) + 2x (k)FU(k) (6) such that GU(k) ≤ λ + Λx(k), where G ∈ R q×Nm , λ ∈ R q , and Λ ∈ R q×n are the problem constraints. Also, the following new matrices are defined in Equation (6) where S is the N-steps state reachability matrix, T is the N-steps free evolution matrix and Q ∈ R Nn×Nn , R ∈ R Nm×Nm are block-diagonal matrices When constraints are not active, the QP Eqaution (6) has the unconstrained optimal solution U * (k) which is obtained by solving the Karush-Kuhn-Tucker optimality conditions [21]: where the superscript * denotes the optimal values. The MPC command at step k then picks up only the first move such that where Γ = [I m 0 . . . 0] with the m by m identity matrix I m . Now, we consider a predesigned linear controller in a state feedback form Then, the objective is to find weight matrices P, Q, and R in Equation (1) such that holds, where K is given. In brief, once the weight matrices P, Q, and R are found that satisfy Equation (14) given K, then the MPC with the cost function Equation (1) behaves as the linear controller of the form Equation (13), when the constraints are not active. The weight matrices P, Q, and R used in H and F in Equations (8)- (10) can be calculated using a numerical optimizer such that the error norm of the matrix ΓH −1 F + K is minimized from Equation (14). Cairano and Bemporad [17] reformulated Equation (14) into the following convex problem: min . . .
with a constraint that the Hessian matrix H = R + S QS is positive definite. It is noted that in Equation (16) the elements of K are arbitrary except for the first one because only the first element is picked up for the MPC command, which was not explicitly stated in Cairano and Bemporad [17]. In this paper, the optimization problem Equation (14) is tackled by solving it using a numerical optimizer YALMIP [22] to obtain optimal weight matrices P, Q, and R. Also, the MPC signal is easily generated using the standard MPC block available in the Model Predictive Control Toolbox of MathWorks [23], and hence, there is no need to implement a custom optimizer to solve the MPC problems in real time.
In the MPC Toolbox the KWIK algorithm [24] is used to solve the QP problem and a prediction model is given by a linear function of the current state and the control/excitation input increments [23]. In the next section, as a special case when N = 2 we derive an analytical solution for P, Q, and R.

An Analytical Solution When N = 2
The optimization problem Eqaution (14) is not trivial to solve because of the existence of the inverse matrix H −1 and the non-invertibility of the matrix Γ. However, for systems with a one-step prediction (N = 2) and a single control input (m = 1), it is possible to derive an analytical solution that immediately generates optimal weight matrices P, Q, and R. A short prediction horizon is sometimes useful because it requires less computational loads and the control performance is independent of the prediction horizon in an unconstrained application. Furthermore, in the area of power electronics and drive, almost all applications use a very short horizon due to the requirement for high speed [25,26]. Since this analytical solution enables immediate calculations for the weight matrices, it is extremely useful in applications in which the weight matrices need to be updated in real time.
When N = 2, the matrices S, T , Q, and R in Equations (9) and (10) are simplified as Here, R is a scalar if a single control input (i.e., m = 1) is considered. Accordingly, the matrices H and F in Equation (8) become Then, the inverse of the 2 × 2 matrix H is directly calculated as where the determinant D is given by When m = 1, we have Γ = [1 0] in Equation (12) and should find the matrices Q, R, and P for given A, B, and K that satisfy the following equation: For brevity, consider the following 1 × 2 gain matrix K (which can be easily extended to general 1 × n cases): Then, since any matrices Q, R, and P satisfying Equation (21) can be solutions, let us first simplify Equation (21) by choosing P such that B P = 0. General solutions to B P = 0 are given by where (B ) + is the Moore-Penrose generalized inverse of B and Ξ is an arbitrary 2 × 2 matrix. However, for the use of MPC formulation, P must be positive definite so we select the matrix Ξ as Ξ = γ(I − (B ) + B ) with a positive constant γ such that Although Equation (24) provides an exact solution to B P = 0, the matrix (I − (B ) + B ) is always rank deficient and P in Equation (24) is always positive semi-definite, not positive definite. It is for this reason that one additional term is added to Equation (24): where ρ is a small positive number such that B P is sufficiently close to the zero matrix.
Once the matrix P is obtained by Equation (25), B P = 0 (approximately) holds and the determinant D becomes Then, the main Equation (21) to be solved for Q and R is simplified as where Φ A + BK and Ψ B A PA(A + BK). For the sake of simplicity, let us assume Q to be diagonal and decompose each matrix as and R is a scalar because only one control input is considered. Substituting Equation (28) into Equation (27) yields the following equation: There are three unknowns (R, Q 1 , Q 2 ) and two equations, and hence, there exist infinite number of solutions. For brevity, let us assume Q 2 = ηQ 1 with some constant η and Equation (29) becomes Then, the solution to Equation (30) is easily found as There are two requirements that must be satisfied such that the obtained matrices P, Q, and R are valid solutions to Equation (14): the invertibility of the coefficient matrix in Equation (30) and the positive definiteness of the Hessian matrix H = R + S QS. First, the invertibility of the coefficient matrix in Equation (30) is guaranteed when its determinant is nonzero such that we assume Q 1 = ηQ 2 instead and follow the same procedure with the assumption B 1 ( Hence, we obtain two conditions for the requirement H > 0 as It is noted that −B A PAB is a negative constant since P is positive definite. Then, two free parameters η and γ can be determined that satisfy the two requirements Equations (33) and (35). In general, by picking up a large γ and making −B A PAB negative with a large magnitude, we can easily find η. Hence, we have obtained the explicit solutions for P, Q, and R that are given in Equations (25) and (31), respectively.

Results: Application to a WEC
The MPC matching problem handled in the previous sections can be applied to a WEC control problem, assuming that a successful feedback controller is already in hand. First, we solve the MPC matching problem and verify the solutions via numerical simulations with the estimated WEC model. Next, constraints on the actuator force and the output are considered.
A model-scale heaving point absorber designed at Sandia National Laboratories (for details see [27]) was considered for the case study discussed herein. This device was designed to be 1/17th a full-scale equivalent WEC. Figure 1 and Table 1 show a diagram of the WEC device and its relevant physical parameters.  The WEC device is modeled in a transfer function form as follows: G o (s) = 0.0002371s 3 + 0.01783s 2 − 0.002049s s 4 + 24.7s 3 + 54.53s 2 + 389.2s + 62.92 , which relates the applied actuator control force (input) and the resulting velocity of the WEC device in vertical motion (output). Note that this model does not include the excitation dynamics, which are not necessary for the control approach consider here. (The excitation force model is used only in the simulation of the controller.) To simplify the problem, let us replace the fourth-order model (36) with the following second-order approximation: The simplified model (37) was obtained by the MATLAB command tfest over the frequency range of interest f ∈ [0.1 1.0] (Hz) assuming that the energy in ocean waves lies within a relatively narrow frequency band (periods of 1 ≤ T ≤ 10 s). Note that, as previously stated, the WEC considered here is a model scale device, which is tested in a wave basin. As such, the frequency range of interest is higher than that of an open ocean device. From the second-order model (37), the natural frequency and the damping ratio are found by ω n = √ 16.2 = 4.025 (rad/s) = 0.641 (Hz) and ζ = 1.38/(2ω n ) = 0.171, respectively. Hence, the resonant frequency is located at ω n 1 − 2ζ 2 = 3.905 (rad/s) = 0.622 (Hz), which is an approximate value because of the existence of the zero in the transfer function. At the resonant frequency, the output is in phase with the input, and for small input frequencies, 90 degrees phase lead occurs while 90 degrees phase lag occurs for large input frequencies. Figure 2 shows the Bode diagram of the original model (36) and the simplified model (37), which shows a good agreement between the two models. Hence, we proceed to MPC design based on the simplified model (37). First, we transform Equation (37) to a state-space equation: where x(t) = [z(t) v(t)] is the state vector, z(t) and v(t) are the vertical position and velocity of the WEC device, respectively. The coefficient matrices of this system are The state Equation (38) is described in continuous time but it must be converted to a discrete-time equation for MPC formulation: where ·(k) is a quantity at the time step k and x(k) = [z(k) v(k)] , and where the sampling interval T s = 0.02 s is assumed. Let us assume that the incoming wave is of JONSWAP type spectrum with a peak period of 3.5 s, a significant wave height of 0.254 m, and a peak enhancement factor of 1. Then, with the plant model Equation (37) or Equation (40), the optimal PI controller that maximizes the energy extraction from the wave is given by where the gain matrix is calculated as which extracts maximum energy from the incoming wave. The tuning of the PI controller is suboptimal but very close to optimum [27].

Numerical Simulations
To examine the performance of this controller, we design an MPC that behaves exactly as the PI controller Equation (42) for three different scenarios (see Table 2 for a summary of the scenarios) and evaluate that controller in numerical simulations.

1.
The system has no constraints. 2.
The system has a nearly hard constraint on its position. 3.
The system has a hard constraint on its control input and a soft constraint on its position.  More specifically, in Scenario 2 both OV.minECR and OV.maxECR values (OV denotes the output variable) for the position used in the MPC Toolbox are set as 2 × 10 −4 , where an equal concern for the relaxation (ECR) value is zero for a completely hard constraint and a larger ECR value means a softer constraint. When constraints are imposed, the value function Equation (1) is rewritten as where µ is the weight for constraints and k (≥ 0) is another free optimization variable to be determined. By default, where W i,j (W = P, Q, R) is the (i, j) component of the matrix W [23]. The constraints Equations (3) and (4) change to where V x min , V x max , V u min , and V u max are the ECR parameters, each of which has the command OV.minECR, OV.maxECR, MV.minECR, and MV.maxECR, respectively in the MPC Toolbox. It is easy to see that when an ECR value is zero, the constraint is a hard one that cannot be violated. The larger ECR value is, the softer the constraint is.
In Scenario 3 both MV.minECR and MV.maxECR values (MV denotes the manipulated variable) are set as 0 for the control input and OV.minECR = OV.maxECR = 0.01 are selected for the position. For all scenarios, the control horizon is equal to the prediction horizon N. It should be mentioned that the parameters are scaled to improve numerical accuracy.

Scenario 1: System with No Constraint
First, we consider an unconstrained case and apply the analytical solution with N = 2. Since there is no constraint in the system, the performance will be independent of the selection of N [23]. The matrix P is calculated by (25) with γ = 10 6 and ρ = 10 −6 : P = 999, 899 −10, 050 −10, 050 101 . (48) Next, the matrices Q and R are calculated by (31) with η = 1: Then, the resulting gain matrix −ΓH −1 F is calculated as which is very close to the original gain K. Also, it is easily checked that the product of P and the matrix B in Equation (41) has the norm of the order of 10 −7 , which validates the accuracy of the analytical solution. Figure 3 compares the time histories of the control forces, positions, power absorptions, and cumulative average power obtained by the PI controller and the MPC. The power absorption was defined as the product of the control force and the velocity of the device at each time. Note that the convention of this paper is to show absorbed power as negative power. Because the obtained gain is very similar to the original one, there are only minor differences between the two controllers.  Figure 3 shows the cumulative average power. Again, we can see that the two controllers perform quite similarly.

Scenario 2: System with Nearly Hard Constraint on the Position
In this scenario, a nearly hard constraint (OV.minECR = OV.maxECR = 2 × 10 −4 ) is imposed on the position of the WEC device and the prediction horizon N = 50 is assumed to anticipate constraint violation early enough to allow corrective action. A 'nearly' hard constraint means that the constraint is quite well satisfied on the whole but a small constraint violation is allowed. In practice, a completely hard constraint should not be imposed on the output because the plant is always subject to disturbance and QP infeasibility is inevitable. We apply the numerical solver 'YALMIP' to yield the following weight matrices: Then, the resulting gain matrix −ΓH −1 F is calculated as which is very close to the original gain K. More explicitly, the norm of the matrix difference ΓH −1 F + K is calculated as 0.052. Now, we impose a nearly hard output constraint such that the position should not exceed 0.1 m in both directions. Figure 4 compares the time histories of the control forces, positions, power absorptions, and cumulative average power obtained by the unconstrained PI controller and the MPC with the output (position) constraint. When the constraint becomes active, the control force has a sudden rise to satisfy the constraint and is released immediately after the constraint turns to inactive. As with almost any system, such behavior is undesirable in a WEC, as it may excite high frequency structural responses and create unnecessary loading. This sudden rise can be relaxed by allowing soft constraints or introducing a longer prediction horizon as will be shown in the next subsection. It is shown that the constraint on the position is satisfied quite accurately. Also, as expected, a constrained system extracts less energy due to the limited motion. However, the reduction in power to satisfy the constraint is quite small: the unconstrained PI controller captures the power of −42.74 W on average for one period of the incoming wave while the constrained MPC captures the power of −39.51 W.

Scenario 3: System with Hard Constraint on the Control Force and Soft Constraint on the Position
As the final scenario, we apply a hard constraint (MV.minECR = MV.maxECR = 0) on the control force and a soft constraint (OV.minECR = OV.maxECR = 0.01) on the position. In general, control force and position constraints cannot always be simultaneously satisfied, and so the position constraint is usually softened for QP feasibility [23]. Bacelli and Ringwood [28] extended the discussion of the force and position constraints for WECs and proposed sufficient conditions for the satisfaction of both constraints for a given hydrodynamic model and set of sea conditions. Figure 5 compares the time histories of the control forces, positions, power absorptions, and cumulative average power obtained by the unconstrained PI controller and the MPC with the hard and soft constraints. Three different MPCs with these constraints are shown; these controllers are summarized in Table 3. In addition to the basic MPC ("MPC (N = 50)"), there are two additional MPCs that attempt to reduce the large spikes observed in the control force. The first of these ("MPC (N = 50, |MV| < 50)") limits the slew rate, that is the rate change of the control force. We can see from Figure

Experimental Tests
The obtained MPC algorithms are now applied to real experiments. The experiments were conducted in the US Navy's Maneuvering and Sea Keeping (MASK) basin, operated by the Naval Surface Warfare Center, Carderock Division (NSWCCD) in Bethesda, Maryland. The data collected are publicly available at https://mhkdr.openei.org. For the details about the test device, hardware setup, and the wave basin, refer to Coe et al. [27]. As before, two controllers are considered: PI controller and MPC that approximates the PI controller. The WEC system is disturbed by the same JONSWAP type wave used for the numerical simulations. Two identical wave trains were run: one in which the PI controller was used and one in which the MPC was used. A sampling interval T s = 0.001 s and a prediction horizon N = 2 are used for the controllers and the commanded gains for the PI controller are selected as In order to check the control performance a metric K FIT was calculated by estimating the K I and K P gains back from the data and comparing them with the actual gains of the PI controller used for the test. This metric is defined by whereK I andK P are the integral and proportional gains, respectively, estimated from the data and K I and K P are the commanded gains shown in Equation (53). The gains were estimated by dividing the measured control force by the measured position and velocity and then picking up the constant optimal values that minimize the fitting error. The estimated gains arẽ The error norm between Equations (55) and (53) is calculated as K − K / K = 0.0044. Also, the metric (54) is calculated as 99.36%, which shows good agreement between the two controllers.
The top-most plot in Figure 6 displays the time history of the control signals generated by the PI controller and the MPC, which verifies that the MPC behaves as the PI controller on the whole. In order to ensure that the two controllers (PI and MPC) were influenced by the same external wave, the wave elevation data collected from the PI controller and MPC experiments are plotted in the frequency domain and compared in Figure 7. In the inset, the wave elevation data in time domain are plotted for reference. The wave elevation was measured by a sensor that was located far enough from the buoy so that any radiation effect from the buoy could be ignored. Figure 7 shows that the two controllers were indeed under the quite similar wave forces. In the middle plot in Figure 6, the time history of the velocities of the WEC device obtained by the PI controller and the MPC is depicted. Since the control signal of the MPC closely emulates the PI control signal and the two controllers were influenced by the same wave force, the two velocities are also similar. The lower plot in Figure 6 displays the time history of the power captured by using the PI controller and the MPC. The average power (P) captured by the PI controller and the MPC is −35.21 W and −34.05 W, respectively, which also shows the performance difference of the two controllers is small (3%). Now, we consider constrained MPCs such that the control signal is saturated with given upper/lower limits (assuming hard constraints). It will also be shown that the power captured by the constrained MPC is reduced only by a small amount when compared with the unconstrained MPC even though a very short prediction horizon N = 2 is used. Let us say that U max is the maximum control magnitude that can be reached by the actuator. Since the maximum control magnitude is found to be about 1200 N from the top-most plot in Figure 6, let us constrain the system by forcing U max to be 900, 600, and 300 N.
As with the previous example, experiments were run with the same wave input; in each of the four cases, the maximum control input was set to U max = [1200 900 600 300] (N). The top-most plot in Figure 8 shows the control forces obtained by these unconstrained and constrained MPCs. It is clear that the control forces successfully stay within the desired range. In the middle plot in Figure 8, the resulting velocities are plotted for the unconstrained and constrained MPCs. With the constrained MPCs, the magnitude of the resulting velocity is greater than the unconstrained case because the saturated control force cannot fully regulate the velocity. The lower plot in Figure 8 displays the time history of the power captured by using the unconstrained and constrained MPCs. As expected, when the constraint is active, less power is captured.  However, Table 4 indicates that the power captured by the constrained MPCs is worse only by a small percentage when compared with the unconstrained MPC. In the Table 4, the average power capture obtained using the unconstrained and constrained MPCs is listed and the power loss in percentage is also presented in the parenthesis compared with the unconstrained case. For example, when the control force is constrained with U max = 600 which corresponds to 50% of the maximum control magnitude obtained when the unconstrained MPC is applied, the power capture is −31.91 W, which is worse only by 6.28% compared with the unconstrained case (−34.05 W) despite the fact that a very short horizon N = 2 is employed.

Conclusions
This paper studies the application of MPC tuning methods via inverse matching to a WEC, such that the MPC behaves exactly like a predesigned linear controller when constraints are inactive. With the PI controller and device model obtained from real experimental data, the MPC behaved exactly as the PI controller when constraints are not active. When constraints are active, the MPC can be designed to follow them to the desired degree. A predesigned causal predictionless controller was easily realized via MPC technique with specially tuned weight matrices. This controller requires only signals/sensors that are readily available such as position/velocity. By incorporating it within MPC frameworks, the ability to incorporate constraints is successfully maintained. In addition to numerical simulations, experimental results show both accuracy and feasibility of this method. Also, all the results were readily implemented by the commercial software package (Model Predictive Control Toolbox of MATLAB).
This approach provides an attractive solution for WEC control, as it exploits the advantages of both a feedback-based controller for impedance matching and MPC for constraint handling. Practically speaking, such a controller can be implemented for a WEC within a gain-scheduling regime in order to deal with nonlinearities (i.e., by designing linear controllers and subsequent MPCs for specific sea states). Future work should include MPC matching for a multi-input, multi-output (MIMO) system and application to a MIMO WEC system where multiple devices in an array are coupled electrically with one another.