Zeros of orthogonal polynomials on the real line

Let p_n(x) be the orthonormal polynomials associated to a measure dμ of compact support in . If Esupp(dμ), we show there is a δ>0 so that for all n, either p_n or p_n+1 has no zeros in (E−δ,E+δ). If E is an isolated point of supp(μ), we show there is a δ so that for all n, either p_n or p_n+1 has at most one zero in (E−δ,E+δ). We provide an example where the zeros of p_n are dense in a gap of supp(dμ).     

Estimates for the first and second derivatives of the Stieltjes polynomials

Let w_λ(x)(1−x²)^λ−1/2 and P_n^(λ) be the ultraspherical polynomials with respect to w_λ(x). Then we denote E_n+1^(λ) the Stieltjes polynomials with respect to w_λ(x) satisfyingIn this paper, we give estimates for the first and second derivatives of the Stieltjes polynomials E_n+1^(λ) and the product E_n+1^(λ)P_n^(λ) by obtaining the asymptotic differential relations. Moreover, using these differential relations we estimate the second derivatives of E_n+1^(λ)(x) and E_n+1^(λ)(x)P_n^(λ)(x) at the zeros of E_n+1^(λ)(x) and the product E_n+1^(λ)(x)P_n^(λ)(x), respectively.     

Local Accuracy for Radial Basis Function Interpolation on Finite Uniform Grids

We consider interpolation on a finite uniform grid by means of one of the radial basis functions (RBF) φ(r)=r^γ for γ>0, γ2 or φ(r)=r^γ ln r for γ2 ₊. For each positive integer N, let h=N⁻¹ and let {x_i: i =1, 2, …, (N+1)^d} be the set of vertices of the uniform grid of mesh-size h on the unit d-dimensional cube [0, 1]^d. Given f: [0, 1]^d→ , let s_h be its unique RBF interpolant at the grid vertices: s_h(x_i)=f(x_i), i=1, 2, …, (N+1)^d. For h→0, we show that the uniform norm of the error f−s_h on a compact subset K of the interior of [0, 1]^d enjoys the same rate of convergence to zero as the error of RBF interpolation on the infinite uniform grid h ^d, provided that f is a data function whose partial derivatives in the interior of [0, 1]^d up to a certain order can be extended to Lipschitz functions on [0, 1]^d.     

Analytic eigendistributions and applications to homogeneous trees

An analytic distribution on is an element, ν, of the dual of the space of analytic functions on K. In particular, ν defines a linear functional on the polynomial ring . In this work, we study the converse problem: given a linear functional on , try to find a minimal set K such that ν extends to an analytic distribution on K. This study was motivated by the desire to generalize a result that allows the representation of functions on a homogeneous tree as integrals of z-harmonic functions oven a certain interval. A function f on a homogeneous tree T of degree q+1 is said to be z-harmonic, if μ₁f=zf, where μ₁ is the nearest neighbor averaging operator. It was proved in [Cohen, Colonna, Adv. Appl. Math. 20 (1998) 253–274] that if |f(v)|MC^|v| for constants M>0 and , then there exist z-harmonic functions k_z such that where I is the interval with endpoints . In the present paper, we study the case when the above exponential growth condition holds with , which necessitates replacing k_z(v) dz with an analytic distribution ν_v satisfying the z-harmonicity condition μ₁ν=zν. We show that to each function on the tree satisfying the above exponential growth condition there corresponds an eigendistribution on an elliptical region containing I as the interval between its foci.     

Generalized fractal dimensions: equivalences and basic properties

Given a positive probability Borel measure μ on , we establish some basic properties of the associated functions τ_μ^±(q) and of the generalized fractal dimensions D_μ^±(q) for . We first give the equivalence of the Hentschel–Procaccia dimensions with the Rényi dimensions and the mean-q dimensions, for q>0. We then use these relations to prove some regularity properties for τ_μ^±(q) and D_μ^±(q); we also provide some estimates for these functions, in particular estimates on their behaviour at ±∞, as well as for the dimensions corresponding to convolution of two measures. We finally present some calculations for specific examples illustrating the different cases met in the article.     

Galerkin solution of a singular integral equation with constant coefficients

Galerkin methods are used to approximate the singular integral equation

with solution φ having weak singularity at the endpoint −1, where a, b≠0 are constants. In this case φ is decomposed as φ(x)=(1−x)^α(1+x)^βu(x), where β=−α, 0<α<1. Jacobi polynomials are used in the discussions. Under the conditions fH^μ[−1,1] and k(t,x)H^μ,μ[−1,1]×[−1,1], 0<μ<1, the error estimate under a weighted L² norm is O(n^−μ). Under the strengthened conditions f^″H^μ[−1,1] and , 2α−<μ<1, the error estimate under maximum norm is proved to be O(n^2α−−μ+), where , >0 is a small enough constant.     

On weighted Lebesgue function type sums

Let I be a finite or infinite interval and dμ a measure on I. Assume that the weight function w(x)>0, w^′(x) exists, and the function w^′(x)/w(x) is non-increasing on I. Denote by ℓ_k's the fundamental polynomials of Lagrange interpolation on a set of nodes x₁<x₂<<x_n in I. The weighted Lebesgue function type sum for 1≤i<j≤n and s≥1 is defined by

In this paper the exact lower bounds of S_n(x) on a “big set” of I and are obtained. Some applications are also given.     

Representation of Reproducing Kernels and the Lebesgue Constants on the Ball

For the weight function (1−x²)^μ−1/2 on the unit ball, a closed formula of the reproducing kernel is modified to include the case −1/2<μ<0. The new formula is used to study the orthogonal projection of the weighted L² space onto the space of polynomials of degree at most n, and it is proved that the uniform norm of the projection operator has the growth rate of n^(d−1)/2 for μ<0, which is the smallest possible growth rate among all projections, while the rate for μ0 is n^μ+(d−1)/2.     

Uniqueness and stability of regional blow-up in a porous-medium equation

We study the blow-up phenomenon for the porous-medium equation in R^N, N1, u_t=Δu^m+u^m, m>1, for nonnegative, compactly supported initial data. A solution u(x,t) to this problem blows-up at a finite time . Our main result asserts that there is a finite number of points x₁,…,x_kR^N, with |x_i−x_j|2R^* for i≠j, such that Here w_*(|x|) is the unique nontrivial, nonnegative compactly supported, radially symmetric solution of the equation in R^N and R^* is the radius of its support. Moreover u(x,t) remains uniformly bounded up to its blow-up time on compact subsets of . The question becomes reduced to that of proving that the ω-limit set in the problem consists of a single point when its initial condition is nonnegative and compactly supported.     

Best Approximation by Normal and Conormal Sets

The aim of the present paper is to develop a theory of best approximation by elements of so-called normal sets and their complements—conormal sets—in the non-negative orthant ^I₊ of a finite-dimensional coordinate space ^I endowed with the max-norm. A normal (respectively, conormal) set arises as the set of all solutions of a system of inequalities f_α(x)0 (αA), x ^I₊ (respectively, f_α(x)0 (αA), x ^I₊), where f_α is an increasing function and A is an arbitrary set of indices. We consider these sets as analogues (in a certain sense) of convex sets, and we use the so-called min-type functions as analogues of linear functions. We show that many results on best approximation by convex and reverse convex sets and corresponding separation theory (but not all of them) have analogues in the case under consideration. At the same time there are no convex analogues for many results related to best approximation by normal sets.     

Classification Theorems for General Orthogonal Polynomials on the Unit Circle

The set of all probability measures σ on the unit circle splits into three disjoint subsets depending on properties of the derived set of {|_n|²dσ}_n0, denoted by Lim(σ). Here {_n}_n0 are orthogonal polynomials in L²(dσ). The first subset is the set of Rakhmanov measures, i.e., of σ with {m}=Lim(σ), m being the normalized (m( )=1) Lebesgue measure on . The second subset Mar( ) consists of Markoff measures, i.e., of σ with mLim(σ), and is in fact the subject of study for the present paper. A measure σ, belongs to Mar( ) iff there are >0 and l>0 such that sup{|a_n+j|:0jl}>, n=0,1,2,…,{a_n} is the Geronimus parameters (=reflectioncoefficients) of σ. We use this equivalence to describe the asymptotic behavior of the zeros of the corresponding orthogonal polynomials (see Theorem G). The third subset consists of σ with {m}Lim(σ). We show that σ is ratio asymptotic iff either σ is a Rakhmanov measure or σ satisfies the López condition (which implies σMar( )). Measures σ satisfying Lim(σ)={ν} (i.e., weakly asymptotic measures) are also classified. Either ν is the sum of equal point masses placed at the roots of zⁿ=λ, λ , n=1,2,…, or ν is the equilibrium measure (with respect to the logarithmic kernel) for the inverse image under an m-preserving endomorphism z→zⁿ, n=1,2,…, of a closed arc J (including J= ) with removed open concentric arc J₀ (including J₀=). Next, weakly asymptotic measures are completely described in terms of their Geronimus parameters. Finally, we obtain explicit formulae for the parameters of the equilibrium measures ν and show that these measures satisfy {ν}=Lim(ν).     

Generating orthogonal matrix polynomials satisfying second order differential equations from a trio of triangular matrices

The method developed in [A.J. Durán, F.A. Grünbaum, Orthogonal matrix polynomials satisfying second order differential equations, Int. Math. Res. Not. 10 (2004) 461–484] led us to consider matrix polynomials that are orthogonal with respect to weight matrices W(t) of the form , , and (1−t)^α(1+t)^βT(t)T^*(t), with T satisfying T^′=(2Bt+A)T, T(0)=I, T^′=(A+B/t)T, T(1)=I, and T^′(t)=(−A/(1−t)+B/(1+t))T, T(0)=I, respectively. Here A and B are in general two non-commuting matrices. We are interested in sequences of orthogonal polynomials (P_n)_n which also satisfy a second order differential equation with differential coefficients that are matrix polynomials F₂, F₁ and F₀ (independent of n) of degrees not bigger than 2, 1 and 0 respectively. To proceed further and find situations where these second order differential equations hold, we only dealt with the case when one of the matrices A or B vanishes.The purpose of this paper is to show a method which allows us to deal with the case when A, B and F₀ are simultaneously triangularizable (but without making any commutativity assumption).     

An approximation for the complex normal probability integral

An analytic approximation to the complex normal probability integral, (x+iy)=(2)^–1/2 _–t ^x exp[–(t–iy)²/2]dt, is given together with a study of the error in the approximation.     

Orthogonal polynomials for exponential weights on , II

Let I=[0,d), where d is finite or infinite. Let W_ρ(x)=x^ρexp(-Q(x)), where and Q is continuous and increasing on I, with limit ∞ at d. We obtain further bounds on the orthonormal polynomials associated with the weight , finer spacing on their zeros, and estimates of their associated fundamental polynomials of Lagrange interpolation. In addition, we obtain weighted Markov and Bernstein inequalities.     

Transformations preserving normality and Wishart-ness

Let R^n×p, (n), Gl(p) and ₊(p) denote respectively the set of n×p matrices, the set of n×n orthogonal matrices, the set of p×p nonsingular matrices and the set of p × p positive definite matrices. In this paper, it is first shown that a bijective and bimeasurable transformation (BBT) g on R^p≡R^p×1 preserving the multivariate normality of N_p(μ, Σ) for fixed μ=μ₁, μ₂ (μ₁≠μ₂) and for all Σ ₊(p) is of the form g(x)=Ax+b a.e. for some (A, b)Gl(p)×R^p. Second, a BBT g on R^n×p preserving the form for certain 's and all Σ ₊(p) is shown to be of the form g(x)=QxA+E a.e. for some (Q, A, E) (n)×Gl(p)×R^n×p. Third, a BBT h on ₊(p) preserving the Wishart-ness of W_p(Σ, m) (m≥p) for all Σ ₊(p) is shown to be of the form h(w)=A′wA a.e. for some AGl(p). Fourth, a BBT k(x, w)=(k₁(x, w), k₂(x, w)) on R^n×p× ₊(p) which preserves the form of for certain 's and all Σ ₊(p) is shown to be of the form k(x, w)=(QxA+E, A′wA) a.e. for some (Q, A, E) (n)×Gl(p)×R^n×p.     

Bounds on Turán determinants

Let μ denote a symmetric probability measure on [−1,1] and let (p_n) be the corresponding orthogonal polynomials normalized such that p_n(1)=1. We prove that the normalized Turán determinant Δ_n(x)/(1−x²), where , is a Turán determinant of order n−1 for orthogonal polynomials with respect to . We use this to prove lower and upper bounds for the normalized Turán determinant in the interval −1<x<1.     

On a problem of H. Shapiro

Let μ be a real measure on the line such that its Poisson integral M(z) converges and satisfies|M(x+iy)|Ae^−cyα, y→+∞,for some constants A,c>0 and 0<α1. We show that for 1/2<α1 the measure μ must have many sign changes on both positive and negative rays. For 0<α1/2 this is true for at least one of the rays, and not always true for both rays. Asymptotical bounds for the number of sign changes are given which are sharp in some sense.     

Quadrature formulae connected to σ-orthogonal polynomials

Let dλ(t) be a given nonnegative measure on the real line , with compact or infinite support, for which all moments exist and are finite, and μ₀>0. Quadrature formulas of Chakalov–Popoviciu type with multiple nodes

where σ=σ_n=(s₁,s₂,…,s_n) is a given sequence of nonnegative integers, are considered. A such quadrature formula has maximum degree of exactness d_max=2∑_ν=1ⁿs_ν+2n−1 if and only if

The proof of the uniqueness of the extremal nodes τ₁,τ₂,…,τ_n was given first by Ghizzetti and Ossicini (Rend. Mat. 6(8) (1975) 1–15). Here, an alternative simple proof of the existence and the uniqueness of such quadrature formulas is presented. In a study of the error term R(f), an influence function is introduced, its relevant properties are investigated, and in certain classes of functions the error estimate is given. A numerically stable iterative procedure, with quadratic convergence, for determining the nodes τ_ν, ν=1,2,…,n, which are the zeros of the corresponding σ-orthogonal polynomial, is presented. Finally, in order to show a numerical efficiency of the proposed procedure, a few numerical examples are included.     

Estimation of dimension functions of band-limited wavelets

The dimension function D_ψ of a band-limited wavelet ψ is bounded by n if its Fourier transform is supported in [−(2ⁿ⁺²/3)π,(2ⁿ⁺²/3)π]. For each and for each , 0<<δ=δ(n), we construct a wavelet ψ with supp

Full-size image

such that D_ψ>n on a set of positive measure, which proves that [−(2ⁿ⁺²/3)π,(2ⁿ⁺²/3)π] is the largest symmetric interval for estimating the dimension function by n. This construction also provides a family of (uncountably many) wavelet sets each consisting of infinite number of intervals.     

Permutation polynomials of the form

Recently, several classes of permutation polynomials of the form (x²+x+δ)^s+x over have been discovered. They are related to Kloosterman sums. In this paper, the permutation behavior of polynomials of the form (x^p−x+δ)^s+L(x) over is investigated, where L(x) is a linearized polynomial with coefficients in . Six classes of permutation polynomials on are derived. Three classes of permutation polynomials over are also presented.