Modelling and Simulation in EngineeringVolume 2011 (2011), Article ID 239743, 6 pagesdoi:10.1155/2011/239743

Research Article

Chaotic-Search-Based Cultural Algorithm for Solving Unconstrained Optimization Problem

Jianjia He1,2 and Fuyuan Xu1,2

¹Business School, University of Shanghai for Science and Technology, Shanghai 200093, China²Center for Supernetworks Research (China), Shanghai 200093, China

Received 16 March 2011; Revised 25 May 2011; Accepted 29 June 2011

Academic Editor: Farouk Yalaoui

Copyright © 2011 Jianjia He and Fuyuan Xu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

For premature convergence and instability of cultural algorithm in solving function optimization problem, based on cultural algorithm and chaos search optimization, a chaos cultural algorithm (CCA) is proposed. The algorithm model consists of a chaos-based population space and a knowledge-storing belief space, uses normative knowledge and situational knowledge for chaos search and chaos perturbation, respectively, effectively avoids premature convergence of cultural algorithm, and overcomes chaos search optimization's sensitivity to initial values and poor efficiency. Test results show that this algorithm is strong in global search and has good performance in searching efficiency, precision, and stability, especially in solving high-dimensional optimization problem.

1. Introduction

Since the found of the bionics in the 1950s, people began to simulate the evolution mechanism of creatures and proposed such evolutionary computation approaches as the genetic programming and immune algorithm to solve the optimization problem. These algorithm researches are mainly focused on the level of natural choice. Many practical cases show that culture can make population evolve and adapt to the environment at certain speed beyond the evolutionary speed which is merely based on the genetic inheritance [1]. In the human society, culture is regarded as the store carrier of information which is widely delivered among and inside the social groups and can be inherited by all the social members so as to effectively guide them to solve problems. Enlightened from this, Reynolds put forward a brand-new evolutionary computation of cultural algorithm [1].

The cultural algorithm simulates the evolution from both the macroscopic and microscopic levels which stimulate and influence each other mutually. This algorithm provides a kind of explicit mechanism for acquiring, preserving, and integrating the solution knowledge so as to improve the search efficiency. Nowadays, the cultural algorithm, as the hotspot of the algorithm research, has been successfully applied to many fields and achieved better results at certain problems than the traditional evolution algorithm [2–4]. As for the function optimization, the cultural algorithm possesses high solution efficiency, but it has such deficiencies as premature convergence, instability, and so on. The chaotic search optimization is characterized with high global search ability, but it is sensitive to the initial values and poor in searching efficiency. Therefore, this study puts forward the chaos cultural algorithm (CCA) by introducing the chaotic search optimization into the cultural algorithm model. Because of taking advantage of the stochastic and ergodic properties of the chaotic search optimization, this algorithm increases the global search ability and avoids the premature convergence; by making use of the knowledge guidance of the cultural algorithm, it improves the search efficiency and precision.

2. The Principle of the Cultural Algorithm

The cultural algorithm framework consists of the population space and the belief space which simulate the evolutionary process from the microscopic and macroscopic levels separately as shown in Figure 1. The former is the solution space of questions, while the latter is used for forming, storing, and delivering knowledge experiences. They are mutually independent and interrelated and exchange information with each other through specific protocol. First of all, the function of accept () delivers the excellent individual experience in the evolution of the population space to the belief space, then the function of update () renews the population experience or knowledge based on the existing and new experiences in the belief space; subsequently, the function of influence () guides the evolution of the population space by means of population experience or knowledge. Furthermore, the objective function of objective () is utilized to evaluate the individual fitness values in the population space. The function of generate () forms the next-generation individuals in light of the population experience. And based on certain rules, the function of select () selects part of the newly generated individuals as their own parent. The pseudocode of the cultural algorithm is illustrated in Algorithm 1 [5].

Begin

;

Initialize Population POP(t);

Initialize Belief Space BLF(t);

Repeat

    Evaluate Population POP(t);

    Update (BLF(t), Accept (POP(t)));

    Variation (POP(t), Influence (BLF(t)));

    t = t + 1;

    Select(POP(t)from POP());

Until termination condition achieved

End

Algorithm 1

Figure 1: The framework of cultural algorithm.

3. The Design of the Chaos Cultural Algorithm

In this study, the chaotic search optimization [6] is embedded into the cultural algorithm model, and then, the chaos cultural algorithm comes into being by using the interval schemata proposed by Chung [7] and Chung and Reynolds [8] to express the knowledge in the belief space. In this way, the global search ability of the chaotic search and the experience knowledge from the iterative process are used to direct the efficient search of the algorithm to improve the search efficiency. In addition, this algorithm adopts the parallel structure of the chaotic search to efficiently reduce its sensitivity to the initial values. Without loss of generality, the issue of the nonlinear functions optimization can be described as follows:

3.1. Chaotic Search Optimization

As a universal nonlinear phenomenon, chaos is of stochastic, ergodic, and regular properties, where ergodicity means that the chaotic serial traverses all the states nonrepeatedly in a certain range. In optimization, ergodicity can be used as a kind of optimization mechanism to effectively avoid the local optimum, just on which the chaotic search optimization is proposed. Bing and Weisun [6] raised a chaotic search algorithm by use of the logistic chaos mapping (see (2)) and the second carrier wave, which attracts widespread attention and causes numerous researches. This algorithm is strong in global search, but it is sensitive to the initial values and poor in the search efficiency. Some papers are aimed to improve the chaotic search optimization method, just like P-Chaos presented by Liang and Gu [9]. It starts parallel computing from groups of initial points and conducts double frequency carrier after searching to some extent to find the optimal solution as soon as possible. But P-chaos can hardly ensure a point is in the neighborhood of real optimal solution after first carrier, which is the foundation of the second carrier. Also, the initial searching range is relatively large, while the precise search is required later on. Last but not least, the searching is always along one side of previous one, equally one-side search failing in global search. As a result, there is still more work in improving the chaotic search optimization

In the above formula, when , the logistic mapping is chaotic and the interval () is the chaos invariant set [10].

The underlying idea of the chaos search optimization is to transform the chaotic state into the optimization variables by means of similar carrier, enlarge the traversing scope of the chaotic motion to the value range of the optimization variable, and then search by using the chaotic variable [6]. The basic steps are as follows.

Step 1. Initialization. Set , , the maximum times of chaotic motion as M. Assign different initial values to the chaotic variable in (2) (the value cannot be the fixed points of the logistic mapping: 0.25, 0.5, 0.75). The chaotic variable is mapped to the value range of the optimization variables according to (3) and the initial optimization variable is obtained. Set ,.

Step 2. Calculate the chaotic variable sequence and the optimization variable sequence.

Step 3. If , then judge: if the inequality holds, set, ; vice versa, keep and invariant and set , and then turn to Step 2. If, end the chaotic search.

After a certain number of iteration, the logistic mapping () distributes a large number of track points (chaotic sequence) in the interval () and its probability density functionis [10]: This formula is the Chebyshev distribution as shown in Figure 2.

Figure 2: Probability density function of chaotic sequence of logistic map.

It can be seen from Figure 2 that the distribution features of logistic mapping chaotic sequence is that it is even at the middle, plenty at both sides, and up to the infinite many at the two ends. By applying this distribution feature and combining with the cultural algorithm superior interval (discussed later), each decision variable of the individuals in the population space can distribute near to the two ends of the superior interval so as to reduce the blindness search and improve the search efficiency.

3.2. The Structure of Belief Space

The structural definition of belief space is very important in the cultural algorithm framework. The primary function is to provide a definite mechanism for acquiring, storing, and delivering knowledge experience. The two kinds of most fundamental knowledge of the cultural algorithm is applied into this paper: situational knowledge and normative knowledge, that is, .

In the algorithm of this study, S is a schedule which stores the set of excellent individuals chosen from the population space at present which are the other individuals’ model, that is, . The data structure is shown in the Figure 3 [11]. Here, every excellent individual involves its own value of decision variable and the responding adaptive value, and In this paper, reads values directly from the scale popsize of the population space. When initializing, the best individual set reads from all the individuals of the initial Population Space.

Figure 3: Data structure of situational knowledge.

There stores sets of variation intervals of decision variables, namely, superior intervals in . These intervals present the most possible range of the optimal solution, which leads the chaotic iteration to arrive at the optimal solution as soon as possible. The data structure is described in Figure 4.

Figure 4: Data structure of normative knowledge.

Therein,and represent the upper and lower boundaries of the th decision variable separately.andstand for the responding fitness values of the upper and lower boundaries, respectively. In the initial situation,andare set as the definition domain intervals, andandas. During the evolutionary process, these four values will adjust themselves continuously in light of the intervals of the excellent individuals.

3.3. Renovation of Belief Space

The situational knowledge in the belief space renovates twice in each generation: one happens after the chaotic search, and the other happens following the chaotic disturbance. The normative knowledge renovates once in each generation.

3.3.1. Renovation of the Situational Knowledge

The approach of renovating the situational knowledge in this paper learns from the tournament method. The specific approach is as follows.(i)After each iteration, combine all the individual sets and the optimal individual sets represented by the situational knowledge into a new set (the scale of the new set in this algorithm is popsize).(ii)Randomly select competitors in the new set (generally). Compare the fitness values of each individual in the new set and the competitors. Record the number of success for each individual. And set the number of success of the optimal fitness individual as the highest so as to make sure the optimal individual stored in the situational knowledge. (iii)According to the number, select the top individuals (here, ) as the new situational knowledge.

3.3.2. Renovation of Normative Knowledge

Suppose the th individual effect the lower boundary of the th decision variable and the th individual effect the upper boundary of the th decision variable, and then, the upper and lower boundaries of and the responding fitness values are separately showed as follows:

3.4. Accept Function

The population space delivers the individual experience to the belief space by the accept function. In this way, it actually provides a group of optimal subsets to the belief space. In the case of optimization, it is usual to select the top several individuals in a linear order by a certain percentage (generally 20%–25%). Also, set the dynamic accept function [12] and fuzzy accept function [13, 14] according to different questions. Because of the application of tournament method to the renovation of the situational knowledge in this study, the processing of the accept function has also been changed to some extent. In the renovation of the situational knowledge, the accept function accepts all the individuals at present; in the renovation of the normative knowledge, the accept function acquires the TOP several of the situational knowledge (here, ).

3.5. Influence Function

In the iterative process, the knowledge in the belief space influences the optimization variables in two ways: (1) The chaotic search at the superior interval represented by the normative knowledge; (2) The chaotic disturbance at the excellent individuals represented by the situational knowledge.

3.5.1. Chaotic Search by Using the Normative Knowledge

The normative knowledge stores the most possible range of the optimal solution which is described with the superior interval. The chaotic variable caused by the logistic map is owned with the characteristics of “many at the two sides, few at the middle” at the interval (). In this paper, based on the responding superior interval, every decision variable of the individuals is the value range of is classified into three intervals: , , . And then, the chaotic variables are randomly mapped to one of the three intervals. By doing this, any decision variable of the new individuals will be mainly distributed near the superior interval so as to direct the chaotic search by using the normative knowledge

3.5.2. Chaotic Disturbance by Using the Situational Knowledge

The situational knowledge is a list storing the excellent individuals at present that is the model of the other individuals so as to be able to direct the generation of the other individuals. The length of the list is just the scale of the population space. In the algorithm of this study, the situational knowledge is the excellent individual experience formed from the chaotic iteration and stores the specific position of the excellent individuals. Therefore, it makes up for the loss of global superiority because of the restriction of the normative knowledge to the chaotic variables, so the situational knowledge possesses strong global guidance ability and plays the escaping role in the searching process.

This study firstly carries on the chaotic search, secondly renovates the situational knowledge and obtains some excellent individuals, then conducts the chaotic disturbance near the excellent individuals which is helpful for picking out better individuals and acquires the global distribution information of the optimal solution at present. As a result, it is effective to avoid being trapped in the local optimal solution. The specific description is as follows:

3.6. Process of the Chaos Cultural Algorithm

In summary, this section describes the process of the proposed chaos cultural algorithm in this study to solve the optimization problem of the nonlinear function as follows.

Step 1. .

Step 2. Initialize the chaotic variable and generate randomly the initial chaotic variable which cannot be the fixed point of the logistic mapping: 0.25, 0.5 and 0.75.

Step 3. Initialize the population space and the relevant parameters.

Step 4. Initialize the situational knowledge and normative knowledge in Belief Space.

Step 5. Evaluate the individual fitness in the population space: if meeting the requirement, terminate the process; if not, continue the process.

Step 6. Renovate the situational knowledge and normative knowledge in the belief space.

Step 7. Iterative search of the chaotic variables. Generate the new chaotic variables according to (2).

Step 8. Under the lead of the normative knowledge in the belief space, conduct the chaotic search according to (7).

Step 9. Evaluate the fitness of individuals in the population space and renovates the situational knowledge in the belief space.

Step 10. Under the guidance of the situational knowledge in the Belief Space, conduct the chaotic disturbance according to (8).

Step 11. .

Step 12. Turn to Step 5 until meeting the terminating conditions.

4. Case Study

The chaos cultural algorithm combines the advantages of both the chaotic search optimization and the cultural algorithm. For convenient comparison, this paper selects 6 functions from [15, 16] to make the optimization design and make comparisons with the parallel chaotic optimization algorithm (Algorithm  1) based on the information exchange of competition-cooperation in [16] and with the cultural algorithm (Algorithm  2) based on the evolution planning in [15]. Both Algorithms  1 and 2 are forefront in their respective fields and able to solve the optimization problem of functions effectively. The algorithm in this study takes the scale of the population space as popsize = 20. Each algorithm runs 30 times independently. The optimization function is as follows [15–18]:

Here, ₁ is the large-variable interval function; ₂ and ₃ are the functions with local extreme points; ₄ and ₅ are high-dimension complex functions with the dimension up to 30. ₆ is a complex function whose optimization point is difficult to obtain by the general methods. The comparative results between the chaos cultural algorithm and the other algorithms are shown in the Table 1 ( is the real optimal value; stands for the average optimal value of the 30 operations; is the average time of the 30 operations with millisecond as the unit; is the superiority ratio, that is, the frequency percentage of obtaining the optimal values).

Table 1: Comparison of CAA and others.

Function

Algorithm  1
Algorithm  2
This algorithm

₁
−1.002
−0.154877
2865
40
−0.526565
1626
53
−1.002
451
100

₂
−1
−0.998745
965
80
−0.996437
445
76
−1
80
100

₃
0
0.148157
1824
73
0.019447
528
90
0
596
93

₄
0
0.284856
6478
0
0
5733
100
0
6547
100

₅
0
68962.5
5156
0
0.014398
10675
56
0
6189
100

₆
0.002288
0.384578
25157
0
0.006874
25345
0
0.002846
25247
10

It is shown in Table 1 that the algorithm adopted in this paper is much better than the other two algorithms. Algorithm  1 takes long time and has difficulty in dealing with the high-dimension and large-variable intervals. Algorithm  2 is easy to be trapped in the local extreme points so that the stability is not high. The chaos cultural algorithm not only overcomes such deficiencies, but also is superior to the other algorithms in search efficiency and precision.

5. Conclusions

In this paper, the authors introduce the chaotic behavior in chaotic search optimization to cultural algorithm and presents a new type of Intelligent optimization algorithms — chaos cultural algorithm, the marriage of chaotic search and cultural algorithm. The chaos cultural algorithm makes use of the ergodic and stochastic properties of the chaotic search optimization so as to increase the global search ability and makes use of knowledge guidance of the cultural algorithm to decrease the frequency of the blindness search to improve the search efficiency and precision and to increase the generality. Based on the parallel structure, this algorithm adopts the chaotic search led by the normative knowledge and the chaotic disturbance led by the situational knowledge, both of which perform interactively. Therefore, this algorithm is effective to reduce the sensitivity to the initial values and enhance the stability and improve the overall performance.

In solving chaos-cultural-algorithm-based unconst-rained optimization problem, the authors come up with phenomenological approach, solution procedures and algorithm flow, which is tested and verified from and experimental viewpoint. The result shows that the chaos cultural algorithm takes the advantages of both chaotic search and cultural algorithm in solving optimization problem to achieve mutual complementation. The algorithm boasts not only higher solving efficiency and precision overall, but also successfully avoiding search converge to local optimal solution prematurely.

Acknowledgments

This work is supported by National Natural Science Foundation of China (Grant no. 71171135), by the Key Discipline of Shanghai (III) (Grant no. S30504), and by the Innovation Fund Project For Graduate Student of Shanghai (Grant no. JWCXSL1001).

References

R. G. Reynolds, “An introduction to cultural algorithms,” in Proceedings of the 3rd Annual Conference on Evolutionary Programming, pp. 131–139, World Scientific Publishing, San Diego, Calif, USA, 1994.

C. A. Coello and R. I. Becerra, “Evolutionary mutilobjective optimization using a cultural algorithm,” in Proceedings of the IEEE Swarm Intelligence Symposium, pp. 6–13, IEEE Press, Indianapolis, Ind, USA, 2003.

S. Saleem and R. Reynolds, “Cultural Algorithms in dynamic environments,” in Proceedings of the 2000 Congress on Evolutionary Computation, pp. 1513–1520, July 2000.

X. H. Yuan and Y. B. Yuan, “Application of cultural algorithm to generation scheduling of hydrothermal systems,” Energy Conversion and Management, vol. 47, no. 15-16, pp. 2192–2201, 2006.

R. G. Reynolds and S. Zhu, “Knowledge-based function optimization using fuzzy cultural algorithms with evolutionary programming,” IEEE Transactions on Systems, Man, and Cybernetics, Part B: Cybernetics, vol. 31, no. 1, pp. 1–18, 2001.

L. Bing and J. Weisun, “Chaos optimization method and its application,” Control Theory and Applications, vol. 14, no. 4, 1997 (Chinese).

C. J. Chung, Knowledge-Based Approaches to Self-adaptation in Cultural Algorithms, Wayne State University, Detroit, Mech, USA, 1997.

C. J. Chung and R. G. Reynolds, “A testbed for solving optimization problems using cultural algorithms,” in Proceeding of the 4th Annual Conference on Evolutionary Programming, pp. 225–236, MIT Press, Cambridge,Mass, USA, 1996.

H. Y. Liang and X. S. Gu, “A novel chaos optimization algorithm based on parallel computing,” Journal of East China University of Science and Technology, vol. 30, no. 4, pp. 450–453, 2004.

D. X. Yang, G. Li, and G. D. Cheng, “Comparative study on chaos optimization algorithm for nonlinear function,” Chinese Journal of Computational Mechanics, vol. 21, no. 3, pp. 257–262, 2004.

C. J. Chung and R. G. Reynolds, “An evolution-based tool for real-valued function optimization using cultural algorithms,” Journal on Artificial Intelligence Tools, vol. 7, no. 3, pp. 239–292, 1998.

S. M. Saleem, Knowledge-Based Solution to Dynamic Optimization Problems Using Cultural Algorithms, Wayne State University, Detroit, Mech, USA, 2001.

S. N. Zhu and R. G. Reynolds, Fuzzy Cultural Algorithms with Evolutionary Programming VII, Springer, Berlin, Germany, 1998.

S. N. ZHU, Fuzzy Cultural Algorithms with Evolutionary Programming for Real-Valued Function Optimization, Wayne State University, Detroit, Mech, USA, 1998.

X. F. Yuan, Y. N. Wang, and L. H. Wu, “Parallel chaotic optimization algorithm based on competitive-cooperative inter-communication,” Control and Decision, vol. 22, no. 9, pp. 1027–1031, 2007 (Chinese).

C. J. Chung and R. G. Reynolds, “Knowledge-based self-adaptation in evolutionary search,” International Journal of Pattern Recognition and Artificial Intelligence, vol. 14, no. 1, pp. 19–33, 2000.

C. A. Coello and R. L. Becerra, “Efficient evolutionary optimization through the use of a cultural algorithm,” Engineering Optimization, vol. 36, no. 2, pp. 219–236, 2004.

A. D. Belegundu, Study of Mathematical Programming Methods for Structural Optimization, University of Iowa, Iowa, USA, 1982.