PROFICIENCY TESTING ACTIVITIES OF FREQUENCY CALIBRATION LABORATORIES IN TAIWAN, 2009

Similar documents
TRACEABILITY SYSTEM OF ROCKWELL HARDNESS C SCALE IN JAPAN

Measurement uncertainty of the sound absorption

Standard BAL-001-0a Real Power Balancing Control Performance

Provläsningsexemplar / Preview TECHNICAL REPORT INTERNATIONAL SPECIAL COMMITTEE ON RADIO INTERFERENCE

Statistical Analysis on Uncertainty for Autocorrelated Measurements and its Applications to Key Comparisons

ANALYSIS OF EXPERIMENTAL ERRORS

Probability, Expectation Value and Uncertainty

Statistics 511 Additional Materials

Accuracy assessment methods and challenges

577. Estimation of surface roughness using high frequency vibrations

Activity 3: Length Measurements with the Four-Sided Meter Stick

x a x a Lecture 2 Series (See Chapter 1 in Boas)

PH 425 Quantum Measurement and Spin Winter SPINS Lab 1

Accuracy of prediction methods for the improvement of impact sound pressure levels using floor coverings

The AMSU Observation Bias Correction and Its Application Retrieval Scheme, and Typhoon Analysis

WHAT IS THE PROBABILITY FUNCTION FOR LARGE TSUNAMI WAVES? ABSTRACT

Chapter 2 Descriptive Statistics

A statistical method to determine sample size to estimate characteristic value of soil parameters

Standard BAL Real Power Balancing Control Performance

FACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING. Lectures

Analysis of Experimental Measurements

CHAPTER 10 INFINITE SEQUENCES AND SERIES

Hypothesis Testing. Evaluation of Performance of Learned h. Issues. Trade-off Between Bias and Variance

The improvement of the volume ratio measurement method in static expansion vacuum system

OBJECTIVES. Chapter 1 INTRODUCTION TO INSTRUMENTATION FUNCTION AND ADVANTAGES INTRODUCTION. At the end of this chapter, students should be able to:

Expected mean in an environmental noise measurement and its related uncertainty

Statistical Fundamentals and Control Charts

SECTION 1.5 : SUMMATION NOTATION + WORK WITH SEQUENCES

CEE 522 Autumn Uncertainty Concepts for Geotechnical Engineering

Free Space Optical Wireless Communications under Turbulence Channel Effect

The target reliability and design working life

PSYCHOLOGICAL RESEARCH (PYC 304-C) Lecture 9

Measures of Spread: Variance and Standard Deviation

Lecture 12: Hypothesis Testing

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +

Kinetics of Complex Reactions

Error & Uncertainty. Error. More on errors. Uncertainty. Page # The error is the difference between a TRUE value, x, and a MEASURED value, x i :

CHAPTER 8 FUNDAMENTAL SAMPLING DISTRIBUTIONS AND DATA DESCRIPTIONS. 8.1 Random Sampling. 8.2 Some Important Statistics

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

Sample Size Estimation in the Proportional Hazards Model for K-sample or Regression Settings Scott S. Emerson, M.D., Ph.D.

1 Review of Probability & Statistics

a. For each block, draw a free body diagram. Identify the source of each force in each free body diagram.

Problem of Technical Losses in Natural Gas Distribution and Transportation Systems of Georgia

6 Sample Size Calculations

A THRESHOLD DENOISING METHOD BASED ON EMD

A GERMAN CALIBRATION GUIDELINE FOR TORQUE WRENCH CALIBRATION DEVICES REVISED

Determining the sample size necessary to pass the tentative final monograph pre-operative skin preparation study requirements

CENTRIFUGAL PUMP SPECIFIC SPEED PRIMER AND THE AFFINITY LAWS Jacques Chaurette p. eng., Fluide Design Inc. November 2004

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

THE SYSTEMATIC AND THE RANDOM. ERRORS - DUE TO ELEMENT TOLERANCES OF ELECTRICAL NETWORKS

Understanding Samples

Analysis of Experimental Data

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,

SNAP Centre Workshop. Basic Algebraic Manipulation

Run-length & Entropy Coding. Redundancy Removal. Sampling. Quantization. Perform inverse operations at the receiver EEE

Sample Size Determination (Two or More Samples)

General IxJ Contingency Tables

Chapter 6 Sampling Distributions

On a Smarandache problem concerning the prime gaps

Drive Technology \ Drive Automation \ System Integration \ Services. Data Sheet. Functional Safety Safety Characteristics for Brake BE..

Oblivious Transfer using Elliptic Curves

AP Calculus Chapter 9: Infinite Series

Simple Random Sampling!

MEASURES OF DISPERSION (VARIABILITY)

Optimally Sparse SVMs

The standard deviation of the mean

Chapter 12 Correlation

STAT 350 Handout 19 Sampling Distribution, Central Limit Theorem (6.6)

Chapter 8: Estimating with Confidence

Rank Modulation with Multiplicity

Control Charts for Mean for Non-Normally Correlated Data

5. A formulae page and two tables are provided at the end of Part A of the examination PART A

Module 5 EMBEDDED WAVELET CODING. Version 2 ECE IIT, Kharagpur

Invariability of Remainder Based Reversible Watermarking

Economics 250 Assignment 1 Suggested Answers. 1. We have the following data set on the lengths (in minutes) of a sample of long-distance phone calls

µ and π p i.e. Point Estimation x And, more generally, the population proportion is approximately equal to a sample proportion

1 Adiabatic and diabatic representations

Big Picture. 5. Data, Estimates, and Models: quantifying the accuracy of estimates.

The axial dispersion model for tubular reactors at steady state can be described by the following equations: dc dz R n cn = 0 (1) (2) 1 d 2 c.

The Pendulum. Purpose

Drive Technology \ Drive Automation \ System Integration \ Services. Data Sheet. Functional Safety Safety Characteristics for BE..

Estimation of a population proportion March 23,

Random Variables, Sampling and Estimation

If, for instance, we were required to test whether the population mean μ could be equal to a certain value μ

Research Article Health Monitoring for a Structure Using Its Nonstationary Vibration

Lecture 5: Parametric Hypothesis Testing: Comparing Means. GENOME 560, Spring 2016 Doug Fowler, GS

The Local Harmonious Chromatic Problem

Infinite Sequences and Series

Section 11.8: Power Series

THE KALMAN FILTER RAUL ROJAS

FIR Filters. Lecture #7 Chapter 5. BME 310 Biomedical Computing - J.Schesser

ECE 8527: Introduction to Machine Learning and Pattern Recognition Midterm # 1. Vaishali Amin Fall, 2015

Data Analysis and Statistical Methods Statistics 651

Some Properties of the Exact and Score Methods for Binomial Proportion and Sample Size Calculation

Statistical Hypothesis Testing. STAT 536: Genetic Statistics. Statistical Hypothesis Testing - Terminology. Hardy-Weinberg Disequilibrium

Math 140 Introductory Statistics

ROLL CUTTING PROBLEMS UNDER STOCHASTIC DEMAND

Principle Of Superposition

Deterministic Model of Multipath Fading for Circular and Parabolic Reflector Patterns

Transcription:

PROFICIENCY TESTING ACTIVITIES OF FREQUENCY CALIBRATION LABORATORIES IN TAIWAN, 009 Huag-Tie Li, Po-Cheg Chag, Jia-Lu Wag, ad Chia-Shu Liao Natioal Stadard Time & Frequecy Lab., TL, Chughwa Telecom Co., Ltd., Taiwa 1, Lae 551, Mi-Tsu Road, Sec. 5, Yag-Mei, Taoyua, Taiwa 36 Tel: +886-3-44-4066, Fax: +886-3-44-5474, E-mail: liht@cht.com.tw Abstract I order to meet the requiremets of ISO 1705 ad the demad of TAF (Taiwa Accreditatio Foudatio) for calibratio iter-laboratory comparisos, the Natioal Time ad Frequecy Stadard Laboratory (TL) [1], playig the role of coordiatig lab, has periodically orgaized the proficiecy testig activities to provide opportuities for capability comparisos of domestic frequecy calibratio laboratories. The latest two activities were performed i 006 [] ad this year (009). The Device uder Test (DUT), together with TL s calibratio system, was trasferred to each participatig laboratory accordig to a predetermied schedule. The DUT was measured by the participatig lab s ad TL s calibratio systems simultaeously. The measured results from both systems were the aalyzed ad compared. There were 1 ad 15 participatig labs joied the activities performed i 006 ad 009, respectively. All of them were TAF-accredited laboratories. As usual, the E value was used to evaluate a lab s capability of calibratig equipmet withi its accredited measuremet ucertaity. I these two activities, the absolute E values of all labs were smaller tha 1, which meas their calibratio ability was completely qualified. I this paper, the related details of these two activity ad the result comparisos betwee activities i 006 ad 009 are illustrated. INTRODUCTION The goal of performig the proficiecy testig activity is to provide a chace for the domestic accredited laboratories to compare their calibratio capability with oe aother. Whe performig o-site evaluatio, a assessmet team is orgaized to examie the techical competece of the labs ad their compliace with the requiremets of ISO/IEC 1705 Geeral requiremets for the competece of testig ad calibratio laboratories. The proficiecy testig results are the importat iformatio for the techical assessmet team to evaluate the ability of the laboratory. I proficiecy testig activities [3,4], the ability of each idividual laboratory is supposed to be determied ad able to achieve its accredited level of measuremet ucertaity. I the case that the performace of the Device uder Test (DUT) has good repeatability ad stability, its referece value 313

will be assiged by the coordiatig lab ad the compared with measuremet results of the participatig labs, as show i Fig 1. Whe the referece value is located withi a lab s ucertaity rage (e.g., lab 1,, ad 3 i Fig.1), we say that this lab works withi its capability ad its accredited level of measuremet ucertaity is suitable for the DUT beig calibrated. Ref Value 1 3 =Ucertaity Rage N N =Results of LAB N Fig. 1. Measuremet results ad ucertaity rages of participatig labs. However, the characteristic of the DUT i the time/frequecy area (such as a frequecy oscillator) is ostatioary, so the referece value may chage gradually. Therefore, together with the DUT, a calibratio system of the coordiatig lab is usually trasferred to each participat laboratory to get the relative referece value simultaeously with the measuremet. I other words, the DUT is measured by the participat s ad the coordiatig lab s calibratio systems at the same time. EVALUATION OF RESULTS I the proficiecy test activities amog calibratio laboratories, the E value is adopted to idicate how well labs are withi their particular measuremet ucertaity, takig accout of the measuremet ucertaity of the referece value. E stads for Error ormalized ad is defied as: Xi X ref E (1) U ( X ) U ( X ) 95 i 95 ref where X i is the participat laboratory s result X ref is the coordiatig laboratory s result U 95 (X i ) is the participatig laboratory s reported ucertaity (95%) U 95 (X ref ) is the coordiatig laboratory s reported ucertaity (95%) E 1 idicates that the result ad the referece value are i agreemet E 1 idicates that the result is differet from the referece value. Note that the calculatio of E value does ot ited to idicate which lab s result is closest to the referece value, as high-level calibratio labs may have their E values similar to that of some other labs with both larger measuremet errors ad ucertaities. Besides, oe should be aware that the referece value itself has a measuremet ucertaity. The coordiatig lab should have the capability to give a better measuremet ucertaity tha the participatig lab s; otherwise, it will be difficult to evaluate each lab s performace. Cosequetly, the E 1 limit really oly represets the cutoff, below which it is likely that the result is acceptable ad above which it is 314

ulikely that the result is acceptable. So whe cosiderig ay result with E greater tha 1, all factors should be evaluated to see if there is a systematic bias that is cosistetly positive or cosistetly egative which causes the problem. PROFICIENCY TESTING ACTIVITY TL, playig the role of coordiatig lab, has periodically orgaized the proficiecy testig activities every 3 years sice 003. The latest two activities were performed i 006 ad 009, with HP866A [5] ad HP3350A frequecy sythesizers as the DUT, respectively. Both HP866A ad HP3350A ca offer 5/10 MHz output sigals, so ay laboratory with a calibratio capability for either 5 MHz or 10 MHz ca joi these activities. The agig rate of HP866A is smaller tha 5.0 10-10 /day, which is lower tha that of HP3350A, which is about 5.5 10-9 /day. Sice the properties of a oscillator s sigal are o-statioary, TL s calibratio system, cosistig of a 5071A cesium clock, a SR60 time-iterval couter, ad a recordig computer alog with the DUT, were set to each participatig lab sequetially. The DUT s referece value ca the be obtaied durig the measuremet process i each lab. The above metioed 5071A cesium clock is TL s portable frequecy referece, which is traceable to the atioal frequecy stadard. To make sure that the HP5071A ca satisfy the requiremet of traceability, the performace of this Cs clock was measured i TL before ad after every measuremet trip to the participatig labs. The block diagram of the related system for the proficiecy testig activities is show i Fig.. 5/10 MHz Participat Laboratory s Calibratio System DUT (HP866a / HP3350a) Power Splitter 5/10 MHz Recordig Computer TIC (SR60) Cs Clock (5071a) TL s Calibratio System Fig.. The system arragemet i proficiecy testig activities. The umber of participatig labs that joied the proficiecy testig activities was 1 i 006 ad icreased to 15 i 009, but the same procedure was followed. Before the measuremet trips to the participatig labs i each activity, a opeig meetig was held i advace to discuss all details about the related activities. A cosesus referrig to the schedule for trasferrig the DUT ad TL s calibratio system, the measuremet method, the ucertaity evaluatio, the result expressio, etc. was reached. Sice 006, TL has recommeded that the participatig labs iclude the best measuremet capability of the idividual lab, the DUT s agig rate ad temperature effects, etc. i their ucertaity budgets. 315

Note that, i geeral calibratio services, a ucertaity budget usually icluded the DUT s agig to assure the validity of the measuremet result durig the period that the DUT retured to its lab till the ext time the DUT was set for calibratio. However, it is t ecessary to cout the agig for a log period (for example, 1 year) i proficiecy testig activities, because the DUT oly serves as a itermediate for proficiecy testig comparisos ad the measuremet time i each participatig labs is ot loger tha 1 hour. Therefore, we aouced DUT s agig rate for 1 day to each participatig lab for their evaluatio of measuremet ucertaity. While fiishig the measuremet, each participatig lab was required to sed its test report, raw measuremet data, ad ucertaity budget to TL withi 3 days. After the measuremet trips to the participatig labs, collectig the measuremet data ad evaluatig the capability of each laboratory, a meetig was held for fial discussio ad providig the chace for faceto-face commuicatio with the participat labs. RESULTS ANALYSIS AND CONCLUSION The proficiecy testig results of 1 participatig labs i 006 ad 15 participatig labs i this year are show i Tables 1 ad, respectively. X i, X ref, U 95 (X i ) ad U 95 (X ref ) stad for the correspodig measuremet results (frequecy accuracy) ad reported ucertaities from participatig labs ad TL. Usig these four parameters, the E value of each participatig lab could be obtaied. Cosiderig a driftig frequecy stadard may ifluece the measuremet result, it is reasoable for participat labs to iclude the latest calibrated accuracy of their idividual frequecy stadards i the X i measuremets, or separately specify both of them i each participat lab s test report. Results of the latter are adopted i Table 1 ad the metioed calibrated accuracy ΔX i ca be used for the E correctio. The corrected E value could be obtaied usig equatio as below: ( X i X i ) X ref E. () U ( X ) U ( X ) 95 i 95 ref Besides, adjustmet of X ref is egligible because the calibrated accuracy of TL s portable frequecy referece ca reach 1.5 10-13, about 4 orders lower tha the magitude of X ref itself. For most of the participats i Table 1, the corrected E are smaller tha or equal to the origial oes. This shows that the E correctio is meaigful, which gives more reasoable results for iter-laboratory comparisos. Equatio () ca also be cosidered as the corrected measuremet differece betwee a participatig lab ad TL, divided by their combied ucertaity. I geeral, U 95 (X i ) is the domiat term i the combied ucertaity, so the corrected 316 E is a good idex to show how well a participatig lab is withi its accredited capability. All the corrected E of 1 participatig labs i 006 ad 15 participatig labs this year are smaller tha 1, which idicates that the calibratio capabilities of all the participats are qualified i these two activities. Moreover, oe may foud that most of the corrected E values i Table are smaller tha those i Table 1. This is caused by the performaces of differet DUTs adopted i these two activities. The agig rate ad temperature coefficiet of a HP3350A (DUT of 009) are larger tha those of a HP866A (DUT

of 006), so the combied ucertaities are commoly larger ad the correspodet E values become smaller. Geerally speakig, a DUT with smaller agig ad temperature coefficiet would be more suitable for the proficiecy testig applicatio, ad the capability of the participatig labs ca be more properly evaluated. REFERENCES [1] H. T. Li, C. S. Liao, S. Y. Li, W. H. Tseg, C. C. Li, P. C. Chag, D. Wag, F. D. Chu, ad H. Shyu, Highlights of the 008 activities of the Natioal Time ad Frequecy Stadard Laboratory of the Telecommuicatio Laboratories, CHT CO. LTD., Taiwa, i Proceedigs of the Asia-Pacific Workshop o Time ad Frequecy (ATF), 30 October-1 November 008, Jakarta, Idoesia. [] P. C. Chag, H. T. Li, ad C. S. Liao, Proficiecy testig activities of frequecy calibratio laboratories i Taiwa, 006, i Proceedigs of the Asia-Pacific Workshop o Time ad Frequecy (ATF), December 006, New Delhi, Idia. [3] NATA Proficiecy Testig Traiig Course. [4] Guide to the Expressio of Ucertaity i Measuremet ( ISO), 1995. [5] HP 866A Sythesized Sigal Geerator Operatig ad Services Maual, Sec. 1 (Hewlett-Packard), pp. 4-9. Table 1. Proficiecy testig results of 1 participatig labs. Report No. X i X ref ΔX i U 95 (X i ) U 95 (X ref ) E E (corrected) TL-PT-01-3.0343E-09-4.30E-09 -.6E-10 1.19E-09 3.76E-10 1.01 0.81 TL-PT-0-3.8000E-09-3.70E-09 6.4E-1 4.10E-10 4.68E-10 0.16 0.15 TL-PT-03.5507E-09 3.00E-09 3.E-10 7.3E-10 3.76E-10 0.55 0.16 TL-PT-04-1.0000E-08-4.30E-09 3.3E-11 5.8E-09 6.83E-10 0.97 0.97 TL-PT-05-3.0800E-09 -.90E-09-1.9E-13 6.18E-09 3.08E-10 0.03 0.03 TL-PT-06-4.550E-09-5.50E-09-1.0E-09 3.16E-09 3.76E-10 0.31 0.01 TL-PT-07-4.4500E-08-1.40E-09 9.3E-09 8.96E-08 4.68E-10 0.48 0.38 TL-PT-08-4.6000E-09-4.90E-09 1.3E-10 9.00E-10 3.76E-10 0.31 0.44 TL-PT-09 -.3400E-09-1.60E-09 1.4E-10 9.13E-10 3.38E-10 0.76 0.6 TL-PT-10-3.6330E-09-3.30E-09 1.5E-1 5.05E-10 3.76E-10 0.53 0.53 TL-PT-11-5.6900E-09-4.74E-09-7.6E-11 9.0E-08 3.76E-10 0.01 0.01 TL-PT-1-3.0000E-09 -.10E-09-3.E-10 3.30E-07 3.76E-10 0.01 0.01 317

Table. Proficiecy testig results of 15 participatig labs. Report No. X i X ref ΔX i U 95 (X i ) U 95 (X ref ) E (corrected) TL-98FMPT-01-4.918564E-07-4.9345E-07-1.70E-13 7.571E-08 6.93E-08 0.005 TL-98FMPT-0-3.177987E-07-3.17600E-07-1.40E-11 1.500E-07 6.93E-08 0.001 TL-98FMPT-03-3.841748E-07-4.3000E-07 3.0E-08 1.500E-06 6.93E-08 0.005 TL-98FMPT-04-4.11597E-07-4.11000E-07 -.40E-10 4.900E-07 6.93E-08 0.001 TL-98FMPT-05-5.097840E-07-5.10000E-07 1.45E-10 4.00E-07 6.93E-08 0.001 TL-98FMPT-06-3.6069E-07-3.0000E-07 3.70E-10 1.400E-07 6.93E-08 0.041 TL-98FMPT-07-4.194414E-07-4.1000E-07 8.10E-11 3.300E-07 6.93E-08 0.004 TL-98FMPT-08 1.11796E-08 1.11000E-08-1.50E-13 4.300E-08 6.93E-08 0.001 TL-98FMPT-09-5.55941E-07-5.5000E-07-5.00E-10 1.400E-07 6.93E-08 0.001 TL-98FMPT-10-4.964631E-07-4.96185E-07 1.50E-10 1.00E-06 6.93E-08 0.001 TL-98FMPT-11-3.686349E-07-3.68000E-07-1.0E-09 1.390E-07 6.93E-08 0.004 TL-98FMPT-1-3.43065E-07-3.493E-07-1.00E-1 1.390E-07 6.93E-08 0.00 TL-98FMPT-13-3.316300E-07-3.37000E-07 6.50E-09.00E-07 6.93E-08 0.005 TL-98FMPT-14-3.99086E-07-3.30000E-07 9.90E-11 6.80E-08 6.93E-08 0.001 TL-98FMPT-15-3.6563E-07-3.94400E-07 6.00E-10 6.80E-08 6.93E-08 0.9 318