CN110011946A - 支持快速译码的增强型可获得发射分集的正交空间调制方法 - Google Patents
支持快速译码的增强型可获得发射分集的正交空间调制方法 Download PDFInfo
- Publication number
- CN110011946A CN110011946A CN201910238613.1A CN201910238613A CN110011946A CN 110011946 A CN110011946 A CN 110011946A CN 201910238613 A CN201910238613 A CN 201910238613A CN 110011946 A CN110011946 A CN 110011946A
- Authority
- CN
- China
- Prior art keywords
- matrix
- block
- disperse
- symbol
- indicate
- 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.)
- Granted
Links
- 238000013459 approach Methods 0.000 title claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims abstract description 133
- 230000004913 activation Effects 0.000 claims abstract description 11
- 238000000034 method Methods 0.000 claims description 65
- 238000005259 measurement Methods 0.000 claims description 18
- 238000001514 detection method Methods 0.000 claims description 12
- 230000000717 retained effect Effects 0.000 claims description 12
- 238000007476 Maximum Likelihood Methods 0.000 claims description 11
- 238000013461 design Methods 0.000 claims description 10
- 230000009471 action Effects 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 7
- 238000012217 deletion Methods 0.000 claims description 5
- 230000037430 deletion Effects 0.000 claims description 5
- 101000859864 Rattus norvegicus Gamma-crystallin E Proteins 0.000 claims description 4
- 238000013507 mapping Methods 0.000 claims description 4
- 229910052705 radium Inorganic materials 0.000 claims description 4
- HCWPIIXVSYCSAN-UHFFFAOYSA-N radium atom Chemical compound [Ra] HCWPIIXVSYCSAN-UHFFFAOYSA-N 0.000 claims description 4
- 238000009825 accumulation Methods 0.000 claims description 3
- 239000000654 additive Substances 0.000 claims description 3
- 230000000996 additive effect Effects 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 238000009826 distribution Methods 0.000 claims description 3
- 238000007689 inspection Methods 0.000 claims description 2
- 230000017105 transposition Effects 0.000 claims description 2
- 230000009466 transformation Effects 0.000 claims 1
- 239000006185 dispersion Substances 0.000 abstract description 3
- 238000004088 simulation Methods 0.000 abstract 1
- 238000001228 spectrum Methods 0.000 description 23
- 230000005540 biological transmission Effects 0.000 description 13
- 230000008901 benefit Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 238000005457 optimization Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000010276 construction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 230000003213 activating effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000005562 fading Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/0413—MIMO systems
- H04B7/0456—Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L1/00—Arrangements for detecting or preventing errors in the information received
- H04L1/004—Arrangements for detecting or preventing errors in the information received by using forward error control
- H04L1/0045—Arrangements at the receiver end
- H04L1/0054—Maximum-likelihood or sequential decoding, e.g. Viterbi, Fano, ZJ algorithms
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L1/00—Arrangements for detecting or preventing errors in the information received
- H04L1/02—Arrangements for detecting or preventing errors in the information received by diversity reception
- H04L1/06—Arrangements for detecting or preventing errors in the information received by diversity reception using space diversity
- H04L1/0618—Space-time coding
- H04L1/0637—Properties of the code
- H04L1/0668—Orthogonal systems, e.g. using Alamouti codes
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
- H04L25/03006—Arrangements for removing intersymbol interference
- H04L25/03178—Arrangements involving sequence estimation techniques
- H04L25/03203—Trellis search techniques
- H04L25/03242—Methods involving sphere decoding
Abstract
本发明公开了一种支持快速译码的增强型可获得发射分集的正交空间调制方法,在每组集合的Q个弥散矩阵中每次激活P个弥散矩阵分别用来发送P个符号的实部和虚部。关于EDA‑QSM方案的解码,在对所有的2Q个分散矩阵进行合理的排序后,EDA‑QSM具有分块正交特性,这使得本发明支持快速译码。利用这种分块正交特性,本发明由于SCMB‑SD算法仅需执行一次球形搜索就可以同时估计出两组激活的序号向量和P个符号,因而具有很低的解码复杂度。仿真结果表明:在不同的系统配置下,EDA‑QSM较现有的其他空间调制方案具有更好的误码性能。
Description
【技术领域】
本发明属于多天线无线通信系统中的发射分集传输技术领域,涉及一种支持快速译码的增强型可获得发射分集的正交空间调制方法。
【背景技术】
空间调制(Spatial Modulation,SM)和空移键控调制(Space Shift Keying,SSK)利用空间维度来传输信息,较传统的MIMO方案可以获得额外的频谱效率。因此近年来,SM和SSK作为一种新颖的MIMO传输技术受到了广泛的关注。然而,SM和SSK在每次传输时只激活一根天线,所以它们不能获得发射分集,只能依赖接收分集来对抗信道衰落。
针对SM和SSK无法获得发射分集的缺陷,学者们提出了多种解决方法。例如,文献将SM的概念推广到空间和时间的维度,进而提出了可获得发射分集的空时移键控调制(Space-Time Shift Keying,STSK)方法。然而STSK的传输速率随着传输时隙数的增加而线性减少,且其最优的散射矩阵集需要用计算机做最优搜索。为了进一步提高STSK的频谱效率,在文献中,Sugiura等人通过在一个GSTSK信号传输时隙内激活多个散射矩阵,提出了广义空时移键控调制方案(Generalized Space-time Shift Keying,GSTSK)。文献将空时编码和SM结合起来提出了空时分组码空间调制方案(Space-Time Block Coded SpatialModulation,STBC-SM)。利用Alamouti空时编码的正交性,该方案可以实现低复杂度的最大似然译码(Maximum Likelihood,ML)。但是在STBC-SM方案中,为了取得二阶的发射分集,需要对旋转角度进行优化,同时空间维度调制所提供的频谱效率较低。为了提高STBC-SM方案的频谱效率,文献提出了一种基于循环结构的STBC-SM方法(STBC-CSM)。尽管STBC-CSM较STBC-SM系统的频谱效率有所提高,但需要优化的角度个数也相应地增多了。显然文献中的散射矩阵的最优搜索和文献中的角度优化都增加了MIMO系统的设计复杂度。随后,文献通过引入空间星座(SC)矩阵的概念,提出了一种高速率的正交STBC-SM方案,称为SM-OSTBC。SM-OSTBC方案可以取得二阶的发射分集而不需要做任何最优搜索和角度的优化。不过,SM-OSTBC方法仅仅适用于偶数根发射天线和射频(RF)链路的MIMO系统,同时发射端至少需要配置4根射频链路。在最新的研究文献中,文献在SM-OSTBC的设计基础上,提出了一种基于双分集的空间调制方法(DT-SM),该方案以二重空时发射分集(DSTTD)码为基本编码矩阵,用SC矩阵与DSTTD相乘来构造DT-SM的编码矩阵,DT-SM可以获得比STBC-SM以及SM-OSTBC更高的频谱效率和更好的误码性能,但是为了获得二阶分集,其SC矩阵中依然有多个旋转角需要进行最优搜索。此外,文献提出了一种可获得发射分集的正交空间调制(DA-QSM)方法,尽管该方案无需作任何的参数优化或者弥散矩阵搜索,但是在某些参数配置下无法保证获得二阶分集且不支持快速译码。
【发明内容】
本发明的目的在于针对现有的SM方案中存在的天线数配置不够灵活、需要进行参数搜索或者矩阵优化等问题,提供一种支持快速译码的增强型可获得发射分集的正交空间调制方法
为达到上述目的,本发明采用以下技术方案予以实现:
支持快速译码的增强型可获得发射分集的正交空间调制方法,包括以下步骤:
步骤1:设计两组弥散矩阵集合A和B,分别存储在发射端;所述弥散矩阵集合A和B中分别有Q个弥散矩阵;
步骤2:信息比特进入发射机,将所有的信息比特分为三部分,其中第一部分用于在弥散矩阵集合A中选择P个矩阵,第二部分用于在弥散矩阵集合B中选择P个矩阵,第三部分用于映射成符号;
步骤3:将从集合A中选择出来的弥散矩阵与符号的实部相乘,将从集合B中选择出来的弥散矩阵与符号的虚部相乘,再将两个相乘的结果相加后发送到无线信道中;
步骤4:在接收端,对弥散矩阵集合A和B以及实符号进行排序,构建状态码字集合,并对接收信号做实数化处理;
步骤5:再采用状态码字匹配的逐块球形译码方法对输入的信息比特进行检测。
本发明进一步的改进在于:
所述步骤1中,设计两组弥散矩阵集合A和B的具体方法如下:
令参数ek表示K阶单位阵IK的第k列,则弥散矩阵集合和中的2Q个弥散矩阵设计为:
其中,表示Kronecker积,2×2维矩阵Ci和Di是SSB码的八个弥散矩阵,定义为:
其中,a=c=1,d=-jb;每一个弥散矩阵集合中包含有Q=2nT个弥散矩阵,且每个弥散矩阵满足和q=1,…,Q。
所述步骤2的具体方法如下:
在T个符号间隔内,总共B=2log2f(Q,P)+Plog2M个比特进入发射机,其中表示从Q中取P的组合数,表示对x向下取整,且为2的幂次方,串并变换后这些比特数被分为三部分:
第一个部分的B1=log2f(Q,P)个比特用来从集合的Q个弥散矩阵中激活P个,用序号向量km=[k(1),k(2),…,k(P)]来表示集合A中被激活的P个弥散矩阵的序号,m=1,…,f(Q,P);
第二个部分的B1=log2f(Q,P)个比特用来从集合的Q个弥散矩阵中激活P个,用序号向量ln=[l(1),l(2),…,l(P)]来表示集合B中被激活的P个弥散矩阵的序号,n=1,…,f(Q,P);
第三个部分的B2=Plog2M个比特被调制为P个M-PSK/QAM符号其中和分别表示符号sp的实部和虚部。
所述步骤3的具体方法如下:
用激活的第p个弥散矩阵Ak(p)和Bl(p)分别对第p个符号sp的实部和虚部进行调制,再将所有的2P项按照如下方式进行线性组合,得到nT×T维的EDA-QSM发射信号:
其中,符号间隔T=2。
所述步骤4中,对接收信号做实数化处理的具体方法如下:
当发射端发送(3)式中的nT×2维EDA-QSM信号时,则nR×2维接收信号表示为:
其中,H和V分别是nR×nT维和nR×2维的信道矩阵和加性高斯噪声矩阵,H和V中的每个元素分别服从和)分布,N0为噪声方差;对(4)式两边同时作按列拉直运算vec(·)和实数化运算(·),得到等价的实的接收信号:
对于复向量x=[x1,x2,…,xn]T,运算(·)表示其中,是与序号向量km=[k(1),k(2),…,k(P)]和ln=[l(1),l(2),…,l(P)]相对应的包含有2P个实符号的符号向量,若P=Q,即所有的2Q个弥散矩阵都被同时激活的话,则 是4nR×4nT维等价实信道矩阵:
其中,是信道矩阵H的实数化矩阵,其中的每个元素按照的形式扩展,因而的维数是2nR×2nT,4nT×2Q维矩阵G是EDA-QSM信号的生成矩阵:
G=[vec(A1),vec(B1),vec(A2),vec(B2)…vec(AQ),vec(BQ)] (7)
对等价信道矩阵作QR分解,得其中矩阵Q是一个标准正交阵,矩阵R是一个2Q×2Q维的上三角矩阵,将代入(5)式后,式子两边再同时左乘Q的转置,即QT,由此得EDA-QSM的最大似然译码度量为:
其中y=QTvec(Y);因此,EDA-QSM方案的ML译码器表示为:
所述步骤4中对弥散矩阵集合A和B进行排序的方法如下:
按照如下的顺序把2Q个弥散矩阵分为nT组:
上三角阵R具有如下的分块正交结构:
其中,每个子块Eij表示其中的元素为任意值的4×4维非零矩阵,i=1,…,Γ-1,j=2,…,Γ;每个子块Dγ是一个4×4维的对角阵,γ=1,…,Γ,Γ是R中子块Dγ的个数。
所述步骤4中实符号的排序方法如下:
步骤4-1:定义一个与km和ln相对应的2×Q维矩阵矩阵的第(1,k(p))个和第(2,l(p))个元素分别为和
步骤4-2:对矩阵作按列拉直操作,得到2Q×1维信号向量从中看出2P个实符号的排列顺序;
步骤4-3:对每个实符号向量定义2Q×1维的状态码字cm,n;令和cm,n(i)分别表示和cm,n的第i个元素,若是非零值,则cm,n(i)=1,反之,若的值为零,则cm,n(i)=0;把所有的f2(Q,P)个状态码字放到集合中,状态码字集合表示为:
所述步骤5具体包括逐块球形检测和状态码字匹配,具体如下:
步骤5-1:逐块球形检测
将(8)式中的符号向量用x表示,ML度量表示为对(9)式应用SD算法时,假定第l层与R中的第γ个块Dγ相对应,则第l层的当前欧式度量表示为:
式中:yl和xl分别表示向量y和x的第l个元素,r(l,l)表示R的第(l,l)个元素,r(l,l+1:L)和x(l+1:L)分别表示向量[r(l,l+1),…,r(l,L)]和[xl+1,…,xL]T,如果第l层与R中的最后一个块DΓ相对应时,则ξl=0,zl=yl;
令表示第γ个块的当前欧式度量,与第γ个块相对应的四层分别为l=4γ-3,l=4γ-2,l=4γ-1和l=4γ层,这样表示为:
令表示第γ个块的累积欧式度量,γ=1,…,Γ,则有:
则:
假定球形检测的初始半径为d,球形译码从第Γ个块逐块向第一个块搜索,在第γ个块,检查是否满足条件:
若上述条件满足则继续向下一个块搜索,直到所有的分支都被检查过;
为了尽可能快地找到(8)式的最优解,在每一层对所有地节点进行排序;令Ω表示一个正方形M-QAM星座,ΩR表示Ω的实部,则ΩR的大小为由于x中共有2P层符号是非零的,剩余(L-2P)层符号为零,因而将第l层的星座定义为在第γ个块,该块所对应的4层符号分别估计为:
xl=round(zl/r(l,l)),(l=4γ-3,4γ-2,4γ-1,4γ) (17)
式中函数round(x)表示为取离x最近的整数;若round(zl/r(l,l))∈Φl,则将round(zl/r(l,l))作为Φl的第一个元素,若则令xl等于Φl的第一个元素;第l层的上述排序方法用函数来执行;
利用经过排序后的四层星座[Φ4γ-3;Φ4γ-2;Φ4γ-1;Φ4γ;],得到第γ个块所对应的所有可能的个分支,得到第γ个块的所有个分支的操作用函数Vecγ=GenVec(Φ(4γ-3:4γ))来执行;个分支一方面用(16)式中的条件进行修剪;
步骤5-2:状态码字匹配
在第γ个块时,首先计算每个分支x(4γ-3:4γ)的状态向量:
b=state(x(4γ-3:4γ)) (18)
令b(4-t)表示b的第t个元素,t=3,2,1,0;若x4γ-t为非零值,则b(4-t)=1,若x4γ-t=0,则b(4-t)=0;
将第γ个块的所有个分支的状态与状态码字集合中与第γ个块所对应的4层状态码字进行对比,把凡是状态能够匹配上的分支保留下来作为第γ个块的候选分支candγ,令Nγ表示candγ中的分支数;
状态码字集合中共有f2(Q,P)个合法的状态码字,根据符号向量x的当前判决值,删除掉中的那些不能与x的当前判决值匹配上的合法状态码字,凡是与x的当前判决值的状态能够匹配上的合法状态码字依然保留在集合中,把不能匹配上的状态码字删除后扔到另一个集合中;
以上的逐块球形搜索和码字匹配方法相互配合工作,直到所有的Γ个块中的所有候选分支全部都被检查完;因此所提出的SCMB-SD算法能够找到(9)式的最优解;最终,保留在集合中的码字即为最优状态码字copt,从copt能够逆映射出两个序号向量km=[k(1),k(2),…,k(P)]和ln=[l(1),l(2),…,l(P)],保留在x中的即为(9)式的最优解xopt。
选择候选分支的操作用函数实现,只有候选分支candγ中的分支在球形搜索的时候才进行检查;
删除合法状态码字的操作用函数执行。
当在检查candγ中的第i个候选分支时,若该候选分支的欧式度量满足条件,就从中把能够与该候选分支匹配上的状态码字从中再召回到集合中去;该召回的操作用函数执行;
以上的逐块球形搜索和码字匹配方法相互配合工作,直到所有的Γ个块中的所有候选分支全部都被检查完;保留在集合中的码字即为最优状态码字copt,从copt能够逆映射出两个序号向量km=[k(1),k(2),…,k(P)]和ln=[l(1),l(2),…,l(P)],保留在x中的为最优解xopt。
与现有技术相比,本发明具有以下有益效果:
由于EDA-QSM的特殊结构和弥散矩阵集合的构造方法,使得本发明中的EDA-QSM具有如下的一些优点:1、适合于任意偶数根发射天线nT,适合于1~nT之间的任意根激活天线;2、由于EDA-QSM方案中每个弥散矩阵集合中弥散矩阵的数量较传统方法翻倍,因而EDA-QSM获得了非常显著的频谱效率提升;3、无需作任何的参数或者矩阵优化,EDA-QSM方案具有永不消失的行列式(NVD)特性,该特性可以保证EDA-QSM方案获得二阶发射分集;4、当所有的2Q个弥散矩阵按照一定的顺序排列后,EDA-QSM方案在编码结构上具有分块正交特性,因而支持快速译码,具有很低的译码复杂度;5、EDA-QSM的分块正交特性使得同一个块中所对应的4个符号之间相互独立,因而可以实现逐块的球形检测,每次并行地检测一个块中的4个符号,同时与EDA-QSM的状态码字进行匹配来提高球形搜索的速度。由于SCMB-SD算法仅需执行一次球形搜索就可以同时估计出两组激活的序号向量和P个符号,因而具有很低的解码复杂度。为了表述简单,将具有nT根发射天线、nR根接收天线和每个集合中激活P个弥散矩阵的EDA-QSM方案表示为EDAQSM(nT,nR,P)。
【附图说明】
图1为本发明EDA-QSM传输方案发射端结构框图;
图2为本发明EDA-QSM、DT-SM与DA-QSM随着发射天线数变化时频谱效率的比较图;
图3为本发明在发射天线数分别为8、16和24时,EDA-QSM与DA-QSM随着每个集合中激活的弥散矩阵数变化时频谱效率的比较图;
图4为本发明EDA-QSM、DT-SM、DA-QSM与STBC-SM在频谱效率为4bits/s/Hz时BER性能的比较图;
图5为本发明EDA-QSM、DT-SM、DA-QSM、STBC-SM和SM-OSTBC在频谱效率为6bits/s/Hz时BER性能的比较;
图6为本发明EDA-QSM、DT-SM和DA-QSM在频谱效率为8bits/s/Hz时BER性能的比较图。
【具体实施方式】
为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分的实施例,不是全部的实施例,而并非要限制本发明公开的范围。此外,在以下说明中,省略了对公知结构和技术的描述,以避免不必要的混淆本发明公开的概念。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
在附图中示出了根据本发明公开实施例的各种结构示意图。这些图并非是按比例绘制的,其中为了清楚表达的目的,放大了某些细节,并且可能省略了某些细节。图中所示出的各种区域、层的形状及它们之间的相对大小、位置关系仅是示例性的,实际中可能由于制造公差或技术限制而有所偏差,并且本领域技术人员根据实际所需可以另外设计具有不同形状、大小、相对位置的区域/层。
本发明公开的上下文中,当将一层/元件称作位于另一层/元件“上”时,该层/元件可以直接位于该另一层/元件上,或者它们之间可以存在居中层/元件。另外,如果在一种朝向中一层/元件位于另一层/元件“上”,那么当调转朝向时,该层/元件可以位于该另一层/元件“下”。
需要说明的是,本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
下面结合附图对本发明做进一步详细描述:
参见图1,本发明支持快速译码的增强型可获得发射分集的正交空间调制(Enhanced Diversity-Achieving Quadrature Spatial Modulation,EDA-QSM)传输方案采用了两组弥散矩阵(Dispersion Matrices,DM)集合,其中一组用于将待传输符号的实部在空-时两个维度上进行扩展,另一组用于将待传输符号的虚部在空-时两个维度上进行扩展。EDA-QSM算法的发射端实现框图如图1所示。本发明中的EDA-QSM方案较现有的其他可获得发射分集的SM传输方案的显著优势是其中无需对任何参数做最优搜索,且发射天线数及激活天线数更为灵活(发射天线数可以是大于3的任意偶数,激活天线数目任意),此外,较现有方案具有明显的频谱效率提升和误码(BER)性能改善。本发明对频谱效率的提升和BER性能的提升效果见图2~图6。
考虑一个具有nT个发射天线和nR个接收天线的MIMO系统,EDA-QSM传输算法及其接收端的快速译码方法,包括以下步骤:
步骤1:设计两组弥散矩阵集合A和B,分别存储在发射端;所述弥散矩阵集合A和B中分别有Q个弥散矩阵;
为了避免弥散矩阵构造过程中复杂的最优搜索,在EDA-QSM方案中采用了对SSB码的八个弥散矩阵作Kronecke r积扩展的方法来构造弥散矩阵。令参数ek表示K阶单位阵IK的第k列,则弥散矩阵集合和中的2Q个弥散矩阵可以设计为:
其中表示Kronecker积,2×2维矩阵Ci和Di i=1,…,4是SSB码的八个弥散矩阵,它们的定义为:
其中a=c=1,d=-jb。从(2)式中可以看出,每一个弥散矩阵集合中包含有Q=2nT个弥散矩阵,且很容易验证每个弥散矩阵满足和(q=1,…,Q)。
以nT=4为例,集合中的Q=2nT=8个弥散矩阵为:
集合中的Q=2nT=8个弥散矩阵为:
步骤2:信息比特进入发射机,将所有的信息比特分为三部分,其中第一部分用于在弥散矩阵集合A中选择P个矩阵,第二部分用于在弥散矩阵集合B中选择P个矩阵,第三部分用于映射成符号;
如图1所示,发射机结构:假定有两组nT×T维的弥散矩阵(Dispersion Matrices,DM)集合和预先设置在发射端,每个集合中包含有Q个弥散矩阵。EDA-QSM的实现需要如下两步:
第一步:在T个符号间隔内,总共B=2log2f(Q,P)+Plog2M个比特进入发射机,其中表示从Q中取P的组合数,表示对x向下取整,且为2的幂次方)。串并变换(S/P)后这些比特数被分为三部分:第一个部分的B1=log2f(Q,P)个比特用来从第一组DM集合的Q个弥散矩阵中激活P个,用序号向量km=[k(1),k(2),…,k(P)](m=1,…,f(Q,P))来表示集合A中被激活的P个弥散矩阵的序号;第二个部分的B1=log2f(Q,P)个比特用来从第二组DM集合的Q个弥散矩阵中激活P个,用序号向量ln=[l(1),l(2),…,l(P)](n=1,…,f(Q,P))来表示集合B中被激活的P个弥散矩阵的序号;第三个部分的B2=Plog2M个比特被调制为P个M-PSK/QAM符号(p=1,...,P),其中和分别表示符号sp的实部和虚部;
EDA-QSM的主要特征
1)EDA-QSM的频谱效率
EDA-QSM方案在T=2个符号周期内发送了B=2log2f(Q,P)+Plog2M个比特,因而EDA-QSM方案的频谱效率为
上述频谱效率明显高于现有的其他可以获得发射分集的空间调制方案。说明书附图2和图3中分别给出了本发明中的EDA-QSM方案与最新的DT-SM方案以及DA-QSM方案的频谱效率随着发射天线数nT和激活的弥散矩阵数P变化时的频谱效率的比较。从图2和图3中可见,EDA-QSM方案的频谱效率具有明显的优势。
2)满足NVD特征
根据空时编码的设计准则,对于任意两个不同的空时码字X和编码增益定义为
经过分析,任意两个EDA-QSM码字X和之间的最小编码增益始终大于零,即δmin(X,X)>0,因而EDA-QSM方案具有永不消失的行列式(Non-Vanishing Determinant,NVD)特性,因此,EDA-QSM在无需作任何参数优化及矩阵搜索的情况下可以保证获得二阶的发射分集。
3)具有分块正交结构
所提出的EDA-QSM方案具有分块正交结构,该结构使得EDA-QSM支持快速译码。在EDA-QSM传输的接收端,对接收信号作按列拉直及实数化处理后,再对等价的信道矩阵作QR分解,得到的上三角矩阵R则具有如下的形式:
其中,每个子块Eij(i=1,…,Γ-1,j=2,…,Γ)表示其中的元素为任意值的4×4维非零矩阵;每个子块Dγ(γ=1,…,Γ)是一个4×4维的对角阵,Γ是R中子块Dγ的个数,因而与每个子块Dγ所对应的4个符号是相互独立的。
步骤3:将从集合A中选择出来的弥散矩阵与符号的实部相乘,将从集合B中选择出来的弥散矩阵与符号的虚部相乘,再将两个相乘的结果相加后发送到无线信道中;
用激活的第p个弥散矩阵Ak(p)和Bl(p)分别对第p个符号sp的实部和虚部进行调制,再将所有的2P项按照如下方式进行线性组合即得到了一个nT×T维的EDA-QSM发射信号:
与现有的其它空时空间调制算法一致,本发明中采用的符号间隔T=2。
步骤4:在接收端,对弥散矩阵集合A和B以及实符号进行排序,构建状态码字集合,并对接收信号做实数化处理;
1)对接收信号做实数化处理
当发射端发送(1)式中的nT×2维EDA-QSM信号时,则nR×2维接收信号可以表示为:
其中H和V分别是nR×nT维和nR×2维的信道矩阵和加性高斯噪声矩阵,H和V中的每个元素分别服从和)分布,N0为噪声方差。对(3)式两边同时作按列拉直运算vec(·)和实数化运算(·)(对于复向量x=[x1,x2,…,xn]T,运算(·)表示),可得一个等价的实的接收信号:
其中是与序号向量km=[k(1),k(2),…,k(P)]和ln=[l(1),l(2),…,l(P)]相对应的包含有2P个实符号的符号向量,若P=Q,即所有的2Q个弥散矩阵都被同时激活的话则是4nR×4nT维等价实信道矩阵,它具有如下形式:
其中是信道矩阵H的实数化矩阵,其中的每个元素按照的形式扩展,因而的维数是2nR×2nT,4nT×2Q维矩阵G是EDA-QSM信号的生成矩阵:
G=[vec(A1),vec(B1),vec(A2),vec(B2)…vec(AQ),vec(BQ)] (6)
对等价信道矩阵作QR分解,可得其中矩阵Q是一个标准正交阵,矩阵R是一个2Q×2Q维的上三角矩阵,将代入(4)式后,式子两边再同时左乘Q的转置,即QT,由此可得EDA-QSM的最大似然(ML)译码度量为:
其中y=QTvec(Y)。因此,EDA-QSM方案的ML译码器可以表示为:
在上三角矩阵R的右上方若出现零元素,则该编码方案可以支持快速译码,若右上方全部都是非零元素的话则该编码不能支持快速译码。
2)弥散矩阵的排序方法
按照如下的顺序把2Q个弥散矩阵分为nT组:
由于对SSB码的弥散矩阵采用了(2)式中的Kronecker积扩展方法,很容易验证,(9)式中每个组的4个弥散矩阵之间两两满足斜对称特性,而属于不同组的弥散矩阵之间不满足斜对称特性,因而若在生成矩阵G中按照(9)中的顺序对2Q个弥散矩阵排列时,得到的上三角阵R就会具有如下的分块正交结构:
其中,每个子块Eij(i=1,…,Γ-1,j=2,…,Γ)表示其中的元素为任意值的4×4维非零矩阵;每个子块Dγ(γ=1,…,Γ)是一个4×4维的对角阵,Γ是R中子块Dγ的个数,这样与每个子块Dγ所对应的4个符号之间是相互独立的。需要注意的是,如果改变G中弥散矩阵的排列方法,上三角阵R中的分块对角结构将不再出现,因此在后文中弥散矩阵都采用(9)式中的排列方法。
3)实符号的排序方法
若P=Q,即所有的2Q个弥散矩阵都被同时激活的话则与(9)式中弥散矩阵的排列顺序相对应,2Q个实符号的排列顺序为若P<Q,则2P个实符号具有与(9)式不同的排列顺序,这是因为是用Ak(p),而不是Ap来调制的,同样,是用Bl(p),而不是Bp,来调制的。以下给出在P<Q时2P个实符号的排列方法。
在EDAQSM(nT,nR,P)方案中,在发送端,假定两个序号向量km=[k(1),k(2),…,k(P)]和ln=[l(1),l(2),…,l(P)](m,n=1,…,f(Q,P))分别用来从集合和中激活P个弥散矩阵。2P个实符号的排列方法可以由以下三步实现:
第一步:定义一个与km和ln相对应的2×Q维矩阵矩阵的第(1,k(p))个和第(2,l(p))个元素分别为和
第二步:对矩阵作按列拉直操作,得到一个2Q×1维信号向量从中就可以看出2P个实符号的排列顺序。
第三步:为了更为清楚地描述2Q个弥散矩阵的激活状态,对每个实符号向量定义一个2Q×1维的状态码字cm,n。令和cm,n(i)分别表示和cm,n的第i个元素,若是非零值,则cm,n(i)=1,反之,若的值为零,则cm,n(i)=0。很显然,状态码字cm,n与两个序号向量km和ln之间存在一一对应的关系,进一步,把所有的f2(Q,P)个状态码字放到一个集合中,因此状态码字集合可以表示为:
以EDAQSM(4,nR,2)为例,在表1中列出了EDAQSM(4,nR,2)的符号排列顺序和对应的状态码字,由于EDAQSM(4,nR,2)共有f2(Q,P)=f2(8,2)=256个状态码字,为节省空间,表1中只列出了10个状态码字。
表1.EDAQSM(4,nR,2)的状态码字
步骤5:再采用状态码字匹配的逐块球形(SCMB-SD)译码方法对输入的信息比特进行检测。
1)逐块球形检测
由于上三角阵R的对角线上包含了Γ个块Dγ(γ=1,…,Γ),每个块Dγ是一个4×4维的对角阵,因而与每个块相对应的4个实符号之间是相互正交且独立的。利用这一特性,可以对这4个符号实现并行检测,因此在。应用球形译码(Sphere Decoding,SD)算法时,并不是逐层进行搜索,而是逐块进行搜索,每个块中对4个符号并行搜索。这里为表达简便,将(7)式中的符号向量用x表示,因此(7)式中的ML度量可以表示为对(8)式应用SD算法时,假定第l层与R中的第γ个块Dγ相对应,则第l层的当前欧式度量可以表示为:
式中:yl和xl分别表示向量y和x的第l个元素,r(l,l)表示R的第(l,l)个元素,r(l,l+1:L)和x(l+1:L)分别表示向量[r(l,l+1),…,r(l,L)]和[xl+1,…,xL]T,如果第l层与R中的最后一个块DΓ相对应时,则ξl=0,zl=yl。
令表示第γ个块的当前欧式度量,与第γ个块相对应的四层分别为l=4γ-3,l=4γ-2,l=4γ-1和l=4γ层,这样可以表示为:
令表示第γ个块的累积欧式度量,则有:
很显然,
假定球形检测的初始半径为d,球形译码从第Γ个块逐块向第一个块搜索,在第γ个块,检查是否满足条件:
若上述条件满足则继续向下一个块搜索,直到所有的分支都被检查过。
为了尽可能快地找到(8)式的最优解,在每一层对所有地节点进行排序。令Ω表示一个正方形M-QAM星座,ΩR表示Ω的实部,则ΩR的大小为由于x中共有2P层符号是非零的,剩余(L-2P)层符号为零,因而将第l层的星座定义为在第γ个块,该块所对应的4层符号分别估计为:
xl=round(zl/r(l,l)),(l=4γ-3,4γ-2,4γ-1,4γ) (15)
式中函数round(x)表示为取离x最近的整数。若round(zl/r(l,l))∈Φl,则将round(zl/r(l,l))作为Φl的第一个元素,若则令xl等于Φl的第一个元素。第l层的上述排序方法用函数[xl,ml CurL,Φl]=OpeEachLayer(l,zl,r(l,l),Φl)来执行。
利用经过排序后的四层星座[Φ4γ-3;Φ4γ-2;Φ4γ-1;Φ4γ;],可以得到第γ个块所对应的所有可能的个分支,得到第γ个块的所有个分支的操作用函数Vecγ=GenVec(Φ(4γ-3:4γ))来执行。个分支一方面可以用(14)式中的条件进行修剪,另外,还可以用下面的状态码字匹配方法来进行快速地和有效地修剪。
2)状态码字匹配方法
在第γ个块时,首先计算每个分支x(4γ-3:4γ)的状态向量:
b=state(x(4γ-3:4γ)) (16)
令b(4-t)表示b的第t(t=3,2,1,0)个元素,若x4γ-t为非零值,则b(4-t)=1,若x4γ-t=0,则b(4-t)=0。随后,将第γ个块的所有个分支的状态与状态码字集合中与第γ个块所对应的4层状态码字进行对比,把凡是状态能够匹配上的分支保留下来作为第γ个块的候选分支candγ,令Nγ表示candγ中的分支数。上述选择候选分支的操作可以用函数来实现,只有候选分支candγ中的分支在球形搜索的时候才进行检查,没有匹配上的分支不检查。这样就可以大大减少球形搜索的次数。
状态码字集合中共有f2(Q,P)个合法的状态码字,这个数量在天线数较多的时候是非常大的,为了加速球形搜索的速度,根据符号向量x的当前判决值,删除掉中的那些不能与x的当前判决值匹配上的合法状态码字,凡是与x的当前判决值的状态能够匹配上的合法状态码字依然保留在集合中,把不能匹配上的状态码字删除后扔到另一个集合中。上述删除合法状态码字的操作可以用函数来执行。采用这种删除合法状态码字的方式,随着检测不断向下一个块进行(即γ值的减小),每个块中的候选分支candγ的个数Nγ将显著减小。
然而,上述删除合法状态码字的做法也面临着一些风险,即在第γ个块的所有分支还没有被检查完的时候最优状态码字copt就被删掉了。为了避免这种情况出现,当在检查candγ中的第i个候选分支(用candγ(:,i)表示)时,若该候选分支的欧式度量满足条件(14),就从中把可以与该候选分支匹配上的状态码字从中再召回到集合中去。该召回的操作用函数来执行。
以上的逐块球形搜索和码字匹配方法相互配合工作,直到所有的Γ个块中的所有候选分支全部都被检查完。因此所提出的SCMB-SD算法能够找到(8)的最优解。最终,保留在集合中的码字即为最优状态码字copt,从copt可以逆映射出两个序号向量km=[k(1),k(2),…,k(P)]和ln=[l(1),l(2),…,l(P)],保留在x中的即为(8)的最优解xopt。
具体地,针对EDA-QSM传输的SCMB-SD检测算法可以由以下方法实现。(算法中的符号说明:1L×1表示L×1维列向量,且该列向量的每个元素都是1;sum(x)表示对向量x的所有元素求和;表示空集;size(R,2)表示求矩阵R的列数)
输入参数:y,R,ΩR,球形初始半径d0
输出:(8)式的最优解xopt,最优状态码字copt
1)初始化:令L=size(R,2),Γ=L/4,
2)令d=d0,γ=Γ,nγ=1;
3)For t=4 to 1,令l=4(γ-1)+t,zl=yl,执行
4)计算然后执行Vecγ=GenVec(Φ(4γ-3:4γ))和
5)如果Tγ<d且则执行否则,跳转至Step 9);
6)如果γ=1,则跳转至Step 8),否则,For t=4 to 1,令l=4(γ-2)+t,计算ξl=r(l,4γ-3:L)x(4γ-3:L)和zl=yl-ξl,然后执行
7)计算和设γ=γ-1,然后执行Vecγ=GenVec(Φ(4γ-3:4γ))和设置nγ=1,跳转至Step 5);
8)设置d=Tγ,保存xopt=x和
9)如果则跳转至Step 11);否则,设置nγ=nγ+1。若nγ>Nγ,则设γ=γ+1,否则跳转至Step 11);
10)如果γ≤Γ,则跳转至Step 9),否则跳转至Step 14);
11)For t=4 to 1,令l=4(γ-1)+t,如果xl≠candγ(t,nγ),计算
12)计算令x(4γ-3:4γ)=candγ(:,nγ)。如果γ≠Γ,则计算否则
13)如果Tγ<d,则执行如果则跳转至Step 5),否则跳转至Step 9);
14)如果(即没有找到解),则增加d0,然后跳转至Step 2),否则终止程序。
仿真实验
对所提出的EDA-QSM算法的误码性能进行蒙特卡洛仿真并与现有的方案进行比较。在所有的仿真图中横轴表示每个接收天线处的信噪比(SNR),纵轴为误比特率(BER),且所有仿真中接收天线的个数均设为nR=4,性能比较均是在BER值为10-5时所做的。
图4给出了EDA-QSM、DT-SM、DA-QSM与STBC-SM在频谱效率为4bits/s/Hz时的BER比较。由图中可见,DT-SM(4,4,4)和DA-QSM(4,4,2)几乎具有相同的BER性能,它们的性能都好于采用了nT=4个发射天线的STBC-SM方案。而EDAQSM(4,4,1)比DT-SM(4,4,4)和DA-QSM(4,4,2)都具有明显的性能优势,它比DT-SM(4,4,4)和DA-QSM(4,4,2)具有大约0.8dB的性能增益。图4中同时给出了EDAQSM(4,4,1)的理论BER曲线,可以看出,EDA-QSM的理论BER曲线和仿真所得的BER曲线在高SNR时完全相符。
图5中比较了几种方案在频谱效率为6bits/s/Hz时的BER比较。从图5可见,在6bits/s/Hz时,DA-QSM(5,4,4)具有比采用nT=8个发射天线的STBC-SM方案和SM-OSTBC C(4,4,4)都具有更好的性能。EDAQSM(4,4,2)和DT-SM(4,4,4)在高SNR时几乎具有相同的BER性能,而前者激活的发射天线数更少。
图6中比较了EDA-QSM和DT-SM以及DA-QSM方案在频谱效率为8bits/s/Hz时的BER性能。从图中可见,DT-SM(4,4,4)的性能优于DA-QSM(4,4,4),而EDAQSM(4,4,3)的性能明显优于DT-SM(4,4,4),且比DT-SM(4,4,4)具有大约0.8dB的性能增益。这是因为尽管EDAQSM(4,4,3)只传输了3个符号,少于DT-SM(4,4,4)传输的符号数,而EDAQSM(4,4,3)却具有f2(Q,P)=f2(8,3)=1024个序号向量的组合。而且从图6中还可以看出,EDAQSM(4,4,8)的性能也明显优于DT-SM(4,4,4),且有有大约2dB的性能增益。这里尽管EDAQSM(4,4,8)的序号向量组合f2(Q,P)=f2(8,8)=1,但是它能够传输8个符号。图6中同时给出了EDAQSM(4,4,8)在采用ML译码和SCMB-SD译码时的BER性能,可以看出,两种检测器的BER性能完全相同,但是SCMB-SD译码具有非常显著的译码复杂度的降低。
综上,由图4~图6的仿真实验可以看出,EDA-QSM方案较现有的几种典型的空时SM传输方案如STBC-SM、SM-OSTBC以及最新提出的DA-QSM和DT-SM方案均具有明显的性能优势。
以上内容仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明权利要求书的保护范围之内。
Claims (10)
1.支持快速译码的增强型可获得发射分集的正交空间调制方法,其特征在于,包括以下步骤:
步骤1:设计两组弥散矩阵集合A和B,分别存储在发射端;所述弥散矩阵集合A和B中分别有Q个弥散矩阵;
步骤2:信息比特进入发射机,将所有的信息比特分为三部分,其中第一部分用于在弥散矩阵集合A中选择P个矩阵,第二部分用于在弥散矩阵集合B中选择P个矩阵,第三部分用于映射成符号;
步骤3:将从集合A中选择出来的弥散矩阵与符号的实部相乘,将从集合B中选择出来的弥散矩阵与符号的虚部相乘,再将两个相乘的结果相加后发送到无线信道中;
步骤4:在接收端,对弥散矩阵集合A和B以及实符号进行排序,构建状态码字集合,并对接收信号做实数化处理;
步骤5:再采用状态码字匹配的逐块球形译码方法对输入的信息比特进行检测。
2.根据权利要求1所述的支持快速译码的增强型可获得发射分集的正交空间调制方法,其特征在于,所述步骤1中,设计两组弥散矩阵集合A和B的具体方法如下:
在一个具有nT个发射天线和nR个接收天线的MIMO系统中,令参数ek表示K阶单位阵IK的第k列,则弥散矩阵集合和中的2Q个弥散矩阵设计为:
其中,表示Kronecker积,2×2维矩阵Ci和Di是SSB码的八个弥散矩阵,定义为:
其中,a=c=1,d=-jb;每一个弥散矩阵集合中包含有Q=2nT个弥散矩阵,且每个弥散矩阵满足和q=1,…,Q。
3.根据权利要求1所述的支持快速译码的增强型可获得发射分集的正交空间调制方法,其特征在于,所述步骤2的具体方法如下:
在T个符号间隔内,总共B=2log2 f(Q,P)+Plog2 M个比特进入发射机,其中 表示从Q中取P的组合数,表示对x向下取整,且为2的幂次方,串并变换后这些比特数被分为三部分:
第一个部分的B1=log2 f(Q,P)个比特用来从集合的Q个弥散矩阵中激活P个,用序号向量km=[k(1),k(2),…,k(P)]来表示集合A中被激活的P个弥散矩阵的序号,m=1,…,f(Q,P);
第二个部分的B1=log2 f(Q,P)个比特用来从集合的Q个弥散矩阵中激活P个,用序号向量ln=[l(1),l(2),…,l(P)]来表示集合B中被激活的P个弥散矩阵的序号,n=1,…,f(Q,P);
第三个部分的B2=Plog2 M个比特被调制为P个M-PSK/QAM符号其中和分别表示符号sp的实部和虚部。
4.根据权利要求3所述的支持快速译码的增强型可获得发射分集的正交空间调制方法,其特征在于,所述步骤3的具体方法如下:
用激活的第p个弥散矩阵Ak(p)和Bl(p)分别对第p个符号sp的实部和虚部进行调制,再将所有的2P项按照如下方式进行线性组合,得到nT×T维的EDA-QSM发射信号:
其中,符号间隔T=2。
5.根据权利要求4所述的支持快速译码的增强型可获得发射分集的正交空间调制方法,其特征在于,所述步骤4中,对接收信号做实数化处理的具体方法如下:
当发射端发送(3)式中的nT×2维EDA-QSM信号时,则nR×2维接收信号表示为:
其中,H和V分别是nR×nT维和nR×2维的信道矩阵和加性高斯噪声矩阵,H和V中的每个元素分别服从和分布,N0为噪声方差;对(4)式两边同时作按列拉直运算vec(·)和实数化运算(·),得到等价的实的接收信号:
对于复向量x=[x1,x2,…,xn]T,运算(·)表示其中,是与序号向量km=[k(1),k(2),…,k(P)]和ln=[l(1),l(2),…,l(P)]相对应的包含有2P个实符号的符号向量,若P=Q,即所有的2Q个弥散矩阵都被同时激活的话,则 是4nR×4nT维等价实信道矩阵:
其中,是信道矩阵H的实数化矩阵,其中的每个元素按照的形式扩展,因而的维数是2nR×2nT,4nT×2Q维矩阵G是EDA-QSM信号的生成矩阵:
G=[vec(A1),vec(B1),vec(A2),vec(B2)…vec(AQ),vec(BQ)] (7)
对等价信道矩阵作QR分解,得其中矩阵Q是一个标准正交阵,矩阵R是一个2Q×2Q维的上三角矩阵,将代入(5)式后,式子两边再同时左乘Q的转置,即QT,由此得EDA-QSM的最大似然译码度量为:
其中y=QTvec(Y);因此,EDA-QSM方案的ML译码器表示为:
6.根据权利要求5所述的支持快速译码的增强型可获得发射分集的正交空间调制方法,其特征在于,所述步骤4中对弥散矩阵集合A和B进行排序的方法如下:
按照如下的顺序把2Q个弥散矩阵分为nT组:
上三角阵R具有如下的分块正交结构:
其中,每个子块Eij表示其中的元素为任意值的4×4维非零矩阵,i=1,…,Γ-1,j=2,…,Γ;每个子块Dγ是一个4×4维的对角阵,γ=1,…,Γ,Γ是R中子块Dγ的个数。
7.根据权利要求6所述的支持快速译码的增强型可获得发射分集的正交空间调制方法,其特征在于,所述步骤4中实符号的排序方法如下:
步骤4-1:定义一个与km和ln相对应的2×Q维矩阵矩阵的第(1,k(p))个和第(2,l(p))个元素分别为和
步骤4-2:对矩阵作按列拉直操作,得到2Q×1维信号向量从中看出2P个实符号的排列顺序;
步骤4-3:对每个实符号向量定义2Q×1维的状态码字cm,n;令和cm,n(i)分别表示和cm,n的第i个元素,若是非零值,则cm,n(i)=1,反之,若的值为零,则cm,n(i)=0;把所有的f2(Q,P)个状态码字放到集合中,状态码字集合表示为:
8.根据权利要求7所述的支持快速译码的增强型可获得发射分集的正交空间调制方法,其特征在于,所述步骤5具体包括逐块球形检测和状态码字匹配,具体如下:
步骤5-1:逐块球形检测
将(8)式中的符号向量用x表示,ML度量表示为对(9)式应用SD算法时,假定第l层与R中的第γ个块Dγ相对应,则第l层的当前欧式度量表示为:
式中:yl和xl分别表示向量y和x的第l个元素,r(l,l)表示R的第(l,l)个元素,r(l,l+1:L)和x(l+1:L)分别表示向量[r(l,l+1),…,r(l,L)]和[xl+1,…,xL]T,如果第l层与R中的最后一个块DΓ相对应时,则ξl=0,zl=yl;
令表示第γ个块的当前欧式度量,与第γ个块相对应的四层分别为l=4γ-3,l=4γ-2,l=4γ-1和l=4γ层,这样表示为:
令表示第γ个块的累积欧式度量,γ=1,…,Γ,则有:
则:
假定球形检测的初始半径为d,球形译码从第Γ个块逐块向第一个块搜索,在第γ个块,检查是否满足条件:
若上述条件满足则继续向下一个块搜索,直到所有的分支都被检查过;
为了尽可能快地找到(8)式的最优解,在每一层对所有地节点进行排序;令Ω表示一个正方形M-QAM星座,ΩR表示Ω的实部,则ΩR的大小为由于x中共有2P层符号是非零的,剩余(L-2P)层符号为零,因而将第l层的星座定义为在第γ个块,该块所对应的4层符号分别估计为:
xl=round(zl/r(l,l)),l=4γ-3,4γ-2,4γ-1,4γ (17)
式中函数round(x)表示为取离x最近的整数;若round(zl/r(l,l))∈Φl,则将round(zl/r(l,l))作为Φl的第一个元素,若则令xl等于Φl的第一个元素;第l层的上述排序方法用函数来执行;
利用经过排序后的四层星座[Φ4γ-3;Φ4γ-2;Φ4γ-1;Φ4γ;],得到第γ个块所对应的所有可能的个分支,得到第γ个块的所有个分支的操作用函数Vecγ=GenVec(Φ(4γ-3:4γ))来执行;个分支一方面用(16)式中的条件进行修剪;
步骤5-2:状态码字匹配
在第γ个块时,首先计算每个分支x(4γ-3:4γ)的状态向量:
b=state(x(4γ-3:4γ)) (18)
令b(4-t)表示b的第t个元素,t=3,2,1,0;若x4γ-t为非零值,则b(4-t)=1,若x4γ-t=0,则b(4-t)=0;
将第γ个块的所有个分支的状态与状态码字集合中与第γ个块所对应的4层状态码字进行对比,把凡是状态能够匹配上的分支保留下来作为第γ个块的候选分支candγ,令Nγ表示candγ中的分支数;
状态码字集合中共有f2(Q,P)个合法的状态码字,根据符号向量x的当前判决值,删除掉中的那些不能与x的当前判决值匹配上的合法状态码字,凡是与x的当前判决值的状态能够匹配上的合法状态码字依然保留在集合中,把不能匹配上的状态码字删除后扔到另一个集合中;
以上的逐块球形搜索和码字匹配方法相互配合工作,直到所有的Γ个块中的所有候选分支全部都被检查完;因此所提出的SCMB-SD算法能够找到(9)式的最优解;最终,保留在集合中的码字即为最优状态码字copt,从copt能够逆映射出两个序号向量km=[k(1),k(2),…,k(P)]和ln=[l(1),l(2),…,l(P)],保留在x中的即为(9)式的最优解xopt。
9.根据权利要求8所述的支持快速译码的增强型可获得发射分集的正交空间调制方法,其特征在于,选择候选分支的操作用函数实现,只有候选分支candγ中的分支在球形搜索的时候才进行检查;
删除合法状态码字的操作用函数执行。
10.根据权利要求8或9所述的支持快速译码的增强型可获得发射分集的正交空间调制方法,其特征在于,当在检查candγ中的第i个候选分支时,若该候选分支的欧式度量满足条件,就从中把能够与该候选分支匹配上的状态码字从中再召回到集合中去;该召回的操作用函数执行;
以上的逐块球形搜索和码字匹配方法相互配合工作,直到所有的Γ个块中的所有候选分支全部都被检查完;保留在集合中的码字即为最优状态码字copt,从copt能够逆映射出两个序号向量km=[k(1),k(2),…,k(P)]和ln=[l(1),l(2),…,l(P)],保留在x中的为最优解xopt。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910238613.1A CN110011946B (zh) | 2019-03-27 | 2019-03-27 | 支持快速译码的增强型可获得发射分集的正交空间调制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910238613.1A CN110011946B (zh) | 2019-03-27 | 2019-03-27 | 支持快速译码的增强型可获得发射分集的正交空间调制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110011946A true CN110011946A (zh) | 2019-07-12 |
CN110011946B CN110011946B (zh) | 2020-03-31 |
Family
ID=67168400
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910238613.1A Expired - Fee Related CN110011946B (zh) | 2019-03-27 | 2019-03-27 | 支持快速译码的增强型可获得发射分集的正交空间调制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110011946B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110855328A (zh) * | 2019-10-25 | 2020-02-28 | 西安交通大学 | 一种基于天线分组的差分空间调制方法、设备及存储介质 |
CN113965987A (zh) * | 2020-07-20 | 2022-01-21 | 维沃移动通信有限公司 | Ssb搜索信息确定方法、装置及通信设备 |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11846738B2 (en) | 2019-04-23 | 2023-12-19 | Cerium Laboratories Llc | Radiation detection systems and methods |
EP4334755A2 (en) * | 2021-05-03 | 2024-03-13 | Cerium Laboratories LLC | Radiation detection systems and methods |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101026435A (zh) * | 2006-02-24 | 2007-08-29 | 中国科学院上海微系统与信息技术研究所 | 通信系统中低复杂度的极大似然检测方法及装置 |
CN101557269A (zh) * | 2009-05-18 | 2009-10-14 | 北京天碁科技有限公司 | 一种基于超大规模集成电路的球形译码检测方法 |
CN102651727A (zh) * | 2012-05-14 | 2012-08-29 | 电子科技大学 | 用于多根天线的sm-ofdm系统中的非相干检测方法 |
CN109327287A (zh) * | 2018-09-10 | 2019-02-12 | 西安交通大学 | 一种采用堆叠式Alamouti编码映射的空间调制方法 |
-
2019
- 2019-03-27 CN CN201910238613.1A patent/CN110011946B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101026435A (zh) * | 2006-02-24 | 2007-08-29 | 中国科学院上海微系统与信息技术研究所 | 通信系统中低复杂度的极大似然检测方法及装置 |
CN101557269A (zh) * | 2009-05-18 | 2009-10-14 | 北京天碁科技有限公司 | 一种基于超大规模集成电路的球形译码检测方法 |
CN102651727A (zh) * | 2012-05-14 | 2012-08-29 | 电子科技大学 | 用于多根天线的sm-ofdm系统中的非相干检测方法 |
CN109327287A (zh) * | 2018-09-10 | 2019-02-12 | 西安交通大学 | 一种采用堆叠式Alamouti编码映射的空间调制方法 |
Non-Patent Citations (3)
Title |
---|
MINH-TUAN LE: "Spatially Modulated Orthogonal Space-Time Block Codes with Non-Vanishing Determinants", 《 IEEE TRANSACTIONS ON COMMUNICATIONS》 * |
SHINYA SUGIURA: ""Reduced-Complexity Iterative-Detection-Aided Generalized Space-Time Shift Keying", 《IEEE TRANSACTIONS ON VEHICULAR TECHNOLOGY》 * |
王磊: "Diversity-Achieving Quadrature Spatial Modulation", 《IEEE TRANSACTIONS ON VEHICULAR TECHNOLOGY》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110855328A (zh) * | 2019-10-25 | 2020-02-28 | 西安交通大学 | 一种基于天线分组的差分空间调制方法、设备及存储介质 |
CN110855328B (zh) * | 2019-10-25 | 2021-01-19 | 西安交通大学 | 一种基于天线分组的差分空间调制方法、设备及存储介质 |
CN113965987A (zh) * | 2020-07-20 | 2022-01-21 | 维沃移动通信有限公司 | Ssb搜索信息确定方法、装置及通信设备 |
CN113965987B (zh) * | 2020-07-20 | 2024-02-09 | 维沃移动通信有限公司 | Ssb搜索信息确定方法、装置及通信设备 |
Also Published As
Publication number | Publication date |
---|---|
CN110011946B (zh) | 2020-03-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110011946A (zh) | 支持快速译码的增强型可获得发射分集的正交空间调制方法 | |
CN101039137B (zh) | Mimo-ofdm系统基于码本搜索减少预编码反馈比特数的方法及装置 | |
CN105554865B (zh) | 一种基于stbc的mimo-scma系统下行链路设计方法 | |
Seifi et al. | Media-based MIMO: A new frontier in wireless communications | |
CN105306175A (zh) | 基于v-blast编码方式的mimo-scma系统上行链路构架方法 | |
CN109327287A (zh) | 一种采用堆叠式Alamouti编码映射的空间调制方法 | |
CN107835068A (zh) | 一种具有发射分集的低复杂度正交空间调制球形译码检测算法 | |
CN106788626A (zh) | 一种能够获得二阶发射分集的改进正交空间调制传输方法 | |
Adejumobi et al. | A study of spatial media-based modulation using RF mirrors | |
CN1455524A (zh) | 用于多发射多接收(mimo)天线阵列系统的空间并行交织方法 | |
CN100361419C (zh) | 一种空时网格码的构造方法 | |
Mushtaq et al. | Novel construction methods of quaternion orthogonal designs based on complex orthogonal designs | |
Gharavi-Alkhansari et al. | Fast antenna subset selection in wireless MIMO systems | |
Gutierrez et al. | Generalized phase spatial shift keying modulation for MIMO channels | |
Patel et al. | Performance Analysis & implementation of different modulation techniques in Almouti MIMO scheme with Rayleigh channel | |
Bale et al. | Computer design of super-orthogonal space-time trellis codes | |
CN109962749A (zh) | 重叠复用系统及其处理方法和装置、存储介质、处理器 | |
Liu et al. | Error probability for orthogonal space-time block code diversity system using rectangular QAM transmission over Rayleigh fading channels | |
CN109842463A (zh) | 数据传输方法、装置、存储介质及处理器 | |
Krasicki | Labelling diversity for MIMO systems | |
Zhuang | Effects of erroneous CSI on the performance of multiple-transmit antenna selection | |
Sharma et al. | Performance of MIMO space-time coded wireless communication systems | |
Shan et al. | Signal constellations for differential unitary space-time modulation with multiple transmit antennas | |
Al-Ghadhban et al. | Iterative joint and interference nulling/cancellation decoding algorithms for multigroup space time trellis coded systems | |
Li et al. | Space-time trellis codes with linear transformation for fast fading channels |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200331 |