DSP Unit I Material
DSP Unit I Material
UNIT-I
Discrete Time Signals and Systems
Discrete Time Fourier Transform
x(t)
0 t
y(t)
t
0
x(n)
012 34 5 6 78 9………………………………………… n
y(n)
012 34 5 6 78 9………………………………………… n
Amplitude restricted version of discrete time signals are called digital signals, for which the
different number of amplitudes are restricted to finite number (Two). y(n) is a digital signal and all
digital signals are discrete time signals. Digital signals can be obtained from discrete time signals by
quantization mechanism.
Examples:
➢ x(t) = 2Cos(3t ) + 3Sin(2t )
➢ x ( n) = 2 n
➢ y(n) = 2Cos(3n − 4) Discrete time signals
➢ z(n) = 2e jn / 3
➢ x(n) = {1,1,0,1,0,1,1,1,0,0,1}
➢ y(n) = {1,−1,1,−1,1,1,−1,−1,1} Discrete time or Digital signals
➢ z(n) = {1,1,2,1,2,1,1,1,2,2,1}
(n)
n
.…-3 -2 -1 0 1 2 3….
u(n)
n
…..... -3 -2 -1 0 1 2 3 4 5 6 7…………
Unit step signal u(n) is the sum of a train of unit sample sequences
u(n) = (n) + (n – 1) + (n – 2) + (n – 3) + ……….+ (n – k)+ ……
n
u(n) = δ(n - k) or u(n) = δ(k)
k=0 k = −
Unit sample sequence (n) is difference between u(n) and u(n-1)
(n) = u(n) – u(n – 1)
A sequence x(n) can be represented using unit sample sequence (n)
x(n) = x(k) δ(n - k)
k = −
n ; n0
r(n) = nu(n) =
0 ; n0
r(n)
7
6
5
4
3
2
1
n
……. -4 -3 -2-1 0 1 2 3 4 5 6 7…………
x(n) = a nu(n),0 a 1
x(n)
n
……… -4 -3 -2 -1 0 1 2 3 4 5 6 7…………
y(n) = a n ,0 a 1
y(n)
n
…… .. -9-8-7-6 -5 -4 -3 -2-1 0 1 2 3 4 5 6 7…………
n
……..-8 -7 -6 -5 -4 -3 -2 -1 0 1 2 3 ……….
y(n) = an , a 1
y(n)
n
…… .. -9-8-7-6 -5 -4 -3 -2-1 0 1 2 3 4 5 6 7 8…………
n
… …. -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7…………
n
y(n) = a , a 1
x(n)
n
… …. -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7…………
n
0
Functional representation of unit sample sequence x(n) = (n)
1 ; n=0
δ(n) =
0 ; n0
Tabular representation of unit sample sequence x(n) = (n)
n … … -3 -2 -1 0 1 2 3 … …
(n) 0 0 0 0 0 1 0 0 0 0 0
n
0 1 2 3 4 5 6 7…………
Functional representation of unit step sequence x(n) = u(n)
1 ; n0
u(n) =
0 ; n0
x(n - 2) = 6,− 2,3,1,4,5,−1,7 , it is shifted to right by 2 units.
x(n + 3) = 6,−2,3,1,4,5,− 1 ,7 , it is shifted to left by 3 units.
Example-2:
Given discrete time signal, x(n) = u(n)
u(n)
n
0 1 2 3 4 5 6 7…………
u(n-2)
n
0 1 2 3 4 5 6 7 8 9……
u(n+2)
1
n
-2 -1 0 1 2 3 4 5 6
7…………
Prof.B.Ramesh Reddy, Dept. of ECE, Coordinator-IQAC, LBRCE, Mylavaram. Page No. 9
(B)Time Scaling Operation:
If the time scaling operation is applied on a discrete time signal x(n), then the signal is compressed
or expanded in time axes without changing its amplitude. It is represented with x(an), where ‘a’ is
time scaling parameter. If a>1, then it is compressed and if 0<a<1, then it is expanded signal.
Example:
Given discrete time signal, x(n) = 1,2,3, 4,5,6,7,8,9
x(2n) = 2, 4,6,8 , it is compressed signal.
x(n/2) = 1,0,2,0,3,0, 4,0,5,0,6,0,7,0,8,0,9 , it is expanded signal.
(C)Time Reversal or Folding Operation:
Time reversal signal can be obtained by interchanging left hand side and right hand side samples
with respect to vertical axes or y-axes. If the given discrete time signal is x(n), then its time reversal
form is represented with x(-n) and the operation is called folding or mirror image or time reversal.
Example:
Given discrete time signal, x(n) = 1,2,3, 4,5,6,7,8,9
Then, time reversal signal x(-n) = 9,8,7,6,5, 4,3,2,1
(D)Amplitude Scaling Operation:
If the amplitude scaling operation is applied on a discrete time signal x(n), then the signal amplitude
may increase or decrease without changing its duration. It is represented with Ax(n), where ‘A’ is
amplitude scaling parameter.
Example:
Given discrete time signal, x(n) = 1,2,3, 4,5,6,7,8,9
y(n) = 2x(n) = 2,4,6, 8 ,10,12,14,16,18
1
z(n) = x(n) = 0.5,1,1.5, 2,2.5,3,3.5,4,4.5
2
(c)y(n)+ z(n) = 27,21,18,12,6, 0,6,12,18,21,27
(E)Convolution Operation:
Convolution is an operation, which is used in almost all signal processing applications to analyze
signals and systems in both the time and frequency domain. Convolution is a special operation,
which includes four different operations, namely
➢ Folding,
➢ Shifting,
➢ Multiplication and
➢ Summation in the case of discrete time signals or
Integration in the case of continues time signals.
(a)Graphical method
➢ Draw the graphical representation of given sequences x1(m) and x2(m).
➢ Take the folding form of x2(m) to get x2(-m).
➢ Shift the folding sequence x2(-m) in different cases to get x2(n-m).
➢ Finally apply convolution formula to get the convoluted sequence.
x(n) = x1(n) x 2 (n) = x (m) x (n − m)
m = −
1 2
7
n = 0 x(0) = x (m) x (−m) = 1 5 = 5
m=0
1 2
7
n = 1 x(1) = x (m) x (1− m) = 1 6 + 2 5 = 6 + 10 = 16
m=0
1 2
7
n = 2 x(2) = x (m) x (2 − m) = 1 7 + 2 6 + 3 5 = 7 + 12 + 15 = 34
m=0
1 2
7
n = 3 x(3) = x (m) x (3 − m) = 1 8 + 2 7 + 3 6 + 4 5 = 8 + 14 + 18 + 20 = 60
m=0
1 2
7
n = 4 x(4) = x (m) x (4 − m) = 1 9 + 2 8 + 3 7 + 4 6 = 9 + 16 + 21 + 24 = 70
m=0
1 2
7
n = 5 x(5) = x (m) x (5 − m) = 2 9 + 3 8 + 4 7 = 18 + 24 + 28 = 70
m=0
1 2
7
n = 6 x(6) = x (m) x (6 − m) = 3 9 + 4 8 = 27 + 32 = 59
m=0
1 2
7
n = 7 x(7) = x (m) x (7 − m) = 4 9 = 36
m=0
1 2
Convoluted Sequence x(n) = x1(n) x2(n) = {5, 16, 34, 60, 70, 70, 59, 36}
x2(m)
x2(-m)
Folding:
Shifting:
Case –1: x2(1-m)
-9-8-7-6-5-4-3-2-1 0 1 2 3 4 5 6 7 8 9 m
m -4 -3 -2 -1 0 1 2 3 4
Given x1(m) 1 2 3 4
Sequences x2(m) 5 6 7 8 9
Folding x2(-m) 9 8 7 6 5 0 0 0
x2(1-m) 9 8 7 6 5 0 0
x2(2-m) 9 8 7 6 5 0
x2(3-m) 9 8 7 6 5
Shifting x2(4-m) 9 8 7 6
x2(5-m) 0 9 8 7
x2(6-m) 0 0 9 8
x2(7-m) 0 0 0 9 Sum
x1(m) x2(-m) 5 0 0 0 = 5
x1(m) x2(1-m) 6 10 0 0 = 16
x1(m) x2(2-m) 7 12 15 0 = 34
x1(m) x2(3-m) 8 14 18 20 = 60
Multiplication
x1(m) x2(4-m) 9 16 21 24 = 70
x1(m) x2(5-m) 0 18 24 28 = 70
x1(m) x2(6-m) 0 0 27 32 = 59
x1(m) x2(7-m) 0 0 0 36 = 36
Convoluted Sequence x(n) = x1(n) x2(n) = {5, 16, 34, 60, 70, 70, 59, 36}
5 6 7 8 9 0 0 0
0 5 6 7 8 9 0 0
0 0 5 6 7 8 9 0
0 0 0 5 6 7 8 9
= 1 2 3 4 0 0 0 0
9 0 0 0 5 6 7 8
8 9 0 0 0 5 6 7
7 8 9 0 0 0 5 6
6 7 8 9 0 0 0 5
= 1 5 1 6 + 2 5 1 7 + 2 6 + 3 5 1 8 + 2 7 + 3 6 + 4 5
1 9 + 2 8 + 3 7 + 4 6 2 9 + 3 8 + 4 7 3 9 + 4 8 4 9
= 5 6 + 10 7 + 12 + 15 8 + 14 + 18 + 20
9 + 16 + 21 + 24 18 + 24 + 28 27 + 32 36
= 5 16 34 60 70 70 59 36
Convoluted Sequence x(n) = x1(n) x2(n) = {5, 16, 34, 60, 70, 70, 59, 36}
Example-2:
Compute the convoluted sequence x(n) = x1(n) x2(n), if (i) b=a. (i) b≠a
Given x1(n)=anu(n) and x2(n)=bnu(n).
x(n) = x1(n) x2(n)
= x (m) x (n − m)
m = −
1 2
= a
m = −
m
u(m) b n-m u(n-m) u(m) = 1, m 0 and
n
u(n-m) = 1, n – m 0 or m n
= a mb n-m
m =0 u(m)u(n-m) = 1, 0 m n
m
a
n
= bn
m =0 b
= (n + 1)a n , n 0
Case-2: if b≠a
x(n) = x1(n) x2(n)
m
a n
= b n
m =0 b
n +1
a
1−
= bn b
a
1−
b
n +1 n +1
b −a
= ,n 0
b−a
Examples:
➢ x(n) = Cos(0.125n) is an even signal.
➢ y(n) = Sin(0.125n) is an odd signal.
➢ z(n) = Cos(0.125n) + Sin(0.125n) is neither even nor odd signal.
➢ Determine the even and odd parts of following signals
j n
(1)x(n)=2n (2) y(n) = 3e 5 (3) z(n) = 1,2,3, 4,5,6,7,8,9
Problem-1:
Given discrete time signal x(n)=2n
x(n) + x( −n) 3n + 3− n
Even part of the signal x e (n) = =
2 2
x(n) − x( −n) 3n − 3− n
Odd part of the signal x o (n) = =
2 2
Problem-2:
j n
Given discrete time signal y(n) = 3e 5 = 3Cos n + j 3Sin n
5 5
−j n
5
y(-n) = 3e = 3Cos n − j3Sin n
5 5
y(n) + y( −n)
Even part of the signal ye (n) = = 3Cos n
2 5
y(n) − y(−n)
Odd part of the signal yo (n) = = j3Sin n
2 5
z(-n) = 9,8,7,6,5, 4,3,2,1
Even part of the signal
z(n) + z( −n) 1
ze (n) = = 9,8,8,8,8, 8 ,8,8,8,8,9 = 4.5,4,4,4,4, 4,4,4,4,4,4.5
2 2
Odd part of the signal
z(n) − z( −n) 1
zo (n) = = − 9,−8,−6,−4,−2, 0,2,4,6,8,9 = − 4.5,−4,−3,−2,−1, 0,1,2,3,4,4.5
2 2
2
Example-1: x(n) = 9Cos n + 3
9
If x(n) is periodic with a period of N samples, then
2 2 2 2
x(n + N) = 9Cos (n + N ) + 3 = 9Cos n+ N + 3 = 9Cos n + 3 + 2k = x(n)
9 9 9 9
2
N = 2k N = 9k
9
Possible integer values of N = 9,18,27,…. for integer values of k = 1,2,3….
Hence, given signal is periodic with a period N=9 samples.
5
Example-2: x(n) = 9Cos n + 3
9
If x(n) is periodic with a period of N samples, then
5 5 5 5
x(n + N) = 9Cos (n + N ) + 3 = 9Cos n + N + 3 = 9Cos n + 3 + 2k = x(n)
9 9 9 9
5 18
N = 2k N = k
9 5
Possible integer values of N = 18,36,54,…. for integer values of k = 5,10,15….
Hence, given signal is periodic with a period N=18 samples.
5
Example-3: x(n) = Sin n +
9
If x(n) is periodic with a period of N samples, then
5 5 5 5
x(n + N) = Sin (n + N ) + = Sin n + N + = Sin n + + 2k = x(n)
9 9 9 9
5 18
N = 2k N = k , Integer value of N is not possible for integer values of k.
9 5
Hence, given signal is aperiodic or non-periodic.
Example-5: x(n) = 3Sin n2
8
If x(n) is periodic with a period of N samples, then
x(n + N) = 3Sin (n + N )2 = 3Sin n2 + N 2 + Nn = 9Cos n2 + 2k + 2k n = x(n)
8 8 8 4 8 1 2
N 2 = 2k1 N 2 = 16k1 N = 4 k1
8
Possible integer values of N = 4,8,12,16… for integer values of k1 = 1,4,9,16…. and
N = 2k2 N = 8k2
4
Possible integer values of N = 8,16… for integer values of k 2 = 1,2….
N=Common & Min (N of k1,N of k2)=8, hence, given signal is periodic with a period N=8 samples.
Example-6: What is the time period of a signal x(n), which is the sum of 3 periodic signals with
periods N1=2, N2=3 and N3=5.
Possible integer values if k = 1,2,3,4,5,6,7,8,9,….
N1k=2,4,6,8,10,12,14,16,18,20,22,24,26,28,30,32,34……
N2k=3,6,9,12,15,18,21,24,27,30,33………
N3k=5,10,15,20,25,30,35,40,45…….
N=Common & Min (N1,N2,N3)=30, hence, given signal is periodic with a period N=30 samples.
2 2
1 N
E=
n = −
x(n) & P = Lt
N→ 2N + 1
x(n)
n= −N
n
1
Example-1: x(n) = 3 u (n)
2
Total Energy:
2
1
n 2
1
n 1 4
E= x(n) = 3 u (n) = 9 = 9 = 9 = 12 jouls
2 n=0 4
1−
4
n = − n = −
1 3
Average Power:
n 2 n
2
1 N 1 N 1 1 N 1
P = Lt x(n) = Lt 3 u (n) = Lt 9
N→ 2N + 1 N→ 2N + 1 N→ 2N + 1
n = −N n = −N 2 n = 0 4
= Lt
9 1 − 4
1
N +1
( )
12
= Lt ( ) N +1
= 12 (1 − 0) = 12 = 0
1 − 1 4
N→ 2N + 1
1 − 14 N→ 2N + 1
Given x(n) is energy signal because the total energy is finite and average power is zero.
Average Power:
2
1 N 1 N 1 N N+1 1
P = Lt x(n) = Lt u(n) = Lt 1 == Lt = Watts
2
Given x(n) is power signal because the average power is finite and total energy is infinity.
Example-4: Determine the energy and power of x(n) = n(− 1)n , n = 1,2,3 and y(n) = x(n + 7k )
k = −
Given x(n) = n(− 1)n , n = 1,2,3 = 0,−1,2,−3 =
2
Average Power:
( )
2
1 N 1 14 14
Px = Lt x(n) = Lt − 1 + 2 + − 3 = Lt = =0
2 2 2
Given x(n) is energy signal because the total energy is finite and average power is zero.
Given y(n) = x(n + 7k ) = .. + x(n − 7) + x(n) + x(n + 7) + .. periodic with a period of 7 samples
k = −
y(n) = .......... ..., −1,2,−3,0,0,0, 0,−1,2,−3,0,0,0,0,−1,2,−3,0,0,0,0,−1,2,−3,.......
2
n = −
Average Power:
1 N
2 0 + − 1 2 + 2 2 + − 3 2 + 0 + 0 + 0 14
Py = Lt y(n) = = =2
N→ 2N + 1
n= −N 7 7
Given y(n) is power signal because the average power is finite and total energy is infinity.
Examples:
1 ; n=0
➢ Digital Impulse Signal or Unit Sample Sequence, δ(n) =
0 ; n0
1 ; n0
➢ Unit Step Signal, u(n) =
0 ; n0
n ; n0
➢ Unit Ramp Signal, r(n) = nu(n) =
0 ; n0
n n
➢ Double Exponential Signals, x(n) = a ,0 a 1 and y(n) = a , a 1
Signals whose characteristics are random in nature are called nondeterministic signals or random
signals, where the mathematical representation is not possible, for example noise signal.
Example:
N(t)
Continuous Time or
Input signal, x(t) Analog System output signal, y(t)
➢ Thermal stability of continuous time or analog systems is poor because of all analog
components are temperature sensitive.
➢ Continuous time or analog systems are non programmable and static in nature.
➢ Continuous time or analog systems are described by differential equation, which involves
only differentials.
Examples:
(a) A simple RC high pass filter acts as a differentiator, where the output y(t) is the
differentiation of input x(t)
d
y(t) = [ x(t) ]
dt
(b) A simple RC low pass filter acts as an integrator, where the output y(t) is the integration
of input x(t)
d
y(t) = x(t) x(t) = [y(t)]
dt
Examples:
➢ y(n) = 2 x(n) + 3 x(n – 1) + 4 x(n – 2)
➢ y(n) = 2 x(n) + 3 x(n – 1) + 4 y(n – 2).
➢ y(n) = 2 x(n) + n x(n – 1) + 4 y(n – 2).
➢ y(n) = 2 x(n2) + 3 x(n + 1) + 4 y(n – 2).
➢ y(n) = 2 x(n/2) + 3 x(n – 1) + 4 y(n – 2).
h(n) = 2,3,4
1
Example-2: Determine the impulse response of a discrete time system y(n) − y (n − 1) = x( n) by
2
assuming zero initial conditions.
Put, x(n) = δ(n) and y(n) = h(n)
1
h(n) = h(n − 1) + (n)
2
1 1
n = −1 h( − 1 ) = h(−2) + (−1) = (0) + 0 = 0
2 2
1 1
n = 0 h( 0 ) = h(−1) + (0) = (0) + 1 = 1
2 2
1 1 1
n = 1 h(1 ) = h(0) + (1) = (1) + 0 =
2 2 2
11 2
1 1
n = 2 h( 2 ) = h(1) + (2) = + 0 =
2 2 2 2
1 1 1 2 1
3
n = 3 h( 3 ) = h(2) + (3) = + 0 =
2 22 2
…..
…..
n n
1 1
h(n) = , n 0 or h(n) = u (n)
2 2
Example-1: Determine the impulse response and response of the discrete time system
y(n) − ay(n − 1) = x(n) by assuming zero initial conditions with input x(n)=bnu(n).
Put, x(n) = δ(n) and y(n) = h(n) in h(n) = ah(n − 1) + (n)
n = −1 h( − 1 ) = ah(−2) + (−1) = a(0) + 0 = 0
n = 0 h( 0 ) = ah(−1) + (0) = a (0) + 1 = 1
n = 1 h(1 ) = ah(0) + (1) = a(1) + 0 = a
Case-1: if b=a
n
y(n) = x(n) h(n) = a n (1) = a n( 1 + 1 + 1 + ........1,Add,n + 1,times) = (n + 1 )a n ,n 0, or, (n + 1 )a nu(n)
m
m=0
Case-2: if b≠a
n +1
b
1−
a n +1 − b n +1 a n +1 − b n +1
m
n
b a
y(n) = x(n) h(n) = a n = a n = , n 0, or, u ( n)
m=0 a 1−
a a−b a−b
b
Response of a discrete time system by applying input with zero initial conditions is called zero state
response or forced response.
() ()
n n
1 1
y (n) = A r n + B r n = A + B
h 1 2 2 3
Use initial conditions y(-1)=y(-2)=1 to evaluate constants A and B
n n
1 1 5 1
y (n) = y (n) = A + B = y (n − 1) − y (n − 2)
h 2 3 6 6
0 0
1 1 5 1
n = 0 y(0) = y (0) = A + B = y(−1) − y (−2)
h 2 3 6 6
5 1 4 2
y (0) = A + B = − = = 3 A + 3B = 2 − − − − − −(1)
6 6 6 3
1 1
1 1 5 1
n = 1 y(1) = y (1) = A + B = y(0) − y(−1)
h 2 3 6 6
A B 5 2 1 10 − 3 7 7
y (1) = + = − = = 3 A + 2 B = − − − − − (2)
2 3 63 6 18 18 3
7 1 2 2 1
Solve equations (1) and (2) 3B − 2 B = 2 − B = − and A = − B = + = 1
3 3 3 3 3
n n
1 11
Natural Response, yn (n) = − , n −2
2 3 3
(b)Forced Response: Response of the system by applying input with zero initial conditions
n
1
Particular solution of given system, y (n) = k u (n).
p 4
n n
1 1
Find the constant ‘k’ by substituting x(n) = u (n) & y(n) = y (n) = k u (n) in given system
4 p 4
10 8 3 − 10 + 8 1
k 1 − + = 5 k = 5 k = 5 k = 15
3 3 3 3
n nn
1 1 1
Forced response of given system, y (n) = y (n) + y (n) = A + B + 15 , n 0
f h p 2 3 4
Find constants A & B with zero initial conditions
n nn 5 n n −1
1 1 1 1 1 1
y(n) = A + B + 15 = y(n − 1) − y(n − 2) + u (n) + u (n − 1)
2 3 4 6 6 4 4
0 0 5
0 0 −1
1 1 1 1 1 1
n = 0 y(0) = A + B + 15 = y(−1) − y(−2) + u (0) + u (−1)
2 3 4 6 6 4 4
5 1
y (0) = A + B + 15 = (0) − (0) + 1 + 0 = 1 A + B = −14 − − − − − − − (1)
6 6
1 1 1 1 0
1 1 1 5 1 1 1
n = 1 y(1) = A + B + 15 = y(0) − y(−1) + u (1) + u (0)
2 3 4 6 6 4 4
A B 15 5 1 1 A B 15 5 1
y (1) = + + = (1) − (0) + + 1 + = − + + + 1
2 3 4 6 6 4 2 3 4 6 4
15 6 5 6 1 6 45 3 3 − 45
3A + 2B = − + + + 6 = − + 5 + + 6 = 11 + = 11 − 21 = −10 − − − (2)
4 6 4 2 2 2
Solve equations (1) and (2)3A+2(-14-A)=-103A-2A=28-10A=18 & B=-14-A=-14-18=-32
n n n
1 1 1
Forced Response, y (n) = 18 − 32 + 15 , n 0
f 2 3 4
(c)Response of the System: Sum of natural and forced response
n n n n n
1 11 1 1 1
y(n) = yn (n) + y (n) = − + 18 − 32 + 15 , n 0
f 2 3 3 2 3 4
n n n
1 1 1 1
y(n) = 19 − + 32 + 15 , n 0
2 3 3 4
n n n
1 97 1 1
y(n) = 19 − + 15 , n 0
2 3 3 4
() ()
n n
1 1
y (n) = A r n + B r n = A + B
h 1 2 2 3
Use initial conditions y(-1)=y(-2)=1 to evaluate constants A and B
n n
1 1 5 1
y (n) = y (n) = A + B = y (n − 1) − y (n − 2)
h 2 3 6 6
0 0
1 1 5 1
n = 0 y(0) = y (0) = A + B = y(−1) − y (−2)
h 2 3 6 6
5 1 4 2
y (0) = A + B = − = = 3 A + 3B = 2 − − − − − −(1)
6 6 6 3
1 1
1 1 5 1
n = 1 y(1) = y (1) = A + B = y(0) − y(−1)
h 2 3 6 6
A B 5 2 1 10 − 3 7 7
y (1) = + = − = = 3 A + 2 B = − − − − − (2)
2 3 63 6 18 18 3
7 1 2 2 1
Solve equations (1) and (2) 3B − 2 B = 2 − B = − and A = − B = + = 1
3 3 3 3 3
n n
1 11
Natural Response, yn (n) = − , n −2
2 3 3
(b)Forced Response: Response of the system by applying input with zero initial conditions
n
1
Particular solution of given system, y (n) = kn u (n).
p 2
n n
1 1
Find the constant ‘k’ by substituting x(n) = u (n) & y(n) = y (n) = kn u (n) in given
2 p 2
system
rn − rn −1 + rn − 2 = 0
2 1
3 9
2 1
r2 − r + = 0
3 9
2
1
r − = 0
3
1
r = r1 = r2 =
3
Homogeneous solution of given system
n
1
y (n) = ( A + nB)(r )n = ( A + nB)
h 3
Use initial conditions y(-1)=y(-2)=1 to evaluate constants A and B
n
1 2
y (n) = y (n) = ( A + nB) = y (n − 1) − y (n − 2)
1
h 3 3 9
0
1 2 1
n = 0 y (0) = y (0) = ( A + (0) B) = y (−1) − y (−2)
h 3 3 9
2 1 6 −1 5
y (0) = A = − = = − − − − − −(1)
3 9 9 9
1
1 2 1
n = 1 y (1) = y (1) = ( A + B) = y (0) − y (−1)
h 3 3 9
A + B 2 5 1 10 − 3 7 7
y(1) = = − = = A + B = − − − − − ( 2)
3 39 9 27 27 9
7 5 2
Solve equations (1) and (2) B = − =
9 9 9
n n n+2
5 2 1 1 1
Natural Response, yn (n) = + n = (5 + 2n ) = (5 + 2n )
1
, n −2
9 9 3 9 3 3
r n − 2r n − 1 + 4r n − 2 = 0 r 2 − 2r + 4 = 0
2 4 − 16 2 − 12 2 j 2 3
r = = = =1 j 3
2 2 2
r1 = 1 + j 3 & r2 = 1 − j 3 Where
r = 1+ 3 = 4 = 2
Homogeneous solution of given system
= Tan−1 ( 3 ) =
y (n) = ( ACosn + BSinn )(r )n 3
h
n n n
y (n) = ACos + BSin (2)
h 3 3
Use initial conditions y(-1)=y(-2)=1 to evaluate constants A and B
n n n
y(n) = y (n) = ACos + BSin (2) = 2 y(n − 1) − 4 y(n − 2)
h 3 3
n = 1 y(1) = ACos + BSin (2)1 = 2 y(0) − 4 y(−1)
3 3
1 3
A + B 2 = 2(−2) − 4
2 2
( )
− 2 + B 3 = −4 − 4 = −8
B 3 = 2 − 8 = −6
6
B=− = −2 3
3
Natural Response,
n n
y (n) = − 2Cos − 2 3Sin (2)n ; n −2
n 3 3
n n
y (n) = −2 Cos + 3Sin (2)n ; n −2
n 3 3
n n
y (n) = − Cos + 3Sin (2)n + 1; n −2
n 3 3
Where T[•] is transform operator and the relation between input x(n) and output y(n) of a discrete
time system is represented with y(n) = T[ x(n)], it shows the output y(n) is the transformation of
input x(n).
Various properties or classification of discrete time systems are given below
➢ Linear and Nonlinear Systems
➢ Shift Invariant and Variant Systems
➢ Static and Dynamic Systems
➢ Causal and Noncausal Systems
➢ Stable and Unstable Systems
Apply x1(n) as input to the system T[•] and observe the output, take it as y1(n)
Apply x2(n) as input to the system T[•] and observe the output, take it as y2(n)
Example-1: Test the discrete time system y(n) = T[ x(n)] = 2x(n)+3 for linearity.
Given discrete time system y(n) = T[ x(n)] = 2x(n)+3
Apply x1(n) as input to the system y(n) = T[ x(n)] = 2x(n)+3 and observe the output, take it as y1(n)
y1(n) = T[ x1(n)] = 2x1(n)+3
Apply x2(n) as input to the system y(n) = T[ x(n)] = 2x(n)+3 and observe the output, take it as y2(n)
y2(n) = T[ x2(n)] = 2x2(n)+3
ay1(n)+by2(n)=a(2x1(n)+3)+b(2x2(n)+3)= 2ax1(n)+3a+2bx2(n)+3b=2ax1(n)+2bx2(n)+3a+3b---(1)
Apply ax1(n)+bx2(n) as input to the system y(n) = T[ x(n)] = 2x(n)+3 and observe the output
T[ax1(n)+bx2(n)]=2(ax1(n)+bx2(n))+3=2ax1(n)+2bx2(n)+3---(2)
Compare equations (1) and (2) ay1(n)+by2(n)≠T[ax1(n)+bx2(n)] Given system is nonlinear.
Example-2: Test the discrete time system y(n) = T[ x(n)] = 2x(n)+3x(n-1) for linearity.
Given discrete time system y(n) = T[ x(n)] = 2x(n)+3x(n-1)
Apply x1(n) as input to the system y(n) = T[ x(n)] = 2x(n)+ 3x(n-1) and observe the output, take it
as y1(n) y1(n) = T[ x1(n)] = 2x1(n)+3x1(n-1)
Apply x2(n) as input to the system y(n) = T[ x(n)] = 2x(n)+3x(n-1) and observe the output, take it as
y2(n) y2(n) = T[ x2(n)] = 2x2(n)+3x2(n-1)
ay1(n)+by2(n)=a(2x1(n)+3x1(n-1))+b(2x2(n)+ 3x2(n-1))= 2ax1(n)+3ax1(n-1)+2bx2(n)+3bx2(n-1)
ay1(n)+by2(n) =2ax1(n)+2bx2(n)+3ax1(n-1)+3bx2(n-1)---(1)
Apply ax1(n)+bx2(n) as input to the system y(n) = T[ x(n)] = 2x(n)+3x(n-1) and observe the output
T[ax1(n)+bx2(n)]=2(ax1(n)+bx2(n))+3(ax1(n-1)+bx2(n-1))
T[ax1(n)+bx2(n)]=2ax1(n)+2bx2(n)+3ax1(n-1)+3b x2(n-1)---(2)
Compare equations (1) and (2) ay1(n)+by2(n)=T[ax1(n)+bx2(n)] Given system is linear.
Apply delayed form x(n), i.e x(n-n0) as input to the system T[•] and observe the output, if the output
is y(n-n0), then the system is shift invariant otherwise the system is shift variant.
Condition for a shift invariant system is y(n-n0)=T[x(n-n0)] and it is also known as time invariant.
Example-1: Test the discrete time system y(n) = T[ x(n)] = 2x(n)+3 for time invariance.
Example-2: Test the discrete time system y(n) = T[ x(n)] = nx(n)+3 for time invariance.
Static systems are those for which the present output y(n) depends on only present input x(n). Static
systems are also known as memory less systems.
Examples:
➢ y(n) = T[ x(n)] = nx(n)
➢ y(n) = T[ x(n)] = 2x(n)+3x2(n)
➢ y(n) = T[ x(n)] = n2x(n)
n
➢ y(n) = T[x(n)] = +2x(n)
x(n)
➢ y(n) = T[ x(n)] = 2x(n)
Dynamic systems are those for which the present output y(n) depends on past inputs x(n-1), x(n-2),.
and/or past outputs y(n-1), y(n-2),… and/or future inputs x(n+1), x(n+2),…. and/or future outputs
y(n+1), y(n+2),….. Dynamic systems are also known as memory systems.
Examples:
➢ y(n) = T[ x(n)] = 2x(n)+3x(n-1)
➢ y(n) = T[ x(n)] = 2x(n)+3x2(n-2)
➢ y(n) = T[ x(n)] = n2x(n+1)
n
➢ y(n) = T[x(n)] = +2y(n-3)
x(n)
➢ y(n) = T[ x(n)] = 2x(n-5)
N M
➢ y(n) = T[x(n)] = - aK y(n − k) + bK x(n − k)
k =1 k =0
➢ y(n) = T[x(n)] = b
k =0
K x(n − k)
Causal systems are those for which the present output y(n) depends on only present input x(n)
and/or past inputs x(n-1), x(n-2),…. and/or past outputs y(n-1), y(n-2),… but does not depends on
future inputs x(n+1), x(n+2),…. and/or future outputs y(n+1), y(n+2),….. All static systems are
causal.
Non-causal systems are those for which the present output y(n) depends on future inputs x(n+1),
x(n+2),…. and/or future outputs y(n+1), y(n+2).
Examples:
➢ y(n) = T[ x(n)] = 2x(n)+3x2(n)
➢ y(n) = T[ x(n)] = n2x(n)+3 Causal and Static (memory less)
➢ y(n) = T[x(n)] = b
k = −4
K x(n − k)
➢ h(n)=2n u(-n)
Noncausal Systems
➢ h(n)=(1/2)n u(-n)
Apply bounded signal x(n) as input to the system T[•] and observe the output y(n), if the output is
bounded, then the system is stable otherwise the system is unstable. Stable systems are called
Bounded Input Bounded Output (BIBO) systems.
i.e., if x(n) < , then y(n) < , for all values of n.
Examples:
➢ y(n) = T[ x(n)] = 2x(n)+3
➢ y(n) = T[ x(n)] = 2x(-n-3)+3x2(n)
➢ y(n) = T[ x(n)] = 3Cos[x(n)]
N M Stable Systems
➢ y(n) = T[x(n)] = - k =1
aK y(n + k) + bK x(n − k)
k =0
➢ y(n) = T[x(n)] = b
k = −4
K x(n − k)
Example-1: Test the discrete time system h(n)=(1/2)n u(n) for stability.
n n 2
1 1 1 1 1
n = −
h(n) =
n = −
u (n) = = 1 + + + ....... =
2 n=0 2 2 2 1− 1
= 2 = Finite
2
Given causal system is stable, because the impulse response is absolutely summable.
Example-2: Test the discrete time system h(n)=2n u(n) for stability.
n = −
h(n) = 2n u(n) = 2n = 1 + 2 + 4 + 8 + ..... = Infinity
n = − n=0
Given causal system is unstable, because the impulse response is not absolutely summable.
n = −
n
2 2 1− 1
= 2 = Finite
2
Given noncausal system is stable, because the impulse response is absolutely summable.
Example-4: Test the discrete time system h(n)=(1/2)n u(-n) for stability.
0
h(n)
n = −
= (1/ 2) u(−n) = (1/ 2)
n = −
n
n = −
n
= 1 + 2 + 4 + 8 + ..... = Infinity
Given noncausal system is unstable, because the impulse response is not absolutely summable.
Fourier Transform of a discrete time signal is called Discrete Time Fourier Transform. The DTFT
of a signal x(n) is represented with X(ej) and it can be computed from the formula.
DTFT[ x(n) ] = X(e jω ) = x(n) e
n= −
- jω n
In general, the frequency domain X(ej) is in complex form and it can be expressed as
X(e jω ) = XR (e jω ) + j XI (e jω )
Where,
XR (e jω ) : Real part of X (e jω )
XI (e jω ) : Imaginary part of X (e jω )
Magnitude of X(ej) is called the magnitude spectrum and it can be computed from the formula
X(e jω ) = [ XR (e jω ) ] 2 + [ XI (e jω ) ] 2
Phase of X(ej) is called the phase spectrum and it can be computed from the formula
XI (e jω )
X (e ) = Tan
jω −1
jω
R
X (e )
Discrete time signal x(n) can be computed from the frequency domain X (e jω ) is called Inverse
Discrete Time Fourier Transform (IDTFT).
π
1
x(n) = IDTFT[X(e j )] =
2π − π
X(ej ) e jω n dω
Discrete time signal x(n) and the frequency domain X (e jω ) are Discrete Time Fourier
Transformable pairs.
⎯DT
⎯⎯FT
→
x(n) X(ej )
⎯ ⎯⎯
IDT FT
( )
X(ejω ) = an u(n)e- jω n = ae - jω
n
( ) 2
= 1 + ae - jω + ae - jω + .......... . =
1
1 − ae - jω
n = − n=0
1 1
X(e jω ) = =
1 − a(Cos − jSin ) 1 − aCos + jaSin
1 1 1
X(e jω ) = = =
1 − aCos + jaSin (1 − aCos ) + (aSin )
2 2
1 + a Cos − 2aCos + a 2 Sin2
2 2
1
X(e jω ) =
1 + a − 2aCos
2
X (e jω )
Phase spectrum X (e jω ) = Tan −1 I jω
R
X (e )
0 aSin −1 aSin
X(e jω ) = (1 + j 0) − (1 − aCos + jaSin ) = Tan−1 − Tan−1 = −Tan
1 1 − aCos 1 − aCos
Example-2: Determine the DTFT of a sequence x(n)=(1/2)n u(n), hence obtain the magnitude and
phase spectrum.
DTFT[ x(n) ] = X(e jω ) = x(n) e
n= −
- jω n
n n 2
1 1 - jω 1 - jω 1 - jω
X(e ) =
jω
n = −
u(n)e
2
- jω n
=
n=0
e = 1 + e + e + .......... .
2 2 2
1 2 2 2
X(e jω ) = = = =
1 − e-jω 2 − e
1 -jω
2 − (Cos − jSin ) 2 − Cos + jSin
2
2 2 2 2
X(e jω ) = = = =
2 − Cos + jSin (2 − Cos ) + (Sin )
2 2
4 + Cos − 4Cos + Sin
2 2
5 − 4Cos
X (e jω )
Phase spectrum X (e jω ) = Tan −1 I jω
X R (e )
0 Sin −1 Sin
X(e jω ) = (2 + j 0) − (2 − Cos + jSin ) = Tan−1 − Tan−1 = −Tan
2 2 − Cos 2 − Cos
-n 2
0 0
1 jω 1 jω 1 jω
X(e ) =
jω
a
n = −
n
u(-n)e - jω n
= ae
n = −
n - jω n
=
n = −
e = 1 + e + e + .......... .
a a a
1 a a a
X(e jω ) = = = =
1 − e jω a − e
1 jω
a − (Cos + jSin ) a − Cos − jSin
a
a a a
X(e jω ) = = =
a − Cos − jSin (a − Cos ) + (− Sin )
2 2
a + Cos − 2aCos + Sin2
2 2
a
X(e jω ) =
1 + a − 2aCos
2
X (e jω )
Phase spectrum X (e jω ) = Tan −1 I jω
X R (e )
0 − Sin −1 Sin
X(e jω ) = (a + j 0) − (a − Cos − jSin ) = Tan−1 − Tan−1 = Tan
a a − Cos a − Cos
Example-4: Determine the DTFT of a sequence x(n)=2n u(-n), hence obtain the magnitude and
phase spectrum.
DTFT[ x(n) ] = X(e jω ) = x(n) e
n= −
- jω n
-n 2
0 0
1 jω 1 jω 1 jω
X(e ) =
jω
2
n = −
n
u(-n)e - jω n
= 2e
n = −
n - jω n
=
n = −
e = 1 + e + e + .......... .
2 2 2
1 2 2 2
X(e jω ) = = = =
1 − e jω 2 − e
1 jω
2 − (Cos + jSin ) 2 − Cos − jSin
2
2 2 2 2
X(e jω ) = = = =
2 − Cos − jSin (2 − Cos ) + (− Sin )
2 2
4 + Cos − 4Cos + Sin
2 2
5 − 4Cos
X (e ) jω
Phase spectrum X (e jω ) = Tan −1 I jω
R
X (e )
0 − Sin −1 Sin
X(e jω ) = (2 + j 0) − (2 − Cos − jSin ) = Tan−1 − Tan−1 = Tan
2 2 − Cos 2 − Cos
(ae ) + (ae-jω )
−1 −
a − ne-jωn + a ne-jωn = jω − n n
X(e jω ) = a e-jωn =
n
n = − n = − n =0 n = −1 n =0
X (e ) = ae + (ae
jω jω jω 2
) + .......... . + 1 + ae -jω
+ (ae
-jω 2
) + .......... .
2
X (e ) = −1 + 1 + ae + ae
jω jω
( jω 2
) + .......... . + 1 + ae -jω 1
+ e-jω + .......... .
2
1 1 (1 − ae-jω + 1 − ae jω ) 2 − a(e jω + e− jω )
X(e ) = −1 +
jω
+ = −1 + = −1 +
1 − ae jω 1 − ae-jω (1 − ae jω )(1 − ae-jω ) 1 + a 2 − a(e jω + e − jω )
2 − 2aCos − (1 + a 2 − 2aCos ) + 2 − 2aCos 1− a2
X(e jω ) = −1 + = =
1 + a − 2aCos
2
1 + a − 2aCos
2
1 + a − 2aCos
2
n
1
Example-6: Determine the DTFT of a sequence x(n) = , hence obtain the magnitude and
2
phase spectrum.
DTFT[ x(n) ] = X(ejω ) = x(n) e
n = −
- jω n
n −1 −n n −n n
1 -jωj 1 -jωj 1 -jωj − 1 jω
1 -jω
X(e jω ) =
n = −
e = e + e = e + e
2 n = − 2 n =0 2 n = −1 2 n =0 2
2 2
1 1 1 1
X (e jω ) = e jω + e jω + .......... . + 1 + e-jω + e-jω + .......... .
2 2 2 2
2 2
1 1 1 1
X (e jω ) = −1 + 1 + e jω + e jω + .......... . + 1 + e-jω + e-jω + .......... .
2 2 2 2
1 1 2 2 2(2 − e-jω + 2 − e jω )
X(e jω ) = −1 + + = −1 + + = −1 +
1 1
1 − e jω 1 − e-jω 2 − e jω 2 − e-jω (2 − e jω )(2 − e-jω )
2 2
2( 4 − 2Cos ) 8 − 4Cos − 5 + 4Cos + 8 − 4Cos 3
X(e jω ) = −1 + = −1 + = =
(4 − 4Cos + 1) 5 − 4Cos 5 − 4Cos 5 − 4Cos
DTFT[ (n) ] = X(e ) = jω
(n) e
n = −
- jω n
= ..... + (0)e- j0 + ..... = 1 1 = 1
1 2 jSin(n / N ) Sin(n / N )
x(n) = =
2π jn n
j e jw + e − jw e j 2 w + e− j 2 w e j 3w + e− j 3w
X(e ) = 1 + +2 +3
2 2 2
1 1 3 3
X(ej ) = 1 + e jw + e − jw + e j 2 w + e − j 2 w + e j 3 w + e − j 3 w
2 2 2 2
3 j 3w 1 1 3
X(ej ) = e + e j 2 w + e jw + 1 + e − jw + e − j 2 w + e − j 3w
2 2 2 2
x(n) e
n = −
-jwn
= x(−3)e j 3w + x(−2)e j 2 w + x(−1)e jw + x(0) + x(1)e− jw + x(2)e− j 2 w + x(3)e− j 3w
3 1 1 3
x(n) = ,1, , 1 , ,1,
2 2 2 2
Example-10: Determine the DTFT of a periodic sequence x(n) = 1 .
x(n)
n
……. -7 -6 -5 -4 -3-2 -1 0 1 2 3 4 5 6 7…………
DTFT[ x(n) ] = X(e jω ) = 2π (ω − 2πm)
m = −
(A)Linear Property:
If x1(n), x2(n) are two discrete time signals and DTFT[x1(n) ] = X1(ej), DTFT[x2(n) ] = X2(ej),
then DTFT[a x1(n) + b x2(n) ] = a X1(ej) + b X2(ej) is called linear property of DTFT.
Proof:
From the basic definition of DTFT
DTFT[ x(n) ] = x(n) e
n = −
- jω n
= a X1 (e ) + b X 2 (e )
jω jω
DTFT[x(n)] = DTFT [a.b nu (n) + b.a nu (n)] = aDTFT [b nu (n)] + bDTFT [a nu (n)]
X(e jω ) = x(n) e
n = −
- jω n
Replace ω by ω + 2π
X(e j(ω+ 2 ) ) = x(n) e
n = −
- j(ω + 2 ) n
= x(n) e
n = −
- jω n
e - j 2 n
= x(n) e
n = −
- jω n
(1)
= x(n) e
n = −
- jω n
= X(e j ω )
Example: Determine the DTFT of a real and even sequence x(n) = 3,2, 1 ,2,3 , hence obtain the
period of X(e j ).
2
DTFT[ x(n) ] = x(n) e
n = −
- jω n
= x(n) e
n = −2
- jω n
= x(-2) e j2ω x(-1) e jω + x(0) e- j0 + x(1) e- jω + x(2) e- j2ω
X(e j ( ω + 2 ) ) = X(e jω )
=e - jω no
X(e jω )
DTFT[x(n)] = DTFT [a.b nu (n) + b.a nu (n)] = aDTFT [b nu (n)] + bDTFT [a nu (n)]
a + b − (a 2 + b 2 )e-jω
Y(e jω ) = e-jω
1 − (a + b)e-jω + abe-j 2ω
Y(e jω
)=
(a + b)e-jω − (a 2 + b2 )e-j 2ω
1 − (a + b)e-jω + abe-j 2ω
jωon
Replace x(n) by e x(n)
DTFT[ e jω o n x(n) ] = e
n = −
jω o n
x(n) e - jω n
= x(n) e
n = −
jω o n
e - jω n
= x(n) e
n = −
- j(ω - ω o) n
= X(e j(ω - ω o) )
Let – n = m
= x(m) e
m = −
- jω (-n)
Then n = - m
−
= x(m) e jω n Range of summation
m=
LL (m) = –(– ) =
UL (m) = –( ) = –
= x(m) e
m = −
- j (-ω ) n
= X (e - jω )
= x * (n) [e
n = −
jω n
]*
= [ x(n)
n = −
e jω n ]*
*
= x(n) e jω n
n = −
*
= x(n) e- j(-ω ) n
n = −
*
= X( e- jω )
n = −
= [X(e )] *
- jω
= X * (e- jω )
Note : Apply time reversal property on DTFT[ x*(n) ] = X*(e – j )DTFT[ x*(–n) ] = X*(e j
)
Example-1: Determine the DTFT of a real sequence x(n) = 1, 2,3,4
2
DTFT[ x(n) ] = x(n) e
n = −
- jω n
= x(n) e
n = −1
- jω n
= x(-1) e jω + x(0) e- j0 + x(1) e- jω + x(2) e- j2ω
Example-2: Determine the DTFT of a real and even sequence x(n) = 3,2, 1 ,2,3
2
DTFT[ x(n) ] = x(n) e
n = −
- jω n
= x(n) e
n = −2
- jω n
= x(-2) e j2ω x(-1) e jω + x(0) e- j0 + x(1) e- jω + x(2) e- j2ω
Example-3: Determine the DTFT of a real and odd sequence x(n) = 3,2, 0 ,−2,−3
2
DTFT[ x(n) ] =
n = −
x(n) e- jω n = x(n) e
n = −2
- jω n
= x(-2) e j2ω x(-1) e jω + x(0) e- j0 + x(1) e- jω + x(2) e- j2ω
If x(n) is a discrete time signal and DTFT[ x(n) ] = X(ej), then DTFT[nx(n)] = j
d
dω
X(e j ω ) is
DTFT[ x(n) ] = x(n) e
n = −
- jω n
X(ejω ) = x(n) e
n = −
- jω n
d d - jω n
X(ejω ) = x(n) e
d n = − d
= x(n) (-jn)e
n = −
- jω n
= − j nx(n) e- jω n
n = −
= − jDTFT [nx(n)]
DTFT [nx(n)] = j
d
d
X(ej ω )
DTFT[ x(n) ] = x(n) e
n = −
-jω n
Apply convolution formula x1(n)*x2(n) = x (m)x (n-m)
m = -
1 2
= x (m)x (n-m) e
n = − m = -
1 2
-jω n
= x1(m) e -jω m X 2(e jω ) = X 2(e jω ) x1(m) e -jω m
m = − m = −
Proof:
2 *
n = −
x(n) = x(n)
n = −
x * (n) = x(n)
n = −
[ x(n) ]
*
1 π 1 π
= x(n) X(e jω ) e jω n
dω = x(n) X * (e jω ) e - jω n dω
n = − 2π − π n = − 2π − π
π π
jω
( )
1 1
=
2π −π
X * (e )
n = −
x(n) e -jωn
dω =
2π − π
X * (e jω ) X(e jω ) dω
2
1
=
2π −
X(ejω ) dω
Example:
π 2
(i) DTFT[ x(n) ] = X (e- jω ) = x(n) e
n = −
- jω n
X(e jω ) w−0 = x(n) = 1 + 4 + 3 + 2 + 5 + 6 + 7 + 0 + 9 = 37
n = −
π
1
(ii) x(n) = IDTFT[X(e )] = X(e
jω jω
) e jωn dω
2 −π
π π
1
Put n = 0 x( 0 ) = X(e jω )dω X(e )dω = 2x(0) = 2 (3) = 6
jω
2 −π −π
π 2 2
= 2 (12 + 42 + 32 + 22 + 52 + 62 + 7 2 + 02 + 92 ) = 442
X(e ) jω
dω = 2
n = −
x(n)
−π
The ratio between output y(n) to input x(n) in frequency domain representation is called Transfer
function or System function or Frequency response of discrete LSI System and it is represented
with H(ej).
DTFT [ y (n)] Y(e jω )
H(e jω ) = =
DTFT [ x(n)] X(e jω )
In general, the frequency response H(ej) is in complex form and it can be expressed as
H(e jω ) = HR (e jω ) + j HI (e jω )
Where,
HR (e jω ) : Real part of H (e jω )
HI (e jω ) : Imaginary part of H (e jω )
Magnitude of H(ej) is called the magnitude response and it can be computed from the formula
Phase of H(ej) is called the phase response and it can be computed from the formula
H (e jω )
H (e jω ) = Tan −1 I jω
HR (e )
Impulse or Unit Sample Response:
Response of a discrete LSI system with an input of impulse or unit sample sequence is called is
called impulse response or unit sample response and it is represented with h(n).
Unit sample response h(n) can be obtained from the frequency response H(ej) by using Inverse
Discrete Time Fourier Transform (IDTFT).
Response of discrete LSI system y(n) can be obtained from the frequency domain by using Inverse
Discrete Time Fourier Transform (IDTFT).
Example-1: Evaluate (a) Frequency Response (b) Magnitude Response and (c) Phase Response
1 1
(d)Unit Sample Response of a discrete LSI system y(n) − y (n − 2) = x(n) + x(n − 1) .
4 2
(a)Frequency Response
1 1
Given discrete LSI system y(n) − y (n − 2) = x(n) + x(n − 1)
4 2
Apply DTFT
1 1
DTFT y(n) − y(n − 2) = DTFT x(n) + x(n − 1)
4 2
1 1
Y ( e j ) − e − j 2 Y ( e j ) = X ( e j ) + e − j X ( e j )
4 2
1 1
Y (e j )1 − e− j 2 = X (e j )1 + e− j
4 2
1 − j 1
Y ( e j ) 1 + 2 e 1 + e − j 1
= = 2 =
j
X (e ) 1 − 1 e − j 2 1 − j 1 − j 1 − 1 e − j
1 + e 1 − e
4 2 2 2
1 2 2 2
H(e jω ) = = = =
1 − e-jω 2 − e
1 -jω
2 − (Cos − jSin ) 2 − Cos + jSin
2
(b)Magnitude Response
2 2 2 2
H(e jω ) = = = =
2 − Cos + jSin (2 − Cos ) + (Sin )
2 2
4 + Cos − 4Cos + Sin
2 2
5 − 4Cos
2
H(e jω ) =
5 − 4Cos
Sin
H(e jω ) = −Tan−1
2 − Cos
(d)Unit Sample Response
1 1 n
h(n) = IDTFT [ H(e )] = IDTFT
jω
= u ( n)
1 − 1 e-jω 2
2
1 1
Example-2: Evaluate the response of a discrete LSI system y(n) − y (n − 2) = x(n) + x(n − 1) with
4 2
n
1
an input of x(n) = u (n)
3
1 1
Given discrete LSI system y(n) − y (n − 2) = x(n) + x(n − 1)
4 2
Apply DTFT
1 1
DTFT y(n) − y(n − 2) = DTFT x(n) + x(n − 1)
4 2
1 1
Y ( e j ) − e − j 2 Y ( e j ) = X ( e j ) + e − j X ( e j )
4 2
1 1
Y (e j )1 − e− j 2 = X (e j )1 + e− j
4 2
1 1
j 1 + e − j 1 + e − j
Y (e ) 1
j
= H ( e j ) = 2 = 2 =
1 − e − j 2 1 + e − j 1 − e − j 1 − e − j
X (e ) 1 1 1 1
4 2 2 2
1 n 1
Y (e ) = X (e ) H (e ), X (e ) = DTFT [ x(n)] = DTFT u (n) =
j j j j
3 1
1 − e− j
3
1 1 −2 3
Y (e ) =
j
= +
1 − 1 e − j 1 − 1 e − j 1 − 1 e − j 1 − 1 e − j
3 2 3 2
n n
ramp sequence.
3. Determine the convoluted signal y(n)=x(n)*h(n), given that x(n)=anu(n) and h(n)=bnu(n),
if (i) a≠b (ii) a=b.
4. Determine the convoluted sequence y(n) = x1 (n) * x2 (n) , where x1 (n) = {1,2, 3,4} and
6. Analyze the following signals for causal, non-causal, bounded, un-bounded (i) x1(n)=(1/2)n
(ii) x2(n)=(2)n (iii) x3(n)=(1/2)nu(n) (iv) x4(n)=(2)nu(n) (v) x5(n)=(1/2)nu(-n-1)
(v) x6(n)=(2)nu(-n-1).
7. Analyze the following discrete time signals for periodic and aperiodic. If periodic, then
compute the period.(i)x(n)=4Cos(nπ/4)+6Sin(nπ/6) (ii) y(n)=4Cos(n/4)+6Sin(n/6).
8. Analyze the following discrete time signals for even and odd. If fails to satisfy even and odd
properties, then evaluate even and odd parts.(i) x1 (n) = {1,2, 3,1,2} (ii) x2 (n) = {−1,−2, 3,1,2} .
10. Solve the difference equation y(n)-ay(n-1)=x(n) to obtain the impulse response, hence obtain
the system response for an input x(n)=bnu(n), where 0<a<1 and 0<b<1.
12. Determine (a) DTFT of x(n) (b) Magnitude spectrum and (c)Phase spectrum.
13. Determine the DTFT of x(n), Magnitude spectrum and Phase spectrum for the signals
n
1
(i)x(n)=(1/2)nu(n) (ii)x(n)=(2)nu(-n-1) (iii) x(n) = .
2
π
jω
14. Apply properties of DTFT and Evaluate (i) X(ej) at = 0 (ii)
−π
X(e ) dω (iii)
2
15. Apply properties and Evaluate (a)Frequency response (b)Magnitude response and (c)Phase
response of a system having LCCDE y(n)-0.5y(n-1)=x(n)+0.5x(n-1).
16. Determine (a)Unit sample response (b)Unit step response and (c)Response of a causal
system for an input x(n)=(1/2)nu(n).Given LCCDE y(n)-0.25y(n-1)=x(n).
17. Find the response of the system for an input x(n)=(1/4)nu(n) having LCCDE y(n)-5/6 y(n-
1)+1/6 y(n-2)=x(n)+x(n-1). Given initial conditions y(-1)=y(-2)=1.
18. For the system y(n)-5/6 y(n-1)+1/6 y(n-2)=x(n)+x(n-1), find(a)Natural response for given
initial conditions y(-1)=y(-2)=1 (b)Forced response for given input x(n)=(1/2)nu(n) and
(c)Response of the system for an input x(n)=(1/2)nu(n) and initial conditions y(-1)=y(-2)=1.
19. Find the natural response of the system having LCCDE y(n)-2 y(n-1)+4 y(n-2)=x(n)+x(n-1)
for given initial conditions y(-1)=y(-2)=1.
20. Analyze the following systems for causality and stability(a)y(n)-0.5y(n-1)=x(n) with zero
initial conditions(b)y(n)-0.5y(n-1)=x(n) with initial condition y(-1)=1(c)y(n)-2y(n-1)=x(n)
with zero initial conditions(d)y(n)-2y(n-1)=x(n) with initial condition y(-1)=1
(1)δ(n)=u(n)-u(n-1) (2)δ(n)=u(n+1)-u(n-1)
1,3,4
n
(3) u (n) = (k ) (4) u (n) = (n − k )
k = − k =0
a.Time shifting b.Time scaling c.Folding d.Amplitude scaling e.Time shifting 3-c,d
4-a,b
4-a,d