Quantization Error - Finite Word Length Effect Unit 4
Quantization Error - Finite Word Length Effect Unit 4
Fixed point and floating point number representation - ADC - quantization - truncation and rounding -
quantization noise - input / output quantization - coefficient quantization error - product quantization error -
overflow error - limit cycle oscillations due to product quantization and summation - scaling to prevent
overflow.
Finite Word length Effects:
• In the design of FIR Filters, The filter coefficients are determined by the system transfer functions. These
filter co-efficient are quantized/truncated while implementing DSP System because of finite length registers.
• Only Finite numbers of bits are used to perform arithmetic operations. Typical word length is 16 bits, 24
bits, 32 bits etc.
• This finite word length introduces an error which can affect the performance of the DSP system.
• The main errors are
1. Input quantization error
2. Co-efficient quantization error
3. Overflow & round off error (Product Quantization error)
• The effect of error introduced by a signal process depend upon number of factors including the
1. Type of arithmetic
2. Quality of input signal
3. Type of algorithm implemented
1. Input quantization error
• The conversion of continuous-time input signal into digital value produces an error which is known as
input quantization error.
• This error arises due to the representation of the input signal by a fixed number of digits in A/D conversion
process.
2. Co-efficient quantization error
• The filter coefficients are compared to infinite precision. If they are quantized the frequency response of
the resulting filter may differ from the desired frequency response.
i.e poles of the desired filter may change leading to instability.
3. Product Quantization error
• It arises at the output of the multiplier
• When a ‘b’ bit data is multiplied with another ‘b’ bit coefficient the product (‘2b’ bits) should be stored
in ‘b’ bits register. The multiplier Output must be rounded or truncated to ‘b’ bits. This known as overflow
and round off error.
*****************************************************************************************
Types of number representation:
There are two common forms that are used to represent the numbers in a digital or any other digital hardware.
1
MAHARAJA AGRASEN INSTITUTE OF TECHNOLOGY
Department of Electronics and Communication Engineering
0.275 * 2 ➔0.55 ➔0
0.55 * 2 ➔1.10 ➔1
0.10 * 2 ➔0.20 ➔0
0.20 * 2 ➔0.40 ➔0
0.40 * 2 ➔0.80 ➔0
0.80 * 2 ➔1.60 ➔1
0.60 * 2 ➔1.20 ➔1
(30.275)10 = (11110.01000110)2 0.20 * 2 ➔0.40 ➔0
3
MAHARAJA AGRASEN INSTITUTE OF TECHNOLOGY
Department of Electronics and Communication Engineering
4
MAHARAJA AGRASEN INSTITUTE OF TECHNOLOGY
Department of Electronics and Communication Engineering
5
MAHARAJA AGRASEN INSTITUTE OF TECHNOLOGY
Department of Electronics and Communication Engineering
Example
Given X 1 = 3.5 *10 −12 , X 2 = 4.75 *106 . Find the product X 1 X 2
X=(3.5 X 4.75) 10(-12+6)
= (16.625)10-6 → in decimal
In binary: (1.5)10 X (1.25)10 = (210.75) X (210.625)
= 2001 X 0.1100 X 2001 X 0.1010
= 2010 X 0.01111
Addition and subtraction:
• Here the exponent of a smaller number is adjusted until it matches the exponent of a larger number.
• Then, the mantissa are added or subtracted
• The resulting representation is rescaled so that its mantissa lies in the range 0.5 to 1.
• Eg: Add (3.0)10 & (0.125)10
(3.0)10 =2010 X 0.1100 = r e1 X M 1
(0.125)10 =2000 X 0.0010 = r e2 X M 2
Now adjust e2 Such that e1=e2
(0.125)10 =2010 X 0.0000100
Addition →2010 (0.110000 + 0.0000100) →2010 X 0.110010
Subraction →2010 X 1.001101
6
MAHARAJA AGRASEN INSTITUTE OF TECHNOLOGY
Department of Electronics and Communication Engineering
7
MAHARAJA AGRASEN INSTITUTE OF TECHNOLOGY
Department of Electronics and Communication Engineering
8
MAHARAJA AGRASEN INSTITUTE OF TECHNOLOGY
Department of Electronics and Communication Engineering
• Let us assume a sinusoidal signal varying between +1 & -1 having a dynamic range 2
• ADC employs (b+1) bits including sign bit. In this case, the number of levels available for quantizing
x(n) is 2b+1.
• The interval between the successive levels is
2
q = b +1 = 2 −b
2
Where q → quantization step size
-3
If b=3 bits, then q=2 =0.125
***************************************************************************************
Quantization Noise power:
Input Quantization error:
*Derive the equation for quantization noise power (or) Steady state Input Noise Power.
9
MAHARAJA AGRASEN INSTITUTE OF TECHNOLOGY
Department of Electronics and Communication Engineering
2
E e 2 (n ) = e 2 (n ) de
1
q q
−
2
q
2
E e 2 (n ) = e (n )de ------------------------------->(4)
1 2
q q
−
2
E e(n ) = 0
E 2 e(n ) = 0 ------------------------------------------->(5)
Substituting (4) and (5) in (1)
q
2
e (n )de − 0
1 2
q q
e2 =
−
2
q
1 e 3 2
=
q 3 −q
2
1 q q
3 3
= − −
3q 2 2
1 q 3 q 3
= − −
3q 8 8
1 q 3 q 3
= +
3q 8 8
1 2q 3
=
3q 8
q2
e2 = ------------------------------------------------->(6)
12
1
In general, b
= 2 −b = q -------------------------------------------->(7)
2
2
=
(2 )−b 2
e
12
−2 b
2
e2 =
----------------------------------------------->(8)
12
Equation (8) is known as the steady state noise power due to input quantization.
R
q= b → in two’s complement representation.
2
R
q= → in sign magnitude (or) one’s complement representation.
2 −1
b
After quantization, we have noise power e2 as input noise power. Therefore, Output noise power of system is given
by
eo2 = e2 h 2 (n ) ------------------------------------>(9)
n =0
1
2j H (Z )H Z −1
= e2
dZ
z
( )
Where the closed contour integration is evaluated using the method of residue by taking only the poles that lie
inside the unit circle.
Z transform of h(n), H (Z ) = h(n )z − n
n =0
Z transform of h2(n) = Z[h2(n)] = h 2 (n )z − n = h(n )h(n )z − n -------------------->(10)
n =0 n =0
h(n ) = H (Z )Z n −1 dZ -------------------------------->(11)
1
By Inverse Z transform,
2j
Substituting (11) in (10)
h (n )z H (Z )Z n −1 dZ h(n )z − n
1
2 −n
=
n =0 n = 0 2j
=
1
H (Z ) h(n )Z −1 dZ
2j n =0
dZ
h (n) = 2j H (Z ) h(n)Z
2 1 −1
Z −n
n =0 n =0
=
1
2j H (Z ) ( )
h(n ) Z − n
−1
Z −1 dZ
n =0
h (n ) = 2j H (Z )H (Z ) Z ----------------------------------->(12)
2 1 dZ −1
n =0
H(z) =
(1 − a )z H(z −1 ) =
(1 − a )z −1
Given
( z − a) (z −1
− a)
−1
Substituting H(z) and H(z ) in equation (1), we have
e 2 (1 − a )z (1 − a )z −1 −1 2 (1 − a )2 dz
e2 =
2j c ( z − a ) (z − a )
−1
z dz = e
2j c ( z − a ) (z −1 − a )
z −1
2 1
= e residue of H(z) H(z −1 ) z −1at z = a + residue of H(z) H(z −1 ) z −1at z =
a
= e ( z − a )
2 (1 − a ) z −1
2
+
(z − a )(z −1 − a )
0
(1 − a )2 2 (1 − a )
= e −1 = e
2
(z − a )
(1 + a )
2 −2 b
Where, e =
2
12
**************************************************************************************
Find the steady state variance of the noise in the output due to quantization of input for the first order
filter. [Apr/May’11] [Nov/Dec-2016]
y ( n) = ay ( n − 1) + x( n)
Solution:
The impulse response for the above filter is given by h(n) = a n u (n)
2 = e2 h 2 (n)
k =0
= 2
e a
k =0
2n
= 1 + a 2 + a 4 + ....
2
e
1
= e2
1 − a2
−2 b
2 1
= (or )
12 1 − a 2
1 z −1
2 = e2
2 j
c ( z − a)( z −1 − a )
dz
z −1
residue of −1
at z = a
2 ( z − a)( z − a)
= e
z −1
+ residue of −1
at z = 1 / a
( z − a)( z − a)
z −1
= e2 ( z − a)
( z − a )( z −1 − a ) z = a
a −1 1
= e2 −1
= e2
a −a 1 − a2
***************************************************************************************
The output of the A/D converter is applied to a digital filter with the system function
H (Z ) =
0.45Z
Z − 0.72
Find the output noise power of the digital filter, when the input signal is quantized to 7 bits.
Given:
H (Z ) =
0.45Z
Z − 0.72
Solution:
0.45Z −1
H (Z )H Z ( )Z −1 −1
=
0.45Z
−1
Z − 0.72 Z − 0.72
Z −1
0.45 2 Z −1
=
(Z − 0.72) 1 − 0.72
Z
0.2025Z −1
=
(Z − 0.72) 1 − 0.72Z
Z
0.2025Z −1 Z
=
(Z − 0.72) Z − 1
0.72
− 0.28125
=
(Z − 0.72)(Z − 1.3889)
Now the poles of H(Z)H(Z-1)Z-1 are p1=0.72 , p2=1.3889
Output noise power due to input quantization
1
eo2 = e2 ( )
H (Z )H Z −1 Z −1 dZ
2j
( )
N
= e2 Re s H (Z )H Z −1 Z −1
i =1 z = pi
( )
N
= e2 Re s H (Z )H Z −1 Z −1
i =1 z = pi
Where p1,p2,…..pn are the poles of H(Z)H(Z-1) Z-1 that lies inside the unit circle in z-plane.
− 0.28125
eo2 = e2 (Z − 0.72)
(Z − 0.72)(Z − 1.3889) Z =0.72
− 0.28125
= e2
0.72 − 1.3889
= 0.4205 e2
***************************************************************************************
1 1
Consider the transfer function H ( z ) = H 1 ( z ) H 2 ( z ) where H 1 ( z ) = −1
and H 2 ( z ) =
1 − a1 z 1 − a 2 z −1
Find the output round off noise power. Assume 1 = 0.5 and 2 = 0.6 and find output round off noise
power.
Solution:
The round off noise model for H ( z ) = H 1 ( z ) H 2 ( z ) is given by,
From the realization we can find that the noise transfer function seen by noise source e1(n) is H(z), where,
1
H ( z) =
(1 − a1 z )(1 − a2 z −1 ) − − − − − − − − − − − − − − − − − − − − − − − (1)
−1
0 2 = 012 + 02 2 − − − − − − − − − − − − − − − − − − − − − − − − − − − −(3)
1
012 =
2j c H ( z ) H ( z −1 ) z −1dz
2 1 1 1 1 1
= e
2j c 1 − a1 z 1 − a 2 z 1 − a1 z 1 − a 2 z
−1 −1
z −1dz
2 1 1
= e of residue of H ( z ) H ( z −1 ) z −1 at po1es z = a1 , z = a 2 , z = and z =
a1 a2
If a1 and a2 are less than the poles z=1/a1 and z=1/a2 lies outside of the circle z = 1. So, the residue of H(z) H(z-1) z-
1
at z=1/a1 and z=1/a2 are zero. Consequently we have,
012 = of residue of H(z)H(z −1 ) z −1 at po1es z = a1 , z = a 2
z −1 z −1
= ( z − a1 ) + ( − )
(1 − a1 z −1 )(1 − a2 z −1 )(1 − a1 z )(1 − a2 z ) z=a1 z a 2
(1 − a1 z −1 )(1 − a2 z −1 )(1 − a1 z )(1 − a2 z ) z=a2
2 1 1
= e +
a2
2
( ) a2
1 − 1 − a 2 (1 − a1a 2 ) 1 − (1 − a1a 2 ) 1 − a 2
2
( )
a1 a1
a1 1 1 a2 1 1
012 = e 2 . . + . . − − − − − − − − − − − − − −( 4)
a1 − a 2 1 − a1
2
1 − a 1 a 2 a 2 − a1 1 − a 2
2
1 − a1 a 2
In the same way,
e2
02 2 =
2j c H 2 ( z ) H 2 ( z −1 ) z −1dz
e2 1 1
=
2j c 1 − a2 z 1 − a2 z
−1
z −1dz
z −1
= e ( z − a2 )
(1 − a2 z −1 )(1 − a2 z ) z=a2
2
z −1
= e (z − a 2 z −1 )
(1 − a2 z )(1 − a2 z ) z=a2
2
−1
2 1
= e 2
− − − − − − − − − − − − − − − − − − − − − − − − − − − − − − − − − − − − − −(5)
1 − a 2
1 a1 1 1 a2 1 1
02 = e2 + . . + . .
1 − a 2
2
a1 − a 2 1 − a1 1 − a1a 2 a 2 − a1 1 − a 2 1 − a1a 2
2 2
1
= e +
(
a1 1 − a 2 − a 2 1 − a1
2 2 2
) ( )
( )( )
2
1 − a2 1 − a1 1 − a 2 (1 − a1a 2 )(a1 − a 2 )
2 2 2
2
= e
1
+
(a1 − a2 )(1 + a1a2 )
1 − a 2
2
( )( )
1 − a1 1 − a 2 (1 − a1a 2 )(a1 − a 2 )
2 2
=
1
2 −2 b
+
(1 + a1a2 )
12
1 − a 2
2
( )(
1 − a1 1 − a 2 (1 − a1a 2 )
2 2
)
The steady state noise power for a1 = 0.5, a 2 = 0.6 is given by
2 −2 b 1 1 + (0.5)(0.6)
= +
12 1 − (0.6) 2
( )(
1 − (0.5) 1 − (0.6) (1 − 0.6(0.5))
2 2
)
2 −2 b
= (5.4315)
12
*************************************************************************************
1
Draw the quantization noise model for a second order system H ( z ) = and find the
1 − 2r cos z −1 + r 2 z −2
steady state output noise variance.
Solution:
Given:
1
H ( z) =
1 − 2r cos z −1 + r 2 z −2
The quantization noise model is,
we know, 0 = 01 + 02
2 2 2
2 −2 b
sin 2 ( n + 1)
0 2 = 2.
12
r 2n
n =0 sin 2
−2 b
=
2 −2 b 1 (1 + r )2 (1 − cos 2 )
6 2 sin 2 (1 − r 2 )(1 − 2r 2 cos 2 + r 4 )
=
2 −2 b (1 + r ) 2
6 (1 − r )(1 − 2r 2 cos 2 + r 4 )
2
******************************************************************************
Co-efficient quantization error
• We know that the IIR Filter is characterized by the system function
M
b z k
−k
H (Z ) = k =0
N
1 + ak z −k
k =1
• After quantizing ,
M
b k q z −k
H (Z )q = k =0
N
1 + a k q z − k
k =1
Where ak q = a k + a k
bk q = bk + bk
• The quantization of filter coefficients alters the positions of the poles and zeros in z-plane.
1. If the poles of desired filter lie close to the unit circle, then the quantized filter poles may lie outside
the unit circle leading into instability of filter.
2. Deviation in poles and zeros also lead to deviation in frequency response.
***************************************************************************************
Consider a second order IIR filter with H ( z ) = 1.0 find the effect on quantization
(1 − 0.5 z −1 )(1 − 0.45 z −1 )
on pole locations of the given system function in direct form and in cascade form. Take b=3bits.
[Apr/May-10] [Nov/Dec-11]
Solution:
Given that,
1.0
H ( z) =
(1 − 0.5 z )(1 − 0.45 z −1 )
−1
1
H ( z) = −1 −1
z ( z − 0.5 z ) z −1 ( z − 0.5)
z2
=
( z − 0.5)( z − 0.45)
The roots of the denominator of H(z) are the original poles of H(z). let the original poles of H(z) be p1 and
p2.
Here p1=0.5 and p2=0.45
Direct form I:
1.0
H ( z) =
(1 − 0.5 z )(1 − 0.45 z −1 )
−1
1
H ( z) = −1
1 − 0.5 z − 0.45 z −1 + 0.225 z − 2
1
= −1
1 − 0.95 z + 0.225 z −2
Let us quantize the coefficients by truncation.
Convert to Truncate to Convert to
.9510 .11112 .1112 .87510
Binary 3-bits decimal
Convert to Truncate to Convert to
.22510 .00112 .0012 .12510
Binary 3-bits decimal
Let H ( z ) be the transfer function of the IIR system after quantizing the coefficients.
1
H ( z) = −1
1 − 0.875 z + 0.125 z −2
Y ( z) 1
let H ( z ) = = −1
X ( z ) 1 − 0.875 z + 0.125 z −2
On cross multiplying the above equation we get,
Y ( z ) − 0.875 z −1Y ( z ) + 0.125 z −2Y ( z ) = X ( z )
Y ( z ) = X ( z ) + 0.875 z −1Y ( z ) − 0.125 z −2Y ( z )
Cascade form:
Given that
1.0
H ( z) = −1
(1 − 0.5 z )(1 − 0.45 z −1 )
In cascade realization the system can be realized as cascade of first order sections.
H(z)=H1(z)+H2(z)
Where,
1 1
H1 (z) = −1
and H 2 (z) =
1 − 0.5 z 1 − 0.45 z −1
Let us quantize the coefficients of H1(z) and H2(z) by truncation.
Convert to Truncate to Convert to
.510 .10002 .1002 .510
Binary 3-bits decimal
Convert to Convert to Convert to
.4510 .01112 .0112 .37510
Binary 3-bits decimal
let , H1 (z) and H 2 (z) be the transfer function of the first-order sections after quantizing the coefficients.
1
H1 (z) =
1 − 0.5 z −1
1
H 2 (z) =
1 − 0.375 z −1
Y (z) 1
let , H1 (z) = 1 =
X (z) 1 − 0.5 z −1
Y1 (z) − 0.5 z −1 Y1 (z) = X(z)
Y1 (z) = X (z) + 0.5 z −1 Y1 (z)
Y (z) 1
let, H 2 (z) = =
Y1 (z) 1 − 0.375 z −1
on cross multiplying the above equation we get ,
Y (z) − 0.375 z −1 Y (z) = Y1 (z)
Y (z) = Y1 (z) + 0.375 z −1Y(z)
***************************************************************************************
Round off effects and overflow in digital filter:
*Explain in detail about round off effects in digital filters.
• The presence of one or more quantizer in the realization of a digital filter results in a non-linear device.
i.e. recursive digital filter may exhibit undesirable oscillations in its output
• In the finite arithmetic operations, some registers may overflow if the input signal level becomes large.
• These overflow represents non-linear distortion leading to limit cycle oscillations
• There are two types of limit cycle oscillations which includes
1. Zero input limit cycle oscillations (Low amplitude compared to overflow limit cycle oscillations)
2. Over flow limit cycle oscillations.
Zero input limit cycle oscillations
• The arithmetic operations produces oscillations even when the input is zero or some non zero constant
values. Such oscillations are called zero input limit cycle oscillations.
Overflow limit cycle oscillations
• The limit cycle occurs due to the overflow of adder is known as overflow limit cycle oscillations.
Dead Band:
The limit cycle occurs as a result of quantization effect in multiplication. The amplitude of the output
during a limit cycle is confined to a range of values called the dead band of the filter.
2 −b
y(n − 1) 2
(1 − a )
Consider a first order filter
y (n ) = ay(n − 1) + x(n ); n >0
After rounding the product
y q (n ) = Qa y(n − 1) + x(n );
The round off error
2 −b 2 −b
− er
2 2
where, er →difference between the quantized value and the actual value.
2 −b
Qay(n − 1) − ay(n − 1)
2
The dead band of the filter for the limit cycle oscillations are
y (n − 1) a 0
Qay(n − 1) =
− y (n − 1) a 0
2−b
y(n − 1) − a y(n − 1)
2
2−b
y(n − 1) (1 − a )
2
2 −b
Dead band of the filter, y(n − 1) 2
(1 − a )
**************************************************************************************
Problem: Consider a 1st order FIR system equation y(n) = x(n) + ay(n − 1) with
0.875 , n = 0
x ( n) =
0 , otherwise
Find the limit cycle effect and the dead band. Assume b=4 and a=0.95. (Nov/Dec-12)(Nov/Dec-15)
[May/June-2016]
Solution:
Given:
0.875 , n = 0
x(n) =
0 , otherwise
−b
2 2 −4
Dead band = = = 0.625
2(1 − a ) 2(1 − 0.95 )
y (n) = x(n) + 0.95 y (n − 1)
Qay(n − 1)
n x(n) y (n − 1) ay (n − 1) y ( n ) = x( n ) + Qay (n − 1)
(round off to 4-bits)
0 0.875 0 0 0.0000 y(0)=0.875
0.875 * 0.95
= (0.1101)2
1 0 0.875 = (0.83125)10 y(1)=0.8125
= 0.8125
= (0.11010)2
0.8125 * 0.95
= (0.1100)2
2 0 0.8125 = (0.77187)10 y (2) = 0.75
= 0.75
= (0.110001)2
0.75 * 0.95
= (0.1011)2
3 0 0.75 = (0.7125)10 y (3) = 0.6875
= 0.6875
= (0.1011011)2
0.6875 * 0.95
= (0.1010)2
4 0 0.6875 = (0.653125)10 y (4) = 0.625
= 0.625
= (0.101001)2
0.625 * 0.95
= (0.1010)2
5 0 0.625 = (0.59375)10 y (5) = 0.625
= 0.625
= (0.10011)2
0.625 * 0.95
= (0.1010)2
6 0 0.625 = (0.59375)10 y (6) = 0.625
= 0.625
= (0.10011)2
Conclusion:
The dead band of the filter is 0.625. When n 5 the output remains constant at 0.625 causing
limit cycle oscillations.
***************************************************************************************
Overflow Limit cycle oscillations:
*What are called overflow oscillations? How it can be prevented?
• We know that the limit cycle oscillation is caused by rounding the result of multiplication.
• The limit cycle occurs due to the overflow of adder is known as overflow limit cycle oscillations.\
• Several types of limit cycle oscillations are caused by addition, which makes the filter output oscilate
between maximum and minimum amplitudes.
• Let us consider 2 positive numbers n1 & n2
n1=0.111→7/8
n2=0.110→6/8
n1 + n2=1.101→-5/8 in sign magnitude form.
The sum is wrongly interpreted as a negative number.
• The transfer characteristics of an saturation adder is shown in fig below
where n → The input to the adder
f(n) → The corresponding output
• Saturation arithmetic eliminates limit cycles due to overflow, but it causes undeniable signal distortion
due to the non linearity of the clipper.
• In order to limit the amount of non linear distortion, it is important to scale input signal and unit sample
response between input and any internal summing node in the system to avoid overflow.
From figure
W ( z) S0 S
H ' ( z) = = −1 −2
= 0
X ( z ) 1 + a1 z + a2 z D( z )
S0 X ( z )
W ( z) = = S0 S( z ) x(Z)
D( z )
1
Where S ( z ) =
D( z )
we have
S
w(n) = 0 S (e j )X (e j )(e jn ) d
2
S2 2
w(n) 2 = 0 2 S (e j )X (e j )(e jn ) d
2
Using Schwartz inequality
w(n) 2 S02 2 S (e j ) d 2 X (e j ) d
2 2
n =0 2 2 S (e j ) d
if z = e then dz = je j d
j
which gives
dz
d =
jz
1
w(n) 2 S02 x 2 (n)
2
2 j
c S ( z ) z −1d z
n =0
1
w(n) 2 S02 x 2 (n)
2 j
c S ( z ) S ( z −1 ) z −1d z
n =0
w2 (n) x 2 (n) when
n =0
1
2 j
S02 c S ( z ) S ( z −1 ) d z = 1
1
=
1 z −1dz
2 j D( z ) D( z −1 )
c
1
S02 =
I
Where I=
Note:
• Because of the process of scaling, the overflow is eliminated. Here so is the scaling factor for the first
stage.
1
• Scaling factor for the second stage = S01 and it is given by S 012 = 2
S0 I 2
1 H 1 (Z )H 1 (Z −1 )Z −1
2j c D2 (Z )D2 (Z −1 )
Where I 2 = dZ
********************************************************************************************
0.25 + 0.7 Z −1
For the given transfer function, H (Z ) = , find scaling factor so as to avoid overflow in
1 − 0.5Z −1
the adder ‘1’ of the filter.
Given:
D(Z) = 1-0.5Z-1
D(Z-1) = 1-0.5 Z
Solution:
1 1 dZ
I =
2j c D(Z )D Z −1
( )
Z
1 1 dZ
= ( )
2j c 1 − 0.5Z (1 − 0.5Z ) Z
−1
1 Z 1 dZ
=
2j c (Z − 0.5) (1 − 0.5Z ) Z
Z 1
Residue of +0
(Z − 0.5) (1 − 0.5Z ) Z =0.5
I=1.3333
1
S0 =
I
1
S0 =
1.333
= 0.866
***************************************************************************************
Consider the recursive filter shown in fig. The input x(n) has a range of values of ±100V, represented
by 8 bits. Compute the variance of output due to A/D conversion process. (6)
Solution:
Given the range is ±100V
The difference equation of the system is given by y ( n) = 0.8 y ( n − 1) + x( n) , whose impulse response h(n)
can be obtained as
h(n) = (0.8) n u (n)
range of the signal
quantization step size =
No. of quantization levels
200
= 8
2
= 0.78125
Variance of the error signal
q2
e2 =
12
(0.78125) 2
=
12
e = 0.05086
2
Variance of output
2 = e2 h 2 (n)
n =0
= (0.05086) (0.8) 2 n
n =0
0.05086
= = 0.14128
1 − (0.8) 2
***************************************************************************************
The input to the system y(n)=0.999y(n-1)+x(n) is applied to an ADC. What is the power produced by
the quantization noise at the output of the filter if the input is quantized to a) 8 bits b)16 bits. May-07
Solution:
y(n)=0.999y(n-1)+x(n)
Taking z-transform on both sides
Y(z)=0.999z-1Y(z)+X(z)
Y ( z) 1
H ( z) = =
X ( z ) 1 − 0.999 z −1
z z −1
H ( z ) H ( z −1 ) z −1 = ( )( −1 ) z −1
z − 0.999 z − 0.999
z −1
=
1
( z − 0.999)(−0.999)( z − )
0.999
−0.001
=
( z − 0.999)( z − 0.001)
Where p1,p2,……pN are poles of H(z)H(z-1)z-1, that lies inside the unit circle in z-plane.
0.001
eoi2 = e2 ( z − 0.999)( )
( z − 0.999)( z − 0.001) z =0.999
= e2 500.25
b) b+1=16 bits
22(15)
2 = (500.25) = 3.882 10−8
12
***************************************************************************************
Find the effect of coefficient quantization on pole locations of the given second order IIR system, when
it is realized in direct form I and in cascade form. Assume a word length of 4 bits through truncation.
1
H ( z) = −1
1 − 0.9 z + 0.2 z −2
Solution:
Direct form I
Let b=4 bits including a sign bit
(0.9)10 = (0.111001...) 2
Integer part
0.9 2
1.8
1
0.8 2
1.6
1
0.6 2
1
1.2
0.2 2
0.4
0
0.4 2
0.8
0
0.8 2
1.6
1
z1=0.695
z2=0.178
If we compare the Poles of H(z) and H ( z ) we can observe that the poles of H ( z ) deviate very much from
the original poles.
Cascade form
1
H ( z) = −1
1 − 0.5 z (1 − 0.4 z −1 )
(0.5)10 = (0.1000) 2
After truncation we get
(0.100)2=(0.5)10
After truncation we get
(0.011)2=(0.375)10
0.4 2
(0.4)10 =
0.8
0
0.8 2
1.6
0
0.6 2
1.2
1
0.2 2
0.4
1
0.4 2
0.8
0
(0.4)10=(0.01100…)2
The system function after coefficient quantization is
1
H ( z) = −1
(1 − 0.5 z )(1 − 0.375 z −1 )
N
eoi
2
= e2 Res H ( z ) H ( z −1 ) z −1
i =1 z = pi
N
eoi
2
= e2 ( z − pi ) H ( z ) H ( z −1 ) z −1
i =1 z = pi
−0.0331
eoi
2
= e2 ( z − 0.68) *
( z − 0.68)( z − 1.4706) z = 0.68
−0.0331
eoi
2
= e2 * = 0.0419 e2
(0.68 − 1.4706)
eoi
2
= 0.0419*1.2716*10−4
eoi
2
= 5.328*10−6
1. Sampling
• Sampling is the conversion of a continuous- tome signal into a discrete-time signal obtained by taking the
samples of continuous-time signal at discrete instants.
• Thus if xa(t) is the input to the sampler, the output is xa(nT)≡x(n), where T is called the sampling interval.
2. Quantisation
• The process of converting a discrete-time continuous amplitude signal into digital signal is called
quantization.
• The value of each signal sample is represented by a value selected from a finite set of possible values.
• The difference between the unquantised sample x(n) and the quantized output xq(n) is called the
quantization error or quantization noise.
eq(n)= xq(n)-x(n)
• To eliminate the excess bits either discard them by the process of truncation or discard them by rounding
the resulting number by the process of rounding.
• The values allowed in the digital signals are called the quantization levels
• The distance ∆ between two successive quantization levels is called the quantization step size or resolution.
• The quality of the output of the A/D converter is measured by the signal-to-quantization noise ratio.
3. Coding
• In the coding process, each discrete value xq(n) is represented by a b-bit binary sequence.