CN112929128B - 一种基于置信度传播的mimo检测方法、装置 - Google Patents

一种基于置信度传播的mimo检测方法、装置 Download PDF

Info

Publication number
CN112929128B
CN112929128B CN202110148962.1A CN202110148962A CN112929128B CN 112929128 B CN112929128 B CN 112929128B CN 202110148962 A CN202110148962 A CN 202110148962A CN 112929128 B CN112929128 B CN 112929128B
Authority
CN
China
Prior art keywords
value
symbol
node
values
observation
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
CN202110148962.1A
Other languages
English (en)
Other versions
CN112929128A (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.)
Network Communication and Security Zijinshan Laboratory
Original Assignee
Network Communication and Security Zijinshan Laboratory
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 Network Communication and Security Zijinshan Laboratory filed Critical Network Communication and Security Zijinshan Laboratory
Priority to CN202110148962.1A priority Critical patent/CN112929128B/zh
Publication of CN112929128A publication Critical patent/CN112929128A/zh
Application granted granted Critical
Publication of CN112929128B publication Critical patent/CN112929128B/zh
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/0045Arrangements at the receiver end
    • H04L1/0054Maximum-likelihood or sequential decoding, e.g. Viterbi, Fano, ZJ algorithms
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Artificial Intelligence (AREA)
  • Radio Transmission System (AREA)

Abstract

本发明公开了一种基于置信度传播的MIMO检测方法、装置,涉及通信领域,其中MIMO检测方法包括:获取MIMO的系统参数;观察节点根据所述系统参数以及各个符号节点的先验概率计算各个接收符号的后验概率的LLR值,并将其传回给各个符号节点;根据观察节点反馈回来的信息,计算符号节点的LLR值;根据观察节点的后验概率的LLR值得到先验概率的LLR值;根据所述先验概率的LLR值计算出符号节点要发给各个观察节点的后验概率值,并开始新一次迭代;迭代完成后,根据需求输出符号;该基于硬件平台的MIMO BP算法的实现方式,采用线性插值的方式进行近似计算,避免了e指数查表需要存储多张表以及地址读取冲突问题,减少了硬件资源的损耗。

Description

一种基于置信度传播的MIMO检测方法、装置
技术领域
本发明涉及通信领域,具体而言,涉及一种MIMO检测方法、装置。
背景技术
随着直播、高清电话会议,物联网通信等高速率高容量低延迟的应用场景的需要,5G比4G具有更高速率和更高的网络容量。为了满足人类的现实需求,另一方面解决频率资源短缺的问题,多输入多输出(MIMO)技术成为5G的关键技术之一。MIMO技术是在发送端和接收端都部署多根天线进行无线通信的技术,主要的研究方向包括空间复用,空间分集,预编码、信号检测等领域。所以MIMO检测算法的好坏决定着整个MIMO系统的整体表现。
现有的MIMO检测算法性能最好的是最大似然检测算法,但是其复杂度会随着系统的天线数和调制阶数指数级增长,其他的线性检测算法虽然复杂度要低于最大似然检测,但是通常是以付出较大的性能损失作为代价的。而另一种非线性检测算法,基于置信度传播(BP)的检测算法在性能上可以逼近最大似然检测,而在复杂度上可以低于最大似然检测。BP检测算法由于其优秀的性能,得到广泛的关注,但是由于其在进行概率更新时涉及到指数运算与除法运算,同样具有一定的复杂度。在硬件实现过程中,指数运算采用查表的方式时,首先要存储的数据要从0到正无穷,这是一个很大的数据量存储量,同时,由于每次都只能读一个地址,为了防止地址冲突,进而需要同时存储多张e指数表,对于硬件资源是一个巨大的挑战。
发明内容
为了解决上述问题,本发明提供一种复杂度低、硬件资源损耗小而且能够避免了e指数查表需要同时存储多张表以及地址读取冲突问题的MIMO检测方法以及装置。
为了实现上述目的,本发明一方面提供一种基于置信度传播的MIMO检测方法,包括以下步骤:
获取MIMO的系统参数;所述系统参数包括发送天线数,接收天线数,系统的调制阶数,接收天线所接收到的发射天线经过信道传输和噪声干扰的接收信号、对应的信道信息以及噪声方差;
观察节点根据所述系统参数以及各个符号节点的先验概率计算各个接收符号的后验概率的LLR值,并将其传回给各个符号节点;
根据观察节点的后验概率的LLR值得到先验概率的LLR值;
根据所述先验概率的LLR值计算出符号节点要发给各个观察节点的先验概率值,并开始新一轮的先验概率的迭代;
当达到最大迭代次数时,根据需求输出符号;
其中所述根据所述先验概率的LLR值计算出符号节点要发给各个观察节点的先验概率值包括:
对于αi,j(s),从
Figure GDA0003854009850000021
选出最大值αmax,然后[αmaxi,j(sk)]乘以常数log2e,得到一个新的数A=u+v,v为小数部分数值,u为整数部分数值,将1-0.5v右移y位,得到
Figure GDA0003854009850000022
的值;
将得到的
Figure GDA0003854009850000023
的值进行累加得到F,根据F最高位相应的2的幂次,确定w的值,将F右移w位得到m的值,则根据ln(F)=ln2·(w+m-1)完成对数运算;
maxi,j(sk)]的值与ln(F)求和,再乘以常数log2e得到一个新的数值A′=u′+v′,v′为小数部分数值,u′为整数部分数值,将1-0.5v′右移u′位,即可得到相对应的先验概率值;
其中,αi,j(s)表示第i个符号节点到第j个观察节点的先验信息的LLR值;
Figure GDA0003854009850000024
s是发送符号值的集合;系统调制阶数为log2Q;F表示所有
Figure GDA0003854009850000025
的值累加得到的数值;w表示F最高位相应的2的幂次的数值;m表示将F右移w位得到的数值。
可选的,所述观察节点根据符号节点的先验概率计算各个接收符号的后验概率的LLR值,并将其传回给各个符号节点,具体包括:
计算观察节点j所收信息的均值
Figure GDA0003854009850000026
1≤j≤2N和方差
Figure GDA0003854009850000027
Figure GDA0003854009850000028
1≤j≤2N;
计算其他符号节点和信道带来的噪声干扰的均值μj,i=μj-hj,isTpi,j,1≤i≤2M,1≤j≤2N和方差
Figure GDA0003854009850000029
1≤i≤2M,1≤j≤2N;
计算观察节点的后验概率的LLR值,
Figure GDA00038540098500000210
式中,M为发送天线数;N为接收天线数;hj,i为信道矩阵H的内部元素;i表示第i个符号节点;j表示第j个观察节点;
Figure GDA00038540098500000211
s是发送符号值的集合;系统调制阶数为log2Q;(*)T表示对矩阵做转置;pi,j为第l次迭代中符号节点i到观察节点j传递的先验概率;σ2为信道噪声;rj为天线接收的信号Y=HX+σ2的内部元素。
可选的,计算过程中的指数运算采用以下方法:
对于输入的值xi,乘以一个常数log2e,得到xilog2e,将其整数部分和小数部分分离,并将小数部分右移一位,1减去右移一位后的小数值再右移整数位即得到对应的指数运算的值。
可选的,计算过程中的对数运算采用以下方法:
对于输入的值xi,找到其转换成二进制表示的整数最高位1的位置w,将xi右移w位得到m,再与w-1相加,最后乘以常数ln2完成对数运算。
另一方面,本发明还提供一种基于置信度传播的MIMO检测装置,包括:
获取单元,用于获取MIMO的系统参数;所述系统参数包括发送天线数,接收天线数,系统的调制阶数,接收天线所接收到的发射天线经过信道传输和噪声干扰的接收信号、对应的信道信息以及噪声方差;
第一计算单元,用于观察节点根据所述系统参数以及各个符号节点的先验概率计算各个接收符号的后验概率的LLR值,并将其传回给各个符号节点;
第二计算单元,用于根据观察节点的后验概率的LLR值得到先验概率的LLR值;
第三计算单元,用于根据所述先验概率的LLR值计算出符号节点要发给各个观察节点的先验概率值,并开始新一轮的先验概率的迭代;
选取单元,用于当达到最大迭代次数时,根据需求输出符号;
其中所述根据所述先验概率的LLR值计算出符号节点要发给各个观察节点的先验概率值包括:
对于αi,j(s),从
Figure GDA0003854009850000031
选出最大值αmax,然后[αmaxi,j(sk)]乘以常数log2e,得到一个新的数A=u+v,v为小数部分数值,u为整数部分数值,将1-0.5v右移u位,得到
Figure GDA0003854009850000032
的值;
将得到的
Figure GDA0003854009850000033
的值进行累加得到F,根据F最高位相应的2的幂次,确定w的值,将F右移w位得到m的值,则根据ln(F)=ln2·(w+m-1)完成对数运算;
maxi,j(sk)]的值与lnF)求和,再乘以常数log2e得到一个新的数值A′=u′+v′,v′为小数部分数值,u′为整数部分数值,将1-0.5v′右移u′位,即可得到相对应的先验概率值;
其中,αi,j(s)表示第i个符号节点到第j个观察节点的先验信息的LLR值;
Figure GDA0003854009850000034
s是发送符号值的集合;系统调制阶数为log2Q;F表示所有
Figure GDA0003854009850000041
的值累加得到的数值;w表示F最高位相应的2的幂次的数值;m表示将F右移w位得到的数值。
在上述MIMO检测装置中,可选的,所述第一计算单元包括LLR计算单元,所述LLR计算单元用于计算观察节点的LLR值,具体包括:
观察节点接收到的信息rj减去均值uj,i,信道hj,i左移一位与rjj,i相乘,得到的结果2hj,i(rjj,i),再与(s-s1)做乘法;信道hj,i与自身做乘法,再与s2-s1 2做乘法,得到
Figure GDA0003854009850000042
两者做一次减法得到值
Figure GDA0003854009850000043
除以
Figure GDA0003854009850000044
得到观察节点的LLR值;
其中,rj为天线接收的信号Y=HX+σ2的内部元素,hj,i为信道矩阵H的内部元素;i表示第i个符号节点;j表示第j个观察节点,σ2为信道噪声。
在上述MIMO检测装置中,可选的,所述第三计算单元包括指数计算,所述指数计算单元用于对整个计算过程中的指数进行运算,具体的运算步骤如下:
对于输入的值xi,乘以一个常数log2e,得到x1log2e将其整数部分和小数部分分离,并将小数部分右移一位,1减去右移一位后的小数值再右移整数位即得到对应的指数运算的值。
在上述MIMO检测装置中,可选的,所述第三计算单元包括对数计算,所述对数计算单元用于对整个计算过程中的对数进行运算,具体的运算步骤如下:
对于输入的值xi,找到其转换成二进制表示的整数最高位1的位置w,将xi右移w位得到m,再与w-1相加,最后乘以常数ln2完成对数运算。
本发明对于MIMO置信度传播检测算法中的概率更新计算公式,采用线性插值的方式进行近似计算,避免了e指数查表需要同时存储多张表以及地址读取冲突问题,使得在确保算法性能的前提下,减少硬件存储单元和计算单元的损耗。
附图说明
图1是本发明中一种基于置信度传播的MIMO检测方法的流程图;
图2是本发明中一种基于置信度传播的MIMO检测装置的结构图;
图3是本发明中一种基于置信度传播的MIMO检测装置中均值计算模块的结构图;
图4是本发明中一种基于置信度传播的MIMO检测装置中方差计算模块的示意图;
图5是本发明中一种基于置信度传播的MIMO检测装置中LLR计算模块的示意图;
图6是本发明中一种基于置信度传播的MIMO检测装置中概率计算模块的结构图;
图7是本发明中一种基于置信度传播的MIMO检测装置中指数计算模块的示意图;
图8是本发明中一种基于置信度传播的MIMO检测装置中对数计算模块的示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
参照图1,本实施例提供一种基于置信度传播的MIMO检测方法,包括以下步骤:
S10:获取MIMO的系统参数;
具体的,系统参数包括发送天线数,接收天线数,系统的调制阶数,接收天线所接收到的发射天线经过信道传输和噪声干扰的接收信号、对应的信道信息以及噪声方差。
S20:观察节点根据所述系统参数以及各个符号节点的先验概率计算各个观察节点的后验概率的LLR值,并将其传回给各个符号节点;发送天线数M,接收天线数N,信道矩阵的维度XN×M,将MIMO系统从复数域扩展到实数域:
X1×2M=[realX1×M,imagX1×M]=[x1,x2,…,x2M]T
Y1×2N=[realX1×N,imagX1×N]=[y1,y2,…,y2N]T
Figure GDA0003854009850000051
hj,i为信道矩阵H的内部元素,1≤i≤2M,1≤j≤2N;
rj为天线接收的信号Y=HX+σ2,的内部元素,1≤j≤2N;
σ2为信道噪声。
相关的计算步骤如下:
步骤一:系统的调制阶数为log2Q;s是发送符号值的集合,
Figure GDA0003854009850000052
假设第l次迭代中符号节点i到观察节点j传递的先验概率
Figure GDA0003854009850000053
包含
Figure GDA0003854009850000054
个概率值。初始化先验概率
Figure GDA0003854009850000055
则对应s中取各个元素的选取
Figure GDA0003854009850000056
1≤i≤2M,1≤j≤2N,
Figure GDA0003854009850000057
步骤二:计算观察节点j所收信息的均值
Figure GDA0003854009850000058
1≤j≤2N和方差
Figure GDA0003854009850000059
Figure GDA00038540098500000510
1≤j≤2N
步骤三:计算其他符号节点和信道带来的噪声干扰的均值μj,i=μj-hj,isTpi,j,1≤i≤2M,1≤j≤2N
和方差
Figure GDA0003854009850000061
1≤i≤2M,1≤j≤2N
建立高斯模型
Figure GDA0003854009850000062
1≤i≤2M,1≤j≤2N。
需要说明的是,上述均值计算的过程如下,以观察节点j为例,与j相关的符号节点有2M个,如图2所示,对于每个符号节点可能的符号值以及其对应的概率做乘法,送入累加器,得到的结果乘以该符号节点i和观察节点j之间的信道信息hj,i,进行累加,对于累加的结果减去各自对应的信息,即高斯模型所需的均值信息。
另外,上述方差计算的过程如下,以观察节点j为例,与j相关的符号节点有2M个,对于每个符号节点可能的符号值的平方以及其对应的概率做乘法,送入累加器。对于每个符号节点可能的符号值以及其对应的概率做乘法,送入累加器,再做平方,得到的结果乘以该符号节点i和观察节点j之间的信道hj,i,进行累加,对于累加的结果减去各自对应的信息再加上信道带来的噪声方差即高斯模型所需的方差信息。
步骤四:计算观察节点的LLR值;
具体的,由以下公式计算得出,
Figure GDA0003854009850000063
S30:根据观察节点的后验概率的LLR值得到先验概率的LLR值;
具体的,对于第i个符号节点,将2N个观察节点反馈回来的信息进行综合,得到
Figure GDA0003854009850000064
为第i号天线上符号的LLR信息,进一步得到先验概率的LLR值αi,j(s),
Figure GDA0003854009850000065
S40:根据所述先验概率的LLR值计算出符号节点要发给各个观察节点的先验概率值,并开始新一轮的先验概率的迭代;
具体的,分析用于计算的αi,j(s)的维度,在具体计算时,以单根天线为例分三步计算相应的概率值:
利用αi,j(s)计算
Figure GDA0003854009850000066
的先验概率值,首先从
Figure GDA0003854009850000067
选出最大值αmax,然后[αmaxi,j(sk)]乘以常数log2e,得到一个新的数A=u+v,v为小数部分数值,u为整数部分数值,将1-0.5v右移u位,得到
Figure GDA0003854009850000068
的值;
将得到的
Figure GDA0003854009850000071
的值进行累加得到F,根据F最高位相应的2的幂次,确定w的值,将F右移w位得到m的值,则根据ln(F)=ln2·(w+m-1)完成对数运算;
maxi,j(sk)]的值与ln(F)求和,再乘以常数log2e得到一个新的数值A′=u′+v′,v′为小数部分数值,u′为整数部分数值,将1-0.5v′右移u′位,即可得到相对应的先验概率值。
需要说明的是,上述F表示所有
Figure GDA0003854009850000072
的值累加得到的数值;w表示F最高位相应的2的幂次的数值;m表示将F右移w位得到的数值。
另外需要说明的是,在上述概率计算的过程中,指数运算采用以下方法:对于输入的值xi,乘以一个常数log2e,得到xilog2e,将其整数部分和小数部分分离,并将小数部分右移一位,1减去右移一位后的小数值再右移整数位即得到对应的指数运算的值。
另外,在上述概率计算的过程中,指数运算采用以下方法:对于输入的值xi,找到其转换成二进制表示的整数最高位1的位置w,将xi右移w位得到m,再与w-1相加,最后乘以常数ln2完成对数运算。
S50:当达到最大迭代次数时,根据需求输出符号。
具体的,在完成迭代后,输出最大后验信息γmax=max(γ)对应s中的符号即为检测结果。
由以上步骤可知,该MIMO BP检测方法在复杂度上低于最大似然检测。相较于MMSE检测算法,避免了该算法中涉及到的矩阵求逆运算,较大程度上减少了硬件资源损耗。对于MIMO BP检测算法中的概率更新计算公式,采用线性插值的方式进行近似计算,避免了e指数查表需要同时存储多张表以及地址读取冲突问题。
参照图2,本实施例还提供本发明提供一种基于置信度传播的MIMO检测装置,包括:
获取单元100,用于获取MIMO的系统参数;具体的获取过程在上述MIMO检测方法的步骤S10中已经详细阐述,在此不再赘述。
第一计算单元200,用于观察节点根据所述系统参数以及各个符号节点的先验概率计算各个接收符号的后验概率的LLR值,并将其传回给各个符号节点;具体的计算过程在上述MIMO检测方法的步骤S20中已经详细阐述,在此不再赘述。
需要说明的是,该第一计算模块中涉及到均值计算是利用均值计算模块实现的,如图3所示,以观察节点j为例,与j相关的符号节点有2M个,对于每个符号节点可能的符号值以及其对应的概率做乘法,送入累加器,得到的结果乘以该符号节点i和观察节点j之间的信道信息hj,i,进行累加,对于累加的结果减去各自对应的信息,即高斯模型所得的均值信息。
同样的,第一计算模块中涉及到方差计算是利用方差计算模块实现的,方差计算模块如图4所示,以观察节点j为例,与j相关的符号节点有2M个,对于每个符号节点可能的符号值的平方以及其对应的概率做乘法,送入累加器。对于每个符号节点可能的符号值以及其对应的概率做乘法,送入累加器,再做平方,得到的结果乘以该符号节点i和观察节点j之间的信道hj,i,进行累加,对于累加的结果减去各自对应的信息,即高斯模型所得的方差信息。
第二计算单元300,用于根据观察节点的后验概率的LLR值得到先验概率的LLR值;具体的计算过程在上述MIMO检测方法的步骤S40中已经详细阐述,在此不再赘述。
需要说明的是,该过程中涉及到的LLR值计算是利用LLR计算模块实现的,LLR计算模块如图5所示,观察节点接收到的信息rj减去均值uj,i,信道hj,i左移一位与rjj,i相乘,得到的结果2hj,i(rjj,i),再与(s-s1)做乘法;信道hj,i与自身做乘法,再与s2-s1 2做乘法,得到
Figure GDA0003854009850000081
两者做一次减法得到值
Figure GDA0003854009850000082
除以
Figure GDA0003854009850000083
得到观察节点的LLR值。
第三计算单元400,用于根据所述先验概率的LLR值计算出符号节点要发给各个观察节点的先验概率值,并开始新一轮的先验概率的迭代;具体的计算过程在上述MIMO检测方法的步骤S40中已经详细阐述,在此不再赘述。
具体的概率计算模块如图6所示:分析用于计算的αi,j(s)的维度,在具体计算时,以单根天线为例分三步计算相应的概率值:
利用αi,j(s)计算
Figure GDA0003854009850000084
的先验概率值,首先从
Figure GDA0003854009850000085
选出最大值αmax,然后[αmaxi,j(sk)]乘以常数log2e,得到一个新的数A=u+v,v为小数部分数值,u为整数部分数值,将1-0.5v右移u位,得到
Figure GDA0003854009850000086
的值;
将得到的
Figure GDA0003854009850000087
的值进行累加得到F,根据F最高位相应的2的幂次,确定w的值,将F右移w位得到m的值,则根据ln(F)=ln2·(w+m-1)完成对数运算;
maxi,j(sk)]的值与ln(F)求和,再乘以常数log2e得到一个新的数值A′=u′+v′,v′为小数部分数值,u′为整数部分数值,将1-0.5v′右移u′位,即可得到相对应的先验概率值。
上述概率计算的过程中涉及到的指数计算利用指数计算模块实现,具体的指数计算模块如图7所示:对于输入的值xi,乘以一个常数log2e,得到xilog2e,将其整数部分和小数部分分离,并将小数部分右移一位,1减去右移一位后的小数值再右移整数位即得到对应的指数运算的值。
上述概率计算的过程中涉及到的对数计算利用对数计算模块实现,对数计算模块如图8所示:对于输入的值xi,找到其转换成二进制表示的整数最高位1的位置w,将xi右移w位得到m,再与w-1相加,最后乘以常数ln2完成对数运算。
选取单元500,迭代完成后,根据需求输出符号。具体的选取过程在上述MIMO检测方法的步骤S50中已经详细阐述,在此不再赘述。
另外,本发明实施例还提供一种计算机可读存储介质,其中,该计算机可读存储介质可存储有程序,该程序执行时包括上述方法实施例中记载的任何MIMO检测方法的部分或全部步骤。
另外,在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储器中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储器中,包括若干指令用以使得一台计算机设备(可为个人计算机、服务器或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储器包括:U盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、移动硬盘、磁碟或者光盘等各种可以存储程序代码的介质。
本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成,该程序可以存储于一计算机可读存储器中,存储器可以包括:闪存盘、只读存储器(英文:Read-Only Memory,简称:ROM)、随机存取器(英文:Random Access Memory,简称:RAM)、磁盘或光盘等。
以上参照附图描述了根据本发明的实施例的用于实现服务链的方法的示例性流程图。应指出的是,以上描述中包括的大量细节仅是对本发明的示例性说明,而不是对本发明的限制。在本发明的其他实施例中,该方法可具有更多、更少或不同的步骤,且各步骤之间的顺序、包含、功能等关系可以与所描述和图示的不同。

Claims (9)

1.一种基于置信度传播的MIMO信号检测方法,其特征在于,包括:
获取MIMO的系统参数;所述系统参数包括发送天线数,接收天线数,系统的调制阶数,接收天线所接收到的发射天线经过信道传输和噪声干扰的接收信号、对应的信道信息以及噪声方差;
观察节点根据所述系统参数以及各个符号节点的先验概率计算各个接收符号的后验概率的LLR值,并将其传回给各个符号节点;
根据观察节点的后验概率的LLR值得到先验概率的LLR值;
根据所述先验概率的LLR值计算出符号节点要发给各个观察节点的先验概率值,并开始新一轮的先验概率的迭代;
当达到最大迭代次数时,根据需求输出符号;
其中所述根据所述先验概率的LLR值计算出符号节点要发给各个观察节点的先验概率值包括:
对于αi,j(s),从
Figure FDA0003943533460000011
选出最大值αmax,然后[αmaxi,j(sk)]乘以常数log2e,得到一个新的数A=u+v,v为小数部分数值,u为整数部分数值,将1-0.5v右移u位,得到
Figure FDA0003943533460000012
的值;
将得到的
Figure FDA0003943533460000013
的值进行累加得到F,根据F最高位相应的2的幂次,确定w的值,将F右移w位得到m的值,则根据ln(F)=ln2·(w+m-1)完成对数运算;
maxi,j(sk)]的值与ln(F)求和,再乘以常数log2e得到一个新的数值A′=u′+v′,v′为小数部分数值,u′为整数部分数值,将1-0.5v′右移u′位,即可得到相对应的先验概率值;其中,αi,j(s)表示第i个符号节点到第j个观察节点的先验信息的LLR值;
Figure FDA0003943533460000014
s是发送符号值的集合;系统调制阶数为log2Q;F表示所有
Figure FDA0003943533460000015
的值累加得到的数值;w表示F最高位相应的2的幂次的数值;m表示将F右移w位得到的数值。
2.根据权利要求1所述的MIMO信号检测方法,其特征在于,所述观察节点根据符号节点的先验概率计算各个接收符号的后验概率的LLR值,并将其传回给各个符号节点,具体包括:计算观察节点j所收信息的均值
Figure FDA0003943533460000016
和方差
Figure FDA0003943533460000017
Figure FDA0003943533460000018
计算其他符号节点和信道带来的噪声干扰的均值μj,i=μj-hj,isTpi,j,1≤i≤2M,1≤j≤2N和方差
Figure FDA0003943533460000021
计算观察节点的后验概率的LLR值,
Figure FDA0003943533460000022
式中,M为发送天线数;N为接收天线数;hj,i为信道矩阵H的内部元素;i表示第i个符号节点;j表示第j个观察节点;
Figure FDA0003943533460000023
s是发送符号值的集合;系统调制阶数为log2Q;(*)T表示对矩阵做转置;pi,j为第l次迭代中符号节点i到观察节点j传递的先验概率;σ2为信道噪声;rj为天线接收的信号Y=HX+σ2的内部元素。
3.根据权利要求1所述的MIMO信号检测方法,其特征在于,计算过程中的指数运算采用以下方法:
对于输入的值xi,乘以一个常数log2e,得到xilog2e,将其整数部分和小数部分分离,并将小数部分右移一位,1减去右移一位后的小数值再右移整数位即得到对应的指数运算的值。
4.根据权利要求1所述的MIMO信号检测方法,其特征在于,计算过程中的对数运算采用以下方法:
对于输入的值xi,找到其转换成二进制表示的整数最高位1的位置w,将xi右移w位得到m,再与w-1相加,最后乘以常数ln2完成对数运算。
5.一种基于置信度传播的MIMO信号检测装置,其特征在于,包括:
获取单元,用于获取MIMO的系统参数;所述系统参数包括发送天线数,接收天线数,系统的调制阶数,接收天线所接收到的发射天线经过信道传输和噪声干扰的接收信号、对应的信道信息以及噪声方差;
第一计算单元,用于观察节点根据所述系统参数以及各个符号节点的先验概率计算各个接收符号的后验概率的LLR值,并将其传回给各个符号节点;
第二计算单元,用于根据观察节点的后验概率的LLR值得到先验概率的LLR值;
第三计算单元,用于根据所述先验概率的LLR值计算出符号节点要发给各个观察节点的先验概率值,并开始新一轮的先验概率的迭代;
选取单元,用于当达到最大迭代次数时,根据需求输出符号;
其中所述根据所述先验概率的LLR值计算出符号节点要发给各个观察节点的先验概率值包括:
对于αi,j(s),从
Figure FDA0003943533460000024
选出最大值αmax,然后[αmaxi,j(sk)]乘以常数log2e,得到一个新的数A=u+v,v为小数部分数值,u为整数部分数值,将1-0.5v右移u位,得到
Figure FDA0003943533460000031
的值;
将得到的
Figure FDA0003943533460000032
的值进行累加得到F,根据F最高位相应的2的幂次,确定w的值,将F右移w位得到m的值,则根据ln(F)=ln2·(w+m-1)完成对数运算;
maxi,j(sk)]的值与ln(F)求和,再乘以常数log2e得到一个新的数值A′=u′+v′,v′为小数部分数值,u′为整数部分数值,将1-0.5v′右移u′位,即可得到相对应的先验概率值;其中,αi,j(s)表示第i个符号节点到第j个观察节点的先验信息的LLR值;
Figure FDA0003943533460000033
s是发送符号值的集合;系统调制阶数为log2Q;F表示所有
Figure FDA0003943533460000034
的值累加得到的数值;w表示F最高位相应的2的幂次的数值;m表示将F右移w位得到的数值。
6.根据权利要求5所述的MIMO信号检测装置,其特征在于,所述第一计算单元包括LLR计算单元,所述LLR计算单元用于计算观察节点的LLR值,具体包括:
观察节点接收到的信息rj减去均值uj,i,信道hj,i左移一位与rjj,i相乘,得到的结果2hj,i(rjj,i),再与(s-s1)做乘法;信道hj,i与自身做乘法,再与s2-s1 2做乘法,得到
Figure FDA0003943533460000035
两者做一次减法得到值
Figure FDA0003943533460000036
除以
Figure FDA0003943533460000037
得到观察节点的LLR值;
其中,rj为天线接收的信号Y=HX+σ2的内部元素,hj,i为信道矩阵H的内部元素;i表示第i个符号节点;j表示第j个观察节点,σ2为信道噪声。
7.根据权利要求5所述的MIMO信号检测装置,其特征在于,所述第三计算单元包括指数计算单元,所述指数计算单元用于对整个计算过程中的指数进行运算,具体的运算步骤如下:对于输入的值xi,乘以一个常数log2e,得到xilog2e将其整数部分和小数部分分离,并将小数部分右移一位,1减去右移一位后的小数值再右移整数位即得到对应的指数运算的值。
8.根据权利要求5所述的MIMO信号检测装置,其特征在于,所述第三计算单元包括对数计算单元,所述对数计算单元用于对整个计算过程中的对数进行运算,具体的运算步骤如下:对于输入的值xi,找到其转换成二进制表示的整数最高位1的位置w,将xi右移w位得到m,再与w-1相加,最后乘以常数ln2完成对数运算。
9.一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如权利要求1至4任一项所述的一种基于置信度传播的MIMO信号检测方法的步骤。
CN202110148962.1A 2021-02-03 2021-02-03 一种基于置信度传播的mimo检测方法、装置 Active CN112929128B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110148962.1A CN112929128B (zh) 2021-02-03 2021-02-03 一种基于置信度传播的mimo检测方法、装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110148962.1A CN112929128B (zh) 2021-02-03 2021-02-03 一种基于置信度传播的mimo检测方法、装置

Publications (2)

Publication Number Publication Date
CN112929128A CN112929128A (zh) 2021-06-08
CN112929128B true CN112929128B (zh) 2023-02-28

Family

ID=76169674

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110148962.1A Active CN112929128B (zh) 2021-02-03 2021-02-03 一种基于置信度传播的mimo检测方法、装置

Country Status (1)

Country Link
CN (1) CN112929128B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114301545B (zh) * 2021-12-03 2024-04-19 网络通信与安全紫金山实验室 信号检测方法、装置、电子设备及存储介质
CN115378524A (zh) * 2022-08-09 2022-11-22 网络通信与安全紫金山实验室 一种优化的期望传播检测方法及信号检测装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103781074A (zh) * 2014-02-08 2014-05-07 北京极科极客科技有限公司 无线局域网安全接入方法及设备
WO2016020611A1 (fr) * 2014-08-05 2016-02-11 Orange Procede et dispositif de reception de signaux
CN106209189A (zh) * 2016-07-19 2016-12-07 清华大学 大规模mimo系统中信号检测装置及方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107241167B (zh) * 2017-06-29 2020-02-18 东南大学 一种基于大规模mimo系统bp检测的改进方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103781074A (zh) * 2014-02-08 2014-05-07 北京极科极客科技有限公司 无线局域网安全接入方法及设备
WO2016020611A1 (fr) * 2014-08-05 2016-02-11 Orange Procede et dispositif de reception de signaux
CN106209189A (zh) * 2016-07-19 2016-12-07 清华大学 大规模mimo系统中信号检测装置及方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
网络空间安全体系与关键技术;罗军舟等;《中国科学:信息科学》;20160820(第08期);全文 *

Also Published As

Publication number Publication date
CN112929128A (zh) 2021-06-08

Similar Documents

Publication Publication Date Title
CN112929128B (zh) 一种基于置信度传播的mimo检测方法、装置
US7457367B2 (en) Detector and method for estimating data probability in a multi-channel receiver
JP6009717B2 (ja) 低複雑度受信機および低密度シグネチャ変調のための方法
CN111224905B (zh) 一种大规模物联网中基于卷积残差网络的多用户检测方法
JP2009527182A (ja) データ信号を処理する方法、データ処理部、及びコンピュータプログラム製品
JP5723371B2 (ja) Ofdm通信システムにおいて通信路を推定するためのプロセス、およびそのための受信機
KR20230133306A (ko) 셀프리 mimo의 다중 사용자 검출 방법 및 장치
CN115037340B (zh) 信号检测方法、装置、电子设备及存储介质
CN113630160B (zh) 一种大规模mimo检测方法、装置、设备和存储介质
CN113055318B (zh) 一种信道估计方法
CN113364535B (zh) 数学形式多输入多输出检测方法、系统、设备及存储介质
WO2023237182A1 (en) Radio receiver with multi-stage equalization
CN100369390C (zh) 多天线阵列接收无线传输的方法
JP6102250B2 (ja) 有効チャネルおよびフィードバック情報を決定する装置および方法
EP3859609A1 (en) Finite resolution decomposition of a matrix for low-complexity and energy-efficient matrix-vector multiplication
Niazadeh et al. An alternating minimization method for sparse channel estimation
CN108075866B (zh) 基于重复传输系统的解调方法和装置
CN115706684B (zh) 信道参数的估计方法、装置及存储介质
CN112490688B (zh) 一种确定方法、设备及计算机存储介质
JP6259214B2 (ja) Mimo受信装置
CN112583748B (zh) 一种信道估计方法、装置及电子设备
US20230179315A1 (en) Method for Disseminating Scaling Information and Application Thereof in VLSI Implementation of Fixed-Point FFT
WO2022222896A1 (zh) 一种发送、接收方法、装置、电子设备和存储介质
KR101343913B1 (ko) 공간 다중화를 지원하는 다중 입출력 시스템에서 병렬 심볼제거 방법 및 장치
EP3378168B1 (en) Apparatus and method for deriving a submatrix

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