MCM-based Uncertainty Evaluations practical aspects and critical issues

Similar documents
Investigation of Uncertainty Sources in the Determination of Beta Emitting Tritium in the UAL

Traceability and uncertainty for phase measurements

Our focus will be on linear systems. A system is linear if it obeys the principle of superposition and homogenity, i.e.

BAYESIAN CURVE FITTING USING PIECEWISE POLYNOMIALS. Dariusz Biskup

Training Course Textbook Emma Woolliams Andreas Hueni Javier Gorroño. Intermediate Uncertainty Analysis

Physics 123. Exam #1. October 11, 2006

Watermarking 2D-Vector Data for Geographical Information Systems

1 Definition of Rademacher Complexity

Guidelines on the Estimation of Uncertainty in Hardness Measurements

Linear Regression Model

AE/ME 339. K. M. Isaac. 8/31/2004 topic4: Implicit method, Stability, ADI method. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept.

ASYMMETRIC TRAFFIC ASSIGNMENT WITH FLOW RESPONSIVE SIGNAL CONTROL IN AN URBAN NETWORK

Some basic statistics and curve fitting techniques

Statistics for Managers Using Microsoft Excel/SPSS Chapter 13 The Simple Linear Regression Model and Correlation

Lecture 3. Camera Models 2 & Camera Calibration. Professor Silvio Savarese Computational Vision and Geometry Lab. 13- Jan- 15.

04 - Treaps. Dr. Alexander Souza

UNCERTAINTY OF AIRCRAFT NOISE MEASUREMENTS: EVALUATION FOR AN AIRCRAFT NOISE MONITORING NETWORK

Several generation methods of multinomial distributed random number Tian Lei 1, a,linxihe 1,b,Zhigang Zhang 1,c

STRATEGIES FOR BUILDING AN AC-DC TRANSFER SCALE

Statistics Chapter 4

Denote the function derivatives f(x) in given points. x a b. Using relationships (1.2), polynomials (1.1) are written in the form

Optimal Marketing Strategies for a Customer Data Intermediary. Technical Appendix

Fall 2012 Analysis of Experimental Measurements B. Eisenstein/rev. S. Errede. ) with a symmetric Pcovariance matrix of the y( x ) measurements V

EQUATION CHAPTER 1 SECTION 1STRAIN IN A CONTINUOUS MEDIUM

Failure of Assumptions

Final Report on the International Comparison of Luminous Responsivity CCPR-K3.b

System in Weibull Distribution

Linear Multiple Regression Model of High Performance Liquid Chromatography

Description of the Force Method Procedure. Indeterminate Analysis Force Method 1. Force Method con t. Force Method con t

PARETO OPTIMAL ROBUST FEEDBACK LINEARIZATION CONTROL OF A NONLINEAR SYSTEM WITH PARAMETRIC UNCERTAINTIES

Uncertainty in measurements of power and energy on power networks

The Bellman Equation

Applied Mathematics Letters

AN ANALYSIS OF A FRACTAL KINETICS CURVE OF SAVAGEAU

CHAPTER 10 ROTATIONAL MOTION

total If no external forces act, the total linear momentum of the system is conserved. This occurs in collisions and explosions.

Least Squares Fitting of Data

Computational and Statistical Learning theory Assignment 4

EURAMET.M.D-S2 Final Report Final report

Lab 2e Thermal System Response and Effective Heat Transfer Coefficient

NON-SYNCHRONOUS TILTING PAD BEARING CHARACTERISTICS

ESS 265 Spring Quarter 2005 Time Series Analysis: Error Analysis

MATHEMATICAL SIMULATION OF THE MAGNETIC FIELD OCCURRED BY THE ARMATURE REACTION OF THE SYNCHRONOUS MACHINE

On the number of regions in an m-dimensional space cut by n hyperplanes

Variance of Estimates in Dynamic Data Reconciliation

Demonstrate Teaching FEM, Particle Dynamics, and Simple MD Codes Using Matlab & Mathematica

AGC Introduction

IMPROVEMENT OF THE REALIZATION OF FORCES BETWEEN 2 MN AND 5 MN AT PTB THE NEW 5 MN FORCE STANDARD MACHINE

LECTURE :FACTOR ANALYSIS

Primary Velocity Distribution in Open Channels with Different Vegetation Layout - Experiment and Numerical Simulation -

Linear classification models: Perceptron. CS534-Machine learning

Systematic Error Illustration of Bias. Sources of Systematic Errors. Effects of Systematic Errors 9/23/2009. Instrument Errors Method Errors Personal

Output-Feedback Model Predictive Control with Constraints Improve Robustness Using Past Information

Prof. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model

The Impact of the Earth s Movement through the Space on Measuring the Velocity of Light

EXACT TRAVELLING WAVE SOLUTIONS FOR THREE NONLINEAR EVOLUTION EQUATIONS BY A BERNOULLI SUB-ODE METHOD

8/25/17. Data Modeling. Data Modeling. Data Modeling. Patrice Koehl Department of Biological Sciences National University of Singapore

Statistics for Economics & Business

Correlation and Regression. Correlation 9.1. Correlation. Chapter 9

A Differential Evaluation Markov Chain Monte Carlo algorithm for Bayesian Model Updating M. Sherri a, I. Boulkaibet b, T. Marwala b, M. I.

8.3 Divide & Conquer for tridiagonal A

A DIMENSIONALITY REDUCTION ALGORITHM OF HYPER SPECTRAL IMAGE BASED ON FRACT ANALYSIS

Analytical Chemistry Calibration Curve Handout

Positron Lifetime Spectroscopy

Research Note NONLINEAR ANALYSIS OF SEMI-RIGID FRAMES WITH RIGID END SECTIONS * S. SEREN AKAVCI **

Numerical Simulation of Melting Process in Single Screw Extruder with Vibration Force Field

Integral Transforms and Dual Integral Equations to Solve Heat Equation with Mixed Conditions

1 Matrix representations of canonical matrices

Microrelief measurements for white-light interferometer with adaptive algorithm interferogram processing

Solutions to selected problems from homework 1.

APPENDIX 2 FITTING A STRAIGHT LINE TO OBSERVATIONS

Finite Vector Space Representations Ross Bannister Data Assimilation Research Centre, Reading, UK Last updated: 2nd August 2003

EMISSION MEASUREMENTS IN DUAL FUELED INTERNAL COMBUSTION ENGINE TESTS

Excess Error, Approximation Error, and Estimation Error

VERIFICATION OF FE MODELS FOR MODEL UPDATING

PHYS 1443 Section 002 Lecture #20

CLOSED-FORM CHARACTERIZATION OF THE CHANNEL CAPACITY OF MULTI-BRANCH MAXIMAL RATIO COMBINING OVER CORRELATED NAKAGAMI FADING CHANNELS

Quantum Particle Motion in Physical Space

Message modification, neutral bits and boomerangs

THE CURRENT BALANCE Physics 258/259

The Analysis and the Performance Simulation of the Capacity of Bit-interleaved Coded Modulation System

IET (2015) ISSN

Invertible Linear Transforms Implemented by Integer Mapping. (Published in Chinese, Science in China, Series E, April 2000, V30, N2, pp )

Spectral Recomposition in Stratigraphic Interpretation*

Formal solvers of the RT equation

The Electromechanical Behavior of a Micro-Beam Driven by Traveling Electrostatic Force

UNCERTAINTY ON SIGNAL PARAMETER ESTIMATION IN FREQUENCY DOMAIN

XII.3 The EM (Expectation-Maximization) Algorithm

Basic Business Statistics, 10/e

PHYS 342L NOTES ON ANALYZING DATA. Spring Semester 2002

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Rectilinear motion. Lecture 2: Kinematics of Particles. External motion is known, find force. External forces are known, find motion

Sampling Theory MODULE V LECTURE - 17 RATIO AND PRODUCT METHODS OF ESTIMATION

/ n ) are compared. The logic is: if the two

Global Sensitivity. Tuesday 20 th February, 2018

ˆ (0.10 m) E ( N m /C ) 36 ˆj ( j C m)

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

Statistical analysis of Accelerated life testing under Weibull distribution based on fuzzy theory

Physics 3A: Linear Momentum. Physics 3A: Linear Momentum. Physics 3A: Linear Momentum. Physics 3A: Linear Momentum

1.3 Hence, calculate a formula for the force required to break the bond (i.e. the maximum value of F)

Transcription:

C-based Uncertanty Evalatons practcal aspects and crtcal sses H. Hatjea, B. van Dorp,. orel and P.H.J. Schellekens Endhoven Unversty of Technology

Contents Introdcton Standard ncertanty bdget de wthot dervatves onte-carlo ncertanty evalaton achne eleent effects Exaples of vrtal easreent systes

Introdcton: What s a vrtal nstrent? Conventonal easrng nstrent Vrtal easrng nstrent d =.577 d =.577 =.75 A vrtal nstrent calclates ts ncertanty tself

The standard ncertanty bdget accordng to EAL and GU: =,..); = ) Qantty Vale of qantty Standard ncertanty ) Senstvty coeffcent c = Contrbton to standard-ncertanty n c ) = ) c ) ) c ) : : : : : N N N ) c N N ) total: ) + ) = +.. N

Uncertanty bdget wthot explct partal dervatves = Qantty.. ± ),..); =...) Vale of qantty Standard ncertanty Contrbton to standard-ncertanty n ) ) ) = +r ),,,..., N ) -,,... N ) ) =, +r ),..., N ) -,... N ) : : : : : N N N ) N =,,..., N +r N )) -,... N ) total ) = + +.. N Vary each paraeter separately wth ts ncertanty? ± = + r r = ± rando nber)

onte Carlo Uncertanty Bdget ) ± + = = r r,,...),....); Vary all paraeter wth ts ncertanty at a te? ± = + r r = ± rando nber) Or r = rando nber wth <r>=, s r = Slaton nber k Estaton of standard-ncertanty ),..., )) ),..., N N N N r r + + = ),..., )) ),..., N N N N r r + + = : K ),..., )) ),..., N N KN N K K r r + + = K.. K + + + =

Characterstcs of onte Carlo Uncertanty Bdget: dfferent dstrbtons can be cobned hgher order ters are taken nto accont nknown systeatc errors can be slated correlatons can be slated can be sed for very coplcated probles not too any slatons needed N=5? s s =%)

Inflence of flterng effects Calclate Paraeters e.g Ra =.86 µ Calclate standard ncertanty for any paraeter P = P nonal P + wth P nonal e.g Ra =.86 ±. µ Calclate Paraeters e.g Ra =.76 µ Exaple how to slate/evalate probe effect

Inflence of teperatre effects o L = f L, α, T, δa, δt ) = L + α δ T + T C) δα) L s length at C L s easred length T s average teperatre of object and scale a s average expanson coeffcent δα s dfference n expanson coeffcent δt s teperatre dfference between object and scale ethod to randose: calclate L back), randose paraeters, calclate new L

-Densonal Geoetrcal devatons Scale devatons are never rando: take correlatons nto accont ) ) ) ) j j j j x x j j x x x AC x x AC x x x x + = + = AC s atocorrelaton fncton dx x x L AC L + = δ δ δ δ ) ) )

-Densonal Geoetrcal devatons analytcal calclaton for easred length 3 devaton n µ - 3 4 5 6 7 8 9 poston n easred devatons along a easreent lne of a C 4 ncertanty n µ 3 L δ AC δ ) = x) x + δ L δ ) dx 3 4 5 6 7 8 9 easred length n Standard ncertanty n a length easreent along ths axs

-Densonal Geoetrcal devatons: How to randose the: Atocorrelaton fncton st be constant 4 evdent possbltes : org nal sgnal flpped left rght error n croeter.5 -.5 error n croeter.5 -.5-3 - 3 reversed sgn flpped left rght and reversed sgn error n croeter.5 -.5 error n croeter.5 -.5-3 poston n - 3 poston n

-Densonal Geoetrcal devatons: How to randose the: Alternatve for perodc sgnal: Take forer transfor randose phases transfor back exaple: rondness dagra alternatve: rotate dagra over rando angle

-Densonal Geoetrcal devatons: How to randose the: Proble for non-contnos sgnal Dscontnty gves hgher haroncs

lt-densonal Geoetrcal devatons Lnear, straghtness and rotatonal errors can all be slated lkewse antanng atocorrelaton) Correlaton between axes: sqareness: slated separately

Exaple: Vrtal roghness tester No onte-carlo: jst ncertanty bdget Qantty vale standard ncertanty contrbton to ncertanty n Ra 8,3 n contrbton to ncertanty n RS 4 µ x-axs ) %,4 µ calbraton z-axs ) %,8 n calbraton λ c,5 %,5 n,4 µ λ s,5 µ %,43 n,3 µ F N 5%, n,3 n rads µ 5%, n, µ step sze,5 µ %, µ nose 3 n RS, n, n total s),6 n,6 µ

Exaple: Srface plate easreent onte-carlo: sperpose rando ncertanty on easreents, re-calclate paraeter exaple: ncertanty n easred heghts

Exaple: Vrtal Fzea nterferoeter Randoze reference srface: flp, rotate, negate

Vrtal C Calbrated/easred devatons n -D s T y y T y x y T x y x T x x

Vrtal C easreent strategy: easre real probed objects/dsplaceents: ths ncldes software correcton easre sall and larger wavelengths For each lnearty/straghtness: 4 µ wth µ steps.6 wth 4 µ steps 3 wth µ steps.5 wth steps)

Vrtal C Straghtness easreent: cerac straghtedge

Vrtal C Straghtness easreent: cerac straghtedge.5 error xty, easred, stepsze..5 error xty, easred, stepsze.6 error n croeter error n croeter -.5.5..5. -.5 5 5 5 3 poston n poston n Short wavelengths Long wavelengths

Vrtal C Lnearty easreent: laser-nterferoeter step gage

Vrtal C Lnearty easreent: laser-nterferoeter step gage xtx error, easred, stepsze. xtx error, easred, stepsze.6 error n croeter.5 -.5 error n croeter.5 -.5 -.5..5. - 5 5 5 3 poston n poston n Short wavelengths Long wavelengths

Vrtal C Note short-wavelength yty-devatons nflence on slghtly salgned straghtedge n easreent of xty.8 straghtness devaton n µ.6.4. -. -.4 -.6 -.8 5 5 5 3 x posto n n

Vrtal C Synthesse rando sgnal fro 3 easreent ranges.8 exaple for xtx.6 error n croeter.4. -...4.6.8.. poston n

Vrtal C Slate and calclaton: synthetc data generaton wth own progra calclaton of geoetrcal paraeters easy ones: own progra coplcated ones: Zess Uess va backdoor, not easy) 8 9 x axs 8 4 6 8 y axs

Vrtal C Reslts for rng gages, as a fncton of nber of easred pont.5 daeter.5 daeter 5.45 odelled easred.45 odelled easred standard ncertanty, croeters.4.35.3.5..5. standard ncertanty, croeters.4.35.3.5..5..5.5 4 6 6 3 4 6 6 3 nber of easreent ponts nber of easreent ponts

Vrtal C Unsolved probles: traceablty: what to do wth ncertanty n calbrated objects repeatablty: enters n all calbratons, bt only once n real easreents dstncton between rando, nknown systeatc, known systeatc vanshes n ths approach alas of short-wavelength errors n long-wavelenght calbratons what to do wth ncertanty-bas n for easreents