Content
Workshop DYNAMIC2016
9th International Workshop on Analysis of Dynamic Measurements
Nov. 910 2016, PTB Berlin, Germany
Dynamic measurements can be found in many areas of metrology and industry, such as, for instance, in applications where mechanical quantities, electrical pulses or temperature curves are measured.
A quantity is called dynamic when its value at one time instant depends on its values at previous time instants.
That is, in contrast to static measurements where a single value or a (small) set of values is measured, dynamic measurements consider continuous functions of time. Since the analysis of dynamic measurements requires different approaches than the analysis of static measurements this part of metrology is often called "Dynamic Metrology".
The mathematical modeling of dynamic measurements typically utilizes methodologies and concepts from digital signal processing. In the language of metrology a signal denotes a dynamic quantity, and a system a measurement device whose input and/or output are signals. The output signal of a system is thus the indication value of the measurement device for a corresponding input signal.
The typical scenario in a dynamic measurement is a timedependent input signal, a linear timeinvariant (LTI) measurement system and a corresponding timedependent ouptut signal. The linearity of the system refers to the superposition of input signals, and the timeinvariance ensures that the system itself does not change over time. Mathematically, the relation between input signal and output signal is in this case given by a convolution
$$ x(t) = \int y(\tau) h(\taut) d\tau , $$
with $y(t)$ the input signal (measurand), $h(t)$ the system's impulse response, and $x(t)$ the system's output signal. The measurand is then the input signal $x(t)$ of the measuring system, and its estimation requires deconvolution.
Statistical modeling of static quantities is realized by univariate or multivariate random variables. The evaluation of uncertainty is then based on probability density functions. In a mathematical sense, the extension of this concept to dynamic measurements requires the utilization of stochastic processes as a model for the uncertain knowledge of the continuous functions $x(t)$ and $y(t)$. That is, the dynamic quantity $x(t)$ is considered as a continuous function in time and uncertain knowledge about its values is expressed by a stochastic process $X_t$ with continuous trajectories. The working group has developed a corresponding framework for the evaluation of uncertainties which is consistent with the framework of the GUM for static measurements.
However, instead of the continuous function $x(t)$, typically its discretization is considered in the application. Therefore, the continuous output signal of the measurement system is stored on a computer as a vector, whereby $\mathbf{x}=(x_1,\ldots,x_N)$ with $x_k=x(t_k)$. Statistical modeling of the uncertain knowledge of the values of such a vector and its propagation can then, in principle, be carried out by using the concepts and framework of GUM Supplement 2. However, this is possible only for measurements on closed time intervals. Thus, it is not possible to treat running (ontime) measurements without making further assumptions.
A characteristic property of a dynamic measurement is that the output signal is not proportional to the input signal. This is due to socalled dynamic effects caused by the imperfect dynamic behavior of the measuring system. For instance, accelerometers typically show a resonance behavior, and hence, for a measured acceleration with a certain frequency content the output signal of the accelerometer shows a significant "ringing".
The aim in the analysis of dynamic measurements is the compensation for timedependent errors, such as, ringing, phase deviations and other effects. In contrast to static measurements, this cannot be accomplished by scaling and shifting the output signal, since it only takes into account the static characteristics of the measuring system. Instead, for linear timeinvariant systems (LTI) a socalled deconvolution has to be carried out. This allows the compensation of dynamic effects and, thereby, in principle, the reconstruction of the values of the dynamic measurand.
The analysis of dynamic measurements in the case of linear timeinvariant (LTI) systems is typically carried out by application of a suitable digital filter. The design of such a filter is based on the available knowledge about the measurement system and aims at compensating its dynamic effects. As illustrated in the example in Fig. 2, the frequency response of the compensation filter is the reciprocal of the system's frequency response up to a certain frequency. Thus, the prerequisite for the design of a compensation filter is a dynamic calibration of the measurement device in a suitable frequency range.
The literal meaning of "dynamic" relates solely to timedependent quantities. However, from a mathematical perspective the definition of a dynamic measurement can be extended to other independent quantities than time. This includes quantities whose value depends on frequency, spatial coordinates, wavelength, etc. The extended definition is reasonable, because the mathematical treatment of such measurements does not depend on the physical interpretation of the independent quantity.
A quantity is a dynamic quantity if its value depends on another, independent quantity. A measurement is dynamic when at least one of the involved quantities is dynamic.
The extended definition of a dynamic measurement contains a broad spectrum of metrological applications. Typical examples are measurements of mechanical quantities as, for example, force, torque and pressure. Further examples are oscilloscope measurements for the characterization of high speed electronics, hydrophone measurements for the characterization of medical ultrasound devices or the spectral characterization of radiation sources. The extended definition also includes, for example, the spectral characterization of luminous sources, spectral color measurements or cameraaided temperature measurements. The applications range from single sensor measurements up to large sensor networks.
In many areas of application, dynamic properties and dynamic errors have been neglected so far. Instead, ruleofthumb correction methods were applied and larger uncertainties assigned to the measurement result. However, in recent years the demand for more precise dynamic measurements with smaller uncertainties has increased steadily.
Research
The PTB Working Group 8.42 has been carrying out research in the field of dynamic measurements for more than 10 years and its scientific work covers almost all areas of dynamic metrology. Publications of the working group contain, for instance, methods for the statistical analysis of dynamic calibration, design of digital deconvolution filters for estimating the value of the measurand, GUMcompliant evaluation of dynamic measurement uncertainty, and efficient implementation of GUM Monte Carlo for the application of digital filters. Moreover, the working group works closely together with industrial partners to promote developments in dynamic metrology.
The mathematical and statistical treatment of dynamic measurements requires different approaches and tools than the analysis of static measurements. Thus, the areas of research in the working group are determined by the challenges in dynamic metrology.
Propagation of uncertainties
In general, the evaluation and propagation of measurement uncertainties for discretized dynamic quantities can be carried out by application of the GUM framework. However, there are many mathematical and practical challenges which require specific developments and research. For instance, in practice the uncertainty associated with a static quantity is often determined by repeated measurements. Therefore, for univariate quantities a rather small number of measurements are sufficient. For multivariate quantities, however, the necessary number of measurements increases with the dimension of the quantity. Discretized dynamic quantities are typically very highdimensional data sets, with a typical dynamic measurement consisting of more than one thousand time instants. The evaluation of uncertainty by means of repeated measurements is thus not possible. To this end, parametric approaches have to be determined. Corresponding methods can be found in the field of time series analysis, but their application in metrology requires significant further developments.
Estimating the measurand
In most cases, estimation of the measurand in dynamic measurements requires deconvolution. However, this is a mathematically illposed inverse problem. That is, it requires some kind of regularization in order to obtain reasonable uncertainties. To this end, a typical approach in signal processing is the application of a suitable lowpass filter to suppress undesired high frequency components. In fact, many classical concepts of deconvolution such as Tikhonov regularization or Wiener deconvolution can be interpreted as a successive application of the reciprocal system response and a lowpass filter. However, this is only possible by taking into account prior knowledge about the measurand, which is currently not considered in the GUM and its supplements. The type of prior knowledge can be, for instance, a parametric model or an upper bound in the frequency domain. In any case, the lowpass filter causes a systematic deviation in the estimation result. A reduction in the systematic deviation is always attended by an increase in the noisedependent variance in the result of deconvolution and vice verse. For metrological applications, these systematic errors have to be considered in the uncertainty budget. However, so far no general guidelines or harmonized treatment of these uncertainty contributions is available.
Practical challenges
The treatment of dynamic measurements leads to a number of practical challenges for metrologists. One of the currently most urging challenges is that of transferring the measurement result. Owing to the high dimensionality of discretized dynamic measurements, the associated uncertainties are also highdimensional. For instance, the transfer of the calibration results for the impulse response of a sampling oscilloscope requires transferring a covariance matrix of the dimension of at least $1000\times 1000$. This is not possible with the current practice of paperbased calibration certificates.
Dynamic calibration
A large number of guidelines and standards are available for carrying out and taking advantage of static calibration. The analysis of dynamic measurements, though, is based on a dynamic calibration of the measuring system. The realization of dynamic calibration requires completely different measurement approaches and mathematical tools than static calibration. In some areas first attempts have been made, but in general this constitutes a significant challenge for future research.
For the treatment of the abovementioned challenges, the working group has developed mathematical methodologies for the
 GUMcompliant evaluation of dynamic calibration experiments,
 design of digital deconvolution filters for the compensation of dynamic errors,
 GUMcompliant evaluation of uncertainties for the application of digital filters,
 efficient implementation of GUM Monte Carlo for the application of digital filters,
 iterative estimation of spectral power distributions,
 regularized deconvolution, including GUMcompliant uncertainty evaluation.
Software
In order to facilitate the application of the methods developed in the working group, suitable methods are available as free software implementations. For questions, remarks and suggestions, please contact Sascha Eichstädt.
 PyDynamic  Software for the analysis of dynamic measurements
As part of the EMPIR project 14SIP08 Dynamic the PTB working group 8.42 together with the National Physical Laboratory (UK) develops a comprehensive Python software package, which contains a large variety of methods for the analysis of dynamic measurements.
 Monte Carlo for dynamic measurements
The propagation of measurement uncertainties for dynamic measurements using the Monte Carlo propagation of distributions requires efficient implementation in order to achieve a high degree of accuracy, since the computer memory requirement for a standard implementation is very high. In the working group 8.42 a MATLAB software package has been developed to carry out Monte Carlo for dynamic measurements with high accuracy on standard desktop computers.
 related publication
S. Eichstädt, A. Link, P. Harris and C. Elster (2012). Efficient implementation of a Monte Carlo method for uncertainty evaluation in dynamic measurements. Metrologia 49, 401410
(doi:10.1088/00261394/49/3/401).
 Spectral deconvolution with RichardsonLucy
The correction of deviations in spectra measured with a spectrometer is often necessary in order to obtain accurate results. The classical approach for such a correction is based on a method by Stearns. However, it has been demonstrated that the RichardsonLucy method can produce much better results here. Therefore, the PTB Working Group has written a software implementation of an adapted RichardsonLucy method with an automatic stopping rule. This software includes MATLAB code, a Python code as well as a graphical user interface written in Python.
 related publication
S. Eichstädt F. Schmähling G. Wübbeler, K. Anhalt, L. Bünger, K. Krüger and C. Elster (2013). Comparison of the RichardsonLucy method and a classical approach for spectrometer bandpass correction. Metrologia 50, 107118 (doi: 10.1088/00261394/50/2/107)
 Evaluation of uncertainties in the application of the DFT
The Fourier transform and its counterpart for discrete time signals, the DFT, are common tools in measurement science and application. Although almost every scientific software package offers readytouse implementations of the DFT, the propagation of uncertainties in line with GUM is typically neglected. This is of particular importance in dynamic metrology, when input estimation is carried out by deconvolution in the frequency domain. To this end, we present the new opensource software tool GUM2DFT, which utilizes closed formulas for the efficient propagation of uncertainties for the application of the DFT, inverse DFT and input estimation in the frequency domain. It handles different frequency domain representations, accounts for autocorrelation and takes advantage of the symmetry inherent in the DFT result for realvalued time domain signals.
 related publication
S. Eichstädt and V. Wilkens "GUM2DFT  A software tool for uncertainty evaluation of transient signals in the frequency domain". Meas. Sci. Technol. 27(5), 055001, 2016 (doi: 10.1088/09570233/27/5/055001)
Workshops
Together with the National Physical Laboratory (UK) and the Laboratoire national de métrologie et d'essais (France), the working group is organizing the workshop series "Analysis of Dynamic Measurements" as part of EURAMET TC1078.
 "Signal processing awareness seminar", NPL, UK, 2006
 "Analysis of dynamic measurements", PTB, Germany, 2007
 "Analysis of dynamic measurements", NPL, UK, 2008
 Session TC21 Dynamical Measurements at IMEKO XIX World Congress, Portugal, 2009
 "5th workshop on the analysis of dynamic measurements", SP, Sweden, 2010
 "6th workshop on the analysis of dynamic measurements", Chalmers University, Sweden 2011
 "7th workshop on the analysis of dynamic measurements", LNE, France, 2012
 "8th workshop on the analysis of dynamic measurements" INRIM, Italy, 2014
 "9th International workshop on the analysis of dynamic measurements", PTB, Germany, 2016
Projects
running projects
 20142017) EMRP ENG63 "Sensor network metrology for the determination of electrical grid characteristics" (project website: GridSens)
 (20152018) MNPQ "Reconstruction of spatially resolved color spectra from continuous line scan camera measurements"
 (20152018) EMPIR SIP "Standards and software to maximise end user takeup of NMI calibrations of dynamic force, torque and pressure sensors"
 (since 2008) EURAMET TCIM 1078 "Development of methods for the evaluation of uncertainty in dynamic measurements" abgeschlossene Projekte
completed projects
 (20112014) EMRP IND09 "Traceable Dynamic Measurement of Mechanical Quantities"
 PTBMitteilungen 2015, Issue 2 "Traceable Dynamic Measurement of Mechanical Quantities"
 (20072008) "Modellgestützte Kalibrierung von Beschleunigungsaufnehmern" (Abschlussbericht)
Publications
• 
S. Eichstädt, C. Elster, I. M. Smith and T. J. Esward
J. Sens. Sens. Syst., 6
97105,
2017.
[DOI: 10.5194/jsss6972017]

• 
S. Eichstädt, N. Makarava and C. Elster
Measurement Science and Technology, 27(12),
125009,
2016.

• 
M. Dierl, T. Eckhard, B. Frei, M. Klammer, S. Eichstädt and C. Elster
Journal of the Optical Society of America A, 33(7),
13701376,
2016.
[DOI: 10.1364/JOSAA.33.001370]

• 
S. Eichstädt, V. Wilkens, A. Dienstfrey, P. Hale, B. Hughes and C. Jarvis
Metrologia, 53(4),
2016.
[DOI: 10.1088/00261394/53/4/S125]

• 
S. Eichstädt and V. Wilkens
Measurement Science and Technology, 27(5),
055001,
2016.

• 
S. Eichstädt and C. Elster
tm  Technisches Messen, 83(2),
6670,
2016.
[DOI: 10.1515/teme20150098]

• 
M. Kobusch, S. Eichstädt, L. Klaus and T. Bruns
ACTA IMEKO, 4(2),
4551,
2015.
(Open Access)

• 
L. Klaus, B. Arendacká, M. Kobusch and T. Bruns
ACTA IMEKO, 3(1),
16,
2015.
(Open Access)

• 
C. Matthews, F. Pennecchi, S. Eichstädt, A. Malengo, T. Esward, I. M. Smith, C. Elster, A. Knott, F. Arrhén and A. Lakka
Metrologia, 51(3),
326338,
2014.
[DOI: 10.1088/00261394/51/3/326]

• 
S. Eichstädt, B. Arendacká, A. Link and C. Elster
EPJ Web of Conferences, 77(3),
2014.
[DOI: 10.1051/epjconf/20147700003]

• 
S. Eichstädt and C. Elster
Journal of Physics: Conference Series, 490(1),
012230,
2014.

• 
B. Arendacká, A. Täubner, S. Eichstädt, T. Bruns and C. Elster
Measurement Science Review, 14(2),
5261,
2014.
[DOI: 10.2478/msr20140009]

• 
S. Eichstädt, F. Schmähling, G. Wübbeler, K. Anhalt, L. Bünger, U. Krüger and C. Elster
Metrologia, 50(2),
107118,
2013.
[DOI: 10.1088/00261394/50/2/107]

• 
S. Eichstädt, A. Link, P. M. Harris and C. Elster
Efficient implementation of a Monte Carlo method for uncertainty evaluation in dynamic measurements.
Metrologia, 49(3),
401,
2012.
[DOI: 10.1088/00261394/49/3/401]

• 
T. Bruns, A. Link and A. Täubner
Metrologia, 49(1),
2731,
2012.
[DOI: 10.1088/00261394/49/1/005]

• 
H. Füser, S. Eichstädt, K. Baaske, C. Elster, K. Kuhlmann, R. Judaschke, K. Pierz and M. Bieler
Measurement Science and Technology, 23(2),
025201,
2012.

• 
S. Eichstädt and C. Elster
In F. Pavese, M. Bär, J.R. Filtz, A. B. Forbes, L. Pendrill, K. Shirono, editor, Volume Advanced Mathematical & Computational Tools in Metrology and Testing IX of Series on Advances in Mathematics for Applied Sciences
Chapter 16, page 126135
Publisher: World Scientific New Jersey,
84 edition
, 2012

• 
T. J. Esward, C. Matthews, S. Downes, A. Knott, S. Eichstädt and C. Elster
In F. Pavese, M. Bär, J.R. Filtz, A. B. Forbes, L. Pendrill, K. Shirono, editor, Volume Advanced Mathematical & Computational Tools in Metrology and Testing IX of Series on Advances in Mathematics for Applied Sciences
Chapter 19, page 143151
Publisher: World Scientific New Jersey,
84 edition
, 2012

• 
S. Eichstädt
PhD Thesis
, 2012

• 
S. Eichstädt, A. Link and C. Elster
Sensors, 10(8),
762131,
2010.
[DOI: 10.3390/s100807621]

• 
S. Eichstädt, A. Link, T. Bruns and C. Elster
Measurement, 43(5),
708713,
2010.

• 
S. Eichstädt, C. Elster, T. J. Esward and J. P. Hessling
Metrologia, 47(5),
522533,
2010.
[DOI: 10.1088/00261394/47/5/003]

• 
A. Link and C. Elster
Measurement Science and Technology, 20(5),
055104,
2009.

• 
G. Wübbeler, A. Link, T. Bruns and C. Elster
In F. Pavese, M. Bär, J.M. Limares, C. Perruchet, N.F. Zhang, editor, Volume Advanced Mathematical & Computational Tools in Metrology VIIIof Series on Advances in Mathematics for Applied Sciences
Chapter 52, page 369374
Publisher: World Scientific New Jersey,
78 edition
, 2009

• 
C. Elster and A. Link
In F. Pavese, M. Bär, J.M. Limares, C. Perruchet, N.F. Zhang, editor, Volume Advanced Mathematical & Computational Tools in Metrology VIIIof Series on Advances in Mathematics for Applied Sciences
Chapter 13, page 8089
Publisher: World Scientific New Jersey,
78 edition
, 2009

• 
C. Elster and A. Link
Metrologia, 45(4),
464473,
2008.
[DOI: 10.1088/00261394/45/4/013]

• 
A. Link, A. Täubner, W. Wabinski, T. Bruns and C. Elster
Measurement, 40(910),
928935,
2007.

• 
C. Elster, A. Link and T. Bruns
Measurement Science and Technology, 18(12),
36823687,
2007.
[DOI: 10.1088/09570233/18/12/002]

• 
H.J. von Martens, C. Elster, A. Link, A. Täubner and T. Bruns
Metrologia, 43(1A),
0900209002,
2006.

• 
A. Link, A. Täubner, W. Wabinski, T. Bruns and C. Elster
Measurement Science and Technology, 17(7),
18881894,
2006.
[DOI: 10.1088/09570233/17/7/030]

• 
A. Link, M. Kobusch, T. Bruns and C. Elster
tm  Technisches Messen, 73(12),
675683,
2006.
[DOI: 10.1524/teme.2006.73.12.675]

• 
H.J. von Martens, W. Wabinski, A. Link, H.J. Schlaak, A. Täubner and U. Göbel
Technisches Messen, 72
141152,
2005.
[DOI: 10.1524/teme.72.3.141.60281]

• 
A. Link, W. Wabinski and H.J. von Martens
tm  Technisches Messen, 72(32005),
153160,
2005.
[DOI: 10.1524/teme.72.3.153.60277]

• 
A. Link and H.J. Martens
Measurement, 35(1),
191199,
2004.

• 
A. Link, W. Wabinski and H.J. Martens
Proceedings of SPIE Volume 5503
page 580587
, 2004

• 
C. Elster
Measurement Science and Technology, 11(9),
1359,
2000.
[DOI: 10.1088/09570233/11/9/315]

• 
H.J. von Martens, A. Täubner, W. Wabinski, A. Link and H.J. Schlaak
Measurement, 28(1),
320,
2000.

• 
J. Gerhardt and H.J. Schlaak
tm  Technisches Messen, 6(1),
274282,
2000.
[DOI: 10.1524/teme.2000.67.6.274]

• 
A. Link and H.J. Martens
Shock. Vib., 7(1),
101112,
2000.
[DOI: 10.1155/2000/531718]

• 
A. Link, W. Wabinski, A. Pohl and H.J. von Martens
Proceedings of SPIE
page 126136
, 2000
