Abstract
An improved antenna array pattern synthesis method using fast Fourier transform is proposed, which can be effectively applied to the synthesis of large planar arrays with periodic structure. Theoretical and simulative analyses show that the original FFT method has a low convergence rate and the converged solution can hardly fully meet the requirements of the desired pattern. A scaling factor is introduced to the original method. By choosing a proper value for the scaling factor, the convergence rate can be greatly improved and the final solution is able to fully meet the expectations. Simulation results are given to demonstrate the effectiveness of the proposed algorithm.
1. Introduction
In order to solve complex antenna pattern synthesis problems, various methods using various optimization algorithms have been developed. In [1], a quadratic program is formed for arbitrary array pattern synthesis. In [2], a convex optimization problem [3] is formulated for pattern synthesis subject to arbitrary upper bounds. For certain cases, the convex programming problem can be reduced to a linear programming problem [4]. In [5] an effective hybrid optimization method is proposed for footprint pattern synthesis of very large planar antenna arrays [6]. The methods mentioned above all adopt conventional optimization algorithms [7]. Global optimization algorithms such as genetic algorithms [8–10] and particle swarm optimization algorithms [11–13] have also been successfully applied in pattern synthesis problems. To approximate the desired pattern for an array, we need to discretize the angular space. The bigger the elements number is, the greater the required discrete density needs to be. Usually, the computational complexity would grow greatly as the elements number increases. As a consequence, normal synthesis techniques using local or global optimization algorithms are usually not suitable for large planar arrays.
As we know, fast Fourier transforms (FFT) are able to quickly compute the radiation pattern of an array with periodic structure. Once the number of FFT points is specified, the computation time is barely affected by the element number. In [14], an FFT method suitable for large planar arrays is proposed. The operation is very straightforward, which mainly involves direct and inverse fast Fourier transforms. In [15], a modified iterative FFT technique is proposed for leaky-wave antenna pattern synthesis. In [16], FFT is used for the pattern synthesis of nonuniform antenna arrays. The iterative FFT method is very efficient as shown in [14], and many examples are presented, but why the method is effective has not been explained.
In this paper, both theoretical and simulative analyses of the FFT method are presented. It is found that, though effective, it is a slow-convergence method and can hardly converge to the optimum solution. Based on the analyses, we introduce a scaling factor and the performance can be greatly improved.
2. Original FFT Method for Pattern Synthesis
First, we establish the relationship between a certain point in the FFT result and the corresponding angle of the radiation pattern. Consider a planar array with elements arranging in a rectangular grid and spacing and between rows and columns. Assume the element pattern is isotropic. The array factor is given bywhere is the complex excitation of the th element and is the wavelength. If the coordinates are used, the array pattern can be written asPerforming points 2D inverse fast Fourier transform (IFFT) on the excitations, we getIf we want to represent the array factor using , then the coordinates and are related bywhere and are integers, making sure both sides of the equations have the same value range. Suppose that both and are even numbers. DefineCombining (4) and (5) and considering the value ranges of and , we haveSo, the relation between the array factor and the IFFT of the array excitation isFinally, the visible space is given byNote that the indices change in (5) is in fact the fftshift operation in Matlab.
Once the corresponding relationship is established, the procedure of the FFT method for pattern synthesis is given as follows.(1)Specify dimension of the array , the initial excitation , the FFT points , and the maximum iteration times .(2)Perform IFFT on the excitation of the th iteration and obtain the array factor .(3)Extract the amplitude and phase of .(4)Compare with the desired pattern . If the computed pattern fully meets the requirements or the maximum iteration time is reached, terminate the procedure; otherwise go to the next step.(5)Obtain the new pattern by replacing the undesired with as follows: where is a set containing all the points where the pattern is undesired. For example, if is a point within the side lobe region and , then it means that the side lobe level exceeds the desired level and that .(6)Perform points 2-D FFT on and choose the first points as the initial excitation for the next iteration . Constraints can be easily made for amplitude-only or phase-only synthesis.(7)Go to step .
The procedure is also illustrated in Figure 1.
3. Algorithm Analysis
The procedure is simple yet effective. In this section, the method is carefully examined. Without loss of generality, take an -element linear array for example. Suppose that in the th iteration the excitation sequence is which has been extended to points by zero padding. is the number of the FFT points. After comparing the obtained pattern with the desired pattern, we can obtain the error pattern . Now we will find the relation between the error pattern in the next iteration and . It is difficult to precisely compute through the above assumption. Let be the corresponding error excitation for . Thenwhere contains points. Consider the elements in If the set only contains points representing the side lobe region, then the magnitude values of the elements in are very small, as well as . So the variation in from is very little. Then we can use as the desired pattern in the th iteration and obtain the error pattern
For complex weighting using both amplitude and phase, is given aswhere is a sequence by setting the first elements in to zero. SoAs a result, we havewhere denotes the -norm.
It is seen from (15) that the method has the ability to converge. However, as increases the difference between and decreases. We can conclude that the method has a lower converge rate for larger FFT points. For amplitude-only synthesis, the analysis is similar and (15) still holds. For phase-only analysis, the excitation of the th iteration iswhere is the given amplitude. Since the magnitude value of can be very small compared to , can be approximated bySo the convergence properties are similar to amplitude-only or complex tapering.
To verify our analysis, we will give an example. Consider a linear array with 60 elements spacing by half wavelength. For the amplitude-only synthesis, the desired pattern has a maximum side lobe level of −40 dB. For the phase-only case, the maximum side lobe level is −18 dB. We use the uniform taper as the initial excitation and set the maximum iteration number to 5000 and 10000 for the two cases, respectively. Figure 2 shows the convergence properties of the norm of the error pattern. We can see that the method has a low convergence rate, as it converges logarithmically in the latter part of the iteration. We can also see that the bigger the FFT points number is, the lower the method converges. Figure 2 also illustrates that although the norm of the error pattern can converge to a very low value, it does not reach zero. It means that the method can hardly synthesize a pattern that fully meets the requirements of the desired pattern.
(a)
(b)
4. Improved FFT Method for Pattern Synthesis
For planar arrays, the elements in the error pattern are given byFirst we divide into and , which are the sets representing the side lobe region and the main lobe region, respectively. When , and the value tends to converge towards zero more slowly as the procedure iterates. So, to change the slowly varying characteristics, we propose the following updating equation for the error pattern:where is a scaling factor within the range . The smaller the value of is, the more intense the excitation changes. In (19), when , the error pattern remains unchanged, since for main lobe shaping synthesis, we need to approximate to a certain error level, whereas in the side lobe region, only is required. So, the final updating equation replacing (9) is given by
To analyze the convergence properties of the proposed method, we still use the prior examples and consider a linear array with 60 elements spacing by half wavelength. The maximum iteration times are 5000 and 10000 as before. As illustrated in Figures 3 and 4, although the curves are not smooth, the norm of the error pattern can converge to zero for all the cases. In both figures, we can also see that the iteration times of the case are much less than that of . Figures 5 and 6 show the influence of value on the maximum iteration times and the array directivity. It is shown that as grows near one, the iterations times increase rapidly but the directivity only changes a little. So we can say that the proposed method can give a good performance by setting to zero. The radiation patterns obtained using are shown in Figure 7, where is measured from one end of the linear array to the other end.
(a)
(b)
(a)
(b)
(a)
(b)
(a)
(b)
(a)
(b)
5. Example for Planar Array Pattern Synthesis
Consider a planar array with elements spacing by in both directions. The desired pattern has a maximum side lobe level of −40 dB and two −60 dB notches at rectangular sectors and . The procedure suggested in [14] is adopted. First the points FFT is used. At the following two phases, and points FFT are performed, with the previous results being the initial excitations for the next phase. At each phase, the maximum iteration time is 1000, and the scaling factor is set to zero. Figure 8 gives the convergence property of the number of undesired radiation points. It shows that, after the previous two phases, the final phase converges very quickly.
The 1D radiation patterns containing the two notch sectors are shown in Figure 9. Figure 10 presents the 2D radiation pattern, where the visible space is determined by (8). The normalized amplitude and phase of the excitations are given in Figure 11.
6. Conclusion
The iterative FFT method is capable of synthesizing a large planar array. In this paper, the method is validated by theoretical and simulative analyses. But the original method has a low convergence rate, and the synthesized results are usually unable to achieve the optimal solution. A scaling factor is introduced to form an improved method, which can avoid the drawbacks of the original method. Analysis and simulation results showed the effectiveness of the improved method.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.