Uncertainty of Stability Variances Based on Finite Differences
|
|
- Reginald Carpenter
- 5 years ago
- Views:
Transcription
1 Uncertainty of Stability Variances Based on Finite Differences C. A. Greenhall* Jet Propulsion Laboratory California Institute of Technology W. J. Riley S ymme t r icom, Inc. Draft 009 Abstract This work gives an algorithm for computing the degrees of freedom of estimators of Allan and Hadamard variances, including their modified versions. A consistent approach is used throughout. Introduction This work gives an algorithm by which one can compute error bars for measurements of frequency stability variances in the presence of powerlaw phase noises. These stability variances fall into two categories: unmodified variances, which use dth differences of phase samples for d = (Allan) or (Hadamard), and modified variances, which use dth differences of averaged phase samples. The corresponding sampling functions that act on phase and frequency are shown in [l]. Each variance is defined as a scaling factor times the expected value of the squared differences. Unbiased estimates of this variance are computed from available phase data by taking time averages of the squared differences. The usual choices for the estimation stride (the time step) are the sample period 70 and the averaging period 7, a multiple of 70. These give respectively the overlapped estimator (OE) and nonoverlapped estimator (NOE) of the stability variance (although nonoverlapped is a misnomer; there is always some overlap between summands). The new algorithm, which computes the equivalent degrees of freedom () of a variance estimator, can replace several ones currently in use with a single, complete, and consistent method. Specifically, this algorithm covers the OE and NOE of the unmodified and modified Allan and Hadamard variances for all common noise types (4 5 a 5 ) at all applicable sample sizes and averaging factors. Previously, for the NOE of Allan and Hadamard variances, sigma confidence limits were generally set by scaling the measured deviation by a noisedependent factor K,/a, where M is the number of summands [][6]. For the OE of Allan variance, empirical approximations [7] were generally used along with chisquared statistics; nonempirical methods for both Allan variance estimators were also published [][4]. Although an computation for the OE of modified Allan (and time) variance was first given in [5], the simpler approach of [6] was generally used; an alternative unpublished approach using the Hadamard formed the basis of the new The work of this author was performed at the Jet Propulsion Laboratory, California Institute of Technology ( JPL ), under a contract with the National Aeronautics and Space Administration. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise, does not constitute or imply its endorsement by the United States Government or JPL.
2 algorithm. As an example of the application of the new algorithm, the Stable program for frequency stability analysis [7] has adopted it (since Version.4) instead of the multiple algorithms previously used for setting confidence limits and error bars. We wish to maintain a clear distinction between a stability variance and its estimators; in particular, we treat both the OE and the NOE of each variance. Even though the OE usually has lower uncertainty than the NOE, the NOE is convenient when phase data are processed in real time or read sequentially from a file. In fact, the TSC 50A Time Interval Analyzer [8], in its averaging mode, computes an NOE of modified Allan variance. The goal of the new algorithm is not closed formulas but fast numerical results whose accuracy is adequate for the purpose at hand. All the calculations are based on the same theoretical principles, with no empirical formulas. For each T, one must choose a dominant phasenoise power law, S, ( f ).fa where a E {,,0,,,, 4) (white PM to random run FM); see [9] for a method of powerlaw identification. The phase noise is assumed to be approximately bandlimited to the Nyquist frequency / (70). Not covered are effects of trend removal, drift removal for Allan variance in particular; the dth phase differences are assumed to have zero mean. One can use Hadamard variance to obtain stability results that are invariant to linear frequency drift. Special longterm stability statistics, such as total Allan variance [0], total Hadamard variance [9], and Theol [], are not covered; these require their own treatments. N Theory of operation Although the presentation of the algorithm is selfcontained, a brief account of the theory behind it may be helpful. The algorithm s output is the equivalent degrees of freedom () of an unbiased estimator V of a stability variance u = EV. Define It has been observed empirically (but not systematically) that ( v / u )V has approximately a xz distribution. Thus, having computed v and observed V, one can obtain confidence intervals of form vv/xz 5 cr 5 v V / z ~from x$ levels < x [7]. The model for phase x ( t )is the 70difference of a pure powerlaw process: x ( t )= ATOW ( t )7 () where w ( t )is a continuoustime process with spectral density for all f > 0, and A is the backward difference operator. Then S, ( f ) is asymptotically proportional to fa as f 0 and has approximate bandwidth / (70); this is the first reason for using w ( t ). Now let ( t )= A$A ( t ), () d where E = TO or T. For E = TO we have z ( t )= A,z ( t ),which leads to an unmodified variance. For E =7 =m~0 we observe from () that A,w ( t )= w ( t ) w ( t TO) = m A,w ( t TO) = m (t TO) = me ( t ), n=o n=o where 5 ( t ) is a discretetime average of m samples of x. In this case, z ( t ) = ma$ ( t ),which leads to a modified variance [6]; this is the second reason for using w ( t ). In either case we
3 assume that z ( t )is a stationary zeromean Gaussian process with autocovariance (ACV) function S, ( t ) = Ez (U t ) z (u). Ignoring the conventional scaling factors, we define the stability variance and its estimator by where the stride S is 70for the OE and 7 for the NOE. The number of terms M depends on the estimator type and the number of data. We have EV = a = s, (0). Then cov [ z ( t ),z (u)]= s; (u t ), and The definition ()) after substitution of ( 5 ), simplifies to The ACV sz ( t )is obtained from () by applying a difference operator of order d generalized autocovariance (GACV) s, ( t )of the powerlaw process w (t) []: sz d ( t )= (ATAT) A,A,s, to the (t). The function s, ( t )is given below [6] for each a. Algorithm for calculation Our purpose is to obtain practical numerical approximations of (6). We give two versions of the algorithm: the simplified version merely truncates the sum in the exact formula; the full version maintains the number of summation terms below a presassigned threshold and avoids catastrophic roundoff errors. They have the same inputs, output, function definitions, and initial step. Some explanation of the approximations is given in Appendix A. Because the results are invariant to time scaling, we may set 7 =, 70= l/m. All arithmetic is to be carried out in doubleprecision floating point. Operations that give signed integers are the floor function x (greatest integer that is 5 x) and ceiling function.r = xj (least integer that is x).. Q Inputs = frequency noise exponent a =,, 0,,,, 4 Noise type = WHPM, FLPM, WHFM, FLFM, RWFM, FWFM, RRFM d = order of phase difference d = : firstdifference variance (included for completeness) d = : Allan variance d = : Hadamard variance Restriction: a d > m = averaging factor 7 / 7 0, positive integer
4 F = filter factor F = : modified variance F = m: unmodified variance S = stride factor (estimator stride = T / S ) S = : nonoverlapped estimator S = m: overlapped estimator N = number of phase data with sample period. TO Output = equivalent degrees of freedom of the variance estimator. Constant and function definitions Set an integer constant J,, (used only in the full version); we suggest J,, = 00. The formal arguments of the following functions have the same names as the input arguments of the main algorithm.. Define the function s, ( t,a ) as follows: a,s ( t, a ) ~tl tn It t4in ~tl [ti5 t6n It 0 lt 4 lt7. (7) The entries with In It( must evaluate to 0 when t = 0.. Define the function s, ( t,00, a ) =,s ( t,a ), 4 5 a Define the function s, ( t,f, a, d ) = ( A I A ~ ) (~t,f s,a ), d =,, ; (9) that is (with dependence on F and a suppressed on the right sides), sz ( t,f, a, ) = 6 s x ( t ) 4sZ ( t ) 4 s X (t SZ ) S, (t ) sx (t ), ( t,f, Q, ) = Osx ( t ) 5sZ ( t ) 5sX (t ) 6s, (t ) 6sX(t ) S, (t ) S, ( t 4. Define the function ( BasicSum(J,M,S,F,a,d)=s;(O,F,a,d) J j=l ( 4 G s, &) ($,F, a, s; ),F,a,d d).
5 .4 Initial steps for both versions. Compute M, the number of summands in the estimator, as follows: if N L, otherwise there are not enough data.. Let J = min ( M,(d )S). Main procedure, simplified version.5 This is just one step: _ BasicSum ( J,M, S, F, a,d ). sz (0, F, a, d ) M () To check the effect of the truncation, one can also,try a larger value of J, say min ( M,6s). Main procedure, full version.6 M Let r =. S There are four cases. The calculations use coefficients from three numerical tables..6. Case. Modified variances: F =, all a This case also applies to unmodified variances when F = m =. If J 5 Jm, _ BasicSum (J,M, S,,a,d ) sz (0,, a, d ) M Else if J > J,, and r d, take ao,al from Table ; then r Jmax Else let m' = (not necessarily an integer); then r _ BasicSum (J", sz (0,, a,d ) Jmm.6. Jmax, m',,a, d ) Case. Unmodified variances, WHFM to RRFM: F = m, a If J I Jmax If m ( d ) 5 J,, then let m' = m else let m' = 00. Then sz (0, m', a, d ) M BasicSum ( J,M, S, m', a, d ) 5 I0
6 Else if J > Jmax and T d, take ao,al from Table ; then r Jmax Else let m' = (not necessarily an integer); then r SP BasicSum (Jmax, Jmax, m', 00, Q, d ) (O,OO, a, d) Jm,.6. Case. Unmodified variances, FLPM: F = m, If J 5 Jmax I _ s? (0, m,,d ) M Q = BasicSum (J,M, S,m,,d ) Remark: For this case, m must be less than about lo6 to avoid roundoff error. Else if J > Jmax and r d, take ao,al from Table (Q = ),bo, bl from Table ; then Jmax (not necessarily an integer); then Else let m' = r _ (bo bl In m)' Jmax BasicSum (J",.6.4 Case 4. Unmodified variances, WHPM: F = m, Jm,, Q m', m',,d ) = This calculation is exact, and can be expressed in closed form. In these formulas, binomial coefficient Let K = IfKSd r.. n! IC! (n I C )!. r Else M where a0 = also given in Table ( a = ). 6 d ( denotes the
7 .7 Tables Table. Coefficients for modified variances d=l d= d= I I Table. Coefficients for unmodified variances d=l d= d= a o $ a z 5. B 4 a a a0 a % z Table. Coefficients for logarithmic denominator, unmodified variances, FLFM ( a = ) d=l bo bl d= bo d= bl bo bl 40 Examples First, we must point out that the new values differ somewhat from older ones because of our choice of phase noise model. Previously (for a < ) the phase was generally assumed to have a pure f a spectrum; here the phase is modeled as the first difference of an process. For example, consider overlapped Allan variance, white FM, 05 phase data. The old results are from [4] (close to those of [7]). 7/70 old new The old and new results are in practical agreement at the larger values of r, where the results matter more. By allowing this mild discrepancy, we were able to make the algorithm simpler and more uniform. Figures and (at the end of the paper) show examples of, computed by the new algorithm, as a function of noise type and of variance type with other parameters fixed. 7
8 5 Conclusion The stability variances considered here can all be put into the same form, namely, the meansquare average of the output of a finitedifference filter acting, not on the phase samples, but on their cumulative sums. With this insight, we have been able to make an algorithm for computing the equivalent degrees of freedom of variance measurements. It covers all the commonly used variances and estimators, and then some. There is now a single unified source for these uncertainty calculations instead of the many papers that each cover only one variance and perhaps only one estimator of it. Complete pseudocode for the new algorithm is given here; software is available on request []. A Appendix: explanat ion of approximations As is, ( 6 ) is unfit for numerical computation. We find empirically that sp ( t )tends rapidly to zero as t increases beyond d. For the accuracy needed here (a few percent), there is no point in allowing j / S to be more than d. Indeed, for sufficiently large t the calculation of sp ( t )blows up from roundoff error, even in double precision, because linear combinations of large sw values are taken to get small s, values. At the very least, one should truncate the sum at j = (d )S, as in the simplified version of the algorithm. we do the The full version of the algorithm uses the following general strategy. If J 5 J,, there are two cases. First, if A4 (d )S then S = m summation ( 6 ). When J > J,, J,,/ (d ) >>. We truncate the sum at ( d )S and approximate it by an integral; this gives dl vd r 0 = (ao r (l;)s:(t)dt?), where These coefficients can be evaluated in advance. Second, if M < (d )S then we do another summation in which J is reduced from M to Jmax and S is reduced proportionately from m. The extra term for j = J in BasicSum makes the sum a trapezoidal approximation to the integral, whether or not the sum is truncated. This method works straightforwardly for Case ; indeed, in this setting the modified variances are simpler than the unmodified ones. In Case, when m is large we compute s, ( t ) using the limiting form s, ( t,m), which is actually sl ( t ). This means that we are treating z ( t ) as w' ( t ), the process w ( t )being differentiable in the meansquare sense. The most troublesome case is the overlapped estimators of the unmodified variances for flicker PM. As S = m 00, s, ( t ) approaches a function with logarithmic singularities. The factor bo bl l n m is an asymptotic form of s, (0). It would be possible (though inconvenient) to add another largem subcase as in Case, but one does not expect flicker PM to be the dominant noise type when m is large. Case 4 is constructed by knowing that the phase samples are accurately uncorrelated when w ( t ) is a Wiener process. The simplified computation () is correct but wasteful because s, ( t ) is a linear combination of hatshaped peaks of width /m. 8
9 References F. Vernotte, 00, Application of the moment condition to noise simulation and to stability analysis, IEEE Transactions on Ultrasonics, Ferroelectrics,and Frequency Control, 49, The unmodified and modified Hadamard variances are called the Picinbono and pulsar variances. P. Lesage and C. Audoin, 97, Characterization of frequency stability: uncertainty due to the finite number of measurements, IEEE Transactions on Instrumentation and Measurement, IM, 576. P. Lesage and C. Audoin, 974, Correction to Characterization of frequency stability: uncertainty due to the finite number of measurements, IEEE Transactions on Instrumentation and Measurement, IM, 0. P. Lesage and C. Audoin, 975, Comments on Characterization of frequency stability: uncertainty due to the number of measurements, IEEE Transactions on Instrumentation and Measurement, IM4, 86 P. Lesage and C. Audoin, 976, Correction to characterization of frequency stability: uncertainty due to the number of measurements, IEEE Transactions on Instrumentation and Measurement., IM5, 70. P. Lesage and C. Audoin, 977, Estimation of the twosample variance with limited number of data, st Annual Symposium on Frequency Control, 8. D. Howe, D. Allan, and J. Barnes, 98, Properties of signal sources and measurement methods, 5th Annual Symposium on Frequency Control, 47. Reprinted in [8], TN460; see also [9], 5. In the flicker PM expression on p. 7, the argument of exp should be raised to the / power; in the flicker FM expression, N should be squared. These corrections have been made in [lo], 4. D. Sullivan, D. Allan, D. Howe, and F. Walls (editors), 990, Characterization of Clocks and Oscillators, NIST Technical Note 7, U. S. Department of Commerce, National Institute of Standards and Technology. Online at go./ timefreq/gener al/tnl7/ S. Stein, 985, Frequency and timetheir measurement and characterization, in E. Gerber and A. Ballato (editors), Precision Frequency Control, Academic Press, New York. Reprinted in [8], TN60. [lo] IEEE Std 9999, IEEE Standard Definitions of Physical Quantities for Fundamental Frequency and Time MetrologyRandom Instabilities. [ll] K. Yoshimura, 978, Characterization of frequency stability: uncertainty due to the autocorrelation of the frequency fluctuations, IEEE Transactions on Instrumentation and Measurement, IM7, 978, 7. [] K. Yoshimura, 989, Degrees of freedom of the estimate of the twosample variance in the continuous sampling method, IEEE Transactions on Instrumentation and Measurement, 8, [] C. Greenhall, 98, A structure function representation theorem with applications to frequency stability estimation, IEEE Transactions on Instrumentation and Measurement, IM,
10 [4] C. Greenhall, 99, Recipes for degrees of freedom of frequency stability estimators, IEEE Transactions on Instrumentation and Measurement, 40, [5] T. Walter, 994, Characterizing frequency stability: a continuous powerlaw model with discrete sampling, IEEE Transactions on Instrumentation and Measurement, 4, [I6 C. Greenhall, 997, The thirddifference approach to modified Allan variance, IEEE Transactions on Instrumentation and Measurement, 46, [7] Stable, Software for frequency stability analysis, Hamilton Technical Services, 95 Woodbury St., S. Hamilton, MA 098, [8] Timing Solutions Corporation, 4775 Walnut St. Suite lb, Boulder, CO , www. timing.com [9] D. Howe, R. Beard, C. Greenhall, F. Vernotte, and W. Riley, 00, A total estimator of the Hadamard function used for GPS operations, Proceedings of the nd Annual Precise Time and Time Interval (PTTI) Meeting, Nov. 000, Reston, VA, [0] D. Howe, 000, The total deviation approach to longterm characterization of frequency stability, IEEE Transactions on Ultrasonics, Ferroelectriqand Frequency Control, 47, 00. [] D. Howe and T. Peppler, 00, Very longterm frequency stability: estimation using a specialpurpose statistic, Proceedings of the 00 IEEE International Frequency Control Symposium, 8. [] Upon request by the second author will send an electronic copy of C source code that implements this algorithm, along with a small Microsoft Windows (r) program that computes the for various cases. 0
11 Overlapping ADEV EDF for N=lOO I 0 00 Averaging Factor, m Figure. Edf vs. averaging factor with powerlaw noise type as a parameter: Allan variance, overlapped estimator, 00 phase samples. 00 Overlapping EDF for N=00 and W FM Noise Averaging Factor, m Figure. Edf vs. averaging factor for three stability variances: overlapped estimator, white FM, 00 phase samples.
UNCERTAINTY OF STABILITY VARIANCES BASED ON FINITE DIFFERENCES
UNCERTAINTY OF STABILITY VARIANCES BASED ON FINITE DIFFERENCES C. A. Greenhall W. J. Riley Jet Propulsion Laboratory Symmetricom, Inc. California Institute of Technology Abstract This work gives an algorithm
More informationA SIMPLE ALGORITHM FOR APPROXIMATING CONFIDENCE ON THE MODIFIED ALLAN VARIANCE AND THE TIME VARIANCE*
A SIMPLE ALGORITHM FOR APPROXIMATING CONFIDENCE ON THE MODIFIED ALLAN VARIANCE AND THE TIME VARIANCE* Marc A. Weiss Time & Frequency Div., 847.5 National Institute of Standards and Technology 325 Broadway,
More informationDiscrete Simulation of Power Law Noise
Discrete Simulation of Power Law Noise Neil Ashby 1,2 1 University of Colorado, Boulder, CO 80309-0390 USA 2 National Institute of Standards and Technology, Boulder, CO 80305 USA ashby@boulder.nist.gov
More informationA SIMPLE ALGORITHM FOR APPROXIMATING CONFIDENCE ON THE MODIFIED ALLAN VARIANCE AND THE TIME VARIANCE*
A SIMPLE ALGORITHM FOR APPROXIMATING CONFIDENCE ON THE MODIFIED ALLAN VARIANCE AND THE TIME VARIANCE* Marc A. Weiss Time & Frequency Div., 847.5 National Institute of Standards and Technology 325 Broadway,
More informationDEGREES OF FREEDOM AND THREE-CORNERED HATS
33rd Annual Precise Time and Time Interval (PTTI) Meeting DEGREES OF FREEDOM AND THREE-CORNERED HATS Christopher R. Ekstrom and Paul A. Koppang U.S. Naval Observatory, Washington, D.C. Abstract The three-cornered
More informationA KALMAN FILTER CLOCK ALGORITHM FOR USE IN THE PRESENCE OF FLICKER FREQUENCY MODULATION NOISE
A KALMAN FILTER CLOCK ALGORITHM FOR USE IN THE PRESENCE OF FLICKER FREQUENCY MODULATION NOISE J. A. Davis, C. A. Greenhall *, and P. W. Stacey National Physical Laboratory Queens Road, Teddington, Middlesex,
More informationVARIANCES BASED ON DATA WITH DEAD TIME BETWEEN THE MEASUREMENTS *
VARIANCES BASED ON DATA WITH DEAD TIME BETWEEN THE MEASUREMENTS * James A. Barnes Austron Incorporated Boulder, Colorado 80301 David W. Allan Time and Frequency Division National Bureau of Standards Boulder,
More informationUncertainties of drift coefficients and extrapolation errors: Application to clock error prediction
Vernotte et al - E 1319 - Final Version - Metrologia - February 1, 21 1 Uncertainties of drift coefficients and extrapolation errors: Application to clock error prediction F Vernotte 1, J Delporte 2, M
More informationMTIE AND TDEV ANALYSIS OF UNEVENLY SPACED TIME SERIES DATA AND ITS APPLICATION TO TELECOMMUNICATIONS SYNCHRONIZATION MEASUREMENTS
MTIE AND TDEV ANALYSIS OF UNEVENLY SPACED TIME SERIES DATA AND ITS APPLICATION TO TELECOMMUNICATIONS SYNCHRONIZATION MEASUREMENTS Mingfu Li, Hsin-Min Peng, and Chia-Shu Liao Telecommunication Labs., Chunghwa
More informationTHE LONG-TERM STABILITY OF THE U.S. NAVAL OBSERVATORY S MASERS
THE LONG-TERM STABILITY OF THE U.S. NAVAL OBSERVATORY S MASERS Demetrios Matsakis, Paul Koppang Time Service Department U.S. Naval Observatory Washington, DC, USA and R. Michael Garvey Symmetricom, Inc.
More informationUse of the Autocorrelation Function for Frequency Stability Analysis
Use of the Autocorrelation Function for Frequency Stability Analysis W.J. Riley, Hamilton Technical Services Introduction This paper describes the use of the autocorrelation function (ACF) as a complement
More informationCONFIDENCE ON THE THREE-POINT ESTIMATOR OF FREQUENCY DRIFT*
CONFIDENCE ON THE THREE-POINT ESTIMATOR OF FREQUENCY DRIFT* Marc A. Weiss and Christine Hackman Time and Frequency Division National Institute of Standards and Technology 325 Broadway Boulder, CO 80303
More informationINVESTIGATING THE BIASES IN ALLAN AND HADAMARD VARIANCES AS MEASURES OF M TH ORDER RANDOM STABILITY
43 rd Annual Precise Time and Time Interval (PTTI) Systems and Applications Meeting INVESTIGATING THE BIASES IN ALLAN AND HADAMARD VARIANCES AS MEASURES OF M TH ORDER RANDOM STABILITY Victor R. Reinhardt
More informationA KALMAN FILTER FOR ATOMIC CLOCKS AND TIMESCALES
33rdAnnual Precise Time and Time Interval (PTTI) Meeting A KALMAN FILTER FOR ATOMIC CLOCKS AND TIMESCALES Lee A. Breakiron U.S. Naval Observatory Washington, DC 2392542,USA Abstract In a test of whether
More informationSTUDIES OF NPL S CLOCK ENSEMBLE ALGORITHM
43 rd Annual Precise Time and Time Interval (PTTI) Systems and Applications Meeting STUDIES OF NPL S CLOCK ENSEMBLE ALGORITHM Setnam L. Shemar, John A. Davis, and Peter B. Whibberley National Physical
More informationTHE LONG-TERM STABILITY OF THE U.S. NAVAL OBSERVATORY S MASERS
THE LONG-TERM STABILITY OF THE U.S. NAVAL OBSERVATORY S MASERS Demetrios Matsakis, Paul Koppang Time Service Department U.S. Naval Observatory Washington, DC, USA and R. Michael Garvey Symmetricom, Inc.
More informationNoise analysis and identification in MEMS sensors, Allan, Time, Hadamard, Overlapping, Modified, Total variance
DT0064 Design tip Noise analysis and identification in MEMS sensors, Allan, Time, Hadamard, Overlapping, Modified, Total variance By Andrea Vitali Main components LSM6DS3H LSM6DSM/LSM6DSL STEVAL-STLKT01V1
More informationA KALMAN FILTER FOR ATOMIC CLOCKS AND TIMESCALES
33rdAnnual Precise Time and Time Interval (PTTI) Meeting A KALMAN FILTER FOR ATOMIC CLOCKS AND TIMESCALES Lee A. Breakiron U.S. Naval Observatory Washington, DC 2392-542,USA Abstract In a test of whether
More informationLee A. Breakiron U.S. Naval Observatory Washington, DC ABSTRACT
A COMPARATIVE STUDY OF CLOCK RATE AND DRIFT ESTIMATION Lee A. Breakiron U.S. Naval Observatory Washington, DC 20392 ABSTRACT Five different methods of drift determination and four different methods of
More informationALLAN VARIANCE ESTIMATED BY PHASE NOISE MEASUREMENTS
ALLAN VARIANCE ESTIMATED BY PHASE NOISE MEASUREMENTS P. C. Chang, H. M. Peng, and S. Y. Lin National Standard Time & Frequency Lab., TL, Taiwan 1, Lane 551, Min-Tsu Road, Sec. 5, Yang-Mei, Taoyuan, Taiwan
More informationKALMAN PLUS WEIGHTS: A TIME SCALE ALGORITHM*
33rdAnnual Precise Time and Time nterval (PlTl) Meeting KALMAN PLUS WEGHTS: A TME SCALE ALGORTHM* Charles A. Greenhall Jet Propulsion Laboratory, California nstitute of Technology 48 Oak Grove Drive, MS
More informationcovariance function, 174 probability structure of; Yule-Walker equations, 174 Moving average process, fluctuations, 5-6, 175 probability structure of
Index* The Statistical Analysis of Time Series by T. W. Anderson Copyright 1971 John Wiley & Sons, Inc. Aliasing, 387-388 Autoregressive {continued) Amplitude, 4, 94 case of first-order, 174 Associated
More informationHOW TO DEAL WITH FFT SAMPLING INFLUENCES ON ADEV CALCULATIONS
HOW TO DEAL WITH FFT SAMPLING INFLUENCES ON ADEV CALCULATIONS Po-Ceng Cang National Standard Time & Frequency Lab., TL, Taiwan 1, Lane 551, Min-Tsu Road, Sec. 5, Yang-Mei, Taoyuan, Taiwan 36 Tel: 886 3
More informationEssentials of expressing measurement uncertainty
Essentials of expressing measurement uncertainty This is a brief summary of the method of evaluating and expressing uncertainty in measurement adopted widely by U.S. industry, companies in other countries,
More informationModeling and Estimation of Stationary and Non-stationary Noises of Rubidium Atomic Clock
ISSN : 48-96, Vol. 4, Issue 7( Version ), July 014, pp.44-49 RESEARCH ARTICLE OPEN ACCESS Modeling and Estimation of Stationary and Non-stationary Noises of Rubidium Atomic Clock 1 Deepak Mishra, Alak
More informationAN ADAPTIVE ALGORITHM TO EVALUATE CLOCK PERFORMANCE IN REAL TIME*
AN ADAPTIVE ALGORITHM TO EVALUATE CLOCK PERFORMANCE IN REAL TIME* Dr. James A. Barnes Austron Boulder, Co. Abstract Kalman filters and ARIMA models provide optimum control and evaluation techniques (in
More information6.435, System Identification
System Identification 6.435 SET 3 Nonparametric Identification Munther A. Dahleh 1 Nonparametric Methods for System ID Time domain methods Impulse response Step response Correlation analysis / time Frequency
More informationStudy of a Non-Intrusrve blectron
Study of a Non-Intrusrve blectron b n i F -910s Diagnosti &m~i Radius LA-OP--77-3/37 Thomas J. T. Kwan, Barbara G. DeVolder, John C. Goldstein, and Charles M. Snell Los Alamos National Laboratory, Los
More informationKALMAN PLUS WEIGHTS: A TIME SCALE ALGORITHM*
33rdAnnual Precise Time and Time nterval (PlTl) Meeting KALMAN PLUS WEGHTS: A TME SCALE ALGORTHM* Charles A. Greenhall Jet Propulsion Laboratory, California nstitute of Technology 48 Oak Grove Drive, MS
More informationStability Variances: A filter Approach.
Stability Variances: A filter Approach. Alaa Madissi, François Vernotte and Emeric De Clercq Abstract We analyze the Allan Variance estimator as the combination of Discrete-Time linear filters. We apply
More informationNON-LINEAR PARAMETER ESTIMATION USING VOLTERRA AND WIENER THEORIES
Journal of Sound and Vibration (1999) 221(5), 85 821 Article No. jsvi.1998.1984, available online at http://www.idealibrary.com on NON-LINEAR PARAMETER ESTIMATION USING VOLTERRA AND WIENER THEORIES Department
More informationLEAST-SQUARES ANALYSIS OF TWO-WAY SATELLITE TIME AND FREQUENCY TRANSFER MEASUREMENTS
33rdAnnual Precise Time and Time Interval (Z'lTl) Meeting LEAST-SQUARES ANALYSIS OF TWO-WAY SATELLITE TIME AND FREQUENCY TRANSFER MEASUREMENTS J. A. Davis, P. M. Harris, S. L. Shemar, and M. G. Cox National
More informationTHE RANGE COVERED BY A CLOCK ERROR IN THE CASE OF WHITE FM
! 90th Annual Precise Time and Time Interval (PTTI) Meeting THE RANGE COVERED BY A CLOCK ERROR IN THE CASE OF WHITE FM Patrizia Tavella* and Daniela MeoO JCIstituto Elettrotecnico Nazionale G. Ferraris
More informationWhy Correlation Matters in Cost Estimating
Why Correlation Matters in Cost Estimating Stephen A. Book The Aerospace Corporation P.O. Box 92957 Los Angeles, CA 90009-29597 (310) 336-8655 stephen.a.book@aero.org 32nd Annual DoD Cost Analysis Symposium
More informationModeling Laser and e-beam Generated Plasma-Plume Experiments Using LASNEX
UCRL-ID-136726 Modeling Laser and e-beam Generated Plasma-Plume Experiments Using LASNEX D.D.-M. Ho December 1,1999 US. Department of Energy Approved for public release; further dissemination unlimited
More information0 0'0 2S ~~ Employment category
Analyze Phase 331 60000 50000 40000 30000 20000 10000 O~----,------.------,------,,------,------.------,----- N = 227 136 27 41 32 5 ' V~ 00 0' 00 00 i-.~ fl' ~G ~~ ~O~ ()0 -S 0 -S ~~ 0 ~~ 0 ~G d> ~0~
More informationA note on a Marčenko-Pastur type theorem for time series. Jianfeng. Yao
A note on a Marčenko-Pastur type theorem for time series Jianfeng Yao Workshop on High-dimensional statistics The University of Hong Kong, October 2011 Overview 1 High-dimensional data and the sample covariance
More informationSTAT 520: Forecasting and Time Series. David B. Hitchcock University of South Carolina Department of Statistics
David B. University of South Carolina Department of Statistics What are Time Series Data? Time series data are collected sequentially over time. Some common examples include: 1. Meteorological data (temperatures,
More informationWavelet Methods for Time Series Analysis. Part IV: Wavelet-Based Decorrelation of Time Series
Wavelet Methods for Time Series Analysis Part IV: Wavelet-Based Decorrelation of Time Series DWT well-suited for decorrelating certain time series, including ones generated from a fractionally differenced
More informationMultiresolution Models of Time Series
Multiresolution Models of Time Series Andrea Tamoni (Bocconi University ) 2011 Tamoni Multiresolution Models of Time Series 1/ 16 General Framework Time-scale decomposition General Framework Begin with
More informationFog Monitor 100 (FM 100) Extinction Module. Operator Manual
Particle Analysis and Display System (PADS): Fog Monitor 100 (FM 100) Extinction Module Operator Manual DOC-0217 Rev A-1 PADS 2.7.3, FM 100 Extinction Module 2.7.0 5710 Flatiron Parkway, Unit B Boulder,
More informationSTUDIES OF THE UNBIASED FIR FILTER FOR THE TIME ERROR MODEL IN APPLICATIONS TO GPS-BASED TIMEKEEPING 1
STUDIES OF THE UNBIASED FIR FILTER FOR THE TIME ERROR MODEL IN APPLICATIONS TO GPS-BASED TIMEKEEPING 1 Yu. Shmaliy and O. Ibarra-Manzano Guanajuato University, FIMEE, Mexico E-mail: shmaliy@salamanca.ugto.mx
More informationADSP ADSP ADSP ADSP. Advanced Digital Signal Processing (18-792) Spring Fall Semester, Department of Electrical and Computer Engineering
Advanced Digital Signal rocessing (18-792) Spring Fall Semester, 201 2012 Department of Electrical and Computer Engineering ROBLEM SET 8 Issued: 10/26/18 Due: 11/2/18 Note: This problem set is due Friday,
More informationTowards an Optimal Noise Versus Resolution Trade-off in Wind Scatterometry
Towards an Optimal Noise Versus Resolution Trade-off in Wind Scatterometry Brent Williams Jet Propulsion Lab, California Institute of Technology IOWVST Meeting Utrecht Netherlands June 12, 2012 Copyright
More informationThe Modified Allan Variance as Time-Domain Analysis Tool for Estimating the Hurst Parameter of Long-Range Dependent Traffic
The Modified Allan Variance as Time-Domain Analysis Tool for Estimating the urst Parameter of Long-Range Dependent Traffic Stefano Bregni, Senior Member, IEEE, Luca Primerano Politecnico di Milano, Dept.
More information5 Generation and Dissemination of Time and Frequency Standard
5 Generation and Dissemination of Time and Frequency Standard 5-1 Algorithm of Ensemble Atomic Time HANADO Yuko, IMAE Michito, AIDA Masanori, HOSOKAWA Mizuhiko, ITO Hiroyuki, NAKAGAWA Fumimaru, and SHIMIZU
More informationPAijpam.eu ADAPTIVE K-S TESTS FOR WHITE NOISE IN THE FREQUENCY DOMAIN Hossein Arsham University of Baltimore Baltimore, MD 21201, USA
International Journal of Pure and Applied Mathematics Volume 82 No. 4 2013, 521-529 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu doi: http://dx.doi.org/10.12732/ijpam.v82i4.2
More informationTHE FUTURE MODEL OF TA (TL)
THE FUTURE MODEL OF TA (TL) Shinn-Yan Lin, Hsin-Ming Peng, Weng-Hung Tseng, Hung-Tien Lin, and Chia-Shu Liao National Standard Time and Frequency Lab., TL, Chunghwa Telecom Co. Ltd. No. 12 Lane 551, Min-Tsu
More informationUPGRADED CALIBRATIONS OF THE THOMSON SYSTEM AT DIII D
GA A23440 UPGRADED CALIBRATIONS OF THE THOMSON SYSTEM AT DIII D by B. BRAY, C. HSIEH, T.N. CARLSTROM, and C.C. MAKARIOU AUGUST 2000 DISCLAIMER This report was prepared as an account of work sponsored by
More informationSignals and Spectra - Review
Signals and Spectra - Review SIGNALS DETERMINISTIC No uncertainty w.r.t. the value of a signal at any time Modeled by mathematical epressions RANDOM some degree of uncertainty before the signal occurs
More information2. SPECTRAL ANALYSIS APPLIED TO STOCHASTIC PROCESSES
2. SPECTRAL ANALYSIS APPLIED TO STOCHASTIC PROCESSES 2.0 THEOREM OF WIENER- KHINTCHINE An important technique in the study of deterministic signals consists in using harmonic functions to gain the spectral
More informationAnalysis of Geodetic Time Series Using Allan Variances
Universität Stuttgart Geodätisches Institut Analysis of Geodetic Time Series Using Allan Variances Studienarbeit im Studiengang Geodäsie und Geoinformatik an der Universität Stuttgart Thomas Friederichs
More informationARIMA Modelling and Forecasting
ARIMA Modelling and Forecasting Economic time series often appear nonstationary, because of trends, seasonal patterns, cycles, etc. However, the differences may appear stationary. Δx t x t x t 1 (first
More informationStochastic Processes: I. consider bowl of worms model for oscilloscope experiment:
Stochastic Processes: I consider bowl of worms model for oscilloscope experiment: SAPAscope 2.0 / 0 1 RESET SAPA2e 22, 23 II 1 stochastic process is: Stochastic Processes: II informally: bowl + drawing
More informationTHE GPS TOOLKIT: OPEN SOURCE CLOCK TOOLS
THE GPS TOOLKIT: OPEN SOURCE CLOCK TOOLS Timothy J. H. Craddock, Richard John Broderick, Colin P. Petersen, and Alex Hu Applied Research Laboratories: The University of Texas at Austin Space and Geophysics
More informationSYSM 6303: Quantitative Introduction to Risk and Uncertainty in Business Lecture 4: Fitting Data to Distributions
SYSM 6303: Quantitative Introduction to Risk and Uncertainty in Business Lecture 4: Fitting Data to Distributions M. Vidyasagar Cecil & Ida Green Chair The University of Texas at Dallas Email: M.Vidyasagar@utdallas.edu
More informationLQSAlamos National Laboratory
DEVELOPMENT AND APPLICATION OF ADVANCED ONEPOINT TURBULENCE MODELS Charles G. Speziale January 1996 Technical Report No. AM96002 DISCLAIMER This report was prepared as an account of work sponsored by an
More informationRCA Analysis of the Polar Codes and the use of Feedback to aid Polarization at Short Blocklengths
RCA Analysis of the Polar Codes and the use of Feedback to aid Polarization at Short Blocklengths Kasra Vakilinia, Dariush Divsalar*, and Richard D. Wesel Department of Electrical Engineering, University
More informationIf we want to analyze experimental or simulated data we might encounter the following tasks:
Chapter 1 Introduction If we want to analyze experimental or simulated data we might encounter the following tasks: Characterization of the source of the signal and diagnosis Studying dependencies Prediction
More informationProbability and Statistics for Final Year Engineering Students
Probability and Statistics for Final Year Engineering Students By Yoni Nazarathy, Last Updated: May 24, 2011. Lecture 6p: Spectral Density, Passing Random Processes through LTI Systems, Filtering Terms
More informationThe Identification of ARIMA Models
APPENDIX 4 The Identification of ARIMA Models As we have established in a previous lecture, there is a one-to-one correspondence between the parameters of an ARMA(p, q) model, including the variance of
More informationTIME SERIES AND FORECASTING. Luca Gambetti UAB, Barcelona GSE Master in Macroeconomic Policy and Financial Markets
TIME SERIES AND FORECASTING Luca Gambetti UAB, Barcelona GSE 2014-2015 Master in Macroeconomic Policy and Financial Markets 1 Contacts Prof.: Luca Gambetti Office: B3-1130 Edifici B Office hours: email:
More informationRESPONSE SPECTRUM METHOD FOR ESTIMATION OF PEAK FLOOR ACCELERATION DEMAND
RESPONSE SPECTRUM METHOD FOR ESTIMATION OF PEAK FLOOR ACCELERATION DEMAND Shahram Taghavi 1 and Eduardo Miranda 2 1 Senior catastrophe risk modeler, Risk Management Solutions, CA, USA 2 Associate Professor,
More informationAlamos National Laboratory is operated by the University of California for the United States Department of Energy under contract W-7405-ENG-36
LA- UR "6 Ea Alamos National Laboratory is operated by the University of California for the United States Department of Energy under contract W-7405-ENG-36 TITLE: AUTHOR(S): SUBMITTED TO: SPATIAL AND ANGULAR
More informationTime series models in the Frequency domain. The power spectrum, Spectral analysis
ime series models in the Frequency domain he power spectrum, Spectral analysis Relationship between the periodogram and the autocorrelations = + = ( ) ( ˆ α ˆ ) β I Yt cos t + Yt sin t t= t= ( ( ) ) cosλ
More informationRESONANCE INTERFERENCE AND ABSOLUTE CROSS SECTIONS IN NEAR-THRESHOLD ELECTRON-IMPACT EXCITATION OF MULTICHARGED IONS
RESONANCE INTERFERENCE AND ABSOLUTE CROSS SECTIONS IN NEAR-THRESHOLD ELECTRON-IMPACT EXCITATION OF MULTICHARGED IONS M. E. BANNISTER Physics Division, Oak Ridge National Labomtory, Oak Ridge, TN 37831-6372,USA
More informationThe Kalman Filter. Data Assimilation & Inverse Problems from Weather Forecasting to Neuroscience. Sarah Dance
The Kalman Filter Data Assimilation & Inverse Problems from Weather Forecasting to Neuroscience Sarah Dance School of Mathematical and Physical Sciences, University of Reading s.l.dance@reading.ac.uk July
More informationGamma process model for time-dependent structural reliability analysis
Gamma process model for time-dependent structural reliability analysis M.D. Pandey Department of Civil Engineering, University of Waterloo, Waterloo, Ontario, Canada J.M. van Noortwijk HKV Consultants,
More informationEconometrics II - EXAM Answer each question in separate sheets in three hours
Econometrics II - EXAM Answer each question in separate sheets in three hours. Let u and u be jointly Gaussian and independent of z in all the equations. a Investigate the identification of the following
More informationSystem Identification Approach Applied to Drift Estimation.
Instrumentation and Measurement Technology Conference - IMTC 007 Warsaw, Poland, May 1-3, 007 System Identification Approach Applied to Drift Estimation. Frans Verbeyst 1,, Rik Pintelon 1, Yves Rolain
More informationFinancial Econometrics and Quantitative Risk Managenent Return Properties
Financial Econometrics and Quantitative Risk Managenent Return Properties Eric Zivot Updated: April 1, 2013 Lecture Outline Course introduction Return definitions Empirical properties of returns Reading
More informationESTIMATION OF TIME-DOMAIN FREQUENCY STABILITY BASED ON PHASE NOISE MEASUREMENT
ESTIMATION OF TIME-DOMAIN FREQUENCY STABILITY BASED ON PHASE NOISE MEASUREMENT P. C. Cang, H. M. Peng, and S. Y. Lin National Standard Time & Frequenc Laborator, TL, Taiwan, Lane 55, Min-Tsu Road, Sec.
More informationPlutonium 239 Equivalency Calculations
LLNL-TR-490356 Plutonium 239 Equivalency Calculations J. Wen July 7, 2011 Disclaimer This document was prepared as an account of work sponsored by an agency of the United States government. Neither the
More informationOn Exponential Decay and the Riemann Hypothesis
On Exponential Decay and the Riemann Hypothesis JEFFREY N. COOK ABSTRACT. A Riemann operator is constructed in which sequential elements are removed from a decaying set by means of prime factorization,
More informationSOLVING POWER AND OPTIMAL POWER FLOW PROBLEMS IN THE PRESENCE OF UNCERTAINTY BY AFFINE ARITHMETIC
SOLVING POWER AND OPTIMAL POWER FLOW PROBLEMS IN THE PRESENCE OF UNCERTAINTY BY AFFINE ARITHMETIC Alfredo Vaccaro RTSI 2015 - September 16-18, 2015, Torino, Italy RESEARCH MOTIVATIONS Power Flow (PF) and
More informationI I FINAL, 01 Jun 8.4 to 31 May TITLE AND SUBTITLE 5 * _- N, '. ', -;
R AD-A237 850 E........ I N 11111IIIII U 1 1I!til II II... 1. AGENCY USE ONLY Leave 'VanK) I2. REPORT DATE 3 REPORT TYPE AND " - - I I FINAL, 01 Jun 8.4 to 31 May 88 4. TITLE AND SUBTITLE 5 * _- N, '.
More informationStatistical Methods in Particle Physics
Statistical Methods in Particle Physics Lecture 11 January 7, 2013 Silvia Masciocchi, GSI Darmstadt s.masciocchi@gsi.de Winter Semester 2012 / 13 Outline How to communicate the statistical uncertainty
More informationProf. Dr. Roland Füss Lecture Series in Applied Econometrics Summer Term Introduction to Time Series Analysis
Introduction to Time Series Analysis 1 Contents: I. Basics of Time Series Analysis... 4 I.1 Stationarity... 5 I.2 Autocorrelation Function... 9 I.3 Partial Autocorrelation Function (PACF)... 14 I.4 Transformation
More informationV003 How Reliable Is Statistical Wavelet Estimation?
V003 How Reliable Is Statistical Wavelet Estimation? J.A. Edgar* (BG Group) & M. van der Baan (University of Alberta) SUMMARY Well logs are often used for the estimation of seismic wavelets. The phase
More informationA source model for ISDN packet data traffic *
1 A source model for ISDN packet data traffic * Kavitha Chandra and Charles Thompson Center for Advanced Computation University of Massachusetts Lowell, Lowell MA 01854 * Proceedings of the 28th Annual
More informationTHE OPTICAL SCIENCES COMPANY P.O. Box 446, Placentia, California Phone: (714) Pulsed Laser Backscatter. Glenn A. Tyler.
THE OPTICAL SCIENCES COMPANY P.O. Box 446, Placentia, California 92670 Phone: (714) 524-3622 Report No. TR-400 (v \P Pulsed Laser Backscatter Glenn A. Tyler and David L. Fried November 1980 WBTmiüTIÖlf
More informationFinancial Econometrics
Financial Econometrics Nonlinear time series analysis Gerald P. Dwyer Trinity College, Dublin January 2016 Outline 1 Nonlinearity Does nonlinearity matter? Nonlinear models Tests for nonlinearity Forecasting
More informationX t = a t + r t, (7.1)
Chapter 7 State Space Models 71 Introduction State Space models, developed over the past 10 20 years, are alternative models for time series They include both the ARIMA models of Chapters 3 6 and the Classical
More informationECE295, Data Assimila0on and Inverse Problems, Spring 2015
ECE295, Data Assimila0on and Inverse Problems, Spring 2015 1 April, Intro; Linear discrete Inverse problems (Aster Ch 1 and 2) Slides 8 April, SVD (Aster ch 2 and 3) Slides 15 April, RegularizaFon (ch
More informationStochastic Processes. A stochastic process is a function of two variables:
Stochastic Processes Stochastic: from Greek stochastikos, proceeding by guesswork, literally, skillful in aiming. A stochastic process is simply a collection of random variables labelled by some parameter:
More informationIndex I-1. in one variable, solution set of, 474 solving by factoring, 473 cubic function definition, 394 graphs of, 394 x-intercepts on, 474
Index A Absolute value explanation of, 40, 81 82 of slope of lines, 453 addition applications involving, 43 associative law for, 506 508, 570 commutative law for, 238, 505 509, 570 English phrases for,
More informationA statistic for describing pulsar and clock stabilities
Astron. Astrophys. 326, 924 928 (1997) ASTRONOMY AND ASTROPHYSICS A statistic for describing pulsar and clock stabilities Demetrios N. Matsakis 1, J. H. Taylor 2, and T. Marshall Eubanks 1 1 U.S. Naval
More informationOn Exponential Decay and the Riemann Hypothesis
On Exponential Decay and the Riemann Hypothesis JEFFREY N. COOK ABSTRACT. A Riemann operator is constructed in which sequential elements are removed from a decaying set by means of prime factorization,
More informationNAVAL POSTGRADUATE SCHOOL Monterey, California THESIS MAXIMIZING THE STABILITY OF AN ENSEMBLE OF CLOCKS. Juan J. Ruiz
NAVAL POSTGRADUATE SCHOOL Monterey, California THESIS MAXIMIZING THE STABILITY OF AN ENSEMBLE OF CLOCKS by Juan J. Ruiz Thesis Advisor: Second Reader: September 3 Alan Washburn Paul Sánchez Approved for
More informationLecture 14: Strain Examples. GEOS 655 Tectonic Geodesy Jeff Freymueller
Lecture 14: Strain Examples GEOS 655 Tectonic Geodesy Jeff Freymueller A Worked Example Consider this case of pure shear deformation, and two vectors dx 1 and dx 2. How do they rotate? We ll look at vector
More informationApproximating the Conway-Maxwell-Poisson normalizing constant
Filomat 30:4 016, 953 960 DOI 10.98/FIL1604953S Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Approximating the Conway-Maxwell-Poisson
More informationSpectral representations and ergodic theorems for stationary stochastic processes
AMS 263 Stochastic Processes (Fall 2005) Instructor: Athanasios Kottas Spectral representations and ergodic theorems for stationary stochastic processes Stationary stochastic processes Theory and methods
More informationStatistics. Lent Term 2015 Prof. Mark Thomson. 2: The Gaussian Limit
Statistics Lent Term 2015 Prof. Mark Thomson Lecture 2 : The Gaussian Limit Prof. M.A. Thomson Lent Term 2015 29 Lecture Lecture Lecture Lecture 1: Back to basics Introduction, Probability distribution
More informationComputational Data Analysis!
12.714 Computational Data Analysis! Alan Chave (alan@whoi.edu)! Thomas Herring (tah@mit.edu),! http://geoweb.mit.edu/~tah/12.714! Introduction to Spectral Analysis! Topics Today! Aspects of Time series
More informationA NEW REALIZATION OF TERRESTRIAL TIME
CCTF/04-17 A NEW REALIZATION OF TERRESTRIAL TIME G. Petit BIPM, Pavillon de Breteuil, 92312 Sèvres Cedex France- e-mail: gpetit@bipm.org ABSTRACT Terrestrial Time TT is a time coordinate in a geocentric
More informationUNIFORMLY MOST POWERFUL CYCLIC PERMUTATION INVARIANT DETECTION FOR DISCRETE-TIME SIGNALS
UNIFORMLY MOST POWERFUL CYCLIC PERMUTATION INVARIANT DETECTION FOR DISCRETE-TIME SIGNALS F. C. Nicolls and G. de Jager Department of Electrical Engineering, University of Cape Town Rondebosch 77, South
More informationHAPPY BIRTHDAY CHARLES
HAPPY BIRTHDAY CHARLES MY TALK IS TITLED: Charles Stein's Research Involving Fixed Sample Optimality, Apart from Multivariate Normal Minimax Shrinkage [aka: Everything Else that Charles Wrote] Lawrence
More informationUncertainty due to Finite Resolution Measurements
Uncertainty due to Finite Resolution Measurements S.D. Phillips, B. Tolman, T.W. Estler National Institute of Standards and Technology Gaithersburg, MD 899 Steven.Phillips@NIST.gov Abstract We investigate
More informationFast and accurate estimation and adjustment of a local
86 ieee transactions on ultrasonics ferroelectrics and frequency control vol. 53 no. 5 may 006 An Unbiased FIR Filter for TIE Model of a Local Clock in Applications to GPS-Based Timekeeping Yuriy S. Shmaliy
More information(4) How do you develop an optimal signal detection technique from the knowledge of
Signal and Noise in Global Warming Detection Final Report Gerald R. North, Principal Investigator Climate System Research Program, Texas A&M University, Texas 1. Program Objectives The specific objectives
More information