Singularity, Wielandt’s lemma and singular values

In this study, some upper and lower bounds for singular values of a general complex matrix are investigated, according to singularity and Wielandt’s lemma of matrices. Especially, some relationships between the singular values of the matrix A and its block norm matrix are established. Based on these relationships, one may obtain the effective estimates for the singular values of large matrices by using the lower dimension norm matrices. In addition, a small error in Piazza (2002) [G. Piazza, T. Politi, An upper bound for the condition number of a matrix in spectral norm, J. Comput. Appl. Math. 143 (1) (2002) 141-144] is also corrected. Some numerical experiments on saddle point problems show that these results are simple and sharp under suitable conditions.     

Stopping criteria for the Ando-Li-Mathias and Bini-Meini-Poloni geometric means

We provide an upper bound for the number of iterations necessary to achieve a desired level of accuracy for the Ando-Li-Mathias [Linear Algebra Appl. 385 (2004) 305-334] and Bini-Meini-Poloni [Math. Comput. 79 (2010) 437-452] symmetrization procedures for computing the geometric mean of n positive definite matrices, where accuracy is measured by the spectral norm and the Thompson metric on the convex cone of positive definite matrices. It is shown that the upper bound for the number of iterations depends only on the diameter of the set of n matrices and the desired convergence tolerance. A striking result is that the upper bound decreases as n increases on any bounded region of positive definite matrices.     

Bounds for norms of the matrix inverse and the smallest singular value

In the first part, we obtain two easily calculable lower bounds for ‖A^-1‖, where ‖·‖ is an arbitrary matrix norm, in the case when A is an M-matrix, using first row sums and then column sums. Using those results, we obtain the characterization of M-matrices whose inverses are stochastic matrices. With different approach, we give another easily calculable lower bounds for ‖A^-1‖_∞ and ‖A^-1‖₁ in the case when A is an M-matrix. In the second part, using the results from the first part, we obtain our main result, an easily calculable upper bound for ‖A^-1‖₁ in the case when A is an SDD matrix, thus improving the known bound. All mentioned norm bounds can be used for bounding the smallest singular value of a matrix.     

The optimal perturbation bounds of the Moore-Penrose inverse under the Frobenius norm

We obtain the optimal perturbation bounds of the Moore-Penrose inverse under the Frobenius norm by using Singular Value Decomposition, which improved the results in the earlier paper [P.-Å. Wedin, Perturbation theory for pseudo-inverses, BIT 13 (1973) 217-232]. In addition, a perturbation bound of the Moore-Penrose inverse under the Frobenius norm in the case of the multiplicative perturbation model is also given.     

Lie algebras of infinitesimal norm isometries

Given a norm on a finite dimensional vector space V, we may consider the group of all linear automorphisms which preserve it. The Lie algebra of this group is a Lie subalgebra of the endomorphism algebra of V having two properties: (1) it is the Lie algebra of a compact subgroup, and (2) it is “saturated” in a sence made precise below. We show that any Lie subalgebra satisfying these conditions is the Lie algebra of the group of linear automorphisms preserving some norm. There is an appendix on elementary Lie group theory.     

On the spectral radius and the spectral norm of Hadamard products of nonnegative matrices

We prove the spectral radius inequality ρ(A₁°A₂°?°A_k)?ρ(A₁A₂?A_k) for nonnegative matrices using the ideas of Horn and Zhang. We obtain the inequality ‖A°B‖?ρ(A^TB) for nonnegative matrices, which improves Schur’s classical inequality ‖A°B‖?‖A‖‖B‖, where ‖·‖ denotes the spectral norm. We also give counterexamples to two conjectures about the Hadamard product.     

On monotonicity of the Lanczos approximation to the matrix exponential

We prove strictly monotonic error decrease in the Euclidian norm of the Krylov subspace approximation of exp(A)φ, where φ and A are respectively a vector and a symmetric matrix. In addition, we show that the norm of the approximate solution grows strictly monotonically with the subspace dimension.     

Path of quasi-means as a geodesic

As a generalization of the Hiai-Petz geometries, we discuss two types of them where the geodesics are the quasi-arithmetic means and the quasi-geometric means respectively. Each derivative of such a geodesic might determine a new relative operator entropy. Also in these cases, the Finsler metric can be induced by each unitarily invariant norm. If the norm is strictly convex, then the geodesic is the shortest. We also give examples of the shortest paths which are not the geodesics when the Finsler metrics are induced by the Ky Fan k-norms.     

Criteria for the stability of the finiteness property and for the uniqueness of Barabanov norms

A set of matrices is said to have the finiteness property if the maximal rate of exponential growth of long products of matrices drawn from that set is realised by a periodic product. The extent to which the finiteness property is prevalent among finite sets of matrices is the subject of ongoing research. In this article, we give a condition on a finite irreducible set of matrices which guarantees that the finiteness property holds not only for that set, but also for all sufficiently nearby sets of equal cardinality. We also prove a theorem giving conditions under which the Barabanov norm associated to a finite irreducible set of matrices is unique up to multiplication by a scalar, and show that in certain cases these conditions are also persistent under small perturbations.     

Procrustes problems and associated approximation problems for matrices with k-involutory symmetries

We say that a matrix R∈C^n×n is k-involutary if its minimal polynomial is x^k-1 for some k?2, so R^k-1=R^-1 and the eigenvalues of R are 1,ζ,ζ²,…,ζ^k-1, where ζ=e^2πi/k. Let α,μ∈{0,1,…,k-1}. If R∈C^m×m, A∈C^m×n, S∈C^n×n and R and S are k-involutory, we say that A is (R,S,μ)-symmetric if RAS^-1=ζ^μA, and A is (R,S,α,μ)-symmetric if RAS^-α=ζ^μA.Let L be the class of m×n(R,S,μ)-symmetric matrices or the class of m×n(R,S,α,μ)-symmetric matrices. Given X∈C^n×t and B∈C^m×t, we characterize the matrices A in L that minimize ‖AX-B‖ (Frobenius norm), and, given an arbitrary W∈C^m×n, we find the unique matrix A∈L that minimizes both ‖AX-B‖ and ‖A-W‖. We also obtain necessary and sufficient conditions for existence of A∈L such that AX=B, and, assuming that the conditions are satisfied, characterize the set of all such A.