<?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="decimation">

  <name>Decimation</name>

  <metadata>
  <md:version>2.10</md:version>
  <md:created>2002/01/10</md:created>
  <md:revised>2003/01/15</md:revised>
  <md:authorlist>
      <md:author id="schniter">
      <md:firstname>Phil</md:firstname>
      
      <md:surname>Schniter</md:surname>
      <md:email>schniter@ee.eng.ohio-state.edu</md:email>
    </md:author>
  </md:authorlist>

  <md:maintainerlist>
    <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="emaloney">
      <md:firstname>Erin</md:firstname>
      
      <md:surname>Maloney</md:surname>
      <md:email>emaloney@rice.edu</md:email>
    </md:maintainer>
    <md:maintainer id="schniter">
      <md:firstname>Phil</md:firstname>
      
      <md:surname>Schniter</md:surname>
      <md:email>schniter@ee.eng.ohio-state.edu</md:email>
    </md:maintainer>
  </md:maintainerlist>
  
  <md:keywordlist>
    <md:keyword>aliasing</md:keyword>
    <md:keyword>anti-aliasing</md:keyword>
    <md:keyword>decimation</md:keyword>
    <md:keyword>decimator</md:keyword>
    <md:keyword>downsampling</md:keyword>
  </md:keywordlist>

  <md:abstract>Introduction to decimation.</md:abstract>
</metadata>



  <content>

    <para id="para1">
Decimation is the process of filtering and downsampling a signal 
to decrease its effective sampling rate, as illustrated in <cnxn target="figfix" strength="9"/>. The 
filtering is employed to prevent aliasing that might otherwise 
result from downsampling.
    </para>


    <figure id="figfix">
      <media type="image/png" src="m10445fig1.png"/>
    </figure>

    <para id="para2">
To be more specific, say that

<m:math display="block">
  <m:apply>
  <m:eq/>
    <m:apply>
	<m:ci type="fn">
	  <m:msub>
	    <m:mi>x</m:mi>
	    <m:mi>c</m:mi>
	  </m:msub>
	</m:ci>
	<m:ci>t</m:ci>
    </m:apply>
    <m:apply>
      <m:plus/>
        <m:apply>
	    <m:ci type="fn">
	      <m:msub>
	        <m:mi>x</m:mi>
	        <m:mi>l</m:mi>
	      </m:msub>
	    </m:ci>
	    <m:ci>t</m:ci>
        </m:apply>   
        <m:apply>
	    <m:ci type="fn">
	      <m:msub>
	        <m:mi>x</m:mi>
	        <m:mi>b</m:mi>
	      </m:msub>
	    </m:ci>
	    <m:ci>t</m:ci>
        </m:apply> 
      </m:apply>
    </m:apply>
  </m:math>

where <m:math>
         <m:apply>
	    <m:ci type="fn">
	      <m:msub>
	        <m:mi>x</m:mi>
	        <m:mi>l</m:mi>
	      </m:msub>
	    </m:ci>
	    <m:ci>t</m:ci>
        </m:apply> 
      </m:math>
is a lowpass component bandlimited to
<m:math><m:apply><m:divide/><m:cn>1</m:cn>
<m:apply><m:times/><m:cn>2</m:cn><m:ci>M</m:ci>
<m:ci>T</m:ci></m:apply></m:apply></m:math> Hz and
<m:math>
         <m:apply>
	    <m:ci type="fn">
	      <m:msub>
	        <m:mi>x</m:mi>
	        <m:mi>b</m:mi>
	      </m:msub>
	    </m:ci>
	    <m:ci>t</m:ci>
        </m:apply> 
      </m:math>
is a bandpass component with energy between
<m:math><m:apply><m:divide/><m:cn>1</m:cn>
<m:apply><m:times/><m:cn>2</m:cn><m:ci>M</m:ci>
<m:ci>T</m:ci></m:apply></m:apply></m:math> and
<m:math><m:apply><m:times/><m:apply><m:divide/><m:cn>1</m:cn>
<m:apply><m:times/><m:cn>2</m:cn>
<m:ci>T</m:ci></m:apply></m:apply><m:ci>Hz</m:ci></m:apply></m:math>.
If sampling
<m:math>
         <m:apply>
	    <m:ci type="fn">
	      <m:msub>
	        <m:mi>x</m:mi>
	        <m:mi>c</m:mi>
	      </m:msub>
	    </m:ci>
	    <m:ci>t</m:ci>
        </m:apply> 
      </m:math>
with interval <m:math><m:ci>T</m:ci></m:math>
yields an unaliased discrete representation
<m:math><m:apply><m:ci type="fn" class="discrete">x</m:ci>
<m:ci>m</m:ci></m:apply></m:math>, then decimating
<m:math><m:apply><m:ci type="fn" class="discrete">x</m:ci>
<m:ci>m</m:ci></m:apply></m:math> by a factor 
<m:math><m:ci>M</m:ci></m:math> will yield
<m:math><m:apply><m:ci type="fn" class="discrete">y</m:ci>
<m:ci>n</m:ci></m:apply></m:math>, an unaliased 
<m:math><m:apply><m:times/><m:ci>M</m:ci>
<m:ci>T</m:ci></m:apply></m:math>-sampled
representation of lowpass component
<m:math>
         <m:apply>
	    <m:ci type="fn">
	      <m:msub>
	        <m:mi>x</m:mi>
	        <m:mi>l</m:mi>
	      </m:msub>
	    </m:ci>
	    <m:ci>t</m:ci>
        </m:apply> 
      </m:math>.

</para>


    <para id="para3">

We offer the following justification of the previously described 
decimation procedure. From the sampling theorem, we have

	  <m:math display="block">
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:ci type="fn">X</m:ci>
		  <m:apply>
		    <m:exp/>
		    <m:apply>
		      <m:times/>
		      <m:imaginaryi/>
		      <m:ci>ω</m:ci>
		    </m:apply>
		  </m:apply>
		</m:apply>
		<m:apply>
             <m:plus/>
             <m:apply>
		  <m:times/>
		  <m:apply>
		    <m:divide/>
		    <m:cn>1</m:cn>
		    <m:ci>T</m:ci>
		  </m:apply>
		  <m:apply>
		    <m:sum/>
		    <m:bvar><m:ci>k</m:ci></m:bvar>
		    <m:condition><m:ci>k</m:ci></m:condition>
		    <m:apply>
		      <m:ci type="fn">
			<m:msub>
			  <m:mi>X</m:mi>
			  <m:mi>l</m:mi>
			</m:msub>
		      </m:ci>
		      <m:apply>
                  <m:times/>
                  <m:imaginaryi/>
                  <m:apply>
			<m:divide/>
			<m:apply>
			  <m:minus/>
			  <m:ci>ω</m:ci>
			  <m:apply>
			    <m:times/>
			    <m:cn>2</m:cn>
			    <m:pi/>
			    <m:ci>k</m:ci>
			  </m:apply>
			</m:apply>
			<m:ci>T</m:ci>
		      </m:apply>
                 </m:apply>
		    </m:apply>
		  </m:apply>
		</m:apply>
            <m:apply>
		  <m:times/>
		  <m:apply>
		    <m:divide/>
		    <m:cn>1</m:cn>
		    <m:ci>T</m:ci>
		  </m:apply>
		  <m:apply>
		    <m:sum/>
		    <m:bvar><m:ci>k</m:ci></m:bvar>
		    <m:condition><m:ci>k</m:ci></m:condition>
		    <m:apply>
		      <m:ci type="fn">
			<m:msub>
			  <m:mi>X</m:mi>
			  <m:mi>b</m:mi>
			</m:msub>
		      </m:ci>
		      <m:apply>
                  <m:times/>
                  <m:imaginaryi/>
                  <m:apply>
			<m:divide/>
			<m:apply>
			  <m:minus/>
			  <m:ci>ω</m:ci>
			  <m:apply>
			    <m:times/>
			    <m:cn>2</m:cn>
			    <m:pi/>
			    <m:ci>k</m:ci>
			  </m:apply>
			</m:apply>
			<m:ci>T</m:ci>
		      </m:apply>
                 </m:apply>
		    </m:apply>
		  </m:apply>
		</m:apply>
            </m:apply>
            </m:apply>
	  </m:math>

    </para>

  <para id="para4">

The bandpass component 
<m:math>
         <m:apply>
	    <m:ci type="fn">
	      <m:msub>
	        <m:mi>X</m:mi>
	        <m:mi>b</m:mi>
	      </m:msub>
	    </m:ci>
          <m:apply>
          <m:times/>
            <m:imaginaryi/>
            <m:ci>Ω</m:ci>
          </m:apply>
        </m:apply> 
      </m:math>
is the removed by
<m:math>
  <m:apply>
  <m:divide/>
    <m:pi/>
    <m:ci>M</m:ci>
  </m:apply>
</m:math>-lowpass
filtering, giving

	  <m:math display="block">
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:ci type="fn">V</m:ci>
		  <m:apply>
		    <m:exp/>
		    <m:apply>
		      <m:times/>
		      <m:imaginaryi/>
		      <m:ci>ω</m:ci>
		    </m:apply>
		  </m:apply>
		</m:apply>
		<m:apply>
		  <m:times/>
		  <m:apply>
		    <m:divide/>
		    <m:cn>1</m:cn>
		    <m:ci>T</m:ci>
		  </m:apply>
		  <m:apply>
		    <m:sum/>
		    <m:bvar><m:ci>k</m:ci></m:bvar>
		    <m:condition><m:ci>k</m:ci></m:condition>
		    <m:apply>
		      <m:ci type="fn">
			<m:msub>
			  <m:mi>X</m:mi>
			  <m:mi>l</m:mi>
			</m:msub>
		      </m:ci>
		      <m:apply>
                  <m:times/>
                  <m:imaginaryi/>
                  <m:apply>
			<m:divide/>
			<m:apply>
			  <m:minus/>
			  <m:ci>ω</m:ci>
			  <m:apply>
			    <m:times/>
			    <m:cn>2</m:cn>
			    <m:pi/>
			    <m:ci>k</m:ci>
			  </m:apply>
			</m:apply>
			<m:ci>T</m:ci>
		      </m:apply>
                 </m:apply>
		    </m:apply>
		  </m:apply>
		</m:apply>
	      </m:apply>
	  </m:math>

Finally, downsampling yields

<equation id="eq1">
	  <m:math display="block">
	      <m:apply>
		<m:eq/>
		<m:apply>
		  <m:ci type="fn">Y</m:ci>
		  <m:apply>
		    <m:exp/>
		    <m:apply>
		      <m:times/>
		      <m:imaginaryi/>
		      <m:ci>ω</m:ci>
		    </m:apply>
		  </m:apply>
		</m:apply>
		<m:apply>
		  <m:times/>
		  <m:apply>
		    <m:divide/>
		    <m:cn>1</m:cn>
                <m:apply>
                  <m:times/>
		      <m:ci>M</m:ci>
                  <m:ci>T</m:ci>
                </m:apply>
		  </m:apply>
		  <m:apply>
		    <m:sum/>
                <m:bvar><m:ci>p</m:ci></m:bvar>
		    <m:lowlimit>
                      <m:cn>0</m:cn>
                </m:lowlimit>
                <m:uplimit>
                  <m:apply>
                    <m:minus/>
                      <m:ci>M</m:ci>
                      <m:cn>1</m:cn>
                  </m:apply>
                </m:uplimit>
		  <m:apply>
		    <m:sum/>
		    <m:bvar><m:ci>k</m:ci></m:bvar>
		    <m:condition><m:ci>k</m:ci></m:condition>
		    <m:apply>
		      <m:ci type="fn">
			<m:msub>
			  <m:mi>X</m:mi>
			  <m:mi>l</m:mi>
			</m:msub>
		      </m:ci>
		      <m:apply>
                  <m:times/>
                  <m:imaginaryi/>
                  <m:apply>
			<m:divide/>
			<m:apply>
			  <m:minus/>
                      <m:apply>
                       <m:divide/>
                         <m:apply>
                           <m:minus/>
              	           <m:ci>ω</m:ci>
                             <m:apply>
                              <m:times/>
                                <m:cn>2</m:cn>
                                <m:pi/>
                                <m:ci>p</m:ci>
                             </m:apply>              
                         </m:apply>
                         <m:ci>M</m:ci>
                      </m:apply>
			  <m:apply>
			    <m:times/>
			    <m:cn>2</m:cn>
			    <m:pi/>
			    <m:ci>k</m:ci>
			  </m:apply>
			</m:apply>
			<m:ci>T</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:times/>
		      <m:ci>M</m:ci>
                  <m:ci>T</m:ci>
                </m:apply>
		  </m:apply>
		  <m:apply>
		    <m:sum/>
                <m:bvar><m:ci>p</m:ci></m:bvar>
		    <m:lowlimit>
                      <m:cn>0</m:cn>
                </m:lowlimit>
                <m:uplimit>
                  <m:apply>
                    <m:minus/>
                      <m:ci>M</m:ci>
                      <m:cn>1</m:cn>
                  </m:apply>
                </m:uplimit>
		  <m:apply>
		    <m:sum/>
		    <m:bvar><m:ci>k</m:ci></m:bvar>
		    <m:condition><m:ci>k</m:ci></m:condition>
		    <m:apply>
		      <m:ci type="fn">
			<m:msub>
			  <m:mi>X</m:mi>
			  <m:mi>l</m:mi>
			</m:msub>
		      </m:ci>
		      <m:apply>
                  <m:times/>
                  <m:imaginaryi/>
                  <m:apply>
			<m:divide/>
			<m:apply>
			  <m:minus/>
                      <m:ci>ω</m:ci>
                      <m:apply>
                      <m:times/>
                        <m:apply>
                        <m:times/>
                          <m:cn>2</m:cn>
                          <m:pi/>
                        </m:apply>
                        <m:apply>
                        <m:plus/>
                          <m:apply>
                          <m:times/>
                            <m:ci>k</m:ci>
                            <m:ci>M</m:ci>
                          </m:apply>     
                          <m:ci>p</m:ci>         
                        </m:apply>
                      </m:apply>
			</m:apply>
                  <m:apply>
                  <m:times/>
			  <m:ci>M</m:ci>
			  <m:ci>T</m:ci>
		      </m:apply>
                 </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:times/>
		      <m:ci>M</m:ci>
                  <m:ci>T</m:ci>
                </m:apply>
		  </m:apply>
		  <m:apply>
		    <m:sum/>
		    <m:bvar><m:ci>l</m:ci></m:bvar>
		    <m:condition><m:ci>l</m:ci></m:condition>
		    <m:apply>
		      <m:ci type="fn">
			<m:msub>
			  <m:mi>X</m:mi>
			  <m:mi>l</m:mi>
			</m:msub>
		      </m:ci>
		      <m:apply>
                  <m:times/>
                  <m:imaginaryi/>
                  <m:apply>
			<m:divide/>
			<m:apply>
			  <m:minus/>
                      <m:ci>ω</m:ci>
                      <m:apply>
                      <m:times/>
                          <m:cn>2</m:cn>
                          <m:pi/>
                          <m:ci>l</m:ci>
                      </m:apply>
			</m:apply>
                  <m:apply>
                  <m:times/>
			  <m:ci>M</m:ci>
			  <m:ci>T</m:ci>
		      </m:apply>
                 </m:apply>
		    </m:apply>
		  </m:apply>
		</m:apply>
            </m:apply>
	      </m:apply>
	  </m:math>
</equation>

which is clearly a <m:math><m:apply><m:times/><m:ci>M</m:ci>
<m:ci>T</m:ci></m:apply></m:math>-sampled version of
<m:math>
         <m:apply>
	    <m:ci type="fn">
	      <m:msub>
	        <m:mi>x</m:mi>
	        <m:mi>l</m:mi>
	      </m:msub>
	    </m:ci>
	    <m:ci>t</m:ci>
        </m:apply> 
      </m:math>.
A frequency-domain illustration for
<m:math><m:apply><m:eq/><m:ci>M</m:ci><m:cn>2</m:cn></m:apply></m:math>
appears in <cnxn target="fig3" strength="9"/>.

    </para>


    <figure id="fig3">
     <media type="image/png" src="m10445fig2.png"/>
    </figure>


  </content>
</document>
