Optimization of a 12 . 5 GHzMicrostrip Antenna Array Using Taguchi ’ s Method

This paper presents Taguchi’s optimization method implemented in the design of a single feed (without any matching network) 5× 5 microstrip antenna array operating around 12.5 GHz. The proposed optimization method is statistical and is widely used for quality assurance in many fields such as mechanical and chemical production, consumer electronics, services; however it has been underused in the field of electromagnetics. It allows optimization of multiparameter, multitarget complex designs in a very short time in conjunction with advanced simulation tools. The proposed antenna has been fully evaluated under Taghuchi’s and PSO’s optimization methods, and the experimental results show total Gain of 15 dB, and good matching with S11 better than −20 dB, in the frequency range 12.3 to 12.8 GHz.


Introduction
The demands of modern telecommunications increase hence, more complex antenna array configurations are required.Optimization and synthesis of microwave antenna arrays and other complex electromagnetic-based structures typically deal with multiobjective functions that are highly nonlinear, discontinuous, and have a large number of optimization parameters.In addition, due to strong mutual coupling and other propagation effects, they require a simultaneous optimization of the design parameters.
The above optimization methods can be divided into two categories: global and local techniques.Global techniques such as EP, GA, PSO, and ANN are capable of handling multidimensional, discontinuous, and nondifferentiable objective functions with many potential local maxima while they are largely independent of initial conditions.However, a main drawback is that the convergence rate is slow [10].
In contrast, the main advantage of local techniques such as the gradient-based algorithms is that the solution converges rapidly.However, local techniques work well only for a small number of continuous parameters highly depending on the starting point or the initial guessing while they reacting are relatively poorly in case of discontinuities in solution spaces [10].
In order to bridge the weak points of these two techniques a new statistical optimization method is necessary.A way of achieving target optimization and reduction in variation around the target is to apply the Design of Experiment (DOE) technique [11].DOE is a powerful statistical technique for improving product or process design as well as for solving production problems.A standardized version of the DOE has been introduced by Taguchi, an easy to learn and apply technique for design optimization and production problem investigation [12,13].Taguchi's optimization technique can handle multidimensional, discontinuous, and International Journal of Antennas and Propagation nondifferentiable objective functions with many potential local maxima whilst converges rapidly to the optimum result but within a well-defined area [14].
The goal of Taguchi's method is the prediction of the parameters that cause variability (noise) and the optimum value of each parameter (control), in order to achieve robustness of the device under test.The Taguchi approach has the great advantage that only a small fraction of all possible parameter-level combinations is tested, based on the theory of highly fractional orthogonal designs (orthogonal arrays).
Previously published work in the field of antenna array using Taguchi's technique had as optimization goal either pattern or input matching criteria.Work on radiation patterns includes optimization of the spiral antenna gain [15] as well as optimizing linear arrays for nulls at specific directions, suitable pattern shapes, and suppressed side lobes [16].Work on input matching includes two-layer wideband patch antenna [17], ultrabroadband zigzag logperiodic antenna [18], and a CPW slot antenna [11].This paper presents, for the first time, the optimization of a simple microstrip planar 5 × 5 antenna array configuration with interconnected elements in order to achieve good matching (S11 < −10 dB) without using of any matching networks and high gain, suitable for the 12.5 GHz broadcasting satellite service (BSS) frequency bands.
The exact dimensions of the array are precisely calculated by applying a Taguchi's optimization algorithm in conjunction with an electromagnetic solver FEKO [19].The objective of Taguchi's algorithm is the minimization of a particular mathematical "fitness function," which is suitably determined according to the criteria mentioned above.
Furthermore, a statistical analysis of variance (ANOVA) of the optimized dimensions of the antenna is applied to the final iteration.The goal of this analysis is the investigation of the significance of each parameter within their final limits in order to fabricate a robust device.

Taguchi Optimization Method
The optimization procedure begins after some initial consideration as regards the initial conditions setting as well as the selection of a proper OA and an appropriate expression of the fitness function (FF).The selection of an OA depends on the number of input parameters and the number of levels for each parameter.The FF is a particular mathematical function and is developed according to the nature of the problem and the optimization goals.
After a simple analysis, the simulation results serve as objective functions for optimization and data analysis, and an optimum combination of the parameter values can be obtained.The log functions of the outputs, named by Taguchi as Signal-to-Noise ratios (S/N), are used for the prediction of the optimum result.It can be demonstrated via statistics that although the number of experiments are dramatically reduced, the optimum result obtained through the orthogonal array usage is very close to that obtained by making use of the full factorial approach.
When the Taguchi method is implemented at the design level and the efforts are focused on the optimization of the control values, the experiments can be replaced with simulations.
In order to achieve as high convergence with the goal as possible, successive implementations of the method have to be applied.Under this procedure, the optimum results of the last iteration serve as central values for the next, reducing each time with a predefined factor the level difference of each parameter.The procedure terminates when the level difference becomes negligible, and maximum available accuracy has been reached.
The procedural steps in detail are shown below.
(1) Consideration of the problem that must be solved.
(2) Extraction of the FF and definition of the optimum goal (minimum, nominal, or maximum).
(3) Definition of the main parameters and their estimated (center) values.
(4) Definition of the levels limit for each parameter within ±10% of the center values.In order to describe the nonlinear effect so as to gradually minimize each iteration level's difference, an odd number of levels must be used for each input parameter.
(5) Definition of the maximum resolution of the parameters.
(6) Design of experiment (DOE) using Taguchi's suggested orthogonal arrays OAn(mk) in order to minimize the effect of any erroneous assumptions that have been made due to effects considered negligible, which consist of (i) n rows (number of experiments), (ii) k columns (number of parameters), and (iii) m levels (on which each parameter will vary).
(7) Simulation using FEKO according to the selected OA.
(8) Evaluation of the compliance of the FF for each combination of the levels of parameters based on the simulation results.
(9) Computation of the mean value of the fitness functions of the experiment (10) Computation of the mean value for each level of each parameter (Example: for the parameter A when level is 1, add the values of all corresponding FF and compute their mean value).
(11) Consideration of the optimum level for each parameter depending on the Y mi and the nature of the goal (minimum, nominal, or maximum). ( The predicted value might not be the optimum because the OA is a fractional factorial design, but, nevertheless, it shows the direction of the optimization.During the next iterations, as the gap between the mean and optimum predicted value becomes smaller, the possibility that the optimum predicted value to be the real optimum value rises significantly. (13) Definition of the reducing percentage (RP) of the initial difference between the levels of the parameters.
The RP depends on the nature of the problem and can be high for simple cases with only one optimum condition or low for more complex situations.( 14) Creation of new level differences (LDs) by multiplying the RP with the initial level of the parameters (15) Creation of new levels for the next iteration by adding the estimated optimum levels of the parameters of the 1th iteration with the LD i .
(16) The procedure stops when the LD i reaches the limits of the allowed resolution of the parameters.
(17) Decision, on which parameters are significant based on their final levels, using statistical analyses of variance (ANOVA).
(18) Construction of the antenna taking special care of the significant parameters as they will be responsible for any variability in the antenna characteristics.

Design of the Antenna
A 5 × 5 microstrip planar array is considered for the study.The array is tuned to work for the 12.5 GHz Broadcasting Satellite Service (BSS) frequency bands.This antenna array consists of equally spaced rectangular patches joined together with microstrip lines and is fed through a simple microstrip line without any matching network, as shown in Figure 1.
In this paper, Taguchi's method will be used to optimize the geometrical dimensions of the antenna so as to simultaneously achieve good matching, low side lobe level, and high gain in the operating frequency range.The substrate chosen for the antenna implementation is the microwave laminate RT/duroid 5880 from Rogers with er = 2.2, tan δ = 0.0009, and thickness of 1.575 mm.
Thus, a weighted FF has been chosen with the form where the weight factors (1000, 250) indicate the significance of the targets, and factor 15 indicates the optimum result (goal + 50%) for input matching and gain.There is also a dead zone in the pattern cut ϕ j(0 ≤ ϕ j ≤ 2π), θ j(6 • ≤ θ j ≤ 14 • ), and θ j(−14 • ≤ θ j ≤ −6 • ) where the results are ignored as they are irrelevant to gain or SLL results.It has to be mentioned that the computation takes into account the case of the first two factors being negative, that is, |S11 (dB) | > 15 and Gain (dB) > 15 and keeps them to zero otherwise; the FF could be driven to false results, discriminating towards one or the other characteristic.
In the design of an array antenna, the most important design parameters are usually the number of elements, spacing between the elements, excitation (amplitude and phase), half-power beamwidth, directivity, and side lobe level  [20,21].Hence, the initial center dimensions for each parameter can be defined as shown below.
The free air wavelength at 12.5 GHz is λ 0 =24 mm.Considering a ± 20% deviation from the center values, initial levels can be created, as shown in Table 1.For a realistic approach in the antenna design, the resolution for all dimensions is set to 1 μm.For an OA with 5 parameters of 3 levels for each parameter, a configuration with at least Taguchi suggests two solutions (i) the OA 18 (37, 2) that can handle up to 7 parameters with 3 levels each and one with 2 levels in an array of 18 rows.
(ii) the OA 27 (313) that can handle up to 13 parameters with 3 levels each in an array of 27 rows.
Since achieving the desired gain and obtaining good matching is a complex task, it was decided to employ a larger OA to increase the possible combinations.Thus, the second solution appears more promising (Table 2) as it has 50% more rows for each OA and offers higher confidence level at the expense of some extracomputation time (11.5 h instead of 7.5 h for this case study).
According to the complexity of the design, a reducing percentage of RP = 0.1 has been assigned.The procedure terminates when the LD is below 1 μm (minimum resolution) after twenty-five optimization cycles (27 each), finally reaching 675 iterations.
FEKO is a software suite built for the analysis of a wide range of electromagnetic problems.Applications include EMC analysis, antenna design, microstrip antennas and circuits, dielectric media, and scattering analysis.
The core of the software package is based on the method of moments (MoMs).The MoM is a full wave solution of Maxwell's integral equations in the frequency domain.The active area of the microstrip array antenna design of this case study consists of 25 metallic, equal in size, squared patches separated by equal connector lines, and a single microstrip feeding line as shown in Figure 1.The calculation of the antenna elements as well as their simulation has been carried out in EDITFEKO module, which is capable of simplifying the process of generating and editing * .pre(PREFEKO input) files as this is a standard text editor that can use ASCII text files with customized functionality.For such a complicated iterative design, the editor EDITFEKO is chosen to simplify the process of generating and editing * .pre(PREFEKO input) files as this is a standard text editor that can use ASCII text files with customized functionality.After meshing the active area of the antenna has been separated into 944 triangles, as illustrated in Figure 2.
To model the dielectric area under the metallic patches, FEKO uses Green's function formulation which implements 2D infinite planes with finite thickness (h = 1.575 mm, er = 2.2, tan δ = 0.0009) to handle the dielectric layer.Under the dielectric layer, real ground layer is modeled using reflection coefficient approximation to complete the microstrip antenna design.Finally POSTFEKO module is used to visually confirm the correctness of the model before starting a potentially time demanding FEKO simulation.It is also used to display the results, of the FEKO simulation, either in the form of position in a 3D window or in arbitrary 2D graphs.
Taking into account the above considerations in combination with FEKO simulation results the Taguchi's optimization algorithm for the 5 × 5 antenna array of Figure 1 is implemented in C++.The use of C++ code offers the ability to alter the antenna's parameters since all the parameters can be entered in a text format and get processed by the program.

Results
For comparison reasons, the parameters of the antenna will also be optimized using the well-known Particle Swarm Optimization (PSO) method.The particles and the FF of the PSO are set to be the same as those of Taguchi's method.Therefore, the optimization performance of both optimizers could be compared to each other according to the total number of iterations required towards the design goal.
For a fare comparison, 675 iterations are also conducted by PSO choosing the best FF every after 27 iterations to compare it with that of Taguchi's method.The optimization procedure graph, shown in Figure 3, presents the curves of the best FF values for every optimization cycle, as they converged under Taguchi and PSO processes.The results show that Taguchi's method excels PSO (some 1 dB less as both algorithms are configured as minimizers), when a well-defined narrow searching area is under consideration.
The results for optimum dimensions extracted after 25 cycles of iterations (less than 4 hours of processing time) for the 12.5 GHz single feed planar antenna array are illustrated in Table 3.
In order to fabricate a good performance and robust design taking into account the fabrication tolerances, the selected OA 27 (35) with the final results are examined under   a statistical analysis of variance (ANOVA) using MINITAB 14 software package.
To execute this analysis, Table II, which includes the levels status of OA 27 (35) corresponding to the antenna parameters, is transferred to a Minitab worksheet where the five parameters model has been assigned to the five columns of the parameters (A, B, C, D, E).The results of the final iteration of the Taguchi's optimization procedure (20•Log10 (FF) values) have been assigned to the next column as response and have been named as 12.5 GHz antenna array.
According to the results which are illustrated in Table 4, only two of the parameters (B.connector length and C. connector width) are significant in the final iteration with a predefined significance level α > 0.005, and they are responsible for the 83.15% of the total variability of the mean value.Special care should be taken during the fabrication procedure of the antenna to keep these two parameters close to the optimum values in order to achieve a robust device.
The analysis of variance table gives, for each term in the model, the degrees of freedom (DF), the sequential sums of squares (Seq SSs), the adjusted (partial) sums of squares (Adj SSs), the adjusted means squares (Adj MSs), and the F-statistic from the adjusted means squares and its P-value, which shows the significance of each parameter within its levels under a predefined level of significance.The same results from ANOVA are presented graphically with the main effects plot (fitted means (dB)) in Figure 4.
Taking into account the above precautions, the antenna array with its optimum dimensions has been fabricated on RT/duroid 5880 microwave laminate, as shown in Figure 5.
Simulations have been carried out using FEKO software package, and the results show high gain ≈ 15 dB, as shown in Figures 6 and 7 for planes A and B, high SLL ≈ 13 dB in B-plane, and very good matching S11 ≈ −7 dB for the central frequency of 12.5 GHz.A small asymmetry has been observed in the pattern of plane B (Figure 7), due to asymmetric current distributions as shown in Figure 4.An array design with a feeding arrangement at the centre element will eliminate this problem.Also, a higher gain closer to 18-20 dB was not achieved as a nonuniform illumination resulted from this simple array configuration.The antenna array gain has been simulated as well as measured in the frequency range 12 to 13 GHz, presenting excellent results in the pattern cut ϕ j(0 ≤ ϕ j ≤ 2π), θ j = 0 • (Gain > 10 dB), as illustrated in Figure 8.
The measurements of S11 have been carried out using a network analyzer (Agilent Network Analyzer N5230A) and shown in Figure 9.A comparison between the simulations and measurements results of S11 has been performed, which shows good agreement, as illustrated in Figure 10.

Conclusions
Analytical description of the Taguchi statistical optimization method for electromagnetic applications has been presented and implemented in the design and fabrication of a 12.5 GHz single feed 5 × 5 planar microstrip antenna array.Following this method, complicated optimization goals such as matching and high directivity and low SLL have been accomplished and verified by the agreement between simulations (FEKO) and measurements (Agilent Network Analyzer N5230A).Besides, the antenna array gain has been measured in the frequency range 12 to 13 GHz indicating high gain between 12.3 and 12.8 GHz, in good agreement with theory and simulation results.Finally, the efficiency measured was very satisfactory throughout the frequency range tested.It has to be mentioned that no feeding networks are used in this configuration.In order to create a robust construction, a statistical analysis of variance (MINITAB) has been performed and the most significant parameters of the antenna design have been discovered.This work is still ongoing, and an extended version of Taguchi's optimization method is under consideration with a wider consideration area for each parameter, and the results will be presented in the near future.

( 12 )
Prediction of the optimum value of the experiment's FF, based on the 20•Log10 values of Y and the Y mi .(The conversion is essential in order to avoid negative values especially at the beginning, when the differences between Y and Y mi are high)

Figure 2 :Figure 3 :
Figure 2: FEKO-meshed design of the array antenna with current distribution.

FF mean (dB) of 12 . 5 Figure 4 :
Figure 4: Main effects plot of the antenna's parameters variance at the achieved optimum level.

Table 1 :
The initial levels for the optimization procedure.

Table 2 :
The OA 27 53as configured with the parameters levels and the responses of the final FF.

Table 3 :
The optimized values of the parameters.

Table 4 :
The analysis of variance (ANOVA) table.