Dependence of simple eigenpairs to differentiable ...aalexan3/articles/eigen.pdf · Dependence of...

2
Dependence of simple eigenpairs to differentiable perturbations Alen Alexanderian * Abstract Let A : R n R n×n be a differentiable matrix valued function. Suppose A(0) has a real simple eigenvalue λ0. We show here that under for ξ sufficiently close to the origin, A(ξ) has a real eigenvalue λ(ξ), and that λ(ξ) depends differentiably on ξ.A related result regarding the eigenvectors is also considered. 1 Introduction Consider a matrix valued function A : R n R n×n , A = A(ξ) with ξ R n . Here we show that if A(0) has a real simple eigenvalue λ 0 and that A depends differentiably on ξ, then there exists a neighborhood U R n of ξ = 0, such that for each ξ U , A(ξ) has a real eigenvalue λ(ξ) which depends differentiably on ξ. Moreover, under the same assumption, we have a unit eigenvector v(ξ) R n for each ξ U corresponding to λ(ξ), and v depends differentiably on ξ in U . The idea is mainly based on that of [1] with a slight modification of considering R as the base field. 2 Implicit Function Theorem Here we state the following special case of the Implicit Function Theorem (see e.g. [3, 2]) which is used in this note. Theorem 2.1. Let D R n × R be an open set and F : D R a C 1 function. Assume that (x 0 ,y 0 ) D is such that F (x 0 ,y 0 )=0, and ∂y F (x 0 ,y) y=y0 6=0. Then, there exists an open set U R n containing x 0 and an open set V R containing y 0 with U × V D, and a unique C 1 function h : U V such that, h(x 0 )= y 0 and F (x,h(x)) = 0, x U. 3 Main results Here we present the main results we are concerned with in this note. The idea behind the proof of the below Theorem is from [1]. The only difference here is that we restrict the base field to reals and consider real eigenvalues. * The University of Texas at Austin, USA. E-mail: [email protected] Last revised: November 12, 2013

Transcript of Dependence of simple eigenpairs to differentiable ...aalexan3/articles/eigen.pdf · Dependence of...

Page 1: Dependence of simple eigenpairs to differentiable ...aalexan3/articles/eigen.pdf · Dependence of simple eigenpairs to differentiable perturbations ... Jerrold E. Marsden and Michael

Dependence of simple eigenpairs to differentiable

perturbations

Alen Alexanderian∗

Abstract

Let A : Rn → Rn×n be a differentiable matrix valued function. Suppose A(0) hasa real simple eigenvalue λ0. We show here that under for ξ sufficiently close to theorigin, A(ξ) has a real eigenvalue λ(ξ), and that λ(ξ) depends differentiably on ξ. Arelated result regarding the eigenvectors is also considered.

1 Introduction

Consider a matrix valued function A : Rn → Rn×n, A = A(ξ) with ξ ∈ Rn. Here weshow that if A(0) has a real simple eigenvalue λ0 and that A depends differentiably onξ, then there exists a neighborhood U ⊂ Rn of ξ = 0, such that for each ξ ∈ U , A(ξ) hasa real eigenvalue λ(ξ) which depends differentiably on ξ. Moreover, under the sameassumption, we have a unit eigenvector v(ξ) ∈ Rn for each ξ ∈ U corresponding to λ(ξ),and v depends differentiably on ξ in U . The idea is mainly based on that of [1] with aslight modification of considering R as the base field.

2 Implicit Function Theorem

Here we state the following special case of the Implicit Function Theorem (seee.g. [3, 2]) which is used in this note.Theorem 2.1. Let D ⊂ Rn × R be an open set and F : D → R a C1 function. Assumethat (x0, y0) ∈ D is such that

F (x0, y0) = 0, and∂

∂yF (x0, y)

∣∣∣y=y0

6= 0.

Then, there exists an open set U ⊂ Rn containing x0 and an open set V ⊂ R containingy0 with U × V ⊂ D, and a unique C1 function h : U → V such that,

h(x0) = y0 and F (x, h(x)) = 0, ∀x ∈ U.

3 Main results

Here we present the main results we are concerned with in this note. The ideabehind the proof of the below Theorem is from [1]. The only difference here is that werestrict the base field to reals and consider real eigenvalues.

∗The University of Texas at Austin, USA. E-mail: [email protected] revised: November 12, 2013

Page 2: Dependence of simple eigenpairs to differentiable ...aalexan3/articles/eigen.pdf · Dependence of simple eigenpairs to differentiable perturbations ... Jerrold E. Marsden and Michael

Differentiable perturbations of simple eigenpairs

Theorem 3.1. Let A : Rn → Rn×n be a differentiable matrix valued function. SupposeA(0) has a real simple eigenvalue λ0. Then, there exists a neighborhood U ⊂ Rn ofthe origin, such that for each ξ ∈ U , A(ξ) has a real eigenvalue λ(ξ). Moreover, λ(ξ)depends differentiably on ξ.

Proof. Consider the function p : Rn ×R→ R given by,

p(ξ, y) = det(A(ξ)− yI).

Note that by the assumptions on λ0, we have p(0, λ0) = 0. Moreover, we note thatp(0, y) = (y − λ0)r(y), with r(y) a polynomial of degree n − 1. Since λ0 is a simpleeigenvalue, i.e. a simple root of p(0, y), it follows that ∂

∂yp(0, y) |y=λ0= r(λ0) 6= 0.

Thus, we have the following:

p(0, λ0) = 0,∂p

∂y(0, y)

∣∣y=λ0

6= 0.

Therefore, by the Implicit Function Theorem, there exists a neighborhood U ⊂ Rn of0, a neighborhood V ⊂ R of λ0, and a C1 function λ : U → V such that λ(0) = λ0 andp(ξ, λ(ξ)) = 0 for every ξ ∈ U . That is

det(A(ξ)− λ(ξ)I) = 0, ∀ ξ ∈ U,

so that λ(ξ) is an eigenvalue of A(ξ) with λ(ξ) ∈ V ⊂ R.

Let A(ξ) be as in the above theorem. For ξ = 0, we have the real simple eigenvalueλ0. Let v0 ∈ Rn be a unit eigenvector corresponding to λ0, and note that choice of v0is unique up to a sign, as the eigenspace corresponding to λ0 is one-dimensional. Wemake the choice of v0 be fixed such that v0 makes a positive angle with the vector e =

(1, 1, . . . , 1) ∈ Rn. Note that e could have been chosen to be any other vector in Rn, thechoice here was for specificity. For each ξ ∈ U , where U is as in the above Theorem, wehave a simple real eigenvalue λ(ξ) ∈ Rn; the choice of the corresponding eigenvectorcan be made v(ξ) can be made unique by first normalizing it, and then fixing the anglebetween v(ξ) and v0 to be in (−π/2, π/2] as follows. Let v(ξ) be a unit eigenvectorcorresponding to λ(ξ) if they are perpendicular, rotating v if necessary, we make theangle between v and v0 to be π/2; otherwise, we multiply v(ξ) by sign

(v(ξ) · v0

). We

call this unique choice of v(ξ) the positively oriented eigenvector corresponding to λξ.

Using Theorem 8 in [1, p. 102], we immediately get the following:Theorem 3.2. Let A(ξ) be as in the previous theorem with λ(ξ) the correspondingeigenvalue for ξ ∈ U as stated there. For each ξ ∈ U , let the v(ξ) ∈ Rn be the cor-responding positively oriented eigenvector. Then, v(ξ) depends differentiably on ξ inU .

References

[1] Peter D. Lax. Linear Algebra. Wiley-Interscience, 1996.

[2] Jerrold E. Marsden and Michael J. Hoffman. Elementary classical analysis. W.H. Freeman andCompany, 2nd edition, 1993.

[3] Walter Rudin. Principles of Mathematical Analysis. McGraw-Hill Book Company, 2nd edition,1964.

2