Refinement Method for Weighting Scheme of Fully Spatial Beamformer

So far, a wideband spatial beamformer has been proposed. This kind of beamformer has a major contribution as its weighting coefficients are real valued in which they can be simply realized by attenuators or amplifiers. However, so far, the range of attenuation or amplification is relatively large which is not practical for hardware realization. Therefore, this paper proposes a concept to reduce the range of weighting coefficients hence, the hardware realization becomes practical. In this paper, a full prototype of wideband spatial beamformer is constructed to reflect the true beamforming performance of the proposed refinement method. Its radiation patterns obtained from simulation and measurement are compared. As a result, we can reduce the attenuation or amplification range while some radiation characteristic is remained.


Introduction
In recent years, wireless communication systems have grown rapidly and contained several technologies operating in different frequencies, for example, WiFi, WiMAX, Bluetooth, and broadband WLAN.From these reasons, the systems are demanded for better coverage, wider frequency range, higher quality of service, and more capacity [1].So far, smart antenna systems have been envisaged to be the solution.This is because the systems can form one beam towards a desired direction and create nulls or sidelobes towards interference directions.As a result, the greatly improving system performance and also saving energy can be obtained [2].The smart antenna systems usually consist of antenna array and a suitable signal-processing unit adjusting the weighting coefficients at individual antenna elements to ease the effect of interference signal [3].The signal processing works according to the utilized algorithms.So far, the beamforming algorithms for smart antennas in wireless communications have concerned only narrowband operation.However, in order to support high data-rate transmission, new efforts are currently required to deal with smart antennas in wide frequency range [4].From literatures [5][6][7], wideband smart antenna technology can be classified into three categories.The first category utilizes space-time signal processing, so called spatiotemporal beamformer.These systems consist of array antennas and tapped-delay line at each branch of the array to deal with the received signal in time domain [8].However, the number of tapped-delay line increases with the bandwidth of utilized frequency.Next category relies on space-frequency signal processing.This is an alternative approach to perform wideband beam-formation without the use of tapped-delay lines.In this method the received signal will be decomposed into non-overlapping narrowband component using band-pass filter.However, for wideband signal, the systems require a large number of frequency filters [9,10].The last category utilizes only spatial signal processing, so-called wideband spatial beamformer.This is a new method for wideband smart antennas.In the part of signal processing, weighting coefficients are calculated using two-dimensional Inverse Discrete Fourier Transform (IDFT) technique applied to the required radiation pattern.The obtained weighting coefficients become real values which can be simply realized by attenuators or amplifiers [11][12][13].The main advantage of this wideband spatial beamformer is that its design does not require phase shifters, delay circuits, or International Journal of Antennas and Propagation e l e m e n t s frequency filters.This is of considerable practical advantage over the other types of beamformer employing complex weighting coefficients [14].However, the range between maximum and minimum values of weighting coefficients (attenuation or amplification) is relatively wide.This is considerably impractical.Therefore, this paper proposes the refinement method for fully spatial beamformer.The aim of the proposed method is to reduce the range of overall weighting coefficients.As a result, the implementation of this kind of beamformer is practical.A full prototype of fully spatial beamformer is developed in order to validate the proposed method.The remainder of this paper is as follows.After brief introduction showing the background and motivation of the paper, the concept of beam formation in wide range of operating frequency band using only spatial signal processing is described in Section 2. Section 3 proposes the method to make the weights more practical, so-called refinement method.A full prototype of fully spatial beamformer is detailed and constructed in Section 4.Then, its beamforming performance is revealed with the use of proposed refinement method in Section 5. Finally, Section 6 concludes the article.

Wideband Beamforming Algorithm Using
Fully Spatial Signal Processing The relationship between M and N is [15].The radiation pattern of the array as a function of angle and frequency can be given by where f is the frequency variable, c is the speed of signal, and G( f , φ) is a function of frequency-angle-dependent gain of each antenna element.Note that signal is incident or transmitted in azimuth direction φ and θ = 90 • .When φ and θ are the azimuth and the elevation angles, respectively.
In order to determinate the weighting coefficients w m1,m2 in (1), a modified IDFT is applied to H as shown in (2) (2) We define two auxiliary functions as Where N u1 and N u2 are number of sampling points in u 1 − u 2 plane in which N u1 ∼ = 2N 1 and N u2 ∼ = 2N 2 .The H(u 1 , u 2 ) is the form of Sinc function to obtain maximum gain in the desired direction [16] which can be given as where φ 0 is the direction at which the main beam to be pointed to which is assumed to be prior known [14], α is the number related to the main beam first-null width, r l and r h are defined by the following expressions: r l = ( f l /c)d and r h = ( f h /c)d, respectively, when d = d 2 1 sin 2 φ 0 + d 2 2 co s 2 φ 0 , f l is lowest frequency and f h is highest frequency.Referring to the previous work presented in [14], null steering implementation is difficult to perform using IDF method.This is because the obtained weighting coefficients become real valued.Therefore, degree of freedom is relatively limited comparing the ones employing complex weighting coefficients.
To give some simulation examples, we consider 4 × 4array operating over frequencies from 1.9 to 2.4 GHz and choose planar monopole having omnidirectional radiation to be an antenna element.From (1) to (3), we calculate weighting coefficients when the main beam is required to be pointed to φ 0 = 45 • .The outcomes of beam formation are shown in Figures 2 and 3.
Figures 2 and 3 show the radiation pattern of the 4 × 4 spatial beamformer in which the objective function as shown in (3) is assumed using Sinc function with the mainbeam pointed to 45 • off the array broadside direction.These figures indicate that the main beams are correctly directed towards the desired direction, and the patterns are alike throughout the desired frequency band, from 1.9 to 2.4 GHz.Also, it can be observed that as frequency increases, the array's sidelobe levels increase, and the null locations slightly change with respect to the operating frequency.However, from the results, the validity of the wideband spatial beamforming concept can be confirmed.

Chebyshev Function.
In the previously section, we consider the Sinc function being as an objective function H(u 1 , u 2 ).As a result, the sidelobe level (SLL) is relatively high.Thus, in this section, the objective function is created, using Chebyshev function instead as shown in (4) [17,18].We consider where From ( 6) direction of arrival on Substituting ( 9) into ( 6) yields or In this section, it is easily done by replacing the objective function H(u 1 , u 2 ) appeared in (3) by the one appeared in (4).Therefore, SLL is expected to be lower comparing to the ones shown in Figures 2 and 3.The simulation results for this assumption are shown in Figures 4 and 5.
Comparing between Figures 4 and 5 and Figures 2 and  3 reveals that the sidelobe level when objective function is assumed using Chebyshev polynomial is lower than using Sinc function for all frequencies.In Figure 5, the null location slightly changes when frequency increases as the ones appeared in Figure 3.

Refinement Method for Weighting Scheme of Fully Spatial Beamformer
Regarding the previous section, it has shown that we can obtain a value of weighting coefficient of smart antenna by taking IDFT to the desired radiation pattern.One example of the obtained weighting coefficients using both Sinc and Chebyshev functions is shown in Table 1.As we can see they are all real valued in which we can simply utilize attenuators or amplifiers to be weights.However, as we can notice in Table 1, they are not integer numbers, and also the range between the maximum and minimum values is relatively wide.This is considerably not practical for hardware implementation.Therefore, the modification of calculated weighting coefficients is needed while maintaining some important characteristic of radiation pattern.
In this section, we describe the proposed refinement method for weighting scheme of fully spatial beamformer in order to make it practical.This method is implemented after initially calculating the weighting coefficients using IDFT as mentioned before.Some significant properties of radiation behavior are maintained as described below.

Refinement Algorithm.
The main proposal of this algorithm is to decrease the operation range of weighting coefficients as an example shown in Table 1.This algorithm is required because simply squeezing the mentioned range affects some importance radiating characteristics such as main-beam's direction and sidelobe levels.However, obtaining the integer numbers in weighting coefficients is byproduct when performing refinement method.
The weighting coefficient from IDFT method detailed in Section 2.1 is now defined as where N is number of array antennas.Usually, the value of W o is general number, then we round the weights obtained from IDFT method to integer-valued weights (W r ) as Then, the required mainbeam direction (φ 0 ), average level of minor lobe (MLL 0 ) and width of main beam (δ 0 ) and weighting coefficient (W r ) are given.Next, we carry out the following steps.
(1) Determine the maximum of W r , that is where k = 1, 2, 3 . . . is order of iteration.
(2) Decrease maximum of weighting coefficient (w (k) max) by decreasing a step size, where Δw > 0. This results in integer number for weighting coefficients.We have where j = 1, 2, 3, . . . .Then we replace w (k) max, j with w (k) max in a set of W r .We consider (3) Multiply W r by steering vector ψ, and we can get output as follows: where (p, q) is index of antenna elements.(4) Now, we can obtain the radiation pattern from (18).Then we can check some significant parameters such as mainbeam direction (φ (k) r ), average minor-lobe level (MLL (k) r ), and beamwidth (δ (k) r ) according to the following conditions: where w 0 is desired maximum of weighing coefficients, Δφ is deviation of main beam angle, ΔMLL is deviation of average minor-lobe level, Δδ is deviation of beamwidth.
(5) Repeat steps 2 to 4 until the set of weighting coefficient is zero according to (19) to ( 22), and then we can obtain a set of new solution (W r ) for the refinement method.
This proposed algorithm can run in both offline or real time processing depending on the purpose of users.In case of ruing in real-time, we need a suitable DSP board to receive the important information from DFT method then perform 5 steps refinement as shown above.Afterwards, the DSP board feedbacks the weighting coefficients to control the attenuators or amplifiers.

Simulation Results.
In this section we show the efficiency of the proposed refinement method using our own developed program in MATLAB.The utilized antennas are arranged in 4 × 4 lattice, and the operating frequency is given from 1.9 to 2.4 GHz.We also assume the decreasing step size of weighting coefficient of 1, desired maximum of weight value of 32 dB, and steering angle from 0 to 90 degrees.Some parameters are given in simulations for 6 cases as shown in Table 2. Table 2 shows 6 cases for the proposed refinement method.Some parameters, for example, maximum weighing coefficients, deviation of mainbeam angle (Δφ), deviation of average minor lobe level (ΔMLL), and deviation of beamwidth (Δδ) are differently given in each case.For case I, we set the conditions not to have any error in the parameters mentioned earlier.From running some simulations, we have found that 67 directions or cases out of 91 directions, 0 • -90 • , are succeeded for refinement.This means that there are only 67 directions out of 91 directions in which we can maintain no error in mainbeam direction, minor lobe level and beamwidth.For case II, we allow the slight error in mainbeam direction, minor lobe level and beamwidth as shown in Table 2.We have found that the number of success cases increases from 67 to 71 directions.Also, higher error in mainbeam direction, minor lobe level, and beamwidth is allowed as shown in Table 2 for cases III, IV, V, and VI.As we can see, the number of success cases increase as we allow more error in mainbeam direction, minor-lobe level and beamwidth.For case VI, 100% success can be achieved when we allow error in mainbeam direction of 4 • , minor lobe level of 4 • and beamwidth of 5 • .For this case, the maximum weighting coefficients can be decreased from 54 dB to 32 dB.Table 3 shows 16 values of weighting coefficients comparing between IDFT method and proposed refinement method when the desired mainbeam direction is given at 45 • .For the refinement method, we set w 0 = 18 dB, Δφ = 0 • , ΔMLL = 3 dB and Δδ = 0 • .As we can see in Table 3, we can decrease maximum weighting coefficient from 22.78 dB to 18 dB which mean the range of attenuation or amplification decreases.
Figure 6 shows beamforming performance for case shown in Table 3 at frequency 2.15 GHz comparing between using IDFT method and proposed refinement method.The obtained result shows that after the refinement method is done, the mainbeam direction is still pointed to 45 • , and there is no error in beamwidth.Also, the error in average minor lobe level is within 3 dB.Figure 8: Radiation pattern at 45 • using proposed refinement method from frequencies 1.9 to 2.4 GHz.
Figures 7 and 8 show beamforming performance of the beamformer throughout the designated band using IDFT method and proposed refinement method, respectively.As we can see, the beamforming performance for both cases is similar while we can decrease the range of attenuation or amplification up to 4.78 dB (22.78-18 dB).

Practical Realizations
The full prototype employing 4 × 4 array antennas for the proposed antenna systems is designed for the operating frequencies from 1.9 to 2.4 GHz as shown in Figure 9.The prototype is a digitally controlled analog beamforming system, and it consists of three functional blocks: antenna array, weighting network, and control devices.A block diagram of smart antenna systems starts from the received incident signals at

Single Element.
A printed monopole antenna as shown in Figure 10 is chosen to be the antenna element for this paper.The reason is that this type of antenna is compact enough to be arranged to form the array having interelement spacing of half-wavelength.As appeared in Figure 10, the feed line is 25-mm in length, and ground plane is sized by 25 × 40 mm 2 .The fabricated antenna is shown in Figure 11 using single-layer FR-4-based board having its thickness of 1.6 mm.This antenna was designed to be able to operate in 1.9-2.4GHz band. Figure 12 shows the measured return loss of the antenna.As we can see, the antenna provides return loss lower than −10 dB within the designated band.Also, its radiation is omnidirectional over the designated band.

Antenna Array.
Here we present the design and development of a 4 × 4-array antenna which is capable of beam steering in the azimuth direction over an increased operational bandwidth.The design focuses on the frequency band from 1.9 to 2.4 GHz.The photograph of constructed antenna array is presented in Figure 13 when the antenna elements are spaced by half-wavelength at 3 GHz (50 mm), which is higher than the upper frequency of the assumed band (2.4 GHz).The use of this higher frequency (3 GHz) is required in the beamforming algorithm to avoid the edge effect when finding the weighting coefficients.

Weighting Network.
The weighting networks consist of four major components: digital attenuators, digital switching, phase shifters, and power combiners.Each component operates very well over wide range of frequency covering frequencies from 1.9 to 2.4 GHz.The detail of each component is shown as follows.
(1) The utilized attenuator is a 50-Ω RF digital step attenuator offering an attenuation range up to 31.5 dB with 0.5 dB step controlled by 6 bit serial interface.The operating frequency covers from DC to 2.4 GHz.
(2) The digital switching is 50-Ω high isolate SPDT RF switch designed for wireless application, covering a broad frequency range from DC up to 3 GHz with low insertion loss.
(3) Phase shifters using FR-4 substrate were designed from CST microwave studio.These components are utilized to 180 • -shift the phase of signal in order to produce the minus sign for some weighting coefficients.
(4) The power combiners are separated into two different types which are 2 : 1 operating in frequency range from 1.6 to 3.3 GHz and 16 : 1 operating for frequencies from 1.8 to 2.6 GHz.

Control Devices.
Twelve ATMEGA328 microprocessors are utilized to digital attenuator and digital switching.The reason to choose this type of control device is that it is low of cost but be able to meet the requirement in this paper.

Spatial Beamformer.
The prototype of 4 × 4 spatial beamformer is constructed to test the validation of the proposed method as shown in Figure 14.The power supply is located on the bottom of the shelf.This supports needed power to microcontrollers, digital attenuators, and digital switches.The control devices as mentioned in last section are located at the 2nd shelf followed by the weighting network which is placed on the 3rd shelf as we can see in Figure 14.
The top of the shelf is where we put the 4 × 4 antenna array as mentioned in last section.Also, the power combiners are placed on the top shelf and pass the output signal to the computer.Figure 18: Measured radiation pattern using refinement method for frequencies from 1.9 to 2.4 GHz, when the main beam direction is 7 • .Figure 20: Simulated radiation pattern using refinement method for frequencies from 1.9 to 2.4 GHz, when the main beam direction is −34 • .

Experimental Results
1.9 to 2.4 GHz.All the process starting from giving desired mainbeam direction performs automatically.However to show some example, the desired directions to be chosen are 7 directions: 2   last section.The performance is compared between employing IDFT method and proposed refinement method.There are 3 cases for this section as follows.Please note that the operating frequencies are assumed from 1.9 to 2.4 GHz.
Case A. the beamforming performance using the proposed refinement method is tested.In this case, the mainbeam direction is chosen at 7 Case B. For this case, the mainbeam is supposed to be pointed to −34 • .Also the parameters for refinement are as follows: w 0 = 4 dB, Δφ = 2 • , ΔMLL = −15 dB, and Δδ = 0 • .After the refinement process is done, we have found that the range of attenuation or amplification can be reduced from 30 dB to 4 dB.This is considerable practical.
Figure 19 shows the beamforming performance obtained from simulation and measurement for Case B at frequency 2.15 GHz.As we can see, after refinement method is done, the error in mainbeam direction is 2 • , and the error in minor lobe level is within −15 dB.Also, there is no error in beamwidth.This error occurred because we allow some error in the refinement method.Comparing with the results obtained from Case A, we can obtain higher range reduction of weighting coefficients while expecting some error in mainbeam direction, minor-lobe level and beamwidth.
Figures 20 and 21 reveal the beamforming performance for Case B throughout the designated band, both from simulation and measurement, respectively.We can observe more error in mainbeam direction, average minor-lobe level and beamwidth compared with the ones presented in Figures 17 and 18.As expected, the ones obtained from measurement are likely to have higher minor-lobe level.Case C. For this case, the beamforming performance using proposed refinement method comparative with IDFT method is pointed out.The parameters for Case C are as follows: w 0 = 4 dB, Δφ = 2 • , ΔMLL = −11 dB, and Δδ = 0 • .Please note that the mainbeam direction for this case is changed to 2 • .As a result, this case allows us to reduce the range of weighting coefficients from 8 dB to 4 dB.
Figure 22 shows the beamforming performance at 2.15 GHz using IDFT method for both simulation and measurement.Figures 23 and 24 show radiation pattern throughout the designated band, 1.9 to 2.4 GHz, obtained from simulation and measurement, respectively.As we can see, its beamforming behavior looks stable throughout the band.However, the ones obtained from measurement have slightly higher minor lobe level.
Figures 25 to 27 show the performance similar to the ones from Figures 22 to 24 but the proposed refinement method is taken place after we obtain the weighting coefficients from IDFT method.Its beamforming at 2.15 GHz is shown in Figure 25 while Figures 26 and 27 show it performance throughout the designated band via simulation and measurement, respectively.As we can see, the beamforming behavior looks similar to the cases using IDFT method with narrower range of weighting coefficients.

Conclusion
This paper has presented the method to reduce the range of weighting coefficients for fully spatial beamformer, socalled refinement method.This is considerably useful as the weighting coefficients for this kind of beamformer are real valued; therefore, narrower range of weighting coefficients is advantage.As a result, the choice of attenuators or amplifiers becomes practical.This refinement method is taken place after weighting coefficients are calculated using IDFT method.This proposed method can be succeeded in every mainbeam direction when error in some parameters is allowed.Those parameters are mainbeam direction, average minor lobe level, and beamwidth.The more errors are allowed, the more success cases can be achieved.A full prototype of fully spatial beamformer was constructed to test the validation of the proposed refinement method.The results in terms of radiation pattern obtained from simulation and measurement are compared.The obtained results reveal that we can reduce the range of attenuation or amplification using the proposed refinement method while some significant radiation behaviors are remained.

Figure 1 :Figure 2 :
Figure 1: Layout of array configuration for fully spatial beamformer.

Figure 3 :
Figure 3: Simulated radiation pattern of 4 × 4 beamformer for frequencies from 1.9 to 2.4 GHz, when the main beam direction is 45 • .

Figure 4 :Figure 5 :
Figure 4: Simulated radiation pattern of 4 × 4 beamformer plotted at 2.15 GHz, when the main beam direction is 45 • , and the objective function is Chebyshev.

Figure 6 :
Figure 6: The comparison of radiation patterns at 45 • using IDFT method and Refinement method at frequency 2.15 GHz.

Figure 9 :
Figure 9: Diagram of smart antenna system.

Figure 10 :
Figure 10: Printed monopole antenna structure at (a) front view and (b) back view.

Figure 15 :Figure 16 :
Figure 15: Testing the constructed prototype in an anechoic chamber.

Figure 17 :
Figure17: Simulated radiation pattern using refinement method for frequencies from 1.9 to 2.4 GHz, when the main beam direction is 7 • .

5. 1 .Figure 19 :
Figure19: Simulated and measured radiation pattern using refinement method for frequency 2.15 GHz, when the main beam direction is 34 • .

Figure 21 :Figure 22 :
Figure21: Measured radiation pattern using refinement method for frequencies from 1.9 to 2.4 GHz, when the main beam direction is −34 • .

Table 1 :
Weighting coefficients when the required radiation pattern is created using Sinc and Chebyshev functions.

Table 2 :
Parameters given in simulation for refinement algorithm of 6 cases.

Table 3 :
The comparison value of attenuator factors between IDFT method and refinement method.
• ,5.2.Experimental Results.In this section we validate the beamforming performance of the proposed refinement method through the prototype measurement described in