This is a national stage filing under section 371 of PCT/EP2014/075267 filed Nov. 21, 2014, which claims the benefit of French Application No. 1362091, filed on Dec. 4, 2013, which application is hereby incorporated herein by reference.
Embodiments of the invention are related to a method and device for compensating bandwidth mismatches of time interleaved analog-to-digital converters.
Certain present-day applications require analog-to-digital converters having a high sampling frequency and a high resolution. Among these applications are cable TV transmissions, optical communications and satellite communications in which one or more modulated signals are transmitted simultaneously on a wide frequency band (for example several GHz wide).
Moreover, the use of increasingly complex modulations requires a digitizing of the signals using a large number of bits.
Single-converter architectures have difficulty in meeting the required performance levels at an acceptable cost in terms of consumption and of space requirement.
Moreover, structures having several time-interleaved analog-to digital converters (TIADCs: Time-Interleaved Analog to Digital Converters) are being imposed in order to meet this requirement. In such a time-interleaved structure, the time-interleaved analog-to-digital converters respectively carry out time-shifted analog-to-digital conversions of an analog signal. In other words, if the structure comprises M converters, the latter successively sample, each in their turn, the signal at a frequency equal to Fs/M, where Fs is the overall sampling frequency of the structure.
However, the disadvantage of this type of structure is that mismatches, even slight, between the converters create parasitic lines or frequency bands that can be situated in the frequency zone of the useful signal. Moreover, in the case, for example, of television signals containing several channels, these parasitic bands or lines can interfere with some of these channels.
These mismatches can have various causes, such as for example timing skews and possibly static gains and/or offsets that are different between the converters of the structure.
Another cause of mismatch is the inequality between the values of the bandwidths of the different converters of the structure.
It is in fact known to those skilled in the art, for example from the article by Tsung-Heng Tsai and others entitiled “Bandwidth Mismatch and its Correction in Time-Interleaved Analog-to-Digital Converters,” IEEE Transactions on Circuits and Systems—II: Express Briefs, vol. 53, no. 10, October 2006, that an analog-to-digital converter can be represented in the first order by a first order low-pass filter and by an ideal sampler. The bandwidth of the converter is then defined by the cut-off frequency or the time constant of this first order low-pass filter. It is then said that there is bandwidth mismatch of the different converters of the structure when the different cut-off frequencies or different time constants of the different first order low-pass filters are different.
In practice, the cut-off frequency of these filters is high, typically of the order of one GHz or of about ten GHz, in order to prevent a cutting off of the useful signal. Such a difference in cut-off frequencies causes frequency bands or parasitic tones in the sampled signal. Moreover, as the input signal becomes higher and close to the cut-off frequency of the filters, the power of the parasitic frequency bands becomes greater with respect to the useful signal. This is particularly critical for analog input signals having high frequencies, typically of the order of a few tens of GHz, such as those found for example in optical communications.
The abovementioned article proposes a solution for correcting the bandwidth mismatch based on the hypothesis that the different time constants of the low-pass filters of the converters are already known, which is a limited solution.
The article by Shahzad Saleem and others, entitled “Adaptative Blind Background Calibration of Polynomial-Represented Frequency Response Mismatches in a Two-Channel Time-Interleaved ADC,” IEEE Transactions on Circuits and Systems—I: Regular Papers, vol. 58, no. 6, June 2011, deals with bandwidth mismatches of any order greater than or equal to 1.
In this second article, the solution described for dealing with these mismatches is valid only for an input signal having strict constraints in terms of frequency spectrum and bandwidth. The frequency spectrum must thus notably exhibit “holes” and the proposed solution aims to minimize the power of the parasitic frequency noise, appearing in these holes because of the mismatches, by way of iterative and retroactive corrections (that is to say carried out upstream of the converters). Here again, such a solution is limited and is moreover complex.
Embodiments of the invention notably but not exclusively apply to the correction of the mismatch between the different bandwidths of several so called “time-interleaved” analog-to-digital converters, that is to say respectively carrying out time-shifted analog-to-digital conversions of an analog signal in order to increase the overall sampling frequency.
Embodiments of the invention relate to the processing of a sampled signal and, more particularly, the correction of that sampled signal notably in order to correct the bandwidth mismatch of several analog-to-digital converters of an analog-to-digital conversion system delivering that sampled signal.
According to one implementation and embodiment, it is proposed to compensate for the bandwidth mismatch of a structure of time-interleaved converters in a “forward” manner, that is to say without retroactive correction of the samples to be corrected based on information situated downstream of the point of reception of those samples, and to do so in a single iteration, that is to say without a feedback loop, and therefore without risk of divergence, and in a “blind” manner, that is to say by acting on the signal actually delivered by the converters, whatever that signal may be, without it being necessary to provide a calibration phase using a known signal at the input of the converters.
According to one aspect, there is proposed a method of direct compensation of the bandwidth mismatch of M time-interleaved converters respectively delivering M original trains of M original samples, M being greater than two, each converter being considered in the first order as comprising a first order low-pass filter. The method comprises a process of digital filtering of the M original trains delivering M filtered trains of filtered samples. This filtering has a transfer function substantially equal to the product of the transfer function of a reference low-pass filter and the transfer function of a derivative filter multiplied by the reference time constant of the reference low-pass filter and, for each original train of original samples.
An estimation process delivers an estimated difference between the time constant of the first order low-pass filter associated with the corresponding converter and the reference time constant, preferably normalized by the reference time constant, and comprising a first generation of a first item of differentiated power information relative to the original train and a second generation of a second item of information of differentiated power relative to at least one of the M filtered trains.
A process of correction of the original samples of the original train uses the filtered samples of the corresponding filtered train and the corresponding estimated difference, in such a way as to deliver a corrected train of corrected samples. The M corrected trains are then considered as coming from the corresponding M converters, the first order low-pass filters of which all have the reference time constant as their time constant.
The method according to this first aspect thus makes it possible to carry out the compensation of bandwidths of the converters, and notably the estimation of the differences of time constants, in a direct (forward) manner, that is to say without retroactive correction of the original samples from an item of information located downstream of the point of reception of the original samples to be corrected. This therefore requires only a single iteration, that it to say no feedback loop and without risk of divergence or any necessity of providing a calibration phase. Because of this, the hardware embodiment is thereby simplified and the estimation of the mismatch is faster.
The term “compensation” must of course be understood to be a compensation that is not necessarily perfect because, in practice, the correction can result in low-pass filters all having the reference time constant to within an estimation error.
Moreover, the concept of “differentiated” power information signifies that in fact this differentiated power information comprises two different power indications associated with the same train, which will make it possible to eliminate the static gain of the corresponding converter in the determination of the estimated difference between the time constant of the first order low-pass filter associated with this converter and the reference time constant.
According to one implementation, the process of filtering comprises a filtering of the original trains of original samples by a digital filter, for example a finite impulse response filter, the transfer function of which is substantially equal to jωτref/(1+jωτref), in which τref denotes the time constant of the first order reference low-pass filter.
According to one implementation, the first generation comprises a first filtering of the corresponding original train with a first filter and a determination of a first indication of power of the original train filtered by this first filter, and a second filtering of the corresponding original train within a second filter, different from the first filter, and a determination of a second indication of power of the original train filtered by this second filter, the first differentiated power information comprising the first and second power indications, and the second generation comprises a third filtering of the at least one of the M trains filtered with the first filter and a determination of a third indication of power of the at least one of the M filtered trains filtered by this first filter, and a fourth filtering of the at least one of the M trains filtered with the second filter and a determination of a fourth indication of power of the at least one of the M filtered trains filtered by this second filter, the second differentiated power information comprising the third and fourth power indications.
This second differentiated power information could, as a variant, be an item of averaged power information obtained for example by averaging the third and fourth power indications associated with some or all of the filtered trains.
According to one implementation, the process of estimation furthermore comprises a determination of a first averaged power indication obtained by an average of the M first power indications and a determination of a second averaged power indication obtained by an average of the M second power indications, and the estimated difference is obtained for the corresponding original train from the first corresponding power indication, from the second corresponding power indication, from the third power indication, from the fourth power indication, from the first averaged power indication and from the second averaged power indication.
According to one implementation, the correction process comprises, for each original train, an addition to each original sample of the product of the homologous filtered sample of the filtered train corresponding to this original train and the corresponding estimated difference.
According to another aspect, a device can be used for the direct compensation of the bandwidth mismatch of M time-interleaved converters, M being greater than two, each converter is considered in the first order as comprising a first order low-pass filter. An input receives M original trains of original samples respectively coming from the M converters. Output delivers M corrected trains of corrected samples. A processor comprises a digital filter configured to carry out a filtering process on the M original trains in such a way as to deliver M corresponding filtered trains of filtered samples. This filter has a transfer function substantially equal to the product of the transfer function of a reference low-pass filter and the transfer function of a derivative filter multiplied by the reference time constant of the reference low-pass filter.
An estimator is configured to carry out an estimation process delivering, for each original train, an estimated difference between the time constant of the first order low-pass filter associated with the corresponding converter and the reference time constant, preferably normalized by the reference constant, and comprising a first generation of a first item of differential differentiated power information relative to each original train and a second generation of a second item of differentiated power information relative to at least one of the M filtered trains. A correction circuit is configured to carry out a process of correction of the original samples of each original train using the filtered samples of the corresponding filtered trains and the corresponding estimated difference, in such a way as to deliver a corresponding corrected train of corrected samples, the M corrected trains then being considered as having come from the corresponding M converters, the first order low-pass filters of which all have the reference time constant as their time constant.
According to one embodiment, the filter comprises a filter whose transfer function is substantially equal to jωτref/(1+jωτref), where τref denotes the time constant of the reference first order low-pass filter.
This filter can for example be a finite impulse response filter.
According to one embodiment, the estimator comprises a first stage comprising a first filter intended to filter the corresponding original train and a first module configured to determine a first power indication of the original train filtered by the first filter, and a second filter, different from the first filter, intended to filter the corresponding original train and a second module configured to determine a second power indication of the original train filtered by the second filter. The first differentiated power information comprises the first and second power indications. A second stage comprises the first filter intended to filter the at least one of the M filtered trains. A third module is configured to determine a third power indication of the at least one of the M filtered trains filtered by the first filter. The second filter is intended to filter the at least one of the M filtered trains. A fourth module is configured to determine a fourth power indication of the at least one of the M filtered trains filtered by this second filter. The second differentiated power information comprises the third and fourth power indications. According to one embodiment, the estimator comprises a third stage configured to determine a first averaged power indication obtained by an average of the M first power indications and a second averaged power indication obtained by an average of the M second power indications. The estimator also comprises a calculator configured to calculate the estimated difference for the corresponding original train, from the corresponding first power indication, from the second corresponding power indication, from the third power indication, from the fourth power indication, from the first averaged power indication and from the second averaged power indication.
According to one embodiment, the correction circuit is configured, for each original train, to add to each original sample the product of the homologous filtered sample of the filtered train corresponding to this original train and the corresponding estimated difference.
According to another aspect, there is proposed an integrated circuit comprising a compensation device such as defined above.
Other advantages and features of the invention will become apparent on examination of the detailed description of implementations and embodiments, that are in no way limiting, and of the appended drawings in which:
In
The analog-to-digital converter ADCm is controlled by the clock signal clkm.
As shown more particularly in
In this embodiment, each converter ADCm delivers an original train of original samples {tilde over (x)}m[k].
As shown diagrammatically in
The transfer function Bm(ω) of the low-pass filter FPBm is defined by the following formula (1):
in which ωm denotes the cut-off pulsatance of the filter.
The time constant τm of the filter is equal to 1/ωm.
The M analog-to-digital converters of the conversion system SCV do not generally all have the same time constant and consequently the same bandwidth. This mismatch of bandwidths, or of time constants, adds parasitic frequency bands to the frequency spectrum of the output signal.
There is therefore proposed a method for processing the sampled signal, the purpose of which is to compensate for this bandwidth or time constant mismatch.
This processing method is used within a processing device DIS the M input terminals BE0-BEM-1 of which receive the M trains of samples {tilde over (x)}m[k], where m varies from 0 to M−1.
In this respect, the device DIS comprises processor MT, which delivers at its output corrected trains of corrected samples {circumflex over (x)}m[k], where m varies from 0 to M−1.
Moreover, these M corrected trains are then considered as having come from the corresponding M converters the first order low-pass filters of which all have the same time constant which, as will be seen in more detail below, is a reference time constant of a reference first order low-pass filter.
Reference will now be made more particularly to
As shown in
These filter MFL has a transfer function substantially equal to the product of the transfer function of a reference low-pass filter and the transfer function of a derivative filter multiplied by the reference time constant τref of the reference low-pass filter.
As is well known by those skilled in the art, a first order low-pass filter can be represented by a capacitive-resistive network. Knowing the hardware structure of each of the M an analog-to-digital converters ADCm, it is therefore possible, notably by simulation, to determine the statistical distribution of the time constants τm of the first order low-pass filters associated with these converters.
Even though it is possible to take an arbitrary value for the reference time constant τref, it is preferable, in order to minimize the degradation of the corrected sampled signal, to choose a value substantially equal to the average of the distribution of the time constants obtained by simulation as the reference time constant τref.
The processor MT also comprises an estimator MEST, which is configured to carry out an estimation process delivering an estimated difference {circumflex over (Δ)}τm between the time constant τm of the first order low-pass filter associated with the corresponding converter ADCm and the reference time constant τref.
This estimated difference is in this case the estimation of the relative or normalized difference:
As will be seen in more detail below, this estimation processing comprises a first generation of a first item of differentiated power information Pm,1, Pm2 relative to each original train {tilde over (x)}m[k].
The estimation processing comprises moreover a second generation of a second item of differentiated power information Pf,1, Pf,2 relative to at least one of the M filtered trains ym(k).
In the example described, the filtered train used is the train yM-1[k]. This being so, it would have been entirely possible to use any one of these filtered trains. It would also have been possible to use some or all of the filtered trains for determining this second differentiated power information, for example by averaging the power indications Pf,1, Pf,2 associated with some or all of the filtered trains. However, it was chosen to use only one single filtered train for reasons of simplification and of space requirement.
As will be seen in more detail below, the concept of “differentiated” power information signifies that this differentiated power information in fact comprises two different power indications associated with the same train, which will make it possible to eliminate the static gain of the corresponding converter in the calculation of {circumflex over (Δ)}τm.
The processor MT comprises moreover a correction circuit MCOR comprising, as will be seen below, multipliers and adders configured to carry out a process of correction of the original samples {tilde over (x)}m[k] so as to deliver a corrected train of corresponding corrected samples {circumflex over (x)}m[k]
This correction processing uses the filtered samples of the corresponding filtered train ym[k] and the corresponding estimated difference {circumflex over (Δ)}τm.
Moreover, as mentioned above, the M corrected trains are then considered as having come from the corresponding M converters, the first order low-pass filters of which FPBm all have the reference time constant τref as their time constant.
A compensation of the bandwidth mismatch of the different analog-to-digital converters ADCm has therefore been carried out here. The term “compensation” should of course be understood to be a compensation that is not necessarily perfect because, in practice, the correction can result in low-pass filters all having the reference time constant to within an estimation error.
As long as the differences {circumflex over (Δ)}τm have not been estimated, the samples are for example corrected with zero values.
Similarly, the estimation can be carried out just once, the set of differences {circumflex over (Δ)}τm then remaining valid throughout the period of operation of the device.
However, the estimation can be repeated several times in order to update the set of differences {circumflex over (Δ)}τm in order to take account notably of the variation of operating temperature for example.
Reference will now be made more particularly to
The filter MFL comprises a filter, the transfer function of which is close to that defined by the formula (2):
This filter is produced, in the example described in
In the example described here, the original trains of samples arrive in parallel and are “serialized” by a multiplexer MUX (M samples of the sampled signal respectively coming from the M successive converters).
The filtered samples are then “re-parallelized” in a demultiplexer DMUX so as to deliver the M filtered trains of filtered samples.
This embodiment using a serializing of the original samples of the different original trains and then a re-parallelizing of the filtered trains is only one possible example of embodiment. In fact, it would be entirely possible to envisage carrying out the filtering processing in a parallel manner.
The frequency response C(ω) of this filter FIR, corresponding to the transfer function defined below by the formula (4), can then be expressed by the following formula (3):
By way of non-limiting example, for a reference time constant τref value equal to 0.45/Fs, it is possible to use a filter FIR having 17 coefficients (K=17), the values of which are given in an appendix.
If
If
The first stage ET1 also comprises a second filter FH2, different from the first filter FH1, intended to filter each original train {tilde over (x)}m[k] and a second module MD2 configured to determine a second power indication Pm,2 of the original train {tilde over (x)}m[k] filtered by the second filter FH2. These two power indications Pm,1 and Pm,2 form, as mentioned above, the first item of differentiated power information.
Several embodiments are of course possible. Thus, either provision is made for processing the original samples in parallel and consequently M assemblies similar to the one shown in
Even though it is possible to use any two filters FH1, FH2, insofar as they are different, that is to say that they have different transfer functions, it is preferable that the transfer functions of these filters should have as few zeros as possible.
Thus, by way of indication, it will be possible to choose the following transfer function H1 for the filter FH1:
H1(z)=z+z−1 (4)
and the following transfer function H2(z) for the filter FH2:
H2(z)=z−z−1 (5)
The first module MD1 comprises moreover a block 61 configured to square the value of each sample, that is to say, to determine its power.
The module MD1 comprises moreover a block 63 intended to calculate, over L samples, the average power of the train so as to deliver the power indication Pm,1.
This power indication Pm,1 is therefore an average power indication in this case.
The average calculation carried out in the block 63 can be of any type, for example a conventional average or a moving average. Moreover, the number L can for example be of the order of 1 million.
The structure of the block 62 and of the block 64 of the second module MD2 connected downstream of the filter FH2 is similar to that of the blocks 61 and 63.
The power indication Pm,2 delivered by the second module MD2 is therefore also an average power indication in this case.
The second stage ET2 (
More precisely, the filtered train yM-1[k] is filtered by the two filters FH1 and FH2 and the modules MD3 and MD4 comprising the blocks 71, 73 and 72, 74, deliver the second item of differentiated power information comprising a third power indication Pf,1 and a fourth power indication Pf,2.
These different power indications Pm,1, Pm,2, Pf,1, Pf,2 will be used in the block BLC of the estimator MEST for determining, for each original train {tilde over (x)}m[k], the estimated time constant difference {circumflex over (Δ)}τm.
More precisely, as shown in
Similarly, the third stage ET3 generates a second averaged power indication
More precisely,
The calculating module MDC then determines (
The fact that the power indications Pm,1 and Pm,2 as well as the power indications Pf,2 and Pf,1 were obtained by two filters having different transfer functions makes it possible for the estimation {circumflex over (Δ)}τm not to be interfered with by the static gain of the converter ADCm.
Similarly, the use of the power indications Pf,2 and Pf,1 makes it possible to calculate this estimated difference {circumflex over (Δ)}τm in a single iteration.
If
In each of the
The input signal is a multi-tone signal having several useful tones TNi.
It can be seen that, before compensation, the sampled signal delivered by the converters has parasitic frequency bands in addition to these tones TNi (
After compensation by the device DIS, these parasitic bands have been eliminated (
It has been seen above that the converters could have mismatches of various types such as for example timing skew mismatches, and also different static gains or different offsets.
The French patent application No. 1 351 447 describes an example of compensation of the offset mismatch.
The French patent application No. 1 351 445 describes an example of compensation of the gain mismatch.
The French patent application No. 1 351 448 describes an example of the timing skew mismatch.
The processor MT of the device DIS that has just been described can easily be incorporated in such a correction system.
Thus, by way of example, it would be possible firstly to compensate for the offsets, then to compensate for the bandwidth mismatch, then to compensate for the static gain mismatch and then, finally, to compensate for the timing skew mismatch.
Number | Date | Country | Kind |
---|---|---|---|
13 62091 | Dec 2013 | FR | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/EP2014/075267 | 11/21/2014 | WO | 00 |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2015/082233 | 6/11/2015 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
7782235 | Velazquez | Aug 2010 | B1 |
7916050 | Mujica | Mar 2011 | B1 |
7916051 | Sestok | Mar 2011 | B1 |
8094050 | Sestok | Jan 2012 | B2 |
8749410 | Van Der Goes | Jun 2014 | B1 |
8890739 | Lewis | Nov 2014 | B2 |
20070262895 | Stein | Nov 2007 | A1 |
20100192027 | Johansson | Jul 2010 | A1 |
20110234435 | Woodward | Sep 2011 | A1 |
20130015990 | Payne | Jan 2013 | A1 |
20130141261 | Johancsik | Jun 2013 | A1 |
20140009318 | Taluja | Jan 2014 | A1 |
20150263749 | Ramakrishnan | Sep 2015 | A1 |
20150303934 | Chiu | Oct 2015 | A1 |
20160079994 | Lee | Mar 2016 | A1 |
20160149582 | Ragab | May 2016 | A1 |
Entry |
---|
Ghanem, F. et al., “A Background Extraction Technique for Bandwidth Mismatch Error in a Two-Channel Time-Interleaved ADC,” IEEE 10th International New Circuits and Systems Conference (NEWCAS), Jun. 17, 2012, pp. 353-356. |
Saleem, S. et al., “Adaptive Blind Background Calibration of Polynomial-Represented Frequency Response Mismatches in a Two-Channel Time-Interleaved ADC,” IEEE Transactions on Circuits and Systems I: Regular Papers, vol. 58, No. 6, Jun. 2011, pp. 1300-1310. |
Tsai, T-H. et al., “Bandwidth Mismatch and Its Correction in Time-Interleaved Analog-to-Digital Converters,” IEEE Transactions on Circuits and Systems II: Express Briefs, vol. 53, No. 10, Oct. 2006, pp. 1133-1137. |
Number | Date | Country | |
---|---|---|---|
20170026052 A1 | Jan 2017 | US |