Next Article in Journal
A Through-Life Cost Analysis Model to Support Investment Decision-Making in Concentrated Solar Power Projects
Previous Article in Journal
Modeling the Charging Behaviors for Electric Vehicles Based on Ternary Symmetric Kernel Density Estimation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Monte Carlo-Based Procedure for Determining the Maximum Energy at the Output of Accelerometers

by
Krzysztof Tomczyk
Faculty of Electrical and Computer Engineering, Cracow University of Technology, Warszawska 24, 31-155 Krakow, Poland
Energies 2020, 13(7), 1552; https://doi.org/10.3390/en13071552
Submission received: 13 February 2020 / Revised: 16 March 2020 / Accepted: 26 March 2020 / Published: 27 March 2020
(This article belongs to the Section I: Energy Fundamentals and Conversion)

Abstract

:
The solutions presented in this paper can be the basis for mutual comparison of different types of accelerometers produced by competing companies. An application of a procedure based on the Monte Carlo method to determine the maximum energy at the output of accelerometers is discussed here. The fixed-point algorithm controlled by the Monte Carlo method is used to determine this energy. This algorithm can only be used for the time-invariant and linear measurement systems. Hence, the accelerometer nonlinearities are not considered here. The mathematical models of the accelerometer and the special filter, represented by the relevant transfer functions, are the basis for the above procedure. Testing results of the voltage-mode accelerometer of type DJB A/1800/V are presented here as an example of an implementation of the solutions proposed. Calculation of the energy was executed in Mathcad 14 program with the built-in Programming Toolbar. The value of the maximum output energy determined for a specified time interval corresponds to the maximum integral-square error of the accelerometer. Such maximum energy can be a comparative ratio just like the accuracy class in the case of instruments used for the static measurements. Hence, the main analytical and technical contributions of this paper concern the development of theoretical procedures and the presentation of their application on the example of a real type of accelerometer.

1. Introduction

The basic condition for mutual comparison of different types of accelerometers [1,2,3] is that the accelerometers under test have the same frequency bandwidth and it is necessary to develop the model of special filter [4] which is the reference for such an analysis. The frequency bandwidth of the filter must be the same as the operation frequency range of the compared accelerometers.
The reference of energies at the outputs of the accelerometer and the filter is necessary because, from the point of view of the practical implementation of measurements, only the value of energy for the operation range of the accelerometer, not for the full frequency bandwidth, is interesting [5]. Such comparison is carried out by subtracting the energy value at the outputs of both accelerometers. The difference between these energies is defined below as the reference energy and, for the purpose of determining its maximum value, it is convenient to apply the algorithm discussed in Section 3.
The reference energy of an accelerometer corresponds to the integral-square error [6] of the sensor related to this type of error at the output of the special filter. Developing the procedure determines the maximum value of the integral-square error (maximum reference energy) and seems to be the justified proposal for the comparison of different accelerometers.
The basis for the implementation of the task described above is the mathematical model of the accelerometer determined through its parametric identification [7,8,9,10,11,12], as well as the mathematical model of the special filter. The model of the accelerometer is defined by the second-order transfer function [1,2,3], whereas the model of the filter should be a high-order (a range of 10 to 15 order is sufficient). Determination of the parameters of the accelerometer model should be carried out on the basis of measurement points of both components of the frequency response (the amplitude and phase) in accordance with the procedures included in the relevant international guides [7]. In this way (e.g., by using the weighted least squares procedure (WLSP)), the accelerometer parameters and associated uncertainties can be obtained [3,13,14]. Additionally, when we want to determine the associated uncertainties, it is necessary to use the Monte Carlo procedure (MCP) [15,16,17,18], whose number of trials can be assumed in advance according to the guidelines included in the JCGM (Joint Committee for Guides in Metrology)—supplement 1 [13]. The model of the filter is expected to meet the so-called non-distortion transformation [19], according to which the amplitude-frequency and phase-frequency responses are flat and linear, respectively, within the frequency bandwidth of interest [20]. However, these assumptions are ideal from a practical point of view. In practice, they are not possible to meet, and one can only strive to attain them [20], for example, by using the above-mentioned high-order filter. The degree of achievement for the approximation therefore depends on the type of filter used and has an effect on the value of the reference energy (i.e., the higher the order of the approximation, the lower value of this energy). Hence, there is the obvious necessity to use the same model of the filter within the considered comparison of the different types of accelerometers.
The maximum value of the reference energy can be determined by using the fixed-point algorithm (FPA) [5,21,22,23] with its extension by the MCP [13,18]. The MCP can be implemented based on the same assumptions as in the case of the procedure intended for determining the uncertainties associated with the parameters of the accelerometer model by using the WLSP. A method based on a system of complex convolution equations could also be applied instead of the above-mentioned FPA [4]. However, the possibilities of its application are limited due to the numerical difficulties of solving a system including more than twenty-five equations [24]. Therefore, for the needs of this paper, the first type of approach is chosen. This choice is also supported by the fact that by employing the FPA, the maximum reference energy can be obtained quickly. The FPA also allows determination of the critical case of the input signal of the accelerometer that is limited both in magnitude and in time interval [25], and which produces this energy. Both the number of reversals of such a signal and the times at which this reversal occurs are determined. It should be emphasized here that any other signal contained in these limits generates less energy than that produced by the critical case of the magnitude-limited signal [4].
The FPA determines, by an application of the MCP, the numerical value of the maximum reference energy in an iterative way for the ranges of variability of the model parameter of the accelerometer by the values of associated uncertainties [26,27]. Taking into account the uniform distribution of parameters in the above ranges, the Monte Carlo simulation (MCS) should therefore be based on a uniform pseudo-random number generator [28].
The assessment of the impact of the uncertainties associated with the parameters of the accelerometer on the value of the absolute error is presented in studies by Tomczyk [26,27]. However, this assessment does not apply to full ranges of parameter variability on the values of uncertainties associated with them. Additionally, the cited assessment concerns only the values of parameters increased or decreased by the values of the uncertainties associated with them. Therefore, it was not necessary to use the MCP in such a case.
The above-mentioned issues are developed in this paper using the MCP. The block diagrams with an exact description of both the proposed MCP as well as the algorithm for determining the maximum reference energy at the output of accelerometers are discussed here. This maximum reference energy can be the basis for mutual comparison [4] of accelerometers with the same bandwidth but produced by different companies. This comparison relates to the maximum integral-square error, which reflects an accuracy of the accelerometer. It can therefore help one choose the proper accelerometer from many types at different prices.
The practical solutions presented here are applied to the testing of the DJB A/1800/V type of accelerometer [29].

2. Accelerometer Modeling

The model of the accelerometer is defined by using the transfer function,
K ( s ) = S a ω 0 2 s 2 + 2 β ω 0 s + ω 0 2 ,
where S a and β are the accelerometer voltage sensitivity and damping factor, respectively, while ω 0 = 2 π f 0 , where f 0 denotes the non-damped natural frequency.
Parametric identification of the accelerometer model is directed towards determining the parameters of the transfer function given by Equation (1) together with associated uncertainties [8,9]. The first step in this subject is the determination of the measurement points for both frequency responses. This is achieved by applying the back-to-back method commonly used in the measurement practice. The block diagram for this type of identification method is shown in Figure 1. The reference accelerometer is used here as the comparison signal source, while the data acquisition card together with the LabVIEW software are employed for the control of the shaker and recording of the measurement data.
The sensitivity S a of the tested accelerometer is determined on the basis of the relation,
S a = S r V a V r ,
where S r , V a , and V r are the sensitivity of the reference accelerometer, the voltage output of the system for signal conditioning at the output of the accelerometer, and the voltage output of the reference, respectively.
Based on the measurement points of the frequency responses of the accelerometer, the estimates S ˜ a , f ˜ 0 , and β ˜ are associated with the parameters S a , f 0 , and β which occur in the transfer function given by Equation (1), as well as the associated uncertainties u ( S ˜ a ) , u ( f ˜ 0 ) , and u ( β ˜ ) , which are determined by applying the WLSP [3,13]. This procedure employs the MCP which utilizes the Box–Muller and the Wichmann–Hill pseudo-random number generator with normal and uniform distributions, respectively. The minimum number M of Monte Carlo trials (MCts) should be determined based on the relation,
M > 10 4 1 p ,
where p denotes the assumed coverage probability [13].

3. Algorithm for Determining the Maximum Energy (ME)

The optimization algorithm for determining the ME is based on the integral equation,
0 T x ( v ) Φ ( τ , v ) d v = ξ · x ( τ ) ,
where
Φ ( τ , v ) = Φ a ( τ , v ) Φ f ( τ , v ) = 0 T k a ( t τ ) k a ( t v ) d t + 0 T k f ( t τ ) k f ( t v ) d t ,   t ( 0 , T ) ,
while x ( t ) is any signal at the input of the accelerometer, ξ is the Lagrange multiplier [5], k a ( t ) and k f ( t ) are the impulse responses of the accelerometer and the filter, respectively, and T denotes the time of accelerometer testing. The impulse responses k a ( t ) and k f ( t ) are determined as the inverse Laplace transforms of the corresponding transfer functions K a ( s ) and K f ( s ) . In turn, the transfer function K a ( s ) is obtained by way of the parametric identification of the accelerometer, while the transfer function K f ( s ) is convenient to be defined by the low-pass Butterworth filter given by,
K f ( s ) = S a i = 1 L ( s 2 π f c e j ( 2 i + L 1 ) π 2 L ) .
In Equation (6), L is the order of the filter, f c is the filter’s cut-off frequency, which is equal to the frequency bandwidth of the accelerometer, and j denotes the imaginary number.
The solution to Equation (4) determines the signal x m ( t ) that yields the maximum value of the coefficient ξ for which the ratio is satisfied,
ξ ( x , y ) = E ( y ) E ( x )
where E ( x ) and E ( y ) are the energies of any input and output signals x ( t ) and y ( t ) , respectively. The error y ( t ) is defined by,
y ( t ) = y a ( t ) y f ( t ) ,
where
y a ( t ) = 0 t k a ( t τ ) x ( τ ) d τ  
and
y f ( t ) = 0 t k f ( t τ ) x ( τ ) d τ  
are the output signals of the accelerometer and filter, respectively.
Both the signal x ( t ) as well as the signal x m ( t ) , maximizing the coefficient ξ , are limited in time T and magnitude a . Thus, we can write,
x ( t ) a , x m ( t ) a .
The denominator and nominator of Equation (7) are represented by,
E ( x ) = 0 T [ x ( t ) ] 2 d t ,
and
E ( y ) = E ( y a ) E ( y f ) = 0 T [ y a ( t ) ] 2 d t 0 T [ y f ( t ) ] 2 d t = 0 T [ y ( t ) ] 2 d t .
Substituting Equations (9) and (10) into Equation (13) and taking into account the relationship between the impulse responses,
k ( t ) = k a ( t ) k f ( t ) = 1 [ K a ( s ) ] 1 [ K f ( s ) ] ,
we finally have
E ( y ) = 0 T [ 0 t k ( t τ ) x ( τ ) d τ ] 2 d t .
When the signal x ( t ) is limited in magnitude M (i.e., it has a rectangular shape), then we can consider the condition given by,
a · sgn [ 0 T x ( v ) Φ ( τ , v ) d v ] = ξ x ( τ )
where sgn denotes the signum operation.
Hence, the signal maximizing the energy E ( y ) is represented by,
x m ( τ ) = a · sgn [ 0 T x ( v ) Φ ( τ , v ) d v ] ξ max
and ξ max is the maximum value of ξ ( x , y ) defined by Equation (7).
The procedure for determining the input signal that maximizes the energy E ( y ) is based on the FPA [5] and includes the following stages:
  • Calculate the function Φ ( τ , v ) based on Equation (5);
  • Determine the initial input signal x 0 ( τ ) , where ξ 0 is the initial Lagrange multiplier, assumed in advance by the formula,
    x 0 ( τ ) = a · sgn [ Φ ( τ , v ) ] ξ 0
  • Determine the maximizing signal on the basis of the iterative algorithm, given by,
    x i + 1 ( τ ) = a · sgn [ 0 T x i ( v ) Φ ( τ , v ) d v ] ξ i + 1 for i = 0 ,   1 ,   2 ,   ,   J 1 ,
    where J denotes the number of assumed iterations.
    It should be emphasized that for each iteration step i , the Lagrange multiplier ξ i + 1 must be determined based on Equation (7).
  • Calculate the energy E ( y ) by using Equation (15).
The signal x m ( t ) that yields the maximum value of the quotient ξ ( x , y ) is the one that maximizes the energy E ( y ) . This ME denoted by E max ( y ) is obtained by substituting the signal x m ( t ) into Equation (15). Then, we have
E max ( y ) = 0 T [ 0 t k ( t τ ) x m ( τ ) d τ ] 2 d t
The energy E max ( y ) is the accelerometer output energy, referred to as the output energy of filter, and is valid for the accelerometer frequency bandwidth.
Figure 2 shows the block diagram of the algorithm described above applied for determining the ME at the output of accelerometers. The switch S has the position 1 for i = 0 , and it has the position 2 when i > 0 . During the particular iterations i = 1, 2, …, J 1 of the algorithm above, the energies E ( y i ) and E ( y i 1 ) are compared with each other, and the value of the current ME and the corresponding input signal limited in magnitude are stored in the memory. This signal can be then processed by using specialized IT tools, e.g., neural network [30].
The final solutions of the ME and corresponding signal are delivered for i = J 1 . It should be noted that this algorithm converges quickly because the maximum number J of iterations does not exceed 30.
The signal producing the ME has the property that any other signal contained in its limitations can only yield less energy than the energy value caused by the maximizing signal [4]. Figure 3 shows examples of the maximizing signal x m ( t ) limited simultaneously in time and magnitude and three other signals x 1 ( t ) , x 2 ( t ) , and x 3 ( t ) contained in its limitations denoted by ± a (magnitude) and T (time).
The limited signal x m ( t ) in Figure 3 has seven time switches denoted by t 1 , t 2 , …, t 7 . The number of these switches and the times at which they occur are determined by the FPA, and they constitute the shape of the signal that maximizes the energy at the output of the accelerometer.

4. Monte Carlo-Based Procedure for Determining the ME

On the basis of obtaining both estimates S ˜ a , f ˜ 0 , and β ˜ of the accelerometer model and the associated uncertainties u ( S ˜ a ) , u ( f ˜ 0 ) , and u ( β ˜ ) , the values of the parameters denoted by S ˜ a max , f ˜ 0 max , and β ˜ max that produce the energy E max ( y ) should be determined. This can be carried out for the rectangular ranges shown in Figure 4 (the lines slanting to the right), by applying the MCP with the number of trials determined using Equation (3), where:
r min ( α ) = α u ( α ) r max ( α ) = α + u ( α ) r ( α ) = r max ( α ) r min ( α )
and α corresponds to the estimates S ˜ a , f ˜ 0 , and β , ˜ respectively, in the case of three ranges shown in Figure 4.
The trials of the MCP executed in three ranges were carried out by employing a uniform distribution. Hence, it is recommended to use the Wichmann–Hill pseudo-random number generator [28], according to the guidelines included in the JCGM supplement 1 [13]. The horizontal ranges of a uniform distribution in Figure 4 are constrained by the values of uncertainty associated with the parameter estimates of the accelerometer, while the vertical ranges of this distribution are contained in the range 0 ,   γ , where γ = 1 / r ( α ) is inversely proportional to the range of the parameter estimate variability.
Figure 5 shows the block diagram of the MCP for determining the E max ( y ) . This procedure includes three main sub-blocks:
  • Back-to-back procedure for identification and determination of the accelerometer parameters and associated uncertainties.
  • Selection of the special filter which is the reference for calculation of the energy at the output of accelerometer.
  • Execution of the FPA.
The input parameter for the MCP is the number of MCts assumed in advance according to Equation (3), while the input parameters for the FPA are: number of iterations J , value of time T , quantization step size Δ s of the calculation, values of the signal limit a , and initial Lagrange multiplier ξ 0 .

5. Results

The results of testing the accelerometer of the type DJB A/1800/V are presented below. Typical parameters for the accelerometer, taken from the corresponding datasheet [29], are as follows:
voltage sensitivity S a (±10%): 1.02 V/(m/s2),
frequency response (accelerometer bandwidth): 0.2 Hz–1 kHz ±5%,
maximum continuous g level: 4903 m/s2,
resonant frequency: 4 kHz.
The B&K8305 charge mode accelerometer was applied as the reference. The parameters of this accelerometer, contained in the datasheet [31], are:
charge sensitivity (±10%): 1.08 pC/g,
amplitude (±10%) and phase (±1°) range: 0.2 Hz–1 kHz.
The results of the back-to-back identification are obtained based on the measurement of the amplitude and phase responses carried out for 42 frequencies within the range 30 Hz–6 kHz. This range was selected in such a way as to exceed the value of the resonant frequency since it has a significant influence on the value of the estimates f ˜ 0 and β ˜ .
The measurement points of both frequency responses, obtained for the above frequencies (the first column denoted by F), are contained in the second and third columns (denoted by A and Φ ) of Table 1.
The measurement points above are registered by the NI-6221 data acquisition card with a sampling rate equal to 100 kS/s and the computer programs developed in the LabVIEW software.
The expanded uncertainties u(A) and u( Φ ) which are contained in the fourth and fifth columns of Table 1 were determined based on Table 1 included in the ISO 16063-21:2003(E) standard (7). These uncertainties are valid for the coverage factor k = 2 .
The accelerometer model responses (solid lines) obtained by utilizing the WLSP are shown in Figure 6. The covariance matrices involved in the basis kernel of the WLSP used for calculating the estimates of the accelerometer model parameters and applied for determining the uncertainties associated with these estimates are determined by using the MCP. The number M of MCts, which is equal to 2 × 10 5 , was taken as the minimum acceptable value obtained based on Equation (3). The coverage probability p equal to 0.95 was assumed in advance.
The expanded uncertainties, u(A) and u( Φ ), are marked in Figure 6 by the solid vertical lines within the measurement points of both responses. The graphs of both responses are presented by using a logarithmic scale due to the wide range of considered frequencies.
The estimates S ˜ a , f ˜ 0 , and β ˜ of the parameters of the accelerometer model and the associated uncertainties u ( S ˜ a ) , u ( f ˜ 0 ) , and u ( β ˜ ) determined by using the WLSP are tabulated in Table 2.
The model responses shown in Figure 6 are obtained by relative approximation by using the functions A ( f ) and φ ( f ) for the parameters tabulated in Table 2, where
A ( f ) = S ˜ a [ 1 ( f / f ˜ 0 ) 2 ] 2 + 4 β ˜ 2 ( f / f ˜ 0 ) 2
and
φ ( f ) = arctan 2 β ˜ ( f / f ˜ 0 ) 1 ( f / f ˜ 0 ) 2
The estimates and associated uncertainties from Table 2 were utilized to determine the parameters S ˜ a max , f ˜ 0 max , and β ˜ max giving the maximum value of the reference energy determined by the FPA controlled by the MCP. The input parameters for the FPA were selected as: J = 30, T = 2 ms, Δ s = 2 μs, and ξ 0 = 1, and the value of signal limit a = 1.037 V was assumed as equal to the estimate of the voltage sensitivity S ˜ a of the considered accelerometer. The time T was obtained from the steady state of the impulse response k ( t ) calculated based on Equation (14). This steady state was assumed to be equal to 5% of the maximum deviation (peak response) from zero of this impulse response. The tenth-order ( L = 10) analog Butterworth filter was adopted and is described by Equation (6). The filter cut-off frequency f c is equal to the frequency response of the considered accelerometer and is 1 kHz. The values of both the gain factor a and the signal limit a are equal to 1.03 V (i.e., the value of the estimate S ˜ a of the accelerometer voltage sensitivity). The upper limit of the filter order results from the computational limitation of the MathCad program which was applied for calculation of the impulse response in the way of determining the inverse Laplace transform on the basis of the relevant transfer function K f ( s ) .
Figure 7 shows the impulse responses k ( t ) , where the maximum deviations of these responses from the steady state value are equal to 1764. The number of MCts for the procedure controlling the FPA was 2 × 10 5 , similar to the WLSP. Based on parameters and uncertainties given in Table 2, the following values were obtained: u min ( S ˜ a ) = 1.034 V/(m/s2); u max ( S ˜ a ) = 1.040 V/(m/s2); u ( S ˜ a ) = 0.016 V/(m/s2); u min ( f ˜ 0 ) = 4067 Hz; u max ( f ˜ 0 ) = 4079 Hz; u ( f ˜ 0 ) = 12 Hz; u min ( β ˜ ) = 0.3178; u max ( β ˜ ) = 0.3186; and u ( β ˜ ) = 8 for the ranges shown in Figure 4.
The values of the parameters S ˜ a max = 1.038 V/(m/s2), f ˜ 0 max = 4070 Hz, and β ˜ max = 0.3178, which give the maximum reference energy E ( y ) max = 5.862 mV·m2, are obtained by applying the MCP for T = 2 ms, m = 127,548, and J = 15. The obtained uncertainty u ( E max ) of the MCP is equal to 0.008 mV·m2.
The relationship between the maximum reference energy E max ( y ) and the time T of the considered accelerometer is linear [5] for the steady state of the impulse response k ( t ) (i.e., for T > 2 ms). Thus, it is possible to determine the functional relationship between this energy and any time T higher than 2 ms. Hence, the value of energy E max ( y ) should be determined for any time higher than 2 ms, and then the relevant formula can be easily determined based on the linear regression. For example, if T = 10 ms, then we have E max ( y ) = 34.115 mV·m2 and u ( E max ) = 0.012 mV·m2, which are obtained for m = 93,515 and J = 9. Based on the two values of the energy E max ( y ) and the values of time T , we have:
E ( y , T ) max = 3531.0 · T 1.2 .
Figure 8 shows the signal x m ( t ) with two reversals producing the E max ( y ) at the output of the considered accelerometer obtained for T = 2 ms.
The reversals of the signal x m ( t ) are 0.68 ms and 1.49 ms.

6. Conclusions

The paper presents a proposal for the application of the FPA controlled by the MCP to determine the maximum reference energy at the output of accelerometers. Block diagrams together with a detailed description of the proposed solutions are presented.
The reference energy of the accelerometer should be understood as its output energy decreased by the output energy of the special filter (modeled by the low-pass type) with the cut-off frequency corresponding to the upper limit of the accelerometer frequency response. The lower limit of this response is neglected in this paper due to its low value (0.2 Hz).
Additionally, the MCP was used for parametric identification of both the parameters associated with the transfer function of the accelerometer and the associated uncertainties. Then, the value of the maximum reference energy was determined along with the associated uncertainties at the output of the accelerometer of the type DJB A/1800/V, considered as an example of the accelerometer. This was performed by checking the full range of the parameter variability of the accelerometer model by value-associated uncertainties. This is possible to perform only by using the MCP based on the pseudo-random number generator with a uniform distribution. This necessity results from the fact that the energy above can be obtained for the values of parameters contained within the variability ranges and not, for example, at their ends or in their middles. Only in this way it is possible to accurately determine the maximum reference energy, which can be the basis for mutual comparison of different types of accelerometers. The presentation of the relevant procedures together with a practical example giving the possibility of such comparison is the main contribution of this paper.

Funding

The presented research was conducted at the Faculty of Electrical and Computer Engineering, Cracow University of Technology, and was financially supported by the Ministry of Science and Higher Education, Republic of Poland (grant no. E-31/2020).

Conflicts of Interest

The authors declare no conflicts of interest.

Glossary

WLSPweighted least squares procedure
FPAfixed-point algorithm
MCPMonte Carlo procedure
MCSMonte Carlo simulation
MCtMonte Carlo trial
JCGMJoint Committee for Guides in Metrology
MEmaximum energy

References

  1. Yu, J.C.; Lan, C.B. System modeling and robust design of microaccelerometer using piezoelectric thin film. In Proceedings of the 1999 IEEE/SICE/RSJ. International Conference on Multisensor Fusion and Integration for Intelligent Systems, Taipei, Taiwan, 18 August 1999; pp. 99–104. [Google Scholar]
  2. Sharapov, V. Piezoceramic Sensors; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  3. Link, A.; Täbner, A.; Wabinski, W.; Bruns, T.; Elster, C. Modelling accelerometers for transient signals using calibration measurement upon sinusoidal excitation. Measurement 2007, 40, 928–935. [Google Scholar] [CrossRef]
  4. Layer, E. Theoretical principles for establishing a hierarchy of dynamic accuracy with the integral-square-error as an example. IEEE Trans. Instrum. Meas. 1997, 46, 1178–1182. [Google Scholar] [CrossRef]
  5. Honig, M.L.; Steiglitz, K. Maximizing the output energy of a linear channel with a time and amplitude limited input. IEEE Trans. Inform. Theory 1992, 38, 1041–1052. [Google Scholar] [CrossRef] [Green Version]
  6. Layer, E.; Tomczyk, K. Determination of Non-Standard Input Signal Maximizing the Absolute Error. Metrol. Meas. Syst. 2009, 17, 199–208. [Google Scholar]
  7. ISO 16063-21: 2003. Methods for the calibration of vibration and shock transducers—Part 21: Vibration calibration by comparison to a reference transducer. Available online: https://www.iso.org/standard/27053.html (accessed on 25 March 2020).
  8. Kollar, I. On frequency-domain identification of linear systems. IEEE Trans. Instrum. Meas. 1993, 42, 2–6. [Google Scholar] [CrossRef]
  9. Juang, J.N. Applied System Identification; Prentice-Hall: Upper Saddle River, NJ, USA, 1994. [Google Scholar]
  10. Guillaume, P. Frequency response measurements of multivariable systems using nonlinear averaging techniques. IEEE Trans. Instrum. Meas. 1998, 47, 796–800. [Google Scholar] [CrossRef]
  11. Zhao, S.; Wang, F.; Xu, H.; Zhu, J. Multi-frequency identification method in signal processing. Digit. Signal Process. 2009, 19, 555–566. [Google Scholar] [CrossRef]
  12. Isermann, R.; Münchhof, M. Identification of Dynamic Systems; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  13. BIPM, IEC, IFCC, ILAC, ISO, IUPAP, OIML, Guide to the expression of uncertainty in measurement. Suplement 1—Propagation of distributions using a Monte Carlo method. 2008. Available online: https://www.bipm.org/utils/common/documents/jcgm/JCGM_101_2008_E.pdf (accessed on 25 March 2020).
  14. BIPM, IEC, IFCC, ILAC, ISO, IUPAP, OIML, Guide to the expression of uncertainty in measurement. Supplement 2—Extension to any number of output quantities. 2011. Available online: https://www.bipm.org/utils/common/documents/jcgm/JCGM_102_2011_E.pdf (accessed on 25 March 2020).
  15. Guimarães Couto, P.R.; Carreteiro Damasceno, J.; Pinheiro de Oliveira, S. Monte Carlo simulations applied to uncertainty in measurement. In Theory and Applications of Monte Carlo Simulations; Chan, V., Ed.; IntechOpen: London, UK, 2013; pp. 27–51. [Google Scholar]
  16. Harris, P.M.; Cox, M.G. On a Monte Carlo method for measurement uncertainty evaluation and its implementation. Metrologia 2014, 51, 176–182. [Google Scholar] [CrossRef]
  17. Tokarska, M.; Gniotek, K. Assessment of measurement uncertainty using Monte Carlo method based on STATISTICA. Prz. Elektrotechniczny 2010, 86, 43–46. [Google Scholar]
  18. Palenčár, R.; Sopkuliak, P.; Palenčár, J.; Ďuriš, S.; Suroviak, E.; Halaj, M. Application of Monte Carlo method for evaluation of uncertainties of ITS-90 by standard platinum resistance thermometer. Meas. Sci. Rev. 2017, 17, 108–116. [Google Scholar] [CrossRef] [Green Version]
  19. Tomczyk, K. Levenberg-Marquardt algorithm for optimization of mathematical models according to minimax objective function of measurement systems. Metrol. Meas. Syst. 2009, 16, 599–606. [Google Scholar]
  20. Bimbi Junior, S.; de Toledo Fleury, A.; Ruas, R.; Chaves de Oliveira, V. Distortion reduction in FIR filters by approximation through window factor on function coefficients. Set IJBE 2016, 2, 13–22. [Google Scholar] [CrossRef]
  21. Wyner, A. Spectra of bounded functions. In Open Problems in Communication and Computation; Cover, T.M., Gopinath, B., Eds.; Springer: New York, NY, USA, 1987; pp. 46–48. [Google Scholar]
  22. Elia, M.; Taricco, G.; Viterbo, E. Optimal energy transfer in band-limited communication channels. IEEE Trans. Inf. Theory 1999, 45, 2020–2029. [Google Scholar] [CrossRef] [Green Version]
  23. Long, G.; Nelakanti, G. Iteration methods for Fredholm integral equations of the second kind. Comput. Math. Appl. 2006, 53, 886–894. [Google Scholar] [CrossRef] [Green Version]
  24. Tomczyk, K. Special signals in the calibration of systems for measuring dynamic quantities. Measurement 2014, 49, 148–152. [Google Scholar] [CrossRef]
  25. Rutland, N.K. The principle of matching: Practical conditions for systems with inputs restricted in magnitude and rate of change. IEEE Trans. Autom. Control 1994, 39, 550–553. [Google Scholar] [CrossRef]
  26. Tomczyk, K. Influence of Monte Carlo generations applied for modelling of measuring instruments on maximum distance error. IEEE Trans. Instrum. Meas. Control 2019, 41, 74–84. [Google Scholar] [CrossRef]
  27. Tomczyk, K. Impact of uncertainties in accelerometer modeling on the maximum values of absolute dynamic error. Measurement 2016, 80, 71–78. [Google Scholar] [CrossRef]
  28. Wichmann, B.A.; Hill, I.D. Generating good pseudo-random numbers. Comput. Stat. Data Anal. 2006, 51, 1614–1622. [Google Scholar] [CrossRef]
  29. DJB Instruments, A/1800/T, A/1800/TC, A/1800/VTC Micro g Piezoelectric IEPE Accelerometer. 2016. Available online: http://www.wibroakustyka.com.pl/dokumenty/A1800_june_2016.pdf (accessed on 25 March 2020).
  30. Dudzik, M.; Jagiello, A.; Drapik, S.; Prusak, J. The selected real tramway substation overload analysis using the optimal structure of an artificial neural network. In Proceedings of the 2018 International Symposium on Power Electronics, Electrical Drives, Automation and Motion (SPEEDAM), Amalfi, Italy, 20–22 June 2018. [Google Scholar]
  31. Brüel & Kjær, Product Data: Reference Standard Accelerometer Types 8305 and 8305-001. 2018. Available online: https://www.bksv.com/-/media/literature/Product-Data/bp2052.ashx (accessed on 25 March 2020).
Figure 1. Block diagram of system for back-to-back identification of accelerometer.
Figure 1. Block diagram of system for back-to-back identification of accelerometer.
Energies 13 01552 g001
Figure 2. Block diagram of algorithm intended for determining the maximum energy (ME).
Figure 2. Block diagram of algorithm intended for determining the maximum energy (ME).
Energies 13 01552 g002
Figure 3. Examples of the limited signal and other signals contained in its limitation.
Figure 3. Examples of the limited signal and other signals contained in its limitation.
Energies 13 01552 g003
Figure 4. Range of the Monte Carlo procedure (MCP) for the estimates associated with the accelerometer model.
Figure 4. Range of the Monte Carlo procedure (MCP) for the estimates associated with the accelerometer model.
Energies 13 01552 g004
Figure 5. Block diagram of Monte Carlo-based procedure for determining the ME at the output of the accelerometer.
Figure 5. Block diagram of Monte Carlo-based procedure for determining the ME at the output of the accelerometer.
Energies 13 01552 g005
Figure 6. Amplitude and phase responses with associated uncertainties.
Figure 6. Amplitude and phase responses with associated uncertainties.
Energies 13 01552 g006
Figure 7. Impulse responses k(t) calculated based on Equation (14).
Figure 7. Impulse responses k(t) calculated based on Equation (14).
Energies 13 01552 g007
Figure 8. Signal x m ( t ) obtained for the considered accelerometer.
Figure 8. Signal x m ( t ) obtained for the considered accelerometer.
Energies 13 01552 g008
Table 1. Measurement points of the frequency responses together with expanded uncertainties (k = 2) for the accelerometer DJB A/1800/V.
Table 1. Measurement points of the frequency responses together with expanded uncertainties (k = 2) for the accelerometer DJB A/1800/V.
F
[Hz]
A
[V/(m/s2)]
Φ [deg.]u(A)
[V/(m/s2)]
u ( Φ )
[deg.]
F
[Hz]
A
[V/(m/s2)]
Φ
[deg.]
u(A)
[V/(m/s2)]
u ( Φ )
[deg.]
301.03−0.10.011.020001.26−20.40.031.0
401.04−0.40.011.022001.25−19.50.041.0
501.00−0.30.011.024001.38−30.00.041.0
601.02−0.50.011.026001.52−40.50.051.0
701.03−0.40.011.028001.66−53.70.051.0
800.98−0.70.011.030001.65−38.10.051.0
901.04−0.70.011.032001.77−41.60.051.0
1000.96−0.90.011.034001.93−53.80.061.0
2000.95−2.30.011.036002.12−58.20.061.0
3001.05−2.00.011.038002.05−73.10.061.0
4001.01−3.20.011.040002.20−93.20.071.0
5000.96−2.40.011.042001.81−108.90.051.0
6000.98−6.30.011.044001.79−99.90.051.0
7001.03−3.00.011.046001.46−126.20.041.0
8001.04−6.00.011.048001.40−111.20.041.0
9000.97−2.40.011.050001.25−117.60.041.0
10001.05−3.50.011.052001.07−138.00.031.0
12001.16−14.90.021.054000.94−140.00.031.0
14001.06−18.70.021.056000.89−143.00.031.0
16001.15−9.60.021.058000.81−154.40.021.0
18001.27−8.20.031.060000.96−158.50.031.0
Table 2. Estimated parameters and associated uncertainties obtained by using the weighted least squares procedure (WLSP).
Table 2. Estimated parameters and associated uncertainties obtained by using the weighted least squares procedure (WLSP).
S ˜ a
[V/(m/s2)]
f ˜ 0
[Hz]
β ˜
[–]
u ( S ˜ a )
[V/(m/s2)]
u ( f ˜ 0 )
[Hz]
u ( β ˜ )
[–]
1.037 4073 0.3182 0.00360.0004

Share and Cite

MDPI and ACS Style

Tomczyk, K. Monte Carlo-Based Procedure for Determining the Maximum Energy at the Output of Accelerometers. Energies 2020, 13, 1552. https://doi.org/10.3390/en13071552

AMA Style

Tomczyk K. Monte Carlo-Based Procedure for Determining the Maximum Energy at the Output of Accelerometers. Energies. 2020; 13(7):1552. https://doi.org/10.3390/en13071552

Chicago/Turabian Style

Tomczyk, Krzysztof. 2020. "Monte Carlo-Based Procedure for Determining the Maximum Energy at the Output of Accelerometers" Energies 13, no. 7: 1552. https://doi.org/10.3390/en13071552

APA Style

Tomczyk, K. (2020). Monte Carlo-Based Procedure for Determining the Maximum Energy at the Output of Accelerometers. Energies, 13(7), 1552. https://doi.org/10.3390/en13071552

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