<?xml version="1.0" encoding="utf-8" standalone="no"?>
<!DOCTYPE document PUBLIC "-//CNX//DTD CNXML 0.5 plus MathML//EN" "http://cnx.rice.edu/cnxml/0.5/DTD/cnxml_mathml.dtd">
<document xmlns="http://cnx.rice.edu/cnxml" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:bib="http://bibtexml.sf.net/" id="m2116">
  
  <name xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">Eigenvalue Decomposition</name>
  
  <metadata xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">
  <md:version xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">2.9</md:version>
  <md:created xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">2001/03/23</md:created>
  <md:revised xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">2002/10/24</md:revised>
  <md:authorlist xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">
    <md:author xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="aca">
      <md:firstname xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">Thanos</md:firstname>
      
      <md:surname xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">Antoulas</md:surname>
      <md:email xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">aca@rice.edu</md:email>
    </md:author>
    <md:author xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="jps">
      <md:firstname xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">John</md:firstname>
      <md:othername xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">Paul</md:othername>
      <md:surname xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">Slavinsky</md:surname>
      <md:email xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">jps@alumni.rice.edu</md:email>
    </md:author>
  </md:authorlist>

  <md:maintainerlist xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">
    <md:maintainer xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="rha">
      <md:firstname xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">Roy</md:firstname>
      
      <md:surname xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">Ha</md:surname>
      <md:email xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">rha@rice.edu</md:email>
    </md:maintainer>
    <md:maintainer xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="lizychan">
      <md:firstname xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">Elizabeth</md:firstname>
      
      <md:surname xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">Chan</md:surname>
      <md:email xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">lizychan@rice.edu</md:email>
    </md:maintainer>
    <md:maintainer xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="aca">
      <md:firstname xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">Thanos</md:firstname>
      
      <md:surname xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">Antoulas</md:surname>
      <md:email xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">aca@rice.edu</md:email>
    </md:maintainer>
    <md:maintainer xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="jps">
      <md:firstname xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">John</md:firstname>
      <md:othername xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">Paul</md:othername>
      <md:surname xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">Slavinsky</md:surname>
      <md:email xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">jps@alumni.rice.edu</md:email>
    </md:maintainer>
  </md:maintainerlist>
  
  <md:keywordlist xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">
    <md:keyword xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">decomposition</md:keyword>
    <md:keyword xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">eigenvalue</md:keyword>
    <md:keyword xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">eigenvector</md:keyword>
    <md:keyword xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">multiplicity</md:keyword>
  </md:keywordlist>

  <md:abstract xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">(Blank Abstract)</md:abstract>
</metadata>


  <content xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">
    <para xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="p0">
      When we apply a matrix to a vector (i.e. multiply them
      together), the vector is transformed.  An interesting question
      to ask ourselves is whether there are any particular
      combinations of such a matrix and vector whose result is a new
      vector that is proportional to the original vector.  In math
      terminology, this question can be posed as follows: if we have a
      matrix <m:math><m:ci type="matrix">A</m:ci></m:math>: 
      <m:math>
<m:apply><m:mo>→</m:mo>
	<m:apply>
	  <m:power/>
	  <m:ci>ℝ</m:ci>
	  <m:ci>n</m:ci>
	</m:apply>
	<m:apply>
	  <m:power/>
	  <m:ci>ℝ</m:ci>
	  <m:ci>n</m:ci>
	</m:apply>
</m:apply>
      </m:math>,
      does there exist a vector
      <m:math>
	<m:apply>
	<m:in/>
	  <m:ci type="vector">x</m:ci>
	  <m:apply>
	    <m:power/>
	    <m:ci>ℝ</m:ci>
	    <m:ci>n</m:ci>
	  </m:apply>
	</m:apply>
      </m:math>
      and a scalar
      <m:math>
	<m:apply>
	  <m:in/>
	  <m:ci>λ</m:ci>
	  <m:complexes/>
	</m:apply>
      </m:math>
      such that
      <m:math>
	<m:apply>
	  <m:eq/>
	  <m:apply>
	    <m:times/>
	    <m:ci type="matrix">A</m:ci>
	    <m:ci type="vector">x</m:ci>
	  </m:apply>
	  <m:apply>
	    <m:times/>
	    <m:ci>λ</m:ci>
	    <m:ci type="vector">x</m:ci>
	  </m:apply>
	</m:apply>
      </m:math>?
      If so, then the complexity of
      <m:math>
	<m:apply>
	  <m:times/>
	  <m:ci type="matrix">A</m:ci>
	  <m:ci type="vector">x</m:ci>
	</m:apply>
      </m:math>
      is reduced.  It no longer must be thought of as a matrix
      multiplication; instead, applying <m:math><m:ci type="matrix">A</m:ci></m:math> to <m:math><m:ci type="vector">x</m:ci></m:math>
      has the simple effect of linearly scaling <m:math><m:ci type="vector">x</m:ci></m:math> by some scalar
      factor <m:math><m:ci>λ</m:ci></m:math>.</para>
    
    <para xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="p1">In this situation, where
      <m:math>
	<m:apply>
	  <m:eq/>
	  <m:apply>
	    <m:times/>
	    <m:ci type="matrix">A</m:ci>
	    <m:ci type="vector">x</m:ci>
	  </m:apply>
	  <m:apply>
	    <m:times/>
	    <m:ci>λ</m:ci>
	    <m:ci type="vector">x</m:ci>
	  </m:apply>
	</m:apply>
      </m:math>,
      <m:math><m:ci>λ</m:ci></m:math> is known as an eigenvalue and <m:math><m:ci type="vector">x</m:ci></m:math>
      is its associated eigenvector.  For a certain matrix, each one
      of its eigenvectors is associated with a particular (though not
      necessarily unique) eigenvalue. The word "eigen" is German and
      means "same"; this is appropriate because the vector <m:math><m:ci type="vector">x</m:ci></m:math> 
      after the matrix multiplication is the same as the original
      vector <m:math><m:ci type="vector">x</m:ci></m:math>, 
      except for the scaling factor.  The following two examples give actual possible values for the matrices, vectors, and values discussed in general terms above.</para>
    
    <equation xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="eq1">
      <m:math mode="block">
	<m:apply>
	  <m:eq/>
	  <m:apply>
	    <m:times/>
	    <m:matrix>
	      <m:matrixrow><m:cn>1</m:cn><m:cn>-1</m:cn></m:matrixrow>
	      <m:matrixrow><m:cn>-1</m:cn><m:cn>1</m:cn></m:matrixrow>
	    </m:matrix>
	    <m:vector>
	      <m:cn>1</m:cn>
	      <m:cn>1</m:cn>
	    </m:vector>
	  </m:apply>
	  <m:apply>
	    <m:times/>
	    <m:cn>0</m:cn>
	    <m:vector>
	      <m:cn>1</m:cn>
	      <m:cn>1</m:cn>
	    </m:vector>
	  </m:apply>
	</m:apply>
      </m:math>
    </equation>
    
    <para xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="paraid">Here,
      <m:math>
	<m:vector>
	  <m:cn>1</m:cn>
	  <m:cn>1</m:cn>
	</m:vector>
      </m:math>
      is the eigenvector and <m:math><m:cn>0</m:cn></m:math>
      is its associated eigenvalue.</para>
    
    <equation xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="eq2">
      <m:math mode="block">
	<m:apply>
	  <m:eq/>
	  <m:apply>
	    <m:times/>
	    <m:matrix>
	      <m:matrixrow><m:cn>2</m:cn><m:cn>1</m:cn></m:matrixrow>
	      <m:matrixrow><m:cn>1</m:cn><m:cn>2</m:cn></m:matrixrow>
	    </m:matrix>
	    <m:vector>
	      <m:cn>1</m:cn>
	      <m:cn>1</m:cn>
	    </m:vector>
	  </m:apply>
	  <m:apply>
	    <m:times/>
	    <m:cn>3</m:cn>
	    <m:vector>
	      <m:cn>1</m:cn>
	      <m:cn>1</m:cn>
	    </m:vector>	
	  </m:apply>
	</m:apply>
      </m:math>
    </equation>
    
    <para xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="paraid2">In this second example,
      <m:math>
	<m:vector>
	  <m:cn>1</m:cn>
	  <m:cn>1</m:cn>
	</m:vector>
      </m:math>
is again the eigenvector but the eigenvalue is now <m:math><m:cn>3</m:cn></m:math>.</para>
    
    <para xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="p4">Now we'd like to develop a method of finding the eigenvalues and eigenvectors of a matrix.  We start with what is basically the defining equation behind this whole idea:</para>
    
    <equation xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="eq3">
      <m:math>
	<m:apply>
	  <m:eq/>
	  <m:apply>
	    <m:times/>
	    <m:ci type="matrix">A</m:ci>
	    <m:ci type="vector">x</m:ci>
	  </m:apply>
	  <m:apply>
	    <m:times/>
	    <m:ci>λ</m:ci>
	    <m:ci type="vector">x</m:ci>
	  </m:apply>
	</m:apply>
      </m:math>
    </equation>
    
    <para xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="p5">Next, we move the
      <m:math mode="block">
	<m:apply>
	<m:times/>
	  <m:ci>λ</m:ci>
	  <m:ci type="vector">x</m:ci>
	</m:apply>
      </m:math>
      term to the left-hand side and factor:</para>
    
    <equation xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="eq4">
      <m:math mode="block">
	<m:apply>
	  <m:eq/>
	  <m:apply>
	    <m:times/>
	    <m:apply>
	      <m:minus/>
	      <m:ci type="matrix">A</m:ci>
	      <m:apply>
		<m:times/>
		<m:ci>λ</m:ci>
		<m:ci type="matrix">I</m:ci>
	      </m:apply>
	    </m:apply>
	    <m:ci type="vector">x</m:ci>
	  </m:apply>
	  <m:cn>0</m:cn>
	</m:apply>
      </m:math>
    </equation>
    
    <para xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="p6">Here's the important rule to remember:  there exists
      <m:math>
	<m:apply>
	  <m:neq/>
	  <m:ci type="vector">x</m:ci>
	  <m:cn>0</m:cn>
	</m:apply>
      </m:math>
      satisfying the equation if and only if
      <m:math>
	<m:apply>
	  <m:eq/>
	  <m:apply>
	    <m:determinant/>
	    <m:apply>
	      <m:minus/>
	      <m:ci type="matrix">A</m:ci>
	      <m:apply>
		<m:times/>
		<m:ci>λ</m:ci>
		<m:ci type="matrix">I</m:ci>
	      </m:apply>
	    </m:apply>
	  </m:apply>
	  <m:cn>0</m:cn>
	</m:apply>
      </m:math>.
      So, to find the eigenvalues, we need to solve this determinant equation.  </para>
    
    <example xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="ex1"> 
      
      <para xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="p7">Given the matrix <m:math><m:ci type="matrix">A</m:ci></m:math>, 
	solve for <m:math><m:ci>λ</m:ci></m:math> in 
	<m:math>
	  <m:apply>
	    <m:eq/>
	    <m:apply>
	      <m:determinant/>
	      <m:apply>
		<m:minus/>
		<m:ci type="matrix">A</m:ci>
		<m:apply>
		  <m:times/>
		  <m:ci>λ</m:ci>
		  <m:ci type="matrix">I</m:ci>
		</m:apply>
	      </m:apply>
	    </m:apply>
	    <m:cn>0</m:cn>
	  </m:apply>
	</m:math>.</para>
      
      <equation xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="eq5">
 	<m:math mode="block">
	  <m:apply>
	    <m:eq/>
	    <m:ci type="matrix">A</m:ci>
	    <m:matrix>
	      <m:matrixrow><m:cn>2</m:cn><m:cn>1</m:cn></m:matrixrow>
	      <m:matrixrow><m:cn>1</m:cn><m:cn>2</m:cn></m:matrixrow>
	    </m:matrix>
	  </m:apply>
	</m:math>
      </equation>
      
      <equation xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="eq6">
 	<m:math mode="block">
	  <m:apply>
	    <m:eq/>
	    <m:apply>
	      <m:determinant/>
	      <m:apply>
		<m:minus/>
		<m:ci type="matrix">A</m:ci>
		<m:apply>
		  <m:times/>
		  <m:ci>λ</m:ci>
		  <m:ci type="matrix">I</m:ci>
		</m:apply>
	      </m:apply>
	    </m:apply>
	    <m:apply>
	      <m:determinant/>
	      <m:matrix>
		<m:matrixrow>
		  <m:apply>
		    <m:minus/>
		    <m:cn>2</m:cn>
		    <m:ci>λ</m:ci>
		  </m:apply>
		  <m:cn>1</m:cn>
		</m:matrixrow>
		<m:matrixrow>
		  <m:cn>1</m:cn>
		  <m:apply>
		    <m:minus/>
		    <m:cn>2</m:cn>
		    <m:ci>λ</m:ci>
		  </m:apply>
		</m:matrixrow>
	      </m:matrix>
	    </m:apply>
	    <m:apply>
	      <m:minus/>
	      <m:apply>
		<m:power/>
		<m:apply>
		  <m:minus/>
		  <m:ci>λ</m:ci>
		  <m:cn>2</m:cn>
		</m:apply>
		<m:cn>2</m:cn>
	      </m:apply>
	      <m:cn>1</m:cn>
	    </m:apply>
	    <m:apply>
	      <m:plus/>
	      <m:apply>
		<m:minus/>
		<m:apply>
		  <m:power/>
		  <m:ci>λ</m:ci>
		  <m:cn>2</m:cn>
		</m:apply>
		<m:apply>
		  <m:times/>
		  <m:cn>4</m:cn>
		  <m:ci>λ</m:ci>
		</m:apply>
	      </m:apply>
	      <m:cn>3</m:cn>
	    </m:apply>
	    <m:cn>0</m:cn>
	  </m:apply>
	</m:math>
      </equation>
      
      <equation xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="eq7">
 	<m:math mode="block">
	  <m:apply>
	    <m:eq/>
	    <m:ci>λ</m:ci>
	    <m:set>
	      <m:cn>3</m:cn>
	      <m:cn>1</m:cn>
	    </m:set>
	  </m:apply>
	</m:math>
      </equation>
      
    </example>
    
    <para xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="p8">After finding the eigenvalues, we need to find the
      associated eigenvectors.  Looking at <cnxn xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" strength="9" target="eq4">the defining equation</cnxn>, we see that the
      eigenvector <m:math><m:ci type="vector">x</m:ci></m:math>
      is annihilated by the matrix
      <m:math>
	<m:apply>
	  <m:minus/>
	  <m:ci type="matrix">A</m:ci>
	  <m:apply>
	    <m:times/>
	    <m:ci>λ</m:ci>
	    <m:ci type="matrix">I</m:ci>
	  </m:apply>
	</m:apply>
      </m:math>.
      So to solve for the eigenvectors, we simply find the kernel (nullspace) of
      <m:math>
	<m:apply>
	  <m:minus/>
	  <m:ci type="matrix">A</m:ci>
	  <m:apply>
	    <m:times/>
	    <m:ci>λ</m:ci>
	    <m:ci type="matrix">I</m:ci>
	  </m:apply>
	</m:apply>
      </m:math>
      using the two eigenvalues we just calculated.  If we did this for the example above, we'd find that the eigenvector associated with
      <m:math>
	<m:apply>
	  <m:eq/>
	  <m:ci>λ</m:ci>
	  <m:cn>3</m:cn>
	</m:apply>
      </m:math>
      is
      <m:math>
	<m:vector>
	  <m:cn>1</m:cn><m:cn>1</m:cn>
	</m:vector>
      </m:math>
      and the eigenvector associated with
      <m:math>
	<m:apply>
	  <m:eq/>
	  <m:ci>λ</m:ci>
	  <m:cn>1</m:cn>
	</m:apply>
      </m:math>
      is
      <m:math>
	<m:vector>
	  <m:cn>1</m:cn><m:cn>-1</m:cn>
	</m:vector>
      </m:math>.</para>
    
    <para xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="p9">You may be wondering why eigenvalue decomposition is useful.  It seems at first glance that it is only helpful in determining the effect a matrix has on a certain small subset of possible vectors (the eigenvectors).  However, the benefits become clear when you think about how many other vectors can be looked at from an eigenvalue perspective by decomposing them into components along the available eigenvectors.  For instance, in the above example, let's say we wanted to apply
      <m:math><m:ci type="matrix">A</m:ci></m:math>
      to the vector
      <m:math>
	<m:vector>
	  <m:cn>2</m:cn><m:cn>0</m:cn>
	</m:vector>
	
      </m:math>.
      Instead of doing the matrix multiply (admittedly not too difficult in this case), the vector
      <m:math>
	<m:vector>
	  <m:cn>2</m:cn><m:cn>0</m:cn>
	</m:vector>
	
      </m:math>
      could be split into components in the direction of the eigenvalues: </para>

    <equation xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="eq8">
      <m:math mode="block">
	<m:apply>
	  <m:eq/>
	  <m:vector>
	    <m:cn>2</m:cn>
	    <m:cn>0</m:cn>
	  </m:vector>
	  <m:apply>
	    <m:plus/>
	    <m:vector>
	      <m:cn>1</m:cn>
	      <m:cn>1</m:cn>
	    </m:vector>
	    <m:vector>
	      <m:cn>1</m:cn>
	      <m:cn>-1</m:cn>
	    </m:vector>
	  </m:apply>
	</m:apply>
      </m:math>
    </equation>
    
    <para xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="p10">Now, each of these components could be scaled by the appropriate eigenvalue and then added back together to form the net result.</para>
    
    <section xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="s2">
      <name xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">Multiplicity</name>    
      
      <para xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="p11">Once we have determined the eigenvalues of a particular matrix, we can start to discuss them in terms of their multiplicity.  There are two types of eigenvalue multiplicity: algebraic multiplicity and geometric multiplicity.</para>
      
      <definition xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="d1">
	<term xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">Algebraic Multiplicity</term>
	<meaning xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">The number of repetitions of a certain eigenvalue.  If, for a certain matrix,
	  <m:math>
	    <m:apply>
	      <m:eq/>
	      <m:ci>λ</m:ci>
	      <m:set>
		<m:cn>3</m:cn>
		<m:cn>3</m:cn>
		<m:cn>4</m:cn>
	      </m:set>
	    </m:apply>
	  </m:math>,
	  then the algebraic multiplicity of <m:math><m:cn>3</m:cn></m:math> would be <m:math><m:cn>2</m:cn></m:math>
	  (as it appears twice) and the algebraic multiplicity of <m:math><m:cn>4</m:cn></m:math>
	  would be <m:math><m:cn>1</m:cn></m:math> 
	  (as it appears once).  This type of multiplicity is normally
	  represented by the Greek letter <m:math><m:ci>α</m:ci></m:math>, 
	  where
	  <m:math>
	    <m:apply>
	      <m:ci type="fn">α</m:ci>
	      <m:ci><m:msub><m:mi>λ</m:mi><m:mi>i</m:mi></m:msub></m:ci>
	    </m:apply>
	  </m:math>
represents the algebraic multiplicity of
	  <m:math><m:ci><m:msub><m:mi>λ</m:mi><m:mi>i</m:mi></m:msub></m:ci></m:math>.</meaning>
      </definition>
      
      <definition xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="d2">
	<term xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">Geometric Multiplicity</term>
	<meaning xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">A particular eigenvalue's geometric multiplicity is defined as the dimension of the nullspace of
	  <m:math>
	    <m:apply>
	      <m:minus/>
	      <m:apply>
		<m:times/>
		<m:ci>λ</m:ci>
		<m:ci type="matrix">I</m:ci>
	      </m:apply>
	      <m:ci type="matrix">A</m:ci>
	    </m:apply>
	  </m:math>.
This type of multiplicity is normally represented by the Greek letter
	  <m:math><m:ci>γ</m:ci></m:math>,  where
	  <m:math>
	    <m:apply>
	      <m:ci type="fn">γ</m:ci>
	      <m:ci><m:msub><m:mi>λ</m:mi><m:mi>i</m:mi></m:msub></m:ci>
	    </m:apply>
	  </m:math>
	  represents the geometric multiplicity of
	  <m:math><m:ci><m:msub><m:mi>λ</m:mi><m:mi>i</m:mi></m:msub></m:ci></m:math>.</meaning>
      </definition>
      
    </section>
    
    <section xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="s3">
      <name xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">Helpful Facts</name>
      
      <para xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="p12">Here are some helpful facts about certain special cases of matrices.</para>
      
      <section xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="s3_1">
	<name xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">Rank</name>
	
	<para xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="p13">A matrix
	  <m:math><m:ci type="matrix">A</m:ci></m:math>
	  is full rank if
	  <m:math>
	    <m:apply>
	      <m:neq/>
	      <m:apply>
		<m:determinant/>
		<m:ci type="matrix">A</m:ci>
	      </m:apply>
	      <m:cn>0</m:cn>
	    </m:apply>
	  </m:math>.
	  However, if
	  <m:math>
	    <m:apply>
	      <m:eq/>
	      <m:ci>λ</m:ci>
	      <m:cn>0</m:cn>
	    </m:apply>
	  </m:math>
	  then
	  <m:math>
	    <m:apply>
	      <m:eq/>
	      <m:apply>
		<m:determinant/>
		<m:apply>
		  <m:minus/>
		  <m:apply>
		    <m:times/>
		    <m:ci>λ</m:ci>
		    <m:ci type="matrix">I</m:ci>
		  </m:apply>
		  <m:ci type="matrix">A</m:ci>
		</m:apply>
	      </m:apply>
	      <m:cn>0</m:cn>
	    </m:apply>
	  </m:math>.
	  This tells us that
	  <m:math>
	    <m:apply>
	      <m:eq/>
	      <m:apply>
		<m:determinant/>
		<m:ci type="matrix">A</m:ci>
	      </m:apply>
	      <m:cn>0</m:cn>
	    </m:apply>
	  </m:math>.
	  Therefore, if a matrix has at least one eigenvalue equal to
	  <m:math><m:cn>0</m:cn></m:math>, then it cannot have full rank.  Specifically, for an
	  <m:math><m:ci>n</m:ci></m:math>-dimensional square matrix: </para>
	
	<list xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="l1">
	  <item xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">When one eigenvalue equals <m:math><m:cn>0</m:cn></m:math>, 
	    <m:math>
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:ci type="fn">rank</m:ci>
		  <m:ci type="matrix">A</m:ci>
		</m:apply>
		<m:apply>
		  <m:minus/>
		  <m:ci>n</m:ci>
		  <m:cn>1</m:cn>
		</m:apply>
	      </m:apply>
	    </m:math>
	  </item>
	  <item xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">When multiple eigenvalues equal <m:math><m:cn>0</m:cn></m:math> 
	    <m:math>
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:ci type="fn">rank</m:ci>
		  <m:ci type="matrix">A</m:ci>
		</m:apply>
		<m:apply>
		  <m:minus/>
		  <m:ci>n</m:ci>
		  <m:apply>
		    <m:ci type="function">γ</m:ci>
		    <m:cn>0</m:cn>
		  </m:apply>
		</m:apply>
	      </m:apply>
	    </m:math>.
	    This property holds even if there are other non-zero eigenvalues</item>
	</list>
      </section>
      
      <section xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="s3_2">
	<name xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">Symmetric Matrices</name>

	<para xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="p14">A symmetric matrix is one whose transpose is equal to itself
	  (<m:math>
	    <m:apply>
	      <m:eq/>
	      <m:ci type="matrix">A</m:ci>
	      <m:apply>
		<m:transpose/>
		<m:ci type="matrix">A</m:ci>
	      </m:apply>
	    </m:apply>
	  </m:math>).
	  These matrices (represented by A 
	  below) have the following properties:</para>

	<list xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="l2" type="enumerated">
	  <item xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">Its eigenvalues are real.</item>
	  <item xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">Its eigenvectors are orthogonal.</item>
	  <item xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">They are always diagonalizable.</item>
	</list>

      </section>
      
    </section>
  </content>
</document>
