Simulate path-following control using adaptive model predictive controller
Model Predictive Control Toolbox / Automated Driving
The Path Following Control System block simulates a path-following control (PFC) system that keeps an ego vehicle traveling along the center of a straight or curved road while tracking a set velocity and maintaining a safe distance from a lead vehicle. To do so, the controller adjusts both the longitudinal acceleration and front steering angle of the ego vehicle. The block computes optimal control actions while satisfying safe distance, velocity, acceleration, and steering angle constraints using adaptive model predictive control (MPC).
This block combines the capabilities of the Lane Keeping Assist System and Adaptive Cruise Control System blocks into a single controller.
To customize your controller, for example to use advanced MPC features or modify controller initial conditions, click Create PFC subsystem.
Set velocity
— Ego vehicle velocity setpointEgo vehicle velocity setpoint in m/s. When there is no lead vehicle, the controller tracks this velocity.
Time gap
— Safe time gapSafe time gap in seconds between the lead vehicle and the ego vehicle. This time gap is used to calculate the minimum safe following distance constraint. For more information, see Safe Following Distance.
Relative distance
— Distance between lead vehicle and ego vehicleDistance in meters between lead vehicle and ego vehicle. To calculate this signal, subtract the ego vehicle position from the lead vehicle position.
Relative velocity
— Velocity difference between lead vehicle and ego vehicleVelocity difference in meters per second between lead vehicle and ego vehicle. To calculate this signal, subtract the ego vehicle velocity from the lead vehicle velocity.
Longitudinal velocity
— Ego vehicle velocityEgo vehicle velocity in m/s.
Curvature
— Road curvatureRoad curvature, specified as 1/R, where R is the radius of the curve in meters.
The road curvature is:
Positive when the road curves toward the positive Y axis of the global coordinate system.
Negative when the road curves toward the negative Y axis of the global coordinate system.
Zero for a straight road.
The controller models the road curvature as a measured disturbance with previewing. You can specify the curvature as a:
Scalar signal — Specify the curvature for the current control interval. The controller uses this curvature value across the prediction horizon.
Vector signal with length less than or equal to the Prediction Horizon — Specify the current and predicted curvature values across the prediction horizon. If the length of the vector is less than the prediction horizon, then the controller uses the final curvature value in the vector for the remainder of the prediction horizon.
Lateral deviation
— Ego vehicle lateral deviationEgo vehicle lateral deviation in meters from the centerline of the lane.
Relative yaw angle
— Angle from lane centerlineEgo vehicle longitudinal axis angle in radians from the centerline of the lane.
Minimum longitudinal acceleration
— Minimum ego vehicle accelerationMinimum ego vehicle longitudinal acceleration constraint in m/s2. Use this input port when the minimum acceleration varies at run time.
To enable this port, select Use external source for the Minimum longitudinal acceleration parameter.
Maximum longitudinal acceleration
— Maximum ego vehicle accelerationMaximum ego vehicle longitudinal acceleration constraint in m/s2. Use this input port when the maximum acceleration varies at run time.
To enable this port, select Use external source for the Maximum longitudinal acceleration parameter.
Minimum steering angle
— Minimum front steering angleMinimum front steering angle constraint in radians. Use this input port when the minimum steering angle varies at run time.
To enable this port, select Use external source for the Minimum steering angle parameter.
Maximum steering angle
— Maximum front steering angleMaximum front steering angle constraint in radians. Use this input port when the maximum steering angle varies at run time.
To enable this port, select Use external source for the Maximum steering angle parameter.
Enable optimization
— Controller optimization enable signalController optimization enable signal. When this signal is:
Nonzero, the controller performs optimization calculations and generates the Longitudinal acceleration and Steering angle control signals.
Zero, the controller does not perform optimization calculations. In this case, the Longitudinal acceleration and Steering angle output signals remain at the values they had when the optimization was disabled. The controller continues to update its internal state estimates.
To enable this port, select the Use external signal to enable or disable optimization parameter.
External control signal
— Control signals applied to ego vehicleActual control signals applied to the ego vehicle. The first element of this signal is the longitudinal acceleration in m/s2, and the second element is the steering angle in radians. The controller uses these signals to estimate the ego vehicle model states. Use this input port when the control signals applied to the ego vehicle do not match the optimal control signals computed by the model predictive controller. This mismatch can occur when, for example:
The Path Following Control System is not the active controller. Maintaining an accurate state estimate when the controller is not active prevents bumps in the control signals when the controller becomes active.
The steering or acceleration actuator fails and does not provide the correct control signal to the ego vehicle.
To enable this port, select the Use external control signal for bumpless transfer between PFC and other controllers parameter.
Vehicle dynamics matrix A
— State matrix of ego vehicle predictive modelState matrix of ego vehicle predictive model. The number of rows in the state matrix corresponds to the number of states in the predictive model. This matrix must be square.
The ego vehicle predictive model defined by Vehicle dynamics matrix A, Vehicle dynamics matrix B, and Vehicle dynamics matrix C must be minimal.
To enable this port, select the Use vehicle model parameter.
Vehicle dynamics matrix B
— Input-to-state matrix of ego vehicle predictive modelInput-to-state matrix of ego vehicle predictive model. The number of rows in this signal must match the number of rows in Vehicle dynamics matrix A.
The ego vehicle predictive model defined by Vehicle dynamics matrix A, Vehicle dynamics matrix B, and Vehicle dynamics matrix C must be minimal.
To enable this port, select the Use vehicle model parameter.
Vehicle dynamics matrix C
— State-to-output matrix of ego vehicle predictive modelState-to-output matrix of ego vehicle predictive model. The number of columns in this signal must match the number of rows in Vehicle dynamics matrix A.
The ego vehicle predictive model defined by Vehicle dynamics matrix A, Vehicle dynamics matrix B, and Vehicle dynamics matrix C must be minimal.
To enable this port, select the Use vehicle model parameter.
Longitudinal acceleration
— Acceleration control signalAcceleration control signal in m/s2 generated by the controller.
Steering angle
— Front steering angle control signalFront steering angle control signal in radians generated by the controller. The front steering angle is the angle of the front tires from the longitudinal axis of the vehicle. The steering angle is positive towards the positive lateral axis of the ego vehicle.
Use vehicle parameters
— Define ego vehicle model using vehicle propertieson
(default) | off
Select this parameter to define the ego vehicle model used by the MPC controller by specifying properties of the ego vehicle. The ego vehicle model is the linear model from the longitudinal acceleration and front steering angle to the longitudinal velocity, lateral velocity, and yaw angle rate.
To define the vehicle model, specify the following block parameters:
Total mass
Yaw moment of inertia
Longitudinal distance from center of gravity to front tires
Longitudinal distance from center of gravity to rear tires
Cornering stiffness of front tires
Cornering stiffness of rear tires
Longitudinal acceleration tracking time constant
For more information on the ego vehicle model, see Ego Vehicle Predictive Model
Selecting this parameter clears the Use vehicle model parameter.
Use vehicle model
— Define ego vehicle model using state-space matricesoff
(default) | on
Select this parameter to define the state-space matrices of the ego vehicle model used by the MPC controller. The ego vehicle model is the linear model from the longitudinal acceleration and front steering angle to the longitudinal velocity, lateral velocity, and yaw angle rate.
To define the initial internal model, specify the A, B, and C state-space matrices. The internal model must be a minimal realization with no direct feedthrough, and the dimensions of A, B, and C must be consistent.
Typically, the ego vehicle model is velocity-dependent, and therefore, it varies over time. To update the internal model at run time, use the Vehicle dynamics A, Vehicle dynamics B, and Vehicle dynamics C input ports.
For more information on the ego vehicle model, see Ego Vehicle Predictive Model
Selecting this parameter clears the Use vehicle parameters parameter.
Total mass
— Ego vehicle mass1575
(default) | positive scalarEgo vehicle mass in kg.
To enable this parameter, select the Use vehicle parameters parameter.
Yaw moment of inertia
— Moment of inertia about the ego vehicle vertical axis2875
(default) | positive scalarMoment of inertia about the ego vehicle vertical axis in mNs2.
To enable this parameter, select the Use vehicle parameters parameter.
Longitudinal distance from center of gravity to front tires
— Distance from the ego vehicle center of mass to its front tires1.2
(default) | positive scalarDistance from the ego vehicle center of mass to its front tires in meters, measured along the longitudinal axis of the vehicle.
To enable this parameter, select the Use vehicle parameters parameter.
Longitudinal distance from center of gravity to rear tires
— Distance from the ego vehicle center of mass to its rear tires1.6
(default) | positive scalarDistance from the ego vehicle center of mass to its rear tires in meters, measured along the longitudinal axis of the vehicle.
To enable this parameter, select the Use vehicle parameters parameter.
Cornering stiffness of front tires
— Front tire stiffness19000
(default) | positive scalarFront tire stiffness in N/rad, defined as the relationship between the side force on the front tires and the angle of the tires to the longitudinal axis of the vehicle.
To enable this parameter, select the Use vehicle parameters parameter.
Cornering stiffness of rear tires
— Rear tire stiffness33000
(default) | positive scalarRear tire stiffness in N/rad, defined as the relationship between the side force on the rear tires and the angle of the tires to the longitudinal axis of the vehicle.
To enable this parameter, select the Use vehicle parameters parameter.
Longitudinal acceleration tracking time constant
— Time constant for acceleration tracking0.5
(default) | positive scalarTime constant for tracking longitudinal acceleration, specified in seconds.
To enable this parameter, select the Use vehicle parameters parameter.
A
— Initial state matrix of ego vehicle predictive modelInitial state matrix of ego vehicle predictive model. The number of rows in the state matrix corresponds to the number of states in the predictive model. This matrix must be square.
The initial ego vehicle predictive model defined by A, B, and C must be minimal.
Typically, the ego vehicle model varies over time. To update the state matrix at run time, use the Vehicle dynamics A input port.
To enable this parameter, select the Use vehicle model parameter.
B
— Initial input-to-state matrix of ego vehicle predictive modelInitial input-to-state matrix of ego vehicle predictive model. The number of rows in this parameter must match the number of rows in A.
The initial ego vehicle predictive model defined by A, B, and C must be minimal.
Typically, the ego vehicle model varies over time. To update the input-to-state matrix at run time, use the Vehicle dynamics B input port.
To enable this parameter, select the Use vehicle model parameter.
C
— Initial state-to-output matrix of ego vehicle predictive modelInitial state-to-output matrix of ego vehicle predictive model. The number of columns in this parameter must match the number of rows in A.
The initial ego vehicle predictive model defined by A, B, and C must be minimal.
Typically, the ego vehicle model varies over time. To update the state-to-output matrix at run time, use the Vehicle dynamics C input port.
To enable this parameter, select the Use vehicle model parameter.
Initial longitudinal velocity
— Initial velocity of the ego vehicle model15
(default) | nonnegative scalarInitial velocity of the ego vehicle model in m/s, which can differ from the actual ego vehicle initial velocity.
This value is used to configure the initial conditions of the model predictive controller. For more information, see Initial Conditions.
Note
A very small initial velocity, for example eps
, can produce a
nonminimal realization for the controller plant model, causing an error. To prevent
this error, set the initial velocity to a larger value, for example
1e-3
.
Transport lag between model inputs and outputs
— Total transport lag in ego vehicle model0
(default) | nonnegative scalarTotal transport lag, τ, in the ego vehicle model in seconds. This lag includes actuator, sensor, and communication lags. For each input-output channel, the transport lag model is:
Maintain safe distance between lead vehicle and ego vehicle
— Enable spacing controlon
(default) | off
To configure the safe following distance, set the Default spacing parameter. For more information on the safe following distance used by the controller, see Safe Following Distance.
Default spacing
— Minimum spacing to lead vehicle10
(default) | nonnegative scalarMinimum spacing in meters between the lead vehicle and the ego vehicle. This value corresponds to the target relative distance between the ego and lead vehicles when the ego vehicle velocity is zero.
This value is used to calculate the:
Minimum safe following distance. For more information, see Safe Following Distance.
Controller initial conditions. For more information, see Initial Conditions.
To enable this parameter, select the Maintain safe distance between lead vehicle and ego vehicle parameter.
Minimum steering angle
— Minimum front steering angle-0.26
(default) | scalar between -pi/2
and pi/2
Minimum front steering angle constraint in radians.
If the minimum steering angle varies over time, add the Minimum steering angle input port to the block by selecting Use external source.
This parameter must be less than the Maximum steering angle parameter.
Maximum steering angle
— Maximum front steering angle0.26
(default) | scalar between -pi/2
and pi/2
Maximum front steering angle constraint in radians.
If the maximum steering angle varies over time, add the Maximum steering angle input port to the block by selecting Use external source.
This parameter must be greater than the Minimum steering angle parameter.
Minimum longitudinal acceleration
— Minimum ego vehicle acceleration-3
(default) | scalarMinimum ego vehicle longitudinal acceleration constraint in m/s2.
If the minimum acceleration varies over time, add the Minimum longitudinal acceleration input port to the block by selecting Use external source.
Maximum longitudinal acceleration
— Maximum ego vehicle acceleration2
(default) | scalarMaximum ego vehicle longitudinal acceleration constraint in m/s2.
If the maximum acceleration varies over time, add the Maximum longitudinal acceleration input port to the block by selecting Use external source.
Sample time
— Controller sample time0.1
(default) | positive scalarController sample time in seconds.
Prediction horizon
— Controller prediction horizon10
(default) | positive integerController prediction horizon steps. The controller prediction time is the product of the sample time and the prediction horizon.
Control horizon
— Controller control horizon3
(default) | positive integer | vector of positive integersController control horizon, specified as one of the following:
Positive integer less than or equal to the Prediction horizon parameter. In this case, the controller computes m free control moves occurring at times k through k+m-1, and holds the controller output constant for the remaining prediction horizon steps from k+m through k+p-1. Here, k is the current control interval.
Vector of positive integers, [m1, m2, …], where the sum of the integers equals the Prediction horizon parameter. In this case, the controller computes M blocks of free moves, where M is the length of the control horizon vector. The first free move applies to times k through k+m1-1, the second free move applies from time k+m1 through k+m1+m2-1, and so on. Using block moves can improve the robustness of your controller.
Weight on velocity tracking
— Tuning weight for longitudinal velocity tracking0.1
(default) | positive scalarTuning weight for longitudinal velocity tracking. To produce smaller velocity-tracking errors, increase this weight.
Weight on lateral error
— Tuning weight for lateral error1
(default) | positive scalarTuning weight for lateral error. To produce smaller lateral errors, increase this weight.
Weight on change of longitudinal acceleration
— Tuning weight for change in longitudinal acceleration0.1
(default) | positive scalarTuning weight for changes in longitudinal acceleration. To produce less-aggressive vehicle acceleration, increase this weight.
Weight on change of steering angle
— Tuning weight for change in steering angle0.1
(default) | positive scalarTuning weight for changes in steering angle. To produce less-aggressive steering angle changes, increase this weight.
Use suboptimal solution
— Apply suboptimal solution after specified number of iterationsoff
(default) | on
Configure the controller to apply a suboptimal solution after a specified maximum number of iterations, which guarantees the worst-case execution time for your controller.
For more information, see Suboptimal QP Solution.
After selecting this parameter, specify the Maximum iteration number parameter.
Maximum iteration number
— Maximum optimization iterations10
(default) | positive integerMaximum number of controller optimization iterations.
To enable this parameter, select the Use suboptimal solution parameter.
Use external signal to enable or disable optimization
— Add port for enabling optimizationoff
(default) | on
To add the Enable optimization input port to the block, select this parameter.
Use external signal for bumpless transfer between PFC and other controllers
— Add external control signal input portoff
(default) | on
To add the External control signal input port to the block, select this parameter.
Create PFC subsystem
— Create custom controllerGenerate a custom PFC subsystem, which you can modify for your application. The configuration data for the custom controller is exported to the MATLAB® workspace as a structure.
You can modify the custom controller subsystem to:
Modify default MPC settings or use advanced MPC features.
Modify the default controller initial conditions.
Use different application settings, such as a custom safe following distance definition.
The default ego vehicle predictive model for path-following control is the combination of two state-space models, one for adaptive cruise control and one for lane keeping.
The predictive state-space model for adaptive cruise control is:
Here, τ is the Longitudinal acceleration tracking time constant parameter.
The input to this model is the longitudinal acceleration in m/s2, and the output is the longitudinal velocity in meters per second.
The predictive state-space model for lane keeping is:
Here:
VX is the longitudinal velocity of the car. At the start of the simulation, this velocity is equal to the Initial condition for longitudinal velocity parameter. At run time, this velocity is equal to the Longitudinal velocity input signal.
m is the Total mass parameter.
IZ is the Yaw moment of inertia parameter.
LF is the Longitudinal distance from center of gravity to front tires parameter.
LR is the Longitudinal distance from center of gravity to rear tires parameter.
CF is the Cornering stiffness of front tires parameter.
CR is the Cornering stiffness of rear tires parameter.
The input to this model is the steering angle in radians. The outputs are the lateral velocity in meters per second and yaw angle rate in radians per second.
The Path Following Control System block combines these models as follows:
The inputs to this combined model are the longitudinal acceleration in m/s2 and steering angle in radians. The outputs are the longitudinal velocity in meters per second, lateral velocity in meters per second, and yaw angle rate in radians per second.
The controller creates its internal predictive model by augmenting the ego vehicle dynamic model. The augmented model includes the road curvature as a measured disturbance input signal.
To define a different ego vehicle predictive model, select the Use vehicle model parameter, and specify the initial state-space model. Then, specify the run-time values of the state-space matrices using the Vehicle dynamics A, Vehicle dynamics B, and Vehicle dynamics C input signals.
When the Maintain safe distance between lead vehicle and ego vehicle parameter is selected, the model predictive controller computes the safe following distance constraint; that is, the minimum relative distance between the lead and ego vehicle, as:
Here:
DS is the Default spacing parameter.
GT is the Time gap input signal.
VE is the Longitudinal velocity input signal.
To define a different safe following distance constraint, create a custom path-following control system by, on the Block tab, clicking Create PFC subsystem.
By default, the model predictive controller assumes the following initial conditions for the ego vehicle:
Longitudinal velocity is equal to the Initial longitudinal velocity parameter.
Longitudinal acceleration is zero.
Lateral velocity is zero.
Steering angle is zero.
Yaw angle rate is zero.
When the Maintain safe distance between lead vehicle and ego vehicle parameter is selected, the controller assumes the following additional initial conditions:
The lead vehicle longitudinal velocity is equal to the Initial longitudinal velocity parameter.
Relative distance between the lead vehicle and ego vehicle is:
Here:
DS is the Default spacing parameter.
GT is the time gap and is assumed to
be 1.4
.
VE is the Initial longitudinal velocity parameter.
If the initial conditions in your model do not match these conditions, the Steering angle and Longitudinal acceleration outputs can exhibit initial bumps at the start of the simulation.
To modify the controller initial conditions to match your simulation, create a custom path-following control system by, on the Block tab, clicking Create PFC subsystem.