Polarized Uniform Linear Array System: Beam Radiation Pattern, Beamforming Diversity Order, and Channel Capacity

There have been many studies regarding antenna polarization; however, there have been few publications on the analysis of the channel capacity for polarized antenna systems using the beamforming technique. According to Chung et al., the channel capacity is determined by the density of scatterers and the transmission power, which is obtained based on the assumption that scatterers are uniformly distributed on a 3D spherical scattering model. However, it contradicts the practical scenario, where scatterers may not be uniformly distributed under outdoor environment, and lacks the consideration of fading channel gain. In this study, we derive the channel capacity of polarized uniform linear array (PULA) systems using the beamforming technique in a practical scattering environment. The results show that, for PULA systems, the channel capacity, which is boosted by beamforming diversity, can be determined using the channel gain, beam radiation pattern, and beamforming diversity order (BDO), where the BDO is dependent on the antenna characteristics and array configurations.


Introduction
Recently, there has been a gradual increase in the demand for polarized antenna systems, especially for 5G networks.This is mainly because antenna polarization is an important parameter employed for the design of space-limited wireless devices [1,2].Several techniques such as space-time diversity, multiplexing, and array processing can be exploited in order to boost the system throughput for polarized antennas.Thus, many published research articles as well as ongoing projects have focused on antenna polarization [3][4][5].However, there have been few reports on the analysis of the channel capacity for polarized antenna systems using the beamforming technique.In [6], the channel capacity is known to be determined by the density of scatterers and the transmission power instead of other system parameters.This is based on the compact uniform linear array (ULA) and assumes that scatterers are uniformly distributed on a 3D spherical scattering model.However, in practical scenarios, the scatterers may not be uniformly distributed, especially in outdoor environments.In addition, [6] does not consider the fading channel, which is important when obtaining the channel capacity in [7,8].In this study, we obtain the channel capacity of polarized ULA (PULA) systems using the beamforming technique in a practical scattering environment.Our results indicate that, for PULA systems, the channel capacity, which is boosted by beamforming diversity, is determined mainly by the channel gain, beam radiation pattern, and beamforming diversity order (BDO).
In this paper, the beam radiation pattern (including the beamwidth, the number of lobes, and the radiation gain of the beam) of conventional ULA systems is described based on the antenna characteristics and array configurations (i.e., dipole antennas with different lengths, the number of array elements, and the array element spacing).And we characterize the beam radiation pattern for PULA systems in both 2D and 3D forms, where the dipole antenna-based PULA is shown to have a higher beam radiation gain compared with isotropic antenna-based systems.The BDO, which is defined with respect to the number of beam lobes generated by the transceiver, is obtained under different transmission scenarios by considering the line-of-sight (LOS) component Boresight direction and scattering environments.In addition, it is dependent on the antenna characteristics and array configurations and has been shown to be linearly proportional to the number of array elements and array element spacing.
The remainder of this paper is organized as follows.In Section 2, we analyze the conventional ULA, while the beam radiation patterns of the PULA in both 2D and 3D forms are presented in Section 3. In Section 4, we introduce the BDO concept, and in Section 5 we discuss the channel capacity of PULA systems.Finally, we conclude the paper in Section 6.

Feature Analysis of Conventional ULA
2.1.Array Factor.Figure 1 illustrates a conventional ULA, where  0 denotes the signal incidence angle shifted from the boresight direction (the antenna boresight is the axis vertical to the orientation of the array alignment [9]),  represents the array element spacing, and  is the number of array elements.The array factor (), which represents the far-field radiation of an isotropic ULA system [10], can be expressed in terms of the inner product of the weighting vector and array propagation vector as follows: The weighting vector is given as and the array propagation vector is given as where  denotes the signal wavelength.If we use the weighting vector to adjust a beam with a specific angle (shifted from the boresight direction) of  0 and use the array spacing factor  to describe the array configuration, that is,  = /, the array factor can be further expressed as Based on the derivation in Appendix A, the normalized () is obtained as Because −2 ≤ sin  − sin  0 ≤ 2 and  ≥ 2 are often used to avoid the array aliasing problem [11], () then approximates to the sinc function as (5) 2.2.Off-Boresight Angle.As discussed in [9], the beamwidth is increased significantly when the beam steers to an angle far off the boresight direction, such as the case when  0 reaches 90 ∘ .Note that the beamwidth generally refers to the width of the main-lobe of a beam.
This illustrates that the beamwidth is inversely proportional to the number of array elements  and array element spacing  but is directly proportional to the off-boresight angle  0 .
Figure 2 shows HPBW examples of isotropic ULA systems that were obtained by varying  0 for three different angles, that is,  0 = 0 ∘ , 45 ∘ , and 90 ∘ .The figure shows that the HPBW extends about four times when  0 changes from 0 ∘ to 90 ∘ .Note that, in this illustration, we set  = 8 and  = 2 and the boresight direction along the horizontal axis.To focus the beamwidth when there is a large off-boresight angle effect, [9] proposes a sensor delay line-based array system, which can effectively prevent the larger off-boresight angle effect by increasing the dimensions of the array system.

Number of Lobes for a
Beam.We previously discussed the beamwidth, which is affected by the number of array elements, the array element spacing, and the off-boresight angle.In this section, we analyze the relationship between these factors and the number of lobes present in a beam.In the case where relays from different directions are used to contribute space diversity in an effort to boost the capacity of a ULA system [12], all of the lobes of a beam need to be considered as the BDO, which is eventually used to derive the channel capacity.
According to (4), there are zero values at sin  = ±( ⋅ /)+sin  0 , where  = 0, . . ., −1 and  ∈ (0, ).Therefore, the lobe-widths of the ULA beam are given by Because each lobe has two zero points and there are total 2 −1 zero points at  ∈ (0, ), we can find (2 −1 − 1) lobes.Due to the symmetrical beam pattern of ULAs, the overall number of lobes is given as (4 −1 − 2), which is unrelated to the off-boresight angle  0 .Figure 3 depicts the number of lobes when  is incremented with  = 2, where the number of lobes can be easily observed as 6, 14, 30, and 62 when  equals 4, 8, 16, and 32, respectively.Based on the above analysis, we can make the following proposition.

Proposition 1. The beamwidth and number of lobes of a beam are relevant to the array configuration. The beamwidth is inversely proportional to the number of array elements and the array element spacing but is directly proportional to the offboresight angle. Meanwhile, the number of lobes of a beam is proportional to the number of array elements and the array
element spacing but has no relationship with the off-boresight angle.

Beam Radiation Pattern of PULA
In this section, we discuss the beam pattern associated with the antenna characteristics, such as the use of dipole antennas, where the beamwidth and beam lobe gain are different from the case that uses the isotropic ULA.

Cross-Array Beam Generation by PULA Configuration.
The dipole antenna or doublet is the most widely used class of antenna, and the radiation pattern analysis throughout this paper is based mainly on the dipole antenna.Figure 4 demonstrates a dipole antenna with a length of .(, , , ) denotes the far-field radiation pattern of the dipole, where  is the angle of the dipole shifted from the -axis and  is the factor used to determine the dipole length; that is,  = .The half-wavelength dipole is defined when  = 1/2.Figure 5 shows our proposed orthogonal triple-polarized dipole ULA, where an array and branch (&) multiple antenna configuration is depicted.In other words, the PULA system is composed of several paralleled array elements, and there are three orthogonal colocated branches/dipoles at each array element.Three branch sets of    1 ,    2 , and    3 are implemented along the -axis, -axis, and -axis, respectively.
Because three colocated branches at each array element are fixed with zero spacing, as shown in Figure 5, which makes the beamwidth scale up to 360 ∘ according to Proposition 1, the beams need to be generated by using corresponding cross-array branches rather than the colocated branches at each array element.Therefore, we can obtain three orthogonal beams generated by a PULA as follows: (i) The beam varied in the - plane is generated by the set of branches of    1 , where  = 1, 2, . . ., .
(ii) The beam varied in the - plane is generated by the set of branches of    2 .
(iii) The beam varied in the - plane is generated by the set of branches of    3 .

PULA Beam Radiation
Pattern in 2D Form.First, the PULA radiation pattern is analyzed in 2D form, where only two orthogonal colocated branches (vertical and horizontal) are implemented at each array element.By fixing the vertical branch along with the vertical axis when  = 0, we have the following: the vertical dipole radiation pattern  V (, ) in 2D form (where  = 0) is presented as Figure 6 depicts the dipole radiation pattern of the vertical branch with  different values.We observe that the dipole does not radiate in the longitudinal direction of the antenna structure and consequently has a higher gain compared with that of an isotropic antenna.In addition, the radiation pattern has a higher gain and narrower width as  increased, where the case of  = 1 achieves the lowest angle of radiation and the highest horizontal range.Note that when  > 1, the radiation pattern of the dipole does not follow the current shape in Figure 6 [13].
The radiation pattern of the colocated dipole can be easily obtained based on (8) as where it is shifted by  degrees and the horizontal dipole pattern is obtained as  = 90 ∘ .The PULA 2D beam radiation pattern is expressed by multiplying the array factor and the radiation pattern of a single polarized antenna.Let  V (, ,  0,V ) and  ℎ (, ,  0,ℎ ) denote the beam radiation patterns generated by the sets of vertical dipoles and horizontal dipoles, respectively.We have where  0,V and  0,ℎ are the off-boresight angles corresponding to the sets of vertical branches and horizontal branches, respectively.Figure 7 is drawn when  = 4 and  = 2.We observe that the beam radiation patterns generated by the sets of vertical dipoles and horizontal dipoles have a 90 ∘ rotation.The beamforming gain can be extended by increasing the dipole length.However, the trade-off regarding the antenna size needs to be considered in the antenna design.The 2D beam radiation pattern for an isotropic ULA is also given in Figure 7 for comparison purposes, where it has a larger sidelobe and smaller main-lobe pattern compared with the dipole-based array system [14].This is mainly because the dipole does not radiate in the longitudinal direction of the antenna structure and maintains a higher radiation gain compared with that of the isotropic antenna.

PULA Beam Radiation
Pattern in 3D Form.In [3], the 3D dipole radiation pattern is found to have a donut shape.In addition, the 3D radiation pattern of the dual-polarized halfwavelength dipole was investigated in [4,15] with different inclination angles of , for example, 0 ∘ , 45 ∘ , and 90 ∘ .In this paper, we further obtain the 3D radiation pattern for a triple-polarized system, where three colocated dipoles ( 1 ,  2 , and  3 ) are deployed along the -axis, -axis, and -axis, respectively.The formulae for the radiation pattern derived based on (31) in [4] for these three colocated dipoles are given as Based on (11) the radiation patterns of triple-polarized half-wavelength dipoles are plotted in Figure 8, where the three patterns affect each other and the donut radiation shape of the dipole antenna is reshaped due to the presence of colocated dipoles.

Beamforming Diversity Order (BDO)
In the relevant studies, for example, [16,17], the interference reduction by beamforming is studied, where the sidelobe radiation of a beam is suppressed via the leaky wave antenna and particle swarm optimization in order to decrease the sidelobe interferences caused at neighbor cells.However, in this paper, we focus on the theoretical analysis under a perfect scenario without any interference.The radiation of sidelobes can also contribute to the SNR that enhances the International Journal of Antennas and Propagation capacity, because the BDO is increased with the additional number of sidelobes.Figure 10 describes a radio transmission scenario when beamforming technology is applied at both the transmitter (Tx) and receiver (Rx).In this case, we define the signal-transmission cutting surface (STCS), which is the surface holds the signal from the Tx to Rx. to Figure 10(a), three are simultaneously held by the at Tx (STCS Tx ), and three lobes are simultaneously held by the STCS at the Rx (STCS Rx ).And we define the BDO as the multiplication of the number of lobes cut by STCS Tx and STCS Rx .For example, in the case of Figure 10(a), the BDO is 9. On the other hand, the BDO may decrease because of the large off-boresight angle effect, as illustrated in Figure 10(b); for example, the BDO is one where only one lobe is held by the STCS Tx and STCS Rx .
Practically, the scatterers are often distributed into clusters [18,19].A group of local scatterers within a cluster have relatively close coordinates on the sphere representing the incidence signal scattered by the roofs, walls, trees, windows, doors, and so on from a correlated direction.Meanwhile, global scatterers among different clusters have independent coordinates, which denote the incidence signal scattered from uncorrelated directions.Given a cluster of local scatterers, for the th scatterer with a polar coordinate of   (  ,   ), the signal response after scattering is given as where  Tx (  ) is the signal propagation delay from the Tx to   ,  Rx (  ) is the signal propagation delay from   to the Rx, and   denotes the signal carrier frequency.  represents the channel gain and exp(⋅) gives the signal phase rotation after scattering.Because the local scatterers in this cluster are very close to each other, we have where  0 represents the coordinates of the center of gravity for the cluster.Then, ℎ 0 () can be used to represent the scattered signal as follows: When the scatterers in different clusters are not close to each other, ( Tx ( 0, ) +  Rx ( 0, )) is independent of ( Tx ( 0,+1 ) +  Rx ( 0,+1 )) for  = 1, . . ., , where  0, is the coordinate of the center of gravity for the th cluster.Then, the received signal power can be denoted as The fully scattered environment refers to the scattered signal coming from all directions; that is,  min = − and  max = , and dense urban and indoor scenarios are usually studied under such an environment.The non-fully-scattered environment is realized when  min > − and  max < .Most rural and remote scenarios, where wide open spaces and flat terrain fulfill the LOS conditions, are studied in this environment [20].
We assume the scenario where the Tx, which steers its main-lobe toward Rx, is located far from Rx. Consequently, the BDO at the Tx side can be assumed to be one, and the overall BDO of the system is determined only by the Rx.In addition, suppose several clusters of local scatterers are located on a sphere that encircles the Rx, as illustrated in Figure 11; we can summarize the system BDO as given by Table 1.From Table 1, we note that BDO equals the number of lobes generated by the Rx in a fully scattered scenario.In the case where the LOS component exists in a non-fullyscattered environment, the main-lobe of the Tx covers half of the surface of the sphere around the Rx.Therefore, the BDO is varied from 2 −1 − 1 to 4 −1 − 2. For the case of Moreover, the fact that each lobe has two zero points was overlooked in [3].

Channel Capacity by PULA with Beamforming
For a fully scattered environment under the assumption that there is no interference, [6] derives the channel capacity for an array system as follows: where  =   /(8(2/) 2  2  ). refers to the density of the scatterers located on the 3D spherical model,   is the number of scatterers,   is the radius of the scattering sphere,  is the transmission power, and  0 indicates the additive noise power density.From (17), we observe that the channel capacity is only related to the scatterer density and transmission power as opposed to other system parameters.Equation (7) concludes in two ways.At first, it describes that the capacity is regardless of the system bandwidth, where a high throughput can be achieved even by a negligibly small bandwidth.And second, it shows that the capacity is irrespective of the antenna array configuration.However, the conclusions in [6] were obtained using a ULA compact antenna with  ≫ 2, and based on Table 1 the diversity order thus decreases significantly.In [6], we also see that ( 17) is employed when the scatterers are uniformly distributed on the sphere.This is not the case practically, and the scatterers' distribution may follow the example shown in Figure 11.In addition, the fading channel gain, which is usually used to derive the channel capacity, was not considered by (7).Therefore, we provide a new method to obtain the channel capacity for PULA systems, and this method is based mainly on the fading channel gain, beam radiation pattern, and BDO.
In [7,8], the channel capacity with Rx diversity is given by where ℎ denotes a gain of fading channel.Here, we define the beam lobe signal to noise power ratio (LSNR) based on the beam radiation pattern as where  denotes the lobe index of the beam, θ (can be found in Appendix B) holds the condition for (, ) having peak values, and   refers to the added noise power at the th lobe.Based on (19), the effect of the off-boresight angle can be disregarded according to Proposition 2.

Proposition 2. The off-boresight angle does not affect the lobe power but affects the lobe-width (proved by Appendix B).
The channel capacity for a PULA system under the Rayleigh fading channel can be given as where  = 6 is the degree-of-freedom of the triplepolarized antenna in [3].ℎ  denotes a gain of fading channel corresponding to the th lobe.Then, the capacity upper bound is found by Jensen's inequality for concave functions as follows: Rayleigh ≤ 2 where the first term contains the capacity contributed by the main-lobe that steers the LOS component ( is the Rician -factor).The second term represents the capacity obtained by this main-lobe for the NLOS component, while the third term indicates the capacity from other lobes of the NLOS component.Note that (21) and ( 22) apply to the scenario involving a fully scattered environment.In the case of a nonfully-scattered environment, the PULA capacity is decreased due to the decrement of BDO.
Figure 12 gives the normalized channel capacity for PULA with beamforming.From Figure 12, we observe that, by extending the dipole length (e.g.,  = 1), we can effectively increase the capacity of the PULA system.This is because the radiation gain increases as  increases, where the case of  = 1 achieves the highest gain, as discussed in Section 3. The channel capacity is also increased when more array elements are implemented (e.g.,  = 8), and this is mainly due to the increment of the BDO.Compared with the capacity increment obtained by extending the dipole length, the capacity increment due to the BDO is small because the array factor is normalized by the number of array elements in (4).Finally, Figure 12 depicts the channel capacity under different channel models, where the AWGN has a higher capacity than in the case of fading channels.In addition, the Rician fading channel with  = 9 dB can achieve a higher channel capacity compared with the Rayleigh fading channel due to the presence of the LOS component.

Conclusions
In this paper, we derive the channel capacity for PULA systems using the beamforming technique.Our results show that the channel capacity is determined mainly by three system indexes, that is, the fading channel gain, beam radiation pattern, and BDO, where the BDO is dependent on the antenna characteristics and array configurations.In contrast to conventional works, we also consider the fading channel in a practical scattering environment, and the numerical results confirm the validity of the channel capacity obtained in this paper.

Figure 1 :
Figure 1: Configuration of the conventional ULA.

Figure 2 :
Figure 2: Illustration of the off-boresight angle effect.

Figure 3 :Figure 4 :
Figure 3: The number of beam lobes varies with different  values.

Figure 7 :
Figure 7: Beam radiation patterns of the PULA in 2D form.

Figure 8 :
Figure 8: Radiation patterns of the triple-polarized half-wavelength dipole in 3D form.

Figure 9 :
Figure 9: Beam radiation patterns of PULA in 3D form.

Figure 12 :
Figure 12: The normalized channel capacity for PULA with beamforming.

Table 1 :
BDO values depending on scenarios.