Abstract
Typical approaches to designing quantum genetic algorithms are based on a concept of a qubit, a two-level quantum system. But many-valued quantum logic is more perspective from the point of view of the computational power. This paper proposes a quantum genetic algorithm based on a three-level quantum system in order to accelerate evolutionary process. Simulation using a set of standard test functions proves that the given algorithm is more effective and precise than the conventional quantum genetic algorithm.
1. Introduction
Evolutionary algorithms are one of the most widely used methods for solving the global optimization or search problem [1]. In its classical implementation, the smallest unit of information is a bit, a structure which may be in one of the following states: 0 and 1. It is important that a classical bit does not allow any intermediate states.
Quantum genetic algorithm (QGA) is a new evolutionary algorithm based on a combination of quantum computation and conventional genetic algorithms technology [2, 3]. This algorithm can be applied to the same set of problems the conventional genetic algorithm is used for, but it allows significantly accelerating the evolutionary process through quantum parallelization and entanglement of the quantum state. Probabilistic mechanism of the quantum computations combined with the evolutionary algorithm provides a global search for a solution with a rapid convergence and a small value of population. These algorithms have demonstrated their effectiveness for solving combinatorial and functional optimization problems, mechanical engineering optimization problems, image processing, and many others [4–8].
QGA main concepts were first proposed by Narayanan and Moore [9]. The basic unit of information used in a quantum computations is a qubit, a quantum system, which may be in the basis state or the basis state. Quantum nature of the qubit lies in the superposition principle, under which the qubit may be in any state which is a linear combination of basis states:with a normalization constraint:
As a consequence of the superposition principle, qubit state space is incomparably greater than state space of a classical bit. Information, stored in and amplitudes, is the quantum part of the information. Practically, та are numbers that specify the probability amplitudes of a qubit being in the states and , respectively.
The outcome of quantum measurement is a qubit in one of the basis states. However, it is important that the result of the measurement is not deterministic, like it is in the classical calculation, but probabilistic.
QGA uses a matrix representation for practical implementation:
Qubit (1) state in this representation can be defined as
A structured set containing qubits represents a quantum chromosome. State vector of such chromosome is a superposition state of basis states of the register , :
All the information about a qubits system is represented by state vector . The only thing that can be done to such a system is transforming the initial state vector into a new state by applying a quantum gate. So quantum genetic algorithm is a state transition process from some starting state to an ending state using a quantum gate algorithm. The information about the problem solution is only contained by the vector of the system ending state .
The main disadvantage appearing during quantum state rotation is the need to use a table to look up the rotating angle, which limits the universality of the search. Fixed rotating angle has a negative impact on the search speed, which is the reason for QGA sometimes being implemented as an adaptive process for determining the rotating angle. To enhance the local searching ability and to get out of the local optimal solutions, GDA can be extended with conventional operators traditionally used in the classical genetic algorithm, such as quantum mutation operation and quantum disaster operation [10].
2. Quantum Genetic Algorithm Based on Qutrits (QGA_3)
Qudit is a structure with several, more than two, states, which can also be used to encode a chromosome. A qudit is basically a unit of quantum information, which may be in any of states or in any superposition of those.
Ternary quantum logic is the simplest type of the many-valued logic. The basic unit of memory is called a qutrit. It has three basis states, , , and . The state of the qutrit can be represented as a superposition in the form of a linear combination:with a normalization constraint:
Matrix representation:
For the practical implementation, qutrit (6) state can be represented as
A system containing qutrits has basis states (as opposed to states for binary logic). When switching from binary to many-valued logic, we get an exponential increase in the number of basis states, which leads to an increase in the algorithm performance for the same search accuracy.
Quantum chromosome length is determined by the desired search accuracy , search area , and the number of quantum system basis states:where is a number of quantum system basis states.
With , search accuracy , and search area , the length of the quantum chromosome has to be at least 21 qubits, while with only 14 qutrits are required.
2.1. Qutrit Encoding
Matrix representation of a quantum chromosome is a structure which consists of qutrits.
| ⋯ | ⋯ | ⋯ | ||||||
| ⋯ | ⋯ | ⋯ | ||||||
| ⋯ | ⋯ | ⋯ |
Here determine the state of th qutrit, аnd qutrits determine one individual in the population. The initial state of the qutrits holds no information about the system state, so the easiest way to initialize the base population is by setting all probability state amplitudes to be equal to one another [2]. This means that at the end of the initialization each qutrit is in the state
2.2. Observation of Genes
Classical information about the problem’s solution is located in the ending state vector and can be retrieved as a result of quantum observation. Qutrit in one of the basis states, which is obtained as a result of quantum observation, is a classical representation of a quantum chromosome.
Based on the approach described in [11], the following pseudocode for quantum -qutrit chromosome state observation can be suggested (Algorithm 1). The product of its work provides the ground state (0, 1 or 2) with a probability of , and , respectively.
| ||||||||||||||||||||
The algorithm leads to the quantum chromosome being transformed into its classical representation in the ternary numbering system.
| ⋯ | ⋯ | ⋯ | ||||||
| ⋯ | ⋯ | ⋯ | ||||||
| ⋯ | ⋯ | ⋯ | ||||||
| Observation | ||||||||
| 0 | 2 | 0 | 1 | 1 | 2 | 1 | 0 | 0 |
2.3. Quantum Rotating Gates
All information about the problem and the algorithm for solving it is contained in the quantum gate, so its algorithm is the key issue of any quantum genetic algorithm construction.
The pseudocode for the genetic operator updating the qutrit state (quantum gate of the system evolution) as a result of probability amplitudes , , and modification is presented in Algorithm 2.
| ||||||||||||||||||||||
The algorithm’s concept can be illustrated by the image given below. The vertical lines set the probability amplitudes , , and . If the classical representation of the population best individual has, for example, a 2 on the position , all the amplitudes, except for , are multiplied by the factor . Its value lies in a range between 0 and 1 and is determined empirically. As a result, only the probability amplitude , which corresponds to the population best individual on the previous evolution stage, will increase. At the same time, probability amplitudes and will decrease.
It is important, that the described quantum gate algorithm eliminates the necessity to use a lookup table, which exists in the traditional QGA.
Initialized population () includes all the possible solutions with the same probability. This means that the QGA process starts with a random search. As the evolution progresses in time, the distribution undergoes changes (Figure 1).
The mechanisms for the local search with are illustrated in Figure 2. They demonstrate the role of parameter. Its fixed value has a bad impact on the convergence speed. So, it is desirable to use an adaptive quantum gate operator. The listed data shows that it can influence the individual in three possible ways:(1)The fitness of the individual is improved, but it does not become the best;(2)The fitness of the individual is improved, and it becomes the best;(3)The fitness of the individual is worsened.
This means that the bigger the number of individuals located in the area of the best individual , the more effective the process of the local convergence.
To sum up, QGA starts with the global search and switches to the local search automatically due to the change in the structural characteristics of the population individuals’ distribution.
2.4. Quantum Disaster Operation
Quantum population evolution may cause its falling into local optimal solution. This question is particularly important when the search area is large (big number of parameters) or when the fitness function has a complex topology around the global optimum.
Taking into consideration the small number of the quantum population, which is usually not bigger than individuals, getting out of the local optimal solution can be achieved by a quantum disaster operation [12]. It can be summarized as initialization of some population individuals, except for the best, with an initial state. Such operation “virtually” enlarges the size of the population and extends the area for optimal value search.
In the simplest case, quantum disaster operation can be implemented as Algorithm 3.
| ||||||||||||||||||
3. Simulation Test
QGA_3 based on a qutrits system is implemented using C++ programming language. The tests have been performed on the Intel Celeron CPU G1840 2.80 GHz, 4.0 GB RAM, Windows 7 Professional environment.
The majority of the quantum algorithms use similar primary ideas and approaches, which are based on the qubit representation of the population. For example, some algorithms based on the quantum computation ideas are presented in [13]. The studies of the quantum algorithms and the possibilities of their modification and improvement are carried out in [14, 15]. Other potential approaches also include the option of the transition to qutrits representation of the population and can be considered as directions for further research. Considering this, the only algorithms compared by effectiveness and performance were QGA and the suggested QGA_3.
The analysis and performance evaluation have been performed on a number of test functions [16]:
De Jong Function
The function has one global minimum solution:
Coldstein Price Function
The function has one global minimum solution:
Schaffer Function
The function has one global minimum solution:
Monopole and Six-Peak Camelback Function
The function has one global maximum solution:
Dual-Pole and Six-Peak Camelback Function
The function has two global minimum solutions:
Multipeak Positive Function
The function has one global maximum solution:
Ackley Function
The function has one global minimum solution:
Rastrigin Function
The function has one global minimum solution:
Optimization and Results. For QGA and QGA_3 simulation testing the following typical parameters were used: population size , total generations of iteration , and precision .
The Choice of Parameter. The value of parameter determines the evolution of the quantum system and has a key role in the algorithm. It determines the balance between the global (quantum observation operation) and local (quantum gate) search for the found solution specification.
To choose the optimal value we have investigated its influence on the behavior of the main parameters which characterize the quantum system evolution, such as the average fitness of the best individual and the standard deviation of function value over 1000 algorithm runs. Its characteristic behavior for the example of function optimization is presented in Figures 3 and 4. The quantum disaster operation was not applied in this case.
The research shows that if , the determinative operation in the population evolution is the quantum observation. Global search and the best individual’s approximation to the optimal value take place because of its probabilistic nature. At the same time, the quantum rotation angle is greater than and is too big to perform a local search.
As value increases, the rotation angle approaches its optimal value ~. It corresponds to ~: with this parameter value the fitness function and standard deviation reach the minimum value (see Figures 3 and 4).
Quantum system behavior for is critically dependent on its evolution time (see Figure 5). Indeed, with system simply does not have enough time to relax to the ending state, which corresponds to the fitness function optimal value. Enlarging the value of evolution time parameter practically simply causes the algorithm running time to increase, with no optimal value improvement.
The role of the quantum disaster operation can be illustrated by Figure 6. The results of the function optimization are given for parameters count . Taking into consideration the complex topology around the optimum, the size of the population .
Based on the given analysis it has been established that for QGA_3 without the quantum disaster operation is the optimal value. Adding the quantum disaster operation to the QGA_3 algorithm allows producing the best result with .
The optimization results of the – test functions with a standard set of input parameters are listed in Table 1.
In Table 1, denotes the function value of optimum; is the average value of the best individual fitness function; sd is standard deviation of fitness function; is average elapsed time need for the algorithm to run.
A search for and functions global minimum has been carried out separately, using parameter count and a standard set of other algorithm parameters (see Table 2). The results reveal that QGA_3 is more, but not enough, effective, especially in regard to function
The result can be significantly improved by enlarging the population size at least to and adding a quantum disaster operation, with the other parameters remaining unchanged. The results of experiments after the described improvements to the algorithm are listed in Table 3.
4. Conclusions
In this paper a new quantum genetic algorithm QGA_3 has been proposed. The algorithm is based on the idea of using qutrits for the representation of the population individuals. A simple algorithm, which does not require a lookup table, has been used for constructing the evolution operator which is responsible for the quantum system state updates.
The experiments have been performed using a typical set of test functions. On the basis of the experimental results we have concluded that QGA_3 allows decreasing the computational time while being more productive and having better convergence, compared to the conventional QGA.
We have also obtained that, in a case with a lot of parameters to be optimized, the QGA_3 efficacy can be improved by adding a quantum disaster operation. In the same time, the increase in algorithm operating time does not exceed 4%.
Further promising investigations include the research of -valued quantum logic units, qudits and their possible applications for representation of the population individuals. This is expected to be a subject of separate studies.
Conflicts of Interest
The author declares that there are no conflicts of interest regarding the publication of this paper.