WO2015096099A1 - Procédé et appareil d'estimation d'un angle d'arrivée, et dispositif électronique - Google Patents

Procédé et appareil d'estimation d'un angle d'arrivée, et dispositif électronique Download PDF

Info

Publication number
WO2015096099A1
WO2015096099A1 PCT/CN2013/090581 CN2013090581W WO2015096099A1 WO 2015096099 A1 WO2015096099 A1 WO 2015096099A1 CN 2013090581 W CN2013090581 W CN 2013090581W WO 2015096099 A1 WO2015096099 A1 WO 2015096099A1
Authority
WO
WIPO (PCT)
Prior art keywords
matrix
antenna array
arrival
redundant dictionary
measurement results
Prior art date
Application number
PCT/CN2013/090581
Other languages
English (en)
Chinese (zh)
Inventor
刘劲楠
王悦
Original Assignee
华为技术有限公司
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 华为技术有限公司 filed Critical 华为技术有限公司
Priority to PCT/CN2013/090581 priority Critical patent/WO2015096099A1/fr
Priority to CN201380002828.4A priority patent/CN104937856B/zh
Publication of WO2015096099A1 publication Critical patent/WO2015096099A1/fr

Links

Classifications

    • 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
    • H04B7/0837—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station using pre-detection combining
    • H04B7/0842—Weighted combining
    • H04B7/086—Weighted combining using weights depending on external parameters, e.g. direction of arrival [DOA], predetermined weights or beamforming

Definitions

  • the present invention relates to the field of communications, and in particular, to a method, apparatus, and electronic device for estimating an angle of arrival. Background technique
  • the base station can measure the terminal by measuring the DOA of the terminal signal (ie, the uplink signal), and can also perform the weight selection of the downlink beamforming (Downlink Beam Forming) according to the DO A.
  • the terminal can also perform weight selection of Uplink Beam Forming by measuring the DOA of the base station signal (ie, the downlink signal).
  • the invention provides a method, a device and an electronic device for estimating an angle of arrival, which are used to solve the technical problem that a storage redundancy dictionary existing in the prior art requires a large storage space.
  • the present invention provides an electronic device, which is applied to an Orthogonal Frequency Division Multiplexing (OFDM) communication system, and includes: a processor, configured to obtain an antenna array measurement result of M array elements; and based on the antenna array measurement result and The redundancy dictionary estimates the angle of arrival of the K wireless signals, wherein the redundancy dictionary is a partial discrete Fourier transform matrix, where M and K are both positive integers.
  • OFDM Orthogonal Frequency Division Multiplexing
  • the processor is specifically configured to: obtain the M array elements respectively An antenna array measurement result of the array elements of each dimension; for each of the array elements, estimating K arrival angles based on the antenna array measurement results and the redundancy dictionary; The K arrival angles of the array elements of the dimensions, and the K arrival angles of the M array elements are estimated.
  • the redundancy dictionary is obtained by: uniformly dividing a range of trigonometric function values into N grids; A direction vector of each of the N grids, and the redundant dictionary is composed of N of the direction vectors.
  • the direction vector is determined according to the following formula:
  • the processor is specifically configured to:
  • the measurement result and the redundancy dictionary are obtained by using the compressed sensing algorithm to obtain a sparse signal matrix, wherein the compressed sensing algorithm includes at least: ⁇ calculating a transposed right multiplied first matrix of the redundant dictionary by using a fast Fourier algorithm Operation and/or redundancy dictionary right multiplication of the second matrix operation, the first matrix and the second matrix being a matrix associated with the antenna array measurement results; estimating the angle of arrival based on the sparse signal matrix .
  • the processor is further configured to:: based on the antenna array measurement result, the redundancy dictionary, and one of the redundancy dictionary The order derivative matrix is obtained by using the compressed sensing algorithm to obtain a sparse signal matrix.
  • the elements of the first derivative matrix are determined according to the following formula:
  • the processor is specifically configured to: expand the first matrix into N Dimension matrix; performing an inverse fast Fourier transform on the N-dimensional matrix and performing fast Fourier shift.
  • the processor is specifically configured to: perform fast fast on the second matrix Inverse transform, obtaining a third matrix; determining an even term in the first M term in the third matrix as an even term of the operation result of the redundant dictionary right multiplied by the second matrix, and The inverted value of the odd term in the first M term in the third matrix is determined to be an odd term of the operation result.
  • the electronic device further includes: an antenna array, configured to perform measurement according to Configuring to perform received signal measurement to obtain the antenna array measurement result, wherein the measurement configuration includes at least a combination of one or more of a sparsity degree, a snapshot number, and a number of snapshot bits.
  • the electronic device further includes: a receiver, configured to receive The antenna array measures the result, and sends the antenna array measurement result to the processor.
  • the electronic device is specifically a terminal or a base station.
  • the present invention provides an apparatus for estimating an angle of arrival, which is applied to an Orthogonal Frequency Division Multiplexing (OFDM) communication system, and includes: an antenna array measurement result receiving unit, configured to obtain an antenna array measurement result of M array elements; And an angle estimating unit, configured to estimate an angle of arrival of the K wireless signals based on the antenna array measurement result and the redundancy dictionary, wherein the redundancy dictionary is a partial discrete Fourier transform matrix, and M and K are positive integers.
  • OFDM Orthogonal Frequency Division Multiplexing
  • the angle of arrival estimating unit specifically includes: an antenna array measurement result input subunit An antenna array measurement result for respectively obtaining an array element of each of the M array elements; a single dimension angle of arrival calculation subunit, configured for the array element of each dimension, based on the antenna array measurement Results and the redundant dictionary, estimating K angles of arrival; multi-dimensional arrival An angle calculation subunit, configured to calculate the K arrival angles of the M array elements based on K arrival angles of the array elements of each of the dimensions.
  • the redundancy dictionary is obtained by: uniformly dividing a range of trigonometric function values into N grids; A direction vector of each of the N grids, and the redundant dictionary is composed of N of the direction vectors.
  • the direction vector is determined according to the following formula:
  • N is an integer power of 2.
  • the angle of arrival estimation unit specifically includes: a sparse signal matrix calculation a subunit, configured to obtain a sparse signal matrix by using the compressed sensing algorithm based on the antenna array measurement result and the redundancy dictionary, where the compressed sensing algorithm includes at least: calculating a redundancy by using a fast Fourier algorithm Transposition of the dictionary by the operation of the first matrix and/or operation of the redundancy dictionary by the operation of the second matrix.
  • the first matrix and the second matrix are matrices related to the measurement results of the antenna array; the estimation subunit And for estimating the angle of arrival based on the sparse signal matrix.
  • the sparse signal matrix calculation subunit is specifically configured to: determine, according to the antenna array measurement result, the redundancy dictionary, and the redundancy The first derivative matrix of the dictionary is obtained by using the compressed sensing algorithm to obtain a sparse signal matrix.
  • the elements in the first derivative matrix are determined according to the following formula: B( jj ii) : j 2;r (jj ⁇ ⁇ c j ⁇ ( jj -i)d-2ii N)
  • jj is the row index of the first derivative matrix
  • ii is the column index of the first derivative matrix
  • the sparse signal matrix calculation subunit includes: an extension subunit, The first matrix is expanded into a dimensional matrix; an inverse fast Fourier transform sub-unit is configured to perform an inverse fast Fourier transform on the one-dimensional matrix, and perform a fast Fourier transform shift.
  • the sparse signal matrix calculation subunit specifically includes: inverse fast Fourier transform a subunit, configured to perform an inverse fast Fourier transform on the second matrix, to obtain a third matrix, and a determining subunit, configured to determine an even term in the preceding item in the third matrix as The redundancy dictionary is right-multiplied by the even term of the operation result of the second matrix, and the inverted value of the odd term in the preceding term in the third matrix is determined as an odd term of the operation result.
  • the antenna array measurement result receiving unit is specifically configured to: The measurement configuration performs a received signal measurement to obtain the antenna array measurement result, wherein the measurement configuration includes at least a combination of one or more of a sparsity degree, a snapshot number, and a number of snapshot bits.
  • the antenna array measurement result receiving unit is specifically configured to: receive The antenna array measurement result transmitted by the peer device.
  • the present invention provides a method for estimating an angle of arrival, which is applied to an Orthogonal Frequency Division Multiplexing (OFDM) communication system, comprising: obtaining an antenna array measurement result of two array elements; based on the antenna array measurement result and redundancy a dictionary, estimating an angle of arrival of the plurality of wireless signals, wherein the redundancy dictionary is a partial discrete Fourier transform matrix, and ⁇ and ⁇ are positive integers.
  • the antenna array measurement result of obtaining the M array elements is specifically as follows: An antenna array measurement result of the array element of each of the M array elements is obtained.
  • the estimating the angle of arrival of the K wireless signals based on the antenna array measurement result and the redundancy dictionary specifically: An array element of one dimension, estimating K arrival angles based on the antenna array measurement result and the redundancy dictionary; calculating the M array elements based on K arrival angles of the array elements of each dimension Said K arrival angles.
  • the redundancy dictionary is obtained by: uniformly dividing a range of trigonometric function values into N grids; A direction vector of each of the N grids, and the redundant dictionary is composed of N of the direction vectors.
  • the direction vector is determined according to the following formula:
  • the performing, based on the antenna array measurement result and the redundancy dictionary, Estimating the angle of arrival of the K wireless signals comprising: obtaining the sparse signal matrix by using the compressed sensing algorithm based on the antenna array measurement result and the redundancy dictionary, wherein the compressed sensing algorithm includes at least: The Fourier algorithm calculates the operation of the transposed right-of-first matrix of the redundancy dictionary and/or the operation of the redundant dictionary by the second matrix. The first matrix and the second matrix are measured with the antenna array. a correlation matrix; estimating the angle of arrival based on the sparse signal matrix.
  • the obtaining, by using the compressed sensing algorithm, a sparse letter based on the antenna array measurement result and the redundancy dictionary The number matrix is specifically: based on the antenna array measurement result, the redundancy dictionary, and a first derivative matrix of the redundancy dictionary, the compressed sensing algorithm is used to obtain a sparse signal matrix.
  • the elements of the first derivative matrix are determined according to the following formula:
  • jj is the number of rows of the first derivative matrix
  • ii is the number of columns of the first derivative matrix
  • the transposed matrix of the redundant dictionary is multiplied by the first
  • the matrix specifically includes: expanding the first matrix into a one-dimensional matrix; performing an inverse fast Fourier transform on the one-dimensional matrix, and performing fast Fourier shift.
  • the method includes: performing an inverse fast Fourier transform on the second matrix to obtain a third matrix; determining an even term in the preceding term in the third matrix as the redundancy dictionary, and multiplying the redundant matrix An even term of the operation result of the two matrices, and an inverse value of the odd term in the preceding term in the third matrix is determined as an odd term of the operation result.
  • the method before the obtaining the array measurement flaw, the method further The method includes: obtaining the antenna array measurement result according to the measurement configuration, wherein the measurement configuration includes at least a combination of one or more of a sparsity degree, a snapshot number, and a number of snapshot bits.
  • the obtaining the array measurement result is specifically: receiving by the opposite end The antenna array measurement result sent by the device.
  • the present invention provides an electronic device, which is applied to an Orthogonal Frequency Division Multiplexing (OFDM) communication system
  • the electronic device includes: a processor, configured to obtain antenna array measurement results of M array elements and measure scar and redundancy based on the antenna array a dictionary, estimating an angle of arrival of the K wireless signals, wherein the redundant dictionary used herein is a partial discrete Fourier transform matrix, that is, the redundant dictionary has a discrete Fourier transform factor, then, the complex
  • the redundant dictionary is stored without consuming additional storage space, so the storage redundancy dictionary existing in the prior art is effectively solved.
  • the technical problem of large storage space saves system resources.
  • FIG. 1 is a functional block diagram of an electronic device according to an embodiment of the present invention.
  • FIGS. 2A-2B are functional block diagrams of an electronic device according to another embodiment of the present invention.
  • FIG. 3 is a schematic structural diagram of a uniform linear array according to an embodiment of the present invention.
  • FIG. 4 is a schematic structural view of an L-shaped array according to an embodiment of the present invention.
  • FIG. 5 is a schematic diagram showing a relationship between a spatial angle and a direction angle according to an embodiment of the present invention
  • FIG. 6 is a functional block diagram of an apparatus for estimating an angle of arrival in an embodiment of the present invention.
  • FIG. 7 is a flow chart of a method of estimating an angle of arrival in an embodiment of the present invention. detailed description
  • the embodiment of the present application solves the technical problem that a storage redundancy dictionary needs a large storage space in the prior art by providing a method, a device, and an electronic device for estimating an angle of arrival.
  • the present invention provides an apparatus for estimating an angle of arrival, which is applied to an Orthogonal Frequency Division Multiplexing (OFDM) communication system, the apparatus comprising: an angle of arrival estimation unit for obtaining antenna array measurement results of M array elements and based on an antenna array Measurement results and redundant dictionary, estimating the angle of arrival of K wireless signals,
  • the redundant dictionary used here is a partial discrete Fourier transform matrix, that is to say, the redundant dictionary has a discrete Fourier transform factor, and then the multiplexed frequency division multiplexing communication system can be multiplexed.
  • a part of the complete discrete Fourier transform matrix does not need to consume additional storage space to store the redundant dictionary. Therefore, the technical problem of requiring a large storage space for the storage redundancy dictionary existing in the prior art is effectively solved, and the storage is saved. space.
  • the present invention provides an electronic device, which may be any node in an Orthogonal Frequency Division Multiplexing (OFDM) communication system, such as a terminal, a base station, and of course, other wireless communication devices.
  • OFDM Orthogonal Frequency Division Multiplexing
  • FIG. 1 is a functional block diagram of an electronic device for estimating an angle of arrival according to an embodiment of the present invention.
  • the electronic device includes: a processor 10, configured to obtain an antenna array measurement result of the M array elements; and an angle of arrival of the K wireless signals based on the antenna array measurement result and the redundancy dictionary, where the redundancy dictionary is partially discrete A Fourier transform matrix, where M and K are both positive integers.
  • the electronic device may further include a transceiver 20 and a memory 30. If the electronic device is a terminal, the display unit, the WIFI module, the I/O interface, and the like may also be included.
  • the above antenna array measurement result is obtained by measuring the wireless signal received by the antenna array, and then there are the following two cases.
  • the antenna array measurement result is obtained by the electronic device itself.
  • the electronic device further includes: an antenna array 11a, configured to perform a received signal measurement according to the measurement configuration to obtain an antenna array measurement result.
  • the antenna array 11a and the processor 10 are disposed in the same device. Then, after receiving the wireless signal transmitted by the signal source, the antenna array 11a measures the signal according to the measurement configuration, and generates an antenna array measurement result, and then The result is sent to the processor 10.
  • the antenna array 11a may only have the function of receiving signals, and when the antenna array 11a receives the wireless signal, it sends it The measurement processor is sent to the electronic device to generate an antenna array measurement result, and the result is sent to the processor 10.
  • the antenna array 11a may be a one-dimensional antenna array, such as a linear array of uniform hooks, or a multi-dimensional antenna array, such as a two-dimensional L-shaped array, which is not specifically limited in the present application.
  • the above measurement configuration includes at least a combination of one or more of sparsity, number of snapshots, and number of snapshot bits per snapshot.
  • the above-mentioned measurement configuration can be determined by the electronic device without being sent by other devices.
  • the antenna array measurement result is measured and sent by the peer device.
  • the electronic device further includes: a receiver lib, configured to receive the antenna array measurement result, and measure the antenna array. The result is sent to the processor.
  • the antenna array and the processor 10 are disposed in different devices, that is, the antenna device is disposed on the opposite device, and after the antenna array measures the wireless signal received by the antenna array and obtains the antenna array measurement result, The measurement result is sent to the electronic device and received by the receiver lib, which in turn is sent to the processor 10.
  • the measurement configuration is sent by the base station to the terminal or negotiated with the base station, and after the terminal completes the measurement and generates the antenna array measurement score, the terminal It is sent to the base station, received by the receiver lib on the base station and sent to the processor 10.
  • the base station may transmit the measurement configuration determined by the terminal or negotiate with the terminal, and send the antenna array measurement result to the receiver lib on the terminal, and It is forwarded to the processor 10 by the receiver lib.
  • the antenna array is a uniform linear array as an example.
  • FIG 3 is a schematic diagram of the structure of a uniform linear array.
  • the antenna array has M array elements, the distance between the two elements is d, and the angle of the wireless signal to each element is (9, where ⁇ is the angle between the direction of arrival of the signal and the direction of the antenna array.
  • the redundant dictionary is a partial discrete Fourier transform matrix, where the discrete Fourier transform factor is a Fourier transform factor in OFDM in the multiplexed electronic device, so that the electronic device can construct redundancy offline. Dictionary, without the need to store the Fourier factor in the redundant dictionary, It is only necessary to call the stored Fourier factor of OFDM at the time of use, saving system resources.
  • the range of trigonometric values is evenly divided into N grids.
  • the redundant dictionary described above can be represented by the sine function of the angle.
  • the sine function value is at -90. ⁇ 90.
  • a person skilled in the art can set a redundancy dictionary by a cosine function or a redundant dictionary by a sine function, which is not specifically limited in the present application.
  • a redundant dictionary is constructed by a cosine function as an example.
  • the direction vector of each of the N grids is calculated, and the N direction vectors form a redundancy dictionary.
  • the direction vector is first obtained for the angle corresponding to each of the N grids, wherein the direction vector of the i-th grid is as shown in the formula (1).
  • N is an integer power of 2
  • FFT Fast Fourier Transformation
  • the processor 10 begins to estimate the angle of arrival.
  • the processor 10 can obtain a sparse signal matrix using a compressed sensing algorithm based on the antenna array measurement results and the redundancy dictionary.
  • a pessimistic algorithm-based compressed sensing algorithm is used to input the redundant dictionary A, the array measurement result Y, that is, the received signal matrix, and the sparsity K, that is, the number of angles to be estimated, and the recovered sparse signal matrix is calculated.
  • X where the set of locations of non-zero elements in X is T.
  • SOMP Simultaneous orthogonal Match Pursuit
  • Step 2 Iterate the calculation to get the X that satisfies the condition.
  • Estimate the value update that is, find the pseudo inverse matrix of A k and multiply Y by right to obtain the estimated value X.
  • the transposed matrix of A appearing in step a is multiplied by r, that is, the first matrix, and A is right multiplied by X in step f, that is, the second matrix, at N
  • r that is, the first matrix
  • X that is, the second matrix
  • the transpose matrix of A is multiplied by r
  • the FFT inverse transform can be performed on the second matrix X to obtain a third matrix, that is, an FFT inverse transform matrix of X. Then, an even term in the first M term in the third matrix is obtained. And determining an even term of the operation result of the redundancy dictionary A by the second matrix X, and determining an inverse value of the odd term in the first M term in the third matrix as an odd term of the operation result.
  • the angle of arrival can be estimated based on the sparse signal matrix.
  • the non-zero element row of the sparse signal matrix obtained by the first step is subscripted as i, and its angle of arrival is ar CC0S (l- 2i / N).
  • the positions of some non-zero elements in the recovered sparse signal matrix may not be on the grid, then, in order to obtain the positions of these non-grid elements, when performing the compressed sensing algorithm SBI-SVD, Sparsity Bayesian Inference-Singular Value Decomposition, Off- Grid Direction of Arrival Estimation Using Sparse Bayesian Inferenc, IEEE Transaction on Signal Processing, vol, 61, No l, Janary, 2013), or Simultaneous Orthogonal Match Pursuit-Last Squared, An alternating descent algorithm for the off-grid DOA estimation problem with sparsity constraints, EUSIPCO, 2012) to obtain a sparse signal matrix.
  • the first derivative matrix of the redundant dictionary is used to obtain the sparse signal matrix.
  • elements in the first-order inverse matrix of the redundancy dictionary may be as shown in equation (5).
  • jj is the number of rows of the first derivative matrix
  • ⁇ is the number of columns of the first derivative matrix
  • the angle of arrival is arccos (l-2). (i+diag(A)) /N).
  • the M array elements may include array elements having at least two dimensions, that is, the antenna array has at least two dimensions.
  • the antenna array is a two-dimensional array, such as an L-shaped antenna array.
  • FIG. 4 is a schematic structural diagram of an L-shaped antenna array.
  • the processor 10 is configured to obtain antenna array measurement results of the array elements of each of the M array elements respectively; for each array element, the antenna array is used to measure the flaws and redundancy.
  • the remainder dictionary estimated K arrival angles; based on the K arrival angles of the array elements of each dimension, the K arrival angles of the M array elements are estimated.
  • the processor 10 obtains a first dimension, such as an antenna array measurement result of the array element in the horizontal direction, and a second dimension, such as an antenna measurement result of the vertical direction array element, where the antenna array of each dimension
  • a first dimension such as an antenna array measurement result of the array element in the horizontal direction
  • a second dimension such as an antenna measurement result of the vertical direction array element
  • the device with the antenna array divides the antenna array into two independent one-dimensional homogenous linear arrays, and estimates the spatial angles ⁇ and ⁇ of each uniform linear array respectively, as shown in Fig. 5, and Fig. 5 is a spatial angle. Schematic diagram of the relationship with the direction angle.
  • sparseness indicates that the antenna receives the signal and the spatial angles ⁇ , ⁇ .
  • the vertical direction that is, the antenna on the y-axis
  • Y x is a received signal matrix representing the array elements on the X-axis
  • Y y is the received signal matrix of the array elements on the y-axis
  • E x , E y are the noise matrices of the signals on the X-axis and y-axis array elements, respectively.
  • the algorithm of compressed sensing to estimate cos a, cos ⁇ according to Y x , A, K, and Y y , A, K respectively. That is, the method for estimating the angle of arrival by using the one-dimensional antenna array for the antenna arrays in two dimensions, that is, using the methods in one or more of the above embodiments to obtain ⁇ ⁇ and X y respectively , then, according to X x Position of non-zero data in i x , get coso
  • the processor 10 is based on ⁇ ⁇ , ⁇ , ⁇ , ⁇ , and Y y , A, B , K , respectively.
  • ⁇ is the first derivative matrix of ⁇
  • COSy 9 l_2 (i y +diag(Ay))/N is obtained.
  • both i x and i y can be vectors.
  • an angle matching algorithm such as fitting to implement an inverse trigonometric function, is used to calculate the angle of arrival.
  • the electronic device in the one or more embodiments is applied to an Orthogonal Frequency Division Multiplexing (OFDM) communication system
  • the electronic device includes: a processor, configured to obtain an antenna array measurement result of the M array elements and Estimating the angle of arrival of the K wireless signals based on the antenna array measurement results and the redundancy dictionary
  • the redundant dictionary used herein is a partial discrete Fourier transform matrix, that is, The redundant dictionary has a discrete Fourier transform factor, and then a part of the complete discrete Fourier transform matrix in the orthogonal frequency division multiplexing communication system can be multiplexed without consuming additional storage space for storing the redundant dictionary. Therefore, the technical problem that the storage redundancy dictionary existing in the prior art requires a large storage space is effectively solved, and system resources are saved.
  • the present invention also provides an apparatus for estimating an angle of arrival.
  • the apparatus is applied to an OFDM communication system, and the apparatus can be disposed on any node of the communication system, as shown in FIG. 6, FIG. A functional block diagram of the apparatus for estimating the angle of arrival in this embodiment.
  • the device includes: an antenna array measurement result receiving unit 61, configured to obtain antenna array measurement results of M array elements; an arrival angle estimation unit 62, configured to estimate arrival of K wireless signals based on antenna array measurement results and a redundancy dictionary An angle, where the redundancy dictionary is a partial discrete Fourier transform matrix, and M and K are positive integers.
  • the antenna array measurement result receiving unit 61 is specifically configured to: obtain the antenna array measurement result by performing the received signal measurement according to the measurement configuration, where the foregoing measurement configuration includes at least sparseness One or more of the number of the metrics, the number of snapshots, and the number of the Snapshots may include other parameters, which are not specifically limited in this application.
  • the antenna array measurement result receiving unit 61 is specifically configured to: receive the antenna array measurement result sent by the peer device.
  • the antenna array may be a one-dimensional antenna array or a multi-dimensional antenna array.
  • the angle of arrival estimating unit 62 may include: an antenna array.
  • the measurement result input subunit is used to respectively obtain antenna array measurement results of the array elements of each of the M array elements; the single dimension angle of arrival calculation subunit is used for the array elements of each dimension, based on the antenna array measurement result And a redundant dictionary, estimating K arrival angles; a multi-dimensional arrival angle calculation sub-unit for calculating K arrival angles of M array elements based on K arrival angles of the array elements of each dimension.
  • the redundant dictionary is obtained by: uniformly dividing a range of trigonometric values into N meshes; calculating a direction vector of each of the N meshes, and composing a redundancy dictionary by the N direction vectors.
  • the direction vector is determined according to the formula (1), then the redundant dictionary is as public Formula (2) is shown.
  • the angle-of-arrival estimation unit 62 may include: a sparse signal matrix calculation sub-unit, configured to obtain a sparse signal matrix by using a compressed sensing algorithm based on the antenna array measurement result and the redundancy dictionary, wherein the compressed sensing algorithm includes at least: The Fourier algorithm calculates the operation of the transposed right-of-first matrix of the redundant dictionary and/or the operation of the redundant dictionary by the second matrix.
  • the first matrix and the second matrix are matrices related to the antenna array measurement results; Unit, used to estimate the angle of arrival based on the sparse signal matrix.
  • the non-zero element is not on the grid
  • the sparse signal matrix calculation sub-unit is specifically configured to: based on the antenna array measurement result, the redundancy dictionary, and the first derivative matrix of the redundancy dictionary, using a compressed sensing algorithm Obtain a sparse signal matrix.
  • the elements in the first derivative matrix are as shown in equation (5).
  • the sparse signal matrix calculation subunit specifically includes: an extended subunit for expanding the first matrix into an N-dimensional matrix; and an inverse fast Fourier transform sub-unit for performing an inverse fast Fourier transform on the N-dimensional matrix And perform a fast Fourier transform shift.
  • the sparse signal matrix calculation subunit specifically includes: an inverse fast Fourier transform subunit, configured to perform an inverse fast Fourier transform on the second matrix to obtain a third matrix; and determine a subunit for using the third matrix
  • the even term in the first M term in the middle is determined as the even term of the operation result of the redundant dictionary by the second matrix, and the inverse value of the odd term in the first M term in the third matrix is determined as the operation result Odd items.
  • the present invention provides a method for estimating an angle of arrival, which is applied to any node in an Orthogonal Frequency Division Multiplexing (OFDM) communication system, and may be applied to a terminal or a base station, which is not specifically limited in this application. .
  • OFDM Orthogonal Frequency Division Multiplexing
  • S101 Obtain an antenna array measurement score of M array elements
  • S102 Estimating the angle of arrival of the K wireless signals based on the antenna array measurement result and the redundancy dictionary, wherein the redundancy dictionary is a partial discrete Fourier transform matrix, and M and K are positive integers.
  • S101 can have two embodiments. First, the antenna array measurement result is obtained by the node through the antenna array set by the node. Then, the S101 may be: performing the received signal measurement according to the measurement configuration to obtain the antenna array measurement result, where the measurement configuration includes at least the sparsity degree, A combination of one or more of the number of snapshots and the number of snapshot bits per snapshot. Second, the antenna array measurement result is sent by the peer device provided with the antenna array. Then, the S101 may be: receiving the antenna array measurement result sent by the peer device.
  • S101 may be: obtaining antenna array measurement results of array elements of each of M array elements respectively; based on antenna array measurement results and redundancy
  • the remainder dictionary, estimating the angle of arrival of the K wireless signals specifically includes: estimating the K arrival angles based on the antenna array measurement results and the redundancy dictionary for each dimension element; K arrivals based on the array elements of each dimension Angle, calculate the K arrival angles of M elements.
  • the redundancy dictionary in one or more of the above embodiments may be obtained by: uniformly dividing a range of trigonometric values into N grids; calculating a direction vector of each of the N grids, and The N direction vectors form a redundant dictionary.
  • the direction vector is as shown in the formula (1)
  • the corresponding redundancy dictionary is as shown in the formula (2).
  • S102 may include: obtaining a sparse signal matrix by using a compressed sensing algorithm based on the antenna array measurement result and the redundancy dictionary, wherein the compressed sensing algorithm includes at least: using a fast Fourier algorithm to calculate a transposed right multiplication of the redundancy dictionary The operation of a matrix and/or the operation of the redundancy dictionary by the second matrix.
  • the first matrix and the second matrix are matrices related to the antenna array measurement results; based on the sparse signal matrix, the angle of arrival is estimated.
  • S102 may be: obtaining a sparse signal matrix by using a compressed sensing algorithm based on the antenna array measurement result, the redundancy dictionary, and the first derivative matrix of the redundancy dictionary. .
  • the elements of the first derivative matrix are as shown in equation (6).
  • the first matrix is first expanded into an N-dimensional matrix, and then the N-dimensional matrix is further advanced.
  • the fast Fourier transform is performed and fast Fourier shift is performed.
  • the second matrix is inversely transformed by the second dictionary
  • the third matrix is obtained by first performing the inverse fast Fourier transform on the second matrix. Determining an even term of the first M term in the third matrix as an even term of the operation result of the redundancy dictionary right by the second matrix, and inverting the inverse of the odd term in the first M term in the third matrix, Determine the odd number of entries as the result of the operation.
  • the electronic device in the one or more embodiments is applied to an Orthogonal Frequency Division Multiplexing (OFDM) communication system
  • the electronic device includes: a processor, configured to obtain an antenna array measurement result of the array element and estimate an angle of arrival of the plurality of wireless signals based on the antenna array measurement result and the redundancy dictionary, where the redundant dictionary used herein is a partial discrete Fourier
  • the leaf transformation matrix that is to say, the redundant dictionary has a discrete Fourier transform factor, then a part of the complete discrete Fourier transform matrix in the orthogonal frequency division multiplexing communication system can be multiplexed without consuming additional
  • the storage space is used to store the redundancy dictionary. Therefore, the technical problem that the storage redundancy dictionary existing in the prior art requires a large storage space is effectively solved, and system resources are saved.
  • embodiments of the present invention can be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment, or a combination of software and hardware. Moreover, the invention can be embodied in the form of one or more computer program products embodied on a computer usable storage medium (including but not limited to disk storage, CD-ROM, optical storage, etc.) in which computer usable program code is embodied.
  • a computer usable storage medium including but not limited to disk storage, CD-ROM, optical storage, etc.
  • the present invention has been described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (system), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flowchart illustrations and/or FIG.
  • These computer program instructions can be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing device to produce a machine for the execution of instructions for execution by a processor of a computer or other programmable data processing device.
  • the computer program instructions can also be stored in a computer readable memory that can direct a computer or other programmable data processing device to operate in a particular manner, such that the instructions stored in the computer readable memory produce an article of manufacture comprising the instruction device.
  • the apparatus implements the functions specified in one or more blocks of a flow or a flow and/or block diagram of the flowchart.
  • These computer program instructions can also be loaded onto a computer or other programmable data processing device such that a series of operational steps are performed on a computer or other programmable device to produce computer-implemented processing for execution on a computer or other programmable device.
  • the instructions provide steps for implementing the functions specified in one or more of the flow or in a block or blocks of a flow diagram.

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Variable-Direction Aerials And Aerial Arrays (AREA)
  • Radio Transmission System (AREA)

Abstract

La présente invention concerne un procédé et un appareil d'estimation d'un angle d'arrivée, et un dispositif électronique. L'invention a pour objectif de résoudre le problème technique lié, dans l'état de la technique, au fait qu'un espace de stockage plus important est requis pour stocker un dictionnaire redondant. Le dispositif électronique, utilisé dans un système de communications à multiplexage par répartition orthogonale de la fréquence, comprend un processeur, utilisé pour obtenir des résultats de mesurage de réseau d'antennes de M éléments de réseau, et estimer des angles d'arrivée de K signaux radio d'après les résultats de mesurage de réseau d'antennes et un dictionnaire redondant, le dictionnaire redondant étant une matrice de transformée de Fourier partiellement discrète, M et K étant des nombres entiers positifs.
PCT/CN2013/090581 2013-12-26 2013-12-26 Procédé et appareil d'estimation d'un angle d'arrivée, et dispositif électronique WO2015096099A1 (fr)

Priority Applications (2)

Application Number Priority Date Filing Date Title
PCT/CN2013/090581 WO2015096099A1 (fr) 2013-12-26 2013-12-26 Procédé et appareil d'estimation d'un angle d'arrivée, et dispositif électronique
CN201380002828.4A CN104937856B (zh) 2013-12-26 2013-12-26 一种估计到达角的方法、装置及电子设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/CN2013/090581 WO2015096099A1 (fr) 2013-12-26 2013-12-26 Procédé et appareil d'estimation d'un angle d'arrivée, et dispositif électronique

Publications (1)

Publication Number Publication Date
WO2015096099A1 true WO2015096099A1 (fr) 2015-07-02

Family

ID=53477371

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2013/090581 WO2015096099A1 (fr) 2013-12-26 2013-12-26 Procédé et appareil d'estimation d'un angle d'arrivée, et dispositif électronique

Country Status (2)

Country Link
CN (1) CN104937856B (fr)
WO (1) WO2015096099A1 (fr)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018076072A1 (fr) 2016-10-28 2018-05-03 Macquarie University Estimation de direction d'arrivée
CN108957390A (zh) * 2018-07-09 2018-12-07 东南大学 一种存在互耦时基于稀疏贝叶斯理论的到达角估计方法
CN109039403A (zh) * 2018-09-14 2018-12-18 中国人民解放军空军预警学院 大规模mimo系统中基于冗余字典的下行链路信道估计方法
US20210033693A1 (en) * 2018-03-19 2021-02-04 King Abdullah University Of Science And Technology Ultrasound based air-writing system and method

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112740058B (zh) * 2019-12-26 2022-04-22 华为技术有限公司 一种到达角的测量方法和装置
CN111431575B (zh) * 2020-03-04 2022-02-11 哈尔滨工程大学 基于常规波束形成的来波方向稀疏重构方法
CN111551902B (zh) * 2020-06-02 2022-11-01 电子科技大学 基于压缩感知技术的fmcw雷达天线缺损时采集信号恢复方法
CN113253194B (zh) * 2021-04-21 2022-07-08 中国电子科技集团公司第二十九研究所 一种基于稀疏表示的宽频段到达角度和极化联合测量方法
CN116450993B (zh) * 2023-04-24 2023-12-12 哈尔滨工业大学 多测量向量卫星数据的处理方法、电子设备及存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070164902A1 (en) * 2005-12-02 2007-07-19 Samsung Electronics Co., Ltd. Smart antenna beamforming device in communication system and method thereof
CN101232316A (zh) * 2008-01-17 2008-07-30 华南理工大学 一种基于移动wimax的波达方向估计方法
CN101860387A (zh) * 2009-04-09 2010-10-13 中兴通讯股份有限公司 一种波达角估计的方法和系统
CN102608565A (zh) * 2012-03-23 2012-07-25 哈尔滨工程大学 一种基于均匀圆阵列的波达方向估计方法
CN103454615A (zh) * 2013-03-01 2013-12-18 西安电子科技大学 稳健盲信号参数估计方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4309110B2 (ja) * 2002-09-27 2009-08-05 パナソニック株式会社 適応アンテナ無線通信装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070164902A1 (en) * 2005-12-02 2007-07-19 Samsung Electronics Co., Ltd. Smart antenna beamforming device in communication system and method thereof
CN101232316A (zh) * 2008-01-17 2008-07-30 华南理工大学 一种基于移动wimax的波达方向估计方法
CN101860387A (zh) * 2009-04-09 2010-10-13 中兴通讯股份有限公司 一种波达角估计的方法和系统
CN102608565A (zh) * 2012-03-23 2012-07-25 哈尔滨工程大学 一种基于均匀圆阵列的波达方向估计方法
CN103454615A (zh) * 2013-03-01 2013-12-18 西安电子科技大学 稳健盲信号参数估计方法

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018076072A1 (fr) 2016-10-28 2018-05-03 Macquarie University Estimation de direction d'arrivée
EP3532858A4 (fr) * 2016-10-28 2020-10-07 MacQuarie University Estimation de direction d'arrivée
US11194005B2 (en) 2016-10-28 2021-12-07 Macquarie University Direction of arrival estimation
IL266212B1 (en) * 2016-10-28 2023-06-01 Univ Macquarie Direction for arrival assessment
US11782118B2 (en) 2016-10-28 2023-10-10 Macquarie University Direction of arrival estimation
US20210033693A1 (en) * 2018-03-19 2021-02-04 King Abdullah University Of Science And Technology Ultrasound based air-writing system and method
US11630518B2 (en) * 2018-03-19 2023-04-18 King Abdullah University Of Science And Technology Ultrasound based air-writing system and method
CN108957390A (zh) * 2018-07-09 2018-12-07 东南大学 一种存在互耦时基于稀疏贝叶斯理论的到达角估计方法
CN108957390B (zh) * 2018-07-09 2022-03-18 东南大学 一种存在互耦时基于稀疏贝叶斯理论的到达角估计方法
CN109039403A (zh) * 2018-09-14 2018-12-18 中国人民解放军空军预警学院 大规模mimo系统中基于冗余字典的下行链路信道估计方法

Also Published As

Publication number Publication date
CN104937856B (zh) 2018-06-05
CN104937856A (zh) 2015-09-23

Similar Documents

Publication Publication Date Title
WO2015096099A1 (fr) Procédé et appareil d'estimation d'un angle d'arrivée, et dispositif électronique
Yan et al. Fast DOA estimation based on a split subspace decomposition on the array covariance matrix
Wen et al. Direction finding in MIMO radar with large antenna arrays and nonorthogonal waveforms
Xue et al. IAA spectral estimation: fast implementation using the Gohberg–Semencul factorization
CN107589399B (zh) 基于多采样虚拟信号奇异值分解的互质阵列波达方向估计方法
CN106569172B (zh) 二维doa估计方法
CN107870315B (zh) 一种利用迭代相位补偿技术估计任意阵列波达方向方法
Liao et al. Direction-of-arrival estimation in subarrays-based linear sparse arrays with gain/phase uncertainties
CN107544051A (zh) 嵌套阵列基于k‑r子空间的波达方向估计方法
WO2016023308A1 (fr) Procédé et appareil de localisation pour une source de décharge partielle d'un transformateur
Xia Joint diagonalization based 2D-DOD and 2D-DOA estimation for bistatic MIMO radar
Yang et al. A variational Bayesian strategy for solving the DOA estimation problem in sparse array
Mao et al. Low complexity DOA estimation using AMP with unitary transformation and iterative refinement
El Badawy et al. Localizing unsynchronized sensors with unknown sources
Pan et al. DOA estimation with planar array via spatial finite rate of innovation reconstruction
Yan et al. Computationally efficient direction of arrival estimation with unknown number of signals
Cai et al. An improved expanding and shift scheme for the construction of fourth-order difference co-arrays
He et al. Efficient underwater two-dimensional coherent source localization with linear vector-hydrophone array
CN114375038B (zh) 定位方法、装置、基站、存储介质和计算机程序产品
Kristensson et al. Modified IQML and weighted subspace fitting without eigendecomposition
Tang et al. PSWF-based decoupled atomic norm minimization for DOD and DOA estimation in MIMO radar with arbitrary linear arrays
CN110138693B (zh) 结构化压缩感知信道估计方法及装置、存储介质和终端
Chen et al. A parallel linear solver for multilevel Toeplitz systems with possibly several right-hand sides
Xu et al. Fast Two‐Dimensional Direction‐of‐Arrival Estimation of Multiple Signals in Coprime Planar Array
Lv et al. Joint DOA and frequency estimation based on spatio-temporal co-prime sampling

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 13900379

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 13900379

Country of ref document: EP

Kind code of ref document: A1