DoA Prediction Based Beamforming with Low Training Overhead for Highly-Mobile UAV Communication with Cellular Networks

: In supporting communications with unmanned aerial vehicles (UAVs) as aerial user equipments (aUEs) in cellular systems, the current beamforming schemes based on channel state estimation are facing severe challenges from the pilot contamination effect, especially in 5G and future networks where the cell size becomes small and the user density is high. Beamforming schemes based on signal direction of arrival (DoA) are regarded as a highly promising alternative to solve this problem. However, to achieve optimal performance for DoA-based beamforming, the error to DoA estimation during pilot signal intervals, caused by the high mobility of UAVs, must be addressed. In the meantime, the training overheads of traditional DoA estimation algorithms must be reduced to save the bandwidth for data communication. This paper investigates uplink beamforming performance enhancement based on signal DoA estimation to support UAV-cellular network communication. We propose a novel DoA estimation algorithm to predict angle variations during the intervals, which achieves high precision even when UAVs are at high mobility. The prediction process requires no pilot signals and enables timely adjustment of the steering vector when calculating the beamforming weight vector. The proposed algorithm contributes to the realisation of a beamforming scheme with real-time steering vector updates, which simultaneously maintains high beamforming gains and low training overheads. Simulation results show that, compared with the conventional DoA-based beamforming scheme, the proposed method yields more accurate DoA estimation output and higher gains. Furthermore, simulation experiments also suggests that applying the proposed scheme can reduce up to 100 pilot signal transmissions per second.


Introduction
With the rapid development of UAVs, these devices are revolutionizing a large quantity of civil application fields, including industrial construction, products inspection, safety monitoring, agriculture, rescue and search, bird-eye viewing and so on [1,2]. Some common functions of UAVs in these applications are photographing and video recording from the air, and transfering high-resolution photos or videos requires a robust and fast communication link between the UAV and the ground server. It is thus useful to incorporate UAVs into the cellular system such as aUEs, due to the large  The methods introduced above cannot maintain highly accurate estimations for fast-moving UAVs without incurring a high system training overhead, because they focus on a single time's connection via beamforming, without giving enough attention to the realistic mechanism of the cellular system. In these schemes, the typical solution to the UAV's high mobility problem is to shorten the pilot signal interval (t p2 − t p1 ), so that the DoA information is updated frequently to avoid large decrease in gain. However, this solution also requires a higher number of pilot signals to be transmitted for each user in the network. The work in [14,20] have reduced the required number of pilot signals in each beamforming training, however, these methods still heavily rely on training sequences because each DoA estimation at least requires one pilot signal. In future cellular networks, the number of users is expected to be increased by a factor on the order of tens [22], and if each individual user requires a large number of pilot signals, the total training overhead will be greatly increased and occupy an excessive amount of the bandwidth originally allocated for data transmission. Furthermore, most real-world deployment of the cellular system will set a minimum interval for transmitting pilot signals, which restricts the frequency for DoA estimation. In [23,24], position-aided beamforming methods are proposed for high-speed train communications. Since the track of the train is known to the BS processor, beamforming can be optimised by casting the highest gain to the train's antenna at any moment, without the assistance of pilot signals. However, these methods are not applicable in UAV communications, since most UAVs do not have pre-defined trajectories. Even if trajectories of UAVs are known in some application fields, such tracks are far less precise than the ground rail way case.
In this paper, we propose an algorithm for evaluating the DoA of a UAV by utilizing position information obtained from the UAV's flight status information (FSI). The algorithm incorporates a prediction method for tracking the angle variations during intervals between the reception of FSI update messages. Since the prediction process does not require the transmission of pilot signals, it can be executed frequently without increasing the training overhead. In this way, a steady and accurate output of DoA estimations is obtained. The proposed algorithm enables the realization of a novel beamforming scheme. As illustrated in the lower half of Figure 1, instead of using the DoA calculated at t p1 , this scheme allows the weight vector to be obtained using up-to-date DoA values generated via prediction at times between t p1 and t p2 . Thus, high gains can be maintained for highly-mobile UAVs. Furthermore, instead of requiring a shorter (t p2 − t p1 ), the proposed scheme can simply generate more predictions to maintain accurate DoA information, resulting in low training overhead.
Simulation experiments carried out to verify the proposed method are presented. Compared with the conventional DoA-based beamforming scheme using an improved variant of the MUSIC algorithm [17], the proposed method yields more accurate DoA estimations and higher gains. The performance margin is especially remarkable when the UAV is flying fast or the DoA estimation interval is large. Besides, simulation experiments also suggests that applying the proposed scheme can reduce up to 100 pilot signal transmissions per second.
The rest of the paper is organized as follows: Section 2 describes the sources of DoA errors and their impact on beamforming methods. Section 3 introduces the proposed algorithm and the method for the real-time steering vector update beamforming scheme. Simulation experiments and results are given in Section 4 to verify their performance. Section 5 concludes this paper.
Notations: We employ boldface for vectors and bold capital lettes for matrices; the superscript (·) T denotes the transpose of a vector or matrix; the superscript (·) H denotes the Hermitian transpose of a vector or matrix; the superscript (·) −1 denotes the inverse of a matrix; I is the identity matrix; exp(x) denotes exponentia; | · | denotes absolute value; the superscript (·) • denotes angle in degrees. The acronym list is shown in Table 1.

Problem Formulation
We consider a UAV-cellular communication scenario in which a UAV, as an aUE, transfers data over the uplink, while other UAVs simultaneously communicating with their serving BSs introduce interference with the target received uplink signal. The UAV's altitude is assumed to be above 40 m, in which case, according to 3GPP Release 15, it has a LoS probability of 1 with the BS. To improve the received signal-to-interference-plus-noise ratio (SINR), the BS applies beamforming based on DoA estimation by first estimating the angle of arrival of the target signal and then executing beamforming.
Here, we illustrate the sources of DoA estimation errors due to the UAV's mobility and the influence of these errors on the performance of beamforming algorithms. A brief discussion on such error has already been provided in Section 1. Here, we present a general introduction to DoA estimation errors. In a typical UAV-cellular communication scenario, multiple factors, including the angular resolution of the algorithms, the signal-to-noise ratio (SNR) and the angular spread effect, can affect the error range of DoA estimation algorithms. Moreover, additional DoA errors are caused by the UAV's mobility during the time interval between two DOA estimations. Suppose that two consecutive DoA estimations are performed at time instants t p n−1 and t p n , with an interval ∆t = t p n − t p n−1 . During ∆t, the actual DoA of the incoming signal continuously varies, while the estimated DoA remains the same. Thus, at any time instant t m , t p n−1 ≤ t m ≤ t p n , a DoA estimation error exists, which is equal to the difference between the estimated angle and the actual angle: where e(t m ) denotes the estimation error at t m and η(t m ) and γ e are the corresponding actual and estimated DoA values, respectively. In addition to a longer DoA estimation interval, a higher UAV speed or a reduced cell radius in the 5G system will also cause this DoA estimation error to increase [20,25]. The gain achieved with classical nonrobust beamforming decreases dramatically with even a minor DoA error [26,27]. Although robust beamforming methods have been developed to suppress such sensitivity, the optimal gains of all these algorithms degrade to various extents as the DoA error increases. We take the widely studied diagonal loading (DL) algorithm for robust beamforming as an example to illustrate this effect. In this algorithm, the beamforming weight vector, w, is obtained as follows: where I L is an identity matrix; s(θ) denotes the steering vector of the target signal and θ is the estimated DoA value; R is an L × L array correlation matrix, the elements of which represent the correlations between the received signals of the array antennas. The coefficient σ, known as the DL factor, modifies R by enhancing the effect of noise, forcing the beamformer to mitigate noise and thus weakening its ability to suppress interference. In this way, the beamformer is made less sensitive to DoA errors. Over a span of many years and up through recent times, a large number of algorithms have been proposed for mapping the DL factor with a certain DoA error range to achieve optimal performance [26,[28][29][30]. However, simply finding the most suitable DL factor will not result in optimal gains since the DoA error range places an upper limit on the achievable gain. A plot of the SINR output of a DL beamformer versus the DL factor is shown in Figure 3, where the input SNR is 10 dB and the number of antenna elements is 8. The blue, red and orange curves represent the outputs for DoA estimation errors of 1 • , 3 • , and 5 • , respectively. In this figure, 3.93 dB and 3.54 dB differences in the maximum gains are observed between the blue and red curves and the red and orange curves, respectively, suggesting that a large DoA error inevitably leads to a low gain. Thus, precise DoA estimation is crucial for proper beamforming performance.

The Proposed DoA Estimation Algorithm and Beamforming Scheme
The proposed algorithm derives real-time DoA information of the UAV from its FSI, and then the calculated DoA values are supplied to the beamformer and used to obtain the optimum beamforming weight vector. Considering the size and power consumption limitation of most commercial UAVs, we assume the target UAV carries a single omni-directional antenna.

UAV FSI
Positions are sourced from the UAV FSI, which incorporates the UAV navigation data. By using a Kalman filter to integrate Global Positioning System (GPS) and on-board sensors' motion information, the UAV's navigator maintains more accurate positions and motion vectors (velocity and acceleration) than simply using GPS systems [31,32]. UAV's FSI is regulated to be transmitted to the ground controllers periodically, for safety concerns. This telemetry message is exchanged at a minimal frequency of 4-5 Hz [1]. Although the BS could command the UAV to report more frequently, we will demonstrate in Section 4 that this is not an actual need. FSI is a reliable information resource-as specified in the 3GPP standard, it has high quality of service (QoS) priority (up to 10 −3 Packet Error Loss Rate) [33]. The first step to calculate the signal DoA consists of extracting flight status data from the telemetry message, which can be realized by creating an interface between the coordination information stack and the beamforming processor.
To summarize, utilizing the FSI has several benefits. First, UAV's position is deducted from processing UAV's on-board GPS and the inertial measurement unit (IMU) data, which obtains higher accuracy than using the GPS sensor alone, and reduces the computational load on the BS end. Second, FSI is regulated to be transmitted for any type of UAVs, so that the proposed algorithm is applicable for all civil applications. Third, by using FSI, which already exists in the Command&Control Links between UAV and cellular networks, we reduce the use of training sequences and thus reduce the training overhead.

DoA Estimation and Prediction
The DoA of the incoming signal is defined as the angle between the signal's propagation direction and the receiver's antenna array. To calculate DoA, the position vector from FSI needs to be first transformed to a common reference which incorporates the BS's position vector. The universal local north-east-down (NED) frame is used to this end, where the z-axis is perpendicular to the earth surface, and x and y axes are both tangent to the earth surface. A GPS position is typically expressed in Geodetic frame as ρ = (ω, ϕ, h) T , where h is the height above the earth surface; ω and ϕ are longitude and latitude, respectively. Set the original UAV and BS's GPS positions as ρ ug and ρ bg , respectively. Therefore, To perform the transformation, their equivalent positions in Earth-Centered-Earth-Fixed (ECEF) frame, p ue and p be , are calculated first [34], where p ue and p be are the UAV and the BS's position vectors, respectively. r e is the earth radius, and e is a constant. Then, the UAV's position in our defined local NED frame, p ul , is calculated as: where T en denotes the transformation matrix from the ECEF coordinates to the local NED coordinates, After incorporating the UAV's position, the DoA estimation is obtained as follows: where (x bl , y bl , z bl ) is the BS's position. As the BS is on the origin, Equation (7) can be simplified as follows: In the prediction algorithm, position variations based on estimated UAV's movements are first calculated, and then transformed into the new estimated DoA angles. Assume a total number of K 0 predictions of the DoA values are performed between two consecutive DoA estimations via FSI, and the n-th calculated UAV's position is denoted as p n = (α n , β n , λ n ) T , 0 ≤ n ≤ K 0 . Then, denote the n-th calculated UAV's velocity and the acceleration as v n = (u n , v n , w n ) T , 0 ≤ n ≤ K 0 , and µ = (µ x , µ y , µ z ) T , respectively. The status transition matrix is then given as follows: where τ is the time step of the prediction process. Here, v 0 and µ are deducted from the original velocity v b and acceleration vectors µ b which are read from FSI in UAV's vehicle carried NED frame, where T bn is the transformation matrix from the vehicle carried NED frame to the local NED frame, where, φ, δ and ψ denote the Yaw, Pitch, and Roll rotational coordinates, respectively [34]. After the new position is predicted, DoA could be updated using Equation (8). Figure 4 illustrates how the prediction algorithm maintains a stable DoA estimation output. Suppose a linear array lies on the x axis and a DoA estimation is executed when the UAV is at the position pr n , with the actual DoA to be η n . By applying the proposed algorithm the BS calculates the assumed DoA γ n , using the UAV's position from the FSI message, pe n . After a certain time, the UAV reaches a new position pr n+1 , and the BS gets a new assumed position of it, pe n+1 . Now the actual and assumed DoA values become η n+1 and γ n+1 , respectively. Thus, the DoA error becomes η n+1 − γ n+1 , whereas on the estimation point it was η n − γ n . It can be seen that compared with Figure 2, the DoA error has not shown notable variations with the movement of the UAV. We use the curves c r and c e to denote the actual and the predicted tracks, respectively. Since the IMU keeps highly accurate UAV velocity and acceleration information during a short time, the position error remains in a limited range during the whole prediction process. c r can be approximated by c e with a vector: = c e + (pr n − pe n ).
Although the prediction process does not obtain a significantly more accurate position information, it keeps the DoA error free from cumulative increasing without adding to the training overhead.

Beamforming Scheme with Real-Time Steering Vector Updates
By using the proposed method, the beamforming scheme is also radically different. The traditional only updates steering vector relying on pilots, now most estimations are done without processing traning sequences, and so steerging vector updates in real-time. In the beamforming process, the steering vector s(θ) is used to adjust the direction of the main beam toward the target. Furthermore, it is derived directly from the DoA value. Thus, based on the proposed DoA estimation method, we propose a beamforming scheme that is suitable for providing service to UAVs with fast and stable uplink data transmission. First, the relation between the DoA value and the steering vector is introduced. Consider an array of L omnidirectional elements. When the DoA information is acquired, the steering vector is calculated as, where t d is the arrival time difference of the transmitted signal between array elements, given by: where r is the distance between the array elements and c is the speed of light. In our proposed beamforming scheme, beamforming is applied based on DoA estimations. This scheme consists of two stages during uplink data transmission. In Stage 1, the processor receives the FSI from the telemetry messages and calculates the DoA and the beamforming weight vector; then, it proceeds to Stage 2. In Stage 2, the processor continues to update the DoA information and the steering vector via the prediction process. When new FSI is obtained, Stage 2 terminates, and the processor returns to Stage 1. Through this process, the steering vector is updated in real-time so that the main beam can be constantly adjusted to track the target UAV, thus providing high gains for the signals of highly-mobile UAVs. The proposed scheme is summarized in Algorithm 1: 3. Calculate the steering vector using Equation (14) and then calculate the beamforming weight vector; 4. Set k = 1, then go to Stage 2; Stage 2 ; if k < K P then 1. Perform DoA variation prediction using Equations (9)-(12) and (8) ; 2. Update the steering vector using Equation (14), then update the beamforming weight vector; 3. Set k = k + 1; else Go to Stage 1; end end

Simulation Experiments
This section presents experiments conducted to verify the proposed DoA estimation algorithm and beamforming scheme. For this purpose, we present a performance comparison with a baseline beamforming scheme using the widely studied MUSIC algorithm. Here, a recently proposed improved variant of the MUSIC algorithm [17] is used, which offers higher accuracy and lower computational complexity than the conventional MUSIC algorithm. In all experiments, the UAVs and the BS are assumed to be in the defined local NDE frame. On the BS side, there is an antenna array with 16 elements lying on the z axis. The BS height is 15 m above the ground. The spacing between the array elements is half of the wavelength, which is 0.165 m in this case. Two UAV models in addition to the target UAV are include, which fly in circles and produce uplink signals that act as interference with the receiver. The speed values of the UAVs are selected in accordance with those applied in 3GPP Release 15. The maximum speed is 160 km/h, and a moderate value of 60 km/h is chosen. The target UAV's transmission power is 23 dBm, and the transmitted signal-to-interference ratios (SIRs) are 1 dB and 2 dB. The noise power is −95 dBm. The carrier frequency is 0.9 GHz, and the speed of light is assumed to be 3.0 × 10 8 m/s. The wireless air-ground channel model is built in accordance with the specifications of 3GPP Release 15 [3]. The small-scale fading parameters of the channel are determined in accordance with the third alternative model for fast fading in that document. The number of clusters contained in a signal is 15 [35]. The sampling interval of the improved MUSIC algorithm is 0.03255 × 10 −9 s, and the number of snapshots is 3000. The symbol length is 66.7 × 10 −9 s.
The receiving power Pr (in mW) of each UAV is expressed as follows: where P t and PL denote the transmission power (in dBm) and the path loss (in dB), respectively. The path loss is a function of the 3-dimensional (3-D) distance between the UAV and the BS and the UAV's altitude [3]: where d 3D and h are the 3-D distance and the UAV's height, respectively, and f c is the carrier frequency. The settings of the interfering UAVs are identical in all experiments. Throughout the rest of this section, 'UAV' refers to the target UAV unless otherwise stated. The simulation parameters are summarized in Table 2.

Performance Comparison for High and Low UAV Speed Scenarios
The purpose of this group of experiments is to compare the DoA estimation and beamforming performance of the proposed and baseline methods. Here, two scenarios are studied to consider the high and moderate mobilities of the UAV, for which the UAV's speed is set to 160 km/h and 60 km/h, respectively. Horizontally, the trajectory of the target UAV is a curve with a radius of 80 m. We set the initial distance between the UAV and the BS as 100 m. For both beamforming schemes, we adopt two widely studied beamformers, a nonrobust method and a robust method, namely, the linearly constrained minimum variance (LCMV) beamformer and the DL beamformer, respectively. Because the typical transmission frequency of FSI messages is 4-5 Hz, the pilot signal interval is set to 200 ms for the proposed algorithm. The prediction process has low complexity and thus can be executed at a high frequency. However, given the limited change in the position of the UAV within a short time duration, the prediction process is executed every 1 ms in our experiments. The improved MUSIC algorithm considered for comparison uses reference signals to perform DoA estimation; considering that the reference signals usually have a higher transmission frequency than the telemetry messages, the DoA estimation interval for the traditional beamforming scheme is set to 50 ms. The errors of both DoA estimation methods and of the beamforming output are examined every 10 ms. Each time the averaged absolute error output of 20 times DoA estimation results are obtained. The widely adopted GPS multivariate Gaussian positioning error model is applied for the proposed method. Both the vertical and horizontal errors have means of zero. Based on a review of popular commercial products [36][37][38], the standard variances of the vertical and horizontal errors are selected to be √ 3 m and √ 2 m, respectively. We compare the computational complexity of both DoA estimation algorithms using simulation experiments in Matlab. The general setting is as follows: considering the realistic execution frequency of both algorithms, the proposed algorithm is performed once and the prediction is executed once per 1ms; the improved MUSIC algorithm is performed four times. The algorithm execution time is obtained by averaging 100 independent executions. The average time duration in executing the conventional MUSIC algorithm is 0.0540 s, while the proposed method is 0.0029 s. The MUSIC algorithm is approximately 19 times higher in computation time than the proposed method. Therefore, the proposed method is of significantly lower computational complexity.
The DoA estimation errors for the high-speed and low-speed scenarios are shown in Figure 5a,b, respectively. From these figures, it can be seen that the proposed method has obtained better DoA estimation performance. In Figure 5a, we have marked the pilot signal intervals for the proposed algorithm during the period of 400-600 ms and for the improved MUSIC algorithm during the period of 450-500 ms. For the conventional method, it can be seen that during 450-500 ms, the UAV's mobility causes the DoA estimation errors to increase. These errors are larger in Figure 5a than in Figure 5b due to the higher speed of the UAV. The proposed prediction method, on the other hand, continues tracking the variations in the DoA value, causing the errors to remain low during the DoA estimation interval. Besides during the intervals, the proposed algorithm is also more accurate on DoA estimations via pilot signals, e.g., at 600 ms. This is because it explores LoS signals directly, thus avoiding the destructive influence of the angular spread effect. It should be noted that if no LoS paths exist, the proposed method can be less accurate. However, according to 3GPP Release 15, a major proportion of the UAVs incorporated into the cellular system are maintaining LoS connections with the BSs, and therefore the proposed algorithm can be found useful in many civil UAV application scenarios.
The timely updating of the DoA estimations directly contributes to the performance of beamforming, as seen in Figure 6c,d for the two different speed scenarios. The red and blue lines represent the outputs of the DL beamformer combined with the proposed algorithm and the improved MUSIC algorithm, respectively, and the black and yellow lines represent the outputs of the LCMV beamformer combined with the proposed algorithm and the improved MUSIC algorithm, respectively. We summarize the statistics of these results in Table 3, where the average gains and standard variances of the DL and LCMV results are provided. In both figures, the proposed scheme using the DL beamformer exhibits the highest performance, with average gains of 21.92 dB for the scenario in Figure 6c and 21.95 dB for the scenario in Figure 7b, whereas the performance of the conventional scheme using the LCMV beamformer is the worst, with corresponding average gains of −2.71 dB and −0.88 dB, respectively. It is also observed that the DL beamformer achieves much higher gains than the LCMV beamformer. However, it should be noted that its computational complexity is also much higher. For both beamformers, there are no significant margins between the low-speed and high-speed cases in terms of the average SINR or variance, thus demonstrating the stability of the proposed beamforming scheme.

Performance Comparison for Varying Pilot Signal Intervals
In this subsection, we examine the beamforming performance of the proposed and conventional schemes considering varying pilot signal intervals. For both schemes, we consider three experiments based on scenarios with interval values of 10 ms, 50 ms and 200 ms. The UAV is flying along a curve with 80 m radius, and the speed is fixed at 160 km/h. The initial distance between the UAV and the BS is 100 m. The SINR outputs of the LCMV and DL beamformers are shown in Figure 7c,d, respectively.  It is observed that the proposed scheme maintains higher performance for both beamformers. Furthermore, the stability of the proposed method is convincingly demonstrated here since the three cases show very similar DoA estimation and beamforming performance, for both the LCMV and DL beamformers. Thus, we conclude that an FSI exchange frequency of 5 Hz is sufficient to achieve high and stable beamforming performance. In this way, the proposed algorithm can save significant training overhead. Accordingly, we compare the training overheads of the two schemes. Based on the previous analysis, we assume that the proposed method receives telemetry messages every 200 ms. Given that the telemetry information already exists on the Command&Control link for UAV-cellular communications, the corresponding overhead is incurred even if the FSI is not applied for DoA estimation; therefore, the traditional method will also require the consumption of a minimum of 5 Hz for telemetry messages. Thus, the proposed method can achieve overhead reductions of 100, 20 and 5 pilot signals per second relative to the conventional scheme when using pilot signal intervals of 10 ms, 50 ms and 200 ms, respectively.
We also note that the pilot signal transmission interval mainly affects the conventional scheme by magnifying the accumulated error of the improved MUSIC algorithm. In the traditional scheme, the DoA information cannot be updated until a new pilot signal is received and processed. Thus, this scheme suffers from large errors when the DoA estimation interval is 200 ms, and the gain also decreases dramatically.

Performance Evaluation for the Random Mobility Scenario
From the previous results, it can be concluded that DoA estimation is crucial for beamforming performance. When all other conditions remain the same, the beamformer's SINR output is inversely proportional to the DoA error. However, the routes of UAVs are more complex than those of most traditional ground vehicles that are served in cellular networks. The varying trajectories and flexible headings of UAVs lead to constant variation in the DoA information and therefore may cause turbulent fluctuations in the output of DoA estimation algorithms.
In this subsection, we use a UAV random walk model to further verify the suitability of the proposed algorithm for UAVs with high mobility. In the model, the initial distance between the UAV and the BS is 490 m. The amplitude of the UAV velocity is constant at 160 km/h, whereas its direction varies every 300 ms. The azimuth and polar angles of the velocity are chosen randomly between 0 • and 360 • and between −60 • and 60 • , respectively. The DoA estimation interval is 200 ms. The random walk process lasts 6000 ms and is executed 20 times at each point, and the average absolute errors are presented in Figure 8. It can be seen that the proposed algorithm maintains accurate DoA estimations throughout the whole process. The overall stability is also remarkable: the average absolute errors remain between 0.08 • and 0.36 • . Compared with the errors reported in the previous subsections, the output is more accurate here. This is because of the longer distance between the UAV and the BS, which causes the same position errors to correspond to smaller DoA estimation errors. Thus, the proposed DoA estimation algorithm maintains more accurate DoA estimations over longer distances.

Conclusions
In this paper, we investigate high performance beamforming for communications between UAVs and cellular networks. The UAV's signal is estimated by a proposed DoA estimation algorithm and its variations are also predicted by using position information. The proposed algorithm maintains precise DoA estimations, and is of low complexity. Meanwhile, by using DoA information from the Command&Control Links instead of from pilot signals, the proposed algorithm incurs much less training overhead to the network. The algorithm is applicable to beamforming, based on which a beamforming scheme with real-time steering vector updates is also proposed. As proved by simulation results, by applying the proposed beamforming scheme higher DoA accuracy and more stable beamforming gains are obtained compared with traditional beamforming methods. In the meantime, training overhead is proved to be reduced significantly.