CN115801009B - Method for compensating time offset error of TIADC parallel acquisition system - Google Patents
Method for compensating time offset error of TIADC parallel acquisition system Download PDFInfo
- Publication number
- CN115801009B CN115801009B CN202310044449.7A CN202310044449A CN115801009B CN 115801009 B CN115801009 B CN 115801009B CN 202310044449 A CN202310044449 A CN 202310044449A CN 115801009 B CN115801009 B CN 115801009B
- Authority
- CN
- China
- Prior art keywords
- output
- channel
- ref
- differentiator
- sub
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Abstract
The invention discloses a method for compensating time offset errors of a TIADC parallel acquisition system, and belongs to the field of high-speed high-precision analog-digital conversion. The method comprises the following steps: calculating to obtain the output of the sub-channel band time mismatch error; calibrating and compensating the output of the sub-channel with time mismatch error; the differentiator is used for deriving and the output of the differentiator is corrected. The invention mainly aims at realizing the differentiator by using a five-point value differential formula by the compensation circuit, and greatly reducing the area required by the differentiator compared with the direct FIR filter for deriving; after the differentiator is adjusted by the differentiator correction module, the performance of the compensated signal is obviously better than that of the compensated signal.
Description
Technical Field
The invention relates to the technical field of high-speed high-precision analog-digital conversion, in particular to a method for compensating time offset errors of a TIADC parallel acquisition system.
Background
A time interleaved ADC is an architecture that loops through a set of M sub-ADCs, so the overall throughput is M times the sampling rate of each sub-ADC. Theoretically, any type or configuration of ADC may be used as a sub-ADC for TIADC. The principle structure of TIADC is shown in figure 1, and compared with a single-channel ADC, the architecture can keep the same precision as that of a sub-ADC and further improve the sampling rate. However, the timing of the channels of TIADC is shown in fig. 2, and various mismatch errors exist between the interleaved channels, which severely degrade the performance of the TI-ADC, wherein the time mismatch errors are most difficult to eliminate.
Assume that the sampling period of the whole system isEach sub-ADC sampling period is +.>Digital output of analog input signal after TIADC sampling quantization>The method comprises the following steps:
wherein the method comprises the steps ofFor the output of sub-ADCs, +.>For TIADC output, +.>Is the sampling period of TIADC. For a pair ofPerforming discrete fourier transform can obtain the expression in the frequency domain as:
for the output frequency response of TIADC, +.>For the output frequency response of the ith sub-channel ADC, M is TIADC channel number, +.>Is the sampling angular frequency. If sampling time mismatch exists, sampling of TIADC is uneven, and the time mismatch error of the ith channel is assumed to be +.>Then the digital output of TIADC is:
for the error coefficient of the ith channel, usually +.>. For->Performing discrete fourier transform can obtain the expression in the frequency domain as:
when the input signal is a sinusoidal signal:
to output the angular frequency +.>Is an impulse function. As can be seen from the formula, the time mismatch error presents redundant harmonics on the frequency spectrum, resulting in a substantial reduction in performance of TIADC, and the frequency point where the time mismatch error occurs on the frequency spectrum is located:
Disclosure of Invention
The invention aims to provide a method for compensating time offset errors of a TIADC parallel acquisition system, which aims to solve the problems in the background technology.
In order to solve the technical problems, the invention provides a method for compensating time offset errors of a TIADC parallel acquisition system, which comprises the following steps:
calculating to obtain the output of the sub-channel band time mismatch error;
calibrating and compensating the output of the sub-channel with time mismatch error;
the differentiator is used for deriving and the output of the differentiator is corrected.
In one embodiment, calculating the output of the sub-channel band time mismatch error includes:
switching the 0 channel output to the reference channel output, and estimatingTime mismatch error of sub-channels;
Switching the 0 channel output from the reference channel back to the 0 channel output after convergence, estimatesTime mismatch error of sub-channel->;
Due to,/>The method comprises the steps of carrying out a first treatment on the surface of the The output of the sub-channel band time mismatch error expressed by the Taylor series expansion formula is as follows:
for the ideal output of the ith sub-channel, k is constant,/->To the k-th power of the time mismatch coefficient of the i-th channel,is the k-th derivative of the output of the ith channel.
In one embodiment, the output of the sub-channel with time mismatch error after calibration compensation is:
to 3 th-order of the time mismatch coefficient of the ith channel,/->Is the 3 rd order derivative of the output of the i-th channel.
In one embodiment, the differentiator calculates the first derivative value using a first order numerical differentiation formula based on a lagrangian interpolation polynomial of a five-point formula:
for the first order of TIADC output signal, < >>、/>、/>、/>Four sampling points adjacent to each other.
In one embodiment, modifying the output of the differentiator comprises:
let the original first derivative value be set to:
in obtaining the time mismatch error between the reference channel and the 0 channelAnd->Later, the relative error between the two channels is obtained>;
The same inputs are infused for the 0 channel and the reference channel, and then the outputs of the two channels are respectively:
for the output of the reference channel->Output of the 0 th channel; expansion to the first order with taylor series yields:
a derivative value for an ideal sampling point; since the differential value of the ideal point cannot be known, the actual output differential value is used>Approximate replacement->Finally expressed as:
because it is assumed that there is a coefficient to be corrected between the output of the differentiator and the actual derivative valueBy extracting correction factors->Expressed by the formula:
after obtaining the coefficient to be corrected, the final compensation is performed on the output of the differentiatorThe correction of the differentiator can be achieved.
The method for compensating the time offset error of the TIADC parallel acquisition system has the following beneficial effects:
(1) The invention mainly aims at realizing the differentiator by using a five-point value differential formula by the compensation circuit, and greatly reducing the area required by the differentiator compared with the direct FIR filter for deriving: only four registers, two multipliers and one adder are needed for one differentiator;
(2) The key path is shorter, and the synthesis is easier to realize higher frequency;
(3) The simulation shows that the structure designed by the traditional five-point value differential formula has the same error as the corrected differentiator structure, and the same frequency point can be improved by more than 30 dB.
Drawings
Fig. 1 is a schematic illustration of TIADC.
Fig. 2 is a timing diagram of the channels of TIADC.
Fig. 3 is a schematic diagram of a TIADC time mismatch error calibration circuit.
Fig. 4 is a schematic diagram of a calibration module configuration.
Fig. 5 is a schematic diagram of the differentiator modification block.
Fig. 6 is a schematic diagram of a differential structure after correction of coefficients.
FIG. 7 is a schematic diagram of a preferred embodiment of the present inventionSFDR before frequency bin calibration.
Detailed Description
The method for compensating the time offset error of the TIADC parallel acquisition system provided by the invention is further described in detail below with reference to the accompanying drawings and specific embodiments. The advantages and features of the present invention will become more apparent from the following description. It should be noted that the drawings are in a very simplified form and are all to a non-precise scale, merely for convenience and clarity in aiding in the description of embodiments of the invention.
The invention provides a method for compensating time offset errors of a TIADC parallel acquisition system, which improves the quality of sampling output signals of the system by compensating time mismatch errors and output errors of a differentiator under the condition that the offset and gain errors of the whole TIADC parallel acquisition system are assumed to be zero.
As shown in fig. 3, the compensation structure of the present invention needs to combine with the TIADC feedforward type time mismatch error estimation structure, in this embodiment, taking four channels as an example, firstly, 0 channel output is switched to reference channel output, and time mismatch errors of four sub-channels are estimated. Switching the 0 channel output from the reference channel back to the 0 channel output after convergence, estimating the time mismatch error of the four sub-channels>。
Typically the time mismatch error of the sub-channels is much smaller than 1, i.e,/>. Therefore, the Taylor series expansion formula is used for representing the output +.>:
K is a constant (k=1, 2, …) for the ideal output of the ith sub-channel,>k-th power of the time mismatch coefficient of the ith channel,/-th power of>Is the k-th derivative of the output of the ith channel. The compensation structure compensates according to the calibration module structure of fig. 4 after obtaining the output of the sub-channel band time mismatch error according to the conventional high-order cascade structure, and the final output result is:
to 3 th-order of the time mismatch coefficient of the ith channel,/->For the output of the ith channel, +.>Is the 3 rd order derivative of the output of the i-th channel. The calibrated output removes the first and second order error terms because the performance of the final output is related to the estimation of the time mismatch error on the one hand and the accuracy of the derivative with the differentiator on the other hand.
The differentiator usually adopts a direct type FIR filter structure, and the order is designed to be more than 30 orders to meet the performance requirement, and the differentiator adopts a first order numerical differential formula of a Lagrange interpolation polynomial based on a five-point formula to calculate the first derivative value of a TIADC output signal:
for the first order of TIADC output signal, < >>、/>、/>、/>Four sampling points adjacent to each other. Since there is a large deviation between the error derived by the five-point value differential formula and the actual derivative value, the output of the differentiator is corrected again, and the structure is as shown in fig. 6, and the original first derivative value is set as:
is a correction coefficient. A coefficient exists between the derivative value obtained by the five-point value differentiation method and the actual derivative value ∈>Representing the deviation from the actual derivative value, an additional differentiator modification block is introduced as shown in fig. 5.
In obtaining the time mismatch error between the reference channel and the 0 channelAnd->Later, the relative error between the two channels can be obtained +.>:
The same inputs are infused for the reference channel and the 0 channel, and then the outputs of the two channels are respectively:
for the output of the reference channel->Is the output of channel 0. Expansion to the first order with the taylor series can result in:
is the derivative value of the ideal sampling point. Since the differential value of the ideal point cannot be known, the actual output differential value is used>Approximate replacement->Finally expressed as:
because it is assumed that there is a coefficient to be corrected between the output of the differentiator and the actual derivative valueCorrection coefficient is extracted by the structure of fig. 5 +.>Expressed by the formula: />
after the coefficients to be corrected are obtained, as shown in FIG. 6, the output of the differentiator is finally compensatedThe correction of the differentiator can be realized, and the signal performance after final compensation can be better improved.
The invention compensates the time mismatch errors existing among multiple channels in TIADC by using a Taylor series expansion formula, and the performance is mainly dependent on the time mismatch errorsDifference of differenceAnd the accuracy of the differentiator output, the invention mainly improves the performance of the signal after compensation by improving the accuracy of the differentiator output.
After the differentiator is adjusted by the differentiator correction module with the structure of fig. 5, the performance of the compensated signal is obviously better than that of the compensated signal. Compared with the direct FIR type filter for deriving, the current structure can have better performance under the same time mismatch error and the same compensation structure by comparing the performance of the 33-order direct FIR type filter. The invention is illustrated with a 4-channel structure that is suitable for TIADC in a multi-channel structure. Based on the above calibration scheme, the four-channel TIADC is set with time mismatch errors of [ -1 mill, -2 mill, -3 mill, -4 millrespectively]The reference channel is set without time mismatch error, inAnd (5) calibrating the frequency points. From fig. 7 and 8, it can be seen that the SFDR of the output signal after calibration is increased from 46.62dB to 102.11dB, and the calibration effect is very obvious.
The above description is only illustrative of the preferred embodiments of the present invention and is not intended to limit the scope of the present invention, and any alterations and modifications made by those skilled in the art based on the above disclosure shall fall within the scope of the appended claims.
Claims (1)
1. A method for compensating for time offset errors in a TIADC parallel acquisition system, comprising:
calculating to obtain the output of the sub-channel band time mismatch error;
calibrating and compensating the output of the sub-channel with time mismatch error;
the differentiator is adopted for deriving, and the output of the differentiator is corrected;
the output of the sub-channel band time mismatch error obtained by calculation comprises the following steps:
switching the 0 channel output to the reference channel output to estimate x sub-channelsTime mismatch error r of tracks ref ,r 1 ,...,r x ;
Switching the 0-channel output from the reference channel back to the 0-channel output after convergence, estimating the time mismatch error r of the x sub-channels 0 ,r 1 ,...,r x ;
Due to r i <<1, i=0, 1,2,; the output of the sub-channel band time mismatch error expressed by the Taylor series expansion formula is as follows:
x i (n) is the ideal output of the ith sub-channel, k is a constant,k-th power of the time mismatch coefficient of the ith channel,/-th power of>A k-order derivative of the output of the ith channel;
the output after calibration compensation for the output of the sub-channel band time mismatch error is:
to the 3 rd power of the time mismatch coefficient of the ith channel, x '' i (n) is a 3-order derivative of the output of the ith channel;
the differentiator calculates a first derivative value using a first order numerical differential equation of a lagrangian interpolation polynomial based on a five-point equation:
y' (n) is a first-order derivative of the TIADC output signal, and y (n+2), y (n+1), y (n-1) and y (n-2) are four sampling points which are adjacent to each other front and back;
correcting the output of the differentiator includes:
the original first derivative value is set as:
correction coefficient r derivator Representing the deviation from the actual derivative value;
in obtaining the time mismatch error r of the reference channel and the 0 channel ref And r 0 Later, a relative error r between the two channels is obtained 0-ref =r 0 -r ref ;
The same inputs are infused for the 0 channel and the reference channel, and then the outputs of the two channels are respectively:
y ref (n)=x(n)
y 0 (n)=x(n+r 0-ref )
y ref (n) is the output of the reference channel, y 0 (n) is the output of the 0 th channel; expansion to the first order with taylor series yields:
y 0 (n)=x(n+r 0-ref )≈x(n)+r 0-ref x(n)
x' (n) is the derivative value of the ideal sampling point; since the differential value of the ideal point cannot be known, the actual differential value y 'is used' 0 (n) approximately replaces x' (n), and is finally expressed as:
y 0 (n)≈x(n)+r 0-ref y 0 ′(n)
because there is a coefficient r to be corrected between the output of the differentiator and the actual derivative value derivator By extracting correction coefficient g derivator Expressed by the formula:
g derivator is the reciprocal of the correction coefficient; because:
y 0 (n)-r 0-ref r derivator y 0 (n)=y 0 (n)-r 0-ref y 0 (n)+r 0-ref (r derivator -1)y 0 (n)
=x(n)+r 0-ref (r derivator -1)y 0 (n)
binding y ref (n) =x (n) gives:
y ref (n)-[y 0 (n)-r 0-ref r derivator y 0 (n)]=r 0-ref (r derivator -1)y 0 (n)
output g derivator Finally expressed as:
after obtaining the coefficient to be corrected, g is compensated at the final output of the differentiator derivator The correction of the differentiator can be achieved.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310044449.7A CN115801009B (en) | 2023-01-30 | 2023-01-30 | Method for compensating time offset error of TIADC parallel acquisition system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310044449.7A CN115801009B (en) | 2023-01-30 | 2023-01-30 | Method for compensating time offset error of TIADC parallel acquisition system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115801009A CN115801009A (en) | 2023-03-14 |
CN115801009B true CN115801009B (en) | 2023-05-09 |
Family
ID=85429104
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310044449.7A Active CN115801009B (en) | 2023-01-30 | 2023-01-30 | Method for compensating time offset error of TIADC parallel acquisition system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115801009B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116594593B (en) * | 2023-07-13 | 2023-09-26 | 上海芯炽科技集团有限公司 | Digital differentiator based on Lechassen extrapolation and optimized structure thereof |
CN116781079A (en) * | 2023-08-22 | 2023-09-19 | 上海芯炽科技集团有限公司 | TIADC time mismatch error calibration circuit based on reference channel |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8009070B2 (en) * | 2007-06-21 | 2011-08-30 | Signal Processing Devices Sweden Ab | Compensation of mismatch errors in a time-interleaved analog-to-digital converter |
JP5286420B2 (en) * | 2009-09-30 | 2013-09-11 | 株式会社日立製作所 | Analog-digital converter and semiconductor integrated circuit device using the same |
CN104038226B (en) * | 2014-06-25 | 2018-06-05 | 华为技术有限公司 | The time-interleaved analog-digital converter of multichannel |
US9531398B2 (en) * | 2015-03-11 | 2016-12-27 | Analog Devices, Inc. | Limiting aging effects in analog differential circuits |
CN108111169B (en) * | 2018-01-03 | 2021-04-20 | 中山大学 | Combined correction method for linear mismatch and nonlinear mismatch of four-channel TIADC |
CN108494402B (en) * | 2018-03-14 | 2021-07-27 | 东南大学 | TIADC system error estimation and compensation method based on sine fitting |
CN113346902B (en) * | 2021-06-16 | 2023-04-07 | 合肥工业大学 | All-digital calibration structure based on TIADC composite output and calibration method thereof |
-
2023
- 2023-01-30 CN CN202310044449.7A patent/CN115801009B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN115801009A (en) | 2023-03-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115801009B (en) | Method for compensating time offset error of TIADC parallel acquisition system | |
CN108471313B (en) | Digital-analog mixed signal-based TIADC system calibration method | |
US20170117914A1 (en) | Method and apparatus for providing digital background calibration for mismatches in m-channel time-interleved adcs (ti-adcs) | |
US8519875B2 (en) | System and method for background calibration of time interleaved analog to digital converters | |
US5061925A (en) | Phase equalization system for a digital-to-analog converter utilizing separate digital and analog sections | |
US8836550B2 (en) | System and method for digitally correcting mismatches in multipath ADCs | |
CN102075464A (en) | Joint estimation and real-time correction method for channel error of TIADC system | |
CN109361390B (en) | Inter-channel sampling time error correction module and method for time-interleaved ADC | |
CN102291141A (en) | Time-interleaved split ADC (Analog-to-Digital Converter) calibration structure without redundant channel and adaptive calibration method thereof | |
CN112751564A (en) | Sampling clock phase mismatch error estimation method and device | |
CN111064469A (en) | Method for correcting TIADC sampling time mismatch error based on adjacent channel autocorrelation function | |
CN104393872A (en) | Sampling time error correction method for multi-channel parallel analog-digital converter (ADC) system | |
US7808407B2 (en) | Sub-channel distortion mitigation in parallel digital systems | |
Centurelli et al. | Improved digital background calibration of time-interleaved pipeline A/D converters | |
US8307248B2 (en) | Method and a system for estimating errors introduced in a time-interleaved analog-to-digital converter system | |
Li et al. | A background correlation-based timing skew estimation method for time-interleaved ADCs | |
US9748967B1 (en) | Periodic signal averaging with a time interleaving analog to digital converter | |
CN115021754A (en) | TIADC sampling time mismatch digital correction method and system | |
CN116594593B (en) | Digital differentiator based on Lechassen extrapolation and optimized structure thereof | |
Khan et al. | A fully digital background calibration technique for M-channel time-interleaved ADCs | |
CN113517890A (en) | Extraction method for sampling time mismatch of time-interleaved ADC (analog to digital converter) | |
CN116781079A (en) | TIADC time mismatch error calibration circuit based on reference channel | |
CN110401445B (en) | TIADC system sampling time mismatch digital background calibration method with low hardware overhead | |
CN113346902B (en) | All-digital calibration structure based on TIADC composite output and calibration method thereof | |
Ta et al. | Fully digital background calibration technique for channel mismatches in TIADCs |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |