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.
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.
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.
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:
where x(t)=cos(2π⋅f0t) and fs>2⋅f0. Two cases are defined as:
case 1: fs=fs1→x1[n]=cos(θ1n)
case 2: fs=1/3⋅fs1→x2[n]=cos(θ2n)
In the equation θ2=α⋅θ1, what is the value of α?
Given the following situation:
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:
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:
How do we have to choose fs in order to obtain no aliasing, thus y(t)=x(t)?
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
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<f≤fs2
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 t⋅fs.
Interpolation equation:y(t)=∞∑n=−∞x[n]p(t−n⋅Ts)
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=n⋅Ts, if samples are taken at a rate fs=1/Ts, that is greater than 2fmax.