Boundary effects on convergence rates for Tikhonov regularization

We consider the Tikhonov regularizer f_λ of a smooth function f ε H^2m[0, 1], defined as the solution (see [1]) to We prove that if f^(j)(0) = f^(j)(1) = 0, J = m, …, k < 2m − 1, then ¦f − f_λ¦_j² Rλ^{(2k − 2j + 3)/2m}, J = 0, …, m. A detailed analysis is given of the effect of the boundary on convergence rates.     

The approximate identity kernels of product type for the Walsh system

At present there are only a few approximate identity kernels for the Walsh system, for example, the p^N-truncated Dirichlet kernel D_{p^N − 1}(t) = ∑_{j = 0}^{pN − 1} w_j(t) [6]; the Abel-Poisson kernel λ_γ(t) = ∑_{k = 0}^∞ γ^kw_k(t) [3], and so on. In [6], Zheng has introduced a new kind of approximate identity kernels for the Walsh system—the kernels of product type. In the present paper we discuss the approximation properties of such product type kernels. Estimates of their moments as well as a direct approximation theorem are obtained. Then, to establish an inverse approximation theorem, we need the p-adic derivative of product type kernels and we estimate this derivative in L¹-norm.     

A Remez-Type Inequality for Non-dense Müntz Spaces with Explicit Bound

LetΛ :=(λ_k)^∞_k=0be a sequence of distinct nonnegative real numbers withλ₀ :=0 and ∑^∞_k=1 1/λ_k<∞. Let(0, 1) and(0, 1−) be fixed. An earlier work of the authors shows that [formula]is finite. In this paper an explicit upper bound forC(Λ, , ) is given. In the special caseλ_k :=k^α,α>1, our bounds are essentially sharp.     

Markov-Type Inequalities for Products of Müntz Polynomials

Let Λ(λ_j)^∞_j=0 be a sequence of distinct real numbers. The span of {x^λ₀, x^λ₁, …, x^λ_n} over is denoted by M_n(Λ)span{x^λ₀, x^λ₁, …, x^λ_n}. Elements of M_n(Λ) are called Müntz polynomials. The principal result of this paper is the following Markov-type inequality for products of Müntz polynomials. T 2.1. LetΛ(λ_j)^∞_j=0andΓ(γ_j)^∞_j=0be increasing sequences of nonnegative real numbers. Let

Full-size image

Then

Full-size image

18(n+m+1)(λ_n+γ_m).In particular ,

Full-size image

Under some necessary extra assumptions, an analog of the above Markov-type inequality is extended to the cases when the factor x is dropped, and when the interval [0, 1] is replaced by [a, b](0, ∞).     

Tractability of Multivariate Integration for Weighted Korobov Classes

We study the worst-case error of multivariate integration in weighted Korobov classes of periodic functions of d coordinates. This class is defined in terms of weights γ_j which moderate the behavior of functions with respect to successive coordinates. We study two classes of quadrature rules. They are quasi-Monte Carlo rules which use n function values and in which all quadrature weights are 1/n and rules for which all quadrature weights are non-negative. Tractability for these two classes of quadrature rules means that the minimal number of function values needed to guarantee error in the worst-case setting is bounded by a polynomial in d and ⁻¹. Strong tractability means that the bound does not depend on d and depends polynomially on ⁻¹. We prove that strong tractability holds iff ∑^∞_j=1 γ_j<∞, and tractability holds iff lim sup_d→∞ ∑^d_j=1 γ_j/log d<∞. Furthermore, strong tractability or tractability results are achieved by the relatively small class of lattice rules. We also prove that if ∑^∞_j=1 γ^1/α_j<∞, where α measures the decay of Fourier coefficients in the weighted Korobov class, then for d1, n prime and δ>0 there exist lattice rules that satisfy an error bound independent of d and of order n^−α/2+δ. This is almost the best possible result, since the order n^−α/2 cannot be improved upon even for d=1. A corresponding result is deduced for weighted non-periodic Sobolev spaces: if ∑^∞_j=1 γ^1/2_j<∞, then for d1, n prime and δ>0 there exist shifted lattice rules that satisfy an error bound independent of d and of order n^−1+δ. We also check how the randomized error of the (classical) Monte Carlo algorithm depends on d for weighted Korobov classes. It turns out that Monte Carlo is strongly tractable iff ∑^∞_j=1 log γ_j<∞ and tractable iff lim sup_d→∞ ∑^d_j=1 log γ_j/log d<∞. Hence, in particular, for γ_j=1 we have the usual Korobov space in which integration is intractable for the two classes of quadrature rules in the worst-case setting, whereas Monte Carlo is strongly tractable in the randomized setting.     

Estimates of Freud-Christoffel functions for some weights with the whole real line as support

Upper and lower bounds for generalized Christoffel functions, called Freud-Christoffel functions, are obtained. These have the form λ_n,p(W,j,x) = infPW_Lp(R)/|P^(j)(X)| where the infimum is taken over all polynomials P(x) of degree at most n − 1. The upper and lower bounds for λ_n,p(W,j,x) are obtained for all 0 < p ∞ and J = 0, 1, 2, 3,… for weights W(x) = exp(−Q(x)), where, among other things, Q(x) is bounded in [− A, A], and Q″ is continuous in β(−A, A) for some A > 0. For p = ∞, the lower bounds give a simple proof of local and global Markov-Bernstein inequalities. For p = 2, the results remove some restrictions on Q in Freud's work. The weights considered include W(x) = exp(− ¦x¦^α/2), α > 0, and W(x) = exp(− exp(¦x¦)), > 0.     

On Positive and Copositive Polynomial and Spline Approximation inLp[−1, 1], 0<p<∞

For a functionfL_p[−1, 1], 0<p<∞, with finitely many sign changes, we construct a sequence of polynomialsP_nΠ_nwhich are copositive withfand such that f−P_npCω(f, (n+1)⁻¹)_p, whereω(f, t)_pdenotes the Ditzian–Totik modulus of continuity inL_pmetric. It was shown by S. P. Zhou that this estimate is exact in the sense that if f has at least one sign change, thenωcannot be replaced byω²if 1<p<∞. In fact, we show that even for positive approximation and all 0<p<∞ the same conclusion is true. Also, some results for (co)positive spline approximation, exact in the same sense, are obtained.     

On the Uniform Approximation of a Class of Analytic Functions by Bruwier Series

For a class of analytic functions f(z) defined by Laplace–Stieltjes integrals the uniform convergence on compact subsets of the complex plane of the Bruwier series (B-series) ∑^∞_n=0 λ_n(f) , λ_n(f)=f⁽ⁿ⁾(nc)+cf⁽ⁿ⁺¹⁾(nc), generated by f(z) and the uniform approximation of the generating function f(z) by its B-series in cones |arg z|< is shown.     

Rate of Convergence of Positive Series

We investigate the rate of convergence of series of the form

The Bernstein Constant and Polynomial Interpolation at the Chebyshev Nodes

In this paper, we establish new asymptotic relations for the errors of approximation in L_p[−1,1], 0<p∞, of x^λ, λ>0, by the Lagrange interpolation polynomials at the Chebyshev nodes of the first and second kind. As a corollary, we show that the Bernstein constant

Full-size image

is finite for λ>0 and .     

Best Approximation with Linear Constraints: A Model Example

We consider best approximation in L_p( ), 1 ≤ p ≤ ∞, by means of entire functions y of exponential type subject to additional constraints Γ_j(y) = 0, j = 1, ..., K. Here Γ_j are (unbounded) linear functionals of the form Γ_j(y) = Dⁿy(s_j) − ∑ a_kD^ky(s_j) where s_j are fixed points.     

On the sequence of powers of a stochastic matrix with large exponent

We consider the class of primitive stochastic n×n matrices A, whose exponent is at least (n²−2n+2)/2+2. It is known that for such an A, the associated directed graph has cycles of just two different lengths, say k and j with k>j, and that there is an α between 0 and 1 such that the characteristic polynomial of A is λⁿ−αλ^n−j−(1−α)λ^n−k. In this paper, we prove that for any mn, if α1/2, then A^m+k−A^m_∞A^m−1w^T_∞, where 1 is the all-ones vector and w^T is the left-Perron vector for A, normalized so that w^T1=1. We also prove that if jn/2, n31 and , then A^m+j−A^m_∞A^m−1w^T_∞ for all sufficiently large m. Both of these results lead to lower bounds on the rate of convergence of the sequence A^m.     

Remarks Concerning Completeness of Translates in Function Spaces

This note explains how to translate the author's old result on cyclic vectors of the multiple shift operator into the language of completeness theorems for integer translates. This translation, together with those results, turns out to be a source for many completeness theorems. In particular, there follows the existence of functions f whose positive integer translates f(x−k), where k ₊ are complete in the spaces C^l₀( ), L^p( ), W^l_p( ), 2<p<∞, l=0, 1, …, as well as in their weighted and/or vector-valued analogues.     

Équations du type Monge-Ampère sur les variétés Riemanniennes compactes, I

Let (V_n, g) be a C^∞ compact Riemannian manifold. For a suitable function on V_n, let us consider the change of metric: g′ = g + Hess(), and the function, as a ratio of two determinants, M() = ¦g′¦ ¦g¦⁻¹. Using the method of continuity, we first solve in C^∞ the problem: Log M() = λ + ƒ, λ > 0, ƒ ε C^∞. Then, under weak hypothesis on F, we solve the general equation: Log M() = F(P, ), F in C^∞(V_n × ¦α, β¦), using a method of iteration. Our study gives rise to an interesting a priori estimate on ¦¦, which does not occur in the complex case. This estimate should enable us to solve the equation above when λ 0, providing we can overcome difficulties related to the invertibility of the linearised operator. This open question will be treated in our next article.     

Strong limit theorems for quasi-orthogonal random fields

This is a systematic and unified treatment of a variety of seemingly different strong limit problems. The main emphasis is laid on the study of the a.s. behavior of the rectangular means ζ_mn = 1/(λ₁(m) λ₂(n)) Σ_i=1^m Σ_k=1ⁿ X_ik as either max{m, n} → ∞ or min{m, n} → ∞. Here {X_ik: i, k ≥ 1} is an orthogonal or merely quasi-orthogonal random field, whereas {λ₁(m): m ≥ 1} and {λ₂(n): n ≥ 1} are nondecreasing sequences of positive numbers subject to certain growth conditions. The method applied provides the rate of convergence, as well. The sufficient conditions obtained are shown to be the best possible in general. Results on double subsequences and 1-parameter limit theorems are also included.     

On Generalized Hermite–Fejér Interpolation of Lagrange Type on the Chebyshev Nodes

For fC[−1, 1], let H_m, n(f, x) denote the (0, 1, …,anbsp;m) Hermite–Fejér (HF) interpolation polynomial of f based on the Chebyshev nodes. That is, H_m, n(f, x) is the polynomial of least degree which interpolates f(x) and has its first m derivatives vanish at each of the zeros of the nth Chebyshev polynomial of the first kind. In this paper a precise pointwise estimate for the approximation error |H_2m, n(f, x)−f(x)| is developed, and an equiconvergence result for Lagrange and (0, 1, …, 2m) HF interpolation on the Chebyshev nodes is obtained. This equiconvergence result is then used to show that a rational interpolatory process, obtained by combining the divergent Lagrange and (0, 1, …, 2m) HF interpolation methods on the Chebyshev nodes, is convergent for all fC[−1, 1].     

The Asymptotic Behavior of a Family of Sequences via Tauberian Theorems

We study the asymptotic behavior of a family of sequences defined by the following nonlinear induction relation c₀ = 1 and c_n ∑^k_{j = 1} r_jc_[n/mj] + ∑^k_{j = k + 1} r_jc_{[(n + 1)^1/mj] − 1} for n ≥ 1, where the r_j are real positive numbers and m_j are integers greater than or equal to 2. Depending on the fact that ∑^k_{j = 1} r_j is greater or lower than 1, we prove that c_n/n^α or c_n/(ln n)^α goes to some finite limit for some explicit α. Our study is based on Tauberian theorems and extends a result of Erdös et al.     

Near-best multivariate approximation by Fourier series, Chebyshev series and Chebyshev interpolation

A set of results concerning goodness of approximation and convergence in norm is given for L_∞ and L₁ approximation of multivariate functions on hypercubes. Firstly the trigonometric polynomial formed by taking a partial sum of a multivariate Fourier series and the algebraic polynomials formed either by taking a partial sum of a multivariate Chebyshev series of the first kind or by interpolating at a tensor product of Chebyshev polynomial zeros are all shown to be near-best L_∞ approximations. Secondly the trigonometric and algebraic polynomials formed by taking, respectively, a partial sum of a multivariate Fourier series and a partial sum of a multivariate Chebyshev series of the second kind are both shown to be hear-best L₁ approximations. In all the cases considered, the relative distance of a near-best approximation from a corresponding best approximation is shown to be at most of the order of Π log n_j, where n_j (j = 1, 2,…, N) are the respective degrees of approximation in the N individual variables. Moreover, convergence in the relevant norm is established for all the sequences of near-best approximations under consideration, subject to appropriate restrictions on the function space.     

Approximation of elements of a generalized Orlicz sequence space by nonlinear, singular kernels

Let l be a generalized Orlicz sequence space generated by a modular (x) = ∑_{i − 0}^∞ _i(¦t_i¦), X = (t_i), with s-convex functions _i, 0 < s 1, and let K_w,j: R₊ → R₊ for j=0,1,2,…, w ε Wwhere is an abstract set of indices. Assuming certain singularity assumptions on the nonlinear kernel K_w,j and setting T_wx = ((T_wx)_i)_{i = 0}^∞, with (T_wx)_i = ∑_{j = 0}ⁱ K_{w,i − j}(¦t_j¦) for x = (t_j), convergence results: T_wx → x in l are obtained (both modular convergence and norm convergence), with respect to a filter of subsets of the set .     

Uniform Estimates for a Class of Evolution Equations

L^∞ estimates are derived for the oscillatory integral ∫⁺₀^∞e^{−i(xλ + (1/m) tλm)}a(λ) dλ, where 2 ≤ m and (x, t) × ₊. The amplitude a(λ) can be oscillatory, e.g., a(λ) = e^{it (λ)} with (λ) a polynomial of degree ≤ m − 1, or it can be of polynomial type, e.g., a(λ) = (1 + λ)^k with 0 ≤ k ≤ (m − 2). The estimates are applied to the study of solutions of certain linear pseudodifferential equations, of the generalized Schrödinger or Airy type, and of associated semilinear equations.