Fourier Transforms3
Fourier Transforms3
Fourier Transforms3
FOURIER TRANSFORMS
1. Introduction
1.1 Definition of the transform and spectrum
Definition: Consider a signal v(t), where t ∈ (−∞, ∞). Let ω be a real number.
The Fourier transform, V (ω), of the signal v(t) is defined by the integral
Z ∞
V (ω) = v(t)e−jωt dt (10.1)
−∞
4
This integral exists whenever
515
516 MATHEMATICS: THE LANGUAGE OF ECE
Z ∞
|v(t)|dt < ∞
−∞
Note that the variable of integration is real, unlike the case of the inverse
Laplace transform. We have seen how to calculate such integrals already, in
Chapter 9.
The Fourier transform can be defined to directly represent frequency in Hertz
by the transformation pair
Z ∞
V (F ) = v(t)e−j2πF t dt,
−∞
Z ∞
v(t) = V (F )ej2πF t dF.
−∞
i.e., the magnitude spectrum is an even function of ω and the phase spectrum
is an odd function of ω. In the sequel, when we plot the spectrum, most of the
time we will only show positive values of ω.
The power spectrum of an aperiodic signal is defined in a manner analogous
to that of a periodic one in Chapter 9, Section 3.3.
Definition: The power spectrum, Ps (ω), of a signal s(t) with Fourier transform
S(ω), is defined as
As with the Fourier series case (see Equation 9.72, page 478), the power
spectrum can also be expressed in decibels:
Therefore,
Equation 10.5 says that the Fourier transform can be found from the Laplace
transform by the substitution s = jω. Inversely, the Laplace transform can be
found from the Fourier transform by the substitution ω = s/j.
Example 10.1. The Laplace transform of the function v(t) = eat u(t) was
found to be
1 In
Chapter 8, we denoted the Laplace transform of v(t) as V (s). We change the notation here to avoid
confusion, since we use V (ω) to denote the Fourier transform of v(t).
518 MATHEMATICS: THE LANGUAGE OF ECE
1
Lv (s) =
s−a
(see Table 1, page 404 in Chapter 8). Therefore, its Fourier transform is
1
V (ω) =
jω − a
————————————————————————————- 4
The Fourier series coefficients Sp (k) can be found from the Fourier transform
S(ω), as follows:
Z T /2
1 2π
Sp (k) = sp (t)e−j T kt
dt
T −T /2
Z
1 T /2 2π
= s(t)e−j T kt dt
T −T /2
Z
1 ∞ 2π
= s(t)e−j T kt dt
T −∞
µ ¶
1 2πk
= S (10.6)
T T
Therefore,
µ ¶
1 2πk
Sp (k) = S (10.7)
T T
Equation 10.7 says that the Fourier series Sp (k) can be found from the Fourier
transform S(ω) by dividing it by T and substituting ω = 2πk T . This is basically
Fourier transforms 519
Z ∞
V (ω) = δ(t − t0 )e−jωt dt
−∞
−jωt0
= e (10.8)
Then
Z
1 ∞
v(t) = V (ω)ejωt dω
2π −∞
Z Z
1 ∞ jωt 1 ∞
= δ(ω − ω0 )e dω + δ(ω + ω0 )ejωt dω
2 −∞ 2 −∞
1 jω0 t 1 −jω0 t
= e + e = cos(ω0 t)
2 2
From the uniqueness property (discussed in Section 3), we recognize V (ω)
as the Fourier transform of the (periodic) cos(ω0 t) signal.
————————————————————————————- 4
We have
520 MATHEMATICS: THE LANGUAGE OF ECE
Z ∞
Pτ (ω) = pτ (t)e−jωt dt
−∞
Z τ /2
= e−jωt dt
−τ /2
1 ¯τ /2
¯
= e−jωt ¯
−jω −τ /2
1 ³ −jωτ /2 ´
= e − e−jω(−τ /2)
−jω
1 ³ jωτ /2 ´
= e − e−jωτ /2
jω
2 sin(ωτ /2)
= sin(ωτ /2) = τ (10.11)
ω ωτ /2
So, finally,
sin(ωτ /2)
Pτ (ω) = τ (10.12)
ωτ /2
where T > 0 is a known parameter. This signal is equal to one period of the
periodic signal in Example 9.11.
The Fourier transform of s(t) can be found as follows:
Z ∞
S(ω) = s(t)e−jωt dt
−∞
Fourier transforms 521
0.8
Magnitude
0.6
0.4
0.2
0
0 1 2 3 4 5 6 7 8 9 10
frequency (Hz)
3.5
2.5
Phase
1.5
0.5
0
0 1 2 3 4 5 6 7 8 9 10
frequency (Hz)
Z T /2
= te−jωt dt
−T /2
1 ¯T /2
−jωt ¯
= 2
e ((−jω)t − 1) ¯
(−jω) −T /2
1 −jωt ¯
¯T /2
= 2
e (jωt + 1)¯
ω −T /2
µ ¶ µ ¶
1 −jω T T 1 jω T T
= e 2 jω + 1 − 2 e 2 −jω + 1
ω2 2 ω 2
µ ¶h i h i
1 T T
−jω 2 T
jω 2 1 −jω T2 jω T2
= jω e + e + e − e
ω2 2 ω2
jT cos(ωT /2) 2j sin(ωT /2)
= − (10.14)
ω ω2
The plot of this spectrum for T = 2 is shown in Figure 10.2. Note again that
the plot is versus frequency F in Hertz, not angular frequency ω.
522 MATHEMATICS: THE LANGUAGE OF ECE
0.8
Magnitude
0.6
0.4
0.2
0
0 1 2 3 4 5 6 7 8 9 10
frequency (Hz)
1
Phase
−1
−2
0 1 2 3 4 5 6 7 8 9 10
frequency (Hz)
————————————————————————————- 4
*Example 10.6. The Fourier series of the periodic ramp, one period of which
is defined as
sp (t) = t, t ∈ [−T /2, T /2)
can be evaluated now using Equation 10.7. We have, from Equation 10.14:
· ¸¯
1 jT cos(ωT /2) 2j sin(ωT /2) ¯¯
Sp (k) = − ¯ 2πk
T ω ω2 ω= T
³ ´ ³ ´
2πk T 2πk T
1 jT cos T 2 2j sin T 2
= 2πk
− ³ ´2
T T
2πk
T
jT cos(πk)
=
2πk
Fourier transforms 523
S(ω) = T Sp (k)|k= T ω
2π
¯
j ¯
= T· cos(kπ)¯¯
kω0 k= T2πω
µ ¶
Tω j
= T · T ω 2π cos π
2π T
2π
µ ¶
jT ωT
= cos
ω 2
and we see that we cannot recover the term with the sin(ωT /2) in Equation
10.14.
————————————————————————————- 4
Table 10.2 summarizes the Fourier transforms of some common signals.
Table 10.1. Fourier transforms of common signals.
v(t) V (ω)
δ(t) 1
1 1
sin(ω0 t) 2j
δ(ω − ω0 ) − 2j
δ(ω + ω0 )
1 1
cos(ω0 t) 2
δ(ω − ω0 ) + 2
δ(ω + ω0 )
at 1
e u(t) jω−a
at ω0
e sin(ω0 t)u(t) 2
(jω−a)2 +ω0
jω−a
eat cos(ω0 t)u(t) 2
(jω−a)2 +ω0
3. Properties
The Fourier transform, like the Laplace transform and Fourier series, pos-
sesses similar properties.
3.1 Uniqueness
For each given function v(t), there exists one and only one function V (ω),
defined via Equation 10.1. Inversely, given a function V (ω), there exists one
524 MATHEMATICS: THE LANGUAGE OF ECE
and only one function v(t), defined via Equation 10.2, that has the function
V (ω) as its Fourier transform.
3.2 Linearity
The following theorem is the analogue of theorems 8.1 and 9.1; it can be
proven in an entirely similar fashion:
Theorem 10.1 Consider given constants ci (real or complex valued) and
functions vi (t) with known Fourier transforms Vi (ω). Let
X
v(t) = ci vi (t) (10.15)
i
3.3 Time-shifting
Let t0 be a given real-valued constant, positive or negative (unlike the Laplace
transform). The analogue of theorems 8.2 and 9.2 is:
Theorem 10.2 Consider a signal v(t), with known Fourier transform V (ω).
Let y(t) = v(t − t0 ), where t0 is a given real-valued constant. The Fourier
transform of the time-shifted signal is given by
Z ∞ Z ∞
4 −jωt
Y (ω) = y(t)e dt = v(t − t0 )e−jωt dt
−∞ −∞
Z ∞
= v(t − t0 )e−jωt0 ejωt0 e−jωt dt
−∞
Z ∞
−jωt0
= e v(t − t0 )e−jω(t−t0 ) dt
−∞
Z ∞
= e−jωt0 v(t)e−jωt dt = e−jωt0 V (ω)
−∞
Fourier transforms 525
t0 = −θ/ω0
From Equation 10.17 and Equation 10.9, page 519, we can now write
3.4 Time-reversal
Consider a signal v(t), with known Fourier transform V (ω).
Theorem 10.3 Let y(t) = v(−t), be the time-reversed version of v(t). The
Fourier transform of the time-reversed signal is given by
Z ∞
4
Y (ω) = y(t)e−jωt dt
−∞
Z ∞
= v(−t)e−jωt dt
−∞
526 MATHEMATICS: THE LANGUAGE OF ECE
Z −∞
= v(u)e−jω(−u) (−du) (10.21)
∞
Z −∞
= − v(u)e−j(−ω)u du
∞
Z ∞
= v(u)e−j(−ω)u du = (10.22)
−∞
= V (−ω)
4
3.5 Time-scaling
Consider a signal v(t), with known Fourier transform V (ω). Let y(t) =
v(at), where a 6= 0 is a given constant. The Fourier transform can be calculated
as follows. Consider first the case a > 0. We have by definition
Z ∞
4
Y (ω) = y(t)e−jωt dt
−∞
Z ∞
= v(at)e−jωt dt
−∞
Z ∞ Z ∞
u du 1 ω
= v(u)e−jω a = v(u)e−j a u du (10.23)
−∞ a a −∞
1
= V (ω/a) (10.24)
a
In Equation 10.23, we have made the substitution of variables u = at.
The case where a < 0 can be derived from Equations 10.20 and 10.24, by
observing that we can write a = −|a|. Then the signal
y1 (t) = v(|a|t)
1
has Fourier transform Y1 (ω) = |a| V (ω/|a|), and the signal
y(t) = y1 (−t)
1 1
has Fourier transform Y (ω) = |a| V (−ω/|a|) = |a| V (ω/a). Therefore,
Theorem 10.4 Consider a signal v(t), with known Fourier transform V (ω).
Let y(t) = v(at), where a 6= 0 is a given constant. The Fourier transform is
given by
1
Y (ω) = V (ω/a) (10.25)
|a|
4
Fourier transforms 527
3.6 Derivatives
The Fourier transform of the n-th derivative of a signal can also be related
to the Fourier transform of the signal itself. The analogue of theorems 9.5 and
8.4 is the following.
Theorem 10.5 Consider a signal v(t) with known Fourier transform V (ω).
The Fourier transform of yk (t) = dk v(t)/dtk , the k-th order derivative of the
signal v(t), is given by
3.7 Examples
*Example 10.8.
Determine the Fourier transform of the signal s(t) in Equation 10.27
0, t < −1,
−1, −1 ≤ t < 0,
s(t) = (10.27)
2, 0 ≤ t ≤ 2,
0, t > 2.
using (a) the direct integration approach and (b) the properties of the Fourier
transform. The signal is shown in Figure 10.3.
Note that when applying the properties, we can determine this transform:
1 using the linearity/time-shifting properties and the transform of (two) cen-
tered pulses pτ (t) (of width τ ):
½
1, −τ /2 ≤ t ≤ τ /2,
pτ (t) = (10.28)
0, otherwise.
One pulse, p1 (t), of width 1, would cover the region −1 ≤ t < 0; the other,
p2 (t), of width 2, would cover the region 0 ≤ t ≤ 2.
2 using the linearity/time-scaling properties and the transform of one cen-
tered pulse p(t) = p1 (t) with unit width. The transform of the pulses was
calculated in Example 10.4.
Approach 1. From Example 10.4, we have
sin(ωτ /2)
Pτ (ω) = τ (10.29)
ωτ /2
528 MATHEMATICS: THE LANGUAGE OF ECE
2.5
1.5
1
s(t)
0.5
−0.5
−1
−1.5
−2
−5 −4 −3 −2 −1 0 1 2 3 4 5
time t
We can write the signal in Equation 10.27 using the two simpler signals in
Equation 10.28 as follows:
Equation 10.30 simply says that the signal s(t) consists of two pulses: one
of width τ = 1, shifted to the left by 0.5 time units and another of width τ = 2,
shifted to the right by 1 time unit. A decomposition like this would avoid using
a time scaling property along with time shifting.
Therefore (applying the linearity and time-shifting properties of the Fourier
transform) and using Equation 10.29:
y(t) = p1 (t/2)
stretches the unit pulse from −1 to +1. We can then write
s(t) = −p1 (t + 0.5) + 2y(t − 1) = −p1 (t + 0.5) + 2p1 (1/2(t − 1)) (10.32)
Then
sin(ω/2)
P1 (ω) = (10.34)
ω/2
and from the time-scaling property
sin(ω)
Y (ω) = 2P1 (2ω) = 2 (10.35)
ω
sin(ω/2) sin(ω)
S(ω) = −ej0.5ω + 2e−jω 2 (10.36)
ω/2 ω
which (of course) agrees with the function S(ω) we derived in Equation 10.31.
Approach 3. From its definition, we have
Z ∞ Z 0 Z 2
−jωt −jωt
S(ω) = s(t)e dt = (−1)e dt + (2)e−jωt dt
−∞ −1 0
Z 0 Z 2
= − e−jωt dt + 2 e−jωt
dt
−1 0
¯0 ¯2
e−jωt ¯¯ e−jωt ¯¯
= − ¯ +2 ¯
−jω ¯−1 −jω ¯0
530 MATHEMATICS: THE LANGUAGE OF ECE
1 ¯0 1 −jωt ¯¯2
¯
= e−jωt ¯ −2 e ¯
jω −1 jω 0
1 h i h
1 −j2ω i
= 1 − ejω − 2 e −1
jω jω
1 h i
= 1 − ejω + 2 − 2e−j2ω
jω
1 h i
= 3 − ejω − 2e−j2ω (10.37)
jω
Equations 10.37 and 10.31 look different. However, using Euler’s identities
in Equation 10.31, we get
1 jω/2 1 jω
j0.5ω 2j
[e − e−jω/2 ] −jω 2j
[e − e−jω ]
S(ω) = −e + 4e
ω/2 ω
1 h ³ ´i
= −ej0.5ω [ejω/2 − e−jω/2 ] + 2e−jω ejω − e−jω
jω
1 h i
= −ejω + 1 + 2 − 2e−j2ω
jω
1 h i
= 3 − ejω − 2e−j2ω (10.38)
jω
which is the same as Equation 10.37.
Figure 10.4 depicts the spectrum of this signal. Note that both negative and
positive frequencies are shown. From the graph, the zeros of the magnitude
spectrum occur at frequencies
F = 1, 2, 3, 4 Hz.
Can you prove that the frequency F = k, where k is an integer, will also yield
a zero, by using Equation 10.37?
————————————————————————————- 4
Write a Matlab script to plot the spectrum of the signal for a frequency range
of 0–10 Hz, with a resolution of 0.001 Hz.
The spectrum consists of two plots, one for the magnitude of S(ω) and one
for the phase of S(ω). A Matlab script to plot the spectrum is the following:
Fourier transforms 531
Magnitude spectrum
4
3
Magnitude
0
−4 −3 −2 −1 0 1 2 3 4
Frequency (Hz)
Phase spectrum
4
2
Phase
−2
−4
−4 −3 −2 −1 0 1 2 3 4
Frequency (Hz)
£ ¤
Figure 10.4. The spectrum of S(ω) = 1
jω
3 − ejω − 2e−j2ω .
S(1) = 1 ;
subplot(2,1,1)
plot( freq_range, magnitude )
xlabel(’frequency (Hz)’)
ylabel(’Magnitude’)
grid on
subplot(2,1,2)
plot( freq_range, phase )
xlabel(’frequency (Hz)’)
ylabel(’Phase’)
grid on
——————————————————————————–
The spectrum produced by this script is shown in Figure 10.5. Note that
the horizontal axis is frequency in Hz, not in radians. When comparing this
spectrum to the one shown in Figure 10.1, notice that the magnitude parts are
equal. The phase parts differ by a factor of 2ω, as you can easily explain by
using the time-shifting property.
————————————————————————————- 4
Example 10.10. For the Fourier transform of the previous example, evaluate
the value of S(ω) for ω = 0 and ω = π/2. Express your results in exponential
form. 4
We have
sin(0)
S(0) = e−j2·0 =1
0
sin(π/2) 1 2
S(π/2) = e−j2π/2 = 2e−jπ = e−jπ
π/2 π π
————————————————————————————- 4
*Example 10.11. Find the Fourier transform of the signal, v(t), in Figure
10.6, by using time scaled and shifted versions of the centered unit pulse signal
p1 (t) in Equation 10.10, page 519, Example 10.4. 4
Fourier transforms 533
0.8
Magnitude
0.6
0.4
0.2
0
0 1 2 3 4 5 6 7 8 9 10
frequency (Hz)
2
Phase
−2
−4
0 1 2 3 4 5 6 7 8 9 10
frequency (Hz)
y1 (t) = p1 (t/2)
y2 (t) = p1 (t/6)
y1 (t) is a centered pulse of width 2 and y2 (t) is a centered pulse of width 6.
Then, from the time-scaling property of the Fourier transform we obtain
534 MATHEMATICS: THE LANGUAGE OF ECE
0
v(t)
−1
−2
−5 −4 −3 −2 −1 0 1 2 3 4 5
t
sin(2ω/2) sin(ω)
Y1 (ω) = 2P1 (2ω) = 2 =2 (10.40)
2ω/2 ω
sin(6ω/2) sin(3ω)
Y2 (ω) = 6P1 (6ω) = 6 =6 (10.41)
6ω/2 3ω
Note that we keep the transform in the form of the sinc function sin(x)/x.
From Figure 10.6, we can write
v(t) = y2 (t + 1) − 2y1 (t − 3)
from which and the time-shifting property of the Fourier transform we get:
sin(3ω) sin(ω)
V (ω) = 6ejω − 4e−3jω (10.42)
3ω ω
Figure 10.7 depicts the spectrum of V (ω). (It looks almost artistic.) Note
that the magnitude spectrum is zero for angular frequencies ω = π, 2π, 3π, . . .,
Fourier transforms 535
or for frequencies f = 0.5, 1, 1.5, 2, . . . Note also how these frequencies are
related to the width (2) of the pulse: they are multiples of 1/pulse width. (This
relationship will be important in later classes.)
Magnitude spectrum
8
6
Magnitude
0
−1.5 −1 −0.5 0 0.5 1 1.5
Frequency f (Hz)
Phase spectrum
4
2
Phase
−2
−4
−1.5 −1 −0.5 0 0.5 1 1.5
Frequency f (Hz)
————————————————————————————- 4
*Example 10.12. Find the Fourier transform of the signal, v(t), in Figure
10.6, by using time scaled and shifted versions of the centered unit pulse signal
p2 (t) in Equation 10.10, example 10.4. 4
Inspecting Figure 10.6, we can write
v(t) = p2 (t + 3) + p2 (t + 1) + p2 (t − 1) − 2p2 (t − 3)
The Fourier transform of the centered unit pulse signal p2 (t) is given by
sin(ω)
P2 (ω) = 2 .
ω
Then, from the linearity and time-shifting properties of the Fourier transform,
we can write
³ ´
V (ω) = ej2π3ω + ej2πω + e−j2πω − 2e−j2π3ω P2 (ω)
536 MATHEMATICS: THE LANGUAGE OF ECE
³ ´ sin(ω)
= ej2π3ω + ej2πω + e−j2πω − 2e−j2π3ω 2
ω
(10.43)
Equations 10.43 and 10.42 appear to be different. Can you reconcile them?
————————————————————————————- 4
Note that the driving force does not have to be a periodic function. We
are interested again in finding vp (t), the particular solution of the equation.
Therefore, throughout this section we will assume that the initial conditions are
0.
As with the Laplace transform approach, the solution methodology is to:
Taking the Fourier transform of both sides of Equation 10.44, and applying
Theorem 10.5, we get
1
Vp (ω) = · Vs (ω) (10.46)
(jω)k + ak−1 (jω)k−1 + · · · + a1 (jω) + a0
Fourier transforms 537
The particular solution can now be found as the following theorem summa-
rizes:
Theorem 10.6 Consider the differential equation 10.44. Let Vs (ω) denote
the Fourier transform of the driving force. Then the particular solution is given
by
Z ∞
Vs (ω)ejωt
vp (t) = dt (10.47)
−∞ (jω)k + ak−1 (jω)k−1 + · · · + a1 (jω) + a0
Note that we will rarely compute the inverse transform but will use the
properties of the transform and tables similar to the Laplace transform to obtain
a complete description of vp (t).
As with the Fourier series case,
Definition: The function
4 1
H(ω) = (10.48)
(jω)k + ak−1 (jω)k−1 + · · · + a1 (jω) + a0
is called the transfer function of the system represented by differential equation
10.44. Note that, in Equation 10.44, there are no derivatives of the driving
force, hence the numerator in Equation 10.48 is equal to 1. Compare Equations
10.48 and 8.58.
4
Set k = 3 in Equation 10.46. From Equations 10.49 and 10.50, the Fourier
transform of the particular solution is given by:
Vs (ω)
Vp (ω) =
(jω)3 + 2(jω)2 + 3(jω) + 4
60πδ(ω − 2π100) + 60πδ(ω + 2π100)
=
−jω 3 − 2ω 2 + 3jω + 4
60πδ(ω − 2π100) + 60πδ(ω + 2π100)
=
(4 − 2ω 2 ) + jω(3 − ω 2 )
60π
= δ(ω − 2π100)
[4 − 2(2π100) ] + j(2π100)[3 − (2π100)2 ]
2
60π
+ δ(ω + 2π100)
[4 − 2(−2π100) ] + j(−2π100)[3 − (−2π100)2 ]
2
(10.51)
60π
= δ(ω − 2π100)
[4 − 2(2π100)2 ] + j[6π100 − (2π100)3 ]
60π
+ δ(ω + 2π100)
[4 − 2(2π100) ] − j[6π100 − (2π100)3 ]
2
60π
= p δ(ω − 2π100)
[4 − 2(2π100) ] + [6π100 − (2π100)3 ]2 ejθ
2 2
60π
+p δ(ω + 2π100)
[4 − 2(2π100)2 ]2 + [6π100 − (2π100)3 ]2 e−jθ
(10.52)
In deriving Equation 10.51, we have used the property of the δ(t) function
in Equation 2.17, page 61.
From Equation 10.19, page 525, the Fourier transform Vp (ω) corresponds to
a cosine time signal, with amplitude equal to
60
p
[4 − 2(2π100) ] + [6π100 − (2π100)3 ]2
2 2
and phase
à !
(2π100)[3 − (2π100)2 ]
θ = arctan ≈ −π/2
4 − 2(2π100)2
Therefore,
60
vp (t) = p cos(2π100t − π/2)
[4 − 2(2π100)2 ]2+ [6π100 − (2π100)3 ]2
60
= p sin(2π100t)
[4 − 2(2π100)2 ]2 + [6π100 − (2π100)3 ]2
Fourier transforms 539
or
fourier(f)
The command
ifourier(F)
syms w
F = 1/(1+j*w) ;
ifourier(F)
——————————————————————————–
returns the expression e−t u(t).
————————————————————————————- 4
Skillset to be developed
Determine analytically the Fourier transform of a signal using the direct
integration approach.
Determine analytically the signal v(t) in the time domain, from knowledge
of its Fourier transform V (ω).
Write Matlab scripts to plot the frequency and power spectrum of a signal.
Problems
Calculation of the Fourier transform
10.1. Using direct integration, find the Fourier transform of the signal
10.2. Using direct integration, find the Fourier transform of the signal
*10.3. Show that the Fourier transform of the triangular pulse signal
½ |t|
s(t) = 1− T , −T ≤ t ≤ T ,
0, otherwise.
is given by
µ ¶2
sin(ωT /2)
S(ω) = T
ωT /2
10.4. Consider a signal s(t) with Fourier transform S(ω). Let ω0 be a real
number. Define a new (complex-valued) signal
*10.5. Consider p(t), the unit pulse signal in Equation 2.12, page 59. Let
t0 = 10.
Determine P (ω) using time shifting and scaling properties, and the results
of Example 10.4.
10.8. Consider the signal h3 (t) shown in Figure 2.16, page 81. Calculate its
Fourier transform.
10.9. Consider the signal h4 (t), defined in Equation 2.41, page 82. Calculate
its Fourier transform.
10.10. Consider the signal h5 (t), defined in Equation 2.43, page 82. Calcu-
late its Fourier transform.
10.11. Consider the centered unit pulse p5 (t), in Example 10.4. Calculate
the Fourier transform of the signal y(t) = 3p5 (t/2).
*10.12. Consider the centered unit pulse p5 (t), in Example 10.4. Calculate
the Fourier transform of the signal y(t) = 3p5 (t/2 + 3).
10.13. Consider the signal v(t), shown in Figure 8.3, page 428.
(a) Calculate the Laplace transform V (s).
(b) Calculate the Fourier transform, V (ω) from the Laplace transform V (s).
sin(2ω) sin(ω)
S(ω) = 8ejω + 2e−jω
2ω ω
Find the Fourier transform of the time-scaled and shifted signal y(t) = s(2t−3).
*10.15. Using the formula for the Fourier transform of the centered unit pulse
in Equation 10.10, page 519, together with time shifting and scaling properties,
find the Fourier transform of the signal defined by
0, t<0
2, 0≤t≤1
s(t) =
−1, 1 < t ≤ 3
0, 3<t
————————————————————————————–
544 MATHEMATICS: THE LANGUAGE OF ECE
*10.16. The Fourier series for the periodic signal, s1 (t), whose single period
is defined by ½
sin(2πt/T ), 0 ≤ t ≤ T /2
s1 (t) =
0, T /2 < t < T
where T = 1/60 seconds, is denoted by S1 (k).
(a) Let a single period of the periodic signal s2 (t), be defined by
½
sin(2πt/T ), 0 ≤ t ≤ T /2
s2 (t) =
0, T /2 < t < T
where T = 1/120 seconds. Find the Fourier series, S2 (k), in terms of S1 (k).
Hint: Use the relationship between the Fourier series and Fourier transform to
obtain both series and look for the similarities.
(b) Let a single period of the periodic signal s3 (t), be defined by
cos(2πt/T ), −T /4 ≤ t ≤ T /4
s3 (t) = 0, −T /2 < t < −T /4
0, T /4 < t < T /2
where T = 1/120 seconds. Find the Fourier series, S3 (k), in terms of S1 (k).
*10.17. The rectified sinusoids are very important in many ECE applications,
particularly in communications and power conversion. It is important that you
understand the relationships among the various forms, half-wave and full-wave
rectified signals, in both the time and frequency domains.
(a) Compute, by direct integration, the Fourier transform of the function,
s1 (t), defined by
0, t ≤ −T /4
s1 (t) = cos(2πt/T ), −T /4 ≤ t ≤ T /4
0, T /4 < t
(b) Assume T = 1/60; plot the magnitude and phase spectra from zero to
10KHz. Plot the power spectrum in decibels as a function of frequency in Hertz.
Recall that the power spectrum in dB is defined by PdB (ω) = 20 log10 (|S(ω)|).
(c) The periodic extension of s1 (t) with a period of T is written
0, −T /2 < t ≤ −T /4
s1p (t) = cos(2πt/T ), −T /4 ≤ t ≤ T /4
0, T /4 < t ≤ T /2
Use the relationship of the Fourier transform and Fourier series to obtain the
Fourier series coefficients of s1p (t).
Fourier transforms 545
(d) Use the information found above to obtain the Fourier series coefficients
of s2p (t) = | cos(2πt/T )|. This is the full-wave rectified sinusoid. Pay partic-
ular attention to the fundamental frequencies of s1p (t) and s2p (t).
(e) Use the information found above and Fourier properties (time scaling and
shifting) to obtain the Fourier series coefficients of s3p (t) = | sin(2π10t/T )|.
————————————————————————————–
*10.18. Let pτ (t) be the unit centered pulse of width τ , in Equation 10.10
of Example 10.4.
(a) Find the power spectrum of a clock generator signal in Equation 2.21,
page 63. For this problem, you will periodically extend the single pulse with
1
a period T = 100 seconds, and set τ = T /2. The power spectrum does not
depend on the phase of the signal, so you may use the symmetric pulse to make
computation easy.
(b) Plot the power spectrum of the clock generator signal in part (a) in deci-
bels.
(a) Compute the Fourier transform of s3 (t). Hint: Use Euler’s identity.
(b) For the case where α = 50, f = 200, and θ = π/4, plot s3 (t) from 0
to 125 ms. Using subplot, plot the magnitude spectrum in decibels from 0 to 2
KHz.
(c) At what frequency does the peak of the spectrum occur? Is this what you
would expect? What parameter(s) would you change to make the peak more
distinct (sharper)? less distinct (broader)?
*10.20. Two time-domain signals, s1(t) and s2(t), and their Fourier trans-
forms, SA(F ) and SB(F ), are shown in Figures 10.8 and 10.9. Using your
knowledge of the location of zeros of the spectra of rectangular pulses, match
the spectra to the corresponding time domain signal. The zeros in the spectra
can be identified as the deep negative dips. Using standard computation, it is
unlikely that an exact zero will be produced.
————————————————————————————–
546 MATHEMATICS: THE LANGUAGE OF ECE
time domain
2
1.5
1
s1(t)
0.5
−0.5
−1
0 1 2 3 4 5 6 7 8
1.5
1
s2(t)
0.5
−0.5
−1
0 1 2 3 4 5 6 7 8
time in seconds
(a) Compute the transfer function of the system represented by the differential
equation.
(b) Compute the output (solution) of the system when the input is
s(t) = cos(5t)
(c) Compute the output (solution) of the system when the input is
s(t) = cos(20t)
Frequency domain
60
40
20
SA(F)
−20
−40
0 5 10 15 20 25
50
SB(F)
−50
0 5 10 15 20 25
frequency in hertz
d2 v(t) dv(t)
2
+ 30 + 900v(t) = 1800vs (t),
dt dt
use the Fourier transform to find the particular solution vp (t), if the input, vs (t)
is defined by
vs (t) = cos(30t)
d2 v(t) R dv(t) 1 1
2
+ + v(t) = vs (t) (10.54)
dt L dt LC LC
where L = 100µH, C = 10µF and vs (t) = 5 cos(2π5000t). Use Fourier
transforms to:
(a) find the particular solution vp (t) when R = 0.6Ω.
(b) find the particular solution vp (t) when R = 12Ω.
548 MATHEMATICS: THE LANGUAGE OF ECE
————————————————————————————–
Does this function represent a Butterworth filter? Hint: you must show that
Equation 10.56 can be cast in the form of Equation 10.55. To do that, find
the conjugate, H ∗ (ω), of the given transfer function and compute
Fourier transforms 549
6 Consider the circuit shown in Figure 10.10. The input is the current source;
the output is the voltage across the resistor. Write the transfer function.
1F
i(t) 1Ω v(t)
7 Show that this circuit implements a Butterworth filter of order 1. Show that
the cutoff frequency is ωc = 1 Hz.
8 Consider the circuit shown in Figure 10.11. The input is the current source;
the output is the voltage across the capacitor. Write the transfer function.
550 MATHEMATICS: THE LANGUAGE OF ECE
1H
1F 1F
i(t) 1Ω v(t)
9 Show that this circuit implements a Butterworth filter of order 3. Show that
the cutoff frequency is ωc = 1 Hz. Show that the transfer function for this
filter is given by
1
|H(ω)| =
(jω)3 + 2(jω)2 + 2(jω) + 1
Fourier transforms 551
4 Suppose that any signal power lower than -30 dB is practically zero. Which
signal consumes more bandwidth?
2 Suppose that any signal power lower than -30 dB is practically zero. Is the
OOK code more expensive, in terms of bandwidth, than the signals in the
previous problem?
1 Determine the magnitude spectra of the BPSK signals in Equations 2.62 and
2.63, page 107. Plot the two spectra.
2 Suppose that any signal power lower than -30 dB is practically zero. Is the
BPSK code more expensive, in terms of bandwidth, than the signals in the
previous two problems?
1 Determine the magnitude spectra of the FSK signals in equations 2.64 and
2.65, page 109. Plot the two spectra.
2 Suppose that any signal power lower than -30 dB is practically zero. Is the
FSK code more expensive, in terms of bandwidth, than the signals in the
previous three problems?
552 MATHEMATICS: THE LANGUAGE OF ECE
*10.30. (ECE402) Binary Phase Shift Keying (BPSK) using raised cosine
signals. Consider the Binary Phase Shift Keying (BPSK) code in problem 2.48.
A variation of this code uses the signal
½A
2 [1 + cos(2πfc t)], |t| < 0.5fc ,
c
c1 (t) = (10.57)
0, |t| ≥ 0.5fc .
to transmit bit 1, while bit 0 is represented by the signal
The duration of the signals c1 (t) and c0 (t) is T seconds, resulting in a rate
R = 1/T bps. Figure 10.12 shows how the ASCII letter V (with ASCII code
1010110) would be transmitted using this code, over a 1 bps channel. Note that
T = 1 sec for such a channel. Compare the signal waveforms to that of the
BPSK code, in Figure 2.26, page 108.
1.5
1 0 1 0 1 1 0
1
m(t)
0.5
0
T 2T 3T 4T 5T 6T 7T
−0.5
−1
−1 0 1 2 3 4 5 6 7 8
1 0 1 0 1 1 0
1
s(t)
−1
−2
−1 0 1 2 3 4 5 6 7 8 9
time (sec)
Figure 10.12. The raised cosine Binary Phase Shift Keying (BPSK) code.
The reason for even considering use of a “complicated” signal like that in
Equation 10.57 is the bandwidth savings it provides, compared to other BPSK
signals (see question 6).
2 Calculate and plot the magnitude spectrum of the signal used to send bit 1.
3 Calculate and plot the magnitude spectrum of the signal used to send bit 0.
Are they the same? Why?
4 Calculate and plot the magnitude spectrum of the signal used to send bit 1,
for the signal used in the code of Problem 2.48.
5 Calculate and plot the magnitude spectrum of the signal used to send bit 0,
for the signal used in the code of Problem 2.48. Are they the same? Why?
6 Is this BPSK code more or less expensive, in terms of bandwidth, than the
code in Problem 2.48? Suppose that any signal power lower than -30 dB is
practically zero.
7 In yet another BPSK variation, bit 1 is represented by the signal
½
Ac , |t| < 0.5fc ,
c1 (t) = (10.59)
0, |t| ≥ 0.5fc .
while bit 0 is represented by the signal
Calculate and plot the magnitude spectrum of the signals c1 (t) and c0 (t).
8 Suppose that any signal power lower than -30 dB is practically zero. Using
the spectrum plots or the formulas you have calculated, show that this BPSK
code is 5 times more expensive, in terms of bandwidth, than the raised cosine
BPSK code.
————————————————————————————–