On extremal perturbations of semi-Fredholm operators

Given a bounded linear operatorA in an infinite dimensional Banach space and a compact subset of a connected component of its semi-Fredholm domain, we construct a finite rank operatorF such that –A+F is bounded below (or surjective) for each ,F ²=0 and rankF=max min{dimN(–A), codimR(–A)}, if ind(–A)0 (or ind(–A)0, respectively) for each .     

Reconstruction of a rational nonsquare matrix function from local data

In many problems the local zero-pole structure (i.e. locations of zeros and poles together with their orders) of a scalar rational functionw is a key piece of structure. Knowledge of the order of the pole or zero of the rational functionw at the point is equivalent to knowledge of the -module (where is the space of rational functions analytic at ). For the more intricate case of a rationalp×m matrix functionW, we consider the structure of the module as the appropriate analogue of zero-pole structure (location of zeros and poles together with directional information), where is the set of column vectors of heightm with entries equal to rational functions which are analytic at . Modules of the form in turn can be explicitly parametrized in terms of a collection of matrices (C ,A ,B ,B , ) together with a certain row-reduced(p–m)×m matrix polynomialP(z) (which is independent of ) which satisfy certain normalization and consistency conditions. We therefore define the collection (C ,A ,Z ,B , ,P(z)) to be the local spectral data set of the rational matrix functionW at . We discuss the direct problem of how to compute the local spectral data explicitly from a realizationW(z)=D+C(z–A) ^–1 B forW and solve the inverse problem of classifying which collections (C ,A ,Z ,B , ,P(z)) satisfying the local consistency and normalization conditions arise as the local spectral data sets of some rational matrix functionW. Earlier work in the literature handles the case whereW is square with nonzero determinant.     

Factorization of a selfadjoint nonanalytic operator function II. Single spectral zone

We consider a selfadjoint and smooth enough operator-valued functionL() on the segment [a, b]. LetL(a)0,L(b)0, and there exist two positive numbers and such that the inequality |(L()f, f)|< ([a, b] f=1) implies the inequality (L'()f, f)>. Then the functionL() admits a factorizationL()=M()(I-Z) whereM() is a continuous and invertible on [a, b] operator-valued function, and operatorZ is similar to a selfadjoint one. This result was obtained in the first part of the present paper [10] under a stronge conditionL()0 ( [a,b]). For analytic functionL() the result of this paper was obtained in [13].     

Some generalizations of the classical moment problem

The article is devoted to two generalizations of the classical power moment problem, namely: 1) instead of representing the moment sequence by ⁿ , a representation by polynomialsP _n (), ¹, connected with a Jacobi matrix, appears; 2) in the representation, instead of ⁿ , the expression ⁿ figures, where is a real generalized function (i.e., we investigate some infinite-dimensional moment problem).The work is partially supported by the DFG, Project 436 UKR 113/39/0 and by the CRDF, Project UM1-2090.     

A sandwich theorem and Hyers—Ulam stability of affine functions

It is shown that two real functionsf andg, defined on a real intervalI, satisfy the inequalitiesf(x + (1 – )y) g(x) + (1 – )g(y) andg(x + (1 – )y) f(x) + (1 – )f(y) for allx, y I and [0, 1], iff there exists an affine functionh: I such thatf h g. As a consequence we obtain a stability result of Hyers—Ulam type for affine functions.     

J-pseudo-spectral andJ-inner-pseudo-outer factorizations for matrix polynomials

For a comonic polynomialL() and a selfadjoint invertible matrixJ the following two factorization problems are considered: firstly, we parametrize all comonic polynomialsR() such that . Secondly, if it exists, we give theJ-innerpseudo-outer factorizationL()=()R(), where () isJ-inner andR() is a comonic pseudo-outer polynomial. We shall also consider these problems with additional restrictions on the pole structure and/or zero structure ofR(). The analysis of these problems is based on the solution of a general inverse spectral problem for rational matrix functions, which consists of finding the set of rational matrix functions for which two given pairs are extensions of their pole and zero pair, respectively.The work of this author was supported by the USA-Israel Binational Science Foundation (BSF) Grant no. 9400271.     

The Padé-Rayleigh-Ritz method for solving large hermitian eigenproblems

We make use of the Padé approximants and the Krylov sequencex, Ax,,...,A ^m–1 x in the projection methods to compute a few Ritz values of a large hermitian matrixA of ordern. This process consists in approaching the poles ofR _x()=((I–A)^–1 x,x), the mean value of the resolvant ofA, by those of [m–1/m]_Rx(), where [m–1/m]_Rx() is the Padé approximant of orderm of the functionR _x(). This is equivalent to approaching some eigenvalues ofA by the roots of the polynomial of degreem of the denominator of [m–1/m]_Rx(). This projection method, called the Padé-Rayleigh-Ritz (PRR) method, provides a simple way to determine the minimum polynomial ofx in the Krylov subspace methods for the symmetrical case. The numerical stability of the PRR method can be ensured if the projection subspacem is sufficiently small. The mainly expensive portion of this method is its projection phase, which is composed of the matrix-vector multiplications and, consequently, is well suited for parallel computing. This is also true when the matrices are sparse, as recently demonstrated, especially on massively parallel machines. This paper points out a relationship between the PRR and Lanczos methods and presents a theoretical comparison between them with regard to stability and parallelism. We then try to justify the use of this method under some assumptions.     

Isolated point of spectrum ofP-hyponormal, log-hyponormal operators

LetT B(H) be a bounded linear operator on a complex Hilbert spaceH. Let ₀ (T) be an isolated point of (T) and let be the Riesz idempotent for ₀. In this paper, we prove that ifT isp-hyponormal or log-hyponormal, thenE is self-adjoint andE H=ker(H–₀)=ker(H–₀ ^*.This research was supported by Grant-in-Aid Research 1 No. 12640187.     

Minimum Divergence Estimators Based on Grouped Data

The paper considers statistical models with real-valued observations i.i.d. by F(x, ₀) from a family of distribution functions (F(x, ); ), R ^s , s 1. For random quantizations defined by sample quantiles (F _n ^–1 (₁),, F _n ^–1 ( _m–1)) of arbitrary fixed orders 0 < ₁ < _m-1 < 1, there are studied estimators _,n of ₀ which minimize -divergences of the theoretical and empirical probabilities. Under an appropriate regularity, all these estimators are shown to be as efficient (first order, in the sense of Rao) as the MLE in the model quantified nonrandomly by (F ^–1 (₁,₀),, F ^–1 ( _m–1, ₀)). Moreover, the Fisher information matrix I _m (₀, ) of the latter model with the equidistant orders = ( _j = j/m : 1 j m – 1) arbitrarily closely approximates the Fisher information J(₀) of the original model when m is appropriately large. Thus the random binning by a large number of quantiles of equidistant orders leads to appropriate estimates of the above considered type.     

Compact perturbation of definite type spectra of self-adjoint quadratic operator pencils

Self-adjoint quadratic operator pencilsL()= ² A + B + C with a noninvertible leading operatorA are considered. In particular, a characterization of the spectral points of positive and of negative type ofL is given, and their behavior under a compact perturbation is studied. These results are applied to a pencil arising in magnetohydrodynamics.     

A variational principle for eigenvalues of pencils of Hermitian matrices

LetH=(A, B) be a pair of HermitianN×N matrices. A complex number is an eigenvalue ofH ifdet(A–B)=0 (we include = ifdetB=0). For nonsingularH (i.e., for which some is not an eigenvalue), we show precisely which eigenvalues can be characterized as _k ⁺ =sup{inf{*A:*B=1,S},SS _k},S _k being the set of subspaces of C ^N of codimensionk–1.Dedicated to the memory of our friend and colleague Branko NajmanResearch supported by NSERC of Canada and the I.W.Killam FoundationProfessor Najman died suddenly while this work was at its final stage. His research was supported by the Ministry of Science of CroatiaResearch supported by NSERC of Canada     

Shift invariant spaces inL 2

If is a complex, separable Hilbert space, letL ² () denote theL ²-space of functions defined on the unit circle and having values in . The bilateral shift onL ²() is the operator (U f)()=f(). A Hilbert spaceH iscontractively contained in the Hilbert spaceK ifHK and the inclusion mapH–K is a contraction. We describe the structure of those Hilbert spaces, contractively contained inL ²(), that are carried into themselves contractively byU . We also do this for the subcase of those spaces which are carried into themselves unitarily byU .     

Diagonals of positive semigroups

LetE be a Dedekind complete complex Banach lattice and letD denote the diagonal projection from the spaceL _r (E) onto the centerZ(E) ofE. Let {T(t)} _t0 be a positive strongly continuous semigroup of linear operators with generatorA. The first main result is that if the spectral bounds(A) equals to zero, then the functionD(T(t)) is a center valuedp-function. The second main result is that if for >0 the diagonalD(R(, A)) of the resolvent operatorR(, A) is strictly positive, then (D(R(, A))) ^–1 is a center valued Bernstein function. As an application of these results it follows that the order limit lim₀D(R(,A)) exists inZ(E) and equals the order limit lim _m D((R(, A)) ^m ) for any >0.     

Eigenvalues of integral operators onL 2(I) given by analytic kernels

Let {_n} _n=0 be the eigenvalue sequence of a symmetric Hilbert-Schmidt operator onL ₂(I). WhenI is an open interval, a necessary condition for {_n} _n=0 to be in the sequence space is obtained. WhenI is a closed bounded interval, sufficient conditions for {_n} _n=0 to be in the sequence space ^– are obtained.     

Conditions for eigenvectors of a selfadjoint operator-function to form a basis

Consider a functionL() defined on an interval of the real axis whose values are linear bounded selfadjoint operators in a Hilbert spaceH. A point ₀ and a vector ₀ H( ₀ 0) are called eigenvalue and eigenvector ofL() ifL() ifL(₀) ₀ = 0. Supposing that the functionL() can be represented as an absolutely convergent Fourier integral, the interval is sufficiently small and the derivative will be positive at some point, it has been proved that all the eigenvectors of the operator-functionL() corresponding to the eigenvalues from the interval form an unconditional basis in the subspace spanned by them.     

Separation of roots of matrix and operator polynomials

In a complex Hilbert spaceX for an arbitrary operator polynomialL() ( C) of degreem the following theorem is proved. If the equation (L()x, x)=0 hasm distinct roots at every pointx X, x=1, then there existm pairwise disjoint connected sets in C such that each set contains a root at everyx. The minimal distance between the roots is separated from zero under the same assumption on the discriminant and the leading coefficient of that equation.     

On eigenvalues of perturbed quadratic matrix polynomials

Among other results, it is shown that ifC andK are arbitrary complexn×n matrices and if det( ₀ ² I₀ C+K)=0 for some ₀0 (resp. ₀=0), then the Newton diagram of the polynomialt(, ) = det(² I+(1+)C+K expanded in (–₀) and , has at least a point on or below the linex+y=b (resp. has no expanded in (–₀) and , has at least a point on or below the of 0 as an eigenvalue of ₀ ² I+₀ C+K. These are extensions of similar results deu to H. Langer, B. Najman, and K. Veseli proved for diagonable matricesC, and shed light on the eigenvalues of the perturbed quadratic matrix polynomials. Our proofs are independent and seem to be simpler     

Extremal feedback perturbations for families of closed operators

Let T- S, be a family of not necessarily bounded semi-Fredholm operators, where T and S are operators acting between Banach spaces X and Y, and where S is bounded with D(S) D(T). For compact sets , as well as for certain open sets , we investigate existence and minimal rank of bounded feedback perturbations of the form F=BE such that min.ind (T-S+F)=0 for all . Here B is a given operator from a linear space Z to Y and E is some operator from X to Z.We give a simple characterization of that situation, when such regularizing feedback perturbations exist and show that for compact sets the minimal rank never exceeds max { min.ind (T-S) }+1. Moreover, an example shows that the minimal rank, in fact, may increase from max {...} to max {...}+1, if the given B enforces a certain structure of the feedbachk perturbation F.However, the minimal rank is equal to max { min.ind (T-S) }, if is an open set such that min.ind (T-S) already vanishes for all but finitely many points . We illustrate this result by applying it to the stabilization of certain infinite-dimensional dynamical systems in Hilbert space.     

Statistical topography. I. Fractal dimension of coastlines and number-area rule for Islands

Summary Statistical topography involves the study of the geometrical properties of the iso-sets (contour lines or surfaces) of a random potential (r). Previous work [1,2] has addressed coastlines on a random relief (x, y) that possess a single characteristic spatial scale with topography belonging to the universality class of the random percolation problem. In the present paper this previous analytical approach is extended to the case of a multiscale random function with a power spectrum of scales, ^H , in a wide range of wavelengths, ₀ < < _m . It is found that the pattern of the coastline differs significantly from that of a monoscale landscape provided that –3/4 <H < 1, with the case –3/4 <H < 0 corresponding to the long-range correlated percolation and 0 <H < 1 to the fractional Brownian relief. The expression for the fractal dimension of an individual coastline is derived,D _h = (10 – 3H)/7, the maximum value of whichD _h = 7/4, corresponds to the monoscale relief. The distribution functionF(a) of level lines over their sizea is calculated:F(a) a ^–4(1-H)/7, for ₀ a _m . A comparison of the theoretical results with geographical data is presented.