1. Introduction
With the development of beamforming technology and unmanned aerial vehicle swarm (UAVs) technology, the antenna array elements equipping the UAVs are combined into a distributed beamforming system through collaboration among the swarm to achieve directional high gain signal transmission. Therefore, using beamforming technology to improve the combat capability of UAVs in a fierce confrontation environment is a growing trend [
1]. It is obvious that taking UAVs as an array has many advantages, such as using low-cost UAVs to achieve high gain signal transmission, which is also known as the beam pattern, but the targets of the array are often located in the near-field, and the motion error [
2,
3,
4,
5] of a single UAV as an array element will deteriorate the beam pattern formed by the array. At present, the array position error calibration is mainly founded on the calibration method of three or more known radiant sources, based on far-field signals, or the rotating array calibration method, based on a single known as a radiant source [
6,
7,
8], but the latter is, in essence, also a multi-radiant source calibration method exchanging time for space. Reference [
2] proposes using three simultaneously existing auxiliary sources to calibrate the position error of the array, but this method can only use the Taylor first-order expansion of the array steering vector with error when the position error has a very small disturbance, which means the position error cannot be corrected by this method when it is too large. In array beamforming, the beam cannot be formed when an array element position error exists, but as long as the phase error caused by the position error is compensated for, the beam can be formed effectively. It is easier to solve a one-dimensional phase error than to solve two- or three- dimensional position errors, and the phase error can be obtained by using a single known radiant source as a reference source [
9,
10,
11,
12]. As for the solution to phase error, based on the information of the reference source, the phase error can be estimated using the phase gradient (PG) value of the received signal between the array elements [
13,
14,
15,
16,
17], but this method requires a high signal-to-noise ratio (SNR), and under the condition of low SNR, the deviation of the phase error estimation value from the real value is large. Phase error can also be solved by iteration. The authors of [
18,
19,
20,
21] take the entropy of the image as the objective function and estimate the phase error by using the Newton iterative method, but this method presents problems, such as the selection of the initial value affecting the algorithm time, causing the algorithm to fall into local optimization. When the expected main-lobe direction of beamforming is the same as the reference source direction, the phase error in main-lobe direction can be compensated for by the estimated phase error compensation function, but when it is not the same as the reference source direction, the mismatch of the phase error compensation function occurs, which will affect the beamforming pattern. Therefore, it is necessary to study the influence of phase error compensation function mismatch, which is the same as the division of the area effectively compensated by the estimated phase error compensation function. The authors of [
22] propose that when the near-field radial compensation filter is used to filter out interference sources, the beamforming pattern will deteriorate when the sources are far away from the array, but the study does not specifically analyze the pattern deterioration level. The work in [
23] proposes a method to form a multi-focus antenna array in the near-field, with specified amplitude and phase conditions for the multi-focus problem of the near-field array, but this study does not involve the analysis of the mismatch of the single-focus filter.
Based on the above analysis, firstly, a near-field array beamforming model with UAV position error is constructed in
Section 2, and the model is approximately simplified by the Taylor expansion of the phase difference function of the near-field signal. Secondly, the improved Newton maximum entropy (IN-ME) algorithm is used to estimate and correct the phase error in
Section 3. Under the condition of low SNR, taking the single known radiant source as a reference source, the initial value of the phase error is estimated through the PG information of the received array signal, and then, taking the maximum entropy function of the array synthetic power at the reference source as the objective function, the Newton iterative solution to the phase error is carried out. After the phase error is compensated, beamforming is carried out. Thirdly, when phase error estimation based on proposed IN-ME algorithm is used to compensate for beamforming, if the array forms beam scanning, the phase compensation function will be mismatched. Based on the phase error compensation function of the IN-ME algorithm, the effectively compensated area which can be compensated for to effectively form the beam pattern is further divided in
Section 4, laying a foundation for the subsequent spatial region division and UAV path planning. The simulation results in
Section 5 show the effectiveness of the algorithm, and the conclusions are given in
Section 6. The adverse influence of UAV position errors on near-field array beamforming is effectively suppressed, and the effectively compensated area, with phase compensation function mismatch, is divided.
3. IN-ME Phase Error Compensation Algorithm Based on Reference Source
Only by compensating for phase error can the array effectively form a beam pattern. We propose the IN-ME algorithm, which, compared to Newton maximum entropy (N-ME) algorithm, can obtain higher convergent value within a shorter period, as well as form a lower sidelobe-level beam pattern with a higher SNR.
3.1. Phase Error Estimation Based on Newton Maximum Entropy Algorithm
Entropy is a measurement of information uncertainty. The maximum entropy principle is a criterion of probabilistic model learning, and when the probability distribution of random variables obeys uniform distribution, the entropy reaches its highest level. As for beamforming, we expect that after phase error compensation, the synthetic powers in the main-lobe direction of the array’s received signals in all snapshots increase, which becomes a multi-objective optimization problem. If we regard the ratio of synthetic power in the main-lobe direction of the array’s received signals in one snapshot to the sum of the synthetic powers in all snapshots as the probability distribution, compared to the maximum entropy principle, it can be seen as that when this ratio obeys uniform distribution, the synthetic power in every snapshot can increase, which means that the phase error compensation is effective for the signals in every snapshot. Therefore, the objective function should be the maximum entropy of that ratio.
Assuming that there is a reference source in the target scene, if the array beamforming is expected to form a main-lobe at the source, the weighted value of the array can be obtained as
. When the snapshots are certain,
, representing the total ideal synthetic power of all snapshots at the reference point, is a constant, and can be obtained by
where,
SS is the total number of snapshots, and
is the synthesized signal at the reference source after ideally weighting and summing the received signals in the
ssth snapshot. The synthetic power in the main-lobe direction of the array’s received signals in one snapshot can be obtained by
where,
is the received signal in the
ssth snapshot of the array element, and
is the phase error of the
nth array element to be estimated. Because
is a constant, the entropy objective function of
can be written as
The entropy objective function
is the function of the phase
to be estimated. Thus, the phase estimation based on the maximum entropy can be expressed as
Theorem 2. ([
25])
Suppose that ,
which is the second derivation of
,
is continuous in an open neighborhood of , and that and is negative definite. Then is a strict local maximizer of .
Proof. Because , which is also known as the Hessian, is a continuous and negative definite at , we can choose a radius so that remains negative definite for all x in the open ball . Taking any nonzero vector with , we have , and so , where, for some . Since , we have , and therefore , giving the results. □
By combining Theorem 1 and Theorem 2, we can use the Newton method for Formula (20). There is an iterative solution of
as
where, superscript
indicates the
lth iteration. To solve Equation (21), the first and second derivation of
against
need to be calculated. The first derivative expression can be expressed as
The second derivative expression can be expressed as
where,
where,
represents the real part operation.
In practical cases, the phase errors in signals received by different array elements are relatively independent from each other [
26]. Therefore, the phase error of each array element is searched separately, and in Formula (23), only the diagonal elements of the Hessian, which is also the second derivative
, of the entropy with respect to phase error are derived, while the off-diagonal elements
are ignored.
By taking Equations (24) and (25) into Equation (22), the analytical formula of the first derivative is obtained
By taking Equation (24) to Equation (28) into Equation (23), the analytical formula of the first derivative is obtained
After performing phase error calibration for Equation (16), the beamforming with error calibration is
where,
represents the Hadamard product of the matrix, and
.
3.2. Initial Value Estimation of Phase Error Based on PG
According to the principle of Newton’s maximum entropy algorithm, the selection of the initial value of the phase error will directly affect the iteration efficiency. If the initial value of the phase error is randomly selected, such as directly setting the initial value of phase error to zero, which is the N-ME algorithm, it may not only greatly slow down the convergence speed of the algorithm, but also make the algorithm fall into the local optimal solution. Therefore, a method based on a reference source to select the initial value of the phase error is proposed to improve the N-ME algorithm.
When there are noises and array position errors, after weighting the received signals, the initial received signals of the array without phase error calibration can be written as
where,
represents conjugation,
represents the initial output signal of the
nth array element, and
represents the initial phase error of the
nth array element caused by noise and array position error.
Ideally, after weighting the received signals, the phase error caused by different distances is fully compensated, and the PG between the output signals of the array elements is constant. For discrete sequences, the phase difference between the output signals of adjacent array elements is the PG. We define the correlation sequence of the array as
The PG between adjacent elements is estimated as
By taking the phase of the first array element as a reference, the initial value of phase error to be compensated for each array element is estimated by calculating the cumulative sum between adjacent array elements. There is
where,
represents the cumulative sum.
The IN-ME phase error compensation algorithm, based on the reference source, is shown in Algorithm 1.
Algorithm 1 The IN-ME phase error compensation algorithm based on reference source. |
>Input: received signals of array , reference source Output: beamform after phase error calibration Step 1: Calculate weighted signal, weights , and weighted signal ; Step 2: Calculate the correlation sequence of the array, ; Step 3: Estimate initial value of phase error based on PG: 1. Estimate the PG between adjacent elements, ; 2. Estimate initial value of phase error, ; Step 4: Improved Newton Maximum Entropy Algorithm: 1. Calculate the array’s synthetic signal at the reference source, ; 2. Form the entropy objective function, ; 3. Calculate the first and second derivatives of to , , ; 4. Use Newton iterative method for optimization, ; Step 5: Beamform after error calibration, . |
4. Effectively Compensated Area
The phase error compensation based on phase error compensation function of the IN-ME algorithm is aimed at the situation in which the expected main-lobe of beamforming is just at the reference source . When the expected main-lobe of beamforming is not at the reference source, there is a mismatch between the phase error compensation function and the actual phase error. Thus, it is necessary to discuss the effectively compensated area based on the phase error compensation function of the IN-ME algorithm.
Assuming that the position errors of the array elements are independent of each other and meet the normal distribution with the mean value of 0 and the variance of
, which is, for the
nth UAV element,
and
, and the two-dimensional joint distribution of the two also meets the normal distribution, according to the “
rule”, the elements can be regarded approximately as that they are located in a circle with the radius centered on the ideal array element position, as shown in
Figure 2. When the expected main-lobe of beamforming is at the source
, but the reference source of the IN-ME algorithm is
, the effectively compensated area based on the phase error compensation function of the IN-ME algorithm can be approximately solved by geometric methods.
Obviously, when the actual position of the array element is on the circle, the phase error is larger than when it is in the circle. Point
A is set as a moving point on the circle to represent the actual position of the
nth array element. When the two sources,
and
, have relative phase errors
, it is considered that the compensation is effective [
27,
28]. There is
where,
and
are respectively the
axis and
axis coordinates of point
A.As shown in
Figure 2, crossing point
A, the plumb line
AF, which is perpendicular to
PO, is made, and the point of intersection is point
F. The plumb line
AE, which is perpendicular to
BO, is made, and the point of intersection is point
E. We let
,
, and
. When
and
, the electromagnetic wave from
or
to any point on the circle can be approximately simplified to the plane wave, and we have
and
. Then Equation (36) is approximately simplified as
By using sum-to-product addition formulas, there is
. Then there is
As shown in
Figure 2, the perpendicular line of the
angular bisector which is made through the center
O of the circle intersects with the circle at
and
. Equation (38) shows that the expected point
A which can get max(
fit) approximately locates at
or
, according to
The area of
that can be compensated based on
for the
nth array element can be obtained by
After finding out the upper and lower bounds of the area of
the phase error for each array element can be effectively compensated, and the area can be divided by finding the intersection of the area of all array elements. As shown in
Figure 3, the connecting line from the
nth array element to
is marked as line
, which is the angular bisector of
, and the line towards the negative angle direction is line
, while the line towards the positive angle direction is line
. The intersection point of line
and line
is marked as point
. The intersection point of line
and line
is marked as point
. The line
parallel to line
is started from point
, and the connecting line between
and
is marked as line
. The intersection point of
and
is marked as point
, and line
parallel to
is started from
. The connecting line between
and
is marked as line
. Then the effectively compensated area of the array for
based on
is within the area encircled by four sides of
g,
e,
d, and
f.
As shown in
Figure 3, the problem of finding the effectively compensated area is transformed into the problem of finding the intersection point of the lines, and the equations of the black line cluster
can be expressed as
The equation of the green line cluster
can be expressed as
The equation of the blue line cluster
can be expressed as
Let
,
, and
, and coordinates of
H,
D and
G can be obtained by
The analytic expressions of the red four sides
g,
e,
d, and
f are
The area encircled by the four edges shown in Equation (45) is the effectively compensated area of that can be compensated based on IN-ME algorithm.
By estimating and compensating the phase errors, the phase errors of the received signals caused by the position of the array elements in the array are calibrated so that the antenna array of UAVs can effectively carry out beamforming. The analysis of the effectively compensated area based on the reference source can provide a reference for the target area division and path planning of UAVs.
6. Conclusions
In order to solve the problems of UAV near-field beamforming, UAV position error compensation, and the effectively compensated area of compensation, this paper firstly constructs a near-field array beamforming model based on element position error, and approximately simplifies the model using the Taylor expansion of the near-field signal phase difference function. Moreover, an IN-ME algorithm is proposed to estimate and compensate for the phase error. Based on the known reference source signal, the initial value of the phase error is estimated through the PG information of the array’s received signal, and then the array signal entropy objective function is established to iteratively optimize the phase error value using the Newton iteration method, and the beam pattern is formed after the phase error compensation. Thirdly, when the beam scanning leads to the phase compensation function mismatch, based on the phase error compensation function of the IN-ME algorithm, the effectively compensated area is further divided, which lays a foundation for the subsequent area region division and UAV path planning. Finally, the validity of the conclusion is verified by simulation. By using the IN-ME algorithm proposed in this paper, the influence of a single UAV’s position error on array beamforming is effectively suppressed, and the effectively compensated area is divided because of phase compensation function mismatch, which provides a theoretical basis for the UAV beamforming application in wireless communication, electronic reconnaissance, and jamming. As the signal wavelength becomes shorter, the adverse influence of the UAV’s position error increases, and the beam pattern will deteriorate sharply, which is worthy of further research.