A new approach to Euler splines

Starting from the exponential Euler polynomials discussed by Euler in “Institutions Calculi Differentialis,” Vol. II, 1755, the author introduced in “Linear operators and approximation,” Vol. 20, 1972, the so-called exponential Euler splines. Here we describe a new approach to these splines. Let t be a constant such that t=|t|e^iα, −π<α<π,t≠0,t≠1.. Let S₁(x:t) be the cardinal linear spline such that S₁(v:t) = t^v for all v ε Z. Starting from S₁(x:t) it is shown that we obtain all higher degree exponential Euler splines recursively by the averaging operation . Here S_n(x:t) is a cardinal spline of degree n if n is odd, while is a cardinal spline if n is even. It is shown that they have the properties S_n(v:t) = t^v for v ε Z.     

Permanence and global stability in a Lotka-Volterra predator-prey system with delays

Consider the permanence and global asymptotic stability of models governed by the following Lotka-Volterra-type system:

, with initial conditions

x_i(t) = φ_i(t) ≥ o, t ≤ t₀, and φ_i(t₀) > 0. 1 ≤ i ≤n

. We define x₀(t) = x_n+1(t)≡0 and suppose that φ_i(t), 1 ≤ i ≤ n, are bounded continuous functions on [t₀, + ∞) and γ_i, α_i, c_i > 0,γ_i,j ≥ 0, for all relevant i,j.Extending a technique of Saito, Hara and Ma[1] for n = 2 to the above system for n ≥ 2, we offer sufficient conditions for permanence and global asymptotic stability of the solutions which improve the well-known result of Gopalsamy.     

An analytic approximation to the cardinal functions of Gaussian radial basis functions on an infinite lattice

Gaussian radial basis functions (RBFs) have been very useful in computer graphics and for numerical solutions of partial differential equations where these RBFs are defined, on a grid with uniform spacing h, as translates of the “master” function (x;α,h)≡exp(-[α²/h²]x²) where α is a user-choosable constant. Unfortunately, computing the coefficients of (x-jh;α,h) requires solving a linear system with a dense matrix. It would be much more efficient to rearrange the basis functions into the equivalent “Lagrangian” or “cardinal” basis because the interpolation matrix in the new basis is the identity matrix; the cardinal basis C_j(x;α,h) is defined by the set of linear combinations of the Gaussians such that C_j(kh)=1 when k=j and C_j(kh)=0 for all integers . We show that the cardinal functions for the uniform grid are C_j(x;h,α)=C(x/h-j;α) where C(X;α)≈(α²/π)sin(πX)/sinh(α²X). The relative error is only about 4exp(-2π²/α²) as demonstrated by the explicit second order approximation. It has long been known that the error in a series of Gaussian RBFs does not converge to zero for fixed α as h→0, but only to an “error saturation” proportional to exp(-π²/α²). Because the error in our approximation to the master cardinal function C(X;α) is the square of the error saturation, there is no penalty for using our new approximations to obtain matrix-free interpolating RBF approximations to an arbitrary function f(x). The master cardinal function on a uniform grid in d dimensions is just the direct product of the one-dimensional cardinal functions. Thus in two dimensions . We show that the matrix-free interpolation can be extended to non-uniform grids by a smooth change of coordinates.     

Properties of Hida processes on . 2. Prediction and interpolation problems for processes on

If E is an ordered set, we study the processes Y_t, t E, for which the vectorial spaces _t generated by all the conditional expectations E(Y_sβ _t) for s ≥ t have finite dimensions d(t) ≤ N. ( _t is some convenient filtration.) We first develop a geometrical approach in the general situation and give a “Goursat's representation” Y_t = Σf_i(t)M_i(t), where the M_i(t) are martingales. We then restrict us to the cases E = or E = ² and give representations of the processes by the mean of stochastic integrals of “Goursat's kernels.” The special case when Y_t is the solution of a differential equation is considered.     

Unbounded divergence of simple quadrature formulas

Given a sequence of real or complex coefficients c_i and a sequence of distinct nodes t_i in a compact interval T, we prove the divergence and the unbounded divergence on superdense sets in the space C(T) of the simple quadrature formulas ∝_Tx(t)du(t) = Q_n(x) + R_n(x) and ∝_Tw(t)x(t)dt = Q_n(x) + R_n(x), where Q_n(x)=∑_i=1^m_n c_ix(t_i), ε C(T).The divergence (not certainly unbounded) for at most one continuous function of the first simple quadrature formula, with m_n = n and u(t) = t, was established by P. J. Davis in 1953.     

Instability of Certain Nonlinear Stochastic Differential Equations

We prove existence and uniqueness of the solution X^ε_t of the SDE, X^ε_t = εB_t + ∫^t₀u^{q −1} _ε(s, X^ε_t) ds, where X^ε_t is a one-dimensional process and u_ε(t, x) the density of X^ε_t (ε > 0, q > 1). We show that the closure of (X^ε_t; 0 ≤ t ≤ 1) with respect to Hölder norm, when ε goes to 0, is a.s. equal to an explicit family of continuous functions. We obtain similar results, considering SDE′s where the drift coefficient is equal to ± sgn(x) u(t, x).     

Persistence, contractivity and global stability in logistic equations with piecewise constant delays

We establish sufficient conditions for the persistence and the contractivity of solutions and the global asymptotic stability for the positive equilibrium N^*=1/(a+∑_i=0^mb_i) of the following differential equation with piecewise constant arguments:

where r(t) is a nonnegative continuous function on [0,+∞), r(t)0, ∑_i=0^mb_i>0, b_i0, i=0,1,2,…,m, and a+∑_i=0^mb_i>0. These new conditions depend on a,b₀ and ∑_i=1^mb_i, and hence these are other type conditions than those given by So and Yu (Hokkaido Math. J. 24 (1995) 269–286) and others. In particular, in the case m=0 and r(t)≡r>0, we offer necessary and sufficient conditions for the persistence and contractivity of solutions. We also investigate the following differential equation with nonlinear delay terms:

where r(t) is a nonnegative continuous function on [0,+∞), r(t)0, 1−ax−g(x,x,…,x)=0 has a unique solution x^*>0 and g(x₀,x₁,…,x_m)C¹[(0,+∞)×(0,+∞)××(0,+∞)].     

Stepanov-like pseudo-almost periodic mild solutions to perturbed nonautonomous evolution equations with infinite delay

In this paper, we prove the invariance of Stepanov-like pseudo-almost periodic functions under bounded linear operators. Furthermore, we obtain existence and uniqueness theorems of pseudo-almost periodic mild solutions to evolution equations u^′(t)=A(t)u(t)+h(t) and on , assuming that A(t) satisfy “Acquistapace–Terreni” conditions, that the evolution family generated by A(t) has exponential dichotomy, that R(λ₀,A()) is almost periodic, that B,C(t,s)_t≥s are bounded linear operators, that f is Lipschitz with respect to the second argument uniformly in the first argument and that h, f, F are Stepanov-like pseudo-almost periodic for p>1 and continuous. To illustrate our abstract result, a concrete example is given.     

Comonotone Jackson's Inequality

Let 2s points y_i=−πy_2s<…<y₁<π be given. Using these points, we define the points y_i for all integer indices i by the equality y_i=y_i+2s+2π. We shall write fΔ⁽¹⁾(Y) if f is a 2π-periodic continuous function and f does not decrease on [y_i, y_i−1], if i is odd; and f does not increase on [y_i, y_i−1], if i is even. In this article the following Theorem 1—the comonotone analogue of Jackson's inequality—is proved. 1. If fΔ⁽¹⁾(Y), then for each nonnegative integer n there is a trigonometric polynomial τ_n(x) of order n such that τ_nΔ⁽¹⁾(Y), and |f(x)−π_n(x)|c(s) ω(f; 1/(n+1)), x , where ω(f; t) is the modulus of continuity of f, c(s)=const. Depending only on s.     

Stability conditions for systems of linear nonautonomous delay differential equations

Systems of linear nonautonomous delay differential equations are considered which are of the form y′_i(t) = ∑_{k = 1}ⁿ ∝₀^T b_ik(t, s) y_k(t − s) dη_ik(s) − c_i(t) y_i(t), where I = 1,…, n. Sufficient conditions are derived for both the asymptotic stability and the instability of the zero solution. The main result is found by a monotone technique using elementary methods only. Moreover, additional criteria are obtained by using the method of Lyapunov functionals.     

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 .     

On sharp Burkholder–Rosenthal-type inequalities for infinite-degree U-statistics

In this paper, we present a method that allows one to obtain a number of sharp inequalities for expectations of functions of infinite-degree U-statistics. Using the approach, we prove, in particular, the following result: Let D be the class of functions f :R₊→R₊ such that the function f(x+z)−f(x) is concave in xR₊ for all zR₊. Then the following estimate holds: for all fD and all U-statistics ∑_1i1<<ilnY_i1,…,il(X_i1,…,X_il) with nonnegative kernels Y_i1,…,il :R^l→R₊, 1i_kn; i_r≠i_s, r≠s; k,r,s=1,…,l; l=0,…,m, in independent r.v.'s X₁,…,X_n. Similar inequality holds for sums of decoupled U-statistics. The class D is quite wide and includes all nonnegative twice differentiable functions f such that the function f″(x) is nonincreasing in x>0, and, in particular, the power functions f(x)=x^t, 1<t2; the power functions multiplied by logarithm f(x)= (x+x₀)^t ln(x+x₀), 1<t<2, x₀max(e^{(3t2−6t+2)/(t(t−1)(2−t))},1); and the entropy-type functions f(x)=(x+x₀)ln(x+x₀), x₀1. As an application of the results, we determine the best constants in Burkholder–Rosenthal-type inequalities for sums of U-statistics and prove new decoupling inequalities for those objects. The results obtained in the paper are, to our knowledge, the first known results on the best constants in sharp moment estimates for U-statistics of a general type.     

Uniform approximation by some Hermite interpolating splines

In this paper some upper bound for the error ∥ s-f ∥_∞ is given, where f ε C¹[a,b], but s is a so-called Hermite spline interpolant (HSI) of degree 2q ?1 such that f(x_i) = s(x_i), f′(rmx_i) = s′(x_i), s^(j) (x_i) = 0 (i = 0, 1, …, n; j = 2, 3, …, q ?1; n > 0, q > 0) and the knots x_i are such that a = x₀ < x₁ < … < x_n = b. Necessary and sufficient conditions for the existence of convex HSI are given and upper error bound for approximation of the function fε C¹[a, b] by convex HSI is also given.     

Compact causal data interpolation

Consider a Hilbert space equipped with a time-structure, i.e., a resolution E of the identity on defined on subsets of some linearly ordered set Λ. For which x and y in is it possible to find a causal (time respecting) compact operator T, so that Tx = y? When T is required to be a Hilbert-Schmidt operator and (Λ, E) is sufficiently regular, this question is answered in terms of the “time-densities” of x and y. The condition is that the integral ∝_gLμ_x({s t})⁻¹ dμ_y(t) should be finite, where μ_x and μ_y are the measures on Λ given by μ_x(Ω) = ¦|E(Ω)x¦|² and μ_y(Ω) = ¦|E(Ω)y¦|². Further a solution is given for the related problem of minimizing the sum of ¦|Tx − y¦|² and the squared Hilbert-Schmidt norm ¦|R¦|₂² of T.     

On the Rate of Almost Everywhere Convergence of Certain Classical Integral Means

The main purpose of this article is to establish nearly optimal results concerning the rate of almost everywhere convergence of the Gauss–Weierstrass, Abel–Poisson, and Bochner–Riesz means of the one-dimensional Fourier integral. A typical result for these means is the following: If the function f belongs to the Besov spaceB^s_{p, p}, 1<p<∞, 0<s<1, thenT_mtf (x)−f(x)=o_x(t^s) a.e. ast→0+.     

Schur functions and the invariant polynomials characterizing U(n) tensor operators

We give a direct formulation of the invariant polynomials _μG_q⁽ⁿ⁾(, Δ_i,;, x_{i,i + 1},) characterizing U(n) tensor operators p, q, …, q, 0, …, 0 in terms of the symmetric functions S_λ known as Schur functions. To this end, we show after the change of variables Δ_i = γ_i − δ_i and x_{i, i + 1} = δ_i − δ_{i + 1} that_μG_q⁽ⁿ⁾(,Δ_i;, x_{i, i + 1},) becomes an integral linear combination of products of Schur functions S_α(, γ_i,) · S_β(, δ_i,) in the variables {γ₁,…, γ_n} and {δ₁,…, δ_n}, respectively. That is, we give a direct proof that _μG_q⁽ⁿ⁾(,Δ_i,;, x_{i, i + 1},) is a bisymmetric polynomial with integer coefficients in the variables {γ₁,…, γ_n} and {δ₁,…, δ_n}. By making further use of basic properties of Schur functions such as the Littlewood-Richardson rule, we prove several remarkable new symmetries for the yet more general bisymmetric polynomials _μ^mG_q⁽ⁿ⁾(γ₁,…, γ_n; δ₁,…, δ_m). These new symmetries enable us to give an explicit formula for both _μ^mG₁⁽ⁿ⁾(γ; δ) and ₁G₂⁽ⁿ⁾(γ; δ). In addition, we describe both algebraic and numerical integration methods for deriving general polynomial formulas for _μ^mG_q⁽ⁿ⁾(γ; δ).     

Admissibility of the natural estimator of the mean of a Gaussian process

Let {X(t): t [a, b]} be a Gaussian process with mean μ L₂[a, b] and continuous covariance K(s, t). When estimating μ under the loss ∫_a^b ( (t)−μ(t))² dt the natural estimator X is admissible if K is unknown. If K is known, X is minimax with risk ∫_a^b K(t, t) dt and admissible if and only if the three by three matrix whose entries are K(t_i, t_j) has a determinant which vanishes identically in t_i [a, b], i = 1, 2, 3.     

I.i.d. representations for the bivariate product limit estimators and the bootstrap versions

Let F(s, t) = P(X > s, Y > t) be the bivariate survival function which is subject to random censoring. Let be the bivariate product limit estimator (PL-estimator) by Campbell and Földes (1982, Proceedings International Colloquium on Non-parametric Statistical Inference, Budapest 1980, North-Holland, Amsterdam). In this paper, it was shown that

Full-size image

, where {ζ_i(s, t)} is i.i.d. mean zero process and R_n(s, t) is of the order O((n⁻¹log n)^3/4) a.s. uniformly on compact sets. Weak convergence of the process {n⁻¹ Σ_{i = 1}ⁿ ζ_i(s, t)} to a two-dimensional-time Gaussian process is shown. The covariance structure of the limiting Gaussian process is also given. Corresponding results are also derived for the bootstrap estimators. The result can be extended to the multivariate cases and are extensions of the univariate case of Lo and Singh (1986, Probab. Theory Relat. Fields, 71, 455–465). The estimator is also modified so that the modified estimator is closer to the true survival function than in supnorm.     

On nonautonomous linear systems of differential and difference equations with -symmetric coefficient matrices

Let denote the set of continuous n×n matrices on an interval . We say that is a nontrivial k-involution if where ζ=e^-2πi/k, d₀+d₁++d_k-1=n, and with . We say that is R-symmetric if R(t)A(t)R^-1(t)=A(t), , and we show that if A is R-symmetric then solving x^′=A(t)x or x^′=A(t)x+f(t) reduces to solving k independent d_ℓ×d_ℓ systems, 0ℓk-1. We consider the asymptotic behavior of the solutions in the case where . Finally, we sketch analogous results for linear systems of difference equations.     

Kolmogorov and Linear Widths of Weighted Sobolev-Type Classes on a Finite Interval, II

Let I be a finite interval, r and ρ(t)=dist{t, ∂I}, tI. Denote by Δ^s₊W^r_p, α, 0α<∞, the class of functions x on I with the seminorm x^(r)ρ^α_Lp1 for which Δ^s_τx, τ>0, is nonnegative on I. We obtain two-sided estimates of the Kolmogorov widths d_n(Δ^s₊W^r_p, α)_Lq and of the linear widths d_n(Δ^s₊W^r_p, α)^lin_Lq, s=0, 1, …, r+1.