Sampling Theorem Course - Notes - 7824 - Part2 PDF
Sampling Theorem Course - Notes - 7824 - Part2 PDF
Sampling Theorem Course - Notes - 7824 - Part2 PDF
Andrew W. H. House
11 May 2004
Obviously in this case, Xc (ω) = 0 for |ω| > ωm , so the signal xc (t) is band-limited.
Let’s choose a value for ωs so that ωs > 2ωm and therefore
2π 2π
>2·
Ts Tm
or in other words Tm > 2Ts or Ts < 12 Tm and then we can plot as below.
ox o x o
ox ox xo ox o x o xo
ox o x ox x ox
ox xo
ox xo xo x o x ox o
xo
x x
o
x o
o x ox
o
o
x x ox
x
o
o
xo o x ox
x xo x ox
Let’s assume that there are two sets of samples take, those marked with ‘x’, and those
marked with ‘o’. We will also assume that in each case, ωs > 2ωm so our sampling rate is
high enough.
Clearly, each set of samples will have different numerical values, since they record different
points on the original CT signal. Yet the sampling theorem suggests that the CT signal can
be reconstructed from either of the sets of samples.
How can this be? It seems counter-intuitive that different sets of samples will reproduce the
same signal. But that is what happens because the samples are used to recover the lowest
frequency signal. That basically means that there is only one possible signal that could have
produced a given set of samples, so long as the assumptions of the sampling theorem were
true.
3 Sampling Processes
We have seen, at least in a qualitative sense, how samples can represent a CT signal provided
that we follow the rules of the sampling theorem. At this point, we should consider how we
can get samples – after wall, without being able to sample real-world CT signals, there’s not
much point to this course. That is, we need to know how to get samples from a CT signal.
xc(t)
? xd(n)
There are many methods to get samples, some of which are covered in the following subsec-
tions.
The sample-and-hold circuit (sometimes called zero-order hold) converts the held analog
values to digital values during the hold interval Ts , as show in the following diagram.
Essentially the circuit latches values of the analog signal during the hold intervals, producing
a series of stepped samples.
This type of sampling multiplies a CT signal xc (t) with a pulse signal p(t) to get a series of
pulses that map certain portions of the original CT signal.
In the pulse train, note that the area of each pulse is 1 – that is, its height can be considered
to be ∆, and its width 1/∆.
When the signal is modulated (multiplied) with the pulses, we get xp (t) = xc (t) · p(t) shown
below.
Each pulse of the modulated signal has an area that is approximately equal to the value of
the original signal at the middle of the pulse width, which is also the value of the DT sample
we would ideally have. This is illustrated in the above diagram.
Each modulated pulse above can be integrated to get its area, as follows.
Z nTs +δ
xd (n) = xp (t)dt ≈ xc (nTs )
nTs −δ
Thus, the area of each modulated pulse is roughly the value of the sample. So this is a good
approximation, but is not exactly what we want.
But, what happens when the width of the pulse narrows, as ∆1 → 0? Our pulse train
p(t) becomes an impulse train, with each pulse being a Dirac delta, δ(t). This suggests a
theoretical approach to sampling, which we explore in the next section.
a.k.a. train to samples. Note that the above is exactly equal, not approximately as with
PAM sampling before.
Therefore, this is our model of ideal sampling, which is the kind of sampling that the sampling
theorem means.
Model:
xc(t) xp(t) Convert impulse xd(n)
train to samples
p(t)
This model above represents the basic concept of a C/D converter, illustrated below.
Ts
Obviously not. In this case, ωs = 2ωm which does not satisfy the require-
ments of the sampling theorem, which requires strictly greater than, not
greater than or equal. Therefore, we cannot guarantee recovery of xc (t).
3. What if Xc (ω) = 0 for |ω| > 15000π?
Here, ωm = 15000π.
Is (ωs = 20000π) > (2ωm = 30000π)?
Obviously not. Therefore, we cannot guarantee recovery of xc (t).
4. What if |Xc (ω)| = 0 for |ω| > 5000π?
This question isn’t quite so obvious. Can we figure out, from the given
information, whether the signal is band-limited? And if so, can we figure
out whether the sampling was fast enough?
Consider that since Xc (ω) = |Xc (ω)| ejΘ(ω) , we can be sure that Xc (ω) = 0
for |ω| > 5000π since in that range, Xc (ω) = 0 · ejΘ(ω) = 0.
So signal is band-limited, and ωm = 5000π.
Is (ωs = 20000π) > (2ωm = 10000π)?
Yes. Therefore, we can guarantee recovery of xc (t).
5. What if <e{Xc (ω)} = 0 for |ω| > 5000π?
Again, this is not an obvious question. Consider this representation of
Xc (ω).
Xc (ω) = <e{Xc (ω)} + j=m{Xc (ω)}
Here, we know that the real part is zero outside of the giving band limits
but we know nothing about the imaginary part. Thus, we can’t say for sure
whether it will be zero or non-zero.
Since we don’t know for sure whether Xc (ω) = 0 for |ω| > 5000π, we cannot
guarantee recovery.
6. What if Xc (ω) ∗ Xc (ω) = 0 for |ω| > 15000π?
This is another tricky situation. What can we figure out about Xc (ω)?
Thinking back to previous signals courses, recall that
Xc (ω) ∗ Xc (ω)
l Fourier Transform
2πxc (t)xc (t)
In this case, the band limits of the convolved transforms of Xc (ω) are twice
those of Xc (ω) itself.
So Xc (ω) = 0 for |ω| > 7500π, so ωm = 7500π.
Is (ωs = 20000π) > (2ωm = 15000π)?
Yes. Therefore, we can guarantee recovery of xc (t).
Additionally, for self-practice, you could identify the range of values of ωs for
each scenario that would allow the signal to be recovered.
The above is true provided that xc (t) satisfies some conditions (such as being
finite energy) or if xc (t) is periodic with a non-standard Fourier transform.
So, for our frequency analysis, let us assume our sampling function is
+∞
X
p(t) = δ(t − nTs )
n=−∞
x(t) · y(t)
l Fourier Transform
1
2π
X(ω) ∗ Y (ω)
so
We know that p(t) is periodic, with period Ts , so P (ω) is a non-standard transform based
on its Fourier series.
So, based on our knowledge of p(t) and Fourier transforms and series,
+∞
X 2π
P (ω) = 2π ak δ(ω − kωs ) where ωs =
k=−∞
Ts
which gives
+∞
1 X
Xp (ω) = Xc (ω − kωs )
Ts k=−∞
With the Xc (ω) shown above, Xp (ω) would be as illustrated below, since it is made of shifted
and scaled replicas of Xc (ω):
1 1 1 1 1
Xp (ω) = . . .+ Xc (ω +2ωs )+ Xc (ω +ωs )+ Xc (ω)+ Xc (ω −ωs )+ Xc (ω −2ωs )+. . .
Ts Ts Ts Ts Ts
We have non-overlapping replicas of Xc (ω) so that we can recover Xc (ω) from Xp (ω) by
defining a reconstruction filter, Hr (ω), defined and illustrated as follows.
−ωs ωs
Ts 2
<ω< 2
Hr (ω) =
0 |ω| > ω2s
Then, if the recovered signal Xr (ω) = Xp (ω) · Hr (ω) then Xr (ω) = Xc (ω) and therefore
xr (t) = xc (t). Thus xc (t) is recovered from its samples!
The reconstruction occurs because it ignores all but one of the replicas (the lowest frequency
one) and undoes the scaling of Xc (ω) that occurs in Xp (ω) for that replica.
So this is how the reconstruction works when the sampling theorem conditions are met.
What happens if the conditions aren’t met? We’ll do two examples to explore that scenario.
With ωs 6> 2ωm , Xp (ω) would be as shown below. Note that the replicas overlap,
and thus add, in places. The added version is shown by a dotted line where it
differs from the non-adding parts.
The above Xr (ω) is clearly not the same as Xc (ω), and so xr (t) 6= xc (t) and the
original signal is not recovered. This is called aliasing due to undersampling.
Also, as an item of note, in this case, if ωs = 2ωm exactly, we still would have
been able to recover Xc (ω). But that is only because, in this case, Xc (ω) = 0
at the band-limits, and thus there is no interference. That’s why the sampling
theorem can only guarantee recovery for ωs > 2ωm , since ωs = 2ωm depends on
the values of Xc (ω).
In the prior example, since the sampling rate was too low to meet the conditions of the
sampling theorem, we had overlapping replicas which added together and thus obscured the
original value. This is called aliasing.
In this case, even if weird stuff didn’t happen due to the undefined boundary
of the reconstruction filter, the replicas clearly cancel each other out. Thus,
the recovered signal will have no frequency components. xr (t) will be a DC
signal, which is clearly wrong.
We see aliasing again in this example. Aliasing can be avoided by following the
requirements of the sampling theorem.
We’ll conclude with an example covering ideal sampling, before looking at the concept of
aliasing in a bit more detail.
for each of the following signals, plot Xp (ω) which is the Fourier transform of
xp (t) = xc (t)·p(t), then determine what happens when that signal is reconstructed
by the ideal reconstruction filter Hr (ω) with cut-off frequency ωs /2, as defined
below.
Ts |ω| < ω2s
Hr (ω) =
0 |ω| > ω2s
sin 4000πt
1. Case 1: xc (t) = πt
Here, we can use the CT Fourier Transform tables in the textbooks to figure
out what Xc (ω) will be.
sin 4000πt
xc (t) = πt
l Fourier Transform
1 |ω| < 4000π
Xc (ω) =
0 |ω| > 4000π
This means that the ideal reconstruction filter will successfully recover the
original signal. In other words, Xr (ω) = Xp (ω) · Hr (ω) = Xc (ω).
This means that the ideal reconstruction filter will successfully recover the
original signal. In other words, Xr (ω) = Xp (ω) · Hr (ω) = Xc (ω).
More examples involving aliasing will follow in the next section.