Title: Summarized by
1Unit - I
- Summarized by
- Neetesh Purohit
- Lecturer, IIIT,
- Allahabad, UP, India
- http//profile.iiita.ac.in/np/
2Classification of signals
- Multichannel
- (signals are generated by multiple
sensors/sources which are measuring the same
parameter. Example measurement of ground
acceleration a few kilometers from epicenter of
an earthquake, ECG etc) - Multidimentional
- (The value of signal is a function of M
independent variables.)
3- Continuous time Vs Discrete time (may not be
discrete valued) -
- Continuous valued Vs discrete valued (quantized
continuous valued discrete time signal) - Deterministic Vs Random
- (Deterministic, if signal can be uniquely
described by an explicit mathematical expression
e. g. sine, cosine, line, parabola etc. Random
signal examples noise, unsynchronized digital
signal etc )
4Signal Types and Properties
- Energy Vs Power signals
- Energy of x(n) E ?-8ltnlt8x(n)2
- If E is finite, x(n) is called energy signal.
- Power of x(n) P lim(N tend to 8) 1/(2N1)
-
?-NltnltNx(n)2 - If E is finite P0 but if E is infinite P may or
may not be finite. (prove it?) - If P is finite then x(n) is called power signal.
5- Periodic Vs Aperiodic signals
- x(nN) x(n) for all n
- If no value of N satisfy it the signal is
aperiodic. - A sinusoid x(n) Asin2pfn is periodic if fk/N
- Power of periodic signals
- (1/N)?0ltnltN-1x(n)2
6- Symmetric (even) Vs Antisymmetric (odd)
- x(n) x(-n) symmetric (even)
- x(n) -x(-n) antisymmetric (odd)
- Any arbitrary signal can be expressed as sum of
two signal components, even and odd. - xe(n) ½x(n) x(-n) xo(n) ½x(n) - x(-n)
7Basic sequences and operations
- Unit sample/discrete time impulse/impulse
- dn 1(if n0) 0 (otherwise)
- - Any sequence can be represented as a sum of
scaled, delayed impulses e.g. -
- if pn a-2,a-1, a0, a1,a2, .. then
- pn .. a-2 dn2 a-1dn1 a0 dn a1
dn-1 .. - More generally,
- xn ? xkdn-k
8- Unit Step
- un 1(ngt0) 0 (nlt0)
- un ? dn-k (0ltklt8)
- dn un un-1
- Exponential
- xn Aan
- Sinusoidal
- xn Acos(?0nf) for all n
9Spectra of signals
Frequency spacing ?0 !
10Sampling
- A signal x(t) can be sampled by multiplying it
with a periodic pulse train s(t). Thus the
sampled signal - xs(t) x(t).s(t)
- Being periodic pulse train s(t) can be expanded
using Fourier series, thus - s(t) c0 ?2cn.cos(n?st)
- (cn will be a sinc function thus tend to zero as
t tend to infinity) - if X(f) and x(t) are Fourier Transform pair then
spectra of xs(t) can be obtained using Modulation
theorem, as- - Xs(f) c0X(f) c1X(ffc) X(f-fc)
- c2X(f2fc)X(f-2fc)
.
11If sampling frequency fs is less than 2W (W is
maximum frequency component in X(f)), the
sampling will introduce an unrecoverable
distortion (Aliasing).
Sampling theorem fs gt 2W (fs2w called Nyquist
rate)
12Antialiasing filters
- Ideal filters (instant cutoff) are unrealizable.
- Thus sampling at nyquist rate will result in
aliasing. - The sound of an aliased signal is very annoying
thus either - - x(t) should be first filtered than sampled,
prefiltering antialiasing filters. - x(t) should be first sampled then distorted part
should be filtered out, postfiltering
antialiasing filter.
- Use of antialiase filter degrades voice quality,
but always better than the quality of distorted
signal. - It also filtered out wideband noise, thus
improves quality. - It can be compensated to some extend by
oversampling (it helps in reconstruction) - Telephone standard original voice (0-10Khz)
prealiase filter (0-3.3Khz), oversampling 8Khz.
13Reconstruction
- If sampled signal does not have aliasing then
theoretically same signal can be reconstructed
using ideal filters. - The impulse response of ideal filters is a sinc
pulse. These sinc pulses of various sampling
instants adds up to interpolate the intermediate
portion of signal in time domain.
Ideal filters are unrealizable but as the order
of filter increase its performance improves.
14Uniform Quantization
- dependent variable (e.g. voltage) continuous ?
discrete - (continuous-valued vs. discrete-valued signal)
- It introduces permanent errors, measured in terms
of quantization noise. - x(t) is normalized to - 1 V, sampled and applied
to q-level quantizer. - Step size 2/q
15Quantization noise
- Quantization error ak xq(kTs) - x(kTs)
- ak a uniformly distributed random variable in
the range 1/q to -1/q irrespective of
instantaneous value of x(kTs). - The PDF of ak can be evaluated as q/2. (?(-1/q to
1/q) k dx 1 thus kq/2)
16- Ea ?(-1/q to 1/q) a(q/2) da 0
- Variance Q-noise power (NQ) Ea2 - Ea
- ?(-1/q to 1/q) a2(q/2) da 1/3q2
- In terms of step size S 1/q (-1/q) 2/q
- The quantization noise (NQ) S2/12
17Source Coding of quantized samples
- Fixed length coding (q2n n bits are required
for each level) - Variable length coding
- Comma code
- (each word will start with 0 and one extra 1
at the end. first code 0) - Tree code
- (no code word appears as prefix in another
codeword, first code 0) - Shannon Fano
- ( Bi partitioning till last two elements. 0
in upper/lower part and 1 in lower/upper part)
- Huffman
- (adding two least symbol probabilities and
rearrangement till two elements, back tracing for
code.) - nth extension
- (form a group by combining n consecutive
symbols then code it.) - Lempel Ziv
- (Table formation for compressing binary data)
18Channel Coding
- Systematic generation and addition of redundant
bits to the message bits, at the transmitter,
such that transmission errors can be minimized by
performing an inverse operation, at the receiver. - Example triple repetition codes
- Logic 0 - 000
- Logic 1 - 111
- It has the ability to correct one bit error and
detect two bit error.
19The Concept of frequency
- For Continuous - time sinusoidal signal x(t)
- x(tT) x(t) T is time period 1/F (freq.)
- Increasing F results in an increase in rate of
oscillations (more periods are included in unit
time) - Two sinusoids x1(t) and x2(t) with distinct freq
F1 and F2 are themselves distinct.
20- For Discrete time sinusoids
- x(n) Acos(?nf) ? 2pf
- f (cycles per sample) F/Fs Fs is the sampling
freq - Periodic, if f is a rational number.
- for periodicity, x(nN) x(n) for all n i.e.
- Acos(2pf(nN)f) Acos(2pfnf)
- This relation is true if and only if there exists
an integer K such that 2pfN 2pk i.e. fk/N - If k and N are relatively prime then N is called
the fundamental period of xn. - A small change in freq can results in large
change in period i.e. f1 31/60 implies N160 but
f230/60 implies N2 2.
21- Two sinusoids x1n and x2n whose frequencies
are separated by an integer multiple of 2p are
identical (more precisely, indistinguishable) - Acos((? 2p)nf) Acos(2pn (?nf))
Acos(?nf) - Sinusoids with plt ? lt p i.e. -1/2ltflt1/2 are
distinct. If it appears that a sinusoid has f
outside this range then definitely in the above
range its identical sinusoid does exists.
22- The highest rate of oscillation in a discrete
time sinusoid is attained when ? p (or - p) or
equivalently f1/2 or -1/2. - As ? varies from 0, p/8, p/4, p/2, p then f
?/2p will increases as 0, 1/16, 1/8, ¼, ½ and N
8, 16, 8, 4, 2 - When ? gt p e.g. 3p/2 then f ¾ gt ½ which should
be represented by its identical sinusoid cos(2p-
?)n.
23Classification of systems
- Static (memoryless) Vs Dynamic (with memory)
- (static, if o/p depends upon present i/p only
e.g. y(n) nx(n)bx3(n). Dynamic, if depends upon
past or present e.g. y(n)x(n-1)3x(n)) - Time-invariant Vs Time Variant systems
- (if i/p, x(n) results in y(n) then x(n-k) must
results in y(n-k), otherwise system is time
variant. Ex TI, y(n)x(n)-x(n-1) TV, y(n)
x(n)cos(wn))
24- Linear Vs Nonlinear
- (linear, must satisfy superposition principle,
otherwise nonlinear.) - Causal Vs Noncausal
- (o/p should depend upon present and past i/p but
not on future i/p x(n1).., otherwise
noncausal. If signal is first recorded then
offline processing is done then only non causal
systems are possible to implement) - Stable Vs Unstable
- (BIBO stable if every bounded i/p produces a
bounded o/p, otherwise unstable.)
25Convolution sum
- y(n) Tx(n)
- T represents the operation performed by the
system - If x(n) is an impulse then y(n) h(n) Td(n)
- y(n) T?(k -8 to 8) x(k) d(n-k)
- T..x(-1)d(n1)x(0) d(n)x(1)d(n-1)..
- x(-1)Td(n1)x(0)Td(n)x(1)Td(n-1)..
- By applying superposition property of linear
systems - ? (k -8 to 8) x(k)Td(n-k)
- Apply time-invariance property
- ? (k -8 to 8) x(k)h(n-k) x(n)h(n)
26Graphical evaluation of convolution
- Fold seq h(n) to get h(-k).
- Multiply x(k) with h(-k) and add the results to
get y(0). - Shift h(-k) right to get h(1-k) and repeat step-2
to get y(1). - Repeat above step to get y(2), y(3)..
- Shift h(-k) left to get h(-1-k) and repeat step-2
to get y(-1). - Repeat above step to get y(-2), y(-3)..
27Properties of convolution operation
- Commutative law
- y(n) x(n)h(n) h(n)x(n)
- The role of x(n) and h(n) can be interchanged
(either seq can be folded) - Distributive law
- y(n) x(n)h1(n)h2(n) x(n)h1(n)x(n)h2(n
) - Overall h(n) of two parallel systems equals the
sum of individual h(n)s of subsystems - Associative law
- y(n) x(n)h1(n)h2(n) x(n)h1(n)h2(n)
- x(n)h2(n)h1(n) (using commutative law)
- The overall h(n) is convolution of h(n)s of
cascade subsystems - Order of intermediate subsystems can be
interchanged
28Condition for Causality
- Using convolution sum,
- y(n) ?(k -8 to 8) h(k)x(n-k)
- ?(k 0 to 8) h(k)x(n-k) ?(k -8 to -1)
h(k)x(n-k) - h(0)x(n)h(1)x(n-1).. h(-1)x(n1)
.. - The first sum involves present and past values of
x(n). - The second sum involves future values of inputs
thus all terms of this sum should be zero for
causality. - It is guaranteed if, h(n)0 for nlt0
- It is necessary and sufficient condition for
causality.
29Condition for BIBO Stability
- y(n) ?(k -8 to 8) h(k)x(n-k)
- lt ?(k -8 to 8) h(k)x(n-k)
- (absolute value of sum is always less than or
equal to sum of absolute values) - x(n) is called bounded when there exists a finite
constant Mx such that x(n)ltMxlt8 - If i/p is bounded then,
- y(n)lt?(k -8 to 8)h(k)Mx Mx?(k -8 to 8)
h(k)
30- Applying same definition, y(n) is bounded if
y(n)lt8. - Thus, Mx?(k -8 to 8) h(k)lt8
- As Mx is constant thus sufficient condition for
BIBO stability is Sh ?(k -8 to 8) h(k)lt8 - Is this necessary condition as well?
- (can be proved by showing that there exists a
bounded i/p for which the o/p is unbounded if Sh
is not bounded)
31- If x(n) h(-n)/h(-n) then it will be a bounded
seq for any value of n. (even if h(-n) is
unbounded for some values of n then also x(n)
will be bounded due to division by its absolute
value. ) - Substitute in convolution formula, y(0) for above
sequence will be Sh thus it is unbounded if, Sh
is unbounded. - Sh is bounded if h(n) approaches zero as n
approaches infinity. - If a limited duration x(n) is applied to a BIBO
stable system then it will produce an transient
O/P which will die out to zero at steady state.
(prove it?)
32FIR systems
- This classification depends upon system
characteristics the impulse response. - If h(n)0 for nlt0 and ngtM then causal FIR system
(finite number of symbols in h(n)) - For a causal FIR system
- y(n) ?(k0 to M-1) h(k)x(n-k)
h(0)x(n)h(1)x(n-1) ----- h(M)x(n-M1)
33- To form the o/p at any instant nn0 the most
recent M input samples are required. - System must have a memory size M to remember past
M Input symbols. - This type of realization is called non recursive
realization of systems. - The FIR system acts as a window of size M thus
the impulse response of FIR system is also called
a window
34IIR systems
- The h(n) has infinite number of symbols e.g.
h(n) anu(n). - y(n) ?(k0 to 8) h(k)x(n-k)
h(0)x(n)h(1)x(n-1) h(2)x(n-2) ----- up to 8 - To form the o/p at any instant of time nn0 the
system must remember all the previous input
samples.
35- To store all previous values of I/P, it will
require infinite memory space, which is
practically impossible. - Is it possible to realize an IIR system?
(fortunately , Yes.)
36Recursive systems
- The classification recursive and non recursive
systems is based on method of implementation. - An FIR system can be implemented by either method
but an IIR system can only be implemented by
recursive method. - In recursive systems the output at any instant of
time nn0 will depend upon one or more previous
values of O/Ps. e.g. - y(n) ay(n-1) bx(n)
37Iterative method
- y(0) ay(-1)x(0)
- y(1) ay(0)x(1) a2y(-1)ax(0)x(1)
- y(2) a3y(-1)a2x(0)ax(1)x(2)
- On generalizing,
- y(n) an1y(-1) ?(k0 to n) akx(n-k)
- y(n) yzi(n) yzs(n)
- The O/P of the system can also be represented as
a summation of two terms, one depends upon
applied input yzs(n) and the other upon initial
condition yzi(n).
38Direct method
- Indirect method uses z-transform.
- Obtain characteristic eq by substituting y(n)
(a)n and x(n)0 in the given equation. - Find all roots a1, a2, a3, ----- except a00
- The homogeneous sol will be
- yh(n) c1a1n c2a2n, c3a3n -----------
?(k1toN)ck(ak)n
39- Assume particular solution as per given input. 0
if impulse, ku(n) if u(n), kMn if AMn etc. - The precaution the assumed particular solution
should not exists in homogeneous solution. - Substitute it, along with given input, in given
equation and solve for k by choosing a value of
n such that no term should vanish.
40- y(n) yh(n)yp(n) find values of y(0), y(1),
y(2), .. - Compute y(0), y(1), . from given eq and equate
these values to solve for unknown c1, c2, c3, . - Substitute the values of c1, c2, c3, . and
given initial conditions in y(n) yh(n)yp(n) to
get the total solution.
41Zero input response
- If input is then particular solution will also be
zero. - Now evaluate c1, c2, c3, with y(n) yh(n) and
given equation. - Substitute these values in y(n) yh(n) to get
zero input response.
42Zero state response
- The total solution is y(n) yh(n)yp(n)
- Substitute all initial conditions to be zero and
find c1, c2, c3, - Substitute these values in above equation to get
zero state response
43When O/P depends upon current as well as past
values of I/P
- first the o/p should be calculated only for x(n)
by setting all previous values zero. - Now using the linearity property other o/ps
should be evaluated i.e. y(n-1) replace each n
into n-1 in the y(n) etc. - The overall o/p will be the sum of these o/ps.
- The final expression should always be in terms of
unit impulse or unit step function.
44LCCDE
- y(n) - ?(k1 to N) ak y(n-k) ?(j0 to M) bj
x(n-j) - y(n) yh(n)yp(n) and y(n) yzi(n)yzs(n)
- y(n) an1y(-1) ?(k0 to n) akx(n-k) for 1st
order system - Just observe that at any instant nn0 the O/P
depends upon all previous values of I/P. Thus a
recursive system described by LCCDE is an IIR
system. However the converse is not true. - If all ak are zero, in above LCCDE, it will
describe a non-recursive y(n) ?(j1 to M) bj
x(n-j) thus FIR system.
45- y(n) an1y(-1) ?(k0 to n) akx(n-k) for 1st
order system - The steady state response yss(n) Lim(n?8) y(n)
- yss(n) 0 a is less than 1 1/(1-a) sum of GP
assuming unit step input - yp(n)
- The transient response the part of y(n) which
has been vanished at steady state - Other methods for evaluating steady state and
transient response are also available.
46Impulse response of LCCDE systems
- For 1st order systems yzs(n) ?(k0 to n)
akx(n-k) - Compare it with y(n)?(k -8 to 8) h(k) x(n-k)
- h(n) anu(n) with following conditions
- h(n) must be causal sequence so that lower limit
k0 and k-8 will produce same result. As if klt0
then h(-ve) 0.
47- x(n) should also be causal so that the upper
limit kn and k8 should produce same result. As
if kgtn then x(n-k) x(-ve)0 - The recursive system described by LCCDE is a
causal IIR system with causal input.
48- If i/p is an impulse then particular solution is
zero thus yzs(n) will only depend upon
homogeneous solution of y(n) e.g. for 1st order
system -
- yzs(n) ?(k0 to n) ak?(n-k) ?(k0 to 8)
akx(n-k) anu(n) h(n) - If input is an impulse then
- yzs(n) yh(n) ?(k1toN) ck(ak)n h(n)
49LCCDE system stability
- h(n) ?(k1toN)ck(ak)n
- For stability ?(n0 to 8) h(n) lt 8
- ?(n0 to 8) ?(k1toN)ck(ak)n lt 8
- ?(n0 to 8) ?(k1toN) ck (ak)n lt 8 using
inequality - ?(k1toN)ck ?(n0 to 8) (ak)n lt 8 reversing
order of summation
50- All cks are finite thus above condition is
satisfied if - ?(n0 to 8) (ak)n lt 8 for all k.
- This summation will converge if all aks are less
than unity (infinite GP sum) - Necessary and sufficient condition for stability
of LCCDE systems is, - all roots of characteristic equation should be
less than unity
51Realization of LCCDE systems
- Let we have to realize following system
- y(n) -a1y(n-1)b0x(n)b1x(n-1)
- There are two methods
- Direct form I
- Direct form II (Canonic form)
- Direct form II can be obtained using
associative property i.e. the order of two sub
systems of a system can be interchanged