US4071842A

US4071842A - Apparatus for analog to digital conversion

Info

Publication number: US4071842A
Application number: US05/756,399
Authority: US
Inventors: Stuart Keene Tewksbury
Original assignee: Bell Telephone Laboratories Inc
Current assignee: AT&T Corp
Priority date: 1975-08-28
Filing date: 1977-01-03
Publication date: 1978-01-31
Anticipated expiration: 1995-01-31
Also published as: JPS5227345A; US4107669A; US4017849A; BE845505A; FR2322487A1; FR2322487B1; DE2638534A1; CA1112366A; NL7609570A; JPS6056335B2; GB1548044A; IT1066933B

Abstract

An input analog signal is converted to a digital output signal by an oversampled predictive DPCM coder which includes an n stage delay line in the feedback loop. The n delay line outputs are weighted by coefficients a_i . . . a_n selected according to the relationship ##EQU1## AND THEN SUMMED. Alternatively, the feedback loop may comprise a chain of n integrators arranged so that the signal fed back to the comparator is the sum of single, double, triple...n order integration. By so doing, the coder attenuates the signal power at the quantizer input while the attenuator coefficients are independent of the input signal statistics.

A similar technique may also be applied in an oversampled error feedback coder, which includes a feedback loop having an n stage delay line. Here again, the delay line outputs are weighted in accordance with the above relationship. Alternatively, a chain of n integrators may be used in the coder input, and an identical chain employed in the feedback loop. By so doing, the coder attenuates the coding noise power in the signal band while the coder design is rendered independent of the quantizing error statistics.

Description

This application is a division of application Ser. No. 608,524 filed Aug. 29, 1975.

BACKGROUND OF THE INVENTION

1. Field of the Invention

The present invention relates generally to conversion of analog signals to digital signals and, more particularly to both predictive and error feedback coders that sample the analog signal at generally high rates and provide a differential pulse code modulation (DPCM) output.

2. Description of the Prior Art

Sophisticated digital integrated circuit technologies have made digital signal processors which perform a variety of functions a practical alternative to conventional analog systems. The trend toward digital processing is accompanied by many advantages, including precisely predictable performance, the ability to share hardware over several separate channels without interchannel distortion, and the flexibility of programming one processor to perform a variety of functions.

Since digital processors must typically interface with continuous signals, analog/digital and digital/analog converters thus become increasingly important system components. Many such converters have been developed, including predictive coders which reduce the dynamic range of the input signal applied to the quantizer within the coder, or error feedback coders which shape the spectral distribution of the quantizing error so as to reduce in-band noise. Both types of coders may make use of tapped delay lines in the feedback path, with the n outputs of the delay line being fed through n separate attenuators having coefficients a₁ . . . a_n, and then summed. With this arrangement, various techniques are available for choosing appropriate attenuator coefficients. In each technique, it is desirable to minimize the number of required quantizaton levels, so that the number of bits needed to adequately represent the signal being encoded is reduced. Certain of the prior art techniques for coefficient optimization have been developed primarily for transmission systems. However, in such systems, oversampling, as desired in the present invention, cannot be used, since it results in an undesirably high transmission rate. Accordingly, one subject of the present invention is to optimize the coefficient selection technique in cases where oversampling is acceptable.

Other prior art techniques known in the field of video transmission have utilized moderate oversampling in conjunction with single or double integration feedback networks. However, these systems have also been limited in the degree of quantizer simplification achieved. Accordingly, another object of the present invention is to retain the signal independent feedback network as in video coders, while achieving a much greater reduction in the number of required quantization levels.

SUMMARY OF THE INVENTION

Each of the foergoing and additional objects are achieved in accordance with the principles of the instant invention by providing, in both predictive and error feedback oversampled coders which include delay lines having n output stages, a series of n attenuators having coefficients selected according to the relationship ##EQU2## Alternatively, in the case of predictive coders, the feedback loop may include a chain of n integrators arranged so that the signal fed back to the comparator is the sum of single, double, triple . . . n order integration. In the case of error feedback coders, serial chains of n integrators may be employed both in the coder input line, and in the feedback loop. By virtue of the advantageous coder arrangement of the present invention, the number of levels required in the quantizer is significantly reduced, while the feedback network is maintained independent of the input signal.

BRIEF DESCRIPTION OF THE DRAWING

The invention will be more clearly understood from a consideration of the following detailed description when read in light of the accompanying drawing in which:

FIG. 1 is a block diagram of a prior art predictive coder;

FIG. 2 is a block diagram of one realization of the feedback filter of FIG. 1;

FIG. 3 is a graph of both the mean square error signal and the reduction in the number of required quantizer levels in the coder of FIGS. 1 and 2, when the coefficients in accordance with the invention are used;

FIG. 4 is a block diagram of a predictive coder constructed in accordance with the principles of the present invention;

FIG. 5 is an alternative arrangement for the coder of FIG. 4;

FIG. 6 is a graph of the transfer function of the coders of FIGS. 4 and 5 as a function of frequency;

FIG. 7 is a block diagram of a prior art error feedback coder;

FIG. 8 is a block diagram of an error feedback coder constructed in accordance with the principles of the instant invention;

FIG. 9 is a block diagram of a third order integrator that may be used in the coder of FIG. 8; and

FIG. 10 is a graph showing the reduction in the number of quantization levels achievable using the coder of FIG. 8.

DETAILED DESCRIPTION

As shown in FIG. 1, a prior art predictive coder typically includes a forward loop containing an analog to digital (A/D) converter such as quantizer 101 and a feedback loop which contains a feedback filter 102, to be discussed in more detail hereinafter. The output of feedback filter 102 is fed back to its input via line 103 and additive circuit 104, which combines the output with the output of quantizer 101. The feedback filter 102 output is also supplied to subtractive circuit 105, thus completng the feedback circuit to the quantizer.

In operation, the quantizer 101 input signal e_n is the difference between the input signal s_n and the output f_n of filter 102. Quantizer 101, assumed to have equal quantization levels, has an output E_n which approximates e_n by the nearest discrete quantization level, i.e.,

E.sub.n = e.sub.n + q.sub.n                                (1)

where q_n is the quantization error. From the above relationships, it is seen that the output of circuit 104 is given by:

X.sub.n = E.sub.n + f.sub.n                                (2)

and the output of circuit 105 is given by:

e.sub.n = s.sub.n - f.sub.n                                (3)

so that, from

equations

1, 2 and 3:

X.sub.n = (e.sub.n +f.sub.n) + q.sub.n = s.sub.n +q.sub.n. (4)

If the characteristis of feedback filter 102 are appropriately chosen, its output f_n is a good prediction of the next input signal sample s_n. As a result, e_n which represents the extrapolation error of the coder, is quantized in quantizer 101, instead of the input signal s_n, as would be the case in an A/D converter not using predictive feedback. Specifically, if s_p and e_p represent the peak values of s_n and e_n, respectively, and it is assumed that s_p and e_p are proportional to their mean square values σ₅ ² and σ_e ², the number of quantization levels required in quantizer 101 is reduced by a factor

θ = √σ.sup.2.sub.e / √σ.sup.2.sub.s (5)

relative to the number required in nonpredictive coders.

One known manner of arranging filter 102 to produce the above results is the use of the delay line having a series of n stages such as

stages

201, 202, 203 and 204 as shown in FIG. 2. connected to the output of each stage 201 through 204 are n attenuators such as

attenuators

211, 212, 213 and 214 each having a coefficient a₁, a₂, a₃ . . . a_n. The outputs of attenuators 211 through 214 are summed in

additive circuits

220, 221, 222 so that the Z-transform K(Z) of the entire filter is given by: ##EQU3##

When the delay line arrangement of FIG. 2 is used in the coder of FIG. 1, the mean square value of th quantizer input signal, neglecting quantizing errors, is given by: ##EQU4## wherein ψ_m, the autocorrelation function of the input signal, is given by: ##EQU5## The standard prior art procedure for minimizng σ² _e is to simultaneously solve the N linear equations:

∂σ.sup.2.sub.e /∂ a.sub.i = 0 for i = 1.2 . . . N .                                                   (9)

while solution of equation (9) yields the filter coefficients a_i, the coefficients are functions of and depend upon th input sample correlation ψ_e. However, such dependence is not desirable in many applications, since the input signals are nonstationary. Accordingly, an alternative filter design approach is needed.

In accordance with the present invention, signal independent attenuator coefficients a_i are computed according to the relationship ##EQU6## This relationship is arrived at by first considering that the input sample correlations ψ_e and the two-sided power spectral density G_s (f) of the input signal are related by: ##EQU7## where f_s is the sample rate and f_o is the input signal bandwidth. Using equation (11) in equation (7) σ² _e may be expanded as a power series in (f_o /f_s)², since it is assumed that f_s >>2f_o due to oversampling. Such expansion yields the general form ##EQU8## where the β_l depend only on the filter coefficients and N. By choosing the filter coefficients in accordance with equation (10), β_l = 0 for l ≦ N-1 in equation (12), leaving the lowest order nonvanishing term in equation (12) as ##EQU9## Combining equations (13) and (5), it is seen that by using the coefficients selected in accordance with equation (10), the number of quantization levels required is significantly reduced. For example, as shown in FIG. 3, for a three stage delay line (N=3) and for a sampling rate f_s of 20 times that of the input signal bandwidth f_o, the number of quantization levels required in quantizer 101 is reduced by more than a factor of 64, as compared to the case of a quantizer without feedback. For the same ratio of f_s /f_o, the reduction of levels for 4 and 5 stage delay lines are greater than 256 and 512, respectively. Alternatively stated, it is seen from FIG. 3 that by the advantageous choice of attenuator coefficients in accordance with the present invention, σ² _e decreases by 6N dB per octave increase in sampling rate, for an N stage delay line; quantization levels are accordingly decreased by a factor of 2N for each doubling of the sampling rate.

Using the results of equation (10) in equation (6), it will be seen that, in Z transform notation, the transfer characteristic of filter 102 is given by:

K(Z) =  1-(1-Z.sup.-1).sup.N.                              (14)

also, it will be seen that the transfer function T(Z) between the output and input of filter 102, taking into account feedback line 103, is: ##EQU10## Substituting equation (14) in (15), the output/input transfer function is thus: ##EQU11## Since equation (16) may be written in the form ##EQU12## and since the transfer function of a single integrator is given by 1/(1-Z^-1), a coder using analog prediction rather than the predictive filter of FIG. 2, may be constructed in accordance with the invention as shown in FIG. 4. As before, the forward loop contains a quantizer 401 which is arranged to sample the analog input signal applied thereto at a rate f_s much greater than the input signal bandwidth f_o, and to provide at its output a digital signal E_n, representative of said input signal. The digital signal is fed, in the feedback loop, first to a digital to analog converter 402, and then to a serial chain of M analog integrators, such as

integrators

410, 411, and 412. The outputs of

integrators

410, 411, and 412 are connected via

lines

420, 421 and 422 to an adding circuit 430, such that the adder output consists of the sum of single, double, triple . . . M order integration. The feedback loop is completed by connecting the output of adding circuit 430, which output is a good prediction of the input signal, to the negative input terminal of a comparator or differencing circuit 440, the positive input to which is the analog input signal.

Inspection of the circuit of FIG. 4 reveals that the transfer function between the output of converter 402 and the output of adder 430 is given, as desired, by equation (17). The initial Z^-1 term is built into the time response of the integrators. The use of D/A converter 402 prior to the integrators 410 through 412 in the feedback loop of course requires that the integrators be designed to operate on analog signals. This arrangement is, however, preferred over the use of digital integrators followed by a digital to analog conversion, since a very high degree of converter precision would then be needed. However, when the coder is configured as shown in FIG. 4, converter 402 may be constructed with a degree of precision only matching that of quantizer 401, which as explained previously, is advantageously quite simple.

An alternative arrangement to the coder of FIG. 4 is shown in FIG. 5. In this embodiment, quantizer 501, converter 502, and differencing circuit 540 are all equivalent to their FIG. 4 counterparts.

Adder circuits

550 and 551 are disposed after each

integrator

510 and 511 in the serial chain of M integrators, with the exception only of the last integrator 512. The output of converter 502 is simultaneously applied to the input of the first integrator 510 in the chain, as well as to each

adder circuit

550 and 551. As will be seen by inspection, the output of the last integrator 512 is the sum of first, second, third . . . M order integration, as was the case for the coder of FIG. 4.

The advantages of the predictive coders of FIGS. 4 and 5 may be still better appreciated by reference to the graph of FIG. 6, in which the coder transfer function (the ratio of the coder output to input, in dB is plotted as a function of the radio f/f_s of the input signal frequency to the sampling frequency. For a second order coder in accordance with the invention, i.e., a coder having a feedback loop in which the output is the sum of single plus double integration, curve 601 shows the coder transfer function, which is given by:

[h(f)].sup.2 = 4(1-cos2πf/f.sub.s).sup.2 .              (18)

By way of comparison, curve 602 shows the transfer function of a prior art coder having two serially connected integrators in the feedback loop. The transfer function in this case is given by: ##EQU13##

As can be seen from the graph, superior attenuation is provided by the coders in accordance with the invention, in the ranges where f/f_s ≲0.2, i.e., where there is oversampling. This is true notwithstanding the use, in prior art apparatus, of intentionally "leaky" integrators, which tend to reduce but not eliminate entirely the peak at the transfer function pole.

The invention described above with respect to predictive coders is also applicable to error feedback coders, such as prior art coders of the type depicted in FIG. 7. As shown, the forward loop includes an A/D converter or quantizer 701, and the feedback loop contains a feedback filter 702. However, in this arrangement, the input to quantizer 701, which is a composite signal x_n comprising the difference between the input signal s_n and a feedback signal f_n, is subtracted from the output thereof by a differencing circuit 704, to form a quantizing error signal q_n. The error signal q_n is filtered by filter 702 to form the aforesaid feedback signal, which is then subtracted from the input signal in a second differencing circuit 705. In this arrangement, assuming that filter 702 is an N tap delay line with coefficients a₁ . . . a_n, the quantizer 701 output signal X_n is the nearest discrete quantization level to the quantizer input sample x_n, wherein ##EQU14## and

X.sub.n = x.sub.n + q.sub.n .                              (21)

Since, as stated previously, equation (6) is the transfer function of filter 702, equation (20) may be written in Z - transform rotation as

X(Z) = S(Z) + [1-H(Z)]·q(Z) .                     (22)

equation (22) shows that the spectral density of the quantizing noise is shaped by the response of the Z-transform 1-H(Z), and thus that appropriate choice of the coefficients a_i, assuming f_s >>f_o, will minimize the amount of in-band noise power by placing most of the noise out-of-band where it may be eliminated using a digital low pass filter. By so doing, the quantization level separation can be increased relative to that of an A/D converter without feedback; this, in turn, permits a desired decrease in the number of quantization levels required.

Minimization of noise in the band between -f_o and f_o is generally accomplished, using prior art techniques, by first recognizing that the noise power Q(f_o) is given by ##EQU15## where G_q (f) is the two-sided power density function of the individual quantizing errors q_n, which, it is assumed, are not correlated with the input samples, and where ##EQU16## Expanding equation (23) using equation (24) gives ##EQU17## where ##EQU18## Filter coefficients which minimize Q(f_o) may then be obtained by solving the simultaneous equations

∂Q(f.sub.o)/∂a.sub.k = 0 K = 1,2 . . . N, (27)

which yields a solution

B =  Γ.sup.-1 . γ                              (28)

where B and γ are vectors with components a_k and β_k, respectively, and Γ is an NXM matrix with components

Γ.sub.j,k = β.sub.j-k .                         (29)

This solution, however, is often undesirable in that the coefficients obtained from equation (29) depend on the noise spectral density G_q (f).

In accordance with the invention, spectrum independent coefficients are obtained in a manner similar to that described above for predictive coders, namely by expanding equation (25) as a power series in (f_o /f_s)² and choosing the a_k such that coefficients of terms up to the N^th term vanish. The resulting coefficients are again given by ##EQU19## and the transfer function of an N-tap delay line with coefficients according to equation (30) is again given by equation (14).

Realization of an error feedback pulse code modulation encoder utilizing analog integrators in lieu of an N tap delay line in the feedback loop is accomplished as shown in FIG. 8. As in FIG. 7, the forward loop includes a quantizer 801 which is arranged to sample the composite signal X(Z) applied thereto at a rate f_s greater than the Nyquist rate 2f_o. However, the input signal is applied to quantizer 801 via a first N^th order integrator 806, followed by a combinatorial circuit 805. Circuit 805 is arranged to algebraically combine (i.e., add or subtract in a desired manner) the input signals applied thereto. In the feedback loop, the output of quantizer 801 is applied to a digital to analog converter 802, which may be constructed with the same simplicity as quantizer 801, to form a direct feedback signal f_n. This signal is applied via line 804 to a positive input of circuit 805, and also to the input of a second N^th order integrator 803, the output of which is applied to a negative input of circuit 805. N^th order integrators 803 and 806 are simply serial chains of N analog integrators. For example, a third order integrator would simply include

integrators

901, 902 and 903, as shown in FIG. 9.

The equivalence between the coders of FIGS. 7 and 8 when the former advantageously utilizes the coefficients of equation (30) may be appreciated by rewriting equation (22) as

X(Z) = S(Z) + G(Z) ·q(Z)                          (31)

wherein

G(Z) .tbd. 1-H(Z) .                                        (32)

bearing in mind that the transfer function of a multistage delay line with coefficients in accordance with equation (30) is

H(Z) = 1-(1-Z.sup.-1).sup.N,                               (14)

then

G(Z) = (1-Z.sup.-1).sup.N .

accordingly, the transfer function of N^th order integrators 803 amnd 806 is given by 1/G(Z), since, as stated previously, the transfer function of an individual integrator is 1/(1-Z^-1). Applying the foreging to the coder FIG. 8, it will be seen that for an input signal S(Z), and an output signal X(Z), the output of integrator 806 is given by S(Z)/G(Z), the output of integrator 803 is given by X(Z)/G(Z), and the signal on line 804 is the output signal X(Z). The output of circuit 805 is thus ##EQU20## and the output of quantizer 801 is thus ##EQU21## Equation (34) when simpified, results in equation (31), as desired.

The advantages of the coder of FIG. 8 can be appreciated by reference to FIG. 10, which shows the reduction in the number of quantization levels achievable by use of the present invention. For a sampling rate f_s of 20 times the highest signal frequency f_o, and for a three stage integrator, a reduction in the number of levels by a factor of 2⁸ is attained.

While the above descriptions of predictive and error feedback coders indicate that performance improves as the number N of delay line stages (or integrators) increases, this theoretical improvement is limited by experimental implementations which indicate that use of feedback networks in which N>5 is not practical. First, the loop gain G is nominally unity but is bounded, for stability reasons by

2.sup.N /(2.sup.N +1) < G < 2.sup.N /(2.sup.N -1)          (35)

for an N-tap predictor. Thus, as N gets larger, the stability limits are narrowed. Furthermore, loop delay is added by the nonzero propagation delays of the A/D and D/A converters and must be compensated by, for example, feedforward techniques in the cascaded integrators of FIG. 5. Neglecting such compensation can also lead to coder instability. Therefore, the simplification of the quantizer is offset by an increasing sensitivity of the coder stability to the analog feedback network as N increases.

In addition to increased stability problems as N becomes large, the quantizer input signal power contains a contribution, neglected above, from the nonzero quantizing errors. Letting σ² _e be the mean square power of the predictive coder's difference signal e_n, σ² _s be the mean square power of the error feedback coder's input signal power and assuming a uniform level separation δ in the quantizer (random quantizing errors and noise power δ² /12), the mean square value σ_in ² of the quantizer input signal is ##EQU22## for predictive coders and ##EQU23## for error feedback coders. In predictive coders, σ² _e decreases with N increasing while δ₂ remains fixed, setting a lower limit on δ_in ² as N increases. With error feedback coders, δ increases with increasing N while σ_s ² remains fixed, again placing a lower limit on σ_in ².

For these reasons, a practical limit on feedback filter complexity in an oversampled coder realization is likely to be N ≲ 5.

An example of the advantages of the present invention is illustrated by the A/D interface unit which is presently employed to quantize the frequency division multiplexed signals in the digital processor described in Communication Technology, Vol. COM-19, No. 6, December 1971 pp. 1050-59. This unit presently quantizes a 108 KHz bandwidth signal using a 13 bit uniform PCM code. However, utilizing a 4 stage predictive coder constructed in accordance with the invention, with a 2 MHz sample rate, only an 8-level (3 bit) quantizer is required.

While not shown in the drawing, remote decoding of the outputs of both the predictive and error feedback coders of the present invention is quite simply attained. For predictive encoders, the output/input relationship for the entire coder is given by

O(Z) = (1-Z.sup.-1).sup.N [I(Z) + q(Z)],                   (38)

where O(Z) is the output signal, I(Z) is the input signal and q(Z) is the quantizing error. Accordingly, the remote decoder may comprise a serial chain of N integrators, which, as stated previously, has a transfer function given by 1/(1-Z^-1)^N. In the case of error feedback encoders, the output/input relationship is given by

O(Z) = I(Z) + (I-Z.sup.- 1).sup.N q(Z) .                   (39)

accordingly, a low pass filter with a stop band chosen to compensate for the increasing noise spectral density at higher frequencies, may conveniently be used for decoding.

Various modifications and adaptations of the present invention will be readily apparent to those skilled in the art. For this reason, it is intended the invention be limited only by the appended claims.

Claims

What is claimed is:

1. An error feedback pulse code modulation encoder for converting an analog input signal to a digital output signal representative thereof, comprising

a first differencing circuit for combining said input signal and an error feedback signal to generate a composite signal,

quantizing means responsive to said first differencing circuit for generating said digital output signal by sampling said composite signal and converting said samples to digital form,

a second differencing circuit for combining said composite signal and said digital output signal to generate a quantizing error signal, and

means for generating said error feedback signal from said quantizing error signal including

1. a delay line having stages 1 through n for receiving said quantizing error signal, each of said stages having an output tap,

2. n attenuators having attenuation coefficients a₁ through a_n connected to respective ones of said output taps, and

3. means for algebraically summing the outputs of said attenuators to form said error feedback signal,

wherein said attenuation coefficients are selected in accordance with the relationship ##EQU24##

2. The invention defined in claim 1 wherein said analog signal has a bandwidth f_o and said quantizer is arranged to sample said error signal at a rate of at least 5 f_o.

3. An error feedback pulse code modulation encoder arranged to convert an analog input signal to a digital output signal representative thereof, comprising

a quantizer for sampling a composite analog signal applied thereto at a rate f_s to generate said digital output signal,

means for generating a quantization error signal from the difference between digital output signal and said composite signal,

a filter for receiving said quantization error signal and for generating a feedback signal therefrom, and

means for generating said composite signal from the difference between said analog input signal and said feedback signal,

wherein said filter has a transfer function K(Z) in the Z domain given by

K(Z) = 1-(1-Z.sup.-1).sup.N

and N is an integer between 2 and 5.