CN113259278A - Millimeter wave multi-antenna channel estimation method and device - Google Patents

Millimeter wave multi-antenna channel estimation method and device Download PDF

Info

Publication number
CN113259278A
CN113259278A CN202110550050.7A CN202110550050A CN113259278A CN 113259278 A CN113259278 A CN 113259278A CN 202110550050 A CN202110550050 A CN 202110550050A CN 113259278 A CN113259278 A CN 113259278A
Authority
CN
China
Prior art keywords
channel estimation
sub
matrix
result
iteration
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.)
Granted
Application number
CN202110550050.7A
Other languages
Chinese (zh)
Other versions
CN113259278B (en
Inventor
漆渊
张纪焱
邓皓戈
左辰宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing University of Posts and Telecommunications
Original Assignee
Beijing University of Posts and Telecommunications
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 Beijing University of Posts and Telecommunications filed Critical Beijing University of Posts and Telecommunications
Priority to CN202110550050.7A priority Critical patent/CN113259278B/en
Publication of CN113259278A publication Critical patent/CN113259278A/en
Application granted granted Critical
Publication of CN113259278B publication Critical patent/CN113259278B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202Channel estimation
    • H04L25/0204Channel estimation of multiple channels
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202Channel estimation
    • H04L25/024Channel estimation channel estimation algorithms
    • H04L25/0242Channel estimation channel estimation algorithms using matrix methods

Abstract

The application discloses a millimeter wave multi-antenna channel estimation method, which comprises the following steps: splitting a matrix for channel estimation into m operation sub-matrices, wherein the operation sub-matrices are: and the ith sub-vector of the measurement vector of the observation signal is equal to the product of the ith first sub-matrix split by the first matrix and the channel estimation result to be estimated, i is a natural number which is more than or equal to 1 and less than or equal to m, the channel estimation sub-results of the m operation sub-matrices are respectively calculated in parallel to obtain m channel estimation sub-results, and the m channel estimation sub-results are averaged to obtain the channel estimation result to be estimated. The invention realizes distributed processing based on computing resources and is beneficial to reducing the running time of channel estimation. Through iterative convergence of each operation submatrix, the error between the channel estimation value and the true value is reduced, the reliability of channel estimation is improved, and the running time of the channel estimation is further reduced.

Description

Millimeter wave multi-antenna channel estimation method and device
Technical Field
The invention relates to the field of wireless communication, in particular to a millimeter wave multi-antenna channel estimation method.
Background
The channel estimation is an essential link for a wireless communication link, is particularly important for a large-scale multi-antenna system, and is a precondition for fully utilizing the multi-antenna gain of the system. In a large-scale multi-antenna system, operations such as precoding, signal detection, resource allocation and the like all require accurate Channel State Information (CSI).
At present, millimeter wave large-scale multi-antenna has definite application prospect in future mobile communication system, and channel estimation is a key link in millimeter wave large-scale multi-antenna scene. However, due to the increase of the antenna size, how to obtain accurate CSI in real time is a great challenge.
In millimeter wave multi-antenna channels, channel estimation simultaneously faces the challenge of channel sparsity, and the channels need to be described according to the sparsity characteristics of millimeter wave channels, limited arrival angles, departure angles and path gains.
Under non-sparse conditions, the channel estimation algorithm mainly includes a Least Square (LS) method and a Minimum Mean Square Error (MMSE) method. Due to the very complicated channel characteristics, the realization of these methods often needs to use long pilot sequences and large-scale channel coefficients, which is not practical for practical millimeter-wave large-scale multiple-input multiple-output (MIMO) systems.
In the millimeter wave massive MIMO system, since signal transmission is distributed on a part of paths, channels are sparse in time and angle, so that most of sparse channel estimation algorithms at present are based on a Compressive Sensing (CS) principle. In millimeter wave channel estimation based on compressed sensing, a multistage process can be used to avoid exhaustive search, but due to limited transmission power, the use of a wide beam is limited, and the method is difficult to apply to an actual channel. Common millimeter wave multi-antenna channel estimation is roughly of two types: one is a non-Bayesian based algorithm, such as Orthogonal Matching Pursuit (OMP), which is an iterative algorithm that finds a sub-optimal solution, and the other is a Bayesian based algorithm that applies an estimation technique to identify a desired Sparse solution, such as Sparse Bayesian Learning (SBL) and Bayesian Compressive Sensing (BCS), using Bayesian methods and appropriate statistical assumptions.
The existing channel estimation algorithm is generally a serial algorithm, because a millimeter wave large-scale antenna is provided with a large-scale antenna array, a measurement matrix and a complete dictionary matrix related to channel estimation have quite large scales, and simultaneously, because the complexity of matrix operation increases exponentially along with the increase of the matrix scale, a large amount of time is consumed by utilizing the serial algorithm, and the time complexity is high.
Disclosure of Invention
The invention provides a millimeter wave multi-antenna channel estimation method, which is used for reducing the running time of channel estimation.
The millimeter wave multi-antenna channel estimation method provided by the invention is realized as follows:
splitting a first matrix for channel estimation into m operation submatrices, wherein the operation submatrices are: the ith sub-vector of the measurement vector of the observation signal is equal to the product of the ith first sub-matrix split by the first matrix and the channel estimation result to be estimated, i is a natural number which is more than or equal to 1 and less than or equal to m,
respectively calculating the channel estimation sub-results of the m operation sub-matrixes in parallel to obtain m channel estimation sub-results,
and averaging the m channel estimation sub-results to obtain a channel estimation result to be estimated.
The present invention also provides a millimeter wave multi-antenna channel estimation apparatus, comprising,
a decomposition unit, configured to split the first matrix for channel estimation into m operator matrices, where the operator matrices are: the ith sub-vector of the measurement vector of the observation signal is equal to the product of the ith first sub-matrix split by the first matrix and the channel estimation result to be estimated, i is a natural number which is more than or equal to 1 and less than or equal to m,
m first calculation units, each for calculating the channel estimation sub-results of the m operation sub-matrices in parallel to obtain m channel estimation sub-results,
and the second calculation unit is used for averaging the m channel estimation sub-results to obtain a channel estimation result to be estimated.
According to the millimeter wave multi-antenna channel estimation method provided by the invention, a large-scale matrix for channel estimation is divided into a plurality of small-scale operation sub-matrixes, so that distributed processing can be carried out according to computing resources, channel estimation sub-results are respectively obtained through each operation sub-matrix in parallel, and the running time of channel estimation is favorably reduced. Furthermore, through iterative convergence of each operation submatrix, the error between the channel estimation value and the true value is reduced, the reliability of channel estimation is improved, and the running time of the channel estimation is further reduced. The method and the device can be suitable for various application scenes by adjusting the number of the operation submatrices, and are high in flexibility.
Drawings
Fig. 1 is a schematic flow chart of a millimeter wave multi-antenna channel estimation method.
Fig. 2 is a schematic flow chart of solving channel estimation in a parallel manner.
Fig. 3 is a diagram illustrating convergence characteristics in solving for channel estimation.
Fig. 4 is a diagram illustrating the convergence of channel estimation with different noise powers.
Fig. 5 is a schematic diagram of a probability density curve of error distribution at different iteration numbers.
Fig. 6 is a schematic diagram of the millimeter wave multi-antenna channel estimation apparatus according to the present application.
Detailed Description
For the purpose of making the objects, technical means and advantages of the present application more apparent, the present application will be described in further detail with reference to the accompanying drawings.
The millimeter wave multi-antenna channel estimation method divides a large-scale matrix for channel estimation into a plurality of small-scale operation sub-matrixes, completes operation on each operation sub-matrix in parallel through distributed computation, and carries out iterative processing on channel estimation sub-results of each operation sub-matrix, gradually converges, and reduces errors until the error is similar to a correct solution.
Referring to fig. 1, fig. 1 is a schematic flow chart of a millimeter wave multi-antenna channel estimation method. The method comprises the steps of (1) carrying out,
step 101, splitting a first matrix for channel estimation into m operator matrixes, wherein the operator matrixes are: the ith sub-vector of the measurement vector of the observation signal is equal to the product of the ith first sub-matrix split by the first matrix and the channel estimation result to be estimated, i is a natural number which is more than or equal to 1 and less than or equal to m,
step 102, calculating the channel estimation sub-results of the m operation sub-matrices in parallel to obtain m channel estimation sub-results,
and 103, averaging the m channel estimation sub-results to obtain a channel estimation result to be estimated.
The channel estimation method provided by the application is based on distributed computation, the reliability of a computation result is guaranteed, the channel estimation speed can be effectively improved, and the number and scale of operation sub-matrixes can be adjusted according to computation resources, so that the method is suitable for different application scenarios and high in flexibility.
The following detailed description is provided to facilitate understanding of the invention.
The signal matrix X is transmitted from a transmitting end, passes through a precoding matrix F, is transmitted through a channel matrix H, reaches a receiving end, first passes through a combining matrix W, and is received by the receiving end to obtain an observation matrix Y. During the channel transmission, the signal will be superimposed with white gaussian noise Q. Thus, the observation matrix Y can be expressed as:
Figure BDA0003075074550000031
wherein the content of the first and second substances,
Figure BDA0003075074550000032
c is a set of plural numbers, NrFor the number of receiving antennas, NtFor the number of transmitting antennas, MrTo receive a measurement sequence, MtFor transmitting measurement sequencesColumn, signal matrix X is of scale MtWith power P. The noise matrix Q is a complex noise matrix with a power of PnThe real and imaginary matrix scales are both Mr×MtThe noise matrix is randomly generated.
After vectorizing the signal, obtaining an observation signal expression after vectorization:
Figure BDA0003075074550000033
wherein A isDFor a dictionary matrix, each column represents the Kronecker product:
Figure BDA0003075074550000041
Figure BDA0003075074550000042
to know
Figure BDA0003075074550000043
Respectively the kth and jth angles, h, corresponding to the dictionary matrixbIs the channel gain vector under the dictionary matrix; and nQ is a noise vector after the complex noise matrix is vectorized.
Defining a first matrix
Figure BDA0003075074550000044
The above formula can be further rewritten as:
yv=Ahb+nQ
thus, the channel estimation problem can be translated into a channel estimation by the first matrix A and the measurement vector yvTo estimate hbThe problem is solved by the matrix of (1). Since the channel estimation is affected by noise, the above equation is rewritten as:
Figure BDA0003075074550000045
in order to conveniently adopt distributed parallel solving, large-scale matrix operation is carried out
Figure BDA0003075074550000046
Splitting into m sub-matrix operations, wherein the number of m can be determined according to computing resources, and any operation sub-matrix i is
Figure BDA0003075074550000047
Each operator matrix can be calculated by a respective calculation unit, and i is a natural number equal to or greater than 1 and equal to or less than m.
It should be understood that the computing unit may be for physical devices, or may be a computing task distributed in each physical device or in the same physical device.
Thus, the large-scale matrix operation is decomposed into parallel operations of m operator matrices, wherein the operator matrices are: i-th subvector y of the measurement vector of the observed signaliEqual to the ith first sub-matrix A divided by the first matrixiAnd the product of the channel estimation result to be estimated. Expressed mathematically as:
Figure BDA0003075074550000048
referring to fig. 2, fig. 2 is a schematic flow chart of solving channel estimation in a parallel manner. The method comprises the steps of (1) carrying out,
step 201, any computing unit solves the channel estimation sub-result of the first iteration of the operator sub-matrix i through the operator sub-matrix i, and records the channel estimation sub-result as
Figure BDA0003075074550000049
Thus, the channel estimation sub-results of m operation sub-matrixes are obtained in total,
step 202, the m channel estimation sub-results are summed and averaged to obtain the channel estimation result of the first iteration
Figure BDA00030750745500000410
Expressed mathematically as:
Figure BDA0003075074550000051
step 203, for any operation sub-matrix i, calculating the channel estimation sub-result of the current iteration of the operation sub-matrix in a manner that the channel estimation sub-result calculated by the current iteration of the operation sub-matrix is converged by using the channel estimation sub-result calculated by the operation sub-matrix and the channel estimation result of the previous iteration,
one of the embodiments of the method for converging the channel estimation result calculated by the current iteration of the operation submatrix is as follows:
updating the difference between the channel estimation sub-result calculated last time of the operation sub-matrix and the channel estimation result of the last iteration by using iteration update parameters to obtain an updated value; accumulating the update value on the basis of the channel estimation sub-result calculated last time by the operation sub-matrix; wherein, the update parameter is the product of the conjugate of the first parameter and the conjugate of the second parameter, the first parameter is related to the maximum eigenvalue and the minimum eigenvalue of the set second matrix, and the second parameter is related to the first sub-matrix.
Thus, for the channel estimation sub-result of the i-th iteration calculated by the operation sub-matrix i, the following formula can be expressed as follows:
Figure BDA0003075074550000052
wherein, the first parameter used for iterative update is gamma, and the second parameter is the first sub-matrix A in the operation sub-matrixiZero space projection of (1), noted as betai
Figure BDA0003075074550000053
I is an identity matrix and is a matrix of the identity,
Figure BDA0003075074550000054
is an operator matrix iThe calculated channel estimation sub-result of the l-th iteration,
Figure BDA0003075074550000055
as a result of the channel estimation of the l-th iteration,
l is more than or equal to 1 and less than or equal to L which is a set iteration threshold.
And step 204, carrying out weighted summation on the channel estimation sub-results of the current iteration calculated by the m operation sub-matrixes i to obtain the channel estimation result of the current iteration.
One of the weighted summation modes is as follows: taking the third parameter for weighted iteration updating as a weighting coefficient, weighting the average value of the channel estimation sub-results of the current iteration calculated by the m operation sub-matrixes i, accumulating the weighting result and the channel estimation result of the previous iteration,
preferably, the channel estimation result of the previous iteration is weighted by using the remaining weights, and the weighted channel estimation result of the previous iteration and the weighted result are accumulated.
The third parameter is related to a maximum eigenvalue and a minimum eigenvalue of the second matrix.
In this way, the channel estimation sub-results of the l iteration calculated by the m operation sub-matrices i are subjected to weighted summation to obtain the channel estimation result of the l iteration, which can be expressed by a mathematical expression:
Figure BDA0003075074550000061
wherein the third parameter is η.
The first parameter and the third parameter are obtained by the following method:
setting a second matrix expressed by the following mathematical formula:
Figure BDA0003075074550000062
obtaining the maximum eigenvalue mu of the matrix according to the second matrixmaxAnd minimum eigenvalue μminAnd constructing an equation set:
Figure BDA0003075074550000063
and solving the equation set to obtain a first parameter and a third parameter. Wherein the third parameter is a real number and the first parameter is within 0 to 2.
And step 205, accumulating the current iteration times, returning to step 203 until the accumulated iteration times reach a set iteration threshold, and outputting a final channel estimation result.
Referring to fig. 3 to 5, fig. 3 to 5 are schematic diagrams illustrating the performance of the preliminary demonstration method through numerical simulation. The number of F matrixes and W matrixes involved in the simulation is less than 4; number of transmitted measurement sequences MtThe number of received measurement sequences Mr is 8 respectively; number of transmitting antennas NtAnd number of receiving antennas NrIs 8; when the calculation units are divided, 8 calculation units are adopted for parallel calculation; the transmission angle is selected to be 0 degrees and 90 degrees, and the receiving angle is selected to be 45 degrees and 135 degrees.
The complex noise matrix Q and the first matrix A are random matrices, so that in order to avoid the accidental of experiments, simulation results are drawn by adopting a method of calculating and averaging for many times, and the accidental errors of the experiments are reduced as much as possible. The attached simulation graphs are drawn by adopting a method of calculating 500 times of averaging.
According to the simulation diagram, the problem of large-scale multi-antenna channel estimation can be effectively solved through lower iteration times. By increasing the number of iterations, the estimated value can be gradually converged toward the true value. When the number of iterations increases to a certain extent (e.g., l > 15), the channel estimate has substantially converged, after which the number of iterations no longer plays a decisive role in reducing the error between the channel estimate and the true value, which error will mainly depend on the magnitude of the noise power.
The invention can be flexibly changed aiming at antenna arrays and channels with different scales. It is mainly characterized in that: 1. and (3) dividing the operation submatrix 2 and selecting parameters. When the size of the transmitting signal matrix and the measuring signal matrix is large, various division options are available, for example, a matrix of 16 × 16 is selected, and the matrix can be divided into 8 calculation units for operation, or can be selected to be divided into 4 machines for operation. For different application scenarios, the requirements for algorithm precision and convergence speed are different, so that different division modes can be selected for channel estimation. The same reasoning is also adopted when the parameters are selected, the selected parameters are different in the matrixes with different scales, and the convergence speed of the matrixes can be changed by changing the parameters.
Referring to fig. 6, fig. 6 is a schematic diagram of the millimeter wave multi-antenna channel estimation apparatus according to the present invention. The device comprises a plurality of devices which are connected with each other,
a decomposition unit, configured to split the first matrix for channel estimation into m operator matrices, where the operator matrices are: the ith sub-vector of the measurement vector of the observation signal is equal to the product of the ith first sub-matrix split by the first matrix and the channel estimation result to be estimated, i is a natural number which is more than or equal to 1 and less than or equal to m,
m first calculation units, each for calculating the channel estimation sub-results of the m operation sub-matrices in parallel to obtain m channel estimation sub-results,
and the second calculation unit is used for averaging the m channel estimation sub-results to obtain a channel estimation result to be estimated.
The first computing unit may include a first computing unit,
the first calculation module is used for calculating the channel estimation sub-result of the current iteration of the operation sub-matrix according to the mode of converging the channel estimation sub-result of the current iteration of the operation sub-matrix by using the channel estimation sub-result of the last iteration of the operation sub-matrix and the channel estimation result of the last iteration;
and the iteration number control module is used for controlling the iteration number of each first calculation module.
The second calculation unit comprises a second calculation unit,
and the second calculation module is used for carrying out weighted summation on the channel estimation sub-results of the current iteration of the m operation sub-matrixes to obtain the channel estimation result of the current iteration.
In this document, relational terms such as first and second, and the like may be used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Also, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other identical elements in a process, method, article, or apparatus that comprises the element.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.

Claims (10)

1. A millimeter wave multi-antenna channel estimation method is characterized by comprising the following steps:
splitting a first matrix for channel estimation into m operation submatrices, wherein the operation submatrices are: the ith sub-vector of the measurement vector of the observation signal is equal to the product of the ith first sub-matrix split by the first matrix and the channel estimation result to be estimated, i is a natural number which is more than or equal to 1 and less than or equal to m,
respectively calculating the channel estimation sub-results of the m operation sub-matrixes in parallel to obtain m channel estimation sub-results,
and averaging the m channel estimation sub-results to obtain a channel estimation result to be estimated.
2. The method of claim 1, wherein said calculating the channel estimation sub-results of the m operation sub-matrices in parallel, respectively, to obtain m channel estimation sub-results, comprises:
for any one of the operator matrices i,
and calculating the channel estimation sub-result of the current iteration of the operation sub-matrix according to the mode of converging the channel estimation sub-result of the current iteration of the operation sub-matrix by using the channel estimation sub-result of the last iteration of the operation sub-matrix and the channel estimation result of the last iteration.
3. The method as claimed in claim 2, wherein said averaging the m channel estimation sub-results to obtain the channel estimation result to be estimated comprises:
weighting and summing the channel estimation sub-results of the current iteration of the m operation sub-matrixes to obtain the channel estimation result of the current iteration,
and repeatedly and respectively calculating the channel estimation sub-result of the next iteration of the m operation sub-matrixes and the channel estimation result of the next iteration until the set iteration time threshold is reached.
4. The method as claimed in claim 3, wherein said calculating the channel estimation sub-result of the current iteration of the operator matrix according to the way of converging the channel estimation sub-result of the current iteration of the operator matrix by using the channel estimation sub-result of the last iteration of the operator matrix and the channel estimation result of the last iteration of the operator matrix comprises:
updating the difference between the channel estimation sub-result of the current iteration of the operation sub-matrix i and the channel estimation result of the current iteration by using iteration update parameters to obtain an updated value;
accumulating the updated value on the basis of the channel estimation sub-result of the current iteration of the operation sub-matrix i;
the weighted summation of the channel estimation sub-results of the current iteration calculated by the m operation sub-matrices i to obtain the channel estimation result of the current iteration includes:
taking the third parameter for updating the weighted iteration as a weighting coefficient, weighting the average value of the channel estimation sub-results of the current iteration of the m operation sub-matrixes to obtain a first weighting result,
weighting the channel estimation result of the last iteration by the residual weights to obtain a second weighting result,
and accumulating the first weighting result and the second weighting result to obtain a channel estimation result of the iteration.
5. The method as claimed in claim 4, wherein said updating the difference between the sub-result of channel estimation of the current iteration of said operator matrix and the result of channel estimation of the current iteration by using an iteration update parameter to obtain an updated value comprises:
multiplying the conjugate of the first parameter for iterative update and the conjugate of the second parameter corresponding to the operator sub-matrix by the difference between the channel estimation sub-result of the current iteration of the operator sub-matrix and the channel estimation result of the current iteration,
wherein the content of the first and second substances,
the first parameter is related to the maximum eigenvalue and the minimum eigenvalue of the set second matrix,
the second parameter is associated with a first sub-matrix of the operator sub-matrices.
6. The method of claim 5 wherein the second parameter is a null-space projection of a first one of the operator submatrices,
the first parameter and the third parameter are obtained through the following method:
constructing a second matrix by using a first sub-matrix in the operation sub-matrices, wherein the second matrix is as follows:
Figure FDA0003075074540000021
wherein A isiIs the first sub-matrix of the operator sub-matrix i,
Figure FDA0003075074540000022
is the first in the operator matrix iA transposed matrix of the sub-matrices, m being the number of the operation sub-matrices;
obtaining the maximum eigenvalue and the minimum eigenvalue of the matrix according to the second matrix,
and constructing the following equation set by using the maximum characteristic value and the minimum characteristic value:
Figure FDA0003075074540000023
wherein, mumaxIs the maximum eigenvalue, muminIs the minimum eigenvalue, η is the third parameter, γ is the first parameter,
and solving to obtain a first parameter and a third parameter according to the equation set.
7. The method of claim 6, wherein the zero-space projection of the first one of the operator submatrices is:
Figure FDA0003075074540000024
wherein, I is an identity matrix,
the first matrix is:
Figure FDA0003075074540000025
wherein, P is the power of the transmitted signal,
Figure FDA0003075074540000026
is a transposed matrix of the precoding matrix at the transmitting end,
Figure FDA0003075074540000027
combining the conjugates of the matrices for the transmitting end, ADIs a dictionary matrix.
8. A millimeter wave multi-antenna channel estimation apparatus, comprising,
a decomposition unit, configured to split the first matrix for channel estimation into m operator matrices, where the operator matrices are: the ith sub-vector of the measurement vector of the observation signal is equal to the product of the ith first sub-matrix split by the first matrix and the channel estimation result to be estimated, i is a natural number which is more than or equal to 1 and less than or equal to m,
m first calculation units, each for calculating the channel estimation sub-results of the m operation sub-matrices in parallel to obtain m channel estimation sub-results,
and the second calculation unit is used for averaging the m channel estimation sub-results to obtain a channel estimation result to be estimated.
9. The apparatus of claim 8, wherein the first computing unit comprises,
the first calculation module is used for calculating the channel estimation sub-result of the current iteration of the operation sub-matrix according to the mode of converging the channel estimation sub-result of the current iteration of the operation sub-matrix by using the channel estimation sub-result of the last iteration of the operation sub-matrix and the channel estimation result of the last iteration;
and the iteration number control module is used for controlling the iteration number of each first calculation module.
10. The apparatus of claim 8, wherein the second computing unit comprises,
and the second calculation module is used for carrying out weighted summation on the channel estimation sub-results of the current iteration of the m operation sub-matrixes to obtain the channel estimation result of the current iteration.
CN202110550050.7A 2021-05-20 2021-05-20 Millimeter wave multi-antenna channel estimation method and device Active CN113259278B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110550050.7A CN113259278B (en) 2021-05-20 2021-05-20 Millimeter wave multi-antenna channel estimation method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110550050.7A CN113259278B (en) 2021-05-20 2021-05-20 Millimeter wave multi-antenna channel estimation method and device

Publications (2)

Publication Number Publication Date
CN113259278A true CN113259278A (en) 2021-08-13
CN113259278B CN113259278B (en) 2022-07-19

Family

ID=77182913

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110550050.7A Active CN113259278B (en) 2021-05-20 2021-05-20 Millimeter wave multi-antenna channel estimation method and device

Country Status (1)

Country Link
CN (1) CN113259278B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115118317A (en) * 2022-05-27 2022-09-27 成都中科微信息技术研究院有限公司 Iterative precoding multi-stream method, medium and device suitable for millimeter waves

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1809043A (en) * 2005-01-17 2006-07-26 都科摩(北京)通信技术研究中心有限公司 Frequency domain sub-space channel estimation device and method, receiver and signal receiving method
US20150295733A1 (en) * 2014-04-09 2015-10-15 Altiostar Networks, Inc. Sparse Ordered Iterative Group Multi-Antenna Channel Estimation
CN109560841A (en) * 2018-12-13 2019-04-02 东北大学 Extensive mimo system channel estimation methods based on improved distributed compression perception algorithm

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1809043A (en) * 2005-01-17 2006-07-26 都科摩(北京)通信技术研究中心有限公司 Frequency domain sub-space channel estimation device and method, receiver and signal receiving method
US20150295733A1 (en) * 2014-04-09 2015-10-15 Altiostar Networks, Inc. Sparse Ordered Iterative Group Multi-Antenna Channel Estimation
CN109560841A (en) * 2018-12-13 2019-04-02 东北大学 Extensive mimo system channel estimation methods based on improved distributed compression perception algorithm

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
ALAM ZAIB等: "Distributed Channel Estimation and Pilot Contamination Analysis for Massive MIMO-OFDM Systems", 《IEEE TRANSACTIONS ON COMMUNICATIONS》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115118317A (en) * 2022-05-27 2022-09-27 成都中科微信息技术研究院有限公司 Iterative precoding multi-stream method, medium and device suitable for millimeter waves
CN115118317B (en) * 2022-05-27 2023-10-03 成都中科微信息技术研究院有限公司 Iterative precoding multi-stream method, medium and device suitable for millimeter waves

Also Published As

Publication number Publication date
CN113259278B (en) 2022-07-19

Similar Documents

Publication Publication Date Title
CN111181619B (en) Millimeter wave hybrid beam forming design method based on deep reinforcement learning
KR100931901B1 (en) Beamforming Method and System of Transmitter for Reduced Complexity Multiple Input Multiple Output Transceiver
CN110824415B (en) Sparse direction-of-arrival angle estimation method based on multiple-sending and multiple-receiving array
CN104977558A (en) Distributed source center direction-of-arrival estimation method based on Bayesian compressed perception
CN113411110A (en) Millimeter wave communication beam training method based on deep reinforcement learning
CN104698430A (en) Virtual array antenna based angle estimation method for improving accuracy
CN113259278B (en) Millimeter wave multi-antenna channel estimation method and device
Ajorloo et al. Antenna placement in a compressive sensing-based colocated MIMO radar
CN103378888A (en) Downlink beam forming method and downlink beam forming equipment
CN102158267B (en) Method and system for processing signals
CN114204971A (en) Iterative aggregation beamforming design and user equipment selection method
CN114239644A (en) Non-grid block sparse angle estimation method based on array element cross coupling optimization
CN106716855B (en) Method of setting cooperation parameter of communication system, computer program thereof, information storage device thereof, and communication system including a plurality of node devices
CN108896967B (en) Method and device for detecting distance extended target based on clutter covariance matrix estimation
CN113422627B (en) Hybrid precoding method and device based on adaptive gradient back propagation
CN115361258A (en) Large-scale MIMO sparse channel estimation method and related equipment
CN111965593B (en) Method and device for determining direction of arrival of signal
CN111555785B (en) Signal precoding method and device and electronic equipment
CN111800176B (en) Equivalent channel matrix generation method, precoding method, device, equipment and medium
KR100965100B1 (en) Adaptive beam former and method of the same
KR101599190B1 (en) MIMO Communication Method and System using the Block Circulant Jacket Matrices
CN114629533A (en) Information geometry method and system for large-scale MIMO channel estimation
KR20120089820A (en) Beam design method and algorithm for time-varying multiuser mimo interference channels transmitting multiple data streams
US11552827B2 (en) Signal estimation apparatus, signal estimation method and program recording medium
EP3726739A1 (en) Memory-assisted radio frequency beam training for mimo channels

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
GR01 Patent grant
GR01 Patent grant