Unit 2
Unit 2
FOURIER TRANSFORM
2.1.1.Introduction:
We can represent any periodic function x (t), over the entire interval (-∞<t<∞) as
a discrete sum of exponential functions. We can also represent any non-periodic
function in terms of exponential functions over any finite interval ( . For
the purpose of frequency analysis, we need to represent every type of driving function
in terms of exponential functions not over the finite interval but over the entire interval
(-∞, ∞).
Let us consider a function x(t) as shown in fig (1). We derive to represent this function
as a sum of exponential functions over the entire interval (-∞, ∞). For the purpose we
construct a new periodic function with period T as shown in fig (2). The period T
is made large enough so that there is no overlap between the pulses of the shape of x(t).
This new function is a periodic function and consequently can be represented with
an exponential fourier series.
Fig: 1 Fig: 2
In the limit if we let T becomes infinite, then the pulses in the periodic function repeat
after an infinite interval. Hence in the limit T→∞, and x(t) are identical i.e.
Thus the fourier series representing over the entire interval will also represent x(t)
over the entire interval if we let T = ∞ in this series.
The term represents the amplitude of the component of frequency . We shall now
assume that T becomes very large. As T becomes larger, (the fundamental frequency)
becomes smaller and the spectrum becomes denser.
As seen from equation (1) when T→∞, the magnitude of each component becomes
infinite similarly small but now there are also an infinite number of frequency
components. The spectrum exists for every value of ω, and is no longer discrete but a
continuous function of ω.
To illustrate this,
T ----------------------> (3)
Then
= ---------------------> (4)
∴ ------------------------> (5)
Let us try a graphical interpretation of equation (6), which represents a discrete sum or
a sum of discrete frequency components. Actually the quantity is complex
in general. However we shall assume that the quantity is real. The figure (3)
shows such a plot of this quantity as a function of ω. This function only at discrete
values of ω at ω1, ω2, ω3, _ _ _ _ _ _, , _ _ _ _ _ etc. Where
Fig: 3
The curve now is a continuous function of ω and is given by . Also at T→∞ the
function and equation (6) and (5) becomes.
Equation (7) and (8) are usually referred as the Fourier Transform pair. Equation (8) is
known as the direct Fourier Transform of x(t) and equation (7) is known as Inverse
Fourier transform of X(ω).
= t>0
x (t) =
0 t<0
1 t≥0
Because u(t) =
0 t<0
=
=
= =
F[ =
Now = =
X(ω) = =j
θ (ω) =
=-
t>0
x(t) =
t<0
FT [x(t)] = X(ω) =
= t
= t
= +
X(ω) =
∴ and θ(ω) = 0
The signal and its spectrum is shown in the fig 5 (a) and 5(b) respectively.
2.1.3.3.
∴ FT [x(t)] = X(ω) =
= t
X(ω) =
a > 0. X [t
∴ =t
The spectral density function X(ω) may be expressed in terms of magnitude and phase-
density functions.
X(ω) = =
2.1.3.4. A Gate function
0 |t|>
Fig: 6(a)
X(ω) =
=A
=A -------------> (16)
Fig: 6(b)
2.1.4.Existence of Fourier Transform:
x(t) =
and X( ) =
For the function x(t) to be Fourier transformable, x(t) should satisfy Dirichlet‟s
conditions given below…
( i.)The function x(t) should be a single valued with a finite number of maxima and
minima and a finite number of discontinuities in any finite time interval
a.) The original time function can be uniquely recovered from it.
b.) It has a property analogous to common logarithm that helps in evaluating the
convolution integral.
2.1.5.Singularity Function:
Consider a unit step voltage applied to a capacitor as shown in figure 7. The current i
through the capacitor is given by i = C
If the same voltage in ramp signal as shown in figure 8. The current through the
capacitor would be a pulse as shown in figure 9.
Fig: 8 Fig: 9
The solution to an ideal unit step voltage does not exist but it is possible to obtain a
solution in the limit by assuming a non ideal source as shown in figure side by and then
letting a go to zero in the limit. Let us designate the un idealized voltage function by
ua (t). In the limit when a goes to zero, this voltage function becomes a unit step
function. The derivation of the function ua(t) is a rectangular pulse of height and
width a. As „a‟ varies, the pulse shape varies but the area of the pulse remains constant.
In the limit when a goes to zero, the height of the pulse goes to infinity and the width of
the pulse is zero. The area of the pulse however, remains unity. The figure 10 shows
that sequence of such pulses as „a‟ varies.
As „a‟ tends to zero the resultant pulse or function is called as unit impulse function.
Therefore we define the unit impulse function as the derivative of unit step function.
The unit impulse function is denoted by δ(t).Therefore we have
-------(17)
But the function is a rectangular pulse of height and width a. This can be
described as
In the limit as „a‟ goes to zero, the function δ (t) assumes the form of a pulse of infinite
height and zero width. The area of the pulse however remains unity. The function δ (t)
exists only at t =0 and has a zero value of all t ≠ 0. This fact can be expressed as.
δ (t) = 0 t ≠ 0
Equation (19) may be taken as a definition of the unit impulse function. This description
of an impulse function is more general and represents the impulse function as a limiting
process without any regard to the actual shape of the pulse; i.e. it is not necessary to
define the impulse function as a limiting form of Gaussian pulse, a triangular pulse, an
doubled or single sided exponential pulse, or any other pulse forms as shown in fig 11.
Fig: 11
The function δ (t- is a unit impulse function which exists only at t= and is zero
everywhere else.
Hence multiplication of any function x(t) by an impulse δ (t- also yields an impulse
of strength x( ) at t= .
i.e.
This also follows
Here x(t) =
∴F [x(t)] = F [ =
F[ =
Let x(t) = A
In a gate function if z→∞, the gate function tends to a constant function A. The Fourier
Transform of a constant A is therefore the Fourier Transform of a gate function as
z→∞.
We have F
∴ F [A] =
= 2πA
= 2πAδ ( )
∴ F [x(t)] = F [ =
=
=
=π
=π
1 t>0
sgn (t) =
-1 t<0
=
F [sgn (t)] =
= =
F [sgn (t)] =
∴ sgn (t) ↔
u (t) =
∴F [u (t)] =
F [u (t)] =
The spectral density function contains an impulse at = 0. Thus the function u(t)
contains a large d.c component as expected. In addition also has other frequency
component. The function u(t) is not a constant. Since it is zero for t<0, and there is an
abrupt discontinuity at t=0 giving rise to other frequency components.
For a true d.c signal, x(t) is constant for -∞<t<∞ . Such a signal indeed has no other
components except dc (i.e. = 0).
= +j
∴F [ =F[ +j
=π
F[ =2π
The Fourier Transform of is therefore a single impulse of strength 2 π at = .
The signal is not a real function of time, and hence it has a spectrum which exists at
= alone.
2.1.13.The Fourier Transform of a Periodic Function: -
Strictly speaking, the Fourier Transform of a periodic function does not exist, since it
fails to satisfy the condition of absolute integrability. For a periodic function x(t):
x(t) = where
F [x(t)] = F
F [x(t)] = 2π
(∵
The above relation states that the spectral density function or the Fourier Transform of a
periodic signal consists of impulses located at the harmonic frequencies of the signal
and the strength of each impulse is the same as times the value of the corresponding
coefficient in the exponential Fourier series.
=
This is a periodic function with periodic T as shown in figure below.
Fig: 18
Where
From the sampling (shifting) property of impulse function, the above equation becomes.
∴F
F
The above relation is very significant. It states that the Fourier Transform of a
unit impulse train of period T is also a train of impulses of strength and separated by
radians ( Therefore, the impulse train function is its own transform.
The conversion formulae for time domain to frequency domain and vice versa
are known as Fourier Transform pair.
∴ X( ) =
And x(t) =
1. Symmetry Property: -
∴2πx(t) =
2πx(-t) =
Replacing t by .
2πX(- ) =
We get 2πX(- ) =
2. Linearity Property: -
If
3. Scaling Property: -
∴F
The function x(at) compressed in the time scale by factor a. Similarly a function
represents a function X( ) expanded in the frequency scale by the same factor „a‟. The
scaling property therefore states that compression in the time domain is equivalent to
expansion in the frequency domain and vice versa.
If x(t) ↔ X( ) then x ↔ X( )
Proof: F
Let then dt = dx
∴F
This theorem states that if a function is shifted in the time domain by the seconds, then
its magnitude spectrum X( ) remains unchanged out the phase spectrum is changed by
an amount -
Proof: F
=
=X
Similarly
Thus the process of modulation translates the frequency spectrum by amount this
result is also known as “modulation theorem”.
i.e. X( ) =
=> =
The above result are valid only if exists, i.e. if is absolutely integrable. This is
possible only if
i.e.
X( ) = F [x(t)] =
Similarly
x(t) =
x(t) =
If and then
i.e.
Proof: - F
∴F
If then
i.e.
LAPLACE TRANSFORM: (LT)
2.2 INTRODUCTION:
► Laplace Transform can be used to analyze a large class of continuous signals that are
not absolutely integrable, for these signals Fourier Series and Fourier Tranform are not
possible.
►The prime rate of the Laplace Transform in engineering is transient and stability
analysis of casual LTI systems described by differential equations.
S = σ + jω
= eσt (cosωt+jsinωt)
= +j ----------------------> (1)
Real Imaginary
Consider an LTI system
Input x(t)=eSt
= h(t)*x(t)
= eSt
Substitute S = σ + jω
= |H (σ + jω) |eσt cos (ωt +∅ (σ+jω)) + j |H (σ + jω) |eσt sin (ωt +∅ (σ+jω)) --------> (2)
The system changes the amplitude of the input by |H (σ + jω) | and shifts the phase of
the sinusoidal components by ∅ (σ + jω). The system does not change the damping
factor σ of the input or sinusoidal frequency ω.
H (S) =
H (σ + jω) =
= -----------------> (4)
H (jω) =
h (t)e-σt = ejωt dω
h (t) = ejωt..eσtdω
= e(σ+jω)tdω
Substitute S = σ + jω => ds =
ω→∞
S = σ + j∞ h (t) = * =
ω → -∞
S = σ - j∞
h (t) → is inverse Laplace Transform
For arbitrary signal x (t)
=> X (S) = x (t) =
Fig: 20
X(s) =
1 t≥0
Since u(t)= => x(s) =
0 t<0
X(s) = =
X(s) = -
X(s) = - X(s) = =
X(s) = +
To converge x(s) the II term becomes zero. For that (s+a)<0 => => 0
Fig: 22
Similarly
LT
e-2tu(t) ↔ Re(s)>-2
LT
e-3tu(t)↔ Re(s)>-3
Pole-zero plot:
LT
-e-2tu(-t) ↔ Re(s)<-2
LT
e-3tu(-t)↔ Re(s)<-3
X(s) = - = ROC:Re(s)<-3
Fig: 24
EX: 5 Determine the Laplace Transform of the two sided signal x(t) = e -a|t| a>0.
Solution:
eat t<0
x(t) = e-a|t|=
e-at t≥0
eatu(-t) ↔ Re(s)<+a
X(s) = - =
e-a|t| ↔
ROC: -a<Re(s)<a
eat 0<t<T
x(t) =
0 otherwise
X(s) = = =
= X(s) = [1- ]
x(-a) = =T
Solution:
X(s) = =
= = =1
X(s) = = =
X(s) =
s>0 =>term I becomes zero.
X(s) = Re(s)>0 LT
u(t) ↔ ; Re(s)>0
X(s) = = =
X(s) = = -
LT[x(t)] =
e-atu(t) ↔ Re(s)>-a
↔ Re(s)>0
LT
↔ Re(s)>0
X(s) = =
X(s) = =
x(t) =
e-atu(t)↔ ; Re(s)>-a
↔ Re(s)>0
↔ Re(s)>0
X(s) = =
1- Knowledge of ROC can be used to determine a unique ILT. The most common form
of other knowledge is causality, stability or existence of the Fourier Transform.
ILT → is obtained by comparing the pole locations to the jω axis. If pole is right of jω
axis the right sided ILT. If pole is lies on the left of jω axis, then the left – sided Inverse
transform is chosen.
Properties of the Laplace Transform:
The properties of the bilateral Laplace Transform are similar to those of the Fourier
transform. The derivations of many of these properties are analogous to those of the
corresponding properties of the Fourier transform.
x2(t) LT
↔ X2(s) with ROC = R2
LT
ax1(t) + bx2(t) ↔ aX1(s) + bX2(s)
=a
= aX1(s) + bX2(s)
Proof: LT [x (t – t0)] =
dτ = dt t = τ + t0
τ → ∞ as t → -∞ and τ → ∞ as t → ∞
LT [x (t) ]= e-st dt
= dt
= X(s – s0).
The ROC associated with X(s-s0) is that of X(s), shifted by R . Note that if X(s) has a
pole or zero at s=a, then X(s-s0) has a pole or zero at (s-s0) = a i.e. s = (a + s0).
IV Time Scaling
τ → ∞ as t = ∞ =
=
= X
=
L{x (at)} = X
=L-1[S X (s)]
∴ ↔ S X(s)
↔ Sn X (s)
Proof: X (s) =
L [-t x (t)] =
t x (t) ↔
x (t) ↔ (-1)n
The convolution property states that convolution in the time domain corresponds to
multiplication in the S-Domain.
L [x1(t) * x2(t)] = dt
=
Inter changing the order of Integration.
= .
Using the time shift Property the term in the bracket becomes X2(s) e-sτ
=x2(s) dτ = X1(s)X2(s)
x (t) ↔ X (s)
x* (t) ↔ X* (s*)
L [x*(t)] = = = [X(s*)]*
= X*(s*)
IX Multiplication property:
= x2(t)e-stdt
Interchanging the order of Integration and noting that X1(s1) does not depend on t1,then
we have
= ds1
Using the shifting in the S-domain property the bracketed term is X2(s-s1).
LT
x1 (t) x2(t) ↔ [X1(s) * X2(s)]
Proof: If a signal x (t) is convolved with a unit step function u(t). We get
x (t) * u (t) = dτ
=0 t<τ
x (t) * u (t) =
The convolution of a signal with a unit step function is the same as the cumulative integral
of the signal.
= x (t) * u (t)
L =L
u (t) ↔
= X (s) u (s)
= X (s)
The Laplace transform of a casual signal can be determined from the knowledge of the
Laplace transform of its first cycle. Consider periodic signal x (t) with period T o.
Fig:26
Let x1(t), x2(t), x3(t) be the signals representing the 1st, 2nd, 3rd cycles of a periodic signal
x (t),therefore x(t) can be written as,
= X1(s)
=X1(s) = X1(s) )n
X (s) =
Fig: 27
y (s) = X (s) H (s)
H (s) =
The impulse response h (t) is defined as the output of an LTI system due to a unit
impulse signal input x (t) applied at time t=0.
h (t) ↔ H (s)
Where H (s) is the system function or transfer function. The step response s(t) is defined
as the output of an LTI system due to a unit step signal x(t) = u(t)
The step response of an LTI system with Impulse response h(t) is given by
= h (t) *u(t)
=
u(t – τ) = 1 t>τ
s (t) = = 0 t<τ
L [s (t)] = L
= (Using Integration property.)
In many applications we are considering the causal signals. In causal system we are
interested in the behavior of the system t ≥ 0.
In causal system we do not need to consider the ROC. All are right sided sequences.
X (s) =
The ULT and BLT are same for signals that are zero for t <0.
Differentiation in time domain:
x (t) ↔
Integrating by parts
x (t) ↔ x (t) +s
= x (t)
= sX (s) – x ( )
General form:
= ______ x (t)
The initial and final value theorem allows us to determine the initial value of x (t), x (0 +),
and the final value from X(s). These theorems are most often used to evaluate either the
initial or final values of a system output without explicitly determining the entire time
response of the system.
x (o+) =
x (∞) =
Proof
LT = = s X(s) – x (0)
lim s → ∞ = lim s → ∞
lim s → 0 = lim s → 0
= lim s → 0
x(∞) = lim s → 0