Data Acquisition Part 2 pptx

25 273 0
Data Acquisition Part 2 pptx

Đang tải... (xem toàn văn)

Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống

Thông tin tài liệu

Data Acquisition 16 3. As a consequence, if one wants to add some WGN to increase performance by averaging, the choice is dictated by the number of samples that may be averaged. This is clearly suggested by the intersecting continuous lines in Fig. 18, and better illustrated by Fig. 19, in which ENOB is plotted as a function of n σ for fixed N . It is clear, for example, that for 4N = it is convenient 0.2 LSB n σ ≈ , etc. Quite surprisingly, the very frequent choice 0.5 n σ = is optimal only for N of the order of 13 2 . 0 2 4 6 8 10 12 14 16 18 7.6 7.8 8 8.2 8.4 8.6 X: 18 Y: 8.274 σ n =0.05 LSB log 2 N b e simulations approx 1 approx 2 Fig. 14. ENOB of an 8-bit linear DAS with input WGN ( 0.05 LSB n σ = ), as a function of the number N of the averaged samples. Noise, Averaging, and Dithering in Data Acquisition Systems 17 0 2 4 6 8 10 12 14 16 18 7 8 9 10 11 12 13 14 15 16 17 X: 18 Y: 10.92 σ n =0.3LSB log 2 N b e simulations approx 1 approx 2 Fig. 15. ENOB of an 8-bit linear DAS with input WGN ( 0.3 LSB n σ = ), as a function of the number N of the averaged samples. Data Acquisition 18 0 2 4 6 8 10 12 14 16 18 7 8 9 10 11 12 13 14 15 16 X: 18 Y: 15.26 σ n =0.5LSB log 2 N b e simulations approx 1 approx 2 Fig. 16. ENOB of an 8-bit linear DAS with input WGN ( 0.5 LSB n σ = ), as a function of the number N of the averaged samples. Noise, Averaging, and Dithering in Data Acquisition Systems 19 0 2 4 6 8 10 12 14 16 18 11.5 12 12.5 13 X: 18 Y: 12.59 σ n =0.1LSB log 2 N b e simulations approx 1 approx 2 Fig. 17. ENOB of a 12-bit linear DAS with input WGN ( 0.1 LSB n σ = ), as a function of the number N of the averaged samples. Data Acquisition 20 Fig. 18. Variation in the ENOB (with respect to the nominal resolution b) as a function of the number N of the averaged samples, for different values of input WGN ( n σ = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6 LSB). The figure compares the approximation given by (24) (approx. 1) with expression (25), in which the approximation (17) of ()g ⋅ is used (approx. 2). n σ [LSB] b Δ [bit] 0 0 0.1 0.59 0.2 1.50 0.3 2.92 0.4 4.92 0.5 7.48 Tab. 2. Maximum (asymptotic) increase of ENOB attainable by averaging, for given levels n σ of input WGN. Noise, Averaging, and Dithering in Data Acquisition Systems 21 Fig. 19. ENOB increase as a function of the input noise n σ , for fixed values of the number N of averaged samples. The maxima of the curves, and the typical values 0.4 LSB n σ = and 0.5 LSB n σ = are highlighted. 7. Conclusions The chapter examines the overall effect, in terms of effective resolution, of input noise and output averaging in linear DAS. The analysis applies to both the cases of unwanted system noise, and of noise purposely added to increase the performance (non-subtractive dithering). After a brief discussion of the ENOB figure of merit, the equations to determine the ENOB in various situations are derived and validated by simulations. The results clarify the nature of the acquisition error in presence of noise – in terms of “dithered quantization error” q d e and “randomized quantization error” q r e – and can be used, for example, to choose the optimal level of input noise in a non-subtractive dithering scheme. The choice is demonstrated to be non-trivial, even if quite simple with the use of the proper equations. In particular, the very common choice 0.5 LSB n σ = is demonstrated to be suboptimal in most practical cases. A very important warning is that the presented analysis is limited to the case of perfectly linear DAS, and is not applicable in the common case of meaningful nonlinearity error affecting the DAS. The case of non-subtractive dithering in nonlinear DAS can be analyzed with means similar to those presented in this chapter. In particular, the optimal levels of Data Acquisition 22 noise for nonlinear DAS are considerably higher than those derived for linear DAS [AGLS07]. This is, however, the subject of a possible future extended version of the chapter. 8. Acknowledgements The authors wish to thank prof. Mario Savino for helpful discussions and suggestions. 9. References [AD09] L. Angrisani and M. D’Arco. Modeling timing jitter effects in digital-to-analog converters. IEEE Trans. Instrum. Meas., 58(2):330–336, 2009. [AGLS07] F. Attivissimo, N. Giaquinto, A. M. L. Lanzolla, and M. Savino. Effects of midpoint linearization and nonsubtractive dithering in A/D converters. Measurement, 40(5):537–544, June 2007. [AGS04] F. Attivissimo, N. Giaquinto, and M. Savino. Uncertainty evaluation in dithered A/D converters. In Proc. of IMEKO TC7 Symposium, pages 121–124, St. Petersburg, Russia, June 2004. [AGS08] F. Attivissimo, N. Giaquinto, and M. Savino. Uncertainty evaluation in dithered ADC-based instruments. Measurement, 41(4):364–370, May 2008. [AH98] O. Aumala and J. Holub. Dithering design for measurement of slowly varying signals. Measurement, 23(4):271–276, June 1998. [BDR05] E. Balestrieri, P. Daponte, and S. Rapuano. A state of the art on ADC error compensation methods. IEEE Trans. Instrum. Meas., 54(4):1388–1394, 2005. [CP94] P. Carbone and D. Petri. Effect of additive dither on the resolution of ideal quantizers. IEEE Trans. Instrum. Meas., 43(3):389 –396, June 1994. [GT97] N. Giaquinto and A. Trotta. Fast and accurate ADC testing via an enhanced sine wave fitting algorithm. IEEE Trans. Instrum. Meas., 46(4):1020–1025, August 1997. [IEE94] IEEE Standards Board. IEEE Standard 1057 for Digitizing Waveform Recorders. IEEE Press, New York, NY, December 1994. [IEE00] IEEE Standards Board. IEEE Standard 1241 for Terminology and Test Methods for Analog-to-Digital Converters . IEEE Press, New York, NY, December 2000. [KB05] I. Kollár and J. J. Blair. Improved determination of the best fitting sine wave in ADC testing. IEEE Trans. Instrum. Meas., 54:1978–1983, October 2005. [Nat97] National Instruments, Inc. PCI-1200 User Manual, January 1997. [Nat05] National Instruments, Inc. PXI-5922 Data Sheet, 2005. [Nat07] National Instruments, Inc. DAQ E-Series User Manual, February 2007. [Sch64] L. Schuchman. Dither signals and their effect on quantization noise. IEEE Trans. Comm. Tech. , 12(4):162–165, December 1964. [SO05] R. Skartlien and L. Oyehaug. Quantization error and resolution in ensemble averaged data with noise. IEEE Trans. Instrum. Meas., 54(3):1303 – 1312, June 2005. [WK08] B. Widrow and I. Kollár. Quantization Noise: Roundoff Error in Digital Computation, Signal Processing, Control, and Communications . Cambridge University Press, Cambridge, UK, 2008. [WLVW00] R. A. Wannamaker, S. P. Lipshitz, J. Vanderkooy, and J. N. Wright. A theory of nonsubtractive dither. IEEE Trans. Signal Process., 48(2):499–516, 2000. 2 Bandpass Sampling for Data Acquisition Systems Leopoldo Angrisani 1 and Michele Vadursi 2 1 University of Naples Federico II, Department of Computer Science and Control Systems 2 University of Naples “Parthenope”, Department of Technologies Italy 1. Introduction A number of modern measurement instruments employed in different application fields consist of an analogue front-end, a data acquisition section, and a processing section. A key role is played by the data acquisition section, which is mandated to the digitization of the input signal, according to a specific sample rate (Corcoran, 1999). The choice of the sample rate is connected to the optimal use of the resources of the data acquisition system (DAS). This is particularly true for modern communication systems, which operate at very high frequencies. The higher the sample rate, in fact, the shorter the observation interval and, consequently, the worse the frequency resolution allowed by the DAS memory buffer. So, the sample rate has to be chosen high enough to avoid aliasing, but at the same time, an unnecessarily high sample rate does not allow for an optimal exploitation of the DAS resources. As well known, the sample rate must be correctly chosen to avoid aliasing, which can seriously affect the accuracy of measurement results. The sampling theorem, in fact, affirms that a band-limited signal can be alias-free sampled at a rate f s greater than twice its highest frequency f max (Shannon, 1949). As regards bandpass signals, which are characterized by a low ratio of bandwidth to carrier frequency and are peculiar to many digital communication systems, a much less strict condition applies. In particular, bandpass signals can be alias-free sampled at a rate f s greater than twice their bandwidth B (Kohlenberg, 1953). It is worth noting, however, that this is only a necessary condition. It is indeed possible to alias-free sample bandpass signals at a rate fs much lower than 2f max, but such rate has to be chosen very carefully; it has been shown in (Brown, 1980; Vaughan et al., 1991; De Paula & Pieper, 1992; Tseng, 2002) that aliasing can occur if fs is chosen outside certain ranges. Moreover, particular attention has to be paid, as bandpass sampling can imply a degradation of the signal-to-noise ratio (Vaughan et al., 1991). Some recent papers have also focused on frequency shifting induced by bandpass sampling in more detail (Angrisani et al. 2004; Diez et al., 2005), providing analytical relations for establishing the final central frequency of the discrete-time signal, which digital receivers need to know (Akos et al., 1999) and determining the minimum admissible value of fs that is submultiple of a fixed sample rate (Betta et al., 2009). Sampling a bandpass signal at a rate lower than twice its highest frequency f max is referred to as bandpass sampling. Bandpass sampling is relevant in several fields of application, such Data Acquisition 24 as optics (Gaskell, 1978), communications (Waters & Jarrett, 1982), radar (Jackson & Matthewson, 1986) and sonar investigations (Grace & Pitt, 1968). It is also the core of the receiver of software-defined radio (SDR) systems (Akos et al., 1999; Latiri et al., 2006). Although the theory of bandpass sampling is now well-established and the choice of sample rate is very important for processing and measurement, at the current state of the art it seems that digital instruments that automatically select the best f s , on the basis of specific optimization strategies, are not available on the market. A possible criterion for choosing the optimal value of f s within the admissible alias-free ranges was introduced some years ago (Angrisani et al., 2004). An iterative algorithm was proposed, which selects the minimum alias-free sample rate that places the spectral replica at the normalized frequency requested by the user. The algorithm, however, cannot be profitably applied to any DAS. Two conditions have, in fact, to be met: (i) the sample rate can be set with unlimited resolution, and (ii) the sample clock has to be very stable. Failing to comply with such ideal conditions may result in an undesired and unpredictable frequency shifting and possible aliasing. More recently, a comprehensive analysis of the effects that the sample clock instability and the time-base finite resolution have on the optimal sample rate and, consequently, on the central frequency of the spectral replicas was developed (Angrisani & Vadursi, 2008). On the basis of its outcomes, the authors also presented an automatic method for selecting the optimal value of f s , according to the aforementioned criterion. The method includes both sample clock accuracy and time-base resolution among input parameters, and is suitable for practical applications on any DAS, no matter its sample clock characteristics. Specifically, the method provides the minimum f s that locates the spectrum of the discrete-time signal at the normalized central frequency required by the user, given the signal bandwidth B, a possible guard band B g , and original carrier frequency f c . Information on the possible deviation from expected central frequency, as an effect of DAS non-idealities, is also made available. In fact, the proposed method is extremely practical, since (i) it can be profitably applied no matter what the time-base resolution of the DAS is, and (ii) it takes into account the instability of the sample clock to face unpredictable frequency shifting and the consequent possible uncontrolled aliasing. A number of tests are carried out to assess the performance of the method in correctly locating the spectral replica at the desired central frequency, while granting no superposition of the replicas. Some tests are, in particular, mandated to highlight the effects of DAS non-idealities on the frequency shifting and consequent unexpected aliasing. This chapter is organized as follows. The theory of bandpass sampling will be presented in Section 2, along with analytical relations for establishing the final central frequency of the discrete-time signal and details and explicative figures on the frequency shifting resulting from the bandpass sampling and on the effects of the sample rate choice in terms of possible aliasing. Section 2 also presents the analysis of the effects that the sample clock instability and the time-base finite resolution which was first introduced in (Angrisani & Vadursi, 2008). Section 3 presents the proposed algorithm for the automatic selection of the sample rate given the user’s input, and shows the results of experiments conducted on real signals. 2. Analysis of the effects of bandpass sampling with a non-ideal data acquisition system Let s(t) be a generic bandpass signal, characterized by a bandwidth B and a central frequency f c . As well known, the spectrum of the discrete-time version of s(t) consists of an infinite set of replicas of the spectrum of s(t), centered at frequencies Bandpass Sampling for Data Acquisition Systems 25 f λ,ν = λ f c + ν f s (1) where ν ∈ Z and λ ∈ {-1;1}. The situation is depicted in Fig.1 with regard to positive frequencies of magnitude spectrum. Fig. 1. Typical amplitude spectrum of (a) a bandpass signal s(t) and (b) its sampled version; f s is the sample rate. Only the positive portion of the frequency axis is considered. Replicas of the 'positive' spectrum (red triangles in Fig.1) are centered at f 1,ν , whereas those peculiar to the 'negative' one (white triangles in Fig.1) are centered at f −1,ν . It can be shown that only two replicas are centered in the interval [0, f s ], respectively at frequencies f λ1,ν1 = f c mod f s (2) and f λ2,ν2 = f s –( f c mod f s ) (3) where mod denotes the modulo operation. The condition to be met in order to avoid aliasing is ( ) 22 s g g f BB BB f ∗ −+ + << (4) where f* is the minimum f λ1,ν1 and f λ2,ν2 . Inequality (4) implies the following condition on (f c mod f s ): ( ) mod , , 22 2 2 sg gsg cs s fBB BB f BB B ff f ⎛⎞ −+ +++ ⎛⎞ ⎜⎟ ∈∪− ⎜⎟ ⎜⎟ ⎜⎟ ⎝⎠ ⎝⎠ (5) [...]... 11.811 02 12. 00000 12. 09677 12. 29508 12. 39669 12. 60504 λ 1 -1 1 -1 1 -1 1 -1 f* [MHz] 3.84615 3.90 625 3.93700 4.00000 4.0 322 5 4.09836 4.1 322 3 4 .20 168 f1,old [MHz] 1. 920 00 1. 920 00 1. 920 00 1. 920 00 1. 920 00 1. 920 00 1. 920 00 1. 920 00 f2,old [MHz] 3.84 923 3.93938 3.98551 4.08000 4. 128 39 4 .22 754 4 .27 835 4.3 825 2 f1 [MHz] 2. 09585 2. 09859 2. 09581 2. 098 62 2.09577 2. 09865 2. 09573 2. 09868 f2 [MHz] 3.67133 3.7 628 5 3.80760... 1 :2: 5 rule 38 Data Acquisition Bg [MHz] 1.16 16.30435 5.43485 0. 12 24.19355 8.06455 0. 12 12. 6050 4 .20 00 0.15 1.16 16.3043 5.4333 0 .20 24 .1935 8.0635 0 .20 8.438 82 2.10970 0.09 1.16 10.81081 2. 7 027 4 0.19 3.84 16 .26 016 4.06496 0.19 0 Δf = 100 Hz 0.10 0 p=4 4 .20 160 3.84 Δf = 10 Hz 12. 60504 0 Δf = 100 Hz Δf* [%] 3.84 p=3 f* [MHz] 0 Δf = 10 Hz fs [MHz] 8.4388 2. 1108 0 .27 1.16 10.8108 2. 70 32 0.33 3.84 0 .21 ... Δf = 100 Hz Δf* [%] 3.84 p=3 f* [MHz] 0 Δf = 10 Hz fs [MHz] 8.4388 2. 1108 0 .27 1.16 10.8108 2. 70 32 0.33 3.84 0 .21 1.16 21 . 621 62 2.7 027 4 0.18 33.05785 4.1 322 5 0 .24 17.1674 2. 1454 0 .25 1.16 21 . 621 6 2. 70 32 0 .25 3.84 Δf = 100 Hz 0 .25 2. 14598 0 p=8 4.06 62 17.16738 3.84 Δf = 10 Hz 16 .26 02 0 33.0579 4.1315 0.39 Table 5 Frequency-domain results related to a QAM signal with bandwidth B = 3.84 MHz and a carrier... Δf ⎞ , −χM ⎟ ⎜− 2 ⎝ ⎠ Δf ⎤ ⎛1 ⎞⎡ ⎜ 2 − ν ⎟ ⎢ − χ M fs − 2 ( 1 − χ M )⎥ ⎝ ⎠⎣ ⎦ B ⎛ Δf ⎞ , χM ⎟ ⎜− ⎝ 2 ⎠ Δf ⎤ ⎛1 ⎞⎡ ( 1 + χ M )⎥ ⎜ −ν ⎟ ⎢ χ M fs − 2 2 ⎠⎣ ⎦ C ⎛ Δf ⎞ , χM ⎟ ⎜ ⎝ 2 ⎠ Δf ⎤ ⎛1 ⎞⎡ ( 1 + χ M )⎥ ⎜ −ν ⎟ ⎢ χ M fs + 2 2 ⎠⎣ ⎦ D ⎛ Δf ⎞ , −χM ⎟ ⎜ ⎝ 2 ⎠ (15) Δf ⎤ ⎛1 ⎞⎡ ⎜ 2 − ν ⎟ ⎢ − χ M fs + 2 ( 1 − χ M )⎥ ⎝ ⎠⎣ ⎦ Table 2 Values assumed by g2(ε, χ) in the vertices of its domain D 2. 1 Replica of the... Acquisition ⎧ ⎡ B + Bg ⎤ Δf p +ν ⎪ fs > ( 1 − χ M )⎥ ⎢ 1 − p ν χM ⎢ 2 2 ⎥ ⎪ ⎣ ⎦ ⎪ ⎨ 1⎞ ⎛ B + Bg + Δf ⎜ ν − ⎟ ( 1 + χ M ) ⎪ 2 ⎝ ⎪f* > p ⎪ 2 + p ( 1 + χ M ) − 2 p ν χ M ⎩ (22 ) when λ = -1 (negative replica) Once the user has entered the desired value of p, the algorithm provides the lowest fs that verifies (20 ) and (21 ) (or (20 ) and (22 ), if λ = -1), given the bandwidth and the central frequency of the... minimized by Δf ⎞ ⎛ ⎜ f s − 2 ⎟ ( 1 − χ M ) − B + Bg ⎡ Δf ( 1 − χ M ) ⎤ ⎝ ⎠ +ν ⎢ χ M fs + ⎥ 2 2 ⎢ ⎥ ⎣ ⎦ ( ) (14) 28 Data Acquisition According to (13) and (14), in the most restrictive case the condition (9) can be rewritten as ⎧ ⎪ ⎪f*< ⎪ ⎨ ⎪ ⎪f* > ⎪ ⎩ Δf ⎞ ⎛ ⎜ f s − 2 ⎟ ( 1 − χ M ) − B + Bg ⎡ Δf ( 1 − χ M ) ⎤ ⎝ ⎠ + ν ⎢ χ M fs + ⎥ 2 2 ⎢ ⎥ ⎣ ⎦ B + Bg Δf ⎡ ⎤ −ν ⎢ χ M fs + ( 1 + χ M )⎥ 2 2 ⎣ ⎦ ( ) Vertex Coordinates... ε ( 1 + χ ) ⎤ ⎣ ⎦ (10) 27 Bandpass Sampling for Data Acquisition Systems ⎛1 ⎞ g2 ( ε , χ ) = ⎜ − ν ⎟ ⎡ χ f s + ε ( 1 + χ ) ⎤ ⎣ ⎦ 2 ⎠ (11) To find the pair {ε, χ} that maximizes g1(ε, χ) in the domain D = [-Δf /2, Δf /2] x [χM, χM] let us first null the partial derivatives of g1 with respect to variables ε and χ : ⎧ ⎪ −ν ( 1 + χ ) = 0 ⎨ ⎪−ν ( f s + ε ) = 0 ⎩ ( 12) The system ( 12) has no solutions in... A ⎛ Δf ⎞ , −χM ⎟ ⎜− ⎝ 2 ⎠ Δf ⎡ ⎤ −ν ⎢ − χ M f s − ( 1 − χ M )⎥ 2 ⎣ ⎦ B ⎛ Δf ⎞ , χM ⎟ ⎜− ⎝ 2 ⎠ Δf ⎡ ⎤ −ν ⎢ χ M f s − ( 1 + χ M )⎥ 2 ⎣ ⎦ C ⎛ Δf ⎞ , χM ⎟ ⎜ 2 ⎝ ⎠ Δf ⎡ ⎤ −ν ⎢ χ M f s + ( 1 + χ M )⎥ 2 ⎣ ⎦ D ⎛ Δf ⎞ , −χM ⎟ ⎜ 2 ⎝ ⎠ Δf ⎡ ⎤ −ν ⎢ − χ M f s + ( 1 − χ M )⎥ 2 ⎣ ⎦ Table 1 Values assumed by g1(ε, χ) in the vertices of its domain D B + Bg 2 Δf ⎡ −ν ⎢ χ M fs + ( 1 + χ M )⎤ ⎥ 2 ⎣ ⎦ (13) With regard... sample rate Δf 100 Hz 10 Hz 1 Hz 10-8 998.3 740 .29 704.369 10-6 1003.4 745.85 710 .25 3 5.47·10-6 1 024 .1 770.78 736.445 10-5 1045.7 797.64 763.857 10-4 1494.9 1319 .26 129 9.107 3.54·10-4 29 15.8 28 19.14 28 19.143 χM Table 4 Optimal sample rate, expressed in kilohertz, for a bandpass signal characterized by a bandwidth B = 140 kHz and a carrier frequency fc = 595. 121 MHz, as a function of different values of... similarly shown that g2(ε, χ) assumes its minimum on one of the vertices of the domain D The four alternatives are enlisted in Table 2 Being ν < 0, the vertex C can be discarded, because g2(C) is sum of all positive terms Moreover, as χM (fs - Δf /2) > 0, g2(A) < g2(B) Finally, posing g2(A) < g2(D) implies -Δf (1 - χM) < 0, which is always true in actual situations In conclusion, g2 assumes its minimum . 12. 00000 -1 4.00000 1. 920 00 4.08000 2. 098 62 3.90350 12. 09677 1 4.0 322 5 1. 920 00 4. 128 39 2. 09577 3.95046 12. 29508 -1 4.09836 1. 920 00 4 .22 754 2. 09865 4.05106 12. 39669 1 4.1 322 3 1. 920 00 4 .27 835. [MHz] f 2, old [MHz] f 1 [MHz] f 2 [MHz] 11.53846 1 3.84615 1. 920 00 3.84 923 2. 09585 3.67133 11.71875 -1 3.90 625 1. 920 00 3.93938 2. 09859 3.7 628 5 11.811 02 1 3.93700 1. 920 00 3.98551 2. 09581. Dithering in Data Acquisition Systems 19 0 2 4 6 8 10 12 14 16 18 11.5 12 12. 5 13 X: 18 Y: 12. 59 σ n =0.1LSB log 2 N b e simulations approx 1 approx 2

Ngày đăng: 21/06/2014, 01:20

Tài liệu cùng người dùng

  • Đang cập nhật ...

Tài liệu liên quan