Low rank perturbation of regular matrix polynomials

Let A(λ) be a complex regular matrix polynomial of degree ? with g elementary divisors corresponding to the finite eigenvalue λ₀. We show that for most complex matrix polynomials B(λ) with degree at most ? satisfying rank the perturbed polynomial (A+B)(λ) has exactly elementary divisors corresponding to λ₀, and we determine their degrees. If does not exceed the number of λ₀-elementary divisors of A(λ) with degree greater than 1, then the λ₀-elementary divisors of (A+B)(λ) are the elementary divisors of A(λ) corresponding to λ₀ with smallest degree, together with rank(B(λ)-B(λ₀)) linear λ₀-elementary divisors. Otherwise, the degree of all the λ₀-elementary divisors of (A+B)(λ) is one. This behavior happens for any matrix polynomial B(λ) except those in a proper algebraic submanifold in the set of matrix polynomials of degree at most ?. If A(λ) has an infinite eigenvalue, the corresponding result follows from considering the zero eigenvalue of the perturbed dual polynomial.     

First order spectral perturbation theory of square singular matrix polynomials

We develop first order eigenvalue expansions of one-parametric perturbations of square singular matrix polynomials. Although the eigenvalues of a singular matrix polynomial P(λ) are not continuous functions of the entries of the coefficients of the polynomial, we show that for most perturbations they are indeed continuous. Given an eigenvalue λ₀ of P(λ) we prove that, for generic perturbations M(λ) of degree at most the degree of P(λ), the eigenvalues of P(λ)+?M(λ) admit covergent series expansions near λ₀ and we describe the first order term of these expansions in terms of M(λ₀) and certain particular bases of the left and right null spaces of P(λ₀). In the important case of λ₀ being a semisimple eigenvalue of P(λ) any bases of the left and right null spaces of P(λ₀) can be used, and the first order term of the eigenvalue expansions takes a simple form. In this situation we also obtain the limit vector of the associated eigenvector expansions.     

Deflating quadratic matrix polynomials with structure preserving transformations

Given a pair of distinct eigenvalues (λ₁,λ₂) of an n×n quadratic matrix polynomial Q(λ) with nonsingular leading coefficient and their corresponding eigenvectors, we show how to transform Q(λ) into a quadratic of the form having the same eigenvalue s as Q(λ), with Q_d(λ) an (n-1)×(n-1) quadratic matrix polynomial and q(λ) a scalar quadratic polynomial with roots λ₁ and λ₂. This block diagonalization cannot be achieved by a similarity transformation applied directly to Q(λ) unless the eigenvectors corresponding to λ₁ and λ₂ are parallel. We identify conditions under which we can construct a family of 2n×2n elementary similarity transformations that (a) are rank-two modifications of the identity matrix, (b) act on linearizations of Q(λ), (c) preserve the block structure of a large class of block symmetric linearizations of Q(λ), thereby defining new quadratic matrix polynomials Q₁(λ) that have the same eigenvalue s as Q(λ), (d) yield quadratics Q₁(λ) with the property that their eigenvectors associated with λ₁ and λ₂ are parallel and hence can subsequently be deflated by a similarity applied directly to Q₁(λ). This is the first attempt at building elementary transformations that preserve the block structure of widely used linearizations and which have a specific action.     

The eigenvalue problem for linear and affine iterated function systems

The eigenvalue problem for a linear function L centers on solving the eigen-equation . This paper generalizes the eigenvalue problem from a single linear function to an iterated function system F consisting of possibly an infinite number of linear or affine functions. The eigen-equation becomes F(X)=λX, where λ>0 is real, X is a compact set, and F(X)=?_f∈Ff(X). The main result is that an irreducible, linear iterated function system F has a unique eigenvalue λ equal to the joint spectral radius of the functions in F and a corresponding eigenset S that is centrally symmetric, star-shaped, and full dimensional. Results of Barabanov and of Dranisnikov-Konyagin-Protasov on the joint spectral radius follow as corollaries.     

Two results on homogeneous Hessian nilpotent polynomials

Let z=(z₁,…,z_n) and , the Laplace operator. A formal power series P(z) is said to be Hessian Nilpotent (HN) if its Hessian matrix is nilpotent. In recent developments in [M. de Bondt, A. van den Essen, A reduction of the Jacobian conjecture to the symmetric case, Proc. Amer. Math. Soc. 133 (8) (2005) 2201-2205. [MR2138860]; G. Meng, Legendre transform, Hessian conjecture and tree formula, Appl. Math. Lett. 19 (6) (2006) 503-510. [MR2170971]. See also math-ph/0308035; W. Zhao, Hessian nilpotent polynomials and the Jacobian conjecture, Trans. Amer. Math. Soc. 359 (2007) 249-274. [MR2247890]. See also math.CV/0409534], the Jacobian conjecture has been reduced to the following so-called vanishing conjecture (VC) of HN polynomials: for any homogeneous HN polynomialP(z) (of degreed=4), we haveΔ^mP^m+1(z)=0for anym?0. In this paper, we first show that the VC holds for any homogeneous HN polynomial P(z) provided that the projective subvarieties Z_P and Z_σ2 of CPⁿ⁻¹ determined by the principal ideals generated by P(z) and , respectively, intersect only at regular points of Z_P. Consequently, the Jacobian conjecture holds for the symmetric polynomial maps F=z−∇P with P(z) HN if F has no non-zero fixed point w∈Cⁿ with . Secondly, we show that the VC holds for a HN formal power series P(z) if and only if, for any polynomial f(z), Δ^m(f(z)P(z)^m)=0 when m?0.     

Singular Sturm-Liouville problems whose coefficients depend rationally on the eigenvalue parameter

Let −Dω(·,z)D+q be a differential operator in L²(0,∞) whose leading coefficient contains the eigenvalue parameter z. For the case that ω(·,z) has the particular form

     

Piecewise linear maps, Liapunov exponents and entropy

Let be a class of piecewise linear maps associated with a transition matrix A. In this paper, we prove that if fA,x∈LA, then the Liapunov exponent λ(x) of fA,x is equal to a measure theoretic entropy hmA,x of fA,x, where mA,x is a Markov measure associated with A and x. The Liapunov exponent and the entropy are computable by solving an eigenvalue problem and can be explicitly calculated when the transition matrix A is symmetric. Moreover, we also show that maxxλ(x)=maxxhmA,x=log(λ₁), where λ₁ is the maximal eigenvalue of A.     

Matrix refinement type equations

Taking advantage of perpetuities and the asymptotic behavior of products of random matrices we obtain the direct form of the Fourier transform of an L¹-solution of the following random matrix refinement type equation

f(x)=_∫Ω|detL(ω)|C(ω)f(L(ω)x-M(ω))P(dω),     

Polynomial dynamical systems and the Korteweg—de Vries equation

We explicitly construct polynomial vector fields L_k, k = 0, 1, 2, 3, 4, 6, on the complex linear space C⁶ with coordinates X = (x₂, x₃, x₄) and Z = (z₄, z₅, z₆). The fields L_k are linearly independent outside their discriminant variety Δ ? C⁶ and are tangent to this variety. We describe a polynomial Lie algebra of the fields L_k and the structure of the polynomial ring C[X,Z] as a graded module with two generators x₂ and z₄ over this algebra. The fields L₁ and L₃ commute. Any polynomial P(X,Z) ∈ C[X,Z] determines a hyperelliptic function P(X,Z)(u₁, u₃) of genus 2, where u₁ and u₃ are the coordinates of trajectories of the fields L₁ and L₃. The function 2x₂(u₁, u₃) is a two-zone solution of the Korteweg–de Vries hierarchy, and ?z₄(u₁, u₃)/?u₁ = ?x₂(u₁, u₃)/?u₃.     

The distance from a matrix polynomial to matrix polynomials with a prescribed multiple eigenvalue

For a matrix polynomial P(λ) and a given complex number μ, we introduce a (spectral norm) distance from P(λ) to the matrix polynomials that have μ as an eigenvalue of geometric multiplicity at least κ, and a distance from P(λ) to the matrix polynomials that have μ as a multiple eigenvalue. Then we compute the first distance and obtain bounds for the second one, constructing associated perturbations of P(λ).     

Palindromic companion forms for matrix polynomials of odd degree

The standard way to solve polynomial eigenvalue problems P(λ)x=0 is to convert the matrix polynomial P(λ) into a matrix pencil that preserves its spectral information — a process known as linearization. When P(λ) is palindromic, the eigenvalues, elementary divisors, and minimal indices of P(λ) have certain symmetries that can be lost when using the classical first and second Frobenius companion linearizations for numerical computations, since these linearizations do not preserve the palindromic structure. Recently new families of pencils have been introduced with the goal of finding linearizations that retain whatever structure the original P(λ) might possess, with particular attention to the preservation of palindromic structure. However, no general construction of palindromic linearizations valid for all palindromic polynomials has as yet been achieved. In this paper we present a family of linearizations for odd degree polynomials P(λ) which are palindromic whenever P(λ) is, and which are valid for all palindromic polynomials of odd degree. We illustrate our construction with several examples. In addition, we establish a simple way to recover the minimal indices of the polynomial from those of the linearizations in the new family.     

On polynomials orthogonal to all powers of a given polynomial on a segment

In this paper we investigate the following “polynomial moment problem”: for a given complex polynomial P(z) and distinct a,b∈C to describe polynomials q(z) orthogonal to all powers of P(z) on [a,b]. We show that for given P(z), q(z) the condition that q(z) is orthogonal to all powers of P(z) is equivalent to the condition that branches of the algebraic function Q(P⁻¹(z)), where , satisfy a certain system of linear equations over Z. On this base we provide the solution of the polynomial moment problem for wide classes of polynomials. In particular, we give the complete solution for polynomials of degree less than 10.     

On an eigenvalue problem involving variable exponent growth conditions

We study the problem in Ω, u=0 on ∂Ω, where Ω is a bounded domain in R^N, is a continuous function and λ and ε are two positive constants. We prove that for any ε>0 each λ∈(0,λ₁) is an eigenvalue of the above problem, where λ₁ is the principal eigenvalue of the Laplace operator on Ω. Moreover, for each eigenvalue λ∈(0,λ₁) it corresponds a unique eigenfunction. The proofs will be based on the Banach fixed point theorem combined with adequate variational techniques.     

Modular metric spaces, I: Basic concepts

The notion of a modular is introduced as follows. A (metric) modular on a set X is a function w:(0,∞)×X×X→[0,∞] satisfying, for all x,y,z∈X, the following three properties: x=y if and only if w(λ,x,y)=0 for all λ>0; w(λ,x,y)=w(λ,y,x) for all λ>0; w(λ+μ,x,y)≤w(λ,x,z)+w(μ,y,z) for all λ,μ>0. We show that, given x₀∈X, the set X_w={x∈X:lim_λ→∞w(λ,x,x₀)=0} is a metric space with metric , called a modular space. The modular w is said to be convex if (λ,x,y)?λw(λ,x,y) is also a modular on X. In this case X_w coincides with the set of all x∈X such that w(λ,x,x₀)<∞ for some λ=λ(x)>0 and is metrizable by . Moreover, if or , then ; otherwise, the reverse inequalities hold. We develop the theory of metric spaces, generated by modulars, and extend the results by H. Nakano, J. Musielak, W. Orlicz, Ph. Turpin and others for modulars on linear spaces.     

Exponential stability of wave equations with potential and indefinite damping

First, we consider the linear wave equation utt−uxx+a(x)ut+b(x)u=0 on a bounded interval (0,L)⊂R. The damping function a is allowed to change its sign. If is positive and the spectrum of the operator (_∂xx−b) is negative, exponential stability is proved for small . Explicit estimates of the decay rate ω are given in terms of and the largest eigenvalue of (_∂xx−b). Second, we show the existence of a global, small, smooth solution of the corresponding nonlinear wave equation utt−σx(ux)+a(x)ut+b(x)u=0, if, additionally, the negative part of a is small enough compared with ω.     

Exposants de Lyapounov pour opérateurs de Schrödinger discrètes quasi-périodiquesLyapounov exponents for discrete,quasi-periodic Schrödinger operators

We consider quasi-periodic Schrödinger operators H on $\text{Z}$ of the form H=H_λ,x,ω=λv(x+nω)δ_n,n′+Δ where v is a non-constant real analytic function on the d-torus $\text{T}{}^{\mathrm{d}}\text{(d?1)}$ and Δ denotes the discrete lattice Laplacian on $\text{Z}$. Denote by L_ω(E) the Lyapounov exponent, considered as function of the energy E and the rotation vector $\text{ω∈}\text{T}{}^{\mathrm{d}}$. It is shown that for |λ|>λ₀(v), there is the uniform minoration $\text{L}{}_{\mathrm{\omega }}\text{(E)>}\text{1}\text{2}\text{log}\text{|λ|}$ for all E and ω. For all λ and ω, L_ω(E) is a continuous function of E. Moreover, L_ω(E) is jointly continuous in (ω,E), at any point $\text{(ω}{}_0\text{,E}{}_0\text{)∈}\text{T}{}^{\mathrm{d}}\text{×}\text{R}$ such that k·ω₀≠0 for all $\text{k∈}\text{Z}{}^{\mathrm{d}}\text{?{0}}$. To cite this article: J. Bourgain, C. R. Acad. Sci. Paris, Ser. I 335 (2002) 529–531.     

Calculation of the eigenvalues of Schrödinger equations by an extension of Hill's method

The eigenfunctions of the one dimensional Schrödinger equation Ψ″ + [E ? V(x)]Ψ=0, where V(x) is a polynomial, are represented by expansions of the form $\text{∑}\text{k}\text{=0}\text{∞}{}^{\mathrm{<mtext>c</mtext>}}\text{k}{}^{\mathrm{?}}\text{k}\text{(ω,}\text{x}\text{)}$. The functions ?_k (ω, x) are chosen in such a way that recurrence relations hold for the coefficients c_k: examples treated are D_k(ωx) (Weber-Hermite functions), exp (?ωx²)x^k, exp (?cx^q)D_k(ωx). From these recurrence relations, one considers an infinite bandmatrix whose finite square sections permit to solve approximately the original eigenproblem. It is then shown how a good choice of the parameter ω may reduce dramatically the complexity of the computations, by a theoretical study of the relation holding between the error on an eigenvalue, the order of the matrix, and the value of ω. The paper contains tables with 10 significant figures of the 30 first eigenvalues corresponding to V(x) = x^2m, m = 2(1)7, and the 6 first eigenvalues corresponding to V(x) = x² + λx¹⁰ and x² + λx¹², λ = .01(.01).1(.1)1(1)10(10)100.     

ω-hypoelliptic differential operators of constant strength

We study ω-hypoelliptic differential operators of constant strength. We show that any operator with constant strength and coefficients in which is homogeneous ω-hypoelliptic is also σ-hypoelliptic for any weight function σ=O(ω). We also present a sufficient condition in order to ensure that a differential operator admits a parametrix and, as a consequence, we obtain some conditions on the weights (ω,σ) to conclude that, for any operator P(x,D) with constant strength, the σ-hypoellipticity of the frozen operator P(x₀,D) implies the ω-hypoellipticity of P(x,D). This requires the use of pseudodifferential operators.     

The transcendental meromorphic solutions of a certain type of nonlinear differential equations

In this paper, we study the differential equations of the following form w²+R(z)²(w(k))=Q(z), where R(z), Q(z) are nonzero rational functions. We proved the following three conclusions: (1) If either P(z) or Q(z) is a nonconstant polynomial or k is an even integer, then the differential equation w²+P²(z)²(w(k))=Q(z) has no transcendental meromorphic solution; if P(z), Q(z) are constants and k is an odd integer, then the differential equation has only transcendental meromorphic solutions of the form f(z)=acos(bz+c). (2) If either P(z) or Q(z) is a nonconstant polynomial or k>1, then the differential equation w²+(z−z₀)P²(z)²(w(k))=Q(z) has no transcendental meromorphic solution, furthermore the differential equation w²+A(z−z₀)²(w^′)=B, where A, B are nonzero constants, has only transcendental meromorphic solutions of the form , where a, b are constants such that Ab²=1, a²=B. (3) If the differential equation , where P is a nonconstant polynomial and Q is a nonzero rational function, has a transcendental meromorphic solution, then k is an odd integer and Q is a polynomial. Furthermore, if k=1, then Q(z)≡C (constant) and the solution is of the form f(z)=Bcosq(z), where B is a constant such that B²=C and q^′(z)=±P(z).     

The inverse of a matrix polynomial

An explicit representation is obtained for P(z)^?1 when P(z) is a complex n×n matrix polynomial in z whose coefficient of the highest power of z is the identity matrix. The representation is a sum of terms involving negative powers of z?λ for each λ such that P(λ) is singular. The coefficients of these terms are generated by sequences u_k, v_k of 1×n and n×1 vectors, respectively, which satisfy u₁≠0, v₁≠0, ∑^k?1_h=0(1?h!)u_k?hP^(h)(λ)=0, ∑^k?1_h=0(1?h!)P^(h)(λ)v_k?h=0, and certain orthogonality relations. In more general cases, including that when P(z) is analytic at λ but not necessarily a polynomial, the terms in the representation involving negative powers of z?λ provide the principal part of the Laurent expansion for P(z)^?1 in a punctured neighborhood of z=λ.