<?xml version="1.0" encoding="utf-8" standalone="no"?>
<!-- This module has two big mathml problems -->
<!-- Namely, the use of presentation for * -->
<!-- and the mover bar. -->
<!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="m2122">
  
  <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/">Controllability and Observability Grammians</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.7</md:version>
  <md:created xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">2001/04/09</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="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/">controllability</md:keyword>
    <md:keyword xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">observability</md:keyword>
    <md:keyword xmlns:m="http://www.w3.org/1998/Math/MathML" xmlns:bib="http://bibtexml.sf.net/">grammian</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/">
    <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="s1">
      <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/">Controllability Grammian</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="p0">The finite controllability grammian at time
	<m:math><m:apply><m:lt/><m:ci>t</m:ci><m:infinity/></m:apply></m:math>
	is defined as follows.  </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>
	  <m:apply>
	    <m:eq/>
	    <m:apply>
	      <m:ci type="fn">P</m:ci>
	      <m:ci>t</m:ci>
	    </m:apply>
	    <m:apply>
	      <m:int/>
	      <m:bvar><m:ci>τ</m:ci></m:bvar>
	      <m:uplimit><m:ci>t</m:ci></m:uplimit>
	      <m:lowlimit><m:cn>0</m:cn></m:lowlimit>
	      <m:apply>
		<m:times/>
		<m:apply>
		  <m:exp/>
		  <m:apply>
		    <m:times/>
		    <m:ci type="matrix">A</m:ci>
		    <m:ci>τ</m:ci>
		  </m:apply>
		</m:apply>
		<m:ci type="matrix">B</m:ci>
		<m:ci><m:msup><m:mi>B</m:mi><m:mi>*</m:mi></m:msup></m:ci>
		<m:apply>
		  <m:exp/>
		  <m:apply>
		    <m:times/>
		    <m:ci><m:msup><m:mi>A</m:mi><m:mi>*</m:mi></m:msup></m:ci>
		    <m:ci>τ</m:ci>
		  </m:apply>
		</m:apply>
	      </m:apply>
	    </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="p1">This grammian has two important properties.  First,
	<m:math>
	  <m:apply>
	    <m:eq/>
	    <m:apply>
	      <m:ci type="fn">P</m:ci>
	      <m:ci>t</m:ci>
	    </m:apply>

	    <m:apply>
	      <m:geq/>
	      <m:apply>
		<m:ci type="fn"><m:msup><m:mi>P</m:mi><m:mi>*</m:mi></m:msup></m:ci>
		<m:ci>t</m:ci>
	      </m:apply>
	      <m:cn>0</m:cn>
	    </m:apply>
	  </m:apply>
	</m:math>.

	Secondly, the columns of
	<m:math><m:apply><m:ci type="fn">P</m:ci><m:ci>t</m:ci></m:apply></m:math>
	span the controllable space, i.e.
	
	<m:math>
	<m:apply>
	    <m:eq/>
	    <m:apply>
	      <m:ci type="fn">im</m:ci>
	      <m:apply>
		<m:ci type="fn">P</m:ci>
		<m:ci>t</m:ci>
	      </m:apply>
	    </m:apply>
	    <m:apply>
	      <m:ci type="fn">im</m:ci>
	      <m:apply>
		<m:ci type="fn">C</m:ci>
		<m:ci>A</m:ci>
		<m:ci>B</m:ci>
	      </m:apply>
	    </m:apply>
	  </m:apply>
	</m:math>
	
	It can be shown that the state defined by <m:math><m:ci type="matrix">A</m:ci></m:math> and <m:math><m:ci type="matrix">B</m:ci></m:math> 
	is controllable if, and only if,
	<m:math><m:apply><m:ci type="fn">P</m:ci><m:ci>t</m:ci></m:apply></m:math>
	is positive definite for some 
<m:math><m:apply><m:gt/><m:ci>t</m:ci><m:cn>0</m:cn></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="p2">Using the controllability grammian, we can determine how to most efficiently take a system from the zero state to a certain state
	<m:math><m:ci><m:mover accent="true"><m:mi>x</m:mi><m:mo>―</m:mo></m:mover></m:ci></m:math>.
	Given that
	<m:math><m:ci><m:mover accent="true"><m:mi>x</m:mi><m:mo>―</m:mo></m:mover></m:ci></m:math>
	is in the controllable space, there exists <m:math><m:ci>ξ</m:ci></m:math> such that
	
	<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>
	    <m:apply>
	      <m:eq/>
	      <m:ci><m:mover accent="true"><m:mi>x</m:mi><m:mo>―</m:mo></m:mover></m:ci>
	      <m:apply>
		<m:times/>
		<m:apply>
		  <m:ci type="fn">P</m:ci>
		  <m:ci><m:mover accent="true"><m:mi>T</m:mi><m:mo>―</m:mo></m:mover></m:ci>
		</m:apply>
		<m:ci>ξ</m:ci>
	      </m:apply>
	    </m:apply>
	  </m:math>
	</equation>
	
	for some
	<m:math><m:apply><m:gt/><m:ci><m:mover accent="true"><m:mi>T</m:mi><m:mo>―</m:mo></m:mover></m:ci><m:cn>0</m:cn></m:apply></m:math>.   In this case, the minimum energy input required to move the system from zero to
	<m:math><m:ci><m:mover accent="true"><m:mi>x</m:mi><m:mo>―</m:mo></m:mover></m:ci></m:math>
	is
	
	<m:math>
	  <m:apply>
	    <m:eq/>
	    <m:ci><m:mover accent="true"><m:mi>u</m:mi><m:mo>―</m:mo></m:mover></m:ci>
	    <m:apply>
	      <m:times/>
	      <m:ci><m:msup><m:mi>B</m:mi><m:mi>*</m:mi></m:msup></m:ci>
	      <m:apply>
		<m:exp/>
		<m:apply>
		  <m:times/>
		  <m:ci><m:msup><m:mi>A</m:mi><m:mi>*</m:mi></m:msup></m:ci>
		  <m:apply>
		    <m:minus/>
		    <m:ci><m:mover accent="true"><m:mi>T</m:mi><m:mo>―</m:mo></m:mover></m:ci>
		    <m:ci>t</m:ci>
		  </m:apply>
		</m:apply>
	      </m:apply>
	      <m:ci><m:mover accent="true"><m:mi>ξ</m:mi><m:mo>―</m:mo></m:mover></m:ci>
	    </m:apply>
	  </m:apply>
	</m:math>
	
	If the controllability matrix is invertible, we can use <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/" target="eq2" strength="9">the relation equation between ξ and certain state</cnxn> to put
	<m:math><m:ci><m:mover accent="true"><m:mi>u</m:mi><m:mo>―</m:mo></m:mover></m:ci></m:math>
in terms of
	<m:math><m:ci><m:mover accent="true"><m:mi>x</m:mi><m:mo>―</m:mo></m:mover></m:ci></m:math>:
	
	<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:ci>
		<m:mover accent="true">
		  <m:mi>u</m:mi><m:mo>―</m:mo>
		</m:mover>
	      </m:ci>
	      
	      <m:apply>
		<m:times/>
		<m:ci><m:msup><m:mi>B</m:mi><m:mi>*</m:mi></m:msup></m:ci>
		<m:apply>
		  <m:exp/>
		  <m:apply>
		    <m:times/>
		    <m:ci><m:msup><m:mi>A</m:mi><m:mi>*</m:mi></m:msup></m:ci>
		    <m:apply>
		      <m:minus/>
		      <m:ci><m:mover accent="true"><m:mi>T</m:mi><m:mo>―</m:mo></m:mover></m:ci>
		      <m:ci>t</m:ci>
		    </m:apply>
		  </m:apply>
		</m:apply>
		<m:apply>
		  <m:inverse/>
		  <m:apply>
		    <m:ci type="fn">P</m:ci>
		    <m:ci><m:mover accent="true"><m:mi>T</m:mi><m:mo>―</m:mo></m:mover></m:ci>
		  </m:apply>
		</m:apply>
		<m:ci><m:mover accent="true"><m:mi>x</m:mi><m:mo>―</m:mo></m:mover></m:ci>
	      </m:apply>
	    </m:apply>
	  </m:math>	
	</equation>
	
      </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="p3">
In general, this minimal energy is exactly equal to
	<m:math>
	  <m:apply>
	    <m:times/>
	    <m:apply>
	      <m:power/>
	      <m:ci><m:mover accent="true"><m:mi>ξ</m:mi><m:mo>―</m:mo></m:mover></m:ci>
	      <m:ci>*</m:ci>
	    </m:apply>
	    <m:apply>
	      <m:ci type="fn">P</m:ci>
	      <m:ci><m:mover accent="true"><m:mi>T</m:mi><m:mo>―</m:mo></m:mover></m:ci>
	    </m:apply>
	    <m:ci><m:mover accent="true"><m:mi>ξ</m:mi><m:mo>―</m:mo></m:mover></m:ci>
	</m:apply>
	</m:math>.
	If the system is controllable, then this formula becomes
	
	<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_2">
	  <m:math>
	    <m:apply>
	      <m:eq/>
	      <m:apply>
		<m:ci type="fn">Energy</m:ci>
		<m:ci><m:mover accent="true"><m:mi>u</m:mi><m:mo>―</m:mo></m:mover></m:ci>
	      </m:apply>
	      <m:apply>
		<m:times/>
		<m:apply>
		  <m:power/>
		  <m:ci><m:mover><m:mi>x</m:mi><m:mo>―</m:mo></m:mover></m:ci>
		  <m:ci>*</m:ci>
		</m:apply>
		<m:apply>
		  <m:inverse/>
		  <m:apply>
		    <m:ci type="fn">P</m:ci>
		    <m:ci><m:mover><m:mi>T</m:mi><m:mo>―</m:mo></m:mover></m:ci>
		  </m:apply>
		</m:apply>
		<m:ci><m:mover><m:mi>x</m:mi><m:mo>―</m:mo></m:mover></m:ci>
	      </m:apply>		
	    </m:apply>
	  </m:math>
	</equation>
	
	If you don't want to start at the zero state, the formulas above can still be applied for taking a system at state
	<m:math><m:ci><m:msub><m:mi>x</m:mi><m:mn>1</m:mn></m:msub></m:ci></m:math>
	to a state
	<m:math><m:ci><m:msub><m:mi>x</m:mi><m:mn>2</m:mn></m:msub></m:ci></m:math>.
	This holds even if
	<m:math><m:ci><m:msub><m:mi>x</m:mi><m:mn>1</m:mn></m:msub></m:ci></m:math>
	and
	<m:math><m:ci><m:msub><m:mi>x</m:mi><m:mn>2</m:mn></m:msub></m:ci></m:math>
	are not controllable; in this case, all that is necessary is for
	<m:math>
	  <m:apply>
	    <m:minus/>
	    <m:ci><m:msub><m:mi>x</m:mi><m:mn>2</m:mn></m:msub></m:ci>
	    <m:ci><m:msub><m:mi>x</m:mi><m:mn>1</m:mn></m:msub></m:ci>
	  </m:apply>
	</m:math>
	to be in the controllable space.  (This makes sense if you think of
	<m:math><m:ci><m:msub><m:mi>x</m:mi><m:mn>1</m:mn></m:msub></m:ci></m:math>
	as being the zero state and
	<m:math><m:ci><m:msub><m:mi>x</m:mi><m:mn>2</m:mn></m:msub></m:ci></m:math>
	as being the general state we are trying to reach; it is the exact analog of the previous case.  Using
<m:math><m:ci><m:msub><m:mi>x</m:mi><m:mn>1</m:mn></m:msub></m:ci></m:math>
	and
	<m:math><m:ci><m:msub><m:mi>x</m:mi><m:mn>2</m:mn></m:msub></m:ci></m:math>
is just like using <m:math><m:cn>0</m:cn></m:math> and <m:math><m:ci>x</m:ci></m:math> 
	with an appropriate offset.)</para>
      
    </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="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/">Observability Grammian</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="p4">The finite observability grammian at time
	<m:math><m:apply><m:lt/><m:ci>t</m:ci><m:infinity/></m:apply></m:math>
is defined as
	
	<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>
	    <m:apply>
	      <m:eq/>
	      <m:apply>
		<m:ci type="fn">Q</m:ci>
		<m:ci>t</m:ci>
	      </m:apply>
	      <m:apply>
		<m:int/>
		<m:bvar><m:ci>τ</m:ci></m:bvar>
		<m:uplimit><m:ci>t</m:ci></m:uplimit>
		<m:lowlimit><m:cn>0</m:cn></m:lowlimit>
		<m:apply>
		  <m:times/>
		  <m:apply>
		    <m:exp/>
		    <m:apply>
		      <m:times/>
		      <m:ci><m:msup><m:mi>A</m:mi><m:mi>*</m:mi></m:msup></m:ci>
		      <m:ci>τ</m:ci>
		    </m:apply>
		  </m:apply>
		  <m:ci><m:msup><m:mi>C</m:mi><m:mi>*</m:mi></m:msup></m:ci>
		  <m:ci>C</m:ci>
		  <m:apply>
		    <m:exp/>
		    <m:apply>
		      <m:times/>
		      <m:ci type="matrix">A</m:ci>
		      <m:ci>τ</m:ci>
		    </m:apply>
		  </m:apply>
		</m:apply>
	      </m:apply>
	    </m:apply>
	  </m:math>
	</equation>
	
	Parallel to the finite controllability grammian, the kernel of finite observability grammian is equal to the kernel of the observability matrix.  (This relationship holds for positive time only.)
	
	<m:math>
	  <m:apply>
	    <m:eq/>
	    <m:apply>
	      <m:ci type="fn">ker</m:ci>
	      <m:apply>
		<m:ci type="fn">Q</m:ci>
		<m:ci>t</m:ci>
	      </m:apply>
	    </m:apply>
	    <m:apply>
	      <m:ci type="fn">ker</m:ci>
	      <m:apply>
		<m:ci type="fn">O</m:ci>
		<m:ci type="matrix">C</m:ci>
		<m:ci type="matrix">A</m:ci>
	      </m:apply>
	    </m:apply>
	  </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="p5">Using this grammian, we can find an expression for
	the energy of the output <m:math><m:ci>y</m:ci></m:math> at time <m:math><m:ci>T</m:ci></m:math>
	caused by the system's initial state <m:math><m:ci>x</m:ci></m:math>:
	
	<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>
	    <m:apply>
	      <m:eq/>
	      <m:apply>
		<m:ci type="fn">Energy</m:ci>
		<m:ci>y</m:ci>
	      </m:apply>
	      <m:apply>
		<m:times/>
		<m:ci><m:msup><m:mi>x</m:mi><m:mi>*</m:mi></m:msup></m:ci>
		<m:apply>
		  <m:ci type="fn">Q</m:ci>
		  <m:ci>T</m:ci>
		</m:apply>
		<m:ci>x</m:ci>
	      </m:apply>
	    </m:apply>
	  </m:math>
	</equation>

      </para>
      
    </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/">Infinite Grammians</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="p6">Consider a continuous-time linear system defined,
	as per normal, by the matrices <m:math><m:ci type="matrix">A</m:ci></m:math>, <m:math><m:ci type="matrix">B</m:ci></m:math>, <m:math><m:ci type="matrix">C</m:ci></m:math>, and<m:math><m:ci type="matrix">D</m:ci></m:math>.
	Assuming that this system is stable (i.e. all of its eigenvalues have negative real parts), both the controllability and observability grammians are defined for
	<m:math><m:apply><m:eq/><m:ci>t</m:ci><m:infinity/></m:apply></m:math>.
	
	<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>
	    <m:apply>
	      <m:eq/>
	      <m:ci>P</m:ci>
	      <m:apply>
		<m:int/>
		<m:bvar><m:ci>τ</m:ci></m:bvar>
		<m:uplimit><m:infinity/></m:uplimit>
		<m:lowlimit><m:cn>0</m:cn></m:lowlimit>
		<m:apply>
		  <m:times/>
		  <m:apply>
		    <m:exp/>
		    <m:apply>
		      <m:times/>
		      <m:ci type="matrix">A</m:ci>
		      <m:ci>τ</m:ci>
		    </m:apply>
		  </m:apply>
		  <m:ci type="matrix">B</m:ci>
		  <m:ci><m:msup><m:mi>B</m:mi><m:mi>*</m:mi></m:msup></m:ci>
		  <m:apply>
		    <m:exp/>
		    <m:apply>
		      <m:times/>
		      <m:ci><m:msup><m:mi>A</m:mi><m:mi>*</m:mi></m:msup></m:ci>
		      <m:ci>τ</m:ci>
		    </m:apply>
		  </m:apply>
		</m:apply>
	      </m:apply>
	    </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>
	    <m:apply>
	      <m:eq/>
	      <m:ci>Q</m:ci>
	      <m:apply>
		<m:int/>
		<m:bvar><m:ci>τ</m:ci></m:bvar>
		<m:uplimit><m:infinity/></m:uplimit>
		<m:lowlimit><m:cn>0</m:cn></m:lowlimit>
		<m:apply>
		  <m:times/>
		  <m:apply>
		    <m:exp/>
		    <m:apply>
		      <m:times/>
		      <m:ci><m:msup><m:mi>A</m:mi><m:mi>*</m:mi></m:msup></m:ci>
		      <m:ci>τ</m:ci>
		    </m:apply>
		  </m:apply>
		  <m:ci><m:msup><m:mi>C</m:mi><m:mi>*</m:mi></m:msup></m:ci>
		  <m:ci type="matrix">C</m:ci>
		  <m:apply>
		    <m:exp/>
		    <m:apply>
		      <m:times/>
		      <m:ci>A</m:ci>
		      <m:ci>τ</m:ci>
		    </m:apply>
		  </m:apply>
		</m:apply>
	      </m:apply>
	    </m:apply>
	  </m:math>
	</equation>
	
	These are called the infinite controllability and infinite observability grammians, respectively.  These grammians satisfy the linear matrix equations known as the <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/">Lyapunov equations</term>.
	
	<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>
	    <m:apply>
	      <m:eq/>
	      <m:apply>
		<m:plus/>
		<m:apply>
		  <m:times/>
		  <m:ci type="matrix">A</m:ci>
		  <m:ci type="matrix">P</m:ci>
		</m:apply>
		<m:apply>
		  <m:times/>
		  <m:ci type="matrix">P</m:ci>
		  <m:ci><m:msup><m:mi>A</m:mi><m:mi>*</m:mi></m:msup></m:ci>
		</m:apply>
		<m:apply>
		  <m:times/>
		  <m:ci type="matrix">B</m:ci>
		  <m:ci><m:msup><m:mi>B</m:mi><m:mi>*</m:mi></m:msup></m:ci>
		</m:apply>
	      </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="eq9">
	  <m:math>
	    <m:apply>
	      <m:eq/>
	      <m:apply>
		<m:plus/>
		<m:apply>
		  <m:times/>
		  <m:ci><m:msup><m:mi>A</m:mi><m:mi>*</m:mi></m:msup></m:ci>
		  <m:ci type="matrix">Q</m:ci>
		</m:apply>
		<m:apply>
		  <m:times/>
		  <m:ci type="matrix">Q</m:ci>
		  <m:ci type="matrix">A</m:ci>
		</m:apply>		
		<m:apply>
		  <m:times/>
		  <m:ci><m:msup><m:mi>C</m:mi><m:mi>*</m:mi></m:msup></m:ci>
		  <m:ci type="matrix">C</m:ci>
		</m:apply>
	      </m:apply>
	      <m:cn>0</m:cn>
	    </m:apply>
	  </m:math>
	</equation>
	
      </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="p7">In the case of infinite grammians, the equations for minimal energy state transfer and observation energy drop their dependence on time. Assuming stability and complete controllability, the minimal energy required to transfer from zero to state
<m:math><m:ci><m:msub><m:mi>x</m:mi><m:mi>c</m:mi></m:msub></m:ci></m:math>
	is
	
	<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="eq10">
	  <m:math>
	    <m:apply>
	      <m:times/>	
	      <m:apply>
		<m:power/>
		<m:ci><m:msub><m:mi>x</m:mi><m:mi>c</m:mi></m:msub></m:ci>
		<m:ci>*</m:ci>
	      </m:apply>
	      <m:apply>
		<m:inverse/>
		<m:ci>P</m:ci>
	      </m:apply>
	      <m:ci><m:msub><m:mi>x</m:mi><m:mi>c</m:mi></m:msub></m:ci>
	    </m:apply>
	  </m:math>
	</equation>
	
	Similarly, the largest observation energy produced by the state
	<m:math><m:ci><m:msub><m:mi>x</m:mi><m:mi>o</m:mi></m:msub></m:ci></m:math>
	is obtained for an infinite observation interval and is equal to:
	
	<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="eq11">
	  <m:math>
	    <m:apply>
	      <m:times/>	
	      <m:apply>
		<m:power/>
		<m:ci><m:msub><m:mi>x</m:mi><m:mi>o</m:mi></m:msub></m:ci>
		<m:ci>*</m:ci>
	      </m:apply>
	      <m:ci>Q</m:ci>
	      <m:ci><m:msub><m:mi>x</m:mi><m:mi>o</m:mi></m:msub></m:ci>
	    </m:apply>
	  </m:math>
	</equation>
	
      </para>
      
    </section>
    
  </content>
</document>
