Digital Signal Processing
Digital Signal Processing
Digital Signal Processing
(DSP)
SIGNAL PROCESSING
Topics
• Continuous time (CT) and Discrete time (DT) signals
• Periodic and pulse signals
• Energy and power in signals
• Standard CT and DT signals
• Impulse, step, pulse, ramp, sine and exponential
• Analysis of CT signals
• Fourier Series and Fourier Transforms
• Convolution and Correlation
• Analysis of DT signals
• Discrete Fourier Transform of DT signals
• Fast Fourier Transform
• Z Transform
• Digital filtering in time domain
• Linear Filters
• FIR Filters
• IIR Filters
SIGNAL
• A signal is defined as any physical quantity that varies with time, space or any other
independent variable or variables.
• Mathematically, a signal can be represented as a function of one or more independent
variables.
S1(t) = 7t
S2(t) = 18t2.
Here S1(t) and S2(t) represents two signals- one that varies linearly with time t and the
other varies quadratically with t.
• Complex signals are there which cannot be expressed by simple mathematical
equations. (eg: speech signal, ECG , EEG…)
SIGNAL
Analog System
• Analog signals are continuous function of an independent variable, such as time, space
etc.
• It is defined for every instant of independent variable, so the magnitude is continuous
in the specified range for analog signals.
Limitations of DSP:
1. The conversion speed of ADC and the processing speed of signal processors
should be very high to perform all real time processing
2. Signals of high bandwidth require fast sampling rate ADCs and fast processors.
Applications of DSP:
n -2 -1 0 1 2 3
x(n) -1 2 1.5 -0.9 1.4 1.6
4. Sequence representation
0 n
δ1(n)
1;n≥0
Delayed impulse δ1(n) = δ(n – n0) = 1
0;n≠0
n0 n
0 n
The unit step signal is related to digital impulse by the
summation relation Unit step signal
n;n≥0
Ramp signal, ur(n) =
0;n<0
n
Ramp signal
4. Exponential signal
g(n)
a ;n≥0
n
0<a<1
Exponential signal, g(n) =
0;n<0
Exponential signal
Mathematical operations on Discrete Time signals
1. Shifting in time
A signal x(n) may be shifted in time by replacing the independent variable ‘n’ by
‘n-k’, where ‘k’ is an integer. If k is positive integer, the time shift results in a
delay by k units of time. If k is negative integer, the time shift results in an
advance of the signal by mod(k) units in time. The delay results in shifting each
sample of x(n) to right. The advance results in shifting each sample of x(n) to left.
x(n) 3
Eg:-
Let x(n) = 1 ; n = 2 2
1
2; n = 3
3; n = 4 0 1 2 3 4 5 6
Let x1(n) = x(n-2), where x1 is delayed signal
of x(n) x1(n) 3
2
when n = 4; x1(4) = x(4-2) = x(2) = 1 1
when n = 5; x1(5) = x(5-2) = x(3) = 2
when n = 6; x1(6) = x(6-2) = x(4) = 3 0 1 2 3 4 5 6
n
The folding of a signal x(n) is performed by changing the sign of the time base n
in the signal x(n). The folding operation produces a signal x(-n) which is mirror
of x(n) with respect to time origin n = 0.
Eg:- Let x(n) = n; -3 ≤ n ≤ 3. Now folded signal x1(n) = x(-n) = -n; -3 ≤ k ≤ 3
x 1(n)
x (n)
-3 -2 -1 1 2 3
0 1 2 3 n 0 n
-3 -2 -1
x(n) x1(n)
n n
0 1 2 3 4 5 0 1 2 3 4 5
5. Signal (vector) addition
The sum of two signals x1(n) and x2(n) is a signal y(n), whose value at any instant
is equal to the sum of the samples of these two signals at that instant.
i.e. , y(n) = x1(n) + x2(n) ; -∞ < n < ∞
The energy of a signal may be finite or infinite, and can be applied to complex
valued and real-valued signals. If E is finite (0 < E < ∞), then x(n) is called an
energy signal.
There are several common conventions for defining the Fourier transform of a complex-
valued Lebesgue integrable function, x In communications and signal processing, for
instance, it is often the function:
When the independent variable f represents time (with SI unit of seconds), the
transform variable represents ordinary frequency (in hertz). If is Hölder continuous,
then it can be reconstructed from X by the inverse transform:
The fourier transform of a finite energy discrete time signal, x(n) is defined as
The Fourier transform is one of the several mathematical tools that is useful in the
analysis and design of LTI (Linear Time Invariant) system. Another one is the Fourier
Series. These signal representations basically involve the decomposition of the signals in
terms of sinusoidal components. In such a decomposition, the signal is said to be represent
in the frequency domain.
To obtain the Fourier Transform representation, we shall start by finding the Fourier
transformation of sampled Analog signal.
Then ,
ax1(n) + bx2(n) ----------------- aX1(ejω) + bX2(ejω)
Proof:
F{ax1(n) + bx2(n)} = Σ(ax1(n) + bx2(n)) e-jωn
= Σ ax1(n) e-jωn + Σ bx2(n) e-jωn
= a X1(ejω) + b X2(ejω)
2. Time Shifting Property
Then ,
x(n-k) ----------------- e-jωk X(ejω)
Proof:
Let n-k = p
Then ,
x(-n) ----------------- X(e-jω)
Proof:
Let -n = p
F{x(p)} = X (e-jω)
(limit is from p=-infinity to + infinity). Hence proved.
4. Convolution Theorem
Proof:
Then,
Proof:
Results:
The fundamental result in LTI system theory is that any LTI system can be characterized
entirely by a single function called the system's impulse response. The output of the
system is simply the convolution of the input to the system with the system's impulse
response. This method of analysis is often called the time domain point-of-view. The same
result is true of discrete-time linear shift-invariant systems, in which signals are discrete-
time samples, and convolution is defined on sequences.
The Impulse response from a simple audio system. Showing the original impulse, with high
frequencies boosted, then with low frequencies boosted.
Practical applications of Impulse response
• Loudspeakers
A very useful real application that demonstrates this idea was the development of
impulse response loudspeaker testing in the 1980s which led to big improvements in
loudspeaker design. Loudspeakers suffer from phase inaccuracy, a defect unlike
normal measured properties like frequency response.
• Digital Filtering
Impulse response is a very important concept in the design of digital filters for
audio processing, because these differ from 'real' filters in often having a pre-echo,
which the ear is not accustomed to.
• Electronic processing
Impulse response analysis is a major facet of radar, ultrasound imaging, and many
areas of digital signal processing. An interesting example would be broadband
internet connections. Where once it was only possible to get 4 kHz speech signal
over a local telephone wire, or data at 300 bit/s using a modem, it is now
commonplace to pass 2 Mb/s over these same wires, largely because of 'adaptive
equalisation' which processes out the time smearing and echoes on the line.
FIR Systems
In FIR (Finite duration Impulse response) systems, the impulse response consists of
finite number of samples. The convolution formula for FIR system is given by,
From eqn 1 , it can be concluded that the impulse response selects only N samples of the
input signal.
Thus, the system acts as a window that views only the most recent N input
signal samples in forming the output. It neglects all prior input samples. So a FIR
system has a finite memory of length N samples.
IIR Systems
In IIR (Infinite duration Impulse Response) systems, the impulse response has infinite
number of samples. The convolution formula for IIR systems is given by,
Where the first summation is from k=1 to N and the second summation is over k=0 to M
Let Z{y(n)} = Y(z) ; Z{y(n-k)} = z-k Y(z)
Let Z{x(n)} = X(z) ; Z{x(n-k)} = z-k X(z)
The simplest analog IIR filter is an RC filter made up of a single resistor (R)
feeding into a node shared with a single capacitor (C). This filter has an exponential
impulse response characterized by an RC time constant.
The Z-transform, like many other integral transforms, can be defined as either a one-
sided or two-sided transform.
Bilateral Z-transform
The bilateral or two-sided Z-transform of a discrete-time signal x[n] is the function X(z)
defined as