Summarized by - PowerPoint PPT Presentation

1 / 51
About This Presentation
Title:

Summarized by

Description:

(signals are generated by multiple sensors/sources which are measuring the ... of ground acceleration a few kilometers from epicenter of an earthquake, ECG etc) ... – PowerPoint PPT presentation

Number of Views:28
Avg rating:3.0/5.0
Slides: 52
Provided by: stude613
Category:
Tags: ecg | summarized

less

Transcript and Presenter's Notes

Title: Summarized by


1
Unit - I
  • Summarized by
  • Neetesh Purohit
  • Lecturer, IIIT,
  • Allahabad, UP, India
  • http//profile.iiita.ac.in/np/

2
Classification 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 )

4
Signal 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)

7
Basic 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

9
Spectra of signals
Frequency spacing ?0 !
10
Sampling
  • 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)
    .

11
If 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)
12
Antialiasing 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.

13
Reconstruction
  • 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.
14
Uniform 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

15
Quantization 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

17
Source 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)

18
Channel 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.

19
The 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.

23
Classification 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.)

25
Convolution 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)

26
Graphical 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)..

27
Properties 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

28
Condition 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.

29
Condition 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?)

32
FIR 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

34
IIR 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.)

36
Recursive 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)

37
Iterative 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).

38
Direct 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.

41
Zero 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.

42
Zero 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

43
When 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.

44
LCCDE
  • 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.

46
Impulse 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)

49
LCCDE 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

51
Realization 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
Write a Comment
User Comments (0)
About PowerShow.com