Next Article in Journal
Multi-Excitation Infrared Fusion for Impact Evaluation of Aluminium-BFRP/GFRP Hybrid Composites
Next Article in Special Issue
Spin-Orbital Conversion of a Strongly Focused Light Wave with High-Order Cylindrical–Circular Polarization
Previous Article in Journal
Stopping the Mobile Robotic Vehicle at a Defined Distance from the Obstacle by Means of an Infrared Distance Sensor
Previous Article in Special Issue
Minimal Focal Spot Size Measured Based on Intensity and Power Flow
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Peculiarities and Applications of Stochastic Processes with Fractal Properties

by
Oleg Semenovich Amosov
1,* and
Svetlana Gennadievna Amosova
2
1
Laboratory of Intellectual Control Systems and Modeling, V.A. Trapeznikov Institute of Control Sciences of Russian Academy of Sciences, 65 Profsoyuznaya Street, 117997 Moscow, Russia
2
Laboratory of Cyber-Physical Systems, V.A. Trapeznikov Institute of Control Sciences of Russian Academy of Sciences, 65 Profsoyuznaya Street, 117997 Moscow, Russia
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(17), 5960; https://doi.org/10.3390/s21175960
Submission received: 31 July 2021 / Revised: 23 August 2021 / Accepted: 30 August 2021 / Published: 5 September 2021
(This article belongs to the Special Issue State-of-the-Art Optical Sensors Technology in Russia 2021-2022)

Abstract

:
In this paper, the fractal properties of stochastic processes and objects in different areas were specified and investigated. These included: measuring systems and sensors, navigation and motion controls, telecommunication systems and networks, and flaw detection technologies. Additional options that occur through the use of fractality were also indicated and exemplified for each application. Regarding the problems associated with navigation information processing, the following fractal nature processes were identified: errors of inertial sensors based on the microelectromechanical systems called MEMS, in particular gyroscopic drift and accelerometer bias, and; the trajectory movement of mobile objects. With regard to navigation problems specifically, the estimation problem statement and its solution are given by way of the Bayesian approach for processing fractal processes. The modified index of self-similarity for telecommunication series was proposed, and the self-similarity of network traffic based on the R/S method and wavelet analysis was identified. In failure detection, fractality manifested as porosity, wrinkles, surface fractures, and ultrasonic echo signals measured using non-destructive sensors used for rivet compound testing.

1. Introduction

Many processes in natural and artificial systems are known to have fractal structures [1,2]. Currently, the fractal properties of objects and processes continue to be studied in many areas [3]. With the development of higher-precision sensors and improved measuring and diagnostic systems, a need to improve our understanding of fractal structures and system errors has arisen. One goal of contemporary research in this area is to improve the cybersecurity of modern systems. The use of fractality allows mathematical models of processes and objects to be refined, thereby helping to overcome problems of analysis and synthesis to achieve greater competence.
This article examines the features of stochastic process models with fractal properties and their applications in some important areas. These areas include: sensors and measuring systems [4,5], navigation and motion controls [6,7], telecommunication systems [3,8,9], and flaw detection technologies [10,11].
The structure of the paper is as follows: processes and objects of a fractal nature in key areas. The possibility of using fractal properties to improve the quality of problem solving is then by evaluated. Following this, the positive effects of fractality in these contexts is considered.

2. Fractal Process Model

First-order Markov processes have been widely used in practice. However, there are others random processes with short- and long-term memory.
One possible model by which to understand the higher-order Markov process as it applies to memory, is a model of fractional Brownian motion [12,13] (fractional Wiener process—FWP) containing the Hurst parameter.
To identify the features of the fractional structure of the process, fractality indicators such as the Hurst parameter are used. Several methods are used to estimate the Hurst parameter, including: rescaled range analysis (the R/S method), variance change graph analysis, and wavelet analysis [13,14]. According to the R/S method, the Hurst parameter of a time series is calculated by the following formula:
H = ln ( R S ) / ln ( n 2 ) ,
where R is the range of the first n cumulative deviations from the mean, S is the standard deviation of observations, and n is the number of observations.
The FWP is called a random process x ( t ) with Hurst parameter H ( 0 < H < 1 ) if Δ x = x ( t ) x ( t 0 ) has a Gaussian distribution with zero mathematical expectation and variance σ 2 ( t t 0 ) 2 H , where σ is a positive constant, namely its probability measure [1].
P ( Δ x < X ) = 1 2 π σ ( t - t 0 ) H X exp [ 1 2 ( u σ ( t t 0 ) H ) 2 ] d u
at H = 0.5, the FWP is consistent with classical Brownian movement (the classical Wiener process—CWP). The increments Δ x = x ( t ) x ( t 0 ) of fractional Brownian motion are called fractional Gaussian noise, with variance obeying the ratio D [ x ( t ) - x ( t 0 ) ] = σ 2 ( t - t 0 ) 2 H .
At Hurst parameter H: (1) if 0 < H < 0 . 5 , then the process has anti-persistent properties, when an upward trend is replaced by a downward trend or vice versa; (2) if H 0.5 , then the process is consistent with the classical Brownian movement; (3) if 0.5 < H < 1 , then the process has persistent (i.e., trend-resistant) properties.

3. Mathematical Model of Inertial Sensor Errors in Navigation Systems with Respect to the Chaotic Indicator

Due to the development of high-performance navigation and motion control systems, the requirements for estimating the accuracy of true process parameters that determine the qualitative characteristics of these systems have increased. When processing navigation information parameters by way of an inertial navigation system, the requirements for precision features can be met by making system error descriptions precise [4,6].
When processing navigation information, it is important to take into account the characteristics of errors in the output signals of navigation system sensors. The Wiener process is widely used to describe modeling errors of inertial sensors with stochastic features [15]. In particular, such a model is used to describe the zero drift of a gyroscope.
The Wiener process model [16] for an error x ( t ) in continuous time is
x ˙ ( t ) = w ( t ) ,
where
E [ w ( t ) ] = 0 , E [ w ( t ) w ( τ ) ] = q ( t ) δ ( t τ ) ,
where q ( t ) is the intensity of continuous white noise w ( t ) .
The discrete-time state equation with a sampling period T = t i t i 1
x i = x i 1 + w i
The covariance Q for discrete noise w i in case q ( t ) = q = c o n s t
Q = σ 2 = q T .
The process w i can be represented by a model using white noise n i with unit intensity via ratio w i = q T n i .
It should be noted that the process x ( t ) described by a stochastic differential Equation (3) is a Wiener process defined through the white noise w ( t ) . A sequence x i (5) is a sequence with uncorrelated increments and is called a Wiener sequence.
Based on the results of other authors’ experimental analyses, error behavior graphs of inertial sensor outputs suggest that error description processes have a fractal nature.
This is indicated by the results of gyroscopic drift testing presented in [17]. Constant and variable components, correlation coefficients, variance within a drift variable, and white noise were taken into account in describing the gyroscopic drift model [17]. The paper [15] also proposes a mathematical model of sensor error based on Allan variation, which uses unknown parameters. The model characterizes the contribution of the error components of the sensor, corresponding to a linear trend, Wiener process, flicker noise, white noise and quantization noise.
Recently, interest in the class of non-stationary processes that have a fractal nature has deepened [1,2,3,5,6,7,8,9,10,14,18,19].
Models of sensor errors based on fractal Wiener process
Choosing a probability measure (2) for increments of a random process x ˙ ( t ) = w ( t ) , x i = x i 1 + w i , Δ x = x ( t ) x ( t 0 ) = x i x i 1 , we obtain a model with the FWP for an error.
Sensor Error Measurement Model
Based on the obtained results of analyzing gyroscope and accelerometer tests, a new mathematical model of inertial sensor error based on the fractal Wiener process is proposed considering given the chaotic indicator
y i = η i + x i + v i ,
where η i is the trend component, x i is the component expressing entropy of a process, modeled by the FWP in accordance with (2) and dependent on H—the Hurst parameter; v i is the random white noise with variance R i .

3.1. Experimental Confirmation of Fractal Properties of Gyroscope Zero Drift Process

The digital gyroscope STMicroelectronics L3G4200D PmodGYRO (Figure 1a) was examined using a myRIO National Instruments (NI) device utilizing its onboard FPGA, a dual-core programmable controller, and LabVIEW software environment [20]. The L3G4200D is a low-power, three-axis, angular rate sensor able to provide unprecedented, zero bias stability (zero-rate level), and sensitivity stability over temperature and time [21].
Measured range: one of the available value ranges, e.g., ±250 dps, or °/s (grad/s). The sensitivity value is 0.00875°/LSB (LSB—least significant bit).
Figure 2a shows the gyroscope drift at the following parameters: data transfer rate is 100 Hz; sampling time is 1 h with sample spacing T = 10 2 s—360,000 values at a sensitivity of 250 dps. The zero bias (digital zero-rate level) is more than 0.5 °/s.
In our experiment (Figure 2a), a signal is shown at the output of the electronic part of the sensor where the noise component (bias instability) is caused by electronic components that are part of an electrical circuit.
The value of the Hurst parameter H was calculated by the method based on wavelet analysis. We used the standard MatLab wfbmesti function in accordance with previous approaches [22,23,24]. As can be seen from Figure 2b, the value of the H parameter varies by about 0.28, i.e., the process has anti-persistent properties.

3.2. Modeling of Accelerometer Bias

The digital accelerometer ADXL345 [25] (Figure 1b) is examined using the myRIO NI device and the LabVIEW software environment. The accelerometer was placed on a fixed insulated base. The output signals were registered along one of the sensitivity axes, orthogonal to the gravitational vector g.
Figure 3a shows the accelerometer bias at the following parameters: 16 g resolution range; 0.0625 g/LSB scale; 100 Hz data transmission rate; 20 min sampling time with T = 10 2 s—127,077 sample spacing values; and the temperature of 23 °C, L S B = 2 g/210.
Figure 3a shows the signal at the output of the accelerometer, where the noise component of bias instability is in the range of −1 to +1 LSB, which is caused by the electronic part of the sensor [25].
When examining the accelerometer bias for 20 min (Figure 3b), the Hurst parameter fluctuates by about 0.03. Therefore, the process has anti-persistent properties.
In inertial sensors, changes in the Hurst parameter from small values to higher ones can serve as an indicator that errors of inertial systems accumulate due to the departure of the zero bias because of technical malfunctions, temperature influence, and vibration. For different types of inertial sensors (gyroscopes, accelerometers), the Hurst parameter will be different.

4. Navigation and Motion Control

The accuracy of parameter estimation in orientation tracking is mostly determined by the choice of mathematical, moving object, descriptions, and the choice of efficient algorithms for nonlinear filtering. In both cases, when solving the problem of trajectory tracking, the mathematical models based on the classical Wiener process are widely applicable. At the same time, an aggregation in the form of the FWP is used to model stochastic processes with fractal properties [1,18].
Mathematical Models of Object Movement for Trajectory Tracking Problems
The object moves in a horizontal plane with coordinates x and y of Cartesian axes. Let us consider the motion model for one coordinate x; the model for another coordinate, y, is similar.
Let us also consider the object movement models traditionally used to solve orientation tracking problems [16,18,26,27].
Model with Wiener process for acceleration
x = w ( t ) ,
with white noise characteristics w ( t ) .
Model with white noise for acceleration, or Wiener process for speed
x ¨ ( t ) = w ( t ) .
A model with white noise for speed, or a Wiener process for a coordinate, is determined by Equations (3) and (4).
It is essential that all presented models of motion, (3), (8), and (9), are defined through the Wiener process, and that their discrete representations with a period of discretization T = t i t i 1 in the form of x i = Φ i x i 1 + w i , are defined through the Wiener sequence (5). Here, w i is an n-dimensional-centered Gaussian white noise sequence with the covariance matrix Q i . Φ i is the known matrix of a respective dimension.
With this in mind, we propose new models of object motion based on the use of the fractional Wiener process (2) [18,28].
Models of object movement based on the FWP
By choosing the probability measure (2) for increments of a random process
x ˙ ( t ) = w ( t ) , x i = x i 1 + w i , Δ x = x ( t ) x ( t 0 ) = x i x i 1 ,
we obtain a model with the FWP for the coordinate. Proceeding similarly for processes (8) and (9), we obtain a model with the FWP for acceleration and a model with the FWP for speed, respectively.
Measurement model for trajectory tracking
A range and azimuth tracking approach can be used to determine the location of an object.
Measuring range ( ρ i ) and bearing ( α i )
ξ i = [ ξ 1 i ξ 2 i ] = [ ρ i α i ] = [ s 1 ( x i , y i ) s 2 ( x i , y i ) ] + [ Δ ρ i Δ α i ] = [ x i 2 + y i 2 a r c t g ( y i / x i ) ] + [ Δ ρ i Δ α i ] .
To obtain a linear model of measurements, the primary measurements of the polar coordinates of an object, range ( ρ i ), and bearing ( α i ), can be converted into a rectangular coordinate system and are presented in the following form [18,26]
x i * = ρ i cos α i x i + Δ x i ;   y i * = ρ i sin α i y i + Δ x i ,
where Δ x i and Δ y i are random errors of coordinate measurements with the following characteristics
E ( Δ x i ) = E ( Δ y i ) = 0 ; σ Δ x i 2 = E ( Δ x i 2 ) = σ Δ ρ i 2 cos 2 α i + σ Δ α i 2 ρ i 2 sin 2 α i ; σ Δ y i 2 = E ( Δ y i 2 ) = σ Δ ρ i 2 sin 2 α i + σ Δ α i 2 ρ i 2 cos 2 α i ; E ( Δ x i Δ y i ) = R Δ x i Δ y i = ( σ Δ ρ i 2 σ Δ α i 2 ρ i 2 ) sin α i cos α i
Linear measurements will be used for the Kalman Filter (KF). Let us now consider an example of fractal process filtering.

5. Fractal Process Filtering

To clearly identify the features of filtration in the fractal process, we intentionally choose the simplest model to express the FWP. It is necessary to evaluate the coordinate of a fractal Wiener process
x i = x i 1 + w i
using linear measures of type y i = x i + v i , where i = 0 , 1 , are the moments of time and w i , v i are the zero-mean Gaussian white noises independent of each other and of x i , with variances Q = σ w 2 and R = σ v 2 , respectively. Here, the values of the process x 0 at the initial represent the random zero-mean Gaussian variable with covariance at the P x 0 = σ x 0 2 .
Modeling parameters: initial state m x 0 = 0 , σ x 0 = 20 ; σ v = 10 , discretization period T = 1   s . The FWP (11) was modeled using the wfbm function in the MatLab mathematical package, which is based on wavelet decomposition and simulates a fractal Wiener sequence with a given Hurst parameter.
Kalman filtering and a non-recurrent Bayesian algorithm (BA) of optimal linear filtering were used to obtain the required estimate [28].
It is essential that if, for the classical Wiener process, increments Wi are independent, then the FWP does not have this property and is not a Markov process, except in cases where H = 0.5 when the FWP is consistent with the CWP. In such cases, KF only gives the optimal solution for this value of H parameter and is not optimal with other H parameter values.
In general, when 0 < H < 1 , optimal estimations of non-recurrent BA and expected (analysis) covariance matrices for estimation error e i = x i x ˜ i ( y i ) , are determined by means of the following expressions [28,29]:
x ˜ i = x - i + P x i y i ( P y i y i ) 1 [ y i y - i ] , P e i = P x i x i P x i y i P y i y i 1 P y i x i
where P a i b i are covariance matrices.
The selective (valid) r.m.s. of the estimation errors are presented, calculated as
σ ˜ i η 1 L j = 1 L ( e η i ( j ) ) 2 , e η i ( j ) = x ˜ i η ( j ) ( y i ( j ) ) x i ( j ) , η = B A , K F , L = 3000 .
Modeling and estimation of the FWP process was carried out with different Hurst parameters: H: 0.1; 0.5; 0.9 (Figure 4). The number of measurements can be expressed as i = 1.200 ¯ .
As Figure 4a–c show σ i F K is the expected r.m.s of estimation errors corresponding to the error variance of optimal estimation, which is the only element of the analysis matrix of KF covariance errors. σ ˜ i F K is the selective r.m.s of filtering errors when using KF; σ i B A is the expected r.m.s. of filtering errors when using non-recurrent BA; σ ˜ i B A is the selective r.m.s. of filtering errors for BA, and σ ˜ i is the selective r.m.s. of measurement errors.
The results show the coincidence of the expected σ i B A (12) and the valid σ ˜ i B A (13) accuracy characteristics of the presented optimal linear algorithm for non-recurrent estimation. When parameter H = 0.1, higher estimation accuracy can be achieved with the help of the optimal linear algorithm than with KF. When the parameter H = 0.5, both KF and the linear optimal algorithm estimate the random process with equal accuracy. When parameter H = 0.9, for KF, estimation error is increased in comparison to the optimal linear algorithm. A mismatch between the results of analysis and the actual (selective) characteristic for KF is observed at H = 0.1 and H = 0.9, respectively.
The study shows that a higher estimation accuracy can be achieved using non-recurrent algorithms, unlike KF.
In navigation and motion control, the values of the Hurst parameter for object motion models may change due to changes in the state of the environment (e.g., turbulence, weather conditions, natural phenomena, etc.).
For Kalman-type algorithms to be efficient, it is necessary to solve the adaptive problem of estimation by introducing and evaluating an unknown parameter—the fractal parameter H.

6. Information and Telecommunication System

The development and research of mathematical models of network telecommunication traffic allow us to define the boundaries of network performance and optimize the network structure of the system to minimize queueing delays. Traditional network traffic models show pulsations on short time scales but very smooth traffic at larger scales. Regardless, traffic shows variability over a wide range of time scales. When studying network telecommunication traffic, one of the important tasks is to identify anomalies and respond to them early.
Analyses of network telecommunication traffic result in time series (TS) analyses most widely used in traditional methods of statistically analyzing random variable functions. Alongside with traditional methods the signal processing methods based on fractal and wavelet-transformations [5,7,9,10,14,18,19], have seen more widespread use in recent years. The first of these analyses is specifically concerned with the identification of fractal properties that allow time series’ to be attributed to predetermined models. One distinct benefit of these analyses is their ability to reveal features of local structures and the different properties of complex signals that would otherwise be invisible using standard, real-time representations.
Telecommunication network traffic exhibits self-similar properties over a wide range of time scales that are very different to the properties of traditional models. Self-similar models are appropriate for telecommunication network traffic, as they have the capacity to estimate network performance/quality, allocate resources, and ensure quality of service [9].
Therefore, there is a need to develop new algorithms and modify those that exist for the purpose of analyzing the information systems (IS) of telecommunication series.
In order to analyze the telecommunication system, the use of an algorithm consisting of the following steps was proposed:
(1)
Pre-processing of the traffic time series consisting of statistical data sampling with the purpose of generating the TS we are interested in;
(2)
Estimation of fractality indicators using different methods such as R/S methods and wavelet analyses;
(3)
Identification of the mathematical models of TS by approximation to known models, or through synthesizing the structures and parameters of mathematical TS models.
To reduce the possibility of misinterpreting results of the parameter H evaluation, we propose to analyze not only the entire sample, but also the individual blocks of the sample. In this context, we propose to consider the blocks separately with a window offset and/or in a progressive way. Depending on the parameter H estimation method, we will make further calculations for each concerned block(s) separately to determine the expected value of the results found. This modified Hurst parameter, in the author’s opinion, will more adequately reflect the real environmental issues associated with identifying the fractal properties of telecommunication time series.
The result of the analyses proposed here will provide full objective information about the telecommunication time series of IS, with the determined characteristics allowing for the identification of a mathematical TS model. To solve some practical problems, for example, the identification of anomalies in the network traffic that may impact IS data integrity, it is necessary to compare the TS with an identified model and compare the results.
In order to analyze telecommunication traffic, we will consider the time series of an organization’s incoming and outgoing internet traffic (Figure 5). The organization selected will have an average traffic volume of more than 160 Gbytes per day.
Let us consider two time series of incoming and outgoing internet traffic for an organization with an average volume of more than 160 Gbyte a day. Let the first TS be collected over a sampling period of 60 days with sample spacing equal to 1 h, and let the second TS be sampled over 7 days, with sample spacing equal to 1 min. The parameter H can be obtained by charting the dependence of lg ( E ( R ( n ) / S ( n ) ) ) on lg ( n ) and then using the obtained points to match a straight line with a slope of parameter H by the least squares method. Here, we will henceforth consider E the symbol of the expected value.
Figure 6 shows the dependence of the standardized range logarithms on the logarithms of range lengths, and the linear regression with a slope of parameter H according to the R/S-method.
The time series were divided into blocks in a progressive way, with spacing of 100 discrete points, i.e., the first block included 100 values, the second included 200, the third included 300, etc. Each block, in turn, was divided into intervals ranging in size from 10 discrete points to a 20% block length increase in size between adjacent intervals. For each interval, the R/S ratio and the R/S mathematical expectation of the given block were calculated.
While investigating the information and telecommunication systems, the parameter H ≥ 0.8. Therefore, this process has persistent properties. Because of this, we propose to use fractality indicators for the time series model description, along with various traffic indicators (byte intensity per time unit, batches intensity per time unit, average batch size, throughput, loss rate, channel delay, loading, etc.).
In addition to the network traffic model, the fractality indicators can also be used to determine threats to information resources by solving a network classification problem [19]. Thus, the paper [19] proposes the modeling and implementation of a threat protection system with an attack forecast model mode of adaptation and parameter setting based on intelligent technologies in addition to fractal and wavelet analyses.
Thus, in telecommunication systems, the Hurst parameter can serve as an indicator when abnormal influences occur in the system (unauthorized access, computer attacks, etc.).

7. Defect Detection

Today, developments relating to the transition to advanced digital intelligent manufacturing technologies, robotic systems, big data processing systems using artificial intelligence, and machine learning continue to actively develop. In every industry, the use of non-destructive testing methods is able to investigate components’ properties on site or in the field. Non-destructive testing methods appeal because they are much cheaper and faster than traditional, destructive, tests and do not cause any damage to the constructions’ components.
Methods of image recognition using photo- and camcorder-sourced images are used to detect defects [30,31]. Another effective approach is an ultrasonic method of non-destructive testing [32,33], which can be operated using flaw detectors.
Fractality manifests as porosity, wrinkles, and surface fractures, all of which can be detected during organoleptic testing and by ultrasonic echo signals by non-destructive control of rivet joints.
The analysis of non-destructive testing signals was carried out to identify their fractal properties.
Figure 7 represents the observed echo signals during non-destructive testing using an ultrasonic flaw detector (Olympus EPOCH LTC) with an accuracy of measuring rivet length up to 0.01 mm. The accuracy of measuring the echo signal amplitude is up to 0.25%.
The red line is a strobe for detecting defects. Rivet defects differ in the type of echo signal in relation to this strobe.
The Hurst parameter was calculated for 28 implementations of echo signals of an ultrasonic flaw detector for each of the normal groups, and for three groups with defective rivet types. For normal rivet joints, the parameter H was approximately 0.2, whereas for all the defective rivets it was about 0.35. This can serve as a reliable indicator of the presence or absence of a defect in binary classification. As we have shown in [34] in the case of binary classification, 100% accuracy can be achieved using both deep neural networks with an LSTM layer, and by estimating the Hurst parameter H for the same source data.

8. Conclusions

Following the specification and investigation of the fractal properties of processes and objects in different areas, the present study showed that when constructing models of processes and objects in various domains, proper allowances must be made for fractality.
A stochastic modeling process based on the fractal Wiener process was proposed, with provision for the Hurst parameter.
Relevant to the problems of navigation information processing, the following fractal nature processes were identified: errors of inertial sensors based on the microelectromechanical systems called MEMS, in particular the gyroscopic drift and accelerometer biases, and the trajectory movement of mobile objects.
New mathematical models of inertial sensor errors based on the fractal Wiener process were presented.
New models of object motion based on the use of the fractional Wiener process were also proposed. The specific nature of estimating fractal process parameters using Kalman-type and non-recurrent Bayesian algorithms were investigated. Our illustrative examples show that, when estimating the state of the studied processes, fractality must be accounted for.
In addition, we revealed that the traffic in information systems has a fractal nature. A modified index of self-similarity for telecommunication series was proposed, and we identified network traffic self-similarity based on the R/S method and wavelet analysis. At the same time, we proposed to use the fractality indicator when identifying threats to information resource protection systems.
In failure detection, fractality is manifested as porosity, wrinkles, surface fractures, and ultrasonic echo signals through the non-destructive testing of rivet compounds. We recommend using the Hurst parameter when seeking the reliable binary classification of defects in rivet joints.

Author Contributions

Conceptualization, O.S.A. and S.G.A.; methodology, O.S.A.; software, S.G.A.; validation, O.S.A. amd S.G.A.; formal analysis, O.S.A. and S.G.A.; investigation, O.S.A.; resources, S.G.A.; data curation, O.S.A.; writing—original draft preparation, O.S.A. and S.G.A.; visualization, S.G.A.; supervision, O.S.A. and S.G.A.; project administration, O.S.A.; funding acquisition, O.S.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

A study involving humans or animals is absent.

Informed Consent Statement

A study involving humans is absent.

Data Availability Statement

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Crownover, R.M. Introduction to Fractals and Chaos; Jones and Bartlett Publishers, Inc.: Boston, MA, USA; London, UK, 1995. [Google Scholar]
  2. Aguirre, J.; Viana, R.L.; Sanjuan, M.A.F. Fractal structures in nonlinear dynamics. Rev. Mod. Phys. 2009, 81, 333–386. [Google Scholar] [CrossRef]
  3. Nualart, D. Fractional Brownian motion: Stochastic calculus and applications. In Proceedings of the International Congress of Mathematicians, Madrid, Spain, 22–30 August 2006; pp. 1541–1562. [Google Scholar]
  4. Stepanov, O.A.; Motorin, A.V. Performance Criteria for the Identification of Inertial Sensor Error Models. Sensors 2019, 19, 1997. [Google Scholar] [CrossRef] [Green Version]
  5. Yang, X.B.; Jin, X.Q.; Du, Z.M.; Zhu, Y.H. A novel model-based fault detection method for temperature sensor using fractal correlation dimension. Build. Environ. 2011, 46, 970–979. [Google Scholar] [CrossRef]
  6. Yang, L.; Sun, Z.L.; Xu, C.; Chai, X.D. The Application of Software Reliability Based on Fractal Theory in Integrated Navigation System. In Proceedings of the 2013 International Forum on Special Equipments and Engineering Mechanics, Nanjing, China, 27–28 July 2016; pp. 118–122. [Google Scholar]
  7. Han, T.; Yang, Y.X.; Huang, G.W. Vertical scaling optimization algorithm for a fractal interpolation model of time offsets prediction in navigation systems. Appl. Math. Model. 2021, 90, 862–874. [Google Scholar] [CrossRef]
  8. Liu, J. Fractal Network Traffic Analysis with Applications. Ph.D. Thesis, Georgia Institute of Technology, Atlanta, GA, USA, 2006; 115p. [Google Scholar]
  9. Radev, D.; Lokshina, I.; Radeva, S. Modeling and Simulation of Self-Similar Traffic for Wireless IP Networks. In Proceedings of the 2009 Wireless Telecommunications Symposium, Prague, Czech Republic, 22–24 April 2009. [Google Scholar] [CrossRef]
  10. Xianglin, Z.; Jin, S.J. A Fractal-Based Flaw Feature Extraction Method for Ultrasonic Phased Array Nondestructive Testing. In Proceedings of the International Conference on Mechatronics and Automation, Changchun, China, 9–12 August 2009. [Google Scholar] [CrossRef]
  11. Jovančević, I.; Pham, H.H.; Orteu, J.J.; Gilblas, R.; Harvent, J.; Maurice, X.; Brèthes, L. 3D Point Cloud Analysis for Detection and Characterization of Defects on Airplane Exterior Surface. J. Nondestruct. Eval. 2017, 36, 1–17. [Google Scholar] [CrossRef] [Green Version]
  12. Kolmogorov, A.N. Wienersche Spiralen und Einige Interessante Kurven im Hilbertschen Raum. C.R. (Doklady) Acad. URSS (N.S.) 1940, 26, 115–118. [Google Scholar]
  13. Benoit, B.M.; Van Ness, J.W. Fractional Brownian Motions, Fractional Noises and Applications. SIAM Rev. 1968, 10, 422–437. [Google Scholar]
  14. Amosov, O.S.; Amosova, S.G.; Muller, N.V. Identification of potential risks to system security using wavelet analysis, the time-and-frequency distribution indicator of the time series and the correlation analysis of wavelet-spectra. In Proceedings of the 2018 International Multi-Conference on Industrial Engineering and Modern Technologies (FarEastCon), Vladivostok, Russia, 3–4 October 2018. [Google Scholar]
  15. Stepanov, O.A.; Motorin, A.V. Comparison of methods for identifying sensor error models based on Allan variations and nonlinear filtering algorithms. In Proceedings of the XXI St. Petersburg International Conference on Integrated Navigation Systems, St. Petersburg, Russia, 2–7 October 2018; pp. 98–103. [Google Scholar]
  16. Bar-Shalom, Y. Estimation with Applications to Tracking and Navigation; Bar-Shalom, Y., Li, X.-R., Kirubarajan, T., Eds.; John Wiley & Sons: New York, NY, USA, 2001; 558p. [Google Scholar]
  17. Gosteva, N.D.; Litvinenko, Y.A. Research of a mathematical model of care of a two-stage float gyroscope. In Proceedings of the XV Conference of Young scientists “Navigation and Motion Control”, St. Petersburg, Russia, 2–7 October 2018; pp. 117–125. [Google Scholar]
  18. Amosov, O.S.; Baena, S.G. Wavelet Based Filtering of Mobile Object Fractional Trajectory Parameters. In Proceedings of the IEEE International Conference on Control and Automation (IEEE ICCA 2017), Ohrid, Macedonia, 3–6 July 2017; pp. 118–123. [Google Scholar]
  19. Amosov, O.S.; Amosova, S.G.; Magola, D.S. Identification of information recourses threats based on intelligent technologies, fractal and wavelet analysis. In Proceedings of the 2018 IEEE International Conference on Applied System Innovation (ICASI), Chiba, Japan, 13–17 April 2018; pp. 528–531. [Google Scholar]
  20. Doering, E. NI myRIO Project Essentials Guide; Rose-Hulman Institute of Technology; National Technology and Science Press, 2016; p. 249. Available online: https://download.ni.com/evaluation/academic/myRIO_project_essentials_guide__Feb_09_2016___optimized.pdf (accessed on 2 September 2021).
  21. MEMS Motion Sensors: Ultra-Stable Three-Axis Digital Output Gyroscope L3G4200D. Available online: https://www.mouser.com/datasheet/2/389/l3g4200d-954834.pdf (accessed on 2 September 2021).
  22. Abry, P.; Flandrin, P.; Taqqu, M.S.; Veitch, D. Self-similarity and long-range dependence through the wavelet lens. Theory Appl. Long-Range Depend. 2003, 1, 527–556. [Google Scholar]
  23. Bardet, J.-M.; Lang, G.; Oppenheim, G.; Philippe, A.; Stoev, S.; Taqqu, M.S. Semi-parametric estimation of the long-range dependence parameter: A survey. Theory Appl. Long-Range Depend. 2003, 557, 557–577. [Google Scholar]
  24. Flandrin, P. Wavelet analysis and synthesis of fractional Brownian motion. IEEE Trans. Inf. Theory 1992, 38, 910–917. [Google Scholar] [CrossRef]
  25. Digital Accelerometer ADXL345. Available online: https://www.analog.com/media/en/technical-documentation/data-sheets/ADXL345.pdf (accessed on 2 September 2021).
  26. Kuzmin, S.Z. Digital Radar. Introduction to the Theory; KVITS Publishing House: Kiev, Ukraine, 2000; 428p. [Google Scholar]
  27. Amosov, O.S. Markov sequence filtering on the basis of bayesian and neural network approaches and fuzzy logic systems in navigation data processing. J. Comput. Syst. Sci. Int. 2004, 43, 551–559. [Google Scholar]
  28. Amosov, O.S. Peculiarities of stochastic processes with fractal properties and their applications in problems of navigation information processing. In Proceedings of the 25th Saint Petersburg International Conference on Integrated Navigation Systems, ICINS 2018, Saint Petersburg, Russia, 28–30 May 2018; Volume 137790, pp. 1–5. [Google Scholar]
  29. Stepanov, O.A.; Amosov, O.S. Nonrecurrent Linear Estimation and Neural Networks. IFAC Proc. Vol. 2004, 37, 213–218. [Google Scholar] [CrossRef]
  30. Lv, X.; Duan, F.; Jiang, J.J.; Fu, X.; Gan, L. Deep Metallic Surface Defect Detection: The New Benchmark and Detection Network. Sensors 2020, 20, 1562. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Chu, H.-H.; Wang, Z.-Y. A vision-based system for post-welding quality measurement and defect detection. Int. J. Adv. Manuf. Technol. 2016, 86, 3007–3014. [Google Scholar] [CrossRef]
  32. Fan, Y.; Park, U.; Udpa, L.; Ramuhalli, P.; Shih, W.; Stockman, G.C. Automated rivet inspection for aging aircraft with magneto-optic imager. Electromagn. Nondestruct. Eval. (IX) 2005, 25, 185–194. [Google Scholar]
  33. Lingvall, F.; Stepinski, T. Automatic detecting and classifying defects during eddy current inspection of riveted lap-joints. NDT&E Int. 2000, 33, 47–55. [Google Scholar]
  34. Amosov, O.S.; Amosova, S.G.; Iochkov, I.O. Defects Detection and Recognition in Aviation Riveted Joints by Using Ultrasonic Echo Signals of Non-Destructive Testing. In Proceedings of the 17th IFAC Symposium on Information Control Problems in Manufacturing, IFAC INCOM 2021, Budapest, Hungary, 7–9 June 2021; pp. 484–489. [Google Scholar]
Figure 1. Inertial sensors: (a) gyroscope L3G4200D; (b) accelerometer ADXL345.
Figure 1. Inertial sensors: (a) gyroscope L3G4200D; (b) accelerometer ADXL345.
Sensors 21 05960 g001
Figure 2. Gyroscope drift according to X coordinate (a) and the Hurst parameter with the gyroscope drift for 1 h (b).
Figure 2. Gyroscope drift according to X coordinate (a) and the Hurst parameter with the gyroscope drift for 1 h (b).
Sensors 21 05960 g002
Figure 3. Accelerometer bias and drift by x coordinate (a), Hurst parameter (b).
Figure 3. Accelerometer bias and drift by x coordinate (a), Hurst parameter (b).
Sensors 21 05960 g003
Figure 4. R.m.s. of estimation errors with different Hurst parameters (H.). (a) H = 0.1, (b) H = 0.5, (c) H = 0.9.
Figure 4. R.m.s. of estimation errors with different Hurst parameters (H.). (a) H = 0.1, (b) H = 0.5, (c) H = 0.9.
Sensors 21 05960 g004
Figure 5. Scheme of the organization’s telecommunication system with incoming and outgoing Internet traffic.
Figure 5. Scheme of the organization’s telecommunication system with incoming and outgoing Internet traffic.
Sensors 21 05960 g005
Figure 6. Characteristics of information and telecommunication systems: (a,c) show the analysis of 2-month, and weekly, series of intensity, respectively, in bytes; (b,d)—in packages.
Figure 6. Characteristics of information and telecommunication systems: (a,c) show the analysis of 2-month, and weekly, series of intensity, respectively, in bytes; (b,d)—in packages.
Sensors 21 05960 g006
Figure 7. Echo signals of ultrasonic flaw detector during non-destructive testing of rivet joints: (a) is a normal rivet; (b) defect in the middle of the rivet at the 80% undercut; (c) defect in the middle of the rivet at the 50% undercut; (d) defect at the rivet cap base at the 50% undercut.
Figure 7. Echo signals of ultrasonic flaw detector during non-destructive testing of rivet joints: (a) is a normal rivet; (b) defect in the middle of the rivet at the 80% undercut; (c) defect in the middle of the rivet at the 50% undercut; (d) defect at the rivet cap base at the 50% undercut.
Sensors 21 05960 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Amosov, O.S.; Amosova, S.G. Peculiarities and Applications of Stochastic Processes with Fractal Properties. Sensors 2021, 21, 5960. https://doi.org/10.3390/s21175960

AMA Style

Amosov OS, Amosova SG. Peculiarities and Applications of Stochastic Processes with Fractal Properties. Sensors. 2021; 21(17):5960. https://doi.org/10.3390/s21175960

Chicago/Turabian Style

Amosov, Oleg Semenovich, and Svetlana Gennadievna Amosova. 2021. "Peculiarities and Applications of Stochastic Processes with Fractal Properties" Sensors 21, no. 17: 5960. https://doi.org/10.3390/s21175960

APA Style

Amosov, O. S., & Amosova, S. G. (2021). Peculiarities and Applications of Stochastic Processes with Fractal Properties. Sensors, 21(17), 5960. https://doi.org/10.3390/s21175960

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