Total positivity of Narayana matrices

We prove the total positivity of the Narayana triangles of type $A$ and type $B$, and thus affirmatively confirm a conjecture of Chen, Liang and Wang and a conjecture of Pan and Zeng. We also prove the strict total positivity of the Narayana squares of type $A$ and type $B$.     

Total positivity and accurate computations with Gram matrices of Said-Ball bases

In this article, it is proved that Gram matrices of totally positive bases of the space of polynomials of a given degree on a compact interval are totally positive. Conditions to guarantee computations to high relative accuracy with those matrices are also obtained. Furthermore, a fast and accurate algorithm to compute the bidiagonal factorization of Gram matrices of the Said-Ball bases is obtained and used to compute to high relative accuracy their singular values and inverses, as well as the solution of some linear systems associated with these matrices. Numerical examples are included.     

The smallest eigenvalue of large Hankel matrices generated by a deformed Laguerre weight

We study the asymptotic behavior of the smallest eigenvalue, λ_N, of the Hankel (or moments) matrix denoted by , with respect to the weight . An asymptotic expression of the polynomials orthogonal with w(x) is established. Using this, we obtain the specific asymptotic formulas of λ_N in this paper. Applying a parallel numerical algorithm, we get a variety of numerical results of λ_N corresponding to our theoretical calculations.     

An inequality for non-negative matrices

Let A be an n × n matrix with non-negative entries and no entry in (0, 1). We prove that there exist integers r, s with 0 r s 2ⁿ such that A^r A^s. We prove that 2ⁿ cannot be replaced with e√n log n. We also give an application to the theory of formal languages.     

Completely normal matrices

Normal matrices in which all submatrices are normal are said to be completely normal. We characterize this class of matrices, determine the possible inertias of a particular completely normal matrix, and show that real matrices in this class are closed under (general) Schur complementation. We provide explicit formulas for the Moore–Penrose inverse of a completely normal matrix of size at least four. A result on irreducible principally normal matrices is derived as well.     

Chains of dog-ears for completely positive matrices

A necessary and sufficient condition to determine the complete positivity of a matrixwith a particular graph, in dependence of complete positivity of smaller matrices, is given. Under some singularity assumptions, this condition furnishes a characterization for completely positive matrices with a “non-crossing cycle” as associated graph. In particular the characterization holds for singular pentadiagonal matrices.     

Inversion of conjugate-Toeplitz matrices and conjugate-Hankel matrices

The inverses of conjugate-Toeplitz (CT) and conjugate-Hankel (CH) matrices can be expressed by the Gohberg–Heinig type formula. We obtain an explicit inverse formula of CT matrix. Similarly, the formula and the decomposition of the inverse of a CH matrix are provided. Also the stability of the inverse formulas of CT and CH matrices are discussed. Examples are provided to verify the feasibility of the algorithms.     

Approximation of high-dimensional kernel matrices by multilevel circulant matrices

Kernels are important in developing a variety of numerical methods, such as approximation, interpolation, neural networks, machine learning and meshless methods for solving engineering problems. A common problem of these kernel-based methods is to calculate inverses of kernel matrices generated by a kernel function and a set of points. Due to the denseness of these matrices, finding their inverses is computationally costly. To overcome this difficulty, we introduce in this paper an approximation of the kernel matrices by appropriate multilevel circulant matrices so that the fast Fourier transform can be applied to reduce the computational cost. Convergence analysis for the proposed approximation is established based on certain decay properties of the kernels.     

On the solutions of the equation AXB = C under Toeplitz‐like and Hankel matrices constraint

In this paper, we are mainly concerned with 2 types of constrained matrix equation problems of the form AXB=C, the least squares problem and the optimal approximation problem, and we consider several constraint matrices, such as general Toeplitz matrices, upper triangular Toeplitz matrices, lower triangular Toeplitz matrices, symmetric Toeplitz matrices, and Hankel matrices. In the first problem, owing to the special structure of the constraint matrix , we construct special algorithms; necessary and sufficient conditions are obtained about the existence and uniqueness for the solutions. In the second problem, we use von Neumann alternating projection algorithm to obtain the solutions of problem. Then we give 2 numerical examples to demonstrate the effectiveness of the algorithms.     

The question of equivalence for generalized Hausdorff matrices

In the present paper we investigate when Hausdorff matrices and generalized Hausdorff matrices, with the same mass function, are equivalent, as bounded operators on c and ?p.     

Accurate computations of matrices with bidiagonal decomposition using methods for totally positive matrices

A class of sign‐symmetric P‐matrices including all nonsingular totally positive matrices and their inverses as well as tridiagonal nonsingular H‐matrices is presented and analyzed. These matrices present a bidiagonal decomposition that can be used to obtain algorithms to compute with high relative accuracy their singular values, eigenvalues, inverses, or their LDU factorization. Copyright © 2012 John Wiley & Sons, Ltd.     

A Sufficient condition for strict total positivity of a matrix

We establish a sufficient condition for strict total positivity of a matrix In particular, we show that if the (positive) elements of a square matrix grow sufficiently fast as their distance from the diagonal of the matrix increases, then the matrix is strictly totally positive.     

Williamson matrices up to order 59

A recent result of Schmidt has brought Williamson matrices back into the spotlight. In this article, a new algorithm is introduced to search for hard to find Williamson matrices. We find all nonequivalent Williamson matrices of odd order n up to n = 59. It turns out that there are none for n = 35, 47, 53, 59 and it seems that the Turyn class may be the only infinite class of these matrices.      

From random matrices to quasi-periodic Jacobi matrices via orthogonal polynomials

We present an informal review of results on asymptotics of orthogonal polynomials, stressing their spectral aspects and similarity in two cases considered. They are polynomials orthonormal on a finite union of disjoint intervals with respect to the Szegö weight and polynomials orthonormal on with respect to varying weights and having the same union of intervals as the set of oscillations of asymptotics. In both cases we construct double infinite Jacobi matrices with generically quasi-periodic coefficients and show that each of them is an isospectral deformation of another. Related results on asymptotic eigenvalue distribution of a class of random matrices of large size are also shortly discussed.     

Toeplitz approximate inverse preconditioner for banded Toeplitz matrices

In this paper we introduce a new preconditioner for banded Toeplitz matrices, whose inverse is itself a Toeplitz matrix. Given a banded Hermitian positive definite Toeplitz matrixT, we construct a Toepliz matrixM such that the spectrum ofMT is clustered around one; specifically, if the bandwidth ofT is , all but eigenvalues ofMT are exactly one. Thus the preconditioned conjugate gradient method converges in +1 steps which is about half the iterations as required by other preconditioners for Toepliz systems that have been suggested in the literature. This idea has a natural extension to non-banded and non-Hermitian Toeplitz matrices, and to block Toeplitz matrices with Toeplitz blocks which arise in many two dimensional applications in signal processing. Convergence results are given for each scheme, as well as numerical experiments illustrating the good convergence properties of the new preconditioner.Partly supported by a travel fund from the Deutsche Forschungsgemeinschaft.Research supported in part by Oak Ridge Associated Universities grant no. 009707.     

Out-of-band quasiseparable matrices

We study structured matrices which consist of a band part and quasiseparable parts below and upper the band. We extend algorithms known for quasiseparable matrices, i.e. for the case when the band consists of the main diagonal only, to a wider class of matrices. The matrices which we consider may be treated as an usual quasiseparable matrices with larger orders of generators. Hence one can apply the methods developed for usual quasiseparable matrices and obtain various linear complexity O(N) algorithms. However in this case the coefficients in N in the complexity estimates turns out to be quite large. In this paper we use the structure more accurately by division of the matrix into three parts in which the middle part is the band instead of diagonal as it is used for usual quasiseparable matrices. This approach allows to use better the structure of the matrix in order to improve the coefficients in N in the complexity estimates for the algorithms. This method works for algorithms which keep invariant the structure.     

Properties of balanced and perfect matrices

In this paper we define wheel matrices and characterize some properties of matrices that are perfect but not balanced. A consequence of our results is that if a matrixA is perfect then we can construct a polynomial number of submatricesA _I,,A _n ofA such thatA is balanced if and only if all the submatricesA ₁,,A _n ofA are perfect. This implies that if the problem of testing perfection is polynomially solvable, then so is the problem of testing balancedness.Partial support under NSF Grants DMS-8606188, ECS-8800281 and DDM-8800281.     

Set inertia-preserving matrices

The inertia-preservers of several sets of matrices are identified. The sets include: all real matrices, all complex matrices, triangular matrices, real symmetric matrices and Hermitian matrices.     

A spectral equivalence for Jacobi matrices

We use the classical results of Baxter and Golinskii–Ibragimov to prove a new spectral equivalence for Jacobi matrices on . In particular, we consider the class of Jacobi matrices with conditionally summable parameter sequences and find necessary and sufficient conditions on the spectral measure such that and lie in or for s1.