Three Bounds on the Minimal Singular Value: A Comparison

Similar documents
ETNA Kent State University

In 1912, G. Frobenius [6] provided the following upper bound for the spectral radius of nonnegative matrices.

Generalized Numerical Radius Inequalities for Operator Matrices

A note on polar decomposition based Geršgorin-type sets

The Eigenvalue Problem: Perturbation Theory

Computational math: Assignment 1

H-matrix theory vs. eigenvalue localization

CSE 291. Assignment Spectral clustering versus k-means. Out: Wed May 23 Due: Wed Jun 13

AN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES

Some bounds for the spectral radius of the Hadamard product of matrices

THE PERTURBATION BOUND FOR THE SPECTRAL RADIUS OF A NON-NEGATIVE TENSOR

be a Householder matrix. Then prove the followings H = I 2 uut Hu = (I 2 uu u T u )u = u 2 uut u

Lecture 10 - Eigenvalues problem

AN EXTREMUM PROPERTY OF SUMS OF EIGENVALUES' HELMUT WIELANDT

AMS526: Numerical Analysis I (Numerical Linear Algebra for Computational and Data Sciences)

Distribution for the Standard Eigenvalues of Quaternion Matrices

BOUNDS OF MODULUS OF EIGENVALUES BASED ON STEIN EQUATION

First, we review some important facts on the location of eigenvalues of matrices.

2. Linear algebra. matrices and vectors. linear equations. range and nullspace of matrices. function of vectors, gradient and Hessian

Linear Algebra: Characteristic Value Problem

Knowledge Discovery and Data Mining 1 (VO) ( )

1 Number Systems and Errors 1

A note on estimates for the spectral radius of a nonnegative matrix

The residual again. The residual is our method of judging how good a potential solution x! of a system A x = b actually is. We compute. r = b - A x!

Jordan Journal of Mathematics and Statistics (JJMS) 5(3), 2012, pp A NEW ITERATIVE METHOD FOR SOLVING LINEAR SYSTEMS OF EQUATIONS

Review of Some Concepts from Linear Algebra: Part 2

Math Spring 2011 Final Exam

Gershgorin s Circle Theorem for Estimating the Eigenvalues of a Matrix with Known Error Bounds

MAT 1332: CALCULUS FOR LIFE SCIENCES. Contents. 1. Review: Linear Algebra II Vectors and matrices Definition. 1.2.

6.241 Dynamic Systems and Control

Notes on Linear Algebra and Matrix Theory

Characterization of half-radial matrices

Bounds for Levinger s function of nonnegative almost skew-symmetric matrices

A NOTE ON CHECKING REGULARITY OF INTERVAL MATRICES

Trace Inequalities for a Block Hadamard Product

Jim Lambers MAT 610 Summer Session Lecture 2 Notes

UCSD ECE269 Handout #18 Prof. Young-Han Kim Monday, March 19, Final Examination (Total: 130 points)

Graphs and matrices with maximal energy

Introduction to Applied Linear Algebra with MATLAB

The eigenvalues of an n x n complex matrix A are the roots of the n th degree polynomial det(a-λi) = 0. Further the positive square roots of the

Recall the convention that, for us, all vectors are column vectors.

Chapter 3 Transformations

An estimation of the spectral radius of a product of block matrices

MAT Linear Algebra Collection of sample exams

OPTIMAL SCALING FOR P -NORMS AND COMPONENTWISE DISTANCE TO SINGULARITY

S.F. Xu (Department of Mathematics, Peking University, Beijing)

Preliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012

Linear Algebra Massoud Malek

AN ITERATION. In part as motivation, we consider an iteration method for solving a system of linear equations which has the form x Ax = b

Linear Algebra Practice Problems

R, 1 i 1,i 2,...,i m n.

Matrix Algebra, part 2

Inequalities For Singular Values And Traces Of Quaternion Hermitian Matrices

Z-Pencils. November 20, Abstract

AMS526: Numerical Analysis I (Numerical Linear Algebra)

Eigenvalue and Eigenvector Problems

We first repeat some well known facts about condition numbers for normwise and componentwise perturbations. Consider the matrix

Interpolating the arithmetic geometric mean inequality and its operator version

Section 3.9. Matrix Norm

Math 443 Differential Geometry Spring Handout 3: Bilinear and Quadratic Forms This handout should be read just before Chapter 4 of the textbook.

FIXED POINT ITERATIONS

Numerical Linear Algebra Homework Assignment - Week 2

Linear Algebra. Shan-Hung Wu. Department of Computer Science, National Tsing Hua University, Taiwan. Large-Scale ML, Fall 2016

OPERATIONS on large matrices are a cornerstone of

Let x be an approximate solution for Ax = b, e.g., obtained by Gaussian elimination. Let x denote the exact solution. Call. r := b A x.

Here is an example of a block diagonal matrix with Jordan Blocks on the diagonal: J

Short proofs of theorems of Mirsky and Horn on diagonals and eigenvalues of matrices

A Smoothing Newton Method for Solving Absolute Value Equations

Lecture 8: Linear Algebra Background

Exponential bases on rectangles of R d

Math 489AB Exercises for Chapter 1 Fall Section 1.0

In particular, if A is a square matrix and λ is one of its eigenvalues, then we can find a non-zero column vector X with

Comparison results between Jacobi and other iterative methods

1. Linear systems of equations. Chapters 7-8: Linear Algebra. Solution(s) of a linear system of equations (continued)

Numerical Analysis: Solving Systems of Linear Equations

Perron Frobenius Theory

The Miracles of Tropical Spectral Theory

NORMS ON SPACE OF MATRICES

(Mathematical Operations with Arrays) Applied Linear Algebra in Geoscience Using MATLAB

Matrix Vector Products

Compound matrices and some classical inequalities

arxiv: v3 [math.ra] 10 Jun 2016

MATH36001 Perron Frobenius Theory 2015

UNIFYING LEAST SQUARES, TOTAL LEAST SQUARES AND DATA LEAST SQUARES

On Powers of General Tridiagonal Matrices

Matrix Inequalities by Means of Block Matrices 1

11 a 12 a 21 a 11 a 22 a 12 a 21. (C.11) A = The determinant of a product of two matrices is given by AB = A B 1 1 = (C.13) and similarly.

Exponentials of Symmetric Matrices through Tridiagonal Reductions

Computational Methods. Eigenvalues and Singular Values

On the Definition of Two Natural Classes of Scalar Product

Announcements Wednesday, November 7

Throughout these notes we assume V, W are finite dimensional inner product spaces over C.

arxiv: v1 [math.ra] 11 Aug 2014

NOTE ON THE SKEW ENERGY OF ORIENTED GRAPHS. Communicated by Ivan Gutman. 1. Introduction

Permutation transformations of tensors with an application

Conceptual Questions for Review

Spectral Properties of Matrix Polynomials in the Max Algebra

MATH.2720 Introduction to Programming with MATLAB Vector and Matrix Algebra

HW2 - Due 01/30. Each answer must be mathematically justified. Don t forget your name.

Example: Filter output power. maximization. Definition. Eigenvalues, eigenvectors and similarity. Example: Stability of linear systems.

Transcription:

Applied Mathematical Sciences, Vol. 4, 010, no. 64, 3195-3199 Three Bounds on the Minimal Singular Value: A Comparison Kateřina Hlaváčková-Schindler Commission for Scientific Visualization Austrian Academy of Sciences Donau-City Str. 1, A-10 Vienna, Austria katerina.schindler@assoc.oeaw.ac.at Abstract We presented a new lower bound on minimal singular values of real matrices based on Frobenius norm and determinant and showed in [4] that under certain assumptions on the matrix is our estimate sharper than a recent lower bound from Hong and Pan [3]. In this paper we show, under which conditions is our lower bound sharper than two other recent lower bounds for minimal singular values based on a matrix norm and determinant, namely the bound from Piazza and Politi and the bound from Hou-Biao Li et al. Keywords: real non-singular matrix, singular values, lower bound 1 Introduction The singular values or eigenvalues of real matrices are fundamental quantities in matrix analysis and other mathematical fields. The term singular value relates to the distance between a matrix and the set of singular matrices. The singular values are difficult to evaluate in general. In some cases at least the knowledge of their approximate values is valuable. The first bounds for eigenvalues were obtained more than a hundred years ago. The first paper using traces in eigenvalue inequalities was from Schur in 1909 [9]. Possibly the best-known inequality on eigenvalues is from Gerschgorin in 1931 []. Recently, several other lower bounds have been proposed for the smallest singular value of a square matrix, such as Johnson s bound, Brauer-type bound, Li s bound and and Ostrowski- type bound [6, 7, 1, 8, 13]. In this paper we deal only with the bounds on the minimal singular values using a matrix norm or determinant. We compare our bound from [4] with two other known bounds using a matrix norm or determinant.

3196 Kateřina Hlaváčková-Schindler Preliminary Notes Let A be a n n, n matrix with real elements. Let A E =( n i,j=1 a ij 1/ be the Frobenius norm of matrix A. Trace of a n n matrix A denotes tr(a = n i=1 a ii. The spectral norm of the matrix A is A = max 1 i n λ i, where λ i is eigenvalue of A T A.Ifλ 1,...,λ n are the eigenvalues of the matrix A, then deta = λ 1 λ...λ n. Denote the smallest singular value of A by σ n and its largest singular value by σ 1. It holds that A E = n i=1 σi = tr(a T A, where trace tr(a T A= n i=1 σi. For symmetric positive definite matrix A, one can suppose A = 1, i.e. that the matrix A is normalized and σ max = 1. The matrix normalization can be always achieved by multiplying the set of equations Ax = b by a suitable constant or for example by the divisive normalization defined by Weiss [1] or Ng et al. [10]. G. Piazza and T. Politi gave in 00 [11] a lower bound on the minimal singular matrix with positive singular values as σ min (n / A E. (1 Hou-Biao Li et al. gave in 010 in [5] another lower bound on minimal singular bound ( n 1 n 1 σ min. ( A E These two lower bounds on the minimal singular value based on determinant and matrix norm will be compared with our bound given in the next section. 3 Main Results 3.1 Our lower bound In [4] we proved the following theorem. Theorem 3.1 Let A be a nonsingular matrix with singular values σ i so that σ max = σ 1... σ n = σ min and let σ max σ min. Let A E = ( n i,j=1 a ij 1/ be the Frobenius norm of matrix A. (i Then for its minimal and maximal singular values holds 0 < A E 1/ nσ max <σ min. (3 n(1 σ max /n

Bounds on minimal singular value 3197 (ii For σ max =1(supposing that 1 holds 0 < ( /n ( A 1/ n <σ min. n( /n 1 3. Comparison with the lower bound from Piazza and Politi We will show that our estimate of minimal singular value is under certain conditions sharper than the estimate from from Piazza and Politi. Theorem 3. Let for a n n- nonsingular matrix A with n>1 hold for its singular values σ i,σ max = σ 1 =1... σ n = σ min and let σ max σ min and A E 0. If for matrix A holds deta < (n n n (n then < A E 1/ nσ max (4 (n / A E n(1 σ max /n i.e. our estimate is sharper than the estimate from Piazza and Politi. Proof: Formula (4 in power is (n A E < A F nσ max. (5 n(1 σ max /n Since between arithmetical mean A(σi and geometrical mean G(σi ofσ1,...,σ n holds A(σi G(σ i then is X = ( 1 n A E σ max 1. The right side of (5 is ( /n σmax /n. ( 1 n A E σ max ( /n σmax. It will then suffice to prove that Since A(σ i G(σ i then also /n 1 n A and n < 1. (6 A n 1 A 1. (7 n /n The left side of (6 is, which using (7 is smaller or equal to A /n n. The last value is smaller than 1 because of the condition on deta n 4/n n from the statement of the theorem.

3198 Kateřina Hlaváčková-Schindler 3.3 Comparison with the lower bound from Hou-Biao Li et al. We will show that under certain conditions is our lower bound of minimal singular value sharper than the lower bound from Hou-Biao Li et al. Theorem 3.3 Let for a n n- nonsingular matrix A with n>1 hold for its singular values σ i,σ max = σ 1 =1... σ n = σ min and let σ max σ min and A E 0. If between, A E holds (n 1 1 1 n < A E then ( n 1 A (n 1/ < A E nσmax n(1 σ max 1/ /n (8 i.e. our estimate is sharper than the estimate from Hou-Biao Li et al. Proof: Denote P = X. /n = ( 1 n A E σ max, L = ( (n 1 n 1 ( /n σmax /n A. From the proof of the previous theorem we know that X = ( 1 n A E σ max 1. ( /n σmax One is to prove that under the given conditions for A E and is L<P. It is obvious that 1 n = 1 n 1 (1 1. (9 n From the condition (n 1 1 1 n < A E follows (n 1 n < A E which is under the rewriting of by means of (9 n equivalent to (n 1( 1 n n 1 < A E. (10 The (n 1-th power of (10 gives (n 1 n 1 ( n < A (n 1 E, (11 from which is ( n 1 (n 1 L = < /n. (1 A E Since X 1, is also L< /n.x = P which was to be proven. ACKNOWLEDGEMENTS. The author is grateful for the support by the projects of the Czech Academy of Sciences GACR 10/08/0567 and MSMT CR C06001.

Bounds on minimal singular value 3199 References [1] A. Brauer, I.C. Gentry: Bounds for the greatest characteristic root of an irreducible nonnegative matrix. Linear Algebra Appl., 8:105-107, 1974. [] S. Gerschgorin: Über die Abgrenzung der Eigenwerte einer Matrix. Izv. Akad. Nauk. USSR Otd. Fiz.-Mat. Nauk 7, 749754, 1931. [3] Y.P. Hong, C. T. Pan.: A lower bound for the smallest singular value, Linear Algebra Appl., 17:7-3 (199. [4] K. Hlaváčková-Schindler: A new lower bound for the minimal singular value for real non-singular matrices by a matrix norm and determinant, submitted to Applied Mathematical Sciences. [5] Hou-Biao Li, Ting-Zhu Huang and Hong Li: Some new results on determinantal inequalities and applications, Journal of Inequalities and Applications, Vo. 010, Article ID 847357, 16 pages, 010. [6] C.R. Johnson: A Gersgorin type lower bound for the smallest singular value, Linear Algebra and Its Applications Vol. 11, pp.1-7, 1989. [7] C.R. Johnson, T. Szulc: Further lower bounds for the smallest singular value, Linear Algebra and Its Applications, Vol. 7, pp. 169-179, 1998. [8] L. Li: Lower bounds dor the smallest singular value, Computers and Mathematics with Applications, Vol. 41, pp.483-487, 001. [9] L. Mirsky: An Introduction to Linear Algebra, Oxford U.P., 1955. [10] A.Y. Ng, M.I. Jordan, Y. Weiss: On spectral clustering: Analysis and an algorithm. In S. Becker, T. G. Dietterich, and Z. Ghahramani (Eds., Advances in neural information processing systems. Cambridge, MA: MIT Press, 00. [11] G. Piazza, T. Politi: An upper bound for the condition number of a matrix in spectral norm. Journal of Computational and Appl. Math., 143, 141-144, 00. [1] Y. Weiss: Segmentation using eigenvectors: A unifying view. In Proceedings IEEE International Conference on Computer Vision (Vol. 9, pp. 97598. Piscataway, NJ: IEEE Press, 1999. [13] Y. Yamamoto: Equality Conditions for Lower Bounds on the Smallest Singular Value of a Bidiagonal Matrix, Applied Mathematics and Computation, Vol. 00, Issue 1, June 008, pp.54-60. Received: May, 010