3.1. Study Scenario Description
At first, we will consider that the value of the external message is zero and the inclusion loop is not taken into account (in
Figure 2, switch
K is in position B). This is denoted as the
reference scenario because in this case the scheme represents the dynamical system described by the logistic map without any other intervention. The results of this scenario (outputs/visualization points marked on the scheme with 1 and 2) will be used for comparison in a
study scenario, that will consider the dynamical system affected by inclusion.
For the study scenario we consider the external message as a byte session (image or text) and the inclusion loop is taken into account (i.e., the switch works in positions A or B). Thus, the inclusion signal will be formed by summing the discretized chaotic signal with the external message. These resulting values will be included in the dynamical system evolution after scaling them by a given F factor.
We will compare the results of this
study scenario with the previous scenario, the reference one, in the visualization points (1,2,3) marked on the diagram in
Figure 2.
3.2. Study Method
In visualization point 1, the random process associated with the logistic map can be followed; specifically, we follow the trajectories (the particular realizations of the random process). Each trajectory is determined by the initial condition and the
R control parameter value. It is known that this random process is ergodic [
12,
13]. In our study, we consider
R = 4; in this case, the first-order probability law of the random process is described by the probability density function in
Figure 1a. The probability law in
Figure 1a is considered in the stationarity region of the random process (after the transient time has elapsed) [
14]. We follow whether the first-order probability law of the chaotic system is affected by inclusion. Therefore, we analyze the random process at observation point 1 of the scheme in the scenario with inclusion and will refer to the
reference scenario. For this we apply the Kolmogorov–Smirnov test.
The Kolmogorov–Smirnov statistical test verifies the concordance between an experimental distribution law and a theoretical one. A detailed description of the steps of the Kolmogorov–Smirnov test is presented in [
14]. The experimental data set used for the Kolmogorov–Smirnov test consists of the values (
,
,
, …,
) obtained by sampling
N trajectories of the random process of the new system (with message inclusion) at iteration
k. To obtain the trajectories, we considered
N initial conditions
,
,
, …,
; randomly generated according to the uniform distribution law in [0, 1] and
R = 4, iterating the dynamical system until the chosen
k iteration. In this paper, we used
N = 10,000 trajectories,
k = 150, and the statistical significance level of the test, α = 0.05. To choose the
k value, we considered an iteration in the stationary region of the logistic map [
14]. Using the experimental data set sampled from the random process at
k iteration, we check whether the experimental probability law verifies the theoretical law.
The two hypotheses of the Kolmogorov–Smirnov statistical test are:
H0: the experimental data set comes from the same theoretical distribution as in (3), the inclusion does not modify the probability law.
H1: the experimental data set does not come from the same theoretical distribution as in (3), the inclusion affects the probability law.
The maximum absolute deviation between the two distribution functions, theoretical
and experimental
, is calculated. This becomes the
δ test value.
The
δ test value is compared to the
threshold value, according to the chosen
α value:
where
N = data volume (number of trajectories considered) and
α = significance level of the test.
If δ ≤ , H0 hypothesis is accepted.
According to the estimation theory [
16], for
α = 0.05 and using Monte Carlo analysis by resuming the statistical test for
L = 500 times, the range of accepted proportions is [0.93, 0.97].
Note: To determine the range of accepted proportions, we use the confidence interval defined as , where ; L = 500; is point value of the standard gaussian law. In our case = 1.96 for = 0.05. Thus, the accepted region will be ; namely [0.93, 0.97].
In observation point 2, we acquire the data representing the trajectory of the chaotic system after discretization. This trajectory appears as a sequence of bytes, and we represent it as an image. We also display the corresponding histogram.
In observation point 3, the output image of the encryption scheme and its histogram are visually inspected. For illustration, we used black and white images of size 256 × 256 pixels, each pixel being represented by 8 bits. Thus, the corresponding histograms are made in 256 values.
The processing scheme and simulations were implemented using the Matlab R2017b development environment.
3.3. Experimental Results
In order to perform the simulations corresponding to the two scenarios (
reference and
study scenario) we considered for the logistic map the control parameter
R = 4 and the initial condition
= 0.2. For the scaling factor
F in the scheme in
Figure 2 we have chosen the value
F =
based on the study made in
Section 3.4.
In the
reference scenario, the scheme in
Figure 2 is in fact the dynamical system logistic map working according to relation (1) (unaffected by inclusion). Thus, the inclusion loop is not taken into account, the decision switch
K is always in position B (open). We illustrate graphically some elements that will matter in discussing the effect of inclusion in the
study scenario.
By plotting the trajectory of the dynamical system at observation point 1 in
Figure 2, for two different initial conditions and
R = 4, we obtain two distinct trajectories; illustrated in
Figure 3.
Figure 4 shows the theoretical first order distribution function and the experimental distribution function of the random variable sampled at iteration
k = 150, obtained using the method described in
Section 3.2. A Monte Carlo analysis was performed, the Kolmogorov–Smirnov test being repeated 500 times. The obtained percentage of
H0 hypothesis acceptance was 94.8%, which is in the estimation range [0.93, 0.97]. Therefore, we can state with 95% statistical confidence that the experimental data come from the theoretical probability law of the logistic map for
R = 4.
Figure 5 shows the image after discretizing the dynamical system values and its corresponding histogram, obtained at the observation point 2. Analyzing the obtained histogram, it can be noticed that it corresponds to the histogram of the probability density of the logistic function from
Figure 1a. This is a result of the ergodicity of the random process because the image is obtained from a temporal analysis of the measurements at visualization point 2, tracked over time, over a large number of iterations, on a random chosen trajectory.
Next, we analyze the behavior of the dynamical system in the
study scenario by comparison with the
reference scenario. To illustrate the
study scenario, we consider the external message to be an image of size 256 × 256 pixels, with low entropy; this is illustrated in
Figure 6.
The first analysis step is the comparison of the evolution of the dynamical system trajectories in the
study scenario with the
reference scenario at observation point 1, as shown in
Figure 7. The trajectories were obtained using
= 0.2,
R = 4 and
F =
.
It can be observed that the trajectories are different. This shows that inclusion leads to a new system with modified dynamics. Further ,we check whether the properties of the new system follow the statistical behavior of the reference system described by the logistic map.
A first step in verifying that the properties are preserved is to perform the Kolmogorov–Smirnov test in the
study scenario for the data sets obtained at observation point 1, using the method described in
Section 3.2. As can be seen in
Figure 8, the experimental distribution function of the new dynamical system corresponds to the theoretical one. A Monte Carlo analysis was performed, the Kolmogorov–Smirnov test was repeated 500 times. The obtained percentage of acceptance of hypothesis
H0 was 95.4%.
Next, we verify in observation point 2 the image obtained after discretizing the values of the new dynamical system and the corresponding histogram, illustrated in
Figure 9. Analyzing the histogram in
Figure 9 by comparison with the histogram obtained for the
reference scenario in
Figure 5, it can be stated that the experimental probability density of the new dynamical system is also similar to the theoretical probability density of the logistic map.
Adding now observation point 3, for the
study scenario we can illustrate the use of the dynamical system for encryption when the input message is introduced in its dynamics (encryption by inclusion). The result is illustrated in
Figure 10 and shows that a visual transformation of the input image is achieved, but with low diffusion. A similar processing was performed in [
7], the difference in the operation of the encryption scheme in the current paper being the introduction of the decision switch. Hence, the addition of this switch does not influence the performance of the encryption, but is intended to ensure the possibility of inclusion and avoid situations where the system exceeds the range of values of the logistic function [0, 1]. In
Section 4, a cryptographic improvement will be made by adding mixing functions [
7] to the processing scheme in
Figure 2.
In conclusion, the statistical behavior of the new system obtained by message inclusion respects the 1st order probability law of the random process associated to the logistic map. Both the addition of the decision switch in the processing scheme and a good choice of the scaling factor contribute to this fact. In order to determine the impact of the scaling factor on the statistical behavior of the new dynamical system, a study is performed in the next section.
3.4. Determining the Upper Limit of the Magnitude of the Scaling Factor
To include the message in the evolution of the dynamical system without affecting its behavior, message scaling is required. In previous studies [
5,
6] the choice of scaling factor was done empirically, and no upper limit was given for it. We propose to determine the limits of the magnitude of the scaling factor. We expect the magnitude of the scaling factor to be small enough to allow the scheme to operate within the range of definition of the logistic function and not to change its statistical properties. No minimum value is given for the lower limit. As the value of the scaling factor tends to 0, the inclusion of the message no longer influences the evolution of the dynamical system.
For this purpose, we applied the Kolmogorov–Smirnov statistical test using the experimental data from observation point 1 for a set of
F values in the
study scenario where we do not have external message. The experimental data set was obtained under the same conditions described in the presentation of the Kolmogorov–Smirnov test in
Section 3.2. To decide on the results, we performed a Monte Carlo analysis by repeating the Kolmogorov–Smirnov test 500 times for each value of
F. Considering the significance level of the test α = 0.05, the acceptance rate of the test is [0.93, 0.97].
Analyzing the results obtained in
Table 1, it can be stated that the experimentally determined upper limit of the scaling factor is
, since for values less than or equal to this value, the Kolmogorov–Smirnov test had an acceptance proportion of the
H0 hypothesis in the range [0.93, 0.97].
Therefore, considering that the acceptance proportion of the H0 hypothesis is in the range [0.93, 0.97] for values of F ≤ , it can be concluded that simply using the decision switch alone—without a proper choice of the magnitude of the scaling factor—is not sufficient to not disturb the statistical behavior of the dynamical system.
Thus, this study contributes to a correct choice of the scaling factor that, together with the use of the decision switch, contributes to the proper functionality of the encryption scheme based on the logistic map for R = 4.