Lecture-17-18 Discrete Time Fourier Series
Lecture-17-18 Discrete Time Fourier Series
❑ There are some important differences that should be considered before the start
of Fourier series representation of discrete time signals.
❑ In particular, the Fourier series representation of a discrete-time periodic signal is
a finite series, as opposed to the infinite series representation required for
continuous time periodic signals.
❑ Consequently, there is no mathematical issues of convergence (as in CT Fourier
series Convergence).
❑ The following properties of CT exponential are different in the DT case:
➢ CT:- The larger the magnitude of 𝝎𝒐 , the higher is the rate of oscillation in the signal.
➢ CT:- Complex exponential, i.e., 𝒆𝒋𝝎° 𝒕 is periodic for any value of 𝝎𝒐 .
➢ CT:- All harmonically related complex exponentials 𝑒 𝑗𝑘𝝎° 𝑡 , are distinct.
➢ DT:-If 𝜔0 is increased from 0 𝑡𝑜 2𝜋, the rate of oscillations first increases (from 0 to 𝜋) and
then decreases (from 𝜋 to 0) .
➢ DT:- The signal with frequency 𝝎𝟎 is identical to signals with frequencies 𝝎𝟎 ± 𝟐𝝅 , 𝝎𝟎 ±
𝟒𝝅, and so on.
➢ DT:- There are only N distinct periodic exponentials, i.e., k + N n = k n
Lecture-17-18 Revision-Chap#01
❑ The set of all discrete-time complex exponential signals that are periodic with period N
is given by
2
k n = e
jk n
jk0 n
=e N
, k = 0, 1, 2,
2𝜋
❑ All these signals have fundamental frequencies that are multiples of and thus are
𝑁
harmonically related.
❑ Consider the discrete-time complex exponential:
j (0 + 2 ) n
e = e j0 n e j 2 n e j0 n e j 2 n = cos ( 2 n ) + j sin ( 2 n )
❑ This is a consequence of the fact that discrete-time complex exponentials which differ in
frequency by a multiple of 𝟐𝝅 are identical. In fact, the signal with frequency 𝝎𝟎 is
identical to signals with frequencies 𝝎𝟎 ± 𝟐𝝅 , 𝝎𝟎 ± 𝟒𝝅, and so on……..
❑ Specifically, 𝜙0 𝑛 = 𝜙N 𝑛 , 𝜙𝟏 𝑛 = 𝜙N+1 𝑛 , and, in general,
k n = k + rN n , r, N Z 3
Lecture-17-18 Revision-Chap#01
k k k
❑ Therefore, Since the sequences 𝝓𝒌 [𝒏] are distinct only over a range of N
successive values of 𝒌, the summation in eq. (above) need only include terms
over this range.
❑ Thus, the summation is on k, as k varies over a range of N successive integers,
beginning with any value of k. We indicate this by expressing the limits of the
summation as 𝒌 = 𝑵 . That is,
2
x n = akk n =
jk n
ak e jk0 n = ak e N
k = N k = N k = N
Lecture-17-18 Revision-Chap#01
k = N k = N k = N
❑ For example,
➢ 𝒌 could take on the values 𝒌 = 0, 1, … , 𝑁 − 1, or 𝒌 = 1, … , 𝑁 or 𝑘 = 2, 4, … , 𝑁 +
1, and so on …..
➢ In either case, the same set of complex exponential sequences appears in the
summation on the right-hand side of the above equation.
k = N k = N k = N
❑ If 𝒌 is in the range from 0 to 𝑵 − 1,
x n = a11 n + a22 n + + aN N n
❑ If 𝑘 is in the range from −1 to 𝑵 − 2 (write yourself in the same manner as above)
x n e − jr0n = j ( k − r )0 n
ak e
n = N n = N k = N
**A closed-form expression is an expression that can be evaluated/solved in a finite number of steps/operations.
n = N
e
j ( k − r )0 n
=
k = N
cos ( k − n ) 0 n + j
k = N
sin ( k − n ) 0 n
N n = N
❑ This provides a closed-form expression for obtaining the Fourier series
coefficients, and the discrete-time Fourier series pair (synthesis
equation and analysis equation) will be:
1 2
ak = x
n e − jk0 n
, where 0 = Analysis equation
N n = N N
The set of coefficients 𝒂𝒌 are often called the Fourier series coefficients or
the spectral coefficients of 𝒙 𝒏 .
© 1992–2008 R. C. Gonzalez & R. E. Woods
Digital Image Processing, 3rd ed.
Lecture-17-18 Revision-Chap#01
x n =
2j
(
1 j0 n − j0 n
e −e , ) Comparing with synthesis eq., by
inspection that
2 2 1 1
1 j N n 1 − j N n a1 = , a−1 = −
x n = e − e 2j 2j
2j 2j
The remaining coefficients over the
2 2
1 j (1) n 1 j ( −1) n interval of summation <N>, i.e.,
x n = e N
− e N
one time period are zero.
2j 2j
1 1
These coefficients repeat with period N; aN +1 = , aN −1 = −
2j 2j
x n = A cos (0 n + )
www.ImageProcessingPlace.c
1 1 om
a1 = , a−1 = − Chapter 10
2j 2j Image Segmentation 2
N0 = m
2
0
Let 0 = , then
5 2
N0 = m 5. where ‘m’ is the smallest integer
2 that makes N as an integer.
5
The Fourier series coefficients for this example with 𝑁 = 5 are illustrated in Figure
below:
𝟐𝝅
Fourier coefficients for 𝒙 𝒏 = 𝐬𝐢𝐧 5
𝒏 .
The fact that they repeat periodically is indicated. However, only one period is utilized in the synthesis equation.
𝟐𝝅
Fourier coefficients for 𝒙 𝒏 = 𝐬𝐢𝐧 𝟑 𝒏.
5
1 j N n − j N n 3 j N n − j N n 1 j N n + 2 − j N n + 2
2 2 2 2 4 4
x n = 1 + e −e + e +e + e +e
2j 2 2
3 1 j N n 3 1 − j N n 1 j 2 j 2 N n 1 − j 2 − j 2 N n
2 2 2 2
x n = 1 + + e + − e + e e + e e
2 2 j 2 2 j 2 2
j
e 2
= cos + j sin j
a0 = 1 2 2
1 j 2 1
3 1 3 1 a2 = e j
a1 = + − j 2 2 ak = 0, k 0, 1, 2
2 2j 2 2
3 1 3 1
1 j 2
a−2 = e
1
− j or ak = 0, k 2
a−1 = − + j 2 2
2 2j 2 2
e a
www.ImageProcessingPlace.c
om
Chapter k10
Image Segmentation
m ak
ak
Show that 𝒙[𝒏] = 𝑨 𝒄𝒐𝒔(𝑩𝒏 + 𝑪), and specify numerical values for the constants A,
B, and C.
Fact-4: Given:
1 9
x n = 50.
2
10 n =0
Using Parseval’s relation,
1
x n = a
2
=50.
2
k
N N N
ak = 50,
2
n =−1
8
a−1 + a0 + a1 + ak = 50,
2 2 2 2
n=2
8
( 5) + a0 + ( 5 ) + ak = 50,
2 2 2 2
n=2
8
a0 + ak = 0, ak = 0 for k = 2,3,
2 2
,8.
n=2
2
jk n 8
1 jk n
N
a e
N
k
0
= ak e
−1
5
,
j n − j n
= a−1e 5
+ a1e 5
,
j n − j n
= 5e 5
+ 5e 5
,
j 5 n − j 5 n
e +e
= 5 2 ,
2
A = 10,
B= ,
= 10 cos n. 5
5
C=0
© 1992–2008 R. C. Gonzalez & R. E. Woods
Digital Image Processing, 3rd ed.
Lecture-17-18 Revision-Chap#01
x2 n ⎯→
FS
bk
where,
1
a0 = a3 = a1 = 1, and b0 = b1 = b2 = b3 = 1.
2
Using the multiplication property in Table 3.1, determine the Fourier series coefficients
𝒄𝒌 for the signal 𝒈[𝒏] = 𝒙𝟏 [𝒏] 𝒙𝟐 [𝒏].
g n = x1 n x2 n ⎯→
FS
ck = ab l k −l
l = N
3
ck = al bk −l a0bk −0 + a1bk −1 + a2bk − 2 + a3bk −3
l =0
3
ck = al bk −l bk −0 + 2bk −1 + 2bk − 2 + bk −3
l =0
3
ck = al bk −l 1 + 2 + 2 + 1 = 6
l =0
N −Chapter
1 j 2 kn 10
a k = e N
Image Segmentation
n =0
a) Show that
a k = N for k = 0, N , 2 N , 3N ,
a) Show that a k = 0 whenever k is not an integer multiple of N. (use the finite sum
formula)
2
j kn
a k = e N
N
a) Let k = pN , p Z . Then,
om
Chapter 10 a k = e N
n =0
Image Segmentation
N −1 2 N −1 N −1
j pNn
a pN = e N
= e j 2 pn
1 = N
n =0 n =0 n =0
q + N −1 2
c) Let j kn
a k = e N
where q is some arbitrary integer.
n=q
Now, 2 N −1 2
j kq j kn
a k = e N
e N using part ‘b’ a k =0, if k pN , p Z
n =0
26
Lecture-17-18 Revision-Chap#01
27
Lecture-17-18 Revision-Chap#01
➢ This shows that the exponential at 𝝎𝟎 + 𝟐𝝅 is the same as that at frequency 𝝎𝟎 . This is
very different when compared with the continuous-time exponential case, in which the
signals 𝑒 𝑗𝜔0 𝑡 are all distinct for distinct values of 𝝎𝟎 . 28
Lecture-17-18 Revision-Chap#01
= (e ) = ( −1)
j n j n n
e e j = cos ( n ) + j sin ( n )
➢ It shows that the signal oscillates rapidly, changing sign at each point in time.
29
Digital Image Processing, 3rd ed.
Lecture-17-18 Revision-Chap#01
3
1
Some Basic Signals
◼ Periodicity Property of Discrete-Time Complex Exponential
◼ It is also considerable value in discrete-time signal and system analysis to consider
sets of harmonically related periodic exponentials – i.e., periodic exponentials with a
common period N.
2
jk n
k n = e N
, k = 0, 1, 2,
◼ In the continuous-time case, all of the harmonically related complex exponentials
2𝜋
𝑗𝑘 𝑡
𝑒 𝑇 , are distinct. However, this is not the case in discrete time. Specifically,
2
j ( k + N ) n
k + N n = e N
,
2
jk n
k + N n = e N
e j 2 n = k n
◼ This implies that there are only N distinct periodic exponentials. For example,
2 4 2
j n j n j ( N −1) n
0 n = 1, 1 n = e N
, 2 n = e N
, , N −1 n = e N
are all distinct, and any other is identical to one of these (e.g., 𝜙0 𝑛 = 𝜙N 𝑛 ).
Digital Image Processing, 3rd ed.
Lecture-17-18 Revision-Chap#01
Point to Ponder
The low-frequency (that is, slowly varying) discrete-time exponentials have values of 𝜔0 near
0, 2𝜋, and any other even multiple of 𝜋, while the high frequencies (corresponding to rapid
variations) are located near 𝜔0 = ±𝜋 and other odd multiples of ±𝜋 .
e
j0 ( n + N )
= e j0 n Periodicity
e j0 n .e j0 N = e j0 n ,
e j0 N = 1,
j ( 2 m )
e j0 N = e ,
0 N = 2 m,
2
N = m
0
2
x ( t ) = 3cos 4t + x (t ) = e
j ( t −1)
T0 =
3 0
0 = 4 0 =
2 2
T0 = T0 = 2
4 2
➢ Example-02 ➢ Example-04
x ( t ) = cos ( 2t )
2
x ( t ) = 10sin 20 t + + 10
3
1 + cos ( 4t )
x (t ) = 0 = 20
2
2 1
1 1
x ( t ) = + cos ( 4t ) T0 =
2 2 20 10
0 = 4
2
T0 =
4 2