8
$\begingroup$

Chebyshev theory provides a very effective method for approximating continuous real valued functions on the unit interval. Is there something similar for continuous real valued functions on the closed unit disc? I can only find papers which discuss analytic complex valued functions, which is not what I need.

UPDATE: Some basic features of Chebyshev theory are as follows. Suppose we have a continuous function $f$ on $[-1,1]$ (which may be expensive to compute). Let $P_nf$ be the unique polynomial of degree $n$ that agrees with $f$ at the Chebyshev points $\cos(k\pi/n)$ for $0\leq k\leq n$. There is an efficient procedure for computing $P_nf$, which is tolerant of rounding errors. It does not involve any integration: we just need to evaluate $f$ at $n+1$ points. There are various results about the convergence of $P_nf$ to $f$, assuming various hypotheses about the smoothness or real-analyticity of $f$. In typical examples, convergence is rapid. There is an extensive theory surrounding all this, in which the Chebyshev polynomials $T_n(x)=\cos(n\,\cos^{-1}(x))$ play a central role.

It is known that everything works much worse if we use equally spaced points instead of Chebyshev points. This is closely related to the Bernstein approximation system mentioned by Liviu Nicolaescu. That has some attractive theoretical properties, but usually converges quite slowly. Similarly, even if $f$ is analytic, the polynomials $P_nf$ typically converge to $f$ more rapidly than the Taylor approximations.

$\endgroup$
  • 1
    $\begingroup$ The Stone-Weierstrass theorem says that if $A$ is a subalgebra of the algebra $C(X)$ of real continuous functions on a compact metric space $X$, and $A$ separates points, that is for $x \neq y$ there are $f,g \in A$ with $f(x) \neq g(y)$, then the closure of $A$ in the uniform norm is the whole of $C(X)$. So for instance the set of polynomials on the disk is such. What exactly are you looking for? I must mention that I'm in pure math, so maybe I'm missing the point here. $\endgroup$ – Frol Zapolsky Jan 11 '15 at 12:14
  • 1
    $\begingroup$ @FrolZapolsky Stone-Weierstrass is next to useless on its own if one wants concrete approximants by polynomials or estimates of rate of convergence in terms of moduli of continuity/smoothness, etc. Note that the OP seems to be asking for effective approximation schemes, not nonconstructive density results $\endgroup$ – Yemon Choi Jan 11 '15 at 14:01
  • $\begingroup$ Neil: perhaps you should remind people what Chebyshev's theory does for the interval - I have a vague recollection that it deals with best interpolants/approximants of a given degree, but perhaps I am misremembering? $\endgroup$ – Yemon Choi Jan 11 '15 at 14:04
  • 1
    $\begingroup$ Maybe this paper hans.munthe-kaas.no/Chebyshev/Welcome_files/munthekaas09mcp.pdf and the references therein might have the answer you are looking for. Apparently some of the one-dimensional Chebyshev interpolation extends to interpolations on fundamental domains of affine Weyl groups. the above paper explains how. $\endgroup$ – Liviu Nicolaescu Jan 11 '15 at 20:06
  • 1
    $\begingroup$ @LiviuNicolaescu: Thanks. By following the references in the paper you mentioned, I found the book "Orthogonal polynomials of several variables" by Dunkl and Xu. It looks like I could extract an answer from there with some labour. I am still hoping that someone will point me to a more convenient reference. $\endgroup$ – Neil Strickland Jan 11 '15 at 21:03
1
$\begingroup$

If the main interest lies in numerical methods designed to compute specifically with functions defined on a disk, the following recent manuscript may be useful:

https://arxiv.org/pdf/1604.03061.pdf

A few basic (theoretical) differences between Chebyshev polynomial approximations on an interval and on a disk are as follows:

the main feature of Chebyshev approximation by polynomials on an interval is the fact that the family of approximants, at each degree $n$, is spanned by a Chebyshev system, that is a set of continuous functions such that every nonzero linear combination has at most $n$ zeros on the interval. This implies existence and uniqueness of a best approximant $p_{n}^{*}$, and also a characterization of $p_{n}^{*}$ in terms of the error function which alternates in sign at $n+2$ extremal points of the interval. A study of general Chebyshev systems can be found in [2].

When approximating real functions on the unit disk by real polynomials in $x$ and $y$, the property of being a Chebyshev system is obviously lost, and thus also the ``Chebyshev theory''. Thus, different arguments are used:

  • Existence of a best approximant follows from the continuity of the sup norm with respect to the coefficients of the polynomials.
  • In general uniqueness of a best approximant is false. For instance, it is known that a best approximant to the monomial $x^{n}y^{m}$ in the disk by polynomials of total degree less than $n+m$ is unique if and only if $n=0$, $m=0$, or $n=m=1$.
  • Characterization of a best approximant can be given. The notion of extremal signatures is sometimes used, see [4]. Such characterizations can be seen as extensions of the alternation property that holds for the Chebyshev systems. They essentially follow from the Hahn-Banach separation theorem. For instance, the classical Kolmogorov criterion ([4] p.65 or [1] p.10) gives the following characterization of a best uniform approximant from a convex set $M$ in $C(X)$ where $X$ is a compact set:

The function $u_0\in M$ is a best uniform approximant to $f\in C(X)$ from $M$ iff $$\forall u\in M,\quad\inf_{x\in P[\epsilon_0]}\text{Re}(\overline{\epsilon_0(x)}(u(x)-u_0(x)))\leq 0,$$ where $\epsilon_0=f-u_0$ and $P[\epsilon_0]=\{x\in X,~|\epsilon_0(x)|=\|\epsilon_0\|\}$. In the particular case when $M$ is an $n$-dimensional linear subspace of real-valued (resp. complex-valued) functions of $C(X)$, a characterization is as follows:

There are $m$ points $x_i\in P[\epsilon_0]$ and $m$ positive numbers $\theta_i$, $m\leq n+1$ (resp. $m\leq 2n+1$) such that $$\forall u\in M,\quad\sum_{i=1}^m\theta_i\overline{\epsilon_0(x_i)}u(x_i)=0.$$ From a numerical point of view, practical implementations of Chebyshev approximation on a - square - are discussed in [5], considering approximation of a continuous function by expansions in terms of products $T_{n}(x)T_{m}(y)$, $n,m\geq 0$, where $T_{n}$ denotes the univariate Chebyshev polynomials.

Finally, analogs of Chebyshev polynomials in the disk have also been studied in [3].

References

[1] D. Braess, Nonlinear approximation theory. Springer Series in Computational Mathematics, 7. Springer-Verlag, Berlin, 1986.

[2] S. Karlin, W.J. Studden, Chebyshev systems: with applications in analysis and statistics, Wiley, New York, 1966.

[3] I. Moale; F. Peherstorfer, Explicit min-max polynomials on the disc. J. Approx. Theory 163 (2011), no. 6, 707-723.

[4] T. Rivlin, Chebyshev polynomials, From approximation theory to algebra and number theory, Pure and Applied Mathematics, Wiley, Inc., New York, 1990.

[5] A. Sommariva, M. Vianello, R. Zanovello, Adaptive bivariate Chebyshev approximation. Numer. Algorithms 38 (2005), no. 1-3, 79-94.

$\endgroup$
3
$\begingroup$

The uniform limit of a sequence of complex polynomials is a holomorphic function. To get beyond holomorphic functions you need to work with polynomials depending on both $z$ and $\bar{z}$. Check this older MO question which deals with uniform approximations on the square via multivariate Bernstein polynomials.

The disk case can be reduced to the square case as follows. Given a continuous function $f(r,\theta)$ on the disk $r\leq 1$ we can extend it to a continuous function $F(r,\theta)$ in the plane by setting

$$ F(r,\theta):=\begin{cases} f(r, \theta), & r\leq 1,\\ f(1,\theta), & r\geq 1. \end{cases} $$

Next approximate $F$ on $[-1,1]\times [-1,1]$ using the bivariate Bernstein polynomials as detailed in the linked MO post.

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.