Improving the Localization of Eigenvalues for Complex Matrices

Similar documents
Eigenvalue localization for complex matrices

Generalization of Samuelson s inequality and location of eigenvalues

TRACES OF HADAMARD AND KRONECKER PRODUCTS OF MATRICES. 1. Introduction

The Method of Least Squares. To understand least squares fitting of data.

On Random Line Segments in the Unit Square

5.1. The Rayleigh s quotient. Definition 49. Let A = A be a self-adjoint matrix. quotient is the function. R(x) = x,ax, for x = 0.

A Hadamard-type lower bound for symmetric diagonally dominant positive matrices

Bounds for the Extreme Eigenvalues Using the Trace and Determinant

A collocation method for singular integral equations with cosecant kernel via Semi-trigonometric interpolation

Stochastic Matrices in a Finite Field

Some properties of Boubaker polynomials and applications

62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +

Definitions and Theorems. where x are the decision variables. c, b, and a are constant coefficients.

Optimization Methods MIT 2.098/6.255/ Final exam

Symmetric Division Deg Energy of a Graph

International Journal of Mathematical Archive-3(4), 2012, Page: Available online through ISSN

ON POINTWISE BINOMIAL APPROXIMATION

Bangi 43600, Selangor Darul Ehsan, Malaysia (Received 12 February 2010, accepted 21 April 2010)

A Note On The Exponential Of A Matrix Whose Elements Are All 1

k-generalized FIBONACCI NUMBERS CLOSE TO THE FORM 2 a + 3 b + 5 c 1. Introduction

The Choquet Integral with Respect to Fuzzy-Valued Set Functions

Matrix Theory, Math6304 Lecture Notes from November 27, 2012 taken by Charles Mills

Sequences of Definite Integrals, Factorials and Double Factorials

PAijpam.eu ON TENSOR PRODUCT DECOMPOSITION

Lecture #20. n ( x p i )1/p = max

A NOTE ON THE TOTAL LEAST SQUARES FIT TO COPLANAR POINTS

COMMON FIXED POINT THEOREMS VIA w-distance

You may work in pairs or purely individually for this assignment.

Taylor polynomial solution of difference equation with constant coefficients via time scales calculus

Computation of Error Bounds for P-matrix Linear Complementarity Problems

ROSE WONG. f(1) f(n) where L the average value of f(n). In this paper, we will examine averages of several different arithmetic functions.

A class of spectral bounds for Max k-cut

NEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE

The picture in figure 1.1 helps us to see that the area represents the distance traveled. Figure 1: Area represents distance travelled

Lecture 2 Clustering Part II

Yuki Seo. Received May 23, 2010; revised August 15, 2010

Pairs of disjoint q-element subsets far from each other

SOME RELATIONS ON HERMITE MATRIX POLYNOMIALS. Levent Kargin and Veli Kurt

A solid Foundation for q-appell Polynomials

5.1 Review of Singular Value Decomposition (SVD)

A Simplified Binet Formula for k-generalized Fibonacci Numbers

SPECTRUM OF THE DIRECT SUM OF OPERATORS

Numerical Conformal Mapping via a Fredholm Integral Equation using Fourier Method ABSTRACT INTRODUCTION

Period Function of a Lienard Equation

Dec Communication on Applied Mathematics and Computation Vol.32 No.4

arxiv: v1 [math.fa] 3 Apr 2016

PAPER : IIT-JAM 2010

A GENERALIZATION OF THE SYMMETRY BETWEEN COMPLETE AND ELEMENTARY SYMMETRIC FUNCTIONS. Mircea Merca

The value of Banach limits on a certain sequence of all rational numbers in the interval (0,1) Bao Qi Feng

x c the remainder is Pc ().

3.2 Properties of Division 3.3 Zeros of Polynomials 3.4 Complex and Rational Zeros of Polynomials

TEACHER CERTIFICATION STUDY GUIDE

On the Jacobsthal-Lucas Numbers by Matrix Method 1

New Bounds for the Resolvent Energy of Graphs

AMS Mathematics Subject Classification : 40A05, 40A99, 42A10. Key words and phrases : Harmonic series, Fourier series. 1.

ON THE LEHMER CONSTANT OF FINITE CYCLIC GROUPS

Math 61CM - Solutions to homework 3

Numerical Solution of the First-Order Hyperbolic Partial Differential Equation with Point-Wise Advance

A 2nTH ORDER LINEAR DIFFERENCE EQUATION

Introducing a Novel Bivariate Generalized Skew-Symmetric Normal Distribution

The log-behavior of n p(n) and n p(n)/n

A constructive analysis of convex-valued demand correspondence for weakly uniformly rotund and monotonic preference

G. R. Pasha Department of Statistics Bahauddin Zakariya University Multan, Pakistan

Asymptotic distribution of products of sums of independent random variables

It is often useful to approximate complicated functions using simpler ones. We consider the task of approximating a function by a polynomial.

Estimation of Backward Perturbation Bounds For Linear Least Squares Problem

Ma 4121: Introduction to Lebesgue Integration Solutions to Homework Assignment 5

The inverse eigenvalue problem for symmetric doubly stochastic matrices

Appendix to Quicksort Asymptotics

Commutativity in Permutation Groups

Basics of Probability Theory (for Theory of Computation courses)

Solution of Differential Equation from the Transform Technique

Complex Analysis Spring 2001 Homework I Solution

Review Problems 1. ICME and MS&E Refresher Course September 19, 2011 B = C = AB = A = A 2 = A 3... C 2 = C 3 = =

Riesz-Fischer Sequences and Lower Frame Bounds

CLOSED FORM FORMULA FOR THE NUMBER OF RESTRICTED COMPOSITIONS

Proof of Fermat s Last Theorem by Algebra Identities and Linear Algebra

A new error bound for linear complementarity problems for B-matrices

PUTNAM TRAINING INEQUALITIES

Concavity Solutions of Second-Order Differential Equations

Algebra of Least Squares

Lecture 20. Brief Review of Gram-Schmidt and Gauss s Algorithm

A Study on Some Integer Sequences

Linear Regression Demystified

Properties of Fuzzy Length on Fuzzy Set

PRACTICE FINAL/STUDY GUIDE SOLUTIONS

A REMARK ON A PROBLEM OF KLEE

The multiplicative structure of finite field and a construction of LRC

THE ASYMPTOTIC COMPLEXITY OF MATRIX REDUCTION OVER FINITE FIELDS

The Growth of Functions. Theoretical Supplement

On forward improvement iteration for stopping problems

Recurrence Relations

17 Phonons and conduction electrons in solids (Hiroshi Matsuoka)

Bull. Korean Math. Soc. 36 (1999), No. 3, pp. 451{457 THE STRONG CONSISTENCY OF NONLINEAR REGRESSION QUANTILES ESTIMATORS Seung Hoe Choi and Hae Kyung

Topics in Eigen-analysis

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j

Disjoint Systems. Abstract

THE KALMAN FILTER RAUL ROJAS

Poincaré Problem for Nonlinear Elliptic Equations of Second Order in Unbounded Domains

Riemann Sums y = f (x)

Transcription:

Applied Mathematical Scieces, Vol. 5, 011, o. 8, 1857-1864 Improvig the Localizatio of Eigevalues for Complex Matrices P. Sargolzaei 1, R. Rakhshaipur Departmet of Mathematics, Uiversity of Sista ad Baluchesta Zaheda, Ira Abstract I this paper, we give a ew boud for λ i, which improves Schur, Kress, ad Eberlei s iequalities. By usig this estimate, we provide disks ad rectagles which iclude the eigevalues of a complex matrix ad the we compare these regios. Fially, we improve these regios by a strip for matrices with characteristic polyomial icludig real coefficiets. Numerical examples are provided to illustrate the results i this cotributio. Mathematics Subject Classificatio: 65F15, 65B99 Keywords: Eigevalues, localizatio, complex matrices, rectagular regio 1 Itroductio Matrices are a key tool i liear algebra. Matrices fid may applicatios, for example, graph theory uses matrices to keep track of distaces betwee pairs of vertices i a graph. Computer graphics uses matrices to project - dimesioal space oto a -dimesioal scree. Matrix calculus geeralizes classical aalytical otios such as derivatives of fuctios or expoetials to matrices. Due to these, A major brach of umerical aalysis is devoted to the developmet of efficiet algorithms for matrix computatios, a subject that is ceturies old but still a active area of research. Let A be a complex matrix of order with eigevalues λ 1,λ,..., λ. Gersgori doe almost first step for localizatio of eigevalues. He showed each 1 Correspodig Author. E-mail: sargolzaei@hamoo.usb.ac.ir

1858 P. Sargolzaei ad R. Rakhshaipur eigevalue λ of A satisfies at least i oe of the followig disks: λ a ii a ij,,,...,. j=1,j i Wolkowicz i [6, 7] provided upper ad lower bouds for (Rλ i ) ad (Iλ i ) ad the Rojo [, 4] provided regios that are cotai eigevalues. More precisely, Eberlei i [1] showed that λ i A AA A A = k 6 A 1 (A), (1) while Kress [] gave a tighter boud λ i ( A 4 AA A A ) 1/ = k (A). () Note that we use. as the Euclidea orm throughout. The rest of this paper is orgaized as comes ext. I Sectio, based o the available literature, we review a upper boud for λ i, which is tighter tha k (A). I sectio, we oce agai review the disks ad rectagles which cotais of eigevalues, accordig to the literature. The, we improve these regios by a strip for matrix with characteristic polyomial icludig real coefficiets. I fact, the cotributio of this paper will be furished i Sectio.1. The paper eds with several examples i Sectio 4. Prerequisites for λ i Theorem.1. Let A be a complex matrix of order with eigevalues λ 1,λ,..., λ. The λ i k (A), () where k (A) =(( A tra ) AA A A ) 1/ + tra. Proof. See [8]. Lemma.. Let A be a complex matrix of order. The k (A) k (A) k 1 (A). (4) Equality holds o the left of (4) if ad oly if A is ormal or tra=0. Equality holds o the right of (4) if ad oly if A is ormal.

Improvig the localizatio of eigevalues for complex matrices 1859 Proof. See [8]. Disks ad rectagles cotai eigevalues We will use the followig lemma, that is proved i [5]: Lemma.1. If z 1,z,..., z are complex umber, the disk z 1 z i 1 ( z i 1 z i ) give boud for them. Theorem.. Let A be a complex matrix of order with eigevalues λ 1,λ,..., λ. The, all of the eigevalues lie i the disks for j =1,,. λ tra ( 1 )1/ (k j (A) tra )1/, (5) Proof. Let z i = λ i i the lemma.1 ad use (1), () ad (). By lemma., it is trivial that disk (5) for j = is tighter. Later we will give rectagular regios that are cotai eigevalues. Accordig to Schur theorem, there exists a uitary U such that U AU = T = D + N where T is a upper triagular matrix, N is a strictly upper triagular matrix ad D = diag(λ 1,λ,,λ ). Let B = 1 (A + A ) ad C = 1 i (A A ). We compute B, C : B = 1 4 T + T = 1 4 (4 (Rλ i ) + N ), therefore ad similarly B = C = (Rλ i ) + 1 N, (Iλ i ) + 1 N.

1860 P. Sargolzaei ad R. Rakhshaipur From the two above relatios, we obtai B C = (Rλ i ) (Iλ i ) = RtrA. O the other had relatio ad (1), (), () ad for j =1,,, we obtai λ i = (Rλ i ) + (Iλ i ). Now, by usig this (Rλ i ) 1 (k j(a)+rtra ) (6) (Iλ i ) 1 (k j(a) RtrA ). (7) Lemma.. If d 1,d,..., d are real umber such that d 1 +d +...+d =0, the d j 1 d k,j =1,,...,. Proof. If d j = 0, we are doe. Suppose x i = d i d j. For i j defie f(x 1,..., x j 1,x j+1,..., x )=1+ x i. Sice x i =1, therefore, f has miimum i j i j i the poit ( 1,..., 1 ), the we have 1 1 1 f( 1,..., 1 1 )= 1 1+ i j x i = d i d j. Theorem.4. Let A be a complex matrix of order, with eigevalues λ 1,λ,..., λ. The all of the eigevalues lie i the rectagle where ad [ RtrA α, RtrA α =( 1 )1/ ( β =( 1 )1/ ( + α] [ ItrA ItrA β, + β], (Rλ i ) (RtrA) ) 1/ (8) (Iλ i ) (ItrA) ) 1/. (9)

Improvig the localizatio of eigevalues for complex matrices 1861 Proof. The real parts ad imagiary parts of umbers λ 1 tra,..., λ tra, satisfy the coditio of lemma.. Therefore, we apply twice this iequality to obtai (Rλ j RtrA ) 1 (Rλ k RtrA ) = 1 ( (Rλ k ) (RtrA) ) ad (Iλ j ItrA ) 1 1 (Iλ k ItrA ) = 1 ( (Iλ k ) (ItrA) ). The rectagle is icluded i the disk defied (5). Because α + β = ( λ i tra ). Also the vertices of the rectagle are o the boudary disk. Rectagle ca be computed without kowig the eigevalues. Thus, we attai that where ad [ RtrA α j, RtrA α j =( 1 β j =( 1 + α j] [ ItrA for j =1,,. Above rectagle for j = is tighter. β j, ItrA + β j ], (10) )1/ ( 1 (k j(a)+rtra ) (RtrA) ) 1/, (11) )1/ ( 1 (k j(a) RtrA ) (ItrA) ) 1/, (1).1 Mai cotributio: Improved regio I order to cotribute ad provide better regios, we keep goig as follows. Now, we suppose A be a complex matrix with characteristic polyomial havig real coefficiets. Therefore if λ be a eigevalue of A, the λ will be a eigevalue of A too. Thus Iλ i ( 1 (Iλ k ) ) 1/, (1) for i =1,,...,. Take j =, the we obtai Iλ 1 (k (A) tra ). (14) The improved eigevalue localizatio regio will be the itersectio of the rectagular regios with the strip defied above.

186 P. Sargolzaei ad R. Rakhshaipur 4 Numerical examples Example 1. Let A = are i 1+i 0 1+i 1 1 4i. Three rectagles of Theorem.4 k 1 =4.0940, Rλ 1.1675, Iλ 5.66, k =0.774, Rλ 1 1.8950, Iλ 5.0609, k =5.1697, Rλ 1 1.17, Iλ 5.788. Also three disks Theorem. are λ 1 5 i.8954, λ 1 5 i.6, λ 1 5 i.071. We have show three above rectagles i figure (1) ad last rectagle with last disk i Figure (), (eigevalues was show with *). 5 4 1 0 1 1.5 1 0.5 0 0.5 1 1.5.5.5 Figure 1: Rectagles of Example 1

Improvig the localizatio of eigevalues for complex matrices 186 5 4 1 0 1 1 0 1 4 5 Figure : Rectagle ad disk of Example 1 1 0 Example. Let B = i i. For this matrix we have 8i 1 0 k (A) =6.57, Rλ.5546, Iλ.951, ad λ 4.8469. I additio, sice characteristic polyomial of B has real coefficiets, by usig strip (14) we obtai Iλ.856. We have show above rectagle ad disk with strip (14) i Figure (): 5 4 1 0 1 4 5 6 4 0 4 6 Figure : Rectagle ad disk of Example with strip (14)

1864 P. Sargolzaei ad R. Rakhshaipur 5 Cocludig remarks The localizatio of the eigevalues of the complex matrices has studied ad ew improved regios were preseted i this work. For future extesios, we ca suggest that oe ca fid ovel relatios betwee the average ad variace of the eigevalues of the complex matrices to fid the sharper areas cotaiig the eigevalues. Refereces [1] P.J. Eberlei, O measures of o-ormality for matrices, Amer. Math. Moth. 7(1965) 995-996. [] R. Kress, H.L.De Vries, R. Wegma, O oormal matrices, Liear Algebra Appl. 8(1974) 109-10. [] O. Rojo, R. Soto, T. Avila, H. Rojo, Localizatio of eigevalues i elliptic regiocs, Comput. Math. Appl. 9(7) (1995) -11. [4] O. Rojo, R. Soto, H. Rojo, New eigevalue estimates for complex matrices, Comput. Math. Appl. 5() (199) -11. [5] O. Rojo, R. Soto, H. Rojo, T.Y. Tam, Eigevalue localizatio for complex matrices, Proyeccioes, 11(1) (199) 11-19. [6] H. Wolkowicz, G.P.H. Stya, Bouds for eigevalues usig traces, Liear Algebra Appl. 9(1980) 471-506. [7] H. Wolkowicz, G.P.H. Stya, More bouds for eigevalues usig traces, Liear Algebra Appl. 1(1980) 1-17. [8] T-Z Huag, L. Wag, Improvig bouds for eigevalues of complex matrices usig traces, Liear Algebra Appl. 46(007) 841-854. Received: December, 010