International Journal of Antennas and Propagation

Volume 2014 (2014), Article ID 392895, 12 pages

http://dx.doi.org/10.1155/2014/392895

Research Article

Direction Finding Using Multiple Sum and Difference Patterns in 4D Antenna Arrays

Quanjiang Zhu, Shiwen Yang, Ruilin Yao, and Zaiping Nie

School of Electronic Engineering, University of Electronic Science and Technology of China (UESTC), Chengdu 611731, China

Received 4 January 2014; Accepted 19 March 2014; Published 21 May 2014

Academic Editor: Paolo Rocca

Copyright © 2014 Quanjiang Zhu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Traditional monopulse systems used for direction finding usually face the contradiction between high angle precision and wide angle-searching field, and a compromise has to be made. In this paper, the time modulation technique in four-dimensional (4D) antenna array is introduced into the conventional phase-comparison monopulse to form a novel direction-finding system, in which both high angle resolution and wide field-of-view are realized. The full 4D array is divided into two subarrays and the differential evolution (DE) algorithm is used to optimize the time sequence of each subarray to generate multibeams at the center frequency and low sidebands. Then the multibeams of the two subarrays are phase-compared with each other and multiple pairs of sum-difference beams are formed at different sidebands and point to different spatial angles. The proposed direction-finding system covers a large field-of-view of up to ±60° and simultaneously maintains the advantages of monopulse systems, such as high angle precision and low computation complexity. Theoretical analysis and experimental results validate the effectiveness of the proposed system.

1. Introduction

Monopulse systems, which were evolved from sequential lobing, are capable of eliminating the errors caused by amplitude fluctuations of target echoes and increase the data rate, since that angle information can be derived from a single pulse [1]. Monopulse has great advantages, such as very high angular precision and very low computational complexity and has been used in many applications, such as air traffic control, missile tracking, and space antennas. A classic and simple type of monopulse is the phase-comparison monopulse, which uses the measured phase differences of two antennas’ outputs to establish the target bearing. An inherent problem with such monopulse system is that the resulting angle information can be ambiguous when the spacing of the two antennas or arrays is more than half a wavelength, which is often the case for high precision measurement [2, 3]. One of the common ways to solve this ambiguity problem is to use staggered baseline length [3, 4], which is to make another measurement with a different antenna spacing (and thus different ambiguity) and combine the two ambiguous measurements into one unambiguous measurement. Other ambiguity resolution techniques include the use of multiple frequencies, baseline rotation or frequency modulation, and so forth. However, these methods either require extra antenna element, receiver electronics and additional installation space, or require the increased complexity or computation time.

As compared to the aforementioned techniques, the four-dimensional (4D) antenna array provides another method to resolve the multiple ambiguities. The 4D antenna array uses time as the fourth design parameter to alleviate the stringent requirement on traditional antenna arrays operating in the 3D space, which is realized by time-modulating the antenna arrays with high speed RF switches. Since its inception in the 1959 [5], many researchers have paid great attention to this novel technology, especially during the past decade [6–19]. As a result of time modulation in 4D arrays, sidebands are generated at multiples of the time modulation frequency. These sideband signals can either be suppressed or be enhanced through the optimization of the exciting time sequences [7–13]. In [14], the sidebands were purposely utilized and a two-element 4D array was configured to operate as a direction-finding system, in which the sum pattern was formed at the center frequency and the difference pattern was formed at the first sideband. Relevant experiment was then conducted in [15], and the measured results confirmed the theoretical analysis. Some monopulse systems based on the 4D arrays have been conceptually designed and numerically verified [16–19]. However, the contradiction between precise angle estimation and large field-of-view was not addressed in these papers. In [20], the time modulation technique in 4D arrays was firstly proposed as a method of resolving the ambiguities existing in phase-comparison monopulse systems. However, no experiment is performed to validate the idea, and the signal processing units are not addressed.

In this paper, a novel phase-comparison monopulse based on the 4D antenna array is presented. The sidebands of a 4D antenna array are enhanced and utilized not only to generate multiple beams but also to resolve the multiple ambiguities. The proposed technique makes it possible to realize unambiguous high accuracy direction finding over a large field-of-view of more than ±60°. This paper is organized as follows. In Section 2, the conventional phase-comparison monopulse and the 4D arrays are introduced. In Section 3, the novel monopulse system based on a 4D array is presented and its working principle is analyzed. In Section 4, an experimental system composed of 4D arrays is built, and the experimental results and corresponding signal processing steps are presented in detail. Finally, some conclusions are drawn.

2. Theoretical Backgrounds

2.1. Conventional Phase-Comparison Monopulse

There are two types of monopulse systems, namely, amplitude-comparison monopulse and phase-comparison monopulse. The phase-comparison monopulse is defined in terms of receiving beams with different phase centers [1] and is adopted as the type of monopulse in this paper.

As illustrated in Figure 1, an array model with 8 isotropic elements is used to demonstrate the basic principles of phase-comparison monopulse. The full array is divided into two subarrays. When a plane wave impinges on the array with an angle measured from broadside, the received signals or voltages of the two subarrays can be written as

where is the received signal voltage of a single element, is the element space, is the subarray distance, and refers to the free space wave number at the operating frequency. It can be seen that the received signals of the two subarrays have the same amplitude, and the only difference is the phase response. The phase difference of the two signals with respect to the center of the array equals , which is directly related to the direction of coming wave.

Figure 1: A phase-comparison monopulse system based on a conventional 8-element array.

The sum and difference of the two received signals are given as:

The difference-sum ratio between the two received signals, , is defined as

Equation (3) is the basic direction-finding equation, which relates the difference-sum ratio to the bearing angle , with the distance of the two subarray as the parameter. The relationship among the three parameters , , and is presented in Figure 2. As can be seen from Figures 2(a)–2(d), ambiguity resolution is required as long as the subarray distance is greater than half a wavelength. In fact, since tangent and sine functions are both periodic functions, any that satisfies

is a valid solution for (3).

(a)

(b)

(c)

(d)

Figure 2: Difference-sum ratio for different subarray distance.

As the distance between the subarrays becomes larger, the slope of the difference-sum ratio becomes steeper and more zeros and poles are formed. For the instance of (shown in Figure 2(d)), the slope of is the steepest of the four cases (meaning the highest measurement accuracy), and three zeros (°, 0°, and 30°) and four poles (°, −14°, 14°, and 48°) are formed. For any particular value of , four ambiguous bearing angles correspond to it. If the angle-searching range is inside the field from −14° to 14°, it is easy to determine the angle as is monotonic inside this small field. However, it turns to be difficult when the searching range is beyond the field, due to ambiguity problem.

2.2. 4D Antenna Arrays

Unlike conventional arrays operating in 3-dimentional (3D) space, the 4D antenna arrays are formed by introducing a fourth dimension, time, into the array design, which can be realized by time-modulating conventional arrays with RF switches. Figure 3 shows the basic topology of an 8-element 4D antenna array at the receiving mode. The 4D linear array consists of 8 isotropic elements with an equal element space . The received signal of each antenna element is time-modulated by the high speed RF switch and then summed to form the output. The RF switch is controlled by the circuit board programmed for specific time sequences. Compared with the conventional arrays, the RF switches have the effects of amplitude and phase weighting by the control of on-off state of each antenna element [21].

Figure 3: Block diagram of an 8-element 4D antenna array.

For a plane wave with a frequency of and the incident angle of with respect to the broadside of the array, the array factor of the 4D array can be given as

where represents the periodic “on-off” state of the th element with a time period of . For different time modulation schemes, has different forms of expression. In this paper, the time modulation scheme of pulse shifting proposed in [22] is adopted. Accordingly, is given by

where and represent the normalized switch-on instant and switch-on duration for the th element, respectively.

Equation (5) differs from the array factor (1) of a conventional array by the time factor . As is a periodic function with the time modulation period , can be decomposed into a Fourier series, given by [22]

where represents the complex excitation weighting for the th element at the th sideband and is the time modulation frequency. Taking into account (5) and (7), (4) can be rewritten as

and can be viewed as a superposition of multiple patterns at multiple frequencies separated by the time modulation frequency and is analogous to the array factors of conventional arrays.

3. Monopulse Based on 4D Antenna Array

3.1. Basic Theory

In theory, the contradiction between the angle measurement precision and large field-of-view exists in any phase-comparison monopulses. To solve this problem, the monopulse technique is combined with the 4D antenna arrays. 4D antenna arrays have been used to generate multiple beams [12, 23], and, more importantly, these beams are formed at different sidebands and point in different directions. Inspired by this unique feature of 4D arrays, a novel phase-comparison monopulse based on a 4D array is presented, as shown in Figure 4. The 8-element 4D arrays are divided into two 4-element subarrays, which are separated by a certain distance ( is supposed to be ). The two subarrays are excited with the same time sequences, implying that the exciting time sequence of element () is the same as that of element . In this way, the power patterns generated by the two subarrays would be the same except for a phase difference, and the received signals of two subarrays, after summed and subtracted with each other, would produce a similar difference-sum ratio as the one shown in Figure 2(d). In order to maximize the signal-to-noise ratio (SNR) of the received signal, a set of band-pass filters can be mounted between the antenna elements and RF switches to filter out the noises outside the center band.

Figure 4: A phase-comparison monopulse system based on an 8-element 4D antenna array.

The basic principle of the proposed monopulse system can be stated as follows. Firstly, the large field-of-view is divided into several sectors. Secondly, the multiple beams generated by the 4D array are used to cover these angle sectors, with one beam for one angular sector. Thirdly, the received signals of two identical 4D arrays are phase-compared, and multiple sum-difference patterns are formed with one sum-difference pattern for one sector. These multiple sum-difference patterns are formed at different sidebands with different bore-sight axes and cover the field-of-view as a whole.

Similar to the conventional phase-comparison monopulse, the difference-sum ratio of the 8-element 4D array at the th sideband can be expressed as

where denotes the amplitude of the difference signal at th sideband, denotes the amplitude of the sum signal at th sideband, and is the equivalent complex excitation weighting for the th element at the th sideband, calculated by (8). As can be seen, the difference-sum ratio of the 4D array is the same as that of the conventional array, independent of the sideband order , as shown in Figure 5.

Figure 5: Difference-sum ratio of the 4D array where isotropic elements are used. For a signal coming from the th sector, is used.

Different from the conventional phase-comparison monopulse, the angle measurement field for the 4D array is divided into several sectors, and each sector is covered by a pair of sum-difference beams. As shown in Figure 5, the angle bearing information can be derived according to the without ambiguity.

3.2. Multiple Sum-Difference Pattern Synthesis

In order to steer the multiple beams to the directions of zero and pole points shown in Figure 5, the DE algorithm is used to optimize the time sequence of each subarray. The cost function of the DE algorithm is set as

where the superscript denotes the sideband number; and are the calculated and desired maximum sideband level; and are the calculated and predefined beam direction; BW is the calculated beamwidth and is the predefined beamwidth to ensure the −3 dB beam crossover level; and , , and are the corresponding weighting factors for each term. The optimized time sequence is shown in Figure 6, and resulting power patterns of the two subarrays are shown in Figure 7. Note that, in Figure 7, the power patterns of the two subarrays are the same and their phase patterns differ by a phase difference of .

Figure 6: The optimized time sequence for the two subarrays.

Figure 7: The multiple beams generated by subarray 1 and subarray 2 at different sidebands with different pointing directions.

By summing and subtracting the patterns generated by the two subarrays with each other, multiple sum-difference patterns are obtained, as shown in Figure 8. As expected, seven pairs of sum-difference patterns are formed in the direction of °, −30°, −14°, 0°, 14°, 30°, and 48°, with each sum-difference pattern covering a sector of the overall field-of-view and the overall field-of-view is extended to ±60°.

(a)

(b)

(c)

(d)

Figure 8: The sum-difference patterns formed at different sidebands and different pointing directions by phase-comparing the multibeams of subarray 1 and subarray 2 (only the positive sidebands are shown).

3.3. Method for Direction Finding

Suppose that a signal impinges on the 4D array. The signal will be received by the two subarrays. The spectra of the two received signals can be computed using the fast Fourier transform (FFT) algorithm and their spectra response is then examined. The sideband that has the maximum response provides a rough estimation of signal direction. Then sum and subtract the two received signals with each other and compute the spectra of the sum and difference signals. By comparing the amplitude of the difference signal with that of the sum signal, the difference-sum ratio can be obtained, and thus the direction of the signal can be easily determined. A detailed example for direction finding will be presented in the next part.

4. Experimental Setup and Results

4.1. Experimental Setup

To validate the theory presented above, an experiment is conducted in an anechoic chamber. The block diagram of the experimental setup is shown in Figure 9. A horn operating at the frequency of 2.6 GHz is used as the transmitting antenna and is kept in a fixed position. The 4D array composed of an 8-element printed dipole and the feeding network is used as the receiving antenna and is placed on a rotating platform. The feeding network includes 8 SPST switches, two combiners, and a CPLD board that controls the on-off state of the switches. The time sequence shown in Figure 6 is programmed into the CPLD board, which is used to generate the repetition pulse signal to drive the switches. The time modulation frequency is set to be  KHz.

Figure 9: Block diagram of the experiment setup.

The two antennas are separated with a distance of 14 m to ensure that the far-field testing condition is satisfied. A sinusoid signal of 2.6 GHz is sent to the horn antenna and transmitted. The signal is received by the elements, time-modulated by RF switches, and then summed into two channels by the two combiners. The signals of the two channels are sampled by a digital oscilloscope synchronously. The oscilloscope works at a sampling rate of 20 Gpts/s, which is 7 times higher than the highest frequency of the transmitted signal and satisfies the sampling theorem. Ten periods of the received signals can be stored. The stored signals are then exported to the personal computer (PC) and postprocessed in Matlab afterwards.

4.2. Calibration of Difference-Sum Ratio

In the theoretical analysis in Parts II and III, isotropic element patterns are used, and the patterns of the two subarrays are the same except for a phase difference. In experiment, the mutual coupling effects, ground plane, and assembly errors affect the arrays’ patterns, and the patterns of the two subarrays are not the same any more.

For accurate angle measurement, the active element patterns of the dipole array are firstly measured without time modulation. The active element pattern is the radiation pattern of an element when it is located in the array with other elements matched perfectly [24], taking into account the effect of mutual coupling between array elements. The active element pattern in the H plane for the azimuth angle from −90° to 90° of each element is measured, and the patterns of element number 1 to number 4 are shown in Figure 10. As can be seen, these patterns are different due to mutual coupling effect, the ground plane, and so forth.

(a) Element 1,

(b) Element 2,

(c) Element 3,

(d) Element 4,

Figure 10: Measured normalized active element amplitude patterns of subarray 1 [denoted as ].

Based on the measured active element patterns, the sum-difference pattern at each sideband can be calibrated and calculated by

where and denote the sum and difference pattern at the th sideband and denotes the active element pattern of the th element, which is shown in Figure 10. In calculation, both amplitude and phase information are taken into account. Figure 11 shows the calibrated sum-difference pattern at , +1, +2, and +3 sideband. As compared to the patterns in Figure 8 where isotropic patterns are used, the calibrated multiple beams of the 7 pairs of sum-difference patterns still point in the directions of °, −30°, −14°, 0°, 14°, 30°, and 48°, but the patterns change more or less.

(a)

(b)

(c)

(d)

Figure 11: Sum-difference patterns taking into account the measured active element patterns.

The difference-sum ratio can also be calibrated by comparing the amplitude of the difference pattern with that of the sum pattern, given by

Figure 12 shows the calibrated for different sideband signals. As can be seen, the at different sideband are not the same any more. For a signal coming from different directions, the corresponding can be used. For instance, for a signal coming from 30°, should be utilized, since the signal is located in the beamwidth of the sum beam at the +2nd sideband (from 24° to 42°).

(a)

(b)

(c)

(d)

Figure 12: Calibrated difference-sum ratio based on measured active element patterns.

4.3. Direction Determination of Signals

Since the horn as the transmitted antenna is at a fixed position, the 4D array is rotated every 5° from −65° to 65° to simulate the situation where the signal comes from a spatial angle. To extract the angle bearing information of the transmitted signal, the exported signals are processed in the following way, as shown in Figure 13.

Figure 13: Signal processing steps in the 4D array for direction finding.

Firstly, the received signals of subarray 1 or subarray 2 are spectrum analyzed using the FFT algorithm. The sideband that has the maximum response indicates the sector of the field-of-view that the signal comes from and provides a coarse estimation of the arrival direction. The number of the sideband having the maximum amplitude response is designated as . Secondly, the received signals of subarray 1 and subarray 2 are summed up and subtracted with each other. The resulting sum and difference signals are also spectrum analyzed using the FFT algorithm. The magnitude of the th sideband of the difference signal is compared with that of the sum signal, and the difference-sum ratio is obtained. Finally, the bearing angle can be precisely determined from obtained difference-sum ratio , according to the sum-difference patterns of the th sideband shown in Figure 12.

Take the case where the transmitter is in the 30° direction as an example. The waveforms of the two received signals of subarray 1 and subarray 2 are shown in Figure 14(a), and their corresponding spectra are shown in Figure 14(b). In Figure 14(b), the maximum amplitude response appears in the +2nd sideband, which indicates that the wave comes from the sector around ° and the difference-sum ratio should be chosen as a reference. By summing and subtracting the received signals of subarray 1 and subarray 2, the waveforms of the sum and difference signals are obtained, as shown in Figure 15(a). By FFT, the corresponding spectra of the sum and difference signals are obtained and shown in Figure 15(b). By comparing the amplitude of the +2nd sideband of the difference signal with that of the +2nd sideband of the sum signal, the ratio is found to be equal to +0.137. By looking up the difference-sum ratio shown in Figure 12, the ratio +0.137 corresponds to a wave coming from the direction of °; thus the measurement error is 0.5°.

(a)

(b)

Figure 14: The procedure to determine the pair of sum-difference patterns that has the maximum amplitude response and should be used for angle measurement. (a) The waveforms of signals received by the two subarrays at °. (b) The spectra of received signals received at °, which are obtained by FFT performed in PC.

(a)

(b)

Figure 15: The procedure to obtain the difference-sum ratio : (a) waveforms of sum and difference signals at ° and (b) spectra of sum and difference signals at °.

The measurement error for a spatial field of ±65° is shown in Figure 16. For the wave coming from the direction of −60° to 60°, the measurement error is less than 2°. The error is caused mainly by phase imbalance within the two channels and also related to the SNR of the received signal. As a set of suitable BPFs are not available and not used in the experiment, the SNRs of the received signals are not high enough, especially for the signals from the edges of the field-of-view. The SNR of the received signal decreases with the increase of . The lower the SNR is, the greater the measurement error will be. In order to provide higher estimation accuracy, more careful calibrations are also needed.

Figure 16: Measured estimation errors for a spatial field of ±65°.

5. Conclusions and Discussions

The 4D antenna array is demonstrated for the application of direction finding over a wide field-of-view in this paper. The time modulation technique in 4D arrays is used to generate multiple beams at different sidebands and the phase-comparison monopulse is used to generate the sum-difference patterns. As a result, multiple pairs of sum-difference patterns with different pointing directions are formed at different sidebands. According to the theoretical analysis, an experimental system is set up and the angle measurement algorithm is calibrated taking into account the active element patterns. The experimental results show that this novel monopulse system provides unambiguous direction finding over a spatial field-of-view of ±60° with very low computational complexity.

The advantage of the proposed technique is that it can achieve direction finding over a much wider field-of-view by adding a set of RF switches to a conventional array. However, the low SNR performance of the received signals in 4D arrays increases the measurement errors. In order to realize a SNR equivalent to that of a conventional array, a set of BPFs must be added after the elements to filter out the noise outside the center band and avoid the aliasing effect, when the bandwidth of the incoming signal is wider than the modulation frequency. The SNR performance of the 4D array will be reported in the future.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was supported by the Natural Science Foundation of China (Grant no. 61125104) and in part by the 111 Project of China (Grant no. B07046).

References

S. M. Sherman and D. K. Barton, Monopulse Principles and Techniques, Artech House, Boston, Mass, USA, 2011.

W. B. Kendall, “Unambiguous accuracy of an interferometer angle-measuring system,” IEEE Transactions on Space Electronics and Telemetry, vol. 11, no. 2, pp. 62–70, 1965.

L. G. Bullock, G. R. Oeh, and J. J. Sparagna, “An analysis of wide-band microwave monopulse direction-finding techniques,” IEEE Transactions on Aerospace and Electronic Systems, vol. 7, no. 1, pp. 188–203, 1971.

E. Jacobs and E. W. Ralston, “Ambiguity resolution in interferometry,” IEEE Transactions on Aerospace and Electronic Systems, vol. 17, no. 6, pp. 766–780, 1981.

H. E. Shanks and R. W. Bickmore, “Four-dimensional electromagnetic radiators,” Canadian Journal of Physics, vol. 37, no. 3, pp. 263–275, 1959.

W. H. Kummer, A. T. Villeneuve, T. S. Fong, and F. G. Terrio, “Ultra-low sidelobes from time-modulated arrays,” IEEE Transactions on Antennas and Propagation, vol. 11, pp. 633–639, 1963.

S. Yang, Y. B. Gan, and A. Qing, “Sideband suppression in time-modulated linear arrays by the differential evolution algorithm,” IEEE Antennas and Wireless Propagation Letters, vol. 1, pp. 173–175, 2002.

S. Yang, Y. B. Gan, A. Qing, and P. K. Tan, “Design of a uniform amplitude time modulated linear array with optimized time sequences,” IEEE Transactions on Antennas and Propagation, vol. 53, no. 7, pp. 2337–2339, 2005.

S. Pal, S. Das, and A. Basak, “Design of time-modulated linear arrays with a multi-objective optimization approach,” Progress in Electromagnetics Research B, no. 23, pp. 83–107, 2010.

L. Poli, P. Rocca, L. Manica, and A. Massa, “Time modulated planar arrays—analysis and optimisation of the sideband radiations,” IET Microwaves, Antennas and Propagation, vol. 4, no. 9, pp. 1165–1171, 2010.

L. Poli, P. Rocca, and A. Massa, “Sideband radiation reduction exploiting pattern multiplication in directive time-modulated linear arrays,” IET Microwaves, Antennas and Propagation, vol. 6, no. 2, pp. 214–222, 2012.

L. Poli, P. Rocca, G. Oliveri, and A. Massa, “Harmonic beamforming in time-modulated linear arrays,” IEEE Transactions on Antennas and Propagation, vol. 59, no. 7, pp. 2538–2545, 2011.

Y. Tong and A. Tennant, “A two-channel time modulated linear array with adaptive beamforming,” IEEE Transactions on Antennas and Propagation, vol. 60, no. 1, pp. 141–147, 2012.

A. Tennant and B. Chambers, “A two-element time-modulated array with direction-finding properties,” IEEE Antennas and Wireless Propagation Letters, vol. 6, pp. 64–65, 2007.

A. Tennant, “Experimental two-element time-modulated direction finding array,” IEEE Transactions on Antennas and Propagation, vol. 58, no. 3, pp. 986–988, 2010.

J. Fondevila, J. C. Brégains, F. Ares, and E. Moreno, “Application of time modulation in the synthesis of sum and difference patterns by using linear arrays,” Microwave and Optical Technology Letters, vol. 48, no. 5, pp. 829–832, 2006.

P. Rocca, L. Manica, L. Poli, and A. Massa, “Synthesis of compromise sum-difference arrays through time-modulation,” IET Radar, Sonar and Navigation, vol. 3, no. 6, pp. 630–637, 2009.

Y. Chen, S. Yang, and Z. Nie, “Design of a novel monopulse antenna system using the time-modulated antenna arrays,” International Journal of RF and Microwave Computer-Aided Engineering, vol. 20, no. 2, pp. 163–169, 2010.

P. Rocca, L. Manica, L. Poli, and A. Massa, “Synthesis of compromise sum-difference arrays through time-modulation,” IET Radar, Sonar and Navigation, vol. 3, no. 6, pp. 630–637, 2009.

X. Huang, S. Yang, G. Li, and Z. Nie, “A novel application for sum-difference pattern detection of signal direction using time-modulated linear arrays,” in Proceedings of the 18th International Symposium on Intelligent Signal Processing and Communication Systems (ISPACS '10), Chengdu, China, December 2010.

Q. Zhu, S. Yang, R. Yao, M. Huang, and Z. Nie, “Unified time and frequency domain study of time modulated arrays,” IEEE Transactions on Aerospace and Electronic Systems, vol. 61, no. 6, pp. 3069–3076, 2013.

L. Poli, P. Rocca, L. Manica, and A. Massa, “Pattern synthesis in time-modulated linear arrays through pulse shifting,” IET Microwaves, Antennas and Propagation, vol. 4, no. 9, pp. 1157–1164, 2010.

Y. Tong and A. Tennant, “Simultaneous control of sidelobe level and harmonic beam steering in time-modulated linear arrays,” Electronics Letters, vol. 46, no. 3, pp. 200–202, 2010.

D. M. Pozar, “Active element pattern,” IEEE Transactions on Antennas and Propagation, vol. 42, no. 8, pp. 1176–1178, 1994.