<?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="None">
  <name>Solving Linear Constant-Coefficient Difference Equations</name>
  <metadata>
  <md:version>1.3</md:version>
  <md:created>2004/07/20 11:25:27 GMT-5</md:created>
  <md:revised>2007/07/17 21:19:12.978 GMT-5</md:revised>
  <md:authorlist>
      <md:author id="richb">
      <md:firstname>Richard</md:firstname>
      <md:othername>G.</md:othername>
      <md:surname>Baraniuk</md:surname>
      <md:email>richb@rice.edu</md:email>
    </md:author>
  </md:authorlist>

  <md:maintainerlist>
    <md:maintainer id="richb">
      <md:firstname>Richard</md:firstname>
      <md:othername>G.</md:othername>
      <md:surname>Baraniuk</md:surname>
      <md:email>richb@rice.edu</md:email>
    </md:maintainer>
    <md:maintainer id="ahlfing">
      <md:firstname>Robert</md:firstname>
      
      <md:surname>Ahlfinger</md:surname>
      <md:email>ahlfing@rice.edu</md:email>
    </md:maintainer>
  </md:maintainerlist>
  
  <md:keywordlist>
    <md:keyword>Difference Equations</md:keyword>
    <md:keyword>ELEC 301</md:keyword>
  </md:keywordlist>

  <md:abstract>A module concerning the concepts involved in solving linear constant-coefficient difference equations.</md:abstract>
</metadata>

  <content>
    <para id="para1">
      <list id="list1" type="named-list">
	<item><name>Step 1</name>Given the input
	  <m:math>
	    <m:apply>
	      <m:ci type="fn" class="discrete">x</m:ci>
	      <m:ci>n</m:ci>
	    </m:apply>
	  </m:math>, find a solution to
	  <m:math display="block">
	    <m:apply>
	      <m:eq/>
	      <m:apply>
		<m:sum/>
		<m:bvar>
		  <m:ci>k</m:ci>
		</m:bvar>
		<m:uplimit>
		  <m:ci>N</m:ci>
		</m:uplimit>
		<m:lowlimit>
		  <m:cn>0</m:cn>
		</m:lowlimit>
		<m:apply>
		  <m:times/>
		  <m:ci>
		    <m:msub>
		      <m:mi>a</m:mi>
		      <m:mi>k</m:mi>
		    </m:msub>
		  </m:ci>
		  <m:apply>
		    <m:ci type="fn" class="discrete">
		      <m:msub>
			<m:mi>y</m:mi>
			<m:mi>p</m:mi>
		      </m:msub>
		    </m:ci>
		    <m:apply>
		      <m:minus/>
		      <m:ci>n</m:ci>
		      <m:ci>k</m:ci>
		    </m:apply>
		  </m:apply>
		</m:apply>
	      </m:apply>
	      <m:apply>
		<m:apply>
		  <m:sum/>
		  <m:bvar>
		    <m:ci>k</m:ci>
		  </m:bvar>
		  <m:uplimit>
		    <m:ci>M</m:ci>
		  </m:uplimit>
		  <m:lowlimit>
		    <m:cn>0</m:cn>
		  </m:lowlimit>
		  <m:apply>
		    <m:times/>
		    <m:ci>
		      <m:msub>
			<m:mi>b</m:mi>
			<m:mi>k</m:mi>
		      </m:msub>
		    </m:ci>
		    <m:apply>
		      <m:ci type="fn" class="discrete">x</m:ci>
		      <m:apply>
			<m:minus/>
			<m:ci>n</m:ci>
			<m:ci>k</m:ci>
		      </m:apply>
		    </m:apply>
		  </m:apply>
		</m:apply>
	      </m:apply>
	    </m:apply>
	  </m:math>
	  <note>Just any old solution will do!</note>
	  <m:math>
	    <m:apply>
	      <m:ci type="fn" class="discrete">
		<m:msub>
		  <m:mi>y</m:mi>
		  <m:mi>p</m:mi>
		</m:msub>
	      </m:ci>
	      <m:ci>n</m:ci>
	    </m:apply>
	  </m:math> - particular solution.
	</item>
	<item>
	  Solve the homogeneous equation
	  <m:math display="block">
	    <m:apply>
	      <m:eq/>
	      <m:apply>
		<m:sum/>
		<m:bvar>
		  <m:ci>k</m:ci>
		</m:bvar>
		<m:uplimit>
		  <m:ci>N</m:ci>
		</m:uplimit>
		<m:lowlimit>
		  <m:cn>0</m:cn>
		</m:lowlimit>
		<m:apply>
		  <m:times/>
		  <m:ci>
		    <m:msub>
		      <m:mi>a</m:mi>
		      <m:mi>k</m:mi>
		    </m:msub>
		  </m:ci>
		  <m:apply>
		    <m:ci type="fn" class="discrete">
		      <m:msub>
			<m:mi>y</m:mi>
			<m:mi>h</m:mi>
		      </m:msub>
		    </m:ci>
		    <m:apply>
		      <m:minus/>
		      <m:ci>n</m:ci>
		      <m:ci>k</m:ci>
		    </m:apply>
		  </m:apply>
		</m:apply>
	      </m:apply>
	      <m:cn>0</m:cn>
	    </m:apply>
	  </m:math>
	  for
	  <m:math>
	    <m:apply>
	      <m:ci type="fn" class="discrete">
		<m:msub>
		  <m:mi>y</m:mi>
		  <m:mi>h</m:mi>
		</m:msub>
	      </m:ci>
	      <m:ci>n</m:ci>
	    </m:apply>
	  </m:math> - homogeneous solution.
	</item>
	<item>
	  Complete solution given by
	  <m:math display="block">
	    <m:apply>
	      <m:eq/>
	      <m:apply>
		<m:ci type="fn" class="discrete">y</m:ci>
		<m:ci>n</m:ci>
	      </m:apply>
	      <m:apply>
		<m:plus/>
		<m:apply>
		  <m:ci type="fn" class="discrete">
		    <m:msub>
		      <m:mi>y</m:mi>
		      <m:mi>p</m:mi>
		    </m:msub>
		  </m:ci>
		  <m:ci>n</m:ci>
		</m:apply>
		<m:apply>
		  <m:ci type="fn" class="discrete">
		    <m:msub>
		      <m:mi>y</m:mi>
		      <m:mi>h</m:mi>
		    </m:msub>
		  </m:ci>
		  <m:ci>n</m:ci>
		</m:apply>
	      </m:apply>
	    </m:apply>
	  </m:math>
	</item>
      </list>
    </para>   
    <section id="sec1">
      <name>Solving The Homogeneous Equation</name>
      <para id="sec1p1">
	<list id="list2">
	  <item>What does it mean?
	    <figure id="fig1">
	      <media type="application/postscript" src="fig1.eps">
		<media type="image/png" src="fig1.png"/>
	      </media>
	    </figure>
	  </item>
	  <item>
	    Clearly
	    <m:math>
	      <m:apply>
		<m:ci type="fn" class="discrete">
		  <m:msub>
		    <m:mi>y</m:mi>
		    <m:mi>h</m:mi>
		  </m:msub>
		</m:ci>
		<m:ci>n</m:ci>
	      </m:apply>
	    </m:math> depends on the <emphasis>INITIAL
	    CONDITIONS</emphasis> of the system
	    <m:math>
	      <m:ci>T</m:ci>
	    </m:math>.  
	    <list id="list3">
	      <item>Linearity</item> <item>Time-Invariance</item>
	      <item>Causality</item> </list> will each depend on these
	      conditions.
	  </item>
	  <item>
	    In this course, we will emphasize the simplest case, when
	    <m:math>
	      <m:ci>T</m:ci>
	    </m:math>
	    is "initially at rest" with "zero initial conditions."
	    <m:math>
	      <m:apply>
		<m:tendsto/>
	      </m:apply>
	    </m:math> we will get LTI and causal solutions. (although
	    possibly at the expense of <emphasis>stability</emphasis>

	  </item>
	</list>
      </para>
      <example id="ex2">
	<para id="ex2p1">
	  Solve
	  <m:math display="block">
	    <m:apply>
	      <m:eq/>
	      <m:apply>
		<m:minus/>
		<m:apply>
		  <m:ci type="fn" class="discrete">y</m:ci>
		  <m:ci>n</m:ci>
		</m:apply>
		<m:apply>
		  <m:times/>
		  <m:ci>a</m:ci>
		  <m:apply>
		    <m:ci type="fn" class="discrete">y</m:ci>
		    <m:apply>
		      <m:minus/>
		      <m:ci>n</m:ci>
		      <m:cn>1</m:cn>
		    </m:apply>
		  </m:apply>
		</m:apply>
	      </m:apply>
	      <m:apply>
		<m:ci type="fn" class="discrete">x</m:ci>
		<m:ci>n</m:ci>
	      </m:apply>
	    </m:apply>
	  </m:math>
	  where
	  <m:math>
	    <m:apply>
	      <m:lt/>
	      <m:apply>
		<m:abs/>
		<m:ci>a</m:ci>
	      </m:apply>
	      <m:cn>1</m:cn>
	    </m:apply>
	  </m:math>
	  for
	  <m:math>
	    <m:apply>
	      <m:eq/>
	      <m:apply>
		<m:ci type="fn" class="discrete">x</m:ci>
		<m:ci>n</m:ci>
	      </m:apply>
	      <m:apply>
		<m:ci type="fn" class="discrete">δ</m:ci>
		<m:ci>n</m:ci>
	      </m:apply>
	    </m:apply>
	  </m:math>.
	  <list id="list4" type="named-list"><item><name>Step 1: Particular Solution</name>
	      Assume
	      <m:math>
		<m:apply>
		  <m:geq/>
		  <m:ci>n</m:ci>
		  <m:cn>0</m:cn>
		</m:apply>
	      </m:math> and "zero initial conditions"
	      <m:math display="block">
		<m:apply>
		  <m:eq/>
		  <m:apply>
		    <m:ci type="fn" class="discrete">
		      <m:msub>
			<m:mi>y</m:mi>
			<m:mi>p</m:mi>
		      </m:msub>
		    </m:ci>
		    <m:cn>0</m:cn>
		  </m:apply>
		  <m:apply>
		    <m:plus/>
		    <m:apply>
		      <m:apply>
			<m:tendsto/>
			<m:apply>
			  <m:ci type="fn" class="discrete">δ</m:ci>
			  <m:cn>0</m:cn>
			</m:apply>
			<m:cn>1</m:cn>
		      </m:apply>
		    </m:apply>
		    <m:apply>
		      <m:times/>
		      <m:ci>a</m:ci>
		      <m:apply>
			<m:apply>
			  <m:tendsto/>
			  <m:apply>
			    <m:ci type="fn" class="discrete">
			      <m:msub>
				<m:mi>y</m:mi>
				<m:mi>p</m:mi>
			      </m:msub>
			    </m:ci>
			    <m:cn>-1</m:cn>
			  </m:apply>
			  <m:cn>0</m:cn>
			</m:apply>
		      </m:apply>
		    </m:apply>
		  </m:apply>
		  <m:cn>1</m:cn>
		</m:apply>
	      </m:math>
	      <m:math display="block">
		<m:apply>
		  <m:eq/>
		  <m:apply>
		    <m:ci type="fn" class="discrete">
		      <m:msub>
			<m:mi>y</m:mi>
			<m:mi>p</m:mi>
		      </m:msub>
		    </m:ci>
		    <m:cn>1</m:cn>
		  </m:apply>
		  <m:apply>
		    <m:plus/>
		    <m:apply>
		      <m:apply>
			<m:tendsto/>
			<m:apply>
			  <m:ci type="fn" class="discrete">δ</m:ci>
			  <m:cn>1</m:cn>
			</m:apply>
			<m:cn>0</m:cn>
		      </m:apply>
		    </m:apply>
		    <m:apply>
		      <m:times/>
		      <m:ci>a</m:ci>
		      <m:apply>
			<m:apply>
			  <m:tendsto/>
			  <m:apply>
			    <m:ci type="fn" class="discrete">
			      <m:msub>
				<m:mi>y</m:mi>
				<m:mi>p</m:mi>
			      </m:msub>
			    </m:ci>
			    <m:cn>0</m:cn>
			  </m:apply>
			  <m:cn>1</m:cn>
			</m:apply>
		      </m:apply>
		    </m:apply>
		  </m:apply>
		  <m:ci>a</m:ci>
		</m:apply>
	      </m:math>
	      <m:math display="block">
		<m:apply>
		  <m:eq/>
		  <m:apply>
		    <m:ci type="fn" class="discrete">
		      <m:msub>
			<m:mi>y</m:mi>
			<m:mi>p</m:mi>
		      </m:msub>
		    </m:ci>
		    <m:cn>2</m:cn>
		  </m:apply>
		  <m:apply>
		    <m:plus/>
		    <m:apply>
		      <m:apply>
			<m:tendsto/>
			<m:apply>
			  <m:ci type="fn" class="discrete">δ</m:ci>
			  <m:cn>2</m:cn>
			</m:apply>
			<m:cn>0</m:cn>
		      </m:apply>
		    </m:apply>
		    <m:apply>
		      <m:times/>
		      <m:ci>a</m:ci>
		      <m:apply>
			<m:apply>
			  <m:tendsto/>
			  <m:apply>
			    <m:ci type="fn" class="discrete">
			      <m:msub>
				<m:mi>y</m:mi>
				<m:mi>p</m:mi>
			      </m:msub>
			    </m:ci>
			    <m:cn>1</m:cn>
			  </m:apply>
			  <m:ci>a</m:ci>
			</m:apply>
		      </m:apply>
		    </m:apply>
		  </m:apply>
		  <m:apply>
		    <m:power/>
		    <m:ci>a</m:ci>
		    <m:cn>2</m:cn>
		  </m:apply>
		</m:apply>
	      </m:math>
	      <equation id="eq1">
		<m:math>
		  <m:apply>
		    <m:eq/>
		    <m:apply>
		      <m:ci type="fn" class="discrete">
			<m:msub>
			  <m:mi>y</m:mi>
			  <m:mi>p</m:mi>
			</m:msub>
		      </m:ci>
		      <m:ci>n</m:ci>
		    </m:apply>
		    <m:apply>
		      <m:power/>
		      <m:ci>a</m:ci>
		      <m:ci>n</m:ci>
		    </m:apply>
		  </m:apply>
		</m:math>
	      </equation>
	      where
	      <m:math>
		<m:apply>
		  <m:geq/>
		  <m:ci>n</m:ci>
		  <m:cn>0</m:cn>
		</m:apply>
	      </m:math>
	      <figure id="fig2">
		<media type="application/postscript" src="fig2.eps">
		  <media type="image/png" src="fig2.png"/>
		</media>
	      </figure>
	    </item>
	    <item><name>Step 2: Homogeneous Solution</name>
	      If
	      <m:math>
		<m:apply>
		  <m:eq/>
		  <m:apply>
		    <m:ci type="fn" class="discrete">x</m:ci>
		    <m:ci>n</m:ci>
		  </m:apply>
		  <m:cn>0</m:cn>
		</m:apply>
	      </m:math>, then 
	      <m:math>
		<m:apply>
		  <m:eq/>
		  <m:apply>
		    <m:minus/>
		    <m:apply>
		      <m:ci type="fn" class="discrete">
			<m:msub>
			  <m:mi>y</m:mi>
			  <m:mi>h</m:mi>
			</m:msub>
		      </m:ci>
		      <m:ci>n</m:ci>
		    </m:apply>
		    <m:apply>
		      <m:times/>
		      <m:ci>a</m:ci>
		      <m:apply>
			<m:ci type="fn" class="discrete">
			  <m:msub>
			    <m:mi>y</m:mi>
			    <m:mi>h</m:mi>
			  </m:msub>
			</m:ci>
			<m:apply>
			  <m:minus/>
			  <m:ci>n</m:ci>
			  <m:cn>1</m:cn>
			</m:apply>
		      </m:apply>
		    </m:apply>
		  </m:apply>
		  <m:cn>0</m:cn>
		</m:apply>
	      </m:math>
	      <m:math display="block">
		<m:apply>
		  <m:eq/>
		  <m:apply>
		    <m:ci type="fn" class="discrete">
		      <m:msub>
			<m:mi>y</m:mi>
			<m:mi>h</m:mi>
		      </m:msub>
		    </m:ci>
		    <m:ci>n</m:ci>
		  </m:apply>
		  <m:apply>
		    <m:times/>
		    <m:ci>a</m:ci>
		    <m:apply>
		      <m:ci type="fn" class="discrete">
			<m:msub>
			  <m:mi>y</m:mi>
			  <m:mi>h</m:mi>
			</m:msub>
		      </m:ci>
		      <m:apply>
			<m:minus/>
			<m:ci>n</m:ci>
			<m:cn>1</m:cn>
		      </m:apply>
		    </m:apply>
		  </m:apply>		
		</m:apply>		
	      </m:math>
	      A solution is given by
	      <equation id="eq2">
		<m:math>
		  <m:apply>
		    <m:eq/>
		    <m:apply>
		      <m:ci type="fn" class="discrete">
			<m:msub>
			  <m:mi>y</m:mi>
			  <m:mi>h</m:mi>
			</m:msub>
		      </m:ci>
		      <m:ci>n</m:ci>
		    </m:apply>
		    <m:apply>
		      <m:times/>
		      <m:ci>c</m:ci>
		      <m:apply>
			<m:power/>
			<m:ci>a</m:ci>
			<m:ci>n</m:ci>
		      </m:apply>
		    </m:apply>
		  </m:apply>
		</m:math>
	      </equation> for all
	      <m:math>
		<m:ci>n</m:ci>
	      </m:math>.
	      <figure id="fig3">
		<media type="application/postscript" src="fig3.eps">
		  <media type="image/png" src="fig3.png"/>
		</media>
	      </figure>
	    </item>
	    <item><name>Step 3: Reconcile</name>
	      <m:math display="block">
		<m:apply>
		  <m:eq/>
		  <m:apply>
		    <m:ci type="fn" class="discrete">y</m:ci>
		    <m:ci>n</m:ci>
		  </m:apply>
		  <m:apply>
		    <m:plus/>
		    <m:apply>
		      <m:ci type="fn" class="discrete">
			<m:msub>
			  <m:mi>y</m:mi>
			  <m:mi>n</m:mi>
			</m:msub>
		      </m:ci>
		      <m:ci>p</m:ci>
		    </m:apply>
		    <m:apply>
		      <m:ci type="fn" class="discrete">
			<m:msub>
			  <m:mi>y</m:mi>
			  <m:mi>h</m:mi>
			</m:msub>
		      </m:ci>
		      <m:ci>n</m:ci>
		    </m:apply>
		  </m:apply>
		</m:apply>
	      </m:math>
	      <m:math display="block">
		<m:apply>
		  <m:eq/>
		  <m:apply/>
		  <m:apply>
		    <m:plus/>
		    <m:apply>
		      <m:times/>
		      <m:apply>
			<m:power/>
			<m:ci>a</m:ci>
			<m:ci>n</m:ci>
		      </m:apply>
		      <m:apply>
			<m:ci type="fn" class="discrete">u</m:ci>
			<m:ci>n</m:ci>
		      </m:apply>
		    </m:apply>
		    <m:apply>
		      <m:times/>
		      <m:apply>
			<m:ci>c</m:ci>
		      </m:apply>
		      <m:apply>
			<m:power/>
			<m:ci>a</m:ci>
			<m:ci>n</m:ci>
		      </m:apply>
		    </m:apply>
		  </m:apply>
		</m:apply>
	      </m:math>
	      How to pick
	      <m:math>
		<m:ci>c</m:ci>
	      </m:math>?  Need auxilliary conditions.
	      <figure id="fig4">
		<media type="application/postscript" src="fig4.eps">
		  <media type="image/png" src="fig4.png"/>
		</media>
	      </figure>
	      If we desire a <emphasis>causal system</emphasis>, then 
	      <emphasis>
		<m:math>
		  <m:apply>
		    <m:eq/>
		    <m:ci>c</m:ci>
		    <m:cn>0</m:cn>
		  </m:apply>
		</m:math>
	      </emphasis> and
	      <m:math id="math1" display="block">
		<m:apply>
		  <m:eq/>
		  <m:apply>
		    <m:ci type="fn" class="discrete">y</m:ci>
		    <m:ci>n</m:ci>
		  </m:apply>
		  <m:apply>
		    <m:times/>
		    <m:apply>
		      <m:power/>
		      <m:ci>a</m:ci>
		      <m:ci>n</m:ci>
		    </m:apply>
		    <m:apply>
		      <m:ci type="fn" class="discrete">u</m:ci>
		      <m:ci>n</m:ci>
		    </m:apply>
		  </m:apply>
		</m:apply>
	      </m:math>
	      <figure id="fig5">
		<media type="application/postscript" src="fig5.eps">
		  <media type="image/png" src="fig5.png"/>
		</media>
	      </figure>
	      If we desire an <emphasis>anticausal system</emphasis>
	      then choose
	      <m:math>
		<m:apply>
		  <m:eq/>
		  <m:ci>c</m:ci>
		  <m:cn>-1</m:cn>
		</m:apply>
	      </m:math>, so
	      <m:math id="math2" display="block">
		<m:apply>
		  <m:eq/>
		  <m:apply>
		    <m:ci type="fn" class="discrete">y</m:ci>
		    <m:ci>n</m:ci>
		  </m:apply>
		  <m:apply>
		    <m:minus/>
		    <m:apply>
		      <m:times/>
		      <m:apply>
			<m:power/>
			<m:ci>a</m:ci>
			<m:ci>n</m:ci>
		      </m:apply>
		      <m:apply>
			<m:ci type="fn" class="discrete">u</m:ci>
			<m:apply>
			  <m:minus/>
			  <m:ci>n</m:ci>
			</m:apply>
		      </m:apply>
		    </m:apply>
		  </m:apply>
		</m:apply>
	      </m:math>
	      This does not assume "system initially at rest!"
	      <figure id="fig6">
		<media type="application/postscript" src="fig6.eps">
		  <media type="image/png" src="fig6.png"/>
		</media>
	      </figure>
	    </item>
	  </list>
	</para>
      </example>
      <para id="para6">
	<list id="list5" type="enumerated">
	  <name>Notes</name> <item><cnxn target="math1">Solution
	  1</cnxn> was causal and stable.</item> <item><cnxn target="math2">Solution 2</cnxn> was anticausal and
	  unstable.</item>
	</list>
	In general, linearity, time-invariance, and causality of a
	system implemented as a DE will depend on the auxilliary
	conditions.  <note type="Fact">If we assume that the system is
	initially at <term>rest</term> ("zero initial conditions"),
	then it will be <emphasis>LINEAR, TIME-INVARIANT, and
	CAUSAL</emphasis>.
	</note>
      </para>
      <para id="para7">
	<name>Note:</name>
	Setting input
	<m:math>
	  <m:apply>
	    <m:eq/>
	    <m:ci type="vector">x</m:ci>
	    <m:ci>
	      <m:msub>
		<m:mi>δ</m:mi>
		<m:mi>0</m:mi>
	      </m:msub>
	    </m:ci>
	  </m:apply>
	</m:math> impulse and setting initial conditions all 
	<m:math>
	  <m:apply>
	    <m:eq/> <m:apply/> <m:cn>0</m:cn>
	  </m:apply>
	</m:math> and solving for
	<m:math>
	  <m:ci type="vector">
	    <m:msub>
	      <m:mi>y</m:mi>
	      <m:mi>p</m:mi>
	    </m:msub>
	  </m:ci>
	</m:math> yields
	<m:math display="block">
	  <m:apply>
	    <m:eq/>
	    <m:ci type="vector">
	      <m:msub>
		<m:mi>y</m:mi>
		<m:mi>p</m:mi>
	      </m:msub>
	    </m:ci>
	    <m:ci type="vector">h</m:ci>
	  </m:apply>
	</m:math>as the impulse response of this LSI system.
      </para>
      <example id="ex3">
	<name>Frequency Response of a "wire"</name>
	<para id="ex3p1">
	  <figure id="fig7">
	    <media type="application/postscript" src="fig7.eps">
	      <media type="image/png" src="fig7.png"/>
	    </media>
	  </figure>
	  <emphasis>Impulse Response:</emphasis>
	  <figure id="fig8">
	    <media type="application/postscript" src="fig8.eps">
	      <media type="image/png" src="fig8.png"/>
	    </media>
	  </figure>
	  <emphasis>so Frequency Response:</emphasis>
	  <m:math display="block">
	    <m:apply>
	      <m:eq/>
	      <m:ci type="matrix">H</m:ci>
	      <m:apply>
		<m:times/>
		<m:apply>
		  <m:power/>
		  <m:ci>𝔽</m:ci>
		  <m:ci>H</m:ci>
		</m:apply>
		<m:apply>
		  <m:ci class="vector">
		    <m:msub>
		      <m:mi>δ</m:mi>
		      <m:mn>0</m:mn>
		    </m:msub>
		  </m:ci>
		</m:apply>
	      </m:apply>
	      <m:apply>
		<m:times/>
		<m:apply>
		  <m:divide/>
		  <m:cn>1</m:cn>
		  <m:apply>
		    <m:root/>
		    <m:ci>N</m:ci>
		  </m:apply>
		</m:apply>
		<m:apply>
		  <m:sum/>
		  <m:bvar>
		    <m:ci>n</m:ci>
		  </m:bvar>
		  <m:uplimit>
		    <m:apply>
		      <m:minus/>
		      <m:ci>N</m:ci>
		      <m:cn>1</m:cn>
		    </m:apply>
		  </m:uplimit>
		  <m:lowlimit>
		    <m:cn>0</m:cn>
		  </m:lowlimit>
		  <m:apply>
		    <m:times/>
		    <m:apply>
		      <m:ci type="fn" class="discrete">h</m:ci>
		      <m:ci>n</m:ci>
		    </m:apply>
		    <m:apply>
		      <m:exp/>
		      <m:apply>
			<m:minus/>
			<m:apply>
			  <m:times/>
			  <m:imaginary/>
			  <m:apply>
			    <m:divide/>
			    <m:apply>
			      <m:times/>
			      <m:cn>2</m:cn>
			      <m:pi/>
			    </m:apply>  
			    <m:ci>N</m:ci>
			  </m:apply>
			  <m:ci>k</m:ci>
			  <m:ci>n</m:ci>
			</m:apply>
		      </m:apply>
		    </m:apply>
		  </m:apply>
		</m:apply>
	      </m:apply>
	      <m:apply>
		<m:times/>
		<m:apply>
		  <m:divide/>
		  <m:cn>1</m:cn>
		  <m:apply>
		    <m:root/>
		    <m:ci>N</m:ci>
		  </m:apply>
		</m:apply>
		<m:apply>
		  <m:sum/>
		  <m:bvar>
		    <m:ci>n</m:ci>
		  </m:bvar>
		  <m:uplimit>
		    <m:apply>
		      <m:minus/>
		      <m:ci>N</m:ci>
		      <m:cn>1</m:cn>
		    </m:apply>
		  </m:uplimit>
		  <m:lowlimit>
		    <m:cn>0</m:cn>
		  </m:lowlimit>
		  <m:apply>
		    <m:times/>
		    <m:apply>
		      <m:ci type="fn" class="discrete">
			<m:msub>
			  <m:mi>δ</m:mi>
			  <m:mn>0</m:mn>
			</m:msub>
		      </m:ci>
		      <m:ci>n</m:ci>
		    </m:apply>
		    <m:apply>
		      <m:exp/>
		      <m:apply>
			<m:minus/>
			<m:apply>
			  <m:times/>
			  <m:imaginary/>
			  <m:apply>
			    <m:divide/>
			    <m:apply>
			      <m:times/>
			      <m:cn>2</m:cn>
			      <m:pi/>
			    </m:apply>  
			    <m:ci>N</m:ci>
			  </m:apply>
			  <m:ci>k</m:ci>
			  <m:ci>n</m:ci>
			</m:apply>
		      </m:apply>
		    </m:apply>
		  </m:apply>
		</m:apply>
	      </m:apply>
	    </m:apply>
	  </m:math>
	  <m:math>
	    <m:apply>
	      <m:eq/>
	      <m:apply>
		<m:ci type="fn" class="discrete">
		  <m:msub>
		    <m:mi>δ</m:mi>
		    <m:mn>0</m:mn>
		  </m:msub>
		</m:ci>
		<m:ci>n</m:ci>
	      </m:apply>
	      <m:apply>
		<m:piecewise>
		  <m:piece>
		    <m:cn>1</m:cn>
		    <m:apply>
		      <m:eq/>
		      <m:ci>n</m:ci>
		      <m:cn>0</m:cn>
		    </m:apply>
		  </m:piece>
		  <m:otherwise>
		    <m:cn>0</m:cn>
		  </m:otherwise>
		</m:piecewise>
	      </m:apply>
	    </m:apply>
	  </m:math>
	  <m:math display="block">
	    <m:apply>
	      <m:eq/>
	      <m:apply/>
	      <m:apply>
		<m:divide/>
		<m:cn>1</m:cn>
		<m:apply>
		  <m:root/>
		  <m:ci>N</m:ci>
		</m:apply>
	      </m:apply>
	    </m:apply>
	  </m:math>
	  <emphasis>Flat</emphasis>
	  <figure id="fig9">
	    <media type="application/postscript" src="fig9.eps">
	      <media type="image/png" src="fig9.png"/>
	    </media>
	  </figure>
	</para>
      </example>
    </section>
  </content>
</document>
