CN116918274A - Method, device and antenna system for estimating wave beam arrival angle - Google Patents
Method, device and antenna system for estimating wave beam arrival angle Download PDFInfo
- Publication number
- CN116918274A CN116918274A CN202180094736.8A CN202180094736A CN116918274A CN 116918274 A CN116918274 A CN 116918274A CN 202180094736 A CN202180094736 A CN 202180094736A CN 116918274 A CN116918274 A CN 116918274A
- Authority
- CN
- China
- Prior art keywords
- antenna array
- array
- antenna
- received signal
- virtual
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 62
- 238000003491 array Methods 0.000 claims abstract description 225
- 238000012545 processing Methods 0.000 claims abstract description 22
- 239000013598 vector Substances 0.000 claims description 60
- 239000011159 matrix material Substances 0.000 claims description 28
- 238000004891 communication Methods 0.000 claims description 25
- 238000010276 construction Methods 0.000 claims description 12
- 238000004590 computer program Methods 0.000 claims description 11
- 238000012549 training Methods 0.000 claims description 11
- 230000001419 dependent effect Effects 0.000 claims description 8
- 238000012986 modification Methods 0.000 abstract description 5
- 230000004048 modification Effects 0.000 abstract description 5
- 238000013461 design Methods 0.000 description 16
- 238000010586 diagram Methods 0.000 description 12
- 230000008569 process Effects 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 230000001360 synchronised effect Effects 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000005562 fading Methods 0.000 description 2
- 239000013307 optical fiber Substances 0.000 description 2
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
- G01S3/04—Details
- G01S3/12—Means for determining sense of direction, e.g. by combining signals from directional antenna or goniometer search coil with those from non-directional antenna
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
- G01S3/14—Systems for determining direction or deviation from predetermined direction
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/08—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
The application discloses a method, a device and an antenna system for estimating a beam arrival angle. The method comprises the following steps: constructing a set of non-uniform linear virtual antenna arrays based on at least two sets of uniform linear antenna arrays; establishing a receiving signal model of a virtual antenna array according to the relative position information between each group of antenna arrays in at least two groups of uniform linear antenna arrays and the receiving signal model of each group of antenna arrays; and estimating the beam arrival angle of the virtual antenna array based on a received signal model of the virtual antenna array by adopting a spatial super-resolution angle estimation algorithm. The method can obtain the angle estimation performance of an ideal array approaching to the same aperture as the virtual antenna array by only carrying out joint processing on the received signals of at least two groups of uniform linear antenna arrays without any actual modification on the physical structure of the antenna array, and can improve the angle estimation performance of the antenna system on the basis of not increasing the physical complexity of the antenna system.
Description
The present application relates to the field of mobile communications technologies, and in particular, to a method, an apparatus, and an antenna system for estimating an arrival angle of a beam.
In radio positioning or directional technology of an antenna system, estimating a beam arrival angle (angle estimation for short) is a basic technology, that is, processing signals received by an antenna array to realize estimation of a source direction of the signals. The angle estimation performance of the antenna system is limited by the aperture of the antenna array, and at present, if the angle estimation accuracy is to be improved, the physical complexity of the antenna system needs to be increased, for example, an antenna array with more array elements or larger aperture is designed.
Therefore, how to improve the angle estimation performance of the antenna system without increasing the physical complexity of the antenna system is a technical problem to be solved by the application.
Disclosure of Invention
The application provides a method, a device and an antenna system for estimating a beam arrival angle, which are used for improving the angle estimation performance of the antenna system on the basis of not increasing the physical complexity of the antenna system.
In a first aspect, a method of estimating an angle of arrival of a beam is provided, which can be applied to a baseband processing unit BBU or to a chip in the BBU. Taking the BBU as an example, the method may include: constructing a group of non-uniform linear virtual antenna arrays based on at least two groups of uniform linear antenna arrays, wherein the virtual antenna arrays meet far field hypothesis conditions; establishing a receiving signal model of a virtual antenna array according to the relative position information between each group of antenna arrays in at least two groups of uniform linear antenna arrays and the receiving signal model of each group of antenna arrays; and estimating the beam arrival angle of the virtual antenna array based on a received signal model of the virtual antenna array by adopting a spatial super-resolution angle estimation algorithm.
The embodiment of the application constructs a group of non-uniform linear virtual antenna arrays based on at least two groups of uniform linear antenna arrays, the construction process does not carry out any actual modification on the physical structure of the antenna arrays, but only carries out joint processing on the received signals of the at least two groups of uniform linear antenna arrays, thereby obtaining the angle estimation performance of an ideal array which approximates to the same aperture as the virtual antenna arrays, and improving the angle estimation performance of the antenna system on the basis of not increasing the physical complexity of the antenna system.
In one possible design, the BBU may build a received signal model of the virtual antenna array comprising: determining the beam arrival angle relation and the phase relation of each group of antenna arrays according to the relative position information of at least two groups of uniform linear antenna arrays; the method comprises the steps of rewriting a received signal model of each group of antenna arrays according to a phase relation and a beam arrival angle relation, enabling dependent variables in the received signal model of each group of antenna arrays to be represented through the beam arrival angle of a virtual antenna array, and enabling a received vector of each array element in the received signal model of each group of antenna arrays to be represented by taking a received vector of a first array element as a reference, wherein the first array element is any array element in at least two groups of uniformly linear antenna arrays; and determining a receiving signal model of the virtual antenna array based on the rewritten receiving signal models of the antenna arrays of each group.
When the received signal model of each antenna array is rewritten, the dependent variable in the received signal model of each antenna array is expressed by the beam arrival angle of the virtual antenna array, and the received vector of each array element in the received signal model of each antenna array is expressed by taking the received vector of the first array element as a reference, so that the received signal model of the virtual antenna array can be stacked based on the rewritten received signal model of each antenna array. The implementation is simple and the reliability is high.
The following is an example of at least two sets of uniformly linear antenna arrays consisting of a first antenna array and a second antenna array:
in one possible design, the relative position information of at least two sets of uniformly linear antenna arrays includes: array element spacing d in the first antenna array and the second antenna array is half wavelength lambda/2, whereinLambda is wavelength, the number of array elements of the first antenna array is 2N+1, the number of array elements of the second antenna array is 2M+1, and the center point O of the first antenna array 1 With the centre point O of the second antenna array 2 Distance of K (lambda/2), O 1 With O 2 The included angle between the connection line of the first antenna array and the radial direction is phi 1 ，O 1 With O 2 The included angle between the connecting line of the second antenna array and the radial direction of the second antenna array is phi 2 ；
The received signal model of the first antenna array taking the first array element as the reference point in Q symbol periods is as follows:
Y 1 ＝A 1 S+W 1 ；
the second antenna array uses the first array element as the reference point in the first q symbol periods to receive the signal model as follows:
Y 2 ＝A 2 S+W 2 ；
wherein Y is 1 ＝[Y 11 ,Y 12 ,…,Y 1Q ]，Y 2 ＝[Y 21 ,Y 22 ,…,Y 2Q ,]Is a receiving matrix;
Y 1q ＝[y 11 (q),y 12 (q),…,y 1,2N+1 (q)] T ,Y 2q ＝[y 21 (q),y 22 (q),…,y 2,2M+1 (q)] T for the received vector in the q-th symbol period, T represents a transpose operation;
is a guide vector; wherein θ is 1 Angle, θ, for the beam of the first antenna array 2 An angle is reached for the beam of the second antenna array;
S＝[S 1 ,S 2 ,…,S Q ]for Q pilot symbolsForming training sequence vectors;
W 1 ＝[W 11 W 12 ,…,W 1Q ]，W 2 ＝[W 21 W 22 ,…,W 2Q ]is a noise matrix;
W 1q ＝[W 11 (q),W 12 (q),…,W 1,2N+1 (q)] T ；
W 2q ＝[W 21 (q),W 22 (q),…,W 2,2M+1 (q)] T ；
W 1q 、W 2q a noise vector for the q-th symbol period;
the array elements of the virtual antenna array are the array elements of the first antenna array and the second antenna array, and the center of the virtual antenna array is O 1 、O 2 A midpoint of the connection line;
the beam angle of arrival relationship of each group of antenna arrays determined from the relative position information between each group of antenna arrays, comprising:
θ 1 ＝φ 1 +θ； (1)
θ 2 ＝φ 2 -θ； (2)
wherein θ is the beam arrival angle of the virtual antenna array;
rewriting a received signal model of each group of antenna arrays according to the phase relation and the beam arrival angle relation, comprising:
the receiving signal taking the first array element of the first antenna array as the reference point in the q-th symbol period is as follows:
With the centre point O of the first antenna array 1 The position is the reference point, the first antenna array is arranged according to the formula (1)The received signals of the n array elements in the q-th symbol period are rewritten as follows:
where n=1, 2, …,2n+1.
The m-th array element of the second antenna array takes the first array element of the second antenna array as a reference point in the q-th symbol period as follows:
with the centre point O of the second antenna array 2 The position is the reference point, and the received signal of the m-th array element of the second antenna array in the q-th symbol period is rewritten as follows according to the formula (2):
with the central point O of the first antenna element 1 The position is a reference point according to the formulas (1) (2) and O 2 With O 1 And (2) rewrites the received signal of the mth element of the second antenna array in the q-th symbol period to:
wherein m=1, 2, …,2m+1;
determining a received signal model of the virtual antenna array based on the rewritten received signal models of the antenna arrays of each group, comprising:
according to formulas (4) and (7), stacking the received signals of the virtual antenna array into a matrix of (2N+2M+2) x Q dimensions, and obtaining a received signal model of the virtual antenna array as follows:
Y＝[y(1)，y(2)，…，y(Q)]＝AS+W；
wherein,
in another possible design, at least two sets of uniform linear antenna arrays are made up of a first antenna array and a second antenna array;
The relative position information of the at least two sets of uniformly linear antenna arrays includes: the array element interval d in the first antenna array and the second antenna array is half wavelength lambda/2, wherein lambda is wavelength, the number of the array elements of the first antenna array is 2N+1, the number of the array elements of the second antenna array is 2M+1, and the center point O of the first antenna array 1 With the centre point O of the second antenna array 2 Distance of K (lambda/2), O 1 With O 2 The included angle between the connection line of the first antenna array and the radial direction is phi 1 ，O 1 With O 2 The included angle between the connecting line of the second antenna array and the radial direction of the second antenna array is phi 2 ；
The received signal model of the first antenna array over the first q symbol periods is:
Y 1 ＝A 1 S+W 1 ；
the received signal model of the second antenna array over the first q symbol periods is:
Y 2 ＝A 2 S+W 2 ；
wherein Y is 1 ＝[Y 11 ,Y 12 ,…,Y 1Q ]，Y 2 ＝[Y 21 ,Y 22 ,…,Y 2Q ,]Is a receiving matrix;
Y 1q ＝[y 11 (q),y 12 (q),…,y 1,2N+1 (q)] T ,Y 2q ＝[y 21 (q),y 22 (q),…,y 2,2M+1 (q)] T for the received vector in the q-th symbol period, T represents a transpose operation;
is a guide vector; wherein θ is 1 Angle, θ, for the beam of the first antenna array 2 An angle is reached for the beam of the second antenna array;
S＝[S 1 ,S 2 ,…,S Q ]training sequence vectors formed by Q pilot symbols;
W 1 ＝[W 11 W 12 ,…,W 1Q ]，W 2 ＝[W 21 W 22 ,…,W 2Q ]is a noise matrix;
W 1q ＝[W 11 (q),W 12 (q),…,W 1,2N+1 (q)] T ；
W 2q ＝[W 21 (q),W 22 (q),…,W 2,2M+1 (q)] T ；
W 1q 、W 2q a noise vector for the q-th symbol period;
the array elements of the virtual antenna array are the array elements of the first antenna array and the second antenna array, and the center of the virtual antenna array is O 1 、O 2 A midpoint of the connection line;
the received signal model of the virtual antenna array is:
Y＝[y(1)，y(2)，…，y(Q)]＝AS+W；
wherein,
virtual days in the q-th symbol periodEach array element in the line array uses the central point O of the first antenna array 1 The received signal for the reference point satisfies the following relationship:
wherein θ 1 ＝φ 1 +θ，θ 2 ＝φ 2 - θ; n=1, 2, …,2n+1; m=1, 2, …,2m+1; θ is the beam angle of arrival of the virtual antenna array.
The two designs adopt two groups of uniform linear antenna arrays to construct a virtual antenna array, and the implementation is simple and the reliability is high.
In one possible design, the BBU may further calculate an initial estimate of the beam arrival angle of each of the at least two groups of antenna arrays based on the received signal model of each of the at least two groups of antenna arrays; and then determining the search range of the spatial super-resolution angle estimation algorithm according to the initial estimation value of the beam arrival angle of each group of antenna arrays.
For example, the search range may be:
wherein, for an initial estimate of the beam angle of arrival of the first antenna array,is the initial estimated value of the beam arrival angle of the second antenna array, and delta is a preset value.
Therefore, the search range of a spatial super-resolution angle estimation algorithm used by the BBU in estimating the wave beam arrival angle of the virtual antenna array can be reduced, and the precision and efficiency of angle estimation are further improved.
In one possible design, delta is inversely related to the signal-to-noise ratio. In other words, the higher the signal-to-noise ratio, the smaller the Δ, and the lower the signal-to-noise ratio, the higher the Δ.
In this way, the accuracy of the angle estimation can be further improved.
In one possible design, after estimating the beam arrival angle of the virtual antenna array based on the received signal model of the virtual antenna array by using the spatial super-resolution angle estimation algorithm, the BBU may further update the initial final estimated value of the beam arrival angle of each group of antenna arrays according to the beam arrival angle of the virtual antenna array, to obtain the final estimated value of the beam arrival angle of each group of antenna arrays.
Thus, the angle estimation precision of each antenna array can be improved.
In a second aspect, an apparatus for estimating an angle of arrival of a beam is provided, such as a BBU, or a chip in a BBU. The apparatus comprises means for performing the method described in the first aspect or any one of the possible designs of the first aspect.
Illustratively, the apparatus may include: a construction module for constructing a set of non-uniform linear virtual antenna arrays based on at least two sets of uniform linear antenna arrays; establishing a receiving signal model of a virtual antenna array according to the relative position information between each group of antenna arrays in at least two groups of uniform linear antenna arrays and the receiving signal model of each group of antenna arrays; and the estimation module is used for estimating the wave beam arrival angle of the virtual antenna array based on the received signal model of the virtual antenna array by adopting a spatial super-resolution angle estimation algorithm.
In a third aspect, an antenna system is provided, comprising: each remote radio unit comprises a group of uniform linear antenna arrays, and the remote radio units are used for receiving wireless signals; and the baseband processing unit is in communication connection with at least two remote radio units and is used for executing the method in the first aspect or any one of the possible designs of the first aspect. For example, the baseband processing unit may be configured to: constructing a group of non-uniform linear virtual antenna arrays based on at least two groups of uniform linear antenna arrays, wherein the virtual antenna arrays meet far field hypothesis conditions; establishing a receiving signal model of a virtual antenna array according to the relative position information between each group of antenna arrays in at least two groups of uniform linear antenna arrays and the receiving signal model of each group of antenna arrays; and estimating the beam arrival angle of the virtual antenna array based on a received signal model of the virtual antenna array by adopting a spatial super-resolution angle estimation algorithm.
In a fourth aspect, a communication device is provided, comprising a processor and a memory; the memory is used for storing computer execution instructions; the processor is configured to execute computer-executable instructions stored in the memory to cause the communication device to perform the method as described in the first aspect or any one of the possible designs of the first aspect.
In a fifth aspect, a communication device is provided, comprising a processor and an interface circuit; the interface circuit is used for receiving the code instruction and transmitting the code instruction to the processor; the processor executes code instructions to perform the method as described in the first aspect or any of the possible designs of the first aspect.
In a sixth aspect, there is provided a computer readable storage medium for storing instructions which, when executed, cause the method described in the first aspect or any one of the possible designs of the first aspect to be implemented.
In a seventh aspect, a chip is provided, the chip being coupled to a memory for reading and executing program instructions stored in the memory, implementing the method as described in the first aspect or any one of the possible designs of the first aspect.
In an eighth aspect, there is provided a computer program product comprising instructions stored therein, which when run on a computer, cause the computer to perform the method described in the first aspect or any of the possible designs of the first aspect.
Technical effects that may be achieved by the above-described possible designs of any one of the second aspect to the eighth aspect and any one of the second aspect are referred to for description of the technical effects that may be achieved by the above-described first aspect and the corresponding possible designs, and detailed descriptions thereof are not repeated here.
Fig. 1 is a schematic diagram of an antenna system according to an embodiment of the present application;
fig. 2 is a flowchart of a method for estimating an angle of arrival of a beam according to an embodiment of the present application;
fig. 3 is a schematic diagram of another antenna system according to an embodiment of the present application;
FIG. 4A is a graph showing MSE statistics for angle estimation errors;
FIG. 4B is a graph showing MAE statistics for angle estimation errors;
fig. 5 is a schematic structural diagram of an apparatus for estimating an angle of arrival of a beam according to an embodiment of the present application;
fig. 6 is a schematic structural diagram of a communication device according to an embodiment of the present application.
The technical scheme of the embodiment of the application can be applied to various communication systems, such as: the technical solution Of the embodiment Of the present application may be adopted in the fifth generation (5th generation,5G) communication system, the sixth generation (6th generation,6G) communication system, or other future evolution systems, or other various wireless communication systems using a wireless access technology, as long as there is a need to estimate the Angle Of Arrival (AOA) Of a beam in the communication system.
A communication system (e.g., a 5G communication system) may employ the millimeter wave band to provide high-rate wireless communication services. The directional performance of the millimeter wave frequency band is strong, and better positioning performance relative to the microwave frequency band can be provided. Meanwhile, millimeter wave frequency band bandwidth resources are rich, and the method is favorable for improving the performance of a communication system by sufficient space diversity gain. However, millimeter wave communication frequency bands are high and fading is serious, resulting in limited transmission distance and reduced corresponding coverage. In addition, millimeter wave penetration is weak and is easily blocked. In order to overcome the problems of communication quality degradation and coverage reduction caused by shielding, a remote radio technology can be adopted to improve the performance of coverage and communication rate. However, especially in the far-radio scenario, the implementation is more complicated if the accuracy of angle estimation is to be improved by designing more array elements or larger-caliber antenna arrays, etc. The present application may be applied to a communication system employing a millimeter wave band, but it should be understood that the present application is not limited thereto.
Referring to fig. 1, an antenna system according to an embodiment of the present application includes at least two remote radio units (Remote Radio Unit, RRU) and a baseband processing unit (Base Band unit, BBU). It should be understood that only two RRUs are taken as examples in fig. 1, namely RRU1 and RRU2, and the number of actual RRUs may not be limited thereto.
The BBU is generally arranged in a central machine room, the RRU is arranged at a far end (near an antenna), and the BBU and the RRU are connected through optical fibers and perform signal transmission. The BBU is mainly responsible for processing all L1 (Layer 1 ), L2 (Layer 2, layer 2), L3 (Layer 3 ), including high-Layer data transmission, scheduling, baseband signal processing, and the like, and the RRU mainly performs up-down conversion, power amplification, filtering, and conversion of medium radio frequency signals.
As shown in fig. 1, a group of uniform linear antenna arrays are disposed on each RRU. The RRU1 includes a plurality of array elements connected in series with each other, and the intervals between adjacent elements are approximately uniform, and the RRU2 includes a plurality of array elements connected in series with each other, and the intervals between adjacent elements are approximately uniform. The aperture of each group of antenna arrays is the ratio of the physical size of the antenna array to the wavelength, i.e. r= (N-1) d/λ, where N is the number of antenna elements, d is the element spacing, and λ is the wavelength.
It should be appreciated that each set of antenna arrays on the antenna system satisfies the far field assumption. By far field hypothesis, it is meant that the distance of the mobile station from the antenna array is much greater than the aperture of the antenna array. For an antenna array that satisfies the far field assumption, the incoming wave signal amplitudes arriving at each element of the antenna array may be considered equal, with only a difference in phase. The mobile station refers to any device capable of wirelessly communicating with the antenna system, for example, a mobile phone, a computer with a mobile terminal device, a portable, pocket, hand-held, a mobile device built in the computer, a wearable device, a vehicle-mounted terminal device, and the like, and the present application is not limited thereto.
It should be understood that only the antenna array on the RRU is shown in fig. 1, and other elements of the RRU, such as a frequency converter, a power amplifier, a filter, etc., are not shown.
Referring to fig. 2, a method for estimating an arrival angle of a beam according to an embodiment of the present application is taken as an example, where the method is applied to the antenna system shown in fig. 1. The method comprises the following steps:
s201, constructing a group of non-uniform linear virtual antenna arrays by the BBU based on at least two groups of uniform linear antenna arrays, wherein the virtual antenna arrays meet far field hypothesis conditions; and establishing a receiving signal model of the virtual antenna array according to the relative position information between each group of antenna arrays in the at least two groups of uniform linear antenna arrays and the receiving signal model of each group of antenna arrays.
It should be understood that at least two sets of uniformly linear antenna arrays means that each set of at least two sets of antenna arrays is uniformly linear, and that the at least two sets of antenna arrays as a whole may be linear or non-linear.
The BBU constructs a set of non-uniform linear virtual antenna arrays based on at least two sets of uniform linear antenna arrays, comprising: the BBU determines the physical conditions of at least two groups of uniformly linear antenna arrays, and obtains the physical conditions of the virtual antenna arrays based on the physical conditions of each group of uniformly linear antenna arrays. For example, all the array elements of the virtual antenna array are all the actual array elements in the at least two groups of uniformly linear antenna arrays, and the center of the virtual antenna array is the center of the actual at least two groups of uniformly linear antenna arrays. It should be understood that the construction process does not make any changes to the actual physical structure of the at least two sets of uniformly linear antenna arrays.
For example, a set of non-uniform linear virtual antenna arrays is configured with two sets of uniform linear antenna arrays, referring to fig. 3, the antenna array on the RRU1 is a first antenna array, the antenna array on the RRU2 is a second antenna array, and the relative position information of the first antenna array and the second antenna array includes:
1) The array element interval d in the first antenna array is half-wavelength lambda/2, and the array element interval d in the second antenna array is half-wavelength lambda/2, wherein lambda is the wavelength;
2) The number of array elements of the first antenna array is 2N+1, and the number of array elements of the second antenna array is 2M+1; wherein M, N is a positive integer, M is the same as or different from N, and the application is not limited;
3) Center point O of first antenna array 1 With the centre point O of the second antenna array 2 Distance K (lambda/2);
4)、O 1 with O 2 The included angle between the connection line of the first antenna array and the radial direction is phi 1 ，O 1 With O 2 The included angle between the connecting line of the second antenna array and the radial direction of the second antenna array is phi 2 。
Correspondingly, the array elements of the virtual antenna array based on the first antenna array and the second antenna array are the array elements of the first antenna array and the second antenna array, and the center of the virtual antenna array is O 1 、O 2 The midpoint of the connection line. The aperture of the virtual antenna array is the ratio of the physical size of the virtual antenna array (i.e., the overall physical size of the first antenna array and the second antenna array) to the wavelength, i.e., r= (N-1) d/λ, N is the number of antenna elements, d is the array element spacing, and λ is the wavelength.
It should be appreciated that in addition to the virtual antenna array should satisfy far field hypothesis conditions, the uplink of the antenna system that constitutes the virtual antenna array also satisfies any one or more of the following hypothesis conditions: 1) The channel has sparse multipath characteristic, and the direct path power is far greater than the indirect path power, so that when the power is distributed according to the water injection theorem, only the angle of the direct path can be estimated and shaped in the direct path direction, and the indirect path can be ignored and the channel capacity is hardly lost; 2) The channel is a slow fading channel, and the channel is kept unchanged in one frame time, so that the beam forming cannot be performed before the angle estimation, the received signal to noise ratio is low, the signal to noise ratio is improved by accumulating the energy of a plurality of symbols, and the angle can be reliably estimated; 3) Pilot symbol sequences transmitted by the mobile stations are mutually orthogonal, so that interference among users can be restrained, and the BBU receiving model can be equivalent to a single-user single-path model; 4) The signals received by each array element on the RRU can be directly transmitted to the BBU after the down-conversion; 5) The RRUs connected with the same BBU through optical fibers are strictly synchronized, so that signals of the RRUs can be processed in a combined mode after reaching the BBU; 6) The additive noise on each array element on RRU is independent from the received signal, so that the algorithms such as MUSIC/ESPRIT can form a spatial super-resolution algorithm by utilizing the orthogonality of the signal space and the noise space.
After determining the physical conditions of each group of actual antenna arrays and the physical conditions of the virtual antenna arrays, the BBU can obtain the received signal models of each group of antenna arrays and the relative position information among each group of antenna arrays, and further calculate the received signal models of the virtual antenna arrays by combining the received signal models of the actual antenna arrays.
One possible method of estimation is: determining the beam arrival angle relation and the phase relation of each group of antenna arrays according to the relative position information of at least two groups of uniform linear antenna arrays; the method comprises the steps of rewriting a received signal model of each group of antenna arrays according to a phase relation and a beam arrival angle relation, enabling dependent variables in the received signal model of each group of antenna arrays to be represented through the beam arrival angle of a virtual antenna array, and enabling a received vector of each array element in the received signal model of each group of antenna arrays to be represented by taking a received vector of a first array element as a reference, wherein the first array element is any array element in at least two groups of uniformly linear antenna arrays; and determining a receiving signal model of the virtual antenna array based on the rewritten receiving signal models of the antenna arrays of each group.
Illustratively, a set of non-uniform linear virtual antenna arrays is still constructed from the two sets of uniform linear antenna arrays shown in fig. 3:
The received signal model of the first antenna array taking the first array element as the reference point in Q symbol periods is as follows:
Y 1 ＝A 1 S+W 1 ；
the second antenna array uses the first array element as the reference point in the first q symbol periods to receive the signal model as follows:
Y 2 ＝A 2 S+W 2 ；
wherein Y is 1 ＝[Y 11 ,Y 12 ,…,Y 1Q ]，Y 2 ＝[Y 21 ,Y 22 ,…,Y 2Q ,]Is a receiving matrix;
Y 1q ＝[y 11 (q),y 12 (q),…,y 1,2N+1 (q)] T ,Y 2q ＝[y 21 (q),y 22 (q),…,y 2,2M+1 (q)] T for the received vector in the q-th symbol period, T represents a transpose operation;
is a guide vector; wherein θ is 1 Angle, θ, for the beam of the first antenna array 2 An angle is reached for the beam of the second antenna array;
S＝[S 1 ,S 2 ,…,S Q ]training sequence vectors formed by Q pilot symbols;
W 1 ＝[W 11 W 12 ,…,W 1Q ]，W 2 ＝[W 21 W 22 ,…,W 2Q ]is a noise matrix;
W 1q ＝[W 11 (q),W 12 (q),…,W 1,2N+1 (q)] T ；
W 2q ＝[W 21 (q),W 22 (q),…,W 2,2M+1 (q)] T ；
W 1q 、W 2q a noise vector for the q-th symbol period;
the beam angle of arrival relationship of each group of antenna arrays determined from the relative position information between each group of antenna arrays, comprising:
θ 1 ＝φ 1 +θ； (1)
θ 2 ＝φ 2 -θ； (2)
wherein θ is the beam arrival angle of the virtual antenna array;
referring to fig. 3, the array elements of the virtual antenna array are the array elements of the first antenna array and the second antenna array, and total 2m+2n+2 array elements; the center of the virtual antenna array is O 1 、O 2 Midpoint O of the connection.
After the above information is obtained, the BBU may rewrite the received signal model of each group of antenna arrays according to the phase relationship and the beam arrival angle relationship of the first antenna array and the second antenna array, and the rewrite process includes the following processes of formulas (3) to (8):
The receiving signal taking the first array element of the first antenna array as the reference point in the q-th symbol period is as follows:
with the centre point O of the first antenna array 1 The position is a reference point, and according to the formula (1), the received signal of the nth array element of the first antenna array in the qth symbol period is rewritten as follows:
where n=1, 2, …,2n+1;
the m-th array element of the second antenna array takes the first array element of the second antenna array as a reference point in the q-th symbol period as follows:
with the centre point O of the second antenna array 2 The position is the reference point, and the received signal of the m-th array element of the second antenna array in the q-th symbol period is rewritten as follows according to the formula (2):
with the central point O of the first antenna element 1 The position is a reference point according to the formulas (1) (2) and O 2 With O 1 And (2) rewrites the received signal of the mth element of the second antenna array in the q-th symbol period to:
wherein m=1, 2, …,2m+1;
determining a received signal model of the virtual antenna array based on the rewritten received signal models of the antenna arrays of each group, comprising:
according to formulas (4) and (7), stacking the received signals of the virtual antenna array into a matrix of (2N+2M+2) x Q dimensions, and obtaining a received signal model of the virtual antenna array as follows:
Y＝[y(1)，y(2)，…，y(Q)]＝AS+W (8)
Wherein,
it should be understood that the actual physical structure of the virtual antenna array is the physical structure of the first antenna array and the second antenna array, so that the virtual antenna array is nonlinear in nature, but the estimation process performs angle compensation on the received signals of each array element to virtualize the whole first antenna array and the second antenna array into a linear antenna array.
It should be emphasized that the above estimation process is merely an example and not limited thereto, and other estimation methods are also possible in practical applications, so long as the method of estimating the received signal model of the whole (i.e. the virtual antenna array) formed by at least two antenna arrays based on the received signal model of each antenna array of the at least two antenna arrays falls within the protection scope of the present application.
In addition, the actual antenna arrays constructing the virtual antenna array are not limited to two groups, but may be more, but it should be noted that when one virtual antenna array is constructed by three or more groups of antenna arrays, the center points of all the three or more groups of antenna arrays need to be on the same straight line. The idea of constructing virtual antenna arrays for three or more antenna arrays may refer to the method of constructing virtual antenna arrays by two antenna arrays, for example, the dependent variable in the received signal model of each antenna array in three antenna arrays may be represented by the beam arrival angle of the virtual antenna array, and the received vector of each array element in the received signal model of each antenna array in three antenna arrays is represented by the received vector of the same array element, so as to further stack the received signal model of the virtual antenna array, which will not be described herein.
S202, the BBU adopts a spatial super-resolution angle estimation algorithm, and estimates the wave beam arrival angle of the virtual antenna array based on a received signal model of the virtual antenna array.
Optionally, the spatial super-resolution angle estimation algorithm is a MUSIC algorithm or an ESPRIT algorithm, which is not limited in the present application.
Based on the above, the embodiment of the application constructs a group of non-uniform linear virtual antenna arrays based on at least two groups of uniform linear antenna arrays, and the construction process does not carry out any actual modification on the antenna arrays, but carries out joint processing on the received signals of at least two groups of uniform linear antenna arrays, so as to obtain the angle estimation performance of an ideal array which approximates to the same aperture as the virtual antenna arrays, and can improve the angle estimation performance of the antenna system on the basis of not increasing the physical complexity of the antenna system.
The above effect is illustrated by a set of simulation data:
the simulation parameters are set as follows: the number of array elements of the uniform linear array 1 with half wavelength interval is 9, i.e. n=4; the number of array elements of the uniform linear array 2 with half wavelength interval is 11, i.e. m=5; the distance between the central connecting lines of the two arrays is 20 times of half wavelength, namely K=20; the number of beats (pilot symbol number) q=200; θ＝10°，θ 1 ＝35°，θ 2 =50°; the signal to noise ratio varies from-10 dB to 10dB. The angle estimation algorithm uses the ESPRIT algorithm. The angle estimation performance index selects the variance (MSE) of the angle estimation error, as well as the Mean Absolute Error (MAE).
The simulation results are shown in fig. 4A and 4B. Fig. 4A is a schematic diagram of MSE statistics of angle estimation errors, and curves from top to bottom in fig. 4A are respectively: the MSE of the array 1, the MSE of the array 2, the MSE of the virtual array (namely the virtual antenna array formed by the array 1 and the array 2) and the MSE of an ideal array with the same caliber as the virtual array. Fig. 4B is a schematic diagram of MAE statistics of angle estimation error, and curves from top to bottom in fig. 4B are respectively: the MAE of array 1, the MAE of array 2, the MAE of the virtual array, and the MAE of the ideal array with the same caliber as the virtual array.
From simulation results, it can be seen that the angle estimation performance of the virtual array is higher than that of a single array (i.e., array 1 or array 2), both from the MSE and MAE perspective, and approximates that of an ideal array of the same caliber as the virtual array.
Optionally, before estimating the beam arrival angle of the virtual antenna array based on the received signal model of the virtual antenna array, the BBU may further calculate an initial estimated value of the beam arrival angle of each group of antenna arrays based on the received signal model of each group of antenna arrays of at least two groups of uniform linearity, and then determine a search range of a spatial super-resolution angle estimation algorithm used by the BBU in estimating the beam arrival angle of the virtual antenna array according to the initial estimated value of the beam arrival angle of each group of antenna arrays. The subsequent BBU can search the angle estimation value only in the search range when using the spatial super-resolution angle estimation algorithm to perform angle estimation on the virtual antenna array.
For example, set upFor an initial estimate of the beam angle of arrival of the first antenna array,for the initial estimated value of the beam arrival angle of the second antenna array, the search range of the spatial super-resolution angle estimation algorithm used by the BBU in estimating the beam arrival angle of the virtual antenna array may be:
wherein, delta is a preset value. Alternatively, the Δ is inversely related to the signal-to-noise ratio, i.e., the higher the signal-to-noise ratio, the smaller the Δ, and the lower the signal-to-noise ratio, the higher the Δ.
Therefore, the search range of a spatial super-resolution angle estimation algorithm used by the BBU in estimating the wave beam arrival angle of the virtual antenna array can be reduced, and the precision and efficiency of angle estimation are further improved.
Further optionally, after estimating the beam arrival angle of the virtual antenna array based on the received signal model of the virtual antenna array, the BBU may further update the initial final estimated value of the beam arrival angle of each group of antenna arrays according to the beam arrival angle of the virtual antenna array, to obtain a final estimated value of the beam arrival angle of each group of antenna arrays.
Still taking the two sets of uniform linear antenna arrays shown in fig. 3 as an example, a set of non-uniform linear virtual antenna arrays is constructed: assuming θ as the estimated value of the beam arrival angle of the virtual antenna array composed of the first antenna array and the second antenna array, the final estimated value of the beam arrival angle of the first antenna array can be determined according to the above formulas (1) and (2) Is phi 1 +θ, determining a final estimate of the beam angle of arrival of the second antenna arrayIs phi 2 -θ。
Thus, the angle estimation precision of each antenna array can be improved.
Based on the same technical idea, an embodiment of the present application also provides an apparatus for estimating an angle of arrival of a beam, for example, a BBU or a chip in a BBU, which includes modules/units for performing the method steps shown in fig. 1.
Illustratively, referring to FIG. 5, the apparatus comprises:
a construction module 501, configured to construct a set of non-uniform linear virtual antenna arrays based on at least two sets of uniform linear antenna arrays; establishing a receiving signal model of a virtual antenna array according to the relative position information between each group of antenna arrays in at least two groups of uniform linear antenna arrays and the receiving signal model of each group of antenna arrays;
the estimating module 502 is configured to estimate a beam arrival angle of the virtual antenna array based on a received signal model of the virtual antenna array using a spatial super-resolution angle estimation algorithm.
Optionally, the building module 501 is specifically configured to, when building a received signal model of a virtual antenna array according to relative position information between each group of antenna arrays in at least two groups of uniformly linear antenna arrays and a received signal model of each group of antenna arrays:
Determining the beam arrival angle relation and the phase relation of each group of antenna arrays according to the relative position information of at least two groups of uniform linear antenna arrays;
the method comprises the steps of rewriting a received signal model of each group of antenna arrays according to a phase relation and a beam arrival angle relation, enabling dependent variables in the received signal model of each group of antenna arrays to be represented through the beam arrival angle of a virtual antenna array, and enabling a received vector of each array element in the received signal model of each group of antenna arrays to be represented by taking a received vector of a first array element as a reference, wherein the first array element is any array element in at least two groups of uniformly linear antenna arrays;
and determining a receiving signal model of the virtual antenna array based on the rewritten receiving signal models of the antenna arrays of each group.
Optionally, at least two groups of uniform linear antenna arrays are composed of a first antenna array and a second antenna array;
the relative position information of the at least two sets of uniformly linear antenna arrays includes: the array element interval d in the first antenna array and the second antenna array is half wavelength lambda/2, wherein lambda is wavelength, the number of the array elements of the first antenna array is 2N+1, the number of the array elements of the second antenna array is 2M+1, and the center point O of the first antenna array 1 With the centre point O of the second antenna array 2 Distance of K (lambda/2), O 1 With O 2 The included angle between the connection line of the first antenna array and the radial direction is phi 1 ，O 1 With O 2 The included angle between the connecting line of the second antenna array and the radial direction of the second antenna array is phi 2 The method comprises the steps of carrying out a first treatment on the surface of the The array elements of the virtual antenna array are the array elements of the first antenna array and the second antenna array, and the center of the virtual antenna array is O 1 、O 2 A midpoint of the connection line;
the received signal model of the first antenna array taking the first array element as the reference point in Q symbol periods is as follows:
Y 1 ＝A 1 S+W 1 ；
the second antenna array uses the first array element as the reference point in the first q symbol periods to receive the signal model as follows:
Y 2 ＝A 2 S+W 2 ；
wherein Y is 1 ＝[Y 11 ,Y 12 ,…,Y 1Q ]，Y 2 ＝[Y 21 ,Y 22 ,…,Y 2Q ,]Is a receiving matrix;
Y 1q ＝[y 11 (q),y 12 (q),…,y 1,2N+1 (q)] T ,Y 2q ＝[y 21 (q),y 22 (q),…,y 2,2M+1 (q)] T for the received vector in the q-th symbol period, T represents a transpose operation;
is a guide vector; wherein θ is 1 Angle, θ, for the beam of the first antenna array 2 An angle is reached for the beam of the second antenna array;
S＝[S 1 ,S 2 ,…,S Q ]training sequence vectors formed by Q pilot symbols;
W 1 ＝[W 11 W 12 ,…,W 1Q ]，W 2 ＝[W 21 W 22 ,…,W 2Q ]is a noise matrix;
W 1q ＝[W 11 (q),W 12 (q),…,W 1,2N+1 (q)] T ；
W 2q ＝[W 21 (q),W 22 (q),…,W 2,2M+1 (q)] T ；
W 1q 、W 2q a noise vector for the q-th symbol period;
the construction module 501 determines a beam angle of arrival relationship of each group of antenna arrays according to the relative position information between each group of antenna arrays, including:
θ 1 ＝φ 1 +θ； (1)
θ 2 ＝φ 2 -θ； (2)
wherein θ is the beam arrival angle of the virtual antenna array;
The construction module 501 is specifically configured to, when writing the received signal model of each group of antenna arrays according to the phase relationship and the beam arrival angle relationship:
the receiving signal taking the first array element of the first antenna array as the reference point in the q-th symbol period is as follows:
with the centre point O of the first antenna array 1 The position is a reference point, and according to the formula (1), the received signal of the nth array element of the first antenna array in the qth symbol period is rewritten as follows:
where n=1, 2, …,2n+1.
The m-th array element of the second antenna array takes the first array element of the second antenna array as a reference point in the q-th symbol period as follows:
with the centre point O of the second antenna array 2 The position is the reference point, and the received signal of the m-th array element of the second antenna array in the q-th symbol period is rewritten as follows according to the formula (2):
with the central point O of the first antenna element 1 The position is a reference point according to the formulas (1) (2) and O 2 With O 1 And (2) rewrites the received signal of the mth element of the second antenna array in the q-th symbol period to:
wherein m=1, 2, …,2m+1;
the construction module 501 is specifically configured to, when determining the received signal model of the virtual antenna array based on the rewritten received signal models of the antenna arrays of each group:
According to formulas (4) and (7), stacking the received signals of the virtual antenna array into a matrix of (2N+2M+2) x Q dimensions, and obtaining a received signal model of the virtual antenna array as follows:
Y＝[y(1)，y(2)，…，y(Q)]＝AS+W；
wherein,
optionally, at least two groups of uniform linear antenna arrays are composed of a first antenna array and a second antenna array;
the relative position information of the at least two sets of uniformly linear antenna arrays includes: the array element interval d in the first antenna array and the second antenna array is half wavelength lambda/2, wherein lambda is wavelength, and the number of the array elements of the first antenna array is 2N+1The number of array elements of the second antenna array is 2M+1, and the center point O of the first antenna array 1 With the centre point O of the second antenna array 2 Distance of K (lambda/2), O 1 With O 2 The included angle between the connection line of the first antenna array and the radial direction is phi 1 ，O 1 With O 2 The included angle between the connecting line of the second antenna array and the radial direction of the second antenna array is phi 2 The method comprises the steps of carrying out a first treatment on the surface of the The array elements of the virtual antenna array are the array elements of the first antenna array and the second antenna array, and the center of the virtual antenna array is O 1 、O 2 A midpoint of the connection line;
the received signal model of the first antenna array over the first q symbol periods is:
Y 1 ＝A 1 S+W 1 ；
the received signal model of the second antenna array over the first q symbol periods is:
Y 2 ＝A 2 S+W 2 ；
wherein Y is 1 ＝[Y 11 ,Y 12 ,…,Y 1Q ]，Y 2 ＝[Y 21 ,Y 22 ,…,Y 2Q ,]Is a receiving matrix;
Y 1q ＝[y 11 (q),y 12 (q),…,y 1,2N+1 (q)] T ,Y 2q ＝[y 21 (q),y 22 (q),…,y 2,2M+1 (q)] T For the received vector in the q-th symbol period, T represents a transpose operation;
is a guide vector; wherein θ is 1 Angle, θ, for the beam of the first antenna array 2 An angle is reached for the beam of the second antenna array;
S＝[S 1 ,S 2 ,…,S Q ]training sequence vectors formed by Q pilot symbols;
W 1 ＝[W 11 W 12 ,…,W 1Q ]，W 2 ＝[W 21 W 22 ,…,W 2Q ]is a noise matrix;
W 1q ＝[W 11 (q),W 12 (q),…,W 1,2N+1 (q)] T ；
W 2q ＝[W 21 (q),W 22 (q),…,W 2,2M+1 (q)] T ；
W 1q 、W 2q a noise vector for the q-th symbol period;
the received signal model of the virtual antenna array is:
Y＝[y(1)，y(2)，…，y(Q)]＝AS+W；
wherein,
each array element in the virtual antenna array in the q-th symbol period uses the central point O of the first antenna array 1 The received signal for the reference point satisfies the following relationship:
wherein θ 1 ＝φ 1 +θ，θ 2 ＝φ 2 - θ; n=1, 2, …,2n+1; m=1, 2, …,2m+1; θ is the beam-to-virtual antenna arrayAngle of arrival.
Optionally, the estimation module 502 is further configured to:
before estimating the beam arrival angle of the virtual antenna array based on the received signal model of the virtual antenna array, calculating an initial estimated value of the beam arrival angle of each group of antenna arrays based on the received signal model of each group of antenna arrays in at least two groups of uniform linear antenna arrays;
and determining the search range of the spatial super-resolution angle estimation algorithm according to the initial estimation value of the beam arrival angle of each group of antenna arrays.
Alternatively, the search range may be:
Wherein, for an initial estimate of the beam angle of arrival of the first antenna array,is the initial estimated value of the beam arrival angle of the second antenna array, and delta is a preset value.
Alternatively, Δ is inversely related to the signal-to-noise ratio.
Optionally, the estimation module 502 is further configured to:
after estimating the beam arrival angle of the virtual antenna array based on the received signal model of the virtual antenna array by using a spatial super-resolution angle estimation algorithm, updating the initial final estimated value of the beam arrival angle of each group of antenna arrays according to the beam arrival angle of the virtual antenna array to obtain the final estimated value of the beam arrival angle of each group of antenna arrays.
Referring to fig. 6, based on the same technical concept, an embodiment of the present application further provides a communication apparatus including a processor 601 and a memory 602; memory 602 is used to store computer-executable instructions; the processor 601 is configured to execute computer-executable instructions stored in the memory 602 to cause the communication device to perform the method as shown in fig. 2.
The processor 601 and the memory 602 may be coupled through an interface circuit, or may be integrated together, which is not limited herein.
The specific connection medium between the processor 601 and the memory 602 is not limited in the embodiment of the present application. In the embodiment of the present application, the processor 601 and the memory 602 are connected through a bus, the bus is shown by a thick line in fig. 6, and the connection manner between other components is only schematically illustrated, but not limited to. The buses may be classified as address buses, data buses, control buses, etc. For ease of illustration, only one thick line is shown in fig. 6, but not only one bus or one type of bus.
It should be understood that the processors mentioned in the embodiments of the present application may be implemented by hardware or may be implemented by software. When implemented in hardware, the processor may be a logic circuit, an integrated circuit, or the like. When implemented in software, the processor may be a general purpose processor, implemented by reading software code stored in a memory.
By way of example, the processor may be a central processing unit (Central Processing Unit, CPU), but may also be other general purpose processors, digital signal processors (Digital Signal Processor, DSP), application specific integrated circuits (Application Specific Integrated Circuit, ASIC), off-the-shelf programmable gate arrays (Field Programmable Gate Array, FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, or the like. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.
It should be understood that the memory referred to in embodiments of the present application may be volatile memory or nonvolatile memory, or may include both volatile and nonvolatile memory. The nonvolatile Memory may be a Read-Only Memory (ROM), a Programmable ROM (PROM), an Erasable PROM (EPROM), an Electrically Erasable EPROM (EEPROM), or a flash Memory. The volatile memory may be random access memory (Random Access Memory, RAM) which acts as an external cache. By way of example, and not limitation, many forms of RAM are available, such as Static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double Data rate Synchronous DRAM (DDR SDRAM), enhanced SDRAM (ESDRAM), synchronous Link DRAM (SLDRAM), and Direct RAM (DR RAM).
It should be noted that when the processor is a general purpose processor, DSP, ASIC, FPGA or other programmable logic device, discrete gate or transistor logic device, discrete hardware components, the memory (storage module) may be integrated into the processor.
It should be noted that the memory described herein is intended to comprise, without being limited to, these and any other suitable types of memory.
Based on the same technical concept, the embodiment of the application also provides a communication device, which comprises a processor and an interface circuit; the interface circuit is used for receiving the code instruction and transmitting the code instruction to the processor; the processor executes code instructions to perform the method as shown in fig. 2.
Based on the same technical concept, the embodiment of the present application also provides a computer readable storage medium for storing instructions that when executed cause the method shown in fig. 2 to be implemented.
Based on the same technical concept, the embodiment of the application further provides a chip, wherein the chip is coupled with the memory and is used for reading and executing the program instructions stored in the memory to realize the method shown in fig. 2.
Based on the same technical idea, an embodiment of the present application also provides a computer program product comprising instructions, which, when run on a computer, cause the computer to perform the method as shown in fig. 2, are stored in the computer program product.
It will be appreciated by those skilled in the art that embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to the application. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
It will be apparent to those skilled in the art that various modifications and variations can be made to the present application without departing from the scope of the application. Thus, it is intended that the present application also include such modifications and alterations insofar as they come within the scope of the appended claims or the equivalents thereof.
Claims (29)
- A method of estimating an angle of arrival of a beam, the method comprising:constructing a group of non-uniform linear virtual antenna arrays based on at least two groups of uniform linear antenna arrays, wherein the virtual antenna arrays meet far field hypothesis conditions; establishing a receiving signal model of the virtual antenna array according to the relative position information between each group of antenna arrays in the at least two groups of uniform linear antenna arrays and the receiving signal model of each group of antenna arrays;and estimating the wave beam arrival angle of the virtual antenna array based on the received signal model of the virtual antenna array by adopting a spatial super-resolution angle estimation algorithm.
- The method of claim 1, wherein building the received signal model of the virtual antenna array based on the relative position information between each of the at least two sets of uniformly linear antenna arrays and the received signal model of each of the sets of antenna arrays comprises:determining the beam arrival angle relation and the phase relation of each group of antenna arrays according to the relative position information of the at least two groups of uniform linear antenna arrays;rewriting a received signal model of each antenna array according to the phase relation and the beam arrival angle relation, enabling a dependent variable in the received signal model of each antenna array to be represented by the beam arrival angle of the virtual antenna array, and enabling a received vector of each array element in the received signal model of each antenna array to be represented by taking a received vector of a first array element as a reference, wherein the first array element is any array element in the at least two groups of uniformly linear antenna arrays;And determining the receiving signal model of the virtual antenna array based on the rewritten receiving signal models of the antenna arrays of each group.
- The method of claim 2, wherein the at least two sets of uniform linear antenna arrays consist of a first antenna array and a second antenna array;the relative position information of the at least two groups of uniformly linear antenna arrays includes: the array element interval d in the first antenna array and the second antenna array is half wavelength lambda/2, wherein lambda is wavelength, the number of array elements in the first antenna array is 2N+1, the number of array elements in the second antenna array is 2M+1, and the center point O of the first antenna array 1 And the center point O of the second antenna array 2 Distance of K (lambda/2), O 1 With O 2 The included angle between the connection line of the first antenna array and the radial direction is phi 1 ，O 1 With O 2 The included angle between the connecting line of the second antenna array and the radial direction is phi 2 ；The received signal model of the first antenna array taking the first array element as the reference point in Q symbol periods is as follows:Y 1 ＝A 1 S+W 1 ；the second antenna array uses the first array element as the reference point in the first q symbol periods to receive the signal model as follows:Y 2 ＝A 2 S+W 2 ；wherein Y is 1 ＝[Y 11 ,Y 12 ,…,Y 1Q ]，Y 2 ＝[Y 21 ,Y 22 ,…,Y 2Q ,]Is a receiving matrix;Y 1q ＝[y 11 (q),y 12 (q),…,y 1,2N+1 (q)] T ,Y 2q ＝[y 21 (q),y 22 (q),…,y 2,2M+1 (q)] T For the received vector in the q-th symbol period, T represents a transpose operation;is a guide vector; wherein θ is 1 An angle θ for the beam arrival of the first antenna array 2 An angle is reached for a beam of the second antenna array;S＝[S 1 ,S 2 ,…,S Q ]training sequence vectors formed by Q pilot symbols;W 1 ＝[W 11 W 12 ,…,W 1Q ]，W 2 ＝[W 21 W 22 ,…,W 2Q ]is a noise matrix;W 1q ＝[W 11 (q),W 12 (q),…,W 1,2N+1 (q)] T ；W 2q ＝[W 21 (q),W 22 (q),…,W 2,2M+1 (q)] T ；W 1q 、W 2q a noise vector for the q-th symbol period;the array elements of the virtual antenna array are the array elements of the first antenna array and the second antenna array, and the center of the virtual antenna array is O 1 、O 2 A midpoint of the connection line;the beam arrival angle relation of each group of antenna arrays determined according to the relative position information among each group of antenna arrays comprises:θ 1 ＝φ 1 +θ； (1)θ 2 ＝φ 2 -θ； (2)wherein θ is a beam arrival angle of the virtual antenna array;rewriting a received signal model of each group of antenna arrays according to the phase relation and the beam arrival angle relation, including:the receiving signals taking the first array element of the first antenna array as a reference point in the q-th symbol period are as follows:at the center point O of the first antenna array 1 The position is a reference point, and according to formula (1), the received signal of the nth array element of the first antenna array in the qth symbol period is rewritten as follows:Where n=1, 2, …,2n+1;the receiving signals taking the first array element of the second antenna array as a reference point in the q-th symbol period of the m-th array element of the second antenna array are as follows:with the center point O of the second antenna array 2 The position is the reference point, and the received signal of the m-th array element of the second antenna array in the q-th symbol period is rewritten as follows according to the formula (2):with the central point O of the first antenna element 1 The position is a reference point according to the formulas (1) (2) and O 2 With O 1 And (2) rewrites the received signal of the mth element of the second antenna array in the qth symbol period to:wherein m=1, 2, …,2m+1;determining the received signal model of the virtual antenna array based on the rewritten received signal models of the antenna arrays of each group, comprising:according to formulas (4) and (7), stacking the received signals of the virtual antenna array into a matrix of (2N+2M+2) multiplied by Q dimension, and obtaining a received signal model of the virtual antenna array as follows:Y＝[y(1)，y(2)，…，y(Q)]＝AS+W；wherein,
- the method of claim 1 or 2, wherein the at least two sets of uniform linear antenna arrays consist of a first antenna array and a second antenna array;The relative position information of the at least two groups of uniformly linear antenna arrays includes: the array element interval d in the first antenna array and the second antenna array is half wavelength lambda/2, wherein lambda is wavelength, the number of array elements in the first antenna array is 2N+1, the number of array elements in the second antenna array is 2M+1, and the center point O of the first antenna array 1 And the center point O of the second antenna array 2 Distance of K (lambda/2), O 1 With O 2 The included angle between the connection line of the first antenna array and the radial direction is phi 1 ，O 1 With O 2 The included angle between the connecting line of the second antenna array and the radial direction is phi 2 ；The received signal model of the first antenna array in the first q symbol periods is:Y 1 ＝A 1 S+W 1 ；the received signal model of the second antenna array in the first q symbol periods is:Y 2 ＝A 2 S+W 2 ；wherein Y is 1 ＝[Y 11 ,Y 12 ,…,Y 1Q ]，Y 2 ＝[Y 21 ,Y 22 ,…,Y 2Q ,]Is a receiving matrix;Y 1q ＝[y 11 (q),y 12 (q),…,y 1,2N+1 (q)] T ,Y 2q ＝[y 21 (q),y 22 (q),…,y 2,2M+1 (q)] T for the received vector in the q-th symbol period, T represents a transpose operation;is a guide vector; wherein θ is 1 An angle θ for the beam arrival of the first antenna array 2 An angle is reached for a beam of the second antenna array;S＝[S 1 ,S 2 ,…,S Q ]training sequence vectors formed by Q pilot symbols;W 1 ＝[W 11 W 12 ,…,W 1Q ]，W 2 ＝[W 21 W 22 ,…,W 2Q ]is a noise matrix;W 1q ＝[W 11 (q),W 12 (q),…,W 1,2N+1 (q)] T ；W 2q ＝[W 21 (q),W 22 (q),…,W 2,2M+1 (q)] T ；W 1q 、W 2q a noise vector for the q-th symbol period;The array elements of the virtual antenna array are the array elements of the first antenna array and the second antenna array, and the center of the virtual antenna array is O 1 、O 2 A midpoint of the connection line;the received signal model of the virtual antenna array is as follows:Y＝[y(1)，y(2)，…，y(Q)]＝AS+W；wherein,each array element in the virtual antenna array uses the center point O of the first antenna array in the q-th symbol period 1 The received signal for the reference point satisfies the following relationship:wherein θ 1 ＝φ 1 +θ，θ 2 ＝φ 2 - θ; n=1, 2, …,2n+1; m=1, 2, …,2m+1; θ is the beam arrival angle of the virtual antenna array.
- The method of any of claims 1-4, wherein prior to estimating the beam angle of arrival of the virtual antenna array based on the received signal model of the virtual antenna array, the method further comprises:calculating an initial estimated value of a beam arrival angle of each group of antenna arrays based on a received signal model of each group of antenna arrays in the at least two groups of uniform linear antenna arrays;and determining the search range of the spatial super-resolution angle estimation algorithm according to the initial estimation value of the beam arrival angle of each group of antenna arrays.
- The method of claim 5, wherein the search range is:Wherein, for an initial estimate of the beam angle of arrival of the first antenna array,and delta is a preset value which is an initial estimated value of the beam arrival angle of the second antenna array.
- The method of claim 6, wherein Δ is inversely related to signal-to-noise ratio.
- The method of any of claims 1-7, wherein after estimating the beam angle of arrival of the virtual antenna array based on a received signal model of the virtual antenna array using a spatial super-resolution angle estimation algorithm, the method further comprises:and updating the initial final estimated value of the beam arrival angle of each group of antenna arrays according to the beam arrival angle of the virtual antenna arrays to obtain the final estimated value of the beam arrival angle of each group of antenna arrays.
- An apparatus for estimating an angle of arrival of a beam, comprising:a construction module for constructing a set of non-uniform linear virtual antenna arrays based on at least two sets of uniform linear antenna arrays; establishing a receiving signal model of the virtual antenna array according to the relative position information between each group of antenna arrays in the at least two groups of uniform linear antenna arrays and the receiving signal model of each group of antenna arrays;The estimation module is used for estimating the wave beam arrival angle of the virtual antenna array based on the received signal model of the virtual antenna array by adopting a spatial super-resolution angle estimation algorithm.
- The apparatus of claim 9, wherein the building module, when building the received signal model of the virtual antenna array based on the relative position information between each of the at least two groups of uniformly linear antenna arrays and the received signal models of each of the groups of antenna arrays, is specifically configured to:determining the beam arrival angle relation and the phase relation of each group of antenna arrays according to the relative position information of the at least two groups of uniform linear antenna arrays;rewriting a received signal model of each antenna array according to the phase relation and the beam arrival angle relation, enabling a dependent variable in the received signal model of each antenna array to be represented by the beam arrival angle of the virtual antenna array, and enabling a received vector of each array element in the received signal model of each antenna array to be represented by taking a received vector of a first array element as a reference, wherein the first array element is any array element in the at least two groups of uniformly linear antenna arrays;And determining the receiving signal model of the virtual antenna array based on the rewritten receiving signal models of the antenna arrays of each group.
- The apparatus of claim 10, wherein the at least two sets of uniform linear antenna arrays consist of a first antenna array and a second antenna array;the relative position information of the at least two groups of uniformly linear antenna arrays includes: the array element interval d in the first antenna array and the second antenna array is half wavelength lambda/2, wherein lambda is wavelength, the number of array elements in the first antenna array is 2N+1, the number of array elements in the second antenna array is 2M+1, and the center point O of the first antenna array 1 And the center point O of the second antenna array 2 Distance of K (lambda/2), O 1 With O 2 The included angle between the connection line of the first antenna array and the radial direction is phi 1 ，O 1 With O 2 The included angle between the connecting line of the second antenna array and the radial direction is phi 2 ；The received signal model of the first antenna array taking the first array element as the reference point in Q symbol periods is as follows:Y 1 ＝A 1 S+W 1 ；the second antenna array uses the first array element as the reference point in the first q symbol periods to receive the signal model as follows:Y 2 ＝A 2 S+W 2 ；wherein Y is 1 ＝[y 11 ,y 12 ,…,y 1Q ]，Y 2 ＝[y 21 ,y 22 ,…,y 2Q ,]Is a receiving matrix;Y 1q ＝[y 11 (q),y 12 (q),…,y 1,2N+1 (q)],Y 2q ＝[y 21 (q),y 22 (q),…,y 2,2M+1 (q)]Is the received vector in the q-th symbol period;is a guide vector; wherein θ is 1 An angle θ for the beam arrival of the first antenna array 2 An angle is reached for a beam of the second antenna array;S＝[S 1 ,S 2 ,…,S Q ]training sequence vectors formed by Q pilot symbols;W 1 ＝[W 11 W 12 ,…,W 1Q ]，W 2 ＝[W 21 W 22 ,…,W 2Q ]is a noise matrix;W 1q ＝[W 11 (q),W 12 (q),…,W 1,2N+1 (q)]；W 2q ＝[W 21 (q),W 22 (q),…,W 2,2M+1 (q)]；W 1q 、W 2q a noise vector for the q-th symbol period;the array elements of the virtual antenna array are the array elements of the first antenna array and the second antenna array, and the center of the virtual antenna array is O 1 、O 2 A midpoint of the connection line;the construction module determines the beam arrival angle relation of each group of antenna arrays according to the relative position information among the groups of antenna arrays, and comprises the following steps:θ 1 ＝ф 1 +θ； (1)θ 2 ＝ф 2 -θ； (2)wherein θ is a beam arrival angle of the virtual antenna array;the construction module is specifically configured to, when rewriting the received signal model of each group of antenna arrays according to the phase relationship and the beam arrival angle relationship:the receiving signals taking the first array element of the first antenna array as a reference point in the q-th symbol period are as follows:at the center point O of the first antenna array 1 The position is a reference point, and according to formula (1), the received signal of the nth array element of the first antenna array in the qth symbol period is rewritten as follows:Where n=1, 2, …,2n+1;the receiving signals taking the first array element of the second antenna array as a reference point in the q-th symbol period of the m-th array element of the second antenna array are as follows:with the center point O of the second antenna array 2 The position is the reference point, and the received signal of the m-th array element of the second antenna array in the q-th symbol period is rewritten as follows according to the formula (2):with the central point O of the first antenna element 1 The position is a reference point according to the formulas (1) (2) and O 2 With O 1 And (2) rewrites the received signal of the mth element of the second antenna array in the qth symbol period to:wherein m=1, 2, …,2m+1;the construction module is specifically configured to, when determining the received signal model of the virtual antenna array based on the rewritten received signal models of the antenna arrays of each group:according to formulas (4) and (7), stacking the received signals of the virtual antenna array into a matrix of (2N+2M+2) multiplied by Q dimension, and obtaining a received signal model of the virtual antenna array as follows:Y＝[y(1)，y(2)，…，y(Q)]＝AS+W；wherein,
- the apparatus of claim 9 or 10, wherein the at least two sets of uniform linear antenna arrays consist of a first antenna array and a second antenna array;The relative position information of the at least two groups of uniformly linear antenna arrays includes: the array element interval d in the first antenna array and the second antenna array is half wavelength lambda/2, wherein lambda is wavelength, the number of array elements in the first antenna array is 2N+1, the number of array elements in the second antenna array is 2M+1, and the center point O of the first antenna array 1 With the second antenna arrayCenter point O 2 Distance of K (lambda/2), O 1 With O 2 The included angle between the connection line of the first antenna array and the radial direction is phi 1 ，O 1 With O 2 The included angle between the connecting line of the second antenna array and the radial direction is phi 2 ；The received signal model of the first antenna array in the first q symbol periods is:Y 1 ＝A 1 S+W 1 ；the received signal model of the second antenna array in the first q symbol periods is:Y 2 ＝A 2 S+W 2 ；wherein Y is 1 ＝[y 11 ,y 12 ,…,y 1Q ]，Y 2 ＝[y 21 ,y 22 ,…,y 2Q ,]Is a receiving matrix;Y 1q ＝[y 11 (q),y 12 (q),…,y 1,2N+1 (q)],Y 2q ＝[y 21 (q),y 22 (q),…,y 2,2M+1 (q)]is the received vector in the q-th symbol period;is a guide vector; wherein θ is 1 An angle θ for the beam arrival of the first antenna array 2 An angle is reached for a beam of the second antenna array;S＝[S 1 ,S 2 ,…,S Q ]training sequence vectors formed by Q pilot symbols;W 1 ＝[W 11 W 12 ,…,W 1Q ]，W 2 ＝[W 21 W 22 ,…,W 2Q ]is a noise matrix;W 1q ＝[W 11 (q),W 12 (q),…,W 1,2N+1 (q)]；W 2q ＝[W 21 (q),W 22 (q),…,W 2,2M+1 (q)]；W 1q 、W 2q a noise vector for the q-th symbol period;the array elements of the virtual antenna array are the array elements of the first antenna array and the second antenna array, and the center of the virtual antenna array is O 1 、O 2 A midpoint of the connection line;the received signal model of the virtual antenna array is as follows:Y＝[y(1)，y(2)，…，y(Q)]＝AS+W；wherein,each array element in the virtual antenna array uses the center point O of the first antenna array in the q-th symbol period 1 The received signal for the reference point satisfies the following relationship:wherein θ 1 ＝ф 1 +θ，θ 2 ＝ф 2 - θ; n=1, 2, …,2n+1; m=1, 2, …,2m+1; θ is the virtualAngle of beam arrival of the antenna array.
- The apparatus of any of claims 9-12, wherein the estimation module is further to:before estimating the beam arrival angle of the virtual antenna array based on the received signal model of the virtual antenna array, calculating an initial estimated value of the beam arrival angle of each group of antenna arrays based on the received signal model of each group of antenna arrays in the at least two groups of uniformly linear antenna arrays;and determining the search range of the spatial super-resolution angle estimation algorithm according to the initial estimation value of the beam arrival angle of each group of antenna arrays.
- The apparatus of claim 13, wherein the search range is:wherein, for an initial estimate of the beam angle of arrival of the first antenna array,and delta is a preset value which is an initial estimated value of the beam arrival angle of the second antenna array.
- The apparatus of claim 14, wherein Δ is inversely related to signal-to-noise ratio.
- The apparatus of any of claims 9-15, wherein the estimation module is further to:and after estimating the beam arrival angle of the virtual antenna array based on the received signal model of the virtual antenna array by using a spatial super-resolution angle estimation algorithm, updating the initial final estimated value of the beam arrival angle of each group of antenna arrays according to the beam arrival angle of the virtual antenna array to obtain the final estimated value of the beam arrival angle of each group of antenna arrays.
- An antenna system, comprising:each remote radio unit comprises a group of uniform linear antenna arrays, and the remote radio units are used for receiving wireless signals;the baseband processing unit is in communication connection with the at least two remote radio units and is used for: constructing a group of non-uniform linear virtual antenna arrays based on at least two groups of uniform linear antenna arrays, wherein the virtual antenna arrays meet far field hypothesis conditions; establishing a receiving signal model of the virtual antenna array according to the relative position information between each group of antenna arrays in the at least two groups of uniform linear antenna arrays and the receiving signal model of each group of antenna arrays; and estimating the wave beam arrival angle of the virtual antenna array based on the received signal model of the virtual antenna array by adopting a spatial super-resolution angle estimation algorithm.
- The system of claim 17, wherein the baseband processing unit is configured to, when establishing the received signal model of the virtual antenna array based on the relative position information between each of the at least two groups of uniformly linear antenna arrays and the received signal models of each of the groups of antenna arrays:determining the beam arrival angle relation and the phase relation of each group of antenna arrays according to the relative position information of the at least two groups of uniform linear antenna arrays;rewriting a received signal model of each antenna array according to the phase relation and the beam arrival angle relation, enabling a dependent variable in the received signal model of each antenna array to be represented by the beam arrival angle of the virtual antenna array, and enabling a received vector of each array element in the received signal model of each antenna array to be represented by taking a received vector of a first array element as a reference, wherein the first array element is any array element in the at least two groups of uniformly linear antenna arrays;and determining the receiving signal model of the virtual antenna array based on the rewritten receiving signal models of the antenna arrays of each group.
- The system of claim 18, wherein the at least two sets of uniform linear antenna arrays consist of a first antenna array and a second antenna array;the relative position information of the at least two groups of uniformly linear antenna arrays includes: the array element interval d in the first antenna array and the second antenna array is half wavelength lambda/2, wherein lambda is wavelength, the number of array elements in the first antenna array is 2N+1, the number of array elements in the second antenna array is 2M+1, and the center point O of the first antenna array 1 And the center point O of the second antenna array 2 Distance of K (lambda/2), O 1 With O 2 The included angle between the connection line of the first antenna array and the radial direction is phi 1 ，O 1 With O 2 The included angle between the connecting line of the second antenna array and the radial direction is phi 2 ；The received signal model of the first antenna array taking the first array element as the reference point in Q symbol periods is as follows:Y 1 ＝A 1 S+W 1 ；the second antenna array uses the first array element as the reference point in the first q symbol periods to receive the signal model as follows:Y 2 ＝A 2 S+W 2 ；wherein Y is 1 ＝[y 11 ,y 12 ,…,y 1Q ]，Y 2 ＝[y 21 ,y 22 ,…,y 2Q ,]Is a receiving matrix;Y 1q ＝[y 11 (q),y 12 (q),…,y 1,2N+1 (q)],Y 2q ＝[y 21 (q),y 22 (q),…,y 2,2M+1 (q)]is the received vector in the q-th symbol period;is a guide vector; wherein θ is 1 An angle θ for the beam arrival of the first antenna array 2 An angle is reached for a beam of the second antenna array;S＝[S 1 ,S 2 ,…,S Q ]training sequence vectors formed by Q pilot symbols;W 1 ＝[W 11 W 12 ,…,W 1Q ]，W 2 ＝[W 21 W 22 ,…,W 2Q ]is a noise matrix;W 1q ＝[W 11 (q),W 12 (q),…,W 1,2N+1 (q)]；W 2q ＝[W 21 (q),W 22 (q),…,W 2,2M+1 (q)]；W 1q 、W 2q a noise vector for the q-th symbol period;the array elements of the virtual antenna array are the array elements of the first antenna array and the second antenna array, and the center of the virtual antenna array is O 1 、O 2 A midpoint of the connection line;the beam arrival angle relation of each group of antenna arrays determined according to the relative position information among each group of antenna arrays comprises:θ 1 ＝ф 1 +θ； (1)θ 2 ＝ф 2 -θ； (2)wherein θ is a beam arrival angle of the virtual antenna array;the baseband processing unit is specifically configured to, when rewriting the received signal model of each group of antenna arrays according to the phase relationship and the beam arrival angle relationship:the receiving signals taking the first array element of the first antenna array as a reference point in the q-th symbol period are as follows:at the center point O of the first antenna array 1 The position is a reference point, and according to formula (1), the received signal of the nth array element of the first antenna array in the qth symbol period is rewritten as follows:where n=1, 2, …,2n+1;the receiving signals taking the first array element of the second antenna array as a reference point in the q-th symbol period of the m-th array element of the second antenna array are as follows:With the center point O of the second antenna array 2 The position is used as a reference point, and the m-th array element of the second antenna array is received in the q-th symbol period according to the formula (2)The signal is rewritten as:with the central point O of the first antenna element 1 The position is a reference point according to the formulas (1) (2) and O 2 With O 1 And (2) rewrites the received signal of the mth element of the second antenna array in the qth symbol period to:wherein m=1, 2, …,2m+1;the baseband processing unit is specifically configured to, when determining the received signal model of the virtual antenna array based on the rewritten received signal models of the antenna arrays of each group:according to formulas (4) and (7), stacking the received signals of the virtual antenna array into a matrix of (2N+2M+2) multiplied by Q dimension, and obtaining a received signal model of the virtual antenna array as follows:Y＝[y(1)，y(2)，…，y(Q)]＝AS+W；wherein,
- the system of claim 17 or 18, wherein the at least two sets of uniform linear antenna arrays consist of a first antenna array and a second antenna array;the relative position information of the at least two groups of uniformly linear antenna arrays includes: the array element interval d in the first antenna array and the second antenna array is half wavelength lambda/2, wherein lambda is wavelength, the number of array elements in the first antenna array is 2N+1, and the number of array elements in the second antenna array is 2N+1 The number of array elements of the second antenna array is 2M+1, and the center point O of the first antenna array 1 And the center point O of the second antenna array 2 Distance of K (lambda/2), O 1 With O 2 The included angle between the connection line of the first antenna array and the radial direction is phi 1 ，O 1 With O 2 The included angle between the connecting line of the second antenna array and the radial direction is phi 2 ；The received signal model of the first antenna array in the first q symbol periods is:Y 1 ＝A 1 S+W 1 ；the received signal model of the second antenna array in the first q symbol periods is:Y 2 ＝A 2 S+W 2 ；wherein Y is 1 ＝[y 11 ,y 12 ,…,y 1Q ]，Y 2 ＝[y 21 ,y 22 ,…,y 2Q ,]Is a receiving matrix;Y 1q ＝[y 11 (q),y 12 (q),…,y 1,2N+1 (q)],Y 2q ＝[y 21 (q),y 22 (q),…,y 2,2M+1 (q)]is the received vector in the q-th symbol period;is a guide vector; wherein θ is 1 An angle θ for the beam arrival of the first antenna array 2 An angle is reached for a beam of the second antenna array;S＝[S 1 ,S 2 ,…,S Q ]training sequence vectors formed by Q pilot symbols;W 1 ＝[W 11 W 12 ,…,W 1Q ]，W 2 ＝[W 21 W 22 ,…,W 2Q ]is a noise matrix;W 1q ＝[W 11 (q),W 12 (q),…,W 1,2N+1 (q)]；W 2q ＝[W 21 (q),W 22 (q),…,W 2,2M+1 (q)]；W 1q 、W 2q a noise vector for the q-th symbol period;the array elements of the virtual antenna array are the array elements of the first antenna array and the second antenna array, and the center of the virtual antenna array is O 1 、O 2 A midpoint of the connection line;the received signal model of the virtual antenna array is as follows:Y＝[y(1)，y(2)，…，y(Q)]＝AS+W；wherein,each array element in the virtual antenna array uses the center point O of the first antenna array in the q-th symbol period 1 The received signal for the reference point satisfies the following relationship:wherein θ 1 ＝ф 1 +θ，θ 2 ＝ф 2 - θ; n=1, 2, …,2n+1; m=1, 2, …,2m+1; θ is the beam arrival angle of the virtual antenna array.
- The system of any of claims 17-20, wherein the baseband processing unit is further configured to:before estimating the beam arrival angle of the virtual antenna array based on the received signal model of the virtual antenna array, calculating an initial estimated value of the beam arrival angle of each group of antenna arrays based on the received signal model of each group of antenna arrays in the at least two groups of uniformly linear antenna arrays;and determining the search range of the spatial super-resolution angle estimation algorithm according to the initial estimation value of the beam arrival angle of each group of antenna arrays.
- The system of claim 21, wherein the search range is:wherein, for an initial estimate of the beam angle of arrival of the first antenna array,and delta is a preset value which is an initial estimated value of the beam arrival angle of the second antenna array.
- The system of claim 22, wherein Δ is inversely related to signal-to-noise ratio.
- The system of any of claims 17-23, wherein the baseband processing unit is further configured to:And after estimating the beam arrival angle of the virtual antenna array based on the received signal model of the virtual antenna array by using a spatial super-resolution angle estimation algorithm, updating the initial final estimated value of the beam arrival angle of each group of antenna arrays according to the beam arrival angle of the virtual antenna array to obtain the final estimated value of the beam arrival angle of each group of antenna arrays.
- A communication device comprising a processor and a memory; the memory is used for storing computer execution instructions; the processor is configured to execute computer-executable instructions stored in the memory to cause the communication device to perform the method of any one of claims 1 to 8.
- A communication device comprising a processor and an interface circuit; the interface circuit is used for receiving code instructions and transmitting the code instructions to the processor; the processor executes the code instructions to perform the method of any one of claims 1 to 8.
- A computer readable storage medium for storing instructions that, when executed, cause the method of any one of claims 1 to 8 to be implemented.
- A chip, characterized in that the chip is coupled to a memory for reading and executing program instructions stored in the memory, implementing the method according to any of claims 1 to 8.
- A computer program product comprising instructions stored therein, which when run on a computer, cause the computer to perform the method of any of claims 1 to 8.
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PCT/CN2021/078766 WO2022183382A1 (en) | 2021-03-02 | 2021-03-02 | Method and apparatus for estimating angle of arrival of beam, and antenna system |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116918274A true CN116918274A (en) | 2023-10-20 |
Family
ID=83153803
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202180094736.8A Pending CN116918274A (en) | 2021-03-02 | 2021-03-02 | Method, device and antenna system for estimating wave beam arrival angle |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN116918274A (en) |
| WO (1) | WO2022183382A1 (en) |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105912791B (en) * | 2016-04-15 | 2019-03-29 | 山东农业大学 | DOA estimation method based on local search in virtual relatively prime array |
| CN106383335B (en) * | 2016-08-18 | 2019-03-29 | 成都中安频谱科技有限公司 | The localization method and device in the fixed direction finding site near-field interference source of high frequency radio |
| EP3418761A1 (en) * | 2017-06-21 | 2018-12-26 | Technische Universität Graz | Method and system for determining a position |
| CN107422295B (en) * | 2017-08-30 | 2019-09-13 | 浙江大学 | The Wave arrival direction estimating method indicated based on relatively prime array virtual Domain equivalent signal atom norm |
-
2021
- 2021-03-02 CN CN202180094736.8A patent/CN116918274A/en active Pending
- 2021-03-02 WO PCT/CN2021/078766 patent/WO2022183382A1/en active Application Filing
Also Published As
| Publication number | Publication date |
|---|---|
| WO2022183382A1 (en) | 2022-09-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10355762B2 (en) | System and method for adaptive beamforming communication | |
| WO2016202280A1 (en) | Beamforming method and device | |
| CN114900400A (en) | Joint sparse channel estimation method based on intelligent reflector assisted Internet of things | |
| US9647739B2 (en) | Method and device for acquiring channel information | |
| WO2019206157A1 (en) | Method for training downlink beam, network device and terminal device | |
| CN110034804B (en) | Method and apparatus for estimating angle information for wireless communication system | |
| CN111865370A (en) | Method, device and system for determining arrival angle of signal | |
| CN108075811A (en) | For mixing the method for precoding and communication equipment | |
| US9762307B2 (en) | Communication system and communication control method using interference channel matrix | |
| CN114600384A (en) | Channel measurement method and communication device | |
| CN111641440B (en) | Method, device, terminal and storage medium for determining precoding matrix | |
| CN116918274A (en) | Method, device and antenna system for estimating wave beam arrival angle | |
| CN110350990B (en) | Phased array network calibration method, device, equipment and storage medium | |
| WO2018068365A1 (en) | Doppler frequency offset estimation method and device based on millimeter wave mimo system | |
| CN116722896A (en) | Multi-user precoding method, equipment and medium for super-directivity antenna array | |
| CN113938168B (en) | Split beam management method and system for uplink and downlink asymmetric communication MIMO system | |
| US11271615B2 (en) | Communications method and apparatus | |
| CN112020889B (en) | Communication method and device | |
| KR101244714B1 (en) | Smart Antenna System for Rejection of Coherent and Incoherent Interferences | |
| CN110417692B (en) | Uplink channel tracking method and device | |
| JP6491339B2 (en) | Directional beam based open loop multi-stream transmission method and base station | |
| CN112272151B (en) | Channel estimation method and device | |
| WO2022266868A1 (en) | Beam tracking method and related device | |
| Senanayake et al. | Distributed LMMSE estimation in cooperative cellular networks | |
| CN105229933B (en) | Space division multiplexing method and base station |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |