Two linear transformations each tridiagonal with respect to an eigenbasis of the other

Let denote a field, and let V denote a vector space over with finite positive dimension. We consider a pair of linear transformations A:V→V and A^*:V→V satisfying both conditions below:

1. [(i)] There exists a basis for V with respect to which the matrix representing A is diagonal and the matrix representing A^* is irreducible tridiagonal.

2. [(ii)] There exists a basis for V with respect to which the matrix representing A^* is diagonal and the matrix representing A is irreducible tridiagonal.

We call such a pair a Leonard pair on V. Refining this notion a bit, we introduce the concept of a Leonard system. We give a complete classification of Leonard systems. Integral to our proof is the following result. We show that for any Leonard pair A,A^* on V, there exists a sequence of scalars β,γ,γ^*,,^* taken from such that both

where [r,s] means rs−sr. The sequence is uniquely determined by the Leonard pair if the dimension of V is at least 4. We conclude by showing how Leonard systems correspond to q-Racah and related polynomials from the Askey scheme.     

Hypernormal matrices

An nxn matrix A is hypernormal if APA^*=A^*PA for all permutation matrices P. We shall explain how to construct hypernormal matrices.     

Complementarity forms of theorems of Lyapunov and Stein, and related results

The well-known Lyapunov's theorem in matrix theory / continuous dynamical systems asserts that a (complex) square matrix A is positive stable (i.e., all eigenvalues lie in the open right-half plane) if and only if there exists a positive definite matrix X such that AX+XA^* is positive definite. In this paper, we prove a complementarity form of this theorem: A is positive stable if and only if for any Hermitian matrix Q, there exists a positive semidefinite matrix X such that AX+XA^*+Q is positive semidefinite and X[AX+XA^*+Q]=0. By considering cone complementarity problems corresponding to linear transformations of the form I−S, we show that a (complex) matrix A has all eigenvalues in the open unit disk of the complex plane if and only if for every Hermitian matrix Q, there exists a positive semidefinite matrix X such that X−AXA^*+Q is positive semidefinite and X[X−AXA^*+Q]=0. By specializing Q (to −I), we deduce the well known Stein's theorem in discrete linear dynamical systems: A has all eigenvalues in the open unit disk if and only if there exists a positive definite matrix X such that X−AXA^* is positive definite.     

The construction of irreducible normalized representations of the general linear group

A method is described for constructing in an explicit form an irreducible representation T of M_n(F), the set of all n × n matrices over the real or complex field F, satisfying the condition T(A^*)=T^*(A) for all A∈M_n(F).     

Weak stability for matrices

An n×n complex matrix A is called weak stable if there exists a matrix W such that W+W^* is positive definite and such that AW+W^*A^* is positive definite. In this note several characterizations for weak stability of a matrix are given, and conditions (on A) allowing W to be a diagonal matrix are also considered. A consequence of our results here is a characterization for nonsingular M-matrices.     

On real matrices with positive definite symmetric component

Inequalities concerning real square matrices A with positive definite symmetric component A+A^*are derived from certain inertia relations which hold for any complex (not necessarily real) square matrices A with positive definite

A+A^*     

A theory of electromagnetic and gravitational fields

Careful calculations using classical field theory show that if a macroscopic ball with uniform surface charge (say, a billiard ball with 1E6 excess electrons) is released near the surface of the earth, it will almost instantaneously accelerate to relativistic speed and blow a hole in the ground. This absurd prediction is just the macroscopic version of the self-force problem for charged particles [1]. Furthermore, if one attempts to develop from electromagnetism a parallel theory for gravitational [2], the result is the same, self-acceleration.

The basis of the new theory is a measure of energy density for any wave equation [3–5]. Given any solution of any four-vector wave equation in spacetime (for example, the potentials (c^-1φA)=(A⁰,A¹,A²,A³) in electromagnetism), one can form the 16th first order partial derivatives of the vector components, with respect to the time and space variables (ct,x) = (x⁰, x¹, x², x³). The sum of the squares of the 16 terms is a natural energy function [6, p. 283] (satisfying a conservation law . Such energy functions are routinely utilized by mathematicians as Lyapunov functions in the theory of stability of waves with boundary conditions. A Lagrangian using this sum leads to a new energy tensor for electromagnetic and gravitational fields, an alternative to that in [7].     

Doubly stochastic matrices whose powers eventually stop

In this note we characterize doubly stochastic matrices A whose powers A,A²,A³,… eventually stop, i.e., A^p=A^p+1= for some positive integer p. The characterization enables us to determine the set of all such matrices.     

Numerical approximation of the product of the square root of a matrix with a vector

Given an n×n symmetric positive definite matrix A and a vector , two numerical methods for approximating are developed, analyzed, and computationally tested. The first method applies a Newton iteration to a specific nonlinear system to approximate while the second method applies a step-control method to numerically solve a specific initial-value problem to approximate . Assuming that A is first reduced to tridiagonal form, the first method requires O(n²) operations per iteration while the second method requires O(n) operations per iteration. In contrast, numerical methods that first approximate A^1/2 and then compute generally require O(n³) operations per iteration.     

On the fourth power moment of Fourier coefficients of cusp form

Let a(n)be the Fourier coefficients of a holomorphic cusp form of weightκ=2n≥12 for the full modular group and A(x)=∑_(n≤x)a(n).In this paper,we establish an asymptotic formula of the fourth power moment of A(x)and prove that ∫T1A~4(x)dx=3/(64κπ~4)s_4;2()T~(2κ)+O(T~(2κ-δ_4+ε))with δ_4=1/8,which improves the previous result.     

Products of complic cosquares and pseudo-involutory matrices

LetF be a field with (nontrivial) involution (i.e.F-conjugation). A nonsingular matrix Aover Fis called a complic F-cosquare provided A=S^*-1for some matrix Sover Fand is called p.i. (pseudo-involutory) provided A=A^-1 It is shown that Ais a complic F-cosquare iff Ais the product of two p.i. matrices over Fand that det (AA)=1 iff Ais the product of two complic F-cosquares (hence iff A is the product of four p.i. matrices over F). It is conjectured that, except for one obvious case (2 x 2 matrices over the field of order 2), every unimodular matrix A over an arbitrary field Fis a product S¹ST:1T with S1 and Tover FThis conjecture is proved for matricesAof order ≤3.     

Some closed form formulae for the intersection of two special matrices under the star order

In this note we present closed form formulae for the intersection A ∧ B under the star order, for certain pairs of special matrices. In particular, it is shown that for two m×n partial isometries P and QP∧Q=P[I-(I-P^*Q)^†(I-P^*Q)]. Some further possible representations for P∧Q are considered.     

On the inertia of the Lyapunov transform AH+HA* for H >0

Given an arbitrary n×n matrix A with complex entries, we characterize all inertia triples (abc) that are attained by the Lyapunov transform AH+ HA^*, as H varies over the set of all n× n positive definite matrices.     

On the normality of operator products

If A and B^* are hyponormal operators such that AB is compact, then the normality of AB implics that of BA. This is a generalization of an old result of Wiegmann and Kaplansky on products of normal operators.     

A characterization and representation of the generalized inverse AT,S and its applications

This paper presents an explicit expression for the generalized inverse A⁽²⁾_T,S. Based on this, we established the characterization, the representation theorem and the limiting process for A⁽²⁾_T,S. As an application, we estimate the error bound of the iterative method for approximating A⁽²⁾_T,S.     

Numerical method for solving 3D inverse scattering problems

Let q(x) L²(D), D R³ is a bounded domain, q = 0 outside D, q is real-valued. Assume that A(\Gj;\t';,\Gj;,k) A(\Gj;\t';,\Gj), the scattering amplitude, is known for all \Gj;|t',\Gj; S², S² is the unit sphere, an d a fixed k \r>0. These data determine q(x) uniquely and a numerical method is given for computing q(x).     

The growth of laguerre matrix polynomials on bounded intervals

Let A be a matrix in ^r×r such that Re(z) > −1/2 for all the eigenvalues of A and let {π_n^(A,1/2) (x)} be the normalized sequence of Laguerre matrix polynomials associated with A. In this paper, it is proved that π_n^(A,1/2) (x) = O(n^(A)/2ln^r−1(n)) and π_n+1^(A,1/2) (x) − π_n^(A,1/2) (x) = O(n^((A)−1)/2ln^r−1(n)) uniformly on bounded intervals, where (A) = max{Re(z); z eigenvalue of A}.     

Minimal enclosing parallelepiped in 3D

We investigate the problem of finding a minimal volume parallelepiped enclosing a given set of n three-dimensional points. We give two mathematical properties of these parallelepipeds, from which we derive two algorithms of theoretical complexity O(n⁶). Experiments show that in practice our quickest algorithm runs in O(n²) (at least for n10⁵). We also present our application in structural biology.     

A note on k-generalized projections

In this note, we investigate characterizations for k-generalized projections (i.e., A^k = A^*) on Hilbert spaces. The obtained results generalize those for generalized projections on Hilbert spaces in [Hong-Ke Du, Yuan Li, The spectral characterization of generalized projections, Linear Algebra Appl. 400 (2005) 313–318] and those for matrices in [J. Benítez, N. Thome, Characterizations and linear combinations of k-generalized projectors, Linear Algebra Appl. 410 (2005) 150–159].