Picard Iteration for Parameter Estimation in Nonlinear Ordinary Differential Equations
Abstract
We consider the problem of using experimental time-series data for parameter estimation in nonlinear ordinary differential equations, focusing on the case where the data is noisy, sparse, irregularly sampled, includes multiple experiments, and does not directly measure the system state or its time-derivative. To account for such low-quality data, we propose a new framework for gradient-based parameter estimation which uses the Picard operator to reformulate the problem as constrained optimization with infinite-dimensional variables and constraints. We then use the contractive properties of the Picard operator to propose a class of gradient-contractive algorithms and provide conditions under which such algorithms are guaranteed to converge to a local optima. The algorithms are then tested on a battery of models and variety of datasets in order to demonstrate robustness and improvement over alternative approaches.
Cite
@article{arxiv.2412.20216,
title = {Picard Iteration for Parameter Estimation in Nonlinear Ordinary Differential Equations},
author = {Aleksandr Talitckii and Matthew M. Peet},
journal= {arXiv preprint arXiv:2412.20216},
year = {2025}
}