Tracking poles and representing Hankel operators directly from data

Summary We propose and analyse a method of estimating the poles near the unit circleT of a functionG whose values are given at a grid of points onT: we give an algorithm for performing this estimation and prove a convergence theorem. The method is to identify the phase for an estimate by considering the peaks of the absolute value ofG onT, and then to estimate the modulus by seeking a bestL ² fit toG over a small arc by a first order rational function. These pole estimates lead to the construction of a basis ofL ² which is well suited to the numerical representation of the Hankel operator with symbolG and thereby to the numerical solution of the Nehari problem (computing the bestH , analytic, approximation toG relative to theL norm), as analysed in [HY]. We present the results of numerical tests of these algorithms.Partially supported by grants from the AFOSR and NSF     

Some Landau's type inequalities for infinitesimal generators

Summary Lett T(t) be a strongly continuous contraction semigroup on a complex Banach space and letA be its infinitesimal generator. We prove that, forx D(A ³), the following inequalities hold true: Ax³ 243/8 x²A ³ x, A ² x 24 xA ³ x². Ift T(t) is a contraction group (resp. cosine function) we get the analogous but better inequalities with constants 9/8 and 3 (resp. 81/40 and 72/25) instead of 243/8 and 24. We consider also uniformly bounded semigroups, groups and cosine functions.     

A strange pigeon-hole principle

Using Ramsey theory, we establish the following pigeon-hole type principle: From a large number of random variables (functions, vectors, etc.) one can always select two, X and Y, such that P(X < Y) 1/2. We apply the principle for a poset problem.     

The umbral calculus and orthogonal polynomials

We determine all orthogonal polynomials having Boas-Buck generating functions g(t)(xf(t)), where% MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGceaqabeaacqqHOo% qwcaGGOaGaamiDaiaacMcacqGH9aqpruqqYLwySbacfaGaa8hiamaa% BeaaleaacaaIWaaabeaakiaadAeacaqGGaWaaSbaaSqaaiaabgdaae% qaaOGaaeikaiaadggacaGGSaGaa8hiaiaadshacaqGPaGaaeilaiaa% bccacaqGGaGaaeiiaiaadggacqGHGjsUcaaIWaGaaiilaiaa-bcacq% GHsislcaaIXaGaaiilaiaa-bcacqGHsislcaaIYaGaaiilaiablAci% ljaacUdaaeaacqqHOoqwcaGGOaGaamiDaiaacMcacqGH9aqpcaWFGa% WaaSraaSqaaiaaicdaaeqaaOGaamOraiaabccadaWgaaWcbaGaaeOm% aaqabaGccaGGOaWaaSqaaSqaaiaaigdaaeaacaaIZaaaaOGaaiilai% aa-bcadaWcbaWcbaGaaGOmaaqaaiaaiodaaaGccaGGSaGaa8hiaiaa% dshacaGGPaGaa8hiamaaBeaaleaacaaIWaaabeaakiaadAeacaqGGa% WaaSbaaSqaaiaabkdaaeqaaOGaaeikamaaleaaleaacaaIYaaabaGa% aG4maaaakiaacYcacaWFGaWaaSqaaSqaaiaaisdaaeaacaaIZaaaaO% Gaaiilaiaa-bcacaWG0bGaaiykaiaacYcacaWFGaWaaSraaSqaaiaa% icdaaeqaaOGaamOraiaabccadaWgaaWcbaGaaeOmaaqabaGccaGGOa% WaaSqaaSqaaiaaisdaaeaacaaIZaaaaOGaaiilaiaa-bcadaWcbaWc% baGaaGynaaqaaiaaiodaaaGccaGGSaGaa8hiaiaadshacaGGPaGaai% 4oaaqaaiabfI6azjaacIcacaWG0bGaaiykaiabg2da9iaa-bcadaWg% baWcbaGaaGimaaqabaGccaWGgbGaaeiiamaaBaaaleaacaqGZaaabe% aakiaacIcadaWcbaWcbaGaaGymaaqaaiaaisdaaaGccaGGSaGaa8hi% amaaleaaleaacaaIYaaabaGaaGinaaaakiaacYcacaWFGaWaaSqaaS% qaaiaaiodaaeaacaaI0aaaaOGaaiilaiaa-bcacaWG0bGaaiykaiaa% -bcadaWgbaWcbaGaaGimaaqabaGccaWGgbGaaeiiamaaBaaaleaaca% qGZaaabeaakiaabIcadaWcbaWcbaGaaGOmaaqaaiaaisdaaaGccaGG% SaGaa8hiamaaleaaleaacaaIZaaabaGaaGinaaaakiaacYcacaWFGa% WaaSqaaSqaaiaaiwdaaeaacaaI0aaaaOGaaiilaiaa-bcacaWG0bGa% aiykaiaacYcaaeaadaWgbaWcbaGaaGimaaqabaGccaWGgbGaaeiiam% aaBaaaleaacaqGZaaabeaakiaacIcadaWcbaWcbaGaaG4maaqaaiaa% isdaaaGccaGGSaGaa8hiamaaleaaleaacaaI1aaabaGaaGinaaaaki% aacYcacaWFGaWaaSqaaSqaaiaaiAdaaeaacaaI0aaaaOGaaiilaiaa% -bcacaWG0bGaaiykaiaacYcacaGGUaGaa8hiamaaBeaaleaacaaIWa% aabeaakiaadAeacaqGGaWaaSbaaSqaaiaabodaaeqaaOGaaeikamaa% leaaleaacaaI1aaabaGaaGinaaaakiaacYcacaWFGaWaaSqaaSqaai% aaiAdaaeaacaaI0aaaaOGaaiilaiaa-bcadaWcbaWcbaGaaG4naaqa% aiaaisdaaaGccaGGSaGaa8hiaiaadshacaGGPaGaaiOlaaaaaa!C1F3!\[\begin{gathered}\Psi (t) = {}_0F{\text{ }}_{\text{1}} {\text{(}}a, t{\text{), }}a \ne 0, - 1, - 2, \ldots ; \hfill \\\Psi (t) = {}_0F{\text{ }}_{\text{2}} (\tfrac{1}{3}, \tfrac{2}{3}, t) {}_0F{\text{ }}_{\text{2}} {\text{(}}\tfrac{2}{3}, \tfrac{4}{3}, t), {}_0F{\text{ }}_{\text{2}} (\tfrac{4}{3}, \tfrac{5}{3}, t); \hfill \\\Psi (t) = {}_0F{\text{ }}_{\text{3}} (\tfrac{1}{4}, \tfrac{2}{4}, \tfrac{3}{4}, t) {}_0F{\text{ }}_{\text{3}} {\text{(}}\tfrac{2}{4}, \tfrac{3}{4}, \tfrac{5}{4}, t), \hfill \\{}_0F{\text{ }}_{\text{3}} (\tfrac{3}{4}, \tfrac{5}{4}, \tfrac{6}{4}, t),. {}_0F{\text{ }}_{\text{3}} {\text{(}}\tfrac{5}{4}, \tfrac{6}{4}, \tfrac{7}{4}, t). \hfill \\\end{gathered}\]We also determine all Sheffer polynomials which are orthogonal on the unit circle. The formula for the product of polynomials of the Boas-Buck type is obtained.     

Quadratic spline interpolation on uniform meshes

The interpolation problem at uniform mesh points of a quadratic splines(x _i)=f _i,i=0, 1,...,N ands(x ₀)=f₀ is considered. It is known that s–f=O(h ³) and s–f=O(h ²), whereh is the step size, and that these orders cannot be improved. Contrary to recently published results we prove that superconvergence cannot occur for any particular point independent off other than mesh points wheres=f by assumption. Best error bounds for some compound formulae approximatingf _i andf _i ⁽³⁾ are also derived.     

Density estimation and adaptive control of markov processes: Average and discounted criteria

We consider a class of discrete-time Markov control processes with Borel state and action spaces, and ^d i.i.d. disturbances with unknown distribution . Under mild semi-continuity and compactness conditions, and assuming that is absolutely continuous with respect to Lebesgue measure, we establish the existence of adaptive control policies which are (1) optimal for the average-reward criterion, and (2) asymptotically optimal in the discounted case. Our results are obtained by taking advantage of some well-known facts in the theory of density estimation. This approach allows us to avoid restrictive conditions on the state space and/or on the system's transition law imposed in recent works, and on the other hand, it clearly shows the way to other applications of nonparametric (density) estimation to adaptive control.Research partially supported by The Third World Academy of Sciences under Research Grant No. MP 898-152.     

On well-ordered mono-unary algebras

Given a linearly ordered set (A, _R ) and an R-monotone function f: AA, we give a necessary and sufficient condition on A, f, _R , involving generating sets and forbidden subalgebras, for _R to be a well-ordering.Partially supported by Hungarian National Foundation for Scientific Research Grant nr. 1813.     

Error for the modified secant method

New error bounds for the modified secant method are provided. We show that our error estimates are better than the ones already existing in the literature, under similar assumptions.     

Parallel and superfast algorithms for Hankel systems of equations

Summary Utilizing kernel structure properties a unified construction of Hankel matrix inversion algorithms is presented. Three types of algorithms are obtained: 1)O(n ²) complexity Levinson type, 2)O (n) parallel complexity Schur-type, and 3)O(n log² n) complexity asymptotically fast ones. All algorithms work without additional assumption (like strong nonsingularity).     

Computing singular solutions to nonlinear analytic systems

Summary A method to generate an accurate approximation to a singular solution of a system of complex analytic equations is presented. Since manyreal systems extend naturally tocomplex analytic systems, this porvides a method for generating approximations to singular solutions to real systems. Examples include systems of polynomials and systems made up of trigonometric, exponential, and polynomial terms. The theorem on which the method is based is proven using results from several complex variables. No special conditions on the derivatives of the system, such as restrictions on the rank of the Jacobian matrix at the solution, are required. The numerical method itself is developed from techniques of homotopy continuation and 1-dimensional quadrature. A specific implementation is given, and the results of numerical experiments in solving five test problems are presented.