---
title: 'Feedback-MPPI: Fast Sampling-Based MPC'
url: https://www.emergentmind.com/topics/feedback-mppi
type: topic
---

# Feedback-MPPI: Fast Sampling-Based MPC

Searching arXiv for recent papers on Feedback-MPPI and closely related MPPI feedback variants.
Feedback-MPPI, abbreviated F-MPPI, is a sampling-based model predictive controller that augments standard Model Predictive Path Integral control with a local linear feedback law computed from the sensitivity of the MPPI solution to the current state. In the specific formulation introduced in “Feedback-MPPI: Fast Sampling-Based MPC via Rollout Differentiation -- Adios low-level controllers,” the purpose of this augmentation is to permit rapid closed-loop corrections around the current state without requiring full re-optimization at each timestep, thereby retaining MPPI’s flexibility for nonlinear dynamics and non-convex costs while adding a Riccati-like inner-loop feedback mechanism [2506.14855]. In a broader usage within the MPPI literature, “feedback” can also denote the receding-horizon closed-loop character of MPPI itself; however, F-MPPI names a more specific architecture in which rollout differentiation supplies an explicit local gain matrix rather than relying only on repeated horizon replanning [2309.13201].

## 1. Standard MPPI background

F-MPPI is defined relative to the standard MPPI formulation for finite-horizon optimal control. The problem is written as
\[
\begin{aligned}
\theta^\ast = \arg\min \;& \ell_N(x_N) + \sum_{i=0}^{N-1} \ell_i(x_i,u_i) \\
\text{s.t. } & x_0 = \hat x, \\
& x_{i+1} = f(x_i,u_i), \\
& u_i = \pi_\theta(x_i,t_i).
\end{aligned}
\]
Here \(x_i\) is the state, \(u_i\) is the control, \(\hat x\) is the current measured state, \(\ell_i\) are running costs, \(\ell_N\) is a terminal cost, \(\theta\) are sampled control parameters, and \(\pi_\theta\) is a parameterized input trajectory policy. The policy parametrization may use direct control sampling, splines, or other low-dimensional parameterizations [2506.14855].

MPPI samples around a nominal parameter vector \(\bar\theta\),
\[
\theta^k = \bar\theta + \Delta^k, \qquad \theta^k \sim \mathcal N(\bar\theta,\Sigma),
\]
rolls each sample through the nonlinear dynamics, and evaluates
\[
J^k := \ell_N(x_N^k) + \sum_{i=0}^{N-1} \ell_i(x_i^k,u_i^k).
\]
The resulting parameter update is a weighted average of perturbations,
\[
\theta^\ast = \bar\theta + \sum_{k=1}^{K} \omega^k \Delta^k,
\]
with importance weights
\[
\omega^k = \frac{\mu^k}{\sum_{j=1}^{K}\mu^j}, \qquad
\mu^k = \exp\!\left(-\frac{1}{\lambda}(J^k-\rho)\right), \qquad
\rho = \min_k J^k.
\]
The control actually applied is
\[
u^\ast = \pi_{\theta^\ast}(\hat x,t_0).
\]
In standard MPPI, this \(u^\ast\) is applied directly and typically held constant until the next optimization cycle. That coupling between optimization frequency and actuation frequency is the limitation F-MPPI is designed to relax [2506.14855].

## 2. Defining feature: local feedback around the MPPI solution

The central addition in F-MPPI is a local linear approximation of the optimal MPPI policy with respect to the initial state,
\[
K := \frac{\partial u^\ast}{\partial x_0}.
\]
This matrix is used as a feedback gain between MPPI solves:
\[
u_{\text{fb}} = u^\ast + K\,(\hat x - x_{\rm sp}),
\]
where \(x_{\rm sp}\) is a local setpoint associated with the MPPI solution. The formulation does not simply reuse \(x_0\) as the setpoint. Instead, the setpoint is continuously updated by linearly interpolating between \(x_0\) and \(x_1\) so that inner-loop tracking is smoother [2506.14855].

This yields a two-rate architecture. The outer loop runs standard MPPI at a moderate frequency and produces a nominal sequence and control action. The inner loop runs at a higher rate and applies the linear feedback correction continuously until the next MPPI update. The paper’s description is explicit: the expensive stochastic horizon optimization is performed less often, while the high-bandwidth stabilization is handled by the linear feedback correction. This suggests a division of labor in which MPPI remains the mechanism for global nonlinear and non-convex trajectory optimization, whereas the gain \(K\) supplies local corrective authority around the current trajectory [2506.14855].

The intended consequence is high-frequency closed-loop action without re-solving the entire sampling problem every timestep. The paper presents this as especially relevant for fast systems where planning at \(50\) Hz is insufficient but actuation at \(200\)–\(500\) Hz is required [2506.14855].

## 3. Rollout differentiation and gain derivation

The feedback gain is obtained by differentiating the MPPI solution with respect to the initial state \(x_0\). Because
\[
u^\ast = \pi_{\theta^\ast}(x_0,t_0),
\]
the chain rule gives
\[
\frac{\partial u^\ast}{\partial x_0}
=
\frac{\partial \pi}{\partial \theta}\frac{\partial \theta^\ast}{\partial x_0}
+
\frac{\partial \pi}{\partial x_0}.
\]

The derivative of the MPPI parameter update follows from
\[
\theta^\ast = \bar\theta + \sum_{k=1}^{K}\omega^k \Delta^k,
\]
which yields
\[
\frac{\partial \theta^\ast}{\partial x_0}
=
\sum_{k=1}^{K} \Delta^k \frac{\partial \omega^k}{\partial x_0}.
\]
The key sensitivity identity is the derivative of the normalized importance weights:
\[
\frac{\partial \omega^k}{\partial x_0}
=
\frac{\omega^k}{\lambda}
\left(
\frac{\partial J^k}{\partial x_0}
-
\sum_{j=1}^{K}\omega^j \frac{\partial J^j}{\partial x_0}
\right).
\]
Substitution into the chain rule produces the general gain expression
\[
K =
\sum_{k=1}^{K}
\frac{\partial \pi}{\partial \theta}\Delta^k
\frac{\omega^k}{\lambda}
\left(
\frac{\partial J^k}{\partial x_0}
-
\sum_{j=1}^{K}\omega^j \frac{\partial J^j}{\partial x_0}
\right)
+
\frac{\partial \pi}{\partial x_0}.
\]
For direct zero-order sampling with
\[
\theta = (u_0,\dots,u_{N-1}),
\]
the parametrization Jacobian reduces to a selector, and the concrete implementation becomes
\[
K =
\sum_{k=1}^{K}
\Delta_0^k
\frac{\omega^k}{\lambda}
\left(
\frac{\partial J^k}{\partial x_0}
-
\sum_{j=1}^{K}\omega^j \frac{\partial J^j}{\partial x_0}
\right).
\]
The gain therefore depends on rollout cost gradients with respect to the initial state rather than on a separate backward dynamic-programming pass [2506.14855].

The paper characterizes this construction as Riccati-inspired rather than Riccati-derived. In a toy validation on a double integrator, it compares MPPI-derived gains to the infinite-horizon LQR gain obtained from the discrete algebraic Riccati equation with \(Q=\mathrm{diag}(1,1)\) and \(R=1\), reporting
\[
K_{\rm LQR} = (-0.958,\,-1.707).
\]
As the number of samples increases, the MPPI gain converges toward the LQR gain; even with fewer samples, the resulting closed loop remains stable. A plausible implication is that the differentiable sampling law can recover familiar local feedback structure in regimes where linear-quadratic approximations are meaningful, while still being embedded in a nonlinear sampling-based MPC stack [2506.14855].

## 4. Computational structure, differentiability, and control-theoretic interpretation

F-MPPI changes the computational profile of MPPI without replacing the MPPI optimizer itself. The gain computation requires differentiating rollout costs \(\partial J^k/\partial x_0\), but this operation is parallelizable because each rollout is independent. The reported implementation uses JAX, runs on CUDA-compatible GPUs, computes rollouts and gradients in parallel, and exploits automatic differentiation. The paper states that the gain computation increases runtime by about \(40\)–\(70\%\) relative to plain MPPI, depending on horizon and sample count, but remains practical because it does not require a new full optimization at the inner-loop frequency [2506.14855].

The framework also imposes a differentiability requirement. Because the gain uses derivatives of rollout costs, local differentiability is needed. The paper states that this is not a serious limitation in robotics because many models are piecewise smooth. It also notes a more specific caveat for constraints represented by hard indicator penalties: such terms are locally constant and therefore do not affect the gain. A smoother barrier formulation is identified as a better alternative, but this was left for future work [2506.14855].

These points delimit the method’s position within control theory. F-MPPI is not a fully gradient-based nonlinear MPC solver, and it does not explicitly solve a Riccati equation. It is instead a sampling-based MPC method supplemented by a local sensitivity-based linear correction. That distinction matters because the feedback law inherits both the strengths and limitations of the rollout distribution and cost model. This suggests that the controller’s local behavior is strongly shaped by the quality of trajectory-cost derivatives, even though the nominal plan is still produced by importance-weighted stochastic sampling [2506.14855].

## 5. Experimental platforms and reported performance

The paper validates F-MPPI on two robotic platforms: a quadrupedal robot in simulation and a quadrotor in simulation and on hardware. In the quadruped study, the platform is a Unitree Aliengo simulated in MuJoCo with a simplified Single Rigid Body Dynamics model and control inputs given by ground reaction forces. The task is to track desired linear velocity references between \(0\) and \(0.5\) m/s over randomly generated rough terrain while withstanding random disturbances of \(\pm 5\) Nm. The cost is
\[
\ell_{\rm track} = \|x-x_{\rm track}\|^2_{Q_{\rm track}} + \|u-u_{\rm reg}\|^2_{R_{\rm reg}}.
\]
The experiments use \(5000\) samples, horizon \(N=10\), \(\delta_t = 0.02\) s, F-MPPI at \(50\) Hz, and baseline MPPI at \(50\), \(80\), and \(100\) Hz. The reported result is that F-MPPI at \(50\) Hz remains competitive with, and in some metrics better than, MPPI running at \(100\) Hz, while the inner-loop feedback keeps the ground reaction forces updated at \(500\) Hz [2506.14855].

In the quadrotor study, the platform is a MikroKopter quadrotor with onboard compute on an Nvidia Jetson Orin NX, state estimation via Telekyb3 and an Unscented Kalman Filter, and simulation in Gazebo. The task is to reach a goal in \(3\)D and, in simulation, navigate around circular obstacles. Rotor speeds \(w_1,\dots,w_4\) are sampled with cubic splines. The main configuration uses horizon \(N=15\), \(K=800\) samples, F-MPPI planning at \(50\) Hz, and inner-loop updates at \(200\) Hz. The stage cost includes
\[
\ell_{\rm goal} = \|x-x_{\rm goal}\|^2_{Q_r} + \|v\|^2_{Q_v} + \|q-q_0\|^2_{Q_q} + \|\omega\|^2_{Q_\omega},
\]
together with
\[
\ell_u = \|u-u_h\|^2_R,
\]
and obstacle avoidance is represented by
\[
\ell_{\rm obs} = Q_{\rm obs}\,\mathbf{1}_{x\in \chi_{\rm obs}}.
\]

The reported outcomes distinguish simulation and real flight. In simulation, the robot successfully navigates through obstacles and reaches the goal without getting stuck in local minima. In real flight, F-MPPI outperforms standard MPPI even when MPPI runs at the highest feasible hardware frequency of \(66.7\) Hz. The paper reports post-transient RMSE values of \(0.017\) m in \(x\) and \(0.027\) m in \(z\) for F-MPPI, versus \(0.048\) m in \(x\) and \(0.038\) m in \(z\) for standard MPPI, corresponding to a \(64.6\%\) reduction in \(x\)-RMSE and a \(28.9\%\) reduction in \(z\)-RMSE in the experiment figure caption [2506.14855].

## 6. Position within the MPPI literature

The term “feedback” has more than one use in the MPPI literature, and the distinction is important. Some papers explicitly state that standard MPPI is already a feedback controller in the model predictive control sense because the state is re-estimated and the optimization rerun at each control step. That interpretation appears in output-sampled MPPI and in safety-oriented MPPI-DBaS, where the feedback character is identified with receding-horizon closed-loop replanning rather than with an added inner-loop gain law [2309.13201]. A related control-theoretic treatment further describes MPPI as a sample-based receding-horizon control method whose feedback behavior is induced by repeated recomputation from the current measured state, and analyzes deterministic MPPI suboptimality as exploration variance is annealed to zero [2502.20953].

Other recent MPPI variants introduce explicit corrective layers, but of different types. Manifold-Constrained MPPI is described as a two-level feedback control architecture in which latent-space MPPI produces a near-feasible reference and an execution-level quadratic program corrects residual manifold mismatch in a single solve rather than through iterative projection [2605.24813]. RL-guided whole-body MPPI places MPPI on top of a pretrained RL locomotion policy as an add-on feedback controller, with the RL policy used as a sampling prior and MPPI continuously correcting that prior online through modular cost terms [2606.25123]. CSC-MPPI is presented as a constrained, feedback-oriented extension that first shifts sampled controls into the feasible region through a primal-dual gradient method and then uses DBSCAN to select representative control inputs instead of taking a global weighted average [2506.16386].

The resulting landscape can be summarized succinctly:

| Method | Feedback mechanism | Distinctive feature |
|---|---|---|
| F-MPPI | Local linear gain from rollout differentiation | High-rate inner-loop correction [2506.14855] |
| MC-MPPI | Execution-level QP correction | Hard manifold-constraint enforcement [2605.24813] |
| RGB | MPPI as add-on feedback above RL prior | Modular objective correction without retraining [2606.25123] |
| CSC-MPPI | Receding-horizon control with feasible-sample correction and cluster selection | Primal-dual constraint handling and DBSCAN [2506.16386] |

A common misconception is therefore that “Feedback-MPPI” names any closed-loop MPPI controller. The literature shows two usages. One is generic, referring to MPPI’s ordinary receding-horizon feedback structure. The other, and more specific, refers to F-MPPI as a controller that differentiates the MPPI solution to obtain an explicit state-feedback gain [2506.14855].

## 7. Significance, scope, and limitations

F-MPPI’s contribution is best understood as a feedback layer on top of MPPI rather than as a new stochastic optimizer. It keeps the standard MPPI sampling rule, importance weights, and nominal-control update, and it adds a local linear gain derived from the state sensitivity of the sampled solution. This allows high-frequency stabilization and tracking between expensive MPPI updates while preserving the flexibility of MPPI for nonlinear, non-convex, and constrained robotic problems [2506.14855].

The method’s scope is bounded by the assumptions of differentiable rollout costs and the quality of the local linear approximation. Because indicator-style obstacle penalties do not contribute to the gain, the inner-loop correction can be partially blind to hard nonsmooth constraints. This suggests that F-MPPI is particularly well matched to settings where a smooth local correction around the nominal MPPI trajectory is valuable, but where the nominal planner still carries responsibility for global obstacle avoidance and long-horizon nonlinear behavior. The empirical results on a quadruped and a quadrotor indicate improved smoothness, stability, and tracking accuracy under those conditions [2506.14855].

Within the broader MPPI literature, Feedback-MPPI marks a shift from viewing sampling-based MPC purely as repeated open-loop sequence optimization toward treating the sampled solution itself as an object that can be differentiated and turned into a fast corrective controller. That is the sense in which the method bridges sampling-based planning and Riccati-style feedback without abandoning the stochastic rollout framework that defines MPPI [2506.14855].

Source: https://www.emergentmind.com/topics/feedback-mppi