CN102412931A - Serial interference cancellation and detection method and system for multi-antenna bit interleaved coding system - Google Patents

Serial interference cancellation and detection method and system for multi-antenna bit interleaved coding system Download PDF

Info

Publication number
CN102412931A
CN102412931A CN2011103651388A CN201110365138A CN102412931A CN 102412931 A CN102412931 A CN 102412931A CN 2011103651388 A CN2011103651388 A CN 2011103651388A CN 201110365138 A CN201110365138 A CN 201110365138A CN 102412931 A CN102412931 A CN 102412931A
Authority
CN
China
Prior art keywords
data
module
mmse
noise ratio
detection signal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN2011103651388A
Other languages
Chinese (zh)
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.)
Xinyuan Technology (xi'an) Co Ltd
Original Assignee
Xinyuan Technology (xi'an) 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 Xinyuan Technology (xi'an) Co Ltd filed Critical Xinyuan Technology (xi'an) Co Ltd
Priority to CN2011103651388A priority Critical patent/CN102412931A/en
Publication of CN102412931A publication Critical patent/CN102412931A/en
Pending legal-status Critical Current

Links

Images

Abstract

The invention relates to a serial interference cancellation and detection method and system for a multi-antenna bit interleaved coding system. The system disclosed by the invention comprises an MMSE (Minimum Mean Square Error) equalization filter module, a post-detection signal-to-noise ratio computing module, a sequencing module, a data channel section module, a soft demodulation module, a data revaluation module, an interference cancellation module, a residual interference estimation module, a similar MMSE equalization filter module, a similar post-detection signal-to-noise ratio computing module and a judgment module. According to the method and the system, the technical problem that the residual interference of 'eliminated' interference is considered to be zero existing in the conventional serial interference cancellation method is solved. In the invention, a correct flow device is provided for an MMSE serial interference cancellation method, a residual interference estimation method is further provided, and a correct form for estimating the MMSE of data in residual layers and a corresponding computation form of a post-detection signal-to-noise ratio used in the presence of 'residual interference' are provided.

Description

The counteracting serial interference detection method and the system of many antennas bit interweaving encoding system
Technical field
The present invention relates to wireless communication field, relate in particular in the MIMO antenna system (MIMO) Interference Cancellation detection technique multi-layer data.
Background technology
MIMO radio communication system can utilize multiplexing to the space, is sending a multi-layer data with a frequency domain resource simultaneously, thereby is increasing the power system capacity of wireless communication system.Yet the detection of every layer data has great influence to the inhibition ability from the interference of other transmission antennas transmit data to systematic function in its testing process.Least mean-square error (MMSE) method is a kind of receiving terminal detection method commonly used, and the present invention has added soft-decision counteracting serial interference detection method and system on the basis of least mean-square error method, progressively eliminate the interference of each layer data to other layer data.
What the present invention was different with the traditional serial interference cancellation method is, the traditional serial interference cancellation method to the interference of " elimination " often think that its residual interference is zero.Yet " elimination " residual interference of disturbing is very important to playing a part in the equalization filtering of follow-up other layer data and the calculating of detection signal-to-noise ratio thereafter.We know that soft bit LLR is used in decoding in the Bit Interleave coded system, and back detection signal-to-noise ratio provides the weighting to the LLR confidence level in the calculating of LLR, if the calculating of back detection signal-to-noise ratio is incorrect, and the very big loss that can bring systematic function.The present invention not only provides the correct flow process and the system of least mean-square error counteracting serial interference method; The method of estimation of residual interference more is provided; Also provide under the situation that has ' residual interference '; The correct form that rest layers data least mean-square error is estimated, and the form of calculation of corresponding back detection signal-to-noise ratio.
Summary of the invention
For the interference of " eliminations " to that solves that existing counteracting serial interference method exists often thinks that its residual interference is zero technical problem; The present invention provides counteracting serial interference detection method and the detection system based on many antennas bit interweaving encoding system; Be applicable to that transmitting terminal and receiving terminal dispose respectively under the situation of many antennas in the radio communication; The correct flow path device of least mean-square error counteracting serial interference method not only is provided, the method for estimation of residual interference more is provided, also provide under the situation that has " residual interference "; The correct form that rest layers data least mean-square error is estimated, and the form of calculation of corresponding back detection signal-to-noise ratio.
Technical solution of the present invention:
The present invention is based on many antennas bit interweaving encoding system the counteracting serial interference detection method based on system configuration comprise: MMSE equalization filter module, back detection signal-to-noise ratio computing module, order module; Choose the data channel module, soft demodulation module, data revaluation module; Interference cancellation module; The residual interference estimation module, type MMSE equalization filter module, type back detection signal-to-noise ratio computing module.
Based on the counteracting serial interference detection method of many antennas bit interweaving encoding system, its special character is: may further comprise the steps:
1] in the MIMO model, calculates MMSE and detect data
Figure BDA0000109414020000022
and back detection signal-to-noise ratio w according to MMSE equalization filter
Figure BDA0000109414020000021
Wherein: MIMO model: y=Hx+n, transmitting terminal configuration Nt transmit antennas, receiving terminal configuration Nr root reception antenna;
G MMSEBe the MMSE equalization filter, its expression formula is: G MMSE=H H(HH H+ σ n 2I Nr) -1
H is a channel matrix, the element h of H IjRepresent the channel coefficients of j transmit antennas, x=[x to i root reception antenna k] T, k=1,2 ... Nt, expression transmission data, y=[y k] T, k=1,2 ... Nr, expression receives data, n=[n k] T, k=1,2 ... Nr, represent the noise of every reception antenna,
Calculate back detection signal-to-noise ratio w:w=1/ (1-diag (G MMSEH))-1;
2] sort from big to small according to back detection signal-to-noise ratio w, obtain ordered series of numbers m=[m 1, m 2... m Nt], satisfy w m 1 > w m 2 > . . . > w m Nt ;
3] choose maximum back detection signal-to-noise ratio data channel m iAs the current data channel of choosing, make i=1, export corresponding back detection signal-to-noise ratio of the current data channel of choosing and corresponding detection data;
4] to step 3] in the current data channel of choosing carry out soft demodulation process, according to
Figure BDA0000109414020000032
Calculate the LLR of this data channel,
Figure BDA0000109414020000033
Represent the current m that chooses iThe LLR value that data channel k bit is corresponding;
5] according to step 4] in the LLR that obtains the detection data of the current data channel of choosing are carried out revaluation, obtain revaluation and detect data
Figure BDA0000109414020000034
P ( θ | x ^ m i ) = Π m = 1 M 1 2 ( 1 + ( 2 c k , m - 1 ) tanh ( L k m i 2 ) )
Wherein θ representes the constellation point value of modulator approach that this data channel adopts, c K, mThe corresponding m bit value (0/1) of expression constellation point value θ;
6] from receive data, eliminate the current interference of choosing data channel, obtain the reception data model
Figure BDA0000109414020000037
after the interference eliminated
7] the current residual interference energy value of choosing data channel after basis
Figure BDA0000109414020000038
the calculating interference eliminated;
8] utilize step 6] in reception data model after the interference eliminated that obtains, use a type MMSE equalization filter, calculate m I+1Detection data on the circuit-switched data passage:
Figure BDA0000109414020000039
And output
Figure BDA00001094140200000310
Type MMSE equalization filter wherein G MMSE ′ = h m i + 1 H ( H M H M H + q m i h m i h m i H + σ n 2 I Nr ) - 1 ,
H M=[h k], k=1,2 ... Nt, k ≠ m 1, m 2... m i, eliminate the data channel that had detected, the channel matrix of remaining data passage among the expression H;
9] basis
Figure BDA00001094140200000312
Calculate m I+1The back detection signal-to-noise ratio of circuit-switched data passage;
10] judge, if i+1 is less than Nt, i++, execution in step 4]; Otherwise if i+1 equals Nt, each data channel all detects, then finishes.
Step 8] middle G ' MMSEAdopt following method to calculate:
Order A i = ( H M H M H + q m i h m i h m i H + σ n 2 I Nr ) , Class MMSE equalization filter is written as following form:
G MMSE ′ = h m i + 1 H A i - 1
A wherein i -1Obtain through following alternative manner:
A i - 1 = A i - 1 - 1 - A i - 1 - 1 h m i ( A i - 1 - 1 h m i ) H 1 q m i - 1 + h m i H ( A i - 1 - 1 h m i )
A 0 - 1 = ( HH H + σ n 2 I ) - 1 .
Comprise the MMSE equalization filter module that connects successively, back detection signal-to-noise ratio computing module, order module; Choose the data channel module, soft demodulation module, data revaluation module; Interference cancellation module, residual interference estimation module, type MMSE equalization filter module; Detection signal-to-noise ratio computing module and judge module after type
MMSE equalization filter module according to the model used in MIMO MMSE equalization filter
Figure BDA0000109414020000045
calculate MMSE detector data
Figure BDA0000109414020000046
where: MIMO models:
Y=Hx+n, transmitting terminal configuration Nt transmit antennas, receiving terminal configuration Nr root reception antenna;
G MMSEBe the MMSE equalization filter, its expression formula is: G MMSE=H H(HH H+ σ n 2I Nr) -1
H is a channel matrix, the element h of H IjRepresent the channel coefficients of j transmit antennas, x=[x to i root reception antenna k] T, k=1,2 ... Nt, expression transmission data, y=[y k] T, k=1,2 ... Nr, expression receives data, n=[n k] T, k=1,2 ... Nr, represent the noise of every reception antenna,
The MMSE that back detection signal-to-noise ratio computing module is used for calculating according to MMSE equalization filter module detects data
Figure BDA0000109414020000047
With formula w=1/ (1-diag (G MMSEH))-1 calculate back detection signal-to-noise ratio w,
Order module is used for according to back detection signal-to-noise ratio w each data passages being sorted from big to small, obtains ordered series of numbers data m=[m 1, m 2... m Nt], wherein
Figure BDA0000109414020000051
Choose the data channel module and be used to choose maximum back detection signal-to-noise ratio data channel m iAs the current data channel of choosing, make i=1, export corresponding back detection signal-to-noise ratio of the current data channel of choosing and corresponding detection data;
Soft demodulation module is used for basis
Figure BDA0000109414020000052
Calculate the LLR of this data channel,
Figure BDA0000109414020000053
Represent the current m that chooses iThe LLR value that data channel k bit is corresponding;
The LLR that data revaluation module is used for obtaining according to soft demodulation module carries out revaluation to the detection data of the current data channel of choosing, and obtains revaluation detection data
Figure BDA0000109414020000054
revaluation principle to do
Figure BDA0000109414020000055
P ( θ | x ^ m i ) = Π m = 1 M 1 2 ( 1 + ( 2 c k , m - 1 ) Tanh ( L k m i 2 ) ) , Wherein θ representes the constellation point value of modulator approach that this data channel adopts, c K, mThe corresponding m bit value (0/1) of expression constellation point value θ;
Interference cancellation module is used for eliminating the current interference of choosing data channel from receiving data, obtains the reception data model
Figure BDA0000109414020000057
after the interference eliminated
The residual interference estimation module is used for according to
Figure BDA0000109414020000058
the current residual interference energy value of choosing data channel after the calculating interference eliminated;
Class MMSE equalization filter module is utilized the reception data model after the interference eliminated that obtains in the interference cancellation module, calculates m I+1Detection data on the circuit-switched data passage:
Figure BDA0000109414020000059
And output
Figure BDA00001094140200000510
Type MMSE equalization filter wherein G MMSE ′ = h m i + 1 H ( H M H M H + q m i h m i h m i H + σ n 2 I Nr ) - 1 ,
H M=[h k], k=1,2 ... Nt, k ≠ m 1, m 2... m i, eliminate the data channel that had detected, the channel matrix of remaining data passage among the expression H;
The detection signal-to-noise ratio computing module is used for basis after type
Figure BDA00001094140200000512
Calculate m I+1The back detection signal-to-noise ratio of circuit-switched data passage;
Judge module is used to judge whether current selected data passage is last data channel, if the time, detection of end then, if not, then need get back to soft demodulation module.
The advantage that the present invention had:
1, existing serial interference detection mode is often thought detected data perfect rejecting from receive data; Do not consider the influence of residual interference; The present invention has overcome prior art prejudice, and the influence and the corresponding processing method of residual interference have been carried out comprehensive consideration.
2, to the ordering of detection order, the present invention has adopted the back detection signal-to-noise ratio sort method of MMSE, and this method no longer need be calculated extra ordering parameter, has reduced operand, and obtains than other sort method more performance.
3, the present invention uses alternative manner calculated for subsequent class MMSE coefficient of equalizing wave filter to reduce operand, and has considered the influence that residual interference is brought in the alternative manner, and result of calculation is more accurate, thereby brings performance boost.
4, provide with the present invention in equalization methods corresponding accurately the back detection signal-to-noise ratio the complete computation method.
Description of drawings
Fig. 1 is system model figure of the present invention;
Fig. 2 is the inventive method flow chart;
Fig. 3 is a system construction drawing of the present invention.
Embodiment
A is balanced about least mean-square error
The model of multiple input, multiple output wireless communication system can be used model description as shown in Figure 1.At transmitter terminal configuration Nt transmit antennas, receiver end configuration Nr root reception antenna.Be formulated as:
y=Hx+n (1)
Wherein H is a channel matrix, the element h of H IjRepresent the channel coefficients of j transmit antennas to i root reception antenna.This paper uses h at the back iThe i row of expression channel matrix H are promptly represented the channel coefficients vector of i transmit antennas to all reception antennas.
Acting as of linear equalization filter obtains sending the estimated value of data to receiving data through the linear equalization filter, can be formulated as:
x ^ = Gy
The citation form of least mean-square error equalization filter is:
G MMSE=E{xy H}[E{yy H}] -1(2)
In the system model of describing like formula (1) of the present invention, the least mean-square error equalization filtering can be written as:
G MMSE=H H(HH Hn 2I Nr) -1(3)
After supposing that a circuit-switched data passage (i circuit-switched data passage) is eliminated in estimation from the primary reception signal, residue receives data and is expressed as:
Figure BDA0000109414020000072
Its effect is:
Figure BDA0000109414020000074
Figure BDA0000109414020000075
According to formula 2, can get rest layers least mean-square error equalization filter and be:
G MMSE ′ = h m i + 1 H ( H M H M H + q m i h m i h m i H + σ n 2 I Nr ) - 1
H M=[h k], k=1,2 ... Nt, k ≠ m 1, m 2... m iEliminate the data channel that had detected, the channel matrix of remaining data passage among the expression H.
B is about soft bit demodulation
The output LLR of soft demodulation module is defined as:
L k m i = ln Λ = ln Σ θ ∈ { Θ : θ k = 1 } exp ( - ( x m i - θ ) 2 2 σ 2 ) - ln Σ θ ∈ { Θ : θ k = 0 } exp ( - ( x m i - θ ) 2 2 σ 2 )
Be that tentation data is made up of useful data and white Gaussian noise, after receiving data x, judge that x is that 1 probability and x are the natural logrithm value of the ratio of 0 probability.Adopt approximate linear approximation methods can be written as:
L k m i = 1 σ 2 [ min θ ∈ { Θ : θ k = 1 } ( | x m i - θ | 2 ) - min θ ∈ { Θ : θ k = 0 } ( | x m i - θ | 2 ) ] - - - ( 5 )
Wherein Θ representes the set of each constellation point on the planisphere.
Figure BDA0000109414020000083
expression back detection signal-to-noise ratio, the absolute value of the LLR of linear approximation is directly proportional with detection signal-to-noise ratio thereafter.
Method flow diagram embodiment below in conjunction with as shown in Figure 2 among the present invention is made detailed description to the present invention:
A) MMSE equalization filter module
This module is used for the MMSE that receives data is detected, and MMSE testing process described in MMSE filter such as the formula (3) is divided into following three steps:
Step 1: calculate a Nr * Nr inverse of a matrix, definition:
Figure BDA0000109414020000084
This inverse of a matrix also will use for the iteration of matrix inversion in the subsequent calculations class MMSE equalization filtering module as an output of MMSE equalization filtering module.
Step 2: calculate the MMSE equalization filter through matrix multiplication according to formula (3).
Step 3: the MMSE that obtains each circuit-switched data passage according to
Figure BDA0000109414020000091
receives the estimated value of data.
B) back detection signal-to-noise ratio computing module
This module functions is to calculate the back detection signal-to-noise ratio that detects each layer data of back through MMSE, and this signal to noise ratio can be used as the reliability of each layer data, is used for the reliability ordering to the transmission data, thus the order of specified data Interference Cancellation.
The back detection signal-to-noise ratio that detects each layer data of back through MMSE can calculate through formula (6):
w=1/(1-diag(G MMSEH))-1(6)
Wherein the diagonal entry of matrix is got in diag () expression.W will use at the soft timing of separating of constellation; W has represented the confidence level of each layer data simultaneously; Utilize vector w that data to be detected are sorted according to confidence level from big to small; The data high according to confidence level detect earlier, eliminate it can improve the counteracting serial interference method to the method for the interference of other layer data performance then.
C) order module
This module according to from big to small sequence arrangement, obtains the maximum circuit-switched data of w to the output w of back detection signal-to-noise ratio computing module, supposes it is i circuit-switched data passage.
D) choose the data channel module
This module is chosen
Figure BDA0000109414020000092
and is sent into soft demodulation module with w (i) according to the result that order module provides.And output
Figure BDA0000109414020000093
E) soft demodulation module
This module is done soft bit demodulation to the data based detection signal-to-noise ratio thereafter of sending into this module.The module effect can be represented by formula (7):
L k m i = w m i [ min θ ∈ { Θ : θ k = 1 } ( | x m i - θ | 2 ) - min θ ∈ { Θ : θ k = 0 } ( | x m i - θ | 2 ) ] - - - ( 7 )
Wherein
Figure BDA0000109414020000101
sends into the back detection signal-to-noise ratio value corresponding with data of this module.
F) data revaluation module
This module reappraises current circuit-switched data passage based on the principle (Soft Decision) of minimized average residual interference power:
Figure BDA0000109414020000102
Figure BDA0000109414020000103
Wherein
Figure BDA0000109414020000105
can do following derivation:
Assume that the first k? Bit sentenced to a probability
Figure BDA0000109414020000106
sentenced to 0 with probability and the first k? bit of LLR is defined as
L k i = ln P 1 | x r k P 0 | x r k
Have again, this bit judges that non-0 is 1, therefore has
P 1 | x r m + P 0 | x r m = 1
Derivation can get:
P 1 | x r m = e L k m i e L k m i + 1 = 1 2 ( 1 + e L k m i - 1 e L k m i + 1 ) = 1 2 ( 1 + tanh ( L k m i 2 ) )
P 0 | x r m = 1 e L k m i + 1 = 1 2 ( 1 - e L k m i - 1 e L k m i + 1 ) = 1 2 ( 1 - tanh ( L k m i 2 ) )
Therefore, receive x r, be judged to constellation point θ probability P (θ | x r) be:
P ( θ | x ^ m i ) = Π m = 1 M 1 2 ( 1 + ( 2 c k , m - 1 ) tanh ( L k m i 2 ) ) - - - ( 8 )
G) interference cancellation module
This module functions is from receive data, to eliminate the interference that i circuit-switched data passage is introduced other layer data, is formulated as:
Its effect is:
Figure BDA0000109414020000112
Figure BDA0000109414020000113
Figure BDA0000109414020000114
H) residual interference estimation module
The function of this module is after estimating that ' eliminations ' i road is disturbed, and estimates that the i circuit-switched data estimates the residual interference energy of all the other reception data.
Make
Figure BDA0000109414020000115
and be evaluated error, then
Figure BDA0000109414020000116
is the average energy of residual interference.Estimate the residual interference energy through formula (9).
Figure BDA0000109414020000117
Wherein
Figure BDA0000109414020000118
obtained by formula (8).
I) type MMSE equalization filter module
The function of this module is to after the interference of having eliminated part (a road or several roads layer emission data), utilizes the linear minimum mean-squared error equalization filter to carry out equalization filtering to the data of rest layers.The output of this module is the detection order that obtains according to order module, the estimated value of following one deck data to be checked.
Suppose and detected l circuit-switched data passage.Last time detected data i road, earlier need be according to the sequence detection j layer data of order module.Then a type MMSE detection filter device is:
G ′ MMSE = h m i + 1 H ( H M H M H + q m i h m i h m i H + σ n 2 I Nr ) - 1 - - - ( 10 )
H wherein M=[h k], k=1,2 ... Nt, k ≠ m 1, m 2... m i
Order
A 0 = HH H + σ n 2 I Nr = Σ k h k h k H + σ n 2 I Nr
A 1 = Σ k k ≠ m 1 h k h k H + σ n 2 I Nr + ( q m 1 - 1 ) h m 1 h m 1 H
= A 0 + ( q m 1 - 1 ) h m 1 h m 1 H
A 2 = A 1 + ( q m 2 - 1 ) h m 2 h m 2 H
...
Class MMSE equalization filter formula (10) is rewritten as formula (11)
G MMSE ′ = h m i + 1 H A i - 1 - - - ( 11 )
Can get according to matrix inversion lemma
A i - 1 = A i - 1 - 1 - A i - 1 - 1 h m i ( A i - 1 - 1 h m i ) H 1 q m i - 1 + h m i H ( A i - 1 - 1 h m i ) - - - ( 12 )
The inverse matrix iteration that only needs to calculate the space of matrices opened through the vector
Figure BDA0000109414020000128
of matrix inversion last time and previous iteration just can obtain this layer MMSE equalization filtering filter need inverse matrix.
J) type back detection signal-to-noise ratio computing module
The function of this module is to calculate m I+1The back detection signal-to-noise ratio of circuit-switched data passage,
Figure BDA0000109414020000129
Calculate through formula (13).
w m i + 1 = 1 / G MMSE , h m i - 1 - - - ( 13 )
K) judge module
Judge whether current selected data passage is last data channel, if the time, detection of end then, if not, then need get back to soft demodulation module.
The invention provides a system based on method, as shown in Figure 3,
Based on the counteracting serial interference detection system of many antennas bit interweaving encoding system, comprise the MMSE equalization filter module that connects successively, back detection signal-to-noise ratio computing module; Order module is chosen the data channel module, soft demodulation module; Data revaluation module, interference cancellation module, residual interference estimation module; Class MMSE equalization filter module, type back detection signal-to-noise ratio computing module and judge module, MMSE equalization filter module is used at the MIMO model according to the MMSE equalization filter
Figure BDA0000109414020000131
Calculate MMSE and detect data
Figure BDA0000109414020000132
Wherein: MIMO model: y=Hx+n, transmitting terminal configuration Nt transmit antennas, receiving terminal configuration Nr root reception antenna; G MMSEBe the MMSE equalization filter, its expression formula is: G MMSE = H H ( HH H + σ n 2 I Nr ) - 1
H is a channel matrix, the element h of H IjRepresent the channel coefficients of j transmit antennas, x=[x to i root reception antenna k] T, k=1,2 ... Nt, expression transmission data, y=[y k] T, k=1,2 ... Nr, expression receives data, n=[n k] T, k=1,2 ... Nr, represent the noise of every reception antenna,
The MMSE that back detection signal-to-noise ratio computing module is used for calculating according to MMSE equalization filter module detects data
Figure BDA0000109414020000134
With formula w=1/ (1-diag (G MMSEH))-1 calculate back detection signal-to-noise ratio w,
Said order module is used for according to back detection signal-to-noise ratio w each data passages being sorted from big to small, obtains ordered series of numbers data m=[m 1, m 2... m Nt], wherein
Figure BDA0000109414020000135
Choose the data channel module and be used to choose maximum back detection signal-to-noise ratio data channel m iAs the current data channel of choosing, make i=1, export corresponding back detection signal-to-noise ratio of the current data channel of choosing and corresponding detection data;
Soft demodulation module is used for basis
Figure BDA0000109414020000136
Calculate the LLR of this data channel,
Figure BDA0000109414020000137
Represent the current m that chooses iThe LLR value that data channel k bit is corresponding;
The LLR that data revaluation module is used for obtaining according to soft demodulation module carries out revaluation to the detection data of the current data channel of choosing, and obtains revaluation detection data
Figure BDA0000109414020000138
revaluation principle to do
Figure BDA0000109414020000141
P ( θ | x ^ m i ) = Π m = 1 M 1 2 ( 1 + ( 2 c k , m - 1 ) Tanh ( L k m i 2 ) ) , Wherein θ representes the constellation point value of modulator approach that this data channel adopts, c K, mThe corresponding m bit value (0/1) of expression constellation point value θ;
Interference cancellation module is used for eliminating the current interference of choosing data channel from receiving data, obtains the reception data model
Figure BDA0000109414020000143
after the interference eliminated
The residual interference estimation module is used for according to
Figure BDA0000109414020000144
the current residual interference energy value of choosing data channel after the calculating interference eliminated;
Class MMSE equalization filter module is utilized the reception data model after the interference eliminated that obtains in the interference cancellation module, calculates m I+1Detection data on the circuit-switched data passage:
Figure BDA0000109414020000145
And output
Figure BDA0000109414020000146
Type MMSE equalization filter wherein G MMSE ′ = h m i + 1 H ( H M H M H + q m i h m i h m i H + σ n 2 I Nr ) - 1 ,
H M=[h k], k=1,2 ... Nt, k ≠ m 1, m 2... m i, eliminate the data channel that had detected, the channel matrix of remaining data passage among the expression H;
The detection signal-to-noise ratio computing module is used for basis after type
Figure BDA0000109414020000148
Calculate m I+1The back detection signal-to-noise ratio of circuit-switched data passage;
Judge module is used to judge whether current selected data passage is last data channel, if the time, detection of end then, if not, then need get back to soft demodulation module.

Claims (3)

1. based on the counteracting serial interference detection method of many antennas bit interweaving encoding system, it is characterized in that: may further comprise the steps:
1] in the MIMO model, calculates MMSE and detect data
Figure FDA0000109414010000012
and back detection signal-to-noise ratio w according to MMSE equalization filter
Figure FDA0000109414010000011
Wherein: MIMO model: y=Hx+n, transmitting terminal configuration Nt transmit antennas, receiving terminal configuration Nr root reception antenna;
G MMSEBe the MMSE equalization filter, its expression formula is: G MMSE=H H(HH H+ σ n 2I Nr) -1
H is a channel matrix, the element h of H IjRepresent the channel coefficients of j transmit antennas, x=[x to i root reception antenna k] T, k=1,2 ... Nt, expression transmission data, y=[y k] T, k=1,2 ... Nr, expression receives data, n=[n k] T, k=1,2 ... Nr, represent the noise of every reception antenna,
Calculate back detection signal-to-noise ratio w:w=1/ (1-diag (G MMSEH))-1;
2] sort from big to small according to back detection signal-to-noise ratio w, obtain ordered series of numbers m=[m 1, m 2... m Nt], satisfy w m 1 > w m 2 > . . . > w m Nt ;
3] choose maximum back detection signal-to-noise ratio data channel m iAs the current data channel of choosing, make i=1, export corresponding back detection signal-to-noise ratio of the current data channel of choosing and corresponding detection data;
4] to step 3] in the current data channel of choosing carry out soft demodulation process, according to
Figure FDA0000109414010000014
Calculate the LLR of this data channel,
Figure FDA0000109414010000015
Represent the current m that chooses iThe LLR value that data channel k bit is corresponding;
5] according to step 4] in the LLR that obtains the detection data of the current data channel of choosing are carried out revaluation, obtain revaluation and detect data
Figure FDA0000109414010000016
Figure FDA0000109414010000017
P ( θ | x ^ m i ) = Π m = 1 M 1 2 ( 1 + ( 2 c k , m - 1 ) tanh ( L k m i 2 ) )
Wherein θ representes the constellation point value of modulator approach that this data channel adopts, c K, mThe corresponding m bit value (0/1) of expression constellation point value θ;
6] from receive data, eliminate the current interference of choosing data channel, obtain the reception data model
Figure FDA0000109414010000021
after the interference eliminated
7] the current residual interference energy value of choosing data channel after basis
Figure FDA0000109414010000022
the calculating interference eliminated;
8] utilize step 6] in reception data model after the interference eliminated that obtains, use a type MMSE equalization filter, calculate m I+1Detection data on the circuit-switched data passage:
Figure FDA0000109414010000023
And output
Figure FDA0000109414010000024
Type MMSE equalization filter wherein G MMSE ′ = h m i + 1 H ( H M H M H + q m i h m i h m i H + σ n 2 I Nr ) - 1 ,
H M=[h k], k=1,2 ... Nt, k ≠ m 1, m 2... m i, eliminate the data channel that had detected, the channel matrix of remaining data passage among the expression H;
9] basis
Figure FDA0000109414010000026
Calculate m I+1The back detection signal-to-noise ratio of circuit-switched data passage;
10] judge, if i+1 is less than Nt, i++, execution in step 4]; Otherwise if i+1 equals Nt, each data channel all detects, then finishes.
2. the counteracting serial interference detection method based on many antennas bit interweaving encoding system according to claim 1 is characterized in that: step 8] middle G ' MMSEAdopt following method to calculate:
Order A i = ( H M H M H + q m i h m i h m i H + σ n 2 I Nr ) , Class MMSE equalization filter is written as following form:
G MMSE ′ = h m i + 1 H A i - 1
Wherein
Figure FDA0000109414010000029
obtains through following alternative manner:
A i - 1 = A i - 1 - 1 - A i - 1 - 1 h m i ( A i - 1 - 1 h m i ) H 1 q m i - 1 + h m i H ( A i - 1 - 1 h m i )
A 0 - 1 = ( HH H + σ n 2 I ) - 1 .
3. according to the counteracting serial interference detection system based on many antennas bit interweaving encoding system of the said method of claim 1, it is characterized in that: comprise the MMSE equalization filter module that connects successively, back detection signal-to-noise ratio computing module; Order module is chosen the data channel module, soft demodulation module; Data revaluation module, interference cancellation module, residual interference estimation module; Class MMSE equalization filter module, type back detection signal-to-noise ratio computing module and judge module
The MMSE equalization filter module according to the model used in a MIMO MMSE equalizer filter
Figure FDA0000109414010000031
calculate MMSE detector data
Figure FDA0000109414010000032
where: MIMO models:
Y=Hx+n, transmitting terminal configuration Nt transmit antennas, receiving terminal configuration Nr root reception antenna;
G MMSEBe the MMSE equalization filter, its expression formula is: G MMSE=H H(HH H+ σ n 2I Nr) -1
H is a channel matrix, the element h of H IjRepresent the channel coefficients of j transmit antennas, x=[x to i root reception antenna k] T, k=1,2 ... Nt, expression transmission data, y=[y k] T, k=1,2 ... Nr, expression receives data, n=[n k] T, k=1,2 ... Nr, represent the noise of every reception antenna,
The MMSE that said back detection signal-to-noise ratio computing module is used for calculating according to MMSE equalization filter module detects data
Figure FDA0000109414010000033
With formula w=1/ (1-diag (G MMSEH))-1 calculate back detection signal-to-noise ratio w,
Said order module is used for according to back detection signal-to-noise ratio w each data passages being sorted from big to small, obtains ordered series of numbers data m=[m 1, m 2... m Nt], wherein
The said data channel module of choosing is used to choose maximum back detection signal-to-noise ratio data channel m iAs the current data channel of choosing, make i=1, export corresponding back detection signal-to-noise ratio of the current data channel of choosing and corresponding detection data;
Said soft demodulation module is used for basis
Figure FDA0000109414010000035
Calculate the LLR of this data channel,
Figure FDA0000109414010000036
Represent the current m that chooses iThe LLR value that data channel k bit is corresponding;
The LLR that said data revaluation module is used for obtaining according to soft demodulation module carries out revaluation to the detection data of the current data channel of choosing, and obtains revaluation detection data
Figure FDA0000109414010000041
revaluation principle to do
P ( θ | x ^ m i ) = Π m = 1 M 1 2 ( 1 + ( 2 c k , m - 1 ) Tanh ( L k m i 2 ) ) , Wherein θ representes the constellation point value of modulator approach that this data channel adopts, c K, mThe corresponding m bit value (0/1) of expression constellation point value θ;
Said interference cancellation module is used for eliminating the current interference of choosing data channel from receiving data, obtains the reception data model
Figure FDA0000109414010000044
after the interference eliminated
Said residual interference estimation module is used for according to
Figure FDA0000109414010000045
the current residual interference energy value of choosing data channel after the calculating interference eliminated;
Said type of MMSE equalization filter module utilized the reception data model after the interference eliminated that obtains in the interference cancellation module, calculates m I+1Detection data on the circuit-switched data passage:
Figure FDA0000109414010000046
And output
Figure FDA0000109414010000047
Type MMSE equalization filter wherein G MMSE ′ = h m i + 1 H ( H M H M H + q m i h m i h m i H + σ n 2 I Nr ) - 1 ,
H M=[h k], k=1,2 ... Nt, k ≠ m 1, m 2... m i, eliminate the data channel that had detected, the channel matrix of remaining data passage among the expression H;
Said type of back detection signal-to-noise ratio computing module is used for basis
Figure FDA0000109414010000049
Calculate m I+1The back detection signal-to-noise ratio of circuit-switched data passage;
Said judge module is used to judge whether current selected data passage is last data channel, if the time, detection of end then, if not, then need get back to soft demodulation module.
CN2011103651388A 2011-11-17 2011-11-17 Serial interference cancellation and detection method and system for multi-antenna bit interleaved coding system Pending CN102412931A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2011103651388A CN102412931A (en) 2011-11-17 2011-11-17 Serial interference cancellation and detection method and system for multi-antenna bit interleaved coding system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2011103651388A CN102412931A (en) 2011-11-17 2011-11-17 Serial interference cancellation and detection method and system for multi-antenna bit interleaved coding system

Publications (1)

Publication Number Publication Date
CN102412931A true CN102412931A (en) 2012-04-11

Family

ID=45914802

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2011103651388A Pending CN102412931A (en) 2011-11-17 2011-11-17 Serial interference cancellation and detection method and system for multi-antenna bit interleaved coding system

Country Status (1)

Country Link
CN (1) CN102412931A (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104506470A (en) * 2014-12-12 2015-04-08 西安电子科技大学 MMSE-IICPD (minimum mean-square error-iterative interference cancellation parallel detection) algorithm applicable to efficient calculation of symbol detection on parallel transmission system
WO2015161508A1 (en) * 2014-04-25 2015-10-29 华为技术有限公司 Method and device for estimating baseband self-interference channel response
CN105721016A (en) * 2014-12-18 2016-06-29 英特尔公司 Interference cancelation
CN103840918B (en) * 2013-11-22 2017-09-29 华为技术有限公司 Method and node in wireless communication system
CN110868273A (en) * 2019-11-29 2020-03-06 中国电子科技集团公司第五十四研究所 Joint decoding method for multipath received signals
CN112702103A (en) * 2020-11-11 2021-04-23 山东星通易航通信科技有限公司 VDES uplink received signal detection method based on serial interference elimination

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103840918B (en) * 2013-11-22 2017-09-29 华为技术有限公司 Method and node in wireless communication system
US10044531B2 (en) 2014-04-25 2018-08-07 Huawei Technologies Co., Ltd. Method for estimating response of baseband self-interference channel and apparatus
CN106464420A (en) * 2014-04-25 2017-02-22 华为技术有限公司 Method and device for estimating baseband self-interference channel response
WO2015161508A1 (en) * 2014-04-25 2015-10-29 华为技术有限公司 Method and device for estimating baseband self-interference channel response
CN106464420B (en) * 2014-04-25 2019-08-20 华为技术有限公司 Estimate the method and apparatus of base band self-interference channel response
CN104506470A (en) * 2014-12-12 2015-04-08 西安电子科技大学 MMSE-IICPD (minimum mean-square error-iterative interference cancellation parallel detection) algorithm applicable to efficient calculation of symbol detection on parallel transmission system
CN104506470B (en) * 2014-12-12 2018-12-04 西安电子科技大学 A kind of II CPD algorithm of MMSE- suitable for parallel transmission system symbol detection
CN105721016A (en) * 2014-12-18 2016-06-29 英特尔公司 Interference cancelation
US9887716B2 (en) 2014-12-18 2018-02-06 Intel Corporation Interference cancelation
CN105721016B (en) * 2014-12-18 2019-01-22 英特尔公司 Interference is eliminated
CN110868273A (en) * 2019-11-29 2020-03-06 中国电子科技集团公司第五十四研究所 Joint decoding method for multipath received signals
CN110868273B (en) * 2019-11-29 2022-03-04 中国电子科技集团公司第五十四研究所 Joint decoding method for multipath received signals
CN112702103A (en) * 2020-11-11 2021-04-23 山东星通易航通信科技有限公司 VDES uplink received signal detection method based on serial interference elimination

Similar Documents

Publication Publication Date Title
CN101322328B (en) The detection of complexity and the method for decoding and device is reduced for the receiver in communication system
CN102412931A (en) Serial interference cancellation and detection method and system for multi-antenna bit interleaved coding system
US8315343B2 (en) Multi-antenna receiver interference cancellation method and apparatus
RU2303330C1 (en) Method for receiving signal in communication system with several channels for transmitting and receiving
CN101981846B (en) Reception device, reception method, and communication system
CN1983910A (en) Method for detecting signal in multi-antenna digital wireless telecommunication system
CN102624666B (en) Cyclic coding method of multi-channel transceiving orthogonal multi-carrier underwater acoustic communication in sparse channel model
CN101442390A (en) Equilibrium acceptance method and apparatus for Turbo of spatial correlation MIMO
CN104301267B (en) The multistage iteration detection method and device of a kind of mimo wireless communication receiver
US20070206696A1 (en) Method and system for an adaptive VBLAST receiver for wireless multiple input multiple output (MIMO) detection
CN102723975B (en) Signal detection method and device of MIMO (multiple input multiple output) system
CN1943155B (en) Apparatus and method for providing transmitting sequence estimation
CN104506470B (en) A kind of II CPD algorithm of MMSE- suitable for parallel transmission system symbol detection
CN103414534B (en) A kind of generalized spatial modulation system receiver detection method in conjunction with threshold judgement
CN101136896B (en) Frequency domain iteration equalizing method based on fast Fourier transformation
CN112215335B (en) System detection method based on deep learning
CN105656823A (en) Underwater communication Turbo receiving system and underwater communication Turbo receiving method based on minimum bit error rate criterion
CN106130615A (en) The activation antenna of generalized spatial modulation system and modulation symbol combined estimation method
CN100444543C (en) Airspace filter detecting method for multiple antenna radio communication system
CN101582742A (en) Method for detecting iteration of multiple input multiple output (MIMO) system, system thereof and device thereof
CN103746728A (en) Mixed adaptive MIMO receiving and detecting method
CN100499611C (en) Inspection of blank field maximum rear-proving probability in wireless communication system
CN105610484B (en) Extensive MIMO Low-complexity iterative receiver methods
CN102487309B (en) Signal detecting method and device under MIMO (Multiple Input Multiple Output) system
JP5506109B2 (en) Perturbation decoder and decoding method in communication system and apparatus using the same

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20120411