Module 17

Inverse Laplace Transform and Waveform Synthesis

Objective:(i) To describe how to obtain inverse Laplace transform making use of the
knowledge of properties of Laplace Transform and properties of ROC.
(ii) To Apply waveform synthesis to find Laplace forms of certain functions
Introduction :
Inverse Laplace transform maps a function in s-domain back to the time domain. One
application is to convert a system response to an input signal from s-domain back to the time
domain.
Since system analysis is usually easier in s-domain, the process is to convert the
system time domain representation to s-domain (both system and inputs),perform system
analysis in s-domain and then convert back to the time domain representation for the
response. The reason to do this process in this convoluted way is that due to its properties, the
Laplace transform converts the differential equations that describe system behaviour to a
polynomial. Also the convolution operation which describes the system action on the input
signals is converted to a multiplication operation. These two properties make it much easier
to do systems analysis in the s-domain.
Inverse Laplace transform is performed using Partial Fraction Expansion that split up
a complicated fraction into forms that are in the Laplace Transform table.
Description :
Concept of Inverse Laplace Transform
We are aware that the Laplace transform of a continuous signal x(t) is given by
∞

𝑋 𝑠 = 𝑥(𝑡)𝑒 −𝑠𝑡 𝑑𝑡
−∞

Since s=σ+jω

∞ ∞
−(σ+jω)𝑡
𝑋 σ + jω = 𝑥(𝑡)𝑒 𝑑𝑡 = 𝑥(𝑡)𝑒 −σ𝑡 𝑒 −jω𝑡 𝑑𝑡 = ℱ 𝑥(𝑡)𝑒 −σ𝑡
−∞ −∞
∞
1
𝑥 𝑡 𝑒 −σ𝑡 = ℱ −1 𝑋 σ + jω = 𝑋 σ + jω 𝑒 jω𝑡 𝑑ω
2𝜋
−∞

We can recover x(t) from its Laplace transform evaluated for a set of values of s=σ+jωin the
ROC, with σfixed and ω varying from -∞ to +∞. Recovering s(t) from X(s) is done by
changing the variable of integration in the above equation from ωto s and using the fact that σ
is constant, so that ds=jdω.
∞ ∞ σ+j∞
1 1 1
𝑥 𝑡 = 𝑋 𝑠 𝑒 σ𝑡 𝑒 jω𝑡 𝑑ω = 𝑋 𝑠 𝑒 s𝑡 𝑑ω = 𝑋 𝑠 𝑒 s𝑡 𝑑𝑠
2𝜋 2𝜋 2𝜋𝑗
−∞ −∞ σ−𝑗 ∞
The contour of integration in above equation is a straight line in the s-plane corresponding to
all points s satisfying Re{s}=σ. This line is parallel to the jω-axis. Therefore, we can choose
any value of σ such that 𝑋 σ + jω converges.
Partial Fraction Expansion
As we know that the rational form of X(s) can be expanded into partial fractions,
Inverse Laplace transform can be taken according to location of poles and ROC of X(s). The
roots of denominator polynomial, i.e., poles can be simple and real, complex or multiple.
We know that X(s) is expanded in partial fractions as
𝑐0 𝑐1 𝑐2 𝑐𝑛
𝑋 𝑠 = + + + ⋯+
𝑠 − 𝑠0 𝑠 − 𝑠1 𝑠 − 𝑠2 𝑠 − 𝑠𝑛
Here the roots s0,s1,s2,...sn can be real, complex or multiple. Then the values of
k0,k1,k2,...kn constants are calculated accordingly.
In order to find the appropriate time domain function, ROC should be indicated for
the s-domain function. Otherwise we may have multiple time-domain functions based on
different possible ROCs.
Illustration
Example for Real roots:
Problem 1:Find out the partial fraction expansion and hence Inverse Laplace transform of the
𝑠 2 +2𝑠−2
function 𝑋 𝑠 = 𝑠(𝑠+2)(𝑠−3) , ROC: Re{s} > 3

Solution:
𝑠 2 +2𝑠−2
The function 𝑋 𝑠 = 𝑠(𝑠+2)(𝑠−3) can be written as,

𝐴 𝐵 𝐶
𝑋 𝑠 = + +
𝑠 𝑠+2 𝑠−3
The constants calculated are A=1/3, B=-1/5, C=13/15
1/3 1/5 13/15
𝑋 𝑠 = + +
𝑠 𝑠+2 𝑠−3
1 1 1 1 13 1
𝑥 𝑡 = ℒ −1 + ℒ −1 + ℒ −1
3 𝑠 5 𝑠+2 15 𝑠−3
From the given ROC: Re{s}>3, the resultant signal x(t) should be right sided.
1 1 13
Therefore, 𝑥 𝑡 = 3 𝑢 𝑡 + 5 𝑒 −2𝑡 𝑢(𝑡) + 15 𝑒 3𝑡 𝑢(𝑡)

Example for Complex roots:

𝑠 2 +2𝑠+1
Problem 2: Obtain right sided time domain signal for the function 𝑋 𝑠 = (𝑠+2)(𝑠 2 +4)
Solution:
𝑠 2 +2𝑠+1 𝐴 𝐵𝑠+𝐶
We can write the given function as 𝑋 𝑠 = (𝑠+2)(𝑠 2 +4) = 𝑠+2 + 𝑠 2 +22

The constants can be calculated as A=1/8, B=0.874, C=0.5

1 1 0.874𝑠+0.5 1 1 𝑠 2
Therefore, 𝑋 𝑠 = 8 𝑠+2 + =8 + 0.874 + 0.25
𝑠 2 +22 𝑠+2 𝑠 2 +22 𝑠 2 +22

1 1 𝑠 2
Finally 𝑥(𝑡) = 8 ℒ −1 + 0.874ℒ −1 + 0.25ℒ −1
𝑠+2 𝑠 2 +22 𝑠 2 +22
1
i.e., 𝑥(𝑡) = 8 𝑒 −2𝑡 𝑢(𝑡) + 0.874 cos 2𝑡 𝑢(𝑡) + 0.25𝑠𝑖𝑛 2𝑡 𝑢(𝑡)

Example for Multiple roots:

𝑠−2
Problem 3: Find out the inverse Laplace transform of 𝑋 𝑠 = 𝑠(𝑠+1)3 , ROC:Re{s}<-1

Solution:
𝑠−2 𝐴 𝐵 𝐶 𝐷
We can write the given function as 𝑋 𝑠 = 𝑠(𝑠+1)3 = + + (𝑠+1) +
𝑠+1 3 𝑠+1 2 𝑠

The constants can be calculated as A=3, B=2, C=2, D=-2

1 1 1 1
Therefore, 𝑋 𝑠 = 3 3
+2 2
+2 −2
𝑠+1 𝑠+1 (𝑠+1) 𝑠

From the given ROC: Re{s}<-1, the resultant signal x(t) should be left sided.
1 1 1 1
Finally,𝑥 𝑡 = 3ℒ −1 + 2ℒ −1 + 2ℒ −1 − 2ℒ −1
𝑠+1 3 𝑠+1 2 (𝑠+1) 𝑠

From the result of Laplace transform

−𝑡 𝑛 −𝑎𝑡 ℒ 1
𝑒 𝑢 −𝑡 𝑛+1
, 𝑅𝑂𝐶: 𝑅𝑒 𝑠 < −𝑎
𝑛! 𝑠+𝑎

−𝑡 2
i.e.,𝑥 𝑡 = 3 2
𝑒 −𝑡 𝑢 −𝑡 − 2𝑡𝑒 −𝑡 𝑢 −𝑡 − 2𝑒 −𝑡 𝑢 −𝑡 + 2𝑢 −𝑡

Laplace transform using Waveform Synthesis

In waveform synthesis, the unit step function u(t) and other functions serve as building blocks
in constructing other waveforms. Once the waveforms are synthesized in the form of other
functions, Laplace transform is found and simplified.
Illustration
For example, we may describe a pulse waveform in terms of unit step functions. A pulse of
unit amplitude from t=a to t=b can be formed by taking the difference between the two step
functions
i.e.,x(t) = u(t-a) - u(t-b)
1 1
Hence 𝑋 𝑠 = 𝑒 −𝑠𝑎 𝑠 − 𝑒 −𝑠𝑏 𝑠

Examples:
Solved Problems:
3𝑠+7
Problem 1: Find the Inverse Laplace transform of 𝑋(𝑠) = 𝑠 2 −2𝑠−3 ROC: Re{s}>3

Solution:
ℒ 1
We know that 𝑒 −𝑎𝑡 𝑢(𝑡) , 𝑅𝑂𝐶: 𝑅𝑒 𝑠 > −𝑎
𝑠+𝑎

Writing X(s) in the form of partial fraction expansion

𝐴 𝐵
𝑋 𝑠 = +
𝑠−3 𝑠+1
The constants can be calculated as A=4, B=-1
𝐴 𝐵 1 1
Therefore, 𝑋 𝑠 = 𝑠−3 + 𝑠+1 = 4 −
𝑠−3 𝑠+1

From the given ROC: Re{s}>3, the resultant signal x(t) should be right sided.
1 1
i.e., 𝑥 𝑡 = 4ℒ −1 − ℒ −1 = 4𝑒 3𝑡 𝑢 𝑡 − 𝑒 −𝑡 𝑢 𝑡
𝑠−3 𝑠+1

−3
Problem 2: Find the Inverse Laplace transform of 𝑋(𝑠) = (𝑠+2)(𝑠−1) ROC: Re{s}<-2

Solution:
ℒ 1
We know that −𝑒 −𝑎𝑡 𝑢(−𝑡) , 𝑅𝑂𝐶: 𝑅𝑒 𝑠 < −𝑎
𝑠+𝑎

Writing X(s) in the form of partial fraction expansion

𝐴 𝐵
𝑋 𝑠 = +
𝑠+2 𝑠−1
The constants can be calculated as A=1, B=-1
𝐴 𝐵 1 1
Therefore, 𝑋 𝑠 = 𝑠+2 + 𝑠−1 = −
𝑠+2 𝑠−1
From the given ROC: Re{s}<-2, the resultant signal x(t) should be left sided.
1 1
i.e., 𝑥 𝑡 = ℒ −1 − ℒ −1 = −𝑒 −2𝑡 𝑢 −𝑡 + 𝑒 𝑡 𝑢 −𝑡
𝑠+2 𝑠−1
1
Problem 3: Find the Inverse Laplace transform of 𝑋(𝑠) = 𝑠 2 +3𝑠+2 ROC: -2<Re{s}<-1

Solution:
We know that
ℒ 1
𝑒 −𝑎𝑡 𝑢(𝑡) , 𝑅𝑂𝐶: 𝑅𝑒 𝑠 > −𝑎
𝑠+𝑎
ℒ 1
−𝑒 −𝑎𝑡 𝑢(−𝑡) , 𝑅𝑂𝐶: 𝑅𝑒 𝑠 < −𝑎
𝑠+𝑎
Writing X(s) in the form of partial fraction expansion
𝐴 𝐵
𝑋 𝑠 = +
𝑠+1 𝑠+2
The constants can be calculated as A=1, B=-1
𝐴 𝐵 1 1
Therefore, 𝑋 𝑠 = 𝑠+1 + 𝑠+2 = −
𝑠+1 𝑠+2

From the given ROC: -2<Re{s}<-1, the two derived conditions are Re{s}<-1 which suits for
1 1
and Re{s}>-2 which suits for . Therefore, the resultant signal x(t) should be two-
𝑠+1 𝑠+2
sided.
1 1
i.e., 𝑥 𝑡 = ℒ −1 − ℒ −1 = −𝑒 −𝑡 𝑢 −𝑡 − 𝑒 −2𝑡 𝑢 𝑡
𝑠+1 𝑠+2
1
Problem 4:Determine right-sided x(t) if 𝑋(𝑠) = 𝑠 2 (𝑠 2 −𝑎 2 ) using convolution theorem

Solution:
1
The given function is 𝑋 𝑠 = = 𝑋1 𝑠 . 𝑋2 (𝑠)
𝑠 2 𝑠 2 −𝑎 2

1 1 1 1 1
𝑋1 𝑠 = 𝑠 2 and 𝑋2 𝑠 = = 2𝑎 − 𝑠+𝑎 (Done using partial fraction expansion)
𝑠 2 −𝑎 2 𝑠−𝑎

1
𝑥1 𝑡 = ℒ −1 = 𝑡𝑢(𝑡)
𝑠2
1 −1 1 1 1 𝑎𝑡
𝑥2 𝑡 = ℒ − ℒ −1 = 𝑒 − 𝑒 −𝑎𝑡 𝑢(𝑡)
2𝑎 𝑠−𝑎 𝑠+𝑎 2𝑎
ℒ
From the Convolution property𝑥1 𝑡 ∗ 𝑥2 𝑡 𝑋1 𝑠 . 𝑋2 (𝑠)
∞ 1 𝑡
i.e., 𝑥 𝑡 = 𝑥1 𝑡 ∗ 𝑥2 𝑡 = 𝑥
−∞ 1
𝜏 𝑥2 𝑡 − 𝜏 𝑑𝜏 = 2𝑎 0
𝜏 𝑒 𝑎(𝑡−𝜏) − 𝑒 −𝑎(𝑡−𝜏) 𝑑𝜏

1 2𝑡 1
𝑥 𝑡 = − + 2 𝑒 𝑎𝑡 − 𝑒 −𝑎𝑡 ;𝑡 > 0
2𝑎 𝑎 𝑎
 1+𝑒 −2𝑠
Problem 5:Find right sided x(t)if𝑋 𝑠 = 3𝑠 2 +2𝑠

Solution:
1+𝑒 −2𝑠 1 𝑒 −2𝑠 1 1 1 1 1 1
Given𝑋 𝑠 = 3𝑠 2 +2𝑠 = 𝑠(3𝑠+2) + 𝑠(3𝑠+2) = 2 − 𝑠+2/3 + 2 𝑒 −2𝑠 − 𝑠+2/3
𝑠 𝑠

Taking inverse Laplace transform and using time shifting property

1 2 1 2
𝑥 𝑡 = 1 − 𝑒 −3𝑡 𝑢(𝑡) + 1 − 𝑒 −3 𝑡−2
𝑢(𝑡 − 2)
2 2
Problem 6:Given, 𝑥(𝑡) = 𝑒 −𝑡 𝑢(𝑡), find the Inverse Laplace transform of 𝑒 −3𝑠 𝑋(2𝑠)
Solution:
From Time scaling property
ℒ
If 𝑥(𝑡) 𝑋(𝑠)
ℒ 1 𝑠
then 𝑥 𝑎𝑡 𝑋
𝑎 𝑎
𝑡 ℒ 1 𝑠
with a=1/2, 𝑥 1 𝑋 1 = 2𝑋(2𝑠)
2 2 2
1 𝑡 ℒ
Therefore, 2 𝑥 𝑋(2𝑠)
2
ℒ
Time shifting property states that then 𝑥 𝑡 − 𝜏 𝑒 −𝑠𝜏 𝑋 𝑠 . Applying this property
to above equation with τ=3,
1 𝑡 − 3 ℒ −3𝑠
𝑥 𝑒 𝑋(2𝑠)
2 2
Applying the transformation to 𝑥(𝑡) = 𝑒 −𝑡 𝑢(𝑡) given by above equation,
1 𝑡−3 1 𝑡−3 𝑡−3
ℒ −1 𝑒 −3𝑠 𝑋(2𝑠) = 𝑥 = 𝑒− 2 𝑢
2 2 2 2

Problem 7: Find the steady state response of the following system to unit step excitation
𝑠+1
𝐻 𝑠 =
𝑠 2 + 3𝑠 + 2
Solution:
We know that for a linear system outputy(t) =h(t)*x(t)
ℒ 1
As the input is𝑥(𝑡) = 𝑢(𝑡) 𝑠
ℒ 𝐻(𝑠)
Output 𝑦 𝑡 = ℎ 𝑡 ∗ 𝑢 𝑡 𝑌 𝑠 = 𝑠
𝐻(𝑠) 𝑠+1 𝑠+1 1
Therefore, 𝑌 𝑠 = = 𝑠(𝑠 2 +3𝑠+2) = 𝑠(𝑠+1)(𝑠+2) = 𝑠(𝑠+2)
𝑠

Using partial fraction expansion

 1 1 1
𝑌 𝑠 = −
2 𝑠 𝑠+2
1−𝑒 −2𝑡
Therefore, step response 𝑦 𝑡 = 𝑢(𝑡)
2

𝑑
Problem 8: Knowing that 𝑈ℒ𝑇 𝑦(𝑡) = 𝑠𝑌 𝑠 − 𝑦(0−)
𝑑𝑡
𝑑
Solve the differential equation 𝑦 𝑡 + 5𝑦 𝑡 = 𝑥(𝑡), with initial condition y(0+)=-2 and
𝑑𝑡
input x(t)=3e-2tu(t)
Solution:
𝑑
The given differential equation is, 𝑑𝑡 𝑦 𝑡 + 5𝑦 𝑡 = 𝑥(𝑡)

Taking Unilateral Laplace transform of above equation

sY(s)-y(0-)+5Y(s)=X(s)
ℒ 1
Substituting initial condition y(0+)=y(0-)=-2 and input 𝑥 𝑡 = 3𝑒 −2𝑡 𝑢 𝑡 𝑋 𝑠 = 3 𝑠+2
and applying partial fraction expansion
3 2 1 1 2 1 3
𝑌 𝑠 = − = − − = −
(𝑠 + 2)(𝑠 + 5) 𝑠 + 5 𝑠 + 2 𝑠 + 5 𝑠 + 5 𝑠 + 2 𝑠 + 5
Taking inverse Laplace transform of the above equation
𝑦 𝑡 = 𝑒 −2𝑡 𝑢 𝑡 − 3𝑒 −5𝑡 𝑢 𝑡
Problem 9: Obtain the impulse response of the system shown and hence prove that the
system is BIBO stable

Solution:
Let us transform all the elements to their Laplace equivalents assuming zero initial
conditions. The Laplace equivalent circuit is shown below
Using voltage divider formula we can write,
1
𝑠𝐶 1
𝑌 𝑠 = 1 𝑋 𝑠 = 𝑋(𝑠)
𝑅 + 𝑠𝐶 𝑠𝑅𝐶 + 1

Therefore, transfer function

𝑌(𝑠) 1 1 1
𝐻 𝑠 = = =
𝑋(𝑠) 𝑠𝑅𝐶 + 1 𝑅𝐶 𝑠 + 1
𝑅𝐶

Taking inverse Laplace transform of above equation, impulse response is

1 −𝑡
ℎ 𝑡 = 𝑒 𝑅𝐶 𝑢(𝑡)
𝑅𝐶
H(s) has one pole located at s=-1/RC
This pole lies in left half of the s-plane. Hence the system is both causal and stable.
This can also be proved checking the absolute integrability of impulse response
∞
ℎ(𝑡) 𝑑𝑡 < ∞
−∞

𝑡 ∞
∞
1 ∞
−
𝑡 1 𝑒 −𝑅𝐶
ℎ(𝑡) 𝑑𝑡 = 𝑒 𝑅𝐶 𝑑𝑡 = = − 𝑒 −∞ − 𝑒 0 = 1 < ∞
−∞ 𝑅𝐶 0 𝑅𝐶 − 1
𝑅𝐶 0

Hence BIBO stable

Problem 10:Obtain the Laplace transform of the triangular pulse as shown below using
waveform synthesis

Solution:
As shown in the figure above, a ramp of slope „A‟ starting at t=0 is taken as f1(t). The
function f2(t) is the ramp of slope -2A starting at t=1. When we add f1(t) and f2(t) we get the
signal shown in figure (c). Observe that a negative going ramp of slope (A-2A=-A) starts at
t=1. To cancel the part of this ramp after t≥2, a positive going ramp of slope +A in figure (d)
is added. Then we get the required triangular pulse of Figure (e).
With the help of step functions, the ramp functions in the figure above can be expressed as
follows:
f1(t)=A t u(t) { The function u(t)=1 for t>0; It indicates that ramp is present only for t≥0}
f2(t)=-2A (t-1) u(t-1) { The function u(t-1)=1 for t>1; It indicates that ramp is present only for
t≥1}
f3(t)=A(t-2) u(t-2)
Therefore, f(t)=f1(t)+f2(t)+f3(t)=A t u(t)-2A (t-1) u(t-1)+A(t-2) u(t-2)
Since u(t), u(t-1) and u(t-2) have values of „1‟ and they just represent the time shifts and
directions of ramp functions, they can be dropped in this expression. Laplace transform of
above equation becomes
1 𝑒 −𝑠 𝑒 −2𝑠 𝐴
ℒ 𝑓 𝑡 =𝐴 2
− 2𝐴 2
+ 𝐴 2
= 2 1 − 2𝑒 −𝑠 + 𝑒 −2𝑠
𝑠 𝑠 𝑠 𝑠
Assignment:
8
Problem 1:Find the inverse Laplace transform of 𝑋 𝑠 = 4𝑠 2 + 15𝑠 + .Assuming
𝑠+2 2 𝑠+1
signal is causal.
2 1
Problem 2:The transfer function of the system is given as 𝐻 𝑠 = 𝑠+3 + 𝑠−2. Determine the
impulse response if the system is (i) stable (ii) causal
Problem 3:Find the current i(t) in a series RLC circuit as shown in figure below, when a
voltage of 100Volts is switched on across the terminals aa‟ at t=0.
Problem 4: The response h(t) of a linear time invariant system to an impulse δ(t), under
initially relaxed condition is h(t) = e−t+ e−2t. Find the step response of this system.
𝑠
Problem 5: Find the Inverse Laplace transform of 𝐺 𝑠 = (𝑠−3)(𝑠 2 −4𝑠+5) , 𝜎 < 2

Problem 6:Given f (t) and g(t) as show below:

Express g(t) in terms of f(t) and hence find Laplace transform

10𝑠 2
Problem 7: Find the Inverse Laplace transform of 𝐺 𝑠 = (𝑠+1)(𝑠+3), σ > 0

Problem 8:Let the Laplace transform of a function f (t) which exists for t > 0 be F1(s)and the
Laplace transform of its delayed version f (t − τ) be F2(s). Let 𝐹1 ∗ (𝑠)be the complex
𝐹2 𝑠 𝐹1 ∗ (𝑠)
conjugate of F1(s) with the Laplace variable set s = σ + jω.If 𝐺 𝑠 = then Find the
𝐹1 (𝑠) 2
inverse Laplace transform of G(s).
Problem 9:Find the Laplace transform of periodic sawtooth waveform using waveform
synthesis

Problem 10:Consider the LTI system for which we are given the following information:
𝑠+2
𝑋 𝑠 =
𝑠−2
𝑥 𝑡 =0; 𝑡 > 0
2 1
and 𝑦 𝑡 = − 3 𝑒 2𝑡 𝑢 −𝑡 + 3 𝑒 −𝑡 𝑢 𝑡

(a) Determine H(s) and its ROC

(b) Determine h(t)

Simulation:
The command one uses now is ilaplace. One also needs to define the symbols t and s.
Let‟s calculate the inverse of the previous function F(s),
𝑠−5
𝐹 𝑠 =
𝑠(𝑠 + 2)2

>>syms t s
>> F=(s-5)/(s*(s+2)^2);
>>ilaplace(F)
ans =
-5/4+(7/2*t+5/4)*exp(-2*t)
>> simplify(ans)
ans =
-5/4+7/2*t*exp(-2*t)+5/4*exp(-2*t)
>> pretty(ans)
- 5/4 + 7/2 t exp(-2 t) + 5/4 exp(-2 t)
Which corresponds to
𝑓 𝑡 = −1.25 + 3.5𝑡𝑒 −2𝑡 + 1.25𝑒 −2𝑡

Alternatively, one can write

>>ilaplace((s-5)/(s*(s+2)^2))

Here is another example.

10(𝑠 + 2)
𝐹 𝑠 =
𝑠(𝑠 2 + 4𝑠 + 5)

>> F=10*(s+2)/(s*(s^2+4*s+5));
>>ilaplace(F)
ans =
-4*exp(-2*t)*cos(t)+2*exp(-2*t)*sin(t)+4
Another example is shown below
>>syms s ;
>>T = 1;
>>F1 = 1 - exp( -T * s )+ exp(-3*T * s/2);
>>F2 = s^2 ;
>>F = F1 / F2;
>>f = ilaplace( F )

f =

t - heaviside(t - 1)*(t - 1) + heaviside(t - 3/2)*(t - 3/2)

>>pretty( simplify ( f ) )

t - heaviside(t - 1) (t - 1) + heaviside(t - 3/2) (t - 3/2)

>>ezplot( f )
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