• 検索結果がありません。

ETNA Kent State University http://etna.math.kent.edu

N/A
N/A
Protected

Academic year: 2022

シェア "ETNA Kent State University http://etna.math.kent.edu"

Copied!
1
0
0

読み込み中.... (全文を見る)

全文

(1)

Electronic Transactions on Numerical Analysis.

Volume 38, pp. 146-167, 2011.

Copyright2011, Kent State University.

ISSN 1068-9613.

ETNA

Kent State University http://etna.math.kent.edu

ROBUST RATIONAL INTERPOLATION AND LEAST-SQUARES

PEDRO GONNETy, RICARDO PACH ´ONy,ANDLLOYD N. TREFETHENy

Abstract. An efficient and robust algorithm and a Matlab coderatdiskare presented for rational interpolation or linearized least-squares approximation of a function based on its values at points equally spaced on a circle. The use of the singular value decomposition enables the detection and elimination of spurious poles or Froissart doublets that commonly complicate such fits without contributing to the quality of the approximation. As an application, the algorithm leads to a method for the stable computation of certain radial basis function interpolants in the difficult case of smoothness parameter"close to zero.

Key words. Rational interpolation, spurious poles, Froissart doublets, Pad´e approximation, radial basis func- tions, ratdisk, singular value decomposition

AMS subject classifications. 41A20, 41A21, 65D05

Received February 10, 2011. Accepted for publication February 28, 2011. Published online May 18, 2011.

Recommended by L. Reichel. P. G. was supported by Swiss National Science Foundation Individual Support Fel- lowships Nr. PBEZP2-127959 and Nr. PA00P2-134146.

yOxford University Mathematical Institute, 25-29 St Giles, Oxford OX1 3LB, UK (Pedro.Gonnet,[email protected]).

146

参照

関連したドキュメント

Examples are presented for: general dense ma- trices, upper triangular matrices, higher order generator semiseparable matrices, quasiseparable matrices, Givens- vector

In this paper we demonstrate, for the first time, that deflation preconditioning can be applied in communication-avoiding formulations of Lanczos-based Krylov methods such as

SOLVING REGULARIZED LINEAR LEAST-SQUARES PROBLEMS BY THE ALTERNATING DIRECTION METHOD WITH APPLICATIONS TO IMAGE..

It is based on a Petrov-Galerkin process and multistep schemes and consists of building, throughout the iterations, an approximation subspace using the previous computations, where

We then introduced the convection-diffusion control problem and illustrated that, with a suitable stabilization technique (the Local Projection Stabilization), the same saddle

We have seen that under rather natural source condi- tions error estimates in Bregman distances can be extended from the well-known quadratic fitting (Gaussian noise) case to

Mainly, by using the extrapolation method, families of estimates can be derived which are valid for any nonsingular matrix and thus can be used for nonsymmetric problems. In

Using Fourier expansions, the solutions of the resulting systems of ordinary differential equations for the Fourier amplitudes can be written, after truncation, in form of