Parameter estimation of ordinary differential equations

Zhengfeng Li, Michael R. Osborne*, Tania Prvan

*Corresponding author for this work

    Research output: Contribution to journalArticlepeer-review

    64 Citations (Scopus)


    This paper addresses the development of a new algorithm for parameter estimation of ordinary differential equations. Here, we show that (1) the simultaneous approach combined with orthogonal cyclic reduction can be used to reduce the estimation problem to an optimization problem subject to a fixed number of equality constraints without the need for structural information to devise a stable embedding in the case of non-trivial dichotomy and (2) the Newton approximation of the Hessian information of the Lagrangian function of the estimation problem should be used in cases where hypothesized models are incorrect or only a limited amount of sample data is available. A new algorithm is proposed which includes the use of the sequential quadratic programming (SQP) Gauss-Newton approximation but also encompasses the SQP Newton approximation along with tests of when to use this approximation. This composite approach relaxes the restrictions on the SQP Gauss-Newton approximation that the hypothesized model should be correct and the sample data set large enough. This new algorithm has been tested on two standard problems.

    Original languageEnglish
    Pages (from-to)264-285
    Number of pages22
    JournalIMA Journal of Numerical Analysis
    Issue number2
    Publication statusPublished - Apr 2005


    Dive into the research topics of 'Parameter estimation of ordinary differential equations'. Together they form a unique fingerprint.

    Cite this