CN101834827B - Signal detection method and device in multiple-input multiple-output system - Google Patents

Signal detection method and device in multiple-input multiple-output system Download PDF

Info

Publication number
CN101834827B
CN101834827B CN2010101354912A CN201010135491A CN101834827B CN 101834827 B CN101834827 B CN 101834827B CN 2010101354912 A CN2010101354912 A CN 2010101354912A CN 201010135491 A CN201010135491 A CN 201010135491A CN 101834827 B CN101834827 B CN 101834827B
Authority
CN
China
Prior art keywords
layer
mrow
bit
nodes
msub
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.)
Expired - Fee Related
Application number
CN2010101354912A
Other languages
Chinese (zh)
Other versions
CN101834827A (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.)
DATANG LINKTECH INFOSYSTEM Co Ltd
Original Assignee
DATANG LINKTECH INFOSYSTEM Co Ltd
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 DATANG LINKTECH INFOSYSTEM Co Ltd filed Critical DATANG LINKTECH INFOSYSTEM Co Ltd
Priority to CN2010101354912A priority Critical patent/CN101834827B/en
Publication of CN101834827A publication Critical patent/CN101834827A/en
Application granted granted Critical
Publication of CN101834827B publication Critical patent/CN101834827B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

The embodiment of the invention discloses signal detection method and device in a multiple-input multiple-output system. The signal detection method comprises the following steps of: converting a constellation point of a quadrature amplitude modulation (QAM) transmitting signal into a bit vector weight sum form to obtain a bit-level-expressed QAM transmitting signal; converting a channel matrix into a compound channel matrix according to the bit-level-expressed QAM transmitting signal; carrying out QR decomposition on the compound channel matrix to obtain an upper triangular matrix; constructing a bit-by-bit layered tree structure by using the upper triangular matrix; searching the bit-by-bit layered tree structure layer by layer based on a breadth-first M algorithm to obtain a signal candidate set; calculating a measure value of each candidate signal in the signal candidate set; and calculating the posterior information of the transmitting signal by using the measure value. The embodiment of the invention can reduce the realization difficulty of a signal detection process.

Description

Signal detection method and device in multi-input multi-output system
Technical Field
The present application relates to the field of communications and computer technologies, and in particular, to a method and an apparatus for detecting signals in a mimo system.
Background
With the development of wireless communication technology, the MIMO (Multiple-input Multiple-output) technology using Multiple antenna elements fully utilizes space resources and can greatly improve the frequency utilization rate of the system within a limited bandwidth, so that the MIMO technology provides a capacity potential which cannot be provided by a single antenna technology for a user and becomes an effective way for improving the capacity and reliability of the system by using space dimensions. Currently, MIMO technology has become one of hot spot technologies in 3G or 4G.
In the MIMO system, since the performance of the receiving end will ultimately greatly affect the transmission rate, the error rate performance and the system complexity of the whole transmission system, the improvement of the performance of the receiving end becomes a focus of attention of researchers. At present, researchers have proposed applying a sequence decoding method to the signal detection process of the MIMO system. The tree searching method using the breadth-first M algorithm can search to obtain a signal candidate sequence, calculate the metric values of the signal candidate sequence, and calculate the posterior information of the bits corresponding to the transmitted signals by using the metric values to complete the signal detection process. The tree searching method based on the breadth-first M algorithm not only ensures that the complexity of calculation is not changed along with the change of the signal-to-noise ratio and the channel condition, but also ensures that a transmission system has good error rate performance.
However, the inventors found in their studies that: the tree searching method using the breadth-first M algorithm is based on the signal tree of the symbol level for searching, under the condition of higher constellation dimensionality, the calculation complexity of the whole tree searching process is high, and the realization difficulty of the signal detection process is increased finally.
Disclosure of Invention
In order to solve the above technical problem, embodiments of the present application provide a method and an apparatus for signal detection in a mimo system, so as to reduce implementation difficulty in a signal detection process.
The embodiment of the application discloses the following technical scheme:
a method of signal detection in a multiple-input multiple-output system, comprising: converting constellation points of Quadrature Amplitude Modulation (QAM) transmitting signals into a bit vector weighted sum form to obtain QAM transmitting signals represented by a bit level, and converting a channel matrix into a composite channel matrix according to the QAM transmitting signals represented by the bit level; performing QR decomposition on the composite channel matrix to obtain an upper triangular matrix, and constructing a bit-by-bit hierarchical tree structure by using the upper triangular matrix; searching the tree structure which is layered bit by bit on the basis of a breadth-first M algorithm layer by layer to obtain a signal candidate set, and calculating the metric value of each candidate signal in the signal candidate set; and calculating posterior information of the transmitting signal by using the metric value.
A signal detection apparatus in a multiple-input multiple-output system, comprising: the conversion unit is used for converting the constellation points of the Quadrature Amplitude Modulation (QAM) sending signals into a bit vector weighted sum form to obtain QAM sending signals represented by a bit level, and obtaining a composite channel matrix according to the QAM sending signals represented by the bit level; the decomposition unit is used for carrying out QR decomposition on the composite channel matrix to obtain an upper triangular matrix, and a bit-by-bit hierarchical tree structure is constructed by utilizing the upper triangular matrix; the searching unit is used for searching the bit-by-bit hierarchical tree structure layer by layer based on a breadth-first M algorithm to obtain a signal candidate set and calculating the metric value of each candidate signal in the signal candidate set; and the detection unit calculates the posterior information of the transmitting signal by using the metric value.
It can be seen from the above embodiments that the bit-level tree search detection method in the present application can reduce the branch paths that need to be searched for each time, and only 2M node metrics need to be calculated after each layer of nodes are expanded correspondingly. Compared with the existing tree searching method based on the symbol, in each stage of the tree searching, each retention path needs to be expanded into a range of constellation set according to the size of the constellation set
Figure GSA00000066275700021
The branch paths need to be calculated after the corresponding nodes at each layer are expanded
Figure GSA00000066275700022
Measurement of child nodes, and thus complexity and modulation order M of the entire tree search algorithmrIn an exponential relationshipBy the implementation scheme, the complexity of the whole algorithm is reduced, and therefore the implementation difficulty of the signal detection process is reduced.
Drawings
In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present application, and it is obvious for those skilled in the art that other drawings can be obtained according to the drawings without inventive exercise.
FIG. 1 is a flow chart of an embodiment of a method for signal detection in a multiple input multiple output system of the present application;
fig. 2 is a schematic diagram of gray mapping according to the present application;
FIG. 3 is a schematic illustration of a natural mapping according to the present application;
fig. 4 is an expanded search graph of the bit-level tree search algorithm when the retention path of the 16QAM constellation is 4 according to the present application;
FIG. 5 is a block diagram of an embodiment of a signal detection apparatus in a MIMO system according to the present application;
fig. 6 is a block diagram of another embodiment of a signal detection apparatus in a mimo system according to the present application.
Detailed Description
In order to make the aforementioned objects, features and advantages of the present application more comprehensible, embodiments accompanying the present application are described in detail below.
Example one
Please refer to fig. 1, which is a flowchart illustrating a signal detection method in a mimo system according to an embodiment of the present invention, the method includes the following steps:
step 101: converting constellation points of Quadrature Amplitude Modulation (QAM) transmitting signals into a bit vector weighted sum form to obtain QAM transmitting signals represented by a bit level, and converting a channel matrix into a composite channel matrix according to the QAM transmitting signals represented by the bit level;
step 102: performing QR decomposition on the composite channel matrix to obtain an upper triangular matrix, and constructing a bit-by-bit hierarchical tree structure by using the upper triangular matrix;
step 103: and searching the bit-by-bit hierarchical tree structure layer by applying a breadth-first M algorithm to obtain a signal candidate set, and calculating the metric value of each candidate signal in the signal candidate set.
Step 104: and calculating posterior information of the transmitting signal by using the metric value.
The above signal detection process is explained in detail below. When the constellation points of the QAM transmission signals are converted into a form represented by bit vectors to obtain QAM transmission signals represented by the bit vectors, and the channel matrix is converted into a composite channel matrix according to the QAM transmission signals represented by the bit vectors, the following procedure may be performed.
Suppose that one uses NtRoot transmitting antenna and NrIn the MIMO system model of the root receiving antenna, the sending symbol is quadrature amplitude modulation of a square constellation, and the size of the constellation set isWherein M iscFor the number of bits of each transmission symbol in the complex field, the discrete-time system model of the MIMO system in the complex field is:
yc=Hcsc+nc (1)
the system model in the complex domain of equation (1) is equivalent to a system model in the real domain as follows:
y=Hs+n (2)
it can be concluded that, in the real number domain, the vector dimension of the received signal y is NRX 1, vector dimension of the transmitted signal s is NTX 1, the vector dimension of the channel matrix H is NT×NTThe vector dimension of the real Gaussian noise N is NRX 1, mean value ofThe variance matrix is
Figure GSA00000066275700043
Wherein N isRIs the equivalent number of receiving antennas in real number domain, NTIs the equivalent number of transmitting antennas in the real number domain, and NR=2Nr,NT=2Nt
When the MIMO channel is a flat fading Rayleigh channel, N frames per frameTMrThe long bit symbols are:
<math> <mrow> <mi>x</mi> <mo>=</mo> <msup> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>x</mi> <mn>1,1</mn> </msub> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <msub> <mi>x</mi> <mrow> <mn>1</mn> <mo>,</mo> <msub> <mi>M</mi> <mi>r</mi> </msub> </mrow> </msub> </mtd> <mtd> <msub> <mi>x</mi> <mn>2,1</mn> </msub> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <msub> <mi>x</mi> <mrow> <mn>2</mn> <mo>,</mo> <msub> <mi>M</mi> <mi>r</mi> </msub> </mrow> </msub> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <msub> <mi>x</mi> <mrow> <msub> <mi>N</mi> <mi>T</mi> </msub> <mo>,</mo> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <msub> <mi>x</mi> <mrow> <msub> <mi>N</mi> <mi>T</mi> </msub> <mo>,</mo> <msub> <mi>M</mi> <mi>r</mi> </msub> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein M isrNumber of bits in real number domain for each transmitted symbol, and Mr=Mc/2. After serial-to-parallel conversion, the symbols are mapped into a transmission symbol s.
When the system uses gray mapping as shown in fig. 2, then the transmitted symbol s on the kth transmit antennakAnd its corresponding mapping bit symbol xk,j,j=1,…,MrThe relationship between can be represented by:
<math> <mrow> <msub> <mi>s</mi> <mi>k</mi> </msub> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>M</mi> <mi>r</mi> </msub> </munderover> <msub> <mi>w</mi> <mi>i</mi> </msub> <munderover> <mi>&Pi;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>i</mi> </munderover> <msub> <mi>x</mi> <mrow> <mi>k</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow> </math>
an example of an 8-PAM gray mapped constellation is given in fig. 2, where an intermediate parameter u is introducedk,i,uk,iE { +1, -1}, let,
<math> <mrow> <msub> <mi>u</mi> <mrow> <mi>k</mi> <mo>,</mo> <mi>i</mi> </mrow> </msub> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>i</mi> </munderover> <msub> <mi>x</mi> <mrow> <mi>k</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mn>1,2</mn> <mo>,</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>,</mo> <msub> <mi>M</mi> <mi>r</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow> </math>
when using ukTo represent uk,i,i=1,…,MrThe composed vector, denoted by w, corresponds to uk,iBy a weight factor wi,i=1,…,MrWhen the vector is formed, then,
<math> <mrow> <msub> <mi>u</mi> <mi>k</mi> </msub> <mo>=</mo> <msup> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>u</mi> <mrow> <mi>k</mi> <mo>,</mo> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>u</mi> <mrow> <mi>k</mi> <mo>,</mo> <mn>2</mn> </mrow> </msub> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <msub> <mi>u</mi> <mrow> <mi>k</mi> <mo>,</mo> <msub> <mi>M</mi> <mi>r</mi> </msub> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> <mo>,</mo> </mrow> </math> <math> <mrow> <mi>w</mi> <mo>=</mo> <msup> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>w</mi> <mn>1</mn> </msub> </mtd> <mtd> <msub> <mi>w</mi> <mn>2</mn> </msub> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <msub> <mi>w</mi> <msub> <mi>M</mi> <mi>r</mi> </msub> </msub> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow> </math>
thus, skThe vector of (d) is represented as:
sk=wTuk (10)
in the case of M-QAM, the,
the weight factor is:
<math> <mrow> <mi>w</mi> <mo>=</mo> <msqrt> <mfrac> <mn>3</mn> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <mi>M</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </mfrac> </msqrt> <msup> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msqrt> <mi>M</mi> </msqrt> <msup> <mn>2</mn> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>,</mo> </mtd> <mtd> <msqrt> <mi>M</mi> </msqrt> <msup> <mn>2</mn> <mrow> <mo>-</mo> <mn>2</mn> </mrow> </msup> <mo>,</mo> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>,</mo> </mtd> <mtd> <mn>1</mn> <mo>,</mo> </mtd> <mtd> <msqrt> <mi>M</mi> </msqrt> <msup> <mn>2</mn> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mi>j</mi> <mo>,</mo> </mtd> <mtd> <msqrt> <mi>M</mi> </msqrt> <msup> <mn>2</mn> <mrow> <mo>-</mo> <mn>2</mn> </mrow> </msup> <mi>j</mi> <mo>,</mo> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <mi>j</mi> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein j represents an imaginary symbol, i.e.
Figure GSA00000066275700052
M is the number of constellation points. When using a 16QAM constellation, for example, M-16,
Figure GSA00000066275700053
from uk,iCan be found ink,iAnd uk,iAnd uk,i-1The following relations exist between the following components:
xk,i=uk,iuk,i-1,i=2,…,Mr (12)
when i is 1, xk,i=uk,i. By substituting equation (12) into equation (7), the relationship between vectors s and u can be obtained:
<math> <mrow> <mi>s</mi> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>I</mi> <msub> <mi>N</mi> <mi>T</mi> </msub> </msub> <mo>&CircleTimes;</mo> <msup> <mi>w</mi> <mi>T</mi> </msup> <mo>)</mo> </mrow> <mi>u</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein,
Figure GSA00000066275700056
representing the Kronecker direct product.
When the formula (13) is brought into the formula (2), there are
Figure GSA00000066275700057
Converting constellation points of Quadrature Amplitude Modulation (QAM) transmitting signals into a bit vector weighted sum form to obtain QAM transmitting signals u represented by bit levels, and defining a composite channel matrix according to the relation between the QAM transmitting signals s and the bit vectors u:
<math> <mrow> <mi>A</mi> <mo>=</mo> <mi>H</mi> <mrow> <mo>(</mo> <msub> <mi>I</mi> <msub> <mi>N</mi> <mi>T</mi> </msub> </msub> <mo>&CircleTimes;</mo> <msup> <mi>w</mi> <mi>T</mi> </msup> <mo>)</mo> </mrow> </mrow> </math>
<math> <mrow> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>h</mi> <mn>1</mn> </msub> <mo>&CircleTimes;</mo> <msup> <mi>w</mi> <mi>T</mi> </msup> </mtd> <mtd> <msub> <mi>h</mi> <mn>2</mn> </msub> <mo>&CircleTimes;</mo> <msup> <mi>w</mi> <mi>T</mi> </msup> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <msub> <mi>h</mi> <msub> <mi>N</mi> <mi>T</mi> </msub> </msub> <mo>&CircleTimes;</mo> <msup> <mi>w</mi> <mi>T</mi> </msup> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>14</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein h isiRepresenting the ith column of the channel matrix H. Definition of Nc=NTMrDimension of the composite channel matrix A is NR×Nc. The equivalent channel model can be written as:
y=Au+n (15)
in addition, when the system uses natural mapping as shown in fig. 3, then the symbol s on the kth transmit antennakAnd its corresponding mapping bit xk,j,j=1,…,MrThe relationship between can be represented by:
<math> <mrow> <msub> <mi>s</mi> <mi>k</mi> </msub> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>M</mi> <mi>c</mi> </msub> </munderover> <msub> <mi>w</mi> <mi>i</mi> </msub> <msub> <mi>x</mi> <mrow> <mi>k</mi> <mo>,</mo> <mi>i</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>16</mn> <mo>)</mo> </mrow> </mrow> </math>
skthe vector of (d) is represented as: sk=wTxk (17)
Wherein, <math> <mrow> <msub> <mi>x</mi> <mi>k</mi> </msub> <mo>=</mo> <msup> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>x</mi> <mrow> <mi>k</mi> <mo>,</mo> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>x</mi> <mrow> <mi>k</mi> <mo>,</mo> <mn>2</mn> </mrow> </msub> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <msub> <mi>x</mi> <mrow> <mi>k</mi> <mo>,</mo> <msub> <mi>M</mi> <mi>r</mi> </msub> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> <mo>.</mo> </mrow> </math>
when a square QAM constellation is used, the weighting factor w in equation (17) is calculated as equation (11). Then, similar to gray mapping, the composite channel matrix is defined as:
<math> <mrow> <mi>A</mi> <mo>=</mo> <mi>H</mi> <mrow> <mo>(</mo> <msub> <mi>I</mi> <msub> <mi>N</mi> <mi>T</mi> </msub> </msub> <mo>&CircleTimes;</mo> <msup> <mi>w</mi> <mi>T</mi> </msup> <mo>)</mo> </mrow> </mrow> </math>
<math> <mrow> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>h</mi> <mn>1</mn> </msub> <mo>&CircleTimes;</mo> <msup> <mi>w</mi> <mi>T</mi> </msup> </mtd> <mtd> <msub> <mi>h</mi> <mn>2</mn> </msub> <mo>&CircleTimes;</mo> <msup> <mi>w</mi> <mi>T</mi> </msup> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <msub> <mi>h</mi> <msub> <mi>N</mi> <mi>T</mi> </msub> </msub> <mo>&CircleTimes;</mo> <msup> <mi>w</mi> <mi>T</mi> </msup> </mtd> </mtr> </mtable> </mfenced> </mrow> </math>
the dimension of the composite channel matrix A is NR×Nc. The equivalent channel model can be written as:
y=Ax+n (18)
wherein, <math> <mrow> <mi>x</mi> <mo>=</mo> <msup> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msubsup> <mi>x</mi> <mn>1</mn> <mi>T</mi> </msubsup> </mtd> <mtd> <msubsup> <mi>x</mi> <mn>2</mn> <mi>T</mi> </msubsup> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <msubsup> <mi>x</mi> <msub> <mi>N</mi> <mi>T</mi> </msub> <mi>T</mi> </msubsup> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> </mrow> </math>
when QR decomposition is performed on the composite channel matrix by using MMSE criterion to obtain an upper triangular matrix, and a bit-by-bit hierarchical tree structure is constructed by using the upper triangular matrix, the method may be performed as follows.
First, in order to obtain an equivalent full rank matrix, based on the MMSE criterion, an extended composite channel matrix is defined as:
<math> <mrow> <munder> <mi>A</mi> <mo>&OverBar;</mo> </munder> <mover> <mo>=</mo> <mi>&Delta;</mi> </mover> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mi>A</mi> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&sigma;</mi> <mi>n</mi> </msub> <msub> <mi>I</mi> <msub> <mi>N</mi> <mi>c</mi> </msub> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>19</mn> <mo>)</mo> </mrow> </mrow> </math>
secondly, for the expanded composite channel matrixAQR decomposition is carried out to obtainAQR, hereThe column vectors thereof are orthogonal to each other, and
Figure GSA00000066275700064
is an upper triangular matrix.
It should be noted here that, in the signal detection process in the mimo system, the metric is calculated according to the following formula:
<math> <mrow> <mi>&mu;</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <msubsup> <mi>&sigma;</mi> <mi>n</mi> <mn>2</mn> </msubsup> </mrow> </mfrac> <msup> <mrow> <mo>|</mo> <mo>|</mo> <mi>y</mi> <mo>-</mo> <mi>Hs</mi> <mo>|</mo> <mo>|</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msup> <mi>x</mi> <mi>T</mi> </msup> <msub> <mi>L</mi> <mi>A</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>20</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein L isA(x) Representing a priori information corresponding to the bit vector x.
Finally, to get better detection, we apply matrices before QR decompositionAThe columns of (a) are reordered.
For Gray mapping, | | y-Hs | | sweet calculation in formula (20)2And
Figure GSA00000066275700066
and then the measured values are correlated to obtain the gray mapping,
| | y - Hs | | 2 = | | y - Au | | 2
= y H y - y H Au - u H A H y + u H A H Au
<math> <mrow> <mo>=</mo> <msup> <mi>y</mi> <mi>H</mi> </msup> <mi>y</mi> <mo>-</mo> <msup> <mover> <mi>u</mi> <mo>^</mo> </mover> <mi>H</mi> </msup> <msup> <mi>R</mi> <mi>H</mi> </msup> <mi>Ru</mi> <mo>-</mo> <msup> <mi>u</mi> <mi>H</mi> </msup> <msup> <mi>R</mi> <mi>H</mi> </msup> <mi>R</mi> <mover> <mi>u</mi> <mo>^</mo> </mover> <mo>+</mo> <msup> <mi>u</mi> <mi>H</mi> </msup> <mrow> <mo>(</mo> <msup> <mi>R</mi> <mi>H</mi> </msup> <mi>R</mi> <mo>-</mo> <msubsup> <mi>&sigma;</mi> <mi>n</mi> <mn>2</mn> </msubsup> <msub> <mi>I</mi> <msub> <mi>N</mi> <mi>c</mi> </msub> </msub> <mo>)</mo> </mrow> <mi>u</mi> </mrow> </math>
<math> <mrow> <mo>=</mo> <msup> <mover> <mi>u</mi> <mo>^</mo> </mover> <mi>H</mi> </msup> <msup> <mi>R</mi> <mi>H</mi> </msup> <mi>R</mi> <mover> <mi>u</mi> <mo>^</mo> </mover> <mo>-</mo> <msup> <mover> <mi>u</mi> <mo>^</mo> </mover> <mi>H</mi> </msup> <msup> <mi>R</mi> <mi>H</mi> </msup> <mi>Ru</mi> <mo>-</mo> <msup> <mi>u</mi> <mi>H</mi> </msup> <msup> <mi>R</mi> <mi>H</mi> </msup> <mi>R</mi> <mover> <mi>u</mi> <mo>^</mo> </mover> <mo>+</mo> <msup> <mi>u</mi> <mi>H</mi> </msup> <msup> <mi>R</mi> <mi>H</mi> </msup> <mi>Ru</mi> <mo>-</mo> <msubsup> <mi>&sigma;</mi> <mi>n</mi> <mn>2</mn> </msubsup> <msup> <mi>u</mi> <mi>H</mi> </msup> <mi>u</mi> <mo>-</mo> <msup> <mover> <mi>u</mi> <mo>^</mo> </mover> <mi>H</mi> </msup> <msup> <mi>R</mi> <mi>H</mi> </msup> <mi>R</mi> <mover> <mi>u</mi> <mo>^</mo> </mover> <mo>+</mo> <msup> <mi>y</mi> <mi>H</mi> </msup> <mi>y</mi> </mrow> </math>
<math> <mrow> <mo>=</mo> <msup> <mrow> <mo>|</mo> <mo>|</mo> <mi>R</mi> <mrow> <mo>(</mo> <mi>u</mi> <mo>-</mo> <mover> <mi>u</mi> <mo>^</mo> </mover> <mo>)</mo> </mrow> <mo>|</mo> <mo>|</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <msubsup> <mi>&sigma;</mi> <mi>n</mi> <mn>2</mn> </msubsup> <msup> <mi>u</mi> <mi>T</mi> </msup> <mi>u</mi> <mo>-</mo> <msup> <mover> <mi>u</mi> <mo>^</mo> </mover> <mi>H</mi> </msup> <msup> <mi>R</mi> <mi>H</mi> </msup> <mi>R</mi> <mover> <mi>u</mi> <mo>^</mo> </mover> <mo>+</mo> <msup> <mi>y</mi> <mi>H</mi> </msup> <mi>y</mi> </mrow> </math>
<math> <mrow> <mo>=</mo> <msup> <mrow> <mo>|</mo> <mo>|</mo> <mi>R</mi> <mrow> <mo>(</mo> <mi>u</mi> <mo>-</mo> <mover> <mi>u</mi> <mo>^</mo> </mover> <mo>)</mo> </mrow> <mo>|</mo> <mo>|</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <msubsup> <mi>&sigma;</mi> <mi>n</mi> <mn>2</mn> </msubsup> <msup> <mi>u</mi> <mi>T</mi> </msup> <mi>u</mi> <mo>-</mo> <msup> <mi>y</mi> <mi>H</mi> </msup> <mrow> <mo>(</mo> <mi>A</mi> <msup> <mrow> <mo>(</mo> <msup> <mi>A</mi> <mi>H</mi> </msup> <mi>A</mi> <mo>+</mo> <msubsup> <mi>&sigma;</mi> <mi>n</mi> <mn>2</mn> </msubsup> <msub> <mi>I</mi> <msub> <mi>N</mi> <mi>c</mi> </msub> </msub> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msup> <mi>A</mi> <mi>H</mi> </msup> <mo>-</mo> <msub> <mi>I</mi> <msub> <mi>N</mi> <mi>R</mi> </msub> </msub> <mo>)</mo> </mrow> <mi>y</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>21</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein, <math> <mrow> <mover> <mi>u</mi> <mo>^</mo> </mover> <mo>=</mo> <msup> <mrow> <mo>(</mo> <msup> <mi>A</mi> <mi>H</mi> </msup> <mi>A</mi> <mo>+</mo> <msubsup> <mi>&sigma;</mi> <mi>n</mi> <mn>2</mn> </msubsup> <msub> <mi>I</mi> <msub> <mi>N</mi> <mi>c</mi> </msub> </msub> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msup> <mi>A</mi> <mi>H</mi> </msup> <mi>y</mi> <mo>.</mo> </mrow> </math>
order to <math> <mrow> <mi>C</mi> <mo>=</mo> <msubsup> <mi>&sigma;</mi> <mi>n</mi> <mn>2</mn> </msubsup> <msup> <mi>u</mi> <mi>T</mi> </msup> <mi>u</mi> <mo>+</mo> <msup> <mi>y</mi> <mi>H</mi> </msup> <mrow> <mo>(</mo> <mi>A</mi> <msup> <mrow> <mo>(</mo> <msup> <mi>A</mi> <mi>H</mi> </msup> <mi>A</mi> <mo>+</mo> <msubsup> <mi>&sigma;</mi> <mi>n</mi> <mn>2</mn> </msubsup> <msub> <mi>I</mi> <msub> <mi>N</mi> <mi>c</mi> </msub> </msub> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msup> <mi>A</mi> <mi>H</mi> </msup> <mo>-</mo> <msub> <mi>I</mi> <msub> <mi>N</mi> <mi>R</mi> </msub> </msub> <mo>)</mo> </mrow> <mi>y</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>22</mn> <mo>)</mo> </mrow> </mrow> </math>
Since the symbols in the intermediate vector u are all bit symbols, | u is satisfiedk,i|2At 1, it is clear that the size of C has no relation to the choice of the intermediate variable u at this time, and therefore it can be omitted from the metric calculation. Then the metric can be re-expressed as:
<math> <mrow> <mi>&mu;</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <msubsup> <mi>&sigma;</mi> <mi>n</mi> <mn>2</mn> </msubsup> </mrow> </mfrac> <msup> <mrow> <mo>|</mo> <mo>|</mo> <mi>R</mi> <mrow> <mo>(</mo> <mi>u</mi> <mo>-</mo> <mover> <mi>u</mi> <mo>^</mo> </mover> <mo>)</mo> </mrow> <mo>|</mo> <mo>|</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msup> <mi>x</mi> <mi>T</mi> </msup> <msub> <mi>L</mi> <mi>A</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>n</mi> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </math>
<math> <mrow> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>c</mi> </msub> </munderover> <mrow> <mo>(</mo> <mo>-</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <msubsup> <mi>&sigma;</mi> <mi>n</mi> <mn>2</mn> </msubsup> </mrow> </mfrac> <msup> <mrow> <mo>|</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>i</mi> </munderover> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mi>j</mi> </msub> <mo>-</mo> <msub> <mover> <mi>u</mi> <mo>^</mo> </mover> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>|</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msub> <mi>x</mi> <mi>i</mi> </msub> <msub> <mi>L</mi> <mi>A</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>23</mn> <mo>)</mo> </mrow> </mrow> </math>
further, for natural mapping, | | y-Hs | | survival in equation (20)2And
Figure GSA00000066275700074
and then the measured values are correlated to obtain the natural mapping,
<math> <mrow> <mi>&mu;</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <msubsup> <mi>&sigma;</mi> <mi>n</mi> <mn>2</mn> </msubsup> </mrow> </mfrac> <msup> <mrow> <mo>|</mo> <mo>|</mo> <mi>R</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <mover> <mi>x</mi> <mo>^</mo> </mover> <mo>)</mo> </mrow> <mo>|</mo> <mo>|</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msup> <mi>x</mi> <mi>T</mi> </msup> <msub> <mi>L</mi> <mi>A</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> </mrow> </math>
<math> <mrow> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>c</mi> </msub> </munderover> <mrow> <mo>(</mo> <mo>-</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <msubsup> <mi>&sigma;</mi> <mi>n</mi> <mn>2</mn> </msubsup> </mrow> </mfrac> <msup> <mrow> <mo>|</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>i</mi> </munderover> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>j</mi> </msub> <mo>-</mo> <msub> <mover> <mi>x</mi> <mo>^</mo> </mover> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>|</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msub> <mi>x</mi> <mi>i</mi> </msub> <msub> <mi>L</mi> <mi>A</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>24</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein, <math> <mrow> <mover> <mi>x</mi> <mo>^</mo> </mover> <mo>=</mo> <msup> <mrow> <mo>(</mo> <msup> <mi>A</mi> <mi>H</mi> </msup> <mi>A</mi> <mo>+</mo> <msubsup> <mi>&sigma;</mi> <mi>n</mi> <mn>2</mn> </msubsup> <msub> <mi>I</mi> <msub> <mi>N</mi> <mi>c</mi> </msub> </msub> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msup> <mi>A</mi> <mi>H</mi> </msup> <mi>y</mi> <mo>.</mo> </mrow> </math>
the above formula (23) and formula (24) are mathematical expressions of signal metric values under gray mapping and natural mapping, respectively, that is, a tree structure layered bit by bit is searched layer by layer based on a breadth-first M algorithm, a signal candidate set can be obtained, and a mathematical model of a metric value process of each candidate signal in the signal candidate set is calculated.
The method comprises the following steps of searching a tree structure which is layered bit by bit on the basis of a breadth-first M algorithm layer by layer to obtain a signal candidate set, and calculating the metric value of the signal candidate set comprises the following steps: calculating the metric values of 2M nodes on 2M branch paths expanded by M retention nodes of the previous layer in the current layer of the bit-by-bit hierarchical tree structure; searching M nodes with the maximum metric value from the 2M nodes, and taking the M nodes with the maximum metric value as retention nodes of the current layer; and sequentially searching layer by layer, selecting M retention nodes from the 2M nodes of the last layer as a signal candidate set, and calculating the metric value of the signal candidate set.
For example, as shown in fig. 4, at the k-th level of the bit-by-bit hierarchical tree structure, metric values of 8 nodes on 8 branch paths expanded by 4 retention nodes of the k-1 level are calculated, metric values of the 8 nodes are compared, and 4 nodes with the largest metric values are selected as retention nodes of the k-th level. And searching layer by layer in sequence, when the last layer is calculated, selecting 4 retention nodes from 8 nodes of the last layer, taking the 4 retention nodes as a signal candidate set, and calculating the metric values of the 4 retention nodes in the signal candidate set.
In gray mapping, the updated calculation of the metric value can be expressed by the following mathematical expression:
<math> <mrow> <msub> <mi>&mu;</mi> <mn>1</mn> </msub> <mo>=</mo> <mo>-</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <msubsup> <mi>&sigma;</mi> <mi>n</mi> <mn>2</mn> </msubsup> </mrow> </mfrac> <msup> <mrow> <mo>|</mo> <msub> <mi>r</mi> <mn>1,1</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mover> <mi>u</mi> <mo>^</mo> </mover> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>|</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msub> <mi>x</mi> <mn>1</mn> </msub> <msub> <mi>L</mi> <mi>A</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </math>
<math> <mrow> <msub> <mi>&mu;</mi> <mi>i</mi> </msub> <mo>=</mo> <msub> <mi>&mu;</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <msubsup> <mi>&sigma;</mi> <mi>n</mi> <mn>2</mn> </msubsup> </mrow> </mfrac> <msup> <mrow> <mo>|</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>i</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mover> <mi>u</mi> <mo>^</mo> </mover> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>i</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mi>j</mi> </msub> <mo>-</mo> <msub> <mover> <mi>u</mi> <mo>^</mo> </mover> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>|</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msub> <mi>x</mi> <mi>i</mi> </msub> <msub> <mi>L</mi> <mi>A</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>,</mo> <mn>2</mn> <mo>&le;</mo> <mi>i</mi> <mo>&le;</mo> <msub> <mi>N</mi> <mi>c</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>25</mn> <mo>)</mo> </mrow> </mrow> </math>
<math> <mrow> <mi>&mu;</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>&mu;</mi> <msub> <mi>N</mi> <mi>c</mi> </msub> </msub> </mrow> </math>
the process of metric hierarchical computation is seen from equation (25), and after the computation of the metric of each layer is completed, the metric is sorted, and M nodes with the maximum metric are reserved as the reserved nodes for the next expansion.
In the natural mapping, the updated calculation of the metric value can be expressed by the following mathematical expression:
<math> <mrow> <msub> <mi>&mu;</mi> <mn>1</mn> </msub> <mo>=</mo> <mo>-</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <msubsup> <mi>&sigma;</mi> <mi>n</mi> <mn>2</mn> </msubsup> </mrow> </mfrac> <msup> <mrow> <mo>|</mo> <msub> <mi>r</mi> <mn>1,1</mn> </msub> <mrow> <mo>(</mo> <msub> <mrow> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>-</mo> <mover> <mi>x</mi> <mo>^</mo> </mover> </mrow> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>|</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msub> <mi>x</mi> <mn>1</mn> </msub> <msub> <mi>L</mi> <mi>A</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </math>
<math> <mrow> <msub> <mi>&mu;</mi> <mi>i</mi> </msub> <mo>=</mo> <msub> <mi>&mu;</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <msubsup> <mi>&sigma;</mi> <mi>n</mi> <mn>2</mn> </msubsup> </mrow> </mfrac> <msup> <mrow> <mo>|</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>i</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mover> <mi>x</mi> <mo>^</mo> </mover> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>i</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>j</mi> </msub> <mo>-</mo> <msub> <mover> <mi>x</mi> <mo>^</mo> </mover> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>|</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msub> <mi>x</mi> <mi>i</mi> </msub> <msub> <mi>L</mi> <mi>A</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>,</mo> <mn>2</mn> <mo>&le;</mo> <mi>i</mi> <mo>&le;</mo> <msub> <mi>N</mi> <mi>c</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>26</mn> <mo>)</mo> </mrow> </mrow> </math>
<math> <mrow> <mi>&mu;</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>&mu;</mi> <msub> <mi>N</mi> <mi>c</mi> </msub> </msub> </mrow> </math>
in a preferred embodiment, performing a layer-by-layer search on a tree structure layered bit by bit based on a breadth-first M algorithm to obtain a signal candidate set, and calculating a metric value of the signal candidate set includes: in the tree structure of the bit-by-bit layering, acquiring prior information of corresponding bits of all layers from a decoder; the layer corresponding to the node with the prior information larger than the first threshold value is not searched and expanded, the preserved path of the preserved node on the upper layer is expanded by using the symbol of the prior information of the current layer, the node with the prior information not larger than the first threshold value is expanded, and M nodes with the maximum metric value are searched as the preserved nodes of the current layer according to the magnitude of the metric value; and sequentially searching layer by layer, selecting M retention nodes from the last layer as a signal candidate set, and calculating the metric value of the signal candidate set.
For example, in a tree structure of bit-by-bit layering, prior information of corresponding bits of all layers is obtained from a decoder, the prior information of the corresponding bits of all layers is compared with a first threshold, for example, if the prior information of the corresponding bits of the K-th layer is greater than the first threshold, the surviving paths of 4 surviving nodes of the K-1-th layer are extended by using symbols of the prior information of the K-th layer, and if the prior information of the corresponding bits of the K + 1-th layer is not greater than the first threshold, 4 surviving nodes of the K-1-th layer are expanded, metric values of 8 nodes on 8 branch paths expanded by the 4 nodes are calculated, the metric values of the 8 nodes are compared, and the largest 4 nodes of the metric values are selected as the surviving nodes of the K + 1-th layer. And sequentially searching layer by layer, selecting 4 retention nodes from the last layer when the last layer is searched, taking the 4 retention nodes as a signal candidate set, and calculating the metric value of the signal candidate set.
Another preferred embodiment is that, performing layer-by-layer search on a tree structure layered bit by bit based on a breadth-first M algorithm to obtain a signal candidate set, and calculating a metric value of the signal candidate set includes: in the tree structure of the bit-by-bit layering, acquiring prior information of corresponding bits of all layers from a decoder; taking the prior information weighted value of the node of which the prior information is greater than the second threshold value as a metric value, searching M nodes with the maximum metric value from 2M nodes spread by M retention nodes on the upper layer, and taking the M nodes with the maximum metric value as the retention nodes on the current layer; and sequentially searching layer by layer, selecting M retention nodes from the last layer as a signal candidate set, and calculating the metric value of the signal candidate set.
For example, in a tree structure with bit-by-bit hierarchy, the prior information of the corresponding bits of all layers is obtained from the decoder, the prior information of the corresponding bits of all layers is compared with the second threshold respectively, and the weighted value of the prior information of the node with the prior information greater than the second threshold is used as the metric value, e.g., for gray mapping, the weighted value of the prior information in the formula (19)
Figure GSA00000066275700091
I.e. a priori information weighted value, for the natural mapping, in equation (20)
Figure GSA00000066275700092
I.e. the prior information weighted value. And if the prior information of the corresponding bit of the k layer is larger than the second threshold, selecting the prior information weighted value of 8 nodes expanded by the 4 retention nodes in the k-1 layer as the metric value, and selecting the 4 nodes with the maximum metric value as the retention nodes of the k layer. . Searching layer by layer in sequence, when the last layer is searched, selecting 4 retention nodes from the last layer, and searching the 4 retention nodesThe retention node serves as a signal candidate set, and a metric value of the signal candidate set is calculated.
Of course, a combination of the two preferred modes can also be achieved. For example, after obtaining the prior information of the corresponding bits of all layers from the decoder in the bit-by-bit hierarchical tree structure, the prior information of the corresponding bits of all layers is compared with the first threshold value, the layer corresponding to the node of which the prior information is greater than the first threshold value is not searched and expanded, for example, when the prior information of the corresponding bit of the k-th layer is greater than the first threshold value, the retention paths of 4 retention nodes of the k-1-th layer are extended by using the symbol of the prior information of the k-th layer, and the prior information weighted value of the node of which the prior information is greater than the second threshold value is taken as the metric value from the remaining layers. And if the prior information of the corresponding bit of the k layer is larger than the second threshold, selecting the prior information weighted value of 8 nodes expanded by the 4 retention nodes in the k-1 layer as the metric value, and selecting the 4 nodes with the maximum metric value as the retention nodes of the k layer. And for the nodes of which the prior information is not greater than the second threshold, searching M nodes with the maximum metric value from the expanded 2M nodes according to the magnitude of the metric value, and taking the M nodes with the maximum metric value as retention nodes of the current layer. And sequentially searching layer by layer, selecting 4 retention nodes from the last layer when the last layer is searched, taking the 4 retention nodes as a signal candidate set, and calculating the metric value of the signal candidate set.
It should be noted that when gray mapping is used, it can be obtained from the definition of the metric if
Figure GSA00000066275700101
Then we need only compute the metric related to the a priori information. The following relationships between the metrics are:
<math> <mrow> <mfrac> <mn>1</mn> <msubsup> <mi>&sigma;</mi> <mi>n</mi> <mn>2</mn> </msubsup> </mfrac> <msup> <mrow> <mo>|</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>i</mi> </munderover> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mi>j</mi> </msub> <mo>-</mo> <msub> <mover> <mi>u</mi> <mo>^</mo> </mover> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>|</mo> </mrow> <mn>2</mn> </msup> <mo>&le;</mo> <mfrac> <mn>1</mn> <msubsup> <mi>&sigma;</mi> <mi>n</mi> <mn>2</mn> </msubsup> </mfrac> <msup> <mrow> <mo>|</mo> <mo>|</mo> <mi>R</mi> <mrow> <mo>(</mo> <mi>u</mi> <mo>-</mo> <mover> <mi>u</mi> <mo>^</mo> </mover> <mo>)</mo> </mrow> <mo>|</mo> <mo>|</mo> </mrow> <mn>2</mn> </msup> <mo>=</mo> <mfrac> <mn>1</mn> <msubsup> <mi>&sigma;</mi> <mi>n</mi> <mn>2</mn> </msubsup> </mfrac> <mrow> <mo>(</mo> <msup> <mrow> <mo>|</mo> <mo>|</mo> <mi>y</mi> <mo>-</mo> <mi>Hs</mi> <mo>|</mo> <mo>|</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mi>C</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>27</mn> <mo>)</mo> </mrow> </mrow> </math>
||y-Hs||2the expectation of + C may be approximated as σn 2(NR+Nc). When | LA(xi)|>>T1,T1≥NR+NcIn which T is1For the second threshold, only the metric associated with the a priori information may be calculated, and the metric associated with the i-th layer signal may be written as:
<math> <mrow> <msub> <mi>&mu;</mi> <mi>I</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msub> <mi>x</mi> <mi>i</mi> </msub> <msub> <mi>L</mi> <mi>A</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>28</mn> <mo>)</mo> </mrow> </mrow> </math>
as the signal-to-noise ratio increases, the absolute value of the a priori information available from the channel decoder increases. The a priori information about these bits can be considered reliable when the absolute value of the a priori information obtained from the decoder exceeds a set threshold (first threshold). Then we do not compute their metrics for these bits any more, and the depth of the search tree can be reduced significantly. The computational complexity of the algorithm can therefore be reduced as the signal-to-noise ratio increases.
The principles described above are equally applicable to natural mapping.
It should be further noted that the value range of the first threshold is 3-5 times of the second threshold.
It can be seen from the above embodiments that the bit-level tree search detection method in the present application can reduce the branch paths that need to be searched for each time, and only 2M node metrics need to be calculated after each layer of nodes are expanded correspondingly. Compared with the existing tree searching method based on the symbol, in each stage of the tree searching, each retention path needs to be expanded into a range of constellation set according to the size of the constellation setThe branch paths need to be calculated after the corresponding nodes at each layer are expanded
Figure GSA00000066275700105
Measurement of child nodes, and thus complexity and modulation order M of the entire tree search algorithmrIn an exponential relation, the implementation scheme of the method reduces the complexity of the whole algorithm, and therefore the implementation difficulty of the signal detection process is reduced.
EXAMPLE III
Corresponding to the signal detection method in the mimo system, the embodiment of the present application further provides a signal detection apparatus in the mimo system. Please refer to fig. 5, which is a block diagram of an embodiment of a signal detection apparatus in a mimo system according to the present application, the apparatus includes a conversion unit 501, a decomposition unit 502, a search unit 503, and a detection unit 504. The internal structure and connection relationship of the device will be further described below in conjunction with the working principle of the device.
A converting unit 501, configured to convert a constellation point of a quadrature amplitude modulation QAM transmission signal into a bit vector weighted sum form to obtain a bit-level-represented QAM transmission signal, and convert a channel matrix into a composite channel matrix according to the bit-level-represented QAM transmission signal;
a decomposition unit 502, configured to perform QR decomposition on the composite channel matrix to obtain an upper triangular matrix, and construct a bit-by-bit hierarchical tree structure using the upper triangular matrix;
the searching unit 503 is configured to perform layer-by-layer search on the bit-by-bit hierarchical tree structure based on a breadth-first M algorithm to obtain a signal candidate set, and calculate a metric value of each candidate signal in the signal candidate set;
a detecting unit 504, configured to calculate a posteriori information of the transmitted signal using the metric value.
Please refer to fig. 6, which is a block diagram illustrating a signal detection apparatus in a mimo system according to another embodiment of the present application. Wherein, the search unit 503 further includes: a calculation sub-unit 5031, a first search sub-unit 5032 and a selection sub-unit 5033,
a computation subunit 5031, configured to compute, in a current layer of the bit-wise hierarchical tree structure, metric values of 2M nodes on 2M branch paths expanded by M retention nodes of a previous layer;
a first searching subunit 5032, configured to search M nodes with the largest metric value from the 2M nodes, and use the M nodes with the largest metric value as retention nodes of a current layer;
a selecting subunit 5033, configured to sequentially search layer by layer, select M surviving nodes from the 2M nodes in the last layer as a signal candidate set, and calculate a metric value of the signal candidate set.
Preferably, the calculating sub-unit 5031 may be replaced by an obtaining sub-unit, and the first searching sub-unit 5032 may be replaced by a second searching sub-unit, and then the searching unit 503 includes:
an obtaining subunit, configured to obtain, from a decoder, prior information of bits corresponding to all layers in the bit-by-bit hierarchical tree structure;
the second searching subunit is used for not searching and expanding the layer corresponding to the node with the prior information being greater than the first threshold, extending the reserved path of the reserved node in the previous layer by using the symbol of the prior information of the current layer, expanding the node with the prior information being not greater than the first threshold, and searching M nodes with the maximum metric value as the reserved nodes of the current layer according to the magnitude of the metric value;
and the selecting subunit is used for sequentially searching layer by layer, selecting M retention nodes from the last layer as a signal candidate set, and calculating the metric value of the signal candidate set.
Preferably, the calculating sub-unit 5031 may be replaced by an obtaining sub-unit, and the first searching sub-unit 5032 may be replaced by a third searching sub-unit, and then the searching unit 503 includes:
an obtaining subunit, configured to obtain, from a decoder, prior information of bits corresponding to all layers in the bit-by-bit hierarchical tree structure;
a third searching subunit, configured to select, as metric values, M nodes with the largest metric values from the 2M nodes expanded by the M retention nodes on the upper layer, and use the M nodes with the largest metric values as retention nodes on the current layer, where the prior information weighted value of the node whose prior information is greater than the second threshold is used as the metric value;
and the selecting subunit is used for sequentially searching layer by layer, selecting M retention nodes from the last layer as a signal candidate set, and calculating the metric value of the signal candidate set.
The mapping mode of the constellation points of the quadrature amplitude modulation QAM transmission signals comprises Gray mapping or natural mapping.
It can be seen from the above embodiments that the bit-level tree search detection method in the present application can reduce the branch paths that need to be searched for each time, and only 2M node metrics need to be calculated after each layer of nodes are expanded correspondingly. Compared with the existing tree searching method based on the symbol, in each stage of the tree searching, each retention path needs to be expanded into a range of constellation set according to the size of the constellation set
Figure GSA00000066275700121
The branch paths need to be calculated after the corresponding nodes at each layer are expanded
Figure GSA00000066275700122
Measurement of child nodes, and thus complexity and modulation order M of the entire tree search algorithmrIn an exponential relation, the implementation scheme of the method reduces the complexity of the whole algorithm, and therefore the implementation difficulty of the signal detection process is reduced.
It should be noted that, as will be understood by those skilled in the art, all or part of the processes in the methods of the above embodiments may be implemented by a computer program, which may be stored in a computer-readable storage medium, and when executed, may include the processes of the embodiments of the methods described above. The storage medium may be a magnetic disk, an optical disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), or the like.
The foregoing detailed description is directed to a method and apparatus for detecting signals in a mimo system, and the principles and embodiments of the present application are explained with reference to specific embodiments, which are merely used to help understand the method and its core ideas of the present application; meanwhile, for a person skilled in the art, according to the idea of the present application, there may be variations in the specific embodiments and the application scope, and in summary, the content of the present specification should not be construed as a limitation to the present application.

Claims (10)

1. A method for signal detection in a multiple-input multiple-output system, comprising:
converting constellation points of Quadrature Amplitude Modulation (QAM) transmitting signals into a bit vector weighted sum form to obtain QAM transmitting signals represented by a bit level, and converting a channel matrix into a composite channel matrix according to the QAM transmitting signals represented by the bit level;
performing QR decomposition on the composite channel matrix to obtain an upper triangular matrix, and constructing a bit-by-bit hierarchical tree structure by using the upper triangular matrix;
searching the tree structure which is layered bit by bit on the basis of a breadth-first M algorithm layer by layer to obtain a signal candidate set, and calculating the metric value of each candidate signal in the signal candidate set;
and calculating posterior information of the transmitting signal by using the metric value.
2. The method of claim 1, wherein the bit-by-bit hierarchical tree structure is searched layer by layer based on a breadth-first M algorithm to obtain a signal candidate set, and calculating the metric value of the signal candidate set comprises:
calculating the metric values of 2M nodes on 2M branch paths expanded by M retention nodes of the previous layer in the current layer of the bit-by-bit hierarchical tree structure, wherein M is the number of constellation points;
searching M nodes with the maximum metric value from the 2M nodes, and taking the M nodes with the maximum metric value as retention nodes of the current layer;
and sequentially searching layer by layer, selecting M retention nodes from the 2M nodes of the last layer as a signal candidate set, and calculating the metric value of the signal candidate set.
3. The method of claim 1, wherein the bit-by-bit hierarchical tree structure is searched layer by layer based on a breadth-first M algorithm to obtain a signal candidate set, and calculating the metric value of the signal candidate set comprises:
in the tree structure of the bit-by-bit layering, acquiring prior information of corresponding bits of all layers from a decoder;
the method comprises the steps that a layer corresponding to a node with prior information larger than a first threshold value is not searched and expanded, a reserved path of a reserved node in the previous layer is expanded by using a symbol of prior information of the current layer, a node with prior information not larger than the first threshold value is expanded, M nodes with the maximum metric value are searched as reserved nodes of the current layer according to the magnitude of the metric value, and M is the number of constellation points;
and sequentially searching layer by layer, selecting M retention nodes from the last layer as a signal candidate set, and calculating the metric value of the signal candidate set.
4. The method of claim 1, wherein the bit-by-bit hierarchical tree structure is searched layer by layer based on a breadth-first M algorithm to obtain a signal candidate set, and calculating the metric value of the signal candidate set comprises:
in the tree structure of the bit-by-bit layering, acquiring prior information of corresponding bits of all layers from a decoder;
using the prior information weighted value of the bit node with the prior information larger than the second threshold value as a metric value, searching M nodes with the maximum metric value from 2M nodes spread by M retention nodes on the upper layer, using the M nodes with the maximum metric value as the retention nodes on the current layer, wherein M is the number of constellation points;
and sequentially searching layer by layer, selecting M retention nodes from the last layer as a signal candidate set, and calculating the metric value of the signal candidate set.
5. The method according to any of claims 1-4, wherein the mapping manner of the constellation points of the QAM transmission signal comprises Gray mapping or natural mapping.
6. A signal detection apparatus in a multiple-input multiple-output system, comprising:
the conversion unit is used for converting the constellation points of the Quadrature Amplitude Modulation (QAM) sending signals into a bit vector weighted sum form to obtain QAM sending signals represented by a bit level, and obtaining a composite channel matrix according to the QAM sending signals represented by the bit level;
the decomposition unit is used for carrying out QR decomposition on the composite channel matrix to obtain an upper triangular matrix, and a bit-by-bit hierarchical tree structure is constructed by utilizing the upper triangular matrix;
the searching unit is used for searching the bit-by-bit hierarchical tree structure layer by layer based on a breadth-first M algorithm to obtain a signal candidate set and calculating the metric value of each candidate signal in the signal candidate set;
and the detection unit calculates the posterior information of the transmitting signal by using the metric value.
7. The apparatus of claim 6, wherein the search unit comprises:
a calculating subunit, configured to calculate, in a current layer of the bit-by-bit hierarchical tree structure, metric values of 2M nodes on 2M branch paths expanded by M surviving nodes of a previous layer, where M is a number of constellation points;
the first searching subunit is configured to search M nodes with the largest metric value from the 2M nodes, and use the M nodes with the largest metric value as retention nodes of a current layer;
and the selecting subunit is used for sequentially searching layer by layer, selecting M retention nodes from the 2M nodes of the last layer as a signal candidate set, and calculating the metric value of the signal candidate set.
8. The apparatus of claim 6, wherein the search unit comprises:
an obtaining subunit, configured to obtain, from a decoder, prior information of bits corresponding to all layers in the bit-by-bit hierarchical tree structure;
the second searching subunit is used for not performing expansion searching on a layer corresponding to a node of which the prior information is greater than the first threshold, directly extending a retention path of a retention node of a previous layer by using a symbol of the prior information of the current layer, expanding a node of which the prior information is not greater than the first threshold, and searching M nodes with the maximum metric value as retention nodes of the current layer according to the magnitude of the metric value, wherein M is the number of constellation points;
and the selecting subunit is used for sequentially searching layer by layer, selecting M retention nodes from the last layer as a signal candidate set, and calculating the metric value of the signal candidate set.
9. The apparatus of claim 6, wherein the search unit comprises:
an obtaining subunit, configured to obtain, from a decoder, prior information of bits corresponding to all layers in the bit-by-bit hierarchical tree structure;
a third searching subunit, configured to select, from 2M nodes spread by M surviving nodes on an upper layer, M nodes with a largest metric value as surviving nodes on a current layer, where M is a number of constellation points, by using a prior information weighted value of a node whose prior information is greater than a second threshold as the metric value;
and the selecting subunit is used for sequentially searching layer by layer, selecting M retention nodes from the last layer as a signal candidate set, and calculating the metric value of the signal candidate set.
10. The apparatus according to any of claims 6-9, wherein the mapping manner of the constellation points of the QAM transmission signal comprises gray mapping or natural mapping.
CN2010101354912A 2010-03-29 2010-03-29 Signal detection method and device in multiple-input multiple-output system Expired - Fee Related CN101834827B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010101354912A CN101834827B (en) 2010-03-29 2010-03-29 Signal detection method and device in multiple-input multiple-output system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010101354912A CN101834827B (en) 2010-03-29 2010-03-29 Signal detection method and device in multiple-input multiple-output system

Publications (2)

Publication Number Publication Date
CN101834827A CN101834827A (en) 2010-09-15
CN101834827B true CN101834827B (en) 2012-07-18

Family

ID=42718759

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010101354912A Expired - Fee Related CN101834827B (en) 2010-03-29 2010-03-29 Signal detection method and device in multiple-input multiple-output system

Country Status (1)

Country Link
CN (1) CN101834827B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102487309B (en) * 2010-12-01 2014-09-17 北京大学 Signal detecting method and device under MIMO (Multiple Input Multiple Output) system
CN102594467B (en) * 2011-08-15 2014-07-02 上海交通大学 Receiver detection method for wireless multiple input multiple output system
CN102970085B (en) * 2012-11-19 2015-01-14 北京航空航天大学 Signal detecting method
CN108334944B (en) * 2016-12-23 2020-04-17 中科寒武纪科技股份有限公司 Artificial neural network operation device and method
CN114731323B (en) * 2020-11-04 2023-09-12 华为技术有限公司 Detection method and device for Multiple Input Multiple Output (MIMO) system
CN114268411B (en) * 2021-11-05 2024-07-23 网络通信与安全紫金山实验室 Hard output MIMO detection method and system, electronic device and storage medium

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080059869A1 (en) * 2006-09-01 2008-03-06 The Regents Of The University Of California Low cost, high performance error detection and correction
CN101541023A (en) * 2008-03-18 2009-09-23 大唐移动通信设备有限公司 Joint iterative detection decoding method and device thereof

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080059869A1 (en) * 2006-09-01 2008-03-06 The Regents Of The University Of California Low cost, high performance error detection and correction
CN101541023A (en) * 2008-03-18 2009-09-23 大唐移动通信设备有限公司 Joint iterative detection decoding method and device thereof

Also Published As

Publication number Publication date
CN101834827A (en) 2010-09-15

Similar Documents

Publication Publication Date Title
CN101834827B (en) Signal detection method and device in multiple-input multiple-output system
JP6405155B2 (en) Signal processing apparatus, signal processing method, and program
CN101997652B (en) Acceptance detection method and device based on LDPC-MIMO (low density parity check-multiple input multiple output) communication system
CN105721106A (en) Multiuser detection method based on serial strategy for SCMA (Sparse Code Multiple Access) uplink communication system
CN107743056B (en) SCMA (sparse code multiple access) multi-user detection method based on compressed sensing assistance
CN106357312B (en) Lattice about subtract auxiliary breadth First tree search MIMO detection method
US8831128B2 (en) MIMO communication system signal detection method
CN105071843A (en) Large-scale MIMO system low-complexity polynomial expansion matrix inversion method and application thereof
Wang et al. Online LSTM-based channel estimation for HF MIMO SC-FDE system
CN105634568A (en) LLR calculation method based on large-scale MIMO system signal detection
CN101227254A (en) Method for detecting V-BLAST in MIMO system
CN101330361B (en) Method and apparatus for detecting signal of multi-code multi-transmission multi-receiving system
CN107276703B (en) Orthogonal space modulation system detection method adopting compressed sensing technology
CN101958875B (en) Detecting method of high order modulated MIMO system in mobile environment
CN111541472B (en) Low-complexity machine learning assisted robust precoding method and device
CN101640649B (en) Method and device for determining channel prediction factor in channel prediction and channel predictor
CN106357318A (en) Large-scale MIMO (Multiple Input Multiple Output) iterative detection method with adjustable convergence rate
CN110504995A (en) Soft output MIMO detection method based on lattice reduction and K-Best
CN113364535A (en) Method, system, device and storage medium for mathematical form multiple-input multiple-output detection
CN111769975A (en) MIMO system signal detection method and system
CN101640583B (en) Method for transmitting preprocessing
CN110868244A (en) Low-complexity communication signal detection method based on channel puncture
CN114050853B (en) Multi-user MIMO transmission method based on joint non-orthogonal codebook and pre-coding design
CN115001544B (en) Signal detection method based on improved Richardson algorithm
CN114978254B (en) Machine learning auxiliary low-complexity robust precoding algorithm based on OFDM system

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120718

Termination date: 20130329