LTV-MPC Approach for Automated Vehicle Path Following at the Limit of Handling

In this paper, a linear time-varying model predictive controller (LTV-MPC) is proposed for automated vehicle path-following applications. In the field of path following, the application of nonlinear MPCs is becoming more common; however, the major disadvantage of this algorithm is the high computational cost. During this research, the authors propose two methods to reduce the nonlinear terms: one is a novel method to define the path-following problem by transforming the path according to the actual state of the vehicle, while the other one is the application of a successive linearization technique to generate the state–space representation of the vehicle used for state prediction by the MPC. Furthermore, the dynamic effect of the steering system is examined as well by modeling the steering dynamics with a first-order lag. Using the proposed method, the necessary segment of the predefined path is transformed, the linearized model of the vehicle is calculated, and the optimal steering control vector is calculated for a finite horizon at every timestep. The longitudinal dynamics of the vehicle are controlled separately from the lateral dynamics by a PI cruise controller. The performance of the controller is evaluated and the effect of the steering model is examined as well.


Introduction
Collision avoidance systems, advanced driver-assistance systems, and the many other types of automated vehicle functions are becoming more and more popular and have become the most related topics in the field of automotive research. The advantages of automatization of different vehicle functions include the potential to improve road safety, reduce pollutant emissions and traveling times, and eliminate human errors, which are the primary cause of accidents. An automated vehicle can avoid collisions by using the steering and the braking system conventionally; furthermore, a new alternative solution for accident prevention is to force the vehicle into an unstable state, in which a control software can drive the vehicle at a level as high as a professional driver, for example, in case of a drift maneuver [1][2][3]. Each of these studies aims to drive a vehicle at the level of a professional human driver. In [1], an LQ controller is used for steady-state drifting, while in [2], an MPC is applied for drifting in varying road surface conditions, and in [3], the drift is realized by a reinforcement learning algorithm.
According to the structure of the automated vehicles, the motion planning and motion control modules are essential to achieve these goals. In this paper, the authors focus on the motion control layer by proposing an LTV-MPC structure for path-following tasks. Based on the existing research and the experimental tests, the authors identify three features to be considered in a path-following controller at higher vehicle speed, the inclusion of the vehicle dynamics in the control law, and the knowledge of the path ahead of the vehicle. Furthermore, the inclusion of the steering dynamics in the control law is identified as the third required feature, which has a high effect on the performance of the path-following controller. Numerous path-following solutions are summarized in [4][5][6], e.g., geometric controllers, such as Pure Pursuit and Stanley, linear quadratic regulators, model predictive controllers, neural network based controllers, etc. The authors decided to apply MPC for the path-following tasks, as the identified features can be incorporated into an MPC and the MPC handles constraints well.
Three main types of MPC methods are used in the field of path tracking, the linear parameter-varying MPC (LPV-MPC) [7][8][9] the linear time-varying MPC (LTV-MPC) [10][11][12][13], and the nonlinear MPC (NLMPC) [10,[14][15][16] solutions. The LPV-MPC applies a linear vehicle model for state prediction and the structure of that model does not change over time; however, a few of the variables, e.g., the velocity of the vehicle, can change. The model is calculated at each timestep based on the current states. The main advantage of the LPV-MPC method is the low computational cost; however, this technique reaches the limit of applicability when the controlled system leaves the linear region, where the system model is defined. Thus, the LPV-MPC method is unable to handle nonlinear system dynamics in a wide range, e.g., tire dynamics in the nonlinear region of the tires. The change in given parameters can be determined by measurements or estimated by state estimators.
The LTV-MPC [10][11][12]17] can handle the nonlinear behavior of the system by linearizing the nonlinear system at every timestep and estimating the states accordingly. The predicted states are still based on a linear system model; however, these are extremely close to the real states of the system, since the controller calculates the predicted states for a short horizon, e.g., a few hundred ms, in which range the prediction is acceptable. In the LTV-MPC structure, the current vehicle state is the basis of the linearization and an additional transformation in the system matrices and the control input is required. The LTV-MPC method claims a bit more computational cost than the LPV-MPC; however, it still can run on an ECU or a rapid prototyping unit in real-time. When using LPV-MPC or LTV-MPC, usually, a quadratic cost function is formulated, the optimum of which is computationally cheaper to find.
Applying the NLMPC, the evolution matrices are calculated based on a nonlinear system model without linearization. A nonlinear cost function is formulated, the optimum of which is computationally expensive to find, and might limit the applicability of the NLMPC, as discussed in [6].
In this paper, the authors propose an LTV-MPC method for automated vehicle pathfollowing scenarios. The aim of this research is to apply a nonlinear vehicle model for state prediction, handled by a successive linearization technique, where the vehicle model is linearized at every timestep by Jacobi linearization. The applied vehicle model is coupled with a Pacejka tire model, the parameters of which are identified based on measurement results. The application of a nonlinear vehicle model coupled with the Pacejka model is expected to be suitable for handling the high sideslip angles of the tires, where the tires are operating in the nonlinear region-which the LPV-MPC is unable to handle-hence, the vehicle should be able to perform the path-following task at the limit of handling.
In the existing research, if the path is defined in a forward-looking way, it is performed by defining the reference path in the global coordinate system by a series of (x,y) points [10,12,[18][19][20][21]. While this solution is obvious, it incorporates high nonlinearity into the plant model, since the motion of the vehicle needs to be transformed from the vehicle's coordinate system-i.e., the ego frame-to the global frame. A novel method for the definition of the reference path is presented, by which the number of the nonlinear terms involved in the entire control problem is successfully reduced. According to our proposal, the path is transformed from the global coordinate system into the vehicle's coordinate system and the reference states, which need to be realized by the vehicle, are expressed according to the state-space representation of the vehicle plant model used for state prediction by the MPC. In the existing research, the reference path is either not defined in a forward-looking way [7,22], or the forward-looking property is limited to the lateral displacement of the vehicle [8,12]. In this article, the authors define both the lateral displacement and the heading angle as a reference series for a finite horizon, as presented in Section 3.
Some studies considered the effect of the actuator fault [23] and the robustness of the controlled system [24]; in this article, the effect of the steering system is analyzed. The effect of the consideration of the steering dynamics in the vehicle model applied for state prediction is discussed in this paper as well. The steering system is modeled as a first-order lag. Moreover, the performance of the controller is evaluated with and without the inclusion of the steering dynamics in the plant model by applying the control input to the same vehicle model, which includes the steering dynamics.
The contribution of this paper is twofold, (1) The proposed LTV-MPC method for path following includes a novel reference definition method that effectively reduces the nonlinear terms of the path-following problem by transforming the path from the global frame to the ego frame. By this transformation, the nonlinearity that needs to be managed during the Jacobian linearization is significantly reduced as the state vector has two elements less, since the position of the vehicle in the global coordinate system is not considered in the plant model. This approach removes two equations from the model while retaining the lookahead principle of the reference definition. (2) The analysis of the effect of the steering dynamics on the path-following problem are incorporated in an LTV-MPC.
The structure of this paper is as follows. Section 2 presents the vehicle model applied in the MPC for state prediction, regarding the steering dynamics and the vehicle model used for testing the controller. In Section 3, the generation of the reference path is presented. The structure of the proposed MPC and the derivation of the quadratic programming (QP) cost function is presented in Section 4. Section 5 describes the results, while Section 6 provides an analysis, in which the controller is tested on a sine wave path and the effect of the application of linear and nonlinear vehicle models in the state prediction is compared. Finally, Section 7 discusses the concluding remarks and suggestions for further research directions.

Vehicle Modelling
In this section, the applied tire and vehicle models are presented. The authors use different vehicle models for state prediction and for simulation testing. Furthermore, a Pacejka tire model is applied in both models, the parameters of which are identified from measurement data.

Vehicle Model for Testing
A four-wheel vehicle model is used for testing the path-following controller in a simulation environment. The roll and pitch dynamics are neglected in the vehicle model since they have no significant effect on the path-following problem. The model describes the planar dynamics of the vehicle: the angular acceleration around the vertical axis (1), the longitudinal (2), and the lateral acceleration (3).
where r is the yaw rate, I z is the moment of inertia of the vehicle around the axis z, M z is the resultant torque, which rotates the vehicle around axis z, a x is the longitudinal acceleration, m is the mass of the vehicle, F x V is the resultant of the longitudinal forces, acting on the vehicle, V y is the lateral velocity of the vehicle in the ego coordinate system, a y is the lateral acceleration, F y V is the resultant of the lateral forces acting on the vehicle, and V x is the longitudinal velocity in the ego frame.
The resultant torque and forces are calculated by (4), (5), and (6) where a and b are the distances between the center of gravity (C.G.) and the front and the rear axle, respectively, s f and s r are the front and rear track of the vehicle, respectively, δ FL and δ FR are the front-left and front-right steering angles of the individual wheels, F xij and F yij are the longitudinal and lateral forces at the individual wheels, where i marks the front and the rear axles, hence, i = R as rear or i = F as front, and j marks the left and the right wheels, hence, j = L as left or j = R as right. In this model, the authors consider a rear-wheel drive vehicle, with equally distributed traction force on the rear wheels, hence, F xRL and F xRR are equal. Furthermore, the authors assume a front-wheel steering vehicle; hence, the steering angle is solely interpreted at the front wheels.
In Figure 1, the tire sideslip angles and the sideslip angle of the vehicle at C.G. are also shown. The vehicle sideslip angle is calculated by (7). The sideslip angle of an individual wheel is calculated using the velocity vector of the given wheel and the steering angle at the front axle. The velocity vector of the wheels is described by (8), the position vector of the wheels is defined by (9), where the origin is the C.G. of the vehicle, and the tire sideslip angles are described by (10) and (11).
where i = FL ∨ FR ∨ RL ∨ RR and α denotes the tire sideslip angle. The horizontal velocity V z in (8) is considered zero as the model focuses solely on the planar dynamics of the vehicle, as stated previously. A Pacejka tire model [20] is applied to calculate the lateral forces at each wheel. During this research, the vehicle moves at constant speed; thus, the longitudinal forces are small, especially compared to the lateral forces and, hence, solely a lateral tire model is used in this article. The Pacejka tire model is described by Equations (12) and (13) [25].
where B is the stiffness factor, C is the shape factor, D is the peak factor, E is the curvature factor, S v and S h are the vertical and the horizontal shift, respectively, and the α sideslip angle is defined in degrees. The Pacejka tire model and the meaning of the factor variables are described in more detail in [25]. The entire vehicle model, including the tire model, is fitted to a test vehicle used for automated vehicle function tests, such as automated drift, Moose test, and other pathfollowing tasks. Further details of the test vehicle setup can be found in [1].
To determine the characteristics of the tires, a ramp steer maneuver is taken by the test vehicle and the necessary variables are measured with the data acquisition system. The parameters in the Pacejka tire model are fitted to the measurement results, as shown in Figure 2. The characteristics of the rear tires have a greater slope at small sideslip angles, where the tire model is nearly linear. That means the rear tires have greater cornering stiffness than the front tires. This meets the expectations of the authors, since the front wheels of the vehicle are mounted with 245/35 R19 tires, while 265/35 R19 tires are applied at the rear wheels-using the same sidewall height, the wider tire becomes stiffer. During the identification of the tire characteristics, the vertical load of the tires is assumed to be constant.
The identified tire models describe the characteristics of the tires located on the same axle, not the individual wheels; hence, the results in Figure 2 contain the lateral forces generated by the two front tire pairs and the two rear tire pairs, respectively. Thus, when modeling one wheel, the value of the lateral force needs to be halved.

Vehicle Model for State Prediction
The states of the vehicle model applied in the MPC for state prediction are advisable to choose in a way to correspond to the control purpose. In this case, the aim is to drive the vehicle along a predefined path as fast as the vehicle is still able to execute the maneuver. The applied vehicle model shown in Figure 3 is a dynamic bicycle model, coupled with different Pacejka tire models at the front and the rear wheels. Since the bicycle model applies one wheel per axle, the identified Pacejka models can be used without halving the values of the lateral forces of the tire characteristics. The applied vehicle model describes the dynamics of the vehicle with four equations, enhanced by a fifth equation for the steering dynamics. The steering dynamics are considered to take into account the dynamic lag of the steering system, in which way the controller considers that the demanded steering angle will only be realized with a time delay. The dynamics of the steering system are modeled by a first-order lag, which has one tuning parameter, the time constant, to identify the measurements made. The time constant is determined by measurements conducted on the test vehicle. The state vector of the system is chosen as x = [v V y ϕ r δ act ] T , where y is the lateral displacement of the vehicle in the ego frame, ϕ is the heading angle of the vehicle, and δ act is the actual steering angle. The vehicle dynamics are described by (14)- (17), while the steering dynamics by (18). . .
where δ dem is the steering angle required by the MPC and T st is the time constant of the first-order steering model. To calculate the F yf and F yr lateral forces, the identified Pacejka models are applied by determining the parameters of (12) and (13). The sideslip angles are calculated by (19) and (20).
where v yF and v yR are the velocity of the front and the rear wheels and α F and α R are the sideslip angles of the front and the rear wheels, respectively. The tests are conducted without the consideration of the steering dynamics, in which case, the state vector is chosen as (18) is neglected during the state prediction.
The future states of a vehicle using a vehicle model can be predicted for a finite time horizon. In this article, the authors apply successive linearization to continuously generate the state-space representation of the vehicle. The successive linearization allows linearizing the nonlinear vehicle model at the current operating point, calculating the state-space representation accordingly, and making the state prediction based on the linearized system. The linearization is conducted at every timestep; hence, the MPC can use the latest state of the vehicle as a basis of the state prediction.
The linearization of the vehicle model is solved by Jacobian linearization and leads to the continuous-time [A c , B c , C c , D c ] state-space representation of the system (21).
where F is the system of nonlinear equations (14)- (18), x is the state vector, u is the control input, which, in this case, is the demanded steering angle δ dem , while C c and D c are considered constant matrices (22).
Furthermore, matrix B c remains constant during the linearization in the following form (23).
The partial derivatives of the matrices are evaluated at the desired operation point, which is specified by the state vector x o , the time derivative of the state vector .
x o , and the control vector u o and, thereby, the partial derivatives in (21) need to be calculated based on the initial conditions x o and u = u o . The evolution of the continuous-time system can be described by (24) .
where x L and u L are the linearized state of the system and the control input, respectively. The constant terms of (24) are incorporated into K c , with which the continuous-time representation of the state-space system becomes (25) and (26).
where y L is the output of the system. Since the MPC is a discrete-time-control technique, the state-space representation needs to be discretized to get the discrete-time state-space representation of the system where T s is the discretization timestep, which is equal to the timestep value of the MPC. Finally, the discrete-time representation of the system (31) can be used for state prediction by the MPC.

Reference Path Definition
As previously introduced, in numerous studies, the control problem is to minimize the error between the spatial reference path and the position of the vehicle. In this research, the state-space representation of the vehicle contains the position of the vehicle in the global frame, according to the reference definition, where the path is defined by a series of (x,y) points in the global frame. The controlled state is the spatial position of the vehicle-X and Y coordinates in the global frame. The spatial path-following problem is possibly supplemented with the tracking of the heading and the yaw rate. The problem with the formulation when the X and Y coordinates are defined as a reference is that it introduces several nonlinear terms into the vehicle model applied for state prediction by the MPC and increases the dimension of the state vector by two. Therefore, if the control problem is defined to follow spatial points, the vehicle needs to be transformed to the global frame in the state-space representation-this transformation is responsible for higher nonlinearity. Thus, if the increase in the dimension of state vector x by two-displacement in X and Y directions-can be omitted, then the computational requirements of the MPC can be reduced, which increases directly proportionally with the increase in the dimension of state vector x, as presented in [26]. In this article, a novel method is proposed for the reduction in nonlinear terms, which is to transform the path to the vehicle and describe the path-following task in a way that retains the advantageous property of the MPC, which is the knowledge of the path ahead of the vehicle-the knowledge of the reference for a finite horizon. In this article, the lateral displacement and the heading of the vehicle are defined as a reference to be followed, in accordance with the states of the vehicle model used for state prediction. The reference state variables are chosen as they clearly define the relationship between the vehicle and the path, with respect to the lateral terms. To calculate the reference for a finite horizon, firstly, the lateral and angular errors are defined, as shown in Figure 4. The lateral error e is the distance between the C.G. point of the vehicle and the intersection point M on the y axis of the ego frame and the reference path. The angular error α is the angle between the heading of the vehicle and the tangent of the path at the intersection point M. To generate the reference for the lateral displacement state y and for the heading angle state ϕ, the vehicle, the reference path, the ego frame, and the global frame need to be considered, as shown in Figure 5. In Figure 5, the lateral errors and the angle errors are shown for a finite N p horizon, where N p is the prediction horizon of the MPC. If the lateral error values from the ego coordinate system can be seen, these errors can be interpreted as reference lateral displacements, which corresponds to y, the first state of the state-space representation. In a similar way, the angular errors can be interpreted as ϕ reference heading angles. By applying this method for reference generation, the transformation of the vehicle coordinates from the ego frame to the global frame becomes avoidable, resulting in omitting several nonlinear terms from the vehicle model. As stated before, if the vehicle model is not transformed into the global frame, the path needs to be transformed into the ego frame of the vehicle. The transformation is defined by Equations (32)-(34), where γ is the transformation angle, e 1 and α 1 are the lateral and angular errors at the M point, respectively, H is an offset vector, which shifts the entire reference trajectory to the necessary lateral displacement, P i is the (x,y) coordinate of the ith path point described in the global frame, and P ref,i is the coordinate of the reference path point in the ego frame. After the transformation, the lateral error references and the angular error references can be calculated; the lateral error reference is the y-direction component of P ref,i , which determines the series of e 1 , e 2 , . . . e Np , while the angular reference can be calculated considering the tangent of the transformed path at every point where the lateral displacements are defined. After the reference values are calculated, the reference vector can be generated by (35) and (36).
While the original reference path is described by N p points, the reference lateral displacements and the heading angles need to be calculated for every path point. The result is a stacked matrix T with a dimension of 2× N p . Using this method of reference definition, several nonlinear terms are omitted from the vehicle model, the reference is coherent to the states of the vehicle model, the reference lateral displacement is coupled with the first state in the state-space model, and the angle reference is coupled with the third state. Furthermore, the presented reference generation method preserves the forward-looking nature of the problem definition, which is necessary for the path-following task.

MPC Structure
As stated in Section 2, the authors apply a successive linearization technique, similar to [12,14,17], to handle the nonlinear vehicle model. During the successive linearization, the nonlinear vehicle model is linearized at every timestep based on the actual state vector of the vehicle. The resulting state-space representation is applied to predict the future vehicle states for a finite prediction horizon, the length of which is N p . Using an MPC controller, the objective is to minimize the difference between the given reference states in the system and the predicted states in the systems by calculating the optimal control input vector as a result of an optimization process, while meeting a set of constraints. The optimization process requires a cost function to minimize. In this section, the derivation of the cost function is presented. The cost function depends on the tracking error and the amount of the control input. The tracking error e(k) in the k-th timestep is defined as e(k) = y(k) − r(k), where y(k) is the current state of the system and r(k) is the given reference. Then, the evolution of the error can be defined for the k-th, k + 1-th, and k + 2-th timesteps, etc., by (37).
The evolution of the error needs to be calculated for the entire prediction horizon, resulting in the error vector e ∈ R N p N o , where N o is the number of outputs in the system, which is equal to the rows of matrix C d . In this case, N o = 2 corresponds to the lateral displacement and the heading angle references.
where u ∈ R N p N u is the control input vector, i.e., the result of the optimization process, N u is the number of the control inputs, which is equal to the columns of matrix B d -in this case-N u = 1 corresponding to the steering input, As the evolution of the tracking error is calculated, the cost function can be defined as In this article, the authors apply a quadratic cost function (41), the optimum of which is found by the built-in numerical solver in the MATLAB software package, called quadprog.

Substituting (40) into (41) results in
The following task is to aggregate the terms that do not depend on the control input u and aggregate the quadratic and linear terms of (44), which leads to (45).
A more compact form of the cost function is in (46).
When solving optimal control problems, the aim is to minimize the tracking errors and the amount of control input. However, a steady-state error might occur while tracking, since the cost of the steady-state error might be less than the cost of the higher value of the control inputs. To avoid the steady-state error, the cost function is defined by the increments in the current control input, rather than defining the individual control inputs. Applying this solution, the effect of the control input values on the cost function is minimal. The control input vector can be transformed as which can be written as where = D ∈ R N p N u ×N p N u is a lower diagonal matrix built of I ∈ R N u ×N u matrices and zero matrices of the same dimensions. Performing the transformation, the evolution of the tracking error can be written as where u(k − 1) is constant for the entire prediction horizon, thus, it can be included in the constant term = K. The new const function, the solution of which is the optimal series of control input increments, is given in the following form where the matrices are defined as

Results
A reference path is defined for testing the controller shown in Figure 6. The path contains a double-lane-change section, which is considered an evasive maneuver [27], and a U-turn section, which is built by a clothoid segment, marked with A and B in Figure 6. As the vehicle drives along the U-turn section, the curvature of the path is the linear function of the length of the arc. During the simulations, the vehicle is started from the left side of the path, using a 0.2 m offset, which is applied to examine how the vehicle can find the path. The simulations are conducted at 50, 60, and 70 km/h, using vehicle models, excluding and including the steering dynamics; hence, in total, six cases are examined. During the simulations, the identified model of the steering system is applied to the vehicle model in every case; hence, the consideration of the steering dynamics in the prediction model is expected to lead to better results. The value of N p and N c is 10 and the sampling time is 0.05 s.
The results are summarized in Table 1, where e is the lateral error and φ is the orientation error. As the velocity of the vehicle is increased, the average and maximal errors are also increased; furthermore, the application of the steering dynamics can reduce the errors significantly.
Please note that the figures below sometimes do not show the difference between the reference and the real steering angle, which is due to the accurate reference tracking. As shown in Figures 7 and 8, at 50 km/h, the consideration of the steering dynamics can increase the accuracy of the controller and result in a smoother steering command, which increases the stability of the vehicle and provides better ride comfort to the passengers. Furthermore, smaller steering interventions result in smaller sideslip angles.   The result at 60 km/h is shown in Figures 9 and 10. The results show similar behavior to the 50 km/h cases; however, the errors are a little higher, although still in an acceptable range. As shown in Figure 2, the tires become saturated when reaching about 4 deg and 2 deg of sideslip at the front and the rear wheels, respectively.
The tire sideslips reach about 7 deg at the front wheels and 5 deg at the rear wheels, which means the tires operate in a saturated state when the vehicle drives at the ending section of the U turn, resulting in a higher lateral error there.
The results at 70 km/h are shown in Figures 11-14. Figures 12 and 14 show the same results as Figures 11 and 13, respectively, but in a smaller time window to make the results more visible. The vehicle reaches and exceeds the friction limit at the U turn, which leads to high lateral errors-the vehicle is physically not able to follow the path due to its great curvature.    The importance of the 70 km/h case is to show whether the controller can handle the vehicle at the friction limit or increase the steering angle to a high value, applying unnecessarily large sideslip angles, which can not decrease the path-following errors. According to the simulation results, the controller can handle the nonlinear characteristics of the tires and operates stably, even when the sideslip angles reach high values. As shown in Figures 11 and 13, the sideslip angles reach a maximal 8 deg and 5 deg at the front and the rear wheels, respectively. The further increase in the sideslip would not generate a larger lateral force. As the cost function penalizes the value of the steering intervention, the steering angle does not increase. In this research, the authors do not apply constraints on the sideslip angle; however, the controller does not generate too-large values due to the inclusion of the nonlinear Pacejka tire model in the plant model.
In the double-lane-change section, the errors do not show a significant increase compared to the lower-speed scenarios; the controller could drive the vehicle faster in this section, but the higher speed would result in a greater error at the U-turn. Please consider that in Table 1, the high values of the average errors at 70 km/h are caused by the large errors at the U-turn, which has a serious effect when calculating the average errors.
Overall, the inclusion of the steering dynamics in the plant model leads to a more accurate path following, with a smoother steering intervention, while the complexity of the model is not seriously increased. Furthermore, as shown by the fourth subgraphs in Figures 7-14, the inclusion of the steering dynamics in the plant model results in a more accurate steering demand tracking, where the tracking error is decreased in every case.

Evaluation of the Effectiveness of the Controller and Comparison of the Linear/Nonlinear Plant Model
In this section, the effectiveness of the controller is evaluated by testing on a sine wave path. Furthermore, the effect of the vehicle plant model is also evaluated by comparing the results of a linear vehicle model and a nonlinear vehicle model. The reference path is shown in Figure 15. The distance between the cones is 30 m, which means 60 m wavelength, and the lateral peak value of the sine wave is 2.5 m in both directions. Each test scenario presented in this section is conducted using the four-wheel vehicle model presented in Section 2.1 and the plant model contains the steering dynamics in each test case. The results of the scenarios conducted on the sine wave path are summarized in Table 2. Please note that when the linear plant model is applied in the 70 km/h scenario, the results are calculated solely for the time interval of 0-4.2 s, since at 4.2 s, the vehicle left the path. In the first scenario, the controller presented in Section 4 is tested, using the vehicle model presented in Section 2.2 for state prediction. The results at 70 km/h vehicle velocity are shown in Figure 16. The controller drives the vehicle accurately and stably, the maximum sideslip angles are 4 deg, and the small values in the lateral and the angular errors also demonstrate the accuracy of the controller. In the next scenario, the nonlinear vehicle plant model is replaced by a linear one, in which small angle assumptions are applied during the linearization of the model, and the vehicle model is coupled with a linear tire model, while the model of the steering dynamics remains unchanged. The resulting linear vehicle model is described by Equations (51)-(55), . .
where the lateral force is a linear function of the sideslip angle (56).
In (56), c α f and c αr are the cornering stiffness of the front and rear tires, respectively, which are identified from the measurement results (Figure 2), and the state vector remains The first test using the linear plant model is conducted at 60 km/h and the results are shown in Figure 17. Both the lateral and angular errors are large, while the controller is unable to drive the vehicle along the path when the velocity is increased to 70 km/h, as shown in Figure 18. In the 70 km/h scenario, the vehicle is unable to follow the path anymore, resulting in leaving the path. The reason for the poor performance of the linear-model-based control is the neglection of the nonlinearities of the controlled vehicle in the state prediction. Both the nonlinear terms in the vehicle model and the tire model are neglected, which results in a less accurate state prediction; hence, in the 60 km/h scenario, performance degradation of the controller is realized, and at 70 km/h, the controller is unable to drive the vehicle along the path.
The results clearly show that the application of the nonlinear vehicle model coupled with the nonlinear tire model for state prediction has a significant advantage over the linear vehicle model and linear tire model couple. Using the nonlinear model, the controller is able to drive the vehicle stable at higher speeds and more accurately than with the linear model. Furthermore, the effectiveness of the proposed LTV-MPC method is also confirmed in this section by testing the controller in the sine wave path.

Conclusions
In this paper, an LTV-MPC structure is proposed for path-following scenarios. The controller uses a nonlinear vehicle model coupled with a Pacejka tire model for state prediction. Furthermore, the dynamics of the steering system are also incorporated into the prediction model. The reference definition is embedded into the LTV-MPC, corresponding to the states in the plant model, which results in a simplified plant model, removing two equations from the model, which are responsible for the coordinate transformation. The presented reference definition method can reduce the nonlinear terms in the path-following problem definition by transforming the path into the ego frame, resulting in a simplified description of the vehicle dynamics, the nonlinearities of which are easier to handle. Using this method, the transformation from the ego frame to the global frame is practically conducted separately from the vehicle model.
During the simulation tests, the inclusion of the steering dynamics results in a more accurate path-following performance with smaller errors and smoother steering command demanded by the controller. Furthermore, the proposed reference definition is proven to be effective, while the nonlinear terms in the plant model are reduced. The proposed controller structure coupled with the plant model can handle the vehicle at sharp corners when the tires operate in the nonlinear region. The vehicle drives stably during the tests, while the errors remain low. According to the results, the controller is applicable for even emergency scenarios, e.g., for a double-lane change. The effectiveness of the proposed controller is tested in a sine wave path, where the controller is proven to be as accurate as on the double-lane-change path. Furthermore, the effect of the linear and nonlinear plant model is also analyzed and the advantage of applying the nonlinear plant model is clearly proven on the sine wave path.
Regarding future research opportunities, the vehicle plant model and the model of the steering system, as an actuator, could be further detailed, which could result in more accurate state prediction. The dynamics of the steering system could be modeled as a function of the vehicle velocity, the yaw rate, and the lateral acceleration, and the self-aligning torque of the tires could be incorporated into the plant model as well. Furthermore, the proposed reference definition may reduce the computational requirements of an NLMPC, as the plant model has fewer nonlinear terms.