Fractal Tent Map with Application to Surrogate Testing

: Discrete chaotic maps are a mathematical basis for many useful applications. One of the most common is chaos-based pseudorandom number generators (PRNGs), which should be computationally cheap and controllable and possess necessary statistical properties, such as mixing and diffusion. However, chaotic PRNGs have several known shortcomings, e


Introduction
Chaos theory presents a prosperous ground for investigation and has a great interdisciplinary influence.It finds applications in various fields of research and development, concerning both theory and practice, e.g., in physics, engineering, economics, and biology [1].Chaotic maps naturally allow for generating pseudorandom sequences, so they are widely used in designing pseudorandom number generators (PRNGs) for various purposes: cryptography [2], secure communications [3], numerical simulations and statistical modeling, noise generation [4], and image encryption [5,6].Unlike traditional PRNGs, such as the Mersenne Twister, which rely on number-theoretical methods, chaos-based PRNGs harness the inherent unpredictability of chaotic systems to produce number sequences with high complexity and better statistical properties [7].
In recent decades, a lot of studies have been carried out concerning PRNGs based on chaotic maps.For example, in [8], Rezk et al. proposed a reconfigurable FPGA-based chaotic PRNG based on the Lorenz and Lü chaotic systems.The authors of [9] proposed a hardware-based PRNG using the single skew tent map [10].Investigations show the high applicability of the coupled skew tent map [11] and the cross-coupled skew tent map [12] to pseudorandom number generation.The well-known chaotic Bernoulli map was investigated in a work [13], where the authors present the results of running standardized tests and prove that the proposed PRNG based on the Bernoulli map has good statistical properties.In another article [14], Moysis et  a pseudorandom bit generator by extracting multiple bits per iteration from the decimal part of a chaotic map.This is carried out by selecting the decimal part of the state in each iteration and comparing each digit separately to a threshold value.In another study [15], Moysis et al. propose a modification of the classical logistic map by introducing fuzzy triangular numbers.The authors show higher complexity compared to the classic logistic map and showcase several phenomena, like antimonotonicity and crisis.Another interesting approach to implementing PRNGs is discussed in a recent study by Kvitko et al. [16], which describes a PRNG scheme based on the Courbage-Nekorkin neuron and its application to chaotic path planning.The authors of [10] show that coupling polynomial chaotic maps expands the number of chaos parameters.To further increase the randomness and uniformity, a PRNG based on the hyperchaotic map with microcontroller implementation was proposed by Murillo-Escobar et al. [17].
Many state-of-the-art PRNGs use modified chaotic maps to obtain necessary statistical properties.However, these systems face several pitfalls arising from a variety of factors.One of them is falling in short time periods [7].To prevent this from happening, the authors of [7] decided to dynamically change the parameters of the chaotic system while increasing the randomness of the generated sequences.Another common shortcoming is a degeneration phenomenon, when the chaotic behavior of the underlying map deteriorates over time, leading to decreased randomness and increased predictability.Garcia-Bosque et al. [18] proposed a transform arithmetic to improve the properties of chaos-based PRNGs.Another drawback of chaos-based algorithms is that generating pseudorandom sequences using some popular approaches, e.g., using the least significant bits or digits of finite-precision numbers, is susceptible to rounding errors and discretization effects, as described in [19].To mitigate this effect, the authors proposed an algorithm for the calculation of the maximum number of bits that are suitable for pseudorandom sequence generation from the binary representation of floating-point numbers.Another issue with PRNGs based on chaotic maps is having a narrow range of parameters suitable for exhibiting dynamical chaos, which reduces the keyspace in encryption tasks.To expand the keyspace of the designed system, Zhao et al. [20] combined a sequence based on quantum random walks with chaotic system outputs, creating a new system with high randomization, unpredictability, and uniform output distribution.The digital implementation of chaotic maps also has a negative influence on the properties of PRNGs.Many researchers apply their efforts to overcome the degradation of chaotic dynamics while being implemented in hardware.For example, Merah et al. [21] proposed a circuit design based on the logistic chaotic map, which enhanced chaotic map randomness and made the cycle length longer using a special perturbation mechanism.
The known solutions to the abovementioned problems, such as additional periodic perturbation or periodic parameter change, seem rather artificial and overcomplicate the algorithm.They also require introducing several excessive parameters, like perturbation frequency, which need to be handled with care, so they are definitely not a "magic pill" in this regard.Another known way to increase the randomness of the output sequence is using a more complicated right-hand-side function in the basic chaotic map.Fractal functions that possess inherently complex structures are good candidates for this role [22].Apart from famous computationally complex variants like the Weierstrass function, several simpler fractal functions exist, for example, fractal piecewise linear functions.In this study, we introduce a fractal generalization of the tent map and apply it to the pseudorandom number generation and surrogate testing tasks.Noise generators (usually generating white noise) are often used as part of some larger system, such as a neural network or synthetic or surrogate data generators.It is very important to maintain a high quality of generated noise for the good operation of the entire system.Testing the signal with surrogate data is performed to check the hypothesis of signal nonlinearity.In our study, we explicitly show that the use of a PRNG as a noise generator improves the accuracy of the surrogate test for biological signal testing.The experimental work demonstrates that the designed generator outperforms algorithms based on well-known simpler maps and, surprisingly, the MATLAB built-in pseudorandom number generator.
The main contributions of this study can be summarized as follows: 1.
A novel fractal chaotic map is presented, being a generalization of the chaotic tent map with a fractal structure.

2.
A pseudorandom number generator (PRNG) was developed based on the fractal tent map proposed in this paper.The developed generator successfully passes all NIST tests.

3.
The developed PRNG was applied to the surrogate data testing task.The test results show that the use of the PRNG based on the fractal tent map shows better accuracy when testing data sequences for nonlinearity than the standard approach.
The rest of this paper is organized as follows.Section 2 contains a description of the materials and methods used in the study, including a description of the proposed chaotic fractal tent map.Section 3 describes the analysis results, including bifurcation diagrams of the proposed map and LLE plots.The results of NIST tests for the developed PRNG based on the proposed map are given, as well as an example of using the developed PRNG for surrogate testing.Finally, Section 4 concludes the paper.

Materials and Methods
In the context of dynamical systems, a map is a mathematical function that defines how a point in a space (e.g., Banach or Euclidean space) evolves over time.Chaotic maps exhibit dynamical chaos and can be characterized by their sensitivity to initial conditions, meaning that even small changes in the starting point can lead to vastly different outcomes, which is a hallmark of chaos theory.Discrete chaotic maps, also known as iterated maps, apply a function repeatedly to iterate points through discrete time.In this study, we introduce a novel chaotic fractal tent map and investigate its properties and possible applications to pseudorandom number generation.To compare the developed chaotic fractal map with existing solutions, we will also consider the well-known logistic map [23] and the tent map [24].

Logistic and Tent Maps
The logistic map is described by the following recurrent formula: where k is a parameter, and x is a state variable.The interval [0, 1] is covered when k = 4.The tent map can be described as follows: where µ is a positive real constant, and x is a state variable.The interval [0, 1] is covered when µ = 2. Figure 1 shows the graphs of the logistic map function (on the left) and the tent map function (on the right).

Proposed Fractal Tent Map
Let us propose the following modification of the tent map: such that the two initial intervals x n < 0.5 and x n ≥ 0.5 are iteratively divided into subintervals of lengths K i D, where i ∈ Z + , 0 < K ≤ 1, and 0 < D ≤ 0.5.The parameter K stands for the contraction ratio of the interval after each division, D stands for the initial division, and i is an iteration counter, standing for the number of divisions.The construction of the map is illustrated in Figure 2.An iterative procedure for searching for the value of x n+1 is given in pseudo-code in Algorithm 1.
In this algorithm, we first determine where x n lies relative to the value 0.5.Then, the two variables m and l, playing the role of the left and right bounds, respectively, are used

Proposed Fractal Tent Map
Let us propose the following modification of the tent map: such that the two initial intervals x n < 0.5 and x n ≥ 0.5 are iteratively divided into subintervals of lengths K i D, where i ∈ Z + , 0 < K ≤ 1, and 0 < D ≤ 0.5.The parameter K stands for the contraction ratio of the interval after each division, D stands for the initial division, and i is an iteration counter, standing for the number of divisions.The construction of the map is illustrated in Figure 2.An iterative procedure for searching for the value of x n+1 is given in pseudo-code in Algorithm 1. Fractal

Proposed Fractal Tent Map
Let us propose the following modification of the tent map: such that the two initial intervals x n < 0.5 and x n ≥ 0.5 are iteratively divided into subintervals of lengths K i D, where i ∈ Z + , 0 < K ≤ 1, and 0 < D ≤ 0.5.The parameter K stands for the contraction ratio of the interval after each division, D stands for the initial division, and i is an iteration counter, standing for the number of divisions.The construction of the map is illustrated in Figure 2.An iterative procedure for searching for the value of x n+1 is given in pseudo-code in Algorithm 1.
In this algorithm, we first determine where x n lies relative to the value 0.5.Then, the two variables m and l, playing the role of the left and right bounds, respectively, are used to search for an exact interval containing x n , starting from one of the two side intervals of length D. The initial values of m and l are 0 and D, respectively, for the left subinterval, and 1 − D and 1, respectively, for the right subinterval.The iterative modification of these bounds belongs to a class of dichotomic search algorithms, like a renowned binary search.When the interval is found where m < x n ≤ l, the equation of a line is used to estimate the value of x n+1 .This map can be considered a generalization of the tent map, and if ) is a hyperbolic iterated function system with the attractor 1 2 , and F f (x) becomes a fractal function [22].
Let us give some examples of the fractal and non-fractal behavior of F f (x).One may recall that the sum of a geometric progression written in terms of D and Therefore, the fractal behavior of F f (x) can be observed for any pair (D, K) satisfying give non-fractal behavior.Figure 3 illustrates phase portraits of the proposed fractal tent map for various values of (D, K) yielding non-fractal (upper row) and fractal (lower row) behavior.
By definition, a fractal function is a function whose graph has, in general, a nonintegral dimension [22,25].To ensure that the function F f is fractal, we estimate the Minkowski-Bouligand (box-counting) dimension of the proposed function's graph by the formula In this algorithm, we first determine where x n lies relative to the value 0.5.Then, the two variables m and l, playing the role of the left and right bounds, respectively, are used to search for an exact interval containing x n , starting from one of the two side intervals of length D. The initial values of m and l are 0 and D, respectively, for the left subinterval, and 1 − D and 1, respectively, for the right subinterval.The iterative modification of these bounds belongs to a class of dichotomic search algorithms, like a renowned binary search.When the interval is found where m < x n ≤ l, the equation of a line is used to estimate the value of x n+1 .
This map can be considered a generalization of the tent map, and if the system [0, 1], F f (x) is a hyperbolic iterated function system with the attractor 1 2 , and F f (x) becomes a fractal function [22].
Let us give some examples of the fractal and non-fractal behavior of F f (x).One may recall that the sum of a geometric progression written in terms of D and Therefore, the fractal behavior of F f (x) can be observed for any pair (D, K) satisfying give non-fractal behavior.Figure 3 illustrates phase portraits of the proposed fractal tent map for various values of (D, K) yielding non-fractal (upper row) and fractal (lower row) behavior.
By definition, a fractal function is a function whose graph has, in general, a nonintegral dimension [22,25].To ensure that the function F f is fractal, we estimate the Minkowski-Bouligand (box-counting) dimension of the proposed function's graph by the formula where N is the number of balls with the radius ε required to cover the function plot.For the parameters D = 1/3, K = 1/3, the numerical estimation gives the value

Surrogate Data Testing
Although PRNGs are typically used for cryptography purposes, they are also an excellent tool for generating noise, which can be used for simulating various processes.In this paper, we apply the proposed chaotic-map-based PRNG to surrogate data testing as a white noise generator.
The main idea of surrogate testing is to first formulate the null hypothesis H0 that the analyzed time series represents some linear process.Subsequently, surrogate datasets are generated that are consistent with the chosen null hypothesis, and finally, the test statistic for the original time series and the surrogate time series is calculated.The null hypothesis is rejected if the value of the test statistic of the original time series is significantly different from the value of the test statistic of the surrogate data.
Surrogate data testing is a popular method for testing the nonlinearity of time series.The null hypothesis is usually that the data were generated by a Gaussian linear stochastic process with constant coefficients [26].

Surrogate Data Testing
Although PRNGs are typically used for cryptography purposes, they are also an excellent tool for generating noise, which can be used for simulating various processes.In this paper, we apply the proposed chaotic-map-based PRNG to surrogate data testing as a white noise generator.
The main idea of surrogate testing is to first formulate the null hypothesis H0 that the analyzed time series represents some linear process.Subsequently, surrogate datasets are generated that are consistent with the chosen null hypothesis, and finally, the test statistic for the original time series and the surrogate time series is calculated.The null hypothesis is rejected if the value of the test statistic of the original time series is significantly different from the value of the test statistic of the surrogate data.
Surrogate data testing is a popular method for testing the nonlinearity of time series.The null hypothesis is usually that the data were generated by a Gaussian linear stochastic process with constant coefficients [26].
Generation of surrogate data.The surrogate data are generated by a Gaussian linear stochastic process to have the same specified properties as the observed time series, such as the mean, variance, and Fourier spectrum.This can be achieved by applying a Fourier transform to the original signal and then randomly mixing the phases and applying the inverse Fourier transform.The output is a surrogate, which, in terms of the spectrum, dispersion, and average value, repeats the original one but is linear [27].This approach is called Fourier-based surrogates and has many variations: the amplitude-adjusted Fourier transform (AAFT) method [28], iteratively refined surrogates (iAAFT) [29], and iterative multivariate surrogates (iAAFTn) [30].Another way to create surrogate data is to use the autoregressive (AR) model and the moving-average (MA) model.Using the autocorrelation method (Yule-Walker method) [31], an AR model is matched to the input data so that their autocorrelation function coincides.As is known, the autocorrelation function is related to the spectrum by the Wiener-Khinchin theorem [32], which states that the power spectral density of a stationary random process in a broad sense is the Fourier transform of the corresponding autocorrelation function.Thus, surrogates generated using the AR process will also have the same spectrum as the original data.
Finding test statistics.To confirm or refute the null hypothesis about the linearity of the signal, the third moment of the function is calculated, which characterizes the asymmetry of the distribution of function values ( x i ).The unbiased estimate of the third moment of the function is defined as where N is the length of the source data.This test statistic is calculated for all surrogate data and the original signal.
Hypothesis testing.Next, a histogram of the distribution of test statistics for the surrogate data is constructed.Hypothesis H0 is considered confirmed if the value of the test statistic of the original signal is close to the distribution of the test statistic of the surrogate data, and it is rejected if the test statistic of the original signal differs from the values of the test statistics for the surrogate data.

Results
This section presents the results of the research concerning both the chaotic and statistical properties of the proposed chaotic map in comparison with the other maps and PRNGs.

Bifurcation Analysis
To construct a 1D bifurcation diagram, we use a modified version of the density-color bifurcation analysis tool proposed in [33].Our modification consists of calculating the relative density of points for each parameter value (for each vertical line) separately.Also, histograms are calculated only over the entire range of values of the state variable.This approach allows us to mark the ranges of values of the state variable with greater contrast and more clearly identify orbits that are close to periodic.Bifurcation diagrams for the logistic and tent maps are shown in Figures 4 and 5, respectively.
The diagrams shown in Figures 4 and 5 without density coloring are often seen in the literature.In this study, they can serve as a reference in comparison with the proposed chaotic map.
For the chaotic fractal tent map, several bifurcation diagrams were constructed with the following parameters: for K = 1 and D ∈ [0; 1], for K = 0.5 and D ∈ [0; 1], and for K ∈ [0; 1] and D = 0.25, as given in Figure 6.The analysis of these bifurcation diagrams reveals several chaotic properties inherent to the proposed map.From Figure 6a, one can clearly see that, in the case of K = 1, for each D ≤ 0.5, the system's chaotic orbits cover the entire phase space, but it is only at the specific values D = 1/2 i , i ∈ N, that the density of the orbit is uniform in the entire range [0, 1] and reaches its maximum values.From Figure 6b, it is clear that the fractal behavior of the function F f (x) occurring when D = 0.25 yields the greatest density possible, and for D < 0.25, the dynamics degrade.For values slightly greater than D = 0.25, the density remains high.In Figure 6c, we observe the same picture concerning the fractal behavior of the function F f (x) near K = 0.5, which also yields the highest possible, and for K < 0.5, the dynamics also degrade.For values slightly greater than K = 0.5, the density remains high, as in the case of exploring the parameter D, which points out that a small perturbation of the system parameters toward increased values (e.g., related to rounding error) would not affect the dynamics, but a small perturbation of their values downward would prevent the system from exhibiting chaotic behavior, so this should be taken into account when implementing the chaotic fractal tent map in practical applications since the fractal behavior corresponds to the edge of chaotic dynamics.The analysis of these bifurcation diagrams reveals several chaotic properties inherent to the proposed map.From Figure 6a, one can clearly see that, in the case of K = 1, for each D ≤ 0.5, the system's chaotic orbits cover the entire phase space, but it is only at the specific values D = 1/2 i , i ∈ N, that the density of the orbit is uniform in the entire range [0, 1] and reaches its maximum values.From Figure 6b, it is clear that the fractal behavior of the function F f (x) occurring when D = 0.25 yields the greatest density possible, and for D < 0.25, the dynamics degrade.For values slightly greater than D = 0.25, the density remains high.In Figure 6c, we observe the same picture concerning the fractal behavior of the function F f (x) near K = 0.5, which also yields the highest possible, and for K < 0.5, the dynamics also degrade.For values slightly greater than K = 0.5, the density remains high, as in the case of exploring the parameter D, which points out that a small perturbation of the system parameters toward increased values (e.g., related to rounding error) would not affect the dynamics, but a small perturbation of their values downward would prevent the system from exhibiting chaotic behavior, so this should be taken into account when implementing the chaotic fractal tent map in practical applications since the fractal behavior corresponds to the edge of chaotic dynamics.The analysis of these bifurcation diagrams reveals several chaotic properties inherent to the proposed map.From Figure 6a, one can clearly see that, in the case of K = 1, for each D ≤ 0.5, the system's chaotic orbits cover the entire phase space, but it is only at the specific values D = 1/2 i , i ∈ N, that the density of the orbit is uniform in the entire range [0, 1] and reaches its maximum values.From Figure 6b, it is clear that the fractal behavior of the function F f (x) occurring when D = 0.25 yields the greatest density possible, and for D < 0.25, the dynamics degrade.For values slightly greater than D = 0.25, the density remains high.In Figure 6c, we observe the same picture concerning the fractal behavior of the function F f (x) near K = 0.5, which also yields the highest possible, and for K < 0.5, the dynamics also degrade.For values slightly greater than K = 0.5, the density remains high, as in the case of exploring the parameter D, which points out that a small perturbation of the system parameters toward increased values (e.g., related to rounding error) would not affect the dynamics, but a small perturbation of their values downward would prevent the system from exhibiting chaotic behavior, so this should be taken into account when implementing the chaotic fractal tent map in practical applications since the fractal behavior corresponds to the edge of chaotic dynamics.Figure 7 shows graphs of approximate and information entropy for the fractal tent map with the parameters K = 1, D ∈ [0, 1] (Figure 6a). Figure 7 shows graphs of approximate and information entropy for the fractal tent map with the parameters K = 1, D ∈ [0, 1] (Figure 6a).  Figure 8 gives the values of the largest Lyapunov exponents (LLEs) of the considered chaotic maps.One can see an increase in the LLE in Figure 8a, corresponding to the fractal tent map when D tends to zero, which relates to greater mixing properties due to the growing number of linear fragments in the function F f (x).To create a PRNG based on a chaotic map, we used the Box-Muller transformation, which allows one to transform a sample of random numbers from uniformly distributed to normally distributed.Thus, for two independent random samples x 1 and x 2 generated by a chaotic map, the Box-Muller transform will be The selected parameters for chaotic maps are as follows: k = 4 for the logistic map; µ = 2 for the tent map; K = 1/3 and D = 1/3 for the fractal tent map. Figure 8 gives the values of the largest Lyapunov exponents (LLEs) of the considered chaotic maps.One can see an increase in the LLE in Figure 8a, corresponding to the fractal tent map when D tends to zero, which relates to greater mixing properties due to the growing number of linear fragments in the function F f (x). Figure 8b, corresponding to the fractal tent map in the fractal regime, shows the same behavior.The two other panels illustrate LLEs depending on parameters for the logistic and tent maps.Notably, in the generalized tent map, high values of LLEs are much easier to obtain.To create a PRNG based on a chaotic map, we used the Box-Muller which allows one to transform a sample of random numbers from uniform normally distributed.Thus, for two independent random samples x 1 and a chaotic map, the Box-Muller transform will be The selected parameters for chaotic maps are as follows: k = 4 for t µ = 2 for the tent map; K = 1/3 and D = 1/3 for the fractal tent map.To create a PRNG based on a chaotic map, we used the Box-Muller transformation, which allows one to transform a sample of random numbers from uniformly distributed to normally distributed.Thus, for two independent random samples x 1 and x 2 generated by a chaotic map, the Box-Muller transform will be The selected parameters for chaotic maps are as follows: k = 4 for the logistic map; µ = 2 for the tent map; K = 1/3 and D = 1/3 for the fractal tent map.

NIST Tests
To transform a generated value from the floating-point representation into binary, we used the method proposed in [34].First, the generated value is incremented by 10 9 , and then it is rounded toward negative infinity; only the remainder from dividing this number by 2 is considered.The combination of values obtained in this way gives us a sequence suitable for testing with the NIST suite.NIST tests were performed for all studied PRNGs.Each of the PRNGs generated a sequence of 10 6 bits and was run 100 times.The significance level α = 0.01.The results obtained are presented in Table 1.The null hypothesis is that the generated sequence is random.If the value of P value ≥ α, the null hypothesis is accepted (i.e., the sequence is random).According to the NIST test package recommendations, the minimum pass rate for each statistical test, except for the random excursion (REVT) test, is 0.96.The minimum pass rate for the REVT test is 0.94.One can see that all the generators under study produce sequences with random characteristics for all cases considered.A key-space analysis of the proposed schemes is performed in order to check the generation algorithms' security against attacks of a brute-force type.The condition that one is required to satisfy is an algorithm having more than 2 100 different key values, which can be used to form streams of random data, that are different from each other [35,36].By taking into account the usage of the double-precision floating-point data type, with mantissa being represented by 53 bits, one can achieve a rough estimation of viable key spaces for proposed schemes (Table 2): Such key-space sizes were acquired by analyzing the different parameter values that can be used in order to generate data sequences.In the case of a logistic map, one might use various starting x 0 values in the range of [0; 1] and the k parameter in the range of roughly (3.55:4] as an initial parameter for outcome control; for the double-precision floating-point data type, this results in a potential key-space size of 2 64 × 2 9 = 2 61 .In the case of a tent map, one obtains a new parameter µ, which can potentially be chosen from anywhere in the range of (1:2), paired with an already-existing x 0 , resulting in 2 52 × 2 52 = 2 104 potential key variants.As for the fractal tent map, one obtains the same result as for the tent map due to the parameter being available to be chosen from the range of (0:1) and parameter D being dependent on it, thus resulting in a key space of 2 52 × 2 52 = 2 104 .The acquired According to the obtained results, one can state that the bias resistance of the proposed PRNGs might be considered high due to the Chi-square metric being significantly lower than that of a required value.
Figure 9 shows a comparison of Gaussian normal distributions and PRNG distributions of a logistic map, a tent map, and a fractal tent map.It can be seen that the studied PRNGs have normal or nearly normal distributions.
12 of 17 the range of (1:2), paired with an already-existing x 0 , resulting in 2 52 × 2 52 = 2 104 potential key variants.As for the fractal tent map, one obtains the same result as for the tent map due to the parameter being available to be chosen from the range of (0:1) and parameter D being dependent on it, thus resulting in a key space of 2 52 × 2 52 = 2 104 .The acquired results for both the tent map and its fractal modification prove that they are both resistant to brute-force attacks, while the logistic map does not allow one to provide a sufficient number of various keys while using a representation of the double floating-point type.
In order to evaluate the PRNGs' bias resistance, one might use a Chi-square formula.Due to algorithms generating data in a binary format, one can calculate the Chi-square for datasets with 1 degree of freedom.In our case, we suggest the null hypothesis that there is no visible correlation between two datasets acquired for different generation parameters [37].If the Chi-square C 2 value is less than the target value of 3.841, then one can state that the null hypothesis is correct and the data are unbiased (Table 3).According to the obtained results, one can state that the bias resistance of the proposed PRNGs might be considered high due to the Chi-square metric being significantly lower than that of a required value.
Figure 9 shows a comparison of Gaussian normal distributions and PRNG distributions of a logistic map, a tent map, and a fractal tent map.It can be seen that the studied PRNGs have normal or nearly normal distributions.To evaluate the dependence between the distribution of PRNGs and different parameters of the underlying chaotic maps, we introduced the parametric histogram plot (Figure 10).In this plot, each parameter value of the chaotic map is represented by a colored vertical line.The color intensity corresponds to the histogram value.It should be noted that, for enhanced clarity, the histogram values are normalized for each parameter value in Figure 10).
To evaluate the dependence between the distribution of PRNGs and different parameters of the underlying chaotic maps, we introduced the parametric histogram plot (Figure 10).In this plot, each parameter value of the chaotic map is represented by a colored vertical line.The color intensity corresponds to the histogram value.It should be noted that, for enhanced clarity, the histogram values are normalized for each parameter value in Figure 10).Figure 10 illustrates that the PRNG distributions for the tent map and the generalized tent map are similar, closely resembling the normal distribution across a broad range of parameter values.However, there are noticeable differences in the distribution width.Conversely, the logistic map exhibits a distinct distribution pattern, with near-zero values Figure 10 illustrates that the PRNG distributions for the tent map and the generalized tent map are similar, closely resembling the normal distribution across a broad range of parameter values.However, there are noticeable differences in the distribution width.Conversely, the logistic map exhibits a distinct distribution pattern, with near-zero values occurring more frequently.These findings suggest that the characteristics of the PRNG are significantly influenced by the specific values of the chaotic map parameters, enabling precise adjustments to the PRNG's behavior.

Surrogate Testing Results
Let us recall that the null hypothesis was that the data being tested were generated by a linear stochastic process with a normal Gaussian distribution.Surrogate data were generated using an autoregressive process of the order p (AR(p) process), defined as follows: where ε t is a normally distributed random value (usually white noise), and a i denotes the parameters of the autoregressive model.The order p of the AR process is determined using the minimum description length (MDL).
To conduct surrogate testing, we use three types of data: 1.
A linear signal that is generated by the AR process.

2.
A nonlinear chaotic signal generated by the logistic map (for k = 4).

3.
A biological signal that is taken from a dataset of 81 recordings of sperm whale clicks.
It is generally a linear signal, but quite noisy.
The results of surrogate testing are shown in Figures 11-13.The histogram shows the distribution of the test statistic for the surrogate data, and the vertical line shows the value of the test statistic for the signal under study.For example, for the first and third types of signals in Figures 11 and 13, we are examining a linear signal, so the test statistic of the signal under study is expected to fit into the distribution of the test statistic of the surrogate data.As we can see in Figures 11-13, in the exemplified runs of the surrogate testing algorithms, all tests give similar, correct results, but during different runs, the random number generator can spontaneously produce an incorrect distribution, which leads to a spurious test result.occurring more frequently.These findings suggest that the characteristics of the PRNG are significantly influenced by the specific values of the chaotic map parameters, enabling precise adjustments to the PRNG's behavior.

Surrogate Testing Results
Let us recall that the null hypothesis was that the data being tested were generated by a linear stochastic process with a normal Gaussian distribution.Surrogate data were generated using an autoregressive process of the order p (AR(p) process), defined as follows: where ε t is a normally distributed random value (usually white noise), and a i denotes the parameters of the autoregressive model.The order p of the AR process is determined using the minimum description length (MDL).
To conduct surrogate testing, we use three types of data: 1.
A linear signal that is generated by the AR process.

2.
A nonlinear chaotic signal generated by the logistic map (for k = 4).

3.
A biological signal that is taken from a dataset of 81 recordings of sperm whale clicks.
It is generally a linear signal, but quite noisy.
The results of surrogate testing are shown in Figures 11-13.The histogram shows the distribution of the test statistic for the surrogate data, and the vertical line shows the value of the test statistic for the signal under study.For example, for the first and third types of signals in Figures 11 and 13, we are examining a linear signal, so the test statistic of the signal under study is expected to fit into the distribution of the test statistic of the surrogate data.As we can see in Figures 11-13, in the exemplified runs of the surrogate testing algorithms, all tests give similar, correct results, but during different runs, the random number generator can spontaneously produce an incorrect distribution, which leads to a spurious test result.For each test, 150 surrogate datasets of length N = 1000 were generated, and each test for four PRNGs was conducted 100 times.Table 4 shows the results of checking the accuracy of determining the signal nonlinearity in terms of the true positive rate.For linear signals, the test statistics of the original signal should fall within the distribution boundary of the test statistics of the surrogate data, and vice versa for a nonlinear signal.As we can see, the best result was obtained when using both the fractal ( D = 1/3, K = 1/3) and common variants of the tent map.Surprisingly, the worst result is achieved when using a standard MATLAB randn() function.number distribution to the normal distribution after applying the Box-Muller transform.The results of NIST tests conducted for the proposed PRNG also show high randomness.The proposed generator occasionally outperforms the MATLAB built-in randn() function, as well as PRNGs based on the logistic map and the conventional tent map.The PRNG based on the proposed chaotic map can be used in stream encryption, chaos-based communications, and other real-world applications that require pseudorandom sequences with a high level of randomness.Using a chaotic-map-based PRNG for surrogate testing shows room for improving the accuracy of rejecting or confirming a nonlinearity hypothesis compared to a standard MATLAB function for generating normally distributed random values.Table 4 shows that the standard MATLAB function randn() makes more errors in determining a linear signal than the proposed PRNG.Errors have also been observed for the PRNG based on the logistic map when determining the linearity of a biological signal with noise.
The obtained results can inspire further studies of fractal chaotic maps and flows from both theoretical and applied points of view.

Figure 1 .
Figure 1.Graphs of the logistic map function (on the left) and the tent map function (on the right).

Figure 1 .
Figure 1.Graphs of the logistic map function (on the left) and the tent map function (on the right).

3 Figure 3 .
Figure 3. Phase portraits of the fractal tent map in Equation (3) for various values of K and D. Panels (a-c) correspond to non-fractal behavior, and panels (d,e) correspond to fractal behavior.Decreasing D leads to smaller intervals, and decreasing K leads to more rapid interval length decay.Fractal behavior occurs if 1 − K = 2D, K ∈ (0, 1].

Figure 3 .
Figure 3. Phase portraits of the fractal tent map in Equation (3) for various values of K and D. Panels (a-c) correspond to non-fractal behavior, and panels (d,e) correspond to fractal behavior.Decreasing D leads to smaller intervals, and decreasing K leads to more rapid interval length decay.Fractal behavior occurs if 1 − K = 2D, K ∈ (0, 1].

Figure 4 .
Figure 4. Relative density-based colored one-dimensional bifurcation diagram for the logistic map.

Figure 5 .
Figure 5. Relative density-based colored one-dimensional bifurcation diagram for the tent map.

Figure 4 .
Figure 4. Relative density-based colored one-dimensional bifurcation diagram for the logistic map.

Figure 4 .
Figure 4. Relative density-based colored one-dimensional bifurcation diagram for the logistic map.

Figure 5 .
Figure 5. Relative density-based colored one-dimensional bifurcation diagram for the tent map.

Figure 5 .
Figure 5. Relative density-based colored one-dimensional bifurcation diagram for the tent map.

Figure 6 .
Figure 6.Relative density-based colored one-dimensional bifurcation diagram for the generalized tent map.Points K = 1 − 2D, where the condition of the fractal function is satisfied, relate to higher orbit densities.Nevertheless, non-fractal parameter combinations can also produce similar densities at some points.

Figure 6 .
Figure 6.Relative density-based colored one-dimensional bifurcation diagram for the generalized tent map.Points K = 1 − 2D, where the condition of the fractal function is satisfied, relate to higher orbit densities.Nevertheless, non-fractal parameter combinations can also produce similar densities at some points.

Figure 7 .
Figure 7. Approximate and information entropy for fractal tent map with parameter K = 1.0.

Figure 8 .
Figure 8  gives the values of the largest Lyapunov exponents (LLEs) of the considered chaotic maps.One can see an increase in the LLE in Figure8a, corresponding to the fractal tent map when D tends to zero, which relates to greater mixing properties due to the growing number of linear fragments in the function F f (x).Figure8b, corresponding to the fractal tent map in the fractal regime, shows the same behavior.The two other panels illustrate LLEs depending on parameters for the logistic and tent maps.Notably, in the generalized tent map, high values of LLEs are much easier to obtain.

Figure 7 .
Figure 7. Approximate and information entropy for fractal tent map with parameter K = 1.0.

Figure 7 .Figure 8
Figure 7. Approximate and information entropy for fractal tent map with parametFigure8gives the values of the largest Lyapunov exponents (LLEs) of chaotic maps.One can see an increase in the LLE in Figure8a, correspondi tent map when D tends to zero, which relates to greater mixing prope growing number of linear fragments in the function F f (x).Figure8b, co the fractal tent map in the fractal regime, shows the same behavior.The tw illustrate LLEs depending on parameters for the logistic and tent maps.generalized tent map, high values of LLEs are much easier to obtain.
Figure 8  gives the values of the largest Lyapunov exponents (LLEs) of chaotic maps.One can see an increase in the LLE in Figure8a, correspondi tent map when D tends to zero, which relates to greater mixing prope growing number of linear fragments in the function F f (x).Figure8b, co the fractal tent map in the fractal regime, shows the same behavior.The tw illustrate LLEs depending on parameters for the logistic and tent maps.generalized tent map, high values of LLEs are much easier to obtain.

Figure 8 .
Figure 8. Dependence of largest Lyapunov exponents on parameters of investigate

Figure 8 .
Figure 8. Dependence of largest Lyapunov exponents on parameters of investigated chaotic maps.

Figure 9 .
Figure 9.Comparison of distributions of sequences generated by PRNG based on logistic map, fractal chaotic map, and tent map with Gaussian normal distribution.

Figure 9 .
Figure 9.Comparison of distributions of sequences generated by PRNG based on logistic map, fractal chaotic map, and tent map with Gaussian normal distribution.

Figure 10 .
Figure 10.Relative density-based two-dimensional histograms showing the dependence of the PRNGs' distributions on map parameters.

Figure 10 .
Figure 10.Relative density-based two-dimensional histograms showing the dependence of the PRNGs' distributions on map parameters.

Figure 11 .
Figure 11.A histogram of distributionsof test statistics for a linear signal.Figure 11.A histogram of distributionsof statistics for linear signal.

Figure 11 .
Figure 11.A histogram of distributionsof test statistics for a linear signal.Figure 11.A histogram of distributionsof statistics for linear signal.
al. present a simple method for designing

Table 1 .
Results of NIST testing for PRNG based on different maps.

Table 2 .
Key-space estimation for every proposed algorithm.

Table 3 .
Chi-squarecalculation for datasets acquired from PRNG based on different parameter values.