1
Estimation of Delta–Sigma Converter Spectrum
E. Nunzi, P. Carbone, D. Petri,
Dipartimento di Ingegneria Elettronica e dell’Informazione,
Università degli Studi di Perugia,
via G. Duranti 93 – 06125 Perugia, Italy.
Phone: ++39 075 5853634, Fax: ++39 075 5853654, Email: nunzi@diei.unipg.it.
Abstract – Effects of the windowing process,
widely investigated by the scientific literature for
narrow–band components embedded in white noise,
is not sufficiently detailed when signals are corrupted by colored noise. Such a phenomenon can
heavily affect the spectral parameters estimation of
the noisy signal. In this paper effects of the windowing on the output of analog–to–digital converters with ∆Σ topology, which present a spectrally
shaped quantization noise, is analyzed. In particular, the spectral leakage of both narrow– and wide–
band components is investigated and a criterion for
choosing the most appropriate window for any given
modulator resolution is given. The proposed analysis validates the use of the Hanning sequence as the
optimum two term cosine window to be employed
for characterizing low order ∆Σ modulators.
cation systems and the estimation of their spectral
parameters is usually carried out in the frequency–
domain by weighting the output data with the Hanning sequence.
In this paper, the spectral leakage effects on the
∆Σ modulator output are considered and a criterion for choosing the most appropriate window for
a given modulator resolution is given. Such analysis validates the use of the Hanning sequence as
the most suitable window to be employed for accurately estimating parameters of both narrow– and
wide–band components. Moreover, it is analyzed
the effect of the windowing process on the power
spectrum of the modulator quantization error.
Keywords – Delta–Sigma spectrum, spectral
leakage, Delta–Sigma characterization.
The architecture of a generic ∆Σ modulator is
composed by an L–th order loop filter, a quantizer
and a negative feedback loop. In spite of its apparent
simple structure and because of the presence of the
quantizer, the analysis of its performance is rather
complex. In order to simplify this task, the internal
ADC is often modeled as an additive error signal e[·]
uniformly distributed and uncorrelated with the input signal x[·]. By assuming such a linear model, the
output of an L–th order modulator can be written as:
I. I NTRODUCTION
Estimation of spectral figures of merit of analog–
to–digital converters (ADCs), such as signal–
to–random noise ratio (SRN R), signal–to–noise
and distortion ratio (SIN AD ), spurious–free dynamic range (SF DR) or total harmonic distortion (T HD ), is usually carried out by employing
frequency–domain based techniques. In particular,
as described in the standards IEEE 1057 and 1241,
such spectral parameters are calculated from the
spectrum of the converter output sequence. This is
usually estimated using the Discrete Fourier Transform (DFT) based on N acquired samples, that allow the power evaluation of the narrow– an wide–
band components [1] [2].
When non–coherent sampling applies, the finite
number of processed samples induces spectral leakage phenomena which may affect the estimation of
the parameter of interest. In order to reduce such effects, windowing is usually applied to the acquired
data [1]–[4]. While this technique has been widely
investigated for classical Nyquist–rate ADCs, the
analysis of the windowed output spectrum of ∆Σ
modulators, which present a spectrally shaped quantization error sequence, is not sufficiently detailed
by the published scientific literature.
Nevertheless, such converters are commonly employed in digital measurement and telecommuni-
II. ∆Σ OUTPUT
SPECTRUM
yL [n] = x[n − L] + qL [n]
(1)
where qL [·] is the quantization error of the L–th order shaped modulator.
Estimates of SRN R, SIN AD or T HD can be
carried out by stimulating the modulator with a sinusoidal signal and by applying the FFT algorithm to
R data records each of length N . The finite length of
the processed signal and the discretization of the frequency axis, induce an estimation bias on the ADC
output spectrum.
In order to reduce the spectral leakage phenomenon, the modulator output samples are usually multiplied by an appropriate window, w[·]. Sequences commonly employed are those belonging
to the cosine–class because they can easily be calculated. In this paper, w[·] has been normalized to
the square root of the energy value in order to bound
the maximum of the correlation sequence values to
1. By indicating with w[·] the normalized window,
0
0
L=0
−50
maxH (dB)
PqL (f )/P
−50
RECT
L=0
qwL (f ) (dB)
P
−100
L=1
−150
L=2
−200
L=3
−250
L=4
−100
L=1
−150
L=2
−200
HANNING
L=3
−250
−300
(a)
−350
L=4
−300
(a)
−350
−4
10
−3
−2
10
10
−1
−4
10
10
−3
−2
10
f
10
−1
10
f
0
0
−50
maxH (dB)
PqL (f )/P
−50
RECT
L=0
−100
qwL (f ) (dB)
P
−100
−150
L=1
L=2
−200
L=3
−150
−200
L=4
−250
−300
(b)
−350
HANNING
L=0
L=1
L=2
−250
L=3
−300
L=4
(b)
−350
−4
10
−3
−2
10
10
−1
10
f
Power spectral densities of the Hanning–
windowed output of an 8–bit modulator stimulated by a
sinewave with amplitude equal to the internal ADC full–
scale under coherent (a) and non–coherent (b) sampling.
Lines in each figure refer to five different order-shaping L
as indicated by the corresponding labels. The input frequency value has been set to the modulator upper–band
edge by considering an OSR = 8.
Figure 1.
it follows that the analyzed signal, y Lw [n], is:
−4
10
−3
−2
10
10
−1
10
f
Figure 2. Solid lines represent the power spectral densities of the 1–bit (a) and 8–bit (b) L–order–shaped quantization noise, normalized to PmaxH . Bolded lines represent
the sidelobe envelope of the windows indicated by the corresponding labels, and centered on the modulator upper–
band edge when considering an OSR = 8.
The frequency axis has been normalized to the converter sampling rate. The input frequency value has
△
yLw [n] = yL [n]w[n] = (x[n − L] + qL [n])w[n], (2) been set equal to f 0 = 1/(2 · OSR), i.e. on the
modulator upper–band edge. Moreover, the modun = 0, ..., N − 1.
lator output has been windowed by the Hanning sequence, and the periodogram has been carried out on
Simulation results presented in this paper, refer to
R = 50 non–overlapped data records each of length
a modulator output weighted by the Hanning winN = 216 . Fig. 1(a) and (b) refer to coherent and
dow. This is the sequence most commonly emnon–coherent sampling of the input sinewave freployed for estimating spectral parameters of ∆Σ
quency, respectively, both normalized to the maxmodulators, because of its suitable properties in the
imum of the periodogram under the coherent confrequency domain.
△
In particular, such a window presents a high sidedition, i.e. PmaxH = F S 2 /(4 · EN BW ), where
2
lobe envelope decay together with a small mainlobe
−1 2
△
N −1
EN BW = N N
is the
width (i.e. 4 frequency bins) thus optimizing the ren=0 w [n]/
n=0 w[n]
quirements of frequency selectivity and low spectral
equivalent noise bandwidth of the window.
leakage. As a consequence, a high resolution quanIn both figures, the power spectral densities of
tizer is needed in order to evaluate effects of spechigh–order modulators converge to a constant value
tral leakage associated with the narrow–band comin the low–frequency band. When non–coherent
ponent.
sampling applies, such a phenomenon is mainly due
Fig. 1 shows the power spectral densities, estito the spectral leakage of the narrow–band compomated by means of the periodogram, PqLw (f ), of
nent. Thus, in the low–frequency region of Fig. 1(b),
samples of the window sidelobes are graphed. Howan 8–bit ∆Σ converter with an oversampling ratio
ever, when coherency is guaranteed, since zeros of
(OSR) equal to 8, when stimulated by a full–scale
the window spectrum are located on integer fre(F S ) input sinewave, for different modulator order–
quency bin values, such a phenomenon is not due
shaping, as indicated by the corresponding labels.
to the narrow–band component, but only to spectral
leaking of the shaped–quantization noise.
Because of the discrete resolution of the frequency axis, which is equal to 1/N , the periodogram shows such a phenomenon only by employing a high number of samples N.
In the following, the effects of the spectral leakage
of the narrow–band component and of the shaped–
quantization noise are analyzed and discussed separately.
III. S PECTRAL LEAKAGE OF THE
NARROW– BAND COMPONENT
The power spectral density of the shaped noise
qL [·] can easily be calculated by applying the Discrete Time Fourier Transform (DTFT) to its autocorrelation sequence, thus obtaining [6]:
PqL (f ) = σe2 22L sin2L (πf ),
|f | < 0.5,
(3)
where σe2 is the quantization error power of the internal ADC and f is the normalized frequency. The
behavior of (3), normalized to PmaxH , is shown in
Fig. 2(a) and (b) (solid lines) for a 1–bit and an 8–
bit modulator, respectively, for different orders L, as
indicated by the corresponding labels.
To determine the effects of spectral leakage of
the narrow–band component on the modulator output spectrum, the envelope of the sidelobes of the
rectangular and Hanning windows have been also
plotted in Fig. 2(a) and (b) with bolded lines for
N = 216 . The envelopes have been traced only for
frequencies lower than that of the input tone since
effects of the window sidelobes at higher frequencies can be neglected because of the higher power
of the shaped noise. Moreover, the image component of the narrow–band signal can be neglected,
thus providing a simpler analysis.
As a consequence, the behavior of the rectangular and Hanning windows can be approximately described by 1/|f0 N − f N | and 1/|f0 N − f N |3 , respectively [5].
The spectral leakage of the narrow–band component can affect the modulator output spectrum only
if the shaped–noise power spectral density is lower
than the envelope of the used window, especially
when high OSRs are employed. As an example,
Fig. 2(b) shows that, by applying the Hanning window, spectral leakage effects appear only for modulator orders higher than or equal to 2.
For low frequency values, the sidelobe envelopes
of the rectangular and Hanning windows converge,
respectively, to 1/(f0 N ) and to 1/(f0 N )3 . By setting these values equal to (3), the frequency f x over
which the leakage of the narrow–band component
do no affect the spectral estimation, can be easily
determined. In particular, since for small f the condition sin2L (πfx ) ≃ (πfx )2L holds true, for any
given L–th order shaping modulator, windowed by a
rectangular or a Hanning window of length N , such
frequency values can be expressed, respectively, as:
fxH = K 3 · W
fxR = K · W,
(4)
1/L
where K = (1/f0 N )1/2L and W = 1/(2πσe ).
Thus, the spectral leakage of the narrow–band
component can be made arbitrary small by decreasing the quantizer resolution (i.e. reducing σ e ) or the
modulator shaping order or by increasing the number of processed samples N .
Notice that, since the first and the last samples of
the Hanning window are equal to 0, the Hanning
sequence is the two term window which presents
a sidelobe decay equal to 1/f 3 instead of 1/f [5].
Thus, it is the most suitable two term window to
employ for reducing spectral leakage effects on the
modulator output spectrum, at least for order shaping lower than 2. Higher order modulators require
windows with more coefficients.
IV. E FFECT OF WINDOWING ON THE
SHAPED – QUANTIZATION NOISE
The power spectral densities of the ∆Σ shaped
noise, windowed by the rectangular and the Hanning sequence, have been calculated as indicated in
App. A. In particular, the expressions of the power
spectral densities of the noise filtered by a first and
second order modulator, and windowed by a rectangular sequence, result to be respectively equal to:
Pq1 wR (f ) = Pq1 (f ) + 2σe2 cos(2πf )/N , |f | < 0.5
2σe2 (4 cos(2πf )
Pq2 wR (f ) = Pq2 (f ) +
−2 cos(4πf ))/N,
(5)
+
|f | < 0.5.
(6)
Expressions (5) and (6) show that the estimated
power spectral densities of the first and second order
shaped noise converge, in the low frequency region,
to a constant value approximately equal to 2σ e2 /N
and 4σe2 /N , respectively.
By applying the same procedure, the expressions
of power spectral densities of the filtered noises,
windowed by the Hanning sequence, are:
Pq1 wH (f ) = Pq1 (f ) + 2σe2 (1+
−RwH [1] cos(2πf )) , |f | < 0.5
(7)
2σe2
(RwH [2] cos(4πf )+
Pq2 wH (f ) = Pq2 (f ) +
−4RwH [1] cos(2πf ) + 3) , |f | < 0.5
(8)
where Rw H [1] and RwH [2] are constant values that
can be calculated from (A.6).
As in the previous case, (7) and (8) show that for
low frequency values, the power spectral densities
of the noise shaped by a first and second order modulator, converge, respectively, to constant values approximately equal to
KH1 = 2σe2 (1 − RwH [1])
KH2 = 2σe2 (RwH [2] − 4RwH [1] + 3)
(9)
(10)
which can be made arbitrarily small by increasing
the number of acquired samples N .
By substituting in (9) and (10), the values of σ e2
and N employed in Fig. 1, we have K H1 = −148.5
dB and KH2 = −249.5 dB.
Thus, the leakage of the wide–band components,
due to the truncation of the modulator output sequence, affects the estimation of the power of the
modulator output error when a low N is employed,
especially when high OSRs are considered.
V. C ONCLUSIONS
In this paper the windowing process on ∆Σ modulator output has been analyzed. The effect of
the spectral leakage of the narrow–band component and of the spectrally shaped–quantization noise
on the estimated power spectral density has been
discussed. In particular, it has been shown that
the spectral leakage of the narrow– and wide–band
component can be made arbitrarily small by increasing the number of processed samples and by properly choosing a window sequence.
Windows most commonly employed are those
attaining to the two term cosine class since they
present a small mainlobe width and are easy to calculate. The presented analysis has shown that the
Hanning sequence is the two term window which
guarantees the lower estimation bias of spectral parameters, thus validating its use in the characterization of ∆Σ modulators.
A PPENDIX A
D ERIVATION
OF EXPRESSIONS
(5)–(8)
By indicating with F{·} the Fourier Transform
operator, the power spectral density of a windowed
△
L–th order shaped noise q Lw [n] = qL [n]w[n] with
autocorrelation function RqLw [m] is:
PqL w (f ) = F {RqL w [m]} = F {RqL [m]Rw [m]} (A.1)
where RqL [m] represent the autocorrelation function
of the L–th order shaped noise and
Rw [m] =
N
−1
w[n + m]w[n],
(A.2)
n=0
m = −(N − 1), ..., (N − 1)
is the aperiodic correlation sequence of the employed normalized window.
The autocorrelation sequence of a white noise
with zero–mean and variance equal to σ e2 , filtered
by a first and second order modulator, respectively,
result to be [6]:
Rq1 [m] = σe2 (2δ[m] − δ[m − 1] − δ[m + 1])
(A.3)
σe2
Rq2 [m] = (6δ[m] − 4δ[m − 1] − 4δ[m + 1]+
+δ[m − 2] + δ[m + 2])
(A.4)
where δ[·] is the discrete Dirac pulse.
By applying (A.2), the autocorrelation of the normalized rectangular and Hanning windows result to
be, respectively:
|m|
, m = −(N − 1), ..., (N − 1) (A.5)
N
2
1
2
2π
|m| +
cos
|m| (N − |m|) +
RwH [m] = −
3
3N
3N
N
RwR [m] = 1 −
2
+
3N
sin
2π
cos 2π
N
2π −
N2π
1 − cos
N
2 sin
N
m = −(N − 1), ..., (N − 1)
sin
2π|m|
N
,
(A.6)
By substituting (A.3) and (A.4) in (A.1) and considering that for real windows Rw [m] = Rw [−m],
expressions (5)–(8) result.
R EFERENCES
[1] Standard for Digitizing Waveform Recorders, IEEE
Std. 1057, Dec. 1994.
[2] Standard for Terminology and Test Methods for
Analog-to-Digital Converters, IEEE Std 1241, Oct.
2000.
[3] P.Carbone, E.Nunzi, D.Petri, “Windows for ADC
Dynamic Testing via Frequency–Domain Analysis,”
IEEE Trans. Instrum. and Meas. Tech., pp. 1679–
1683, Dec. 2001.
[4] I.Kollár, “Evaluation of Sinewave Tests of ADC’s
from Windowed Data,” Proc. 4–th Intern. Workshop
ADC Modelling and Testing, pp. 64–68, Bordeaux,
France, Sept. 9–10, 1999.
[5] A.H.Nuttall, “Some Windows with Very Good Sidelobe Behavior,” IEEE Trans. on Acoustics, Speech,
and Signal Processing, pp. 84–91, Feb. 1981.
[6] S. R. Norsworthy, R. Schreier, G. C. Temes , “Oversampling Delta-Sigma Data Converters : Theory,
Design, and Simulation”, IEEE Press., 1997.