<?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:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="m10557">
  <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 Matrix Spectral Representation Exercises</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.4</md:version>
  <md:created xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">2002/04/09</md:created>
  <md:revised xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">2002/08/05</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="cox">
      <md:firstname xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">Steven</md:firstname>
      
      <md:surname xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">Cox</md:surname>
      <md:email xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">cox@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="cox">
      <md:firstname xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">Steven</md:firstname>
      
      <md:surname xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">Cox</md:surname>
      <md:email xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">cox@rice.edu</md:email>
    </md:maintainer>
    <md:maintainer xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/" id="jgrab">
      <md:firstname xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">Jacob</md:firstname>
      
      <md:surname xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">Grabczewski</md:surname>
      <md:email xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">jgrab@owlnet.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/">exercise</md:keyword>
    <md:keyword xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">spectral representation</md:keyword>
    <md:keyword xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">symmetric matrix</md:keyword>
  </md:keywordlist>

  <md:abstract xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">A set of exercises after learning about the Spectral Representation, Gram-Schmidt Orthogonalization, and Diagonalization of a Symmetric Matrix.</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="para1">
      The stiffness matrix associated with the <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/" document="m10147" target="simple_swing" strength="9">unstable swing</cnxn> is

      <m:math>
        <m:apply>
          <m:eq/>
	  <m:ci type="matrix">S</m:ci>
	  <m:matrix>
	    <m:matrixrow>
	      <m:cn>1</m:cn>
	      <m:cn>0</m:cn>
	      <m:cn type="integer">-1</m:cn>
	      <m:cn>0</m:cn>
	    </m:matrixrow>
	    <m:matrixrow>
	      <m:cn>0</m:cn>
	      <m:cn>1</m:cn>
	      <m:cn>0</m:cn>
	      <m:cn>0</m:cn>
	    </m:matrixrow>
	    <m:matrixrow>
	      <m:cn type="integer">-1</m:cn>
	      <m:cn>0</m:cn>
	      <m:cn>1</m:cn>
	      <m:cn>0</m:cn>
	    </m:matrixrow>
	    <m:matrixrow>
	      <m:cn>0</m:cn>
	      <m:cn>0</m:cn>
	      <m:cn>0</m:cn>
	      <m:cn>1</m:cn>
	    </m:matrixrow>
	  </m:matrix>
        </m:apply>
      </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="para2">
      <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="para2list1" 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/">
	  Find the three distinct eigenvalues,

	  <m:math>
	    <m:apply>
	      <m:eq/>
              <m:ci><m:msub><m:mi>λ</m:mi><m:mn>1</m:mn></m:msub></m:ci>
              <m:cn>1</m:cn>
	    </m:apply>
	  </m:math>,

	  <m:math>
	    <m:apply>
	      <m:eq/>
              <m:ci><m:msub><m:mi>λ</m:mi><m:mn>2</m:mn></m:msub></m:ci>
              <m:cn>2</m:cn>
	    </m:apply>
	  </m:math>,

	  <m:math>
	    <m:apply>
	      <m:eq/>
              <m:ci><m:msub><m:mi>λ</m:mi><m:mn>3</m:mn></m:msub></m:ci>
              <m:cn>0</m:cn>
	    </m:apply>
	  </m:math>,

	  along with their associated eigenvectors

	  <m:math>
	    <m:ci type="vector">
	      <m:msub>
		<m:mi>e</m:mi>
		<m:mrow>
		  <m:mn>1</m:mn>
		  <m:mo>,</m:mo>
		  <m:mn>1</m:mn>
		</m:mrow>
	      </m:msub>
	    </m:ci>
	  </m:math>,

	  <m:math>
	    <m:ci type="vector">
	      <m:msub>
		<m:mi>e</m:mi>
		<m:mrow>
		  <m:mn>1</m:mn>
		  <m:mo>,</m:mo>
		  <m:mn>2</m:mn>
		</m:mrow>
	      </m:msub>
	    </m:ci>
	  </m:math>,

	  <m:math>
	    <m:ci type="vector">
	      <m:msub>
		<m:mi>e</m:mi>
		<m:mrow>
		  <m:mn>1</m:mn>
		  <m:mo>,</m:mo>
		  <m:mn>3</m:mn>
		</m:mrow>
	      </m:msub>
	    </m:ci>
	  </m:math>,

	  <m:math>
	    <m:ci type="vector">
	      <m:msub>
		<m:mi>e</m:mi>
		<m:mrow>
		  <m:mn>3</m:mn>
		  <m:mo>,</m:mo>
		  <m:mn>1</m:mn>
		</m:mrow>
	      </m:msub>
	    </m:ci>
	  </m:math>, and projection matrices,

	  <m:math>
	    <m:ci type="matrix"><m:msub><m:mi>P</m:mi><m:mn>1</m:mn></m:msub></m:ci>
	  </m:math>,

	  <m:math>
	    <m:ci type="matrix"><m:msub><m:mi>P</m:mi><m:mn>2</m:mn></m:msub></m:ci>
	  </m:math>,

	  <m:math>
	    <m:ci type="matrix"><m:msub><m:mi>P</m:mi><m:mn>3</m:mn></m:msub></m:ci>
	  </m:math>. 

	  What are the respective geometric multiplicities?
	</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/">
	  Use the folk theorem that states "in order to transform a matrix
	  it suffices to transform its eigenvalues" to arrive at a guess for

	  <m:math>
	    <m:apply>
	      <m:power/>
              <m:ci type="matrix">S</m:ci>
              <m:cn type="rational">1<m:sep/>2</m:cn>
	    </m:apply>
	  </m:math>.

	  Show that your guess indeed satisfies

	  <m:math>
	    <m:apply>
	      <m:eq/>
              <m:apply>
                <m:times/>
		<m:apply>
		  <m:power/>
		  <m:ci type="matrix">S</m:ci>
		  <m:cn type="rational">1<m:sep/>2</m:cn>
		</m:apply>
		<m:apply>
		  <m:power/>
		  <m:ci type="matrix">S</m:ci>
		  <m:cn type="rational">1<m:sep/>2</m:cn>
		</m:apply>
              </m:apply>
              <m:ci type="matrix">S</m:ci>
	    </m:apply>
	  </m:math>.

	  Does your 

	  <m:math>
	    <m:apply>
	      <m:power/>
              <m:ci type="matrix">S</m:ci>
              <m:cn type="rational">1<m:sep/>2</m:cn>
	    </m:apply>
	  </m:math>

	  have anything in common with the element-wise square root of 
	  <m:math><m:ci type="matrix">S</m:ci></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/">
	  Show that 

	  <m:math>
	    <m:apply>
	      <m:eq/>
              <m:apply>
                <m:times/>
		<m:reals/>
		<m:ci type="matrix"><m:msub><m:mi>P</m:mi><m:mn>3</m:mn></m:msub></m:ci>
              </m:apply>
              <m:apply>
                <m:times/>
		<m:naturalnumbers/>
		<m:ci type="matrix">S</m:ci>
              </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/">
	  Assemble

	  <m:math display="block">
	    <m:apply>
	      <m:eq/>
              <m:ci type="matrix"><m:msup><m:mi>S</m:mi><m:mo>+</m:mo></m:msup></m:ci>
              <m:apply>
                <m:plus/>
		<m:apply>
		  <m:times/>
		  <m:apply>
		    <m:divide/>
		    <m:cn>1</m:cn>
		    <m:ci><m:msub><m:mi>λ</m:mi><m:mn>1</m:mn></m:msub></m:ci>
		  </m:apply>
		  <m:ci type="matrix"><m:msub><m:mi>P</m:mi><m:mn>1</m:mn></m:msub></m:ci>
		</m:apply>
		<m:apply>
		  <m:times/>
		  <m:apply>
		    <m:divide/>
		    <m:cn>1</m:cn>
		    <m:ci><m:msub><m:mi>λ</m:mi><m:mn>2</m:mn></m:msub></m:ci>
		  </m:apply>
		  <m:ci type="matrix"><m:msub><m:mi>P</m:mi><m:mn>2</m:mn></m:msub></m:ci>
		</m:apply>
              </m:apply>
	    </m:apply>
	  </m:math>

	  and check your result against <code xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">pinv(S)</code> in Matlab.
      </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/">
    Use

    <m:math><m:ci type="matrix"><m:msup><m:mi>S</m:mi><m:mo>+</m:mo></m:msup></m:ci></m:math> to solve

    <m:math>
      <m:apply>
	<m:eq/>
	<m:apply>
	  <m:times/>
	  <m:ci type="matrix">S</m:ci>
	  <m:ci>x</m:ci>
	</m:apply>
	<m:ci type="matrix">f</m:ci>
      </m:apply>
    </m:math> 

    where 

    <m:math>
      <m:apply>
	<m:eq/>
	<m:ci type="matrix">f</m:ci>   
	<m:apply>
	  <m:transpose/>
	  <m:matrix>
	    <m:matrixrow>
	      <m:cn>0</m:cn>
	      <m:cn>1</m:cn>
	      <m:cn>0</m:cn>
	      <m:cn>2</m:cn>
	    </m:matrixrow>
	  </m:matrix>
	</m:apply>
      </m:apply>
    </m:math>

        and carefully draw before and after pictures of the unloaded
        and loaded swing.
      </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/">
        It can be very useful to sketch each of the eigenvectors in
        this fashion. In fact, a movie is the way to go. Please run the Matlab
        truss demo by typing <code xmlns:md="http://cnx.rice.edu/mdml/0.4" xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">truss</code> and view all 12 of the movies.
        Please sketch the 4 eigenvectors of (1) by showing how they
        deform the swing.
      </item>

      </list>

    </para>   
  </content>
  
</document>
