Specral Graph Theory and its Applications September 7, Lecture 2 L G1,2 = 1 1.

Similar documents
An Algorithmist s Toolkit September 15, Lecture 2

Spectral Graph Theory Lecture 3. Fundamental Graphs. Daniel A. Spielman September 5, 2018

Fiedler s Theorems on Nodal Domains

An Algorithmist s Toolkit September 10, Lecture 1

Fiedler s Theorems on Nodal Domains

Spectral Graph Theory Lecture 2. The Laplacian. Daniel A. Spielman September 4, x T M x. ψ i = arg min

Rings, Paths, and Cayley Graphs

Rings, Paths, and Paley Graphs

Lecture Introduction. 2 Brief Recap of Lecture 10. CS-621 Theory Gems October 24, 2012

1 Matrix notation and preliminaries from spectral graph theory

Lecture 12: Introduction to Spectral Graph Theory, Cheeger s inequality

ON THE QUALITY OF SPECTRAL SEPARATORS

Ma/CS 6b Class 23: Eigenvalues in Regular Graphs

Lecture 1. 1 if (i, j) E a i,j = 0 otherwise, l i,j = d i if i = j, and

The Simplest Construction of Expanders

Strongly Regular Graphs, part 1

Lecture 10. Lecturer: Aleksander Mądry Scribes: Mani Bastani Parizi and Christos Kalaitzis

Laplacian Matrices of Graphs: Spectral and Electrical Theory

1 Review: symmetric matrices, their eigenvalues and eigenvectors

Laplacians of Graphs, Spectra and Laplacian polynomials

Laplacian Integral Graphs with Maximum Degree 3

Spectra of Adjacency and Laplacian Matrices

Linear algebra and applications to graphs Part 1

Ring Sums, Bridges and Fundamental Sets

Data Analysis and Manifold Learning Lecture 3: Graphs, Graph Matrices, and Graph Embeddings

Lecture 13: Spectral Graph Theory

Using Laplacian Eigenvalues and Eigenvectors in the Analysis of Frequency Assignment Problems

Laplacians of Graphs, Spectra and Laplacian polynomials

U.C. Berkeley Better-than-Worst-Case Analysis Handout 3 Luca Trevisan May 24, 2018

Algebraic Constructions of Graphs

Effective Resistance and Schur Complements

Graphs, Vectors, and Matrices Daniel A. Spielman Yale University. AMS Josiah Willard Gibbs Lecture January 6, 2016

Laplacian eigenvalues and optimality: II. The Laplacian of a graph. R. A. Bailey and Peter Cameron

Properties of Expanders Expanders as Approximations of the Complete Graph

Walks, Springs, and Resistor Networks

On Hadamard Diagonalizable Graphs

Lecture: Local Spectral Methods (3 of 4) 20 An optimization perspective on local spectral methods

Iterative solvers for linear equations

Lecture 7. 1 Normalized Adjacency and Laplacian Matrices

Quadrature for the Finite Free Convolution

Lectures 15: Cayley Graphs of Abelian Groups

j=1 u 1jv 1j. 1/ 2 Lemma 1. An orthogonal set of vectors must be linearly independent.

Sparsification by Effective Resistance Sampling

Data Analysis and Manifold Learning Lecture 7: Spectral Clustering

Convergence of Random Walks

Current; Forest Tree Theorem; Potential Functions and their Bounds

Physical Metaphors for Graphs

1 Matrix notation and preliminaries from spectral graph theory

Lecture 1 and 2: Random Spanning Trees

Algebra I Fall 2007

Iterative solvers for linear equations

Lecture 1: Systems of linear equations and their solutions

Spectral Theory of Unsigned and Signed Graphs Applications to Graph Clustering. Some Slides

CMPSCI 791BB: Advanced ML: Laplacian Learning

A lower bound for the Laplacian eigenvalues of a graph proof of a conjecture by Guo

Spectral Graph Theory and its Applications. Daniel A. Spielman Dept. of Computer Science Program in Applied Mathematics Yale Unviersity

Iterative solvers for linear equations

Kernels of Directed Graph Laplacians. J. S. Caughman and J.J.P. Veerman

MAT2342 : Introduction to Applied Linear Algebra Mike Newman, fall Projections. introduction

Lecture 14: Random Walks, Local Graph Clustering, Linear Programming

Note on deleting a vertex and weak interlacing of the Laplacian spectrum

Remark 3.2. The cross product only makes sense in R 3.

Introducing the Laplacian

An Algorithmist s Toolkit Nov. 10, Lecture 17

Math 113 Homework 5 Solutions (Starred problems) Solutions by Guanyang Wang, with edits by Tom Church.

5 Quiver Representations

Lecture 2: Linear Algebra Review

Review problems for MA 54, Fall 2004.

LINEAR ALGEBRA BOOT CAMP WEEK 4: THE SPECTRAL THEOREM

18.01 Single Variable Calculus Fall 2006

CS168: The Modern Algorithmic Toolbox Lectures #11 and #12: Spectral Graph Theory

Graph fundamentals. Matrices associated with a graph

Spectral and Electrical Graph Theory. Daniel A. Spielman Dept. of Computer Science Program in Applied Mathematics Yale Unviersity

Diffusion and random walks on graphs

SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS

Introduction. Disclaimer. 1.1 First Things. 1.2 Introduction. Spectral Graph Theory Lecture 1. Daniel A. Spielman September 2, 2015

An Algorithmist s Toolkit Lecture 18

Spectral Analysis of Random Walks

MATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators.

Spectral Algorithms I. Slides based on Spectral Mesh Processing Siggraph 2010 course

CS224W: Social and Information Network Analysis Jure Leskovec, Stanford University

Graph Embeddings and Laplacian Eigenvalues

Petar Maymounkov Problem Set 1 MIT with Prof. Jonathan Kelner, Fall 07

17.1 Directed Graphs, Undirected Graphs, Incidence Matrices, Adjacency Matrices, Weighted Graphs

Math 443/543 Graph Theory Notes 5: Graphs as matrices, spectral graph theory, and PageRank

Lecture 9: Low Rank Approximation

Lecture 1: Graphs, Adjacency Matrices, Graph Laplacian

Lecture 24: Element-wise Sampling of Graphs and Linear Equation Solving. 22 Element-wise Sampling of Graphs and Linear Equation Solving

Lecture 4 Eigenvalue problems

MATH 423 Linear Algebra II Lecture 20: Geometry of linear transformations. Eigenvalues and eigenvectors. Characteristic polynomial.

8.1 Concentration inequality for Gaussian random matrix (cont d)

6.854J / J Advanced Algorithms Fall 2008

(a) If A is a 3 by 4 matrix, what does this tell us about its nullspace? Solution: dim N(A) 1, since rank(a) 3. Ax =

Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science Algorithms For Inference Fall 2014

ORIE 6334 Spectral Graph Theory October 13, Lecture 15

Graph Partitioning Algorithms and Laplacian Eigenvalues

Partial cubes: structures, characterizations, and constructions

CS286.2 Lecture 13: Quantum de Finetti Theorems

6.045: Automata, Computability, and Complexity (GITCS) Class 15 Nancy Lynch

. The following is a 3 3 orthogonal matrix: 2/3 1/3 2/3 2/3 2/3 1/3 1/3 2/3 2/3

Transcription:

Specral Graph Theory and its Applications September 7, 2004 Lecturer: Daniel A. Spielman Lecture 2 2.1 The Laplacian, again We ll begin by redefining the Laplacian. First, let G 1,2 be the graph on two vertices with one edge. We define [ ] def 1 1 L G1,2 = 1 1. Note that x T L G1,2 x = (x 1 x 2 ) 2. (2.1) In general, for the graph with n vertices and just one edge between vertices u and v, we can define the Laplacian similarly. For concreteness, I ll call the graph G u,v and define it by 1 if i = j and i u, v L Gu,v (i, j) def = 1 if i = u and j = v, or vice versa, 0 otherwise. For a graph G with edge set E, we now define L(G) def = (u,v) E L(G u,v ). Many elementary properties of the Laplacian follow from this definition. In particular, we see that L G1,2 has eigenvalues 0 and 2, and so is positive semidefinite, where we recall that a symmetric matrix M is positive semidefinite if all of its eigenvalues are non-negative. Also recall that this is equivalent to x T Mx 0, for all x R n. It follows immediately that the Laplacian of every graph is positive semidefinite. One way to see this is to sum equation (2.1) to get x T L G x = (x u x v ) 2. (u,v) E Remark Since the vertex set really doesn t matter, I actually prefer the notation L(E) where E is a set of edges. Had I used this notation above, it would have eliminated some subscripts. For example, I could have written L{u,v} instead of L Gu,v. We now sketch the proof of one of the most fundamental facts about Laplacians. 2-1

Lecture 2: September 7, 2004 2-2 Lemma 2.1.1. Let G = (V, E) be a connected graph, and let λ 1 λ 2 λ n be the eigenvalues of its Laplacian. Then, λ 2 > 0. Proof. Let x be an eigenvector L G of eigenvalue 0. Then, we have and so L G x = 0, x T L G x = (x u x v ) 2 = 0. (u,v) E Thus, for each pair of vertices (u, v) connected by an edge, we have x u = x v. As the graph is connected, we have x u = x v for all pairs of vertices u, v, which implies that x is some constant times the all 1s vector. Thus, the eigenspace of 0 has dimension 1. This lemma led Fiedler to consider the magnitude of λ 2 as a measure of how well-connected the graph is. In later lectures, we will see just how good a measure it is. Accordingly, we often call λ 2 the Fiedler value of a graph, and v 2 its Fiedler vector. In general, it is pointless to consider unconnected graphs, because their spectra are just the union of the spectra of their components. For example, we have: Corollary 2.1.2. Let G = (V, E) be a graph, and let L be its Laplacian. Then, the multiplicity of 0 as an eigenvalue of L equals the number of connected components of G. 2.2 Some Fundamental Graphs We now examine the eigenvalues and eigenvectors of the Laplacians of some fundamental graphs. In particular, we will examine The complete graph on n vertices, K n, which has edge set {(u, v) : u = v}. The path graph on n vertices, P n, which has edge set {(u, u + 1) : 1 u < n}. The ring graph on n vertices, R n, which has all the edges of the path graph, plus the edge (1, n). The grid graph on n 2 vertices. To define this graph, we will identify vertices with pairs (u 1, u 2 ) where 1 u 1 n and 1 u 2 n. Each node (u 1, u 2 ) has edges to each node that differs by one in just one coordinate. Lemma 2.2.1. The Laplacian of K n has eigenvalue 0 with multiplicity 1 and n with multiplicity n 1.

Lecture 2: September 7, 2004 2-3 Proof. Let L be the Laplacian of K n. Let x be any vector orthogonal to the all 1s vector. Consider the first coordinate of Lx. It will be n 1 times x 1, minus x u = x 1, u>1 as x is orthogonal to 1. Thus, x is an eigenvector of eigenvalue n. Lemma 2.2.2. The Laplacian of R n has eigenvectors x k (u) = sin(2πku/n), and y k (u) = cos(2πku/n), for k n/2. Both of these have eigenvalue 2 2 cos(2πk/n). Note x 0 should be ignored, and y 0 is the all 1s vector. If n is even, then x n/2 should also be ignored. Proof. The best way to see that x k and y k are eigenvectors is to plot the graph on the circle using these vectors as coordinates. That they are eigenvectors is geometrically obvious. To compute the eigenvalue, just consider vertex u = n. Lemma 2.2.3. The Laplacian of P n has the same eigenvalues as R 2n, and eigenvectors for 0 k < n x k (u) = sin(πku/n + π/2n). Proof. This is our first interesting example. We derive the eigenvectors and eigenvalues by treating P n as a quotient of R 2n : we will identify vertex u of P n with both vertices u and 2n u of R 2n. Let z be an eigenvector of R 2n in which z(u) = z(2n u) for all u. I then claim that the first n components of z give an eigenvector of P n. To see why this is true, note that an eigenvector z of P n must satisfy for some λ 2z(u) z(u 1) z(u + 1) = λz(u), for 1 < u < n, and z(1) z(2) = λz(2) z(n) z(n 1) = λz(n). One can now immediately verify that the first n components of such a z satisfy these conditions. The only trickly parts are at the endpoints. Now, to obtain such a vector z, we take which is in the span of x k and y k. z k (u) = sin (2πku/(2n) + π/(2n)), The quotient construction used in this proof is an example of a generally applicable technique.

Lecture 2: September 7, 2004 2-4 2.2.1 Products of Graphs We will now determine the eigenvalues and eigenvectors of grid graphs. We will go through a more general theory in which we relate the spectra of products graphs to the spectra of the graphs of which they are products. Definition 2.2.4. Let G = (V, E) and H = (W, F ) be graphs. Then G H is the graph with vertex set V W and edge set ( ) (v 1, w), (v 2, w) where (v 1, v 2 ) E and ( ) (v, w 1 ), (v, w 2 ) where (w 1, w 2 ) F. Theorem 2.2.5. Let G = (V, E) and H = (W, F ) be graphs with Laplacian eigenvalues λ 1,..., λ n and µ 1,..., µ m, and eigenvectors x 1,..., x n and y 1,..., y m, respectively. Then, for each 1 i n and 1 j m, G H has an eigenvector z of eigenvalue λ i + µ j such that z(v, w) = x i (v)y j (w). Proof. To see that this eigenvector has the propper eigenvalue, let L denote the Laplacian of G H, d v the degree of node v in G, and e w the degree of node w in H. We can then verify that (Lz)(v, w) = (d v + e w )(x i (v)y j (w)) (x i (v 2 )y j (w)) (x i (v)y j (w 2 )) (v,v 2 ) E (w,w 2 ) F = (d v )(x i (v)y j (w)) (x i (v 2 )y j (w)) + (e w )(x i (v)y j (w)) (x i (v)y j (w 2 )) (v,v 2 ) E (w,w 2 ) F = y j (w) d v x i (v) (x i (v 2 )) + x i (v) e w y j (w) (x i (v)y j ) (v,v 2 ) E = y j (w)λ i x i )(v) + x i (v)µ j y j (w) = (λ i + µ j )(x i (v)y j (w)). (w,w 2 ) F Thus, the eigenvalues and eigenvectors of the grid graph are completely determined by the eigenvalues and eigenvectors of the path graph. This completely explains the nice spectral embedding of the n-by-n grid graph: its smallest non-zero eigenvalue has multiplicity two and its eigenspace is spanned by the vectors 1 v 2 and v 2 1. 2.3 Bounding Eigenvalues It is rare to be able to precisely determine the eigenvalues of a graph by purely combinatorial arguments. Usually, such arguments only provide bounds on the eigenvalues. The most commonly used tool for bounding eigenvalues is the Courant-Fischer characterization of eigenvalues:

Lecture 2: September 7, 2004 2-5 Theorem 2.3.1 (Courant-Fischer). Let A be an n-by-n symmetric matrix and let 1 k n. Let λ 1 λ 2 λ n be the eigenvalues of A. Then, x T Ax λ k = min max, (2.2) S of dim k x S x T x and x T Ax λ k = max min, (2.3) S of dim n k 1 x S xt x For example, this allows us to obtain an easy upper bound on λ 2 (P n ). Setting S = {x : x 1}, and applying (2.3), we obtain x T Ax λ 2 (L) min. x 1 x T x This allows us to obtain an easy upper bound on λ 2 (P n ). Consider the vector x such that x(u) = (n + 1) 2u. This vector satisfies x 1, so (x(u) x(u + 1)) 2 λ 2 (P n ) 1 u<n x(u) 2 u 1 u<n 2 2 = i(n + 1 2u) 2 4n n(n 2 /3) 12 =. n 2 So, we can easily obtain an upper bound on λ(p n ) that is of the right order of magnitude. It is more difficult to obtain lower bounds. In the next lecture, we will introduce a technique for proving lower bounds. We will apply the technique to the example of the complete binary tree. The complete binary tree on n = 2 d 1 nodes, B n, is the graph with edges of the form (u, 2u) and (u, 2u + 1) for u < n/2. To upper bound λ 2 (B n ), we construct a vector x as follows. We first set x(1) = 0, x(2) = 1, and x(3) = 1. Then, for every vertex u that we can reach from node 2 without going through node 1, we set x(u) = 1. For all the other nodes, we set x(u) = 1. We then have (u,v) B n (x u x v ) 2 = (x 1 x 2 ) 2 + (x 1 x 3 ) 2 = 2/n. x2 n For now, we leave with some lower bounds on λ n. u u Lemma 2.3.2. Let G = (V, E) be a graph and let vertex w V have degree d. Then, λ max (G) d.

Lecture 2: September 7, 2004 2-6 Proof. By (2.2), Let x be the vector given by Then, we have λ max (G) max xt L(G)x. x x T x { 1 u = w x(u) = 0 otherwise. (u,v) E (x u x v ) 2 d = = d. x 2 u 1 We can improve this lemma slightly. Lemma 2.3.3. Let G = (V, E) be a graph and let vertex w V have degree d. Then, λ max (G) d + 1. Proof. Let x be the vector given by d u = w x(u) = 1 if (u, w) E, 0 otherwise. Then, we have ( u,v ) E ) 2 d(d + 1) 2 = = d + 1. (x u x v x2 d 2 + d u Can we prove a larger lower bound? To find out, let s consider the star graph with n vertices, S n. It has edges of the form (1, i) for 2 i n. To determine its eigenvalues, we apply the following lemma Lemma 2.3.4. Let G = (V, E) be a graph, and let v 1 and v 2 be vertices of degree one that are both connected to another vertex w. Then, the vector x given by 1 u = v 1 x(u) = 1 u = v 2 0 otherwise, is an eigenvector of the Laplacian of G of eigenvalue 1. Proof. One can immediately verify that Lx = x.

Lecture 2: September 7, 2004 2-7 Thus, we see that S n has eigenvectors of eigenvalue 1 pairing vertices 2 and 3, 3 and 4, up to 2 and n. The space spanned by these vectors has dimension n 2. As the vector 1 is also an eigenvector, only one vector may remain. We don t even have to compute it to determine its eigenvalue: the sum of the eigenvalues is the trace of the matrix, which is 2(n 1). The n 1 eigenvectors we have found so far account for n 2 of this sum, so the last eigenvector must have eigenvalue n. This establishes that Lemma 2.3.4 is tight. In fact, this last eigenvector turns out to be the vector constructed in the proof of that lemma.

MIT OpenCourseWare http://ocw.mit.edu 18.409 Topics in Theoretical Computer Science: An Algorithmist's Toolkit Fall 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.