<?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="new0">
  <name>Parameterization of Scaling Coefficients</name>
  <metadata>
  <md:version>2.1</md:version>
  <md:created>2003/05/04</md:created>
  <md:revised>2003/05/05</md:revised>
  <md:authorlist>
      <md:author id="kileen">
      <md:firstname>Kileen</md:firstname>
      
      <md:surname>Cheng</md:surname>
      <md:email>kileen@rice.edu</md:email>
    </md:author>
  </md:authorlist>

  <md:maintainerlist>
    <md:maintainer id="kileen">
      <md:firstname>Kileen</md:firstname>
      
      <md:surname>Cheng</md:surname>
      <md:email>kileen@rice.edu</md:email>
    </md:maintainer>
    <md:maintainer id="jago">
      <md:firstname>Adan</md:firstname>
      
      <md:surname>Galvan</md:surname>
      <md:email>jago@rice.edu</md:email>
    </md:maintainer>
    <md:maintainer id="charlet">
      <md:firstname>Charlet</md:firstname>
      
      <md:surname>Reedstrom</md:surname>
      <md:email>charlet@rice.edu</md:email>
    </md:maintainer>
  </md:maintainerlist>
  
  <md:keywordlist>
    <md:keyword>alpha</md:keyword>
    <md:keyword>beta</md:keyword>
    <md:keyword>Daubechies</md:keyword>
    <md:keyword>degrees of freedom</md:keyword>
    <md:keyword>filter coefficients</md:keyword>
    <md:keyword>Haar</md:keyword>
    <md:keyword>phi</md:keyword>
    <md:keyword>psi</md:keyword>
    <md:keyword>scaling function</md:keyword>
    <md:keyword>wavelet function</md:keyword>
    <md:keyword>wavelet systems</md:keyword>
  </md:keywordlist>

  <md:abstract>Given the conditions our wavelet system needs to satisfy, how do we construct a valid scaling filter?</md:abstract>
</metadata>

  <content>
    <para id="para1">
      The properties of the scaling filter  
      <m:math>
	<m:apply>
	  <m:ci type="fn">h</m:ci>
	  <m:ci>n</m:ci>
	</m:apply>
      </m:math> can be used as criteria in the design of our wavelet
      system.  Given a scaling filter that satisifies our desired
      properties, the scaling and wavelet functions can then be
      calculated.  (see <cnxn document="m11146" strength="8">calculating the scaling function</cnxn>)
    </para>

    <section id="sec1">
      <name>Compact FIR Filters</name>
      <para id="para2">
	A very important class of wavelet systems are those with
	compact support.  These give rise to simple finite impulse
	response (FIR) filters with conventient time-localization
	properties.
      </para>

      <para id="para2b">
	Minimal requirements for these compact FIR filters are:
	<list id="list1" type="enumerated">
	  <item>
	    The length of the scaling filter 
	    <m:math>
	      <m:apply>
		<m:ci type="fn">h</m:ci>
		<m:ci>n</m:ci>
	      </m:apply>
	    </m:math> must be even.
	  </item>
	  <item>
	    <m:math>
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:sum/>
		  <m:bvar><m:ci>n</m:ci></m:bvar>
		  <m:domainofapplication>
		    <m:ci>n</m:ci>
		  </m:domainofapplication>
		  <m:apply>
		    <m:ci type="fn" class="discrete">h</m:ci>
		    <m:ci>n</m:ci>
		  </m:apply>
		</m:apply>
		<m:apply>
		  <m:root/>
		  <m:cn>2</m:cn>
		</m:apply>
	      </m:apply>
	    </m:math>
	  </item>
	  <item>
	    <m:math>
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:sum/>
		  <m:bvar><m:ci>n</m:ci></m:bvar>
		  <m:domainofapplication>
		    <m:ci>n</m:ci>
		  </m:domainofapplication>
		  <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:ci type="fn" class="discrete">h</m:ci>
		      <m:apply>
			<m:minus/>
			<m:ci>n</m:ci>
			<m:apply>
			  <m:times/>
			  <m:cn>2</m:cn>
			  <m:ci>k</m:ci>
			</m:apply>
		      </m:apply>
		    </m:apply>
		  </m:apply>
		</m:apply>
		<m:apply>
		  <m:ci type="fn">δ</m:ci>
		  <m:ci>k</m:ci>
		</m:apply>
	      </m:apply>
	    </m:math>
	  </item>
	</list>

	After the <m:math><m:ci>N</m:ci></m:math> linear constraint of
	(ii) and the 
	<m:math>
	  <m:apply>
	    <m:divide/>
	    <m:ci>N</m:ci>
	    <m:cn>2</m:cn>
	  </m:apply>
	</m:math>
	bilinear constraints of (iii), there are 
	<m:math>
	  <m:apply>
	    <m:minus/>
	    <m:apply>
	      <m:divide/>
	      <m:ci>N</m:ci>
	      <m:cn>2</m:cn>
	    </m:apply>
	    <m:cn>1</m:cn>
	  </m:apply>
	</m:math> remaining degrees of freedom that we can adjust to
	give a valid 
	<m:math>
	  <m:apply>
	    <m:ci type="fn">h</m:ci>
	    <m:ci>n</m:ci>
	  </m:apply>
	</m:math>.
      </para>

      <section id="sec1.1">
	<name>Length-2 Scaling Filter</name>
	<para id="para3">
	  Requirements:

	  <list id="list2" type="enumerated">
	    <item>
	      <m:math>
		<m:apply>
		  <m:eq/>
		  <m:ci>N</m:ci>
		  <m:cn>2</m:cn>
		</m:apply>
	      </m:math>
	    </item>
	    <item>
	      <m:math>
		<m:apply>
		  <m:eq/>
		  <m:apply>
		    <m:plus/>
		    <m:apply>
		      <m:ci type="fn">h</m:ci>
		      <m:cn>0</m:cn>
		    </m:apply>
		    <m:apply>
		      <m:ci type="fn">h</m:ci>
		      <m:cn>1</m:cn>
		    </m:apply>
		  </m:apply>
		  <m:apply>
		    <m:root/>
		    <m:cn>2</m:cn>
		  </m:apply>
		</m:apply>
	      </m:math>
	    </item>
	    <item>
	      <m:math>
		<m:apply>
		  <m:eq/>
		  <m:apply>
		    <m:plus/>
		    <m:apply>
		      <m:power/>
		      <m:apply>
			<m:ci type="fn">h</m:ci>
			<m:cn>0</m:cn>
		      </m:apply>
		      <m:cn>2</m:cn>
		    </m:apply>
		    <m:apply>
		      <m:power/>
		      <m:apply>
			<m:ci type="fn">h</m:ci>
			<m:cn>1</m:cn>
		      </m:apply>
		      <m:cn>2</m:cn>
		    </m:apply>
		  </m:apply>
		  <m:cn>1</m:cn>
		</m:apply>
	      </m:math>
	    </item>
	  </list>

	  <note>
	    Degrees of freedom: 
	    <m:math>
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:minus/>
		  <m:apply>
		    <m:divide/>
		    <m:ci>N</m:ci>
		    <m:cn>2</m:cn>
		  </m:apply>
		  <m:cn>1</m:cn>
		</m:apply>
		<m:apply>
		  <m:factorial/>
		  <m:cn>0</m:cn>
		</m:apply>
	      </m:apply>
	    </m:math>
	  </note>
	  After the minimal requirements for the scaling filter are
	  satisfied, there are no degrees of freedom left to give us
	  flexibility in the design of 
	  <m:math>
	    <m:apply>
	      <m:ci type="fn">h</m:ci>
	      <m:ci>n</m:ci>
	    </m:apply>
	  </m:math>.  There is only one set of possible coefficients:

	  <equation id="eqn1">
	    <m:math>
	      <m:apply>
		<m:eq/>
		<m:ci><m:msub>
		    <m:mi>h</m:mi>
		    <m:mi>D2</m:mi>
		  </m:msub></m:ci>
		<m:list>
		  <m:apply>
		    <m:ci type="fn">h</m:ci>
		    <m:cn>0</m:cn>
		  </m:apply>
		  <m:apply>
		    <m:ci type="fn">h</m:ci>
		    <m:cn>1</m:cn>
		  </m:apply>
		</m:list>
		<m:list>
		  <m:apply>
		    <m:divide/>
		    <m:cn>1</m:cn>
		    <m:apply>
		      <m:root/>
		      <m:cn>2</m:cn>
		    </m:apply>
		  </m:apply>
		  <m:apply>
		    <m:divide/>
		    <m:cn>1</m:cn>
		    <m:apply>
		      <m:root/>
		      <m:cn>2</m:cn>
		    </m:apply>
		  </m:apply>
		</m:list>
	      </m:apply>
	    </m:math>
	  </equation>

	  The length-2 scaling coefficient vector is also known as the
	  Haar or Daubechies-2 coefficients and will be discussed in further
	  length <cnxn document="m11150" strength="8">later</cnxn>.
	</para>
      </section>

      <section id="sec1.2">
	<name>Length-4 Scaling Filter</name>
	<para id="para4">
	  Requirements:

	  <list id="list3" type="enumerated">
	    <item>
	      <m:math>
		<m:apply>
		  <m:eq/>
		  <m:ci>N</m:ci>
		  <m:cn>4</m:cn>
		</m:apply>
	      </m:math>
	    </item>
	    <item>
	      <m:math>
		<m:apply>
		  <m:eq/>
		  <m:apply>
		    <m:plus/>
		    <m:apply>
		      <m:ci type="fn">h</m:ci>
		      <m:cn>0</m:cn>
		    </m:apply>
		    <m:apply>
		      <m:ci type="fn">h</m:ci>
		      <m:cn>1</m:cn>
		    </m:apply>
		    <m:apply>
		      <m:ci type="fn">h</m:ci>
		      <m:cn>2</m:cn>
		    </m:apply>
		    <m:apply>
		      <m:ci type="fn">h</m:ci>
		      <m:cn>3</m:cn>
		    </m:apply>
		  </m:apply>
		  <m:apply>
		    <m:root/>
		    <m:cn>2</m:cn>
		  </m:apply>
		</m:apply>
	      </m:math>
	    </item>
	    <item>
	      <m:math>
		<m:apply>
		  <m:eq/>
		  <m:apply>
		    <m:plus/>
		    <m:apply>
		      <m:power/>
		      <m:apply>
			<m:ci type="fn">h</m:ci>
			<m:cn>0</m:cn>
		      </m:apply>
		      <m:cn>2</m:cn>
		    </m:apply>
		    <m:apply>
		      <m:power/>
		      <m:apply>
			<m:ci type="fn">h</m:ci>
			<m:cn>1</m:cn>
		      </m:apply>
		      <m:cn>2</m:cn>
		    </m:apply>
		    <m:apply>
		      <m:power/>
		      <m:apply>
			<m:ci type="fn">h</m:ci>
			<m:cn>3</m:cn>
		      </m:apply>
		      <m:cn>2</m:cn>
		    </m:apply>
		  </m:apply>
		  <m:cn>1</m:cn>
		</m:apply>
	      </m:math>
	      and 
	      <m:math>
		<m:apply>
		  <m:eq/>
		  <m:apply>
		    <m:plus/>
		    <m:apply>
		      <m:times/>
		      <m:apply>
			<m:ci type="fn">h</m:ci>
			<m:cn>0</m:cn>
		      </m:apply>
		      <m:apply>
			<m:ci type="fn">h</m:ci>
			<m:cn>2</m:cn>
		      </m:apply>
		    </m:apply>
		    <m:apply>
		      <m:times/>
		      <m:apply>
			<m:ci type="fn">h</m:ci>
			<m:cn>1</m:cn>
		      </m:apply>
		      <m:apply>
			<m:ci type="fn">h</m:ci>
			<m:cn>3</m:cn>
		      </m:apply>
		    </m:apply>
		  </m:apply>
		  <m:cn>0</m:cn>
		</m:apply>
	      </m:math>
	    </item>
	  </list>

	  <note>
	    Degrees of freedom: 
	    <m:math>
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:minus/>
		  <m:apply>
		    <m:divide/>
		    <m:ci>N</m:ci>
		    <m:cn>2</m:cn>
		  </m:apply>
		  <m:cn>1</m:cn>
		</m:apply>
		<m:cn>1</m:cn>
	      </m:apply>
	    </m:math>
	  </note>
	  With a length-4 scaling vector, there is still one degree of
	  freedom remaining after the minimal requirements have been
	  satisfied.  Letting <m:math><m:ci>α</m:ci></m:math>
	  represent this degree of freedom parameter, we can formulate
	  scaling filter coefficient equations such that:

	  <equation id="eqn2">
	    <m:math>
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:ci type="fn">h</m:ci>
		  <m:cn>0</m:cn>
		</m:apply>
		<m:apply>
		  <m:divide/>
		  <m:apply>
		    <m:plus/>
		    <m:apply>
		      <m:minus/>
		      <m:cn>1</m:cn>
		      <m:apply>
			<m:cos/>
			<m:ci>α</m:ci>
		      </m:apply>
		    </m:apply>
		    <m:apply>
		      <m:sin/>
		      <m:ci>α</m:ci>
		    </m:apply>
		  </m:apply>
		  <m:apply>
		    <m:times/>
		    <m:cn>2</m:cn>
		    <m:apply>
		      <m:root/>
		      <m:cn>2</m:cn>
		    </m:apply>
		  </m:apply>
		</m:apply>
	      </m:apply>
	    </m:math>
	  </equation>

	  <equation id="eqn3">
	    <m:math>
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:ci type="fn">h</m:ci>
		  <m:cn>1</m:cn>
		</m:apply>
		<m:apply>
		  <m:divide/>
		  <m:apply>
		    <m:plus/>
		    <m:cn>1</m:cn>
		    <m:apply>
		      <m:cos/>
		      <m:ci>α</m:ci>
		    </m:apply>
		    <m:apply>
		      <m:sin/>
		      <m:ci>α</m:ci>
		    </m:apply>
		  </m:apply>
		  <m:apply>
		    <m:times/>
		    <m:cn>2</m:cn>
		    <m:apply>
		      <m:root/>
		      <m:cn>2</m:cn>
		    </m:apply>
		  </m:apply>
		</m:apply>
	      </m:apply>
	    </m:math>
	  </equation>

	  <equation id="eqn4">
	    <m:math>
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:ci type="fn">h</m:ci>
		  <m:cn>2</m:cn>
		</m:apply>
		<m:apply>
		  <m:divide/>
		  <m:apply>
		    <m:minus/>
		    <m:apply>
		      <m:plus/>
		      <m:cn>1</m:cn>
		      <m:apply>
			<m:cos/>
			<m:ci>α</m:ci>
		      </m:apply>
		    </m:apply>
		    <m:apply>
		      <m:sin/>
		      <m:ci>α</m:ci>
		    </m:apply>
		  </m:apply>
		  <m:apply>
		    <m:times/>
		    <m:cn>2</m:cn>
		    <m:apply>
		      <m:root/>
		      <m:cn>2</m:cn>
		    </m:apply>
		  </m:apply>
		</m:apply>
	      </m:apply>
	    </m:math>
	  </equation>

	  <equation id="eqn5">
	    <m:math>
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:ci type="fn">h</m:ci>
		  <m:cn>3</m:cn>
		</m:apply>
		<m:apply>
		  <m:divide/>
		  <m:apply>
		    <m:minus/>
		    <m:apply>
		      <m:minus/>
		      <m:cn>1</m:cn>
		      <m:apply>
			<m:cos/>
			<m:ci>α</m:ci>
		      </m:apply>
		    </m:apply>
		    <m:apply>
		      <m:sin/>
		      <m:ci>α</m:ci>
		    </m:apply>
		  </m:apply>
		  <m:apply>
		    <m:times/>
		    <m:cn>2</m:cn>
		    <m:apply>
		      <m:root/>
		      <m:cn>2</m:cn>
		    </m:apply>
		  </m:apply>
		</m:apply>
	      </m:apply>
	    </m:math>
	  </equation>

	  We can adjust <m:math><m:ci>α</m:ci></m:math> to give
	  us a wavelet system with the desired properties.  However,
	  most values of <m:math><m:ci>α</m:ci></m:math> do not
	  lead to a useful wavelet.  The Daubechies wavelet with
	  filter length 4 arises from 
	  <m:math>
	    <m:apply>
	      <m:eq/>
	      <m:ci>α</m:ci>
	      <m:apply>
		<m:divide/>
		<m:pi/>
		<m:cn>3</m:cn>
	      </m:apply>
	    </m:apply>
	  </m:math>.

	  <equation id="eqn6">
	    <m:math>
	      <m:apply>
		<m:eq/>
		<m:ci><m:msub>
		    <m:mi>h</m:mi>
		    <m:mi>D4</m:mi>
		  </m:msub></m:ci>
		<m:list>
		  <m:apply>
		    <m:divide/>
		    <m:apply>
		      <m:plus/>
		      <m:cn>1</m:cn>
		      <m:apply>
			<m:root/>
			<m:cn>3</m:cn>
		      </m:apply>
		    </m:apply>
		    <m:apply>
		      <m:times/>
		      <m:cn>4</m:cn>
		      <m:apply>
			<m:root/>
			<m:cn>2</m:cn>
		      </m:apply>
		    </m:apply>
		  </m:apply>
		  <m:apply>
		    <m:divide/>
		    <m:apply>
		      <m:plus/>
		      <m:cn>3</m:cn>
		      <m:apply>
			<m:root/>
			<m:cn>3</m:cn>
		      </m:apply>
		    </m:apply>
		    <m:apply>
		      <m:times/>
		      <m:cn>4</m:cn>
		      <m:apply>
			<m:root/>
			<m:cn>2</m:cn>
		      </m:apply>
		    </m:apply>
		  </m:apply>
		  <m:apply>
		    <m:divide/>
		    <m:apply>
		      <m:minus/>
		      <m:cn>3</m:cn>
		      <m:apply>
			<m:root/>
			<m:cn>3</m:cn>
		      </m:apply>
		    </m:apply>
		    <m:apply>
		      <m:times/>
		      <m:cn>4</m:cn>
		      <m:apply>
			<m:root/>
			<m:cn>2</m:cn>
		      </m:apply>
		    </m:apply>
		  </m:apply>
		  <m:apply>
		    <m:divide/>
		    <m:apply>
		      <m:minus/>
		      <m:cn>1</m:cn>
		      <m:apply>
			<m:root/>
			<m:cn>3</m:cn>
		      </m:apply>
		    </m:apply>
		    <m:apply>
		      <m:times/>
		      <m:cn>4</m:cn>
		      <m:apply>
			<m:root/>
			<m:cn>2</m:cn>
		      </m:apply>
		    </m:apply>
		  </m:apply>
		</m:list>
	      </m:apply>
	    </m:math>
	  </equation>

	  Note that for 
	  <m:math>
	    <m:apply>
	      <m:in/>
	      <m:ci>α</m:ci>
	      <m:set>
		<m:cn>0</m:cn>
		<m:apply>
		  <m:divide/>
		  <m:pi/>
		  <m:cn>2</m:cn>
		</m:apply>
		<m:apply>
		  <m:divide/>
		  <m:apply>
		    <m:times/>
		    <m:cn>3</m:cn>
		    <m:pi/>
		  </m:apply>
		  <m:cn>2</m:cn>
		</m:apply>
		<m:pi/>
	      </m:set>
	    </m:apply>
	  </m:math>, we get the length-2 Haar coefficients.
	</para>
      </section>

      <section id="sec1.3">
	<name>Length-6 Scaling Filter</name>
	<para id="para5">
	  Requirements:

	  <list id="list4" type="enumerated">
	    <item>
	      <m:math>
		<m:apply>
		  <m:eq/>
		  <m:ci>N</m:ci>
		  <m:cn>6</m:cn>
		</m:apply>
	      </m:math>
	    </item>
	    <item>
	      <m:math>
		<m:apply>
		  <m:eq/>
		  <m:apply>
		    <m:plus/>
		    <m:apply>
		      <m:ci type="fn">h</m:ci>
		      <m:cn>0</m:cn>
		    </m:apply>
		    <m:apply>
		      <m:ci type="fn">h</m:ci>
		      <m:cn>1</m:cn>
		    </m:apply>
		    <m:apply>
		      <m:ci type="fn">h</m:ci>
		      <m:cn>2</m:cn>
		    </m:apply>
		    <m:apply>
		      <m:ci type="fn">h</m:ci>
		      <m:cn>3</m:cn>
		    </m:apply>
		    <m:apply>
		      <m:ci type="fn">h</m:ci>
		      <m:cn>4</m:cn>
		    </m:apply>
		    <m:apply>
		      <m:ci type="fn">h</m:ci>
		      <m:cn>5</m:cn>
		    </m:apply>
		  </m:apply>
		  <m:apply>
		    <m:root/>
		    <m:cn>2</m:cn>
		  </m:apply>
		</m:apply>
	      </m:math>
	    </item>
	    <item>
	      <m:math>
		<m:apply>
		  <m:eq/>
		  <m:apply>
		    <m:plus/>
		    <m:apply>
		      <m:power/>
		      <m:apply>
			<m:ci type="fn">h</m:ci>
			<m:cn>0</m:cn>
		      </m:apply>
		      <m:cn>2</m:cn>
		    </m:apply>
		    <m:apply>
		      <m:power/>
		      <m:apply>
			<m:ci type="fn">h</m:ci>
			<m:cn>1</m:cn>
		      </m:apply>
		      <m:cn>2</m:cn>
		    </m:apply>
		    <m:apply>
		      <m:power/>
		      <m:apply>
			<m:ci type="fn">h</m:ci>
			<m:cn>3</m:cn>
		      </m:apply>
		      <m:cn>2</m:cn>
		    </m:apply>
		    <m:apply>
		      <m:power/>
		      <m:apply>
			<m:ci type="fn">h</m:ci>
			<m:cn>4</m:cn>
		      </m:apply>
		      <m:cn>2</m:cn>
		    </m:apply>
		    <m:apply>
		      <m:power/>
		      <m:apply>
			<m:ci type="fn">h</m:ci>
			<m:cn>5</m:cn>
		      </m:apply>
		      <m:cn>2</m:cn>
		    </m:apply>
		  </m:apply>
		  <m:cn>1</m:cn>
		</m:apply>
	      </m:math>
	      and 
	      <m:math>
		<m:apply>
		  <m:eq/>
		  <m:apply>
		    <m:plus/>
		    <m:apply>
		      <m:times/>
		      <m:apply>
			<m:ci type="fn">h</m:ci>
			<m:cn>0</m:cn>
		      </m:apply>
		      <m:apply>
			<m:ci type="fn">h</m:ci>
			<m:cn>2</m:cn>
		      </m:apply>
		    </m:apply>
		    <m:apply>
		      <m:times/>
		      <m:apply>
			<m:ci type="fn">h</m:ci>
			<m:cn>1</m:cn>
		      </m:apply>
		      <m:apply>
			<m:ci type="fn">h</m:ci>
			<m:cn>3</m:cn>
		      </m:apply>
		    </m:apply>
		    <m:apply>
		      <m:times/>
		      <m:apply>
			<m:ci type="fn">h</m:ci>
			<m:cn>2</m:cn>
		      </m:apply>
		      <m:apply>
			<m:ci type="fn">h</m:ci>
			<m:cn>4</m:cn>
		      </m:apply>
		    </m:apply>
		    <m:apply>
		      <m:times/>
		      <m:apply>
			<m:ci type="fn">h</m:ci>
			<m:cn>3</m:cn>
		      </m:apply>
		      <m:apply>
			<m:ci type="fn">h</m:ci>
			<m:cn>5</m:cn>
		      </m:apply>
		    </m:apply>
		  </m:apply>
		  <m:cn>0</m:cn>
		</m:apply>
	      </m:math>
	    </item>
	  </list>

	  <note>
	    Degrees of freedom: 
	    <m:math>
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:minus/>
		  <m:apply>
		    <m:divide/>
		    <m:ci>N</m:ci>
		    <m:cn>2</m:cn>
		  </m:apply>
		  <m:cn>1</m:cn>
		</m:apply>
		<m:cn>2</m:cn>
	      </m:apply>
	    </m:math>
	  </note>

	  We now have two degrees of freedom.  Defining our freedom
	  parameters as <m:math><m:ci>α</m:ci></m:math> and
	  <m:math><m:ci>β</m:ci></m:math>, our resulting
	  coefficient vector becomes:

	  <equation id="eqn7">
	    <m:math>
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:ci type="fn">h</m:ci>
		  <m:cn>0</m:cn>
		</m:apply>
		<m:apply>
		  <m:divide/>
		  <m:apply>
		    <m:plus/>
		    <m:apply>
		      <m:times/>
		      <m:apply>
			<m:plus/>
			<m:cn>1</m:cn>
			<m:apply>
			  <m:cos/>
			  <m:ci>α</m:ci>
			</m:apply>
			<m:apply>
			  <m:sin/>
			  <m:ci>α</m:ci>
			</m:apply>
		      </m:apply>
		      <m:apply>
			<m:minus/>
			<m:apply>
			  <m:minus/>
			  <m:cn>1</m:cn>
			  <m:apply>
			    <m:cos/>
			    <m:ci>β</m:ci>
			  </m:apply>
			</m:apply>
			<m:apply>
			  <m:sin/>
			  <m:ci>β</m:ci>
			</m:apply>
		      </m:apply>
		    </m:apply>
		    <m:apply>
		      <m:times/>
		      <m:cn>2</m:cn>
		      <m:apply>
			<m:sin/>
			<m:ci>β</m:ci>
		      </m:apply>
		      <m:apply>
			<m:cos/>
			<m:ci>α</m:ci>
		      </m:apply>
		    </m:apply>
		  </m:apply>
		  <m:apply>
		    <m:times/>
		    <m:cn>4</m:cn>
		    <m:apply>
		      <m:root/>
		      <m:cn>2</m:cn>
		    </m:apply>
		  </m:apply>
		</m:apply>
	      </m:apply>
	    </m:math>
	  </equation>

	  <equation id="eqn8">
	    <m:math>
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:ci type="fn">h</m:ci>
		  <m:cn>1</m:cn>
		</m:apply>
		<m:apply>
		  <m:divide/>
		  <m:apply>
		    <m:plus/>
		    <m:apply>
		      <m:times/>
		      <m:apply>
			<m:plus/>
			<m:cn>1</m:cn>
			<m:apply>
			  <m:cos/>
			  <m:ci>α</m:ci>
			</m:apply>
			<m:apply>
			  <m:sin/>
			  <m:ci>α</m:ci>
			</m:apply>
		      </m:apply>
		      <m:apply>
			<m:minus/>
			<m:apply>
			  <m:plus/>
			  <m:cn>1</m:cn>
			  <m:apply>
			    <m:cos/>
			    <m:ci>β</m:ci>
			  </m:apply>
			</m:apply>
			<m:apply>
			  <m:sin/>
			  <m:ci>β</m:ci>
			</m:apply>
		      </m:apply>
		    </m:apply>
		    <m:apply>
		      <m:times/>
		      <m:cn>2</m:cn>
		      <m:apply>
			<m:sin/>
			<m:ci>β</m:ci>
		      </m:apply>
		      <m:apply>
			<m:cos/>
			<m:ci>α</m:ci>
		      </m:apply>
		    </m:apply>
		  </m:apply>
		  <m:apply>
		    <m:times/>
		    <m:cn>4</m:cn>
		    <m:apply>
		      <m:root/>
		      <m:cn>2</m:cn>
		    </m:apply>
		  </m:apply>
		</m:apply>
	      </m:apply>
	    </m:math>
	  </equation>

	  <equation id="eqn9">
	    <m:math>
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:ci type="fn">h</m:ci>
		  <m:cn>2</m:cn>
		</m:apply>
		<m:apply>
		  <m:divide/>
		  <m:apply>
		    <m:minus/>
		    <m:apply>
		      <m:minus/>
		      <m:apply>
			<m:plus/>
			<m:cn>1</m:cn>
			<m:apply>
			  <m:cos/>
			  <m:ci>α</m:ci>
			</m:apply>
		      </m:apply>
		      <m:apply>
			<m:plus/>
			<m:ci>β</m:ci>
			<m:apply>
			  <m:sin/>
			  <m:ci>α</m:ci>
			</m:apply>
		      </m:apply>
		    </m:apply>
		    <m:ci>β</m:ci>
		  </m:apply>
		  <m:apply>
		    <m:times/>
		    <m:cn>2</m:cn>
		    <m:apply>
		      <m:root/>
		      <m:cn>2</m:cn>
		    </m:apply>
		  </m:apply>
		</m:apply>
	      </m:apply>
	    </m:math>
	  </equation>

	  <equation id="eqn10">
	    <m:math>
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:ci type="fn">h</m:ci>
		  <m:cn>3</m:cn>
		</m:apply>
		<m:apply>
		  <m:divide/>
		  <m:apply>
		    <m:minus/>
		    <m:apply>
		      <m:minus/>
		      <m:apply>
			<m:plus/>
			<m:cn>1</m:cn>
			<m:apply>
			  <m:cos/>
			  <m:ci>α</m:ci>
			</m:apply>
		      </m:apply>
		      <m:apply>
			<m:minus/>
			<m:ci>β</m:ci>
			<m:apply>
			  <m:sin/>
			  <m:ci>α</m:ci>
			</m:apply>
		      </m:apply>
		    </m:apply>
		    <m:ci>β</m:ci>
		  </m:apply>
		  <m:apply>
		    <m:times/>
		    <m:cn>2</m:cn>
		    <m:apply>
		      <m:root/>
		      <m:cn>2</m:cn>
		    </m:apply>
		  </m:apply>
		</m:apply>
	      </m:apply>
	    </m:math>
	  </equation>

	  <equation id="eqn11">
	    <m:math>
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:ci type="fn">h</m:ci>
		  <m:cn>4</m:cn>
		</m:apply>
		<m:apply>
		  <m:minus/>
		  <m:apply>
		    <m:minus/>
		    <m:apply>
		      <m:divide/>
		      <m:cn>1</m:cn>
		      <m:apply>
			<m:root/>
			<m:cn>2</m:cn>
		      </m:apply>
		    </m:apply>
		    <m:apply>
		      <m:ci type="fn">h</m:ci>
		      <m:cn>0</m:cn>
		    </m:apply>
		  </m:apply>
		  <m:apply>
		    <m:ci type="fn">h</m:ci>
		    <m:cn>2</m:cn>
		  </m:apply>
		</m:apply>
	      </m:apply>
	    </m:math>
	  </equation>

	  <equation id="eqn12">
	    <m:math>
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:ci type="fn">h</m:ci>
		  <m:cn>5</m:cn>
		</m:apply>
		<m:apply>
		  <m:minus/>
		  <m:apply>
		    <m:minus/>
		    <m:apply>
		      <m:divide/>
		      <m:cn>1</m:cn>
		      <m:apply>
			<m:root/>
			<m:cn>2</m:cn>
		      </m:apply>
		    </m:apply>
		    <m:apply>
		      <m:ci type="fn">h</m:ci>
		      <m:cn>1</m:cn>
		    </m:apply>
		  </m:apply>
		  <m:apply>
		    <m:ci type="fn">h</m:ci>
		    <m:cn>3</m:cn>
		  </m:apply>
		</m:apply>
	      </m:apply>
	    </m:math>
	  </equation>

	  The length-6 Daubechies wavelet is generated with 
	  <m:math>
	    <m:apply>
	      <m:eq/>
	      <m:ci>α</m:ci>
	      <m:cn>1.35980373244182</m:cn>
	    </m:apply>
	  </m:math>, and 
	  <m:math>
	    <m:apply>
	      <m:eq/>
	      <m:ci>β</m:ci>
	      <m:cn>-0.78210638474440</m:cn>
	    </m:apply>
	  </m:math>.  Note that for 
	  <m:math>
	    <m:apply>
	      <m:eq/>
	      <m:ci>α</m:ci>
	      <m:ci>β</m:ci>
	    </m:apply>
	  </m:math>, we get the length-2 Haar coefficients.  Length-4
	  Daubechies coefficients are found if 
	  <m:math>
	    <m:apply>
	      <m:eq/>
	      <m:ci>α</m:ci>
	      <m:apply>
		<m:divide/>
		<m:pi/>
		<m:cn>3</m:cn>
	      </m:apply>
	    </m:apply>
	  </m:math>, and 
	  <m:math>
	    <m:apply>
	      <m:eq/>
	      <m:ci>β</m:ci>
	      <m:cn>0</m:cn>
	    </m:apply>
	  </m:math>.  Realizing that the values of these parameters is
	  unusual, the formula used to calculate each parameter is
	  given below: 

	  <equation id="eqn13">
	    <m:math>
	      <m:apply>
		<m:eq/>
		<m:ci>α</m:ci>
		<m:apply>
		  <m:arctan/>
		  <m:apply>
		    <m:divide/>
		    <m:apply>
		      <m:divide/>
		      <m:apply>
			<m:plus/>
			<m:apply>
			  <m:minus/>
			  <m:apply>
			    <m:times/>
			    <m:cn>2</m:cn>
			    <m:apply>
			      <m:plus/>
			      <m:apply>
				<m:power/>
				<m:apply>
				  <m:ci type="fn">h</m:ci>
				  <m:cn>0</m:cn>
				</m:apply>
				<m:cn>2</m:cn>
			      </m:apply>
			      <m:apply>
				<m:power/>
				<m:apply>
				  <m:ci type="fn">h</m:ci>
				  <m:cn>1</m:cn>
				</m:apply>
				<m:cn>2</m:cn>
			      </m:apply>
			    </m:apply>
			  </m:apply>
			  <m:cn>1</m:cn>
			</m:apply>
			<m:apply>
			  <m:plus/>
			  <m:apply>
			    <m:ci type="fn">h</m:ci>
			    <m:cn>2</m:cn>
			  </m:apply>
			  <m:apply>
			    <m:ci type="fn">h</m:ci>
			    <m:cn>3</m:cn>
			  </m:apply>
			</m:apply>
		      </m:apply>
		      <m:apply>
			<m:root/>
			<m:cn>2</m:cn>
		      </m:apply>
		    </m:apply>
		    <m:cn>2</m:cn>
		  </m:apply>
		</m:apply>
	      </m:apply>
	    </m:math>
	  </equation>

	  <equation id="eqn14">
	    <m:math>
	      <m:apply>
		<m:eq/>
		<m:ci>β</m:ci>
		<m:apply>
		  <m:minus/>
		  <m:ci>α</m:ci>
		  <m:apply>
		    <m:arctan/>
		    <m:apply>
		      <m:divide/>
		      <m:apply>
			<m:minus/>
			<m:apply>
			  <m:ci type="fn">h</m:ci>
			  <m:cn>2</m:cn>
			</m:apply>
			<m:apply>
			  <m:ci type="fn">h</m:ci>
			  <m:cn>3</m:cn>
			</m:apply>
		      </m:apply>
		      <m:apply>
			<m:plus/>
			<m:apply>
			  <m:ci type="fn">h</m:ci>
			  <m:cn>2</m:cn>
			</m:apply>
			<m:apply>
			  <m:minus/>
			  <m:apply>
			    <m:ci type="fn">h</m:ci>
			    <m:cn>3</m:cn>
			  </m:apply>
			  <m:apply>
			    <m:divide/>
			    <m:cn>1</m:cn>
			    <m:apply>
			      <m:root/>
			      <m:cn>2</m:cn>
			    </m:apply>
			  </m:apply>
			</m:apply>
		      </m:apply>
		    </m:apply>
		  </m:apply>
		</m:apply>
	      </m:apply>
	    </m:math>
	  </equation>

	  Given these equations, we can work both forwards and
	  backwards to determine the filter coefficients of a
	  particular system.  Because
	  <m:math>
	    <m:apply>
	      <m:in/>
	      <m:apply>
		<m:and/>
		<m:ci>α</m:ci>
		<m:ci>β</m:ci>
	      </m:apply>
	      <m:interval>
		<m:apply>
		  <m:minus/>
		  <m:pi/>
		</m:apply>
		<m:pi/>
	      </m:interval>
	    </m:apply>
	  </m:math>, all possible  
	  <m:math>
	    <m:apply>
	      <m:ci type="fn">h</m:ci>
	      <m:ci>n</m:ci>
	    </m:apply>
	  </m:math> can be generated.
	</para>

	<para id="para6">
	  More about Daubechies filters will be discussed <cnxn document="m11150" strength="8">later</cnxn>.
	</para>
      </section>
    </section>
  </content>
  
</document>
