Stable polefinding and rational least-squares fitting via eigenvalues
A common way of finding the poles of a meromorphic function f in a domain, where an explicit expression of f is unknown but f can be evaluated at any given z, is to interpolate f by a rational function pq such that r(γi)=f(γi) at prescribed sample points {γi}Li=1 , and then find the roots of...
| Main Authors: | , |
|---|---|
| Format: | Journal article |
| Language: | English |
| Published: |
Springer
2018
|
| _version_ | 1826257378390048768 |
|---|---|
| author | Ito, S Nakatsukasa, Y |
| author_facet | Ito, S Nakatsukasa, Y |
| author_sort | Ito, S |
| collection | OXFORD |
| description | A common way of finding the poles of a meromorphic function f in a domain, where an explicit expression of f is unknown but f can be evaluated at any given z, is to interpolate f by a rational function pq such that r(γi)=f(γi) at prescribed sample points {γi}Li=1 , and then find the roots of q. This is a two-step process and the type of the rational interpolant needs to be specified by the user. Many other algorithms for polefinding and rational interpolation (or least-squares fitting) have been proposed, but their numerical stability has remained largely unexplored. In this work we describe an algorithm with the following three features: (1) it automatically finds an appropriate type for a rational approximant, thereby allowing the user to input just the function f, (2) it finds the poles via a generalized eigenvalue problem of matrices constructed directly from the sampled values f(γi) in a one-step fashion, and (3) it computes rational approximants p^,q^ in a numerically stable manner, in that (p^+Δp)/(q^+Δq)=f with small Δp,Δq at the sample points, making it the first rational interpolation (or approximation) algorithm with guaranteed numerical stability. Our algorithm executes an implicit change of polynomial basis by the QR factorization, and allows for oversampling combined with least-squares fitting. Through experiments we illustrate the resulting accuracy and stability, which can significantly outperform existing algorithms. |
| first_indexed | 2024-03-06T18:17:15Z |
| format | Journal article |
| id | oxford-uuid:0508eab6-abc2-45f1-ac55-998a51da82fb |
| institution | University of Oxford |
| language | English |
| last_indexed | 2024-03-06T18:17:15Z |
| publishDate | 2018 |
| publisher | Springer |
| record_format | dspace |
| spelling | oxford-uuid:0508eab6-abc2-45f1-ac55-998a51da82fb2022-03-26T08:55:00ZStable polefinding and rational least-squares fitting via eigenvaluesJournal articlehttp://purl.org/coar/resource_type/c_dcae04bcuuid:0508eab6-abc2-45f1-ac55-998a51da82fbEnglishSymplectic Elements at OxfordSpringer2018Ito, SNakatsukasa, YA common way of finding the poles of a meromorphic function f in a domain, where an explicit expression of f is unknown but f can be evaluated at any given z, is to interpolate f by a rational function pq such that r(γi)=f(γi) at prescribed sample points {γi}Li=1 , and then find the roots of q. This is a two-step process and the type of the rational interpolant needs to be specified by the user. Many other algorithms for polefinding and rational interpolation (or least-squares fitting) have been proposed, but their numerical stability has remained largely unexplored. In this work we describe an algorithm with the following three features: (1) it automatically finds an appropriate type for a rational approximant, thereby allowing the user to input just the function f, (2) it finds the poles via a generalized eigenvalue problem of matrices constructed directly from the sampled values f(γi) in a one-step fashion, and (3) it computes rational approximants p^,q^ in a numerically stable manner, in that (p^+Δp)/(q^+Δq)=f with small Δp,Δq at the sample points, making it the first rational interpolation (or approximation) algorithm with guaranteed numerical stability. Our algorithm executes an implicit change of polynomial basis by the QR factorization, and allows for oversampling combined with least-squares fitting. Through experiments we illustrate the resulting accuracy and stability, which can significantly outperform existing algorithms. |
| spellingShingle | Ito, S Nakatsukasa, Y Stable polefinding and rational least-squares fitting via eigenvalues |
| title | Stable polefinding and rational least-squares fitting via eigenvalues |
| title_full | Stable polefinding and rational least-squares fitting via eigenvalues |
| title_fullStr | Stable polefinding and rational least-squares fitting via eigenvalues |
| title_full_unstemmed | Stable polefinding and rational least-squares fitting via eigenvalues |
| title_short | Stable polefinding and rational least-squares fitting via eigenvalues |
| title_sort | stable polefinding and rational least squares fitting via eigenvalues |
| work_keys_str_mv | AT itos stablepolefindingandrationalleastsquaresfittingviaeigenvalues AT nakatsukasay stablepolefindingandrationalleastsquaresfittingviaeigenvalues |