Introduction To Discrete Time Signals & System
Introduction To Discrete Time Signals & System
Introduction To Discrete Time Signals & System
apter 01 Questions
Dec. 2010
Give any five classifications of Discrete time systems with examples
X (t) = sin(480 t) + 3 sin (720 t) is sampled with Fs = 600 times per sec.
(1) What are the frequencies in radians in the resulting DT signal x (n)?
June 2011
If x (n) = { 2, -1, 3, 0, 4 } obtain following:
(i) x (-n)
(ii) x (n-1) (iii) x (n+1) (iv) x (-n+2)(v) x (2n)
Classify following DT System on linearity/ causality and time variance:(i) y (n) = 2x(n) + x(n-1)
(ii) y (n) = x(2n) +2
1. x (n)= [1, 2, 1, 2, 3, 4]
2. Explain whether the following signals are power signal or energy signal.
a.
0.5 n u (n)
b.
A cos (n)
j( /4) n
x (n) = e
b.
y [n] = [8 4 8 4 8 4 8 4]
Signal
A signal is defined as any physical quantity that varies with:
1. Time,
2. Space or
3. Any other independent variable or variables
- Examples of Signals
- Speech
- Music
- Pictures
- Video
- ECG
- Mathematically, we describe a signal as a
function of one or more independent variables
s1(t ) = 5 t
s2( t) = 20 t2
f (x, y) = 3x + 2x y + 10 y2
two independent variables x, and y
t, Time
- An image signal is function of two variables:
x, y (Coordinates of an image)
- These signals are generated by some means:
- Speech signal by vocal cord and air flow.
- Image signals by exposing light- sensitive sensors to light
- The signals are generated by a system
Voice amplifier:
- Amplifier is made up of resistors, capacitors, transistors etc.
- Amplifier takes the analog signal (continuous) from microphone,
amplifies it , and produces analog output signal for speaker
Analog processing system:
- Takes analog input, process it in analog form and
produces analog output
x (t)
Analog
Input signal
t is time
x (n)
Digital
Input signal
n is number
(sample)
y (n)
Digital
Output signal
y (t)
Analog
Output signal
Storage of digital signals is very easy on magnetic media, memory, CD, DVD
pen drive.
Limitations of DSP:
1. The use of ADC and DAC may make the processing
comparatively, slower (Analog systems are faster)
2. The power consumption of digital signal processors may be higher
3. Not suitable for signals with huge bandwidth
(ADC, Sampling rate has to be two times the bandwidth)
Classification of signals
- Signal Processing technique for any signal depends upon the
characteristics of the signal.
Multi-channel and Multidimensional Signals:
MC-Multiple components
of same signal
MD- Multiple inputs
Real valued signal
s1(t) = A sin 3 t
Complex valued signal
s2 (t ) = A e j3 t = A cos3 t + j A sin 3t
Vector representation
The signal generated by:
multiple sources or
multiple sensors
z
Three component
signal
y
0
22
10
12
14
16
18
20
Classification of signals
x2(t) = e - | t |
- < tx<(t)
where
t
0 1 2 3 4 5 6 7 8 9 10
( = 2 f)
x (tn) = e
( instead of x (t) )
(T = Time period)
- There are some discrete time signals which are inherently discrete and
do not require the sampling of continuous signal
i.e.
Accumulating a variable over a period of time
(i.e. number of cars/ hour)
x (n) = 0.8n
x (n) = 0 for n <0
Digital signal:
A discrete-time signal having a set of discrete values.
For digital processing of a signal:
- It must be discrete in time
- Its values also must be discrete
- Digital signal is obtained by:
First - Obtaining a discrete time signal by
sampling an analog signal at discrete instants in time
Then Quantizing the values of discrete time to
a set of discrete values
Quantization:
It is the process of converting continuous-valued signal into
discrete valued signal by simple rounding / truncation process
or by mapping to a set of finite values
(1, 2, 3 and 4)
1, 1, 2, 1, 3, 2, 4, 2, 1}
X(n)
-2 -1 0 1 2 3 4 5 6 7
3. Sequential representation:
x (n) = { . . . 0, 0, 1, 4, 1, 0, 0, . . . }
(4-point sequence
Denoted as (n)
0, for n 0
Graphical representation
of unit impulse
-2 -1 0 1
2 3 4
Denoted as u (n)
1, for n 0
u (n)
0 , for n <0
Graphical representation
of step signal
-2 -1 0 1
2 3 4
denoted as u r (n ) u (n)
n, for n 0
0 , for n <0
-2 -1 0 1
2 3 4
for all n
a>1
x (n) = a n
x(n) = a n
x (n) = a n
If the parameter a is real x (n) is a real signal
- if the parameter a is complex valued,
a r e j
n
xx(n)
=
a
(n ) = r n: e j n
= r n ( cos n + j sin n)
Real part
xg (n) = r n cos n
r = 0.9
and = /10
Real part xg
Imaginary part xi
x (n ) = r n e j n
/2
0
2
-<
or 0 < 2
x (n) = x (-n)
x(-n) = - x(n) or
x (n) = - x (-n)
For an even signal x(1) = x (-1),
For an odd signal x(1) = - x(-1),
- If x (n) is odd then
x (2) = x (-2) . . .
x (2) = - x(-2) . . .
x(0 ) = 0
Examples:
A sine wave is odd,
while a cosine wave is even signal
sin (50) = - sin(-50)
cos (50) = cos ( -50)
x (n)
-4 -3 -2 -1 0 1 2 3 4
even signal
x (n) = x (-n)
x (n)
-4 -3 -2 -1
0 1 2 3 4
odd signal
x (n) = - x (-n)
x(n)
x (0 ) = 0
Any arbitrary signal can be expressed as the sum of the two components:
- One even and the
- Other odd
- The even signal component can be found out by:
x (-n) = {2, 2, 1, 4, 3
2,
1}
Check if
x (n)
f = 0.01/ 2
is defined as E |x (n) | 2
n=-
- As the magnitude square is used for x (n),
The definition is applies to real as well as complex-valued signals
Energy signal
- If Energy E is finite then x (n ) is called an energy signal
- The finite Energy can be called Ex of signal x (n)
N
N+1 n=-N
|x (n)| 2
E Lim
then
as N 1/(2N +1) 0
thus Average power P 0
Thus Power for a finite energy signal is zero
EN
1
0 1 2 3..
E |x (n) | 2 = |u (n) | 2 = 1 + 1 + 1 . . .
=
n=-
0
Since the E is infinite, the unit step signal is not an energy signal
Energy
1
N
The average Power P = lim
-----------
u 2 (n)
N 2N +1
n=0
( summation of u2 (n) is N +1)
1 (N+1)
1 + 1/ N
1
= lim
----------- = lim
----------- = ------N 2N +1
N
2 + 1/N
2
Consequently, the unit step sequence is power signal
and its energy is infinite
|an u (n)|2
EN
for n = 0 to
|an . 1|2
a| < 1
|x (n)|2
Thus E = 1/ 1- a2
if |a2| < 1
Explain whether the following signals are power signal or energy signal.
i 0.5 n u (n)
- First calculate energy
E of x (n) is given by:
of the signal
n=
n=
E = | x (n)| 2 =
| 0.5n u (n)| 2
n=-
n=-
-Since this signal is multiplied by unit step: n is from 0 -
( signal is zero for n > 0)
n=
n=
n=
E=
| (1/2) n | 2 =
(1/2) 2 n =
(1/4) n
n=0
n =0
n=0
n=
{Standard Geometric series
A n 1 + A + A2+ A 3 . . . 1/ (1-A) (A = )}
n=0
Thus E = 1 / (1 -1/4 ) = 1/ = 4/ 3 = 1.333
- Thus E is finite
Since Energy is finite it is energy signal
Dependent variable,
Amplitude
Delay
(k is an integer)
y (n) = x (n - k)
for k = 3, y (0) =x (-3)
X (n ) = { -1, 0, 1, 2, 3 , 4, 4, 4, 4, 4}
X(n-3) = { -1, 0, 1, 2, 3 , 4, 4, 4, 4, 4}
shifting delay
k=3
X( n +2) { -1, 0, 1, 2, 3 , 4, 4, 4, 4, 4}
advancing
k= -2 (+2 becomes origin)
2. Time reversal
Folding/ Reflection about the time origin n = 0:
- Replacing n by n
y (n) = x (-n)
x (n ) =
{2, 2, 2, 0, 1, 2, 3, 4}
x (-n) = {4, 3, 2, 1, 0, 2, 2, 2}
x (-n - 2) = {4, 3, 2, 1, 0, 2, 2, 2}
(Folding and delayed by 2)
3. Time scaling:
- Changing in time scale
- Two types down scaling and up scaling
- Down scaling is represented as:
y (n) = T [x (n)]
Example
y (n) = x (2n)
x (n)
= { -3, -2, -1, 0, 1, 2, 3, 4, 4, 4, 4, 4, 4, 4}
y (n) = x (2n)
= { -2,
0, 2, 4, 4, 4, 4}
x (n) = xa (nT)
Up scaling
y (n) = x (n / 2 )
x (n) = { 5, 4, 3, 2, 1, 2, 3, 4, 5}
y (0 )= x (0/2) = x(0) = 1
y (1)= x (1/2) = x(0.5) = no sample
y (2 )= x (2/2) = x(1) = 2
y (3 )= x (3/2) = x(1.5) = no sample
y (4 )= x (4/2) = x(2) = 3
y (6 )= x (6/2) = x(3) = 4
y (8 )= x (8/2) = x(4) = 45
Thus y (n) is expanded version of x (n) with y (1), y(3), y(5) . . No sample
-8 -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7 8
The sampling rate is increased from 1/T to 2/T
Upscaling
At -2, -1 . 0 = 0
-1, 2 . 0 = 0
0, 2 . 1 = 2
1, 1 . 2 = 2
2, 0 . 3 = 0
or
x (n )
y (n)
where T is transformation
(transformed to )
x (n)=
{ 0, 3, 2, 1, 0, 1 ,2, 3, 0}
x(n ) = { 0, 0, 3, 2, 1, 0, 1 , 2, 3, 0, 0, 0}
Y(0) = 1/3[ x (1) +x (0) + x(-1)] = 2/3
=
{ 0, 3, 3, 3, 2, 1, 2, 3, 3, 3, 0 }
n
y (n) =
k=-
x (n)= { 0, 3, 2, 1, 0, 1 ,2, 3, 0}
y (1) = 4 . x (1),
y (2) = 4 . x (2),
...
Dynamic system
- All the systems that are not static.
- In dynamic systems, the output may depends on the
present as well as
past input signals (Even future signals)
y (n) = x (n) + x (n-2) is a dynamic system
- as for finding y (4) = x (4) + x (2)
past input is required
- Thus the dynamic system requires memory to store the past inputs
Dynamic system is also called memory system
Examples:
y (n) = [x (n) + x (n-12)]
for y (15), x (15) and x (3) are required
y (n) = x (n) . x (n-1)
n
y (n) = x (n -k)
Requires Finite memory
k=0
n
y (n) = x (n - k) Requires Infinite memory
k= -
+ x (n-2)
past
+ x (n-3)
past)
Non-causal system:
Future
x1 (n)
x2 (n)
a1
a2
+ a2 x2 (n)]
a1
AND
x1 (n)
x2 (n)
T
T
+
a2
As y (n) = T [x (n)]
Hence, the system is n . x (n)
For linearity:
Check if T [x1 (n) + x2 (n)] = T [ x1(n) ] + T[x2 (n)] or
y1 (n)
=
y2 (n)
so y1 (n) = T [x1 (n) + x2 (n)]
= n. [x1 (n) + x2 (n)] = n. x1 (n) + n.x2 (n)
y2 (n) = T [ x1(n) ] + T[x2 (n)] = n . x1(n) + n. x2 (n)
Hence y1 (n)
= y2 (n)
x (n)
y (n)
x (n -k)
y (n -k)
x (n)
- Delay input by k
y (n, k) = e x (n- k)
x (n-k)
bk x (n - k)
and
(1)
elementary signal component
(2)
is ck [ yk (n)]
=T [ ck xk (n) ] = ck T [ xk (n)]
k
k
Using additive property of the linear system
y (n) = T [x (n)]
y (n) = ck yk (n)
k
x (n)
where k is the delay of the unit sample sequence
x (3)
k =3
-3 -2 -1 0 1 2 3
x (n) =
x (k) (n -k)
k=-
= summation of an infinite number of unit sample sequence , (n -k)
having amplitude values of x (k)
= 0.5 . 1 = 0.5
( multiply)
(n)
1
x (0) . (n)
0
0.5
n
-3 -2 -1 0 1 2 3
in general, x (n) =
x (k). (n -k)
k=-
Let x (n ) = {2, 4, 0, 3}
If input x (n) is a unit impulse (n) then the output of the system is
known as impulse response h (n)
(n)
h (n)
Convolution sum
Since input signal x (n) can be represented as
weighted sum of discrete impulses
i. e. x (n) =
x (k). (n- k)
(where n is the time index,
k=-
k is a parameter showing
the location of input impulse
x (k). (n -k)
k=-
- If the system is linear:
Let T [ (n -k)]
y (n) =
= h (n, k)
x (k). h (n, k)
k=-
y (n) =
k=-
x (k). h (n - k)
y (n) = x (k). h (n - k)
k=-
(Output y at n = sum of values of input x (n) at n = k,
multiplied by unit impulse response at n -k,
for all values of k
The convolution sum is represented as
y (n) = x (n)
* h (n)
Properties of Convolution
1. Commutative Law : x (n) * h (n) = h (n) * x (n)
k
k
y (n) = x (k). h (n - k)
k=-
h ( k). x (n -k).
k=-
y (n) =
x (k) . h (n - k)
k=-
( Response y (n) at n, is summation of values of x (n) at k,
multiplied by unit impulse response at n - k, for all values of k)
Calculating the value y (n) for time instant n = 0
y(0) =
x (k). h(0 - k) =
x (k). h( - k)
k=-
k=-
( h ( -k) indicates folding )
For n = 1
y(1) =
x (k). h(1 - k)
k=-
=
k=-
x (k). h(- k + 1)
y(2) =
x (k). h(1 - k)
= x (k). h(- k + 2)
k=-
k=-
h(- k + 2) indicates shifting of folded signal h
h (-k +2) Indicates h (-k), the folded signal is delayed by 2 sample
And so on
h ( -k)
h (-k +1)
3 Multiplication: Multiply
x (k) by h (n k),
Range of n and k , for calculating y (n) =
k=-
Range of n
x (k) . h (n - k)
to
Highest value of n in y( n) = highest value of n in x (n) +
highest value n in h (n)
Range of k
The range k will be same as range of n in x (n)
Input Signal
x (n) = { 1, 2, 3, 1}
to (2 + 3 = 5 ) = -1 to + 5
Range of k = range of n in x (n ) = 0 to 3
Determining the response/ output y:
For n = 0: y (0)
h (n) =
{ 1, 2, 1, -1}
=4
( 0
+1
+4
+3
For n = 2: y (2)
y (2) = x (k) . h (-k +2)
y (2) =
+0) = 8
k varies from 0 to 3
x (k)
. h (-k +2 )
{1, 2, 3, 1} .
{ -1, 1, 2, 1}
= {1, 2, 3, 1} .
{ -1, 1, 2, 1}
= { 1x-1, 2x1,
3x2, 1x 1}
=(
+6
-1
+2
+1) = 8
For n = 3: y (3)
y (3) = x (k) . h (-k +3)
k varies from 0 to 3
x (k)
. h (-k +3 )
{1, 2, 3, 1} .
{ 0, -1, 1, 2, 1}
y (3) =
= {1, 2, 3, 1, 0} .
{ 0, -1, 1, 2, 1}
(0
-2
+3
+2 +0 ) = 3
Similarly y (4) = -2
y( 5) = -1
y ( -1) = 1
Entire response of the system { . . . ., 0, 0, 1, 4, 8, 8, 3, -2, -1, 0, 0, . . . }
h (n) = { 1, 1, 1}
y (n) =
x (k). h (n - k)
k=-
Range of n for y (n) = xL + hL = 0 + 0 = 0 to xh + hh = 2 + 3 = 5
= 0 to 5
Range of k : same as the value of n in x (n) = 0 to 3
x (n) 1 1 0 1
h (n)
1 1
1 0 1
1
1
.1
-2
-2
-2
-2
-2
. -2
-3
-3
-3
-3
-3
.-3
.4
Origin
1 -2 +0 +4, -2 -3 +0, -3+4 , 4}
3,
-5,
1,
4 }
Correlation:
- Correlation is used for the comparison of two signals
- It is measure of degree to which two signals are similar
- Cross correlation
Correlation of two separate signals is known as
cross correlation
- Auto correlation
Correlation of a signal with itself is known as
auto correlation
Cross Correlation:
Correlation between two signals x (n) and y (n)
+
rx y ( l ) = x (n) . y (n - l)
where l = 0, 1, 2, 3 . . .
n=-
Index l is time shift (or lag) parameter and
L
the subscripts x, y on cross correlation sequence rx y ( l )
indicates the sequences being correlated
Auto correlation
The correlation of a signal with itself to determine time delay
between the transmitted signal and the received signal
+
rx x (l) =
n=-
x (n) . x (n - l)
where l = 0, 1, 2, 3 . . .
Range of n in x (n) is -2 to +1
Range of n in y (n) is -2 to +1
Range of l:
since the equation y (n- I) should have maximum value at y (n) = y(1)
nl=1
maximum
Starting value of n in summation is -2
-2 l = 1 for n = -2
I = -1 -2 = -3
or starting value of l = -3
in x (n)
The range of n = -2 to +1
The range for l is -3 to 3
1
rx y(l) =
n = -2
x (n) . y (n - l)
where l = -3 to + 3
=2
x (n) = { 1, 1, 0, 1}
- Fold y (n):
y (n) = { 4, -3, -2, 1}
1 1
1 1
-2 -2
-3 -3
4 4
0 1
0 1
0 -2
0 -3
0 4
-1,
-5,
2,
2,
-3,
4}
4}
x2 (n) be {1, 2, 1, 1}
Folding x2 (n)
x2(-n)= { 1, 1, 2, 1}
x1(n)
1, 2, 1, 1
x2(-n)
1 1 2 1 1
1 1 2 1 1
2 2, 4, 2, 2
1 1, 2, 1, 1
{ 1,
+ 3,1+2+2,
+5, 1+1+4+1,
+ 7,
5, ,+2+1,
3, +11}
}
x x(l)
rx xR(I)
= ={ 1,
2 +1,
1++2+2,
Properties of Correlation
1. The result of autocorrelation is maximum
when the signal matches with itself and there is no phase shifting.
2. Auto-correlation is an even function
rx x (l) = rx x (- l)
3. The cross-correlation is not commutative.
That means rx y(l) ry x (l)
y (n) =
x (k). h (n - k)
k=-
h (k) is the impulse response of the system
FIR
If h (k) is of finite duration, the system is called
Finite Impulse Response (FIR) system
IIR
If h (k) is of infinite duration, the system is called
Infinite Impulse Response (IIR) system
END
2. x (t) = sin(480 t) + 3 sin (720 t) is sampled with Fs = 600 times per sec.
(1) What are the frequencies in radians in the resulting DT signal x (n)?
(2) If x (n) is passed through an ideal interpolator,
what is the reconstructed signal?
To find DT signal frequencies sample the CT signal
Put
t = n Ts = n / Fs = n/600
Ts sampling time period
Fs sampling frequency = 1/ Ts
x ( t) = sin ( 480 t ) + 3 sin (720 t)
x [ n]
= sin ( 480 n / 600) + 3 sin (720 n/ 600)
= sin ( 0.8 n ) + 3 sin ( 1.2 n )
1.2 n = (2 0.8) n = -0.8 n
= sin (0.8 n ) + 3 sin (-0.8 n) = - 2 sin (0.8 n)
= - 2 sin (w n)
thus w = 0.8 radians
To find reconstructed signal put n = t . Fs = 600 t
= - 2 sin ( 0.8 . 600 t) = - 2 sin ( 480 t)