7 Properties of Fourier Series and Complex Fourier Spectrum

Module -7
Properties of Fourier Series and Complex Fourier Spectrum.
Objective:To understand the change in Fourier series coefficients due to different signal
operations and to plot complex Fourier spectrum.
Introduction:
The Continuous Time Fourier Series is a good analysis tool for systems with
periodicexcitation. Understanding properties of Fourier series makes the work simple in
calculating the Fourier series coefficients in the case when signals modified by some basic
operations.
Graphical representation of a periodic signal in frequency domain represents Complex
Fourier Spectrum.
Description:
Properties of continuous- time Fourier series
The Fourier series representation possesses a number of important properties
that are useful for various purposes during the transformation of signals from one form to other .
Some of the properties are listed below.
[x1(t) and x2(t)] are two periodic signals with period T and with Fourier series
coefficients Cn and Dn respectively.
1) Linearity property
The linearity property states that, if
𝐹𝑠
x1(t) Cn and x2 (𝑑)
then
proof:
𝐹𝑠
𝐹𝑠
Dn
Ax1(t)+Bx2(t) ACn+BDn
From the definition of Fourier series, we have
1 𝑑 +𝑇
FS[Ax1(t)+Bx2(t)]=𝑇 𝑑 0 [𝐴π‘₯1 𝑑 + 𝐡π‘₯2 𝑑 ]𝑒 βˆ’π‘—π‘› πœ” 0 𝑑 𝑑𝑑
0
1
𝑑 +𝑇
1
=𝐴(𝑇 𝑑 0 π‘₯1 𝑑 𝑒 βˆ’π‘—π‘› πœ” 0 𝑑 𝑑𝑑)+ 𝐡(𝑇
0
=ACn+BDn
𝐹𝑠
or
Ax1(t)+Bx2(t) ACn+BDn
proved
2) Time shifting property
The time shifting property states that,if
𝐹𝑠
x(t) Cn
𝐹𝑠
then
x(t-t0) 𝑒 βˆ’π‘—π‘› πœ” 0 𝑑0 Cn
proof: From the definition of Fourier series, we have
𝑗𝑛 πœ” 0 𝑑
x(t)= ∞
𝑛=βˆ’βˆž 𝐢𝑛 𝑒
𝑗𝑛 πœ” 0 (π‘‘βˆ’π‘‘ 0 )
x(t-t0)= ∞
𝑛=βˆ’βˆž 𝐢𝑛 𝑒
∞
βˆ’π‘—π‘› πœ” 0 𝑑 𝑗𝑛 πœ” 0 𝑑
= 𝑛=βˆ’βˆž [𝐢𝑛 𝑒
]𝑒
-1
βˆ’π‘—π‘› πœ” 0 𝑑
=FS [Cn𝑒
]
x(t-t0)
𝐹𝑠
𝑒 βˆ’π‘—π‘› πœ” 0 𝑑0 Cn
𝑑 0 +𝑇
π‘₯2
𝑑0
𝑑 𝑒 βˆ’π‘—π‘› πœ” 0 𝑑 𝑑𝑑)
Or
3) Time reversal property
The time reversal property states that,if
proved
𝐹𝑠
x(t) Cn
𝐹𝑠
then
x(-t) C-n
proof: From the definition of Fourier series, we have
𝑗𝑛 πœ” 0 𝑑
x(t)= ∞
𝑛=βˆ’βˆž 𝐢𝑛 𝑒
∞
x(-t)= 𝑛=βˆ’βˆž 𝐢𝑛 𝑒 𝑗𝑛 πœ” 0 (βˆ’π‘‘)
substituting n = -p in the right hand side, we get
𝑗 (βˆ’π‘)πœ” 0 (βˆ’π‘‘)
𝑗𝑝 πœ” 0 𝑑
x(-t)= βˆ’βˆž
= ∞
𝑝=∞ πΆβˆ’π‘ 𝑒
𝑝=βˆ’βˆž πΆβˆ’π‘ 𝑒
substituting p = n , we get
𝑗𝑛 πœ” 0 𝑑
x(-t)= ∞
= FS-1[C-n]
𝑛=βˆ’βˆž πΆβˆ’π‘› 𝑒
𝐹𝑠
x(-t) C-n
proved
4) Time scaling property
The time scaling property states that , if
𝐹𝑠
x(t) Cn
𝐹𝑠
x(𝛼t) Cn with Ο‰0→𝛼ω0
From the definition of Fourier series, we have
𝑗𝑛 πœ” 0 𝑑
x(t)= ∞
𝑛=βˆ’βˆž 𝐢𝑛 𝑒
𝑗𝑛 πœ” 0 𝛼𝑑
𝑗𝑛 (πœ” 0 𝛼)𝑑
x(𝛼t)= ∞
= ∞
= FS-1[Cn]
𝑛=βˆ’βˆž 𝐢𝑛 𝑒
𝑛=βˆ’βˆž 𝐢𝑛 𝑒
where Ο‰0→𝛼ω0.
then
proof:
𝐹𝑠
x(𝛼t) Cn with fundamental frequency of 𝛼ω0
proved.
5) Time differential property
The time differential property states that, if
𝐹𝑠
x(t) Cn
𝑑π‘₯ (𝑑) 𝐹𝑠
jnω0Cn
From the definition of Fourier series, we have
𝑗𝑛 πœ” 0 𝑑
x(t)= ∞
𝑛=βˆ’βˆž 𝐢𝑛 𝑒
Differentiating both sides with respect to t, we get
then
proof:
𝑑π‘₯ (𝑑)
𝑑𝑑
=
𝑑𝑑
∞
𝑛=βˆ’βˆž
𝐢𝑛
𝑑(𝑒 𝑗𝑛 πœ” 0 𝑑 )
𝑑𝑑
=
∞
𝑛=βˆ’βˆž
𝑑π‘₯ (𝑑) 𝐹𝑠
𝑑𝑑
𝐢𝑛 (jnΟ‰0 ) 𝑒 𝑗𝑛 πœ” 0 𝑑
𝑗𝑛 πœ” 0 𝑑
= ∞
= FS-1[jnω0 Cn]
𝑛=βˆ’βˆž 𝐢𝑛 (jnΟ‰0 ) 𝑒
jnω0Cn
6) Time integration property
Proved.
The time integration property states that,
𝐹𝑠
if
x(t) Cn
𝑑
π‘₯
βˆ’βˆž
then
𝜏 π‘‘πœ
𝐹𝑆
𝐢𝑛
(if C0 = 0)
𝑗𝑛 πœ” 0
proof:
From the definition of Fourier series, we have
𝑗𝑛 πœ” 0 𝑑
x(t)= ∞
𝑛=βˆ’βˆž 𝐢𝑛 𝑒
𝑑
𝑑
𝑗𝑛 πœ” 0 𝜏
π‘₯ 𝜏 π‘‘πœ = βˆ’βˆž [ ∞
]π‘‘πœ
𝑛=βˆ’βˆž 𝐢𝑛 𝑒
βˆ’βˆž
Interchanging the order of integration and summation, we get
𝑑
π‘₯
βˆ’βˆž
𝜏 π‘‘πœ =
𝑑
∞
𝑗𝑛 πœ” 0 𝜏
π‘‘πœ
𝑛=βˆ’βˆž 𝐢𝑛 βˆ’βˆž 𝑒
=
=
𝑑
π‘₯
βˆ’βˆž
𝜏 π‘‘πœ
𝐹𝑆
∞
𝑛=βˆ’βˆž
𝑒 𝑗𝑛 πœ” 0 𝜏
] π‘‘βˆ’βˆž
𝑗𝑛 πœ” 0
𝐢𝑛
∞
𝑗𝑛 πœ” 0 𝑑
=
𝑛=βˆ’βˆž (𝑗𝑛 πœ” ) 𝑒
0
𝐢𝑛
𝑗𝑛 πœ” 0
𝐢𝑛 [
𝐢
FS-1[𝑗𝑛 πœ”π‘› ]
0
Proved.
, C0 = 0
7) Convolution theorem or property
The convolution theorem or property states that , The Fourier series of the
convolution of two time domain functions x1(t) and x2(t) is equal to the multiplication of their
Fourier series coefficients, i.e.”Convolution of two functions in time domain is equivalent to
multiplication of their Fourier coefficients in frequency domain”.
The convolution property states that,if
𝐹𝑠
x1(t) Cn and x2 (𝑑)
𝐹𝑠
𝐹𝑠
Dn
then
x1(t)*x2(t) TCnDn
proof: From the definition of Fourier series, we have
1 𝑑 +𝑇
1 𝑇
FS[x1(t)*x2(t)]=𝑇 𝑑 0 [x1 (t) βˆ— x2 (t)] 𝑒 βˆ’π‘—π‘› πœ” 0 𝑑 𝑑𝑑 =𝑇 0 [x1 (t) βˆ— x2 (t)] 𝑒 βˆ’π‘—π‘› πœ” 0 𝑑 𝑑𝑑
0
But from the convolution integral for a periodic signal, we have
𝑇
𝑇
x1 (t) βˆ— x2 (t)= 0 x1 (Ο„) βˆ— x2 (t βˆ’ Ο„) π‘‘πœ or x1 (t) βˆ— x2 (t)= 0 x1 (t βˆ’ Ο„) βˆ— x2 (Ο„) π‘‘πœ
Substituting back, we get
1 𝑇 𝑇
FS[x1(t)*x2(t)] =𝑇 0 ( 0 x1 (Ο„) βˆ— x2 (t βˆ’ Ο„) π‘‘πœ) 𝑒 βˆ’π‘—π‘› πœ” 0 𝑑 𝑑𝑑
Interchanging the order of integration, we get
𝑇
1 𝑇
FS[x1(t)*x2(t)] = 0 x1 (Ο„)[ 0 x2 (t βˆ’ Ο„) 𝑒 βˆ’π‘—π‘› πœ” 0 𝑑 𝑑𝑑] π‘‘πœ
𝑇
Substituting t-𝛕 =t' in RHS, we have dt=dt'. Substituting back, we get
β€²
π‘‡βˆ’Ο„
1 𝑇
FS[x1(t)*x2(t)] =𝑇 0 x1 (Ο„)[ βˆ’Ο„ x2 (tβ€²) 𝑒 βˆ’π‘—π‘› πœ” 0 (𝑑 +Ο„) 𝑑𝑑′] π‘‘πœ
1
𝑇
1
𝑇
=T(𝑇 0 x1 (Ο„)𝑒 βˆ’π‘—π‘› πœ” 0 𝑑 𝑑𝛕) (𝑇 0 x1 (tβ€²)𝑒 βˆ’π‘—π‘› πœ” 0 𝑑′ 𝑑tβ€²)
=T(Cn)(Dn)
0 to T or –𝛕 to T-𝛕 or t0 to t0+T will have the same period.
𝐹𝑠
x1(t)*x2(t) TCnDn
Proved.
8) Modulation or Multiplication property
The Modulation or Multiplication property states that,if
𝐹𝑠
x1(t) Cn and x2 (𝑑)
then
proof:
𝐹𝑠
Dn
𝐹𝑠
∞
x1(t)x2(t)
𝑙=βˆ’βˆž 𝐢𝑙 π·π‘›βˆ’π‘™
From the definition of Fourier series, we have
1 𝑑 +𝑇
FS[x1(t)x2(t)]= 𝑇 𝑑 0 x1 t x2 t 𝑒 βˆ’π‘—π‘› πœ” 0 𝑑 𝑑𝑑
0
𝑑 +𝑇
=𝑇 𝑑 0 x1
0
1 𝑑 0 +𝑇
=𝑇 𝑑
x1
0
1
t (
∞
𝑗𝑛 πœ” 0 𝑑
) 𝑒 βˆ’π‘—π‘› πœ” 0 𝑑 𝑑𝑑
𝑙=βˆ’βˆž 𝐢𝑙 𝑒
βˆ’π‘— (π‘›βˆ’π‘™)πœ” 0 𝑑
t ( ∞
) 𝑑𝑑
𝑙=βˆ’βˆž 𝐢𝑙 𝑒
Interchanging the order of integration and summation, we have
1 𝑑 0 +𝑇
FS[x1(t)x2(t)]= ∞
x1 t 𝑒 βˆ’π‘— (π‘›βˆ’π‘™)πœ” 0 𝑑 𝑑𝑑) = ∞
𝑙=βˆ’βˆž 𝐢𝑙 (𝑇 𝑑
𝑙=βˆ’βˆž 𝐢𝑙 π·π‘›βˆ’π‘™
0
x1(t)x2(t)
𝐹𝑠
Proved.
∞
𝑙=βˆ’βˆž 𝐢𝑙 π·π‘›βˆ’π‘™
9) Conjugation and Conjugate symmetry property
𝐹𝑠
If
Then the Conjugation property states that
x*(t)
𝐹𝑠
x(t) Cn
βˆ—
πΆβˆ’π‘›
[for complex x(t)]
proof: Conjugation property
From the definition of Fourier series, we have
1 𝑑 +𝑇
1
FS[x*(t)]= 𝑇 𝑑 0 x βˆ— (t) 𝑒 βˆ’π‘—π‘› πœ” 0 𝑑 𝑑𝑑=(𝑇
1
=(𝑇
x*(t)
𝐹𝑠
0
𝑑 0 +𝑇
x(t) 𝑒 𝑗𝑛 πœ” 0 𝑑 𝑑𝑑)*
𝑑0
𝑑 0 +𝑇
x(t) 𝑒 βˆ’π‘— (βˆ’π‘›)πœ” 0 𝑑 𝑑𝑑)*=[C-n]*
𝑑0
βˆ—
πΆβˆ’π‘›
Proved.
10) Parseval’s Relation or Theorem or property
𝐹𝑠
If
x1(t) Cn
𝐹𝑠
[for complex x1(t)]
And
x2 (𝑑) Dn [for complex x2(t)]
Then, the parsevals relation states that
1 𝑑 0 +𝑇
βˆ—
x1 t x2βˆ— t 𝑑𝑑= ∞
𝑛=βˆ’βˆž Cn 𝐷𝑛 [for complex x1(t) and x2(t)]
𝑇 𝑑0
And parseval’s identity states that
1 𝑑 0 +𝑇
2
⎹x(t)⎸2 𝑑𝑑 = ∞
if x1(t) = x2(t) = x(t)
𝑛=βˆ’βˆž ⎹Cn ⎸
𝑇 𝑑0
Proof: parsevals relation
1 𝑑 +𝑇
1 𝑑 +𝑇
𝑗𝑛 πœ” 0 𝑑
LHS =𝑇 𝑑 0 x1 t x2βˆ— t 𝑑𝑑 = 𝑇 𝑑 0 ( ∞
) π‘₯2βˆ— 𝑑𝑑
𝑛=βˆ’βˆž 𝐢𝑛 𝑒
0
0
Interchanging the order of integration and summation in the RHS , we have
1 𝑑 0 +𝑇
1 𝑑 0 +𝑇 βˆ—
x1 t x2βˆ— t 𝑑𝑑 = ∞
x2 t 𝑒 𝑗𝑛 πœ” 0 𝑑 𝑑𝑑)
𝑛=βˆ’βˆž 𝐢𝑛 (𝑇 𝑑
𝑇 𝑑
0
0
=
∞
𝑛=βˆ’βˆž
1
𝐢𝑛 (𝑇
𝑑 0 +𝑇
x2
𝑑0
t 𝑒 βˆ’π‘—π‘› πœ” 0 𝑑 𝑑𝑑)*=
∞
𝑛=βˆ’βˆž Cn
[Dn ]βˆ—
Parseval’s identity
𝑑 0 +𝑇
𝑇 𝑑0
1
Proved
x1 t
x2βˆ—
t 𝑑𝑑=
∞
𝑛=βˆ’βˆž Cn
π·π‘›βˆ—
If x1(t) = x2(t) = x(t), then the relation changes to
1 𝑑 0 +𝑇
𝑇 𝑑0
x t x βˆ— t 𝑑𝑑=
∞
𝑛=βˆ’βˆž Cn
πΆπ‘›βˆ—
Since ⎹x(t)⎸2 = x t x βˆ— t and ⎹Cn ⎸2 = CnπΆπ‘›βˆ— , substituting these values in above equation, we get
1 𝑑 0 +𝑇
⎹x(t)⎸2
𝑇 𝑑0
𝑑𝑑 =
∞
2
𝑛=βˆ’βˆž ⎹Cn ⎸
Proved.
Complex Fourier Spectrum
The Fourier spectrum of a periodic signal x(t) is a plot of its Fourier coefficients versus
frequency Ο‰. It is in two parts : (a) Amplitude spectrum and (b) phase spectrum. The plot of the
amplitude of Fourier coefficients verses frequency is known as the amplitude spectra, and the
plot of the phase of Fourier coefficients verses frequency is known as phase spectra. The two
plots together are known as Fourier frequency spectra of x(t).This type of representation is also
called frequency domain representation. The Fourier spectrum exists only at discrete frequencies
nΟ‰o, where n=0,1,2,….. Hence it is known as discrete spectrum or line spectrum. The envelope
of the spectrum depends only upon the pulse shape, but not upon the period of repetition.
The trigonometric representation of a periodic signal x(t) contains both sine and cosine
terms with positive and negative amplitude coefficients (an and bn) but with no phase angles.
The cosine representation of a periodic signal contains only positive amplitude
coefficients with phase angle ΞΈn. Therefore, we can plot amplitude spectra(An versus Ο‰)and phase
spectra (ΞΈn versus Ο‰). Since, in this representation, Fourier coefficients exist only for positive
frequencies, this spectra is called single-sided spectra.
The exponential representation of a periodic signal x(t) contains amplitude coefficients
Cn Which are complex. Hence, they can be represented by magnitude and phase. Therefore, we
can plot two spectra, the magnitude spectrum(⎹Cn⎸versus Ο‰) and phase spectrum ( 𝐢𝑛 versus Ο‰).
The spectra can be plotted for both positive and negative frequencies.Henceit is called two-sided
spectra.
The below figure (a) represents the spectrum of a trigonometric Fourier series extending
from 0 to∞,producing a one-sided spectrum as no negative frequencies exist here. The figure (b)
represents the spectrum of a complex exponential Fourier series extending from βˆžπ‘‘π‘œβˆž,producing a two-sided spectrum.
The amplitude spectrum of the exponential Fourier series is symmetrical Fourier series is
symmetrical about the vertical axis. This is true for all periodic functions.
Fig: Complex frequency spectrum for (a) Trigonometric Fourier series and (b) complex
exponential Fourier series.
If Cn is a general complex number, then
Cn = ⎹CnβŽΈπ‘’ 𝑗 πœƒπ‘›
C-n = ⎹CnβŽΈπ‘’ βˆ’π‘— πœƒπ‘›
Cn = ⎹C-n⎸
The magnitude spectrum is symmetrical about the vertical axis passing
through the origin, and the phase spectrum is antisymmetrical about the vertical axis passing
through the origin. So the magnitude spectrum exhibits even symmetry and phase spectrum
exhibits odd symmetry.
When x(t) is real , then C-n=πΆπ‘›βˆ— , the complex conjugate of Cn.
Solved Problems:
Problem 1: Show that the magnitude spectrum of every periodic function is symmetrical about
the vertical axis passing through the origin and the phase spectrum is antisymmetrical about the
vertical axis passing through the origin.
Solution:
The coefficient Cn of exponential Fourier series is given by
1
Cn = 𝑇
T
2
𝑇
βˆ’
2
1
And
C-n = 𝑇
π‘₯ 𝑑 𝑒 βˆ’π‘—π‘› πœ” 0 𝑑 𝑑𝑑)
T
2
𝑇
βˆ’
2
π‘₯ 𝑑 𝑒 𝑗𝑛 πœ” 0 𝑑 𝑑𝑑)
It is evident from these equations that the coefficients Cn and C-n are complex conjugate of each
other , that is
Cn= πΆπ‘›βˆ—
Hence
⎹Cn ⎸=⎹Cβˆ’n ⎸
It, therefore, follows that the magnitude spectrum is symmetrical about the vertical axis
passing through the origin, and hence is an even function of Ο‰.
If Cn is real, then C-n is also real and Cn is equal to C-n. is complex, let
Cn=⎹Cn ⎸ejθ n
C-n=⎹Cn ⎸eβˆ’jΞΈ n
then
The phase of Cn is ΞΈn ;however , the phase of C-n isβˆ’ΞΈn .Hence , it is obvious that the phase
spectrum is antisymmetrical, and the magnitude spectrum is symmetrical about the vertical axis
passing through the origin.
Problem 2: With regard to Fourier series representation, justify the following statement: odd
functions have only sine terms.
Solution: We know that the trigonometric Fourier series of a periodic function x(t) in any
𝑇
𝑇
interval t0 to t0+T or 0 to T or -2 to 2 is given by
∞
𝑛=1(π‘Žπ‘›
x(t)= π‘Ž0 +
T
2
𝑇
βˆ’
2
2
π‘Žπ‘› = 𝑇
and𝑏𝑛 =
T
2
𝑇
βˆ’
2
1
π‘Ž0 =
where
𝑇
cos π‘›πœ”0 𝑑 + 𝑏𝑛 sin 𝑛ω0 t)
π‘₯ 𝑑 𝑑𝑑
π‘₯ 𝑑 cos π‘›πœ”0 𝑑 𝑑𝑑
T
2
𝑇
βˆ’
2
2
𝑇
π‘₯ 𝑑 sin 𝑛ω0 t 𝑑𝑑
For an odd function
x(t) = -x(-t)
Also even part is zero, i.e. xe(t) = 0 and x(t) = x0(t)
1
π‘Ž0 =
𝑇
T
2
𝑇
βˆ’
2
π‘₯ 𝑑 𝑑𝑑 =
1
𝑇
T
2
𝑇
βˆ’
2
π‘₯0 𝑑 𝑑𝑑 = 0
Since the integration of an odd function,over one cycle is always zero.
2
π‘Žπ‘› = 𝑇
T
2
𝑇
βˆ’
2
π‘₯ 𝑑 cos π‘›πœ”0 𝑑 𝑑𝑑 =
2
𝑇
T
2
𝑇
βˆ’
2
π‘₯0 𝑑 cos π‘›πœ”0 𝑑 𝑑𝑑
Here π‘₯0 (t) is odd and cos π‘›πœ”0 𝑑 is even. So π‘₯0 𝑑 cos π‘›πœ”0 𝑑, i.e. the product if an odd function
and even function is odd. So integration over one cycle is zero. Therefore,
π‘Žπ‘› = 0
𝑏𝑛 =
Now,
2
𝑇
T
2
𝑇
βˆ’
2
π‘₯ 𝑑 sin 𝑛ω0 t 𝑑𝑑 =
2
𝑇
T
2
𝑇
βˆ’
2
π‘₯0 𝑑 sin 𝑛ω0 t 𝑑𝑑
Now, π‘₯0 (t) is odd sin 𝑛ω0 t is also odd. Therefore, the product π‘₯0 (t)sin 𝑛ω0 t is even. So we
have to evaluate the integral:
𝑏𝑛 =
2
𝑇
T
2
𝑇
βˆ’
2
π‘₯ 𝑑 sin 𝑛ω0 t 𝑑𝑑 =
x(t)=
4
T
2
𝑇 0
∞
𝑛=1 𝑏𝑛
π‘₯ 𝑑 sin 𝑛ω0 t 𝑑𝑑
sin 𝑛ω0 t
Thus, the Fourier series of odd functions contains only sine terms.
Problem3: With regard to Fourier series representation, justify the following statement: Even
functions have no sine terms.
Solution: We know that the trigonometric Fourier series of a periodic function x(t) in any
𝑇
𝑇
interval t0 to t0+T or 0 to T or -2 to 2 is given by
x(t)= π‘Ž0 +
1
π‘Ž0 =
where
T
2
𝑇
βˆ’
2
2
π‘Žπ‘› = 𝑇
𝑇
∞
𝑛=1(π‘Žπ‘›
T
2
𝑇
βˆ’
2
cos π‘›πœ”0 𝑑 + 𝑏𝑛 sin 𝑛ω0 t)
π‘₯ 𝑑 𝑑𝑑
π‘₯ 𝑑 cos π‘›πœ”0 𝑑 𝑑𝑑
𝑏𝑛 =
and
T
2
𝑇
βˆ’
2
2
𝑇
π‘₯ 𝑑 sin 𝑛ω0 t 𝑑𝑑
For an even function
x(t) = x(-t)
Also odd part is zero, i.e. x(t) = xe(t)
π‘Ž0 =
π‘Žπ‘› =
T
2
𝑇 βˆ’π‘‡
2
T
2 2
𝑇 βˆ’π‘‡
2
1
π‘₯ 𝑑 𝑑𝑑 =
1
𝑇
T
2
𝑇
βˆ’
2
2
π‘₯𝑒 𝑑 𝑑𝑑 = 𝑇
π‘₯ 𝑑 cos π‘›πœ”0 𝑑 𝑑𝑑 =
2
𝑇
T
2
𝑇
βˆ’
2
T
2
0
π‘₯ 𝑑 𝑑𝑑
π‘₯𝑒 𝑑 cos π‘›πœ”0 𝑑 𝑑𝑑
π‘₯𝑒 (t) is even and also cos π‘›πœ”0 𝑑 is even. So, the product π‘₯𝑒 (t)cos π‘›πœ”0 𝑑 is even. So we have to
evaluate the integral:
2
π‘Žπ‘› = 𝑇
𝑏𝑛 =
2
𝑇
T
2
𝑇
βˆ’
2
T
2
𝑇
βˆ’
2
4
π‘₯𝑒 𝑑 cos π‘›πœ”0 𝑑 𝑑𝑑 = 𝑇
π‘₯ 𝑑 sin 𝑛ω0 t 𝑑𝑑 =
2
𝑇
T
2
0
T
2
𝑇
βˆ’
2
π‘₯ 𝑑 cos π‘›πœ”0 𝑑 𝑑𝑑
π‘₯𝑒 𝑑 sin 𝑛ω0 t 𝑑𝑑
Here π‘₯𝑒 (t) is even sin 𝑛ω0 tand is odd. So, the product π‘₯𝑒 𝑑 sin 𝑛ω0 t is odd. Thus, the integral
over one complete cycle is zero.
𝑏𝑛 = 0
x(t)= π‘Ž0 + ∞
𝑛=1 π‘Žπ‘› cos π‘›πœ”0 𝑑
Thus, the Fourier series of even functions contains no sine terms.
Problem 4: With regard to Fourier series representation, justify the following statement:
Functions with half wave symmetry have only odd harmonics.
Solution: We know that the trigonometric Fourier series of a periodic function x(t) in any
𝑇
𝑇
interval t0 to t0+T or 0 to T or -2 to 2 is given by
x(t)= π‘Ž0 +
∞
𝑛=1(π‘Žπ‘›
cos π‘›πœ”0 𝑑 + 𝑏𝑛 sin 𝑛ω0 t)
π‘Ž0 =
where
2
π‘Žπ‘› = 𝑇
T
2
𝑇
βˆ’
2
T
2
𝑇
βˆ’
2
1
𝑇
π‘₯ 𝑑 𝑑𝑑
π‘₯ 𝑑 cos π‘›πœ”0 𝑑 𝑑𝑑
T
2
𝑇
βˆ’
2
2
𝑏𝑛 =
and
𝑇
π‘₯ 𝑑 sin 𝑛ω0 t 𝑑𝑑
For a half wave symmetric function,
𝑇
x(t) = -x(t±2 )
π‘Ž0 = 0
Therefore,
2
π‘Žπ‘› = 𝑇
T
π‘₯
0
𝑑 cos π‘›πœ”0 𝑑 𝑑𝑑
T
2
2
=𝑇[
0
T
π‘₯ 𝑑 cos π‘›πœ”0 𝑑 𝑑𝑑 +
𝑇
2
π‘₯ 𝑑 cos π‘›πœ”0 𝑑 𝑑𝑑 ]
To have the limits of second integration also from 0 to T/2 , change the variable t by t + (T/2) in
the second integration.
T
2
2
π‘Žπ‘› =𝑇 [
0
T
2
π‘₯ 𝑑 cos π‘›πœ”0 𝑑 𝑑𝑑 +
T
T
π‘₯ 𝑑 + 2 cos π‘›πœ”0 (𝑑 + 2 ) 𝑑𝑑 ]
0
𝑇
But
x(t) = -x(t±2 )
T
2
2
π‘Žπ‘› = 𝑇 [
π‘₯ 𝑑 cos π‘›πœ”0 𝑑 𝑑𝑑 +
0
T
2
βˆ’π‘₯ 𝑑 cos(π‘›πœ”0 𝑑 + 𝑛π) 𝑑𝑑 ]
0
2
=𝑇[
0
π‘Žπ‘› =
T
2
0
π‘₯ 𝑑 cos π‘›πœ”0 𝑑 βˆ’ π‘₯ 𝑑 cos(π‘›πœ”0 𝑑 + 𝑛π)] 𝑑𝑑
π‘“π‘œπ‘Ÿ 𝑒𝑣𝑒𝑛 𝑛
T
2
4
𝑇 0
π‘₯ 𝑑 cos π‘›πœ”0 𝑑
π‘“π‘œπ‘Ÿ π‘œπ‘‘π‘‘ 𝑛
Now,
𝑏𝑛 =
2
T
2
π‘₯ 𝑑 sin 𝑛ω0 t 𝑑𝑑
𝑇 0
T
2
2
=𝑇[
0
π‘₯ 𝑑 sin 𝑛ω0 t 𝑑𝑑 +
T
𝑇
2
π‘₯ 𝑑 sin 𝑛ω0 t 𝑑𝑑]
To have the limits of second integration also from 0 to T/2 change the variable t by t + (T/2) in
the second integration.
2
𝑏𝑛 = 𝑇 [
But
T
2
0
T
2
π‘₯ 𝑑 sin 𝑛ω0 t 𝑑𝑑 +
0
𝑇
x(t) = -x(t±2 )
𝑇
T
π‘₯ 𝑑 + 2 sin 𝑛ω0 (t + 2 ) 𝑑𝑑]
2
𝑏𝑛 = 𝑇 [
T
2
0
π‘₯ 𝑑 sin 𝑛ω0 t 𝑑𝑑 +
2
=𝑇[
T
2
0
0
𝑏𝑛 =
4
T
2
0
βˆ’π‘₯ 𝑑 sin( 𝑛ω0 t + nΟ€) 𝑑𝑑]
π‘₯ 𝑑 sin 𝑛ω0 t 𝑑𝑑 -
T
2
0
π‘₯ 𝑑 sin 𝑛ω0 t cos nΟ€) 𝑑𝑑]
π‘“π‘œπ‘Ÿ 𝑒𝑣𝑒𝑛 𝑛
T
2
π‘₯ 𝑑 sin 𝑛ω0 t
π‘“π‘œπ‘Ÿ π‘œπ‘‘π‘‘ 𝑛
Thus, the Fourier series of half wave symmetric function consists of only odd harmonics.
𝑇 0
Problem 5: Draw the complex fourier spectrum of given figure.
Solution: It is a plot between Cn and nω0 = ω. As Cn is real, only amplitude plot is sufficient.
The amplitude versus frequency plot, called the amplitude spectrum.
To draw the frequency spectrum,
x(t) =
2𝐴
πœ‹
+
A βˆ’1 n
∞
𝑛=βˆ’βˆž [
𝑛≠0 Ο€ 2n+1
exponential Fourier series coefficients
2𝐴
C0 = πœ‹
2𝐴
C1= 3πœ‹
2𝐴
C2= 15πœ‹
2𝐴
C3= 35πœ‹
+
βˆ’1 n +1
2nβˆ’1
]ej2nt
Problem 6: For the periodic gate function shown in below figure, plot the magnitude and phase
spectra.
Solution:
Amplitude and phase spectra are shown in below figures respectively.
Example 7: Find the complex exponential Fourier series and the trigonometric Fourier series of
unit impulse train Ξ΄T (t) shown in bellow figure.
Solution:
Let
The periodic waveform shown in above figure with period T can be expressed as
x(t) = δT (t) = ∞
𝑛=βˆ’βˆž 𝛿(𝑑 βˆ’ 𝑛𝑇)
𝑇
t0= - 2
𝑇
So
t0+T = - 2
In one period, only 𝛿(t) exists.
Fundamental frequency Ο‰0 = (2πœ‹/T)
1
Cn =
𝑇
1
=𝑇
𝑑 0 +𝑇
π‘₯(𝑑)𝑒 βˆ’π‘—π‘› πœ” 0 𝑑 𝑑𝑑
𝑑0
T
2
𝑇
βˆ’
2
1
𝛿(t) 𝑒 βˆ’π‘—π‘› πœ” 0 𝑑 𝑑𝑑
𝛿(t)= 1
=0
at t = 0
elsewhere
1
= 𝑇 e0 = 𝑇
1
Cn = 𝑇
The exponential Fourier series is
x(t)=
∞
𝑛=βˆ’βˆž 𝐢𝑛
𝑒 𝑗𝑛 πœ” 0 𝑑 =
1 𝑗𝑛 (2Ο€ )𝑑
∞
T
𝑛=βˆ’βˆž 𝑇 𝑒
1
=𝑇
∞
𝑛=βˆ’βˆž
𝑒 𝑗𝑛
2Ο€
T
𝑑
To get the trigonometric Fourier coeifficients , we have
1
π‘Ž0 = C0 = 𝑇
2
π‘Žπ‘› = 𝑇
T
π‘₯ 𝑑 cos π‘›πœ”0 𝑑 𝑑𝑑
0
2 T
π‘₯ 𝑑 sin 𝑛ω0 t 𝑑𝑑
𝑇 0
=
T
Ξ΄
𝑇 0
2 T
π‘₯
𝑇 0
2
2
𝑑 cos π‘›πœ”0 𝑑 𝑑𝑑 = 𝑇
𝑏𝑛 =
=
𝑑 sin 𝑛ω0 t 𝑑𝑑 = 0
Therefore , the trigonometric Fourier series is
x(t) +π‘Ž0 + ∞
𝑛=1 π‘Žπ‘› cos π‘›πœ”0 𝑑 + 𝑏𝑛 sin 𝑛ω0 t =
The complex Fourier series coefficients are
1
Cn = 𝑇
for all n
The magnitude and phase spectrum are
1
⎹Cn ⎸= βŽΉπ‘‡ ⎸for all n
𝐢𝑛 = 0o
for all n
The frequency spectra are plotted in bellow figure.
1
𝑇
2
+𝑇
∞
𝑛=1 cos π‘›πœ”0 𝑑
Problem 8: Find the Fourier series of the signal x(t)=e-t with T = 1 sec as shown in bellow
figure. Draw its magnitude and phase spectra.
x(t)=e-t
Period T = 1 sec
2πœ‹
2πœ‹
Fundamental frequency Ο‰0 = 𝑇 = 1 =2πœ‹
The exponential Fourier series is
Solution:
Given signal
1
Cn = 𝑇
=
=
x(t)=
𝑇
π‘₯(𝑑)𝑒 βˆ’π‘—π‘› πœ” 0 𝑑
0
1
𝑑𝑑 = 1
1 βˆ’t βˆ’π‘—π‘› 2πœ‹π‘‘
e 𝑒
0
𝑑𝑑
1 βˆ’(1+𝑗𝑛 2πœ‹)𝑑
𝑒 βˆ’(1+𝑗𝑛 2πœ‹ )𝑑 1
𝑒
𝑑𝑑
=
[
]
0
βˆ’(1+𝑗𝑛 2πœ‹) 0
𝑒 βˆ’(1+𝑗𝑛 2πœ‹ ) βˆ’1
1βˆ’π‘’ βˆ’1 𝑒 βˆ’π‘—π‘› 2πœ‹
1βˆ’π‘’ βˆ’1
βˆ’(1+𝑗𝑛 2πœ‹)
∞
𝑛=βˆ’βˆž 𝐢𝑛
=
𝑒 𝑗𝑛 πœ” 0 𝑑 =
=
(1+𝑗𝑛 2πœ‹)
(1+𝑗𝑛 2πœ‹)
1βˆ’π‘’ βˆ’1
∞
𝑗𝑛 2πœ‹π‘‘
𝑛=βˆ’βˆž (1+𝑗𝑛 2πœ‹ ) 𝑒
The spectra are shown in bellow figure.
Problem 9: Obtain the exponential Fourier Series for the wave form shown in below figure.
Also draw the frequency spectrum.
Solution: The periodic waveform shown in fig with a period T= 2Ο€ can be expressed as:
A
0≀t≀π
x(t) =
βˆ’A
Ο€ ≀ t ≀ 2Ο€
Let
t 0 = 0, t 0 +T = 2Ο€
2Ο€
2Ο€
and
Fundamental frequency ⍡0 = T = 2Ο€ =1
Exponential Fourier series
1
C0 = T
Cn =
T
T 0
= 2Ο€
x(t)e
=
Cn =
1
+ 2Ο€
2Ο€
Ο€
βˆ’A dt = 0
dt
Ο€
2Ο€ 0
A
1
=βˆ’
∴
x(t) dt
Ο€
A dt
0
βˆ’jn ⍡0 t
1
1
T
0
1
Aeβˆ’jnt dt + 2Ο€
βˆ’Aeβˆ’jnt dt
(βˆ’1)n βˆ’ 1 βˆ’ 1 βˆ’ (βˆ’1)n = βˆ’j
j2nΟ€
2A
βˆ’j Ο€n
A
2nΟ€
for odd n
0
x(t) = C0 +
2Ο€
Ο€
for even n
∞
jn ⍡0 t
n=βˆ’βˆž Cn e
n≠0
2𝐴
=
A
∞
jnt
n=βˆ’βˆž βˆ’j 2nΟ€ e
2𝐴
2𝐴
C0 = 0, C-1 = C1 = πœ‹ , πΆβˆ’3 = C3 = 3πœ‹ , πΆβˆ’5 = C5= 5πœ‹
The frequency spectrum is shown in the below figure.
Problem10: Find the exponential Fourier series and plot the frequency spectrum for the full
wave rectified sine wave given in below figure
Solution: The waveform shown in fig can be expressed over one period(0 to Ο€) as:
2Ο€
x(t) = A sin ⍡t where ⍡= 2Ο€ =1
because it is part of a sine wave with period = 2Ο€
x(t) = A sin ⍡t 0≀ t ≀ Ο€
The full wave rectified sine wave is periodic with period T = Ο€
Let
t 0 = 0, t 0 +T = 0+Ο€ = Ο€
2Ο€
2Ο€
and
Fundamental frequency ⍡0 = T = Ο€ = 2
The exponential Fourier series is
jn ⍡0 t
j2nt
x(t) = ∞
= ∞
n=βˆ’βˆž Cn e
n=βˆ’βˆž Cn e
1 T
where
Cn = T 0 x(t)eβˆ’jn ⍡0 t dt
Ο€
A Ο€
A sin t eβˆ’j2nt dt = Ο€ 0 sin t eβˆ’j2nt
0
A e j 1βˆ’2n t βˆ’e 0
e βˆ’j 1βˆ’2n t βˆ’e 0
1
=Ο€
= j2Ο€
j 1βˆ’2n
βˆ’
dt
βˆ’j 1βˆ’2n
2A
∴Cn = Ο€(1βˆ’4n 2 )
T
0
1 Ο€
Ο€ 0
1
C0 = T
x(t) dt
A
=
A sin t dt = Ο€ βˆ’ cos 𝑑 πœ‹0 =
The exponential Fourier series is given by
2A
Ο€
x(t) =
2A
∞
j2nt 2A
= Ο€
n=βˆ’βˆž Ο€(1βˆ’4n 2 ) e
2𝐴
2𝐴
2A
+Ο€
∞
n=βˆ’βˆž
n≠0
2𝐴
e j2nt
1βˆ’4n 2
2𝐴
C0 = πœ‹ , C-1 = C1 = 3πœ‹ , πΆβˆ’2 = C2 = βˆ’ 15πœ‹ , πΆβˆ’3 = C3 = βˆ’ 35πœ‹
The frequency spectrum is shown in the below figure.
Assignment Problems:
1. Find the trigonometric Fourier series for the periodic signal x(t) shown in figure below. Also
draw the frequency spectrum.
2. Find the trigonometric Fourier series for the periodic signal x(t) shown in figure below. Also
draw the frequency spectrum.
3. Compute the exponential Fourier series of the signal shown in figure below. Also draw the
frequency spectrum.
4. Find the exponential series of the signal shown in the figure below. Also draw the frequency
spectrum.
5. Find the trigonometric Fourier series of x(t) =t2 over the interval (-1, 1).
6. Find the Fourier series for the periodic signal x(t) =t2 for 0 ≀ 𝑑 ≀ 1, so that it repeats every 1
sec.
7. Determine the time signal corresponding to the magnitude and phase spectra shown in below
figure with ⍡o = πœ‹.
Figure (a) Magnitude and (b) Phase spectra
8. Determine the Fourier series representation for x(t) = 2 sin(2πœ‹π‘‘ βˆ’ 3)+sin6πœ‹π‘‘.
9. Find the complex exponential Fourier series coefficient of the signal x(t) = 3cos 4⍡0 t. Also
draw the frequency spectrum.
10. Find the complex exponential Fourier series coefficient of the signal x(t) = 2cos 3⍡0 t. Also
draw the frequency spectrum.
Simulation:
MATLAB can be used to find and plot the exponential Fourier series of a periodic function.
Using the symbolic math toolbox, we can perform the integrations necessary to find the
coefficients Fn.
Exponential Fourier Series
For the waveform, f(t), shown below:
we can write the exponential Fourier Series expansion as:
∞
𝐹𝑛 𝑒 𝑗𝑛 πœ” π‘œ 𝑑
𝑓 𝑑 =
𝑛 =βˆ’βˆž
Where
The Fn are complex coefficients. Thus, when we plot them we need to plot the magnitude
andphase separately. In Matlab we use the β€œabs” function to find the magnitude and β€œangle” to
findthe phase angle. The following Matlab statements show one method for obtaining these
plots:
Program
clc;clear all; close all;
n=-10:10;
Fn=(1./(j*n*2*pi)).*(4*exp(j*n*2*pi/7)+ 2*exp(-j*n*12*pi/7)-6*exp(j*n*6*pi/7));
Fn(find(n==0))=10/7;
stem(n,abs(Fn))
xlabel('n-->');
ylabel('|Fn|-->');
title('Magnitude spectrum');
figure,
stem(n,angle(Fn))
xlabel('n-->');
ylabel('L Fn -->');
title('Phase spectrum');
Since the value of Fn is different at n = 0, we assign that value using the find command.
Thestatement: find(n==0) will locate the index of the n vector that contains a zero. In this case
thatindex is 6, since the zero is found in the sixth position in n. That value is then used to index
Fn,so the value of Fn(6) is set to 10/7. This corresponds to the value for n = 0. These
commandsresult in the following magnitude and phase plots:
Note that the magnitude spectrum is even and the phase spectrum is odd as expected.
References:
[1] Alan V.Oppenheim, Alan S.Willsky and S.Hamind Nawab, β€œSignals & Systems”, Second
edition, Pearson Education, 8th Indian Reprint, 2005.
[2] M.J.Roberts, β€œSignals and Systems, Analysis using Transform methods and MATLAB”,
Second edition,McGraw-Hill Education,2011
[3] John R Buck, Michael M Daniel and Andrew C.Singer, β€œComputer explorations in Signals
and Systems using MATLAB”,Prentice Hall Signal Processing Series
[4] P Ramakrishna rao, β€œSignals and Systems”, Tata McGraw-Hill, 2008
[5] Tarun Kumar Rawat, β€œSignals and Systems”, Oxford University Press,2011
[6] A.Anand Kumar, β€œSignals and Systems” , PHI Learning Private Limited ,2011