Article,

Vector versions of Prony's algorithm and vector-valued rational approximations

.
Advances in Computational Mathematics, 46 (2): 30 (Mar 19, 2020)
DOI: 10.1007/s10444-020-09751-9

Abstract

Given the scalar sequence \$\backslash\f\_\m\\backslash\^\\backslashinfty \\_\m=0\\$ that satisfiesfm=∑i=1kai$\zeta$im,m=0,1,łdots,\$\$ f\_\m\ = \backslashsum\backslashlimits\_ı=1\^\k\ \a\_ı\\\\backslashzeta\\_ı\^\m\,\backslashquad m=0,1,\backslashldots, \$\$ where \$a\_ı\, \backslashzeta \_ı\\backslashin \backslashmathbb \C\\$ and $\zeta$i are distinct, the algorithm of Prony concerns the determination of the ai and the $\zeta$i from a finite number of the fm. This algorithm is also related to Padé approximants from the infinite power series \$\\backslashsum \^\\backslashinfty \\_\j=0\f\_\j\z^\j\\$. In this work, we discuss ways of extending Prony's algorithm to sequences of vectors \$\\backslash\\backslashboldsymbol \f\\_\m\\backslash\\^\\backslashinfty \\_\m=0\\$ in \$\backslashmathbb \C\^\N\\$ that satisfyfm=∑i=1kai$\zeta$im,m=0,1,łdots,\$\$ \backslashboldsymbol\f\\_\m\ = \backslashsum\backslashlimits\_ı=1\^\k\ \backslashboldsymbol\a\\_ı\ \\backslashzeta\\_ı\^\m\, \backslashquad m=0,1,\backslashldots, \$\$ where \$\backslashboldsymbol \a\\_ı\\backslashin \backslashmathbb \C\^\N\\$ and \$\backslashzeta \_ı\\backslashin \backslashmathbb \C\\$. Two distinct problems arise depending on whether the vectors ai are linearly independent or not. We consider different approaches that enable us to determine the ai and $\zeta$i for these two problems, and develop suitable methods. We concentrate especially on extensions that take into account the possibility of the components of the ai being coupled. One of the applications we consider concerns the case in whichfm=∑i=1rai$\zeta$im,m=0,1,łdots,rlarge,\$\$ \backslashboldsymbol\f\\_\m\ = \backslashsum\backslashlimits\_ı=1\^\r\ \backslashboldsymbol\a\\_ı\ \\backslashzeta\\_ı\^\m\, \backslashquad m=0,1,\backslashldots,\backslashquad r \backslashtext\ large\, \$\$ and we would like to approximate/determine of a number of the pairs ($\zeta$i, ai) for which |$\zeta$i| are largest. We present the related theory and provide numerical examples that confirm this theory. This application can be extended to the more general case in whichfm=∑i=1rpi(m)$\zeta$im,m=0,1,łdots,\$\$ \backslashboldsymbol\f\\_\m\ = \backslashsum\backslashlimits\_ı=1\^\r\ \backslashboldsymbol\p\\_ı\ (m)\\backslashzeta\\_ı\^\m\, \backslashquad m=0,1,\backslashldots, \$\$ where \$\backslashboldsymbol \p\\_ı\(m)\backslashin \backslashmathbb \C\^\N\\$ are some (vector-valued) polynomials in m, and \$\backslashzeta \_ı\\backslashin \backslashmathbb \C\\$ are distinct. Finally, the methods suggested here can be extended to vector sequences in infinite dimensional spaces in a straightforward manner.

Tags

Users

  • @dblp
  • @gdmcbain

Comments and Reviews