<?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="m11465">
  <name>Systems view of sampling and reconstruction</name>
  <metadata>
  <md:version>1.20</md:version>
  <md:created>2003/07/18 02:40:00 GMT-5</md:created>
  <md:revised>2004/03/02 05:07:03.863 US/Central</md:revised>
  <md:authorlist>
    <md:author id="Anders">
      <md:firstname>Anders</md:firstname>
      
      <md:surname>Gjendemsjo</md:surname>
      <md:email>gjendems@NO-SPAM.tele.ntnu.no</md:email>
    </md:author>
  </md:authorlist>

  <md:maintainerlist>
    <md:maintainer id="Anders">
      <md:firstname>Anders</md:firstname>
      
      <md:surname>Gjendemsjo</md:surname>
      <md:email>gjendems@NO-SPAM.tele.ntnu.no</md:email>
    </md:maintainer>
  </md:maintainerlist>
  
  <md:keywordlist>
    <md:keyword>Sampling</md:keyword>
    <md:keyword>Reconstruction</md:keyword>
    <md:keyword>Hold</md:keyword>
    <md:keyword>System</md:keyword>
  </md:keywordlist>

  <md:abstract>An overview of sampling and reconstruction on a system level.</md:abstract>
</metadata>

  <content>
    <section id="s1">
        <name>Ideal reconstruction system</name>
	<para id="s1p1">
	    <cnxn target="f1"/> shows the ideal reconstruction system based
	    on the results of the Sampling theorem <cnxn document="m11423">proof</cnxn>.
	</para>
	<para id="s1p2">
	    <cnxn target="f1"/> consists of a sampling device which produces a time-discrete sequence 
	    <m:math><m:apply><m:ci><m:msub><m:mi>x</m:mi><m:mi>s</m:mi></m:msub></m:ci><m:ci>n</m:ci></m:apply></m:math>.
	    The reconstruction filter, <m:math><m:apply><m:ci>h</m:ci><m:ci>t</m:ci></m:apply></m:math>, is
	    an ideal analog 
	    <link src="http://ccrma-www.stanford.edu/~jos/Interpolation/sinc_function.html">sinc</link> filter, with
	    <m:math>
		<m:apply>
		    <m:eq/>
		    <m:apply>
		    	<m:ci>h</m:ci>
			<m:ci>t</m:ci>
		    </m:apply>
		    <m:apply>
	        	<m:ci type="fn">sinc</m:ci>
			<m:apply>
			    <m:divide/>
			    <m:ci>t</m:ci>
	        	    <m:ci><m:msub><m:mi>T</m:mi><m:mi>s</m:mi></m:msub></m:ci>
			</m:apply>
	      	    </m:apply>
	    	</m:apply>
	    </m:math>. We can't apply the time-discrete sequence
	    <m:math><m:apply><m:ci><m:msub><m:mi>x</m:mi><m:mi>s</m:mi></m:msub></m:ci><m:ci>n</m:ci></m:apply></m:math>
	    directly to the analog filter <m:math><m:apply><m:ci>h</m:ci><m:ci>t</m:ci></m:apply></m:math>.
	    To solve this problem we turn the sequence into an analog signal using <cnxn document="m11450">delta functions</cnxn>.
	    Thus we write
	    <m:math>
		<m:apply>
		    <m:eq/>
		    <m:apply>
		    	<m:ci><m:msub><m:mi>x</m:mi><m:mi>s</m:mi></m:msub></m:ci>
			<m:ci>t</m:ci>
		    </m:apply>
		    <m:apply>
			<m:sum/>
			<m:bvar><m:ci>n</m:ci></m:bvar>
			<m:lowlimit><m:apply><m:minus/><m:infinity/></m:apply></m:lowlimit>
			<m:uplimit><m:infinity/></m:uplimit>
	        	<m:apply>
			    <m:times/>
			    <m:apply>
				<m:ci><m:msub><m:mi>x</m:mi><m:mi>s</m:mi></m:msub></m:ci>
				<m:ci>n</m:ci>
			    </m:apply>
			    <m:apply>
				<m:ci>δ</m:ci>
				<m:apply>
				    <m:minus/>
				    <m:ci>t</m:ci>
				    <m:apply>
					<m:times/>
					<m:ci>n</m:ci>
					<m:ci>T</m:ci>
				    </m:apply>
				</m:apply>
			    </m:apply>
			</m:apply>
	      	    </m:apply>
	    	</m:apply>
	    </m:math>.

	    <figure id="f1">
	        <media type="image/jpg" src="ideal.jpg"/>
		<caption>Ideal reconstruction system</caption>
	    </figure>
	    But when will the system produce an output
	    <m:math>
		<m:apply>
		    <m:eq/>
		    <m:apply>
			<m:ci><m:mover accent="true"><m:mi>x</m:mi><m:mo>̂</m:mo></m:mover></m:ci>
		        <m:ci>t</m:ci>
		    </m:apply>
		    <m:apply>
			<m:ci>x</m:ci>
			<m:ci>t</m:ci>
		    </m:apply>
		</m:apply>
	    </m:math>?
	    According to the <cnxn document="m11419" target="s4">sampling theorem</cnxn> we have
	    <m:math>
		<m:apply>
		    <m:eq/>
		    <m:apply>
			<m:ci><m:mover accent="true"><m:mi>x</m:mi><m:mo>̂</m:mo></m:mover></m:ci>
		        <m:ci>t</m:ci>
		    </m:apply>
		    <m:apply>
			<m:ci>x</m:ci>
			<m:ci>t</m:ci>
		    </m:apply>
		</m:apply>
	    </m:math>
	    when the sampling frequency, <m:math><m:ci><m:msub><m:mi>F</m:mi><m:mi>s</m:mi></m:msub></m:ci></m:math>,
	    is at least twice the highest frequency component of <m:math><m:apply><m:ci>x</m:ci><m:ci>t</m:ci></m:apply></m:math>.


	</para>
    </section>

    <section id="s2">
        <name>Ideal system including anti-aliasing</name>
	<para id="s2p1">
	    To be sure that the reconstructed signal is free of aliasing it is customary to
	    apply a lowpass filter, an <cnxn document="m11443" target="s1ss2p4">anti-aliasing filter</cnxn>, before
	    sampling as shown in <cnxn target="f2"/>.
	    
	    <figure id="f2">
	        <media type="image/jpg" src="sampling_antialias.jpg"> 
                  <!--<param name='height' value='76'/>
                    <param name='width' value='657'/>-->
                </media>
		<caption>
		    Ideal reconstruction system with <cnxn document="m11443" target="s1ss2p4">
		    anti-aliasing filter</cnxn>
		</caption>
	    </figure>
	    Again we ask the question of when the system will produce an output
	    <m:math>
		<m:apply>
		    <m:eq/>
		    <m:apply>
			<m:ci><m:mover accent="true"><m:mi>x</m:mi><m:mo>̂</m:mo></m:mover></m:ci>
			<m:ci>t</m:ci>
		    </m:apply>
		    <m:apply>
		        <m:ci>s</m:ci>
			<m:ci>t</m:ci>
		    </m:apply>
		</m:apply>
	    </m:math>?
	    If the signal is entirely confined within the passband of the lowpass filter we will
	    get perfect reconstruction if <m:math><m:ci><m:msub><m:mi>F</m:mi><m:mi>s</m:mi></m:msub></m:ci></m:math>
	    is high enough.
	</para>
	<para id="s2p2">
	    But if the anti-aliasing filter removes the "higher" frequencies, (which in fact is the job
	    of the anti-aliasing filter), we will <emphasis>never</emphasis> be able
	    to <emphasis>exactly</emphasis> reconstruct the original signal,
	    <m:math><m:apply><m:ci>s</m:ci><m:ci>t</m:ci></m:apply></m:math>. If
	    we sample fast enough we can reconstruct <m:math><m:apply><m:ci>x</m:ci><m:ci>t</m:ci></m:apply></m:math>,
	    which in most cases is satisfying.
	</para>

	<para id="s2p3">
	    The reconstructed signal,
	    <m:math>
		<m:apply>
		    <m:ci><m:mover accent="true"><m:mi>x</m:mi><m:mo>̂</m:mo></m:mover></m:ci>
		    <m:ci>t</m:ci>
		</m:apply>
	    </m:math>, will not have aliased frequencies. This is essential for further use of the signal.
		    
	</para>
    </section>

    <section id="s3">
        <name>Reconstruction with hold operation</name>
	<para id="s3p1">
	    To make our reconstruction system realizable there are many things to look into.
	    Among them are the fact that any practical reconstruction system must input finite length pulses into the
            reconstruction filter. This can be accomplished by the <cnxn document="m11458">hold operation</cnxn>.
	    To alleviate the distortion caused by the hold opeator we apply the output from the hold device
	    to a compensator. The compensation can be as accurate as we wish, this is cost and application consideration.
	    <figure id="f3">
	        <media type="image/jpg" src="sampling_hold.jpg"/>
		<caption>
		    More practical reconstruction system with a <cnxn document="m11458">
		    hold component</cnxn>
		</caption>
	    </figure>
	    By the use of the hold component the reconstruction will not be exact, but as mentioned
	    above we can get as close as we want.
	</para>
    </section>
    <section id="s4">
       <list id="l1" type="inline">
       <item><cnxn document="m11419">Introduction</cnxn></item>
       <item><cnxn document="m11423">Proof</cnxn></item>
       <item><cnxn document="m11443">Illustrations</cnxn></item>
       <item><cnxn document="m11549">Matlab example</cnxn></item>
       <item><cnxn document="m11458">Hold operation</cnxn></item>
       <item><cnxn document="m11448">Aliasing applet</cnxn></item>
       <item><cnxn document="m11442">Exercises</cnxn></item>
       </list>
    </section>

 
  </content>
  
</document>
