Evolutionary Optimization of Case-Based Forecasting Algorithms in Chaotic Environments

The problem of dynamic adaptation of prediction algorithms in chaotic environments based on identification of the situations-analogs in the database of retrospective observations is considered. Under conditions of symmetrical and unsymmetrical chaotic dynamics, traditional computational schemes of precedent prediction turn out to be ineffective. In this regard, a dynamic adaptation of precedent analysis algorithms based on the method of evolutionary modeling is proposed. Implementation of the computational precedent prediction scheme for chaotic processes as well as the evolutionary modeling method are described.


Introduction
The state evolution of a controlled object immersed in an instable medium is described by an additive two-parameter model with a system component (represented by the realization of dynamic chaos) and a random component (represented by a nonstationary stochastic process). Such dynamics are quite typical, for example, for turbulent gas-dynamic and hydrodynamic media [1][2][3][4][5][6]. The presence of a chaotic component in the series of observations significantly hinders the prediction problem, because it contradicts the basic assumption of the probabilistic-statistical paradigm. Symmetry can play an important role in the prediction of nonlinear systems. This implies the repeatability of the experimental conditions under identical conditions. As a result, standard extrapolation computational forecasting schemes based on statistical data analysis are ineffective.
An alternative method based on the precedent data analysis was introduced in [7]. This approach represents the generalization of the nearest neighbor method [8][9][10][11][12] and relies on machine learning techniques and pattern recognition [13][14][15]. The precedent data analysis is based on the hypothesis that if the situation, which is similar to the current situation, has occurred before, the consequences can be used as predictive scenarios. This approach is close to the processes, which takes place in human consciousness, and is applicable to many artificial intelligence tasks.
However, due to the weak recurrence of processes in chaotic environments, the proposed method cannot provide a high reliability of the prediction. In this regard, this work proposes an improved version of precedent analysis, based on dynamic adaptation of the precedent detection algorithm. A modified method of evolutionary modeling proposed in [16] is used as an adaptation mechanism, which has found application in many data analysis problems [17][18][19][20].

Materials and Methods
The time series of observations Y(k), k = 1, · · · , N, formed by the monitoring system of a real physical process, which takes place in an unstable environment, is used as the initial data for the problem to be solved. 2 of 15 Traditional models of direct observations of the dynamic system state can be described by the additive relation of the following form, Y(k) = X(k) + υ(k), k = 1, · · · , N (1) where the system component X(k) is an unknown deterministic process to be identified, and υ(k) is a random process which simulates system noise or measurement errors. Assumptions concerning the probabilistic nature of the noise component are acceptable in many applications. In particular, the typical assumptions are about the stationarity, ergodicity, independence, and Gaussianity of the process υ(k), k = 1, · · · , N. In order to justify the normality of the distribution, references to the central limit theorems of probability theory are usually used [21,22].
If the above-mentioned conditions are fulfilled, then recovery of a system component is reduced to a procedure of sequential estimation of model parametersX(k), k = 1, · · · , N based on well-known methods of statistical data analysis-the least squares method, maximum likelihood method, etc. The mathematical structure of the model is specified sequentially, if necessary. The value ofX(k + τ), k = 1, · · · , N, where τ is the forecast interval, is used as a prognosis of the system state at the moment of time k + τ.
The fundamental differences of observation series (1), which are typical for the considered problem, are the system component X(k), k = 1, · · · , N of the observed process represents a realization of dynamic chaos and usually takes the form of an oscillatory non-periodic process with a set of false local trends of indefinite duration and the noise component υ(k), k = 1, · · · , N of the observed process is a nonstationary heteroscedastic process with a non-degenerate autocorrelation function that is timedependent.
Under these conditions, traditional statistical estimates of model parameters not only lose their effectiveness, but may not be valid at all [23,24].
In this regard, an approach based on the precedent analysis of data was proposed [7]. It is based on the typical approach of epistemology, based on the assumption that identical dynamic states of the system correspond to statistically similar after-effects.
To describe the current dynamic state of the controlled object a sliding observation window is used covering both the current observation and L + 1 previous samples. The size of the observation window L in traditional statistical studies is estimated based on the minimization of the total error squared, including both the standard deviation and the magnitude of the bias due to dynamic errors. However, this approach is unacceptable for chaotic environments due to the non-stationarity of the observational series. The value of the state window should be considered as an option to be specified during model adaptation.
The search for analogous situations is performed by sequentially browsing the retrospective database with a sliding scan window: where N S is the size of the retrospective database used to locate the analogs. During the scanning process, the search of M a analogs with the numbers i * 1 , · · · , i * M a corresponding to the smallest values of the similarity measure is performed. The scanning window with the minimum value of similarity metrics is considered as a precedent. As such similarity measures, distances of the standard deviation It should be noted that in cases (4) and (5), the centered values of observations (CVO) in windows (2) and (3) are used. In addition, in the case of significant heteroscedasticity of the observed process, it is advisable to normalize the state and scanning windows by the estimates of the standard deviation.
In some problems, which are focused on estimation of a local trend, the difference between the coefficients of linear approximation of the observation series in windows (2) and (3) can serve as a measure of similarity: The aftereffect with following precedent window is used as a prediction: i.e., behind the scanning window with the minimum value of the selected metric, The probability of several equal minima is very low. It is only possible at a rough meshing. In this case any minimum can be used as a precedent indicator. To increase the stability of the predicted result, the average value of aftereffects, which follow the obtained windows-analogs, can be used as a prediction A natural generalization of the prediction (9) is a weighted averaging proportional to the degree of relative similarity with respect to the best result (precedent).
It is impossible to obtain analytical relations for estimating the prediction efficiency in the chosen class of observational series due to the infinite variety of possible variations of the chaotic process. In this connection, the results of numerical studies were presented in [7], which showed a relatively low prediction accuracy. In particular, even in the best case for averaging over M a = 3 analogs, the relative accuracy of the prediction over an interval of 200 min fluctuates around 10-25%.
A variant of developing the precedent prediction based on dynamic adaptation of the model parameters and structure at each step of the prediction is proposed. The fullest degree of optimization adaptation is given by the method of full enumeration of values of model parameters in the range of their admissible variations. The total number of variants for such a run is determined by the number of model parameters N M and the number of variants for the value of each of the adapted parameters N Var (j), j = 1, · · · , N M . Obviously, with an increase in the parameters number and the accuracy of their determination, the computation amount increases exponentially. The solution of the problem becomes problematic, even for modern computers.
A possible alternative to the full search is based on a random search in N M space of parameter values. However, this method is characterized by a rather slow convergence and does not guarantee finding the best variant for multidimensional problems.
In the present work, algorithms of the evolutionary modeling method, adapted for specifics of the precedent prediction problem under chaotic dynamics, are used as the main tool for dynamic adaptation.

Implementation of the Computational Precedent Prediction Scheme for Chaotic Processes
As an example of the implementation of a precedent prediction scheme, the problem of predicting a series of observational variations in the normalized density of turbulent flow is considered (Figure 1). The "star" symbol divides the process into two parts. The first part with an interval of 200 days is used as a database of retrospective data, designed to search for analogous situations. The second part with an interval of 100 days is used to test the prediction algorithm.

Implementation of the Computational Precedent Prediction Scheme for Chaotic Processes
As an example of the implementation of a precedent prediction scheme, the problem of predicting a series of observational variations in the normalized density of turbulent flow is considered (Figure 1). The "star" symbol divides the process into two parts. The first part with an interval of 200 days is used as a database of retrospective data, designed to search for analogous situations. The second part with an interval of 100 days is used to test the prediction algorithm. The precedent prediction is implemented according to the following algorithm: 1. At each k-th step of the prediction the current state window ( ) (2) with length of L samples is formed.
3. According to (8), the scanning window corresponding to the minimum value of similarity metrics is found. Subsequent sequence (7) of length in samples is taken as a prediction for the k-th situation.
4. Prediction refinement can be carried out on the basis of (9) by averaging of the aftereffects for the situation-analogs with the lowest values of similarity metrics. Figure 2 shows the results of the program based on the above algorithm in terms of identifying the situations of analogs. The precedent prediction is implemented according to the following algorithm: 1. At each k-th step of the prediction the current state window Y C (k) (2) with length of L samples is formed.
3. According to (8), the scanning window corresponding to the minimum value of similarity metrics is found. Subsequent sequence (7) of length in samples is taken as a prediction for the k-th situation.
4. Prediction refinement can be carried out on the basis of (9) by averaging of M a the aftereffects for the situation-analogs with the lowest values of similarity metrics. Figure 2 shows the results of the program based on the above algorithm in terms of identifying the situations of analogs. In all four diagrams of Figure 2, the same process corresponding to the real situation and the four analog situations with the lowest values of similarity metrics identified during the process of scanning the retrospective database are presented. Figure 3 shows the aftereffects corresponding to the analogs found, which are used as predictions. For comparisons of their efficiency, the same graphs show the dynamics of the real process in the prediction section. In all four diagrams of Figure 2, the same process corresponding to the real situation and the four analog situations with the lowest values of similarity metrics identified during the process of scanning the retrospective database are presented. Figure 3 shows the aftereffects corresponding to the analogs found, which are used as predictions. For comparisons of their efficiency, the same graphs show the dynamics of the real process in the prediction section. In all four diagrams of Figure 2, the same process corresponding to the real situation and the four analog situations with the lowest values of similarity metrics identified during the process of scanning the retrospective database are presented. Figure 3 shows the aftereffects corresponding to the analogs found, which are used as predictions. For comparisons of their efficiency, the same graphs show the dynamics of the real process in the prediction section.  Figure 3 shows the previously mentioned feature of chaotic dynamics, which does not provide a simple linear conclusion that the degree of similarity of the situation-indicator unambiguously corresponds to the degree of similarity of the aftereffect. In particular, it can be shown that the prediction corresponding to the precedent situation and corresponding to the minimum value of the chosen similarity metric is not the best one.
It should be noted that the similarity of the aftereffects for similar situations is not guaranteed for stationary random processes either, but on average this requirement will be fulfilled. The prediction improvement of chaotic processes by averaging aftereffects is not guaranteed and requires additional research.
Nevertheless, the above example shows possibility of application of the precedent analysis for the problem of chaotic process forecasting. The improvement of this approach on the basis of certain methods of structural and parametric adaptation of the model used to form the prediction algorithm is to be discussed.

Preliminary Numerical Studies and Results
First, it is necessary to verify the validity of the precedent analysis in the problem of chaotic dynamics prediction. It is necessary to confirm the deterioration of the prediction accuracy with the increase in the number of the reference window. Otherwise, if the accuracy turns out to be unchanged on average, it will mean that the consequence cannot be used as a forecast.
The average value of the mean square error of the forecast is used as an estimate of the quality of the forecast: Calculations of CVO were performed for two metrics, (4) and (5), used in the process of search for window analogs. Averaging was carried out for N = 100 realizations. Table 1 shows the growth of the CVO prediction (10) by increasing number of the analog window, i.e., when its degree of similarity to the current situation decreases. It is not possible to obtain a uniformly increasing sequence under unstable dynamics of the observed process. Nevertheless, the obtained result allows us to conclude that it is possible to apply the precedent analysis for the class of observational series. It is necessary to search for N A windows in order to obtain a uniformly increasing sequence. A search for N A window-analogs of the current situation is performed for this purpose, the consequences of which are used as a preliminary estimate of the forecast. The final prediction Y(k + 1, k + τ) is made by piecewise averaging of the preliminary estimates. The gain in accuracy as compared to the prediction using the precedent aftereffects is estimated as a relative reduction of the prediction error. The corresponding result for N A = 1, 2, . . . , 10 window-analogs is shown in Figure 4. It shows that the proposed computational scheme makes it possible to significantly increase the prediction accuracy.  One of the important questions without an explicit answer for chaotic dynamics problems is the determination of the best scanning window size L. Obviously, an increase in the window size leads to an increase in the dynamical error due to poorly predictable and rapid changes ( ), = 1, ⋯ , in the chaotic system component of the observed process, while its reduction leads to a substantial increase in the influence of the random component ( ), = 1, ⋯ , in the model (1). The methodology for determining the optimal size L by minimizing the total square of the error used in such cases is also not viable due to the essential non-stationarity of the observed process. The approach for selection of the scanning window size based on numerical studies of large observational arrays is considered.
The problem of estimation the prediction quality for different observation sites and different values of the value L is to be discussed.
In this study, it is not the absolute prediction accuracy which is important, but the relative value that characterizes the change in the prediction quality depending on the growth of the state and scanning windows size. In this regard, only the precedent, i.e., the window with the highest degree of similarity to the current situation, will be used as a prototype.
The range of variation of the counts = [50, 400] is considered. As an indicator of forecast quality, as in the previous case, the average value of cumulative forecast errors ( ) =̃( + 1, + ) − ( + 1, + ) in 100-150 steps of the forecast is used. The results of the program to estimate the dependence of the forecast quality on the size of the scanning window are shown in Figure 5. The best value of the scanning window size lies in the range of change in the counts * = [250, 350]. However, due to the non-stationarity of the observed process, the best value is not a constant and, therefore, it is advisable to include this value in the list of refined parameters for the adaptive computational scheme. One of the important questions without an explicit answer for chaotic dynamics problems is the determination of the best scanning window size L. Obviously, an increase in the window size leads to an increase in the dynamical error due to poorly predictable and rapid changes X(k), k = 1, · · · , N in the chaotic system component of the observed process, while its reduction leads to a substantial increase in the influence of the random component υ(k), k = 1, · · · , N in the model (1). The methodology for determining the optimal size L by minimizing the total square of the error used in such cases is also not viable due to the essential non-stationarity of the observed process. The approach for selection of the scanning window size based on numerical studies of large observational arrays is considered.
The problem of estimation the prediction quality for different observation sites and different values of the value L is to be discussed.
In this study, it is not the absolute prediction accuracy which is important, but the relative value that characterizes the change in the prediction quality depending on the growth of the state and scanning windows size. In this regard, only the precedent, i.e., the window with the highest degree of similarity to the current situation, will be used as a prototype.
The range of variation of the counts dL = [50, 400] is considered. As an indicator of forecast quality, as in the previous case, the average value of cumulative forecast errors in 100-150 steps of the forecast is used. The results of the program to estimate the dependence of the forecast quality on the size of the scanning window are shown in Figure 5. The best value of the scanning window size lies in the range of change in the counts dL * = [250, 350]. However, due to the non-stationarity of the observed process, the best value is not a constant and, therefore, it is advisable to include this value in the list of refined parameters for the adaptive computational scheme.
Another important aspect to be investigated the determination of the feasibility and parameters of the preliminary smoothing of the observational series. According to the previously used computational scheme, search for window analogs was carried out by direct unsmoothed measurements. Unsmoothed measurements were used as prediction, which does not correspond to the traditional concept of prediction, which implies restoration of the system component of the model (1). In this connection, it is expedient to implement the whole process of the precedent prediction on smoothed data, formed by one or another method of dynamic filtering. Then, the search process in the retrospective data will be Another important aspect to be investigated the determination of the feasibility and parameters of the preliminary smoothing of the observational series. According to the previously used computational scheme, search for window analogs was carried out by direct unsmoothed measurements. Unsmoothed measurements were used as prediction, which does not correspond to the traditional concept of prediction, which implies restoration of the system component of the model (1). In this connection, it is expedient to implement the whole process of the precedent prediction on smoothed data, formed by one or another method of dynamic filtering. Then, the search process in the retrospective data will be focused on the system component of the process, which should increase the reliability of intermediate decisions on the choice of the situations-analogs.
In this work, the simplest exponential smoothing was used with a smoothing coefficient = 0.01-0.3. The main problem of sequential smoothing is that it leads to estimation bias caused by a dynamic delay of the smoothed process relative to the initial series of observations. Increasing the filter transfer coefficient leads to decreasing the dynamic bias, but at the same time the degree of smoothing of the observation series decreases. The latter statement is illustrated in Figure 6 for two values of the filter coefficient 1 = 0.01 and 2 = 0.1.
In this case, it is important to know the changes in the quality of the precedent prediction when switching to smoothed observational series. To this end, the problem of prediction quality estimation is solved again, provided that the search for precedent and counterpart situations is performed on the basis of the similarity analysis of smoothed versions of the observed process. Exponential smoothing (11)   In this work, the simplest exponential smoothing was used with a smoothing coefficient α = 0.01-0.3. The main problem of sequential smoothing is that it leads to estimation bias caused by a dynamic delay of the smoothed process relative to the initial series of observations. Increasing the filter transfer coefficient α leads to decreasing the dynamic bias, but at the same time the degree of smoothing of the observation series decreases. The latter statement is illustrated in Figure 6 for two values of the filter coefficient α 1 = 0.01 and α 2 = 0.1. In this case, it is important to know the changes in the quality of the precedent prediction when switching to smoothed observational series. To this end, the problem of prediction quality estimation is solved again, provided that the search for precedent and counterpart situations is performed on the basis of the similarity analysis of smoothed versions of the observed process. Exponential smoothing (11) was performed for transfer coefficients with values of α = 0.01, 0.025, 0.05, 0.1, 0.25, 0.5, 0.75, and 1.0. The degree of smoothness of the process is an inverse value to the value of the filter gain.
The dependence of the average prediction slope (10) from the value of the exponential filter gain is shown in Figure 7. The average estimate of the CVO for the non-smoothed process (dashed line), with changes of a non-systemic nature caused by the non-stationary nature of the observation series is shown.
The accuracy of the precedent prediction can be significantly improved by filtering the noise of the observations. Nevertheless, the degree of smoothing of the observations is limited by the growth of the bias. In the general case, the filter transfer coefficient α, as well as other optional parameters, is advisable to be included in the list of parameters subject to parametric optimization.

Evolutionary Modeling Method Applied to the Problem of Dynamic Adaptation of the Precedent Prediction Algorithm
The technology of evolutionary modeling proposed in [16] and developed in a large number of studies in the field of data analysis [17][18][19][20] is based on the simulation of the Darwinian self-developing process of evolution of biological populations using the principles of variability and selection.
The possibility of algorithmic realization of these principles as applied to the problem under consideration to be considered.
The principle of variability is based on the emergence of relatively small random changes in model parameters during the formation of a set of descendants of i-th generation {M D } i from a set of ancestors {M A } i models. In formalized formulation it is a question of formation and application of i-th generation models' variability and reproduction operator, where N A is the number of ancestor models, N D = N A k B is the number of descendant models in each generation, and k B is the coefficient of multiplication of ancestor models. Combination of sets of ancestor and descendant models forms new generation of mod- subjected to testing, selection, and selection. The number of changes of generations is set beforehand, or determined from the condition of stopping the growth of the efficiency index of the problem to be solved.
The number of modifiable parameter, the type, and value of the modifiable parameter are determined by casting a random variable r ∈ U[0, 1] according to the scheme adopted in the Monte Carlo method.
Three methods of modifying the selected parameter are applied: 1. Natural variability. This type of variability is the most common, occurs with probability P 1 > 0.9 − 0.95, and represents relatively small random fluctuations in model parameters. Natural variability is usually formed as a random variable subject to the normal law a ik ∈ N a i−1,k , s 2 k , where a ik is the value of k-th parameter in i-th model, a i−1,k is the value of the same parameter in the previous generation of models, s k is estimate of the standard deviation (RR) of possible variations a ik .
2. Parametric mutation. It differs from natural variability by low probability of occurrence P 2 < 0.05 − 0.1 and high values of a changeable parameter. It is estimated by drawing a random variable a ik ∈ N a i−1,k , 3s 2 k . 3. Nonparametric mutation. Characterized by unlikely P 2 < 0.05 changes in model structure. In the simplest case, it can consist of a change in the order or dimension of the model.
The problem of variability can be solved in different ways. In nature, in most cases formation of progeny genome is performed by combining genes from genomes of different sex parents. Selection of offspring parameters is performed by random selection of parameters from one or another parent. This approach is used when genetic algorithms of sequential optimization are used.
The principle of selection and selection consists in the sequential selection of the best models for a new generation of parent models. The selection process is carried out on the basis of the selected performance indicator of the selected model E f f M i,k . Values of similarity metrics (4)-(6) can serve as intermediate performance indicators. However, if the metrics themselves are attributed to the number of modifiable parameters, it is obviously incorrect to compare the degree of similarity on the basis of different metrics. In this connection, it is reasonable to choose as the terminal performance indicator an estimate of prediction accuracy based on the chosen variant of the model.
The selection and selection operator for the i-th generation of models can be represented as where E f f (M 1 ) ≥ · · · ≥ E f f (M N a ) ≥ E f f M j , ∀j > N a , N a is number of "surviving" models, which are allowed for further reproduction-modification; N G = N a (1 + k B ) is number of strategies of one generation. If M 0 = M{S 0 , a 0 } is a particular variant of a predictive model given with given structure S 0 and parameters a 0 , and accepted as a base "parent model", then the technology of evolutionary optimization is reduced to the cyclic repetition of the sequence of operators (12) and (13): As selection is performed by the criterion of superiority, strict optimality of the terminal solution is not guaranteed. However, the found solution will be the best one out of the whole set of the random search formed in the process of implementation of the evolutionary technology.

Numerical Studies of Evolutionary Adaptation of the Precedent Data Analysis Algorithm
Implementation of the precedent prediction algorithm described above was supplemented by an adaptation contour implemented on the basis of the evolutionary modeling method. The adaptation was reduced to the cyclic alternation (14) of the operators of variability (14), and selection and selection (13).
Varying type was performed with a probability of 0.9 for small natural variability, 0.05 for parametric mutation, and 0.05 for structural mutation. The vector of parameters [L, Shi f t, α, n A , nMetica] acted as genome, where L is the scan window size, Shi f t is the scan step value, α is exponential filter transfer coefficient (11), n A is the number of analog windows used in prediction formation, and nMetica is the number of similarity metric used. The first three parameters of genome were subjected to natural variability in variation ranges dL = [50, 400], dShi f t = [25,200], dα = [0, 1] and parametric mutation, at which the given variation ranges were tripled. Nonparametric mutation was carried out for parameters n A , changing in the integer range of 1 to 6 and nMetica, with an integer range of variation of 1 to 3. As metrics the relations (4) and (5) and the quadratic metric with normalization of scanning windows and the state of the estimation of the number of observations in these windows were used: where std is the operator for calculating CVO.
The value of the model reproduction coefficient kb was chosen from the range [3 − 5], the number of surviving models in each generation N A = 3 − 4. The given ranges were chosen for empirical considerations, and their further increase does not improve convergence of estimation, but leads to a significant increase in the running time of the program.
As an example, Figure 8 shows the convergence of forecast quality index when solving the precedent prediction problem using evolutionary adaptation of prediction models.
From the above graph, we can conclude that it is possible in principle to significantly improve the accuracy when solving the problem of predicting chaotic processes due to the evolutionary adaptation technology considered in the paper.
It is important to note that the result obtained refers only to a specific realization of a chaotic process; an attempt to transfer the found values of the model parameters to another prediction section almost always leads to a substantial decrease in the accuracy of the predictive estimates. In addition, when assessing the accuracy of the generated result, it essentially depends on the substantive formulation of the problem for the benefit of which the prediction is carried out. A relative increase in the prediction accuracy does not yet mean its sufficiency, for example, for effectively solving the problem of proactive control of the technological process. chaotic process; an attempt to transfer the found values of the model parameters to another prediction section almost always leads to a substantial decrease in the accuracy of the predictive estimates. In addition, when assessing the accuracy of the generated result, it essentially depends on the substantive formulation of the problem for the benefit of which the prediction is carried out. A relative increase in the prediction accuracy does not yet mean its sufficiency, for example, for effectively solving the problem of proactive control of the technological process.

Specific Features of Evolutionary Adaptation of the Precedent Analysis Algorithm for Chaotic Processes
Numerous applications of evolutionary modeling and related genetic algorithms [16][17][18][19][20] enable us to find effective solutions for nonstationary processes with inertia or other stabilizing factors. With respect to turbulent gas-dynamic flows, with low degree of inertia and connectivity, the chaotic component is dominant and high accuracy of extrapolation prediction should not be expected.
A preliminary evaluation of the precedent prediction has shown the possibility in principle of such an approach, but in a chaotic environment, the degree of similarity of the precedent does not guarantee the similarity of the found consequence to the predicted situation. In this regard, unlike the traditional approach based on the nearest-neighbor method, the task consists in finding analogous situations, relying not only on the degree of its similarity according to the chosen proximity criterion, but, above all, on choosing such a model, which would provide the greatest terminal efficiency, evaluated through prediction accuracy.
An alternative solution to the proposed option consists in the a priori formation of a pool of prognostic algorithms, including sets of solving rules, lists of modifiable parameters, and their critical values and ranges of variation. The advantage of this approach is the convenience of its implementation. However, it limits the number of freedom degrees of machine choice; there is no possibility to obtain radically new, unexpected solutions, not foreseen by a programmer or a system analyst. It is clear that with any regularization, any additional set of constraints can close access to unexpected original solutions.
At the same time, the complete removal of constraints in the random formation of predictive algorithms will lead to an infinite number of meaningless computational

Specific Features of Evolutionary Adaptation of the Precedent Analysis Algorithm for Chaotic Processes
Numerous applications of evolutionary modeling and related genetic algorithms [16][17][18][19][20] enable us to find effective solutions for nonstationary processes with inertia or other stabilizing factors. With respect to turbulent gas-dynamic flows, with low degree of inertia and connectivity, the chaotic component is dominant and high accuracy of extrapolation prediction should not be expected.
A preliminary evaluation of the precedent prediction has shown the possibility in principle of such an approach, but in a chaotic environment, the degree of similarity of the precedent does not guarantee the similarity of the found consequence to the predicted situation. In this regard, unlike the traditional approach based on the nearest-neighbor method, the task consists in finding analogous situations, relying not only on the degree of its similarity according to the chosen proximity criterion, but, above all, on choosing such a model, which would provide the greatest terminal efficiency, evaluated through prediction accuracy.
An alternative solution to the proposed option consists in the a priori formation of a pool of prognostic algorithms, including sets of solving rules, lists of modifiable parameters, and their critical values and ranges of variation. The advantage of this approach is the convenience of its implementation. However, it limits the number of freedom degrees of machine choice; there is no possibility to obtain radically new, unexpected solutions, not foreseen by a programmer or a system analyst. It is clear that with any regularization, any additional set of constraints can close access to unexpected original solutions.
At the same time, the complete removal of constraints in the random formation of predictive algorithms will lead to an infinite number of meaningless computational schemes. Waiting for any sensible solution to emerge would require a time commensurate with real biological evolution.
Evolutionary technology, like the whole probabilistic-statistical paradigm, is implicitly oriented toward the comfortable hypothesis of repeatable experiments under unchanging or slowly changing conditions. The transition to non-stationary, and even more chaotic processes, inevitably destroys the statistical basis of predictive technology. To some extent, this also applies to evolutionary modeling. However, in the nature of chaos, as a rule, there are some regularizing effects that reduce the degree of total uncertainty. If evolutionary technology can isolate, at least not explicitly, and use such latent regularities, the task of building effective predictive technologies may be feasible.
In addition, the application of the evolutionary computational scheme will make it possible to answer the question about the fundamental admissibility of one or another class of predictive technologies.

Analysis of the Dynamic and Statistical Characteristics of an Unstable Gas-Dynamic System
The series of observations were obtained during the monitoring of the technological process (TP) which takes place in an unstable gas-dynamic environment. In this case, the instability of the medium is defined as the chaotic nature of the dynamic process which describes the unsymmetrical changes in the phase state parameters (or state vector) of the simulated object. The data used are centered and normalized deviations of measured physical parameters from their values defined by technological regulations of the TP.
As a system component a posteriori conditional mean, obtained on the basis of one of the known algorithms of statistical smoothing, was used; in particular, polynomials of a given order with parameter fitting using the least squares method (LQM). Plots of the system component X(t) of the observation model (1), obtained by MNA fitting of the corresponding noisy curves (upper figures), and plots of the noise component V(t) (lower figures) are presented in Figure 9. In both cases, a polynomial approximation of the 15th degree was used to extract the system component. At the same time, generalized inverse was used in solving the system of normal equations, as the use of approximating polynomials of degree 4 and higher led to degeneracy of the observation matrix.
changing or slowly changing conditions. The transition to non-stationary, and even more chaotic processes, inevitably destroys the statistical basis of predictive technology. To some extent, this also applies to evolutionary modeling. However, in the nature of chaos, as a rule, there are some regularizing effects that reduce the degree of total uncertainty. If evolutionary technology can isolate, at least not explicitly, and use such latent regularities, the task of building effective predictive technologies may be feasible.
In addition, the application of the evolutionary computational scheme will make it possible to answer the question about the fundamental admissibility of one or another class of predictive technologies.

Analysis of the Dynamic and Statistical Characteristics of an Unstable Gas-Dynamic System
The series of observations were obtained during the monitoring of the technological process (TP) which takes place in an unstable gas-dynamic environment. In this case, the instability of the medium is defined as the chaotic nature of the dynamic process which describes the unsymmetrical changes in the phase state parameters (or state vector) of the simulated object. The data used are centered and normalized deviations of measured physical parameters from their values defined by technological regulations of the TP.
As a system component a posteriori conditional mean, obtained on the basis of one of the known algorithms of statistical smoothing, was used; in particular, polynomials of a given order with parameter fitting using the least squares method (LQM). Plots of the system component X(t) of the observation model (1), obtained by MNA fitting of the corresponding noisy curves (upper figures), and plots of the noise component V(t) (lower figures) are presented in Figure 9. In both cases, a polynomial approximation of the 15th degree was used to extract the system component. At the same time, generalized inverse was used in solving the system of normal equations, as the use of approximating polynomials of degree 4 and higher led to degeneracy of the observation matrix. An alternative approach to extract the system component is related to the application of spline approximation. The selection of spline approximation nodes was performed An alternative approach to extract the system component is related to the application of spline approximation. The selection of spline approximation nodes was performed manually at points close to local extrema and inflection points. An example of the formation of a reference curve for a one-day monitoring interval is shown in Figure 10.
The graphs show an oscillating, non-periodic, unsymmetrical process with multiple local trends which do not passes for explicit extrapolation computational schemes. If desired, this observation can be confirmed formally by using the eigenvalues of the linearization matrix in the neighborhood of the stationary point. In the general case, Lyapunov characteristic indices are used to study the behavior in the neighborhood of a random trajectory [25]. The Lyapunov exponents serve as a measure of chaotic behavior: if there are positive exponents, then the behavior of the dynamical system is chaotic.
Statistical analysis of the noise component was carried out based on known methods for testing statistical hypotheses [26,27]. The stationarity of the noise component was checked using the ascending and descending series criterion and the Abbe criterion. The stationarity hypothesis was rejected at the 0.95 confidence level. Pearson's criterion was used for normality checking. Analysis of the Gaussianity of the observations showed a slow convergence to a normal distribution. The main difference is in the weighting of the distribution tips, indicating the presence of anomalous observations. Symmetry 2021, 13, x FOR PEER REVIEW 14 of 16 manually at points close to local extrema and inflection points. An example of the formation of a reference curve for a one-day monitoring interval is shown in Figure 10. Figure 10. System component of a process based on cubic spline approximation (green-experimental data a one-day monitoring, blue-spline approximation).
The graphs show an oscillating, non-periodic, unsymmetrical process with multiple local trends which do not passes for explicit extrapolation computational schemes. If desired, this observation can be confirmed formally by using the eigenvalues of the linearization matrix in the neighborhood of the stationary point. In the general case, Lyapunov characteristic indices are used to study the behavior in the neighborhood of a random trajectory [25]. The Lyapunov exponents serve as a measure of chaotic behavior: if there are positive exponents, then the behavior of the dynamical system is chaotic.
Statistical analysis of the noise component was carried out based on known methods for testing statistical hypotheses [26,27]. The stationarity of the noise component was checked using the ascending and descending series criterion and the Abbe criterion. The stationarity hypothesis was rejected at the 0.95 confidence level. Pearson's criterion was used for normality checking. Analysis of the Gaussianity of the observations showed a slow convergence to a normal distribution. The main difference is in the weighting of the distribution tips, indicating the presence of anomalous observations.

Conclusions
Numerous applications of evolutionary modeling presented in the article allow the drawing of conclusions about the possibility of predicting the state of dynamic systems in unsymmetrical, unstable environments, whose dynamics are described by chaotic processes, based on technology of precedent data analysis. In this case, an increase in the accuracy of prediction can be achieved by means of the computational scheme of dynamic adaptation, based on the well-known method of evolutionary modeling.
In this case, the inertia, which is typical of any physical process, serves as an implicit regulating factor which allows one to reveal the local ordering in a chaotic process. By transition to inertia-free environments, such as the dynamics of quotations of currency instruments in the Forex market, the result is less effective.

Conclusions
Numerous applications of evolutionary modeling presented in the article allow the drawing of conclusions about the possibility of predicting the state of dynamic systems in unsymmetrical, unstable environments, whose dynamics are described by chaotic processes, based on technology of precedent data analysis. In this case, an increase in the accuracy of prediction can be achieved by means of the computational scheme of dynamic adaptation, based on the well-known method of evolutionary modeling.
In this case, the inertia, which is typical of any physical process, serves as an implicit regulating factor which allows one to reveal the local ordering in a chaotic process. By transition to inertia-free environments, such as the dynamics of quotations of currency instruments in the Forex market, the result is less effective.
The presence of "vertical" correlations in series of multidimensional chaotic processes and "horizontal" autocorrelations allowing dynamically optimizing the size of the scanning window should be considered as alternative regularizing factors.