Lecture 5
Sampling and Reconstruction

C.L. Wyatt

2025-07-11

Today’s lecture is a review of material from ECE 2714. We will review sampling theory, the Nyquist rate, practical sampling, and reconstruction techniques.

Recall that hybrid systems are ones with both CT and DT signals. Such systems require the ability to:

Ideal Sampling Theory

The impulse train

xp(t)=n=δ(tnT)x_p(t) = \sum\limits_{n = -\infty}^{\infty} \delta(t-nT)

is a periodic function with period TT. It has a CTFS representation

xp(t)=k=akej2πTktx_p(t) = \sum\limits_{k = -\infty}^{\infty} a_k e^{j\frac{2\pi}{T} k t}

where ak=1Ta_k = \frac{1}{T} are the coefficients.

Its Fourier Transform is then

Xp(jω)=k=δ(ω2πTk)X_p(j\omega) = \sum\limits_{k = -\infty}^{\infty} \delta\left(\omega-\frac{2\pi}{T} k\right)

Given ans arbitrary signal x(t)x(t) we can represent sampling as multiplication by the impulse train

xs(t)=x(t)xp(t)=n=x(t)δ(tnT)=n=x(nT)δ(tnT)x_s(t) = x(t)\cdot x_p(t) = \sum\limits_{n = -\infty}^{\infty} x(t)\cdot\delta(t-nT) = \sum\limits_{n = -\infty}^{\infty} x(nT)\cdot\delta(t-nT)

where we define x(nT)=x[n]x(nT) = x[n].

In the frequency domain, multiplication is convolution

Xs(jω)=12πXp(jω)*X(jω)=1Tn=X(j(ω2πTn))X_s(j\omega) = \frac{1}{2\pi} X_p(j\omega)*X(j\omega) = \frac{1}{T} \sum\limits_{n = -\infty}^{\infty} X\left(j\left(\omega - \frac{2\pi}{T} n \right) \right)

The frequency content of the sampled signal consists of copies of the frequency content of the original signal called images.

Nyquist Criteria

To ensure there is do distortion in the signal we require the images do not overlap, which occurs if T<12BT < \frac{1}{2B}, where BB is the bandwidth of x(t)x(t). This is the Nyquist criteria. In terms of a frequency this sample rate is then ω0=2πT\omega_0 = \tfrac{2\pi}{T} radians-per-second.

If BB is finite then x(t)x(t) is and we can choose TT to meet the Nyquist rate.

If x(t)x(t) is not band-limited, or T>12BT > \frac{1}{2B}, then we get aliasing. To illustrate this suppose the band-limited signal has a magnitude spectrum like the following

"cartoon of the magnitude spectrum of a band-limited signal"

After sampling (y(t)=x1(t)x2(t)y(t) = x_1(t)\cdot x_2(t)) and assuming ω0>4πB\omega_0 > 4\pi B the spectrum of the sampled signal is:

"cartoon of the magnitude spectrum of a band-limited signal after sampling at the Nyquist rate"

If instead ω0<4πB\omega_0 < 4\pi B the images overlap and we get aliasing, where high frequency content gets added to the lower frequency content. This is shown below with the lighter lines showing the images and the heavier line showing their sum.

"cartoon of the magnitude spectrum of a band-limited signal after sampling below the Nyquist rate"

Practical Sampling

While ideal sampling thoery is mathematically simple and intuitive, ideal sampling cannot be implemented in practice because we cannot generate ideal impulses for the impulse train.

We can generate short pulses that approximate the impulse however. This is done through a circuit called a sample and hold.

"schematic of a circuit implementing a sample and hold"

where p(t)p(t) is the pulse train with pulse width P<TP < T p(t)=n=u(tnT)u(tPnT)p(t) = \sum\limits_{n=-\infty}^{\infty} u(t-nT) - u(t-P-nT)

We can model this as well since p(t)=xp(t)*[u(t)u(tp)]p(t) = x_p(t) * \left[u(t) - u(t-p) \right]. The sampled signal is then

xs(t)={xp(t)*[u(t)u(tp)]}x(t)x_s(t) = \left\{ x_p(t) * \left[u(t) - u(t-p) \right] \right\}\cdot x(t) or in the frequency domain Xs(jω)={Xp(jω)[u(t)u(tp)]}*X(jω)X_s(j\omega) = \left\{ X_p(j\omega)\cdot \mathcal{F}\left[u(t) - u(t-p) \right] \right\}* X(j\omega)

To store the output of the sample and hold we need to quantize it into NN bits, e.g. N=8,12,24N = 8, 12, 24 are common choices. The mathematical model for this process is outsie the scope of 3704 (but is covered in 4624 DSP), so we will assume x(nT)x(nT) \in\mathbb{R}.

Is we have no control over the input x(t)x(t) and have to fix a sample time, we need to force x(t)x(t) to be approximately band-limited usng a mow-pass filter called an anti-aliasing filter, whose cut-off frequency is relative to the sampling frequency.

Reconstruction Theory

Given a DT signal x[n]x[n] and a sample spacing TT, we can define a corresponding CT signal as xp(t)=n=x[n]δ(tnT)x_p(t) = \sum\limits_{n = -\infty}^{\infty} x[n] \, \delta(t-nT) the impulse train with each impulse weighted by the DT signal.

CT signal reconstruction can be viewed from two different (but equivalent) perspectives. In the time domain perspective, the CT signal x(t)x(t) corresponding to a DT signal x[n]x[n] can be viewed as interpolation, where the values of the CT signal are equal to the DT signal at intervals of the sample time, i.e. x(nT)=x[n]x(nT) = x[n], and in between the value of x(t)x(t) is interpolated. If the interpolation is of zero-order, the value at x(nT)x(nT) is held constant until x(nT+T)x(nT+T). This is called a zero-order hold, and can mathematically modeled as convolution of the weighted impulse train with a pulse p(t)=u(t)u(tT)p(t) = u(t) - u(t-T) whose width is the sample time, called the interpolation function. y(t)=p(t)*xp(t)y(t) = p(t)*x_p(t)

Practical Reconstruction

As we have seen before we cannot physically represent the impulse train nor the ideal low-pass filter. Thus practical reconstruction uses an approximation of the ideal reconstruction filter by a digital-to-analog converter (DAC), followed by a causal (and thus physically possible) low-pass filter.

Zero-order hold using an R-2R ladder

A zero-order hold DAC can be implemented by a circuit called a resistor ladder. Consider a digital output with NN bits and a reference voltage VrefV_{ref} (for example an 8-bit output port on a micro-controller using CMOS 3.3v logic).

If this port is connected to a resistor network consisting of resistor values RR and 2R2R as follows

Resistor ladder circuit diagram for a sample and hold system. The circuit has a horizontal row of six 2R resistors, each connected to a digital input labeled b0 through b3, bN2, and bN1. Below each 2R resistor is a corresponding R resistor, forming vertical pairs. The right end of the ladder connects to a voltage output labeled positive V-o.

then depending on the bit pattern at the output port VV, the output of the buffer op-amp will be Vo=VrefV2NV_o = V_{ref}\frac{V}{2^N}

If the port value is changed every sample time TT, then the resister ladder and buffer op-amp combine to implement a zero-order hold circuit.

Reconstruction(anti-imaging) filter

The zero-order hold is followed by the reconstruction (anti-imaging) filter which low-pass filters the output and smooths-out the jumps from value to value.

With or without and anti-imaging filter the signal is positive. This can be followed by a summing amplifier to shoft the signal voltage up or down and amplify as needed.

Lab 1

The goal of Lab 1 is the build a sampling and reconstruction system using a protoboard and microcontroller, an implement a DT identity (all-pass) system.