Next Article in Journal
Studies of Fluorescence Lifetimes of Biological Warfare Agents Simulants and Interferers Using the Stroboscopic Method
Previous Article in Journal
Tensile Specimen Circular Grid Pattern and AI-Based Strain Calculation Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Method for Estimating the Resultant Expanded Uncertainty Value Based on Interval Arithmetic

Department of Measurement Science, Electronics and Control, Silesian University of Technology, ul. Akademicka 10, 44-100 Gliwice, Poland
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(16), 7334; https://doi.org/10.3390/app14167334
Submission received: 24 June 2024 / Revised: 16 August 2024 / Accepted: 19 August 2024 / Published: 20 August 2024

Abstract

:

Featured Application

Estimation of expanded uncertainty in case of the analysis of an uncertainty budget containing many components, without using Monte Carlo simulations.

Abstract

The article describes a method for determining the resultant expanded uncertainty value in the case of analyzing an uncertainty budget with many components. It employs interval arithmetic, known for its low calculation complexity compared to other methods in the literature. The article includes a detailed explanation of the proposed method and showcases its applications. Results are compared with those from the Monte Carlo method and assumptions based on the central limit theorem. The relative error in estimating the resulting expanded uncertainty value typically stays within ±5% of the Monte Carlo method’s result.

1. Introduction

In both direct and indirect measurements, calculating the value of the measured quantity y ( t ) is always encumbered with a certain error e y ( t ) . This error is typically delineated as the variance between the actual value of the measured quantity y ˜ ( t ) and the ideal value of this quantity y ˙ ( t ) (achieved when the measurement remains unaffected by any phenomenon, commonly referred to as an error source) [1]:
e y t = y ˜ t y ˙ t .
Through repeated measurements, a specific set of instances of the discussed error signal e y ( t ) is acquired. Parameters like the expected value E [ e y ( t ) ] , variance σ e , y 2 , standard deviation σ e , y , and the probability density function’s shape g e , y ( x ) of the signal’s distribution can depict the signal’s characteristics. It is crucial to mention that when multiple measurements are conducted for various values of the measured quantity, the set of error signal instances e y ( t ) encompasses both A and B type components, detailed in the guide [1].
In measurement practice, besides stating the measured quantity’s obtained value, it is essential to specify a quantity that describes the precision with which the value was determined. The two primary measures commonly utilized for this purpose are the standard uncertainty and the expanded uncertainty [1]. In the classical definition of the error signal, as referenced in Equation (1), under the assumption of a symmetric distribution of instances and a signal with a zero expected value, the standard uncertainty equals the signal’s standard deviation. In this scenario, the expanded uncertainty can be calculated using the equation:
P e ^ y t U y = γ ,
where | e ^ y ( t ) | represents the absolute value of any instance of the error signal e y ( t ) , U y signifies the expanded uncertainty value, and γ = 1 α represents the specified confidence level [1], where γ ( 0 ; 100 ) % . The concept behind the expanded uncertainty metric lies in assuming that if the error signal value is equal to or lower than the expanded uncertainty value, this applies to γ across all signal instances. This metric can also be established by correlating its value with the standard uncertainty value:
U y = c e , y σ e , y ,
where c e , y represents the coverage factor derived from the distribution shape of the analyzed error signal implementations and the chosen confidence level [1]. The characteristics of specific distribution types and the connection between the expansion coefficient values and the provided confidence level are extensively discussed in the study by Horálek [2].
The determination of the expanded uncertainty value, as defined in Equation (2), can also be achieved through experimental means. To do so, successive realizations of the signal outlined in Equation (1) need to be acquired via measurement or simulation experiments. Given a discrete set comprising N instances of the studied error signal, and assuming h e , y ( x ) as a function that specifies the count of error signal instances e y ( i ) with a value of x for i [ 0 ; N 1 ] , where the signal’s mean value is zero, Equation (2) is expressed as:
1 N x = U y U y h e , y x = γ ,
for cases where the probability density function of the analyzed signal is positive and symmetrical with respect to the ordinate axes. A graphical representation of Equations (2) and (4) is shown in Figure 1.
Considering the information above, it is evident that the expanded uncertainty metric effectively quantifies the accuracy of determining the measured quantities’ values. Unlike the standard uncertainty metric, it enables the assessment of the likelihood of obtaining an error signal value equal to or less than the assumed expanded uncertainty value, with the flexibility to freely set the confidence level within the range γ ( 0 ; 100 ) % . Nevertheless, the challenge with this metric lies in the intricate calculations involved in determining its final value, especially when multiple error signals perturb the output quantity.
If the uncertainty budget encompasses N error signals e 0 ( t ) , e 1 ( t ) , …, e N 1 ( t ) , with variances σ 0 2 , σ 1 2 , …, σ N 1 2 correspondingly, the standard uncertainty of the signal e Σ ( t ) = e 0 ( t ) + e 1 ( t ) + + e N 1 ( t ) is defined by the equation cited in [1]:
σ Σ = σ 0 σ 1 σ N 1 T 1 r 0 , 1 r 0 , N 1 r 1 , 0 1 r 1 , N 1 r N 1 , 0 1 σ 0 σ 1 σ N 1 ,
where symbol r i , j denotes subsequent Pearson correlation coefficients, determined on the basis of the equation:
r i , j = r j , i = σ i , j 2 σ i 2 σ j 2 2 σ i σ j ,
wherein the mentioned signals may represent both sources of errors related to type A and type B uncertainties.
When dealing with the expanded uncertainty metric, utilizing Equation (5) is unfeasible. Various methods in the literature offer ways to determine the resulting expanded uncertainty value, such as the distributions propagation method [3,4], the expanded uncertainty combination rule method [5,6], and the fuzzy sets-based method [7]. However, these approaches often rely on analytical techniques, necessitating intricate calculations involving convolution operations. In practical measurements, the Monte Carlo method is frequently employed [8,9], or calculations are simplified by assuming compliance with the central limit theorem conditions, as suggested in the guide [1]. Nonetheless, each of these methods has its specific limitations, as elaborated upon later in the text.
Considering the existing knowledge and challenges identified, the objective is to propose a method for determining the expanded uncertainty value with reduced computational complexity (in comparison to the Monte Carlo method) while maintaining an acceptable level of accuracy. The Monte Carlo method, renowned for its precision, serves as the benchmark method in this study. It is anticipated that the values obtained through the proposed method should deviate from the reference values by no more than ±5%. Additionally, the research will present results obtained through the classical approach [1], which assumes adherence to the central limit theorem conditions. Owing to the complexity associated with the other aforementioned methods, their results are not included in the analysis.
The article is structured into six sections. Section 1 serves as an introduction, outlining the work’s motivation, current knowledge, and key assumptions. Section 2 elaborates on the fundamental assumptions of the proposed method in a general context. Section 3 details the process of determining coherence coefficient values for both uncorrelated and correlated error signals. In Section 4, a simulation experiment is conducted to validate the method’s efficacy, presenting the obtained results. Section 5 showcases practical examples applying the discussed method, including one that addresses an existing measurement setup detailed earlier in the work [10]. Section 6 encapsulates the key findings from the study, highlighting the most important advantages and disadvantages of the method.
A crucial component of the project is the literature review, essential for consolidating information related to the analysis method under discussion. However, the primary focus of the study delves into unexplored territories, encompassing experimental descriptions and method modifications aimed at fulfilling the previously mentioned assumptions.

2. Reductive Interval Arithmetic Method

The reductive interval arithmetic method, rooted in interval arithmetic originally attributed to Moore, R. E. [11], has been a subject of study since 1958. Concurrently, Warmus, M. [12], delved into this area since 1956. It was Moore, R. E. who initially demonstrated the practical application of this method. Interval arithmetic facilitates the examination of interval interactions, enabling the treatment of the expanded uncertainty measure as an interval. Consequently, it can be leveraged to compute the cumulative sum of successive expanded uncertainties.
When adhering to the classical definition of expanded uncertainty, as outlined in Equations (2) and (4), this measure can be represented as an unbiased interval with a radius equivalent to its value. Notably, in the realm of expanded uncertainty measures, the summation of two such intervals does not align with the sum derived from classical interval arithmetic, typically resulting in a smaller real value. Jakubiec, J., proposed a modified interval arithmetic approach, termed “reductive interval arithmetic”, denoting the capability to narrow down the resultant interval [13,14,15,16].
Considering an error signal e Σ ( t ) = e 0 ( t ) + e 1 ( t ) + + e N 1 ( t ) comprising N partial error signals, each associated with expanded uncertainties determined at the same confidence level γ = 1 α as U 0 , U 1 , …, U N 1 , respectively, the calculation of the expanded uncertainty using the reduced interval arithmetic method is expressed in the equation cited in [13]:
U Σ = U 0 U 1 U N 1 T 1 h 0 , 1 h 0 , N 1 h 1 , 0 1 h 1 , N 1 h N 1 , 0 1 U 0 U 1 U N 1 ,
where the symbol h i , j = h j , i signifies the coherence coefficient value for the i-th and j-th error signal pair. These coefficients reflect the correlation between the analyzed error signals and the interdependency of their distribution shapes with the resulting error signal implementation’s distribution shape. Establishing these coefficient values is a pivotal undertaking in the application of the method under discussion, detailed further in the subsequent sections of the study.
The discussed method was proposed and used by its author to analyze the uncertainty budget of measurement chains processing series of measurement data. A feature of the indicated measurement chains is the fact that they process multiple successive implementations of error signals originating from the same source, and, therefore, have identical parameters and properties. Unfortunately, due to further research focused on other areas, subsequent works describing the indicated method have not been published, similarly to its further development. Additionally, in the English-language literature there are very few works available on the method of reductive interval arithmetic. One of its applications can be found in the work [17], however, most of the works devoted to the discussed method were published in Polish. Due to the unique features of the indicated method, the authors of this work decided to continue the research conducted by Jakubiec, J., additionally modifying the original algorithm of its application, thus expanding the area of application of the method.

3. Coherence Factors Values Estimation

Analyzing the case that assumes the existence of only two error signals e a ( t ) and e b ( t ) with expanded uncertainties U a and U b , the Equation (7) takes the form:
U a , b = U a U b T 1 h a , b h b , a 1 U a U b = U a 2 + U b 2 + 2 U a U b h a , b ,
and for the known resultant value of the expanded uncertainty U a , b , resulting from the combination of the uncertainties U a and U b , the value of the coherence coefficient h a , b = h b , a can, therefore, be determined according to the following equation:
h a , b = h b , a = U a , b 2 U a 2 U b 2 2 U a U b .
Determining the coefficient value h a , b in the analyzed scenario involves conducting an experiment using the Monte Carlo method. This process entails iteratively summing the implementation values of the studied error signals, expressed as e ^ a , b ( t ) = e ^ a ( t ) + e ^ b ( t ) , and subsequently calculating the resulting expanded uncertainty U a , b following Equation (4) [1,9]. Alternatively, analytical methods requiring convolution operations can be employed to determine the coherence coefficient value, as outlined in [13]. Notably, the coherence coefficient value obtained is only valid for a signal pair with parameters identical to those of the e a ( t ) and e b ( t ) signals, rendering this method impractical in real-world applications.
The essence of employing the interval arithmetic reduction method is to offer a way to ascertain the coherence coefficient values for any number and parameters of error signals without the need for repetitive Monte Carlo simulations or analytical derivations. The approach involves estimating these coefficients appropriately, a process elucidated in subsequent sections of the study. Given that practical measurement scenarios often lack precise realization values and probability density function insights for the identified error signals, assumptions are made in this paper regarding the unknown nature of these parameters.

3.1. Case of Uncorrelated Error Signals

For uncorrelated error signals, the methodology presented in [18] offers a viable solution. This approach involves initially determining the shape factor value s a , b for a signal pair e a ( t ) and e b ( t ) with specified parameters. Subsequently, this factor is adjusted based on the actual parameters of the studied signals. This technique found practical application in the study by [19], with the correction initially proposed in [18]. The correction primarily accounted for the assumptions of the central limit theorem in computing the resultant expanded uncertainty value. Further refinements, including adjustments to address discrepancies in the expanded uncertainty values of the examined error signals, are elaborated upon in detail.
Assuming two uncorrelated error signals e a ( t ) and e b ( t ) , whose expanded uncertainty values are U = U a = U b , the value of the shape factor s a , b can be determined for the pairs of signals discussed in the form:
s a , b = s b , a = U a , b 2 U a 2 U b 2 2 U a U b = U a , b 2 2 U 2 1 ,
where the expanded uncertainty values for both signals need to be determined at the same confidence level γ . It is worth mentioning that with a sufficiently large dataset (approximately 100,000 implementation values) for the signals in question, the coefficient’s value can be estimated using the Monte Carlo method. In cases where equal expanded uncertainty values are assumed ( U a = U b ), the coefficient’s value is influenced solely by the probability density function shape of the signals under analysis. Alternatively, an approach to determining the shape factor values could involve assuming identical standard uncertainties ( σ = σ a = σ b ). Under this assumption, Equation (10) is expressed as:
s a , b = s b , a = U a , b 2 c a 2 σ a 2 c b 2 σ b 2 2 c a σ a c b σ b = U a , b 2 2 c a c b σ 2 1 ,
However, the values obtained in this way may differ from those obtained by adopting the assumption described for Equation (10). The values of shape factors should be determined for pairs of all unique shapes of error signal distributions occurring in the analyzed measurement chain.
Table 1 presents the acquired shape factor values for typical probability density distribution functions, considering the specified confidence level and the method employed for determining these values. Upon examining the data within the table, it becomes evident that the shape factor values escalate with the rise in confidence level. Additionally, there is a trend where values calculated using Equation (10) are generally higher than those derived from Equation (11). The impact of the method used to ascertain shape factor values on the efficacy of the proposed approach for determining the resultant expanded uncertainty value is explored in subsequent sections of the study.
Upon analyzing Equations (9)–(11), it is observed that the shape coefficients determined through the chosen method align with the coherence coefficients for signals sharing identical parameters to those used in the shape factor calculations. Any alterations in signal parameters or their quantity lead to variations in coherence coefficient values. Hence, a methodology for estimating coherence coefficients based on shape coefficient values and the current signal parameters needs to be proposed.
The initial aspect to address pertains to the central limit theorem. As the number of summed independent random variables increases, the probability density distribution of their sum tends towards a normal distribution [1]. Referencing the outcomes in Table 1, it is apparent that the closer the distributions of a pair of summed signals resemble the normal distribution, the closer the shape coefficient value approaches zero. Consequently, as the quantity of summed signals rises, the subsequent coherence coefficients utilized in Equation (7) tend towards zero. In the work [18], it was proposed to correct the shape coefficients according to the equation:
k i , j = k j , i = U i 2 + U j 2 k = 0 N 1 U k 2 ,
while the application of the described correction allows for obtaining results that are overestimated by approximately 5% in relation to those obtained using the Monte Carlo method (only when assuming equal expanded uncertainty values for analyzed error signals).
The second significant phenomenon concerns the impact of variations in the expanded uncertainty values associated with the analyzed signal pair, for which the coherence coefficient is defined in Equation (9). To explore this effect, an experiment was conducted using the Monte Carlo method, determining the coherence coefficient for a signal pair with distinct expanded uncertainty values U a and U b in accordance with Equation (9), where the ratio U a U b ranged from 1 to 10. The experimental findings for specific confidence levels are depicted in Figure 2.
Generally, the absolute value of the coherence coefficient diminishes as the disparity between the expanded uncertainty values of the signal pair increases. However, this relationship is non-linear, as in some instances, for lower expanded uncertainty ratios, the coherence coefficient values increase before declining with larger disparities. Observing the results presented in Figure 2, it can be noticed that in the case when the summed signals are characterized by distributions close to normal, the dependence of the coherence coefficient value as a function of the ratio of the expanded uncertainty values associated with these signals becomes monotonic and begins to approach a linear function. The discussed dependence becomes strongly non-linear and devoid of monotonicity in the case of distributions that are characterized by a large width (in particular in the case of the distribution in the shape of the letter “u” [2]). The indicated circumstances mean that the exact value of the correction taking into account the occurring phenomena should be determined separately for each pair of signals with the given distributions; however, a universal approach is proposed in the paper, discussed later in the section.
Taking into account the conclusions presented, a universal correction is proposed to address the phenomenon in question:
p a , b = p b , a = min U a , U b max U a , U b ,
Nevertheless, it is essential to note that this correction may not consistently capture the dynamics observed in Figure 2.
Lastly, in scenarios involving uncorrelated error signals, the coherence coefficient value can be determined based on the adjusted shape coefficient value:
h i , j = h j , i = s i , j p i , j k i , j = s i , j min U i , U j max U i , U j U i 2 + U j 2 k = 0 N 1 U k 2 ,
while in the case of the original method [18], these values were determined according to the equation:
h i , j = h j , i = s i , j k i , j = s i , j U i 2 + U j 2 k = 0 N 1 U k 2 .

3.2. Case of Correlated Error Signals

Dealing with correlated error signals presents a significantly more intricate scenario. In such instances, the coherence coefficient value is a product of the convolution of the distributions’ shapes and their correlation coefficient. The existing literature [13,15] delves into analytical techniques for determining coherence coefficient values in correlated scenarios. Nonetheless, these methods pertain to specific signal parameters, necessitating a strategy for estimating these coefficients akin to the approach taken for uncorrelated signals.
However, before delving into estimation methods, it is imperative to establish a clear understanding of how a pair of correlated error signals is defined. In experimental settings, a correlated signal pair e a ( t ) and e b ( t ) is characterized by:
e b t = r a , b e a t + 1 r a , b 2 e c t where r a , b > 0 r a , b e a t 1 r a , b 2 e c t where r a , b < 0 ,
where the signal e c ( t ) is a signal with a given distribution, uncorrelated with the signal e a ( t ) , and r a , b = r b , a is the Pearson correlation coefficient of signals e a ( t ) and e b ( t ) , determined according to the Equation (6). It is crucial to acknowledge that the definition outlined in Equation (16) may not precisely mirror the relationships within the analyzed measurement setup. In such instances, the definition’s formulation should be adjusted to align with the specific case under investigation. Subsequently, the dependencies discussed in the subsequent subsection should be determined to adhere to the established methodology.
For a signal pair governed by the relationship in Equation (16), the coherence coefficient value, contingent on the expanded uncertainty values and correlation coefficient, can be computed using Equation (9). Figure 3 illustrates the coherence coefficient values’ dependency for a signal pair with defined probability density function shapes, assuming equivalent standard uncertainties across the signals. Figure 4 showcases the coherence coefficient values’ correlation with varying expanded uncertainty values and correlation coefficients for signal pairs sharing a similar distribution type. In this analysis, the distribution of e a ( t ) aligns with the first type, while the distribution of e c ( t ) , as previously outlined in Equation (16), corresponds to the second type.
Upon scrutinizing the outcomes of the conducted experiments, it is evident that the coherence coefficient value for a correlated signal pair is contingent on both the standard uncertainty ratio of these signals and the Pearson correlation coefficient. Establishing a relationship akin to that described in Equation (14), allowing for a reasonably accurate estimation of the coherence coefficient utilizing the shape coefficients determined from Equation (11), proves to be exceedingly challenging in this scenario. The discussed relationships are solely explored for signal pairs; complexity escalates with an increasing number of error signals.
It can be noticed that in the case of the results presented in Figure 3, one can draw conclusions similar to those discussed earlier for the results presented in Figure 2. Regardless of the adopted confidence level, the value of the coherence coefficient is directly proportional to the value of the Pearson correlation coefficient of the analyzed signals. The main reason for the non-linearity of the discussed characteristics is the fact that in the case of no correlation (for r a , b = 0 ), the value of the coherence coefficient is nonzero for signal pairs whose shape of the resultant distribution of implementation values deviates from the shape of the normal distribution. As depicted in Figure 3, the relationships become notably more non-linear with higher confidence levels. Additionally, based on the data presented in Figure 4, it can be seen that the absolute value of the coherence coefficient decreases with increase in the ratio of the expanded uncertainty values of the analyzed pair of signals. The same property was demonstrated by the experiments, the results of which are presented in Figure 2.
In conclusion, to implement the discussed method effectively, the initial step involves estimating the resultant expanded uncertainty for all sets of correlated error signals. Various approaches can be employed for this purpose, for example, utilizing the method outlined in the study to estimate coherence coefficients based on the insights from Figure 3 and Figure 4, determining correlation coefficient values analytically per the works of [15], or following the method elucidated in the guide [1]. Subsequent to determining the parameters for the resultant error signal in the correlated signal grouping, Equation (7) can be applied for coherence coefficients estimated according to Equation (14). Notably, in scenarios of complete correlation within the analyzed signals, the coherence coefficient will be 1 for full positive and −1 for full negative correlation, serving as a crucial feature of the method.

4. Verification of the Proposed Analysis Method

To assess the efficacy of the proposed method, a plethora of simulation experiments were conducted. These endeavors aimed to elucidate the conditions warranting the method’s application and to ascertain the typical discrepancies between the results derived from the proposed method and those derived from the Monte Carlo method. Typical discrepancies refer to the range encompassing 95% of the error values in estimating the resultant expanded uncertainty, calculated using the relationship:
δ U , = 100 U U s U s e q : u n c e r r o r ,
where U denotes the estimated expanded uncertainty value following Equation (7) or Equation (3), while U s represents the expanded uncertainty value determined based on Equation (4) for the resultant signal error:
e Σ i = e 0 i + e 1 i + + e N 1 i .
The parameters and quantity of signals e i ( i ) varied depending on the experiment variant. Subsequently, 100,000 successive implementation values of this signal were computed for the resulting error signal e Σ ( i ) . Given the intricate nature of correlated signals, all experiments in this section emphasized the absence of correlation among the analyzed error signals.
A single iteration of the experiment therefore consisted of:
  • Generating 100,000 subsequent implementation values for N error signals e 0 ( i ) , e 1 ( i ) , …, e N 1 ( i ) with selected parameters;
  • Determining the implementation value of the resulting error signal e Σ ( i ) , according to Equation (18);
  • Determining the expanded uncertainty value U s for the signal e Σ ( i ) , according to Equation (4), for the confidence level γ ;
  • Determining the value of the expanded uncertainty U c , 1 for the signal e Σ ( i ) , according to Equation (7), using the shape factor values previously determined according to Equation (10) and coherence coefficients estimated according to Equation (14), for the confidence level γ ;
  • Determining the value of the expanded uncertainty U c , 2 for the signal e Σ ( i ) , according to Equation (7), using the shape factor values previously determined according to Equation (11) and coherence coefficients estimated according to Equation (14), for the confidence level γ ;
  • Determining the value of the expanded uncertainty U c , 3 for the signal e Σ ( i ) , according to Equation (3), assuming that the distribution of the implementation of this signal is close to the normal distribution [1];
  • Determining the value of the relative expanded uncertainty estimation error δ U , 1 for U c , 1 , δ U , 2 for U c , 2 and δ U , 3 for U c , 3 according to Equation (17).
For different variants of the experiment:
  • A selected confidence level of γ was used, constant for each iteration of the experiment implementation;
  • For a single iteration, the number of error signals was drawn from the interval N [ N m i n ; N m a x ] with the same probability of obtaining each possible value;
  • For each error signal e j ( i ) , the expanded uncertainty value was drawn from the interval U j [ U m i n ; U m a x ] with the same probability of obtaining each possible value associated with this signal;
  • The shape of the distribution of the signal implementation value e j ( i ) was each time drawn from a pool of distributions: normal, uniform, triangular, u-shape (sine function distribution), and each type of distribution was equally likely to be drawn.
  • The values N m i n , N m a x , U m i n , and U m a x were constant for all iterations of the implemented variant of the experiment;
Upon deriving the relative error values δ U , for estimating the expanded uncertainty U in each experiment variant, the average, standard deviation, and expanded uncertainty values U and U + were computed for a confidence level of 95%. Table 2 showcases specific results from the experiment variants, while Figure 5 displays histograms depicting the distribution of δ U , values obtained.
Analyzing the presented results, it can be noticed that the use of shape coefficients determined on the basis of Equation (11) causes the obtained uncertainty values to be smaller than those determined according to Equation (10). For selected cases, this allows to obtain more accurate results, while assessing the entire experimental results, it is concluded that it is more beneficial to use the coefficients determined in accordance with Equation (10). In measurement practice, it is a much better approach to slightly overestimate the uncertainty value than to underestimate it. Due to the above, the subsequent results presented will only concern the case of using Equation (10) and will be compared to those obtained by the Monte Carlo method and those obtained assuming the fulfillment of the conditions of the central limit theorem. Figure 6 shows selected experimental results for the indicated calculation methods.
Analyzing the results presented in Figure 6, it can be noticed that the shapes of the histograms obtained for the (a) and (b) methods are almost identical. A significant difference is the fact that for the (a) case, a higher expected value can be observed than for the (b) case. The indicated property means that the estimated values of the resultant expanded uncertainty will be higher for the (a) method. In the case of using the (c) method, a much larger spread of the implementation values of the relative error in estimating the resultant expanded uncertainty can be noticed than in the case of the methods (a) and (b). The discussed property is particularly visible in the case when the analyzed uncertainty budget consists of a small number of error sources and there is no dominant source of error, i.e., in the case when the conditions of the central limit theorem have not been met. In the case of all analyzed methods, the obtained distributions of the implementation of the relative error in the estimation of the resultant value of the expanded uncertainty are asymmetric, which results directly from the features of the distributions of the error signal components, the properties of which were discussed earlier by analyzing Figure 2, Figure 3 and Figure 4.

5. Example of the Proposed Analysis Method Application

To exemplify the practical application of the analysis method in question, two examples are introduced later in this section. The initial example is entirely theoretical and does not pertain to the actual measurement process. The aim of this example is to provide minimal detail and concentrate solely on elucidating the discussed method. The second example was drawn from the established measurement setup previously detailed in the study [10]. It is essential to acknowledge that the results obtained through the proposed method will not be juxtaposed with outcomes from other methods, as this comprehensive comparison was extensively outlined in the preceding section.

5.1. First Example

In the example, it is assumed that the time-varying physical quantity s ( t ) undergoes processing by a measuring transducer, which outputs a representation of this quantity as a time-varying electrical voltage y 1 ( t ) . Subsequently, this voltage is amplified by an amplifier and input as a signal y 2 ( t ) to the A/D converter. The converter then takes consecutive implementations of the input quantity and converts them into discrete samples represented as x ( i ) . The diagram depicting the discussed measurement pathway is illustrated in Figure 7.
It was assumed that the measuring transducer introduces white noise disturbance into the output signal, with the variance transferred to the value x ( i ) being σ w 2 , and the realization values of this signal following a normal distribution. The measurement amplifier introduces an error related to zero shift due to temperature variations into the output quantity, with the variance transferred to the quantity x ( i ) being σ z 2 , and subsequent implementations following a triangular distribution. Subsequently, the A/D converter introduces quantization-related error into the value x ( i ) , with a variance of σ q 2 , and succeeding implementations following a uniform distribution.
Given the outlined assumptions, the vector comprising expanded uncertainties for the specified quantity x ( i ) can be expressed as:
U = U w U z U q = c w σ n c z σ z c q σ q ,
while the coherence matrix takes the form:
h = 1 h w , z h w , q h z , w 1 h z , q h q , w h q , z 1 ,
whereby:
h w , z = h z , w = s n , t min U w , U z max U w , U z U w 2 + U z 2 U w 2 + U z 2 + U q 2 ,
h w , q = h q , w = s n , u min U w , U q max U w , U q U w 2 + U q 2 U w 2 + U z 2 + U q 2 ,
h z , q = h q , z = s t , n min U z , U q max U z , U q U z 2 + U q 2 U w 2 + U z 2 + U q 2 ,
and finally, based on Equation (7):
U Σ = U · h · U T .

5.2. Second Example

In the second example, the analysis method was employed to ascertain the resulting expanded uncertainty value of the measurement chain outlined in the study [10], specifically when the mentioned measurement chain handles a polyharmonic input signal. As per the referenced work, the uncertainty vector for successive error signals is structured as:
U = U x , r U x , d , 1 U x , d , 2 U x , d , N ,
where the number N of the components of the indicated vector depends on the number of analyzed harmonics of the dynamic error signal, the subsequent variances of which are described by the equation:
σ x , d , i 2 = E s , o , i cos φ ˜ y ω s , o , i E s , o , i 2 + E s , o , i sin φ ˜ y ω s , o , i 2 2 ,
φ ˜ y ω 6.26 · 10 13 ω 2 5.73 · 10 7 ω ,
where the amplitude of the selected harmonic
E s , o , i = π 8 2 i 1 2 E s , o ,
depends on the amplitude E s , o of the processed signal [10]. As indicated in the cited work, the harmonics of the dynamic error signal are analyzed whose pulsation meets the condition ( 2 i 1 ) ω x , o < 1 2 ω s for the sampling pulsation ω s = 2 π 48,000   rad . The expanded uncertainty values for the analyzed signals are as follows:
U x , r = c n σ x , r ,
U x , d , i = c d σ x , d , i ,
where σ x , r 2 = 0.13   μ V [10].
Analyzing the case where the pulsation of the signal s ( t ) processed by the discussed measurement chain is ω s , o = 2 π 4000   rad , the vector described in Equation (25) takes the form:
U = U x , r U x , d , 1 U x , d , 2 U x , d , 3 U x , d , 4 ,
while the matrix of coherence coefficients takes the form:
h = 1 h r , d , 1 h r , d , 2 h r , d , 3 h r , d , 4 h d , 1 , r 1 h d , 1 , 2 h d , 1 , 3 h d , 1 , 4 h d , 2 , r h d , 2 , 1 1 h d , 2 , 3 h d , 2 , 4 h d , 3 , r h d , 3 , 1 h d , 3 , 2 1 h d , 3 , 4 h d , 4 , r h d , 4 , 1 h d , 4 , 2 h d , 4 , 3 1 ,
where the subsequent values of its elements are determined following Equation (14), akin to the previous example.
Ultimately, the calculation of the resulting expanded uncertainty value is carried out in accordance with Equation (7), after substituting the values specified in Equations (31) and (32). The same analysis method was also applied in the study [19], although it did not incorporate the correction proposed in Equation (13), leading to less precise results compared to utilizing the suggested method version.

6. Conclusions

The experiments conducted demonstrate that the method for determining the resulting expanded uncertainty yields results nearly identical to those obtained through the Monte Carlo method. This method is applicable for the specified confidence level and exhibits minimal sensitivity to its conditions of use, such as the number of error signals or the range of expanded uncertainty values of these signals. In most scenarios, the results obtained using this method prove more accurate than those derived from the approach outlined in [1]. Notably, the method boasts a straightforward application algorithm and low computational complexity. However, it is essential to highlight that when accepting an error in estimating the resultant expanded uncertainty exceeding approximately ±15%, the method delineated in the guide [1] should be preferred.
A notable drawback of the method lies in the prerequisite of initially determining the shape coefficient values, necessitating knowledge of all probability density function shapes for the analyzed error signals. Despite this, the process is not intricate and allows for analytical, simulation-based, and measurement-based approaches. The determined coefficients can be stored in a microcontroller’s memory for subsequent estimation of coherence coefficient values, considering the current parameters of the analyzed error signals. Challenges arise when dealing with correlated error signals, where it is suggested to analyze correlated signal groups separately to determine the resultant expanded uncertainty value for subsequent consideration in the analysis.
The key achievements of the authors of this work include:
  • Proposing a correction described in Equation (12), thanks to which the results obtained are much more accurate compared to those described so far in the literature [18];
  • Indication of the difference in the effectiveness of the method for different ways of determining the shape coefficient values, described in Equations (10) and (11);
  • Examining the effectiveness of the proposed method depending on the value of the confidence level parameter, number of error signals and their expanded uncertainty range;
  • Comparison of the effectiveness of the method to the classical approach proposed in the guide [1].
In conclusion, the employed analysis method is justified when precise determination of the resultant expanded uncertainty is crucial, irrespective of the number of error signals or the chosen confidence level. Its seamless integration within the measurement system’s software and minimal computational demands make it a practical choice. Notably, alterations in signal parameters entail no additional calculations, underscoring a key advantage of this method.
Future research endeavors could involve refining the method for estimating shape coefficients with greater precision than outlined in this study. Additionally, exploring the analysis of correlated error signals, particularly when the signal count exceeds two, presents a significant area for investigation. One potential approach could involve a strategy akin to that proposed in Equation (14), where a function determining the coherence coefficient value for each pair of distributions individually is derived, or delving into the feasibility of employing artificial neural networks for coherence coefficient estimation.
To facilitate user understanding of the method’s functionality, an illustrative application example has been developed for the GNU Octave (Authors used GNU Octave version 9.2.0, Copyright (C) 1993–2024 The Octave Project Developers) program [20]. This example enables users to replicate the results demonstrated in the study and gain familiarity with the method’s operational algorithm. The example encompasses shape coefficient values for confidence levels γ [ 50 ; 99 ] % and γ N . Users can access the example on GitHub [21].

Author Contributions

Data curation, Ł.D.; formal analysis, M.K. and J.R.; funding acquisition, M.K.; investigation, Ł.D.; methodology, M.K. and Ł.D.; project administration, M.K.; resources, Ł.D.; software, Ł.D.; supervision, J.R.; validation, M.K. and J.R.; writing—original draft, Ł.D.; writing—review and editing, M.K. and J.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by the Polish National Science Centre (NCN), grant number 2022/47/B/ST7/00047, and by the Rector of the Silesian University of Technology, grant number 05/020/RGJ24/0084.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors upon request. Most of the used data can be replicated using example [21].

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
A/AAnalog-to-Analog
A/DAnalog-to-Digital
CLTCentral Limit Theorem
RIAReductive Interval Arithmetic

References

  1. Joint Committee for Guides in Metrology. Evaluation of Measurement Data—Guide to the Expression of Uncertainty in Measurement; JCGM: Sèvres, France, 2008; Available online: https://www.bipm.org/documents/20126/2071204/JCGM_100_2008_E.pdf (accessed on 19 August 2024).
  2. Horálek, V. Analysis of basic probability distributions, their properties and use in determining type B evaluation of measurement uncertainties. Measurement 2013, 46, 16–23. [Google Scholar] [CrossRef]
  3. Koliander, G.; El-Laham, Y.; Djurić, P.M.; Hlawatsch, F. Fusion of probability density functions. Proc. IEEE 2022, 110, 404–453. [Google Scholar] [CrossRef]
  4. Zhang, Z.; Wang, J.; Jiang, C.; Huang, Z.L. A new uncertainty propagation method considering multimodal probability density functions. Struct. Multidiscip. Optim. 2019, 60, 1983–1999. [Google Scholar] [CrossRef]
  5. Dieck, R.H. Measurement Uncertainty: Methods and Applications; ISA: Durham, NC, USA, 2017. [Google Scholar]
  6. Yang, L.; Guo, Y. Combining pre-and post-model information in the uncertainty quantification of non-deterministic models using an extended Bayesian melding approach. Inf. Sci. 2019, 502, 146–163. [Google Scholar] [CrossRef]
  7. Urbanski, M.K.; Wąsowski, J. Fuzzy approach to the theory of measurement inexactness. Measurement 2003, 34, 67–74. [Google Scholar] [CrossRef]
  8. Janssen, H. Monte-Carlo based uncertainty analysis: Sampling efficiency and sampling convergence. Reliab. Eng. Syst. Saf. 2013, 109, 123–132. [Google Scholar] [CrossRef]
  9. Joint Committee for Guides in Metrology. Evaluation of Measurement Data—Propagation of Distributions Using a Monte Carlo Method; JCGM: Sèvres, France, 2008; Available online: https://www.bipm.org/documents/20126/2071204/JCGM_101_2008_E.pdf (accessed on 19 August 2024).
  10. Kampik, M.; Roj, J.; Dróżdż, L. Error Model of a Measurement Chain Containing the Discrete Wavelet Transform Algorithm. Appl. Sci. 2024, 14, 3461. [Google Scholar] [CrossRef]
  11. Moore, R.E.; Kearfott, R.B.; Cloud, M.J. Introduction to Interval Analysis; Other Titles in Applied Mathematics; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 2009. [Google Scholar]
  12. Warmus, M. Calculus of approximations. Bull. l’Académie Pol. Sci. Cl. 3 1956, 4, 253–259. [Google Scholar]
  13. Jakubiec, J. Application of Reductive Interval Arithmetic to Uncertainty Evaluation of Measurement Data Processing Algorithms; Wydawnictwo Politechniki Śląskiej: Gliwice, Poland, 2002. [Google Scholar]
  14. Jakubiec, J.; Konopka, K. Reducing interval arithmetic in dynamic error evaluation. In Proceedings of the XVI IMEKO World Congress, Vienna, Austria, 25–28 September 2000; pp. 100–105. [Google Scholar]
  15. Jakubiec, J. Reductive interval arithmetic application to uncertainty calculation of measurement result burdened correlated errors. Metrol. Meas. Syst. 2003, 10, 137–156. [Google Scholar]
  16. Jakubiec, J.; Konopka, K. Coherence coefficients as uncertainty parameters of error value set. In Proceedings of the IMEKO TC7 Symposium—Measurement Science of the Information Era, Cracow, Poland, 25–27 June 2002; pp. 76–81. [Google Scholar]
  17. Batko, W.; Pawlik, P. Uncertainty evaluation in modelling of acoustic phenomena with uncertain parameters using interval arithmetic. Acta Phys. Pol. A 2012, 121, A-152–A-155. [Google Scholar] [CrossRef]
  18. Jakubiec, J. Błędy i Niepewności Danych w Systemie Pomiarowo-Sterującym; Wydawnictwo Politechniki Śląskiej: Gliwice, Poland, 2010. [Google Scholar]
  19. Kampik, M.; Roj, J.; Dróżdż, L. Estimation of the resultant expanded uncertainty of the output quantities of the measurement chain using the discrete wavelet transform algorithm. Appl. Sci. 2024, 14, 3691. [Google Scholar] [CrossRef]
  20. Eaton, J.W. GNU Octave Official Website. 2024. Available online: https://octave.org (accessed on 19 August 2024).
  21. Dróżdż, L. Reductive Interval Arithmetic Method Application Example on GitHub. 2024. Available online: https://github.com/Kuszki/Octave-Uncertainty-RIA-Example (accessed on 19 August 2024).
Figure 1. Graphical representation of the expanded uncertainty measure U y for the γ = 1 α confidence level assuming a zero expected value and a symmetric distribution of subsequent error signal e y ( i ) implementations.
Figure 1. Graphical representation of the expanded uncertainty measure U y for the γ = 1 α confidence level assuming a zero expected value and a symmetric distribution of subsequent error signal e y ( i ) implementations.
Applsci 14 07334 g001
Figure 2. The value of the coherence coefficient depending on the ratio of the expanded uncertainty values of the analyzed signals (case for uncorrelated signals).
Figure 2. The value of the coherence coefficient depending on the ratio of the expanded uncertainty values of the analyzed signals (case for uncorrelated signals).
Applsci 14 07334 g002
Figure 3. The value of the coherence coefficient depending on the value of the correlation coefficient (case for equal values of the standard deviation of signals).
Figure 3. The value of the coherence coefficient depending on the value of the correlation coefficient (case for equal values of the standard deviation of signals).
Applsci 14 07334 g003
Figure 4. The value of the coherence coefficient depending on the value of the correlation coefficient and the ratio of the standard uncertainty values of the analyzed signals with the same shapes of the probability density function.
Figure 4. The value of the coherence coefficient depending on the value of the correlation coefficient and the ratio of the standard uncertainty values of the analyzed signals with the same shapes of the probability density function.
Applsci 14 07334 g004
Figure 5. Obtained histograms for the implementation of the relative error in estimating the expanded uncertainty value, corresponding to the data listed in Table 2.
Figure 5. Obtained histograms for the implementation of the relative error in estimating the expanded uncertainty value, corresponding to the data listed in Table 2.
Applsci 14 07334 g005
Figure 6. Dependence of the average value of the relative error in estimating the expanded uncertainty value depending on the confidence level for the analyzed method (a) (left side) and assuming the fulfillment of the conditions of the central limit theorem (c) (right side), along with the expanded uncertainty for these quantities, determined for the confidence level γ = 95% (the estimated area is filled with color).
Figure 6. Dependence of the average value of the relative error in estimating the expanded uncertainty value depending on the confidence level for the analyzed method (a) (left side) and assuming the fulfillment of the conditions of the central limit theorem (c) (right side), along with the expanded uncertainty for these quantities, determined for the confidence level γ = 95% (the estimated area is filled with color).
Applsci 14 07334 g006
Figure 7. Diagram of an example measurement chain.
Figure 7. Diagram of an example measurement chain.
Applsci 14 07334 g007
Table 1. Summary of the values of the shape ratio s a , b coefficient for the selected confidence level, where the following symbols denote the distribution: ( n ) normal, ( u ) uniform, ( t ) triangular, ( d ) u-shape (sine function distribution).
Table 1. Summary of the values of the shape ratio s a , b coefficient for the selected confidence level, where the following symbols denote the distribution: ( n ) normal, ( u ) uniform, ( t ) triangular, ( d ) u-shape (sine function distribution).
γ Method s a , b nutd
75% U a = U b  (10)n0.00000−0.07510−0.03966−0.07214
u−0.07510−0.11115−0.11109−0.08076
t−0.03966−0.11109−0.07710−0.10742
d−0.07214−0.08076−0.10742−0.02616
σ a = σ b  (11)n0.00000−0.07926−0.04082−0.07559
u−0.07926−0.11098−0.11099−0.08057
t−0.04082−0.11099−0.07709−0.10640
d−0.07559−0.08057−0.10640−0.02628
85% U a = U b  (10)n0.00000−0.00828−0.029850.06124
u−0.008280.03912−0.026000.14582
t−0.02985−0.02600−0.057330.05345
d0.061240.145820.053450.28045
σ a = σ b  (11)n0.00000−0.00835−0.030490.05931
u−0.008350.03909−0.026190.14242
t−0.03049−0.02619−0.057280.04806
d0.059310.142420.048060.28052
95% U a = U b  (10)n0.000000.131490.014970.28894
u0.131490.335860.177350.52378
t0.014970.177350.040320.35107
d0.288940.523780.351070.71555
σ a = σ b  (11)n0.000000.124730.015060.25356
u0.124730.335780.172190.51584
t0.015060.172190.040110.32181
d0.253560.515840.321810.71555
Table 2. Values of the relative error in estimating the expanded uncertainty value for selected experimental conditions.
Table 2. Values of the relative error in estimating the expanded uncertainty value for selected experimental conditions.
γ N U j (a) Coeficients  s a , b  for case (10)(b) Coeficients  s a , b  for case (11)(c) Assumming CLT is met [1]
U 1 , U 1 , + δ ¯ U , 1 σ δ , 1 U 2 , U 2 , + δ ¯ U , 2 σ δ , 2 U 3 , U 3 , + δ ¯ U , 3 σ δ , 3
75%[2; 5][1; 5]−1.75+3.44+0.811.32−1.85+3.43+0.761.34−7.25+1.09−2.642.16
[1; 10]−1.52+3.39+0.811.30−1.58+3.39+0.761.32−7.77+0.87−2.882.35
[1; 20]−1.61+3.35+0.811.29−1.69+3.34+0.761.31−8.04+0.82−3.042.50
[1; 30]−1.70+3.33+0.811.28−1.79+3.32+0.761.30−8.15+0.93−3.102.56
[2; 9][1; 5]−1.82+3.39+0.721.33−1.87+3.37+0.651.35−5.91+0.84−1.851.77
[1; 10]−1.95+3.42+0.721.34−1.98+3.40+0.661.35−6.45+1.04−2.061.94
[1; 20]−1.70+3.45+0.731.34−1.78+3.43+0.671.36−6.91+0.77−2.182.07
[1; 30]−1.97+3.41+0.721.33−2.07+3.39+0.661.35−6.95+0.77−2.212.10
85%[2; 5][1; 5]−2.05+4.57+1.201.64−2.05+4.44+1.121.61−3.91+1.46−1.221.24
[1; 10]−2.01+4.22+1.091.52−2.01+4.10+1.021.50−4.28+1.88−1.201.38
[1; 20]−1.99+4.08+1.041.48−1.99+3.96+0.981.46−4.52+2.17−1.181.48
[1; 30]−1.94+3.96+1.011.45−1.96+3.84+0.941.43−4.64+2.28−1.181.53
[2; 9][1; 5]−2.01+5.83+1.901.99−2.01+5.70+1.811.96−2.72+0.85−0.930.97
[1; 10]−1.96+5.66+1.791.92−1.96+5.52+1.701.90−3.15+1.24−0.951.07
[1; 20]−2.03+5.51+1.711.88−2.05+5.37+1.631.85−3.36+1.43−0.971.14
[1; 30]−1.89+5.51+1.691.88−1.89+5.37+1.611.85−3.48+1.54−0.971.17
95%[2; 5][1; 5]−2.80+3.07+0.141.34−3.31+2.73−0.291.35−1.12+10.69+3.173.70
[1; 10]−2.96+3.03+0.031.35−3.45+2.71−0.371.35−0.88+13.53+3.714.66
[1; 20]−3.02+2.97−0.031.34−3.49+2.65−0.421.35−0.93+15.36+4.075.29
[1; 30]−3.07+2.97−0.051.35−3.54+2.67−0.431.36−0.89+16.15+4.235.57
[2; 9][1; 5]−2.00+4.09+1.041.49−2.54+3.61+0.541.46−0.90+7.32+2.112.85
[1; 10]−2.16+3.99+0.911.48−2.68+3.55+0.431.44−0.75+8.99+2.473.58
[1; 20]−2.26+3.94+0.841.47−2.76+3.48+0.361.43−0.92+10.23+2.714.08
[1; 30]−2.33+3.94+0.811.48−2.83+3.51+0.341.44−0.97+10.74+2.804.27
Mean * −2.10+3.94+0.891.48−2.30+3.76+0.821.48−3.95+4.71+2.242.56
* in the case of δ ¯ U , , the mean of the absolute values was calculated.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kampik, M.; Roj, J.; Dróżdż, Ł. A Method for Estimating the Resultant Expanded Uncertainty Value Based on Interval Arithmetic. Appl. Sci. 2024, 14, 7334. https://doi.org/10.3390/app14167334

AMA Style

Kampik M, Roj J, Dróżdż Ł. A Method for Estimating the Resultant Expanded Uncertainty Value Based on Interval Arithmetic. Applied Sciences. 2024; 14(16):7334. https://doi.org/10.3390/app14167334

Chicago/Turabian Style

Kampik, Marian, Jerzy Roj, and Łukasz Dróżdż. 2024. "A Method for Estimating the Resultant Expanded Uncertainty Value Based on Interval Arithmetic" Applied Sciences 14, no. 16: 7334. https://doi.org/10.3390/app14167334

APA Style

Kampik, M., Roj, J., & Dróżdż, Ł. (2024). A Method for Estimating the Resultant Expanded Uncertainty Value Based on Interval Arithmetic. Applied Sciences, 14(16), 7334. https://doi.org/10.3390/app14167334

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop