CN115021867B - MIMO-LDPC efficient detection decoding method based on factor graph - Google Patents

MIMO-LDPC efficient detection decoding method based on factor graph Download PDF

Info

Publication number
CN115021867B
CN115021867B CN202210617770.5A CN202210617770A CN115021867B CN 115021867 B CN115021867 B CN 115021867B CN 202210617770 A CN202210617770 A CN 202210617770A CN 115021867 B CN115021867 B CN 115021867B
Authority
CN
China
Prior art keywords
node
unreliable
nodes
bit
variable
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202210617770.5A
Other languages
Chinese (zh)
Other versions
CN115021867A (en
Inventor
费泽松
李欢
郭婧
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN202210617770.5A priority Critical patent/CN115021867B/en
Publication of CN115021867A publication Critical patent/CN115021867A/en
Application granted granted Critical
Publication of CN115021867B publication Critical patent/CN115021867B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L1/00Arrangements for detecting or preventing errors in the information received
    • H04L1/004Arrangements for detecting or preventing errors in the information received by using forward error control
    • H04L1/0056Systems characterized by the type of code used
    • H04L1/0057Block codes
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03MCODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M13/00Coding, decoding or code conversion, for error detection or error correction; Coding theory basic assumptions; Coding bounds; Error probability evaluation methods; Channel models; Simulation or testing of codes
    • H03M13/03Error detection or forward error correction by redundancy in data representation, i.e. code words containing more digits than the source words
    • H03M13/05Error detection or forward error correction by redundancy in data representation, i.e. code words containing more digits than the source words using block codes, i.e. a predetermined number of check bits joined to a predetermined number of information bits
    • H03M13/11Error detection or forward error correction by redundancy in data representation, i.e. code words containing more digits than the source words using block codes, i.e. a predetermined number of check bits joined to a predetermined number of information bits using multiple parity bits
    • H03M13/1102Codes on graphs and decoding on graphs, e.g. low-density parity check [LDPC] codes
    • H03M13/1105Decoding
    • H03M13/1111Soft-decision decoding, e.g. by means of message passing or belief propagation algorithms
    • H03M13/1125Soft-decision decoding, e.g. by means of message passing or belief propagation algorithms using different domains for check node and bit node processing, wherein the different domains include probabilities, likelihood ratios, likelihood differences, log-likelihood ratios or log-likelihood difference pairs
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413MIMO systems
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Probability & Statistics with Applications (AREA)
  • Theoretical Computer Science (AREA)
  • Error Detection And Correction (AREA)

Abstract

The invention provides a factor graph-based MIMO-LDPC efficient detection decoding method, which comprises the following steps: initializing parameters; during the t-th iteration, respectively calculating the mean value and the variance of the probability estimation values transmitted to the unreliable variable nodes by the observation nodes and the symbol probability information transmitted to the unreliable bit nodes by the check nodes; calculating joint probability information of all unreliable variable nodes; obtaining a current detection decoding result according to the joint probability information; if the iteration times reach the maximum, outputting a current detection decoding result, otherwise, if the current detection decoding result is a correct code word, outputting the current result, otherwise, judging a reliable variable node and a reliable bit node, calculating the mean value and the variance transmitted to the observation node by the reliable variable node, calculating the mean value and the variance transmitted to the observation node by the unreliable variable node, and transmitting the information transmitted to the check node by the unreliable bit node to enter the next iteration. The invention can reduce the overall operation complexity while ensuring the detection decoding performance.

Description

MIMO-LDPC efficient detection decoding method based on factor graph
Technical Field
The invention relates to a factor graph-based MIMO-LDPC efficient detection decoding method, and belongs to the technical field of wireless communication.
Background
Multiple-Input Multiple-Output (MIMO) technology has been widely used in 5G, and the technology sets Multiple antennas on a transmitter and a receiver respectively to simultaneously transmit and receive radio signals, so as to improve channel capacity, transmission reliability and spectrum efficiency. Compared with a common single-input single-output system, the MIMO technology can effectively inhibit channel fading caused by multipath effect, and improves the effectiveness and reliability of a communication system. The Low-Density Parity-Check (LDPC) code is used as a long code block channel coding scheme with the error correction capability approaching to the Shannon limit, has the characteristics of good error correction capability, strong burst error resistance capability and the like, and can be well applied to a multi-antenna and multi-carrier system. The combined system of the MIMO technology and the LDPC code effectively improves the reliability, throughput and spectrum efficiency of the wireless communication system, and the research of the system and the receiving end detection decoding algorithm thereof has great significance for the development of 5G and next generation communication technologies.
In recent years, research on detection decoding algorithms in a joint system of MIMO technology and LDPC codes has also been continuously introduced. For a combined system of an MIMO technology and an LDPC code, T.L.Nar et al propose a MIMO-LDPC combined detection and decoding method, namely combining an MIMO detection factor graph and an LDPC decoding factor graph, and performing iterative updating of detection and decoding information in one iteration. Although the joint detection decoding method has excellent performance, the operation complexity is high because the MIMO detection and the LDPC decoding need to be performed simultaneously in one iteration. The QR decomposition is carried out on the channel, so that the times of iterative information updating among nodes in a factor graph are reduced, and the operation complexity of the method is reduced. However, the MIMO-LDPC joint detection decoding method based on the factor graph still has higher operation complexity, and how to reduce the operation complexity while ensuring the detection decoding performance is still an unsolved problem.
Disclosure of Invention
The invention provides a factor graph-based MIMO-LDPC efficient detection decoding method, which can reduce the overall operation complexity while ensuring the detection decoding performance and has the characteristics of low complexity and high throughput.
The invention is realized by the following technical scheme:
a MIMO-LDPC high-efficiency detection decoding method based on a factor graph, wherein the factor graph is provided with observation nodes, reliable variable nodes, unreliable variable nodes, reliable bit nodes, unreliable bit nodes and check nodes, and the method comprises the following steps:
s1, initializing and detecting decoding parameters, wherein the parameters comprise reliable variable node sets
Figure BDA0003673944380000021
Unreliable variable node set U = { x = 1 ,x 2 ,x 3 ,...,x n Get, reliable bit node set->
Figure BDA0003673944380000022
Unreliable bit node set a = { c 1 ,c 2 ,c 3 ,...,c k The reliability threshold T res Detecting the maximum iteration number I of the decoding MAX Iteration counter t =1;
s2, respectively calculating an observation node y in the factor graph during the t iteration j Passing to unreliable variable node x in set U l Mean value of the probability estimates of (1)
Figure BDA0003673944380000023
And variance->
Figure BDA0003673944380000024
And by check node e v To the unreliable bit node c in set A i Log likelihood ratio of->
Figure BDA0003673944380000025
And combines the log-likelihood ratio>
Figure BDA0003673944380000026
Conversion to unreliable variable node x l Symbol probability information of>
Figure BDA0003673944380000027
Wherein, theta r Set s = { theta ] for all possible values after signal modulation 12 ,...,θ q The value in (c) };
s3, according to the mean value
Figure BDA0003673944380000028
Variance +>
Figure BDA0003673944380000029
And symbol probability information->
Figure BDA00036739443800000210
Calculating all unreliable in the set U at the t-th iterationVariable node x l Is combined probability information->
Figure BDA00036739443800000211
Step S4, according to the joint probability information in the step S3
Figure BDA00036739443800000212
Carrying out bit judgment to obtain a current detection decoding result;
s5, judging whether the iteration times reach the maximum iteration times I MAX If yes, outputting the current detection decoding result of the step S4, otherwise, entering the step S6;
s6, judging whether the current detection decoding result is a correct code word or not by using the check matrix of the LDPC code, if so, outputting the current detection decoding result of the S4, otherwise, entering the S7;
s7, making T = T +1, and according to the reliability threshold T res And judging reliable variable nodes in the set U and reliable bit nodes in the set A, calculating the mean value and the variance transmitted to the observation nodes by the reliable variable nodes, moving the reliable variable nodes from the set U to the set V, moving the reliable bit nodes from the set A to the set B, respectively calculating the mean value and the variance transmitted to the observation nodes by the unreliable variable nodes in the updated set U and the information transmitted to the check nodes by the unreliable bit nodes in the updated set A, and entering the step S2.
Further, in the step S2, the mean value is calculated according to the following formula
Figure BDA0003673944380000031
And variance->
Figure BDA0003673944380000032
Figure BDA0003673944380000033
Figure BDA0003673944380000034
Wherein j = 1.. Multidot.m,
Figure BDA0003673944380000035
indicating the j-th symbol, N, received by the receiver x (j) Representing and observing node y j Set of connected variable nodes, N x (j) L represents N x (j) In addition to x l Set of variable nodes outside, h jl' For fading coefficients of a multipath channel>
Figure BDA0003673944380000036
Is a set N x (j) Unreliable variable node x in/l l′ To the observation node y j The average value of (a),
Figure BDA0003673944380000037
is a set N x (j) Unreliable variable node x in/l l′ To observation node y j Variance information of +>
Figure BDA0003673944380000038
Is the receiver noise variance;
by the v-th check node e v To the ith unreliable bit node c in set A i The log likelihood ratio of (d) is:
Figure BDA0003673944380000039
c i ∈N c (v) Wherein v =1 c (v) Representing and checking node e v Set of connected bit nodes, N c (v) I represents N c (v) In addition to c i A set of other bit nodes, based on the number of bit nodes in the set>
Figure BDA00036739443800000310
Set of representations N c (v) Unreliable bit node c in/i i′ To check node e v A log-likelihood ratio of (d);
according to the formula
Figure BDA00036739443800000311
r =1,2, \ 8230;, q will the log-likelihood ratio->
Figure BDA00036739443800000312
Conversion to unreliable variable node x l Is compared with the symbol probability information->
Figure BDA00036739443800000313
Wherein N is c (l) Representing unreliable variable node x in set U l A set of connected bit nodes.
Further, the step S3 specifically includes the following steps:
step S31, calculating variable node x in set U l Information of probability estimation from observation node
Figure BDA0003673944380000041
Is based on the mean value->
Figure BDA0003673944380000042
And variance->
Figure BDA0003673944380000043
The specific calculation method is as follows: />
Figure BDA0003673944380000044
Figure BDA0003673944380000045
Wherein N is y (l) Representing unreliable variable node x in set U l A set of connected observation nodes;
Figure BDA0003673944380000046
represents the complex conjugate of the channel fading coefficient;
step S32, according to the probability information
Figure BDA0003673944380000047
Is based on the mean value->
Figure BDA0003673944380000048
And variance>
Figure BDA0003673944380000049
Calculating specific probability information
Figure BDA00036739443800000410
Figure BDA00036739443800000411
Step S33, calculating the unreliable variable node x in the set U during the t iteration l Joint probability information of
Figure BDA00036739443800000412
Figure BDA00036739443800000413
Wherein Z is a normalization factor of the set s.
Further, the step S4 specifically includes the following steps:
step S41, all unreliable variable nodes x are utilized l Joint probability information of
Figure BDA00036739443800000414
Calculating its corresponding bit log-likelihood ratio>
Figure BDA00036739443800000415
Figure BDA00036739443800000416
Wherein the content of the first and second substances,
Figure BDA00036739443800000417
respectively represent and variable node x l Connected, so that bit node c i Takes a value of 1 and 0 respectively>
Figure BDA00036739443800000418
A probability set;
step S42, utilizing bit log likelihood ratio
Figure BDA0003673944380000051
A decision will be made:
Figure BDA0003673944380000052
i =1, 2.. K, resulting in the current detection decoding result ≥ s>
Figure BDA0003673944380000053
Further, in step S6, the check matrix of the LDPC code is H check If, if
Figure BDA0003673944380000054
If not, the current detection decoding result is incorrect code word, and the step S7 is entered.
Further, in step S7, the process of determining the reliable variable node and the reliable bit node includes:
step S71, calculating the unreliable variable node x in the set U according to the following formula l Mean value of
Figure BDA0003673944380000055
And variance>
Figure BDA0003673944380000056
Step S72,Calculating the unreliable variable node x according to the following formula l Reliability of (2)
Figure BDA0003673944380000057
Figure BDA0003673944380000058
r=1,2,...,q;
Step S73, for all theta in the set S r Calculated by it
Figure BDA0003673944380000059
Are all greater than the reliability threshold T res If the test result is not reliable, judging that the test result is unreliable; if a calculated value is greater than or equal to>
Figure BDA00036739443800000510
Less than a reliability threshold T res Then the unreliable variable node x is determined l Is a reliable variable node, denoted as x l″ Unreliable bit node c connected to the reliable variable node i Is a reliable bit node, denoted as c i″
Further, in the step S7, the reliable variable node x is calculated according to the following formula l″ To the observation node y j Mean value of
Figure BDA00036739443800000511
And variance>
Figure BDA00036739443800000512
y j ∈N y (l),/>
Figure BDA00036739443800000513
y j ∈N y (l)。
Further, in step S7, the unreliable variable node x in the updated set U at the t-th iteration is calculated according to the following formula l To observation node y j Mean value of
Figure BDA00036739443800000514
And variance->
Figure BDA00036739443800000515
Figure BDA00036739443800000516
y j ∈N y (l),/>
Figure BDA00036739443800000517
y j ∈N y (l)。
Further, in step S7, the updated unreliable bit node c in the set a during the t-th iteration is calculated according to the following formula i To check node e v Information of
Figure BDA0003673944380000061
Figure BDA0003673944380000062
e v ∈N e (i)。
The invention has the following beneficial effects:
1. the invention calculates the joint probability information of the unreliable variable nodes according to the mean value and the variance of the probability estimated value of the unreliable variable nodes and the log-likelihood ratio transmitted to the unreliable bit nodes by the check nodes, carries out bit judgment according to the joint probability information to obtain the current decoding result, judges the reliability of the unreliable variable nodes and the unreliable bit variables before the next iteration if the current decoding result is not a correct code word, enables the reliable variable nodes to transmit convergent feedback information to the observation nodes, moves the reliable variable nodes and the reliable bit nodes to other sets, does not update the information of the reliable variable nodes and the reliable bit nodes in the next iteration, reduces the iteration times of detection and decoding, thereby reducing the overall operation complexity while ensuring the detection and decoding performance, accelerating the convergence of the detection and decoding and obtaining higher communication throughput.
Drawings
The present invention will be described in further detail with reference to the accompanying drawings.
FIG. 1 is a flow chart of the present invention.
FIG. 2 is a diagram of a factor graph model according to the present invention.
FIG. 3 is a block diagram of the system components of the present invention.
FIG. 4 is a comparison graph of the performance simulation of the present invention.
FIG. 5 is a complexity comparison graph of the present invention.
Detailed Description
FIG. 2 is a diagram of a factor graph model adopted in the present embodiment, in which the factor graph has an observation node y j Reliable variable node x l″ Unreliable variable node x l Reliable bit node c i″ Unreliable bit node c i And check node e v The specific structure of the factor graph is the prior art. Fig. 3 shows a communication system adopted in this embodiment, which is a multi-antenna receiving end, including serial numbers 1,2 R N of (A) R A receiving antenna, a detection decoder and a signal sink, wherein the symbols received by the receiving end are marked as
Figure BDA0003673944380000071
The output symbol is recorded as ^ after detection and decoding>
Figure BDA0003673944380000072
This embodiment considers a LDPC coded single cell uplink multi-user MIMO system, where the cell has L =3 users, and each user has N t If =2 transmitting antennas, the transmitting end has N in total T =6 transmit antennas; the number of receiving antennas at the base station end is N R And (6). Considering the channel as a combined fading channel, i.e. the channel coefficient matrix +>
Figure BDA0003673944380000073
Each term G obeys a Rayleigh distribution, and xi is a diagonal element->
Figure BDA0003673944380000074
Diagonal matrix of (4), whose diagonal element->
Figure BDA0003673944380000075
Wherein it is present>
Figure BDA0003673944380000076
Obeying the lognormal distribution with the mean value of 0 and the variance of 6, and the distances between the three users and the base station are respectively D 1 =1,D 2 =1.5,D 3 =2 (km), path fading factor v =4. The LDPC codeword length is 960, the code rate is 1/3, and the modulation scheme is QPSK, i.e., q =4.
As shown in fig. 1, the factor graph-based MIMO-LDPC efficient detection decoding method includes the following steps:
s1, initializing detection decoding parameters, wherein the parameters comprise the probability of transmitting variable nodes to observation nodes
Figure BDA0003673944380000077
Mean value information &>
Figure BDA0003673944380000078
Variance information ≥>
Figure BDA0003673944380000079
Bit node passing to check node &>
Figure BDA00036739443800000710
Reliability threshold T res (ii) a Reliable variable node set->
Figure BDA00036739443800000711
Unreliable variable node set U = { x = 1 ,x 2 ,x 3 ,...,x n }; reliable set of bit nodes->
Figure BDA00036739443800000712
Unreliable bit node set a = { c = 1 ,c 2 ,c 3 ,...,c k }; detecting maximum iteration number I of decoding MAX =10, iteration counter t =1;
wherein x is l Representing the l variable node in the factor graph, wherein the value range of l is 1,2,3, n, and n is the total number of symbols sent by a sending end; y is j Representing the jth observation node in the factor graph, wherein the value range of j is 1,2, 3.., m, and m is the total number of symbols received by a receiving end; c. C i Representing the ith bit node in the factor graph, wherein the value range of i is 1,2,3, and k is the bit number generated by the information source; e.g. of the type v Representing the v-th check node in the factor graph, wherein the value range of v is 1,2, 3., w; q is the set of all possible values s = { theta after signal modulation 12 ,...,θ q The number of elements in the Chinese character, s, is the value set of the variable node; in particular, it is possible to use, for example,
Figure BDA00036739443800000713
Figure BDA0003673944380000081
wherein N is T The number of the antennas at the transmitting end is, and in addition, the reliable variable node set V and the unreliable variable node set U are not intersected with each other, that is: />
Figure BDA0003673944380000082
Similarly, reliable bit node set B and non-reliable bit node set A satisfy >>
Figure BDA0003673944380000083
S2, during the t-th iteration, calculating the observed node y in the factor graph respectively j Passing to unreliable variable node x in set U l Mean value of the probability estimates of (1)
Figure BDA0003673944380000084
And variance->
Figure BDA0003673944380000085
And by check node e v To the unreliable bit node c in set A i Is greater than or equal to>
Figure BDA0003673944380000086
And taking the log likelihood ratio->
Figure BDA0003673944380000087
Conversion to unreliable variable node x l Symbol probability information of>
Figure BDA0003673944380000088
Wherein, theta r Set s = { theta ] for all possible values after signal modulation 12 ,...,θ q The value in (c) };
specifically, the observation node y is calculated according to the following formula j To unreliable variable node x in set U l Of the mean of the probability estimates
Figure BDA0003673944380000089
And variance->
Figure BDA00036739443800000810
/>
Figure BDA00036739443800000811
Figure BDA00036739443800000812
Wherein j = 1.. Multidot.m,
Figure BDA00036739443800000813
indicating the j-th symbol, N, received by the receiver x (j) Representing and observing node y j Set of connected variable nodes, N x (j) L represents N x (j) In addition to x l Set of variable nodes outside, h jl' For fading coefficients of a multipath channel>
Figure BDA00036739443800000814
Is a set N x (j) Unreliable variable node x in/l l′ To observation node y j The average value of (a),
Figure BDA00036739443800000815
is a set N x (j) Unreliable variable node x in/l l′ To observation node y j Variance information of +>
Figure BDA00036739443800000816
Is the receiver noise variance;
by the v-th check node e v To the ith unreliable bit node c in set A i The log-likelihood ratio of (a) is:
Figure BDA00036739443800000817
c i ∈N c (v) Wherein v =1 c (v) Representing and checking nodes e v Set of connected bit nodes, N c (v) I represents N c (v) In addition to c i A set of other bit nodes, based on the number of bit nodes in the set>
Figure BDA00036739443800000818
Set of representations N c (v) Unreliable bit node c in/i i′ To check node e v A log-likelihood ratio of;
according to the formula
Figure BDA0003673944380000091
r =1,2, \ 8230;, q will the log-likelihood ratio->
Figure BDA0003673944380000092
Conversion to unreliable variable node x l Is compared with the symbol probability information->
Figure BDA0003673944380000093
Wherein N is c (l) Representing unreliable variable nodes x in the set U l A set of connected bit nodes.
Step S3. According to the mean value
Figure BDA0003673944380000094
Variance->
Figure BDA0003673944380000095
And symbol probability information>
Figure BDA0003673944380000096
Calculating all unreliable variable nodes x in the set U at the t iteration l Is combined probability information->
Figure BDA0003673944380000097
The method specifically comprises the following steps:
step S31, calculating unreliable variable node x in set U l Processing probability estimation information from observation nodes
Figure BDA0003673944380000098
In (d) is based on the mean value>
Figure BDA0003673944380000099
And variance->
Figure BDA00036739443800000910
The specific calculation method is as follows:
Figure BDA00036739443800000911
Figure BDA00036739443800000912
wherein N is y (l) Representing unreliable variable node x in set U l A set of connected observation nodes;
Figure BDA00036739443800000913
represents the complex conjugate of the channel fading coefficient;
step S32, according to the probability information
Figure BDA00036739443800000914
Is based on the mean value->
Figure BDA00036739443800000915
And variance->
Figure BDA00036739443800000916
Calculating concrete probability information
Figure BDA00036739443800000917
Figure BDA00036739443800000918
/>
Step S33, calculating a variable node x in the set U during the t iteration l Joint probability information of
Figure BDA00036739443800000919
Figure BDA00036739443800000920
Wherein Z is a normalization factor of the set s.
Step S4, according to the joint probability information in the step S3
Figure BDA00036739443800000921
Carrying out bit judgment to obtain a current detection decoding result;
the method specifically comprises the following steps:
step S41, all unreliable variable nodes x are utilized l Joint probability information of
Figure BDA0003673944380000101
Calculate its corresponding bit log likelihood ratio->
Figure BDA0003673944380000102
Figure BDA0003673944380000103
Wherein the content of the first and second substances,
Figure BDA0003673944380000104
respectively represent and variable node x l Connected, so that bit node c i Takes a value of 1 and 0 respectively>
Figure BDA0003673944380000105
A probability set;
step S42, utilizing bit log likelihood ratio
Figure BDA0003673944380000106
A decision will be made:
Figure BDA0003673944380000107
i =1, 2.. K, resulting in the current detection decoding result ≥ s>
Figure BDA0003673944380000108
S5, judging whether the iteration times reach the maximum iteration times I MAX =10, if yes, outputting the current detection decoding result of the step S4, otherwise, entering the step S6;
s6, judging whether the current detection decoding result is a correct code word or not by using a check matrix of the LDPC code, if so, outputting the current detection decoding result of the S4, and otherwise, entering the S7;
specifically, the check matrix of the LDPC code is H check If at all
Figure BDA0003673944380000109
If not, the current detection decoding result is incorrect code word, and the step S7 is entered.
S7, enabling T = T +1, and enabling the reliability threshold T to be used according to res Judging reliable variable nodes in the set U and reliable bit nodes in the set A, calculating the mean value and the variance transmitted to observation nodes by the reliable variable nodes, moving reliable node variables from the set U to the set V, moving the reliable bit nodes from the set A to the set B, respectively calculating the mean value and the variance transmitted to the observation nodes by the unreliable variable nodes in the updated set U and the log-likelihood ratio transmitted to the check nodes by the unreliable bit nodes in the updated set A, and entering the step S2;
the method specifically comprises the following steps:
step S71, calculating the unreliable variable node x in the set U according to the following formula l Mean value of
Figure BDA0003673944380000111
And variance
Figure BDA0003673944380000112
Step S72, calculating the unreliable variable node x according to the following formula l Reliability of (2)
Figure BDA0003673944380000113
/>
Figure BDA0003673944380000114
r=1,2,...,q;
Step S73, for all theta in the set S r Calculated by it
Figure BDA0003673944380000115
Are all greater than the reliability threshold T res If the test result is not reliable, judging that the test result is unreliable; if a calculated value is greater than or equal to>
Figure BDA0003673944380000116
Less than a reliability threshold T res Then the unreliable variable node x is determined l Is a reliable variable node, denoted as x l″ An unreliable bit node c connected to the reliable variable node i Is a reliable bit node, denoted as c i″
Step S74, calculating reliable variable node x according to the following formula l″ To observation node y j Mean value of
Figure BDA0003673944380000117
And variance->
Figure BDA0003673944380000118
y j ∈N y (l),/>
Figure BDA0003673944380000119
y j ∈N y (l);
Step S75, calculating the unreliable variable node x in the updated set U during the t iteration according to the following formula l To the observation node y j Mean value of
Figure BDA00036739443800001110
And variance->
Figure BDA00036739443800001111
Figure BDA00036739443800001112
y j ∈N y (l),
Figure BDA00036739443800001113
y j ∈N y (l);
Step S76, calculating the unreliable bit node c in the updated set A during the t iteration according to the following formula i To check node e v Information of (2)
Figure BDA00036739443800001114
Fig. 4 is a diagram showing simulation effect of block error rate when the received signal-to-noise ratio changes in the embodiment, where the abscissa in the diagram is the received signal-to-noise ratio ρ r The calculation formula is
Figure BDA00036739443800001115
Wherein, P t =1 is the symbol normalized power of the transmitting end; the ordinate represents the Block Error rate (BLER). The simulation experiment carries out comparative analysis on the two methods: 1) Using reliability threshold T in detecting decoding process res Judging the nodes, and setting reliable nodes to transmit convergence information to the observation nodes, namely the method; 2) A MIMO-LDPC detection decoding method without reliable node judgment; in addition, in the simulation experiment process, the reliability thresholds of the method are set to be T respectively res =0、T res =1×10 -7 And (6) carrying out simulation.
It can be seen from fig. 4 that when the reliable node feeds back the convergence information, the BLER performance of decoding detection is significantly improved; at T res When the node is not equal to 0, the method has 0.4dB gain compared with a detection decoding method without reliable node judgment; at T res =1×10 -7 Compared with a detection decoding method without reliable node judgment, the method has 1.7dB gain, and the system throughput is effectively improved; meanwhile, as the SNR of the receiving end increases, the BLER gain also gradually increases.
As shown in fig. 5, the method adopts the complexity comparison between different reliability thresholds and the unreliable node determination method, the abscissa is different complexity comparison types, and the ordinate is the percentage of the reduced operation times of the MIMO-LDPC detection decoding method compared with the unreliable node determination in each frame in this embodiment; when the reliable node transmits convergence information to the observation node, the total iteration times of detection decoding are greatly reduced, and the operation complexity of the detection decoding method is further reduced.
The above description is only a preferred embodiment of the present invention, and therefore should not be taken as limiting the scope of the invention, and the equivalent variations and modifications made in the claims and the description of the present invention should be included in the scope of the present invention.

Claims (5)

1. A MIMO-LDPC high-efficiency detection decoding method based on a factor graph is provided, wherein the factor graph comprises observation nodes, reliable variable nodes, unreliable variable nodes, reliable bit nodes, unreliable bit nodes and check nodes, and is characterized in that: the method comprises the following steps:
s1, initializing and detecting decoding parameters, wherein the parameters comprise reliable variable node sets
Figure QLYQS_1
Unreliable variable node set U = { x = 1 ,x 2 ,x 3 ,...,x n Get, reliable bit node set->
Figure QLYQS_2
Unreliable bit node set a = { c 1 ,c 2 ,c 3 ,...,c k }, reliability threshold T res Detecting the maximum iteration number I of decoding MAX Iteration counter t =1;
s2, respectively calculating an observation node y in the factor graph during the t iteration j Passing to unreliable variable node x in set U l Of the mean of the probability estimates
Figure QLYQS_3
And variance->
Figure QLYQS_4
And by check node e v To the unreliable bit node c in set A i Is greater than or equal to>
Figure QLYQS_5
And combines the log-likelihood ratio>
Figure QLYQS_6
Conversion to unreliable variable node x l Symbol probability information of>
Figure QLYQS_7
Wherein, theta r Set s = { theta ] for all possible values after signal modulation 12 ,...,θ q The value in (c) };
s3, according to the mean value
Figure QLYQS_8
Variance->
Figure QLYQS_9
And symbol probability information->
Figure QLYQS_10
Calculating all unreliable variable nodes x in the set U at the t iteration l Is combined probability information->
Figure QLYQS_11
Step S4, according to the joint probability information in the step S3
Figure QLYQS_12
Carrying out bit judgment to obtain a current detection decoding result;
s5, judging whether the iteration times reach the maximum iteration times I MAX If yes, outputting the current detection decoding result of the step S4, otherwise, entering the step S6;
s6, judging whether the current detection decoding result is a correct code word or not by using a check matrix of the LDPC code, if so, outputting the current detection decoding result of the S4, otherwise, entering the S7;
s7, making T = T +1, and according to the reliability threshold T res Judging reliable variable nodes in the set U and reliable bit nodes in the set A, calculating the mean value and the variance transmitted to the observation nodes by the reliable variable nodes, moving the reliable variable nodes from the set U to the set V, moving the reliable bit nodes from the set A to the set B, and respectively calculating the mean value and the variance transmitted to the observation nodes by the unreliable variable nodes in the updated set U and the signal transmitted to the check nodes by the unreliable bit nodes in the updated set AThen, the step S2 is carried out;
the step S3 specifically includes the following steps:
step S31, calculating variable node x in set U l Information of probability estimation from observation node
Figure QLYQS_13
In (d) is based on the mean value>
Figure QLYQS_14
And variance>
Figure QLYQS_15
The specific calculation method is as follows:
Figure QLYQS_16
Figure QLYQS_17
wherein N is y (l) Representing unreliable variable node x in set U l A set of connected observation nodes;
Figure QLYQS_18
represents the complex conjugate of the channel fading coefficient;
step S32, according to the probability information
Figure QLYQS_19
Is based on the mean value->
Figure QLYQS_20
And variance->
Figure QLYQS_21
Calculating specific probability information>
Figure QLYQS_22
/>
Figure QLYQS_23
Step S33, calculating unreliable variable node x in set U during the t iteration l Joint probability information of
Figure QLYQS_24
Figure QLYQS_25
Wherein Z is a normalization factor of the set s;
the step S4 specifically includes the following steps:
step S41, all unreliable variable nodes x are utilized l Joint probability information of
Figure QLYQS_26
Calculate its corresponding bit log likelihood ratio->
Figure QLYQS_27
Figure QLYQS_28
Wherein, the first and the second end of the pipe are connected with each other,
Figure QLYQS_29
respectively represent and variable node x l Connected, so that the bit node c i Takes a value of 1 and 0 respectively>
Figure QLYQS_30
A probability set;
step S42, utilizing bit log likelihood ratio
Figure QLYQS_31
A decision will be made:
Figure QLYQS_32
get the current detection decoding result>
Figure QLYQS_33
In step S7, the process of determining the reliable variable node and the reliable bit node includes:
step S71, calculating the unreliable variable node x in the set U according to the following formula l Mean value of
Figure QLYQS_34
And variance->
Figure QLYQS_35
Figure QLYQS_36
Step S72, calculating the unreliable variable node x according to the following formula l Reliability of (2)
Figure QLYQS_37
Figure QLYQS_38
Step S73, for all theta in the set S r Calculated by it
Figure QLYQS_39
Are all greater than the reliability threshold T res If the test result is not reliable, judging that the test result is unreliable; if a certain calculated->
Figure QLYQS_40
Less than a reliability threshold T res Then the unreliable variable node x is determined l Is a reliable variable node, denoted as x l″ Unreliable bit node c connected to the reliable variable node i Is a reliable bit node, denoted as c i″
2. The factor graph-based MIMO-LDPC efficient detection and decoding method according to claim 1, wherein: in the step S2, the mean value is calculated according to the following formula
Figure QLYQS_41
And variance->
Figure QLYQS_42
Figure QLYQS_43
Figure QLYQS_44
Wherein j = 1.. Multidot.m,
Figure QLYQS_45
indicating the j-th symbol, N, received by the receiver x (j) Representing and observing node y j Set of connected variable nodes, N x (j) L represents N x (j) In addition to x l Set of variable nodes outside, h jl' For the fading coefficient of the multipath channel, is->
Figure QLYQS_46
Is a set N x (j) Unreliable variable node x in/l l′ To the observation node y j The average value of (a),
Figure QLYQS_47
is a set N x (j) Unreliable variable node x in/l l′ To observation node y j Variance information of +>
Figure QLYQS_48
As noise to the receiverVariance;
by the v-th check node e v To the ith unreliable bit node c in set A i The log likelihood ratio of (d) is:
Figure QLYQS_49
wherein v =1 c (v) Representing and checking node e v Set of connected bit nodes, N c (v) I represents N c (v) In addition to c i A set of other bit nodes, based on the number of bit nodes in the set>
Figure QLYQS_50
Set of representations N c (v) Unreliable bit node c in/i i′ To check node e v A log-likelihood ratio of;
according to the formula
Figure QLYQS_51
Pick up the log likelihood ratio>
Figure QLYQS_52
Conversion to unreliable variable node x l Is compared with the symbol probability information->
Figure QLYQS_53
Wherein N is c (l) Representing unreliable variable nodes x in the set U l A set of connected bit nodes.
3. The factor-graph-based MIMO-LDPC efficient detection decoding method according to claim 1 or 2, wherein: in the step S6, the check matrix of the LDPC code is H check If at all
Figure QLYQS_54
If the current detection decoding result is correct code word, the current detection decoding result is output, otherwise, the current detection decoding result is incorrect code word, and step S7 is entered.
4. The factor-graph-based MIMO-LDPC efficient detection decoding method according to claim 1 or 2, wherein: in the step S7, the reliable variable node x is calculated according to the following formula l″ To the observation node y j Mean value of
Figure QLYQS_55
And variance>
Figure QLYQS_56
Figure QLYQS_57
Figure QLYQS_58
5. The factor-graph-based MIMO-LDPC efficient detection decoding method according to claim 1 or 2, wherein: in step S7, the unreliable variable node x in the updated set U during the t-th iteration is calculated according to the following formula l To observation node y j Mean value of
Figure QLYQS_59
And variance->
Figure QLYQS_60
Figure QLYQS_61
/>
CN202210617770.5A 2022-06-01 2022-06-01 MIMO-LDPC efficient detection decoding method based on factor graph Active CN115021867B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210617770.5A CN115021867B (en) 2022-06-01 2022-06-01 MIMO-LDPC efficient detection decoding method based on factor graph

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210617770.5A CN115021867B (en) 2022-06-01 2022-06-01 MIMO-LDPC efficient detection decoding method based on factor graph

Publications (2)

Publication Number Publication Date
CN115021867A CN115021867A (en) 2022-09-06
CN115021867B true CN115021867B (en) 2023-04-14

Family

ID=83072645

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210617770.5A Active CN115021867B (en) 2022-06-01 2022-06-01 MIMO-LDPC efficient detection decoding method based on factor graph

Country Status (1)

Country Link
CN (1) CN115021867B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106330199A (en) * 2016-08-22 2017-01-11 电子科技大学 Factor graph-based sparse code multiple access (SCMA) and low density parity check (LDPC) joint detection decoding algorithm and apparatus
CN107094026A (en) * 2017-04-10 2017-08-25 东南大学 The figure of NB LDPC codings merges detection interpretation method
CN111342935A (en) * 2020-03-04 2020-06-26 西安电子科技大学 High-performance soft-decision decoding method based on QC-LDPC code
WO2022041031A1 (en) * 2020-08-27 2022-03-03 华为技术有限公司 Ldpc code construction method and communication apparatus

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050265387A1 (en) * 2004-06-01 2005-12-01 Khojastepour Mohammad A General code design for the relay channel and factor graph decoding
CN101997639B (en) * 2009-08-10 2014-10-22 中兴通讯股份有限公司 Iterative receiver method of low density parity check-multi-input/output communication system
CN106301517B (en) * 2016-08-10 2018-07-06 清华大学 Based on the satellite multi-beam joint-detection and interpretation method it is expected to propagate and system
CN109889208A (en) * 2019-01-04 2019-06-14 浙江大学 LDPC channel decoding method based on normalization minimum-sum algorithm in NAVDAT

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106330199A (en) * 2016-08-22 2017-01-11 电子科技大学 Factor graph-based sparse code multiple access (SCMA) and low density parity check (LDPC) joint detection decoding algorithm and apparatus
CN107094026A (en) * 2017-04-10 2017-08-25 东南大学 The figure of NB LDPC codings merges detection interpretation method
CN111342935A (en) * 2020-03-04 2020-06-26 西安电子科技大学 High-performance soft-decision decoding method based on QC-LDPC code
WO2022041031A1 (en) * 2020-08-27 2022-03-03 华为技术有限公司 Ldpc code construction method and communication apparatus

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Zesong Fei;Jing Guo.RIS-Assisted Spectrum Sharing Between MIMO Radar and MU-MISO Communication Systems.《 IEEE Wireless Communications Letters》.2021,全文. *
薛英健,项海格.基于低密度校验码编码的MIMO系统方案.电子与信息学报.2004,(10),全文. *

Also Published As

Publication number Publication date
CN115021867A (en) 2022-09-06

Similar Documents

Publication Publication Date Title
US8462893B2 (en) Receiver and method for decoding a receive signal in a wireless communication system
US20070127603A1 (en) Method for the soft bit metric calculation with linear MIMO detection for LDPC codes
US7702050B2 (en) Method and system for an adaptive VBLAST receiver for wireless multiple input multiple output (MIMO) detection
EP2410707A1 (en) Space time coding method in orthogonal network and relay transmission system
CN100505713C (en) MIMO interate detection method of low complexity
CN101997639B (en) Iterative receiver method of low density parity check-multi-input/output communication system
CN102158311B (en) A kind of iteration detection method optimizing serial interference elimination order
CN101442394B (en) Network encode collaboration communication method capable of iteratively decoding
Narasimhan et al. Factor graph based joint detection/decoding for LDPC coded large-MIMO systems
CN101026435A (en) Low-complexity maximum likelihood detecting method and device for communication system
CN101582742A (en) Method for detecting iteration of multiple input multiple output (MIMO) system, system thereof and device thereof
CN104009822A (en) Novel demodulation correcting method for estimating narrowband-containing interference based on non-ideal channel
CN115021867B (en) MIMO-LDPC efficient detection decoding method based on factor graph
CN103746772A (en) Optimization method of demodulator output soft information for LDPC (Low Density Parity Code) code modulation system
CN108512580B (en) Large-scale multi-user MIMO iterative detection method suitable for low-precision quantization
CN110336640A (en) A kind of long system LT code interpretation method of short code
CN101009533A (en) Detection method for MIMO system
CN101005298B (en) Method and system for signal processing
KR101543621B1 (en) Apparatus and method for detecting signal in multiple input multiple output wireless communication system
CN101582743B (en) MIMO detection method and system used for iterative receivers
Benzid et al. LDPC for receive antennas selection in massive MiMo
CN113078983B (en) LLR (LLR) calculation method based on double-Gaussian approximation
CN109245858B (en) Improved joint network-Turbo coding method based on decoding forwarding
US8767896B2 (en) Methods and apparatus for MIMO detection
Zhou et al. Low-complexity sphere decoding for polar-coded MIMO systems

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant