An error analysis for absolute and relative approximation

We consider the vectorial algorithm for finding best polynomial approximationsp P _n to a given functionf C[a, b], with respect to the norm · _s , defined byp – f _s =w ₁ (p – f)+w ₂ (p – f) A bound for the modulus of continuity of the best vectorial approximation operator is given, and using the floating point calculus of J. H. Wilkinson, a bound for the rounding error in the algorithm is derived. For givenf, these estimates provide an indication of the conditioning of the problem, an estimate of the obtainable accuracy, and a practical method for terminating the iteration.This paper was supported in part by the Canadian NCR A-8108, FCAC 74-09 and G.E.T.M.A.Part of this research was done during the first-named author's visit to theB! Chair of Applied Mathematics, University of Athens, Spring term, 1975.     

Error bounds for a near-minimax approximation

It is known that a near minimax polynomial approximation tof C [–1, 1] is provided by a finite carrier projectionM _n fromC[–1, 1] onto the subspace of all polynomials of degree n, such thatM _nf is a weighted least squares approximation off on the set consisting of the extreme points of the Chebyshev polynomialT _{2n + 1}. In this paper, upper bounds for the error f–M _n f are given in terms of divided differences.     

Hankel operators and problems of best approximation of unbounded functions

For each function f, f VMO, there exists a unique function f₀, analytic in the circle and such that f–f₀=_f{gVMO_A}. We define the operator of best approximation (nonlinear) A, Af=f₀, fVMO, In the paper one considers the question of the preservation of a class under the action of the operator i.e. finding the classes X, X VMO, AX X. One investigates the classes X containing unbounded functions. It is proved that if P_X is the space of the symbols of the Hankel operators from a Banach space E of functions into the Hardy space H², then AX X. For E one can take almost any space.Translated from Zapiski Nauchnykh Seminarov Leningradskogo Otdeleniya Matematicheskogo Instituta im. V. A. Steklova AN SSSR, Vol. 141, pp. 5–17, 1985.     

Numerical computation of least constants for the Sobolev inequality

Summary Least constantsc for the well-known Sobolev inequality fcf _{m, G} ,fH ^m (G) are obtained in closed form by a reproducing kernel technique, where the Sobolev spaceH ^m (G) for a domainG in ⁿ is defined as the completion ofC ^m (G) with respect to the Sobolev norm given by , where is the norm ofL ₂ (G) and is the supremum norm onG. Numerical values for the case whereG is the ⁿ are given.     

High levelm-harmonic cardinal B-splines

We generalize the notion ofm-harmonic cardinal B-spline defined in [Rabut, [6c]] to obtain B-splines on an infinite regular grid, which are halfway between elementary B-splines and the Lagrangean cardinal spline function. We give the main properties of these functions: Fourier transform, decay when x , integration,P _k -reproduction (fork<-2m–1) of the associated B-spline approximation, etc. We show that, in some sense, high levelm-harmonic B-splines may be considered as a finer regular approximation of the Dirac distribution than the elementarym-harmonic B-splines are.     

Isomorphismen von rechtseiträumen

Spaces called rectangular spaces were introduced in [5] as incidence spaces (P,G) whose set of linesG is equipped with an equivalence relation and whose set of point pairs P² is equipped with a congruence relation , such that a number of compatibility conditions are satisfied. In this paper we consider isomorphisms, automorphisms, and motions on the rectangular spaces treated in [5]. By an isomorphism of two rectangular spaces (P,G, , ) and (P,G, , ) we mean a bijection of the point setP onto P which maps parallel lines onto parallel lines and congruent points onto congruent points. In the following, we consider only rectangular spaces of characteristic 2 or of dimension two. According to [5] these spaces can be embedded into euclidean spaces. In case (P,G, , ) is a finite dimensional rectangular space, then every congruence preserving bijection ofP onto P is in fact an isomorphism from (P,G, , ) onto (P,G, , ) (see (2.4)). We then concern ourselves with the extension of isomorphisms. Our most important result is the theorem which states that any isomorphism of two rectangular spaces can be uniquely extended to an isomorphism of the associated euclidean spaces (see (3.2)). As a consequence the automorphisms of a rectangular space (P,G, , ) are precisely the restrictions (onP) of the automorphisms of the associated euclidean space which fixP as a whole (see (3.3)). Finally we consider the motions of a rectangular space (P,G, , ). By a motion of(P. G,, ) we mean a bijection ofP which maps lines onto lines, preserves parallelism and satisfies the condition((x), (y)) (x,y) for allx, y P. We show that every motion of a rectangular space can be extended to a motion of the associated euclidean space (see (4.2)). Thus the motions of a rectangular space (P,G, , ) are seen to be the restrictions of the motions of the associated euclidean space which mapP into itself (see (4.3)). This yields an explicit representation of the motions of any rectangular plane (see (4.4)).

---

---

Herrn Professor Burau zum 85. Geburtstag gewidmet     

Computing the Hilbert transform of a Jacobi weight function

We discuss the evaluation of the Hilbert transformf _–1 ¹ (t-)^–1 w(, )(t)dt,–1<<1, of the Jacobi weight functionw(, )(t)=(1–t))(1+t) by analytic and numerical means and also comment on the recursive computation of the quantitiesf _–1 ¹ )(t–)^–1 _n (t;w ^{(, ))} w ^{(, )}(t)dt,n=0, 1, 2, ..., where _n (·;w ^{(, )}) is the Jacobi polynomial of degreen.The work of the first author was supported in part by the National Science Foundation under grant DCR-8320561. The work of the second author was supported by the National Science Foundation under grant DMS-8419086.     

Theorem on a convergence condition in the spaces

For a given -function (u), a condition on a -function (u) is found such that it is necessary and sufficient for the following to hold: if f_n(x) f(x) and f _n (x)M (n=1, 2, ...) where M>0 is an absolute constant, then f _n (x)–f(x)0(n). An analogous condition for convergence in Orlicz spaces is obtained as a corollary.Translated from Matematicheskie Zametki, Vol. 21, No. 5, pp. 615–626, May, 1977.The author thanks V. A. Skvortsov for his constant attention and guidance on this paper.     

The Spectral Theory of Amenable Actions and Invariants of Discrete Groups

Let G denote a semisimple group, a discrete subgroup, B=G/P the Poisson boundary. Regarding invariants of discrete subgroups we prove, in particular, the following:(1) For any -quasi-invariant measure on B, and any probablity measure on , the norm of the operator () on L ²(B,) is equal to (), where is the unitary representation in L ²(X,), and is the regular representation of .(2) In particular this estimate holds when is Lebesgue measure on B, a Patterson–Sullivan measure, or a -stationary measure, and implies explicit lower bounds for the displacement and Margulis number of (w.r.t. a finite generating set), the dimension of the conformal density, the -entropy of the measure, and Lyapunov exponents of .(3) In particular, when G=PSL₂() and is free, the new lower bound of the displacement is somewhat smaller than the Culler–Shalen bound (which requires an additional assumption) and is greater than the standard ball-packing bound.We also prove that ()=_G() for any amenable action of G and L ¹(G), and conversely, give a spectral criterion for amenability of an action of G under certain natural dynamical conditions. In addition, we establish a uniform lower bound for the -entropy of any measure quasi-invariant under the action of a group with property T, and use this fact to construct an interesting class of actions of such groups, related to 'virtual' maximal parabolic subgroups. Most of the results hold in fact in greater generality, and apply for instance when G is any semi-simple algebraic group, or when is any word-hyperbolic group, acting on their Poisson boundary, for example.     

Superstable implicit Runge-Kutta methods for second order initial value problems

Using the well known properties of thes-stage implicit Runge-Kutta methods for first order differential equations, single step methods of arbitrary order can be obtained for the direct integration of the general second order initial value problemsy=f(x, y, y),y(x _o)=y _o,y(x _o)=y _o. These methods when applied to the test equationy+2y+ ² y=0, ,0, +>0, are superstable with the exception of a finite number of isolated values ofh. These methods can be successfully used for solving singular perturbation problems for which f/y and/or f/y are negative and large. Numerical results demonstrate the efficiency of these methods.     

A property of a system of functions close to exponential functions

We consider the system {f _n=xⁿ[l+_n]} in the interval [a,b] (0 a n > 0 and _n(x) such as the condition, we obtain a bound for the coefficients of the polynomial P(x)=#x2211;c_n f _n(x) in terms of P(x)L_p[a,b]. It is found that this bound is not valid without this condition (assuming the other conditions to remain the same).Translated from Matematicheskie Zametki, Vol. 12, No. 1, pp. 29–36, July, 1972.     

Fehlerabschätzungen für Gauß-Quadraturformeln

Summary We consider Gauss quadrature formulaeQ _n ,n, approximating the integral ,w an even weight function. Let be analytic inK _r :={z:|z|<r},r>1, and . The error functionalR _n :=I-Q _n is continuous with respect to |·|_r and the relation , q_2k (x):=x ^2k holds.In this paper estimates for R _n are given. To this end we first derive two new representations of R _n which are essential for our further investigations. The R _n =r ² R _n (), with (x):=1/(r ²-x ²), is estimated in various ways by using the best uniform approximation of in P_2n-1, and also the expansion of with respect to Chebyshe polynomials of the first and second kind. Forw(x)=(1-x ²), =±1/2, R _n is calculated. The asymptotic behaviour, forr1+, of R _n and of the derived error bounds is also discussed. Finally, we compare different error bounds and give numerical examples.

     

Hilbertraum-Approximation und Jacobi-Matrizen

Zusammenfassung SeiH ein Hilbertraum mit Norm und Skalarprodukt (,). Seien _n ,n=0, 1, 2, ... undf Elemente H, s, i natürliche Zahlen und >0 derart dass gilt: 1) _n =1,2)( _f ), _k )=0 für j–k >s, 3) ( _f , _n ) fürj k, 4) (f, _f =0 fürj>i. Seienf _n undf ^* die Projektionen vonf auf die durch ₀,..., _n bezw. ₀,₁,₂,... aufgespannten abgeschlossenen Unterräume. Das Hauptresultat der Arbeit besagt dass für hinreichend kleines KonstantenC>0 und 0<q<1 existieren mit: f _n -f ^*Cq ⁿ . Dieses Resultat steht in enger Beziehung zu gewissen unendlichen MatrizenA=(a _jk ), die charakterisiert sind durch: (*) es existiertm>0 so dassa _jk =0 für |j-k|>m. Das Hauptresultat wird auf unendliche lineare GleichungssystemeAf=0,Af=g angewandt, woA eine Matrix mit der Eigenschaft (*) ist, deren Diagonalen gewissen Wachstumsbedingungen genügen.

---

LetH be a Hilbert space with norm and scalar product (,). Let _n ,n=0, 1, 2, ... andf be elements H, i, s integers and >0 such that: 1) _n =1,2)( _f ), _k )=0 for j–k >s, 3) ( _f , _n ) forj k, 4) (f, _f =0 forj>i. Letf _n andf ^* be the projections off onto the closed subspaces spanned by ₀,..., _n and ₀,₁,₂ .... respectively. The main result says that for sufficiently small there are constantsC>0 and 0<q<1 with f _n -f ^*Cq ⁿ . This result is closely related to certain infinite matricesA=(a _jk ) with the property: (*) there exists anm>0 such thata _jk =0 for |j-k|>m. The result is applied to infinite systems of linear equationsAf=0 andAf=g, whereA is a matrix with property (*), whose diagonals satisfy certain growth conditions.

---

---

Dem Andenken von Professor Eduard Stiefel gewidmet     

On discrete and continuous norms in Paley-Wiener spaces and consequences for exponential frames

It is well known that for certain sequences {t_n}_n the usual L^p norm ·_p in the Paley-Wiener space PW ^p is equivalent to the discrete norm f_p,{tn}:=( _n=– |f(t_n)|^p)^1/p for 1 p = < and f_,{tn}:=sup _n|f(t_n| for p=). We estimate f_p from above by Cf_p, _n and give an explicit value for C depending only on p, , and characteristic parameters of the sequence {t_n}_n. This includes an explicit lower frame bound in a famous theorem of Duffin and Schaeffer.     

Existenzaussagen für gewisse Zweipunkt-Randwertaufgaben

Zusammenfassung Für Randwertaufgaben der Form–u–l ₀ ...u–l ₀ ...u=f(x, u) mitl ₀R,lR,f definiert und stetig auf {a<-x<-b, |u|<} wird eine Existenzaussage gewonnen, fallsf inu linear durch die aufeinanderfolgenden Eigenwerte der zugehörigen linearen Aufgabe beschränkt ist. Zum Beweis betrachtet man die äquivalente Hammersteinsche Integralgleichung mit nichtsymmetrischem Kern. Mit Hilfe des Schauderschen Fixpunktsatzes erhält man für diese Integralgleichung Existenzaussagen, welche Ergebnisse von Dolph verallgemeinern.

---

Summary This note contains an existence theorem for a two-point boundary value problem of the form–u–l ₀ u–l ₀ u=f(x, u) wherel ₀R,l ₀R,f defined and continuous on {a<-xb, |u|<} iff is linear bounded inu by the successive eigenvalues of the corresponding linear problem. To prove this result we consider the equivalent Hammerstein integral equation with non-symmetric kernel. Schauders fixpoint principle supplies existence theorems for integral equations of this type which generalize results of Dolph in some sense.

     

Elementarym-harmonic cardinal B-splines

We generalize the notion of B-spline to the thin plate splines and to otherd-dimensional polyharmonic splines as defined in [Duchon, [3]]; for regular nets, we give the main properties of these B-splines: Fourier transform, decay when x , stability, integration property, links between B-splines of different orders or of different dimensions and in particular link with the polynomial B-splines, approximation using B-splines... We show that, in some sense, B-splines may be considered as a regularized form of the Dirac distribution.     

L p -estimates for the solutions of second order elliptic equations

Summary We investigate the homogeneous Dirichlet problem in H^2,p for a second order elliptic partial differential equation in nondivergence form Lu=f in the case in which the leading coefficients of L belong to H^1,n(), Rⁿ. We prove that if p belongs to a suitable neighbourhood of 2, then the above problem, has a unique solution u satisfying D²u_p Cf_p; furthermore, if f H^k,p, k=1,2, ..., and the coefficients of L satisfy some natural conditions, then the solution satisfies .Lavoro eseguito nell'ambito del gruppi 40% e 60% del M.P.I.     

Best approximation and best unilateral approximation of the kernel of a biharmonic equation and optimal renewal of the values of operators

For the classB _p , 0 < 1, 1p , of 2-periodic functions of the form f(t)=u(,t), whereu (,t) is a biharmonic function in the unit disk, we obtain the exact values of the best approximation and best unilateral approximation of the kernel K(t) of the convolution f= K *g, gl, with respect to the metric of L₁. We also consider the problem of renewal of the values of the convolution operator by using the information about the values of the boundary functions.Translated from Ukrainskii Matematicheskii Zhurnal, Vol.47, No. 11, pp. 1549–1557, November, 1995.     

On quadratic spline interpolation

Using the quadratic spline interpolates(x) fitting the data (x _i,y _i), 0in and satisfying the end conditions_o=y_o, we give formulae approximatingy andy at selected knots by orders up toO(h ⁴).     

L 2-approximation by the translates of a function and related attenuation factors

Summary Let be thek-dimensional subspace spanned by the translates (·–2j/k),j=0, 1, ...,k–1, of a continuous, piecewise smooth, complexvalued, 2-periodic function . For a given functionfL ²(–, ), its least squares approximantS _kf from can be expressed in terms of an orthonormal basis. Iff is continuous,S _kf can be computed via its discrete analogue by fast Fourier transform. The discrete least squares approximant is used to approximate Fourier coefficients, and this complements the works of Gautschi on attenuation factors. Examples of include the space of trigonometric polynomials where is the de la Valleé Poussin kernel, algebraic polynomial splines where is the periodic B-spline, and trigonometric polynomial splines where is the trigonometric B-spline.