# Quantum Bee Colony Optimization and Non-dominated Sorting Quantum Bee Colony Optimization Based Multi-relay Selection Scheme.

1. IntroductionIt is well known that the relay nodes play an import role in conventional cellular networks to help enlarge the coverage of base station or increase the overall cell throughput in 3GPP Long Term Evolution-Advanced (LTE-A) [1]. Relaying is an effective and important technology to overcome the limitation of cell coverage and especially cell-edge users' throughput [2-3]. Relay technology is very important in other cooperative networks, such as ad hoc networks. In order to exploit the advantage of the relay node deployment in the wireless networks, relay selection, power and bandwidth allocation have been investigated in the literature, in which relay selection is the key issue of the radio resource management (RRM) in relaying systems. Most of the relay selection researches are based on certain function of channel state information (CSI), which are physical distance, path loss or end-to-end SNR [4]. In this scenario, the receiver knows all of the CSI between the source and the relay and all of the CSI between the relay and the destination thus chooses one relay according to certain function of CSI [5-6]. However, selecting single relay in wireless networks may result in imbalance of resource utilization, and moreover, the "emergence" diversity gain among multiple relays cannot be achieved. Furthermore, single assisted relay may have the disadvantage of heavy load. In order to solve this problem, [7] has proposed a load-based relay selection scheme. However, single relay selections cannot avoid the fading effect of wireless channels, therefore multi-relay selection schemes are widely researched. It has the ability of both improving the stability of network and maximizing the end-to-end SNR or power efficiency. It is especially useful for wireless relay networks with multiple relays and many complex constraints. However, with the number of relays increasing, the network may have the problems of much more interference and resource crashing. Thus, how to choose a set of suitable relays is very important. As it is well known, energy efficiency or power efficiency is of great importance in green communications, which is also the same in the relay networks. In [8], several relay selection strategies for multi-relay scenarios are proposed, which takes the instantaneous error rate and fast fading channels into consideration. In [9], a novel relay selection scheme considering energy-efficiency is proposed and the relay node which has the best energy efficiency is selected. In [10], a relay ordering based relay selection scheme which considers end-to-end SNR and end-to-end power efficiency is proposed. But the scheme for energy-efficiency only can get a sub-optimal solution, and the results of simulation in [10] illustrate that the solution given by the relay ordering scheme has a large gap compared with exhaustive search scheme.

The relay selection problems in [7-10] are all single objective optimization multi-relay selection problems, which can be modeled as '0-1' optimization problem. Intelligence algorithm can be used to solve the problems, which can get an approximate optimal solution. Some classical intelligence algorithms are widely researched and applied, like particle swarm optimization (PSO) [11]. Quantum genetic algorithm (QGA) is the combination of quantum theory and genetic algorithm therefore has the advantage of faster convergence rate, stronger searching abilities, less computing time. Quantum particle swarm optimization (QPSO) which combines PSO with quantum computing theory is a novel swarm intelligence algorithm, which has a better performance for multi-relay selection problem [12]. Quantum bee colony optimization (QBCO) is a novel swarm intelligence algorithm for solving cognitive radio spectrum allocation problem which is proposed in [13]. Quantum theory has great efficiency and effectiveness in the intelligence algorithm domain, therefore QBCO is designed to solve single objective multi-relay selection problems.

Single objective optimization based multi-relay selection schemes have some limitations. In SNR maximization based multi-relay selection scheme, we can only obtain the solution which has the maximum SNR value. However, it can't reveal the relationship between power consumption and the marginal effect of the SNR. It may be energy inefficiency in high SNR region. Decreasing the transmission power only has little impact on the obtained SNR. Therefore we can decrease the consumed power with the cost of little SNR degradation. In power efficiency maximization based multi-relay selection scheme, we can only obtain the solution which has the maximum power efficiency value. However, such optimization doesn't consider the obtained SNR and data rate. It may give the largest power efficiency with very low obtained SNR. Obviously, such solution can't fulfill the SNR requirement for quality of service (QoS) guarantee in the transmission. In this optimization problem, the optimal power efficiency solution is not "optimal transmission". In order to overcome the disadvantage of single objective optimization based multi-relay selection schemes, we propose the multi-relay selection problems considering multiple objectives (SNR maximization and power consumption minimization or SNR maximization and power efficiency maximization). Multi-objective optimization problems represent an important class of real-world problems. In principle, they are very different from the single objective optimization problems. In single objective optimization, the goal is to obtain the best design or decision, which is usually the global minimum or global maximum on a particular performance indicator depending on the optimization problem of minimization or maximization. In multi-objective optimization, however, there does not exist one solution which is the best with respect to all objectives. Typically such problems involve tradeoffs. In a typical multi-objective optimization problem, there exist a set of solutions which are superior to the rest of solutions in the search space when all objectives are considered but are inferior to other solutions in the space in one or more objectives are considered. The solutions are known as Pareto front solutions or non-dominated solutions. The rest of the solutions are known as dominated solutions.

A number of multi-objective evolutionary algorithms have been proposed in literature, such as classical NSGA [14], NSGA-II [15], SPEA [16] and SPEA2 [17]. [18-23] have done reaserch on multi-objective algorithms based on classical algorithms. In order to improve convergence and diversity of solutions of Pareto front, [18] improves solutions diversity of Pareto front of a well known multi-objective optimization algorithm NSGA-II. [19] proposes a new solution based on multi-objective optimization using the genetic algorithm NSGA-II for security, QoS, and energy efficiency in Wireless Sensor Networks (WSNs). Resource constraints as well as QoS requirements are respected through use of optimal security level based on evolutionary strategy. [20] has researched the multi-objective combinatorial optimization in the design of fiber-based distribution networks. In [21], a multi-objective two-nested genetic algorithm is used to solve clustering homogeneous wireless sensor networks. [22] proposes a new Multi-Objective Optimization Algorithm Based on Non-Dominated Sorting and Bidirectional Local Search (NSBLS) to solve high-dimensional big data. [23] considers many-objective problems and proposes an early-developed and computationally expensive strength Pareto based evolutionary algorithm by introducing an efficient reference direction based density estimator, a new fitness assignment scheme, and a new environmental selection strategy. The majority of multi-objective algorithms make use of the Pareto dominance concept to assign a single fitness value for each individual in the population. This is used to select a set of Pareto front solutions. The diversity of solutions in the Pareto front to cover different tradeoffs of the problem objectives and the distance to the actual front (optimal solutions) are two main issues that should be considered carefully and are affected by the fitness assignment and Pareto front individual selection techniques. Using an external memory (archive) to store non-dominated solutions found during the search process is a common approach to maintain the Pareto front. In multi-objective optimization, it is also important to choose the evolutionary method. Due to the effectiveness of QBCO [13], we choose QBCO as the evolutionary method. Therefore we propose the NSQBCO to solve the multi-objective multi-relay selection problem (SNR maximization and power consumption minimization or SNR maximization and power efficiency maximization).

The rest of the paper is organized as follows. The network model and problem illustration are described in Section 2. The QBCO based single objective multi-relay selection scheme is proposed in Section 3. Section 4 gives the NSQBCO based multi-objective multi-relay selection scheme. Section 5 analyses the simulation. Section 6 concludes the paper.

2. Network model and problem statement

In this section, a cooperative wireless relay system model is considered, which consists of one transmitter for transmission and one receiver for reception and R relays for cooperation as described in Fig. 1.

With relay selection schemes, a set of relays is selected from the R potential relays which maximize SNR or power efficiency proposed later. There is no direct link between the transmitter and the receiver. However, the results can be applied to the case with a direct link straightforwardly.

It is assumed that each relay has only one antenna which can be used for both transmission and reception. Denote [f.sub.i] as CSI from s to the i-th relay and [g.sub.i] as CSI from the i-th relay to y. It is assumed that [f.sub.i] and [g.sub.i] are full known by the i-th relay, all CSI [f.sub.1], [f.sub.2],..., [f.sub.R] and [g.sub.1], [g.sub.2],... , [g.sub.R] are full known by the receiver. Assume that all CSI are normalized independent identical distribution (i.i.d.) Rayleigh random variables. P denotes the transmission power of the transmitter, and [P.sub.i] denotes the transmission power of the i-th relay. Note that power control is not considered in the model, that is to say, the transmitter or relay cannot save power and sponsor the transmissions with better channels. A relay only has two choises, either cooperate or not cooperate at all. A two-step AF protocol is used to forward information without decoding. The transmitter sends the signal [square root of P]z to all the other relay nodes in the first transmission process, where z is the signal to transmit and z is normalized as E[|z|.sup.2] = 1. The i-th relay (if chosen) amplifies the signal received and forwards it to the receiver with power [P.sub.i] in the second transmission process.

Firstly, existing relay selection scheme which selecting one relay based on SNR is reviewed and the cumulative distribution function (CDF) is derived and simulated.

Considering the system model in Fig. 1, the received signal at the i-th relay is

[x.sub.i] = [f.sub.i][square root of P]z + [v.sub.i], (i = 1,..., R) (1)

where [v.sub.i] represents the additive white Gaussian noise (AWGN) with zero-mean and unit-variance at the i-th relay. The relay scales the signal received from the transmitter, then forwards it to the receiver, therefore the received signal at the receiver can be represented by

[mathematical expression not reproducible] (2)

where w represents the AWGN with zero-mean and unit-variance at the receiver. The overall SNR of the transmission is

[mathematical expression not reproducible] (3)

In the relay selection, the relay which has the largest SNR will be chosen to amplify and forward the transmitting data from the source to the destination. Set P = [P.sub.i] = 1 for simplicity, and considering only one relay is selected among R available relays, the CDF of transmission SNR can be derived as

[mathematical expression not reproducible] (4)

where [[gamma].sub.0] represents the SNR value of between the transmitter and the receiver, and [I.sub.1] (*) represents the first order modified Bessel function of the second kind.

Proof: The CDF of F([[gamma].sub.0]) can be derived from single transmitter, single relay and single receiver scenario, which is defined as [F.sub.S,R,D]([[gamma].sub.0]). First, we calculate the CDF of [F.sub.S,R,D] ([[gamma].sub.0]):

[mathematical expression not reproducible]

where X = [[GAMMA].sub.s,R] = [|f|.sup.2] represents SNR from transmitter to relay, Y = [[GAMMA].sub.R,D] = [|g|.sup.2] represents SNR from relay to receiver. Therefore, [P.sub.r] (X [less than or equal to] [[gamma].sub.0]) = 1 -[e[-[gamma].sub.0]], [P.sub.r] (Y [less than or equal to] [[gamma].sub.0]) = 1 -[e[-[gamma].sub.0]] and we use two regions to calculate the integration, that is, [D.sub.1] = {X < [[gamma].sub.0]} U {Y < [[gamma].sub.0]}, [D.sub.2] = D - [D.sub.1]. So the result of the integral can be written as [F.sub.S,R,D] = [F.sub.1] + [F.sub.2].

[mathematical expression not reproducible]

Define [X.sub.1] = X - [[gamma].sub.0],[Y.sub.1] = Y - [[gamma].sub.0]. Since-XY/X+Y+1= [[gamma].sub.0], we can obtain [X.sub.1][Y.sub.1]= [[gamma].sup.0.sub.0] + [[gamma].sub.0]. Therefore,

[mathematical expression not reproducible]

The above integral has the same form with the following:

[mathematical expression not reproducible]

Assume z= [e.sup.t], and place it into the above equations, we can obtain

[mathematical expression not reproducible]

Now it is obvious that the above can be calculated by the integral of [I.sub.1] (x). Therefore,

[mathematical expression not reproducible]

where a = [[gamma].sup.2.sub.0] + [[gamma].sub.0] and b = 1. Use the above result, [F.sub.2]([[gamma].sub.0]) can be written as

[mathematical expression not reproducible]

So

[mathematical expression not reproducible]

For the single relay selection based on SNR, the relay which has the maximum end-to-end SNR is chosen. Therefore, [Y.sub.max] = {[Y.sub.1],[Y.sub.2],..., [Y.sub.R]}. The CDF of [Y.sub.max] can be represented as

[mathematical expression not reproducible]

Considering theory simulation and Monte-Carlo simulation, the CDF varies with SNR which select one best relay in cooperative relay networks based on SNR is simulated. Fig. 2 shows the result when the relay number is 1, 5 and 10 respectively. From Fig. 2, it is obvious that the CDF obtained by theory has almost the same performance with the CDF obtained by

Monte-Carlo, which demonstrates the accuracy of our derivation.

From Fig. 2, we can see that although the diversity gain increases with relay number, however it is rather limited. Fig. 3 presents the simulation of CDF which consider selecting one relay and multiple relays (the best multi-relay selection is obtained through exhaustive search scheme, where all of the [2.sup.R] solutions are calculated and the optimal one is chosen).

From Fig. 3, it is obvious that selecting multiple relays has much more diversity compared with selecting one relay. Besides, if only one relay is selected, it is likely to have heavy load problem. In order to explore more diversity of selecting multiple relays, the multi-relay selection problems and schemes are proposed to select a set of relays from the R potential relays. Since the data to transmit is distributed among multiple relays, the multiple relays can overcome the problem of heavy load transmission of the single relay.

In multi-relay system, the transmitter sends the signal [square root of P]z to the i-th relay in the first transmission process, while in the second transmission process, the i-th relay amplifies its signal received from transmitter by [mathematical expression not reproducible] (therefore the transmission power is [a.sub.i] [.sup.2][P.sub.i]) and then forwards it, where [a.sub.i] represents whether the i-th relay is chosen or not. If [a.sub.i] = 1, the i-th relay is chosen, otherwise [a.sub.i] = 0.

Assume that the relays transmit at the same time, that is to say, synchronization problem is not considered. The angle [[theta].sub.i] = -(arg [f.sub.i] + arg [g.sub.i]) represents the phase of the signal. The received signal by the i-th relay's is defined as |[f.sub.i]|[square root of P]s + [v.sub.i].Thus, the received signal is

y = [square root of P] [R.summation over (i=1)][[a.sub.i]|[f.sub.i][g.sub.i]|[square root of [P.sub.i]]]/[[square root of 1+[|[f.sub.i]|.sup.2]P]]s+[R.summation over (i=1)][[a.sub.i]|[g.sub.i]|[square root of [P.sub.i]]]/[[square root of 1+[|[f.sub.i]|.sup.2]P]][u.sub.i]+w (5)

where w is AWGN at the receiver and [mathematical expression not reproducible] while [v.sub.i] is the AWGN at the i-th relay. w and [v.sub.i] are assumed to be i.i.d. complex Gaussian random variables with zero-mean and unit-variance. It is obvious that [u.sub.i] and [v.sub.i] have the same distribution. The average SNR of the communication system is

[gamma] = P[([R.summation over (i=1)][[a.sub.i]|[f.sub.i][g.sub.i]|[square root of [P.sub.i]]]/[[square root of 1+[|[f.sub.i]|.sup.2]P]]).sup.2]/(1+[R.summation over (i=1)][[a.sup.i.sub.2][|[g.sub.i]|.sup.2][P.sub.i]]/[1+[|[f.sub.i]|.sup.2]P]) (6)

The single objective multi-relay selection optimization problem considering SNR can be written as

[mathematical expression not reproducible] (7)

In wireless networks, energy expenditure is also an emergency problem, so energy-saving communication are widely researched. It is easy to see that the power of all selected relays, [P.sub.total] = [k.summation over (i=1)] [P.sub.i] increases with the selected relays number k. Define the power efficiency as the ratio of [gamma] to all of the transmission powers (the transmitter power and the relay power). The single objective multi-relay selection problem considering power efficiency can be written as

[mathematical expression not reproducible] (8)

Assume that the receiver knows all CSI, this problem is equivalent of solving the problem of the SNR or power efficiency maximization, which is like the problem in [24]. But, here the power control problem is not taken into consideration. Instead, each relay has only two choices: to cooperate with full power or not to cooperate at all. Since every relay has two choices, the problems considering SNR or power efficiency are general 0-1 optimization problems. Exhaustive search scheme has the ability to solve the problem, but the complexity increases exponential with relay number. Therefore QBCO is used to solve the multi-relay selection problems to get a better solution, which will be presented in Section 3.

Since the SNR target increases with the power, the SNR target maximization and the power consumption target minimization are contradictive. Considering the SNR and power consumption simultaneously, multi-objective relay selection problem is proposed, which is in the following

[mathematical expression not reproducible] (9)

Also, the SNR target increases with the power increasing, while the power efficiency may be decrease with the increased power, the SNR target maximization and the power efficiency target maximization are also contradictive, i.e., they cannot get the largest value with the same relay selection scheme. Considering the SNR and power efficiency simultaneously, another multi-objective multi-relay selection problem is proposed, which is in the following

[mathematical expression not reproducible] (10)

Exhaustive search can be used to solve multi-objective multi-relay selection problem (9) and (10), but the complexity is intolerable, that is to say it cannot be used in technology application. In this paper, NSQBCO is proposed to solve the multi-objective multi-relay selection problems, which will be illustrated in Section 4.

3. Single objective multi-relay selection scheme

The relay ordering multi-relay selection schemes proposed in [10] is reviewed and then we propose the single objective multi-relay selection schemes based on QBCO.

The process of relay ordering schemes can be illustrated in the following:

Step1: Order the available relays according to certain functions (Best Worse Channel Selection, which can be written asmin {P|[f.sub.i]|, [P.sub.i]|[g.sub.i]|}, Best Harmonic Mean Selection, which can be written as [([P.sup.-1][|[f.sub.i]|.sup.-2] + [P.sup.-1.sub.i] [|[g.sub.i]|.sup.-2]).sup.-1], SNR-based Selection, which can be written as|[f.sub.i][g.sub.i]|P[P.sub.i]/1+[|[f.sub.i]|.sup.2] P+[|[g.sub.i]|.sup.2][P.sub.i]), therefore get an ordering ([x.sub.1], [x.sub.2],..., [x.sub.R]) of (1,2,..., R). That is to say, if i < j, then relay [x.sub.i] is prior to relay [x.sub.j]. In other words, if relay [x.sub.i] is not selected, relay [x.sub.j] should not be selected either.

Step2: Calculate the SNR, [gamma]([x.sub.1]), [gamma]([x.sub.1], [x.sub.2]),..., [gamma]([x.sub.1], [x.sub.2],..., [x.sub.R]).

Step3: The destination finds the [x.sub.i] such that [gamma]([x.sub.1],[x.sub.2],...,[x.sub.i]) is the largest among the [gamma]([x.sub.1]), ..., [gamma]([x.sub.1], [x.sub.2],..., [x.sub.i]),..., [gamma]([x.sub.1], [x.sub.2],..., [x.sub.R]). The relays which take part in the cooperative communication can be represented as ([x.sub.1], [x.sub.2],..., [x.sub.i]).

For power efficiency problem, [10] modifies the relay ordering schemes, that is instead of selecting the [x.sub.i] which has the largest value of [gamma]([x.sub.1]),..., [gamma]([x.sub.1], [x.sub.2],..., [x.sub.i]) [gamma]([x.sub.1], [x.sub.2],..., [x.sub.R]), choose the smallest [x.sub.i] such that [gamma]([x.sub.1], [x.sub.2],..., [x.sub.i]) > [gamma]([x.sub.1], [x.sub.2],...,[x.sub.i+1]). In other words, we examine the following values [gamma]([x.sub.1]),..., [gamma]([x.sub.1], [x.sub.2],..., [x.sub.i]),..., [gamma]([x.sub.1],[x.sub.2],...,[x.sub.R]), till the SNR values stop increasing. This makes the "worse" relays are not selected, although some overall SNR may be lost.

It has been proposed in [10] that for wireless communication relay networks which has more than 2 relays, there exists no optimal relay ordering. So the schemes proposed in [10] are not global-optimal, that is to say, only a sub-optimal solution is obtained. So we propose the single objective multi-relay selection schemes based on QBCO.

In this paper QBCO is used to solve multi-relay selection problems, which is referred by social behaviour of quantum bees. Quantum employed bees, quantum onlooker bees and quantum scouts bees consist of the colony of quantum bees. They look for food resources (which are represented by quantum position) in an R dimensions space according to its own and its parteners' historical experiences; where R represents the optimization problem's dimension. In QBCO, quantum coding is used to represent the probabilistic state, and the quantum position can be updated by quantum rotation angle (which is defined by certain quantum bee, its local experiences and the whole quantum colony's experiences), which is similar to the bird's flying process. One quantum position can be written by a pair of numbers (a, [beta]), where [|a|.sup.2] + [|[beta]|.sup.2] = 1. The bit position is decided by the quantum position with certain functions, where [|[alpha]|.sup.2] decides the probability that the bit position is in the '0' state and [|[beta]|.sup.2] decides the probability that the bit position is in the '1' state.

The i-th quantum bee's quantum position is

[mathematical expression not reproducible] (11)

where[|[[alpha].sub.ij]|.sup.2] + [[[ss].sub.ij]|.sup.2] = 1, (j = 1,2,...,R), therefore [v.sub.i] can represent [2.sub.R] values simultaneously. For efficiency, set [[alpha].sub.ij] and [[beta].sub.ij] are real numbers and 0 [less than or equal to] [[alpha].sub.ij] [less than or equal to] 1, 0 [less than or equal to] [[beta].sub.ij] [less than or equal to] 1. Therefore [[alpha].sub.ij] =[square root of 1-[[ss].sup.2.sub.ij]], and equation (11) can be simplified as

[v.sub.i] = [[[alpha].sub.i1] [[alpha].sub.i2]... [[alpha].sub.iR]] = [[v.sub.i1] [v.sub.i2]... [v.sub.iR]] (12)

3.1 Evolutionary process of quantum employed bees

The first half of quantum bees in the quantum colony are quantum employed bees. The quantum position is mainly calculated by quantum rotation angle. In QBCO, for simplicity, the i-th (i = 1,2,...,h/2) quantum employed bee's j-th quantum position [v.sub.ij] is updated as

[v.sup.t-1.sub.ij] = abs([v.sup.t.sub.ij] x [cos[theta].sub.t+1.sub.ij] -[square root of [1 - ([v.sup.t.sub.ij]).sup.2]] x sin[[theta].sup.t+1.sub.ij]) (13)

where superscript t is the number of generations (which is also iterations), abs(*) represents the absolute function which makes quantum position in the domain[0, 1], and [[theta].sup.t-1.sub.ij] is the quantum rotation angle calculated through (15) which will be described later.

If [[theta].sup.t+1.sub.ij]=0, according to (13), [v.sup.t+1.sub.ij] = [v.sup.t.sub.ij]. To keep the diversity of the colony of quantum employed bees, if [[theta].sup.t-1.sub.ij] = 0 the quantum position [v.sup.ij] is updated in certain possibility by

[v.sup.t+1.sub.ij] = [square root of 1-[([v.sup.t.sub.ij]).sup.2]] (14)

The quantum colony consists of h /2 quantum employed bees that flies in a space of R dimensions, [x.sub.i] =([x.sub.i1], [x.sub.i2],..., [x.sub.iR]) represents the i-th quantum employed bee's bit position in the space. [v.sub.i] =([v.sub.i1], [v.sub.i2],...,[v.sub.iR]) = [[[alpha].sub.i1] [[alpha].sub.i2]... [[alpha].sub.iR]] represents the i-th quantum employed bee's quantum position and until now the best bit position (the local optimal bit position) of the i-th quantum employed bee is [p.sub.i] = ([p.sub.i1], [p.sub.i2],..., [p.sub.iR]). The global optimal bit position found by the whole quantum bee colony is [p.sub.g] =([p.sub.g1], [p.sub.g2],..., [p.sub.gR]). At each iteration, the quantum rotation angle, quantum position and bit position is updated by

[[theta].sup.t+1.sub.ij] = [e.sub.1]([p.sup.t.sub.ij] - [x.sup.t.sub.ij]) + [e.sub.2]([p.sup.t.sub.gi] - [x.sup.t.sub.ij]) (15)

[mathematical expression not reproducible] (16)

[mathematical expression not reproducible] (17)

where r is a uniform random number in the real domain[0, 1], [c.sub.1] is a constant among [0,1/ R], [[gamma].sup.t+1.sub.ij] [member of] [0,1] is uniform random number, [([v.sup.t+1.sub.ij]).sup.2] defines the selection probability of bit position state in the (t + 1)-th generation. The value of [e.sub.1] and [e.sub.2] represents the relative important degree of [p.sup.t.sub.i] and [p.sup.t.sub.g].

After updating the quantum and bit position, calculate the fitness of each quantum employed bee based on certain function, that is (7) or (8). If the fitness of [x.sup.t+1.sub.i] is better than that of [p.sup.t.sub.i], then update [p.sup.t+1.sub.i] as [p.sup.t+1.sub.i]. If the fitness of [p.sup.t+1.sub.i] is better than that of [p.sup.t.sub.g], then update [p.sup.t+1.sub.g] as [p.sup.t+1.sub.i].

3.2 Evolutionary process of quantum onlooker bees

The second half of quantum bees in the quantum colony are quantum onlooker bees. The quantum onlooker bees' quantum updating process is based on the selected quantum employed bee's quantum position. The selection possibility of the k-th (k = 1,2,..., h/2) quantum employed bee can be calculated by the following equation:

[p.sup.t.sub.k] = U([x.sub.k])/[h/2.summation over (i=1)]U([x.sub.i]) (18)

where U([x.sub.k]) represents the fitness of [x.sub.k], which is [gamma] in (7) or [eta] in (8).

At each iteration, the quantum rotation angles and velocities of the i-th (i = h/2 +1, h/2 +2,...,h) quantum onlooker bee are updated by the following equations, assume that the k-th quantum employed bee is selected as the guidance of the quantum onlooker bee:

[[theta].sup.t+1.sub.ij] = [e.sub.1]([p.sup.t.sub.kj] - [x.sup.t.sub.ij]) + [e.sub.2]([p.sup.t.sub.gj] - [x.sup.t.sub.ij]) (19)

[mathematical expression not reproducible] (20)

[mathematical expression not reproducible] (21)

After updating the quantum and bit position, calculate the fitness of each quantum onlooker bee based on certain function, that is (7) or (8). If the fitness of [x.sup.t+1.sub.i] is better than that of [p.sup.t.sub.i], then update [p.sup.t+1.sub.i] as [x.sup.t+1.sub.i]. If the fitness of [p.sup.t+1.sub.i] is better than that of [p.sup.t+1.sub.g], then update [p.sup.t+1.sub.g] as [p.sup.t+1.sub.i].

3.3 Evolutionary process of quantum scout bees

When the fitness of each quantum employed bees and quantum onlooker bees does not change in limit times, then it becomes a quantum scout bee, which drops the experiences it has owned and has the ability to find new food resources, thus the quantum position is selected randomly from the R dimensions space, while the bit position is generated according to the quantum position.

3.4 Evolutionary process of QBCO

From what we have discussed above, we can see that the proposed QBCO has the advantage of both quantum computing and bee colony optimization. The processes of quantum bee colony optimization for multi-relay selection are shown below:

Step1: Suppose that the receiver knows the CSI [f.sub.1], [f.sub.2],..., [f.sub.R] and [g.sub.1], [g.sub.2],..., [g.sub.R].

Step2: Create an initial quantum bee colony randomly based on quantum coding.

Step3: For all quantum bees, calculate the fitness (i.e., [gamma] or [eta]) for each quantum bee.

Step4: Update each quantum bee's quantum position, bit position, the local optimal position and the global optimal position of the whole quantum bee colony through the evolutionary process of quantum employed bee, quantum onlooker bee and quantum scout bee.

Step 5: If the maximum iteration is reached, stop and output the relay selection result; if not, go to step 3.

4. Multi-objective multi-relay selection scheme

Most of the relay selections in the current literatures only consider one objective, i.e., SNR or power efficiency. Considering two objectives simultaneously, i.e., SNR and power consumption (or SNR and power efficiency), we propose NSQBCO to solve the multi-objective multi-relay selection problems. NSQBCO is based on non-dominated sorting, where the entire population is sorted into various non-dominated levels. This provides the means for selecting the individuals in the better fronts, hence providing the necessary selection pressure to push the population towards the Pareto front. To maintain population diversity, the crowding distance methods adopted by NSGA-II [15] is used, which will be described in the following part.

4.1 Non-dominated sorting and crowding distance

If we want to minimize [f.sub.m](x) (m = 1,..., M), where M is the number of objectives we want to optimize, then for solutions u and v, if for all m = 1, M, [f.sub.m](u) [less than or equal to] [f.sub.m](v), and [there exists]m = 1,... M, [f.sub.m] (u) < [f.sub.m] (v), then define u dominates v, and u is a non-dominated solution, which means for all objectives, solution u is not worse than solution v and at least there exists an objective which solution u is better than solution v. If for all m = 1,..., M, [f.sub.m](u) [greater than or equal to] [f.sub.m](v), and [there exists]m = 1,..., M, [f.sub.m](u) > [f.sub.m](v), then define v dominates u, and v is a non-dominated solution. Otherwise, u and v have no dominating relationship.

The process of non-dominated sorting can be described as follows:

For each solution calculate two entities: 1) domination count [n.sub.p], the number of solutions which dominate p ; 2) [S.sub.p], this set contains all the individuals (each individual is defined as q) that are being dominated by p.

All solutions in the first non-dominated front will have their domination count as zero. Now, for each solution p with [n.sub.p] = 0 , we visit each member q of its set [S.sub.p] and reduce its domination count by one. In doing so, if for any member q the domination count becomes zero, it is put in a separate list Q. These members belong to the second non-dominated front. Now the above procedure is continued with each member of Q and the third front is identified. This process continues until all fronts are identified.

Along with convergence to the Pareto front, it is also desired that the algorithm maintains a good spread of solutions in the obtained set of solutions. We calculate the average distance of two points along each of the objectives. The crowding distance is used to maintain population diversity, and the calculation process will be described in the following.

The crowding-distance computation requires sorting the population according to each objective value in ascending order of magnitude for every front. Therefore, for each objective function, the boundary solutions (solutions with smallest and largest function values) are assigned an infinite distance value. All other intermediate solutions are assigned a distance value equal to the absolute normalized difference in the function values of two adjacent solutions. The calculation is continued with other objective functions. The overall crowding distance value is calculated as the sum of individual distance values corresponding to each objective.

From the description of non-dominated sorting and crowding distance, we can see that the solutions with better front and larger crowding distance are better than others.

4.2 Non-dominated sorting quantum bee colony optimization

The process of NSQBCO uses QBCO proposed in Section 3 as the evolutionary algorithm. The process can be summarized in the following steps:

Step 1: Initialize quantum bee colony S, including the quantum bees' quantum positions and bit positions, then evaluate each quantum bee in the colony. The number of quantum bees in S is recorded as h.

Step 2: Excute non-dominated sorting to obtain non-dominated solutions in S. Calculate the crowding distance and sort the individuals in each front in a descending order. Choose the first h /2 bees as the quantum employed bees, and update the quantum position according to (15)-(17). The global best solution [p.sub.g] is chosen from a specified top part (top 5%) of the sorted S randomly, while the local best solution [p.sub.j] (j = 1,2,..., h / 2) is chosen from the sorted S randomly. Then form a new quantum bee colony S (').

Step 3: Choose the last h /2bees as the quantum onlooker bees, and update the quantum position according to (19)-(21). The global best solution [p.sub.g] is chosen from a specified top part (top 5%) of S randomly, while the selected bee's solution [p.sub.j] (j = h/2 +1, h/2 +2,..., h) is chosen from S randomly. Then form a new quantum employed bee S ('').

Step 4: Combine S , S (') and S ('') thus form a new quantum bee colony. Execute non-dominated sorting and crowding-distance computation and choose the best h quantum bees to form a new quantum bee colony S which will take in the next generation.

Step 5: If it has reached the maximum generation T, then stop and the non-dominated solutions in the S are the Pareto front solutions. Otherwise, go to Step 2 until it has reached the maximum generation.

From the above, we can select the non-dominated solutions in the current bee colony and the parent bee colony and combine them. Then we reject the dominated solutions in the combined bee colony. Through the iteration of the evolutionary process, we can get the non-dominated solutions nearly to the true Pareto front solutions.

4.3 NSQBCO based multi-objective multi-relay selection scheme

According to the above analysis, the processes of NSQBCO based multi-objective multi-relay selection scheme are shown below:

Step 1: Assume the CSI [f.sub.1], [f.sub.2],..., [f.sub.R] and [g.sub.1], [g.sub.2],..., [g.sub.R] are obtained at the receiver before the relay selection process.

Step 2: Using NSQBCO (while one objective is [gamma] and the other is power consumption or one objective is [gamma] and the other is power efficiency) to obtain the Pareto front solutions.

Step 3: The relaying systems choose one solution from the Pareto solutions according to the tradeoff of [gamma] and power consumption or the tradeoff of [gamma] and power efficiency to take part in the cooperative transmission.

5. Simulation results and analysis

In this section, we first show the simulated [gamma] and [eta] of the proposed QBCO based multi-relay selection scheme with relay ordering multi-relay selection schemes, exhaustive search scheme, single relay selection scheme and QPSO scheme proposed in [12]. Then simulation results of NSQBCO based multi-relay selection is presented compared with exhaustive search scheme. In the simulation, all channels and noises at all of the relays and destination are normalized i.i.d. Rayleigh random variables. For QBCO, set the maximal generation to 100, h = 20, [e.sub.1] = 0.06, [e.sub.2] = 0.03,[c.sub.1] = 1/300.

5.1 QBCO based single-objective multi-relay selection scheme

First, 15 relays are adopted in the simulation and they have the same power value [P.sub.i] =0.1*P. Fig. 4 shows the simulation results which [gamma] varise with P. We can see that SNR increases with the power. From Fig. 4(a), we can also see that the three relay ordering multi-relay selection schemes perform almost the same, and multi Best Worst Channel Selection performs the worst, while the multi SNR-based Selection performs the best among the three relay ordering multi-relay selection schemes, but QBCO performs better than all of the relay ordering multi-relay selection schemes, which is the same as exhaustive search. The gap between QBCO and the other schemes is obvious. Also, it is obvious that multi-relay is much more effective than single-relay.

Then set the number of relays to 20, Fig. 4(b) shows the simulation results. From Fig. 4(b), we can see that QBCO perform better than the other relay selection schemes including QPSO scheme proposed in [12], and compared with Fig. 4(a), we can see that when the relay number increases, QBCO can find an optimal solution compared with other algorithms.

Now let we consider the power efficiency problem. Fig. 5 shows the simulation results. Fig. 5(a) considers the case when the number of relays is 15, while Fig. 5(b) considers the case when relay number is 20. Among the three relay ordering multi-relay selection schemes, the

Best Worst scheme performs the worst, while the SNR-based performs the best, which has similar performance with Fig. 4. Our scheme, QBCO, performs better than the other relay selection schemes and has the same performance as exhaustive search when R = 15.

From Fig. 4 and Fig. 5, the differences between the QBCO multi-relay selection scheme and relay ordering multi-relay selection schemes which maximize SNR or power efficiency is obvious. And if the relay number increases, the advantage of the QBCO based multi-relay selection scheme is much more obvious.

5.2 NSQBCO based multi-objective multi-relay selection scheme

Next, consider the proposed multi-objective multi-relay selection scheme. For NSQBCO, set [c.sub.1] = 1/ R, the number of quantum bees h = 20, [e.sub.1] and [e.sub.2] are random numbers among [0,1], the predefined maximum generation is 500 (T=500).

Set P = 2W and [P.sub.i] = 0.1 * P = 0.2W, taking both SNR optimization and power consumed minimization into consideration (9), the performance of all solutions (there are 15 relays in the simulation, so the number of solutions is [2.sup.15]) which is obtained through exhaustive search are plotted in Fig. 6. We can find out that there does not exist one solution which can maximize the SNR while minimize the power consumption simultaneously, that is to say, we have to look for tradeoffs. Fig. 6 also shows that there are a series of solutions that are non-dominated solutions, which are not inferior to other solutions in both of the two objectives.

Fig. 7 shows the tradeoffs between the SNR and power consumption of NSQBCO based multi-relay selection scheme and non-dominated solutions obtained by the exhaustive search (computed by the non-domination sorting of all the possible solutions in Fig. 6). Also, the solutions obtained by the SNR based relay ordering and QBCO scheme for SNR target are presented for comparison. From Fig. 7, NSQBCO based multi-relay selection scheme obtains the same solution as the exhaustive search but cost less time when the number of relays is not very large, which shows the effectiveness of NSQBCO scheme. The solutions obtained by QBCO scheme and SNR based relay ordering scheme are non-dominated solutions, which shows the effectiveness of the proposed single objective QBCO based multi-relay selection scheme. Moreover, we can see that the single-objective optimization can only obtain one solution which maximizes SNR value without considering the power consumed. However, as we can see from Fig. 7, in high SNR region, the power consumed has little effect on SNR. If the power consumed is larger than 4W, the SNR value remains almost constant. In the design of relay networks, we can decrease the consumed power with the cost of little SNR degradation. The non-dominated solutions contain the solution obtained by single objective optimization, which means the multi-objective multi-relay selection scheme has a much wider application range.

When the number of relays increases, the exhaustive search cannot be used due to algorithm complexity. But NSQBCO can still be used to solve multi-objective multi-relay selection schemes. The non-dominated solutions are presented in Fig. 8, in addition to the solutions obtained by QBCO and the SNR based relay ordering schemes. The solution obtained by QBCO scheme for SNR target is still one of the non-dominated solutions, while the solution obtained by the SNR based relay ordering scheme is not, which shows the advantage of the single objective QBCO and the multi-objective NSQBCO relay selection schemes again.

Set P = 4W and [P.sub.i] = 0.1 * P = 0.4W, considering the SNR and power efficiency simultaneously (10), all solutions (the number of solutions is [2.sup.15] = 32768) obtained through exhaustive search are plotted in Fig. 9.

It is obviously that there also does not exist one solution which can maximize the SNR as well as power efficiency, that is to say, we can only get tradeoffs. Fig. 9 also shows that there are a series of solutions which are non-dominated solutions. These solutions are not worse than the other solutions in both objectives. The proposed NSQBCO scheme aims to obtain these solutions.

Fig. 10 shows the tradeoffs between SNR and power efficiency optimization of NSQBCO based multi-relay selection scheme, in addition to the performance of non-dominated solutions obtained by the exhaustive search. The solutions obtained by QBCO scheme for SNR and power efficiency target respectively and the solutions obtained by the SNR and power efficiency based relay ordering schemes are also plotted for comparison.

The NSQBCO based multi-relay selection scheme obtains the same solution as the exhaustive search but cost less time. The solutions obtained by QBCO scheme for SNR and power efficiency target are among the non-dominated solutions, but the solutions obtained by the relay ordering schemes are really much worse, especially for power efficiency target. Moreover, if we only optimize the power efficiency target, we can only obtain the solution which has the maximum power efficiency value, but the SNR value is rather limited. Such optimization doesn't consider the obtained SNR and data rate. Obviously, such solution can't fulfill the SNR requirement for QoS guarantee in the transmission. In this optimization problem, the optimal power efficiency solution is not "optimal transmission". However, through multi-objective multi-relay selection schemes, we can choose one non-dominated solution which doesn't have the largest power efficiency value but guarantee QoS. This shows the wide application range of multi-objective optimizations.

When the number of relays increases, the exhaustive search cannot be used due to algorithm complexity. But the NSQBCO can still solve this problem efficiently. Simulations are shown in Fig. 11 when the relay number is 30.

From Fig. 11, the solutions obtained by QBCO schemes for SNR and power efficiency target respectively and the solutions obtained by the SNR and power efficiency based relay ordering schemes are also presented for comparison. The solutions obtained by QBCO scheme for SNR and power efficiency target are among the non-dominated solutions, but the solutions obtained by the relay ordering schemes are rather poor. This demonstrates that the multi-objective schemes have a wider application field compared with single objective schemes. All these present the advantage of the proposed single objective QBCO based multi-relay selection scheme and the proposed multi-objective NSQBCO based multi-relay selection scheme.

6. Conclusion

This paper has proposed multi-relay selection schemes considering single objective and multi-objective in the cooperative relay networks. Firstly, the single objective optimization problems of the best cooperative relay nodes selection for SNR maximization or power efficiency optimization are solved respectively based on QBCO schemes, and simulation results show that compared with other multi-relay selection schemes in the literature, the proposed schemes have a much better SNR or power efficiency performance. Then, considering SNR maximization and power consumption minimization or SNR maximization and power efficiency maximization simultaneously, this paper has proposed the NSQBCO based multi-objective multi-relay selection schemes, which can obtain the non-dominated solutions. Simulation results show that NSQBCO based multi-relay selection schemes obtain the same Pareto solutions as exhaustive search when the number of relays is not very large. However, when the number of relays is very large, exhaustive search cannot be used due to complexity but NSQBCO based multi-relay selection schemes can still be used to solve the problems. Besides, the solution obtained by QBCO scheme for single objective optimization is included in the non-dominated solutions, which demonstrates the wider application range of NSQBCO based multi-relay selection scheme and the effectiveness of both QBCO and NSQBCO based multi-relay selection schemes.

References

[1] 3GPP TR 36.814, "Further Advancement for E-UTRA Physical Layer Aspects," v 1.5.2, December, 2009.

[2] Laneman, J. Nicholas, D. N. C. Tse, and G. W. Wornell, "Cooperative diversity in wireless networks : Efficient protocols and outage behavior," IEEE Transactions on Information Theory, vol. 50, no. 12, pp. 3062-3080, December, 2004. Article (CrossRef Link)

[3] Nosratinia, A., T. E. Hunter, and A. Hedayat, "Cooperative communication in wireless networks," IEEE Communications Magazine, vol. 42, pp. 68-73, October, 2004. Article (CrossRef Link)

[4] V.Sreng, H.Yanik, D.Falconer, "Relayer Selection Strategies in Cellular Networks with Peer-to-Peer Relaying," in Proc. of 2003 58th IEEE Vehicular Technology Conference, pp.1949-1953, October 4-12, 2003. Article (CrossRef Link)

[5] A. Bletsas, A. Khisti, D. P. Reed, and A. Lippman, "A simple cooperative diversity method based on network path selection," IEEE Journal on Selected Areas in Communications, vol. 24, pp. 659-672, March, 2006. Article (CrossRef Link)

[6] Ding, Zhiguo, H. Dai, and H. V. Poor. "Relay Selection for Cooperative NOMA," 2016. Article (CrossRef Link)

[7] X. Lin and L. Cuthbert, "Load Based Relay Selection Algorithm for Fairness in Relay Based OFDMA Cellular Systems," in Proc. of Wireless Communications and Networking Conference, pp. 1-6, April 5-7, 2009. Article (CrossRef Link)

[8] H. Eghbali, S. Muhaidat, S. A. Hejazi and Y. Ding, "Relay Selection Strategies for Single-Carrier Frequency-Domain Equalization Multiple relay Cooperative Networks," IEEE Transactions on Wireless Communications, vol. 12, no. 5, pp. 2034-2045, May 2013. Article (CrossRef Link)

[9] T. Zhang, S. Zhao, L. Cuthbert and Y. Chen, "Energy-efficient cooperative relay selection scheme in MIMO relay cellular networks," in Proc. of IEEE International Conference on Communication Systems (ICCS), pp. 269-273, Nov, 2010. Article (CrossRef Link)

[10] Y. Jing and H. Jafarkhani, "Single and multiple relay selection schemes and their available divercity orders," IEEE Transactions on Wireless Communication, vol. 8, no. 3, pp. 1414-1423, March 2009. Article (CrossRef Link)

[11] J. Kennedy and R. Eberhart, "Discrete binary version of the particle swarm optimization," in Proc. of IEEE International Conference on Computational Cybernetics and Simulation, pp. 4104-4108, 1997. Article (CrossRef Link)

[12] Cao J, Zhang T, Zeng Z, et al, "Multi-relay selection schemes based on evolutionary algorithm in cooperative relay networks," International Journal of Communication Systems, vol.27, no.4, pp. 571-591, 2014. Article (CrossRef Link)

[13] Hongyuan Gao, Jinlong Cao, "Quantum-inspired bee colony optimization algorithm and its application for cognitive radio spectrum allocation," Journal of Central South University, vol. 43, no. 12, pp. 4743-4749, 2012.

[14] Srinivas N, Deb K, "Mutiobjective optimization using nondominated sorting in genetic algorithms," Evolutionary Computation, vol.2, no.3, pp. 221-248, 1994. Article (CrossRef Link)

[15] Deb K, Pratap A, Agarwal S, Meyarivan T, "A fast and elitist multiobjective genetic algorithm: NSGA-II," IEEE Transactions on Evolutionary Computation, vol.6, no.2, pp. 182-197, 2002. Article (CrossRef Link)

[16] Zitzler, Eckart, and L. Thiele, "Multiobjective evolutionary algorithms: a comparative case study and the strength Pareto approach," IEEE Transactions on Evolutionary Computation, vol.3, no.4, pp. 257-271, 1999. Article (CrossRef Link)

[17] Zitzler E, Laumanns M, Thiele L, "SPEA2: improving the strength Pareto evolutionary algorithm," in Proc. of evolutionary for design, optimization and control with application to an industrial problems (EUROGEN2001), pp.95-100, September 19-21, 2002. Article (CrossRef Link)

[18] Vachhani, Vimal L., V. K. Dabhi, and H. B. Prajapati, "Improving NSGA-II for solving multi objective function optimization problems," in Proc. of International Conference on Computer Communication and Informatics, January 05-07, 2016. Article (CrossRef Link)

[19] Rachedi, Abderrezak, and A. Benslimane, "Multi-objective optimization for Security and QoS adaptation in Wireless Sensor Networks," in Proc. of IEEE International Conference on Communications, May 23-27, 2016. Article (CrossRef Link)

[20] Haidine, Abdelfatteh, "Design of reliable fiber-based distribution networks modeled by multi-objective combinatorial optimization," International Journal of Communication Systems, vol.26, no.10, pp. 1227-1242, 2012. Article (CrossRef Link)

[21] Peiravi A, Mashhadi HR, Hamed Javadi S, "An optimal energy-efficient clustering method in wireless sensor networks using multi-objective genetic algorithm," International Journal of Communication Systems, vol.26, no.1, pp. 114-126, 2013. Article (CrossRef Link)

[22] F. Lin; J. Zeng; j. xiahou; S. Lin; W. Zeng; H. Lv, "Multi-Objective Evolutionary Algorithm Based On Non-Dominated Sorting and Bidirectional Local Search for Big Data," in Proc. of IEEE Transactions on Industrial Informatics, vol.PP, no.99, pp.1-1, 2017. Article (CrossRef Link)

[23] S. Jiang; S. Yang, "A Strength Pareto Evolutionary Algorithm Based on Reference Direction for Multi-objective and Many-objective Optimization," in Proc. of IEEE Transactions on Evolutionary Computation, vol. PP, no.99, pp.1-1,2017. Article (CrossRef Link)

[24] Jing, Yindi, and H. Jafarkhani, "Network beamforming using relays with perfect channel information," Information Theory IEEE Transactions on, vol.55, no.6, pp.2499-2517, 2009. Article (CrossRef Link)

Qiang Ji received his M.S. degree in Pattern Recognition and Intelligent System from Harbin Engineering University, Harbin, Heilongjiang, P.R.China, in April 2005. He is currently pursuing a doctoral degree at the College of Aerospace Science and Engineering in National University of Defense Technology. His current research interests include the design of avionics system, wireless communication networking, relay technology and artificial intelligence algorithms.

Shifeng Zhang received Ph.D. degree in Control Theory and Engineering from National University of Defense Technology in 2000. He was awarded as a Professor and Doctoral Supervisor from the College of Aerospace Science and Engineering, National University of Defense Technology in 2010. His research interests include aircraft overall design, flight dynamics, guidance and control, inertial navigation, and measurement and precision analysis.

HaoGuang Zhao received his M.S. degree in Changchun University of Science and Technology, Changchun, Jilin, P.R.China, in April 2006. He is currently pursuing a doctoral degree at the school of National University of Defense Technology. His current research interests include pattern recognition and its application in future embedded systems, radar target recognizer.

Dr Tiankui Zhang BSc(Eng) PhD is an associate professor of school of information and communication engineering at the Beijing University of Posts and Telecommunications (BUPT). He received his Ph.D. degree in information and communication engineering from BUPT in 2008 and M.S. in communication engineering from BUPT in 2003. His research interests are in wireless communication networking, green wireless networking, signal processing for wireless communications, content centric wireless networks. He had published more than 100 papers including journal papers on IEEE Journal on Selected Areas in Communications, IEEE Transaction on Communications, IEEE Communication Letters etc., and conference papers such as IEEE GLOBECOM, IEEE ICC and so on. Email: zhangtiankui@bupt.edu.cn

Jinlong Cao is a technology engineer in China Telecommunications now. He received his M.S. degree in information and communication engineering from BUPT in 2014. His research interests are in wireless communication networking, relay technology in wireless communications, cognitive radio networks, artificial intelligence algorithms and artificial neural networks. He had published more than 20 papers including journal papers on IEEE Communications Letters, International Journal of Communication Systems, Journal of Systems Engineering and Electronics, Journal of Central South University of Technology etc., and conference papers such as IEEE GLOBECOM, IEEE WPMC and so on.

Qiang Ji (1*), Shifeng Zhang (1), Haoguang Zhao (1) , Tiankui Zhang (2), Jinlong Cao (2)

(1) College of Aerospace Science and Engineering, National University of Defense Technology Hunan, China

[e-mail: heujq@sina.com]

(2) School of Information and Communication Engineering, Beijing University of Posts and Telecommunications Beijing, China

[e-mail: tkzhang@gmail.com]

(*) Corresponding author: Qiang Ji

Received November 8, 2016; revised April 12, 2017; accepted May 25, 2017; published September 30, 2017

doi.org/10.3837/tiis.2017.09.010

Printer friendly Cite/link Email Feedback | |

Author: | Ji, Qiang; Zhang, Shifeng; Zhao, Haoguang; Zhang, Tiankui; Cao, Jinlong |
---|---|

Publication: | KSII Transactions on Internet and Information Systems |

Article Type: | Report |

Date: | Sep 1, 2017 |

Words: | 8967 |

Previous Article: | Phase Noise Self-Cancellation Scheme Based on Orthogonal Polarization for OFDM System. |

Next Article: | The Improved Velocity-based Models for Pedestrian Dynamics. |

Topics: |