CN107809247A - A kind of High Speed High Precision ADC dynamic input-output characteristic curve method for rapidly testing - Google Patents
A kind of High Speed High Precision ADC dynamic input-output characteristic curve method for rapidly testing Download PDFInfo
- Publication number
- CN107809247A CN107809247A CN201710902634.XA CN201710902634A CN107809247A CN 107809247 A CN107809247 A CN 107809247A CN 201710902634 A CN201710902634 A CN 201710902634A CN 107809247 A CN107809247 A CN 107809247A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- mfrac
- formula
- 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.)
- Pending
Links
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
- H03M1/00—Analogue/digital conversion; Digital/analogue conversion
- H03M1/10—Calibration or testing
- H03M1/1071—Measuring or testing
- H03M1/1095—Measuring or testing for ac performance, i.e. dynamic testing
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
The invention discloses a kind of method for rapidly testing of High Speed High Precision ADC dynamic input-output characteristic curve.For High Speed High Precision ADC, this method of testing is modeled using the first kind and Chebyshev polynomial of the second kind to its dynamic input-output characteristic curve simultaneously first, and thus obtains exporting the relation of code value and analog input value.The ADC output digital codes of lower collection are then inputted using high frequency sinusoidal signal based on the model and known pumping signal information establishes one group of matrix equation, and solved using least-square fitting approach, the last quick dynamic input-output characteristic curve for obtaining fitting and corresponding dynamic I NL estimation results.The method of testing avoids tradition and improves the problems such as histogram sampling number is more and related frequency domain estimates strict correlation sampling condition in ADC dynamic input-output characteristic curve method of testings, and these achievements are for realizing that quickly test and checking have highly important practical value to High Speed High Precision ADC.
Description
Technical field
The present invention relates to High Speed High Precision ADC testing field, and in particular to ADC batch testings and checking, for ADC once
Property quick estimating algorithm research etc..
Background technology
Analog-digital converter (Analog-to-Digital Convertor, ADC) is current digital-to-analogue mixed signal sum word letter
Extremely critical comprising modules in number two big systems of processing.It is right recently as the fast development of SOC technologies and communications industry
ADC performance requirement also more and more higher.Wherein production by assembly line due to can reach simultaneously at a high speed, high-precision performance requirement,
Have been widely used for every field, such as radio-frequency technique, multimedia-data procession and automatic test instrument etc..For needing
The ADC manufacturers of this kind of device quality are assessed, the testing time is one of mostly important focus.How fast and effectively to test
High Speed High Precision ADC had both been directly connected to production cycle and its service life of chip, and influenceed the market accreditation of chip indirectly
Degree.For High Speed High Precision ADC chip, its design and processes cost of manufacture originally takes greatly, if surveyed in chip
The testing time short enough can not also be ensured in examination flow, totle drilling cost will greatly increase.Therefore, in order to improve the reliable of product
Property, and consider for testing cost angle, current maximum obstacle is to lack a kind of effective High Speed High Precision ADC fast
Fast method of testing.
Usual ADC performance indications are divided into two classes, and one kind is comprising the static state including differential nonlinearity, integral nonlinearity etc.
Parameter, one kind are comprising the dynamic parameter including SFDR, signal to noise ratio, total harmonic distortion etc..Existing international survey
Test-object standard is to analyze above-mentioned two classes parameter by gathering two groups of different digital codes.Such as 16 High Speed High Precision ADCs, on the one hand,
Under conditions of each output digital code 30 samples of averaged acquisition are ensured, static parameter is tested using histogram method.It is another
Aspect, under the conditions of clock synchronization and correlation sampling is ensured, the test of dynamic parameter is completed using FFT.It is true
On, if from the few dynamic parameter of sample size needed for test static parameter can be estimated indirectly, then test and test entirely
The time needed for static parameter test need not be just consumed in card.
In this context, based on to ADC single data acquisition come test the estimation problem of ADC performance indications early stage by
Gradually attract attention, and up to the present have larger development.Wherein more representative ADC estimating algorithms, are such as based on
DFT frequency domain evaluation method, parametric spectral estimation algorithm and the frequency domain evaluation method based on Bessel function etc., it has been applied to real
The pure static tests of ADC on border.But last decade, as high speed digital oscilloscope and High Speed High Precision ADC circuit are in Electronic Design
In more and more extensive application, the test for High Speed High Precision ADC also increasingly attracts attention.For this kind of ADC, due to it
Application background is different from common ADC, and it more pays close attention to the lower kinematic nonlinearity (Dynamic presented of high-frequency signal input
Nonlinearity).Now, all kinds of evaluation methods having pointed out in early days are because a variety of conditions limit, as do not considered in model construction
Kinematic nonlinearity, ensure strict correlation sampling condition, need single Frequency Estimation accurate etc. enough, height can not be applied to
Fast high-precision adc is fast and accurately tested.
Therefore, being directed to above-mentioned restrictive condition, the present invention studies and improves traditional evaluation method, is ensureing acceptable survey
Realize that the nonlinear quick test of the lower High Speed High Precision ADC of high-frequency signal input has become current mixed signal and surveyed under accuracy of measurement
One of the hot research problem in examination field, checking and test of its achievement in research for analog-digital blended signal device have very heavy
The meaning wanted.
The content of the invention
It is an object of the invention to provide a kind of quick test side of High Speed High Precision ADC dynamic input-output characteristic curve
Method, in the case where ensureing enough estimation precisions, it substantially reduces the sampling number needed for test, and is adopted without strictly relevant
Batten part, greatly shortens the testing time.
To achieve the above object, the technical solution adopted by the present invention is:
A kind of High Speed High Precision ADC dynamic input-output characteristic curve method for rapidly testing, using the first kind and the second class
Chebyshev polynomials model to ADC dynamics input-output characteristic curve to be measured;Then on the basis of reference error optimal estimation
On the ADC of lower collection inputted using sinusoidal signal export digital code and known sinusoidal excitation signal information builds existing dynamic
Input-output characteristic curve model;The unknown parameter in the model is then identified in time domain using least-square fitting approach,
The final dynamic input-output characteristic curve for estimating fitting.
Concretely comprise the following steps:
If ADC to be measured pumping signal is that sinusoidal signal is as follows,
X (t)=V cos (2 π fxt)+C (1)
In formula, V is to input sinusoidal signal amplitude, fxTo input sinusoidal signal frequency, C is input sinusoidal signal offset;
So, the complete model of ADC dynamics input-output characteristic curve should be modeled as:
In formula, NhThe referred to as exponent number of the model, by ckThe sum term of composition is the description to static non linear, and by dkInstitute
The sum term of composition is the description to kinematic nonlinearity, and e (t) represents the summation of all systematic errors;
E (t) influence is put aside in subsequent algorithm estimation process, influences of the e (t) for estimation precision will optimize
Provided in estimation result, then the input-output characteristic curve now considered is actually that the smooth fitting of an approximation is bent
Line;It is more using the first kind and the second class Chebyshev with reference to the non-monodrome form for the dynamic input-output characteristic curve analyzed
The relation of Xiang Shiyu trigonometric functions, such as following two formula,
Tn(cos θ)=cos (n θ) (3)
In formula, TnFor Chebyshev polynomial of the first kind, UnFor Chebyshev polynomial of the second kind;
Mathematic(al) manipulation, which is carried out, with reference to above-mentioned formula (2), (3), (4) obtains the input-output characteristic curve g obtained by estimation1
(x), g2(x),
In formula, tkIt is to correspond to Chebyshev polynomial of the first kind coefficient, ukIt is corresponding Chebyshev polynomial of the second kind system
Number, coefficient c0With coefficient c in formula (2)0Unanimously;
The method estimated using parameter spectral line, the identification of spectral line component directly is carried out to ADC output data in time domain.
For spectrum analysis, it is by ADC output characteristics model modeling:
In formula, fsFor ADC sampling rate, fxFor frequency input signal, N is sampling number;
For frequency input signal unknown in formula (7), f is estimated using iterative methodx;
Then, as shown in formula (14) and formula (15), by the relation of Chebyshev polynomials coefficient and fourier coefficient, directly
Obtain to obtain dynamic input-output characteristic curve g1And g (x)2(x);
In formulaFor first phase deviation;As shown in (16), it can be estimated by once non-same phase harmonic component;
So far, this method of testing is achieved that quick obtaining dynamic input-output characteristic curve g1And g (x)2(x)。
3rd, High Speed High Precision ADC dynamic input-output characteristic curve method for rapidly testing according to claim 2, its
It is characterised by:F is estimated using iterative methodxSpecific implementation process it is as follows:
1) f is obtained using IpDFTxIteration preliminary examination value, it is assumed that its angular frequency w=2 π fx;
2) determineAngular frequency initial valuew0Afterwards, the estimation of primary parameter spectral line is completed using formula (7), so that it is determined that in formula (7)
Each fourier coefficient initial value;
3) the input signal angular frequency for assuming ith iteration is wi.According to Taylor expansion,
cos wtn≈coswitn-tnsinwitnΔwi (8)
sin wtn≈sinwitn+tncoswitnΔwi (9)
In formula,
Δwi=w-wi (10)
Formula (8) and formula (9) are substituted into formula (7), readjusting and simplifying obtains,
Y=Hixi (11)
In formula,
Obviously, sampling number N > > 2N under normal circumstancesh+ 2, therefore, formula (11) meets the spy of overdetermined linear system
Sign, estimate x is obtained by least-square fitting approachi;Iteration is continued thereafter with, until estimate xiMeet set essence
Degree, just obtains every Fourier coefficient and frequency input signal fxEstimate.
The method of testing quantitative relationship of analysis model error on sampling number and model order also in emulation, is obtained
Obtain error optimization estimation;Estimation precision of the user in actual test needed for setting measurement, come then referring to error optimization estimation
Complete ADC to be measured non-linear test.
Beneficial effect:
The advantage of the invention for reasonably utilizing time-domain analysis, and model and estimate using Chebyshev polynomials on the basis of
Mode fit ADC dynamic input-output characteristic curves, its implementation complexity is low, suitable for most of High Speed High Precision ADCs
Non-linear test.
The present invention significantly reduces the sampling number needed for test in the case where ensureing acceptable estimation precision, realizes fast
The purpose of speed test.
Brief description of the drawings
Fig. 1 is the major architectural that ADC is tested in the present invention;
Fig. 2 is ADC dynamics input-output characteristic curve method for rapidly testing broad flow diagram of the present invention;
Fig. 3 is digital-to-analogue mixed signal test system schematic diagram of the present invention;
Fig. 4 is the maximum dynamic I NL estimation errors of the present invention and model order and sampling number graph of a relation;
Fig. 5 is the present invention based on the dynamic I NL Error Graphs 1 for improving histogram test and this method of testing;
Fig. 6 is the present invention based on the dynamic I NL Error Graphs 2 for improving histogram test and this method of testing.
Embodiment
The technical scheme of invention is described in detail below in conjunction with the accompanying drawings.
Technical scheme is described in detail below in conjunction with the accompanying drawings.It should be understood that these embodiments are to be used for
Illustrate the present invention and be not limited to limit the scope of the present invention.The implementation condition used in embodiment can be according to the bar of specific producer
Part does further adjustment, and unreceipted implementation condition is usually the condition in normal experiment.
Embodiment:
Present embodiment describes a kind of method for rapidly testing of High Speed High Precision ADC dynamic input-output characteristic curve, Fig. 1
The major architectural that ADC is tested in the present embodiment is given, the sinusoidal signal that enough purity is produced by high-precision generator is used as
ADC to be measured pumping signal, while to the clock source of the ultralow shake of ADC to be measured offers, when ADC is in stable collection, by surveying
Test system collection ADC output digital codes, are finally based on this analysis of test methods to complete ADC dynamic input-output characteristic curves
Quick estimation.In the present embodiment, ADC to be measured resolution ratio is 14 bits, sampling rate 125MHz.Fig. 2 gives the present invention
ADC dynamic input-output characteristic curve method for rapidly testing broad flow diagrams.
1. the appropriate signal source of selection.High RST source precision is to ensure the key factor that ADC is accurately tested, when certain money of selection
, it is necessary to the signal source of appropriate choice accuracy during ADC to be measured.Generally selected signal source precision 2 ratio at least higher than ADC to be measured
It is more than spy.
2. build digital-to-analogue mixed signal test system.Possesses stability, the test system of high efficiency contributes to ADC's to be measured
Quick test, Fig. 3 give the structure of digital-to-analogue mixed signal test system in the present embodiment, and in whole system, power supply is adopted
0.016%+1.5mV Agilent N6705 are reached with voltage accuracy, test and excitation source device, which uses, possesses 16 bit accuracies
Agilent 33522A, Clock generation module are less than -140dBc/Hz NI PXIe 5450, data using phase noise density
Capture card uses the boards of NI PXIe 6556 of the internal high-speed figure input and output that can support 200MHz, and data processing and deposits
Store up by way of the Labview software platforms of NI companies call Matlab high speed data processing instruments to complete.
3. handled based on the digital code that this method of testing gathers to test system and analyze to obtain the non-linear of ADC to be measured
(dynamic input-output characteristic curve).The specific implementation process of this method of testing is as follows:
In the test for high-frequency input signal, during ADC input-output characteristic curve model constructions, not only to consider pure
The introduced same phase harmonic wave of static non linear, it is also contemplated that the non-same phase harmonic wave that kinematic nonlinearity is introduced.In addition, also wrap
Include the factors of influence such as thermal noise, Gaussian noise, quantizing noise, collectively referred to here in as systematic error.Assuming that ADC to be measured pumping signal
It is as follows for sinusoidal signal,
X (t)=V cos (2 π fxt)+C (1)
In formula, V is to input sinusoidal signal amplitude, fxTo input sinusoidal signal frequency, C is input sinusoidal signal offset.
So, the complete model of ADC dynamics input-output characteristic curve should be modeled as,
In formula, NhThe referred to as exponent number of the model, by ckThe sum term of composition is the description to static non linear, and by dkInstitute
The sum term of composition is the description to kinematic nonlinearity, and e (t) represents the summation of all systematic errors.
E (t) influence is put aside in subsequent algorithm estimation process, influences of the e (t) for estimation precision will optimize
Provided in estimation result.The input-output characteristic curve that so now this method of testing is considered is actually that an approximation is smooth
Matched curve.With reference to the non-monodrome form for the dynamic input-output characteristic curve analyzed, cut using the first kind and the second class
Than the relation of snow husband's multinomial and trigonometric function, such as formula (3) and formula (4),
Tn(cos θ)=cos (n θ) (3)
In formula, TnFor Chebyshev polynomial of the first kind, UnFor Chebyshev polynomial of the second kind.
Convolution (2), formula (3) and formula (4), which carry out mathematic(al) manipulation, can obtain the input-output characteristic curve g obtained by estimation1
(x), g2(x),
In formula, tkIt is to correspond to Chebyshev polynomial of the first kind coefficient, ukIt is corresponding Chebyshev polynomial of the second kind system
Number, coefficient c0With coefficient c in formula (2)0Unanimously.
In order to avoid the spectrum analysis of FFT on frequency domain, that is, the limitation of sampling clock stringent synchronization condition is broken away from, used
The method of parameter spectral line estimation, the identification of spectral line component directly is carried out to ADC output data in time domain.For spectral line point
Analysis, it is by ADC output characteristics model modeling,
In formula, fsFor ADC sampling rate, fxFor frequency input signal, N is sampling number.
For frequency input signal unknown in formula (7), best mode is to estimate f using iterative methodx.It is implemented
Process is as follows,
1) f is obtained using IpDFTxIteration preliminary examination value, it is assumed that its angular frequency w=2 π fx;
2) determineAngular frequency initial valuew0Afterwards, the estimation of primary parameter spectral line is completed using formula (7), so that it is determined that in formula (7)
Each fourier coefficient initial value;
3) the input signal angular frequency for assuming ith iteration is wi.According to Taylor expansion,
cos wtn≈coswitn-tnsinwitnΔwi (8)
sin wtn≈sinwitn+tncoswitnΔwi (9)
In formula,
Δwi=w-wi (10)
Formula (8) and formula (9) are substituted into formula (7), readjusting and simplifying obtains,
Y=Hixi (11)
In formula,
Obviously, sampling number N > > 2N under normal circumstancesh+ 2, therefore, formula (11) meets the spy of overdetermined linear system
Sign, can obtain estimate x by least-square fitting approach (Least Square Fitting, LS)i.Continue thereafter with repeatedly
Generation, until estimate xiMeet set precision, just obtain every Fourier coefficient and frequency input signal fxEstimation
Value.
Then,, can by the relation of Chebyshev polynomials coefficient and fourier coefficient as shown in formula (14) and formula (15)
Directly obtain dynamic input-output characteristic curve g1And g (x)2(x)。
In formulaFor first phase deviation.As shown in (16), it can be estimated by once non-same phase harmonic component.
So far, this method of testing is achieved that quick obtaining dynamic input-output characteristic curve g1And g (x)2(x)。
4. optimize estimation result.ADC quickly estimates that the maximum cost paid is the decline of test result precision.For
This, Fig. 4 gives maximum dynamic I NL estimation errors and the quantitative relationship of model order and sampling number.Using 14 bit A/D C as
Example, being usually only necessary to 8000 sampling numbers, 60 model orders just can greatly optimize the estimation result of this method of testing, thus
Had great practical value in actual test.
5. Fig. 5,6 are given based on the dynamic I NL Error Graphs for improving histogram test and this method of testing.Based on Fig. 5,6
Result, table 1 gives detailed test data contrast.Understand that this patent can be with by above test result figure and test data
The testing time is reduced, while ensures reliable measurement accuracy.
Table 1 is based on the dynamic I NL test results for improving histogram test and this method of testing
The foregoing examples are merely illustrative of the technical concept and features of the invention, its object is to allow the person skilled in the art to be
Present disclosure can be understood and implemented according to this, it is not intended to limit the scope of the present invention.It is all smart according to the present invention
The equivalent transformation or modification that refreshing essence is done, should all be included within the scope of the present invention.
Claims (4)
- A kind of 1. High Speed High Precision ADC dynamic input-output characteristic curve method for rapidly testing, it is characterised in that:Using the first kind ADC dynamics input-output characteristic curve to be measured is modeled with Chebyshev polynomial of the second kind;Then estimate reference error is optimal The ADC output digital codes of lower collection are inputted using sinusoidal signal and known sinusoidal excitation signal information is built on the basis of meter Some dynamic input-output characteristic curve models;Then identified using least-square fitting approach in time domain in the model not Know parameter, finally estimate the dynamic input-output characteristic curve of fitting.
- 2. High Speed High Precision ADC dynamic input-output characteristic curve method for rapidly testing according to claim 1, its feature It is:Concretely comprise the following steps:If ADC to be measured pumping signal is that sinusoidal signal is as follows,X (t)=Vcos (2 π fxt)+C (1)In formula, V is to input sinusoidal signal amplitude, fxTo input sinusoidal signal frequency, C is input sinusoidal signal offset;So, the complete model of ADC dynamics input-output characteristic curve should be modeled as:<mrow> <mi>y</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <msub> <mi>c</mi> <mn>0</mn> </msub> <mn>2</mn> </mfrac> <mo>+</mo> <msubsup> <mo>&Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>h</mi> </msub> </msubsup> <msub> <mi>c</mi> <mi>k</mi> </msub> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mrow> <mo>(</mo> <mi>k</mi> <mn>2</mn> <msub> <mi>&pi;f</mi> <mi>x</mi> </msub> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msubsup> <mo>&Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>h</mi> </msub> </msubsup> <msub> <mi>d</mi> <mi>k</mi> </msub> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mrow> <mo>(</mo> <mi>k</mi> <mn>2</mn> <msub> <mi>&pi;f</mi> <mi>x</mi> </msub> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <mi>e</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>In formula, NhThe referred to as exponent number of the model, by ckThe sum term of composition is the description to static non linear, and by dkFormed Sum term be description to kinematic nonlinearity, e (t) represents the summation of all systematic errors;E (t) influence is put aside in subsequent algorithm estimation process, influences of the e (t) for estimation precision will be estimated in optimization As a result provided in, then the input-output characteristic curve now considered is actually a smooth matched curve of approximation;Knot Close the non-monodrome form of analyzed dynamic input-output characteristic curve, using the first kind and Chebyshev polynomial of the second kind with The relation of trigonometric function, such as following two formula,Tn(cos θ)=cos (n θ) (3)<mrow> <msub> <mi>U</mi> <mi>n</mi> </msub> <mrow> <mo>(</mo> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&theta;</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mrow> <mo>(</mo> <mo>(</mo> <mrow> <mi>n</mi> <mo>+</mo> <mn>1</mn> </mrow> <mo>)</mo> <mi>&theta;</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>sin</mi> <mi>&theta;</mi> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>In formula, TnFor Chebyshev polynomial of the first kind, UnFor Chebyshev polynomial of the second kind;Mathematic(al) manipulation, which is carried out, with reference to above-mentioned formula (2), (3), (4) obtains the input-output characteristic curve g obtained by estimation1(x), g2 (x),<mrow> <msub> <mi>g</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <msub> <mi>c</mi> <mn>0</mn> </msub> <mn>2</mn> </mfrac> <mo>+</mo> <msubsup> <mo>&Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>h</mi> </msub> </msubsup> <msub> <mi>t</mi> <mi>k</mi> </msub> <msub> <mi>T</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>x</mi> <mo>-</mo> <mi>C</mi> </mrow> <mi>V</mi> </mfrac> <mo>)</mo> </mrow> <mo>+</mo> <msqrt> <mrow> <mn>1</mn> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>x</mi> <mo>-</mo> <mi>C</mi> </mrow> <mi>V</mi> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <msubsup> <mo>&Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>h</mi> </msub> </msubsup> <msub> <mi>u</mi> <mi>k</mi> </msub> <msub> <mi>U</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>x</mi> <mo>-</mo> <mi>C</mi> </mrow> <mi>V</mi> </mfrac> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow><mrow> <msub> <mi>g</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <msub> <mi>c</mi> <mn>0</mn> </msub> <mn>2</mn> </mfrac> <mo>+</mo> <msubsup> <mo>&Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>h</mi> </msub> </msubsup> <msub> <mi>t</mi> <mi>k</mi> </msub> <msub> <mi>T</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>x</mi> <mo>-</mo> <mi>C</mi> </mrow> <mi>V</mi> </mfrac> <mo>)</mo> </mrow> <mo>-</mo> <msqrt> <mrow> <mn>1</mn> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>x</mi> <mo>-</mo> <mi>C</mi> </mrow> <mi>V</mi> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <msubsup> <mo>&Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>h</mi> </msub> </msubsup> <msub> <mi>u</mi> <mi>k</mi> </msub> <msub> <mi>U</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>x</mi> <mo>-</mo> <mi>C</mi> </mrow> <mi>V</mi> </mfrac> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>In formula, tkIt is to correspond to Chebyshev polynomial of the first kind coefficient, ukIt is corresponding Chebyshev polynomial of the second kind coefficient, is Number c0With coefficient c in formula (2)0Unanimously;The method estimated using parameter spectral line, the identification of spectral line component directly is carried out to ADC output data in time domain.In order to Spectrum analysis, it is by ADC output characteristics model modeling:<mrow> <mi>y</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <msub> <mi>a</mi> <mn>0</mn> </msub> <mn>2</mn> </mfrac> <mo>+</mo> <msubsup> <mo>&Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>&infin;</mi> </msubsup> <msub> <mi>a</mi> <mi>k</mi> </msub> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mrow> <mo>(</mo> <mi>k</mi> <mfrac> <mrow> <mn>2</mn> <msub> <mi>&pi;f</mi> <mi>x</mi> </msub> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> <mi>n</mi> <mo>)</mo> </mrow> <mo>+</mo> <msubsup> <mo>&Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>&infin;</mi> </msubsup> <msub> <mi>b</mi> <mi>k</mi> </msub> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mrow> <mo>(</mo> <mi>k</mi> <mfrac> <mrow> <mn>2</mn> <msub> <mi>&pi;f</mi> <mi>x</mi> </msub> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> <mi>n</mi> <mo>)</mo> </mrow> <mo>,</mo> <mi>n</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>N</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>In formula, fsFor ADC sampling rate, fxFor frequency input signal, N is sampling number;For frequency input signal unknown in formula (7), f is estimated using iterative methodx;Then, as shown in formula (14) and formula (15), by the relation of Chebyshev polynomials coefficient and fourier coefficient, directly obtain Obtain dynamic input-output characteristic curve g1And g (x)2(x);In formulaFor first phase deviation;As shown in (16), it can be estimated by once non-same phase harmonic component;So far, this method of testing is achieved that quick obtaining dynamic input-output characteristic curve g1And g (x)2(x)。
- 3. High Speed High Precision ADC dynamic input-output characteristic curve method for rapidly testing according to claim 2, its feature It is:F is estimated using iterative methodxSpecific implementation process it is as follows:1) f is obtained using IpDFTxIteration preliminary examination value, it is assumed that its angular frequency w=2 π fx;2) angular frequency initial value w is determined0Afterwards, the estimation of primary parameter spectral line is completed using formula (7), so that it is determined that each Fu in formula (7) Vertical leaf system number initial value;3) the input signal angular frequency for assuming ith iteration is wi.According to Taylor expansion,cos wtn≈cos witn-tnsin witnΔwi (8)sin wtn≈sin witn+tncos witnΔwi (9)In formula,Δwi=w-wi (10)Formula (8) and formula (9) are substituted into formula (7), readjusting and simplifying obtains,Y=Hixi (11)In formula,<mrow> <mtable> <mtr> <mtd> <mrow> <mi>H</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mrow> <mi>cos</mi> <mrow> <mo>(</mo> <msub> <mi>w</mi> <mi>i</mi> </msub> <mo>/</mo> <msub> <mi>f</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mrow> <mi>sin</mi> <mrow> <mo>(</mo> <msub> <mi>w</mi> <mi>i</mi> </msub> <mo>/</mo> <msub> <mi>f</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>sin</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <mn>2</mn> <msub> <mi>w</mi> <mi>i</mi> </msub> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>h</mi> </msub> </munderover> <mrow> <mo>(</mo> <mo>-</mo> <mfrac> <mrow> <msub> <mi>a</mi> <mi>i</mi> </msub> <mi>k</mi> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> <mi>sin</mi> <mo>(</mo> <mfrac> <mrow> <msub> <mi>w</mi> <mi>i</mi> </msub> <mi>k</mi> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> <mo>)</mo> <mo>+</mo> <mfrac> <mrow> <msub> <mi>b</mi> <mi>i</mi> </msub> <mi>k</mi> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> <mi>cos</mi> <mo>(</mo> <mfrac> <mrow> <msub> <mi>w</mi> <mi>i</mi> </msub> <mi>k</mi> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> <mo>)</mo> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mn>2</mn> <msub> <mi>w</mi> <mi>i</mi> </msub> <mo>/</mo> <msub> <mi>f</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mrow> <mi>sin</mi> <mrow> <mo>(</mo> <mn>2</mn> <msub> <mi>w</mi> <mi>i</mi> </msub> <mo>/</mo> <msub> <mi>f</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>sin</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <mn>4</mn> <msub> <mi>w</mi> <mi>i</mi> </msub> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>h</mi> </msub> </munderover> <mrow> <mo>(</mo> <mo>-</mo> <mfrac> <mrow> <mn>2</mn> <msub> <mi>a</mi> <mi>i</mi> </msub> <mi>k</mi> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> <mi>sin</mi> <mo>(</mo> <mfrac> <mrow> <msub> <mi>w</mi> <mi>i</mi> </msub> <mi>k</mi> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> <mo>)</mo> <mo>+</mo> <mfrac> <mrow> <mn>2</mn> <msub> <mi>b</mi> <mi>i</mi> </msub> <mi>k</mi> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> <mi>cos</mi> <mo>(</mo> <mfrac> <mrow> <msub> <mi>w</mi> <mi>i</mi> </msub> <mi>k</mi> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> <mo>)</mo> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mrow> <mi>cos</mi> <mrow> <mo>(</mo> <msub> <mi>Nw</mi> <mi>i</mi> </msub> <mo>/</mo> <msub> <mi>f</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mrow> <mi>sin</mi> <mrow> <mo>(</mo> <msub> <mi>Nw</mi> <mi>i</mi> </msub> <mo>/</mo> <msub> <mi>f</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>sin</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>Nw</mi> <mi>i</mi> </msub> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>h</mi> </msub> </munderover> <mrow> <mo>(</mo> <mo>-</mo> <mfrac> <mrow> <msub> <mi>Na</mi> <mi>i</mi> </msub> <mi>k</mi> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> <mi>sin</mi> <mo>(</mo> <mfrac> <mrow> <msub> <mi>w</mi> <mi>i</mi> </msub> <mi>k</mi> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> <mo>)</mo> <mo>+</mo> <mfrac> <mrow> <msub> <mi>Nb</mi> <mi>i</mi> </msub> <mi>k</mi> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> <mi>cos</mi> <mo>(</mo> <mfrac> <mrow> <msub> <mi>w</mi> <mi>i</mi> </msub> <mi>k</mi> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> <mo>)</mo> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> </mrow><mrow> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>=</mo> <msup> <mrow> <mo>&lsqb;</mo> <mfrac> <msub> <mi>a</mi> <mn>0</mn> </msub> <mn>2</mn> </mfrac> <mo>,</mo> <msub> <mi>a</mi> <mn>1</mn> </msub> <mo>,</mo> <mo>...</mo> <mo>,</mo> <msub> <mi>a</mi> <msub> <mi>N</mi> <mi>h</mi> </msub> </msub> <mo>,</mo> <msub> <mi>b</mi> <mn>1</mn> </msub> <mo>,</mo> <mo>...</mo> <mo>,</mo> <msub> <mi>b</mi> <msub> <mi>N</mi> <mi>h</mi> </msub> </msub> <mo>,</mo> <msub> <mi>&Delta;w</mi> <mi>i</mi> </msub> <mo>&rsqb;</mo> </mrow> <mi>T</mi> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> </mrow>Obviously, sampling number N > > 2N under normal circumstancesh+ 2, therefore, formula (11) meets the feature of overdetermined linear system, passes through Least-square fitting approach obtains estimate xi;Iteration is continued thereafter with, until estimate xiMeet set precision, just obtain Obtain every Fourier coefficient and frequency input signal fxEstimate.
- 4. High Speed High Precision ADC dynamic input-output characteristic curve method for rapidly testing according to claim 1, its feature It is:The method of testing quantitative relationship of analysis model error on sampling number and model order also in emulation, obtain Error optimization is estimated;Estimation precision of the user in actual test needed for setting measurement, has come then referring to error optimization estimation Into ADC to be measured non-linear test.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710902634.XA CN107809247A (en) | 2017-09-29 | 2017-09-29 | A kind of High Speed High Precision ADC dynamic input-output characteristic curve method for rapidly testing |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710902634.XA CN107809247A (en) | 2017-09-29 | 2017-09-29 | A kind of High Speed High Precision ADC dynamic input-output characteristic curve method for rapidly testing |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107809247A true CN107809247A (en) | 2018-03-16 |
Family
ID=61584136
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710902634.XA Pending CN107809247A (en) | 2017-09-29 | 2017-09-29 | A kind of High Speed High Precision ADC dynamic input-output characteristic curve method for rapidly testing |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107809247A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109143034A (en) * | 2018-09-10 | 2019-01-04 | 上海华虹集成电路有限责任公司 | Chip ADC automatic performance test system and method |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103475369A (en) * | 2013-07-04 | 2013-12-25 | 东南大学 | High-precision ADC test method based on signal-source-error one-time correction and identification |
CN103529379A (en) * | 2013-07-04 | 2014-01-22 | 东南大学 | Selection method used for low-resolution signal source in high-precision ADC test |
-
2017
- 2017-09-29 CN CN201710902634.XA patent/CN107809247A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103475369A (en) * | 2013-07-04 | 2013-12-25 | 东南大学 | High-precision ADC test method based on signal-source-error one-time correction and identification |
CN103529379A (en) * | 2013-07-04 | 2014-01-22 | 东南大学 | Selection method used for low-resolution signal source in high-precision ADC test |
Non-Patent Citations (1)
Title |
---|
周凯: ""基于参数谱估计的高速高精度ADC测试算法研究与实现"", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109143034A (en) * | 2018-09-10 | 2019-01-04 | 上海华虹集成电路有限责任公司 | Chip ADC automatic performance test system and method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109444717B (en) | Novel on-chip S parameter error calibration method and device | |
US8310385B2 (en) | Systems and methods for vector-based analog-to-digital converter sequential testing | |
CN103457603B (en) | A kind of method that ADC dynamic parameter is tested based on average frequency spectrum | |
CN104111435B (en) | A kind of test fixture error elimination method | |
CN103475369B (en) | The high-precision ADC test of identification is disposably calibrated based on signal source error | |
CN107181489A (en) | A kind of analog-to-digital conversion calibration method and device | |
CN103529379B (en) | For the choosing method in low-resolution signal source in high-precision adc test | |
Ramos et al. | Recent developments on impedance measurements with DSP-based ellipse-fitting algorithms | |
CN105978643A (en) | Measurement method for time delay test of quadrature phase shift keying modulator | |
CN106063132A (en) | On-chip analog-to-digital converter (adc) linearity test for embedded devices | |
CN105471431A (en) | Method for testing differential linearity error and integral linearity error of analog to digital converter | |
CN107809247A (en) | A kind of High Speed High Precision ADC dynamic input-output characteristic curve method for rapidly testing | |
Silva et al. | Methods for fast characterization of noise and offset in dynamic comparators | |
CN106291102B (en) | A kind of Frequency Standard Comparison device and method | |
CN102324935A (en) | Zero crossing detection-based digital-to-analog converter static testing method and system | |
Jin et al. | Code-density test of analog-to-digital converters using single low-linearity stimulus signal | |
Jin et al. | Testing of precision DAC using low-resolution ADC with wobbling | |
Ireland et al. | A method for using Josephson voltage standards for direct characterization of high performance digitizers to establish AC voltage and current traceability to SI | |
CN107800434A (en) | A kind of fast A/D C method of testings based on parameter extraction | |
CN107786204A (en) | Digital analog converter parameter test system and method | |
CN205945703U (en) | Digital analog converter parameter testing system | |
CN106685423A (en) | Method for testing static parameters of analog-digital converters by aid of sinusoidal waves | |
Stefani et al. | Simple and time-effective procedure for ADC INL estimation | |
Björsell | Modeling analog to digital converters at radio frequency | |
CN111999587A (en) | Iterative vector error correction method for modulation domain network parameter test |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180316 |