The invention relates to high speed analog-to-digital converters (ADC) and, more particularly, to digital equalization of analog-to-digital conversion systems with an ADC that consists of a plurality of time interleaved sub-ADCs.
An increase in the sampling rate of analog-to-digital conversion may be achieved by the use of composite ADCs. A composite ADC contains a number of interleaved sub-ADCs with a common input and a sequential timing. If the number of sub-ADCs equals N, then the resulting conversion rate is N times larger than the rate of one sub-ADC.
Each sub-ADC incorporated in a composite ADC has its own amplitude frequency response and phase frequency response. The misalignment of amplitude and phase frequency responses of different sub-ADCs causes specific signal distortions, with the appearance of spurious frequency components being of prime importance.
The main way to prevent the appearance of the specific distortions in a composite ADC is to use equalization of its output digital signal. There are several patents concerned with digital equalization of a composite ADC output signal, for example U.S. Pat. No. 5,239,299, U.S. Pat. No. 7,408,495, U.S. Patent Application Publication Nos. US 2005/0151679, US 2010/0182174, and others. The equalizer in these patents is an FIR filter (or a set of FIR filters), with the samples coming from each of sub-ADCs being corrected with equalizer coefficients that are calculated from the frequency responses of this sub-ADC.
The ADC equalizer is built usually as a conventional Finite Impulse Response (FIR) filter. The most resource consuming components of FIR filter are multipliers. As the equalizer length L may reach several hundreds of taps, the required number of multipliers becomes the main reason that makes it necessary to use in the equalizer design, more FPGAs and/or FPGAs of bigger size.
It is well known in the art that there is a need for reducing the number of multipliers in FIR filters by using more effective algorithms. There are different approaches to the solution of this problem. One that is most successful and most suitable for digital equalizer design was developed in the works of S. Winograd, Z. Mou and P. Duhamel.
The digital equalizer, like any FIR filter, forms its output sample by calculating a convolution between the input samples and equalizer coefficients. The Winograd-Mou-Duhamel algorithm reduces the number of multiplication in convolution calculation by using, at each clock cycle, some intermediate calculation results obtained in the preceding cycle.
A digital equalizer for an interleaved ADC is a time variant device. The equalizer coefficients that are used at a current clock cycle depend on responses of the correspondent sub-ADC and are different from coefficients used in the previous cycle. For this reason, the Winograd-Mou-Duhamel algorithm, as it is, cannot be directly used to reduce the multipliers number in an ADC digital equalizer.
The present technology provides a digital equalizer for an interleaved ADC that performs equalization of the frequency responses with a reduced number of multipliers.
The foregoing and other objects, features and advantages will be apparent from the following more particular description of the embodiments, as illustrated in the accompanying drawings in which like reference characters refer to the same parts throughout the different views. The drawings are not necessarily to scale, emphasis instead being placed upon illustrating the principles of the embodiments.
A digital equalizer for correction of the frequency responses of an interleaved ADC is described by an assembly of coefficients sets Hm[i]. Here, m is the number of a set and i, 0≦i<L, is the number of a coefficient in the set, L being the length of the equalizer. The operation of a equalizer at any arbitrary time instant is controlled by the set Hm[i], where the number m of the set coincides with the number of the sub-ADC that produced the sample coming at this instant to the equalizer input. The coefficients set Hm[i] for a specific m, is determined from a measured frequency responses of the sub-ADC with the number m. If x[n] are samples at the input of an equalizer, then its output sample y[n] is formed as a convolution of L last input samples and the coefficients Hm[i]:
The sequence of samples x[n] produced by an interleaved ADC has an inherent regularity: the samples that are spaced by M clock cycles (M being the number of sub-ADCs in the interleaved ADC) come from the same sub-ADC. It means that at time instants spaced by M clock cycles the operation of the equalizer is controlled by the same set of coefficients. This fact creates the opportunity to reduce the number of multipliers in the equalizer: when calculating the equalizer output sample y[n] according to the equation (1), it is possible to use some intermediate results that were obtained M clock cycles before during the calculation of the sample y[n−M].
The digital equalizer for correction of the frequency responses of an interleaved ADC according to the present technology, is built with the use of parallel processing: the stream of digital samples created by an interleaved ADC and having samples rate F, is split into N samples sub streams with samples rate F/N each, the internal component units of equalizer handling the samples sub streams simultaneously. The samples rates ratio N is a parameter of the design and determines, in particular, the complexity of the resulting hardware. By increasing the parameter N, the factor of multipliers reduction may be increased.
As it will be clear from following disclosure, the internal structure of the equalizer makes it important that the parameter N equals two to the power of v, where v is an arbitrary integer.
A block diagram of an interleaved ADC with a digital equalizer according to the present technology is shown in
The interleaved ADC 1 converts input analog signal into a stream of digital samples x[n]. This stream comes to the input of demultiplexer 2. The demultiplexer 2 splits the samples stream x[n] into N sub streams, the sample rate of each sub stream being reduced by factor N. N samples appear at the N outputs of the demultiplexer 2 at each period of the reduced samples rate; when the number of a period equals k, these N samples form a group (x[Nk], x[Nk+1], x[Nk+2], . . . , x[Nk+N−1]).
An example of a sequence of input N-groups (N-group being a group of N samples) for N=4 is shown in
N-groups from the outputs of the demultiplexer 2 are applied to the inputs of the samples repositioning unit 3. The samples repositioning unit 3 transposes samples, collecting in one N-group, samples that were produced by the same sub-ADC (see for example
There are discontinuities in the sequence of repositioned N-groups at the outputs of the samples repositioning unit 3. These discontinuities appear repeatedly with an interval of M N-groups. An example of discontinuity may be seen in the
The sequence of input samples x[n] may by broken into segments with a length of NM samples, with the segment with the number r comprising samples with numbers r·NM≦n≦(r+1)·NM−1. The samples repositioning unit 3, while processing the samples from the segment with the number r, forms M N-groups with numbers m, 0≦m≦M−1. All four samples of a quad with the number k belong to a sub stream of input samples that are produced by a sub-ADC with a number m. The lagging N-group with a number m, contains samples
(x[r·NM+m], x[r·NM+m+M], x[r·NM+m+2M], . . . , x[r·NM+m+(N−1)M]),
while the leading N-group with a number m, contains samples
(x[r·NM+m−(N−1)M], x[r·NM+m−(N−2)M], x[r·NM+m−(N−3)M], . . . , x[r·NM+m])
that are ahead of the lagging quad samples by (N−1)M.
(In the example shown in
(x[r·4M+m], x[r·4M+m+M], x[r·4M+m+2M], x[r·4M+m+3M]),
while the leading quad, with a number m, contains samples
(x[r·4M+m−3M], x[r·4M+m−2M], x[r·4M+m−M], x[r·4M+m]),
which are ahead of the lagging quad samples by 3M).
The leading N-group and the lagging N-group pass from the outputs of the samples repositioning unit 3 to the inputs of corresponding PreFIRs transformers 4 and 9. A PreFIRs transformer converts N-group into a set of K samples, with K being equal to three raised to the power of v.
A PreFIRs transformer for N=2 converts two input samples (a[1], a[2]) into three output samples (b[1], b[2], b[3]) in compliance with the next equations:
b[1]=a[1];
b[2]=a[1]+a[2];
b[3]=a[2].
A PreFIRs transformer for N>2 is constructed by an iterative procedure. The procedure is carried out step by step, with the iterative parameter n being equal to two at the first step, and being doubled at each transition to the next step. At each step, parameter k is supposed to correspond to the parameter n (with n being equal to two to the power of v, and k being equal to three to the power v, with v being the same arbitrary integer in both cases). A PreFIRs transformer of order 2n is constructed from two PreFIRs transformers of order n, and k PreFIRs transformers of order 2, in accordance with the next instructions (an order of PreFIRs transformer coincides with the number of its inputs):
Using the instructions presented above, the next relationship is determined between the input and output samples of the PreFIRs transformer of order 4:
b[1]=a[1];
b[2]=a[1]+a[3];
b[3]=a[3];
b[4]=a[3]+a[4];
b[5]=a[1]+a[2]+a[3]+a[4];
b[6]=a[1]+a[2];
b[7]=a[2];
b[8]=a[2]+a[4];
b[9]=a[4].
A block diagram of a PreFIRs transformer for the case of N=4 is shown in the
The sets of K samples from the outputs of PreFIRs transformers 4 and 9 go to the correspondent inputs of the assembly of K double buffer FIR filters 5. The connections inside the assembly of double buffer FIR filters 5 are shown in the
A block diagram of a double buffer FIR filter is shown in the
Each double buffer FIR filter with a number 1, 1≦1≦K possesses coefficients h[1, m, i] that are applied to the corresponding multipliers of this double buffers FIR filter. Here m, 0≦m<M, is the number of sub-ADC that produced current input sample, and i, 0≦i<L/4, is the ordinal number of the coefficient in the set. These coefficients are calculated starting from the specified equalizer coefficients Hm[i] according to an iterative procedure defined by the collection of statements:
h[0, m, i]=Hm[i],
h[1, m, i]=Hm[i]+Hm[i+2M], and
h[2, m, i]=Hm[i+2M],
h2[3r, m, i]=h1[r, m, i], 0≦r<k, 0≦i<L/(2n);
h2[3r+1, m, i]=h1[r, m, i]+h1[r, m, i+M], 0≦r<k, 0≦i<L/(2n);
h2[3r+2, m, i]=h1[r, m, i+M], 0≦r<k, 0≦i<L/(2n).
As an example, the described procedure was applied to a double buffer FIR filter with the parameters N=4 and K=9. The resulting sets of coefficients follow the equations:
double buffer FIR filter #1 h[1, m, i]=Hm[i],
double buffer FIR filter #2 h[2, m, i]=Hm[i]+Hm[i+2M],
double buffer FIR filter #3 h[3, m, i]=Hm[i+2M],
double buffer FIR filter #4 h[4, m, i]=Hm[i+2M]+Hm[i+3M],
double buffer FIR filter #5 h[5, m, i]=Hm[i]+Hm[i+M]+Hm[i +2M]+Hm[i+3M],
double buffer FIR filter #6 h[6, m, i]=Hm[i]+Hm[i+M],
double buffer FIR filter #7 h[7, m, i]=Hm[i+M],
double buffer FIR filter #8 h[8, m, i]=Hm[i+M]+Hm[i+3M],
double buffer FIR filter #9 h[9, m, i]=Hm[i+3M].
The samples from the outputs of the assembly of double buffers FIR filters 5 are applied to the inputs of the PostFIRs transformer 6. A PostFIRs transformer 6 converts K input samples into a set of N samples.
A PostFIRs transformer for N=2 converts three input samples (a[1], a[2], a[3]) into two output samples (b[1], b[2]) in compliance with the equations:
b[1]=a[1]+Da[3];
b[2]=a[2]−a[1]−a[3];
where D means a delay for M clock cycles.
A PostFIRs transformer for N>2 is built by an iterative procedure with n as iteration parameter (n equals two to the power of v, and k equals three to the power of v). A PostFIRs transformer of order 2n is built from k PostFIRs transformers of order 2, and two PostFIRs transformers of order n in accordance with the instructions:
The instructions presented above were used to build a PostFIRs transformer of order 4. The resulting block diagram is shown in
b[1]=a[1]+Da[3]+D(a[8]−a[7]−a[9]),
b[2]=a[6]+Da[4]−a[7]−Da[9],
b[3]=a[2]−a[1]−a[3]+a[7]+Da[9],
b[4]=a[5]−a[4]−a[6]−a[8]+a[7]+a[9]−a[2]+a[1]+a[3]).
The outputs of the PostFIRs transformer 6 are connected to the inputs of the samples sequence restoration unit 7. This unit transforms the input samples
(y[r·NM+m], y[r·NM+m+M], y[r·NM+m+2M], . . . , y[r·NM+m+(N−1)M]),
into a set of output samples (y[4k], y[4k+1], . . . , y[4k+N−1]).
The outputs of the samples sequence restoration unit 7 are connected to the inputs of the multiplexer 8. The multiplexer 8 combines its input samples into an output sequence y[n], with the samples y[n] following the equation (1) presented above in the beginning of the detailed description.
If the interleaved ADC works with the samples rate F samples/s, and the equalizer length is L taps, then an equalizer with a structure of a common FIR filter has to carry out N1=F·L multiplication per second. In
When the samples stream produced by the interleaved ADC is split into two samples sub streams (N=2), the number of multipliers needed to construct a digital equalizer for an interleaved ADC according the present technology is reduced by a factor of N1/N2=4/3=1.33, which is noticeable smaller than the factor N1/N2=1.78 achievable for N=4 (splitting into four samples sub streams). On the other hand, when N>4 and the samples stream produced by the interleaved ADC is split into eight, sixteen or bigger number of samples sub streams, the factor of multipliers reduction increses, but the increased number of additional adders and delay lines makes the needed hardware prohibitively expensive. It seems that the case N=4 is most preferable from a practical use standpoint.
One skilled in the art will realize the invention may be embodied in other specific forms without departing from the spirit or essential characteristics thereof. The foregoing embodiments are therefore to be considered in all respects illustrative rather than limiting of the invention described herein. The scope of the invention is thus indicated by the appended claims, rather than by the foregoing description, and all changes that come within the meaning and range of equivalency of the claims are therefore intended to be embraced therein.
This application claims the benefit of U.S. Provisional Patent Application No. 61/601,360, filed Feb. 21, 2012, the contents of which are incorporated by reference herein in it's entirety.
Number | Date | Country | |
---|---|---|---|
61601360 | Feb 2012 | US |