Novel Method for Optimal Synthesis of 5 G Millimeter Wave Linear Antenna Array

We will propose a useful method for 5G mm wave antenna array synthesis, based on Genetic Algorithm for the synthesis of linear arraywith nonuniform interelement spacing. Our designmethodwas used to obtain the optimal position of the elements in order to get theminimum side lobe level and nulls in desired directions.The simulation results verify that proposedmethod outperforms the previously published methods in terms of suppression side lobe level while maintaining nulls in specified directions.The flexibility of proposed algorithm shows good potential for the antenna array synthesis.


Introduction
In order to satisfy necessary 5G quality-of-standard criteria and meet user mobility, due to the higher path loss at mm wave frequency range, multiple antenna arrays are typically used in outdoor mm wave systems for providing additional gain [1,2].They play an important role in detecting and processing signals arriving from different directions.The performance of the systems depends on the efficient design of the antenna arrays.The desired radiation pattern can be realized by choosing geometrical shape of antenna array, element excitation amplitude, element excitation phase, and element spacing.The objectives in design of the antenna arrays are archiving minimum side lobe level (SLL) and narrow first null beam width (FNBW).Interference issues will become of crucial importance for coexistence of 5G devices, since number of mm wave devices will grow extensively in near future [3].The increasing growth of 5G devices number will prompt the study of array pattern nulling techniques.This problem will reach its culmination in future dense mm wave application, such as for the 60 GHz dense indoor communication, and various outdoor mm wave 5G communication in densely populated urban areas [4].So the goal of antenna array synthesis is to achieve the minimum SLL with narrow FNBW and obtaining narrow or broad nulls in directions of interfering signals.The broad nulls are needed when the direction of arrival of interference may vary slightly with time.
For the linear array geometry, by optimizing the spacing between the elements while keeping the uniform excitation we can suppress side lobe level while preserving the gain of the main beam and can control nulling.
Methods used for the antenna array synthesis can be classified in two categories: deterministic and stochastic.There are several deterministic methods: Schelkunoff Polynomial Method [5], Dolph-Chebyshev Method [6,7], Woodward-Lawson Method [8,9], and Fourier Transform Method [6,8,9].In modern research the most common methods are stochastic, because stochastic methods can almost all be proven to find a global minimum with asymptotic convergence in probability.The biggest advantages of using stochastic methods are their ability in dealing with large number of optimization parameters and avoiding getting stuck in local minima.Some of the methods used for the antenna array synthesis are Evolutionary Algorithms (EA) [10], Genetic Algorithms (GA) [11][12][13], Tabu Search (TS) [11,14], Particle Swarm Optimization (PSO) [15][16][17], Ant Colony Optimization (ACO) [18], Nature-Inspired Cuckoo Search (CS) [19,20], Combination of Global and Local Search [21], and Cat Swarm Optimization [22].All the above-mentioned methods have shown the capability of searching for global solution in electromagnetic optimization problems.
In this paper, we proposed a useful method based on combination of two global searches.The main idea comes from method based on combination of global and local search [21].Changing local search into global search we avoid getting stuck in local minimum and this leads to better exploration of feasible region.Better exploration leads to better solution of optimization problem.Presented algorithm was used to obtain the optimal position of the elements in order to get the minimum side lobe level and nulls in desired directions.We applied for optimization the same radiation pattern requirements of the selected papers.The simulation results reveal that design of antenna arrays using the presented method provides considerable enhancements compared with the synthesis obtained from other published methods.In this approach we will concern arrays of isotropic point sources; this is of great value because the pattern of any antenna can be regarded as being produced by an array of point sources.
Description of antenna array design problem is given in Section 2. In Section 3, a description of the used algorithm is presented.Results are presented in Section 4. Finally, the conclusion of this work is presented in Section 5.

Problem Formulation
The antenna array radiation pattern (RP) may be found according to the pattern multiplication theorem (see, e.g., Ch. 6 from [6,24] and Ch. 3 from [25]): EL (element pattern) is the pattern of the individual array element.
AF (array factor) is a function dependent on the physical placement of antenna elements, amplitude, and phase of excitation.
If we replace each element of the antenna array with an isotopic point source the resulting pattern is the array factor, because the element pattern of isotopic point source is equal to 1.The geometry of the linear antenna array with 2 isotropic elements placed symmetrically along the -axis is given in Figure 1.
In this paper linear antenna array has identical elements with nonuniform interelement spacing and uniform excitation.Array factor in the azimuth plane can be expressed as follows (for more see [6,8,9,[24][25][26]): where  = 2/ is wave number,  is the wavelength,   is excitation amplitude of th element,   is excitation phase of th element,   is location of th element, and  is angle between the line of observer and the source position.
If we assume a uniform excitation of phase as   = 0 for each antenna array element, the array factor can be written in a simple form as follows: In the field of electromagnetism the decibel (dB) is often used as a unit of measurement for array factor.Normalized array factor and normalized array factor in dB are given by

Algorithm
The presented algorithm is applied to achieve the desired radiation pattern with minimum SLL and narrow or broad nulls in specified directions.Our approach uses a combination of two global searches: global search with respect to  and global search with respect to .The global searches were performed using Genetic Algorithm with existing function "ga" from MATLAB software package (R2016a).The global search with respect to  was performed on region [ 11 ,  12 ] with fixed permanent solution   for the cost function (4).With this global search we locate maximum side lobe level.
The cost function for step global search for  following expression was used: First term in ( 6) is employed to minimize the side lobe level between the desired angles  11 and  12 , whereas the second one is for achieving deep nulls in desired directions   .For reducing mutual coupling effects between the elements of the antenna array the following conditions must be satisfied: The global search with respect to  on a small neighbourhood of current solution was performed with fixed   .Lower and upper bounds for each coordinate of vector   are defined using the following expression: Pseudocode of proposed algorithm is given in Algorithm 1.

Simulation Results
In the first experiment we performed the synthesis of 2 = 28 element antenna array in order to get the minimum side lobe level and nulls at  1 = 30 ∘ ,  2 = 32.5 ∘ , and,  3 = 35 ∘ .In this paper we performed experiments for  = 0.05 and 0.1 (with respect to /2) where  11 = 4 ∘ and  12 = 90 ∘ .
The cost function for the algorithm following expression was used: We can see that in this paper cost function is not same as that in paper [23].The main reason why the form of the cost function (9) was taken, rather than form (10), is that the value of the second part of cost function tends towards negative infinity for AF(  ) norm = 0.This case would adversely affect the search in such a way that the search for the second part of (10) would be predominant.In this case, we took normalized International Journal of Antennas and Propagation array factor and therefore it is unlikely that such a case could occur.
The radiation pattern has been shown in Figure 2. As it is seen from Figure 2 obtained results give lower side lobe level then results from paper [23].Region [30 ∘ , 35 ∘ ] is under −40 dB and it is lower than obtained region from published work.The optimal locations with respect to  obtained by using the provided algorithm are shown in Table 1.In Table 2 SLL, FNBW, and nulls depth at  1 = 30 ∘ ,  2 = 32.5 ∘ , and  3 = 35 ∘ are shown.Convergence characteristics in terms of side lobe level versus the number of iterations are shown in Figure 3.It is observed that the cost function value converges to the optimum result quickly.
The observed results were better for  = 0.05.From Table 2, it is seen that the results provide side lobe level lower than −18 dB, with nulls deeper than −100 dB.Minimum side lobe level of the array optimized using presented algorithm showed improvement (4 dB) over the published works GPS [23], ACO [18], and CSO [22] (Table 3).All three nulls are deeper than −100 dB and all nulls are deeper than nulls from published works.The FNBW is slightly higher than previously published works.
In the second experiment, we performed the synthesis of 2 = 20 elements' antenna array with fixed excitation Dolph-Chebyshev amplitudes, which form the side lobe level at −30 dB.Excitation amplitudes have values = [1, 0.97, 0.912, 0.831, 0.731, 0.620, 0.504, 0.391, 0.285, 0.325].Null is located at the position of  1 = 20 ∘ , whereas  11 = 8 ∘ and  12 = 90 ∘ .For the objective function the following relations were taken: Mutual coupling effects were under consideration during the processing.In this experiment, we fixed  = 0.1.The optimal locations with respect to  obtained by using the provided algorithm are shown in Table 4 together with the results obtained in [23].Figure 4 shows a radiation pattern produced by the presented method together with radiation pattern obtained in the work [23].We can notice that the resulting null is at the position of 20 ∘ .In Table 5 SLL, FNBW, and null depth at  1 = 20 ∘ are shown.We can also notice that SLL obtained by presented method is −0.8 dB lower than the examined results obtained in [23].Radiation pattern obtained with the presented method has a deeper null then the published work [23].FNBW has a negligible difference.

Conclusion
In the presented work, we proposed a new algorithm for finding the optimal position of the 5G mm wave antenna array elements which minimize side lobe level and placing nulls in desired directions.In both experiments, suppression of SLL was achieved, while nulls were maintained in specified directions and cost function value of presented method converges to the optimum result quickly.Given algorithm provides better results than the compared methods and shows good potential for the 5G mm wave antenna synthesis.GPS [23] GA [23] Example 1

Figure 2 :
Figure 2: The normalized radiation pattern of 28-element linear array.

Figure 3 :
Figure 3: The convergence characteristics of side lobe level for optimization of a 28-element linear array.

Figure 4 :
Figure 4: The normalized radiation pattern of 20-element linear array.

Table 1 :
Element position of the 28-element linear array (normalized with respect to ).

Table 3 :
Comparison with published results for 2 = 28 elements optimized with respect to position.

Table 4 :
Element position of the 20 element linear array (normalized with respect to ).

Table 5 :
Comparison with published results for 2 = 20 elements optimized with respect to position.