REVIEW 4 major objections 4 minor 19 references
Kernelized offset-free data-driven predictive control for nonlinear systems
T0 review · 4 major / 4 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read The paper proves that a kernelized velocity-form data-driven predictive controller can achieve offset-free tracking for unknown nonlinear systems, with recursive feasibility and Lyapunov stability guarantees under an exactness assumption.
desk verdict Useful data-driven MPC extension with a clean least-squares learning step, but the stability proof has a load-bearing gap that needs repair before the guarantees can be trusted. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the velocity form (incremental) state-space model, where the state is $z_k = \mathrm{col}(y_{k-1}, \Delta x_k)$ and the dynamics matrices $A$, $B$, $C$ depend nonlinearly on scheduling variables. The paper's key move is to parameterize each unknown gradient matrix ($\nabla_x f$, $\nabla_u f$, $\nabla_x h$) as a product of a constant coefficient matrix and a kernel basis matrix built from the data. Because these matrices enter linearly, the representer theorem applies and the coefficients are obtained by least squares. This kernelized velocity model has the same dimension as the analytic model, and the prediction matrices $\Psi$ and $\Gamma$ are formed by products of these kernel matrices, so the MPC problem can be solved with sequential QP by fixing the scheduling variables and iterating.
What would settle it
For the pendulum example, recompute the terminal set $Z_T$ from the kernelized model and test whether $A_{cl}(\rho) Z_T \subseteq Z_T$, $K Z_T \subseteq \Delta U$, and the LMI of Assumption 4.2 hold at a dense grid of points inside $Z_T$; any violation would invalidate the recursive feasibility proof, since the paper does not verify Assumption 4.2 over the whole set.
Extended reading notes
Core claim
The central claim is that the structure of the analytic velocity form can be preserved exactly when functions are replaced by kernel-based representations, because the unknown gradient matrices appear linearly. This turns model learning into two independent least-squares problems: one for the output-gradient coefficients and one for the state/input-gradient coefficients. The learned model has the same state dimension as the analytic velocity model, so the online predictive control problem remains a parameterized nonlinear program solvable by sequential QP. The paper then states Assumption 4.2, a terminal cost and terminal set condition involving a locally stabilizing feedback and an LMI, and proves that under exact kernel representation this yields recursive feasibility and asymptotic Lyapunov stability of the closed loop. Simulations on a discretized pendulum show offset-free tracking for piecewise-constant references under piecewise-constant disturbances, with performance essentially identical to the analytic velocity-form NMPC.
Load-bearing premise
The kernelized velocity model exactly represents the true velocity dynamics; if the finite-sample kernel model is only approximate, the paper's recursive feasibility and asymptotic stability guarantees do not apply.
Editorial extensions
If this is right
- Any nonlinear system whose velocity-form gradient functions lie in the RKHS of the chosen kernel can be controlled offset-free from noiseless data alone, without solving a nonlinear identification problem.
- The terminal cost and terminal set conditions (Assumption 4.2) apply to both analytic and data-driven velocity-form MPC, replacing the more conservative terminal equality constraint used in earlier velocity-form NMPC.
- Because the kernelized model has the same dimension as the analytic model, the online computation cost of vKDPC is comparable to that of model-based velocity-form NMPC; in the pendulum example it averaged 0.0683 s per sample versus 0.0438 s.
- The learning step reduces to solving two least-squares problems, and the paper reports 2.0 s to fit a 2000-sample dataset for the pendulum, so the approach scales linearly with data size in the coefficient matrices.
- When the kernel model is inexact, the guarantees degrade from asymptotic stability to input-to-state stability, and recursive feasibility can be recovered by leaving the initial extended state as a free optimization variable, per the robust framework referenced in the paper.
Reading between the lines
- A natural next test is whether the least-squares kernel model remains accurate when the data are noisy; the paper assumes noiseless data, and the offset-free property in that setting is not established.
- The same structural preservation could be applied to other basis functions, e.g., polynomials or neural features, where the linear-in-parameters property holds, potentially giving offset-free data-driven MPC beyond kernels.
- The proof of recursive feasibility depends on the kernel model being exact; a practical user would need to verify Assumption 4.2's LMI over the entire terminal set, or rely on the robust variant, before deploying vKDPC safely.
- Because the scheduling variables $\rho$ are fixed during each QP step, the convergence of Algorithm 1 to a local optimum of Problem 4.1 is not guaranteed; a reader could test whether the 4\text{--}5 average iterations reach the same cost as the analytic vNMPC on harder nonlinearities.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a data-driven predictive control scheme for nonlinear systems based on a kernelized velocity-form state-space model. The unknown gradients in the velocity model are parameterized with kernel functions and learned by least squares, and the resulting controller is formulated as a nonlinear program solved by sequential quadratic programming. The authors state terminal cost and set conditions intended to guarantee recursive feasibility and Lyapunov stability, and they illustrate the method on a discretized pendulum example, comparing it with an analytic velocity-form NMPC.
Significance. The main idea, preserving the velocity-model structure so that a fixed-dimensional kernelized predictor can be learned by a closed-form least-squares problem, is attractive, and the simulation results are encouraging. The least-squares derivation in Lemma 3.3 is clean for noiseless data, and the authors are explicit about the exactness assumption on the learned model. If the theoretical gaps in Section IV were repaired, the paper would be a useful contribution to data-driven nonlinear MPC. As it stands, however, the recursive feasibility and stability claims are not established as written, and the numerical example does not verify the assumptions under which those claims are proven.
major comments (4)
- [Section IV, Theorem 4.3 proof] The proof of Theorem 4.3 defines the shifted sequence tilde u_{k+1} in absolute control inputs, but the vKDPC problem in (15) optimizes over incremental inputs and the prediction matrices (13)-(14) act on Delta u. The displayed derivation then applies the convolution formula to absolute inputs, which is inconsistent with the model. The correct candidate is the incremental sequence {Delta u*_{1|k}, ..., Delta u*_{N-1|k}, K(z*_{N|k}-r)}. Consequently, recursive feasibility of Problem 4.1 at time k+1 is not established by the proof as written.
- [Section IV, Assumption 4.2 and Theorem 4.3] Assumption 4.2 and the terminal step of the proof of Theorem 4.3 are not stated in consistent coordinates. The terminal constraint is (z_{N|k}-r) in Z_T, so Z_T is a set of deviations, while the proof concludes that tilde z_{N|k+1} = A_cl(rho) z*_{N|k} lies in Z_T. Since A(rho)r = r for the velocity model (2), one must instead show that A_cl(rho)(z*_{N|k}-r) lies in Z_T. The proof omits the subtraction of r and therefore does not establish the terminal constraint. Relatedly, the inclusion Z_T subset Z in Assumption 4.2 is not dimensionally consistent if Z_T is a set of deviations and Z is a set of absolute states; the condition should be {r+e : e in Z_T} subset Z, or Z should be redefined in error coordinates.
- [Section IV, exactness assumption] Before Assumption 4.2, the paper assumes that the kernelized velocity dynamics hat A, hat B, hat C are an exact approximation of the true velocity dynamics. For a finite sample of noiseless data from an unknown system, this exactness assumption is generally false. Theorems 4.3 and 4.4 are proved only for this exact representation, and Remark 4.5 merely defers the inexact case to robust methods in [18]. The simulation in Section V does not verify exactness or the terminal LMI over the computed terminal set, so the theoretical guarantees do not apply to the controller actually implemented in the example. The paper should either provide a quantified approximation bound and adapt the stability analysis, or clearly present the numerical results as purely empirical.
- [Section V and Remark 4.6] In Remark 4.6 and Section V, the terminal ingredients P, K, and Z_T are computed from the linearized dynamics at the reference, i.e., (hat A(x_r,u_r)+hat B(x_r,u_r)K). Assumption 4.2, however, requires the nonlinear, state-dependent LMI and set invariance to hold for all e in Z_T. No argument is given that the computed terminal set satisfies the nonlinear conditions, for example via Lipschitz bounds or a grid verification. This gap is load-bearing because the recursive feasibility and stability theorems rely on Assumption 4.2.
minor comments (4)
- [Lemma 3.3] The full rank hypothesis is ill-posed: the matrices [Ky_1 ... Ky_{s-1}] and [Kx_1 ... Kx_{s-1}] have dimensions (n s) x (s-1) and ((n+m)s) x (s-1), so they cannot have full row rank for s>1; if full column rank is meant, the least-squares solution is not unique and the pseudoinverse formulas give only the minimum-norm minimizer. Please clarify.
- [Theorem 4.4 proof] In the proof of Theorem 4.4, the expression (z*_{N|k}-r)(Q+K^T R K)(z*_{N|k}-r) should be a quadratic form with transposes: (z*_{N|k}-r)^T (Q+K^T R K) (z*_{N|k}-r).
- [Remark 3.2] The claim of equal computational complexity should be qualified: the optimization variables have the same dimensions as in the analytic velocity model, but evaluating the kernel matrices has cost O(s) per entry.
- [General presentation] There are several typos, e.g., 'dyanmics' in Remark 4.6 and 'a offset-free' in the abstract; also, Fig. 2 would benefit from labeled axes and a quantitative validation metric for the kernel model.
Circularity Check
No significant circularity: the kernelized model is learned from data and benchmarked against an independent analytic vNMPC controller; stability claims are conditional on standard terminal assumptions.
full rationale
The derivation chain is self-contained rather than circular. The kernelized velocity model is obtained by applying the representer theorem and solving least-squares problems (Lemma 3.3) on noiseless data, and the multi-step predictor is validated against separate test data (Fig. 1). The central closed-loop claim is an external comparison with an analytic velocity-form NMPC controller on the pendulum benchmark, so the reported offset-free performance is not forced by the fitted coefficients. The recursive-feasibility and stability results are conditional on Assumption 4.2, which postulates standard terminal cost, terminal set, and local stabilizing gain; these are assumptions used in a standard MPC argument, not conclusions derived from the data fit or from the target result itself. The paper does cite prior and same-author work, notably [2] for the velocity-model structure and [5]/[18] for incremental-input offset-free DPC and robust extensions, but those citations are not load-bearing for the claimed kernelized formulation, the least-squares learning problem, or the vNMPC comparison. A correctness gap appears to exist in the printed proof of Theorem 4.3: the shifted sequence is written in absolute inputs while the velocity model uses increments, and the terminal-invariance step does not explicitly carry the reference subtraction in the error coordinate. However, that is a proof-rigor issue rather than a circularity issue; the result does not reduce by construction to its own inputs. No self-definitional reduction, fitted-parameter-renamed-as-prediction, or author-imported uniqueness argument is present. Thus no significant circularity is found.
Assumptions & free parameters
free parameters (4)
- Kernel width sigma =
sqrt(200)
- Kernel coefficient matrices A_alpha, B_alpha, C_alpha =
Not reported numerically
- Terminal cost P and gain K =
Not reported
- Cost weights Q and R and horizon N =
Q=1000 I3, R=10, N=20
assumptions (5)
- standard math The unknown nonlinear functions in the velocity model lie in the finite-dimensional RKHS spanned by kernel sections at the training points.
- domain assumption Training data are noiseless and the kernel matrices have full rank.
- domain assumption The velocity-form representation in [2, Chapter 4] is exact, with the practical approximation that the convex-hull points can be replaced by current measurements.
- ad hoc to paper Assumption 4.2: a single terminal gain K, matrix P, and set Z_T satisfy the nonlinear LMI and set inclusions.
- domain assumption The disturbance in (1) is piecewise constant and does not enter the stability analysis.
Cite this review
Pith. "Pith review of Kernelized offset-free data-driven predictive control for nonlinear systems." pith.science (2026). https://pith.science/paper/NGMSPNN6
@misc{pith2026241118762,
author = {Pith},
title = {Pith review of: Kernelized offset-free data-driven predictive control for nonlinear systems},
year = {2026},
howpublished = {\url{https://pith.science/paper/NGMSPNN6}},
note = {Machine review of arXiv:2411.18762}
}
read the original abstract
This paper presents a kernelized offset-free data-driven predictive control scheme for nonlinear systems. Traditional model-based and data-driven predictive controllers often struggle with inaccurate predictors or persistent disturbances, especially in the case of nonlinear dynamics, leading to tracking offsets and stability issues. To overcome these limitations, we employ kernel methods to parameterize the nonlinear terms of a velocity model, preserving its structure and efficiently learning unknown parameters through a least squares approach. This results in a offset-free data-driven predictive control scheme formulated as a nonlinear program, but solvable via sequential quadratic programming. We provide a framework for analyzing recursive feasibility and stability of the developed method and we demonstrate its effectiveness through simulations on a nonlinear benchmark example.
Figures
Reference graph
Works this paper leans on
-
[13]
Direct data- driven state-feedback control of general nonlinear system s,
C. V erhoek, P . J. W. Koelewijn, S. Haesaert, and R. T ´ oth, “Direct data- driven state-feedback control of general nonlinear system s,” in 2023 62nd IEEE Conference on Decision and Control , 2023, pp. 3688–3693
work page 2023
-
[14]
C. V erhoek and R. T ´ oth, “Kernel-based multi-step predicto rs for data- driven analysis and control of nonlinear systems through th e velocity form,” arXiv preprint arXiv:2408.00688, 2024
work page Pith review arXiv 2024
-
[18]
Koopman Data-Driven Predictive Control with Robust Stability and R ecursive Feasibility Guarantees,
T. de Jong, V . Breschi, M. Schoukens, and M. Lazar, “Koopman Data-Driven Predictive Control with Robust Stability and R ecursive Feasibility Guarantees,” arXiv, 2405.01292, 2024, to appear, CDC 2024
arXiv 2024
-
[1]
Offset-free mp c explained: novelties, subtleties, and applications,
G. Pannocchia, M. Gabiccini, and A. Artoni, “Offset-free mp c explained: novelties, subtleties, and applications,” IF AC-PapersOnLine, vol. 48, no. 23, pp. 342–351, 2015, 5th IFAC Conference on Nonlinear M odel Predictive Control NMPC 2015
work page 2015
-
[2]
P . S. G. Cisneros, Quasi-linear model predictive control: Stability, modelling and implementation . Technical University Hamburg, 2021, phD thesis, https://doi.org/10.15480/882.3574
-
[3]
Spc: Subspace predic tive control,
W. Favoreel, B. D. Moor, and M. Gevers, “Spc: Subspace predic tive control,” IF AC Proceedings V olumes, vol. 32, no. 2, pp. 4004–4009, 1999, 14th IFAC World Congress 1999, Beijing, China
work page 1999
-
[4]
Data-enabled predi ctive control: In the shallows of the deepc,
J. Coulson, J. Lygeros, and F. D¨ orfler, “Data-enabled predi ctive control: In the shallows of the deepc,” in 2019 18th European Control Conference (ECC). IEEE, 2019, pp. 307–312
work page 2019
-
[5]
Offset–free data–driven predic tive control,
M. Lazar and P . V erheijen, “Offset–free data–driven predic tive control,” in 2022 IEEE 61st Conference on Decision and Control (CDC) . IEEE, 2022, pp. 1099–1104
work page 2022
Show all 19 references
-
[6]
Data-Driven quasi- LPV Model Predictive Control Using Koopman Operator Techni ques,
P . S. Cisneros, A. Datar, P . G¨ ottsch, and H. Werner, “Data-Driven quasi- LPV Model Predictive Control Using Koopman Operator Techni ques,” IF AC-PapersOnLine, vol. 53, no. 2, pp. 6062–6068, 2020, 21st IFAC World Congress
2020
-
[7]
Basis-functions nonlinear data-enabled predic tive control: Consistent and computationally efficient formulations,
M. Lazar, “Basis-functions nonlinear data-enabled predic tive control: Consistent and computationally efficient formulations,” i n 2024 Euro- pean Control Conference (ECC) . IEEE, 2024, pp. 888–893
2024
-
[8]
A generalized r epresenter theorem,
B. Sch¨ olkopf, R. Herbrich, and A. J. Smola, “A generalized r epresenter theorem,” in International conference on computational learning theor y. Springer, 2001, pp. 416–426
2001
-
[9]
Learning controllers from d ata via kernel-based interpolation,
Z. Hu, C. De Persis, and P . Tesi, “Learning controllers from d ata via kernel-based interpolation,” in 2023 62nd IEEE Conference on Decision and Control (CDC) , 2023, pp. 8509–8514
2023
-
[10]
Exploring the links betwee n the fundamental lemma and kernel regression,
O. Molodchyk and T. Faulwasser, “Exploring the links betwee n the fundamental lemma and kernel regression,” IEEE Control Systems Letters, vol. 8, pp. 2045–2050, 2024
2024
-
[11]
K PC: Learning-based model predictive control with determinist ic guarantees,
E. T. Maddalena, P . Scharnhorst, Y . Jiang, and C. N. Jones, “K PC: Learning-based model predictive control with determinist ic guarantees,” in Proceedings of the 3rd Conference on Learning for Dynamics a nd Control, vol. 144. PMLR, 2021, pp. 1015–1026
2021
-
[12]
Robust and kernelized data- enabled predictive control for nonlinear systems,
L. Huang, J. Lygeros, and F. D¨ orfler, “Robust and kernelized data- enabled predictive control for nonlinear systems,” IEEE Transactions on Control Systems Technology , vol. 32, no. 2, pp. 611–624, 2024
2024
-
[15]
Positive definite kernels: past, present a nd future,
G. E. Fasshauer, “Positive definite kernels: past, present a nd future,” Dolomites Research Notes on Approximation , vol. 4, pp. 21–63, 2011
2011
-
[16]
V . I. Paulsen and M. Raghupathi, An introduction to the theory of reproducing kernel Hilbert spaces . Cambridge university press, 2016, vol. 152
2016
-
[17]
Universal kernels,
C. A. Micchelli, Y . Xu, and H. Zhang, “Universal kernels,” Journal of Machine Learning Research , vol. 7, no. 95, pp. 2651–2667, 2006
2006
-
[19]
Constrained model predictive control: Stability and optimality,
D. Mayne, J. Rawlings, C. Rao, and P . Scokaert, “Constrained model predictive control: Stability and optimality,” Automatica, vol. 36, no. 6, pp. 789–814, 2000
2000
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.