About the Authors:
Soodabeh Darzi
* E-mail: [email protected]
Affiliation: Center of System and Machine Intelligence, College of Engineering, Universiti Tenaga Nasional, Selangor, Malaysia
Sieh Kiong Tiong
Affiliation: Power Engineering Center, Universiti Tenaga Nasional, Selangor, Malaysia
Mohammad Tariqul Islam
Affiliation: Department of Electrical, Electronic & Systems Engineering, Universiti Kebangsaan Malaysia, Selangor, Malaysia
Hassan Rezai Soleymanpour
Affiliation: Department of Electrical Engineering, Semnan University, Semnan, Iran
Salehin Kibria
Affiliation: Department of Electrical, Electronic & Systems Engineering, Universiti Kebangsaan Malaysia, Selangor, Malaysia
Abstract
An experience oriented-convergence improved gravitational search algorithm (ECGSA) based on two new modifications, searching through the best experiments and using of a dynamic gravitational damping coefficient (α), is introduced in this paper. ECGSA saves its best fitness function evaluations and uses those as the agents’ positions in searching process. In this way, the optimal found trajectories are retained and the search starts from these trajectories, which allow the algorithm to avoid the local optimums. Also, the agents can move faster in search space to obtain better exploration during the first stage of the searching process and they can converge rapidly to the optimal solution at the final stage of the search process by means of the proposed dynamic gravitational damping coefficient. The performance of ECGSA has been evaluated by applying it to eight standard benchmark functions along with six complicated composite test functions. It is also applied to adaptive beamforming problem as a practical issue to improve the weight vectors computed by minimum variance distortionless response (MVDR) beamforming technique. The results of implementation of the proposed algorithm are compared with some well-known heuristic methods and verified the proposed method in both reaching to optimal solutions and robustness.
Figures
Table 11
Fig 1
Fig 2
Fig 3
Table 1
Table 2
Table 3
Table 4
Fig 4
Fig 5
Table 5
Table 6
Table 7
Fig 6
Fig 7
Table 8
Table 9
Fig 8
Table 10
Table 11
Fig 1
Fig 2
Fig 3
Citation: Darzi S, Tiong SK, Tariqul Islam M, Rezai Soleymanpour H, Kibria S (2016) An Experience Oriented-Convergence Improved Gravitational Search Algorithm for Minimum Variance Distortionless Response Beamforming Optimum. PLoS ONE 11(7): e0156749. https://doi.org/10.1371/journal.pone.0156749
Editor: Wen-Bo Du, Beihang University, CHINA
Received: September 29, 2015; Accepted: May 18, 2016; Published: July 11, 2016
Copyright: © 2016 Darzi et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: This study is supported by Universiti Tenaga Nasional. The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
1.Introduction
Adaptive array techniques have been utilized to mitigate multipath fading and high co-channel interference in mobile, wireless, radar communications, satellite and many other similar applications to achieve robust performance and high data transmission rate [1]. One of the known array processing techniques is the minimum variance distortionless response (MVDR) technique, or Capon’s algorithm [2]. This approach is typically applied to cancel unwanted signals (interfering signals) and generate a strong beam towards the desired signal through its computed weight vectors. However, MVDR may have unacceptably low nulling in dynamic and various interference scenarios [3]. Recently, enormous amount of research has been dedicated to applying several optimization approaches for such problems. Thus, approaches such as genetic algorithms (GA) [4, 5], tabu search (TS) [6–8], particle swarm optimization (PSO) [9, 10], ant colony optimization (ACO) [11–13], differential evolution (DE) [14], clonal selection (CS) [15, 16], artificial bee colony (ABC) [17] and artificial immune system (AIS) [18,19] have been suggested and implemented to improve the robust performance, like interference cancelling, in antenna systems.
Recently, GSA was presented as a heuristic optimization algorithm inspired by Newtonian laws of gravitation. It was shown to outperform similar algorithms like PSO and GA for common benchmark functions [20]. Since the inception of GSA, a large number of researchers have introduced various modified and improved gravitational search algorithms based on the idea of memory and social information of PSO, novel strategies to define the agents’ search pattern along with other stochastic strategies [21–26]. The effectiveness of GSA and its version for binary encoded problems (BGSA) [27] in solving a set of nonlinear benchmark functions has been proven. However, very few of these GSA variants have been implemented for the beamforming applications [28].
In this study, we propose an experience oriented-convergence improved gravitational search algorithm (ECGSA) to deal with complicated optimization problems like the MVDR. ECGSA benefits from two modifications: one is to save the best fitness function evaluation of the agents during search process and to treat them as the agents’ effective positions in terms of applying force to other agents. By this modification, no discovered optimal trajectories are lost during the search process and also the agents can search through these optimal trajectories while avoiding local optimums. This is contrary to GSA, in which, the search trajectory is unstable as the agents move stochastically without any feedback to their best experiments. In ECGSA, if any agent finds better solution it updates its best position to its current location. The second modification is to use a dynamic gravitational damping coefficient parameter known as α. The proposed dynamic α is introduced to balance the exploration and exploitation properties of the standard GSA. With a relatively low value of α during first stage of search process, the agents can have larger velocities for a better exploration. On the other hand, the agents can converge to the optimal solution effectively at the final stage of the search process by rapidly increasing the value of α, consequently, reducing the agents velocities.
The most significant difference between stochastic leader gravitational search algorithm (SLGSA) [28] and ECGSA is that ECGSA retains memory of the results from previous iterations, whereas, SLGSA is memory less (i.e. uses results from current iteration only). This results in loss of optimal trajectory in SLGSA. Furthermore, SLGSA uses a stochastic method to determine which agents are allowed to influence other agents. This increases the exploration process by increasing diversity, hence, preventing premature convergence. ECGSA uses the more conventional approach to maintain diversity in the early phase by using high value of gravitational coefficient function, G(t). Hence, ECGSA is fundamentally less stochastic than SLGSA. Thus ECGSA converges with greater precision than SLGSA. As SLGSA has the stochastic method of increasing diversity, it does not use the activation ratio, σ parameter. Overall, the two algorithms provide fundamentally different approaches to balancing the exploration and exploitation aspects of the search process.
In this study, the ECGSA based MVDR has been used in beamforming with uniform linear antenna arrays of 0.5λ spacing between adjacent elements and radiating at 2.3 GHz. The proposed ECGSA is also verified by testing it on some well-known complicated optimization problems. The rest of the paper is organized in the following order: Section 2 presents a brief review of GSA. Our proposed ECGSA is introduced in section 3. The validation of the proposed ECGSA via various functions and the simulation results are reported in section 4. Section 5 presents the basics of the conventional MVDR. Section 6 demonstrates the integration of MVDR with ECGSA along with performance comparison of conventional MVDR with MVDR-PSO, MVDR-GSA and MVDR-ECGSA in different interference scenarios. Finally, Section 7 presents concluding remarks about this investigation.
2.Gravitational Search Algorithm
The GSA algorithm was proposed by Rashedi, et al. [20, 28] as a global optimization method in 2009. The GSA agent’s movements are estimated through their masses. All the GSA search agents gravitate toward the agents with heavier masses due to their gravitational attraction as shown in Fig 1. Agents (masses) with better fitness value in GSA will correspond to larger masses.
[Figure omitted. See PDF.]
Fig 1. Mass acceleration toward the result force in GSA [20].
https://doi.org/10.1371/journal.pone.0156749.g001
For the ith agent of N agents of GSA, the position vector is defined by:(1)where xid represents the position of the ith agent in the dth dimension (decision variable), and n represents the total number of decision variables. Inertia masses are calculated by the fitness evaluation as follows:(2)(3)where fiti(t) represents the fitness value of the ith agent at iteration t; worst(t) and best(t) represent the worst and best fitness of all agents, respectively; and Mi(t) represents the mass value of the ith agent at iteration t. The force operating on mass i from mass j can be calculated as follows:(4)where Mi and Mj represent masses of the ith and jth agents respectively; ε is a significantly small constant that prevents division by zero; Ri,j(t) represents the Euclidean distance between the ith and jth agents defined as follows:(5)and G(t) is the gravitational coefficient at iteration t. The gravitational coefficient G(t) is a monotonically decreasing function of time when it is set to G0 at the starting and will be exponentially damped to examine the accuracy of the algorithm according to the equation below:(6)where α is gravitational damping coefficient which is constant in GSA, t is the current iteration and tmax is the maximum iteration number. The total force affecting the ith agent in the dth dimension is calculated as follows:(7)where kbest is a time dependent function that is able to qualify the exploitation and exploration. It is initialized to the total number of agents (k0) at the start and linearly decreased to one with time [20]. randj is a random number between zero and unity. By the laws of motion, , the acceleration of agent i in dimension d, could be obtained as follows:(8)
The position of agent i depends on the agent velocity and it can be updated using the following equations:(9)(10)where vid(t) represents the velocity of agent i in dimension d at iteration t; and randi is a random number between zero and unity.
3.Experience Oriented-Convergence Improved Gravitational Search Algorithm
In this study, two distinct modifications, using the best found experiments as the agents’ present positions and introducing a dynamic gravitational damping coefficient value, are introduced to GSA as follows:
A. Using the best found experiments of agents as their present positions
The standard GSA does not retain the best fitness function evaluations that the agents find during search process. This is in contradiction of other stochastic search algorithms such as PSO [29], GA [30], DE [31] and simulated annealing (SA) [32]. GSA only relies on the current information of the agents to determine the present search trajectory. This is a weakness of GSA, because the optimal found trajectories may be lost due to updating them by non-optimal new ones. In this paper, we overcome this defect by using the best fitness function evaluations of the agents as their present positions for calculating new movements. To do this, the position updating (10) is modified as follows:(11)
Where represents the unconfirmed position of agent i in dimension d which is calculated from the following equation:(12)
And is the unconfirmed position vector of agent i containing .
Eq (11) simply states that the new movement of an agent is accepted to be its new position, if the movement results in a better solution, otherwise the position of the agent remains unchanged at the present iteration. By this modification no optimal found trajectory is lost, and also the movement of any agent is started from its best experiment at any iteration, which achieves good local search around the optimal found solutions.
B. Dynamic gravitational damping coefficient value
Ideally, a stochastic search based algorithm should focus on exploration and follow it with exploitation. The exploration phase allows the algorithm to gather information about the overall search space. Once the search space is sufficiently explored, the focus can be shifted to exploiting the information gathered by exploration. In GSA, the gravitational coefficient G(t) balances the exploration and exploitation capabilities through increasing/decreasing the magnitude of agent acceleration vector. Higher values of G(t) are needed for relatively large movements at initial stages of search process, while lower values are required for a better convergence to optimal solution at the final phase. The value of G(t) is exponentially decreased during iterations according to (6). The gravitational damping coefficient α adjusts the degree of damping. The more gravitational damping coefficient α is, faster the gravitational coefficient G(t) drops. But finding a value of α to be optimal for both exploration and exploitation is a hard decision. If α is too high, the algorithm suffers from premature convergence. Similarly, if it is too low, the algorithm couldn’t converge to the optimal region at the final iterations. In ECGSA, we propose a gravitational damping coefficient α according to the following equation:(13)where σ is a constant value between 0 and 1, representing the ratio of the first sub function being activated. The variations of the proposed α(t) along with G(t) are shown on Fig 2.
[Figure omitted. See PDF.]
Fig 2. The variations of the proposed α(t) along with G(t).
https://doi.org/10.1371/journal.pone.0156749.g002
According to Fig 2, the value of α is constant and relatively low for a large portion of iterations resulting in good exploration capability. It is linearly increased to a high value at the final iterations resulting in good exploitation and convergence behavior.
According to Fig 3, selecting the value of σ from 0.2 to 0.8 yields better results, although, the algorithm is less sensitive with respect to this parameter.
[Figure omitted. See PDF.]
Fig 3. Performance comparison of different values of σ for multimodal function.
https://doi.org/10.1371/journal.pone.0156749.g003
Effect of the parameter σ on optimization performance was studied on multimodal function F8 and the result is shown in Fig 3. In this function, the special case of ECGSA with only exploration ability (σ = 1) cannot achieve the best solution. The lower values of σ produce significantly better performance than σ = 1 consistently due to exploitative nature of ECGSA with σ≤0.8. In the multimodal case, the rapid convergence of σ = 0 is premature and the algorithm stagnates at local minima. ECGSA with σ = 1 lacks the exploitative ability almost completely, thus it is out performed by ECGSA with 0.2≤σ≤0.8. The interval 0.2≤σ≤0.8 produces a balanced exploration and exploitation ability as is obvious in F8. This case shows that there is negligible variation in performance in that interval, which indicates that the choice of σ is not critical. Only extreme values of σ should be avoided in real world optimization problems in order to maintain robustness of ECGSA.
The proposed dynamic gravitational damping coefficient somehow is similar to the temperature value of the SA algorithm [32] which controls the exploration and exploitation search abilities trade-off of the algorithm, but they are implemented via different equations.
4.Verification of the Proposed ECGSA
For verification, ECGSA is tested on eight representative standard benchmark functions along with six composite test functions. All the test functions are to be minimized and the relevant information can be found in [20] and [33] for the standard benchmark functions and the composite test functions, respectively.
For the proposed ECGSA, N is set to 50, tmax is set to1000, G0 is set to 100, αmin and αmax are set to 10 and 25, respectively, and σ is set to 0.6. The other parameters of ECGSA are set the same as in the standard GSA [20].
The standard benchmark functions are summarized in Tables 1 and 2 for unimodal and multimodal functions, respectively. The results of applying the algorithm to the unimodal and multimodal benchmark functions over 30 runs are reported in Tables 3 and 4, respectively. According to the tables, the consequences are reported for mean, median and best values of the obtained solutions, and the standard deviation value of them. Since GSA was shown to outperform similar algorithms like PSO and GA for the benchmark functions [20], the results are compared with GSA and some of its variants.
[Figure omitted. See PDF.]
Table 1. Unimodal benchmark functions.
https://doi.org/10.1371/journal.pone.0156749.t001
[Figure omitted. See PDF.]
Table 2. Multimodal benchmark functions.
https://doi.org/10.1371/journal.pone.0156749.t002
[Figure omitted. See PDF.]
Table 3. Minimization result of the unimodal benchmark functions in Table 1.
https://doi.org/10.1371/journal.pone.0156749.t003
[Figure omitted. See PDF.]
Table 4. Minimization result of multimodal benchmark functions in Table 2.
https://doi.org/10.1371/journal.pone.0156749.t004
As it can be seen, ECGSA outperforms GSA and its variants for all the benchmark functions tested in this study. The standard deviation of the optimization result from different runs also verifies the robustness of the proposed algorithm.
According to Table 3, the proposed ECGSA can optimize unimodal functions more consistently than the GSA variants. This shows that the exploitation search ability of ECGSA is significantly better than GSA and its variants as the unimodal functions require strong exploitation ability for optimization because of the convexity of such problems.
The convergence characteristic of the proposed algorithm with respect to the standard GSA is shown on Fig 4 for validating the results related to the unimodal functions. The premature convergence of GSA for F3 occurs at about 200 to 300 iterations. But due to setting a low value for α, ECGSA can resist premature convergence. Considering the best fitness function evaluations for improving the optimal search trajectory results in good exploration. This strategy would result in poor convergence speed if αmin is used throughout the optimization process. At 600 iterations, ECGSA optimization of F3 shows a sharp decrease due the start of the exploitation phase (α increasing linearly). The results suggest that a good trade-off has been established in the exploitation and exploration aspects of ECGSA for unimodal functions.
[Figure omitted. See PDF.]
Fig 4. Performance comparison of ECGSA and GSA for F3.
https://doi.org/10.1371/journal.pone.0156749.g004
Table 4 presents the results of optimization of multimodal functions in Table 2 using ECGSA, GSA and its variants. The performance of ECGSA is significantly better than GSA and its variants as illustrated by the mean optimization results in Table 4. Furthermore, ECGSA also has lower standard deviation showing that it is more stable than GSA and its variants. This establishes the superior ability to escape local minima and exploration performance of ECGSA.
Fig 5 shows similar traits to Fig 4. The exploration phase of ECGSA is benefitted by lower value of αmin compared to α of GSA. Thus GSA prematurely converges to suboptimal local minima before 200 iterations. During the exploitation phase, once α begins to increase, the sharp convergence leads to a rapid drop after 600 iterations in Fig 5. Using the best found experiments of agents as their present positions results in the ECGSA curve falling away from the GSA curve in both Figs 4 and 5, even before the α starts to vary dynamically. This clearly indicates that using the best found fitness function evaluations of agents as their present positions results in improved exploration in the early phase. As α starts to vary at 600th iteration, the best fitness function evaluations of agents improved exploration allows greater exploitation by attracting the agents towards their best locations, rather than the current, possibly suboptimal, location. Thus, the convergence process is more stable, allowing the proposed algorithm to optimize multiple complex composite test functions more accurately and efficiently than most well-known heuristic algorithms.
[Figure omitted. See PDF.]
Fig 5. Performance comparison of ECGSA and GSA for F8.
https://doi.org/10.1371/journal.pone.0156749.g005
The Wilcoxon’s Rank Sum test [37] is also used to analyze and to validate the proposed algorithm. Since the underlying probability density function for the algorithm solutions is unknown, techniques like T-test are unusable. Table 5 illustrates the Wilcoxon's rank-sum test between GSA and ECGSA for the benchmark functions. Sufficiently large sample sizes are considered to ensure a viable comparison.
[Figure omitted. See PDF.]
Table 5. Wilcoxon's rank-sum test between GSA and ECGSA for benchmark functions.
https://doi.org/10.1371/journal.pone.0156749.t005
The table shows extremely low p-values for the results indicating that ECGSA can be considered superior to GSA for both unimodal and multimodal benchmark functions within a confidence interval of 99.9% (p-value<0.001).
For more verification, we have tested ECGSA on the composite benchmark functions introduced recently to validate the performance of heuristic algorithms [33]. These functions combine multiple features of various unimodal and multimodal functions. These functions test several properties of the proposed algorithm simultaneously. The results of the algorithm for the composite function optimizations are from 20 runs, for the sake of comparison, as suggested in literature [33].
Table 6 compares the performance of the proposed algorithm versus several well-known heuristic methods. The results clearly indicate that the proposed ECGSA achieves significantly better solutions than all other algorithms considered in Table 6. Only DE marginally outperforms GSA in CF6. However, ECGSA comprehensively outperforms DE for the other 5 composite test functions. ECGSA shows high stability as it exhibits the lowest standard deviation for most of the composite test functions.
[Figure omitted. See PDF.]
Table 6. Composition tests functions.
https://doi.org/10.1371/journal.pone.0156749.t006
Table 7 confirms the superiority of proposed ECGSA over GSA using the Wilcoxon Rank Sum test. The p-values indicate that ECGSA achieves superior results than GSA within a confidence interval of 99.9% for all of the composite functions except CF3.
[Figure omitted. See PDF.]
Table 7. Wilcoxon's rank-sum test between GSA and ECGSA for composition functions.
https://doi.org/10.1371/journal.pone.0156749.t007
5.MVDR Beamforming
MVDR is a minimum output energy (MOE) system based beamforming technique, as explicitly described previously [28]. MVDR keeps a distortionless main lobe response towards the wanted signal while simultaneously minimizing the array output power. The weights of MVDR can be written as [18]:(14)where R is the covariance matrix, H is the Hermitation transpose and a(θ) is the steering vector of the desired signal. Where the steering vector a(θ) is given by:(15)where d is the space between the elements of the antenna, θi is the desired angle, and m is the number of elements. However, typically, low nulling levels towards multiple interference sources is one of the major drawbacks of MVDR beamforming technique [3, 18, 38]. Therefore, many studies have tried to achieve nulls in the interference directions with various solutions [18, 38–39].
6.Model Application
The robustness and efficiency of the proposed method have been evaluated through eight benchmark functions and six composite test functions in the previous section. In this section, the calculation of signal to interference noise ratio (SINR) and applicability of the proposed method to increase SINR of MVDR technique is presented and it is compared with conventional MVDR, MVDR-PSO and MVDR-GSA. Two cases of interference scenarios are considered in this study. Both cases have one user at 0°. The first case contains a single interference source positioned at 30° and the second case considers two interference sources at 30° and 50°. The employment of the ECGSA for MVDR analysis has been carried out using MATLAB®.
6.1 Signal to interference and noise ratio calculation
In this paper, the ECGSA was employed to enhance the MVDR beamforming technique performance by increasing the SINR value of four element array. The MVDR integrated with ECGSA, GSA and PSO method will optimize the SINR via the complex weights, as described in detail previously [28].
In these algorithms, the wmvdr (MVDR weight vector) will be retained for one member in the initial population, whereas, the rest of the population will be initialized randomly as shown in (17). This system will initiate by generating the N agents, which is indicated by WN weight vectors. The weight vectors in every agent contain M number of weight vectors, where M is the number of elements in the array.(16)where population size, i, in ECGSA, GSA and PSO is set to N, and dimension d is equal to number of sensor M as shown in (16).
The weight vectors of the entire population can be illustrated in the matrix format below:(17)
Where
The fitness function, according to literature [28], is the SINR, as shown in (18). Thus, the optimization processes will try to maximize the fitness function and find the corresponding optimal weight vectors as shown in Fig 6.
(18)
Where
[Figure omitted. See PDF.]
Fig 6. The simplified flowchart of ECGSA beamforming.
https://doi.org/10.1371/journal.pone.0156749.g006
6.2 Simulation
In this section, multiple interference scenarios were simulated to validate the proposed approach for real world applications. The parameters of GSA are chosen according to the guidelines and recommendations presented in [20]. These configurations of GSA have also been utilized extensively after the development of GSA [34–36]. In all of the PSO-based simulations, the inertia weight w = 0.7298, and the acceleration constants c1 = c2 = 1.49618 were utilized. These values are chosen based on the common settings in the literature [40–42]. The algorithms are simulated 20 times with maximum number of iteration set to 100, and the best results are recorded.
6.2.1 Case 1: One user one interferences
One interference source at 30° and user at 0° has been assumed in the first case study. Table 8 illustrates the weights corresponding to optimal solutions determined by conventional MVDR, PSO-MVDR, GSA-MVDR, SLGSA-MVDR [28] and ECGSA-MVDR. These weights yield the power response as shown in Fig 7.
[Figure omitted. See PDF.]
Fig 7. Comparison of performance of power response with 100 iterations for user at 0° with interference at 30°.
(a) MVDR, (b) PSO-MVDR, (c) GSA-MVDR, (d) SLGSA-MVDR [28] and (e) ECGSA-MVDR.
https://doi.org/10.1371/journal.pone.0156749.g007
[Figure omitted. See PDF.]
Table 8. Comparison of weight vectors for conventional MVDR, PSO-MVDR, GSA-MVDR, SLGSA-MVDR [28] and ECGSA-MVDR for user at 0° and interferences at 30°.
https://doi.org/10.1371/journal.pone.0156749.t008
Table 9 illustrates improvement of SINR by PSO and GSA, compared to conventional MVDR, are 38.89% and 65.04%, respectively. The SLGSA [28] and ECGSA show the best performance of 72.17% percent improvement. The SINR of 69.99 is the global maximum for this case, which is achieved by both SLGSA-MVDR and ECGSA-MVDR as only one deep null is required. Failure of MVDR, GSA-MVDR and PSO-MVDR to achieve sufficiently deep nulls clearly illustrates their limitations.
[Figure omitted. See PDF.]
Table 9. Comparison of SINR calculation for conventional MVDR, PSO-MVDR, GSA-MVDR, SLGSA-MVDR [28] and ECGSA-MVDR for user at 0° and interference at 30°.
https://doi.org/10.1371/journal.pone.0156749.t009
6.2.2 Case 2: One user two interferences
Two interferences at 30°, 50° and user at 0°has been assumed as the second case study. The weights corresponding to optimal solutions determined by the different beamforming algorithms are shown in Table 10. The resulting power responses from these weights are illustrates in Fig 8.
[Figure omitted. See PDF.]
Fig 8. Comparison of performance of power response with 100 iterations for user at 0°with two interferences at 30° and 50°.
(a) MVDR, (b) PSO-MVDR, (c) GSA-MVDR, (d) SLGSA-MVDR [28] and (e) ECGAS-MVDR.
https://doi.org/10.1371/journal.pone.0156749.g008
[Figure omitted. See PDF.]
Table 10. Comparison of weight vectors for conventional MVDR, PSO-MVDR, GSA-MVDR, SLGSA-MVDR [28] and ECGSA-MVDR for user at 0° and interferences at 30°and 50°.
https://doi.org/10.1371/journal.pone.0156749.t010
Table 11 shows that the improvement of SINR by PSO, GSA and SLGSA [28], compared to conventional MVDR method, are 67.75%, 87.84% and 105.84%, respectively. It can be seen that ECGSA shows the best performance, which is 106.58% percent of improvement. The increase in number of interference sources inherently increases the difficulty of the optimization problem. The superiority of ECGSA-MVDR in case 2 is more apparent than in case 1, especially for SLGSA-MVDR [28]. Thus, ECGSA is more versatile and robust than GSA, PSO and SLGSA.
[Figure omitted. See PDF.]
Table 11. Comparison of SINR calculation for conventional MVDR, PSO-MVDR, GSA-MVDR, SLGSA-MVDR [28] and ECGSA-MVDR for user at 0° and interference at 30°and 50°.
https://doi.org/10.1371/journal.pone.0156749.t011
Overall, the performance disparity between proposed algorithm and SLGSA is expected to increase further with higher number of interference sources. In case 1, as SLGSA-MVDR is able to locate the optimal solution for one interference source, there is no potential for further improvement. As a result, the proposed ECGSA algorithm provides the same result. Subsequently, case 2 was investigated to reveal the minor, yet tangible, SINR improvement by the proposed algorithm over SLGSA-MVDR. The superior performance of ECGSA in case 2, along with the results of the benchmark problems, establishes ECGSA as a suitable interference mitigation algorithm.
7.Conclusions
An experience oriented-convergence improved gravitational search algorithm with dynamic gravitational damping coefficient, α, and searching through the best experiments is presented in this paper. GSA outperforms most known heuristic algorithms but the constant α and unstable search trajectory without any feedback to their best experiments decrease the effectiveness of GSA. Consequently, the proposed dynamic gravitational damping coefficient alters the gravitational constant, resulting in better balance between exploration and exploitation. Using the best found experiments of agents as their present positions provides steadier trajectory, as GSA uses the stochastically varying current position. This modified system exhibits better convergence rate, precision, stability and robustness compared to GSA on fourteen unimodal, multimodal and composite test functions. The proposed algorithm has been applied on MVDR beamforming problems in order to verify the effectiveness of ECGSA in real world applications. Results show superior performance of ECGSA-MVDR in comparison to the conventional MVDR, MVDR-PSO, MVDR-GSA and MVDR-SLGSA due to its ability in achieving higher SINR.
Supporting Information
[Figure omitted. See PDF.]
S1 File. Data for ECGSA result in Tables 3 and 4.
https://doi.org/10.1371/journal.pone.0156749.s001
(XLSX)
S2 File. Data for ECGSA and GSA results in Table 6.
https://doi.org/10.1371/journal.pone.0156749.s002
(XLSX)
Author Contributions
Conceived and designed the experiments: SD SKT. Performed the experiments: MTI. Analyzed the data: SD SK. Contributed reagents/materials/analysis tools: SD HRS. Wrote the paper: SD.
Citation: Darzi S, Tiong SK, Tariqul Islam M, Rezai Soleymanpour H, Kibria S (2016) An Experience Oriented-Convergence Improved Gravitational Search Algorithm for Minimum Variance Distortionless Response Beamforming Optimum. PLoS ONE 11(7): e0156749. https://doi.org/10.1371/journal.pone.0156749
1. Barrettand M., Arnott R., 1994. Adaptive antennas for mobile communications, Electronics and Communication Engineering Journal, 6(4), 203–214.
2. Capon J., 1969. High-resolution frequency-wavenumber spectrum analysis. Institute of Electrical and Electronics Engineers (IEEE), 57 (8), 1408–1418.
3. Dahrouj H., Yu W., 2010. Coordinated beamforming for the multicell multi-antenna wireless system, IEEE transactions on wireless communications, 9(5), 1748–1759.
4. Kiong T. S., Ismail M., Hassan A., 2006. WCDMA forward link capacity improvement by using adaptive antenna with genetic algorithm assisted MDPC beamforming technique, Journal of Applied Sciences, 6(8), 1766–1773.
5. Xu Z., Li H., Liu Q. Z., Li J. Y., 2008. Pattern synthesis of conformal antenna array by the hybrid genetic algorithm, Progress In Electromagnetics Research, 79, 75–90.
6. Kaid Omar, O., Debbat F., Boudghene Stambouli, A., 2012. Null steering beamforming using hybrid algorithm based on honey bees mating optimization and tabu search in adaptive antenna array, Electromagnetic Research C, 32, 65–80.
7. Guney K., Akdagli A., 2001. Null steering of linear antenna arrays using modified tabu search algorithm, Progress In Electromagnetics Research-pier, 33, 167–182.
8. Akdagli A., Guney K., 2003. Shaped-beam pattern synthesis of equally and unequally spaced linear antenna arrays using a modified tabu search algorithm, Microwave and Optical Technology Letters, 36, 16–20.
9. Zaharis Z. D., Yioultsis T. V., 2011. A novel adaptive beamforming technique applied on linear antenna arrays using adaptive mutated Boolean PSO, Electromagnetic Research, 117, 165–179.
10. Khodier M. M., Christodoulou C. G., 2005. Linear array geometry synthesis with minimum side lobe level and null control using particle swarm optimization, IEEE Transactions On Antennas And Propagation, 53(8), 2674–2679.
11. Karaboga N., Guney K., Akdagli A., 2002. Null steering of linear antenna arrays by using modified touring ant colony optimization algorithm, International Journal of RF and Microwave Computer Aided Engineering, 12, 375–383.
12. Karaboga D., Guney K., Akdagli A., 2004. Antenna array pattern nulling by controlling both the amplitude and the phase using modified touring ant colony optimization algorithm, International Journal of Electronics, 91, 241–251.
13. Akdagli A., Guney K., Karaboga D., 2006. Touring ant colony optimization algorithm for shaped-beam pattern synthesis of linear antenna arrays, Electromagnetics Research, 26, 615–628.
14. Yang S. W., Gan Y. B., Qing A. Y., 2004. Antenna-array pattern nulling using a differential evolution algorithm”, International Journal of RF and Microwave Computer Aided Engineering, 14, 57–63.
15. Babayigit B., Akdagli A., Guney K., 2006. A clonal selection algorithm for null synthesizing of linear antenna arrays by amplitude control, Journal of Electromagnetic Waves and Applications, 20, 1007–1020.
16. Guney K., Akdagli A., Babayigit B., 2006. Shaped-beam pattern synthesis of linear antenna arrays with the use of a clonal selection algorithm, Neural Network World, 16, 489–501.
17. Guney K., Onay M., 2007. Amplitude-only pattern nulling of linear antenna arrays with the use of bees algorithm, Progress In Electromagnetics Research, 70, 21–36.
18. Tiong, S. K., Salem, B., Koh, S. P., Sankar, K. P., 2014. Minimum variance distortionless response (MVDR) beamformer with enhanced nulling level control via dynamic mutated artificial immune system, The Scientific World Journal, Article ID 164053.
19. Darzi, S., Tiong, S. K., Islam, M. T., Ismail, M., Kibria, S., Salem, B., 2014. Null steering of adaptive beamforming using linear constraint minimum variance assisted by particle swarm optimization, dynamic mutated artificial immune system, and gravitational search algorithm, The Scientific World Journal, Article ID 724639.
20. Rashedi E., Nezamabadi-pour H., Saryazdi S., 2009. GSA: Gravitational Search Algorithm, Information Sciences, 179, 2232–2248.
21. Mirjalili S., Hashim S. Z. M., Sardroudi H. M., 2012. Training feed forward neural networks using hybrid particle swarm optimization and gravitational search algorithm, Applied Mathematics and Computation, 218(22), 11125–11137.
22. Ashouri M., Hosseini S. M., 2013. Application of new hybrid particle swarm optimization and gravitational search algorithm for non-convex economic load dispatch problem, Journal of Advances in Computer Research, 4(2), 41–51.
23. Li Ch., Zhou J., 2011. Parameters identification of hydraulic turbine governing system using improved gravitational search algorithm, Energy Conversion and Management, 52, 374–381.
24. Khatibiniaa M., Khosravi Sh., 2014. A hybrid approach based on an improved gravitational search algorithm and orthogonal crossover for optimal shape design of concrete gravity dams, Applied Soft Computing, 16, 223–233.
25. Li Ch., Zhou J., Xiao J., Xiao H., 2012, Parameters identification of chaotic system by chaotic gravitational search algorithm, Chaos, Solitons & Fractals, 45, 539–547.
26. Han X. H., Chang X. M., 2012. Chaotic digital secure communication based on a modified gravitational search algorithm filter”, Information Sciences, 208, 14–27.
27. Rashedi E., Nezamabadi-pour H., Saryazdi S., 2010. BGSA: binary gravitational search algorithm, Natural Computing, 9(3), 727–745.
28. Darzi S., Islam M. T., Tiong S. K., Kibria S., Singh M., 2015. Stochastic leader gravitational search algorithm for enhanced adaptive beamforming technique, Plos One |
29. Kennedy, J., Eberhart, R., 1995. Particle swarm optimization, in Proceedings of the IEEE International Conference on Neural Networks, 1942–1948.
30. Holland J. H., 1973. Genetic algorithms and the optimal allocation of trials, SIAM Journal on Computing, 2(2), 88–105.
31. Storn R., Price K., 1997. Differential evolution—a simple and efficient heuristic for global optimization over continuous spaces, Journal of Global Optimization, 11, 341–359.
32. Nouraniy Y., Andresenz B., 1998. A comparison of simulated annealing cooling strategies, J. Phys. A: Math. Gen, 31, 8373–8385.
33. Liang, J., Suganthan, P., Deb, K., 2005. Novel composition test functions for numerical global optimization, in Swarm Intelligence Symposium, Proceedings 2005 IEEE, 68–75.
34. Binjie G., Feng P., 2013,Modified gravitational search algorithm with particle memory ability and its application, International Journal of Innovative Computing, Information and Control, 9(11), 4531–4544.
35. Khajehzadeh M., Taha M.R., Shafie A. E., Eslami M., 2012. A modified gravitational search algorithm for slope stability analysis, Engineering Applications of Artificial Intelligence, 25, 1589–1597.
36. Khajehzadeh M., Taha M. R., Eslami M., 2013. Efficient gravitational search algorithm for optimum design of retaining walls, Structural Engineering and Mechanics, 45(1), 111–127.
37. Derrac J., García S., Molina D., Herrera F., 2011. A practical tutorial on the use of nonparametric statistical tests as a methodology for comparing evolutionary and swarm intelligence algorithms, Swarm and Evolutionary Computation, 1(1), 3–18.
38. Liu F. L., Wang J. K., 2011. Robust mvdr beamformer for nulling level control via multi-parametric quadratic programming, Progress In Electromagnetic s Research C, 20, 239–254.
39. Lorenz R. G., Boyd S. P., 2005. Robust minimum variance beamforming, IEEE Trans. Signal Process, 53(5), 1684–1696.
40. Van Den Bergh, F., 2001. An analysis of particle swarm optimizers, Ph.D. dissertation, Faculty Nat. Agricult. Sci., Univ. Pretoria, Pretoria, South Africa.
41. Shi Y., Eberhart, R., 1998. A modified particle swarm optimizer, IEEE Int. CEC, 69–73.
42. Xue B., Zhang M., Browne W. N., 2013. Particle swarm optimization for feature selection in classification: a multi-objective approach, IEEE Transactions on Cybernetics, 43, (6).
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2016 Darzi et al. This is an open access article distributed under the terms of the Creative Commons Attribution License: http://creativecommons.org/licenses/by/4.0/ (the “License”), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
An experience oriented-convergence improved gravitational search algorithm (ECGSA) based on two new modifications, searching through the best experiments and using of a dynamic gravitational damping coefficient (α), is introduced in this paper. ECGSA saves its best fitness function evaluations and uses those as the agents’ positions in searching process. In this way, the optimal found trajectories are retained and the search starts from these trajectories, which allow the algorithm to avoid the local optimums. Also, the agents can move faster in search space to obtain better exploration during the first stage of the searching process and they can converge rapidly to the optimal solution at the final stage of the search process by means of the proposed dynamic gravitational damping coefficient. The performance of ECGSA has been evaluated by applying it to eight standard benchmark functions along with six complicated composite test functions. It is also applied to adaptive beamforming problem as a practical issue to improve the weight vectors computed by minimum variance distortionless response (MVDR) beamforming technique. The results of implementation of the proposed algorithm are compared with some well-known heuristic methods and verified the proposed method in both reaching to optimal solutions and robustness.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer