3 LTI System
3 LTI System
x (t )
p(t)
1/
t t
–/2 /2
Fig. 2.1 : Arbitrary signal x(t) and practical delta function p(t)
Now, x(t) p(t) is x(0) p(t). Plot of x(0) p(t) is shown in Fig. 2.2(a) and plot of x(0) p(t) is shown in
Fig. 2.2(b). Similarly Fig. 2.2(c) and Fig.2.2 (d) show x(t) p(t – ) i.e. x() p(t – ) and x() p(t – ).
x(0)p( t) x(0)p(t )
x(0)/ x(0)
t t
–/2 /2
(a) (b)
x()p( t – ) x()p( t – )
x( )/ x()
t t
(c) (d)
x ( t ) = [ x(k )
k
p (t k )] ...(1)
t
5 3 0 3 5 7
2 2 2 2 2 2 2
Since, 0 will make input practical delta function into ideal delta function and summation of equation (1)
will change to integration. Thus,
x(t) =
x ( ) ( t ) d ...(3)
The above equation show how any time domain signal can be written in terms of impulse function.
Remember: We can intuitively think that x(t ) as a sum of weighted shifted impulses, where the weight on the impulse
(t ) is x()d
LTI
(t ) h( t)
system
Impulse Impulse
input response
Sy stem
( t ) h( t )
Since, the system is time invariant so.
Sy stem
( t ) h( t )
Since system is linear so if we multiply input with a constant then output also get multiplied by same
constant.
Sy stem
x( ) ( t ) x( ) h( t )
Since, integration is a linear operation applying this operation at input and output side.
Sy stem
x( ) ( t ) d
x ( ) h( t ) d
Sy stem
x ( t ) x ( ) h( t ) d
So output, y(t) = x ( ) h( t ) d
...(4)
Equation (4) is referred as the convolution integral or superposition integral and symbolically
represented as,
y ( t ) = x(t) h(t) ...(5)
Thus, in case of any linear time invariant system we can find output for any arbitrary input x(t) by
convolution of x(t) with impulse response of system h(t).
Procedure that should be followed to find out convolution of two signals:
Suppose, y ( t ) = x ( t ) h( t ) x( ) h( t ) d
Signals x(t) and h(t) are given and we have to find y(t), or we have to compute the above integral. Here
the integration variable is ‘’ and ‘t’ is a constant.
S t e p - 1 : Change variable ‘t’ by ‘’ in x(t) to get x().
S t e p - 2 : From h(t) find h() by changing variable ‘t’ by ‘’ then plot h(–) by taking mirror image of h()
from amplitude axis. To get h(t – ) we have to do the required shifting in h(–).
S t e p - 3 : To find y(t) multiply x() and h(t – ) and find, x ( ) h( t ) d .
x(t) h(t)
1.0
e –t
1.0
t t
0 0
x() h( )
1.0
e–
1.0
0 0
h( –)
1.0
0
Now to find h(t – ) we have to shift h(–). For t > 0 the h(–) will shift right by t units.
h( t – )
1.0
(for t > 0)
t
t 0
x( ) h(t – )
1
–
–
1.0 e
e
× =
t 0 t
t
For t > 0, y(t) = e d
0
= (1 – e–t)
Case-2: For t < 0, plot of x() and h(t – ) are
x( ) h( t – )
1.0
–
e
× = 0
t
(1 e t ), t 0
so, y(t) =
0, t 0
= (1 – e–t) u(t)
Example 2.1
–t –t –t
e e e
t t
1.0
– –
e e e
1.0
e
–1.0
Plot of h(t – ) will be shifted h(–) by ‘t’, if t is positive h(–) shift right by ‘t’ units and if t is negative h(–) shift
left by ‘t’ units. Another way is to use equation to visualise, h(t) = et u(t 1)
These two figures show two cases when t > 1 and t < 1.
h(t – ) h(t – )
e–(t – ) e–(t – )
( t – 1) ( t – 1)
y(t) = x( ) h( t ),
putting x() we get
0
= e h(t ) d e h(t ) d
0
Analysing the plot of h(t – ) we can see that if t < 1 then h(t – ) = 0 for t > 0. So for t < 1
0
y(t) = e
h( t ) d 0
e( t 2)
y(t) = for t 1
2
Now for t > 1, h(t – ) = e–(t – ) for < t – 1
0 t 1
y(t) = e e( t ) dt e
e ( t ) dt
0
0 t 1
t 2 t e t
= e e dt e d = et ( t 1)
0 2
e t
y(t) = e t ( t 1) for t 1
2
Example 2.2
Consider a LTI system with impulse response h(t) and input x(t), the output of system is y(t). The value of
output at t = 5 is
x (t ) h( t)
1.0
2.0
t t
–5 –4 4 5 0 10
Solution 2.2
Since y(t) is output,
y ( t ) = x(t) h(t)
y(t) = x ( ) h ( t ) d
1.0 2.0
–5 –4 4 5 –5 5
So, y(5) = x( ) h(5 ) d
2.0
–5 –4 4 5
So, y(5) = 18
Example 2.3
x( ) x( –)
–T/2 T/2 –T/2 T/2
The plot of x(t – ) will be, for t > 0 the signal x(–) will move right side by t units
So, when t > 0
x(t – )
y(t) = x ( ) x( t ) d
T T
A,
we can seen that, x() = 2 2
t – /2 t + /2
0 else wh ere
T T
A, t t
and x(t – ) = 2 2
0 else wh ere
T /2
y(t) = A2 d
t T / 2
y ( t ) = A2 (T – t)
So, y(t) = A2 (T – t) when t > 0 and t < T
because overlap between x() and x(t – ) exist for t < T and for t > T the overlap of x() and x(t – ) will
be zero and their product will also be zero thus y(t) = 0 for t > T.
Now when t is negative then x(–) will move to left by t units, so plot of x(t – ) for t < 0 is
x(t – )
t – /2 t + /2
T T
A,
Now, x() = 2 2
0 , else wh ere
T T
A, t t
and x(t – ) = 2 2
0 , else wh ere
t T / 2
y(t) = x ( ) x ( t ) d = A2
T / 2
y(t) = A2(t
+ T) for t < 0 and t > – T
because for t < –T overlap between x() and x(t – ) do not exist and y(t) = 0 for t < –T.
2
AT
–T T
Study Notes
Generally in GATE exam whenever we have to find convolution of two signals we will not use the method shown above
because this is a lenghty method, it is better to use laplace transform to do the same. We will learn laplace transform in
Chapter 4. The above method will be helpful in those questions where x(t) is input to an LTI system with impulse response
h(t) and we have to find value of output at a particular time instant as shown in Example 2.2.
x (t ) h( t) y ( t) = h( t ) x (t ) y ( t)
Fig. 2.5
Proof: By definition,
x(t) h(t) = x ( ) h( t ) d
If we perform a change of variables by letting (t – ) = , then = (t – ), d = –d, as –, and
– as . Therefore,
Study Notes
Suppose, y(t) = x(t) h(t),
to find the result we can perform,
If we have signal x(t) simpler than h(t) then we perform h(t) x(t) and if h(t) is simpler then we perform
x(t) h(t).
x(t) h(t) = x ( ) h( t ) d
x ( t) h1 ( t ) h2 ( t ) y (t ) x (t ) h1 ( t ) h2 ( t ) y (t )
Fig. 2.6
h1(t)
h2(t)
Fig. 2.7
Proof: By definition,
x(t) h(t) = x( ) h( t ) d y( t )
x(t) h (t – t0) = x( ) h( t t 0 ) d y( t t 0 )
b
x t
1 b a
x(t) (at + b) = x ( t ) ( t )
a a a
Study Notes
To perform the convolution of a signal with delta function, we should follow these steps
1
• First of all find the area and location of delta function and write delta function as ( t -location) .
area
• Now convolution of signal x(t) with given delta fuction is found by replacing variable ‘t’ in x (t) by (t - location) and
divide expression by area of delta function.
Study Notes
Suppose y(t) = x(t) h(t), signal x(t) is non-zero for TxL < t < TxH and h(t) is non-zero for ThL < t < ThH then y(t) will be non-
zero only for ThL + TxL < t < TxH + ThH ,and width of x(t) and h(t) is Wx = TxHTxL and Wh = ThHThL.
d d d
then, x ( t ) h( t ) = x ( t ) h( t ) y( t ) ...(12)
dt dt dt
Proof : By definition,
y(t) = x ( t ) h( t ) x( ) h( t ) d
d2 d 2 x( t ) d 2 h( t ) dx ( t ) dh( t )
2
y( t ) = h( t ) x ( t )
2
dt dt dt 2 dt dt
x(at) h(at) = x ( a ) h( a( t )) d x( a) h( at a) d
1
d = d , as , and as
a
1 1
Therefore, x(at) h(at) = x ( ) h( at d y( at )
a a
Case-II: a < 0
By definition, x(t) h(t) = x ( ) h( t d y( t )
1
d = d , as – and – as
a
1 1
Therefore, x(at) h(at) = x( ) h( at ) d a y( t )
a
From the above two cases it is evident that,
1
x(at) h(at) = y( at )
a
• x( t ) u ( t ) x( ) d
Example 2.4
y ( t)
3
h( t)
x (t ) 2 2
1
1 1
t t t
–1 0 1 –2 –1 0 –3 –2 –1 0
Example 2.5
Consider an LTI system with input and output related through the following equations:
t
(t )
y(t) = e x( 2) d
A change of variables is performed by letting = ( – 2), which also yields d = d, – as –,
and (t – 2) as t. Therefore,
t2
y(t) = e( t 2 ) x ( ) d e
( t 2 )
u ( t 2 ) x( ) d
Comparing with, y(t) = h( t ) x ( ) d
Example 2.6
Study Notes
The step response s(t) of a continuous-time LTI system is defined as the response of the system when the input is u(t).
In many applications, the step response s(t) is also a useful characterization of the system. The step response s(t) can be
easily determined by equation, that is
t
Thus, the step response s(t) can be obtained by integrating the impulse response h(t). Differentiating with respect to t, we
get
d s( t )
h(t) = s( t )
dt
Thus, the impulse response h(t) can be determined by differentiating the step response s(t).
x ( t) h1 ( t ) h2 ( t ) y (t ) x (t ) h( t ) y (t )
h1( t)
+
x ( t) y (t ) x (t ) h1(t ) ± h2( t) y (t )
–
h2( t)
So, parallel connection can be replaced by single system with impulse response,
h(t) = h1(t) ± h2(t)
x[0], n 0
x[0] [n] =
0, n 0
x[1], n 1
x[1] [n – 1] =
0, n 1
More generally, by including additional shifted, scaled impulses, we can write
x[n] = ... + x[–3] [n + 3] + x[–2] [n + 2] + x[–1] [n + 1] + x[0] [n]
+ x[1] [n – 1] + x[2] [n – 2] + x[3] [n – 3] + .... ...(17)
For any value of n, only one of the terms on the right-hand side of equation (17) is non-zero, and the scaling
associated with that term is precisely x[n]. Writing this summation in a more compact form, we have,
x[n] = x[k ] [n k ]
k
...(18)
This corresponds to the representation of an arbitrary sequence as a linear combination of shifted unit
impulses [n – k], where the weights in this linear combination are x[k].
LTI
[n ] Impulse response
system
h[n ]
Sy stem
x[n ]
m
x[m ] h[n m ]
Equation (19) is referred as the convolution sum or superposition sum and symbolically represented as
y[n] = x[n] h[n] =
m
x[m ] h[n m ] ...(20)
Thus in case of any linear time invariant we can find output for any arbitrary input x[n] by convolution of
x[n] with impulse response of system h[n].
Procedure that should be followed to find out convolution of two signals:
Signal x[n] and h[n] are given and we have to find y[n] or we have to compute above summation. In this
summation variable is ‘m’ and ‘n’ is a constant. So we will follow this procedure.
S t e p - 1 : Change variable ‘n’ to ‘m’ in x[n] and h[n] to get x[m] and h[m].
S t e p - 2 : Now we get h[–m] by inversion of h[m]. To get h[n – m] we have to do shifting in h[–m].
S t e p - 3 : To find y[n] we multiply x[m] and h[n – m] and find
m
x[m ] h[n m ] .
1.0 1.0
0.5 .....
.....
n n
0 1 2
(0.5) m , m 0
So, x[m] = (0.5) m u [m ]
0, m 0
1, m 0
and for h[m] = u [ m ]
0 , m 0
Now, h[–m] = u[–m]
1, n m 0 1, m n
Thus, h[n – m] = u[n – m] = =
0 , n m 0 0, m n
For n 0,
h[n – m]
1.0
..... ...
m
–4 –3 –2 –1 0 1 2 n
Study Note
The discrete-time LTI system which have impulse response h[n] non-zero only for finite value of ‘n’ is called finite impulse
response filter on FIR filter. If impulse h[n] is non-zero for infinite values of ‘n’ then it is called infinite impulse response
filter IIR filter.
Example 2.7
0, n < -5 0, n3
n n
x[n] = 1 and h[n] = 1
≥ , n 3
, n - 5
2
3
Solution 2.7
So we have to find, y[n] = x[n] h[n]
=
m
x[m ] h[n m ]
0, m 5
m
we can define, x[m] = 1
, m 5
2
0, nm 3 0, m n 3
nm nm
and h[n – m] = 1 1
, nm 3 , m n3
3 3
Now, using above definition we can say that, product x[m] h[n – m] will be non zero for –5 m n – 3,
n 3 m nm
1 1
y[n] = x[m ] h[n m ]
2 3
m m 5
3 n 3 3 n 3
n 5
1 1
5 n
1 3 2 2 1 2
= =
3 2 3 3 3 1
1 2
2
5 n n 3
2 1 3
So, for n –2, y[n] = 2 1
3 3 2
and for n –2, y[n] = 0
Example 2.8
m 1, N m N
Since, x[m] = rect
2 N 0 , else where
1, N n m N
and x[n – m] =
0 , else wh ere
1, N n m N n
=
0 , else wh ere
when n 0 that is when n is positive then,
N
y[n] = x[m ] x[n m ] 1 (2 N n )
m N n
and for n > 2N, overlap between x[m] and x[n – m] is 0 so,
y[n] = 0 for n > 2N
Now for n < 0, using above definition we get,
N n
y[n] = x[m] x[n m] 1 2N n
m N
and y[n] = 0 for n < –2N as overlap between x[m] and x[n – m] is zero.
So, y[n] = 2N – m 0 n 2N
y[n] = 2N + n –2N n < 0
Example 2.9
Compute the convolution y(n) = x(n) h(n) of the following pairs of signals:
( a ) x(n) = h(n) = u(n)
( b ) x(n) = (0.8) n u(n) and h(n) = (0.4) n u(n)
( c ) x(n) = h(n) = a n u(n)
( d ) x(n) = u(n – 1) and h(n) = n u(n – 1)
( e ) x(n) = r(n) = nu(n) and h(n) = a –n u(n – 1), where a < 1
Solution 2.9
By definition:
(a) y ( n ) = x(n) h(n)
= x ( k ) h( n k ) u( k ) u ( n k )
k k
The lower limit on the convolution sum simplifies to k = 0 [because u(k) = 0, k < 0], the upper limit
to k = n [because u(n – k) = 0, k > n], and we get
n
y(n) = 1 ( n 1) u ( n ) r( n 1)
k 0
The lower limit on the convolution sum simplifies to k = 0 [because u(k) = 0, k < 0], the upper limit
to k = n [because u(n – k) = 0, k < n], and we get
n n
y(n) = (0.8)k (0.4) n k (0.4) n ( 2)k
k 0 k 0
n
n
1 2 n 1
= (0.4) 2 k (0.4) n (0.4) n ( 2 n 1 1)
k 0 12
The lower limit on the convolution sum simplifies to k = 0 [because u(k) = 0, k < 0], the upper limit
to k = n [because u(n – k) = 0, k > n], and we get
n n
y(n) = ak a n k an 1 ( n 1) a n u ( n )
k 0 k 0
1, n 1
( d ) Given that, x ( n ) = u ( n 1)
0 , n 1
n , n 1
and h ( n ) = n u (n 1)
0, n 1
By definition, y(n) = x ( n ) h( n ) x ( m ) h( n m )
m
1, m 1
Here x ( m ) = u ( m 1)
0 , m 1
n m , n m 1 n m , m n 1
and h(n – m) = n m u ( n m 1)
0 , n m 1 0 , m n 1
using the above definitions we can see that, x(m) h(n – m) will be zero for m < 1 and m > n1, and non
zero for m n 1 and m 1 ,thus n has to be greater than equal to 2
n 1 n 1
and y(n) = x( m ) h( n m ) = a n m
m 1 m 1
n 1 n 2 n 2
n n
= am n ( m 1) n 1 ( m ) 1
m 1 m 0 m 0
n
, n2
Therefore, y(n) = 1
0, n0
n
or, more compactly, y(n) = u ( n 2)
1
n , n 0
( e ) Given that, x ( n ) = n u( n )
0 , n 0
n a n , n 1
and h ( n ) = a u ( n 1)
0 , n 1
By definition, y ( n ) = x ( n ) h( n ) x ( k ) h( n k )
k
k , k 0
Since, x ( k ) = k u( k )
0 , k 0
( n k ) a ( n k ) , n k 1 a ( n k ) , k n 1
and h(n – k) = a u ( n k 1)
0 , n k 1 0 , k n 1
Thus product x(k) h(n – k) will be non zero for k n 1 and k 0 , that is 0 k n – 1 and n 1. Thus
n 1 n 1
y(n) = x ( k ) h( n k ) = k a ( n k ) ;
k 0 k 0
n 1
n a n 1
= a k ak
(1 a )2
[1 n a n 1 ( n 1)a n ]
k 0
a n 1
[1 n a n 1 ( n 1) a n ], n 1
Therefore, y ( n ) = (1 a )2
0; n 1
a n 1
or, more compactly, y(n) = 2
[1 n a n 1 ( n 1) a n ] u ( n 1)
(1 a )
Example 2.10
Consider a LTI system with impulse response h[n] = u[n] and input x[n] then, output y[n] will be _____.
Solution 2.10
Since, y[n] = x[n] h[n]
=
m
x[m ] h[n m ]
1, n m 0
h[n – m] = u[n – m] =
0 , n m 0
1, m 0
h[n – m] =
0 , m 0
n
Thus, y[n] =
m
x [m ]
Example 2.11
Fin d:
( a ) x[n] [n – n 0 ]
( b ) x[n] [an]
( c ) x[n] [an + b]
Solution
n
(a) x[n] [n – n0] =
m
x[m ] [n m n 0 ]
Since, x[n] =
m
x[m ] [n m ]
b
[an + b] = n if b/a is integer and
a
if b/a is not integer then [an + b] = 0.
b
So, if b/a is integer then x[n] [an + b] = x n .
a
Otherwise, x[n] [an + b] = 0.
Study Notes
Suppose, y(n) = x(n) h(n),
to find the result we can perform,
x(n) h(n) = x ( k ) h( n k )
k
or h(n) x(n) = h( k ) x ( n k )
k
If we have signal x(n) simpler than h(n) then we perform h(n) x(n) and if h(n) is simpler then we perform
x(n) h(n).
For example if, h(n) = an u(n) and x(n) = u(n)
then to find y(n) it will be easier to perform,
h(n) x(n) = h( k ) x ( n k )
k
Study Notes
To perform the convolution of a signal with delta function, we should follow these steps
• First of all find the location of delta function if location is an integer then move to step 2 else answer is zero.
• Now convolution of signal x(n) with given delta fuction is found by replacing variable ‘n’ in x(n) by
(n – location), where location is location of the delta function.
REMEMBER If h[n] is impulse response of a LTI system then step response will be h[n] u[n]
n
[n] = h[n ] u [n ] h[ m ]
and from [n] we can get h[n] by
h[n] = [n] – [n – 1]
Sy = y( n ) x( k ) h( n k )
n n k
2.6 Relation between LTI System Properties and the Impulse System
The impulse response completely characterizes the input-output behaviour of an LTI system. Hence, the
properties of the system, such as memory, causality, and stability, are related to the system’s impulse response.
= ... + h(–2) x(n + 2) + h(–1) x(n + 1) + h(0) x(n) + h(1) x(n – 1) + h(2) x(n – 2) + ....
For the system to be memoryless, y(n) must depend only on x(n) and therefore cannot depend on x(n – k) for
k 0. Hence, every term in the above equation must be zero, except h(0) x(n). The condition implies that h(k) = 0 for
k 0, thus a discrete-time LTI system is memoryless if and only if
h(n) = K (n) ...(26)
where K = h(0) is a constant and the convolution sum reduces to the relation,
y(n) = K x(n)
If a discrete-time LTI system has an impulse response h(n) that is not identically zero for n 0, then the
system has memory.
From equation (4), we can deduce similar properties for continuous-time LTI systems with an without
memory. In particular, a continuous-time LTI system is memoryless if h(t) = 0 for t 0, and such a memoryless
LTI system has the form,
y ( t ) = K x(t)
for some constant K and has the impulse response,
h ( t ) = K (t) ...(27)
y(t) = h( ) x ( t ) d
0
x ( ) h( t ) d
= h( k ) x( n k )
k k
h( k ) x ( n k )
If we assume that the input is bounded, or x( n ) Bx , then x(n k ) Bx , and it follows that,
y( n ) Bx
k
h( k )
From the above equation, we can conclude that if the impulse response is absolutely summable, i.e. if
k
h( k ) <
y( t ) = h ( ) x ( t ) d h() x (t d Bx h ( ) d
Therefore, the system is stable if the impulse response is absolutely integrable, i.e., if
h( ) d <
y ( t)
x (t ) h( t) hI( t) x (t )
y ( n)
x (n ) h( n) hI( n) x (n )
The relationship between the impulse response, h(t), and that of the corresponding inverse system hI(t), is
easily derived. The impulse response of the cascade connection in Fig. 2.11 is the convolution of h(t) and hI(t). We
require the output of the cascade to equal the input, or
x(t) [h(t) hI(t)] = x(t)
This requirement implies that,
h(t) hI(t) = (t) ...(32)
Similarly, the impulse response of a discrete-time LTI system, hI(n), must satisfy
h(n) hI(n) = (n) ...(33)
REMEMBER In discrete-time LTI system also we can have series/cascade connection and parallel connection
as that of continuous time LTI system.
h[t ]
x [n] h1[n ] h2[n ] y [n ] x [n] = h1[n ] h2 [n]
y [n ]
h1[n]
+ h[n ]
x [t ] y [n ] x [n ] = h1[n ] h2[n ]
y [n ]
–
h2[n]
Rxy() = x ( t ) y ( t ) dt
x( t ) y ( t ) dt ...(34)
The index is the (time) shift (or lag) parameter and the subscripts xy in the cross-correlation function Rxy()
indicate the signals being correlated. The order of the subscripts, with x proceeding y, indicates the direction in which
one signal is shifted relative to the other.
If we reverse the roles of x(t) and y(t) in equation (36) and therefore reverse the order of the indices xy, we
obtain the cross-correlation function,
Ryx() = y( t ) x ( t ) dt y( t ) x( t ) dt ...(37)
Study Note
Therefore, Rxy(t) is simply the folded version of Ryx(t), where the folding is done w.r.t. t = 0. Hence, Rxy(t) provides the same
information as Ryx(t) w.r.t. the similarly of x(t) to y(t).
The similarities between the computation of the cross-correlation of two signals and the convolution integral
of two signals is apparent. For energy signals, there is a simple mathematical relationship between correlation and
convolution integral.
Rxy() = x() y(–) ...(39)
Proof:
Rxy() = x() y(–)
= x ( t ) y( ( t )) dt x ( t ) y( t ) dt
Remember: For real signals cross correlation is defined as Rxy() = x() y(–) and for non real signals Rxy() = x() y*(–).
Similarly for real signals Ryx() = y() x(–) and for non real signals Ryx() = y() x*(–).
An important special case of correlation of power signals is the correlation between two periodic signals
whose fundamental periods are such that the product of the two signals is also periodic. This will happen any
time the fundamental periods of the two periodic signals have an LCM.
For two signals whose product has a period T, the general for of the correlation function (for real power
signals) is given by,
T /2
1
R xy () = Lim x ( t ) y( t ) dt
T T T / 2
T /2
1
can be replaced by R xy () = x ( t ) y( t ) dt ...(42)
T T / 2
Remember: for real signal auto correlation is Rxx() = x() x(–) and for non real signal Rxx() = x() x*(–)
If x(t) is a power signal, its autocorrelation is
T /2 T /2
1 1
R xx () = Lim x ( t ) x ( t ) dt Lim x ( t ) x ( t ) dt ...(44)
T T T / 2 T T
T / 2
Proof: R xx () = x ( t ) x( t ) dt
R x x ( ) 0 = x( t ) x( t ) dt
2
Rxx(0) = x
( t ) dt Ex
T /2
1
Proof: Rxx() = Lim x ( t ) x ( t ) dt
T T T / 2
T /2
R x x ( ) 0 = Lim 1 x ( t ) x ( t ) dt
T T
T / 2
T /2
1
Rxx(0) = Lim x 2 ( t ) dt Px
T T T / 2
T /2
1
Rxx( – T) = x ( t ) x ( t ( T )) dt
T T / 2
T /2 T /2
1 1
= x ( t ) x ( t T ) dt x ( t ) x ( t ) dt R x x ( )
T T / 2 T T / 2
R xy () = x( t ) y( t ) dt
R x y ( ) = R x y (0) x( t ) y( t ) dt 0
0
The index m is the (time) shift (or lag) parameter and the subscripts xy on the cross-correlation function
Rxy(m) indicate the signals being correlated. The order of the subscripts with x preceeding y, indicates the
direction in which one signal is shifted relative to the other.
If we reverse the roles of x(n) and y(n) equation (52) and therefore reverse the order of the indices xy,
we obtain the cross-correlation function,
R yx( m ) = y( n ) x ( n m ) y( n m ) x ( n )
n n
...(53)
Study Note
Therefore, Rxy(m) is simply the folded version of Ryx(m), where the folding is done w.r.t. m = 0. Hence, Rxy(m)
provides the same information as Ryx(m) w.r.t. the similarity of x(n) to y(n).
The similarities between the computation of the cross-correlation of two signals and the convolution
sum of two sequences is apparent. For energy signals, there is a simple mathematical relationship between
correlation and convolution sum:
R x y ( m ) = x(m) y(–m) ...(55)
Remember: For real signals the crosscorrelation is Rxy(m) = x(m) y(–m) and for non real signals crosscorrelation
is Rxy(m) = x(m) y*(–m). Similarly for real signals the crosscorrelation is Ryx(m) = y(m) x(–m) and for non real
signals crosscorrelation is Ryx(m) = y(m) x*(–m).
Proof:
R x y ( m ) = x(m) y(–m)
= x( n ) y( ( m n )) x( n ) y( n m )
n n
1
= Lim
N 2N 1
x( n m ) y ( n )
n N
1
= Lim x ( n m ) y( n )
N 2N 1 n N
An important special case of correlation of power signals is the correlation between two periodic signals
whose fundamental periods are such that the product of the two signals is also periodic. This will happen any
time the fundamental periods of the two periodic signals have an LCM.
For two signals whose product has a period N, the general form of the correlation function (for real
power signals).
N
1
R x y ( m ) = Lim
N
x (n) y (n m)
2N 1 nN
can be replaced by,
N
1
R xy( m ) = x ( n ) y ( n m)
2 N 1 n N
...(58)
N
1
= Lim
N 2N 1
x( n m ) x( n )
n N
Remember: For real signal autocorrelation is Rxx(m) =x(m) x(–m) and for non real signal autocorrelation will
be R xx(m) =x(m) x*(–m)
Proof: By definition,
Rxx(m) = x( n ) x ( n m )
n
= Rxx(–m)
2. Relation to signal energy and signal power.
Case-1: If x(n) is an energy signal, then
2
Rxx(0) = Ex x
n
( n) ...(62)
Proof: By definition,
Rxx(m) = x( n ) x ( n m )
n
2
R x x ( m ) m 0 R x x (0) =
n
x( n ) x( n ) x
n
( n ) Ex
Proof: By definition,
N
1
Rxx(m) = Lim
N 2N 1
x( n ) x( n m )
n N
N
1
R x x ( m ) m 0 R x x (0) = Lim
N 2N 1
x ( n ) x( n )
n N
N
1
= Lim x 2 ( n ) Px
N 2N 1 n N
Example 2.12
Let h(t) be the triangular pulse shown in Fig. 2.12 (a) and let x(t) be the unit impulse train Fig. 2.12 (b)
expressed as:
x( t ) T ( t ) ( t nT )
n
Determine and sketch y(t) = h(t) x(t) for the following values of T:
(a) T = 3, (b) T = 2, (c) T = 1.5.
h( t) T (t )
1 1
t t
–1 0 1 –2T –T 0 T 2T
(a) (b)
Fig. 2.12
Solution 2.12
y ( t ) = h(t) (t)
= h( t ) ( t nT ) h( t ) ( t nT )
n n
= h ( t nT )
n
y(t)
T=3
1.0
t
–7 –6 –5 –4 –3 –2 –1 0 1 2 3 4 5 6 7
(a)
y(t)
T=2
1.0
t
–7 –6 –5 –4 –3 –2 –1 0 1 2 3 4 5 6 7
(b)
y(t)
T = 1.5
t
–7 –6 –5 –4 –3 –2 –1 0 1 2 3 4 5 6 7
(c)
which is sketched in Fig. 2.13 (c). Note that when T < 2, the triangular pulses are no longer separated
and they overlap.
Example 2.13
The system shown in figure formed by connecting two systems in cascade. The impulse
responses of the systems are given by h 1 (t) and h 2 (t) respectively and
h 1 ( t ) = e –2t u(t);
h 2 ( t ) = 2e –t u(t)
( a ) Find the impulse response h(t) of the overall system shown in figure.
w( t )
x (t ) h1 ( t ) h2( t) y ( t)
(a)
x (t ) h( t ) y ( t)
(b)
1
= 2 e d e2 d 2 1 1
0 0 2
Example 2.14
1
1
t
0 1 2 3
t
0 1 2 3 4 –1
(a) (b)
Example 2.15
(b) From figure, we see that, h(t) 0 for t < 0. Hence, the system is not causal.
h( t )
1/T
t
–T/2 0 T/2
Example 2.16
= k n k u[k ] u [n k ]
k
1 0 k n
Since, u[k] u[n – k] =
0 oth erwise
n n k
we have, y[n] = k n k n ,n0
k 0 k 0
n 1 ( / ) n 1
u [n ]
Thus , y[n] = 1 ( / )
n
( n 1) u[n ]
1
( n 1 n 1 ) u[n ]
or, y[n] =
n ( n 1) u[n ]
(b) y[n] = x k ] h[n k ] k u[k ] ( n k ) u [( n k )]
k k
= n 2 k u[k ] u[k n ]
k
For n 0, we have
1 nk
u[k] u[k – n] =
0 oth erwise
n 2k
n
Thus, we have y[n] = ,n 0
k 0 1 2
For n greater than zero we have
1 nk
u[k] u[k – n] =
0 oth erwise
n 2k n
2 n
y[n] =
k n 1 2
n
= ; n 0
1 2
n
we obtain, y[n] = ; ( All n )
1 2
which is sketched in figure.
y [n]
n
–2 –1 0 1 2 3
Example 2.17
Consider a discrete-time LTI system whose input x[n] and output y[n] are related by,
y[n] = 2 k n x[k + 1]
k=
Is the system causal?
Solution: 2.17
By the definition the impulse response of system is given by
n n n
h[n] = 2 k n [k 1] 2 ( n 1) [k 1] = 2 ( n 1) [k 1]
k k k
From equation we have h[–1] = u[0] = 1 0. Thus, the system is not causal.
Example 2.18
Consider a discrete-time LTI system with impulse response h[n] given by,
h[n] = n u[n]
( a ) Is this system causal?
( b ) Is this system BIBO stable?
Solution: 2.18
( a ) Since h[n] = 0 for n < 0, the system is causal.
( b ) For BIBO stability,
h[k ] = k u[n ]
k k
k 1
=
1
( 1)
k
Example 2.19
Find the impulse response h[n] for each of the causal LTI discrete-time systems satisfying the following
difference equations and indicate whether each system is an FIR or an IIR system.
(a) y[n] = x[n] – 2x[n – 2] + x[n – 3]
( b ) y[n] + 2y[n – 1] = x[n] + x[n – 1]
1
( c ) y[ n] y [ n 2] 2 x[ n ] x[ n 2]
2
Solution 2.19
( a ) By definition,
h[n] = [n] – 2[n – 2] + [n – 3]
or, h[n] = {1, 0, –2, 1}
Since, h[n] has only four terms, the system is an FIR system.
(b) h[n] = –2h[n – 1] + [n] + [n – 1]
Since, the system is causal, h[–1] = 0.
Then, h[0] = –2h[–1] + [0] + [–1]
= [0] = 1
h[1] = –2h[0] + [1] + [0]
= –2 + 1 = –1
h[2] = –2h[1] + [2] + [1]
= –2(2) = –2 2
h[n] = –2h[n – 1] + [n] + [n – 1] = (–1)n 2n – 1
Hence, h[n] = [n] + (–1)n 2n – 1 u[n – 1]
Since, h[n] has infinite terms, the system is an IIR system.
1
(c) h[n] = h[n 2] 2 [n ] [n 2]
2
Since, the system is causal,
h[–2] = h[–1] = 0, Then
1
h[0] = h[2] 2 [0] [2] 2 [0] 2
2
1
h[1] = h[1] 2 [1] [1] 0
2
1 1
h[2] = h[0] 2 [2] [0] (2) 1 0
2 2
1
h[3] = h[1] 2 [3] [1] 0
2
Hence, h[n] = 2[n]
Since, h[n] has only one term, the system is an FIR system.
5 . Let, y(t) = e t u ( t ) ( t 3 k )
k
t
x [n ] h1[n] h2[n ] h2[n] y [n ] ( t )
y(t) = e x ( 2) d
(a)
11
What is the impulse response h(t) for this system?
10 (b) Determine the response of the system when the
8 input x(t) is as shown in figure.
5
4 x ( t)
1 1
n
–1 0 1 2 3 4 5 6 7 1
(b)
t
(a) Find the impulse response h1[n]. –1 2
(b) Find the response of the overall system to the
input, 1 5 . Suppose that the signal,
x[n] = [n] – [n – 1] x( t) = u(t + 0.5) – u(t – 0.5)
is convolved with the signal,
1 1 . The following are the impulse responses of discrete-
time LTI systems. Determine whether each system h(t) = e j0 t
is causal and/or stable. (a) Determine a value of 0 which ensures that,
n
y(0) = 0
1
(a) h [n ] u[n ] where, y(t) = x(t) h(t)
5
(b) Is your answer to the previous part unique?
(b) h[n] = (0.8)n u[n + 2]
n
1 6 . Determine whether each of the following statements
1 concerning LTI systems is true or false. Justify your
(c) h [n ] u [ n ]
2 answers.
(d) h[n] = (5)n u[3 – n] (a) If h(t) is the impulse response of an LTI system
1
n and h(t) is periodic and non-zero, the system is
(e) h [n ] u [n ] (1.01) n u [n 1] unstable.
2
(b) The inverse of a causal LTI system is always
n
1 causal.
(f) h [n ] u [n ] (1.01) n u [1 n ]
2
(c) If h[n ] K for each n, where K is a given
n
1
(g) h [n ] n u[n 1] number, then the LTI system with h[n] as its
3 impulse response is stable.
(d) If a discrete-time LTI system has an impulse
1 3 . The following are the impulse response of continuous
response h[n] of finite duration, the system is
time LTI systems. Determine whether each system
stable.
is causal and/or stable.
(e) If an LTI system is causal, it is stable.
(a) h(t) = e–4t u(t – 2)
(b) h(t) = e–6t u(3 – t) (f) The cascade of a non-causal LTI system with a
(c) h(t) = e–2t u(t + 50) causal one is necessarily non-causal.
(d) h(t) = e2t u(–1 – t) (g) A continuous-time LTI system is stable if and
6 t
only if its step response s(t) is absolutely
(e) h( t ) e integrable – that is, if and only if
(f) h(t) = t e –t u(t)
(g) h(t) = (2e–t – e(t – 100)/100) u(t) s( t ) dt
1 4 . (a) Consider an LTI system with input and output
(h) A discrete-time LTI system is causal if and only
related through the equation,
if its step response s[n] is zero for n < 0.
1 7 . Consider the cascade of two systems shown in figure. Show that the commutativity property does not
The first system, A is known to be LTI. The second hold for these two systems by computing the
system, B is known to be the inverse of system A. Let impulse response of the cascade combination in
y1(t) denote the response of system A to x1(t), and figure (a) and (b), respectively.
let y2(t) denote the response of system A to x2(t).
System System
x [n] y [n ]
LTI y ( t) System A B
x ( t) x (t )
system A B
(a)
(a) What is the response of system B to the input
ay1(t) + by2(t), where a and b are constants? System System
x [n] y [n ]
B A
(b) What is the response of system B to the input
y1(t – )? (b)
1 8 . In the text, we saw that the overall input-output (b) Suppose that we replace system B in each of the
relationship of the cascade of two LTI systems does interconnected systems of figure by the system
not depend on the order in which they are cascaded. with the following relationship between its input
This fact, known as the commutativity property, w[n] and output z[n];
depends on both the linearity and the time z[n] = w[n] + 2
invariance of both systems. In this problem, we Repeat the calculations of part (a) in this case.
illustrate the point.
(a) Consider two discrete-time systems A and B,
where system A is an LTI system with unit sample
response h[n] = (1/2)n u[n]. System B, on the
other hand, is linear but time varying.
Specifically, if the input to system B is w[n], its
output is
z[n] = n w[n]