Digital Signal Processing by S. Salivahanan - Text
Digital Signal Processing by S. Salivahanan - Text
SIGNAL
PROCESSING
lowers Signals and Systems
en S Salivahanan
=) A Vallavaraj
C Gnanapriya
Information contained in this work has been obtained by
Tata McGraw-Hill, from sources believed to be reliable.
However, neither Tata McGraw-Hill nor its authors
guarantee the accuracy or completeness of any information
published herein, and neither Tata McGraw-Hill nor its
authors shall be responsible for any errors, omissions, or
damages arising out of use of this information. This work
is published with the understanding that Tata McGraw-Hill
and its authors are supplying information but are not
attempting to render engineering or other professional
services. If such services are required, the assistance of an
appropriate professional should be sought.
Tata McGraw-Hill
ISBN 0-07-463996-X
Foreword v
Preface vii
Classification of Signals
and Systems
1.1 INTRODUCTION
Signals play a major role in our life. In general, a signal can be a
function of time, distance, position, temperature, pressure, etc., and it
represents some variable of interest associated with a system. For
example, in an electrical system the associated signals are electric
current and voltage. In a mechanical system, the associated signals may
be force, speed, torque, etc. In addition to these, some examples of
signals that we encounter in our daily life are speech, music, picture
and video signals. A signal can be represented in a number of ways.
Most of the signals that we come across are generated naturally.
However, there are some signals that are generated synthetically. In
general, a signal carries information, and the objective of signal
processing is to extract this information.
Signal processing is a method of extracting information from the
signal which in turn depends on the type of signal and the nature of
information it carries. Thus signal processing is concerned with
representing signals in mathematical terms and extracting the
information by carrying out algorithmic operations on the signal.
Mathematically, a signal can be represented in terms of basic functions
in the domain of the original independent variable or it can be
represented in terms of basic functions in a transformed domain.
Similarly, the information contained in the signal can also be extracted
either in the original domain or in the transformed domain.
A system may be defined as an integrated unit composed of diverse,
interacting structures to perform a desired task. The task may vary such
as filtering of noise in a communication receiver, detection of range of a
target in a radar system, or monitoring steam pressure in a boiler. The
function of a system is to process a given input sequence to generate an
output sequence.
2 Digital Signal Processing
= t
0
(a) Continuous-time signal
xr (n)
E UE 2T
EE n oe
(b) Discrete-time signal
Fig. 1.1 Continuous-Time and Discrete-Time Signals
4 Digital Signal Processing
(a) Continuous-time
x(n)
Where
-n, -2sn<0
x[nj=4 n Osn<2
0, otherwise
No=4 and T=0.25
(b) Discrete-time
Solution
(a) x,(t) = sin 15 rt is periodic.
ee P 2n_ 27
The fundamental period is T, = eT
—— = —
= 0.1333333333... seconds
(b) x(t) = sin 20 nt is periodic.
2n_ 20
The fundamental period is T, = w ae 0.1 seconds
Classification of Signals and Systems 7
= 1.41421356... seconds
(d) x4(t) = sin 5rt is periodic.
The fundamental period is T, = =— = 0.4 seconds
Sinc (x)
Fig. 1.3 Typical Examples for (a) Odd Signal and (b) Even Signal
(c) the Sinc (x) Function
The energy signal is one which has finite energy and zero average
power, i.e. x(t) is an energy signal if 0 < E < œ, and P = 0. The power
signal is one which has finite average power and infinite energy, i.e.
0< P < œ, and E = æ. If the signal does not satisfy any of these two
conditions, then it is neither an energy nor a power signal.
J Sode=1 (1.18)
The Equations 1.17 and 1.18 indicate that the area of the impulse
function is unity and this area is confined to an infinitesimal interval on
the t-axis and concentrated at t = 0. The unit impulse function is very
useful in continuous-time system analysis. It is used to generate the
system response providing fundamental information about the system
characteristics. In discrete-time domain, the unit-impulse signal is
called a unit-sample signal. It is defined as
1 n=0
8(n)(n) = {
6, n20 (1.19) J
= fore
Classification of Signals and Systems U1
That is,
24
t, t>0 “ee
A ramp signal starts at t = 0 and increases linearly with time, t. In
discrete-time domain, the unit-ramp signal is defined as
M(t)t) = u(t+2) uf 2)
1
ult+—)-uft-= 1
(1.26)
1.26
The signals u(t+ 3)and u(t- 3)are the unit-step signals shifted by i
units in the time axis towards the left and right, respectively.
Figure 1.4 shows some of the singularity functions. The advantage of
the singularity function is that any arbitrary signal that is made up of
straight line segments can be represented in terms of step and ramp
functions.
1.3.5 Properties of ô (t)
1. f 5) dt=1
4. | MS-A dà = xlt)
5. Slat) = + 8t)
lal
6. x(t) (t — to) = x(ty)
A x(to) clt - to) = x(to)
t2
8. J x(t) &" (t — to) dt = (-1)" x(t)
ti
12 Digital Signal Processing
öt) 5(n)
:
3
2
1
— ke t -3 -2 -1 0 1 2 3 4 an
í.
A
(a) Unit-impulse function
r(t)
Oo 4 2 3 t -3 -2 -1 0 1 2 3 4n
nie
(c) Unit-ramp function
oT -05 0
(d) knee
Fig. 1.4 Singularity Functions (a) Unit-Impulse Function (b) Unit-Step Function
(c) Unit-Ramp Function (d) Unit-Pulse Function
Classification of Signals and Systems 13
Proof
[ x ( t )
8(¢- t)]? = jx(t) 5(t — to) dt + å (to)
ti
LHS = 0.
t
x(t)
t t
-3⁄2 0 -144 01/4 3/4
Fig. E1.2 (a) Fig. El.2 (b)
(b) x(t) = 2M(¢ — 1/4)
Here the signal shown in Fig. E1.2(b) is shifted to the right, with
centre at 1/4. Since a = 1, the signal width is 1 and amplitude is 2.
(c) x(t) = cos(20 nt- 52)
Here the signal x(t) shown in Fig. E1.2(c) is shifted by quarter cycle
to the right.
Fig. El.2(c)
(d) x(t) =r(—0.5¢ + 2) x
Fig. E1.3
Solution
Representation through addition of two unit step functions
The signal x (t) can be obtained by adding both the pulses, i.e.
x(t) = 2[u(t) — u (t — 2)]+[u(t — 3) - u (t - 5))
Representation through multiplication of two unit step functions
x(t) = 2[u (t) u(-t + 2)) + [u (t — 3) u(t + 5))
= 2(u (t) u(2 — t) + u (t — 3) u(5
- t))
ei(2aft+o)
x(t) = Se
2
From this the amplitude spectrum for the signal x(t) consists of two
components of amplitude, viz. A/2 at frequency ‘f’ and A/2 at frequency
‘f’. Similarly, the phase spectrum also consists of two phase
components one at ‘f’ and the other at ‘-/’. The frequency spectrum of
the signal, in this case, is called a double-sided spectrum. The following
example illustrates the single-sided and double-sided frequency spectra
of a signal.
$3 |
q
i
g 10
Amplitude
]
‘J|
$ 2 +
|
1=
-10 te)
Fig. E.1.4 Amplitude and Phase Spectra (a) Single-Sided and (b) Double-Sided
Classification of Signals and Systems \7
x A
Hlaz
(t)+azx(t)]
x(t)
x(t)
aHix(t)] + a2gH(x20)]
x2(h
Solution Let the response of the system to x,(t) be y,(t) and the
response of the system to x(t) be y(t) . Thus, for the input x,(¢), the
describing equation is
dy, (t)
“= + 2y; (t) = x,(t)
i.e.
oth
(9 +y (t)+4= x(t)
and for input x, (t),
be + Yq (t) + 4 = x3 (t)
Shifting
In the case of discrete-time signals, the independent variable is the time,
n. A signal x(n) may be shifted in time, i.e. the signal can be either
advanced in the time axis or delayed in the time axis. The shifted signal
is represented by x(n — k), where k is an integer. If ‘k’ is positive, the
signal is delayed by k units of time and if k is negative, the time shift
results in an advance of signal by k units of time. However, advancing
the signal in the time axis is not possible always. If the signal is available
in a magnetic disk or other storage units, then the signal can be delayed
or advanced as one wishes. But in real time, advancing a signal is not
possible since such an operation involves samples that have not been
generated. As a result, in real-time signal processing applications, the
operation of advancing the time base of the signal is physically
unrealizable.
Folding
This operation is done by replacing the independent variable n by —n.
This results in folding of the signal about the origin, i.e. n = 0. Folding is
also known as the reflection of the signal about the time origin n = 0.
Folding of a signal is done while convoluting the signal with another.
Time scaling
This involves replacing the independent variable n by kn, where k is an
integer. This process is also called as down sampling. If x(n) is the
discrete-time signal obtained by sampling the analog signal, x(t), then
x(n) = x(nT), where T is the sampling period. If time-scaling is done,
then the time-scaled signal, yin] =x (kn) =x(knT). This implies that the
sampling rate is changed from 1/T to 1/kT. This decreases the sampling
rate by a factor of k. Down-sampling operations are discussed in detail
in Chapter 11 of this book. The folding and time scaling operations are
shown in Fig. 1.7(a) and (b).
22 Digital Signal Processing
-8 -6 -4 -2 (0) 2 4 6 8
(a) Unit-step signal
u(n- 4)
TA . = . = .—ei
u(n+ 4)
eer =8 -4
TE LLTI
-2 0 2 4 6 8
L,
(c) Advancing
the unit-step signal by 4 units
x(n)
Ma HOOKS S81 1.08
tet,
-7 -6 -6 -4 -3 -2
-7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7 8H
(a) Folding
x(n)
x(n) ={1, 1,2,1, = 2,0, 1, 3, 1}
x(2=n{1), 2,1,2, 1, 1}
a a
eee
-7 -6 -5 -4 -3 -2 -1 0
|
1 N a a a a s aol ay
(a) Continuous-time
: lon | y(n)
=Hy(x(n))
Jes
—
y(n)
— T7
ys(n) = Ha[ya(n)]
yaln) =Haly(n)]
(b) Discrete-time
twice the input delayed twice, x(n — 2). Let the input sequence be x(n) =
{0,1, 1, 2, 0, 0, 0, ...}. The output sequence for the system as described by
Eq. 1.32 is y(n) = {0, 1, 4, 7, 8, 4, 0, 0, ...}. The block diagram
representation of the system described by Eq.1.32 is shown in Fig. 1.9.
x(n) y(n)
Asn
time). If y(t) is the system response for an input x(t), then the response
of the system when x(t) = &(t) is y(t) = A(t).
The impulse response of a system can be directly obtained from the
solution of the differential or difference equation characterising the
system. The impulse response is also determined by finding out the
output of the system to the rectangular pulse input x(¢) = zIl ($) and
then taking the limit of the resulting system response, y(t) as € > 0. The
unit-impulse function is nothing but the derivative of the unit-step
signal. Therefore, the impulse response of the system can also be
obtained by computing the derivative of the step response of the system.
1.7.3 State-Variable Technique
The state-variable technique provides a convenient formulation
procedure for modelling a multi-input, multi-output system. This
technique also facilitates the determination of the internal behaviour of
the system very easily. The state of a system at time tọ is the minimum
information necessary to completely specify the condition of the system
at time fy and it allows determination of the system outputs at any time
t > to, when inputs upto time ¢ are specified. The state of a system at
time fy is a set of values, at time tọ, of a set of variables. These variables
are called the state variables. The number of state variables is equal to
the order of the system. The state variables are chosen such that they
correspond to physically measurable quantities. It is also convenient to
consider an n-dimensional space in which each coordinate is defined by
one of the state variables
x,, Xo, ...,X,, where n is the order of the system.
This n-dimensional space is called the state space. The state vector is
an n-vector x whose elements are the state variables. The state vector
defines a point in the state space at any time t. As the time changes, the
system state changes and a set of points, which is nothing but the locus
of the tip of the state vector as time progresses, is called a trajectory of
the system.
A linear system of order n with m inputs and k outputs can be
represented by n first-order differential equations and & output
equations as shown below.
dx,
ry = 04, X1 + AygQXqt ... + AyyXp_ + Oy) Uy + Ogu g+... + bim Um
dxs
E = 1 Xy + Ogg Xo + ... + Aon Xn + O91 Uy + Ogg g+... + DamUm
` (1.84)
we
d a= an1 ž1 + Ang Xz+... t brougt
+ Ann Xn + b On)Uy+b ... + Onm Um
+...+6,,,u
Classification of Signals and Systems 27
and
Yı = Cy Xy + Cig Xot... + Cy_ Xp, + yy Uy t+diglot ... + dy, Up
Yq = C21Xy + Cog Xo+ ... + ConXn+ dgy Uy + doggt ... + dom Um
. (1.35)
Yp = CyyXy
+Cho Xot ... + Can Xn +tAyy Uy + Ayg llgt ... + dpmUm
where u;, i = 1, 2, ..., m are the system inputs, x;, i = 1, 2,3, ..., n are
called the state variables and y; i = 1, 2, 3, ..., k are the system outputs.
Equations 1.34 are called the state equations, and Eqs 1.35 are the
output equations. Equations 1.34 and 1.35 together constitute the state-
equation model of the system. Generally, the a’s, b’s, c’s and d’s may be
functions of time. The solution of such a set of time-varying state
equations is very difficult. If the system is assumed to be time-invariant,
then the solution of the state equations can be obtained without much
difficulty.
The state variable representation of a system offers a number of
advantages. The most obvious advantage of this representation is that
multiple-input, multiple-output systems can be easily represented and
analysed. The model is in the time-domain, and one can obtain the
simulation diagram for the equations directly. This is of much use when
computer simulation methods are used to analyse the system. Also, a
compact matrix notation can be used for the state model and using the
laws of linear algebra the state equations can be very easily
manipulated. For example, Eqs 1.34 and 1.35 expressed in a compact
matrix form is shown below. Let us define vectors
xy uy yı
x u
x= tg , L= $ > ya sa (1.36)
Xn Um Yk
and matrices
Oy Gg - + + Om by biz ~~ + bim
G2, An > - + Ban ba ba . . . bom
A=]. Pia sae Sh. BS s E r -& (1.37)
Cu Ci Cin dı d dim
C21 C22 Con dı dz dom
Ce , D=
Ls|Sampler jra
cee BEBie eal
Continuous-time Discretoiüme Discrete-time poni output
continuous-amplitude Continuous-amplitude discrete-amplitude
input signal signal signal
Samples of x(t)
x(t)
0 T 2T 3T 4T 5T 6er t
(a) Samples of x (4)
|Switch (9 ai
x(t) x (i) e
= a g(t)
(b) Modelling a sampler as a switch (c) Mode! of a sampler
>t H-
| |
a ee khae A eee ki eer
T 2T 3T 4T 5T 6T
where
r
C= goetan: (1.42)
B
2
nÒ FC = Eoen (143)
na-am Nao
Thus,
xh
gi= E C,e"r4
i.
where
i"
Gal face etar (1.50)
T 7
g()
1
A
t
-6T -5T -4T-3T -2T-T 0 T 27 3T 4T ST 6T
(a)
Ax)
Fig. 1.14 (a) Impulse Sampling Function (b) Spectrum of the Signal x(t)
(c) Spectrum of Impulse Sampled Signal
Cc n =7E
eet o_1_ "T fe (1.52)
5
Thus C, is same as the sampling frequency f,, for all n. The spectrum
of the impulse sampled signal, x,(t) is given by
The spectra of the signal x(t) and the impulse sampled signal X, (t)
are shown in Figs 1.14 (b) and (c). The effect of impulse sampling is
same as sampling with a train of pulses. However, all the frequency
translated spectra have the same amplitude. The original signal X(f
can be reconstructed from X,(f) using a low-pass filter. Figure 1.15
shows the effect of sampling at a rate lower than the Nyquist rate.
Consider a bandlimited signal x(t), with f, as its highest frequency
content, being sampled at a rate lower than the Nyquist rate, i.e.,
sampling frequency f, < 2f}. This results in overlapping of adjacent
34 Digital Signal Processing
x(f)
-h 0 fn
Xan)
f
-6 -h O h fs-ħfhs fth
(b) Spectrum
of the sampled signal's tor fs > 2fp
>» f
of the low-pass filter will bef, = a. Therefore, the unit impulse response
of an ideal filter for this bandwidth is
h12
ht)=T fe df (1.54)
~f,/2
That is
h(t)
=jax"
—P_ (eit lt _ gg
T
git jx fat -jR fat
xt)= Y xnT)ht-nT)
Using Eq.1.55, we get
(b)
Fig. 1.16 Signal Reconstruction (a) Reconstruction Filter
(b) Time Domain Representation
Quantisation level
[0] E 2T 3T 4T
Fig. 1.17 Quantizing and Encoding
REVIEW QUESTIONS
1.1 What are the major classifications of signals?
1.2 With suitable examples distinguish a deterministic signal from
a random signal.
1.3 What are periodic signals? Give examples.
1.4 Describe the procedure used to determine whether the sum of
two periodic signals is periodic or not.
1.5 Determine which of the following signals are periodic and
determine the fundamental period.
(a) x(t) = 10 sin 25 nt (b) x(t) = 10 sin V5 xt
38 Digital Signal Processing
< t< o0
(c) x(t) = 100 sin (107: - zz) +50 cos{25 xt - z), —20
How are systems classified?
Distinguish static systems from dynamic systems.
What is linear system ?
Determine whether the following systems are linear
(a) we n oun
+ 5y(t) + 2 = x(t) b) 5 —— + y(t) = 5x(t)
lo) ao
o+ y(t) + 5 = 10x(t)
1,22 What isym system?
1.23 What is a causal system? Why are non-causal systems
unrealisable?
1.24 What is BIBO stability?
Classification of Signals and Systems 39
2.1 INTRODUCTION
A signal which is repetitive is a periodic function of time. Any periodic
function of time f(t) can be represented by an infinite series called the
Fourier Series. A function of time f(t) is said to be periodic of period T
if f(t) = f(t + T) for all ¢. For example, the periodic waveforms of
sinusoidal and exponential forms are shown in Fig. 2.1.
4 f(t)
0 Ti2 T t
{c)
Fig. 2.1 Waveforms Representing Periodic Functions
Fourier Analysis of Periodic and Aperiodic Continuous-Time Signals and Systems 41
£ + a3 cos 3t + ...
fH = Jao+ a, cos ot+ a, COS 2M
+ b, Sin Wp t+ by sin 20t + bz sin 30t + ... (2.1)
where @) = 2nf= z, f is the frequency and a’s and b’s are the
coefficients. The Fourier series exists only when the function ft)
satisfies the following three conditions called Dirichlets conditions.
(i) f(t) is well defined and single-valued, except possibly at a finite
number of points, i.e.
f (t) has a finite average value over the period T.
(ii) f(t) must posses only a finite number of discontinuities in the
period T.
(iii) f(t) must have a finite number of positive and negative maxima in
the period T.
Equation 2.1 may be expressed by the Fourier series
Therefore, ff(t) dt = i aT
-T/2 2
T12
Hence, ap= > ffOdt (2.3)
FT -Ti2
42 Digital Signal Processing
T
or, equivalently ag = żj f(t)dt
0
Multiplying both sides of Eq. 2.2 by cos m Wot and integrating, we have
T/2 T/2
J Feos moot dt= > Jao cos moot dt+
-T12 -T12
TIZ 9 T a
J $a, cos nwo t cos moot dt + Yb, sin nwo t cos m wot dt
-Tign=1 -Tign=1
172
Here, — fÍdy cos M Wot dt
=0
-T/2
T/2 a, 72
Jaa, cos n wot cos mog dt= —* J[cos (m + nwo t + cos(m ~ n) wot] dt
-T/2 2 te
|0, form#n
—a,, form=n
2
T/2 T/2
fon sin N@ot cos mogtdi = m fisin (m + n) Wot - sin (m — n)@ot]dt
-T12 -T/2
=0
T12 Ta
Therefore, fro COS NWot dt = —*, form=n
-T/2
T12
Hence, an = 2 fro COS NWy t dt (2.4)
TT
T
or, equivalently a, = Zjf(t) cos n og t dt
0
T w TI2 6
+ Í J a, cos n Wot sin m Wo t dt + Í F bn sin n wot sin m wp t dt
-Tign=1 -Ti2n=1
Fourier Analysis of Periodic and Aperiodic Continuous-Time Signals and Systems 43
T/2
Here, > Í dy sin MWg
£dt = 0
-T12
T12
Í a, COS N Wọ É sin MWg tdt = 0
-T12
T12 0, form#n
¢ sin mota=lT,
fon sin nọ ieman
-T/2 gn
T12
Therefore, fro sin MWg tdt = Zon form=n
-T12
T/2
Hence, b, = 2 JF@)sin nootdt (2.5)
-T12
T
or, equivalently, b„ = Zj f(t)sin n Wotdt
0
The number n = 1, 2, 3, ... gives the values of the harmonic
frequencies.
Symmetry Conditions
(i) If the function f(t) is even, then f-t) = At). For example, cos t, t,
t sin t, are all even. The cosine is an even function, since it may be
expressed as the power series
cos t= 1 — —
frod = affit)at
A 4
f(t) f(t)
(a) (b)
f(t) f(t)
(c) (d)
Fig. 2.2 Waveforms Representing Even Functions
of an odd function is symmetrical about the origin. Iff (t) is odd, ff(t) dt
= 0. The sum of two or more odd functions is an odd function and the
product of two odd functions is an even function.
(a) (b)
fit)
A
Fig. E2.1
Solution Since the given waveform is symmetrical about the
horizontal axis, the average area is zero and hence the d.c. term
a, = 0. In addition, f(t) = f(-t) and so only cosine terms are present,
i.e., b, = 0.
T12
Now, a= 2 fro cos n Wot dt
-T/2
~A, from -T/2<t<-T/4
where f(t)=4+A, from -T/4<t<+T/4
-A, from +T7/4<t<+T/2
Therefore ,
2a T T/4 T/2
a, = Al fc COS NWot)dt + Joos NWotdt + J(-cos NW, t) a
~T/2 -T74 T/4
i -T14 i T14 y T/2
_2A [2] (Saree e [=ne] |
T noo Jar noo Jra no Jra
-24 -sin (= $27) + sin (2287), sin (2907)
NO oT 4 2 4
` jeer : (7) : (**)]
-sin |——— |-sin + sin |——
4 2 4
8A (=27) 4A (752)
= sin |—— |- sin |——
nT 4 NWT 2
When wọ T = 2n, the second term is zero for all integer values of n.
Hence,
8A . (=) 4A. (=)
a, = —— sin |— |= — sin |—
2nn 2 nt 2
ay = 0 (d.c. term)
46 Digital Signal Processing
ay = tA sin(n)=0
4A in(32) 4A
|
A
~ 72 - 7/4 0 7/4 T2
Fig. E2.2
Solution The given waveform for one period can be written as
0, for -T/2<t<-T/4
f(t)= 4A, for -T/4<t<T/4
0, for T/4<t<T/2
For the given waveform, f (—t) = f (t) and hence it is an even function
and has b, = 0.
The value of the d.c.term is
2 T/2
a, = — J £© cos n wot dt
-T/2
Fourier Analvsis of Periodic and Aperiodic Continuous-Time Signals and Systems 47
7 z nof
T -T74
4A
= noT sin (n0 T/4)
2
=0, forn=2, 4,6,
F 372
Fig. E2.3
Solution As the waveform shows no symmetry, the series may
contain both sine and cosine terms. Here, f(t) = A sin Wot
To evaluate ap:
T
=afA sin Wot dt
0
T/2
a [A sin Wot dt
0
48 Digital Signal Processing
art 2A 2A
cos wt]? =ort cos (9T'/2) + 1]
2
a= žfro COS Notdt
at
p 4i Wo
£ cos n Wotdt
a, = g leona N
n(l-n
Age
==r en 2@ot dt
A T72
= T 7l cos 209 t],
2
6, = 2 sinNO, tdt
T/2
=7 JAsin Wy tsin NW,
t dt
0
A à T12
2A | nsin Wgt cos NW, t — sin n Og £ COS Wp t í
T -n? +1 b
Fourier Analysis of Periodic and Aperiodic Continuous-Time Signals and Systems 49
_ 2A [st- sinzeit)"
2 a4
When 0) T = 21, we have b, = A,
Substituting the values of the coefficients in Eq. 2.2, we get
Fig. E2.4
Solution
(i) As the waveform has equal positive and negative area in one
cycle, the average value of aç = 0.
(ii) As f(t) = -f (t), it is an odd function and hence a, = 0
472
and b, = T [iain NWotdt
(iii) Here, f(t + 7/2) = -f (t). Hence it has half-wave odd symmetry
and a, = b, = 0 for n even.
(iv) To find f(t) for the given waveform
ft)-0 _ 0-A
t-0 0-T/4
4A
Therefore,
erefore, f(t) =—T
Pare T
7 <t< >
theregion —<t<—,
Foror th
fit)-A_ A-0
joe FT
4 4 2
4A T 4A
fit)-Az= aC *)- atta
T/4 T/2
_4 4A). 4 4A :
=T J(“A)esinnwg tdt + FJ (Fit 24) sinn oy a
+
8A [= naar)”
—— |M
T =n Oo Tl4
Substituting wp = =. we have
b3 =
8A .
sin
3n =
8A
T- T T+d/2
Fig. E2.5
Solution The periodic function of the Fourier series for the given
pulse train is expressed by
d/2
2Ad
=2 [ae 24
22 giz, = T
T an
2A/snneat)
= 22 2A [si(2904) : (=24)]
= sin -sin |———
T NW Jan NOT 2 2
4A . n@d
—————
81 Dr
NW 2
Ad , 2Ad © sin (nw, d/2)
Hence, f(t)= > +5 »y“aa oe
n=l
ej”
eot _ g7inoot
sin
n Mt =
2j
Substituting these quantities in the expression for the Fourier series
gives
jnwot — JN wot oo Jn@ot _ pinot
po= 4a, + È o,(oe $ [eee]
n=l n=1
172
= [roete at (2.7)
T tye
172
and Cn? 5 Jro [cos n Wp t +j Sin n Wy t] dt
T -T/2
T/2
=2 frien! at (2.8)
T -T/2
- r =| i
with f= cot $ cepet Y cp ef"! (2.9)
nal n=-%
where the values of n are negative in the last term and are included
under the E sign. Also, co may be included under the = sign by using the
value of n = 0. Therefore,
Bg
(a) Find the trigonometric Fourier series of the waveform shown in
Fig. E2.6 and
(b) Determine the exponential Fourier series and hence find a, and b,,
of the trigonometric series and compare the results.
Fig. E2.6
Solution The function of the given waveform for one period can be
written as
fit)= a for -T/2 <t<0O
+A, for 0 < t < T/2
As the waveform is symmetrical about the origin, the function of the
waveform is odd and hence ay = a,, = 0, and
T/2
bes f f (t) sin nœtdt
-772
2 0 T12
= 7- cos
(2 Wy T/2))
noT
When œ = ="
3 | we have
2A
=f
re g cos nt]
0 , ifniseven
4A
bn =) 44 |itn isodd
nn
Here cp = |
$a =0
To evaluate c,,
Since the wave is odd, c, consists of pure imaginary coefficients. From
Eq. 2.7, we have
c, = 17
al f(t)e j2"! dt
0
1 0 ; T/2 f
-+ Í (- A) e270! dt + Í Ae Jnot J
T -T/2 0
o T/2
i A |-1) 1 ei af 1 e inet
E (- jno) -rz L(- Jn ) 0
„A,
s2. Lh +e Jno (712) en ina (T/2)
-
0
T (- cat i }
When œ = ZE we get
Therefore,c, = -j (24) te
for odd n only.
56 Digital Signal Processing
a
Example 27
(a) Find the trigonometric Fourier series of the waveform shown in
Fig. E2.7 and
(b) Determine the exponential Fourier series and hence finda, and b,
of the trigonometric series and compare the results.
A A(t)
54
@gt
0 2n án °
Fig. E2.7
Solution (a) As the waveform is periodic with period 27 in œ t and
continuous for 0 < @t < 2n, with discontinuities at Mot = n (27),
where n = 0, 1, 2, ..., the Dirichlet conditions are satisfied.
To find f(t) for the given waveform of region 0 < Opt < 2n:
pee;aes es S
The equation of the straight line is———
~% %,—%
= (2n, 5), we get
Substituting (x, yı) =(0, 0) and (x, girs
ft)-0 _ 0-5
@)t-0 0-27
T/2
sales
f fe ae
-T12
2r
2 5
On JOn Oo t d (wg t)
a t27?
“@ =
__10 (21)? _ 5
(Qn)? 2
T/2
Using Eq. 2.4, we obtain a, = Z Í f(t) cos nOg t dt
-T12
2z
2 5
=— — |@ t cos noy t d (0y £)
Qn J(4) i g $
1 —
(5
==z tac t sin
si n W td
td (Wo t t)
2r
W o t y t ; non]
= z| - cos n a + -Fsin
2 n n 0
ss%
nt
Combining the average term and the sine-term coefficients, the series
becomes
-Ž
= sin oot-
lo > sin 2a ¢ - > sin 3a ¢ -..
a
58 Digital Signal Processing
_5 < arot
rola r
(b) To determine exponential Fourier series
Here, co = la
3 2
To evaluate c,,:
From Eq. 2.7, we have
i?
7l fde d
= -jnwot t
2n
Cn l | (zante ac t)
“On
5 fein ir 5
= (-jnea,t-v| | =j
ean
(Qn)? Fax ee
Substituting the coefficients c, in Eq. 2.9, the exponential Fourier
series is
, 6 ~j2ogt . 6 ~j@ot 5 . 5 jOgt
(H) =... j —— eitt j eltt g Hy jel
f 4x Jon 2 7 on
. 5 -j2opt
+J — P ias
TIn
By using Eq. 2.11, the trigonometric Fourier series coefficients a, and
6, can be evaluated as
2 (ag/2) ey
Therefore, ae | [f(y dt= T f (f@jae
-T12 -T12
Fourier Analysis of Periodic and Aperiodic Continuous-Time Signals and Systems 59
aE
te
= | f(t)dt
2 T/2
a,= J f(t) cos ng t dt
-T/2
T/2
bn 2 f f(t) sin ngtdt
T -T12
Therefore, substituting all these values, we get
T/2 à
1a I (peor? 2 aedt= =(S2.)
2) + 1+2 È (ito)2 (249
This is the Parseval’s identity.
P== pf iol 2 dt
The Fourier series for the signal f(t) is
f(t) = 5 Cp einot
n= -w
= ob +S cf +3c5
+... (2.15)
| Example 2.8] The complex exponential Fourier representation of a
signal f(t) over the interval (0, T ) is
= — 3 jnnt
fo p 4+ (nn)? :
(a) What is the numerical value of T ?
(b) One of the components f(t) is A cos3nt. Determine the value of A.
(c) Determine the minimum number of terms which must be retained
in the representation off (t) in order to include 99.9% of the energy
in the interval.
Fourier Analysis of Periodic and Aperiodic Continuous-Time Signals and Systems 61
f= 2 Tran
= < 3 jnrt
2n
Hence, T
— m =m ie
ie. T=2
ei Oe
4+(3n)?
a < 3
(c) Total (maximum) erP, = ——, | = 0.669
å ponen E 4+ (nn)?
The power in f(t) is
_=|=]13? +2 3 m
F Md 3
PEET
F A
a.
PET)
f 3a
alPa)
F
4 4+(n) 4+(2n) 4+ (37) 4+(4n)
62 Digital Signal Processing
f= Fc, e/ 7%!
T/2 l
where c, = VT f f(t) ereot at
-T/2
In the limit, for a single pulse, we have
T > =~, @ = 2n/T > dw (a small quantity)
or VT = W/2n > dw/2n
Furthermore, the n” harmonic in the Fourier series is n @ > nda.
Here n must tend to infinity as @) approaches zero, so that the product
is finite, i.e. n Wy > @.
In the limit, the £ sign leads to an integral and we have
_ do J
ae 7 fie -jot dt
a er do
f(t) = mad F(ja) e/
or, equivalently,
E= | fede
E==fJ ros,
L ftJ Fue)
riie)elmt
e/% dw dt
s>
= 1 n ruw]
; T
fOe jot
afao
= — | F(jo) F(-jo)do
is
“35 |FUer® (jo) do
64 Digital Signal Processing
ar e
oma! F(jo)|° do
Operation fW
Time-Scaling f (at)
dF(j@)
Frequency-differentiation (— jt) f (t) da
(Contd.)
Fourier Analysis of Periodic and Aperiodic Continuous-Time Signals and Systems 65
Operation
t
J A@he-vde
Frequency convolution
2.7.1 Linearity
The Fourier transform is a linear operation. Therefore, if
f(t) < Fy (j o)
h(t) e Fo (j o)
then, af; (t) + bf, (t) = aF, (j œ) + bF; (j w)
where a and b are arbitrary constants.
2.7.2 Symmetry
If f(t) e F (jo)
then, F( jt) <= 2nf (—)
Proof
F [f (at)) = f flat) es dt
fat) ka F(Z2)
ja} Xa
We conclude that larger the duration of the time function, smaller is
the bandwidth of its spectrum by the same scaling factor. Conversely,
smaller the duration of the time function, larger is the bandwidth of its
spectrum. This scaling property provides an inverse relationship
between time-duration and bandwidth of a signal i.e. the time-
bandwidth product of an energy signal is a constant.
2.7.4 Convolution
Convolution is a powerful way of characterising the input-output
relationship of time-invariant linear systems. There are two convolution
Fourier Analysis of Periodic and Aperiodic Continuous-Time Signals and Systems 67
theorems, one for the time domain and another for the frequency
domain.
Time Convolution
Proof
= jHof]h(t
-t)e1” alas
=X jo) H( jo)
Hence, the convolution of the signals in the time domain is equal to
the multiplication of their individual Fourier transforms in the
frequency domain.
x(t) = etPC
Fig. E2.9
68 Digital Signal Processing
_aa1
-+—
10 RC
Similarly, the transfer function of the network is
H(jo)= _VjoC _ a
(z + een
1 ) (j@RC +1)
JoC
BSE N 1
~ RC (Jo+ 5)
Re
3 ; ‘ 1 1
Hence, Y¥(j@) = X(j@) H( jo) = RGT 1
(e+e)
ther = 1 [Y n = —1 te _RC
yt)=7 7 [YCj o) RG u(t)
joc
The given input time function is x(t) = te~”®
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
70 Digital Signal Processing
1
Asja +j) ljo- -a= Ga
H . | Peat: a
(b-a) lz+jo3 )e i (b+ 1 ;
3 Yi
ence, Y(jo) = ———
jal
Taking inverse Fourier transform, we get
y(t)=gigle
phen | - D- -bt uw]
When b = a, the partial fraction expansion is invalid. Hence,
z 1
Y( jo) = ——
j (a+ jo)?
ig A 1
2/2]
Using dual of the differentiation property,
etulthe 1
+jo
- .d 1 1
te“ ult) j— ear
: iil (a+ jo)?
Therefore, y(t) = te~™ u(t)
2.7.5 Frequency Convolution
If f(t) <= F( j@) and g(t) @ Gj),
then f(t) g(t) = Z Foe G(jo)
Proof
The inverse transform of [F( jœ) * G( j @)]/2n is
2 = -
p-[Fue) on _(1
*«G(j@) |- t : ere
(4) le J FGwade ju) du dw
1
J FOW |GGo- ju) ei” dodu
(27)? a 00
Putting x = o — u, then
o = x + u and dx = do
Therefore
. . 2 œ Gd
2n
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
74 Digital Signal Processing
F(jo) =2
1 e-20 w+ D
+1
Using the properties of the Fourier transform, write the Fourier
transforms of
t
eft-
l
+D
CPIu09)
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
82 Digital Signal Processing
(t) A
F(jo)=F(f@1 = | fe? at
T/2 ,
Í 1-e°/°! dt
-T/2
1 2, T/2
em dad
1 [e7402 - ej9712])
-jo
Comments
(1) The phase of the amplitude spectrum is exactly the same as that
in the previous case given in Fig. 2.7. :
(2) There is an additional uniform phase shift factor e/°”? which
changes the phase spectrum of the previous case.
obtained.
10
F(jo)= | fOe at
2 f -jot 2
=f 10677” dt = 10)£ |
0 “jo n
A 10(/ -jo
= 20e /° sine ©
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
90 Digital Signal Processing
N
f
EIS
/
/
P,/ \ Ps ase
- Ti2 o 1/2
Fig. 2.10 Triangular Pulse
Equation of line P. A is
2
f(t)p,p, = a
——t+A= A(1+ T t)
f tt PaPa =-4t+A=A(1-
T/2
A
21)
T
2
2 T
Therefore f(t)=A l+at for ~5<tso
Now,
T/2
F (jo)= f reese" dt = | f(t)e4 de
EA -T12
T/2
= ffiee dt + [|oase dt
-T/2
0 T72
f A(1+Ztlevmars f A(1-2t)e
-T/2 T 0 T
T42
=A fel dt+A Jeiet dt
-T?2
0 Ti2
2a jista a- 24 J tet” de
-T72 0
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
94 Digital Signal Processing
Therefore,
xf 2
F [e] = vx ete)
a
F( jo) = FÐ = f Se de=1
Hence, we have the pair 5(¢) = 1.
The frequency spectrum of the impulse function 6(t) shown in
Fig. 2.13 (a) has a constant amplitude and extends over positive and
negative frequencies.
4 A(t) F( ja)
8(t)
— f en ne
te) 0
Fig. 2.13 (a) Impulse Function and its Spectrum
Slope =-b
> t — — > ©
0 b 0 0 g
0 oo
- [eze] [=e]
© L-jo J. L -je Jo
1,1 _2
>
o
= = L Sa E
0 o
(a) (b)
Fig. 2.18 (a) Amplitude and (b) Phase Spectra of the Signum Function
Therefore,F[f(¢)] = Í foei at
0 ‘a
=- Í e@ e J! at +f erg Ie" dt
- 0
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
102 Digital Signal Processing
—32n-24n-16n-8n 0 lón
Ea |
24x 327
40n
eY
fit) = px einwot
acess
17?
where, Ch= > ree dt
-T12
T/2
atT faire -jnoot qt = T
~T/2
tf(t)
t
-37 -27 -T 0 T 2T 3T
=a
n(2)
ja dt
ti
tet r
=Q [4%], a =A?
Since the signal has finite power, it is a power signal and E, = œ.
Fig. E2.24(a)
Solution
Here f(t)=A, for-Tst<0
=-A,for0<t<T
=0, otherwise
- 0 e
F(jo)= Jre dt= O dt+ OT dt
= 0
o y
= fro el dt+ fro ei% dt
-TF 0
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
110 Digital Signal Processing
o T
= alfre? eft do + frei”? ejot 2o)
2n T o
= lei” [=] lei [e]
2 jt jy 2 Jt Jo
2 [se] 2 fest
J
1 : P pia xi Şi ái
=— [e72 Z eie jTt +e jni2 eiT! -e sat |
1 2sin?( =) T*t Tt
—{1-cosTt
[ ] =—_12/ = — sinc? |—
t t 2 2
ae i [ea + geese
2
1 eres jb)t _ ee ii
1 m i |
(a+ jœ- jb) (a+jo+
jb)
a+ jo+jb+a+
jo- jb]
7
oOo,
eee
r
a J (a + jo)? +b?
a+ jo
(a+ jo)? +b?
Ei Eri jeto Jy
+ ieK aina - T aar 2 |
2| (-s(e+0)} [-j(a+o)]? Jo
šH t
2| (- i)? Ca +o)
;
1 1 1
mi =t
2|(-a+0) (a+) z
1a? +@? +2aw +a? +@? -2a0
2 (a? - a2}
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
118 Digital Signal Processing
0 Te
Fig. £2.33
Solution
Fijo)= |fe de
Here f(t) = A sin Wt for 0 < t < = 7/2 where wọ = 2n/T.
= 0, otherwise
T/S
F(j@)= fAsinogt e” dt
0
2 a (ea
m je dt
0 2j
T/2
z A ffo — Hoot dt
2j 4
(Wp - 9) (Oo + @)
A (ei -w)T/2 _ 1)(Wy +0)+ enna - 1)(© - 0)
2 (o-o
2
?)
-A 3 d -i ý
= ae} ara G (et )T/2 +e jlog ae)
Fig. Q2.12
2.13 Obtain the Fourier series for the full wave rectified sine waves
shown in Figs Q2.13 (a) and (b).
A f(t
Fig. Q2.13
ee: S —(at i A
y E ae e" (s+a)
t]
Hence, L{sin wt} = A [ecet — Lle}Po
2j
ak es DA E
eee pe
2ils-jM st+jo} s*+03
Q
H ence, „Lisin % nt}
L{si% t}= Saal
a (3.5)
4. Cosine Function
F(t) = cos Wot
We know that cos Wot = A ia karm
=-
1 1= +
1 C
sEE i
2[s- joo s+j®o s? +0
Hence, Li{cos W t} = z (3.6)
2 s? +03
5. Hyperbolic Sine and Cosine Functions
Lt f(t)= Lt sF(s)
t30* se
Proof We know that
LiF = sle fO — f0)
By taking the limit s — œ on both sides
Ltt O= Lt [sF(s) — f(0)]
So
Lt |fe“ dt= Lt [sF(s)-f(0)]
0 soe
i.e. J soe
Lt (P@e“Ide =0
0
Lt sF(s)- f (0) =0
s%
pr Sa
RC (s+ —
1 i
RC
Taking inverse Laplace transform, we obtain
ölt)
ölt- a)
u(t)
u(t - a)
e“ ult)
n,-at
tec u(t)
n!
cos(wgt) u(t)
t cos(t) u(t)
t sin(@pt) u(t)
e™ sin(wpt) u(t)
(Contd.)
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
$46 Digital Signal Processing
biralD, at ~ Pi) }
Substituting s = p; in the above equation, we get
A= oe d
Similarly , A=
1
Generally, A,= 1a F,(s)|,-p, Where
n = 0, 1, 2,...n-1.
Hís) = Y(s) _
= Laplace transform of output (3.18)
e X(s) Laplace transform of input |4) initial conditions are zero
z
Ps m— f $a + = oO
4,09 A og
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
154 Digital Signal Processing
F(s)= [f(e~™ dt
0
27 n+ DT
z dt +
dt+--+ f fitje
= fr@e-* dt + froe
T nT
°
As f (t) is periodic, the above equation becomes
T T T
= froe dt+e7*? [fOe dt + mte T ffe dt +-
0 o 0
T
=[L+e-8 +078 nt erT +-]ff@e* dt
0
= [1+ eT} (e-*7)? tent (eT) + |e)
T
where F(s)= fre dt
0
Here, F,(s) = £{[u(t)— u (t - T )) f(t)}, which is the transform ofthe first
period of the time function, and {[u (t) — u (t — T )] f @)} has non-zero only
in the first period of f(t).
When we apply the binomial theorem to the bracketed expression, it
becomes 1/(1- e~°7)
Fl) = ——
— T
r froe ita SER
ee
Therefore, Z {f(t)}= soeer Jro e`“ dt
l]
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
158 Digital Signal Processing
Hence, I(s)=en ee
7 (e+)
2
LR
a
s+R/L Ris
E
s+R/L
Taking inverse Laplace transform, we get
R
io- i-e d (3.27)
Impulse Response
For the impulse response, the input excitation is x(t) = 5(t). Hence, the
differential equation becomes
3 R
t=0
x(ĝ On’ C
mm
Fig. 3.4 Series R-C circuit
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
166 Digital Signal Processing
250
I(s)= in
(s? + 625) (s + 2)
Using partial fractions, the above equation can be expanded as
250
I(s)=
(s + 2)(s + j25) (s - 725)
A, A, A;
I(s) =|—>+ deem T
s+2 s+ j25 s-j25
a...
eee nee. a
© (2- j25)(-
j50) (25+ j2)
Ag = (s -j 25) I(s)|, = j25
7 250
(s + 2)(s + J25)}, - jo5
a. mene eee
(2 + j25)(j50) (25- j2)
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
170 Digital Signal Processing
1
Ao 0 =F(s)s7|,-9=—~|
s)s“|s=0 Fl =1
1
A. sh, 4 ee 1 eee
1 ds s+1\,-6 (s +1)”
dente DF()|4.-1= 5] ad
S"ls=-1
Therefore, 3 1 = A _1i 4 1
s“(s+1) s s (s+1
Ha 1 jet 1+e7'
s“(s+1)
x -5 -2s
Therefore, i(t)= L '{U(s)]= £72 aroei
2s“ (s+ 1)
1 z
= slt-1+e |[u(t) - 2u(¢ - 1) + u(t - 2))
Fig. E3.22
Solution Let us consider the switch be at position 2. By applying
Kirchhoff ’s law, we have
di(t) 5
0.2 —— + 4i (t) = 40
a e0
Taking Laplace transform on both sides, we get
|v(t) 12
; C fa T c
_ — j
0 m2 T 372 2T
(a) (b)
Fig. E3.28
Solution The function for the first period of the given waveform is
v(it)=1 for O<tsT/2
=0 forT/2<t<T
1 T
V(s) (s) = ———
1 ent t) dt
jro
1 T12 T
= mlfre“ dt + foe a
l-e 0 T12
i ent?
TIt [-sS |,
1 1—e 87/2
“=| s |
Alternate method to find Laplace transform of the given
periodic waveform
The input periodic pulse train can be represented as
v(t)= u(t)—u(t-T/2) + u(t-T)-u(t-3T/2) + u(t- 2T)- u(t- 5T/2) +...
Its Laplace transform is
Vis)=24 fi e872 4 eT -3572 ge-BT _e-55T/2;.]
s
= -[1- e857? 40 87 (1-087) + oT (e872)
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
182 Digital Signal Processing
= +... 4
Therefore, His) =<——__+—__-
(s+1)(s?+s+1) s+1 s?+s+1
We know that Y(s) = H(s) X(s)
1 s 1
|] =
1 1
EEE
Y (8) = | —— =a
@) (4 zo) +)
(i) |een PE
s(s+1) s s+1
Gi) 4 -1 __1
(s+1)(s+3) 2(s+1) 2(s+3)
Git) s = Z3/7 ,3/7s+1/7
(s+3)(s?+s+1) s+3 s°+s+1
1 s+1
>
:
(iv) —— i = - —— +
(s+1)(s?+s+1) stl s*+54+1
Therefore,Y(s) = 1- — +—I— -—
s+1 2(s+1) 2(s+3)
1 1 3/7 3/7s+1/7
{s+ s?+s+1 s+3 s?+s+l
1 +s+ 1
s+1 s?+s+1
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
190 Digital Signal Processing
25Q
Xo
d : WIA
|
| l }10H
a Šasa ;
} Zsa
Fig. Q3.45
3.46 In the circuit of Fig. Q 3.46, the switch S is closed and steady-state
conditions have been reached. At t = 0, the switch S is opened.
Obtain the expression for the current through the inductor.
20
oe)
Lee a
tov * a a1 pF
Fig. Q3.46
Ans : 5cos1000t
3.47 In the circuit of Fig. Q3.47, the switch S is closed at t = 0 after the
switch is kept open for a long time. Determine the voltage across
the capacitor.
T $A |
i(th=10A t m i| S Xt=0
ee | | +
Fig. Q3.47
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
194 Digital Signal Processing
Table 4.1
mlS) In(2)
(vi) If x(n) is two-sided, and if the circle |z| = rg is in the ROC, then
the ROC will consist of a ring in the z-plane that includes the circle
|z| =r. That is, the ROC includes the intersection of the ROC’s of
the components.
(vii) IfX(z) is rational, then the ROC extends to infinity, i.e. the ROC is
bounded by poles.
(viii) If x(n) is causal, then the ROC includes z = æ.
(ix) If x(n) is anti-causal, then the ROC includes z = 0.
To determine the ROC for the series expressed by the Eq. 4.2, which
is called a two-sided signal z-transform, this equation can be written as
æ -1 ~
The first series, a non-causal sequence, converges for |z| <r , and the
second series, a causal sequence, converges for |z| > rj, resulting
in an annular region of convergence. Then the Eq. 4.2 converges for
rı < |z| < rg. provided r, < rg. The causal, anti-causal and two-sided
signals with their corresponding ROCs are shown in Table 4.2. Some
important commonly used z-transform pairs are given in Table 4.3.
Table 4.2 The Causal, anti-causal and two-sided signals and their ROCs
Signals ROCs
(a) Finite duration signals
Causal
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
202 Digital Signal Processing
íd) x(n) = { 0, 0, 1, 2, 5, 4, 0, 1)
Taking z-transform, we get
X(z) = 27 + 2294 5244425427,
ROC: Entire z-plane except z = 0.
(e) x(n) = &(n), hence X(z) = 1, ROC: Entire z-plane.
(f) x(n) = &n — k), k > 0, hence X(z) = z*, ROC: Entire z-plane except
z=0
(g) x(n) = &(n + k), k > 0, hence X(z) = z*, ROC: Entire z-plane except
z=,
a”, n20
x(n)=
0, n<O
Solution The z-transform for the given x(n) is
< 1
We e know know that
tha 2a n= T i jaļ< 1
Hence, X(z)= l == z
-az z-a
This converges when |az™| < 1 or |z| > |a]. Values of z for which
X(z) = 0 are called zeros of X(z), and values ofz for which X(z) > œ are
called poles of X(z).
Here the poles are at z = a and zeros at z = 0. The region of
convergence is shown in Fig. E 4.3.
Zero at
origin z= 0
Iz|= lal
The ROC ofz™* X(z) is the same as that of X(z) except forz = 0 if k > 0
and z = œ if k <0.
Solution
147
if X(@) = —4—
1~=2z7
14,7 1 1z
Solution Given X(z) = 2 = +—2
l-1
1-=z la
1~=z t
1-=z
2 2 2
1
Therefore, x(n) = Z7! ———_ + =1 z”
1-22 2 1-1,71
)"at + Ha i -)
OHE
acy
T A s T 1 n+1
X(z2):) ———— z?
(z+ 1) (z- 1?
i 2 ý Áz :
Hence F(z)= 2 -— 2 -A , Ag
z (z+D(z-1? (z+1) (z-1) (z-1
= 2 DE a Pg:i
Ag =(z-1 Fa) |201= Cope =9g
d=| z? D 2z- z?
(z + , 3
A, = — re N ==
dz((z+1)j,-4 (z+) emy 4
= [0 + 34 $n |uin)
Alternate Method
X(z)= aeae a
(1+274)(1-274)
A A
"Thetis ‘ach
Equating the numerators, we get
1= A,(1—27!)? + Ag (1+272) (1-271)+ A, (1+ 271)
= A (1-227) +.27*) + Ay (1-277) + Ag (14 277)
Here, A, +A, +A; =1
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
222 Digital Signal Processing
2+3z)
X(z) =
a+27)(14+2
+2
ae
at 27 JQ-2Pid )
i 1 4 l_-ı
1+ 5 z7! + 1 z7?
8
EOE 94.1,-1_
Therefore, X(z) = 2+ 37 87
1,-2, 41,-8
+ 32 z
T D i
Taking inverse z-transform, we get x(n)=| ° 2° 8° 32’
T
(ii) Partial Fraction Expansion Method
24327!
ear dete )(i-2e)
A A A.
= AL+
l+z 1+7 ei
x(n) =
(pt *tta" u(n — 1)
Fig. £4.29
Solution The transfer function of the given circuit in the s-domain
can be expressed as
VRC 1
©) = S URC 341
Hi = = -——
H(z) = —
z-e?
Also, the given input function x(¢) = e`”! may be expressed in the
z-plane as
z
X(z )= peer
Y¥(z) = =Z Z
@-eT) (z-e?)
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
230 Digital Signal Processing
| =2 for 6sns10
| =3 for n>10
4 5
Ans: X(z) =1 +z +z? +23 427447754
2ft 74 rhe 4%) 4
Be Ue Pe.)
4.24 Determine the z-transform of the following sequences
ot
(a) um-4) fae =z) >1
() ôn- 5) Ans: 2°
o (Gju
íd)
1\"
3 u(-n) Ans:
depts
1
3; lal<3
1
92°
(e) 3"u(n - 2) Ans: ———
oe 1-327
4.25 Find the z-transform of the sequence x(n) = na"u(n)
1
4.45 Finindd ththe causal l sisigngnalal x(n) for Xa)X(z) = 1-
l+z`
iL27 +0
r.5oz°
s
Ans: x(n) = 104) cos (= - 71.565°)u(n)
Oe P
f 32* -4z +1
where the ROC is (i) |z| > 1 and (ii) |z| < i using the long
division method.
;
1,4,7, 10,13,
Ans: (a) x(n) = {
++ 14, 11,8, 5, gi
(b) x(n) = { +
4.48 Determine the causal signal x(n) having the z-transform
- 80 (2) |u(n)
4.49 Using (i) the long division method, (ii) partial fraction method
and (iii) residue method, find x(n) and verify the results in each
case for n in the range 0 sn <3.
z+3
(a) Xz) = 7025
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
238 Digital Signal Processing
This expression gives the output response y (n) of the LTI system as a
function of the input signal x (7) and the unit impulse (sample) response
h (n) and is referred to as a convolution sum.
5.1.2 Unit Step Response [u (n)]
The unit step sequence u (n) is defined by
0, n<O
u(n)= P nsô (5.3)
-2 -1 0 1 2 3 4 5 n
(a) (b)
e
Fig. 5.5 (a) The Unit-step Sequence u (n) and
(b) The Shifted Unit-step Sequence u (n - 2).
The step response can be obtained by exciting the input of the system
by a unit-step sequence, i.e, x (n) = u (n). Hence, the output response
y (n) is obtained by using the convolution formula as
yia)= Ð h(k)u(n-k)
k=-
Therefore,
bya (t) a
22 + ay; (t) + bys (t)
2
F [a x, (t) + b x (t)) = a la y; (t)+ b yo(t)]
+ la yi (t) +b y0) [2(ayt)+b ya (t))+ 1
Here aF [x,(t)] + bF [xo(t)] # F [ax, (t) + bx,(t)] and hence the system is
non-linear.
5.2.2 Time- Invariance
A DSP system is said to be time-invariant if the relationship between
the input and output does not change with time. It is mathematically
defined as
ify (n) = F [x (n)], then y (n — k) = F [x (n—k)] = 27 F [x (n)] (5.6)
for all values of k. This is true for all possible excitations. The operator
z™* represents a signal delay of k samples.
|Example 5.9| Check whether the following systems are linear and
time invariant.
(a) F [x (n)) = n [æ (n)?
(b) F [x (n))=a [x (n)? +b
x (n)
Solution
(a) (i) F [x (n)) = n kx (ny
Here, F {x, (n)) = n [xy (n)? and
F ix; (n))= n [x (n)?
Therefore, F [x, (n)] + F [xq (n)] =n Ilx; (n)}? + {x3 (n)}]
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
246 Digital Signal Processing
in|sober|
-in 1- log (1 + œ?)
V¥1+
o?
log(1 + w)
mS
Therefore,
1 log
2
(1 + œ’)
in{ 1 )
= = —lio
log (1+*)2
fi
diys
| =e | Eg
1+0?
i
17 loøg(1+w?)
-a ee 0
Y =
1 =æ
z
as l-az? z-a
Here the pole is atz =a and hence for the system to be stable, ja |< 1.
(b) y(n) =x (n) +e" y(n- 1)
Taking z-transform, we have
Y (z)=X(z) +e" z! Y(z)
Y (2) [1 - e" 27] =X (2)
Y (2z) I 1 he
Therefore, H(z)=
X(z) 1-etz! z-e"
Here the pole is at z = e° and hence |e*| < 1, i.e. a < 0 for stability.
5.2.5 Bounded Input- Bounded Output (BIBO) Stability
Stability is of utmost importance in any system design. There are many
definitions for stability. One of them is BIBO stability. A sequence x(n)
is bounded if there exists a finite M such that |x(n)| < M for all n. Any
system is said to be BIBO stable if and only if every bounded input gives
a bounded output.
For any linear time invariant (LTI) system, the BIBO stability
depends on the impulse response of that system. To obtain the necessary
and sufficient condition for BIBO stability, consider the convolution
property which relates the input and the output of a LTI system as
y(n)= ¥ x(n-k)h
(k)
k=-%
It follows that
Y |A(k)| <2
k=0
È [hlk)| <
k=0
|h (R)| < œ.
k=0
bilap , OSnsm
h (n)= { : (5.16)
0 , otherwise
As this is obviously of finite duration, it represents a FIR system.
He) = $, hn) ei
n=-s0
Therefore |H (e®)| =
GEDI4 GER)4
Ow) == o + arg(ejo _1)_
1) arg (ejo _1)_
1) arg (ejo _33)
Solution
To find the frequency response H (e/“):
Given,y (n)+ iya- 1)=x(n)- x(n- 1)
Taking z-transform, we get
¥(2)+ 521¥@)=X@)-27X()
1o-1). z1
¥@)[1+32 |-x@u zy
H (z) =
Y(z) _ l-z?
a X(z) 441,71
2
l-e
Therefore, the frequency response is, H (e ™) = -jù
i+i¢
2
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
270 Digital Signal Processing
p
Ho LL {Vector magnitude from the i“ zero to the frequency point
i=
q
A {Vector magnitude from the i” pole to the frequency point
nt
onw =z) , —-1<n<15
0 , otherwise
Ans: (a), (c) and (e) are causal and stable.
5.26 For each of the following discrete-time signals, determine
whether or not the system is linear, shift-invariant, causal and
stable.
(a)y (n)=x(n+7), @©)y(n)= (n), Cy (n)=nx (n)
4
(dy m=a+ 2 x(n - k), wis a non-zero constant.
Hj = ——-1___,
1+a,z-° +d9z
by evaluating its poles and restricting them to be inside the unit
circle.
5.31 A causal LTI system is described by the difference equation
y(n)=y(n-1)+y(n-2)+x(n-1)
where x (n) is the input and y (n) is the output.
(i) Find the system function H (z) = za for this system, plot the
In the first case the impulse response A (n ) is folded and shifted, and
x(n) is the excitation signal. In the second case the input signal x(n) is
folded and shifted. Here k(n) acts as the excitation signal.
Commutative Law
Convolution satisfies commutative law, i.e. x(n) » h(n) = h(n) * x(n)
x(n) ho) y(n) =x(n) *h(n) hn) x(n) h(n) = h(n) +x (n)
araia. e]
y(n) yin)
y yin)
t
2) enen E atm) ad
Fig. 6.4 Associative Property
=
17 [60+4+4)=17
1 M ERT
+4j e)? e/*]
ri -4jjee" +4j
x3 (2) = İf60
= £{60-4j(-1)+4j(-D]=15
x4(3) = 260+ (-4j) 2 + 4je792]
= F160 +(-4))(-) +450)
= <[60~ 4-4] =13
Therefore, x,(n) = [ 15, 17, 15, 13 ]
Note: From the above results, we find that the resulting sequences
obtained by both linear convolution and circular convolution have
different values and length. Linear convolution results in an aperiodic
sequence with a length of (2N- 1), i.e. seven in this case, whereas
circular convolution results in a periodic sequence with a length of N,
i.e. four in this case. Circular convolution will produce the same
sequence values as those produced by linear convolution if three zeros
are padded at the end of the two given sequences x ,(n)and x,(n).
Therefore,
y(0)= J x(k)h(- k)
k=-m
y)= J x(k)h(1-k)
k=-œ
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
294 Digital Signal Processing
When n=l
y= J x(k) h1- k)
ka-a
=... + x(—1) A(2) + x(0) ACL +)x(1) A(0) +...
1 ) + ( 1 ) ( 1 ) + 0 . . . = 2
=0+(1)(
When na2
When n=-1
=... +x(-1)
A(1) + x(0) A(-1)
+ x(1) A(-2) +...
=0+(1)(1)+(1)()+0...=2
When n=-2
When n=0
yO) = J h(k)
xi- k)
j
=... + A(0)x(0)
+ A) x(-1) +...
=0+1+0+...=1
When n=1
y= SAR) x(1-k)
how
x(3 - k)
y(3)= J h(k)
kee
=... + A(O) x(3) + ACL) x(2) + A(2) x(1) + A(3) x(0)
+Ah(4)x(-l +...
= 0 + (1) (1) + (a) (1) + (a2) (1) + (a3) (1)+0 +...
=lt+a+a?+a°
Fig. E6.4(b)
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
302 Digital Signal Processing
When n=3
y(3)= J xk)h(3-k)=0
huee
y(4)= J x(k)h(4-k)
ku-
=... + x(0) A(4) + x(1) A(3) + x(2) A(2) +...
= 0 + (1) (0) + (2) (0) + 3 (-1)+0....=-3
These sequence values are plotted in Fig. E 6.5(b).
x(n) A h(n)
Hence the left and right extremes of the convoluted signal y(n) are
calculated as
yp =x, +h, =-2+(-3)=-5
y-=x, +h, =2+0=2
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
306 Digital Signal Processing
By comparing the above equations with that of DFT pair and defining
a sequence x(n) which is identical to x, (n) over a single period, we get
X(k) = Ne,
If a periodic sequence x, (n) is formed by periodically repeating x(n)
every N samples, i.e.
xpin)= $ x(n-1N)
=-~%
Hence,
x’(-n,(mod N)) = x"(N - n) ae x(k)
8. Circular Convolution
9. Circular Correlation
For complex-valued sequences x(n) and y(n),
m $ x(n) e7211™3
n=0
5
= Dy x(n) e717"
nae F A
= 1407" + QeF™ + 2e" 4 Betr 4 Be V5m
=1-—1+ 2(1) + 2(-1) + 3(1)+3-D=0
For
k =4
5
n=0
= 14 eF4™3 4 278m3 4 2e vAn 4 BoV EMS 4 3p V20K/3
= 1+ (-0.5 + j 0.866) + 2(-0.5 -j 0.866) + 2(1)
+ 3(—0.5 + 0.866) + 3(-0.5 -j 0.866)
= -1.5 -j 0.866
Fork=5
= ¥ x(n) a
n=0
as ~j5x/3 -j10x/3 jn ~j20%/3 ~j25n/3
=l+e +2e +2e7™ + 3e +3e
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
318 Digital Signal Processing
(1+ 2(-j)+3(-)+4))
(-2+2j)=-=+j
Cun
Aj
Dje
Sia
3
Given N= 4, x(n) = +>: X(k) ef ™*2 O<n<3
k=0
When
n =0
3
x0)= 1 F, Xhe
4 2
is j
x(1) e x X (k)
ef*#/2
Table 6.1
we ae Teg:
be = a > « Baa-Wnb
-1
Fig. 6.9 Basic Butterfly Flow Graph for the Computation in the DIT FFT Algorithm
In the 8-point DIT FFT flow graph shown in Fig. 6.8, W,°,W,* and W,®
are equal to 1, and hence these scale factors do not actually represent
complex multiplications. Also, since W, Wi, and W equal to —1, they
do not represent a complex multiplication, where there is just a change
in sign. Further, W,', Wit, Wê and W? are j or -j, they need only sign
changes and interchanges of real and imaginary parts, even though they
represent complex multiplications. When N = 24, the number of stages
of computations is L = log, N. Each stage has N complex multiplications
and N complex additions. Therefore, the total number of complex
multiplications and additions in computing all N-DFT samples is equal
to N log, N. Hence, the number of complex multiplications is reduced
from N? to N log 2N.
In the reduced 8-point DIT FFT flow graph shown in Fig. 6.10, there
are actually only four non-trivial complex multiplications corresponding
to these scale factors. When the size of the transform is increased, the
proportion of nontrivial complex multiplications is reduced and N log,N
approximation becomes a little closer.
The reduced flow-graph for 16-point decimation-in-time FFT
algorithm is shown in Fig. 6.11.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
330
S
OL
(0)¥p= —o `> Aa K] o 0z=(0)K
o—¢ pit'z/-8z8's-=(
gM b= l-
o v= (b)x eS
€=(2)x
—o
Digital Signal Processing
0=(2)x
z=(9)x o O!-z2t0=
pty ()x
1=8m
=98
aoo
——oz=(
0=(p)x ;)x
=2@
o y=(€)x . 0=(9)x ©
‘314
s193
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
334 Digital Signal Processing
Using DIT FFT algorithm, we can find X(k) from the given
sequence x(n) as shown in Fig. E 6.17.
Therefore, X(k) = (255, 48.63 + 7166.05, —51 + 7102, -78.63 + j46.05,
—85, —78.63 - j 46.05, -51 —j102, 48.63 — 7166.05}
| Example 6.18 |Given x(n) = {0, 1, 2, 3}, find X(k) using DIT FFT
algorithm.
Solution Given N= 4
- (2%)
Wize tw)
Wf =1 and Wi =e7*? =-j
Using DIT FFT algorithm, we can find X(k) from the given
sequence x(n) as shown in Fig. E 6.18.
Therefore, X(k) = {6, -2 + j2, -2, -2-j 2}
=1
) x(n) Wy" + Wy? +
NADY (N /2 )- 1
= SY x( n +. N/ 2) Wr * (6.24)
n=0
n=0
_ ,2n N
Since, WiN/2* =e 45 2" = cos (nk) —j sin nk = (-1)*, we obtain
(N/2)-1 (N/2)-1
X(k)= J, x(W +CD? J, x(n+.N/2)
wy"
n=0 n=0
(N/2)-1
= > [x(n + D" x (n+ Z)] wy (6.25)
n=0 j 2
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
338
(05 oy*
<
<
ox
(0)
(o)x
A
4-
La
wy qwiogz
Lag
(ixo ae a 7 gs — a (px
(8
Miik
=
Nm
(2)
6Oe
O—-—
(ax> X J l:
Digital Signal Processing
(Ne
()x N ra i
(px HP
X Qa e 05° E
cca
e e
L ee
(g)x° o
WaS Z Wied la o~ (3)x
(9)xO aa2
a a eane
WEA
ax
b= l-
4=(o)x
d-— esz = (0x O--—~¢-
x
gM
Z=(L)¥
= 0—4 —¢+ O— $8-=(r)x
0
89 am iS-
p=(z)x
d óe —— l Q -4 pee
o 1S-=(2)Xx
Z01/+
SO pare
Digital Signal Processing
Q Q O O O (Ex
= +6984- S0'9p/
t9
o us
=(9)x
(2)x
= 82k —p Q so'991 /-€9°8h =(L)x ©
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
346 Digital Signal Processing
X(0) = 6 _ x(0) =0
XN =-2 + 20—
X(2)=-2
Fig. E6.23
Hence, x(n) = {0, 1, 2, 3)
(0)x
= 9€
(0)x
=1
oO
9s9'°6/+p-=(1)x g
(p)x= g
€=
(2)x
(6X
= +p- 9S
Digital Signal Processing
9 9'1/ L=(9)x
vl+y-=(2)x
v-=(b)x
z=(4)x
9s9'1/
o 9 = (g)x
y!-y-=(9)Xx
T
=(€)
y x
Z S J bi i
i= i ; A
O O Q F : c 8 =(Z)x
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
354 Digital Signal Processing
AlsoX;
,(k + 3) = X; (k)
2
X,(k) = F x(2n) W3"
n=0
x(3)=4° 2 X(1)=-3+/5.196
x(2)=3°
gs We
x(5)=6° PIY O AN X(5)=-3-/5.196
(5) wg %(1)=-3
We
Fig. E6.29(b) DIT FFT Flow Diagram for N = 6 = 3.2
|Example 6.30 |Develop the DIT FFT algorithm for decomposing the
DFT for N = 12 and draw the flow diagram.
Solution For N = 12 = 3 . 4, where m, = 3 and N, = 4, Eq. 6.40
becomes
. i ;
Xh) = Yx(8n) Wy" + Yx(Gn + DWG”
n=0 n=0
3
+ ¥x(n +2)wiyrt?*
n=0
3 3
= Y a(n) Wit + Wk Y x(8n+)W;*
n= 0 n=0
3
+ Wa J x(n +2)wy*
n=0
=X (k) + W$,
Xa(k) + WÈ
Xa(k)
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
362 Digital Signal Processing
x(8)=0 Xi2}=0
x(12)=0 p X(3)=0
x(1)=1 i> Xx(4)=0
E A>O x(6)=0
x@)=1 of
; LA2SANII A
LL I
ESS
SP
A NA
4
A SSBKEF Ce.
ASOs
aSLALA
x(13)=1 } SP A0
7
OY» RO VLN
cy ZSA
— PRY
K KY
BROW
OAA See
<>
x(6)=0 Vi SEVA
SA) CO
BE R S
NE x X(9)=2. +j2.6112
'9)=2.. 6112+/2
>
0 >» WX
1
= Dlr) + x(n + 2) WH + x(n + 4) WS") wg
n=0
i
X(3k) = Y [x(n) + x(n + 2) + x(n + 4] WS"
n=0
1
X(3k + 1) = Z [zn +x(n+ 2)W2 +x(n+ a)wijwg wgn
n=0
1
X(3k + 2)= > [x(n)+ x(n+2)Wé +x(n+ 4) we ]we” wink
n=0
X(2)
x(5)
we we We
Fig. E6.32 (b) DIF FFT Flow Diagram for Decomposing the DFT for N = 6 = 2.3
Fig. E6.34(a)
Steps (b), (c) and (d) are described below using the direct
implementation of circular convolution.
Circular convolution of data blocks x,(n) and x(n) with h(n)
padded with (N — 1), i.e. five zeros is given in Fig. E6.34(b) and (c).
CIEE
on Bye oe EE
pat (c)
Fig. E6.34 (b) and (c)
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
374 Digital Signal Processing
or, equivalently,
or, equivalently,
SO or —-NSOS+T
o
Therefore,
O(o)=-
at
where tis a constant phase delay expressed in number of samples. Using
Eq. 7.2,
Dlo)
-1 Im
= tan! ———_-
H(e?”) = - wt
Re H(e’*)
or
M-1
»y h(n)sinan
-1 n=0
ot = tan) o
¥ h(n oson
n=0
or
M-1
2, h(n) sinon
tan t= fy
¥ h(n) cos wn
n=0
Simplifying, we get
M-1
¥ A(n) sin (ot - wn) =0 (7.3)
n=0
and
h(n)
= A(M - 1- n) for
0 < n < M-1 (7.5)
If Eqs 7.4 and 7.5 are satisfied, then the FIR filter will have constant
phase and group delays and thus the phase of the filter will be linear.
The phase and group delays ofthe linear phase FIR filter are equa! and
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
386 Digital Signal Processing
where
and re
type-I design and type-II design. In the Type-I design, the set of
frequency samples includes the sample at frequency œ = 0. In some
cases, it may be desirable to omit the sample at œ = 0 and use some
other set of samples. Such a design procedure is referred to as the Type-
II design.
Type-I Design
The samples are taken at the frequency
27k
mAr =0,1,...,M-1
Ok M’ (7.22 )
The samples of the desired frequency response at these frequencies
are given by
If these numbers are all reai, then these can be considered as the
impulse response coefficients of an FIR filter. This can happen when all
the complex terms appear in complex conjugate pairs, and then all the
terms can be matched by comparing the exponentials. The term Hk)
e/?""*/M should be matched with the term that has the exponential
e/?xnk/M as a factor. The matching terms are then H (k) e/?*"*/™ and
H (M — k) e/?*"(M-0M since 2nn (M-k/M = 2nn — (2nnk/M). These
terms are complex conjugates if H (0) is real and
(i) For M odd:
H (M-k) = H“(k), k = 1, 2,...,(M-1)/2 (7.25)
(ii) For M even:
H(M-k) = Ë (k), k= 1, 2,..., M/2-1 (7.26)
H(M/2) =
The desired frequency response H; (e/®) is chosen such that it satisfies
the Eqs 7.25 and 7.26 for M odd or even, respectively. The filter
coefficients can then be written as
as 2nn
p OSnsM-1 (7.44)
0, otherwise
The window function of a non-causal Hanning window is expressed by
l-n, l<n<
0, otherwise
Therefore, with M = 7,
w(0) = 0.08, w(1) = 0.31, w(2) = 0.77, w(3) = 1,w(4) = 0.77,
w(5) = 0.31, w(6) = 0.08.
The filter coefficients of the resultant filter are then,
h(n) =hy(n).w(n) n=0, 1, 2,3, 4, 5, 6.
Therefore,
h(0) = 0.006, A(1) = — 0.0494, h(2) = 0.1733, h(3) = 0.75,
h(4) = 0.1733, h(5) =- 0.0494 and h(6) = 0.006.
The frequency response is given by
6
H(e/®) = YAM) ejen
n=0
where
f,=0.5(f,+f,) and AF=f,-f, (7.66)
Bandpass FIR Filter
where
AF AF
fa=fa zy fee =fp2+ KE
where
AF AF
fer =fp1+ => fee ie
There are four different cases that result in a linear phase FIR filter,
viz., (i) symmetric unit impulse response and the length of the filter, M
odd, (ii) symmetric unit impulse response and M even, (iii) anti-
symmetric unit impulse response and M odd, and (iv) antisymmetric
unit impulse response and M even. The first case is discussed below in
detail and other cases are listed in Table 7.2.
In the symmetric unit impulse response case, h(n) = H(M - 1- n). The
real-valued frequency response characteristics |H(e/)| = |H,(e/®)|,
given in Eq. 7.14, is
(M -3)
where
Table 7.2 Magnitude Response Functions for Linear Phase FIR Filters
From Table 7.2, it can be seen that the magnitude response function
can be written as given in Eq. 7.77, for the four different cases.
|H(e/®)| = Q(@) Po) (7.77)
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
416 Digital Signal Processing
fa tm(z)
Re(z)
Unit Circle
s-plane
Fig. 8.2 The Mapping of Eq. 8.5 into the z-plane
It can be seen that the mapping of Eq. 8.5 takes the left-half plane
of s-domain into the corresponding points inside the circle of radius 0.5
and centre at z = 0.5, and the right-half of the s-plane is mapped outside
the unit circle. As a result, this mapping results in a stable analog filter
transformed into a stable digital filter; however, as the locations of poles
in the z-domain are confined to smaller frequencies, this design method
can be used only for transforming analog low-pass filters and bandpass
filters having smaller resonant frequencies. Neither a high-pass filter
nor a band reject filter can be realised using this technique.
The forward difference can be substituted for the derivative instead
of the backward difference. This gives,
dy) _ y(nT+T)-y(nT)
dt T
vad va) (8.12)
The transformation formula will be
z-1
s = -—
T (8.13 )
or,
z=1l+sT (8.14)
The mapping of Eq. 8.14 is shown in Fig. 8.3. This results in a worse
situation than the backward difference substitution for the derivative.
When s = j Q, the mapping of these points in the s-domain results in a
straight line in the z-domain with coordinates (Zreai Zimag) = (1,27).
Consequently, stable analog filters do not always map into stable digital
filters.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
Infinite Impulse Response (IIR) Filters 425
jQ Im Z
Unit Circle
Taking T = 1s,
H(z)= u z
1- (0.8187) ( 220) z -1 =
1- 2(0.8187) (- 0.99) z~- + 0.6703 z
That is,
H(@)= 1+ (0.8105) z -1
1+ 16210 z~} + 0.6703 z~?
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
430 Digital Signal Processing
[Hw | 8 ô
|H(jQ)|A 4p
Q, a
Fig. 8.5 Relationship Between @ and Q as Given in Eq. 8.40
2zn +0.1
H(z) = a O
[2 e—»
T (z+)
+oa ] +9
_ __(2/T) (2-1) (z+) +0-1(2+ 1)?
[(2/T) (z - 1) + 0.1(z+ D? +9 (z + 1)?
Substituting T = 0.276 s,
TEE Q E A2 ETy
[cavaz) - 1°" [v09 - 3]
Step (iv) Determination of H, (s).
From Eq. 8.49,
His)= BQ, (N-1/2 B, Q?
8+09Q., py 8° +b S+ Q?
-( By Q. | B, Q2 )
$+)Q, )\ 8? +b,Q.8+0,Q?
From Eq. 8.50,
That is,
2.5467 6.4857
Mers GED (z-1)
RE
zg ar [2auf
- + 2.54675 + 6.4857
Simplifying we get,
H(z) = — 16.5171
a (z + 1)
70.832° + 31.1205z? + 27.2351z + 2.948
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
Infinite Impulse Response (IIR) Filters 445
e? Ch (Q,/Q)
<8, 222, (8.63b)
1+e? CÈ (Q,/9)
When Q = Q, Eq. 8.63b becomes
a= —
2 +e?
Rearranging,
öz
£= Gna (8.64)
8.12 An analog filter has the following system function. Convert this
filter into a digital filter using the impulse invariant technique.
1
H(s) = —~~
ý (s+0.1)7 +9
8.13 Convert the analog filter to a digital filter whose system
function is
H(s) = ——_,
D= rap
8.14 Convert the analog filter to a digital filter whose system
function is
H(s) =(s+0.1)?
-—*8 +36
The digital filter should have a resonant frequency of œ, = 0.2 n.
Use impulse invariant mapping.
8.15 What is bilinear transformation ?
8.16 Compare bilinear transformation with other transformations
based on their stability.
8.17 Obtain the transformation formula for the bilinear
transformation.
8.18 An analog filter has the following system function. Convert this
filter into a digital filter using bilinear transformation.
1
Hi= yaa) +6
8.19 Convert the analog filter to a digital filter whose system
function is
r 1
H(s)
~ (s+2)? (s+1)
using bilinear transformation.
8.20 Convert the analog filter to a digital filter whose system
function is
H(s) =
36
(s +0.1)? +36
The digital filter should have a resonant frequency of œ, = 0.27.
Use bilinear transformation.
8.21 What is meant by frequency warping ? What is the cause of this
effect ?
8.22 Describe Butterworth filters ?
8.23 Comment on the passband and stopband characteristics of
Butterworth filters.
8.24 Describe Chebyshev filters ?
8.25 Describe inverse Chebyshev filters ?
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
Realisation of Digital Linear Systems 455
A7
z
|
i
P
Xz) )
gkbrranchanthy y z z by
noor
Y(z)
ay
a b |
4
(a) (b)
Fig. 9.2 (a) Basic Realisation Block Diagram Representing a First-order
Digital System and
(b) Its Corresponding Signal Flow Graph.
(a)
Fig. E9.2
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
464 Digital Signal Processing
The cascade realisation of the system transfer function is shown
in Fig. E9.4.
Fig. 9.9 Parallel Form Realisation Structure With the Real and Complex
Poles Grouped in Pairs
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
472 Digital Signal Processing
Hia) = ———___4+—___~
(1+4z Tra +r
$ *)( i j 3 a
Solution
Cascade Realisation To obtain the cascade realisation, the
transfer function is broken into a product of two functions as
H(z) = H,(z) H(z)
14427 1
where H,(z)= —— and H,(z)= te
1+=2z7! 1+=2724+=27
2 2 4
The cascade realisation structure for this system function is shown in
Fig. E9.8(a).
x(n)
x(n)
Fig. E9.8(b)
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
498 Digital Signal Processing
The largest error occurs when all the discarded bits are one’s. When
the number x is negative, truncation results in reduction of the
magnitude only. However, because of the negative sign, the resulting
number will be greater than the original number. For example, let the
number be x = — 0.375. That is, in sign magnitude form it is represented
asx = 1011 and after truncation of one bit, Q(x) = 1 01. This is equivalent
to — 0.25 in decimal. But — 0.25 is greater than — 0.375. Therefore, the
truncation error is positive and its range is
O<ep (28-274) (10.3)
The overall range of the truncation error for the sign magnitude
representation is
— (2-8 -2°*) Seps (27-2) (10.4)
(ii) Truncation error for two’s complement representation When the
input number is positive, truncation results in a smaller number, as in
the case of sign magnitude numbers. Hence, the truncation error is
negative and its range is same as that given in Eq. 10.2. If the number is
negative, truncation of the number in two’s complement form results in
a smaller number and the error is negative. Thus the complete range of
the truncation error for the two’s complement representation is
-(2%-2+")<e,<0 (10.5)
äi) Round-off error for sign magnitude and two’s complement
representation The rounding of a binary number involves only the
magnitude of the number and is independent of the type of fixed-point
binary representation. The error due to rounding may be either positive
It has been assumed that the noise resulting from the quantisation
process is a white noise. For this case, we have
Voe (M)=02, and %e(m)=02 (10.21)
where oĉ, is the output noise power (or power of the output error) and
6? is the input noise power. Using Eq. 10.21 in Eq. 10.20 and replacing
the variable & with n,
Y= f HOH ez dz
n=0 2n J Cc
where E (e/“) is the frequency response of a linear phase FIR filter that
has {e(7)} as the impulse response for the first half and the second half
can be obtained using e(n) = e(M — 1 —- n). Thus, the filter with its
coefficients rounded-off can be considered as a parallel connection of the
ideal filter (infinite precision) with a filter whose frequency response is
E (e/%) e12 ™ -2 Since the error e (n) due to rounding of the filter
-B :
|E (e/®)| <
(E +2
E emils (4 -n)a]
J le(n)|
n=0
(10.36)
Eag
[E jjj (e7®)|
< sJe 7 (=> +a) |cosk œ|
k=1
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
512 Digital Signal Processing
0.5
ls (0 <bz< 1) (10.48)
1-Jbl
Table 10.3
Fin-D y (n-2 a, Y (n~ 1)-az 9 n-2)| F (M=Qell| Fn)
in decimal
¥ (-1)20,110 | F (-2)= 1,010 1,010100
¥(O)=1,011 | ¥ (-1)=0,110 10,110011
y (1)=0,110 | ¥ ()=1,011 1,001101
ï (2)=1,010 | ¥ (1)=0,110 10,101100
J (3)=0,110 | F (2)=1,010 14010100
For n = 1, the sum of two products has resulted in a carry bit to the
left of the sign bit that is automatically lost, resulting in a positive
number. The same thing happens forn = 3 and also for other values ofn.
It can be noted from the above table that the output swings between
positive and negative values and the swing of oscillations is also large.
Such limit cycles are referred to as overflow limit cycle oscillations.
The study of limit cycles is important for two reasons. In a
communication environment, when no signal is transmitted, limit cycles
can occur which are extremely undesirable. For example, in a telephone
no one would like to hear unwanted noise when no signal is put in from
the other end. Consequently, when digital filters are used in telephone
exchanges, care must be taken regarding this problem. The: second
reason for studying limit cycles is that this effect can be effectively used
in digital waveform generators. By producing desirable limit cycles in a
reliable manner, these limit cycles can be used as a source in digital
signal processing.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
Multirate Digital Signal Processing 525
rs
|
t=nT
x(t) x(n)
x(t)
0 t
s(t)
sell CRA DR G
J [i PA
T O T 2T 3T 4T 5T 6T t
x(nT)
-T 0 T 2T 37 47 ST 6T nT
F F=LF F
x(n) re | x(t) yim)
Sampling rate
expander/
up sampler
ix(e!®)j
2r e!
Iwe”)
0 nil x 2r w
[Y (01)
T
0 niL x 2n a
Sadie Pa n = multiples of M
0, otherwise
where M = 2.
x(n)
=0, 1, 2, 3, 4, 5,...
y(n)
= 0, 0, 1, 0, 2, 0, 3, 0, 4, 0, 5, 0,...
In general, to obtain the expanded signal y(n) by a factor M,
(M — 1) zeros are inserted between the samples of the original signal
x(n).
The z-transform of the expanded signal is
Y) =X(z™),
M=2.
The input and output signals are shown in Fig. E 11.2.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
536 Digital Signal Processing
te y(m) = (x(mL),
0 otherwise
(g) Modulator
s(n)
Fig. 11.9 Branch Operations in Signal Flow Graphs
node of the branch. External signals enter the input branches and
signals at the output branches are terminal signals. The sum of the
signals entering the node is equal to the sum of the signals leaving the
node. Based on the signal flow graph of Fig. 11.10, the network
equations can be written as follows.
At the input node,
r(n) = x(n) + a,r(n — 1) + agr(n — 2) (11.21)
At the output node,
y(n) = r(n) + birin — 1) + barin - 2) (11.22)
Combining both the equations,
y(n) = x(n) + byx(n — 1) +b, x(n — 2) + a yin — 1) + aQy(n—-2) (11.23)
r(n)
x(n) y(n)
r(n-2)
(11.26)
= 4 E >
an-1 by
Solution
H(z) = Ey(z*) +27) Elz”)
where Ey (2°), E (Z 2) are polyphase components.
_ 1-427
H(z)
1+5271
(1-427!) (1-527)
© (14527?) (1-527)
_ 1-927! +202?
1- 2527
1+2027 ., -9
so tz? CO
1- 2527? 1- 2527?
The polyphase components are
E (z) =
1+ 202%
1-2572 and E2) = ia.
1- 2527?
11.6.1 General Polyphase Framework
The z-transform of an anti-aliasing filter shown in Fig. 11.17 (a) with
impulse response A(n) is given by
Decimation
filter
H(z) = x h(n) 2”
n=0
= A(0) + A(1) 271+ h(2)277 +...
which can be partitioned into M sub-signals where M represents
decimation factor. Hence,
H(z) = A(O) + A(M) z~™ + hM) 2-2 +...
+z AC) + AM + 1) 27! 4+ (QM +1277 +...)
+27 Dih(M-1)+h(Q2M-1)z-™ +...} (11.32)
Equation 11.32 can be written as
M-1 ~
He)= > F himM + ky zr
k=0 m=0
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
550 Digital Signal Processing
x(n) 4
i
i
~ r} ” | m |e á
Fig. 11.18(e)
Polyphase
inplementation
x(n) — > te} ARoz?) |
y(n)
y(n)
ä ao]
0 40 50
a LPF 1
[ioA ni|
400 Hz LPF 2
f(Hz) E AE | f(Hz)
w El ais areal and Two-stage Network he Sones
eW | T ~ him) || fh=5f
> y(m)
Hie]
|H(e?°) |
A Passband Stopband
x ar, 2 Fs,
> je — m
—- — n and
Fp N N
f 2h,
Fig. 11.28 The Comb Filter for N = 5
Ru om 12 x 2000 = 12,000
Ry, p = 30x 20% = 60,000
The overall number of multiplications per second for a three-stage
realisation is given by
Ru, c + Ryu,s + Ru, p= 1,25,500
The number of multiplications per second for a three-stage realisation
is more than that of a two-stage realisation. Hence higher than two-
stage realisation may not lead to an efficient realisation.
Comb Filters
The impulse response of a comb filter (FIR filter) is given by,
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
Multirate Digital Signal Processing 581
H(z) 14 t4 Gy(z)
c ta G,(z) y(n)
t2 G,(z)
REVIEW QUESTIONS
11.1 What is the need for multirate signal processing ?
11.2 Give some examples of multirate digital systems.
11.3 Explain the interpolation process with an example.
11.4 Explain the decimation process with an example.
11.5 Write the input-output relationship for a decimation processing
a factor of five.
11.6 With an example explain the sampling process.
11.7 What is meant by aliasing ?
11.8 How can aliasing be avoided ?
a", n>0
11.9 The signal x(n) is defined by x(n) = {
0, otherwise
(a) Obtain the decimated signal with a factor of three.
(b) Obtain the interpolated signal with a factor of three.
1.10 Explain polyphase decomposition process.
1.11 How can sampling rate be converted by a rational factor M/L ?
1.12 Draw the block diagram of a multistage decimator and
integrator.
1.13 What are the characteristics of a comb filter ?
1.14 Explain with block diagram the general polyphase framework
for decimators and interpolators.
1.15 What is a signal flow graph ?
(Contd.)
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
586 Digital Signal Processing
Variance
Var a(fii=Tair)=|1(FE2E
L) |
_ 1 (sinnfN .
where =Waar f) = 77
sin nf
When N > œ,
v2
E [P xf Tax(f) | Waan(@)d@
-1/2
= Wart (0) T,, A)
=T,,(f)
var [P,,(f lo Tr f)
This is asymptotically unbiased estimate, but not consistent as
variance does not approach zero when N > œ.
The quality factor is,
-GO
v=
which is constant and independent of N specifies the poor quality.
Bartlett Power Spectrum Estimate
Mean
1/2
EL PR" M]= f TOW pon (f- Odo
-1/2
Variance
ver Pi on] =z eo aR |
x 2
r 2
where Waa (f) = Š (rr)
EIP (Fy) Px f)
ny =0 ng=0
N-1 Rai
24 x e712" fi- fa) > ef 2h - fans
4
Ss ny=0 ng =0
~ Nè N-i ’ N-1 |.
$ X eizh + f)ni 5 ei 2h + fh) ng
nysO ng =0
where
r,, (k) is the autocorrelation sequence of x(n)
ry (k) is the crosscorrelation sequence of x(n) and d (n)
If d(n) = x(n),
N
Y, ay (k) ralk -D = ry D1 = 1, 2, -n N
k=1
The above set of equations are called normal equations or the Yule-
Walker equation.
Solve the equations recursively. First consider a predictor of order
one, which is given by
Ny D
a, (1) =
Ty (0)
The least squares error becomes
È o eiki
In matrix form
E"(fi)a=1
where Efi) =11 e771... e/?*Pfy
Minimise the variance o? subjected to the constraint specified, yields
an FIR filter which allows the f, frequency components undistorted.
Other frequency components are attenuated. This yields,
â =TQ EU VE (fi)Ta Ef)
The variance becomes
ite
ee Ss
mn” ETETEA)
The minimum variance power spectrum estimate at frequency f; is
represented in the above equation. By varying the frequency f; from 0 to
0.5, the power spectrum estimate can be obtained. Even if f, changes,
Tz is computed only once. The denominator of o2;,, can be computed
using single DFT. If R, is the estimate of Tx, Rẹ can replace I,.,and the
minimum variance power spectrum estimate of |Capon’s method iis
where
Dy (n) = }, w`’ Xy O
1=0
—estimated crosscorrelation vector
The solution can be obtained as
hy = Rj} (n) Dy (n) (13.39)
Ry (n) and Dy (n) can be computed recursively by
Ry (n)= w Ry (n — 1) + Xy (n) X(n) (13.40)
This is known as the time update equation for Ry (n)
Dy (n)
= w Dy (n - 1) + d (n)Xy (n) (13.41)
Matrix Inversion Lemma
Let A and B be two positive definite M x M matrices, D is a positive
definite N x N matrix and C is an M x N matrix.
A=B'+¢C-D'cT (13.42)
A`! can be obtained from matrix inversion lemma as
A` =B -BC [D + C” B C)’ C7 B (13.43)
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
Adaptive Filters 653
Dw Y wo, GD (13.65)
n-1 ¥
Km- = È w” fan- O bm- 0- V+ fp 1) By 1 0-1)
i=l
n-i
x w) i faa ibm- li-l) + fp_-1 (2) by (n1)
i=1
=w Ky —1(n-1)+ fy 1 (2) bm (0-1) (13.66)
Similarly,
EQ _,()=wE® _, fi, 1) vse
and
E®_\(n)=wE®_,(n-1)+62,_,()
(13.68)
ED (n= F wi Pw
isl
‘ ye [fm 1) +P Mba- i- DY
;
F E urio- i=l
2
= EY? Bae iy. | sae |E® (n-1)
m1 BO mn- D Ewin-D
= EP (n) -ram
K?
(13.70)
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
cin 14
Applications ofDigital Signal
Processing
14.1 INTRODUCTION
Digital signal processing techniques are used in a variety of areas which
include speech, radar, sonar, image, etc. These techniques are applied
in spectral analysis, channel vocoders, homomorphic processing
systems, speech synthesisers, linear prediction systems, analysing the .
signals in radar tracking, etc.
[l
= >
n-an
E x MAMT-rT)eF"? — (14.7)
m=0 r=a-~(m-DN+1
yji w-1 P
Y YLxtal-rT-mNT)ACT
-mNT) een?
m=0 l=0
an Ne eae IT- e]
l=0 m=0 Adt+mNT) eFr)
where
[$]
gil, n)= > x(nT-1IT-mNT)hUT
+ mNT) (14.10)
m=0
i - wavelength
Pulsed Doppler signals can be used to obtain both range and velocity
resolution.
14.3.1 Signal Design
Transmitting narrow pulse provides good range but poor velocity
measurement. A wide pulse of single frequency gives good velocity but
bad range information.
Consider the radar model shown in Fig. 14.16. Let the signal be
generated digitally and transmitted through an analog filter. The
transmitted signal is s(t). The received signal is s(t — t) e?“ -© which
is delayed and frequency shifted. The received signal is passed through
To transmitter From
and antenna receiver
Y
in time and centered at the location of the Dirac. For small a’s, the
transform “zooms-in” to the Dirac with good localisation for very small
scales.
Frequency localisation
Consider the sinc wavelet, i.e. a perfect bandpass filter. Its magnitude
spectrum is 1 for || between n and 2r. Consider a complex sinusoid of
unit magnitude and at frequency @). The highest frequency wavelet that
passes the sinusoid having a scale factor of m/w, (gain of Vx/«,) while
the low frequency wavelet that passes the sinusoid having a scale factor
of 27/0 (gain of /2n/a,).
(viii) Reproducing Kernel
The CWT is a very redundant representation since it is a 2-D expansion
of a 1-D function. Consider the space V of a square integrable function
over the plane (a, b) with respect to da db /a?. Only a subspace H of V
corresponds to wavelet transforms of functions from L?(R).
If a function W (a, b) belongs to H, i.e. it is the wavelet transform of
f(t), then W (a, b) satisfies
1 da db
W (ap, bo) = Z SJE oboa, b) W (a, b) a? (14.34)
where Klao, bo, a, b) = < Wa, bẹ Wa, b > is the reproducing kernel.
al
1
4
Fig. 14.22 Time-Frequency Cells that Correspond to Dyadic Sampling
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
704 Digital Signal Processing
1 a 8 T
Gain
dB
in
——>-
\1 $38
dB
in
Gain
——>
$
$
1 N
radians
in
Phase
—>
(4 0.1 0.2 0.3 0.4 05 06 07 0.8 09 1
(b) Normalised frequency —>
(n,wn)=cheblord(wl,w2,rp,rs);
[b,a)=chebyl(n,rp,wn, ‘high’);
w=0:.01/pi:pi;
[h, omJ =freqz(b,a,w);
m=20*10g10
(abs (h) );
an=angle(h);
subplot (2,1, 1);plot(om/pi,m);
ylabel (‘Gain indB-->’);xlabel(‘(a) Normalised frequency -->’) ;
subplot(2 ,1,2);plot(om/pi,an) ;
xlabel (` (b) Normalised frequency -->');
ylabel (‘Phase in radians -->’);
As an example,
enter the passband ripple... 0.3
enter the stopband ripple... 60
enter the passband freq... 1500
enter the stopband freq... 2000
enter the sampling freq... 9000
The amplitude and phase responses of Chebyshev type - 1 high-pass
digital filter are shown in Fig. 15.24.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
MATLAB Programs T31
an=angle(h);
subplot (2,1, 1);plot(om/pi,m);
ylabel (‘Gain indB-->');xlabel(‘(a) Normalised frequency -->’);
subplot (2 ,1,2);plot(om/pi, an);
xlabel(* (b) Normalised frequency -->’);
ylabel (‘Phase in radians -->’);
As an example,
enter the passband ripple... 0.35
enter the stopband ripple... 35
enter the passband freg... 1500
enter the stopband freq... 2000
enter the sampling freq... 8000
The amplitude and phase responses of Chebyshev type - 2 low-pass
digital filter are shown in Fig. 15.27.
dB
in
Gain
———>
-100 ——— j ee | =) a a ee
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
4 oe =k —ae borar F- a. a ]
o
ee
9
Phase
in
radians
——-> L
a r a S EE S S,N |
orm 0.1 0.2 0.3 0.4 05 06 07 08 09 1
(b) Normalised frequency —->
%HIGH-PASS FILTER
b=firl(n-1,wp, ‘high’,y);
{h,o]=freqz(b,1, 256);
m=20*1log10 (abs (h));
subplot (2,2,2) ;plot(o/pi,m) ;ylabel (‘Gain in dB -->’);
xlabel(* (b) Normalised frequency -->‘);
%BAND-PASS FILTER
wn= [wp ws];
b=firl(n-1l,wn,y);
(h,o]=freqz(b,1,256);
m=20*1lo0g10
(abs (h));
subplot (2,2,3);plot(o/pi,m) ;ylabel(‘GainindB-->’);
xlabel (* (c) Normalised frequency -->’);
%BAND-STOP FILTER
b=firl(n-1,wn, ‘stop’,y);
(h,o]=freqz(b,1,256);
m=20*1log10(abs(h));
subplot (2,2,4);plot(o/pi,m) ;ylabel (‘Gain in dB -->’);
xlabel ( * (d) Normalised frequency -->');
As an example,
enter the passband ripple 0.03
enter the stopband ripple 0.02
enter the passband freq 1800
enter the stopband freq 2400
enter the sampling freq 10000
enter the ripple value(in dBs) 40
The gain responses of low-pass, high-pass, bandpass and bandstop
filters using Chebyshev window are shown in Fig. 15.34.
15.14.5. Hamming Window
Algorithm
1. Get the passband and stopband ripples
2 . Get the passband and stopband edge frequencies
3. Get the sampling frequency
4. Calculate the order of the filter
5. Find the window coefficients using Eq. 7.40
Geis Draw the magnitude and phase responses.
%Program for the design of FIR Low pass, High pass, Band pass
and Bandstop filters using Hamming window
cle;clear all;closeall;
rp=input (‘enter the passband ripple’);
rs=input (‘enter the stopband ripple’);
fp=input (‘enter the passband freq’);
fs=input (‘enter the stopband freq’);
f =input (‘enter the sampling freq’);
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
758 Digital Signal Processing
rsl=filter(h,1l,y);
efor i=1;length (ns),
% rsl(i)=rs6(i);
send
for i=l:lengthi(ns),
rs(i)=rsl(i);
end
subplot (2,2,2);plot(rs) ;title(‘noise signal’);
xlabel (‘Time ——') ;ylabel (‘Amplitude ——’);
%Far end signal
fs1=(552525252525255555555§822222222222225);
trs=sign(rs2);
%Far end signal is digitally modulated and plotted
zl = dmod(fs1, fc, fd, fs,’psk’);
for i =1 :length (ns),
z(i) =z1(i);
end
subplot (2, 2, 3);plot(z);title (‘far-end signal’);
xlabel (‘Time ——;) ;ylabel (‘Amplitude ——’);
%Echo and the far end modulated signal is added in the hybrid
ql=z1 + rs1;
for i=1:length (ns),
a(i)=ql(i);
end
subplot (2,2,4);plot(q) ;title(‘received signal’);
xlabel (‘Time ——’); ylabel (‘Amplitude ——’);
q2=xcorr(q);
+Auto correlation is taken for the near end signal
ar=xcorr(ns);
%cross correction is taken for the near end and far end signal
erd=xcorr(rs,ns);
ll=length(ar) ;j=1;
for i=round(11/2):11,
arl(j)=ar(i)
j=j+1;
end
$Toeplitz matrix is taken for the auto correlated signal
r=toeplitz (arl);
12=length(cr ;j=
d)1;
for isround (12/2):12,
eral (j)=crd(i);
j=j+1;
end
p=crdl';
%Maximum and minimum eigen values are calculated from the
toeplitz matrix
lam=max(eig(r));la=min(eig(r));l=lam/la;
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
a
You have either reached a page that is unavailable for viewing or reached your viewing limit for this
book.
770 Digital Signal Processing
(c) Repeat parts (a) and (b) using the Hamming window
(d) Repeat parts (a) and (b) using the Bartlett window.
15.20 Design an FIR Linear Phase, bandstop filter having the ideal
frequency response
1, for|a|< 2/6
0, for 2/6 <| s 2/3
Hal) = |1, for 213 <|alsa
(a) Determine the coefficient of a 25 tap filter based on the
window method with a rectangular window.
(b) Determine and plot the magnitude and phase response of
the filter.
(c) Repeat parts (a) and (b) using the Hamming window
(d) Repeat parts (a) and (b) using the Bartlett window.
15.21 A digital low-pass filter is required to meeth the following
specfications:
Passband ripple <1 dB
Passband edge 4 KHz
Stopband attenuation 240 dB
Stopband edge 6 KHz
Sample rate 24 KHz
The filter is to be designed by performing a bilinear
transformation on an analog system function. Determine what
order Butterworth, Chebyshev and elliptic analog design must
be used to meet the specifications in the digital
implementation.
15.22 An IIR digital low-pass filter is required to meet the following
specfications
Passband ripple < 0.5 dB
Passband edge 1.2 KHz
Stopband attenuation 240 dB
Stopband edge 2 KHz
Sample rate 8 KHz
Use the design formulas to determine the filter order for
(a) Digital Butterworth filter
(b) Digital Chebyshev filter
(c) Digital elliptic filter
15.23 An analog signal of the form x,t) = a(t) cos(2000 mt) is
bandlimited to the range 900 <F < 1100 Hz. It is used as an
input to the system shown in Fig. Q15.23.
n)
if (M< =0)
{
M = fabs
(M) ;
}
M1 = fabs (20 * log10(M));
if (( F>=edge[1]) && (F< =edge[2]) && (fx[1] ==1))
{M1 =0;}
else if ( (F >=edge[3]) && (F<=edge[4]) && (fx[2] ==1))
{M1 =0;}
else if ((F>=edge[5]) && (F< =edge[6]) && (fx[3] ==1))
{M1 =0;}
if ( M1 > 85)
{M1 = 85;}
x1 =25 +1000 *F;
Y1 = 25 + M1;
X3 = 25 +1000 * F;
Y3=170-15*N;
line(X1,Y1,X2,Y2);
line (X3,Y3,X4,Y4);
}
getch( );
closegraph( );
if (nfilt !=0)
{ goto agl;}
return ;
}
Salient Features
P> Overview of Signals a Systems concepts.