Abstract
A Discrete Artificial Bee Colony (DABC) is presented for joint symbol detection at the receiver in a multidevice SpaceTime Block Code (STBC) MutliInput MultiOutput (MIMO) communication system. Exhaustive search (maximum likelihood detection) for finding an optimal detection has a computational complexity that increases exponentially with the number of mobile devices, transmit antennas per mobile device, and the number of bits per symbol. ABC is a new populationbased, swarmbased Evolutionary Algorithms (EA) presented for multivariable numerical functions and has shown good performance compared to other mainstream EAs for problems in continuous domain. This algorithm simulates the intelligent foraging behavior of honeybee swarms. An enhanced discrete version of the ABC algorithm is presented and applied to the joint symbol detection problem to find a nearly optimal solution in real time. The results of multiple independent simulation runs indicate the effectiveness of DABC with other wellknown algorithms previously proposed for joint symbol detection such as the nearoptimal sphere decoding, minimum mean square error, zero forcing, and semidefinite relaxation, along with other EAs such as genetic algorithm, estimation of distributions algorithm, and the more novel biogeographybased optimization algorithm.
1. Introduction
MultiInput MultiOutput (MIMO) communication systems can offer spatial diversity gains in the fading channels and have significantly higher channel capacity than the SingleInput SingleOutput (SISO) systems for the same total transmission power and bandwidth [1, 2]. The system proposed in this paper comprises of one receiving station and multiple transmitting devices. The receiver’s front end has multiple antennas, and each transmitting device has multiple transmit antennas. Employing the Space Time BlockCode (STBC) is realized to increase the capacity of MIMO systems and consequently improves data throughput and spectral efficiency [3]. Multiantenna systems are widely used because of their ability of dramatically increasing the channel capacity in fading channels [4]. Each transmit device uses an STBC; the receiver side performs the joint signal detection. Such a system is referred to as a multidevice (MD) STBCMIMO system. Generally in an MDSTBCMIMO system, the number of receive antennas is typically smaller than the cumulative number of transmit antennas used by all transmitting devices in the system. An example of MDSTBCMIMO, with a smaller number of antennas at the base station or access point, would be the uplink multiple access communication in cellular systems.
This paper addresses the symbol detection in MDSTBCMIMO systems. The Maximum A Posteriori (MAP) detection, which reduced to the Maximum Likelihood (ML) detection in the case of a priori equally likely symbol blocks, minimizes the probability of detection error, and thus is optimal and is further explained in Section 3. However, a computationally efficient algorithm for achieving MAP or ML detection is not known. Some studies with Sphere Decoding (SD) algorithms exhibit that their expected computational complexity grows polynomially with the problem size up to some value of for the cases of small constellation sizes [5], but it grows exponentially for the cases of large constellation sizes. Also, for some SD algorithms, operation at a low SNR requires inordinately high computation; yet operation at a high SNR is efficient. In fact, [6] shows that even the expected computational complexity of the SD grows exponentially with the problem size in MIMO communication systems. In any case, an algorithm with polynomial growth of expected complexity for all values of the problem size has not yet been found.
Due to the combinatorial nature of the problem, the ML detection is a choice to obtain the optimal solution; yet it has a high computational complexity. Taking advantage of heuristic algorithms, evolutionary algorithms (EAs) more specifically, and their ability to solve optimization problems efficiently facilitates finding an optimal solution with relatively low computational cost. evolutionary algorithms (EAs) are a subset of heuristic algorithms, which is inspired by biological evolution and mutation.
The Artificial Bee Colony (ABC) algorithm is one of the novel EAs that has been introduced by Karaboga in [7]. He presented the ABC algorithm for real (continuous) parameter optimization in unconstrained optimization problems in [8]. ABC is a populationbased, stochastic global optimizer Evolutionary Algorithm. It is based on the theory of foraging bees searching for food sources for their nectar and sharing the information of food sources’ locations to other bees in the hive. This algorithm demonstrates good accuracy and efficiency, compared with other mainstream EAs. In [8], Karaboga and Akay showed that ABC algorithm outperforms other EAs such as Differential Evolution (DE) [9], Particle Swarm Optimization (PSO) [10], and Genetic Algorithm (GA) [11] for numeric problems with multidimensions. In his recent paper [12], Karaboga compared the ABC programming with various genetic programming techniques and crossover methods available in [13] and demonstrated the superiority of ABC over these schemes through simulations. After introducing the ABC algorithm in 2007, some papers are published on the applications of ABC to different optimization problems [14–16], in addition to various realworld applications including filter design [17], image processing [18], control engineering [19], computer science [20, 21], neural networks [22], and even biology [23]. In these studies, ABC outperforms other EAs and is turned into a popular global optimization solver to the continuous optimization problems and applications.
The ABC algorithm discussed in [24] is primarily presented for continuous functions; yet the MDSTBCMIMO detection problem is in the discrete domain. The impressive results of ABC implementation for continuous problems bring up the idea that a discretized versions of ABC can be a potential highefficient lowcomplex solver for discrete or numerical optimization problems. In this paper we introduce a discrete ABC algorithm, in addition to some new features to enhance the overall algorithm’s performance. Although there are a few discretized version of ABC in the literature, we will discuss in Section 4 that their efficiency is the result of their high complexity; thus they cannot be utilized for a realtime symbol detection problem such as MDSTBCMIMO.
Our simulation results show that discrete ABC can meet the best known semioptimal detector (i.e., SD) with less complexity and has better performance than other methods such as Minimum Mean Square Error (MMSE), Zero Forcing (ZF), and SemiDefinite Relaxation (SDR) [25], while it outperforms other EAs such as GA, Estimation of Distributions Algorithm (EDA) [26], and the recently proposed BiogeographyBased Optimization (BBO) [27].
In the rest of this paper, the system model is presented in Section 2. The application of existing symbol detection algorithms is discussed in Section 3. The discrete ABC algorithm and its application to the symbol detection problem are presented in Section 4. The computational complexity comparison of EAs and other solvers is discussed in Section 5. The simulation results are presented in Section 6, and Section 7 contains the conclusion and the future work.
2. System Model
The system consists of mobile devices transmitting signals and one receiver. This system can model the uplink communication of the cellular system. Each mobile device has transmit antennas that apply STBC, whereas the receiver front end has receive antennas. The multiple mobile devices in the proposed systems can cause cochannel interference. An IQmodulation scheme (e.g., QAM, PSK) maps source information into complex numbers. Even if each transmit device employs an orthogonal spacetime code, the absence of coding across different wireless devices cannot guarantee the orthogonality among their signals. In the case of a single mobile device , the wireless device transmits using transmit antennas and communicates with a receiver that has antennas. The number of time slots in the spacetime block code is denoted by . The channel is assumed to be quasistatic; that is, the channel gain remains constant during each time block of data. It is also assumed that the channel gain at each time block is known to the receiver. This assumption is often used in the literature and reasonable if training or pilot signals are used. A complex dimensional matrix represents the MIMO channel, and another complex dimensional matrix represents the input signal in a spacetime block code. The relationship between the input and output signals can be expressed as where is a dimensional complex output matrix and denotes the additive white noise matrix.
Equation (1) describes the relation between the input transmitted signals and the output received signals in terms of complexvalued matrix equation. The relation between input and output of the channel in a system with linear dispersion spacetime coding can be equivalently expressed in terms of a realvalued matrix equation. We now briefly discuss the realvalued matrix equation. The input signal in (1) in the case of the linear dispersion code [28] is denoted by a complexvalued matrix that takes the following form: Here, is the complex number that represents the th symbol, where and correspond to the real and imaginary parts of the symbol, respectively, and indicates the number of symbols conveyed in a spacetime code block. In the IQ constellation diagram, and are discretevalued variables, such that corresponds to a symbol in the constellation diagram. In 4QAM, for example, each of these two variables can take values of , and thus determines one of the four possible symbols arranged in the square grid of , , , and . These symbols can be represented as a dimensional realvalued row vector , whose components are constituted by , , . The real and imaginary parts of matrix ’s components can be arranged as a dimensional realvalued row vector . The relation between and in this new alternative form can be expressed as where is a dimensional realvalued vector representing noise, and is a realvalued matrix derived from the component of matrices , , , and . Equation (1) in the case of multiple wireless devices can be expressed as where is a dimensional complex matrix of the input signal from wireless device and the dimensional complex matrix represents the channel from the th device to the receiver. As a result, (3) can be written as where , is a dimensional realvalued row vector that represents the complex symbols sent from the th wireless device in a spacetime code block. Note that (5) can model the case in which different wireless devices use different code rates and different spacetime codes. The total number of symbols transmitted from all wireless devices in a spacetime code block through all transmit antennas is denoted by .
3. Signal Detection
The ML detection is known to yield the lowest symbol error probability in the case of a priori equally likely symbols. In the case of our problem, the detector at the receiver has to choose from possible sequences of symbols transmitted in a spacetime code block, where is the size of the symbol constellation associated with the modulation scheme. ML detection chooses transmitted symbols that maximize . In the case of additive white Gaussian noise , the ML detection is reduced to choosing the vector from possibilities that has the shortest Euclidean distance , which is expressed as
The ML detection scheme can be implemented by searching through all possible symbol sequences, where and is the size of the symbol constellation. Performing such an exhaustive search to find the minimum Euclidean distance in (6) is computationally inefficient, especially for large . Computational complexity increases exponentially with , the number of bits per symbol, transmit antennas per device, and the number of wireless devices . Highspeed communication requirements demand a lowcomplexity detection scheme. For lowcomplexity nearoptimal detection, in this paper the ABC algorithm is applied to this MDSTBCMIMO detection problem. The MDSTBCMIMO detection problem is converted into a discrete optimization problem that searches the space of symbol combinations. Section 4.2 describes how discrete ABC is applied to the signal detection. In the following sections, the performance of the discrete ABCbased detector is compared with other lowcomplexity suboptimal algorithms such as MMSE, ZF, SDR, and SD.
4. Discrete Artificial Bee Colony
This section presents the discrete Artificial Bee Colony (ABC) algorithm. A general description of ABC is given in the next subsection, followed by the modified discrete ABC algorithm in the subsequent section.
4.1. The Artificial Bee Colony Algorithm
evolutionary algorithms (EAs) have been often used to solve difficult optimization problems. Most of the EAs are inspired by the theory of biological evolution (e.g., selection, crossover, mutation, recombination, and reproduction). The ABC algorithm has been recently presented by Karaboga for real (continuous) parameter optimization in unconstraint optimization problems [24], which is based on a particular intelligent behavior of the honey bee swarms. This algorithm demonstrates good accuracy and efficiency compared to other EAs such as differential evolution (DE) [29], ant colony optimization (ACO) [30], PSO, and GA, for numeric problems with multidimensions [8].
Consider an optimization problem where is a vector and is a constraint set. In the original ABC, each candidate solution is represented by a vector variable of the optimization problem. In the context of evolutionary algorithms, a candidate solution is also referred to as an “individual,” and a group of candidate solutions is referred to as a “population” of individuals. In ABC, each individual (candidate solution to an optimization problem) is analogically considered as a food source position. The fitness value, , of each individual (food source) corresponds to the nectar quality of the food source.
This algorithm imitates the behaviors of the real bees on finding food source locations and sharing the information of food sources to the other bees in the nest. In this algorithm colony bees are classified into three types with certain responsibilities: employed bees, onlooker bees, and scout bees. Employed bees are the bees that have already been assigned to a food source. Each of them saves the food source position and selects another food source in her neighbor and chooses out of two the one that has a better nectar. Then they return to the hive and start to dance based on the quality of the nectar of their associated food source. An onlooker bee watches the dance of employed bees at the hive and selects an employed bee based on the dances observed so that the probability of choosing an employee bee is proportional to the nectar quality of that employee bee. Then the onlooker bee receives the information of the chosen employed bee associated with food source (the food source position and its nectar quality) from her and becomes an employed bee associated with that food source. Since then, the new employed bee (former onlooker) performs the same act as the employed bee in the previous phase; that is, she searches for a new food source in the neighbor of her associated food source for higher nectar quality and saves the best food source and its nectar to her memory. Finally, scout bees are free bees responsible for finding new food sources and evaluating their nectar. As soon as a scout bee finds a food source, she turns into an employed bee.
The algorithm assumes that there is only one employed bee for every food source; thus the number of employed bees is equal to the number of individuals in the population . If there is no improvement in the nectar quality of a food source after certain number of trials, the food source will be abandoned and the employed bee assigned to that food source will become a scout bee that looks for a new food source. A pseudocode of the ABC algorithm is given in Pseudocode 1.

At the first step, ABC generates randomly distributed initial food source positions of the size , whereas each individual solution , is a dimensional vector of numbers. In this step each scout bee that finds a food source location saves the current location in her memory and becomes an employed bee. In the employed bees phase, each employed bee finds a new food source position in the neighborhood of her current associated food source , and if the new food source has a better nectar, she saves the new position to her memory. In the original ABC algorithm, an employed bee locates the new food source positions using the following expression [24]: In this equation is the th component of the th individual of the population, is a randomly selected food source location (different from ), and is a randomly chosen index. is a random real number between that controls the production of a neighbor food source around .
At the beginning of the onlooker bees phase, employed bees share their information about the quality of food sources with onlooker bees. An onlooker bee chooses an employed bee to take the food source information based on the following probability: where is the fitness value of the th solution in the population . After an onlooker bee selects a food source, she becomes an employed bee and locates a new food source in the neighborhood using (8), then she compares its nectar quality with the current food source, and saves the food source position that has a better nectar quality to her memory and returns back to the hive to share this information.
If the number of trials of a food source is not improved through a predetermined number of trials , it will be removed from the population, and the employed bee assigned to that bee becomes a scout that searches for a new food source. Each component of the new food source is randomly selected from , where and are the minimum and maximum of the allowable values in the problem domain. The previous steps are repeated until the termination condition is satisfied, which here is a preset maximum number of generations .
4.2. Discrete ABC Algorithm for Joint Symbol Detection
This subsection introduces a discrete version of the ABC algorithm. The ABC algorithm discussed in the previous subsection is for optimization problems in the continuous domain. However, some of the previous steps have to be modified because the decision variables, which in this case are the transmitted symbols of the MDSTBCMIMO problem, are a set of nonnegative integer numbers. In the discrete ABC (DABC) algorithm, we define a new expression to search in the neighborhood of the current food source position as a replacement to expression (8): where “randint ” returns a random integer number between and (more specifically a random integer number starting from to . Note that if the result of (10) falls beyond the problem’s integer domain, that number is replaced with the closest integer defined within the boundaries. This expression is used in both employed bees and onlooker bees phases. Moreover, the onlooker bees select employed bees with the following probability: which has been observed to increase the efficiency of the algorithm more than other selection methods and the one in (9). The evaluation of (10) is explained in the appendix.
We have applied more enhancements to the algorithm in order to reduce its computational complexity. The most complex section of ABC, and most other EAs, is the fitness function evaluation. Therefore, if the number of these evaluation decreases, the algorithm runs faster. In this version of the discrete ABC, during the employed bees and onlookers bees phases, the algorithm is set to only evaluate those individuals that are modified during the greedy selection process. When the algorithm uses (10) to select a neighbor food source, it may not always return a new food source position due to the stochastic nature of (10). The DABC is set to check whether a food source has been modified prior to proceeding to the fitness function evaluation. In this case, it eliminates a number of fitness evaluations for some individuals that were already evaluated during the previous generations.
The last phase of the DABC is similar to the scout bee phase of the original ABC. In the scout bee phase of DABC, the algorithm selects only one food source that exceeds maximum allowable number of trials to abandon and sends one scout bee to explore new food source positions. The scout bee randomly selects one food source, evaluates its nectar quality, and saves it in her memory. This procedure helps the algorithm to explore the search space more effectively, which is an advanced version of the mutation process in GA and some other EAs that they randomly mutate any individual of the population. A detailed pseudocode for the discrete ABC algorithm is given in Pseudocodes 2, 3, 4, and 5.



We found a few articles in the literature that employ ABC for numerical optimization problem. In [31], Tsai et al. presented an enhanced ABC, which is applied to numerical optimization problems, and called it the interactive ABC (IABC). They have modified the way in which onlooker bees choose a neighboring food source position. IABC incorporates the concept of universal gravitation into the consideration of the affection between employed bees and onlooker bees, and their simulation results demonstrate the high performance of IABC compared with ABC and PSO. Note that the higher performance of their new method comes with the price of a significant increase in the complexity of the algorithm—there is heavy computational load for calculating the gravitation between employed bees in (10) and (11) in [31] for every food source in every algorithm iteration. As mentioned before, an important aspect of EAs’ implementation for wireless communication problems, specifically the MDSTBCMIMO detection, is to have high performance results while maintaining low complexity.
There are other articles that apply ABC to optimization problems with a binary domain. Wang et al. have presented a binary selection method to the ABC algorithm in [32]. Their binary encoding method employs a sigmoid function of velocity as a logical choice for binary selection. However, because they are incorporating (8) with a sigmoid function, their method has a higher complexity than (10) presented in the present paper. Salim et al. introduce a discrete bee algorithm for numerical optimization. However, their algorithm includes binary operators ((8) and (11)–(13) in [33]). They show that their algorithm has better performance than ABC; but if this algorithm is applied to a numerical optimization problem, integertobinary and binarytointeger builtin functions are required to convert all of the integer individuals of each population to binary for food source exploration, and they all have to be converted back into integer for fitness function evaluation after population modification. This procedure has to be done at least two times (during employed bee and onlooker bee phases) for all individuals in all generations, which dramatically increases the algorithm complexity. The algorithm in the present paper has two advantages over the aforementioned papers: it can be implemented for numerical optimization problems, including problems with a binary domain, and it clearly has less complexity compared with algorithms presented in other papers. DABC’s superior performance compared with other mainstream EAs is demonstrated through simulation results in Section 6.
4.3. Application of DABC to MDSTBCMIMO Joint Symbol Detection
The aim of applying DABC to the MDSTBCMIMO symbol detection problem is to minimize the Euclidian distance defined in (6). Therefore, the Euclidian distance in (6) represents the fitness function or nectar quality, and shorter Euclidian distance means better fitness. An individual of the discrete ABC algorithm corresponds to a possible solution to the joint symbol detection problem, that is, a set of conveyed symbols from the transmit devices.
In the MDSTBCMIMO system discussed in this paper, transmitted symbols are chosen from an IQ modulation such as QAM or PSK constellation diagram. We represent each of the possible points in the constellation by a unique integer in the set . The system comprises transmit devices, each device indexed by transmitting QAM symbols in a spacetime code block. Therefore, a DABC individual (a food source location) can be defined as a dimensional () integer row vector where , . For DABC, the integer vector represents the vector [] in expression (6), and the fitness function is translated accordingly. Consequently, a scout bee generates a random vector of integer numbers in as a new individual (new food source location) in line 1 of Pseudocode 2 and line 3 of Pseudocode 5; that is, , .
5. Computational Complexity
A motivation for applying the proposed nearoptimal algorithms to an MDSTBCMIMO problem is their low computational complexity. In this section, the computational complexity of DABC for MDSTBCMIMO symbol detection is compared with that of ZF, MMSE, SD, SDR, EDA, BBO, and GA. The computational complexity of the exhaustive search (an implementation of the ML detector) is , so exhaustive search is usually impractical for realtime operations of symbol detection. A number of suboptimal detection schemes with better computational complexity have been presented in the literature.
The worstcase complexity of SD is exponential, and its expected complexity depends on the problem size and SNR [5]. SD has high complexity of [34] at low SNRs, where . However, it has polynomial complexity, often roughly cubic complexity, at high SNRs [5]. MMSE is one of the suboptimal detectors that involves inverting a matrix, and its computational complexity is [35]. The computational complexity of SDR [36] in each iteration is where stands for the number of transmit antennas.
Typically, the computational complexity of populationbased algorithms is analyzed in terms of the number of fitness function evaluations, which in this paper would be (6). One important reason is that their complexity is highly dependent on their implementation and coding efficiency. The number of function evaluations in, BBO, GA, and EDA, is the same and equal to , where and represent the total number of generations and the population size, respectively [9]. The reason is that in all these algorithms, every individual is evaluated just once during one generation.
In the ABC algorithm however, there is more than one fitness function evaluation for each individual during a generation. During one generation in ABC, during the employed bees phase, each employed bee tests a neighbor food source for its quality; thus the fitness function evaluation has to be run for the whole food source positions once. By the same token, during the onlooker bees phase there are fitness function evaluations for every food source. So the overall number of fitness function evaluations for the algorithm for these two phases would be . In the scout bees phase, the algorithm selects only one individual that exceeds trials to abandon and replaces its employed bee with a scout. Hence the first individual to exceed trials would be in the th algorithm generation. After that, in order to determine a feasible number of function evaluations the worst case has to be considered, in which after the th generation there is a maximum of one abandoned food source in every generation. As a result, the total number of fitness function evaluations for ABC would be as follows: This complexity is higher than the complexity of other aforementioned EAs. However, the complexity of DABC presented in this paper is yet less than (12) because this algorithm does not run the function evaluation procedure for all the individuals in the employed bees and onlooker bees phases due to the conditions in line 7 of Pseudocode 3 and line 9 in Pseudocode 4. This stochastic behavior prevents DABC to have a closedform number of fitness function evaluations.
6. Simulation Results
This section contains the simulation results of the proposed DABCbased detection and its comparison with other detection techniques applied to an MDSTBCMIMO system. The channels are assumed to be quasistatic, and different channels in MDSTBCMIMO are assumed to be independent. In all our simulations, it is assumed that the mobile data is transmitted in a form of 4QAM modulation from all wireless devices . For simulation experiments we assumed that each of the devices transmits the same number of symbols . Therefore, there are symbols conveyed from the transmit devices to the receiver. Each point in the plots of Figures 1, 2, 3, 4, 5, 6, 7, and 8 is a value averaged over multiple independent runs. During each simulation runs, the set of symbols transmitted in a spacetime code block are generated randomly and independently of other simulation run. Also the noise term is generated randomly and independently of other simulation trials.
(a)
(b)
(c)
(a)
(b)
(c)
(a)
(b)
(c)
(a)
(b)
(c)
(a)
(b)
In order to present a fair comparison between the EAs, they are tuned to their best performance, and they are sharing the same number of generations, population size, and initial population. Other EA parameters are kept constant during all simulation runs. A list of parameters set for each algorithm is presented in a table for each simulation result next to its figure. Moreover, GA employs a greedy selection scheme [11], and BBO uses a low complex emigrationbased migration scheme with a constant emigration curve and linear immigration curve [9].
Each point in the plots of Figures 1–8 is a value averaged over multiple independent runs. In each simulation run, the set of transmitted symbols ( in (5)), channel matrices ( in (5)), and noise ( in (5)) are generated randomly and independently of other trials. Hence, in each simulation run the received signal in (6) is set from those randomly generated variables in accordance with (5). Then the algorithms are run to search for the value of the integer vector (which represents as mentioned in Section 4.3) that minimizes . Therefore, the results averaged over different simulation runs are in fact averaged over the different channel and noise realizations and also different realizations of the algorithm’s evolution in the case of probabilistic algorithms such as GA, DABC, EDA, and BBO. This experimental setup compares different algorithms in terms of the averaged performance over different channel and noise realizations.
The simulation results in Figures 1 through 8 show the BER performance comparison between ZF, MMSE, SDR, SD, GA, EDA, BBO, and DABC detectors. The MDSTBCMIMO system configuration, , is set to , , , and for Figures 1, 2, 3, and 4, respectively. The Alamouti spacetime coding [24] is used in Figures 1, 2, and 3; but for Figure 4, a nonorthogonal four transmit antennas configuration is used for each mobile device.
EA’s shared parameters , are set to , , , and in Figures 1, 2, 3, and 4, respectively. For these figures, the total number of symbols transmitted from all users is set to 8, 10, 12, and 14; the algorithms are searching through a search space of , , , and possible solutions, respectively. For each simulation run, the pair of is selected not only to make the EAs’ results close to SD’s, but to choose the smallest possible and to reduce their computational complexity.
In these figures, three EAs (DABC, BBO, and EDA) return the closest results to the SDs. In most of the cases, DABC and SD exactly match together and seem as a united line. Observing these figures shows that the best algorithm that almost always returns the same result as the SD is DABC, followed by the BBO decoder that returns results with about 95% of SD. The third place is for EDA, followed by GA. All EAs outperform other suboptimal detection methods in all the five figures and can meet the optimal result by searching through a much smaller set of individuals by selecting a decent pair of .
Figures 1, 2, and 3 indicate that in the MDSTBCMIMO, DABC has significant better BER performance than ZF, MMSE, SDR, GA, EDA, and BBO, while it closely matches SD. In Figure 3, for example, although DABC and other EAs are searching over approximately individuals over the entire search space of possible solutions, DABC’s BER performance meets the nearoptimal SD. According to this figure, DABC matches SD, where BBO, EDA, and GA require about 0.1 dB, 0.5 dB, and 0.6 dB more SNRs than BBO to achieve BER of , respectively.
The experiment in Figure 4 was performed on a nonorthogonal spacetime code, whereas the experiments in the other figures were performed on the Alamouti code (simple and orthogonal). The total number of symbols transmitted from all users in a space time code block is 12. Hence, using 4QAM, the size of the search space would be 4^{12}, while EAs only evaluate points in the search space, which is significantly smaller than the search space ML has to cover. Similar to other figures, SD and DABC have the best BER performance. In higher SNRs, GA’s performance diminishes notably, while DABC pursues the nearoptimal SD. It can be observed from the figure that BBO, EDA, and GA require 0.3, 0.4, and 1.5 dB less SNR than SD and DABC to achieve BER of 10^{−2}, respectively. DABC perceptibly behaves as the best detection algorithm among other suboptimal detection methods in all the four figures.
From the computational complexity point of view in EAs, finding the optimal pair of is essential in order to minimize the processing power and the required memory. According to the computational complexity order of these algorithms, with a fixed population size , more iterations until termination means more computation. Figure 5 shows the number of iterations required by each detection scheme to achieve a desirable BER. The MIMO system configuration is , it uses the Alamouti STBC and quasistatic channel, and the SNR is fixed to 8 dB. Figure 5 shows that the discrete ABC algorithm with the population size fixed to 100 is the first algorithm that can reach SD’s performance in less than 65 iterations. Other EAs whether cannot reach the SD results, or require much more iterations to converge to the SD results. This improved performance is consistently observed in several other simulations with different system configurations. As a result, not only DABC algorithm outperforms other suboptimal algorithms, it delivers better results than other wellknown EAs, such as GA, BBO, and EDA, and can reach SD.
Figures 6 to 8 show the tradeoff between the population size and the iterations required to achieve a desired BER in GA, BBO, and DABC. The MIMO system configuration is , using the Alamouti STBC and quasistatic channel, and the SNR is set to 8 dB. The detailed system configuration is given in Table 1. This tradeoff is useful from the system design point of view. If a hardware system has high processing capabilities and low memory, then the population size can set lower to get same BER performance and vice versa. (Higher and need more memory.)
Finally, we compare the number of fitness function evaluations between EAs to demonstrate the superiority of DABC over other EAs. As mentioned in Section 5, BBO, EDA, and GA all have the same number of fitness function evaluations equal to . For the parameters of a system described in Figure 2 for instance, these EAs require fitness function evaluations. The original ABC requires , while DABC requires an average of 13,600 evaluations per independent simulation run, which is comparable to other EAs, less than the original ABC and much less than the optimal ML detector with a search space of 4^{10} possible solutions. As a result, DABC would be a significantly considerable choice for joint symbol detection in MDSTBCMIMO systems.
7. Conclusion and Future Work
In this paper, a modified version of the Artificial Bee Colony algorithm is presented for the optimization problems in discrete domain and is applied to a Multi Device (MD) SpaceTime Block Code (STBC) MultiInput MultiOutput (MIMO) system. The enhancements in this algorithm have reduced its complexity, which is much less as compared with optimal ML detector. Thus it is suitable for costeffective highspeed realtime communications. In addition, compared to other evolutionary algorithms like GA, EDA, and BBO, the presented discrete ABC (DABC) detection in MDSTBCMIMO shows a significantly better performance. The proposed algorithm also has consistently better performancecomplexity tradeoffs at low SNRs, in comparison to the existing algorithms. Even at high SNRs, this algorithm has relatively good performancecomplexity tradeoffs. In conclusion, the proposed DABC is suitable for costeffective highspeed realtime communications in MIMO systems.
DABC is a good candidate for solving the same type of computationally complex problems in wireless communication because of its simplistic model, low implementation complexity, and convergence to a nearly optimal solution with a small number of iterations. In this work, our main purpose was to employ a lowcomplexity algorithm that can beat the existing suboptimal detection schemes. To this end, some of the evolutionary procedures of the algorithm have been modified.
Currently, the authors are working on hybridizing DABC with other EAs to take the best aspects of each algorithm for specific applications. The authors are also working on applying DABC to other types of optimization problems, such as multiobjective or constrained optimization problems.
Appendix
Food Source Selection in a Neighborhood in DABC
In (8), is a random number between . We choose the minimum and maximum of this interval and apply it to (8). The minimum is ; hence: The maximum is ; therefore: Putting (A.1) and (A.2) together and using a random integer number generator function that returns a random integer between and , for the integer values of we obtain
Acknowledgment
This work was supported in part by a discovery grant from the National Science and Engineering Research Council (NSERC) of Canada.