CN112490688B - Determination method, determination device and computer storage medium - Google Patents
Determination method, determination device and computer storage medium Download PDFInfo
- Publication number
- CN112490688B CN112490688B CN202011416510.9A CN202011416510A CN112490688B CN 112490688 B CN112490688 B CN 112490688B CN 202011416510 A CN202011416510 A CN 202011416510A CN 112490688 B CN112490688 B CN 112490688B
- Authority
- CN
- China
- Prior art keywords
- matrix
- transformed
- received signal
- equivalent channel
- weight
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 56
- 239000011159 matrix material Substances 0.000 claims abstract description 522
- 230000001131 transforming effect Effects 0.000 claims abstract description 16
- 230000009466 transformation Effects 0.000 claims description 16
- 238000004422 calculation algorithm Methods 0.000 claims description 10
- 238000004891 communication Methods 0.000 claims description 8
- 230000005540 biological transmission Effects 0.000 claims description 7
- 238000010276 construction Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 description 38
- 238000010586 diagram Methods 0.000 description 11
- 238000012545 processing Methods 0.000 description 8
- 238000004590 computer program Methods 0.000 description 7
- 238000000354 decomposition reaction Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 4
- 230000007774 longterm Effects 0.000 description 4
- 238000007792 addition Methods 0.000 description 3
- 230000005291 magnetic effect Effects 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 2
- 238000011426 transformation method Methods 0.000 description 2
- 230000002087 whitening effect Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005562 fading Methods 0.000 description 1
- 230000005294 ferromagnetic effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01Q—ANTENNAS, i.e. RADIO AERIALS
- H01Q21/00—Antenna arrays or systems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Radio Transmission System (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
The embodiment of the application discloses a determination method, which comprises the following steps: after receiving the received signal, obtaining an equivalent channel matrix, transforming the equivalent channel matrix to obtain a transformed matrix, when i is greater than or equal to 1 and less than or equal to Nt, obtaining a matrix composed of an ith column element of the transformed matrix and a front i column element of the transformed matrix, calculating an ith column element of the iterative matrix according to the matrix composed of the ith column element and the front i column element based on an iterative formula of the constructed weight matrix of the received signal, updating i to i+1, returning to obtain the matrix composed of the ith column element of the transformed matrix and the front i column element of the transformed matrix when i is greater than or equal to 1 and less than or equal to Nt, and determining the weight matrix of the received signal according to the iterative matrix when i is greater than Nt so as to estimate the transmitted signal of the transmitting device. The embodiment of the application also provides equipment and a computer storage medium.
Description
Technical Field
The present invention relates to a technique for determining a weight matrix of a received signal in a large-scale antenna technique, and in particular, to a determining method, apparatus, and computer storage medium.
Background
Massive multiple input multiple output (Massive MIMO, large Scale Multi Input Multi Output) is one of the key technologies for long term evolution (LTE, long Term Evolution) and New air interface (NR, new Radio).
MIMO refers to that multiple antennas are used at the transmitting end and the receiving end to transmit and receive signals, and space resources are used to improve the channel capacity of the system and the quality of the received signals without increasing video resources and antenna transmit power.
In a MIMO system, N is assumed at the transmitting end t A root transmitting antenna, a receiving end with N r The received signal can be expressed by the following terms when the channel is quasi-static flat fading in one time slot:
y＝Hx+n (1)
wherein,,
is a channel matrix>
Is the transmit signal vector, ">
Is the received signal vector, ">
Is an additive complex gaussian white noise vector that is uncorrelated with the transmitted signal.
Let the receiving end's estimate of the transmitted signal be expressed as:
wherein,,
is a weight matrix.
Then the estimation error can be expressed as:
the weight matrix can be derived according to minimum mean square error (MMSE, minimum Mean Square Error) criteria as follows:
W MMSE ＝arg min W E[||Wy-x|| 2 ]＝R x H H (HR x H H +R n ) -1 (4)
wherein,,
is the autocorrelation matrix of the transmitted signal vector x, < >>
Is the autocorrelation matrix of the noise vector n.
Assuming that the transmitted signal symbols are independent of each other and equal in energy, their autocorrelation matrix can be expressed as:
assuming that the noise of each receiving antenna is independent of each other, but the noise energy may be unequal, there are:
substituting the above formula (4) can be simplified as:
W MMSE ＝H H (HH H +R) -1 (7)
wherein,,
in practical applications, the receiving end calculates the weight matrix W by Cholesky decomposition or LDLT decomposition to estimate the transmission signal MMSE However, the method has more calculation steps, so that the calculation is more complex, and the operation amount is larger; therefore, the existing weight matrix determination method has the technical problem of large calculated amount.
Disclosure of Invention
The embodiment of the application provides a determining method, determining equipment and a computer storage medium, which can improve the accuracy of video classification.
The technical scheme of the application is realized as follows:
the embodiment of the application provides a determining method, which is applied to receiving equipment, and the number of antennas of transmitting equipment is Nt, and comprises the following steps:
after receiving the received signal, obtaining an equivalent channel matrix;
transforming the equivalent channel matrix to obtain a transformed matrix;
when i is greater than or equal to 1 and less than or equal to Nt, acquiring a matrix formed by an ith column element of the transformed matrix and a front i column element of the transformed matrix; wherein, the initial value of i is 1;
according to a matrix formed by the ith column element and the previous ith column element, calculating the ith column element of the iterative matrix based on an iterative formula of the constructed weight matrix of the received signal;
i is updated to be i+1, and the matrix formed by the ith column element of the transformed matrix and the front ith column element of the transformed matrix is acquired when i is more than or equal to 1 and less than or equal to Nt;
and when i is greater than Nt, determining a weight matrix of the received signal according to the iteration matrix so as to estimate the transmitted signal of the transmitting device.
The embodiment of the application provides a receiving device, the number of antennas of a sending device is Nt, including:
the first acquisition module is used for acquiring an equivalent channel matrix after receiving the received signal;
the transformation module is used for transforming the equivalent channel matrix to obtain a transformed matrix;
the second acquisition module is used for acquiring a matrix formed by the ith column element of the transformed matrix and the front i column element of the transformed matrix when i is more than or equal to 1 and less than or equal to Nt; wherein, the initial value of i is 1;
the iteration module is used for calculating the ith column element of the iteration matrix based on an iteration formula of the constructed weight matrix of the received signal according to the matrix formed by the ith column element and the previous ith column element;
the updating module is used for updating i into i+1, and returning to execute the matrix formed by the ith column element of the transformed matrix and the front i column element of the transformed matrix when i is more than or equal to 1 and less than or equal to Nt;
and the determining module is used for determining a weight matrix of the received signal according to the iteration matrix when i is greater than Nt so as to estimate the transmitted signal of the transmitting device.
The embodiment of the application also provides a receiving device, the number of antennas of the sending device is Nt, and the receiving device includes: a processor and a storage medium storing instructions executable by the processor, the storage medium performing operations in dependence upon the processor through a communication bus, the instructions, when executed by the processor, performing the method of determining of one or more of the embodiments described above.
Embodiments of the present application provide a computer storage medium storing executable instructions that, when executed by one or more processors, perform the determination method described in one or more embodiments above.
The embodiment of the application provides a determining method, a device and a computer storage medium, wherein the method is applied to receiving equipment, the number of antennas of transmitting equipment is Nt, and the method comprises the following steps: after receiving a received signal, acquiring an equivalent channel matrix, transforming the equivalent channel matrix to obtain a transformed matrix, when i is greater than or equal to 1 and less than or equal to Nt, acquiring a matrix formed by an ith column element of the transformed matrix and a front i column element of the transformed matrix, calculating the ith column element of the iterative matrix according to the matrix formed by the ith column element and the front i column element based on an iterative formula of the constructed weight matrix of the received signal, updating i to i+1, and returning to acquire the ith column element of the transformed matrix and the matrix formed by the front i column element of the transformed matrix when i is greater than or equal to 1 and less than or equal to Nt, and determining the weight matrix of the received signal according to the iterative matrix when i is greater than Nt so as to estimate the transmitted signal of the transmitting device; that is, in the embodiment of the present application, after receiving a received signal, the equivalent channel matrix is transformed first, and then elements in each column of the iteration matrix are calculated column by using the transformed matrix and the constructed iteration formula, so as to obtain an iteration matrix, and finally, a weight matrix of the received signal is determined according to the iteration matrix, so that a transmission signal is estimated, so that elements in each column of the iteration formula can be calculated column by using the constructed iteration formula, and the iteration matrix can be calculated by using a parallel calculation mode.
Drawings
FIG. 1 is a flow chart of an alternative determination method according to an embodiment of the present application;
FIG. 2 is a flow chart of an example of an alternative determination method provided by embodiments of the present application;
fig. 3 is a schematic structural diagram of a receiving device according to an embodiment of the present application;
fig. 4 is a schematic structural diagram of a receiving device according to an embodiment of the present application.
Detailed Description
The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application.
Example 1
The embodiment of the application provides a determining method, which is applied to a receiving device, where the number of antennas of the sending device is Nt, and fig. 1 is a schematic flow diagram of an alternative determining method provided in the embodiment of the application, and referring to fig. 1, the determining method may include:
s101: after receiving the received signal, obtaining an equivalent channel matrix;
the existing calculation method for solving the weight matrix based on Cholesky decomposition or LDLT decomposition is more complex in steps and relatively large in operation amount, wherein elements on diagonal lines of the matrix need to be solved when the calculation method is used for solving, then other elements in the row or the column need to be solved, and the dependence causes large calculation amount of the weight matrix.
In order to reduce the calculation amount of the solution weight matrix, the present application provides a determining method, which is applied to a receiving device, wherein the receiving device may be a terminal device or a network device, and embodiments of the present application are not limited in detail herein.
Here, the transmitting device transmits a transmission signal, the receiving device receives a reception signal, and the receiving device receives due to the influence of a network path between the receiving device and the transmitting deviceThe received signal and the transmitted signal have a deviation, and in order to estimate the transmitted signal as much as possible, a weight matrix W of the received signal needs to be determined MMSE The receiving device then uses the received signal and the weight matrix W of the received signal MMSE Being able to estimate the estimated value of the transmitted signal, it is necessary to determine the weight matrix W of the received signal after the received signal is received MMSE 。
To determine the weight matrix W of the received signal MMSE Here, the receiving device first acquires the equivalent channel matrix H, and in order to acquire the equivalent channel matrix H, in an alternative embodiment, S101 may include:
after receiving the received signal, obtaining a channel matrix;
and determining the channel matrix as an equivalent channel matrix.
Specifically, after receiving a received signal, the receiving device acquires a channel matrix, and directly uses the channel matrix as an equivalent channel matrix, and determines a weight matrix based on the equivalent channel matrix; yet another alternative embodiment, S101 may include:
after receiving the received signal, obtaining a channel matrix and a preset matrix;
and determining the product of the channel matrix and the preset matrix as an equivalent channel matrix.
That is, after receiving the reception signal, a channel matrix and a preset matrix are acquired, wherein the preset matrix includes: a beam forming matrix and/or a precoding matrix; then, the product of the channel matrix and the beamforming matrix may be used as an equivalent channel matrix, the product of the channel matrix and the precoding matrix may be used as an equivalent channel matrix, and the product of the channel matrix, the beamforming matrix and the precoding matrix may be used as an equivalent channel matrix.
Taking 5G NR (5G New Radio) and long term evolution (LTE, long Term Evolution) systems as an example, the transmitted data symbols are usually mapped to different layers, and then are transmitted through actual physical antennas after operations such as precoding and beamforming, where the received signals can be expressed as
y＝HP 2 P 1 x+n (8)
Wherein,,
is a channel matrix>
Is a beamforming matrix, < >>
Is a precoding matrix, & gt>
Is the transmit signal vector, ">
Is the received signal vector, ">
Is the noise vector, N p Is the number of antenna ports, N l Is the number of transport layers.
At this time, the equivalent channel matrix may be made to be H e ＝HP 2 P 1 Transforming the received signal to y=h e x+n。
S102: transforming the equivalent channel matrix to obtain a transformed matrix;
after the equivalent channel matrix is obtained, the equivalent channel matrix is required to be transformed to obtain a transformed matrix, wherein the value of each column element of the iteration matrix can be calculated through the transformed matrix, so that the iteration matrix can be calculated, and the weight matrix can be finally determined.
S103: when i is greater than or equal to 1 and less than or equal to Nt, acquiring a matrix formed by an ith column element of the transformed matrix and a front i column element of the transformed matrix;
s104: according to a matrix formed by the ith column element and the former ith column element, calculating the ith column element of the iterative matrix based on an iterative formula of the constructed weight matrix of the received signal;
s105: i is updated to i+1, and the process returns to S103;
s106: when i is greater than Nt, a weight matrix of the received signal is determined according to the iteration matrix to estimate the transmitted signal of the transmitting device.
Specifically, determining an initial value of i, wherein the initial value of i is 1, acquiring a matrix composed of a transformed matrix 1 st column element and a transformed matrix 1 st column element, substituting the matrix composed of the transformed matrix 1 st column element and the transformed matrix 1 st column element into a constructed iteration formula, and calculating to obtain a value of the iteration matrix 1 st column element; and updating i to 2, and as 2 is greater than or equal to 1 and less than Nt, acquiring a matrix composed of the 2 nd column elements of the matrix after transformation and the first 2 column elements of the matrix after transformation, substituting the matrix composed of the 2 nd column elements of the matrix after transformation and the first 2 column elements of the matrix after transformation into an iteration formula to obtain the value of the 2 nd column elements of the iteration matrix, and so on, when i is greater than Nt, obtaining the iteration matrix of Nt columns, determining the weight matrix of the received signal according to the iteration matrix, and estimating the transmitted signal through the received signal and the weight matrix.
Further, in order to obtain an iteration matrix of the value matrix, the equivalent channel matrix is transformed first, where there may be multiple transformation modes, in order to reduce the calculation amount, it is convenient to obtain values of each column element of the iteration matrix in an iteration formula of the weight matrix, and in an alternative embodiment, the following formula is used to transform the equivalent channel matrix to obtain a transformed matrix:
wherein H is an equivalent channel matrix,
is a transformed matrix;
accordingly, the weight matrix of the received signal is determined using the following formula:
wherein,,
as an iterative matrix, W MMSE Is a weight matrix of the received signal.
Specifically, the conjugate transpose of the equivalent channel matrix is calculated, and the conjugate transpose is determined as the transformed matrix, so that an iteration formula is called according to the conjugate transpose of the equivalent channel matrix to obtain an iteration matrix, and the transformed matrix is the conjugate transpose of the equivalent channel matrix, so that the conjugate transpose of the iteration matrix is calculated again, and the conjugate transpose of the iteration matrix is determined as the weight matrix of the received signal.
In order to further reduce the amount of computation and to reduce the complexity of the operation, in an alternative embodiment, the equivalent channel matrix is transformed using the following formula to obtain a transformed matrix:
wherein,,
diagonal matrix of Nr rows and Nr columns, H is equivalent channel matrix, +.>
Is a transformed matrix;
accordingly, the weight matrix of the received signal is determined using the following formula:
wherein,,
as an iterative matrix, W MMSE A weight matrix for the received signal;
wherein Nr is the number of antennas of the receiving device;
that is, here the first diagonal matrix to be preset
And after the iteration matrix is calculated, in the process of determining the weight matrix according to the iteration matrix, corresponding transformation needs to be carried out on the iteration matrix, namely, the product of the iteration matrix and a preset first diagonal matrix is calculated, and the product is determined as the weight matrix.
Here, transforming the equivalent channel matrix may reduce the dimensions of the matrix in the expression of the weight matrix in the calculation process, thereby reducing the calculation amount.
To further reduce the complexity of the calculation, different calculation modes may be selected for different numbers of antennas, and in an alternative embodiment, S102 may include:
when Nr is equal to Nt, the equivalent channel matrix is transformed by adopting the following formula (9) to obtain a transformed matrix;
wherein H is an equivalent channel matrix,
is a transformed matrix;
correspondingly, determining a weight matrix of the received signal by adopting the following formula (10);
wherein,,
as an iterative matrix, W MMSE A weight matrix for the received signal;
that is, when the number of antennas of the receiving apparatus is the same as the number of antennas of the transmitting apparatus, the equivalent channel matrix is a square matrix, and at this time, the transformation method used is to calculate the conjugate transpose matrix of the equivalent channel matrix, and the conjugate transpose matrix of the equivalent channel matrix is still a square matrix, and the number of columns of the square matrix is Nr/Nt columns, so that the calculation amount can be reduced by adopting the transformation method.
When Nr is larger than Nt, the equivalent channel matrix is transformed by adopting the following formula (11) to obtain a transformed matrix;
wherein,,
diagonal matrix of Nr rows and Nr columns, H is equivalent channel matrix, +.>
Is a transformed matrix;
correspondingly, determining a weight matrix of the received signal by adopting the following formula (12);
wherein,,
as an iterative matrix, W MMSE A weight matrix for the received signal;
where Nr is the number of antennas of the receiving device.
For the case that the number of antennas of the receiving device is greater than that of antennas of the transmitting device, calculating the product of the first diagonal matrix and the equivalent channel matrix, wherein the column number of the matrix obtained by the product is Nt column, if a transformation mode of a conjugate transposed matrix is selected, the column number of the matrix after transformation is Nr column, and obviously, when Nr is greater than Nt, the calculation amount can be further corrected by adopting the transformation mode of the product of the first diagonal matrix and the equivalent channel matrix according to the calculation method of an iteration formula.
In order to determine the weight matrix, an iterative formula of the weight matrix needs to be pre-constructed, and in an alternative embodiment, the method further includes:
constructing a matrix M by using an equivalent channel matrix and a preset diagonal matrix;
determining an expression of a generalized inverse matrix of the matrix M;
substituting the expression of the generalized inverse matrix of the matrix M into a preset iterative algorithm to determine an iterative formula of the weight matrix.
Specifically, after determining the equivalent channel matrix, the matrix M may be constructed using the equivalent channel matrix and a preset diagonal matrix, where the preset diagonal matrix is a diagonal matrix of Nr rows and Nr columns; then, determining an expression of a generalized inverse matrix of the matrix M, wherein the expression of the generalized inverse matrix of the matrix M comprises an expression of a weight matrix; then, the expression is substituted into a preset iterative algorithm, for example, a preset Greville iterative algorithm, so that an expression of the weight matrix, that is, an iterative formula of the weight matrix, can be deduced.
In order to further reduce the calculation amount, different methods are adopted to construct the matrix M for different numbers of antennas, and in an alternative embodiment, the matrix M is constructed by using an equivalent channel matrix and a preset diagonal matrix, including:
when Nr is greater than Nt, a matrix M is constructed using the following formula:
wherein,,
for the first diagonal matrix of the preset diagonal matrices,/i>
Is Nt order identity matrix;
accordingly, determining the generalized inverse of the matrix M includes:
transforming the expression of the weight matrix to obtain a transformed expression of the weight matrix;
determining an expression of a generalized inverse matrix of the matrix M;
the generalized inverse matrix of the matrix M contains the transformed expression of the weight matrix.
Specifically, for the case that Nr is greater than Nt, the first row and the first column of the matrix M are the product of the first diagonal matrix and the equivalent channel matrix, and the second row and the first column are Nt-order identity matrices, after the matrix M is constructed, the expression of the weight matrix is transformed to obtain the expression of the weight matrix transformed, and then the generalized inverse matrix of the matrix M is calculated, so that the generalized inverse matrix of the matrix M contains the transformed expression of the weight matrix; in this way, the iteration formula of the weight matrix can be deduced through the generalized inverse matrix expression of the matrix M and a preset iteration algorithm.
For example, the expression of the weight matrix is as in the above formula (7), assuming that
And has
Construction matrix
Wherein->
The conjugate transpose of its matrix M is +.>
The weight matrix may be transformed into:
W MMSE ＝H H (M H M) -1 (14)
whereas the generalized inverse matrix (Moore-Penrose) of matrix M can be expressed as:
order the
The comparison can be obtained:
thus, W can be obtained by solving the generalized inverse of M MMSE The generalized inverse matrix can be calculated by using a Greville iterative algorithm, and the specific formula is as follows:
c k ＝m k -M k-1 d k
wherein m is k Is the kth column of M, M k Is a matrix of the first k columns of M.
Order the
Will->
Substituting the above iterative formula and according to +.>
Is characterized by a diagonal matrix, W MMSE The calculation process of (2) can be simplified as follows:
wherein,,
is->
Column k, < >>
Is->
A matrix of the first k columns of (b).
Thus, an iterative formula of the weight matrix can be obtained.
Additionally, in an alternative embodiment, constructing the matrix M using the equivalent channel matrix and a preset diagonal matrix includes:
when Nr is equal to Nt, the matrix M is constructed using the following formula:
wherein,,
is the second diagonal matrix of the preset diagonal matrices.
Specifically, for the case that Nr is equal to Nt, the first row and the first column of the matrix M are conjugate transpose matrices of the equivalent channel matrix, the second row and the first column are second diagonal matrices, and after the matrix M is constructed, a generalized inverse matrix of the matrix M is calculated, so that the generalized inverse matrix of the matrix M includes an expression of the weight matrix; in this way, the iteration formula of the weight matrix can be deduced through the generalized inverse matrix expression of the matrix M and a preset iteration algorithm.
For example, the expression of the weight matrix is as in the above formula (7), assuming that
And has
Construction matrix
Wherein->
Is N t The order identity matrix, the weight matrix may be transformed into:
comparing generalized inverse matrix calculation formula of M:
order the
The comparison can be obtained:
order the
Will->
Substituting Greville iteration formula, W MMSE The calculation process of (2) can be simplified as follows:
wherein,,
is->
Column k, < >>
Is->
A matrix of the first k columns of (b).
The following examples illustrate the determination of one or more of the embodiments described above.
Fig. 2 is a flowchart of an example of an alternative determining method provided in an embodiment of the present application, as shown in fig. 2, where the determining method may include:
s201: judging whether Nr is equal to Nt; if yes, executing S202, if no, executing S203;
s202: when N is r ＝N t Time, order
S204 is executed;
s203: when N is r ＞N t Time, order
S204 is executed;
s204: n is carried out t And (3) carrying out iterative calculation:
specifically, the transformed equivalent channel matrix, diagonal matrix and identity matrix are substituted into the above formula (12) or the above formula (16) to obtain an iterative matrix
S205: when N is r ＝N t In the time-course of which the first and second contact surfaces,
s206: when N is r ＞N t In the time-course of which the first and second contact surfaces,
in addition, since the weight matrix is used to estimate the transmitted signal, the calculation is also performed
Therefore, in the concrete implementation, usually +.>
The received signal y is whitened to obtain +.>
To reduce the amount of computation. In the following algorithm complexity analysis, for comparative rationality, the statistics are also +.>
Is not the operand of (a)
Is a calculation amount of (a).
Wherein, N is as described above r ＝N t The calculation method is a direct calculation method, and N is as follows r ＞N t Is a whitening calculation method.
Wherein the complex multiplication is converted into 4 real multiplications and 2 real additions, the complex addition is converted into 2 real additions, which can be deduced from the above table 1, when N r ＝N t In this case, the complexity of the direct calculation method is minimal, when N r ＞N t In this case, the whitening calculation method of this example has the smallest complexity, and the calculation complexity of this example is lower than that of the corresponding Cholesky decomposition method.
In the example, a matrix M is respectively constructed by two methods, the calculation of a weight matrix is converted into generalized inverse matrix calculation, and the iterative calculation formula is simplified according to the characteristics of the constructed matrix M; according to the number of the receiving and transmitting antennas, different weight matrix calculation methods, namely a method with lower corresponding calculation complexity, are applied, and the two methods can multiplex a main iterative calculation part through different preprocessing and post-processing operations; according to the collectionThe method for selecting the number of the transmitting antennas to have lower calculation complexity, wherein the calculation complexity is lower than that of the traditional Cholesky decomposition method, and the calculation process can be completely carried out according to a row vector or a column vector (the vector length is N r ) And parallel processing is performed, so that the method is very suitable for parallel calculation to improve the operation efficiency and reduce the operation delay.
Example two
Fig. 3 is a schematic structural diagram of a receiving device provided in an embodiment of the present application, as shown in fig. 3, where the number of antennas of a transmitting device is Nt, including:
a first obtaining module 31, configured to obtain an equivalent channel matrix after receiving a received signal;
a transformation module 32, configured to transform the equivalent channel matrix to obtain a transformed matrix;
a second obtaining module 33, configured to obtain, when i is greater than or equal to 1 and less than or equal to Nt, a matrix formed by an i-th column element of the transformed matrix and a front i-column element of the transformed matrix; wherein, the initial value of i is 1;
the iteration module 34 is configured to calculate, according to a matrix formed by the ith column element and the previous ith column element, the ith column element of the iteration matrix based on an iteration formula of the constructed weight matrix of the received signal;
an updating module 35, configured to update i to i+1, and return to obtain a matrix composed of an i-th column element of the transformed matrix and a front i-column element of the transformed matrix when i is greater than or equal to 1 and less than or equal to Nt;
a determining module 36, configured to determine a weight matrix of the received signal according to the iteration matrix when i is greater than Nt, so as to estimate the transmission signal of the transmitting device.
Optionally, the transformation module 32 is specifically configured to:
transforming the equivalent channel matrix by adopting the following formula to obtain a transformed matrix:
wherein H is an equivalent channel matrix,
is a transformed matrix;
accordingly, the determining module 36 is specifically configured to:
the weight matrix of the received signal is determined using the following formula:
wherein,,
as an iterative matrix, W MMSE Is a weight matrix of the received signal.
Optionally, the transformation module 32 is specifically configured to:
transforming the equivalent channel matrix by adopting the following formula to obtain a transformed matrix:
wherein,,
diagonal matrix of Nr rows and Nr columns, H is equivalent channel matrix, +.>
Is a transformed matrix;
accordingly, the determining module 36 is specifically configured to:
the weight matrix of the received signal is determined using the following formula:
wherein,,
as an iterative matrix, W MMSE A weight matrix for the received signal;
where Nr is the number of antennas of the receiving device.
Optionally, the transformation module 32 is specifically configured to:
when Nr is equal to Nt, the equivalent channel matrix is transformed by adopting the following formula to obtain a transformed matrix:
wherein H is an equivalent channel matrix,
is a transformed matrix;
accordingly, the determining module 36 is specifically configured to:
the weight matrix of the received signal is determined using the following formula:
wherein,,
as an iterative matrix, W MMSE A weight matrix for the received signal;
the transformation module 32 is specifically configured to:
when Nr is greater than Nt, the equivalent channel matrix is transformed by adopting the following formula to obtain a transformed matrix:
wherein,,
diagonal matrix of Nr rows and Nr columns, H being equivalent channel momentArray (S)>
Is a transformed matrix;
accordingly, the determining module 36 is specifically configured to:
the weight matrix of the received signal is determined using the following formula:
wherein,,
as an iterative matrix, W MMSE A weight matrix for the received signal;
where Nr is the number of antennas of the receiving device.
Optionally, the first obtaining module 31 is specifically configured to:
after receiving the received signal, obtaining a channel matrix;
and determining the channel matrix as an equivalent channel matrix.
Optionally, the receiving device further includes:
the construction module is specifically used for:
constructing a matrix M by using an equivalent channel matrix and a preset diagonal matrix; the preset diagonal matrix is a diagonal matrix of Nr rows and Nr columns;
determining an expression of a generalized inverse matrix of the matrix M; wherein the expression of the generalized inverse matrix of the matrix M comprises the expression of the weight matrix;
substituting the expression of the generalized inverse matrix of the matrix M into a preset iterative algorithm to determine an iterative formula of the weight matrix.
Optionally, the constructing module constructs the matrix M by using the equivalent channel matrix and a preset diagonal matrix, including:
when Nr is greater than Nt, a matrix M is constructed using the following formula:
wherein,,
for the first diagonal matrix of the preset diagonal matrices,/i>
Is Nt order identity matrix;
correspondingly, the construction module determines that the generalized inverse matrix of the matrix M comprises:
transforming the expression of the weight matrix to obtain a transformed expression of the weight matrix;
determining an expression of a generalized inverse matrix of the matrix M; the generalized inverse matrix of the matrix M contains the transformed expression of the weight matrix.
Optionally, the constructing module constructs the matrix M by using the equivalent channel matrix and a preset diagonal matrix, including:
when Nr is equal to Nt, the matrix M is constructed using the following formula:
wherein,,
a second diagonal matrix in the preset diagonal matrices;
the second diagonal matrix is an inverse matrix of the first diagonal matrix in the preset diagonal matrices.
Optionally, the first obtaining module 31 is specifically configured to:
after receiving the received signal, obtaining a channel matrix and a preset matrix; the preset matrix comprises the following components: a beam forming matrix and/or a precoding matrix;
and determining the product of the channel matrix and the preset matrix as an equivalent channel matrix.
In practical applications, the first obtaining module 31, the transforming module 32, the second obtaining module 33, the iterating module 34, the updating module 35, the determining module 36 and the constructing module may be implemented by a processor located on the receiving device, specifically, a central processing unit (CPU, central Processing Unit), a microprocessor (MPU, microprocessor Unit), a digital signal processor (DSP, digital Signal Processing) or a field programmable gate array (FPGA, field Programmable Gate Array).
Fig. 4 is a second schematic structural diagram of a receiving device according to an embodiment of the present application, as shown in fig. 4, an embodiment of the present application provides a receiving device 400, including:
a processor 41 and a storage medium 42 storing instructions executable by the processor 41, the storage medium 42 performing operations in dependence on the processor 41 through a communication bus 43, the instructions, when executed by the processor 41, performing the determination method of the above embodiment.
In practical use, the components in the terminal are coupled together via the communication bus 43. It is understood that the communication bus 43 is used to enable connected communication between these components. The communication bus 43 includes a power bus, a control bus, and a status signal bus in addition to the data bus. But for clarity of illustration the various buses are labeled as communication bus 43 in fig. 4.
The present application provides a computer storage medium storing executable instructions that, when executed by one or more processors, perform the method of determining of embodiment one.
The computer readable storage medium may be a magnetic random access Memory (ferromagnetic random access Memory, FRAM), read Only Memory (ROM), programmable Read Only Memory (Programmable Read-Only Memory, PROM), erasable programmable Read Only Memory (Erasable Programmable Read-Only Memory, EPROM), electrically erasable programmable Read Only Memory (Electrically Erasable Programmable Read-Only Memory, EEPROM), flash Memory (Flash Memory), magnetic surface Memory, optical disk, or compact disk Read Only Memory (Compact Disc Read-Only Memory, CD-ROM).
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 a hardware embodiment, a 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, magnetic disk storage, 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 embodiments of 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.
The foregoing description is only of the preferred embodiments of the present application and is not intended to limit the scope of the present application.
Claims (11)
1. A determining method, wherein the method is applied to a receiving device, and the number of antennas of a transmitting device is Nt, and the method comprises:
after receiving the received signal, obtaining an equivalent channel matrix;
transforming the equivalent channel matrix to obtain a transformed matrix;
when i is greater than or equal to 1 and less than or equal to Nt, acquiring a matrix formed by an ith column element of the transformed matrix and a front i column element of the transformed matrix; wherein, the initial value of i is 1;
according to a matrix formed by the ith column element and the previous ith column element, calculating the ith column element of the iterative matrix based on an iterative formula of the constructed weight matrix of the received signal;
i is updated to be i+1, and the matrix formed by the ith column element of the transformed matrix and the front ith column element of the transformed matrix is acquired when i is more than or equal to 1 and less than or equal to Nt;
when i is greater than Nt, determining a weight matrix of the received signal according to the iteration matrix so as to estimate a transmission signal of the transmission device;
wherein the method comprises the following steps:
constructing a matrix M by using the equivalent channel matrix and a preset diagonal matrix; wherein the preset diagonal matrix is a diagonal matrix of Nr rows and Nr columns;
determining an expression of a generalized inverse matrix of the matrix M; wherein the expression of the generalized inverse matrix of the matrix M comprises the expression of the weight matrix;
substituting the expression of the generalized inverse matrix of the matrix M into a preset iterative algorithm, and determining the iterative formula of the weight matrix.
2. The method of claim 1, wherein the equivalent channel matrix is transformed using the formula:
wherein H is an equivalent channel matrix,
is a transformed matrix;
accordingly, the weight matrix of the received signal is determined using the following formula:
wherein,,
for the iterative matrix, W MMSE Is a weight matrix of the received signal.
3. The method of claim 1, wherein the equivalent channel matrix is transformed using the formula:
wherein,,
diagonal matrix of Nr rows and Nr columns, H is equivalent channel matrix, +.>
Is a transformed matrix;
accordingly, the weight matrix of the received signal is determined using the following formula:
wherein,,
for the iterative matrix, W MMSE A weight matrix for the received signal;
and Nr is the number of antennas of the receiving device.
4. A method according to any one of claim 1 to 3, wherein,
when Nr is equal to Nt, transforming the equivalent channel matrix by using the following formula to obtain a transformed matrix:
wherein H is an equivalent channel matrix,
is a transformed matrix;
accordingly, the weight matrix of the received signal is determined using the following formula:
wherein,,
for the iterative matrix, W MMSE A weight matrix for the received signal;
when Nr is larger than Nt, the equivalent channel matrix is transformed by adopting the following formula to obtain a transformed matrix:
wherein,,
diagonal matrix of Nr rows and Nr columns, H is equivalent channel matrix, +.>
Is a transformed matrix;
accordingly, the weight matrix of the received signal is determined using the following formula:
wherein,,
for the iterative matrix, W MMSE A weight matrix for the received signal;
and Nr is the number of antennas of the receiving device.
5. A method according to any one of claims 1 to 3, wherein said obtaining an equivalent channel matrix after receiving a received signal comprises:
after receiving the received signal, obtaining a channel matrix;
and determining the channel matrix as the equivalent channel matrix.
6. The method of claim 5, wherein constructing the matrix M using the equivalent channel matrix and a predetermined diagonal matrix comprises:
when Nr is greater than Nt, a matrix M is constructed using the following formula:
wherein,,
for the first diagonal matrix of the preset diagonal matrices,/a first diagonal matrix of the preset diagonal matrices>
Is Nt order identity matrix;
accordingly, the determining the generalized inverse matrix of the matrix M includes:
transforming the expression of the weight matrix to obtain a transformed expression of the weight matrix;
determining an expression of a generalized inverse matrix of the matrix M; the generalized inverse matrix of the matrix M comprises a transformed expression of the weight matrix.
7. The method of claim 6, wherein constructing the matrix M using the equivalent channel matrix and a predetermined diagonal matrix comprises:
when Nr is equal to Nt, the matrix M is constructed using the following formula:
wherein,,
and the second diagonal matrix is the second diagonal matrix in the preset diagonal matrices.
8. A method according to any one of claims 1 to 3, wherein said obtaining an equivalent channel matrix after receiving a received signal comprises:
after receiving the receiving signals, obtaining a channel matrix and a preset matrix; wherein the preset matrix comprises: a beam forming matrix and/or a precoding matrix;
and determining the product of the channel matrix and the preset matrix as the equivalent channel matrix.
9. A receiving apparatus, characterized in that the number of antennas of a transmitting apparatus is Nt, comprising:
the first acquisition module is used for acquiring an equivalent channel matrix after receiving the received signal;
the transformation module is used for transforming the equivalent channel matrix to obtain a transformed matrix;
the second acquisition module is used for acquiring a matrix formed by the ith column element of the transformed matrix and the front i column element of the transformed matrix when i is more than or equal to 1 and less than or equal to Nt; wherein, the initial value of i is 1;
the iteration module is used for calculating the ith column element of the iteration matrix based on an iteration formula of the constructed weight matrix of the received signal according to the matrix formed by the ith column element and the previous ith column element;
the updating module is used for updating i into i+1, and returning to execute the matrix formed by the ith column element of the transformed matrix and the front i column element of the transformed matrix when i is more than or equal to 1 and less than or equal to Nt;
a determining module, configured to determine a weight matrix of the received signal according to the iteration matrix when i is greater than Nt, so as to estimate a transmission signal of the transmitting device;
wherein the receiving apparatus further comprises a construction module for:
constructing a matrix M by using the equivalent channel matrix and a preset diagonal matrix; wherein the preset diagonal matrix is a diagonal matrix of Nr rows and Nr columns;
determining an expression of a generalized inverse matrix of the matrix M; wherein the expression of the generalized inverse matrix of the matrix M comprises the expression of the weight matrix;
substituting the expression of the generalized inverse matrix of the matrix M into a preset iterative algorithm, and determining the iterative formula of the weight matrix.
10. A receiving apparatus, characterized in that the number of antennas of a transmitting apparatus is Nt, the receiving apparatus comprising: a processor and a storage medium storing instructions executable by the processor, the storage medium performing operations in dependence on the processor through a communication bus, the instructions, when executed by the processor, performing the method of determining of any one of the preceding claims 1 to 8.
11. A computer storage medium storing executable instructions which, when executed by one or more processors, perform the method of determining of any one of claims 1 to 8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011416510.9A CN112490688B (en) | 2020-12-04 | 2020-12-04 | Determination method, determination device and computer storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011416510.9A CN112490688B (en) | 2020-12-04 | 2020-12-04 | Determination method, determination device and computer storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112490688A CN112490688A (en) | 2021-03-12 |
| CN112490688B true CN112490688B (en) | 2023-06-23 |
Family
ID=74939873
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011416510.9A Active CN112490688B (en) | 2020-12-04 | 2020-12-04 | Determination method, determination device and computer storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112490688B (en) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110581815A (en) * | 2019-10-28 | 2019-12-17 | 中国联合网络通信集团有限公司 | Signal detection method and device |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101888286A (en) * | 2009-05-11 | 2010-11-17 | 京信通信系统（中国）有限公司 | Decoding method of MIMO space multiplexing and device thereof |
| EP2474104B1 (en) * | 2009-09-04 | 2015-02-11 | Hitachi, Ltd. | Tomlinson harashima precoding with additional receiver processing in a multi-user multiple-input multiple-output wireless transmission system |
| US9729277B2 (en) * | 2015-12-14 | 2017-08-08 | Tsinghua University | Signal detecting method and device |
| US9467181B2 (en) * | 2015-03-12 | 2016-10-11 | Freescale Semiconductor, Inc. | Radio frequency receiver capable of determining a noise estimate in case of received power unbalanced antennas and method of operating thereof |
| CN108965172B (en) * | 2018-07-10 | 2021-06-08 | 电子科技大学 | Method for estimating uplink channel of large-scale MIMO system with phase noise |
| CN108736938B (en) * | 2018-09-07 | 2021-01-22 | 电子科技大学 | Method for large-scale MIMO uplink channel estimation and data demodulation |
| WO2020158966A1 (en) * | 2019-01-29 | 2020-08-06 | 엘지전자 주식회사 | Method for performing channel estimation for massive mimo channel, and communication device therefor |
| CN111162828B (en) * | 2019-12-12 | 2022-03-25 | 重庆邮电大学 | Low-complexity signal detection method of large-scale MIMO system |
| CN111970034A (en) * | 2020-08-27 | 2020-11-20 | 东南大学 | Distributed phased array MIMO channel acquisition method |
-
2020
- 2020-12-04 CN CN202011416510.9A patent/CN112490688B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110581815A (en) * | 2019-10-28 | 2019-12-17 | 中国联合网络通信集团有限公司 | Signal detection method and device |
Non-Patent Citations (2)
| Title |
|---|
| MediaTek Inc..R1-134451 "System modeling methodology for reduced-complexity ML (R-ML) receivers".3GPP tsg_ran\WG1_RL1.2013,(第TSGR1_74b期),全文. * |
| 依据列相关性优化高斯测量矩阵;边胜琴;徐正光;张利欣;;计算机测量与控制(第11期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112490688A (en) | 2021-03-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106209195B (en) | Channel state information acquisition method, channel state information feedback method and device | |
| TWI622276B (en) | Channel state information acquisition method, channel state information feedback method and device | |
| EP2332300B1 (en) | Method for precoding signals for transmission in wireless mimo system | |
| CN110545131B (en) | Antenna array design method and system under millimeter-wave line-of-sight MIMO channel | |
| TW201635736A (en) | Signal detection method and device | |
| JP2017528960A (en) | Data transmission method and apparatus | |
| CN112929128B (en) | MIMO detection method and device based on confidence propagation | |
| JP6641506B2 (en) | Method and device for transmitting data | |
| CN112490688B (en) | Determination method, determination device and computer storage medium | |
| Yu et al. | Blind identification of multi-channel ARMA models based on second-order statistics | |
| US10250307B2 (en) | Method for setting cooperation parameters of communication system, computer program thereof, information storage means thereof, and communication system including plurality of node devices | |
| CN108964795B (en) | Antenna reciprocity calibration method, device and system | |
| CN113839695B (en) | FDD large-scale MIMO and rate optimal statistical precoding method and device | |
| CN112803976B (en) | Large-scale MIMO precoding method and system and electronic equipment | |
| CN114839604A (en) | Orthogonal waveform design method and system for MIMO radar | |
| CN108366025B (en) | Signal synthesis method and system | |
| WO2017101097A1 (en) | Channel statistical information obtaining method and receiver | |
| CN100369390C (en) | Method of receiving radio transmission by multiple antenna array | |
| CN112019461A (en) | Channel prediction method, wireless communication system, and storage device | |
| Asadi et al. | A fast matrix completion method for index coding | |
| EP4091265B1 (en) | Improved internal data transfer in a multiple-antenna communication system | |
| CN110417692B (en) | Uplink channel tracking method and device | |
| CN115706684B (en) | Channel parameter estimation method, device and storage medium | |
| CN113810096B (en) | Uplink transmission method and system for large-scale MIMO satellite mobile communication | |
| CN113992482B (en) | Signal detection method and system based on 2PPJ iteration |
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 |