A Suboptimal PTS Algorithm Based on Particle Swarm Optimization Technique for PAPR Reduction in OFDM Systems

A suboptimal partial transmit sequence (PTS) based on particle swarm optimization (PSO) algorithm is presented for the low computation complexity and the reduction of the peak-to-average power ratio (PAPR) of an orthogonal frequency division multiplexing (OFDM) system. In general, PTS technique can improve the PAPR statistics of an OFDM system. However, it will come with an exhaustive search over all combinations of allowed phase weighting factors and the search complexity increasing exponentially with the number of subblocks. In this paper, we work around potentially computational intractability; the proposed PSO scheme exploits heuristics to search the optimal combination of phase factors with low complexity. Simulation results show that the new technique can e ﬀ ectively reduce the computation complexity and PAPR reduction.


INTRODUCTION
Orthogonal frequency division multiplexing technique (OFDM) is a multicarrier modulation technology which can decrease the effect of thenoise andinterferences efficiently.Meanwhile, it has many advantages, such as: senior band efficiency and less impact of intersymbol inference.The high peak-to-average power ratio (PAPR) is the main drawback of the OFDM system, in which the OFDM transmitters require expensive linear amplifiers with wide dynamic range.Moreover, the amplifier nonlinearity will cause intermodulation products resulting in unwanted out-of-band power and increased interference.
OFDM is an attractive technique for achieving high bit rate wireless data transmission in frequency-selective fading channels [1].Recently, many schemes of reduction in reductions PAPR have been proposed for OFDM system, as clipping [2] and peak windowing, block coding [3], nonlinear companding transform schemes [4,5], active constellation extension [6], selective mapping [7,8], and partial transmit sequences (PTSs) [9][10][11][12][13][14][15][16][17], which are the most attractive ones due to good system performance and low complexity.Among these methods, PTS scheme is the most efficient approach and a distortionless scheme for PAPR reduction by optimally combining signal subblocks.In PTS technique, the input data block is broken up into disjoint subblocks.The subblocks are multiplied by phase weighting factors and then added together to produce alternative transmit containing the same information.The phase weighting factors, whose amplitude is usually set to 1, are selected such that the resulting PAPR is minimized.The number of allowed phase factors should not be excessively high, to keep the number of required side information bits and the search complexity within a reasonable limit.However, the exhaustive search complexity of the ordinary PTS technique increases exponentially with number of subblocks, so it is practically not realizable for a large number of subblocks.To find out a best weighting factor is a complex and difficult problem.In this paper, we present a novel approach to tackle the PAPR problem to reduce the complexity based on the relationship between the phase weighting factors and the subblock partition schemes.The rest of this paper is organized as follow.In Section 2, definition of PAPR of OFDM system and the principles of PTS techniques are introduced.The particle swarm optimization (PSO) algorithm-based PTS OFDM system is examined in Section 3. The results of simulation are discussed in Section 4 and some conclusions for the proposed scheme are drawn in Section 5.

OFDM systems and peak-to-average power ratio (PAPR)
In an OFDM system with N subcarriers, the complex baseband representation of an OFDM signal is expressed as where T is an input symbol sequence and t stands for a discrete time index.The PAPR of the OFDM signal sequence, defined as the ratio of the maxim power to the average power of the signal, can be expressed by where E[•] denotes the expected value [18].

Optimum partial transmit sequence-OFDM system model
The principle structure of PTS method is shown in Figure 1 as that in [15].In PTS approach, the input data block is partitioned into disjoint subblocks.Each subblock is multiplied by a phase weighting factor, which is obtained by the optimization algorithm to minimize the PAPR value.We define the data block as a vector X = [X 1 X 2 • • • X N ] T , where N denotes the number of subcarriers in OFDM frame.Then, X is partitioned into M disjoint subblocks represented by the vector X i , i = 1, 2, . . ., M such that ( Here, it is assumed that the clusters X i consist of a set of subblocks with equal sizes.Then, the goal of the PTS approach is to form a weighted combination of the M subblocks which is written as where W i , i = 1, 2, . . ., M. The phase weighting factor can be chosen freely within [0,2π).In general, the selection of the phase weighting factors is limited to a set with finite number of elements to reduce the search complexity.After transforming to the time domain, the new time-domain vector becomes These partial sequences are independently rotated by phase weight factors W i , i = 1, 2, . . ., M. The optimal phase weighting factor W i that minimizes the PAPR can be obtained from a comprehensive simulation of all possible 2 M−1 combination.The objective of the PTS technique is to choose a phase weighting vector W = {W 1 , W 2 , . . ., W i } to reduce the PAPR of X , and the optimum parameters for an OFDM symbol can be by The known subblock partitioning can be classified into three categories.The first and the simplest category is called adjacent method which allocates N/M successive symbols to the same subblock.The second category is based on interleaving.In this method, the symbols with distance M  areallocated to the same subblock.The last one is called random partitioning method in which the input symbol sequence is partitioned randomly.The random partitioning is known as to have the best performance in PAPR reduction [16].It is well known that the PAPR performance will be improved as the number of subblocks M is increased for OPTS technique, optimum PAPR can be found after searching 2 M−1 computation if the number of subblock is M. A preset threshold can be used to reduce the computational complexity.We search the PAPR values through phase optimizer and the search is stopped once the PAPR drops bellow the preset threshold.By this way, the computational complexity can be significantly reduced.

PARTICLE SWARM OPTIMIZATION-BASED PTS
Basically, the PSO [19][20][21][22][23][24][25][26][27] technique-based PTS technique described below can be implemented by appropriately changing the optimization for block W in Figure 1.In this context, the population is called a swarm and the individuals are called particles.Resembling the social behavior of a swarm of bees to search the location with the most flowers in a field, the optimization procedure of PSO is based on a population of particles which fly in the solution space with velocity dynamically adjusted according to its own flying experience and the flying experience of the best among the swarm.
Figure 2 shows the flow chart of a PSO algorithm.During the PSO process, each potential solution is represented as a particle with a position vector x, referred to as phase weighting factor and a moving velocity represented as W and v, respectively.Thus for a K-dimensional optimization, the position and velocity of the ith particle can be represented as ), respectively.Each particle has its own best position W P i = (W i,1 , W i,2 , . . ., W i,K ) corresponding to the individual best objective value obtained so far at time t, referred to as pbest.The global best (gbest) particle is denoted by W G =(W g,1 , W g,2 , . . ., W g,K ), which represents the best particle found so far at time t in the entire swarm.The new velocity v i (t + 1) for particle i is updated by where v i (t) is the old velocity of the particle i at time t.Apparent from this equation, the new velocity is related to the old velocity weighted by weight w and also associated to the position of the particle itself and that of the global best one by acceleration factors c 1 and c 2 .The c 1 and c 2 are therefore referred to as the cognitive and social rates, respectively, because they represent the weighting of the acceleration terms that pull the individual particle toward the personal best and global best positions.The inertia weight w in ( 7) is employed to manipulate the impact of the previous history of velocities on the current velocity.Generally, in population-based optimization methods, it is desirable to encourage the individuals to wander through the entire search space, without clustering around the local optima, during the early stage of the optimization.
A suitable value for w(t) provides the desired balance between the global and local exploration ability of the swarm and, consequently, improves the effectiveness of the algorithm.Experimental results suggest that it is preferable to initialize the inertia weight to a large value, giving priority to global exploration of the search space, linear decreasing w(t) so as to obtain refined solutions [20][21][22].For the purpose of intending to simulate the slight unpredictable component of natural swarm behavior, two random functions r 1 and r 2 are applied to independently provide uniform distributed numbers in the range [0, 1] to stochastically vary the relative pull of the personal and global best particles.Based on the updated velocities, new position for particle i is computed according the following equation: The populations of particles are then moved according to the new velocities and locations calculated by ( 7) and ( 8), and tend to cluster together from different directions.Thus, the evaluation of each associate fitness of the new population of particles begins again.The algorithm runs through these processes iteratively until it stops.In this paper, the current positioncan be modified by [24] where w max is the initial weight, w min is the final weight, iter max is maximum number of iterations, and iter is the current iteration number.The procedures of standard PSO can be summarized as follows.
Step 1. Initialize a population of particles with random positions and velocities, where each particle contains K variable.
Step 2. Evaluate the fitness values of all particles, let pbest of each particle and its objective value equal to its current position and objective value, and let gbest and its objective value equal to the position and objective value of the best initial particle.
Step 3. Update the velocity and position of each particle according to (7) and (8).
Step 4. Evaluate the objective values of all particles.
Step 5.For each particle, compare its current objective value with the object value of its pbest.If current value is better, then update pest and its object value with the current position and objective value.Furthermore, determine the best particle of current warm with the best objective values.If the objective value is better than the object value of gbest, then update gbest and its objective value with the position and objective value of the current best particle.Original OFDM PSO Step 6. Termination criteria: if a predefined stopping criterion is met, then output gbest and its objective value; otherwise go back to Step 3.

SIMULATION RESULTS AND DISCUSSIONS
To evaluate and to compare the performance of the suboptimal PTS, numerous computer simulations have been conducted todetermine the PAPR improvements.QPSK modulation is employed with N = 256 subcarriers.The phase weighting factors W = [0, 2π) have been used.In order to generate the complementary cumulative distribution function (CCDF) [18] of the PAPR, 10000 random OFDM frames have been generated.The sampling rates for an accurate PAPR need to be increased by 4 times.The cumulative distribution function (CDF) of the PAPR is one of the most frequently used performance measures for PAPR reduction techniques.The CDF of the amplitude of a signal sample is given by CDF = 1 − exp(PAPR 0 ).In the performance comparison, the parameter of CCDF is defined as In Figure 3, some results of the CCDF of the PAPR are simulated for the OFDM system with 256 subcarriers, in which M = 16 subblockemploying random partition and the phase weight factor W = {±1} M uniformly distributed random variables are used for PTS.As we can see that the CCDF of the PAPR is gradually promoted upon increasing the numbers of generations due to the limited phase weighting factor.As the numbers of generation are increased, the CCDF of the PAPR has been improved.For a generation G n = 40, we can see that the PSO-based PTS technique is capable of attaining a near OPTS technique performance, when P r (PAPR > PAPR 0 ) = 10 −3 .
In Figure 3, we compare the PAPR performance of different numbers of particle generations G n for c 1 = c 2 = 2. Basically, the PAPR performance is improved with G n increasing.However, the degree of improvement is limited when G n is above 40.On the other hand, the computational complexity is increased with G n .Only a slight improvement is attained for increasing G n = 20 to 40.The computational complexity of G n = 40 is double of that of G n = 20.Hence, based on the trade-off between the PAPR reduction and computational complexity, G n = 20 is a suitable choice for our proposed PSO-based PTS technique.
Figure 4 shows the simulated results of the PSO-assisted PTS technique, in comparison against normal OFDM for number of subblocks M. M is one of value in the set {2, 4, 8, 16, 32}.In particular, the PAPR of an OFDM signal exceeds 12 dB for 10 −3 of the possible transmitted OFDM blocks.However, by introducing PTS approach with M = 16 clusters partition with phase factors limited to W = 4, the 10 −3 PAPR reduces to 7.5 dB.In short, new approach can achieve a reduction of PAPR by approximately 3.5 dB at the 10 −3 PAPR.Thus, the performance of the techniques is better for larger M since larger numbers of vectors are searched for larger M in every update of the phase weighting factors.Moreover, it can be observed that probability of very high peak power has been increased significantly if PTS techniques are not used.As the number of subblocks and the set of phase weighting factor are increased, the performance of the PAPR reduction becomes better.However, the processing time gets longer because of much iteration.From Figure 4, as expected, the improvement increases as number of clusters increases.Thus, using the PSO technique, we can obtain better results than presented previously.
The subblock partition for proposed suboptimal method involves dividing the subblocks into multiple disjoint subblocks.Therefore, determining which subblock partition method produces the best performance is important.Figure 5 shows that the subblock partition for proposed suboptimal method involves three dividing subblocks: adjacent method, interleaving method, and pseudorandom method.In the viewpoint of PAPR reduction, pseudorandom subblock partitioning has better performance than others.
In Figure 6, for a fixed number of clusters, the phase weighting factor can be chosen from a larger set of {2, 4, 8, 16}.It is shown that the added degree of freedom in choosing the combining phase weighting factors provides an additional reduction.When the number of phase weighting factor W = 2 and number of subblocks M = 4, PAPR can be reduced about 2.78 dB at 10 −4 from 12 dB to 9.22 dB.When W = 4 and M = 4, at 10 −4 PAPR can be reduced about 4.2 dB from 12 dB to 7.8 dB.As the number of subblocks and the set of phase weighting factor are increased, the performance of the PAPR reduction becomes better.However, the processing time gets longer because of much iteration.
In this section, a threshold is also applied to reduce calculation complexity and is calculated from the CCDF equation, which has given optimal threshold for the number Original OFDM of subblocks as follows.When N subcarriers and M sublocks are assumed, the probability that the PAPR will exceed certain of PAPR 0 is represented [18] as From (11), PAPR 0 ofthreshold ζ, which is satisfied with given probability CCDF, can be represented as which p is pr{ }.The CCDF, P r (PAPR > PAPR 0 ) for M = 8, is shown in Figure 7.The 10 −4 PAPR of an original OFDM frame was 12 dB.OPTS and the PSO-PTS improved on this by 7.6 dB with nearly the same performance for M = 8, while the performance loss with the iteration PTS is 8.4 dB.The iteration number of proposed technique is shown in Table 1.For M = 8, the OPTS technique requires 128 iterations per OFDM frame, while iteration PTS technique requires 16 iterations and the PSO-PTS technique without a threshold requires 88 iterations per OFDM frame.The complexity of iteration PTS is only 12.5% (16/128) of that of the PTS technique.The PSO-PTS technique with a threshold value is exhibited a lower complexity that only requires 23 iterations per OFDM frame.Thus, compared to the OPTS technique, the complexity of the PSO-PTS with threshold is only 18% (23/128 = 0.18).
Figure 8 illustrates some performance of the PTS technique in PAPR reduction for OFDM using PSO with acceleration factors c 1 and c 2 when N = 128, M = 4, and W = 2.It can be seen that when the acceleration factors increases resulted in the PAPR depression increasing.For example, at the level of CCDF being 0.1%, the acceleration factors c 1 = 0.5 and c 2 = 0.5, the PAPR is 8.3 dB, and acceleration factors c 1 = 2 and c 2 = 2, the PAPR is 6.8 dB.By these two examples of the acceleration factors c 1 and c 2 , the improvement in PAPR reduction is about 1.5 dB.Furthermore, we see that the PAPR  reduction of c 1 = 2 and c 2 = 2; and c 1 = 2.5 and c 2 = 2 have similar performance.Hence, after taking the effect of the reduction and the computational complexity into account, c 1 = 2 and c 2 = 2 is a suitable choice for our proposed PSO-based PTS technique.The values of c 1 and c 2 affect the behavior of the swarm in different ways: a bigger c 1 can increase the attraction of W P i for every particle and prevent the particle converging to W G quickly, while a bigger c 2 can decrease the attraction of W P i and prompt the swarm converging to the same W G .

CONCLUSION
In this paper, we analyze the PAPR reduction performance which is derived by using adjacent, interleaved, and random subblock partitioning methods.Random subblock partitioning method has derived the most effective performance, and interleaved subblock partition method has derived the worst.As the number of subblocks is increased, PAPR can be further reduced.Moreover, we formulate the phase weighting factors searching of PTS as a particular combination optimization problem and we apply the PSO technique to search the optimal combination of phase weighting factors for PTS to obtain almost the same PAPR reduction as that of optimal PTS while keeping low complexity.Simulations results show that PSO-based PTS method is an effective method to compromise a better tradeoff between PAPR reduction and computation complexity.By appropriate selection of phase weighting factors according to the required performance and tolerable complexity, the proposed partition scheme can be adaptive to QOS requirement.We illustrated that with this method we can develop algorithms which can achieve better performance-complexity tradeoff than the existing approaches.Additionally, the performance of the proposed method was slightly degraded compared to that of optimum method, PTS.However, the complexity of the proposed method was remarkably lower than that of optimum method.

Figure 1 :
Figure 1: The structure of transmitter with PSO-based PTS scheme.

Figure 3 :
Figure 3: CCDF of PSO technique for different G n when M = 16 and W = 4.

Figure 7 :
Figure 7: Comparison of the PSO-PTS technique with threshold PAPR, iterative PTS, PSO-PTS, and OPTS methods when M = 8.

Figure 8 :
Figure 8: CCDFs comparison of the PSO-based PTS scheme with different combinations of acceleration constants when N = 128, M = 4, and W = 2.

Table 1 :
The computational complexity of the OPTS, IPTS and PSO-PTS techniques with phase weighting factor W = 2.