Basics of sampling and reconstruction

Introduction

Screencast video [⯈]



Module overview

This module will cover one of the most fundamental aspects in signal processing, namely the principle behind converting a continuous-time signal into a discrete-time signal in order to process it with the help of computers.

  1. Sampling of sinusoidal signals - First the concept of sampling is discussed, which represents the conversion from the continuous-time to discrete-time domain. The frequency of the sampled signal is better to be described by a relative frequency, because of the effect of the sampling process. Through the sampling there can also be a loss of spectral information, known as aliasing. This is caused by the uniqueness issue.
  2. Reconstruction of sinusoidal signals [⯈] - Similarly to the sampling, the signal can also be converted back from the discrete-time domain to the continuous-time domain. In order to prevent aliasing the sampling theorem has to be satisfied.
  3. Examples [⯈] - In order to get the reader more acquainted with the sampling procedure, this section includes a screencast video with several examples.



Exercises

In this section several exercises are available, including their answers. The exercises marked in blue are explained by means of more extensive pencast videos.

Video quiz

Given the following situation:
CD converter, question 1.
where x(t)=cos(2πf0t) and fs>2f0. Two cases are defined as:
case 1: fs=fs1x1[n]=cos(θ1n)
case 2: fs=1/3fs1x2[n]=cos(θ2n)
In the equation θ2=αθ1, what is the value of α?






Given the following situation:
CD and DC converter, question 2.
where x(t)=cos(2π50t+π3)sin(2π700tπ3).
What is the minimal sampling rate fs to obtain no aliasing, thus y(t)=x(t).







Given the following situation:
CD and DC converter, question 3.
where
fsi=500 [samples/sec]
fso=400 [samples/sec]
x(t)=cos(2π100t) and y(t)=cos(2πfyt)
What is the value of fy?








Given the following situation:
CD and DC converter including input signal, question 4.
How do we have to choose fs in order to obtain no aliasing, thus y(t)=x(t)?






Exercise bundle

This browser does not support PDFs. Please download the PDF to view it: Download PDF.

Answers

Download the answers here.

Pencast videos [⯈]

The above video player contains a playlist of all pencast videos which can be expanded by clicking the playlist icon in the upper-right corner.



MATLAB lab

Accompanied to this modules are some exercises in MATLAB, which will test your knowledge of the module and will help improve your MATLAB skills.

Lab assignment

This browser does not support PDFs. Please download the PDF to view it: Download PDF.

MATLAB demo [⯈]



Summary

Relation absolute- and relative- frequency:θ=ωTs=2πffs


Frequency Symbol Unit
Absolute frequency f [Hz]
Absolute radial frequency ω [rad/sec]
Relative frequency θ [-]

The spectral representation of a discrete-time signal x[n] is not unique. The relative frequency representation is periodic with period 2π.

Frequency Symbol Fundamental Interval
Relative frequency θ π<θπ
Absolute frequency f fs2<ffs2
Absolute radial frequency ω πfs<ωπfs

An ideal D-to-C convertor chooses the relative frequencies within the FI: π<θπ. This implies that in case a continuous-time absolute frequency is converted with an C-to-D conversion outside the FI, such a frequency will always first be mapped inside the FI and the D-to-C convertor will use this frequency for the conversion of the discrete-time signal into a continuous-time signal. This mapping back into the FI is called Aliasing

Only in case of sampled sinusoidal signals the ideal D-to-C convertor in effect replaces n by tfs.

Interpolation equation:y(t)=n=x[n]p(tnTs)
Sampling theorem: A continuous-time signal x(t) with frequencies no higher than fmax can be reconstructed exactly from its samples x[n]=x(t)t=nTs, if samples are taken at a rate fs=1/Ts, that is greater than 2fmax.