CN102006148B - 基于宽度优先树形搜索的mimo信号检测方法 - Google Patents
基于宽度优先树形搜索的mimo信号检测方法 Download PDFInfo
- Publication number
- CN102006148B CN102006148B CN 201010577316 CN201010577316A CN102006148B CN 102006148 B CN102006148 B CN 102006148B CN 201010577316 CN201010577316 CN 201010577316 CN 201010577316 A CN201010577316 A CN 201010577316A CN 102006148 B CN102006148 B CN 102006148B
- Authority
- CN
- China
- Prior art keywords
- path
- layer
- real number
- signal
- subgroup
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Radio Transmission System (AREA)
Abstract
本发明公开了一种基于“宽度优先”树形搜索的MIMO信号检测方法,主要解决现有MIMO检测方法复杂度高的问题。其实现步骤为(1)***预处理,包括信道矩阵QR分解,预计算路径度量因子;(2)利用“宽度优先”树搜索思想逐层检测信号:2a)利用上一层幸存路径扩展本层路径:利用Schnorr-Euchner枚举法依次确定路径扩展顺序、扩展路径和路径度量;2b)利用归并排序法对扩展路径的路径度量排序以确定幸存路径;(3)若所有层信号检测完成,将最小路径度量对应的幸存路径矢量作为最终检测输出;否则将幸存路径传给下一层检测下一层信号。本发明具有复杂度低,性能损失小的优点,可用于下一代宽带无线通信***MIMO接收机的信号检测。
Description
技术领域
本发明属于无线通信***领域,涉及一种MIMO空间复用***的信号检测方法,适合于下一代宽带无线通信***中的MIMO接收机。
背景技术
MIMO技术充分利用空间的多径特性可以达到扩充信道容量或提高信号可靠度的作用,成为下一代无线通讯的关键技术之一,但是在接收端如何可靠的检测多天线信号成为MIMO技术实现的关键之一,最优的最大似然ML检测算法由于其随天线数及调制阶数指数级增长的复杂度而难以在实际中应用,因此如何设计低复杂度,性能近似ML检测的MIMO检测算法,是一个很有实际意义的工作。
传统的MIMO检测算法主要包括线性与非线性检测部分,线性有迫零ZF、最大均方误差MMSE检测等,非线性包括并行干扰抵消、串行干扰抵消检测等,这些算法复杂度虽低但性能较ML相差较大,因此通常用在一些对误码率要求不高的场合。
近年来人们研究发现MIMO***具有树结构特性,因此可以用树形搜索方法来进行检测。树形搜索算法由于低复杂度却近似ML检测性能的特性得到了广泛的应用,Garret在专利(07782984)”Method Of SphereDecoding With Low Complexity And Good Statistical Output”中提出利用基于“深度优先”的球形检测SD算法,但该类算法的复杂度随信道条件变 化,不利于硬件实现;Reuven在专利(7720169)“Multiple-inputmultiple-output(MIMO)detector incorporating efficient signal point searchand soft information refinement”中,提出了针对2发2收MIMO***的基于“宽度优先”搜索策略的K-best检测方法,但该方法亦不太适合高速并行处理。
基于“宽度优先”搜索策略的K-Best检测算法由于路径度量计算及排序等过程,还是具有较高的复杂度,不方便硬件实现,因此需要进一步降低检测的复杂度,以便于实际应用。
发明内容
本发明目的在于克服上述已有技术的不足,提出一种基于“宽度优先”树形搜索的MIMO信号检测方法,在保证其吞吐率的前提下,降低传统K-Best检测方法的复杂度,以满足未来宽带无线通信***的要求。
本发明是这样实现的:
一.技术原理:
定义:具有NT个发送天线、NR个接收天线的MIMO***信号模型为:
r=Hs+n
其中, 是NR维接收列向量,ri为第i根接收天线的接收信号,上标T表示转置运算; 为NT维发送信号列向量,si为第i根发射天线的发送信号; 为NR维噪声列向量,ni表示第i根接收天线的信号噪声; 是NR×NT维信道传输矩阵,hij表示发送天线j到接收天线i的信道衰落系数,i=1,2,..NR, j=1,2,..NT,NT表示发射天线数,NR表示接收天线数。
本发明涉及的MIMO信号检测即根据所述的接收信号r和信道矩阵H估计发射信号s。
二.技术方案
本发明基于“宽度优先”树形搜索的MIMO信号检测方法,包括:
(1)***预处理步骤:
1a)对接收信号r及信道矩阵H进行实数化处理,得到实数接收向量rRe与实数化矩阵HRe;
1b)对实数化矩阵HRe进行QR分解,确定正交阵Q和上三角阵R,并利用Q阵与实数接收向量rRe计算变换的接收信号;r′=QT·rRe,上标T表示转置运算;
1c)将实数星座点Ω与上三角阵R相乘,得到预计算度量因子阵Γ;
1d)令当前层索引i=2NT,初始化根节点,即将2NT+1层的幸存路径为空,将路径度量设为0,确定每层检测要保留的幸存路径数K=Mc,Mc表示实数星座点的大小;
(2)从第2NT层到第1层逐层检测信号步骤:
2a)利用第i+1层保留的幸存路径扩展第i层节点:扩展前先利用Schnorr-Euchner枚举方法确定路径符号扩展的顺序,再以Schnorr-Euchner枚举结果进行路径扩展;当i=2NT时,扩展出Mc条新路径;其它层每条路径扩展出Mc个分支,共扩展得到K×Mc条新路径,确定每条新路径对应的路径度量值;
2b)利用归并排序方法对步骤2a)中扩展得到的路径按路径度量值的 大小进行排序,确定K个最小度量值,并将这些最小度量值对应的路径放入路径寄存器中;
(3)更新层索引步骤:
检验当前层索引,如果当前层索引i=1,执行步骤(4),否则将路径寄存器中存放的幸存路径及对应的路径度量值输出到第i-1层,更新当前层索引将i减1,返回步骤(2),对第i-1层的信号进行检测,直到当前层索引i=1为止;
(4)输出检测信号步骤:
输出路径寄存器中最小度量值对应的路径矢量作为输出,检测完毕。
本发明与现有技术相比具有如下优点:
1)本发明由于在预处理步骤预先计算星座点集Ω与上三角阵R的乘积即路径度量因子,并利用星座点集Ω的特性简化了计算过程,从而避免了在路径扩展时对路径度量因子的重复计算,有效的节省了路径度量计算的复杂度和处理时间,降低了逐层检测信号的复杂度;
2)本发明由于在路径扩展及路径度量计算前,利用Schnorr-Euchner枚举的方法使得每条路径扩展的分支按升序排列,从而简化了后续的排序过程;
3)本发明由于用归并排序的方法代替传统的冒泡排序方法,能有效降低排序的复杂度,如对于K=8的64QAM调制信号的检测,采用冒泡排序需要执行约64+63+62+…+57=484次比较及交换操作,如果采用本发明中的归并排序方法,仅需执行约36*(4+2+1)=252次比较及交换操作,该方法还采用并行流水结构实现,与冒泡排序相比节省了处理时间。
仿真结果表明,本发明降低复杂度所造成的检测性能损失很小,在保证吞吐率的情况下,能满足下一代宽带无线通信***高速处理和误比特率性能的要求。
附图说明
图1是本发明的检测流程图;
图2是本发明归并排序的子流程图;
图3是对图2中子组两两合并排序的子流程图;
图4是本发明与传统K-Best检测算法的性能仿真比较图。
具体实施方式
下面将结合附图对本发明的具体实现过程进行详细描述。
参照附图1,本发明检测方法的具体实现步骤如下:
步骤一:***预处理步骤:
1a)对接收到的NR维接收信号r及NR×NT维的信道矩阵H进行实数化处理,得到2NR维的实数接收向量rRe与2NR×2NT维的实数矩阵HRe,并确定对应的实数星座点集Ω,对于QPSK调制信号Ω={1,-1};16QAM调制信号Ω={1,3,-1,-3};64QAM调制信号Ω={1,3,5,7,-1,-3,-5,-7};
所述的实数化处理按如下公式进行:
这里,Re(g)表示取实部运算,Im(g)表示取虚部运算;
实数矩阵HRe及实数接收向量rRe分别为:
1b)对实数矩阵HRe进行QR分解,计算2NT×2NT维的上三角阵R及2NT维的变换接收信号r′, ,上标T表示转置运算;采用Givens变换实现QR分解,对QR分解模块输入[HRe|rRe],经过QR分解,输出[R|r′]
其中rij表示上三角阵R第i行第j列元素,i=1,2,..2NT,j=1,2,..2NT,NT表示发送天线数;
1c)执行预计算,即将实数星座点集Ω与上三角阵R相乘,得到路径度量因子阵Γ=Ω·R,这里根据Ω中元素正负对称的特性,只计算Ω正数部分点与R的乘积,如64QAM调制,实数星座点集为Ω={1,3,5,7,-1,-3,-5,-7},只计算子集{1,3,5,7}与三角阵R的乘积,以此简化计算结构,得到2NT×2NT×Mc维的度量因子阵Γ;其中Mc表示实数星座点集Ω的大小,NT表示发送天线数;
1d)令当前层索引i=2NT,初始化根节点,即将2NT+1层的幸存路径为空,将路径度量设为0,确定每层检测要保留的幸存路径数K=Mc,Mc表示实数星座点的大小。
步骤二:从第2NT层到第1层逐层检测信号。
2a):利用i+1层保留的幸存路径扩展第i层的符号节点:
2a1)利用Schnorr-Euchner枚举方法确定路径的符号扩展顺序,
Schnorr-Euchner枚举方法是E.Agrell等在文献“Closet Point Searchin Lattice”中提出的一种快速搜索最近符号点的方法,本发明中采用的Schnorr-Euchner枚举步骤如下:
首先,计算第i层第j条路径符号估计点dj=[r′i+1,j]i/Ri,i,r′i+1,j表示第i+1层第j条幸存路径对应的接收信号矢量,Ri,i表示上三角R阵第i行第i列值,i为当前层索引,j=1,2,..K;
接着,利用符号估计点dj确定符号扩展顺序,
对于64QAM信号Ω={1,3,5,7,-1,-3,-5,-7},第i层第j条路径符号扩展顺序为下面Φ矩阵的第m行所给出的符号顺序;
对于16QAM信号,Ω={1,3,-1,-3},第i层第j条路径符号扩展顺序为下面Φ矩阵的第m行所给出的符号顺序;
对于QPSK信号,当dj>0时,符号扩展顺序为{1,-1},当dj<0时,符号扩展顺序为{-1,1};这里Φ为符号扩展顺序矩阵, 表示根据dj取值确定在符号扩展顺序矩阵Φ中的行索引, 表示取不大于dj的最大整数;按照Schnorr-Euchner枚举确定的符号扩展顺序调整Γi中每一行顺序,得到新的度量因子向量Γ′i,这样在计算路径度量时,每条路径扩展出的分支按路径度量值升序排列,便于后续的排序操作,其中Γi表示度量因子阵Γ第i列;
2a2)计算扩展的路径矢量对应的路径度量值Pi,j,l。
路径度量Pi,j,l的计算通过如下公式进行:
Pi,j,l=Pi+1,j+([r′i+1,j]i-[Γ′i,i]l)2,j∈1,2,..K,l∈1,2,..Mc 5)
上式中,Pi,j,l表示第i层第j条幸存路径第1条分支对应的路径度量,r′i,j,l表示第i层第j条幸存路径第1条分支在去除2NT:(2NT-i)层干扰后的剩余接收向量,Γi表示度量因子阵Γ的第i列,Γ′i表示Γi由2.1)调整顺序后得到的向量,Γ′i,i表示Γ′i的第i行元素,K表示上一层检测得到的幸存路径数,Mc表示实数星座点的大小,[·]n表示矢量的第n个分量, 表示矢量的前n个分量,i=2NT....1;
通过式5)-6)的计算公式,第i层符号扩展时,输入K个幸存路径矢量,输出K个新的路径矢量分组,每个分组有Mc条路径,且该分组中 路径以路径度量Pi,j,l大小按升序排列,即Pi,j,1<Pi,j,2<L<Pi,j,Mc,j=1,2,LK;
2b)利用归并排序对步骤2a)扩展得到的K×Mc条新路径进行排序,得到K个最小路径度量值节点及其所对应的路径矢量。
参照图2,本步骤的归并排序如下:
2b1)将每条路径扩展的分支作为一个子组,则共有K个子组,每个子组中度量值按升序排列;
2b2)将这些升序排列的子组两两合并排序:
首先,设定一个大小为K的最小值列表,定义该最小值列表的长度L=K;
接着,将两个子组奇偶间隔合并为一个队列,定义该队列长度N=2n,n表示子组的长度;
接着,将队列相邻元素进行两两比较,并将队头元素输出到最小值列表的第K-L+1位置处,更新最小值列表长度,将最小值列表长度L减1;删除队头元素,更新队列长度,将队列长度N减1;比较更新后的队列长度N与更新后的最小值列表长度L:如果N>L,则将队尾元素删除,并再次更新队列长度将队列长度N减1;否则不删除队尾元素;对队列更新K次,将最小值列表填满,再将最小值列表中的元素作为本次合并排序的输出,如图3所示;
2b3)重复执行步骤2b2),直到当前只有一个子组为止,再将该子组作为本层的幸存路径输出。
归并排序完成后,从K×Mc条扩展路径中得到K个最小路径度量值 Pi,0,Pi,1,...Pi,K-1,i为当前层索引,将这K个最小路径度量值对应的路径矢量作为第i层的幸存路径存入路径寄存器中。
步骤三:更新层索引:
检验当前层索引,如果当前层索引i=1,执行步骤四;否则将路径寄存器中存放的幸存路径及对应的路径度量值输出到第i-1层,更新当前层索引将i减1,返回步骤二,对下一层信号进行检测,直到当前层索引i=1为止;
步骤四:将路径寄存器中最小路径度量值对应的路径矢量作为输出,检测完毕。
本发明的效果可以通过仿真说明:
1)仿真条件:***使用2发2收MIMO***,信道采用瑞利块衰落信道,调制方式选择为64-QAM。
2)仿真内容与结果:
用本发明的K-Best方法与传统的K-Best方法在上述仿真条件下针对误比特率随信噪比变化的性能进行仿真比较,仿真结果如图4所示。
由图4可见,对64QAM调制方式,本发明与传统的K-Best方法相比,其误比特率性能损失可以忽略,且本发明的实现复杂度仅为传统PD方法的60%左右。
Claims (4)
1.一种基于“宽度优先”树形搜索的MIMO信号检测方法,包括:
(1)***预处理步骤:
1a)对接收信号r及信道矩阵H进行实数化处理,得到实数接收向量rRe与实数化矩阵HRe;
1b)对实数化矩阵HRe进行QR分解,确定正交阵Q和上三角阵R,并利用Q阵与实数接收向量rRe计算变换的接收信号;r′=QT·rRe,上标T表示转置运算;
1c)将实数星座点Ω与上三角阵R相乘,得到预计算度量因子阵Γ;
1d)令当前层索引i=2NT,初始化根节点,即将2NT+1层的幸存路径为空,将路径度量设为0,确定每层检测要保留的幸存路径数K=Mc,Mc表示实数星座点的大小,NT表示发射天线数;
(2)从第2NT层到第1层逐层检测信号步骤:
2a)利用第i+1层保留的幸存路径扩展第i层节点:扩展前先利用Schnorr-Euchner枚举方法确定路径符号扩展的顺序,再以Schnorr-Euchner枚举结果进行路径扩展;当i=2NT时,扩展出Mc条新路径;其它层每条路径扩展出Mc个分支,共扩展得到K×Mc条新路径,确定每条新路径对应的路径度量值;
2b)利用归并排序方法对步骤2a)中扩展得到的路径按路径度量值的大小进行排序,确定K个最小度量值,并将这些最小度量值对应的路径放入路径寄存器中;
(3)更新层索引步骤:
检验当前层索引,如果当前层索引i=1,执行步骤(4),否则将路径寄存器中存放的幸存路径及对应的路径度量值输出到第i-1层,更新当前层索引将i减1,返回步骤(2),对第i-1层的信号进行检测,直到当前层索引i=1为止;
(4)输出检测信号步骤:
将路径寄存器中最小度量值对应的路径矢量作为输出,完成检测。
2.根据权利要求1所述的检测方法,其中步骤2a)所述的以Schnorr-Euchner枚举结果进行路径扩展,是利用Schnorr-Euchner枚举确定的符号排列顺序调整度量因子向量Γi中每一行中的元素顺序,得到新的度量因子向量Γ′i,Γi为度量因子阵Γ的第i列。
4.根据权利要求1所述的检测方法,步骤2b)所述的利用归并排序方法进行排序,按照如下步骤进行:
2b1)将每条路径扩展的分支作为一个子组,则共有K个子组,每个子组中度量值按升序排列;
2b2)将这些升序排列的子组两两合并排序:
首先,设定一个大小为K的最小值列表,定义该最小值列表的长度L=K;
接着,将两个子组奇偶间隔合并为一个队列,定义该队列长度N=2n,n表示子组的长度;
接着,将队列相邻元素进行两两比较,并将队头元素输出到最小值列表的第K-L+1位置处,更新最小值列表长度,将最小值列表长度L减1;删除队头元素,更新队列长度,将队列长度N减1;比较更新后的队列长度N与更新后的最小值列表长度L:如果N>L,则将队尾元素删除,并再次更新队列长度将队列长度N减1;否则不删除队尾元素;对队列更新K次,将最小值列表填满,再将最小值列表中的元素作为本次合并排序的输出;
2b3)重复执行步骤2b2),直到当前只有一个子组为止,再将该子组作为本层的幸存路径输出。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010577316 CN102006148B (zh) | 2010-12-07 | 2010-12-07 | 基于宽度优先树形搜索的mimo信号检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010577316 CN102006148B (zh) | 2010-12-07 | 2010-12-07 | 基于宽度优先树形搜索的mimo信号检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102006148A CN102006148A (zh) | 2011-04-06 |
CN102006148B true CN102006148B (zh) | 2013-04-17 |
Family
ID=43813245
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010577316 Expired - Fee Related CN102006148B (zh) | 2010-12-07 | 2010-12-07 | 基于宽度优先树形搜索的mimo信号检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102006148B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102231641B (zh) * | 2011-07-21 | 2013-08-14 | 西安电子科技大学 | 多输入多输出mimo逐级并行检测方法 |
CN102724160B (zh) * | 2012-05-04 | 2015-03-04 | 西安电子科技大学 | 高阶调制多入多出***中的信号检测方法 |
CN104901910A (zh) * | 2014-03-07 | 2015-09-09 | 电信科学技术研究院 | 一种mimo***的检测方法及装置 |
CN105281814B (zh) * | 2014-07-23 | 2021-03-26 | 锐迪科(重庆)微电子科技有限公司 | 基本奇偶归并网格单元和幸存路径选取构建方法及装置 |
CN104168073B (zh) * | 2014-08-19 | 2016-09-28 | 大唐移动通信设备有限公司 | 一种信号检测方法及装置 |
CN104796239B (zh) * | 2015-01-30 | 2019-03-19 | 苏州恩巨网络有限公司 | 一种mimo无线通信***及信号检测装置和方法 |
CN104932864A (zh) * | 2015-06-25 | 2015-09-23 | 许继电气股份有限公司 | 基于流水线进程的归并排序方法及使用该方法的阀控装置 |
WO2018130292A1 (en) | 2017-01-13 | 2018-07-19 | Huawei Technologies Co., Ltd. | An optimized architecture for a signal decoder |
CN110109059B (zh) * | 2019-03-27 | 2023-04-07 | 西安电子科技大学 | 一种基于深度学习网络的雷达辐射源信号识别方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101431358A (zh) * | 2008-12-19 | 2009-05-13 | 西安电子科技大学 | 基于m-精英进化算法的垂直分层空时信号检测方法 |
CN101582750A (zh) * | 2009-06-02 | 2009-11-18 | 北京天碁科技有限公司 | 一种基于宽度优先搜索的球形译码检测方法 |
-
2010
- 2010-12-07 CN CN 201010577316 patent/CN102006148B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101431358A (zh) * | 2008-12-19 | 2009-05-13 | 西安电子科技大学 | 基于m-精英进化算法的垂直分层空时信号检测方法 |
CN101582750A (zh) * | 2009-06-02 | 2009-11-18 | 北京天碁科技有限公司 | 一种基于宽度优先搜索的球形译码检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102006148A (zh) | 2011-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102006148B (zh) | 基于宽度优先树形搜索的mimo信号检测方法 | |
CN1808959B (zh) | 一种传输数据的方法及通信*** | |
CN101674160B (zh) | 多输入多输出无线通信***信号检测方法及装置 | |
US7660363B2 (en) | Minimum error rate lattice space time codes for wireless communication | |
US8018828B2 (en) | Near ML decoding method based on metric-first search and branch length threshold | |
CN104243069B (zh) | 一种多天线交通通信网络***及信号检测方法 | |
CN101499840B (zh) | 多入多出***的迭代检测方法 | |
CN101662342B (zh) | 一种多输入多输出信号检测方法和装置 | |
CN101427485A (zh) | 复杂性降低的波束受控mimo ofdm*** | |
CN102723975B (zh) | Mimo***的信号检测方法及装置 | |
CN103685090A (zh) | 多输入多输出信道表现的预测方法 | |
CN101383652A (zh) | 一种多输入多输出***的信号检测方法及装置 | |
CN101442390A (zh) | 空间相关MIMO的Turbo均衡接收方法与装置 | |
CN108964725B (zh) | 时变大规模mimo网络中信道参数的稀疏估计方法 | |
JP4854378B2 (ja) | 無線伝送システムおよび無線伝送方法 | |
CN103188703A (zh) | 幸存星座点选择方法和qrm-mld信号检测方法 | |
CN110086743A (zh) | 一种基于差分编码的短突发mimo-ofdm通信***及方法 | |
CN103532889A (zh) | 一种软输出并行堆栈mimo信号检测方法 | |
CN101997652A (zh) | 基于ldpc-mimo通信***的接收检测方法和装置 | |
CN109981151A (zh) | 大规模mimo***中改进的高斯树近似消息传递检测算法 | |
CN102227098A (zh) | 一种多模mimo-scfde自适应传输***频域承载点选取方法 | |
CN103746728A (zh) | 一种混合自适应的mimo接收检测方法 | |
CN102231641B (zh) | 多输入多输出mimo逐级并行检测方法 | |
CN102487309A (zh) | 一种mimo***下的信号检测方法和装置 | |
CN101964667B (zh) | 用于长期演进方案的高效多天线检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130417 Termination date: 20181207 |
|
CF01 | Termination of patent right due to non-payment of annual fee |