Abstract
Differential Evolution (abbreviation for DE) is showing many advantages in solving optimization problems, such as fast convergence, strong robustness, and so on. However, when DE faces a complex target space, the diversity of its population will degenerate in a small scope; even sometimes it is premature to fall into the local minimum. All things contend in beauty in the world; a Shuffled Frog Leaping Algorithm (abbreviation for SFLA) has a strong global ability; unfortunately, its convergence speed is also slow. In order to overcome the shortcoming, this article suggests a Shuffled Frog-leaping Differential Evolution (abbreviation for SFDE) algorithm in a cognitive radio network, which combines Differential Evolution with Shuffled Frog Leaping Algorithm. This proposed method hikes its local searching for a certain number of subgroups, and their individuals join together and share their mutual information among different subgroups, which improves the population diversity and achieves the purpose of fast global search during the whole Differential Evolution. The SFDE is examined by 20 well-known numerical benchmark functions, and those obtained results are compared with four other related algorithms. The experimental simulation in solving the problem of effective throughput optimization for cognitive users shows that the proposed SFDE is effective.
1. Introduction
The convergence of the basic Differential Evolution (abbreviation for DE) is closely related to its control parameters. Many researchers have worked on improving its performance in various ways and developed many variants. In [1], a method of Self-adaptive Differential Evolution algorithm with discrete mutation control parameters (DMPSADE) is proposed. In DMPSADE, each vector has its own mutation, cross parameters, and control strategies. Those original unified control parameters and mutation strategies are divided into each individual, so that the evolutionary granularity becomes smaller and the performance improves. In [2], a fuzzy system is introduced to dynamically adaptive control parameters to improve DE. An improved DE named ADE is proposed in [3], where the author suggests DE with two-level parameter adaptation. The first level of control parameters FP and CRP is used to control the global optimization ability, and the second level of control parameters Fi and CRi corresponds with each individual to improve the local search of the algorithm. The main characteristic of DE is its mutation strategy. Price and Storn proposed more than ten different differential strategies to implement mutation operation such as DE/rand/1/bin, DE/rand/2/bin, DE/best/1/bin, DE/rand-to-best/bin, and so on. Among them, DE/rand/1/bin and DE/best/2/bin are the most widely used and most successful differential strategies. The former can maintain the diversity of the population, and the latter pays more attention to the speed of its convergence.
Reference [4] introduces a fitness Euclidean-distance ratio for multimodal optimization. In [5], a novel mutation strategy, named elite/rand/1, which divides a whole population into two subpopulations on the basis of the fitness and then extracts the maximum information from an elite individual. MDE in [6] includes three new steps: Opposition-Based Learning (OBL), tournament method for mutation, and a single population structure. A parameter adaptive selection can dynamically evaluate their performance during evolution and guide the direction of its development in the next iteration process [7]. Based on the abstract convex theory, a dynamic adaptive differential evolution algorithm (DADE) is proposed in [8]. Shuffled Frog Leaping Algorithm (SFLA) is proposed in [9], which is inspired by the individuals learning from each other during the foraging behaviour of frogs. In [10], a multiphase hybrid algorithm implements a local depth search using SFLA for every cluster and readjusts its global solutions. An improved hybrid bifurcation algorithm based on cloud model is proposed in [11], which focuses on the transformation between qualitative and quantitative computing. Reference [12] proposed crossing and variation frog leaping algorithm (KSFLA), where the individuals of the subpopulations vary from the ranking list before producing new individuals instead of the poor ones. Reference [13] proposed a Self-adaptive Differential Evolution algorithm with Improved Mutation Mode (IMMSADE), which improves the mutation model of DE and introduces some new control parameters.
Based on the above discussion, SFLA has slow the convergence, low efficiency, and precision in the high dimension continuous optimization problems, in that its core strategy is to continuously update the location of the worst individual. But DE does not do anything with the worst individual. A hybrid algorithm based on frog leaping algorithm and Differential Evolution is suggested in this article, which combines the advantages of SFLA and DE, in order to obtain a significant performance on convergence speed and the accuracy.
The rest of this paper is organized as follows. In order to ensure the readability, the related mathematics description of Differential Evolution and Shuffled Frog Leaping Algorithm is summarized in Section 2. Section 3 formally defines models of an improvement in a Shuffled Frog-leaping Differential Evolutionary algorithm (SFDE). Section 4 designs the comparative experiments of well-known public benchmark functions and analyses the related T-test and F-test results for proof of the effective SFDE. The analysis of convergence in SFDE is described in detail in Section 4. Simulation experiments on the problem of effective throughput optimization for cognitive users are given in Section 5. Section 6 concludes this article with a summary and future direction.
2. The Related Work
2.1. About Differential Evolution
DE is a random search based on population and also a heuristic swarm intelligence algorithm. It mainly uses mutation, crossover, and selection operations to perform an intelligent evolutionary search. DE has five operations, which are described as follows.
(1) Initialization Operation. Initializing a population randomly in the search space, the j-th parameter of the i-th individual in the 0-th generation is initialized by formula (1):In formula (1), xUp and xLow are the upper and lower boundary of the j-th parameter of the i-th individual, respectively. denotes a random number, which is limited between 0 to 1.
(2) Mutation Operation. In DE, the most common in mutation operations is DE/rand/1, which selects two individuals at random in a population. Then the vectors are combined with the target vector for mutating into an intermediate vector . The process can be expressed as follows in formula (2):Obviously, the smaller the difference between and is, the smaller the difference between and is. This means that, in the initial stage, the disturbances are large because of the far distance between individuals, and the search range of the algorithm is large. But in the later period of the iteration, the detective range will be a shrink. Several common mutation strategies are shown in the following formulas (3), (4), and (5):① DE/rand/2② DE/current-to-best/1③ DE/best/1Among them, the subscript variables , , , , and represent different individuals and represents the optimal vector in the -th generation.
(3) Correcting Operation. After the mutation operation, the algorithm must ensure that each component (gene) of the intermediate vector satisfies the predefined boundary constraints. If a component (gene) of an intermediate vector exceeds the scope, the vector individual must be corrected. The following formula (6) is a simple but most frequently used correction strategy:
(4) Crossover. Crossover is to mixture the and in the -th generation, producing test vectors. The process can be expressed as follows in formula (7):In formula (7), CR is the cross probability and is a random integer of .
(5) Selecting Operation. By comparing the function values of the target vectors in the test function with the test function values of the corresponding test vectors , which is determined whether the test vectors can be put into the next generation or not, the standard Differential Evolution algorithm employs a greedy algorithm for the selection operation. The process can be expressed as follows in formula (8):
2.2. Shuffled Frog Leaping Algorithm
Suppose that there are frogs in the n-dimensional space catching their food. The basic Shuffled Frog Leaping Algorithm sorts all the individuals in descending order according to the fitness of the objective function and then puts the first frog individual into the first group, the second frog individual into the second group, and so forth. The first frogs entry own group in order, respectively; thus each group contains a frog. Beginning with the frog, the remaining frogs join their own group according to the rule of Round Robin (namely, Mod m).
In each group, represents the current position of the frog. is the optimal location in the group, called the local optimal position. is the worst in the group, called the local worst position.
Let be a minimal objective function, and use formula (9) to calculate the local optimal frog location for the group.Assume that the number of the frogs is , and it is divided into groups according to the grouping operator. The optimal position for each group is the global optimal location, as shown in formula (10).According to the local location update operator, adjusting the position of the worst frog in each group, the concrete adjustment method is as follows.
The frog moving distance is as shown in formula (11):Update the worst frog individuals as follows:In formula (11), is a random number between 0 and 1. in formula (12) represents the maximum distance that the frog can move during its foraging for food. If the positions obtained from formulas (11) and (12) are closer to the food source than the position of the worst frog in the group, the worst frog leaps to this new position. The position of the worst frog is further closer to the food source. Otherwise, is replaced by , and the worst frog continues to search for a new leaping position as follows:If the position obtained from formulas (13) and (14) is still not closer to the food source than the worst frog in the group, or the worst frog individual needs to skip more than the maximum step size, the worst frog leaps to a new position at random. Otherwise, the worst frog leaps to the new position generated from formula (14). SFLA performs an iterative evolution of the frogs within each group via adjusting and updating of the worst frog.
In the evolutionary process of SFLA, each iteration involves the local optimal solution and the global optimal . In the whole population, there exist groups, and therefore there are local optimal solutions, and there is only one global optimal. SFLA uses a grouping operator and a mechanism of individual integration into groups for transmitting messages and sharing information. The grouping operator divides all the frogs in the population into groups, and then the worst frog in each group updates itself local search strategy. When all groups finish searching, a shuffle share information and interact with each other. This mechanism makes SFLA have better global search ability. The procedure of basic SFLA is shown as follows.
Step 1. Initialize the positions of all the frogs in the population at random, to make the frogs randomly distributed in the search space, and initialize the population size , the number of frogs in a group , the number of evolution of the , the maximum step size of the frogs allowed to move , and maximum number of iterations of population .
Step 2. Calculate the fitness of all frogs in the population and sorting them in descending order. is initialized with the best frog position, and the frogs are sorted into groups according to a grouping operator.
Step 3. Adjust the worst frog in the group according to formulas (11) and (12).
Step 4. If the position of the frog in Step 3 is improved, replace the worst frog position in the group with the new position, then skip to Step 6, or readjust the worst frog position of the group according to formulas (13) and (14).
Step 5. If a position of the frog from Step 4 is improved, replace the worst frog position in the group with the new position, then skip to Step 6, or generate a new random position instead of the worst frog position in the current group.
Step 6. Calculate the fitness of all frogs in the group, and sort them in descending order to determine if the internal iteration is over. If not, skip to Step 3.
Step 7. Share shuffle information among all the frogs in the same group.
Step 8. Judge whether the end condition meets or not. If it does, output the optimal solution directly. If not, skip to Step 2 and it continues.
The end condition of SFLA can achieve the maximum number of iterations or meet an error criterion, which may be a combination of the above two. In the case where the optimal value of the objective function and the allowable error are determined in advance, the end condition can be set to satisfy a certain error range. When the difference between the objective value and the optimum is less than the setting error, output the optimal value; otherwise it continues.
3. The Shuffled Frog-Leaping Differential Evolutionary (SFDE)
Although combining a new algorithm with simple assembly can improve the accuracy or convergence of the algorithm to a certain degree, it is still far from the theoretical value of the test function, so how to integrate the two algorithms is the most important problem.
3.1. Integration Strategy between SFLA and DE
There are three main steps in the integration strategy between SFLA and DE. Step : initialize populations. Step : divide the whole population into several groups, and integrate within a group. Step : shuffle the whole population after iterations in each group. We use the idea of dividing the whole population into several subpopulations and considering each subpopulation as a group. So before dividing the whole population into groups, individuals are sorted by the fitness value of the test function, then a new group with the strategy is shown as in Figure 1.

As shown in Figure 1, SFLA helps to adjust the worst individuals in the groups. If the fitness value of the individuals generated from formulas (11) and (12) is improved, the worst one in the original group will be replaced. Otherwise, its position of the group should be readjusted according to formulas (13) and (14). If it is worse than original one; then discard it and randomly generate a new one to replace it. Then we use DE to deal with the new population. The operations include mutation, crossover, and selection.
The evolutionary process of the individual population is shown in Figure 2, where the sequence (1, 2, 3, 4) is the stage number of each group evolution. The worst individual position is adjusted for each subpopulation according to Figure 1, then the operation of the Differential Evolution for all the individuals in the subpopulation, that is, the variation of the evolution and the selection operation. The whole population run mutation and cross operation of Differential Evolution stage, thus resulting in the middle of the vector of individuals (represented in Figure 2 with a red circle of individuals) for the latter part of the selection operation. Select the operation for the generation of intermediate individuals and original target individuals.

In Figure 2, the next generation has the intermediate vector individuals (individual with red virtual coil) and the original target individuals. In SFDE, the mutation strategy uses DE/best/1/bin.
3.2. The Basic Process of Shuffled Frog-Leaping Differential Evolutionary (SFDE)
The basic SFDE can be described in Algorithm 1.
|
In Algorithm 1, T is the number of subpopulation; denotes the best individual in population; is the worst individual in subpopulation; pbesti is the best individual in subpopulation; denotes the current mutated vectors.
4. Experimentation
This section compares the proposed SFDE with twenty classic well-known test functions.
4.1. Experimental Setting and Parameterization
Twenty test functions with 30 high dimensions check the performance of SFDE. The detailed definitions of them are shown in Table 1.
Experimental environment configuration: operation system is Windows 7. Minimum memory is 4G; processor type is Intel (R) Core (TM) i5-2370M @ 2.40GHz; development tools and version are Matlab-R2012a.
To be fair, the initial conditions of each algorithm are consistent. Population size of all Differential Evolution algorithms is set to be 200. The number of subpopulations of SFDE is set to 10. The number of iterations is set according to the complexity of the problem, but the number of iterations of all Differential Evolution algorithms is equal every time.
4.2. Computational Results and Discussion
We record the optimal value obtained at the end of every computing in 30 runs. Then, the average optimal value and the standard deviation of the optimal value are calculated by using the 30 dependently runs. And the standard deviation of the best average value and the optimal value is inputted into the table to facilitate the analysis of the experimental results. In order to highlight the experimental data close to the theoretical value of the function, we test each test function in the table closest to the theoretical value of the experimental data bold. The output images reflect the variation curve of the average optimal value for each test function in 30 separate experiments.
The convergence curves plots of five algorithms for 9 benchmark functions with 30 dimensions are drawn in Figure 3 (Ackle, Michalwicz Function and Schwefel’s Problem 1.2), Figure 4 (Perm Function, Rastrigin, Sphere Function), and Figure 5 (Rotated Hyper-Ellipsoid, Shifted and Zakharov Function, Dixon-Price Function), respectively. The 9 results confirm the good performance of SFDE. The optimal value of the function curve is shown in Figure 3 and the full optimal value of the function is shown in Table 2.

(a) Ackle Function

(b) Michalwicz Function

(c) Schwefel’s Problem 1.2

(a) Rastrigin Function

(b) Perm Function

(c) Sphere Function

(a) Rotated Hyper-Ellipsoid

(b) Zakharov Function

(c) Dixon-Price Function
Schwefel’s Problem 1.2 function’s algorithm is difficult to obtain the theoretical value on 30 dimensions (iterations: 2000 times). From the experimental results, it can be seen that the convergence rate of SFDE is faster than that of other algorithms, and the final accuracy is also better than other algorithms, and even 34 units of magnitude are higher than the accuracy of EPSDE algorithm.
Rotated Hyper-Ellipsoid is a continuous, convex single-peak, which is an extension of the axis-parallel super-ellipsoid function and also becomes a square sum function. From the experimental results, SFDE is better than other algorithms, and with the iteration of the algorithm, the accuracy of the algorithm is getting better.
The Perm Function is characterized by a large number of local optimal value of the single-peak function, and it is a well test function which can test the function of the algorithm. It can be seen from the experimental results that SFDE has the highest precision and the fastest convergence rate for Zakharov Function with 30 dimensions (iterations: 2000 times).
The convergence curves of five algorithms for 9 benchmark functions with 30 dimensions are drawn in Figure 3, respectively. For each benchmark function, there are two curves, the left one is plotted with iterations, and average function values are plotted with NFEs at right. Figures 3(a)–3(c) and 4(a) present multimodal benchmark functions and Figures 4(b)-4(c), 5(a)–5(c) present unimodal benchmark functions. Figure 3(a) is the curves of average optimal of F1. In Figure 3(b), SFDE and CODE almost have the same best value on F2. The convergence curves of compared algorithms for F3 function and F13 function are given in Figures 3(c) and 4(a). Figure 4(b) present convergence curves of test results for F10 function. In Figure 4(c), 5(a)-5(b), plots of SFDE are much better than other variants for F4, F6, and F8 functions. In Figure 5(c), SFDE demonstrates best average and convergence speed for 30 dimensions on F9 function and MDE is the second best on F9.
4.3. Results Analysis
Table 2 shows the comparison between SFDE and classic DE. It can be seen from the above experimental results that SFDE does better in most of the base functions compared with other differential evolutionary algorithms, and the convergence speed is also accelerated.
In Table 2, the eight test functions, , and , are unimodal functions, and the solving accuracy of SFDE is better than that of other algorithms. In the results of the multimodal test function, the individual test function is poor performance, such as and . However, the solving accuracy of SFDE is much higher than that of other algorithms and SFDE is more stable. More importantly, SFDE algorithm itself is very simple and fast and has an absolute advantage.
In order to determine whether SFDE is more effective than other methods, statistical methods need to detect the results of CPU running time; this paper selected the paired T-test and the paired F-test. The T-test is an important test of the mean difference between the two samples, with the aim of checking the system error. The F-test is a significant test of the variance of the two samples, with the aim of testing the accident error. Usually we compare the value of 0.05. If the test result is less than 0.05, the cell will mark the result with ‘+’, instead, the cell will mark the result with ‘-’. The percentage of ‘+’ in all results is the likelihood that SFDE will handle the reliability and advantage of the test problem. Specific test results as shown in Table 3, each column is test function, efficiency, detection methods, and a variety of intelligent optimization.
To compare SFDE with other algorithms, we select advanced variants of DE, namely, JADE, DMPSADE, and IMMSADE, which is introduced in [12]. The comparison results are as shown as follows in Table 4.
Table 4 shows a clear comparison of SFDE-SaDE, SFDE-JADE, SFDE-DMPSADE, and SFDE-IMMASDE. From Table 4, it is clear to see that the convergence speed of SFDE is faster among the considered algorithms for most of test functions.
5. Simulation of Spectrum Throughput Optimization
There are three ways of spectrum sensing: energy detection, matched filter detection, and cyclone-stationary property detection, the most commonly used of which is energy detection. In this way, if the length of each frame is fixed, its sensing time synchronizes with its energy accumulation. At the same time, if authorized users have good protection, the false alarm probability will be smaller and the detection rate will be larger, but the data transmission time will be shortened [14]. That is to say, the idle time of the cognitive user’s usage spectrum becomes shorter, resulting in lower throughput of the cognitive user. Therefore, how to arrange an effective schedule to maximize the throughput of cognitive users (i.e., authorized users have a right of first class service) is a hot research topic.
At present, Monte Carlo is a common method, but its implementation is a complex and time-consuming task [15], and the more important point we have noticed is that it cannot also meet the real-time requirements. It is high time to improve the sensing time by using the improved Differential Evolution algorithm, such as good searching performance, fewer parameters, and faster convergence rate, so as to maximize the system throughput.
5.1. Mathematical Model Based on Energy Detection
Without interfering with the authorized users, it is necessary to determine if they are busy or idle by continuously detecting signals received in a certain frequency band, so as to determine whether these resources are accessible or not, which is a binary hypothesis testing problem, we assume that the authorized user accesses his reserved frequency for periodic , and its disappearance takes as a cycle; the energy detection model is as formula (15):In the above model, is the signal received by the cognitive user. is a Gaussian noise with a mean of 0 and a variance of . is the signal belonged to authorized user, the mean value is 0, and the variance is .
Spectrum sensing concerns two parameters (i.e., detection probability and false alarm probability) [16]. Let be the duration of observation, and while the statistic index denotes as formula (16)where N is the amount of those perceived samples and f is the sampling frequency, it is obvious that . Suppose that probability density is a random variable with distribution of and a given threshold ; the false alarm probability can be given by the following formula (17):where refers to the probability of misjudging an authorized user without his/her perceiving, is the variance of , and is the mean of , respectively.
However, assuming the given and a threshold , the probability density will be a noncentral distribution variable. Then the corresponding detection probability is shown in formula (18).where refers to the probability of correctly judging the existence of authorization and is the signal-to-noise ratio from cognitive users to authorized users. From formulas (17) and (18) combined, we infer that formulas (19) and (20) are correct.When the timely requirement of an authorized user is detected, the related cognitive user needs to quickly give a way to the authorized user. Based on the above formulas (19) and (20), it is obvious that the higher the value of is, the better the related protection is obtained. The lower the is, the higher the reusability of wireless resource is. So a good algorithm should have high and low .
5.2. A Trade-Off between User Throughput and Perceived Time
The time cost of processing one data frame of cognitive radio includes two parts, i.e., frame perceived phase and data transmission phase. When the first perception becomes longer, the detection probability becomes larger, the probability of false alarm becomes smaller, but the cognitive user’s data transmission time will be reduced [17]. Therefore, how to balance the sensing and transmission time under the precondition of authorizing users to operate safely and then improve the throughput of the cognitive radio system is an urgent problem to be solved. In this section, the periodic structure of the radio frame is the basis for research, as shown in Figure 6.

When the authorized user is not in the channel, there is no false alarm probability [18]; the cognitive user’s throughput is ; if an authorized user exists, the cognitive user throughput is . Obviously C0 > C1, and , where Ps is the average cognitive user power and is the noise power. Then , where is the interference degree of the authorized user at the cognitive user acceptance point [19]. It can be found that, in order to protect authorized users, the probability of is close to 1. At this time the system throughput is as formula (21):Among them, , and denote the probability that the first user will not appear and appear, respectively, and . Since C0 > C1 and approximately equal 1, is close to zero. Then formula (21) can simplify the following formula (22):Since the Q is a decreasing function, when the sensing time becomes larger, under the specified detection probability, the false alarm probability will be small, and the data transmission time will be smaller. That is to say, the value of formula (22) may become larger or smaller. The trade-off between perceived user throughput and perceived time is to find the optimal sensing time τ under the protection of the authorized user, so that the system arrives at the maximum effective throughput. Therefore, the optimization problem of looking for the optimal perception time is transformed into formula (23):To solve problem of the effective throughput of spectrum resource in wireless cognition network, we need to find the best service mode for the personalized users. This is a problem of function optimization. This section uses the SFDE to find the optimal perception time.
5.3. Specific Solution Based on SFDE
Specific implementation steps are as follows, which is based on SFDE.
Step initializes a population. The entire range of variables is from 0.001 to 1 into 45 intervals. The upper limit of the i-th segment is bounds and the lower bound is limited to bounds . The specific evaluation process is shown as formula (24):where each interval assumes as an individual and the whole has 45 individual vectors. The whole population divide themselves into the three subpopulations, each of which has 15 individuals. Then each individual is initialized as shown in formula (25):where the 0.0002 is the length of itself individual gene in the population.
Step finds the optimal individual Xlbest in each subpopulation. Each gene vector is brought into fitness formula (22). The Xlbest is based on the fitness and the maximum value of each gene for each individual. With randomly selecting the parameters F and CR for each individual from the parameter candidate pool, we can get the sniffing values F_ad and F_de of the parameter F corresponding to each individual are obtained and use the optimal individual to compute the variation vector.
Step chooses sniffing values CR_ad and CR_de of the parameter CR corresponding to each individual. And we can combine the intermediate vector with the mutation vector of Step .
Step lets the intermediate vector obtained in Step (3) participate in the selection operation, and we can compute the fitness values for each gene in each intermediate vector by formula (22) and compare the fitness value with the target vector. If the fitness is greater than the target vector, the gene will replace the gene corresponding to the target individual. Repeat Steps until all individuals in the population have completed the operation of mutation, crossover, and selection. Then the new population goes into the next generation.
Step repeats Steps until the maximum number of iterations is reached. The optimal sensing time and the corresponding maximum effective throughput are obtained.
5.4. Parameters Setting
In order to prove the effectiveness of the SFDE algorithm in the process of system effective throughput optimization, the simulation carries out and compares with the Monte Carlo method. In order to fully reflect the superiority and contrast integrity of the improved SFDE, the experiment will contrast those brand-new variants, such as JADE, MDE, EPSDE, and CODE. In addition, this experiment will also join the standard particle swarm algorithm to illustrate the Differential Evolution more effectively. The simulation in the hardware configuration for the INTER CORE i7, CPU 3.19GHz, memory 1.98GB running on the computer, the program written to run the software for the MATLAB 7.0. The parameters used in the system simulation experiment are shown in Table 5.
The convergence curve obtained by the simulation is shown in Figure 7. From the effective throughput convergence charts obtained in Figure 5, it can be seen that the optimal sensing time of the SFDE is 3.67ms and 3.7ms for Monte Carlo. The maximum effective throughput of SFDE is 5.2231bps/Hz, while Monte Carlo’s throughput is 5.1947bps/Hz. The SFDE achieves better goat than the Monte Carlo, which is close to the theoretical value of 5.3bps/Hz. For PSO, it is the maximum throughput with 4.5126bps/Hz, much worse than that of SFDE. JADE is second only to Monte Carlo, but its corresponding sensing time is 5.67ms, which is significantly worse. In summary, SFDE in the throughput performance is very good.

It can be seen from the experimental results in Figure 8 that the convergence rate of SFDE in the optimization process is much faster than that of others, especially Monte Carlo. The SFDE spends only 2.1987s on finding the best sensing time, while Monte Carlo spends 52.1423s. So the sensing time of SFDE is 23.7 times faster than Monte Carlo. It is even 5.69 times faster than the EPSDE. Monte Carlo is not dominant in the selection of sensing time; that is, it is not suited to real-time.

From the above experimental results, SFDE increases the system throughput with ensuring the convergence of the premise. Eventually SFDE greatly improve the convergence rate. In other words, SFDE overcomes the shortcomings of Monte Carlo and others that are not suitable for real-time, which makes it more practical in solving the effective throughput optimization of cognitive users.
6. Conclusion
In this paper, in order to optimize the radio time perception, so that the system throughput can be maximized, a Shuffled Frog-leaping Differential Evolution (SFDE) algorithm is put forward, which combines shuffled frog leaping with differential evolution. The experimental studies show that in the same situation the proposed SFDE algorithm improves the existing performance of others. This article simulates the effective throughput of cognitive users in radio and uses the improved SFDE to optimize. Experimental results show that, under the same situation, the SFDE can reach the theoretical value of throughput optimization accuracy, and the convergence speed is also faster.
Data Availability
The experimental Excel data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
Financial supports from the National Natural Science Foundation of China (no. 61572074), the China Scholarship Council for visiting to UK (Grant no. 201706465028), and the 2012 Ladder Plan Project of Beijing Key Lab of Knowledge Engineering for Materials Science (no. Z12110 1002812005) are highly appreciated.