discrete -time signals: time -domain representation
TRANSCRIPT
![Page 1: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/1.jpg)
DiscreteDiscrete--Time Signals:Time Signals:
TimeTime--Domain RepresentationDomain Representation
Tania Stathaki
811b
![Page 2: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/2.jpg)
What is a signal ?
A signal is a function of an independent
variable such as time, distance, position,
temperature, pressure, etc.
![Page 3: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/3.jpg)
For example…
• Electrical Engineering
voltages/currents in a circuit
speech signals
image signals
• Physics
radiation
• Mechanical Engineering
vibration studies
• Astronomy
space photos
![Page 4: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/4.jpg)
or
• BiomedicineEEG, ECG, MRI, X-Rays, Ultrasounds
• Seismologytectonic plate movement, earthquake prediction
• Economicsstock market data
![Page 5: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/5.jpg)
![Page 6: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/6.jpg)
What is DSP?
Mathematical and algorithmic manipulation of discretized and
quantized or naturally digital signals in order to extract the most
relevant and pertinent information that is carried by the signal.
What is a signal?
What is a system?
What is processing?
![Page 7: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/7.jpg)
Signals can be characterized in several ways
Continuous time signals vs. discrete time signals (x(t), x[n]).Temperature in London / signal on a CD-ROM.
Continuous valued signals vs. discrete signals.Amount of current drawn by a device / average scores of TOEFL in a school over years.
–Continuous time and continuous valued : Analog signal.
–Continuous time and discrete valued: Quantized signal.
–Discrete time and continuous valued: Sampled signal.
–Discrete time and discrete values: Digital signal.
Real valued signals vs. complex valued signals.Resident use electric power / industrial use reactive power.
Scalar signals vs. vector valued (multi-channel) signals.Blood pressure signal / 128 channel EEG.
Deterministic vs. random signal: Recorded audio / noise.
One-dimensional vs. two dimensional vs. multidimensional signals.
Speech / still image / video.
![Page 8: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/8.jpg)
![Page 9: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/9.jpg)
Systems
• For our purposes, a DSP system is one that can mathematically
manipulate (e.g., change, record, transmit, transform) digital
signals.
• Furthermore, we are not interested in processing analog signals
either, even though most signals in nature are analog signals.
![Page 10: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/10.jpg)
Various Types of Processing
Modulation and demodulation.
Signal security.
Encryption and decryption.
Multiplexing and de-multiplexing.
Data compression.
Signal de-noising.
Filtering for noise reduction.
Speaker/system identification.
Signal enhancement –equalization.
Audio processing.
Image processing –image de-noising, enhancement,
watermarking.
Reconstruction.
Data analysis and feature extraction.
Frequency/spectral analysis.
![Page 11: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/11.jpg)
Filtering
• By far the most commonly used DSP operationFiltering refers to deliberately changing the frequency content of the signal, typically, by removing certain frequencies from the signals.
For de-noising applications, the (frequency) filter removes those frequencies in the signal that correspond to noise.
In various applications, filtering is used to focus to that part of the spectrum that is of interest, that is, the part that carries the information.
• Typically we have the following types of filtersLow-pass (LPF) –removes high frequencies, and retains (passes) low frequencies.
High-pass (HPF) –removes low frequencies, and retains high frequencies.
Band-pass (BPF) –retains an interval of frequencies within a band, removes others.
Band-stop (BSF) –removes an interval of frequencies within a band, retains others.
Notch filter –removes a specific frequency.
![Page 12: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/12.jpg)
A Common Application: Filtering
![Page 13: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/13.jpg)
Components of a DSP System
![Page 14: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/14.jpg)
Components of a DSP System
![Page 15: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/15.jpg)
Components of a DSP System
![Page 16: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/16.jpg)
Analog-to-Digital-to-Analog…?
• Why not just process the signals in continuous time domain? Isn’t it just a waste of time, money and resources to convert to digital and back to analog?
• Why DSP? We digitally process the signals in discrete domain, because it is
– More flexible, more accurate, easier to mass produce.
– Easier to design.
• System characteristics can easily be changed by programming.
• Any level of accuracy can be obtained by use of appropriate number
of bits.
– More deterministic and reproducible-less sensitive to component values,
etc.
– Many things that cannot be done using analog processors can be done
digitally.
• Allows multiplexing, time sharing, multi-channel processing, adaptive filtering.
• Easy to cascade, no loading effects, signals can be stored indefinitely w/o loss.
• Allows processing of very low frequency signals, which requires unpractical
component values in analog world.
![Page 17: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/17.jpg)
Analog-to-Digital-to-Analog…?
• On the other hand, it can be
– Slower, sampling issues.
– More expensive, increased system complexity,
consumes more power.
• Yet, the advantages far outweigh the
disadvantages. Today, most continuous time
signals are in fact processed in discrete time
using digital signal processors.
![Page 18: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/18.jpg)
Analog-Digital
Examples of analog technology
• photocopiers
• telephones
• audio tapes
• televisions (intensity and color info per scan line)
• VCRs (same as TV)
Examples of digital technology
• Digital computers!
![Page 19: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/19.jpg)
In the next few slides you can see some
real-life signals
![Page 20: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/20.jpg)
Electroencephalogram (EEG) Data
![Page 21: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/21.jpg)
Stock Market Data
![Page 22: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/22.jpg)
Satellite image
Volcano Kamchatka Peninsula, Russia
![Page 23: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/23.jpg)
Satellite image
Volcano in Alaska
![Page 24: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/24.jpg)
Medical Images:
MRI of normal brain
![Page 25: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/25.jpg)
Medical Images:
X-ray knee
![Page 26: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/26.jpg)
Medical Images: Ultrasound
Five-month Foetus (lungs, liver and bowel)
![Page 27: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/27.jpg)
Astronomical images
![Page 28: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/28.jpg)
DiscreteDiscrete--Time Signals:Time Signals:
TimeTime--Domain RepresentationDomain Representation
• Signals represented as sequences of
numbers, called samples
• Sample value of a typical signal or sequence
denoted as x[n] with n being an integer in
the range
• x[n] defined only for integer values of n and
undefined for noninteger values of n
• Discrete-time signal represented by x[n]
∞≤≤∞− n
![Page 29: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/29.jpg)
• Here, n-th sample is given by
• The spacing T is called the sampling interval or
sampling period
• Inverse of sampling interval T, denoted as , is
called the sampling frequency:
DiscreteDiscrete--Time Signals:Time Signals:
TimeTime--Domain RepresentationDomain Representation
),()(][ nTxtxnx anTta == = KK ,1,0,1,2, −−=n
TF1)( −= TFT
![Page 30: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/30.jpg)
DiscreteDiscrete--Time Signals:Time Signals:
TimeTime--Domain RepresentationDomain Representation• Two types of discrete-time signals:
- Sampled-data signals in which samples
are continuous-valued
- Digital signals in which samples are
discrete-valued
• Signals in a practical digital signal
processing system are digital signals
obtained by quantizing the sample values
either by rounding or truncation
![Page 31: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/31.jpg)
256x256 64x64
2 Dimensions
From Continuous to Discrete: Sampling
![Page 32: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/32.jpg)
Discrete (Sampled) and Digital (Quantized) Image
![Page 33: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/33.jpg)
Discrete (Sampled) and Digital (Quantized) Image
![Page 34: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/34.jpg)
256x256 256 levels 256x256 32 levels
Discrete (Sampled) and Digital (Quantized) Image
![Page 35: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/35.jpg)
256x256 256 levels 256x256 2 levels
Discrete (Sampled) and Digital (Quantized) Image
![Page 36: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/36.jpg)
DiscreteDiscrete--Time Signals:Time Signals:
TimeTime--Domain RepresentationDomain Representation
• A discrete-time signal may be a finite-
length or an infinite-length sequence
• Finite-length (also called finite-duration or
finite-extent) sequence is defined only for a
finite time interval:
where and with
• Length or duration of the above finite-
length sequence is
21 NnN ≤≤
1N<∞− ∞<2N 21 NN ≤
112 +−= NNN
![Page 37: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/37.jpg)
DiscreteDiscrete--Time Signals:Time Signals:
TimeTime--Domain RepresentationDomain Representation
• A right-sided sequence x[n] has zero-
valued samples for
• If a right-sided sequence is called a
causal sequence
,01 ≥N
1Nn <
nN1
A right-sided sequence
![Page 38: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/38.jpg)
DiscreteDiscrete--Time Signals:Time Signals:
TimeTime--Domain RepresentationDomain Representation
• A left-sided sequence x[n] has zero-valued
samples for
• If a left-sided sequence is called a
anti-causal sequence
2Nn >
,02 ≤N
2Nn
A left-sided sequence
![Page 39: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/39.jpg)
Operations on SequencesOperations on Sequences
• A single-input, single-output discrete-time
system operates on a sequence, called the
input sequence, according some prescribed
rules and develops another sequence, called
the output sequence, with more desirable
properties
x[n] y[n]
Input sequence Output sequence
Discrete-time
system
![Page 40: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/40.jpg)
Example of an Operation on a Sequence:Example of an Operation on a Sequence:
Noise RemovalNoise Removal
• For example, the input may be a signal
corrupted with additive noise
• A discrete-time system may be designed to
generate an output by removing the noise
component from the input
• In most cases, the operation defining a
particular discrete-time system is composed
of some basic operations
![Page 41: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/41.jpg)
Basic OperationsBasic Operations
• Product (modulation) operation:
– Modulator
• An application is the generation of a finite-length
sequence from an infinite-length sequence by
multiplying the latter with a finite-length sequence
called an window sequence
• Process called windowing
×x[n] y[n]
w[n]][][][ nwnxny ⋅=
![Page 42: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/42.jpg)
Basic OperationsBasic Operations
• Addition operation:
– Adder
• Multiplication operation
– Multiplier
][][][ nwnxny +=
Ax[n] y[n] ][][ nxAny ⋅=
x[n] y[n]
w[n]
+
![Page 43: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/43.jpg)
Basic OperationsBasic Operations
• Time-shifting operation:
where N is an integer
• If N > 0, it is a delay operation
– Unit delay
• If N < 0, it is an advance operation
– Unit advance
][][ Nnxny −=
y[n]x[n] z
1−z y[n]x[n] ][][ 1−= nxny
][][ 1+= nxny
![Page 44: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/44.jpg)
Basic OperationsBasic Operations
• Time-reversal (folding) operation:
• Branching operation: Used to provide
multiple copies of a sequence
][][ nxny −=
x[n] x[n]
x[n]
![Page 45: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/45.jpg)
Combinations of Basic Combinations of Basic
OperationsOperations
• Example -
]3[]2[]1[][][ 4321 −+−+−+= nxnxnxnxny αααα
![Page 46: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/46.jpg)
Sampling Rate AlterationSampling Rate Alteration
• Employed to generate a new sequence y[n]
with a sampling rate higher or lower
than that of the sampling rate of a given
sequence x[n]
• Sampling rate alteration ratio is
• If R > 1, the process called interpolation
• If R < 1, the process called decimation
TF
'TF
T
T
F
FR
'
=
![Page 47: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/47.jpg)
Sampling Rate AlterationSampling Rate Alteration
• In up-sampling by an integer factor L > 1,
equidistant zero-valued samples are
inserted by the up-sampler between each
two consecutive samples of the input
sequence x[n]:
1−L
±±=
=otherwise,0
,2,,0],/[][
LLLnLnxnxu
L][nx ][nxu
![Page 48: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/48.jpg)
Sampling Rate AlterationSampling Rate Alteration
• An example of the up-sampling operation
0 10 20 30 40 50-1
-0.5
0
0.5
1Output sequence up-sampled by 3
Time index n
Amplitude
0 10 20 30 40 50-1
-0.5
0
0.5
1Input Sequence
Time index n
Amplitude
![Page 49: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/49.jpg)
Sampling Rate AlterationSampling Rate Alteration
• In down-sampling by an integer factor
M > 1, every M-th samples of the input
sequence are kept and in-between
samples are removed:
1−M
][][ nMxny =
][nx ][nyM
![Page 50: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/50.jpg)
Sampling Rate AlterationSampling Rate Alteration
• An example of the down-sampling
operation
0 10 20 30 40 50-1
-0.5
0
0.5
1Output sequence down-sampled by 3
Amplitude
Time index n
0 10 20 30 40 50-1
-0.5
0
0.5
1Input Sequence
Time index n
Amplitude
![Page 51: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/51.jpg)
Classification of Sequences Classification of Sequences
Based on SymmetryBased on Symmetry
• Conjugate-symmetric sequence:
If x[n] is real, then it is an even sequence
][*][ nxnx −=
An even sequence
![Page 52: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/52.jpg)
Classification of Sequences Classification of Sequences
Based on SymmetryBased on Symmetry
• Conjugate-antisymmetric sequence:
If x[n] is real, then it is an odd sequence
][*][ nxnx −−=
An odd sequence
![Page 53: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/53.jpg)
Classification of Sequences Classification of Sequences
Based on SymmetryBased on Symmetry
• It follows from the definition that for a
conjugate-symmetric sequence x[n], x[0]
must be a real number
• Likewise, it follows from the definition that
for a conjugate anti-symmetric sequence
y[n], y[0] must be an imaginary number
• From the above, it also follows that for an
odd sequence w[n], w[0] = 0
![Page 54: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/54.jpg)
• Any complex sequence can be expressed as a
sum of its conjugate-symmetric part and its
conjugate-antisymmetric part:
where
Classification of Sequences Classification of Sequences
Based on SymmetryBased on Symmetry
( )][*][][21 nxnxnxcs −+=
( )][*][][21 nxnxnxca −−=
][][][ nxnxnx cacs +=
![Page 55: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/55.jpg)
Classification of Sequences Classification of Sequences
Based on SymmetryBased on Symmetry
• Any real sequence can be expressed as a
sum of its even part and its odd part:
where
][][][ nxnxnx odev +=
( )][][][21 nxnxnxev −+=
( )][][][21 nxnxnxod −−=
![Page 56: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/56.jpg)
Classification of Sequences Classification of Sequences
Based on PeriodicityBased on Periodicity• A sequence satisfying
is called a periodic sequence with a period N where N is
a positive integer and k is any integer
• Smallest value of N satisfying
is called the fundamental period
• A sequence not satisfying the periodicity condition is
called an aperiodic sequence
][~ nx ][~][~ kNnxnx +=
][~][~ kNnxnx +=
![Page 57: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/57.jpg)
Classification of Sequences:Classification of Sequences:
Energy and Power SignalsEnergy and Power Signals
• Total energy of a sequence x[n] is defined by
• An infinite length sequence with finite sample
values may or may not have finite energy
• A finite length sequence with finite sample
values has finite energy
∑=∞
−∞=nnx2
x ][ε
![Page 58: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/58.jpg)
Classification of Sequences:Classification of Sequences:
Energy and Power SignalsEnergy and Power Signals
• The average power of an aperiodic
sequence is defined by
• We define the energy of a sequence x[n]
over a finite interval as
∑=−=
+∞→
K
KnKK
nxP2
12
1
x ][lim
KnK ≤≤−
∑=−=
K
KnKx nx
2
, ][ε
![Page 59: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/59.jpg)
Classification of Sequences:Classification of Sequences:
Energy and Power SignalsEnergy and Power Signals
• The average power of a periodic sequence
with a period N is given by
• The average power of an infinite-length
sequence may be finite or infinite
∑−
==
1
0
21N
nNx nxP ][~
][~ nx
![Page 60: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/60.jpg)
Classification of Sequences:Classification of Sequences:
Energy and Power SignalsEnergy and Power Signals
• Example - Consider the causal sequence
defined by
• Note: x[n] has infinite energy
• Its average power is given by
5.412
)1(9lim19
12
1lim
0
=++
=
∑
+=
∞→=∞→ K
K
KP
K
K
nKx
<≥−=00013
nnnx
n
,,)(][
![Page 61: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/61.jpg)
Classification of Sequences:Classification of Sequences:
Energy and Power SignalsEnergy and Power Signals
• An infinite energy signal with finite average
power is called a power signal
Example - A periodic sequence which has a
finite average power but infinite energy
• A finite energy signal with zero average
power is called an energy signal
![Page 62: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/62.jpg)
Classification of Sequences:Classification of Sequences:
DeterministicDeterministic--StochasticStochastic
![Page 63: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/63.jpg)
Other Types of ClassificationsOther Types of Classifications
• A sequence x[n] is said to be bounded if
• Example - The sequence is a
bounded sequence as
∞<≤ xBnx ][
nnx π= 3.0cos][
13.0cos][ ≤π= nnx
![Page 64: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/64.jpg)
Other Types of ClassificationsOther Types of Classifications
• A sequence x[n] is said to be absolutely
summable if
• Example - The sequence
is an absolutely summable sequence as
∑ ∞<∞
−∞=nnx ][
<≥=00030
nnny
n
,,.][
∞<=−
=∑∞
=428571
301
130
0
..
.n
n
![Page 65: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/65.jpg)
Other Types of ClassificationsOther Types of Classifications
• A sequence x[n] is said to be square-
summable if
• Example - The sequence
is square-summable but not absolutely
summable
∑ ∞<∞
−∞=nnx2][
nn
nh π= 4.0sin][
![Page 66: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/66.jpg)
Basic SequencesBasic Sequences
• Unit sample sequence -
• Unit step sequence -
≠
==
0,0
0,1][
n
nnδ
1
–4 –3 –2 –1 0 1 2 3 4 5 6n
<
≥=
0,0
0,1][
n
nnµ
–4 –3 –2 –1 0 1 2 3 4 5 6
1
n
![Page 67: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/67.jpg)
Basic SequencesBasic Sequences• Real sinusoidal sequence -
where A is the amplitude, is the angular
frequency, and is the phase of x[n]
Example -
)cos(][ φ+ω= nAnx o
oωφ
0 10 20 30 40-2
-1
0
1
2
Time index n
Amplitude
ωo = 0.1
![Page 68: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/68.jpg)
Basic SequencesBasic Sequences• Complex exponential sequence -
where A and are real or complex numbers
• If we write
then we can express
where
,][ nAnx α= ∞<<∞− n
α,)( oo j
eω+σ=α ,φ= jeAA
],[][][)(
nxjnxeeAnx imrenjj oo +== ω+σφ
),cos(][ φ+ω= σneAnx o
nre
o
)sin(][ φ+ω= σneAnx o
nim
o
![Page 69: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/69.jpg)
Basic SequencesBasic Sequences• and of a complex exponential
sequence are real sinusoidal sequences with
constant , growing , and
decaying amplitudes for n > 0
][nxre ][nxim
( )0=σo ( )0>σo
( )0<σo
njnx )exp(][612
1 π+−=
0 10 20 30 40-1
-0.5
0
0.5
1
Time index n
Amplitude
Real part
0 10 20 30 40-1
-0.5
0
0.5
1
Time index n
Amplitude
Imaginary part
![Page 70: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/70.jpg)
Basic SequencesBasic Sequences
• Real exponential sequence -
where A and α are real or complex numbers
,][ nAnx α= ∞<<∞− n
0 5 10 15 20 25 300
10
20
30
40
50
Time index n
Amplitude
α = 1.2
0 5 10 15 20 25 300
5
10
15
20
Time index n
Amplitude
α = 0.9
![Page 71: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/71.jpg)
Basic SequencesBasic Sequences
• Sinusoidal sequence and
complex exponential sequence
are periodic sequences of period N if
where N and r are positive integers
• Smallest value of N satisfying
is the fundamental period of the sequence
• To verify the above fact, consider
)cos( φ+ω nA o
)exp( njB oωrNo π=ω 2
rNo π=ω 2
)cos(][1 φ+ω= nnx o
))(cos(][2 φ++ω= Nnnx o
![Page 72: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/72.jpg)
Basic SequencesBasic Sequences
• Now
which will be equal to
only if
and
• These two conditions are met if and only if
or
))(cos(][2 φ++ω= Nnnx o
NnNn oooo ωφ+ω−ωφ+ω= sin)sin(cos)cos(
][)cos( 1 nxno =φ+ω
0sin =ω No1cos =ω No
rNo π=ω 2rN
o=ω
π2
![Page 73: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/73.jpg)
Basic SequencesBasic Sequences
• If is a noninteger rational number, then
the period will be a multiple of
• Otherwise, the sequence is aperiodic
• Example - is an aperiodic
sequence
oωπ/2
oωπ/2
)3sin(][ φ+= nnx
![Page 74: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/74.jpg)
Basic SequencesBasic Sequences
• Here
• Hence period for r = 0
0=ωo
10
2=
π=
rN
0 10 20 30 400
0.5
1
1.5
2
Time index n
Amplitude
ω0 = 0
![Page 75: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/75.jpg)
Basic SequencesBasic Sequences
• Here
• Hence for r = 1
π=ω 1.0o
201.0
2==
ππr
N
0 10 20 30 40-2
-1
0
1
2
Time index n
Amplitude
ω0 = 0.1π
![Page 76: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/76.jpg)
Basic SequencesBasic Sequences
• Property 1 - Consider and
with and
where k is any positive
integer
• If then x[n] = y[n]
• Thus, x[n] and y[n] are indistinguishable
)exp(][ 1njnx ω=)exp(][ 2njny ω= π<ω≤ 10
)1(22 2 +π<ω≤π kk
,212 kπ+ω=ω
![Page 77: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/77.jpg)
Basic SequencesBasic Sequences
• Property 2 - The frequency of oscillation of
increases as increases from 0
to π, and then decreases as increases from to
• Thus, frequencies in the neighborhood of
are called low frequencies, whereas,
frequencies in the neighborhood of are
called high frequencies
)cos( nA oω oω
oωπ π2
π=ω0=ω
![Page 78: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/78.jpg)
Basic SequencesBasic Sequences
• Because of Property 1, a frequency in
the neighborhood of ω = 2π k is indistinguishable from a frequency
in the neighborhood of ω = 0and a frequency in the neighborhood of
is indistinguishable from a
frequency in the
neighborhood of ω = π
oω
ko π−ω 2
oω
)12( +π−ω ko
)12( +π=ω k
![Page 79: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/79.jpg)
Basic SequencesBasic Sequences
• Frequencies in the neighborhood of ω = 2π kare usually called low frequencies
• Frequencies in the neighborhood of
ω = π (2k+1) are usually called high frequencies
• is a low-
frequency signal
• is a high-
frequency signal
)9.1cos()1.0cos(][1 nnnv π=π=
)2.1cos()8.0cos(][2 nnnv π=π=
![Page 80: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/80.jpg)
Basic SequencesBasic Sequences• An arbitrary sequence can be represented in
the time-domain as a weighted sum of some
basic sequence and its delayed (advanced)
versions
]2[]1[5.1]2[5.0][ −−−++= nnnnx δδδ]6[75.0]4[ −+−+ nn δδ
![Page 81: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/81.jpg)
The Sampling ProcessThe Sampling Process• Often, a discrete-time sequence x[n] is
developed by uniformly sampling a
continuous-time signal as indicated
below
• The relation between the two signals is
)(txa
),()(][ nTxtxnx anTta == =KK ,2,1,0,1,2, −−=n
![Page 82: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/82.jpg)
The Sampling ProcessThe Sampling Process
• Time variable t of is related to the time
variable n of x[n] only at discrete-time
instants given by
with denoting the sampling
frequency and
denoting the sampling angular
frequency
)(txa
TTn
nFnnTt
Ωπ=== 2
nt
TFT /1=
TT Fπ=Ω 2
![Page 83: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/83.jpg)
• Consider the continuous-time signal
• The corresponding discrete-time signal is
where
is the normalized digital angular frequency
of x[n]
The Sampling ProcessThe Sampling Process
)cos()2cos()( φ+Ω=φ+π= tAtfAtx oo
)2
cos()cos(][ φ+ΩΩπ
=φ+Ω= nAnTAnxT
oo
)cos( φ+ω= nA o
ToToo Ω=ΩΩπ=ω /2
Hertz
seconds
radians per second
radians per sample
![Page 84: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/84.jpg)
The Sampling ProcessThe Sampling Process
• If the unit of sampling period T is in
seconds
• The unit of normalized digital angular
frequency is radians/sample
• The unit of normalized analog angular
frequency is radians/second
• The unit of analog frequency is hertz
(Hz)
oω
oΩ
of
![Page 85: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/85.jpg)
The Sampling ProcessThe Sampling Process• The three continuous-time signals
of frequencies 3 Hz, 7 Hz, and 13 Hz, are
sampled at a sampling rate of 10 Hz, i.e.
with T = 0.1 sec. generating the three
sequences
)6cos()(1 ttg π=)14cos()(2 ttg π=)26cos()(3 ttg π=
)6.2cos(][3 nng π=)6.0cos(][1 nng π= )4.1cos(][2 nng π=
![Page 86: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/86.jpg)
The Sampling ProcessThe Sampling Process• Plots of these sequences (shown with circles)
and their parent time functions are shown
below:
• Note that each sequence has exactly the same
sample value for any given n
0 0.2 0.4 0.6 0.8 1-1
-0.5
0
0.5
1
time
Amplitude
![Page 87: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/87.jpg)
The Sampling ProcessThe Sampling Process
• This fact can also be verified by observing that
• As a result, all three sequences are identical
and it is difficult to associate a unique
continuous-time function with each of these
sequences
( ) )6.0cos()6.02(cos)4.1cos(][2 nnnng π=π−π=π=
( ) )6.0cos()6.02(cos)6.2cos(][3 nnnng π=π+π=π=
![Page 88: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/88.jpg)
The Sampling ProcessThe Sampling Process
• The above phenomenon of a continuous-
time signal of higher frequency acquiring
the identity of a sinusoidal sequence of
lower frequency after sampling is called
aliasing
![Page 89: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/89.jpg)
The Sampling ProcessThe Sampling Process
• Since there are an infinite number of
continuous-time signals that can lead to the
same sequence when sampled periodically,
additional conditions need to imposed so
that the sequence can
uniquely represent the parent continuous-
time signal
• In this case, can be fully recovered
from x[n]
)(][ nTxnx a=
)(txa
)(txa
![Page 90: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/90.jpg)
The Sampling ProcessThe Sampling Process
• Example - Determine the discrete-time
signal v[n] obtained by uniformly sampling
at a sampling rate of 200 Hz the continuous-
time signal
• Note: is composed of 5 sinusoidal
signals of frequencies 30 Hz, 150 Hz, 170
Hz, 250 Hz and 330 Hz
)340cos(2)300sin(3)60cos(6)( ttttva π+π+π=)660sin(10)500cos(4 tt π+π+
)(tva
![Page 91: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/91.jpg)
The Sampling ProcessThe Sampling Process
• The sampling period is
• The generated discrete-time signal v[n] is
thus given by
sec005.0200
1 ==T
)7.1cos(2)5.1sin(3)3.0cos(6][ nnnnv π+π+π=
)()( )3.02(cos2)5.02(sin3)3.0cos(6 nnn π−π+π−π+π=
)3.3sin(10)5.2cos(4 nn π+π+
)()( )7.04(sin10)5.02(cos4 nn π−π+π+π+
![Page 92: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/92.jpg)
The Sampling ProcessThe Sampling Process
• Note: v[n] is composed of 3 discrete-time
sinusoidal signals of normalized angular
frequencies: 0.3π, 0.5π, and 0.7π
)5.0cos(4)3.0cos(2)5.0sin(3)3.0cos(6 nnnn π+π+π−π=
)7.0sin(10 nπ−
)7.0sin(10)6435.05.0cos(5)3.0cos(8 nnn π−+π+π=
![Page 93: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/93.jpg)
The Sampling ProcessThe Sampling Process
• Note: An identical discrete-time signal is
also generated by uniformly sampling at a
200-Hz sampling rate the following
continuous-time signals:
)140sin(10)6435.0100cos(5)60cos(8)( ttttaw π−+π+π=
)260sin(10)100cos(4)60cos(2)( ttttag π+π+π=)700sin(3)460cos(6 tt π+π+
![Page 94: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/94.jpg)
The Sampling ProcessThe Sampling Process
• Recall
• Thus if , then the corresponding
normalized digital angular frequency of
the discrete-time signal obtained by
sampling the parent continuous-time
sinusoidal signal will be in the range
• No aliasing
T
oo Ω
Ωπ=ω2
oωoT Ω>Ω 2
π<ω<π−
![Page 95: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/95.jpg)
The Sampling ProcessThe Sampling Process
• On the other hand, if , the
normalized digital angular frequency will
foldover into a lower digital frequency
in the range
because of aliasing
• Hence, to prevent aliasing, the sampling
frequency should be greater than 2
times the frequency of the sinusoidal
signal being sampled
oT Ω<Ω 2
π<ω<π−π⟩ΩΩπ⟨=ω 2/2 Too
TΩoΩ
![Page 96: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/96.jpg)
The Sampling ProcessThe Sampling Process
• Generalization: Consider an arbitrary
continuous-time signal composed of a
weighted sum of a number of sinusoidal
signals
• can be represented uniquely by its
sampled version x[n] if the sampling
frequency is chosen to be greater than 2
times the highest frequency contained in
)(txa
)(txa
TΩ
)(txa
![Page 97: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/97.jpg)
The Sampling ProcessThe Sampling Process
• The condition to be satisfied by the
sampling frequency to prevent aliasing is
called the sampling theorem
• A formal proof of this theorem will be
presented later
![Page 98: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/98.jpg)
DiscreteDiscrete--Time SystemsTime Systems
• A discrete-time system processes a given
input sequence x[n] to generates an output
sequence y[n] with more desirable
properties
• In most applications, the discrete-time
system is a single-input, single-output
system:
SystemtimeDiscrete−
x[n] y[n]
Input sequence Output sequence
![Page 99: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/99.jpg)
DiscreteDiscrete--Time Systems: Time Systems:
ExamplesExamples
• 2-input, 1-output discrete-time systems -
Modulator, adder
• 1-input, 1-output discrete-time systems -
Multiplier, unit delay, unit advance
![Page 100: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/100.jpg)
DiscreteDiscrete--Time Systems: ExamplesTime Systems: Examples
• Accumulator -
• The output y[n] at time instant n is the sum
of the input sample x[n] at time instant n
and the previous output at time
instant which is the sum of all
previous input sample values from to
• The system cumulatively adds, i.e., it
accumulates all input sample values
∑=−∞=
n
xnyl
l][][ ][]1[][][1
nxnynxxn
+−=+∑=−
−∞=ll
]1[ −ny
,1−n
1−n∞−
![Page 101: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/101.jpg)
DiscreteDiscrete--Time Systems:ExamplesTime Systems:Examples
• Accumulator - Input-output relation can
also be written in the form
• The second form is used for a causal input
sequence, in which case is called
the initial condition
∑+∑==
−
−∞=
n
xxny0
1
][][][ll
ll
,][]1[0
∑+−==
n
xyl
l
]1[−y
0≥n
![Page 102: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/102.jpg)
DiscreteDiscrete--Time Systems:ExamplesTime Systems:Examples
• M-point moving-average system -
• Used in smoothing random variations in
data
• An application: Consider
x[n] = s[n] + d[n],
where s[n] is the signal corrupted by a noise
d[n]
∑ −=−
=
1
0
][1][M
k
knxM
ny
![Page 103: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/103.jpg)
d[n] - random signal],)9.0([2][ nnns =
DiscreteDiscrete--Time Systems:ExamplesTime Systems:Examples
0 10 20 30 40 50-2
0
2
4
6
8
Time index n
Amplitude
d[n]s[n]x[n]
0 10 20 30 40 500
1
2
3
4
5
6
7
Time index n
Amplitude
s[n]y[n]
![Page 104: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/104.jpg)
DiscreteDiscrete--Time Systems:ExamplesTime Systems:Examples
• Linear interpolation - Employed to estimate
sample values between pairs of adjacent
sample values of a discrete-time sequence
• Factor-of-4 interpolation
0 1 2
3 4
5 6 7 8 9 10 11 12n
y[n]
![Page 105: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/105.jpg)
DiscreteDiscrete--Time Systems: Time Systems:
ExamplesExamples
• Factor-of-2 interpolator -
• Factor-of-3 interpolator -
( )]1[]1[21][][ ++−+= nxnxnxny uuu
( )]2[]1[31][][ ++−+= nxnxnxny uuu
( )]1[]2[32 ++−+ nxnx uu
![Page 106: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/106.jpg)
DiscreteDiscrete--Time Systems: Time Systems:
ClassificationClassification
• Linear System
• Shift-Invariant System
• Causal System
• Stable System
• Passive and Lossless Systems
![Page 107: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/107.jpg)
Linear DiscreteLinear Discrete--Time SystemsTime Systems
• Definition - If is the output due to an
input and is the output due to an
input then for an input
the output is given by
• Above property must hold for any arbitrary
constants and and for all possible
inputs and
][1 ny
][1 nx
][2 nx
][2 ny
][][][ 21 nxnxnx βα +=
][][][ 21 nynyny βα +=
α ,β][1 nx ][2 nx
![Page 108: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/108.jpg)
AccumulatorAccumulator::
Linear DiscreteLinear Discrete--Time System?Time System?
• Accumulator -
• For an input
the output is
• Hence, the above system is linear
∑=∑=−∞=−∞=
nn
xnyxnyll
ll ][][,][][ 2211
][][][ 21 nxnxnx βα +=
( )∑ +=−∞=
n
xxnyl
ll ][][][ 21 βα
][][][][ 2121 nynyxxnn
βαβα +=∑+∑=−∞=−∞= ll
ll
![Page 109: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/109.jpg)
Causal AccumulatorCausal Accumulator::
Linear DiscreteLinear Discrete--Time System?Time System?
• The outputs and for inputs
and are given by
• The output y[n] for an input
is given by
∑=
+−=n
xyny0111 1
l
l][][][
∑=
+−=n
xyny0222 1
l
l][][][
][ny1 ][ny2][nx2
][nx1
][][ nxnx 21 βα +
∑=
++−=n
xxyny0
211l
ll ])[][(][][ βα
![Page 110: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/110.jpg)
Causal AccumulatorCausal Accumulator cont.:cont.:
Linear DiscreteLinear Discrete--Time System?Time System?
• Now
• Thus if
][][ nyny 21 βα +
)][][( ∑=
+−+n
xy022 1
l
lβ)][][( ∑=
+−=n
xy011 1
l
lα
)][][(])[][( ∑∑==
++−+−=nn
xxyy02
0121 11
ll
ll βαβα
][][][ nynyny 21 βα +=
][][][ 111 21 −+−=− yyy βα
![Page 111: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/111.jpg)
Causal Accumulator Causal Accumulator cont.:cont.:
Linear DiscreteLinear Discrete--Time System?Time System?
• For the causal accumulator to be linear the
condition
must hold for all initial conditions ,
, , and all constants α and β• This condition cannot be satisfied unless the
accumulator is initially at rest with zero
initial condition
• For nonzero initial condition, the system is
nonlinear
][][][ 111 21 −+−=− yyy βα][ 1−y
][ 11 −y ][ 12 −y
![Page 112: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/112.jpg)
A Nonlinear DiscreteA Nonlinear Discrete--Time SystemTime System
• Consider
• Outputs and for inputs
and are given by
][][][][ 112 +−−= nxnxnxny
][ny1
][][][][ 11 11211 +−−= nxnxnxny
][ny2
][][][][ 11 22222 +−−= nxnxnxny
][nx2
][nx1
![Page 113: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/113.jpg)
A Nonlinear DiscreteA Nonlinear Discrete--Time System cont.Time System cont.
• Output y[n] due to an input
is given by
][][ nxnx 21 βα +
221 ][][][ nxnxny βα +=
][][][][ 1111 2121 +++−+−− nxnxnxnx βαβα
][][][ 11 1121
2 +−−= nxnxnxα
][][][ 11 2222
2 +−−+ nxnxnxβ
][][][][][][ 11112 212121 −+−+−−+ nxnxnxnxnxnxαβ
![Page 114: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/114.jpg)
A Nonlinear DiscreteA Nonlinear Discrete--Time System cont.Time System cont.
• On the other hand
• Hence, the system is nonlinear
][][ nyny 21 βα +
][][][ 11 1121 +−−= nxnxnxα
][][][ 11 2222 +−−+ nxnxnxβ][ny≠
![Page 115: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/115.jpg)
Shift (Time)Shift (Time)--Invariant SystemInvariant System
• For a shift-invariant system, if is the
response to an input , then the response to
an input .
is simply
where is any positive or negative integer
• The above relation must hold for any arbitrary
input and its corresponding output
• If n is discrete time, the above property is called
time-invariance property
][ny1][nx1][][ onnxnx −= 1
][][ onnyny −= 1
on
![Page 116: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/116.jpg)
UpUp--SamplerSampler::
ShiftShift--Invariant System?Invariant System?
• Example - Consider the up-sampler with an
input-output relation given by
• For an input the output
is given by
±±==
otherwise,
.....,,,],/[][
0
20 LLnLnxnxu
][][ onnxnx −=1 ][, nx u1
±±==
otherwise,
.....,,,],/[][, 0
2011
LLnLnxnx u
±±=−=
otherwise,
.....,,,],/)[(
0
20 LLnLLnnx o
![Page 117: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/117.jpg)
UpUp--SamplerSampler::
ShiftShift--Invariant System?Invariant System?
• However from the definition of the up-sampler
• Hence, the up-sampler is a time-varying system
][ ou nnx −
±±=−=
otherwise,
.....,,,],/)[(
0
2LnLnnnLnnx oooo
][, nx u1≠
![Page 118: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/118.jpg)
Linear TimeLinear Time--Invariant SystemInvariant System
• Linear Time-Invariant (LTI) System -
A system satisfying both the linearity and
the time-invariance property
• LTI systems are mathematically easy to
analyze and characterize, and consequently,
easy to design
• Highly useful signal processing algorithms
have been developed utilizing this class of
systems over the last several decades
![Page 119: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/119.jpg)
Causal SystemCausal System
• In a causal system, the -th output sample
depends only on input samples x[n]
for and does not depend on input
samples for
• Let and be the responses of a
causal discrete-time system to the inputs
and , respectively
on
onn ≤onn >
][ ony
][ny1 ][ny2
][nx2
][nx1
![Page 120: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/120.jpg)
Causal SystemCausal System
• Then
for n < N
implies also that
for n < N
• For a causal system, changes in output
samples do not precede changes in the input
samples
][][ 21 nxnx =
][][ 21 nyny =
![Page 121: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/121.jpg)
Causal SystemCausal System• Examples of causal systems:
• Examples of noncausal systems:
][][][][][ 321 4321 −+−+−+= nxnxnxnxny αααα
][][][][ 21 210 −+−+= nxbnxbnxbny
][][ 21 21 −+−+ nyanya
][][][ nxnyny +−= 1
])[][(][][ 112
1 ++−+= nxnxnxny uuu
![Page 122: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/122.jpg)
Causal SystemCausal System
• A noncausal system can be implemented as
a causal system by delaying the output by
an appropriate number of samples
• For example a causal implementation of the
factor-of-2 interpolator is given by
])[][(][][ nxnxnxny uuu +−+−= 212
1
![Page 123: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/123.jpg)
Stable SystemStable System
• There are various definitions of stability
• We consider here the bounded-input,
bounded-output (BIBO) stability
• If y[n] is the response to an input x[n] and if
for all values of n
then
for all values of n
xBnx ≤][
yBny ≤][
![Page 124: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/124.jpg)
Stable SystemStable System
• Example - The M-point moving average
filter is BIBO stable:
• For a bounded input we have
∑−
=−=
1
0
1M
kM
knxny ][][
xBnx ≤][
∑∑−
=
−
=−≤−=
1
0
11
0
1M
kM
M
kM
knxknxny ][][][
xxMBMB ≤≤ )(
1
![Page 125: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/125.jpg)
Passive and Lossless SystemsPassive and Lossless Systems
• A discrete-time system is defined to be
passive if, for every finite-energy input x[n],
the output y[n] has, at most, the same energy,
i.e.
• For a lossless system, the above inequality is
satisfied with an equal sign for every input
∞<≤ ∑∑∞
−∞=
∞
−∞= nn
nxny22][][
![Page 126: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/126.jpg)
Passive and Lossless SystemsPassive and Lossless Systems
• Example - Consider the discrete-time
system defined by with N
a positive integer
• Its output energy is given by
• Hence, it is a passive system if and is
a lossless system if
][][ Nnxny −=α
∑∑∞
−∞=
∞
−∞=≤
nn
nxny222][][ α
1≤α1=α
![Page 127: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/127.jpg)
Impulse and Step ResponsesImpulse and Step Responses
• The response of a discrete-time system to a
unit sample sequence δ[n] is called the unit sample response or simply, the
impulse response, and is denoted by h[n]
• The response of a discrete-time system to a
unit step sequence µ[n] is called the unitstep response or simply, the step response,
and is denoted by s[n]
![Page 128: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/128.jpg)
Impulse ResponseImpulse Response
• Example - The impulse response of the
system
is obtained by setting x[n] = δ[n] resulting in
• The impulse response is thus a finite-length
sequence of length 4 given by
][][][][][ 321 4321 −+−+−+= nxnxnxnxny αααα
][][][][][ 321 4321 −+−+−+= nnnnnh δαδαδαδα
,,,][ 4321 αααα↑
=nh
![Page 129: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/129.jpg)
Impulse ResponseImpulse Response
• Example - The impulse response of the
discrete-time accumulator
is obtained by setting x[n] = δ[n] resulting in
∑−∞=
=n
xnyl
l][][
][][][ nnhn
µδ == ∑−∞=l
l
![Page 130: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/130.jpg)
Impulse ResponseImpulse Response
• Example - The impulse response h[n] of
the factor-of-2 interpolator
is obtained by setting and is
given by
• The impulse response is thus a finite-length
sequence of length 3:
])[][(][][ 112
1 ++−+= nxnxnxny uuu
])[][(][][ 112
1 ++−+= nnnnh δδδ
.,.][ 50150↑
=nh
][][ nnxu δ=
![Page 131: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/131.jpg)
TimeTime--Domain Characterization Domain Characterization
of LTI Discreteof LTI Discrete--Time SystemTime System
• Input-Output Relationship -
It can be shown that a consequence of the
linear, time-invariance property is that an
LTI discrete-time system is completely
characterized by its impulse response
• Knowing the impulse response one
can compute the output of the system for
any arbitrary input
![Page 132: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/132.jpg)
TimeTime--Domain Characterization Domain Characterization of LTI Discreteof LTI Discrete--Time SystemTime System
• Let h[n] denote the impulse response of a
LTI discrete-time system
• We compute its output y[n] for the input:
• As the system is linear, we can compute its
outputs for each member of the input
separately and add the individual outputs to
determine y[n]
]5[75.0]2[]1[5.1]2[5.0][ −δ+−δ−−δ++δ= nnnnnx
![Page 133: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/133.jpg)
TimeTime--Domain Characterization Domain Characterization of LTI Discreteof LTI Discrete--Time SystemTime System
• Since the system is time-invariant
input output
]2[]2[ +→+δ nhn
]1[]1[ −→−δ nhn
]2[]2[ −→−δ nhn
]5[]5[ −→−δ nhn
![Page 134: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/134.jpg)
TimeTime--Domain Characterization Domain Characterization of LTI Discreteof LTI Discrete--Time SystemTime System
• Likewise, as the system is linear
• Hence because of the linearity property we
get
]5[75.0]5[75.0 −→−δ nhn
input output
]2[5.0]2[5.0 +→+δ nhn
]2[]2[ −−→−δ− nhn
]1[5.1]1[5.1 −→−δ nhn
][.][.][ 151250 −++= nhnhny
][.][ 57502 −+−− nhnh
![Page 135: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/135.jpg)
TimeTime--Domain Characterization Domain Characterization
of LTI Discreteof LTI Discrete--Time SystemTime System
• Now, any arbitrary input sequence x[n] can
be expressed as a linear combination of
delayed and advanced unit sample
sequences in the form
• The response of the LTI system to an input
will be
∑ −δ=∞
−∞=kknkxnx ][][][
][][ knkx −δ ][][ knhkx −
![Page 136: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/136.jpg)
TimeTime--Domain Characterization Domain Characterization
of LTI Discreteof LTI Discrete--Time SystemTime System
• Hence, the response y[n] to an input
will be
which can be alternately written as
∑ −δ=∞
−∞=kknkxnx ][][][
∑ −=∞
−∞=kknhkxny ][][][
∑∞
−∞=−=
k
khknxny ][][][
![Page 137: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/137.jpg)
Convolution SumConvolution Sum
• The summation
is called the convolution sum of the
sequences x[n] and h[n] and represented
compactly as
∑∑∞
−∞=
∞
−∞=−=−=
kk
nhknxknhkxny ][][][][][
y[n] = x[n] h[n]*
![Page 138: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/138.jpg)
Convolution SumConvolution Sum
• Properties -
• Commutative property:
• Associative property :
• Distributive property :
x[n] h[n] = h[n] x[n]* *
(x[n] h[n]) y[n] = x[n] (h[n] y[n])****
x[n] (h[n] + y[n]) = x[n] h[n] + x[n] y[n]** *
![Page 139: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/139.jpg)
Simple Interconnection Simple Interconnection
SchemesSchemes
• Two simple interconnection schemes are:
• Cascade Connection
• Parallel Connection
![Page 140: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/140.jpg)
Cascade ConnectionCascade Connection
• Impulse response h[n] of the cascade of two
LTI discrete-time systems with impulse
responses and is given by
][nh1][nh2][nh1 ][nh2≡
][][ nhnh 1=][nh2][nh1 *≡
][nh1 ][nh2
][nh2][][ nhnh 1= *
![Page 141: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/141.jpg)
Cascade ConnectionCascade Connection
• Note: The ordering of the systems in the
cascade has no effect on the overall impulse
response because of the commutative
property of convolution
• A cascade connection of two stable systems
is stable
• A cascade connection of two passive
(lossless) systems is passive (lossless)
![Page 142: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/142.jpg)
Cascade ConnectionCascade Connection
• An application is in the development of an
inverse system
• If the cascade connection satisfies the
relation
then the LTI system is said to be the
inverse of and vice-versa
][nh1][nh2
][nh2][1 nh ][nδ=*
![Page 143: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/143.jpg)
Cascade ConnectionCascade Connection• An application of the inverse system
concept is in the recovery of a signal x[n]
from its distorted version appearing at
the output of a transmission channel
• If the impulse response of the channel is
known, then x[n] can be recovered by
designing an inverse system of the channel
][ˆ nx
][nh2][nh1][nx ][nx
channel inverse system][nx
][nh2][nh1 ][nδ=*
![Page 144: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/144.jpg)
Cascade ConnectionCascade Connection
• Example - Consider the discrete-time
accumulator with an impulse response µ[n]• Its inverse system satisfy the condition
• It follows from the above that for
n < 0 and
for
02 =][nh
112 =][h
002 =∑
=
n
hl
l][ 2≥n
][nh2][nµ ][nδ=*
![Page 145: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/145.jpg)
Cascade ConnectionCascade Connection
• Thus the impulse response of the inverse
system of the discrete-time accumulator is
given by
which is called a backward difference
system
]1[][][2 −δ−δ= nnnh
![Page 146: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/146.jpg)
Parallel Connection
• Impulse response h[n] of the parallel
connection of two LTI discrete-time
systems with impulse responses and
is given by
][nh2
][nh1+ ][][ nhnh 1=
][nh2][nh1≡ +
][nh1][nh2
][][][ nhnhnh 21 +=
![Page 147: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/147.jpg)
Simple Interconnection SchemesSimple Interconnection Schemes
• Consider the discrete-time system where
][nh2
][nh1 +
+
][nh4
][nh3
],1[5.0][][1 −δ+δ= nnnh
],1[25.0][5.0][2 −δ−δ= nnnh
],[2][3 nnh δ=
][)5.0(2][4 nnh nµ−=
![Page 148: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/148.jpg)
Simple Interconnection SchemesSimple Interconnection Schemes
• Simplifying the block-diagram we obtain
][nh2
][nh1 +
][][ 43 nhnh +
][nh1 +
])[][(][ 432 nhnhnh +*
≡
![Page 149: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/149.jpg)
Simple Interconnection SchemesSimple Interconnection Schemes
• Overall impulse response h[n] is given by
• Now,
][][][][][ nhnhnhnhnh 42321 ++=
])[][(][][][ nhnhnhnhnh 4321 ++= *
* *
][2])1[][(][][41
21
32 nnnnhnh δ−δ−δ=
]1[][21 −δ−δ= nn
* *
![Page 150: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/150.jpg)
Simple Interconnection SchemesSimple Interconnection Schemes
• Therefore
]1[)(][)( 121
21
21 −µ+µ−= − nn nn
]1[)(][)(21
21 −µ+µ−= nn nn
][][)(21 nnn δ−=δ−=
][][]1[][]1[][][21
21 nnnnnnnh δ=δ−−δ−δ+−δ+δ=
( )][)(2])1[][(][][21
41
21
42 nnnnhnh nµ−−δ−δ=* *
![Page 151: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/151.jpg)
BIBO Stability Condition of an BIBO Stability Condition of an
LTI DiscreteLTI Discrete--Time SystemTime System
• BIBO Stability Condition - A discrete-
time is BIBO stable if the output sequence
y[n] remains bounded for all bounded
input sequence x[n]
• An LTI discrete-time system is BIBO stable
if and only if its impulse response sequence
h[n] is absolutely summable, i.e.
∞<= ∑∞
−∞=nnh ][S
![Page 152: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/152.jpg)
BIBO Stability Condition of an BIBO Stability Condition of an
LTI DiscreteLTI Discrete--Time SystemTime System
• Proof: Assume h[n] is a real sequence
• Since the input sequence x[n] is bounded we
have
• Therefore
∞<≤ xBnx ][
][][][][][ knxkhknxkhnykk
−≤−= ∑∑∞
−∞=
∞
−∞=
SBkhB xk
x =≤ ∑∞
−∞=][
![Page 153: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/153.jpg)
BIBO Stability Condition of an BIBO Stability Condition of an
LTI DiscreteLTI Discrete--Time SystemTime System
• Thus, S < implies
indicating that y[n] is also bounded
• To prove the converse, assume y[n] is
bounded, i.e.,
• Consider the input given by
∞ ∞<≤ yBny ][
yBny ≤][
=−≠−−=0
0
][if,
][if]),[sgn(][
nhK
nhnhnx
![Page 154: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/154.jpg)
BIBO Stability Condition of an BIBO Stability Condition of an
LTI DiscreteLTI Discrete--Time SystemTime System
where sgn(c) = +1 if c > 0 and sgn(c) =
if c < 0 and
• Note: Since , x[n] is obviously
bounded
• For this input, y[n] at n = 0 is
• Therefore, implies S <
1−1≤K
yBny ≤][ ∞
1≤][nx
Skhkhyk
∑∞
−∞=== ][])[sgn(]0[ ∞<≤ yB
![Page 155: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/155.jpg)
• Example - Consider a causal LTI discrete-time system with an impulse response
• For this system
,
• Therefore S < if for which the system is BIBO stable
• If , the system is not BIBO stable
Stability Condition of an LTI Stability Condition of an LTI DiscreteDiscrete--Time SystemTime System
∞ 1<α
1=α
][)(][ nnh nµα=
1<αα
αµα−
=== ∑∑∞
=
∞
−∞= 1
1][
0n
n
n
n nS
![Page 156: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/156.jpg)
Causality Condition of an LTI Causality Condition of an LTI
DiscreteDiscrete--Time SystemTime System
• Let and be two input sequences
with
• The corresponding output samples at
of an LTI system with an impulse response
h[n] are then given by
][nx1 ][nx2
][][ nxnx 21 = onn ≤for
onn =
![Page 157: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/157.jpg)
Causality Condition of an LTI Causality Condition of an LTI
DiscreteDiscrete--Time SystemTime System
∑∑∞
=
∞
−∞=−=−=
0222
ko
koo knxkhknxkhny ][][][][][
∑−
−∞=−+
1
2k
o knxkh ][][
∑∑∞
=
∞
−∞=−=−=
0111
ko
koo knxkhknxkhny ][][][][][
∑−
−∞=−+
1
1k
o knxkh ][][
![Page 158: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/158.jpg)
Causality Condition of an LTI Causality Condition of an LTI
DiscreteDiscrete--Time SystemTime System
• If the LTI system is also causal, then
• As
• This implies
][][ nxnx 21 = onn ≤for
][][ oo nyny 21 =
∑∑∞
=
∞
=−=−
02
01
ko
ko knxkhknxkh ][][][][
∑∑−
−∞=
−
−∞=−=−
1
2
1
1k
ok
o knxkhknxkh ][][][][
![Page 159: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/159.jpg)
Causality Condition of an LTI Causality Condition of an LTI
DiscreteDiscrete--Time SystemTime System
• As for the only way
the condition
will hold if both sums are equal to zero,
which is satisfied if
][][ nxnx 21 ≠ onn >
∑∑−
−∞=
−
−∞=−=−
1
2
1
1k
ok
o knxkhknxkh ][][][][
0=][kh for k < 0
![Page 160: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/160.jpg)
Causality Condition of an LTI Causality Condition of an LTI
DiscreteDiscrete--Time SystemTime System
• An LTI discrete-time system is causal
if and only if its impulse response h[n] is a
causal sequence
• Example - The discrete-time system defined
by
is a causal system as it has a causal impulse
response
]3[]2[]1[][][ 4321 −α+−α+−α+α= nxnxnxnxny
][ 4321 αααα=nh↑
![Page 161: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/161.jpg)
Causality Condition of an LTI Causality Condition of an LTI
DiscreteDiscrete--Time SystemTime System
• Example - The discrete-time accumulator
defined by
is a causal system as it has a causal impulse
response given by
][][][ nnyn
µ=δ= ∑−∞=l
l
][][][ nnhn
µ=δ= ∑−∞=l
l
![Page 162: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/162.jpg)
Causality Condition of an LTI Causality Condition of an LTI
DiscreteDiscrete--Time SystemTime System
• Example - The factor-of-2 interpolator
defined by
is noncausal as it has a noncausal impulse
response given by
( )]1[]1[][][21 ++−+= nxnxnxny uuu
5.015.0][ =nh↑
![Page 163: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/163.jpg)
Causality Condition of an LTI Causality Condition of an LTI
DiscreteDiscrete--Time SystemTime System
• Note: A noncausal LTI discrete-time system
with a finite-length impulse response can
often be realized as a causal system by
inserting an appropriate amount of delay
• For example, a causal version of the factor-
of-2 interpolator is obtained by delaying the
input by one sample period:
( )][]2[]1[][21 nxnxnxny uuu +−+−=
![Page 164: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/164.jpg)
FiniteFinite--Dimensional LTI Dimensional LTI DiscreteDiscrete--Time SystemsTime Systems
• An important subclass of LTI discrete-time
systems is characterized by a linear constant
coefficient difference equation of the form
• x[n] and y[n] are, respectively, the input and
the output of the system
• and are constants characterizing
the system
kd kp
∑∑==
−=−M
kk
N
kk knxpknyd
00
][][
![Page 165: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/165.jpg)
FiniteFinite--Dimensional LTI Dimensional LTI
DiscreteDiscrete--Time SystemsTime Systems
• The order of the system is given by
max(N,M), which is the order of the difference
equation
• It is possible to implement an LTI system
characterized by a constant coefficient
difference equation as here the computation
involves two finite sums of products
![Page 166: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/166.jpg)
FiniteFinite--Dimensional LTI Dimensional LTI DiscreteDiscrete--Time SystemsTime Systems
• If we assume the system to be causal, then
the output y[n] can be recursively computed
using
provided
• y[n] can be computed for all ,
knowing x[n] and the initial conditions
00 ≠d
onn ≥
][,...],2[],1[ Nnynyny ooo −−−
][][][1 01 0
knxd
pkny
d
dny
M
k
kN
k
k −+−−= ∑∑==
![Page 167: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/167.jpg)
Classification of LTI DiscreteClassification of LTI Discrete--Time SystemsTime Systems
Based on Impulse Response Length -
• If the impulse response h[n] is of finite
length, i.e.,
then it is known as a finite impulse
response (FIR) discrete-time system
• The convolution sum description here is
2121 ,andfor0][ NNNnNnnh <><=
∑=
−=2
1
][][][N
Nk
knxkhny
![Page 168: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/168.jpg)
Classification of LTI DiscreteClassification of LTI Discrete--
Time SystemsTime Systems
• The output y[n] of an FIR LTI discrete-time
system can be computed directly from the
convolution sum as it is a finite sum of
products
• Examples of FIR LTI discrete-time systems
are the moving-average system and the
linear interpolators
![Page 169: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/169.jpg)
Classification of LTI DiscreteClassification of LTI Discrete--
Time SystemsTime Systems
• If the impulse response is of infinite length,
then it is known as an infinite impulse
response (IIR) discrete-time system
• The class of IIR systems we are concerned
with in this course are characterized by
linear constant coefficient difference
equations
![Page 170: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/170.jpg)
Classification of LTI DiscreteClassification of LTI Discrete--
Time SystemsTime Systems
• Example - The discrete-time accumulator
defined by
is seen to be an IIR system
][]1[][ nxnyny +−=
![Page 171: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/171.jpg)
Classification of LTI DiscreteClassification of LTI Discrete--
Time SystemsTime Systems
• Example - The familiar numerical
integration formulas that are used to
numerically solve integrals of the form
can be shown to be characterized by linear
constant coefficient difference equations,
and hence, are examples of IIR systems
∫ ττ=t
dxty
0
)()(
![Page 172: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/172.jpg)
Classification of LTI DiscreteClassification of LTI Discrete--
Time SystemsTime Systems
• If we divide the interval of integration into n
equal parts of length T, then the previous
integral can be rewritten as
where we have set t = nT and used the
notation
∫−
ττ+−=nT
Tn
dxTnynTy
)1(
)())1(()(
∫ ττ=nT
dxnTy
0
)()(
![Page 173: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/173.jpg)
Classification of LTI DiscreteClassification of LTI Discrete--
Time SystemsTime Systems
• Using the trapezoidal method we can write
• Hence, a numerical representation of the
definite integral is given by
)())1(()(2
)1(
nTxTnxdx TnT
Tn
+−=ττ∫−
)())1(())1(()(2
nTxTnxTnynTy T +−+−=
![Page 174: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/174.jpg)
Classification of LTI DiscreteClassification of LTI Discrete--
Time SystemsTime Systems
• Let y[n] = y(nT) and x[n] = x(nT)
• Then
reduces to
which is recognized as the difference
equation representation of a first-order IIR
discrete-time system
)())1(())1(()(2
nTxTnxTnynTy T +−+−=
]1[][]1[][2
−++−= nxnxnyny T
![Page 175: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/175.jpg)
Classification of LTI DiscreteClassification of LTI Discrete--
Time SystemsTime Systems
Based on the Output Calculation Process
• Nonrecursive System - Here the output can
be calculated sequentially, knowing only
the present and past input samples
• Recursive System - Here the output
computation involves past output samples in
addition to the present and past input
samples
![Page 176: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/176.jpg)
Classification of LTI DiscreteClassification of LTI Discrete--
Time SystemsTime Systems
Based on the Coefficients -
• Real Discrete-Time System - The impulse
response samples are real valued
• Complex Discrete-Time System - The
impulse response samples are complex
valued
![Page 177: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/177.jpg)
Correlation of SignalsCorrelation of Signals
Definitions
• A measure of similarity between a pair of
energy signals, x[n] and y[n], is given by the
cross-correlation sequence defined by
• The parameter called lag, indicates the
time-shift between the pair of signals
l
][lxyr
...,,,],[][][ 210 ±±=−= ∑∞
−∞=lll
nxy nynxr
![Page 178: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/178.jpg)
Correlation of SignalsCorrelation of Signals
• If y[n] is made the reference signal and we
wish to shift x[n] with respect to y[n], then
the corresponding cross-correlation
sequence is given by
• Thus, is obtained by time-reversing
∑∞−∞= −=
nyx nxnyr ][][][ ll
][][][ ll −=+=∑∞−∞= xym
rmxmy
][lyxr][lxyr
![Page 179: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/179.jpg)
Correlation of SignalsCorrelation of Signals
• The autocorrelation sequence of x[n] is
given by
obtained by setting y[n] = x[n] in the
definition of the cross-correlation sequence
• Note: , the energy
of the signal x[n]
∑∞−∞= −=
nxx nxnxr ][][][ ll
][lxyr
xnxx nxr E∑∞−∞= == ][][ 20
![Page 180: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/180.jpg)
Correlation of SignalsCorrelation of Signals
• From the relation it follows
that implying that is
an even function for real x[n]
• An examination of
reveals that the expression for the cross-
correlation looks quite similar to that of the
linear convolution
][][ ll −= xyyx rr
][][ ll −= xxxx rr ][lxxr
∑∞−∞= −=
nxy nynxr ][][][ ll
![Page 181: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/181.jpg)
Correlation of SignalsCorrelation of Signals
• This similarity is much clearer if we rewrite
the expression for the cross-correlation as
• The cross-correlation of y[n] with the
reference signal x[n] can be computed by
processing x[n] with an LTI discrete-time
system of impulse response ][ ny −
][][)]([][][ llll −=−−=∑∞−∞= yxnynxr
nxy *
][ ny −][nx ][nrxy
![Page 182: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/182.jpg)
Correlation of SignalsCorrelation of Signals
• Likewise, the autocorrelation of x[n] can be
computed by processing x[n] with an LTI
discrete-time system of impulse response
][ nx −][nx ][nrxx
][ nx −
![Page 183: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/183.jpg)
Properties of Autocorrelation and Properties of Autocorrelation and
CrossCross--correlation Sequencescorrelation Sequences
• Consider two finite-energy sequences x[n]
and y[n]
• The energy of the combined sequence
is also finite and
nonnegative, i.e.,
][][ l−+ nynxa
∑∑ ∞−∞=
∞−∞= =−+
nnnxanynxa ][])[][( 222
l
02 2 ≥−+−+ ∑∑ ∞−∞=
∞−∞= nn
nynynxa ][][][ ll
![Page 184: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/184.jpg)
Properties of Autocorrelation and Properties of Autocorrelation and
CrossCross--correlation Sequencescorrelation Sequences
• Thus
where and
• We can rewrite the equation on the previous
slide as
for any finite value of a
00202 ≥++ ][][][ yyxyxx rrara l
00 >= xxxr E][ 00 >= yyyr E][
[ ] 010
01 ≥
arr
rra
yyxy
xyxx
][][
][][
l
l
![Page 185: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/185.jpg)
Properties of Autocorrelation and Properties of Autocorrelation and
CrossCross--correlation Sequencescorrelation Sequences
• Or, in other words, the matrix
is positive semidefinite
•
][][
][][
0
0
yyxy
xyxx
rr
rr
l
l
000 2 ≥− ][][][ lxyyyxx rrr
yxyyxxxy rrr EE=≤ ][][|][| 00l
![Page 186: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/186.jpg)
Properties of Autocorrelation and Properties of Autocorrelation and
CrossCross--correlation Sequencescorrelation Sequences
• The last inequality on the previous slide
provides an upper bound for the cross-
correlation samples
• If we set y[n] = x[n], then the inequality
reduces to
xxxxx rr E]0[|][| =≤l
![Page 187: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/187.jpg)
Properties of Autocorrelation and Properties of Autocorrelation and
CrossCross--correlation Sequencescorrelation Sequences
• Thus, at zero lag ( ), the sample value
of the autocorrelation sequence has its
maximum value
• Now consider the case
where N is an integer and b > 0 is an
arbitrary number
• In this case
0=l
xy b EE 2=
][][ Nnxbny −±=
![Page 188: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/188.jpg)
Properties of Autocorrelation and Properties of Autocorrelation and
CrossCross--correlation Sequencescorrelation Sequences
• Therefore
• Using the above result in
we get
xxyx bb EEEE == 22
yxyyxxxy rrr EE=≤ ][][|][| 00l
][][][ 00 xxxyxx rbrrb ≤≤− l
![Page 189: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/189.jpg)
Correlation Computation Correlation Computation
Using MATLABUsing MATLAB
• The cross-correlation and autocorrelation
sequences can easily be computed using
MATLAB
• Example - Consider the two finite-length
sequences
[ ]244121231 −−=][nx
[ ]321412 −−=][ny
![Page 190: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/190.jpg)
Correlation Computation Correlation Computation
Using MATLABUsing MATLAB
• The cross-correlation sequence
computed using Program 2_7 of text is
plotted below
][nrxy
-4 -2 0 2 4 6 8
-10
0
10
20
30
Lag index
Amplitude
![Page 191: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/191.jpg)
Correlation Computation Correlation Computation
Using MATLABUsing MATLAB
• The autocorrelation sequence
computed using Program 2_7 is shown below
• Note: At zero lag, is the maximum
][lxxr
][0xxr
-5 0 5-20
0
20
40
60
Lag index
Amplitude
![Page 192: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/192.jpg)
Correlation Computation Correlation Computation
Using MATLABUsing MATLAB
• The plot below shows the cross-correlation
of x[n] and for N = 4
• Note: The peak of the cross-correlation is
precisely the value of the delay N
][][ Nnxny −=
-10 -5 0 5-20
0
20
40
60
Lag index
Amplitude
![Page 193: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/193.jpg)
Correlation Computation Correlation Computation Using MATLABUsing MATLAB
• The plot below shows the autocorrelation of
x[n] corrupted with an additive random noise generated using the function randn
• Note: The autocorrelation still exhibits a
peak at zero lag
-5 0 50
20
40
60
80
Lag index
Amplitude
![Page 194: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/194.jpg)
Correlation Computation Correlation Computation
Using MATLABUsing MATLAB
• The autocorrelation and the cross-
correlation can also be computed using the function xcorr
• However, the correlation sequences
generated using this function are the time-
reversed version of those generated using
Programs 2_7 and 2_8
![Page 195: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/195.jpg)
Normalized Forms of Normalized Forms of Correlation Correlation
• Normalized forms of autocorrelation and
cross-correlation are given by
• They are often used for convenience in
comparing and displaying
• Note: and
independent of the range of values of x[n]
and y[n]
][][
][][,
][
][][
000 yyxx
xyxy
xx
xxxx
rr
r
r
r ll
ll == ρρ
1≤|][| lxxρ 1≤|][| lxyρ
![Page 196: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/196.jpg)
Correlation Computation for Correlation Computation for
Power SignalsPower Signals
• The cross-correlation sequence for a pair of
power signals, x[n] and y[n], is defined as
• The autocorrelation sequence of a power
signal x[n] is given by
∑−=∞→
−+
=K
KnKxy nynx
Kr ][][lim][ ll
12
1
∑−=∞→
−+
=K
KnKxx nxnx
Kr ][][lim][ ll
12
1
![Page 197: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/197.jpg)
Correlation Computation for Correlation Computation for
Periodic SignalsPeriodic Signals
• The cross-correlation sequence for a pair of
periodic signals of period N, and ,
is defined as
• The autocorrelation sequence of a periodic
signal of period N is given by
][nx~ ][ny~
][nx~
∑ −= −= 10
1 NnNxy nynxr ][][][ ll
~ ~~~
∑ −= −= 10
1 NnNxx nxnxr ][][][ ll
~ ~~~
![Page 198: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/198.jpg)
Correlation Computation for Correlation Computation for
Periodic SignalsPeriodic Signals
• Note: Both and are also
periodic signals with a period N
• The periodicity property of the
autocorrelation sequence can be exploited to
determine the period of a periodic signal
that may have been corrupted by an additive
random disturbance
][lxyr~~ ][lxxr~~
![Page 199: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/199.jpg)
Correlation Computation for Correlation Computation for
Periodic SignalsPeriodic Signals
• Let be a periodic signal corrupted by
the random noise d[n] resulting in the signal
which is observed for where
][nx~
][][][ ndnxnw += ~
10 −≤≤ MnNM >>
![Page 200: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/200.jpg)
Correlation Computation for Correlation Computation for
Periodic SignalsPeriodic Signals
• The autocorrelation of w[n] is given by
∑ −= −= 10
1 ][][][MnMww nwnwr ll
∑ −= −+−+= 10
1 ])[][])([][(MnM
ndnxndnx ll~~
∑∑ −=
−= −+−= 1
011
01 ][][][][
MnM
MnM
ndndnxnx ll
∑∑ −=
−= −+−+ 1
011
01 ][][][][
MnM
MnM
nxndndnx ll
~ ~
~ ~
][][][][ llll dxxdddxx rrrr +++= ~~~~
![Page 201: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/201.jpg)
Correlation Computation for Correlation Computation for
Periodic SignalsPeriodic Signals
• In the last equation on the previous slide,
is a periodic sequence with a period N and
hence will have peaks at
with the same amplitudes as approachesM
• As and d[n] are not correlated, samples
of cross-correlation sequences and
are likely to be very small relative to the
amplitudes of
...,2,,0 NN=ll
][lxxr ~~
][nx~
][lxxr ~~
][lxdr~ ~ ][ldxr
![Page 202: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/202.jpg)
Correlation Computation for Correlation Computation for
Periodic SignalsPeriodic Signals• The autocorrelation of d[n] will show
a peak at = 0 with other samples having
rapidly decreasing amplitudes with
increasing values of
• Hence, peaks of for > 0 are
essentially due to the peaks of and can
be used to determine whether is a
periodic sequence and also its period N if
the peaks occur at periodic intervals
l
l
||l
][lwwr
][lxxr ~~
][nx~
][lddr
![Page 203: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/203.jpg)
Correlation Computation of a Correlation Computation of a
Periodic Signal Using MATLABPeriodic Signal Using MATLAB
• Example - We determine the period of the
sinusoidal sequence ,
corrupted by an additive
uniformly distributed random noise of
amplitude in the range
• Using Program 2_8 of text we arrive at the
plot of shown on the next slide
)25.0cos(][ nnx =950 ≤≤ n
]5.0,5.0[−
][lwwr
![Page 204: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/204.jpg)
Correlation Computation of a Correlation Computation of a Periodic Signal Using MATLABPeriodic Signal Using MATLAB
• As can be seen from the plot given above,
there is a strong peak at zero lag
• However, there are distinct peaks at lags that are multiples of 8 indicating the period of the sinusoidal sequence to be 8 as expected
-20 -10 0 10 20-60
-40
-20
0
20
40
60
Lag index
Amplitude
![Page 205: Discrete -Time Signals: Time -Domain Representation](https://reader030.vdocument.in/reader030/viewer/2022012409/616a4e0511a7b741a3510ae9/html5/thumbnails/205.jpg)
Correlation Computation of a Correlation Computation of a Periodic Signal Using MATLABPeriodic Signal Using MATLAB
• Figure below shows the plot of
• As can be seen shows a very strong
peak at only zero lag
-20 -10 0 10 20-2
0
2
4
6
8
Lag index
Amplitude
][lddr
][lddr