Publication Cover
Applicable Analysis
An International Journal
Volume 88, 2009 - Issue 1
283
Views
52
CrossRef citations to date
0
Altmetric
Original Articles

Parameter estimation for signals described by differential equations

Pages 29-52 | Received 01 Feb 2008, Accepted 13 Oct 2008, Published online: 21 May 2009
 

Abstract

An estimation method is presented for signals described by linear differential equations whose coefficients are functions of the unknown parameters. Various types of identifiability are defined according to these functions. In the linear case, an estimation algorithm is derived directly from the identifiability conditions, by use of elementary rules of operational calculus. The main steps of the algorithm are first presented through an introductory example of a damped sinusoid estimation. Unlike the modified Prony's method (used as reference), the presented one (1) provides explicit closed-form expressions and (2) remains valid for linear differential equations with non-constant coefficients. These expressions depend only on iterated integrals of the observed signals. Some basic features of the method are analysed and, in particular, we show that a least squares interpretation can be attached to it. Application to the estimation of a chirp signal is also presented with numerical simulations.

AMS Subject Classifications::

Notes

Notes

1. See also Citation8 for an algebraic frequency estimation.

2. The invertibility of the system in the operational domain is not necessary. For example, although x and are proportional, their respective images in the time domain, x and , are not.

3. Note that a fully discrete-time counterpart of the presented method can be found in Citation15–17.

4. In all the simulations, the integrals are numerically computed using the ‘work horse’ trapezoidal method.

5. The Matlab code used for simulating this method was found at http://www.statsci.org

6. In general either d of D is zero.

7. Compare with (Equation24) which describes the true parameter vector Θ. Indeed, since the term Q′ in the second member of (Equation24) depends on Θ, that equation is not computable.

Reprints and Corporate Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

To request a reprint or corporate permissions for this article, please click on the relevant link below:

Academic Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

Obtain permissions instantly via Rightslink by clicking on the button below:

If you are unable to obtain permissions via Rightslink, please complete and submit this Permissions form. For more information, please visit our Permissions help page.