CN105743618B - 一种基于qr分解的混合空间调制天线选择方法 - Google Patents
一种基于qr分解的混合空间调制天线选择方法 Download PDFInfo
- Publication number
- CN105743618B CN105743618B CN201410751818.7A CN201410751818A CN105743618B CN 105743618 B CN105743618 B CN 105743618B CN 201410751818 A CN201410751818 A CN 201410751818A CN 105743618 B CN105743618 B CN 105743618B
- Authority
- CN
- China
- Prior art keywords
- antenna
- antenna combination
- selection method
- sequence numbers
- combination
- 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
- 238000000034 method Methods 0.000 title claims abstract description 15
- 238000002156 mixing Methods 0.000 title abstract 2
- 238000010187 selection method Methods 0.000 claims abstract description 28
- 239000011159 matrix material Substances 0.000 claims abstract description 21
- 238000000354 decomposition reaction Methods 0.000 claims description 17
- 238000010845 search algorithm Methods 0.000 claims description 10
- 238000010586 diagram Methods 0.000 claims description 7
- 230000005540 biological transmission Effects 0.000 abstract description 6
- 238000001228 spectrum Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Landscapes
- Radio Transmission System (AREA)
Abstract
本发明公开了一种基于QR分解的混合空间调制天线选择方法,包括:采用穷搜索算法实现最佳天线组合的选择方法;Φ(l)表示对于给定天线组合l所有可能天线序号的集合;计算天线组合l对应的最小欧式距离,根据所述最小欧式距离以及发送数据的天线序号的集合重新确定最佳天线组合的选择方法;定义一个上三角矩阵使得D(l)中非零元素的最小值等价于天线组合l对应的最小欧式距离;分别对所述上三角矩阵D(l)进行简化处理;得到简化以后的最佳天线组合的选择方法。本发明有益效果:在保证系统性能基本不变的情况下,大大降低了算法复杂度,减小了算法运行时间。
Description
技术领域
本发明属于无线通信调制与编码领域,特别涉及一种MIMO系统中基于QR分解的混合空间调制天线选择方法。
背景技术
虽然传统空间调制(SM)技术完全避免了信道间干扰,不需要发射天线间较高的同步性要求,但它要求发射天线总数必须是2的幂次方。而后提出的扩展的空间调制方案(ESM)虽然允许所有天线在同一时隙内发送数据,但这势必引起天线间的干扰。而混合空间调制(HSM)技术通过不断在SM与ESM之间切换的方式,结合了SM的ESM的优点,避免了其各自的缺点。同时,将天线选择应用到HSM技术上,特别是针对基于欧式距离的混合空间调制天线选择算法HSM-EDAS(Hybrid Spatial Modulation based on Euclidean Distanceoptimized Antenna Selection),大大降低了系统的误比特率。
然而,当发射天线总数较多、星座图较大时,传统的天线选择算法,复杂度普遍较高。比如穷搜索算法,穷搜索算法在接收端进行检测时要对所有的可能的发送天线与所有的星座符号都进行一次遍历,增加了算法复杂度以及算法运行时间。
发明内容
本发明的目的就是为了解决上述问题,提出了一种基于QR分解的混合空间调制天线选择方法,该方法在保证系统性能基本不变的情况下,大大降低了选择算法复杂度,减小了算法运行时间。
为了实现上述目的,本发明采用如下技术方案:
一种基于QR分解的混合空间调制天线选择方法,包括以下步骤:
(1)假设发送端总共有Ns根发送天线,从中选择Nt根发送数据,接收端有Nr根接收天线,则接收信号y=Hx+n,采用穷搜索算法实现最佳天线组合的选择方法;
(2)定义表示发送符号xi、xj中非零元素所在位置的序号的集合,即发送数据的天线序号的集合,Φ(l)表示对于给定天线组合l所有可能天线序号的集合;计算天线组合l对应的最小欧式距离,根据所述最小欧式距离以及发送数据的天线序号的集合重新确定最佳天线组合的选择方法;
(3)定义一个上三角矩阵其中,使得D(l)中非零元素的最小值等价于天线组合l对应的最小欧式距离;
(4)对于上三角矩阵D(l),在i=j和i≠j时,分别对所述上三角矩阵D(l)进行简化处理;得到简化以后的最佳天线组合的选择方法。
所述步骤(1)中采用穷搜索算法实现最佳天线组合的选择方法具体为;
其中,xi=[0,...,s1,0,...,s1,0,...]T,表示发
送端发送符号,Λ是所有可能发送符号的集合,表示天线组
合l对应的Nt列服从独立同分布的标准高斯信道,Γ表示所有种可能的天线组合。
所述步骤(2)中,天线组合l对应的最小欧式距离为:
其中,表示发送符号xi、xj中非零元素所在位置的序号的集合,即发送数据的天线序号的集合,Φ(l)表示对于给定天线组合l所有可能天线序号的集合。s1,s2∈S,S是所有可分解M维正交幅度调制星座图所有符号的集合,表示天线组合l中对应第p或q根发送天线的信道向量。
所述步骤(2)中,根据所述最小欧式距离以及发送数据的天线序号的集合重新确定最佳天线组合的选择方法具体为:
其中,为天线组合l对应的最小欧式距离,xi,xj表示为xi=[0,...,s1,0,...,s1,0,...]T,xj=[0,...,s2,0,...,s2,0,...]T,其中s1,S是M维的正交幅度调制星座图所有符号的集合。
所述步骤(3)中的上三角矩阵D(l)具体为:
其中,表示发送符号xi、xj中非零元素所在位置的序号的集合,即发送数据的天线序号的集合,Φ(l)表示对于给定天线组合l所有可能天线序号的集合。s1,s2∈S,S是所有可分解M维正交幅度调制星座图所有符号的集合,表示天线组合l中对应第p或q根发送天线的信道向量。
所述步骤(4)中,当i=j时,对所述上三角矩阵D(l)进行简化处理的具体方法为:
当i=j时,要求s1≠s2,并且对于QAM调制,任意两个星座点之间的最小距离可表示为此时,公式(4.a)简化为:
其中,表示发送符号xi中非零元素所在位置的序号的集合,即发送数据的天线序号的集合,Φ(l)表示对于给定天线组合l所有可能天线序号的集合;表示天线组合l中对应第p根发送天线的信道向量。
所述步骤(4)中,当i≠j时,对所述上三角矩阵D(l)进行简化处理的具体方法为:将公式(4.b)表示为:
其中对做QR分解得到矩阵Q和R,其
中,QHQ=I|Q|且R是一个(4×4)的上三角矩阵,则公式(4.b)初步简化为:
其中,siI和siQ分别表示si的实部和虚部,且M=N1N2,i=1,2;
由QAM星座点的对称性可知,因此公式(4.b)最终简化为:
其中,一个可分解M维正交幅度调制星座图可以分解为两个N1维和N2维的脉冲幅度调制星座图,M=N1N2,M1和M2分别代表N1维和N2维的脉冲幅度调制星座图x轴以上的部分。
所述步骤(4)中简化以后的最佳天线组合的选择方法为:
其中,
其中,表示发送符号xi中非零元素所在位置的序号的集合,即发送数据的天线序号的集合,Φ(l)表示对于给定天线组合l所有可能天线序号的集合。表示天线组合l中对应第p根发送天线的信道向量,siI和siQ分别表示si的实部和虚部,i=1,2,M1和M2分别代表N1维和N2维的脉冲幅度调制星座图x轴以上的部分。
本发明的有益效果是:
本发明提出的基于QR分解的欧氏距离混合空间调制天线选择算法,避免了当发射天线总数较多、星座图较大时,传统的天线选择算法(比如穷搜索算法)需要对所有天线与符号组合都进行一次遍历的缺点,在保证系统性能基本不变的情况下,大大降低了算法复杂度,减小了算法运行时间。
附图说明
图1是本发明的基于QR分解的混合空间调制天线选择的系统架构图;
图2是频谱效率为3bit/s/Hz,从4根发送天线选择2根进行发送时本发明混合空间调制方案与传统SM、ESM方案的比较示意图;
图3是频谱效率为6bit/s/Hz,从6根发送天线选择4根进行发送时本发明混合空间调制方案与传统SM、ESM方案的比较示意图。
具体实施方式:
下面结合附图与实施例对本发明做进一步说明:
假设发送端总共有Ns根发送天线,从中选择Nt根发送数据,接收端有Nr根接收天线则接收信号y=Hx+n,系统框图如图1,假设信道变化缓慢且频谱效率一定,HSM-EDAS方案信息传输过程为:在每次信息发送前,接收端检测信道状态,分别计算出采用SM和ESM方案发送数据时的最小平方欧式距离,选出最小平方欧式距离较大的方案用来实际发送符号,接收端采用最大似然进行检测。则穷搜索算法实现最佳天线组合的选择方法如下表示
其中表示发送端发送符号,Λ是所有可能发送符号的集合,表示天线组合l对应的Nt列服从独立同分布的标准高斯信道,
Γ表示所有种可能的天线组合。
由上述HSM-EDAS方案中ESM发送符号的特征可知,xi,xj可表示为xi=[0,...,s1,0,...,s1,0,...]T,xj=[0,...,s2,0,...,s2,0,...]T,其中s1,(S是M-QAM所有符号的集合)。我们定义表示发送符号xi、xj中非零元素所在位置的序号的集合(即发送数据的天线序号的集合),Φ(l)表示对于给定天线组合l所有可能天线序号的集合。例如,xi=[+10],xj=[+1+1],则所以,上式(1)可表示为
其天线组合l对应的最小欧式距离为
从公式(2)可以得出计算lselected的复杂度为为了降低基于欧式距离的混合空间调制方案(HSM-EDAS)的复杂度,本发明提出一种简化公式(2)的算法RC-HSM-EDAS。
定义一个上三角矩阵(其中),它的第i行j列元素表示为
则公式(2)等价于
其中,minD(l)表示D(l)中非零元素的最小值。
当i≠j时,上式(4.b)可表示为其中 对做QR分解得到矩阵Q和R,其中,QHQ=I|Q|且R是一个(4×4)的上三角矩阵。则公式(4.b)可重写为:
siI和siQ分别表示si(i=1,2)的实部和虚部且M=N1N2。因此,公式(4)可简化为
考虑到M维正交幅度调制M-QAM星座图的对称性,公式(7)可被进一步简化。
首先,对于公式(7.a),要求s1≠s2。而对于QAM调制,任意两个星座点之间的最小距离可表示为因此,公式(7.a)可改写为:
其次,对于公式(7.b),对所有的星座点都要计算矩阵且s1I,s2I∈N1-
PAM,s1Q,s2Q∈N2-PAM。然而,由QAM星座点的对称性可知因此公式(7.b)等价于
其中,一个可分解M维正交幅度调制星座图可以分解为两个N1维和N2维的脉冲幅度
调制N1-PAM和N2-PAM,M=N1N2,M1和M2分别代表N1维和N2维的脉冲幅度调制星座图x轴以上的
部分。因此,当i=j时D(l)第i行j列的计算复杂度从降低到O(1);当i≠j
时从O(N1N2)降低到
频谱效率为3bit/s/Hz,从4根发送天线中选择2根进行发送时不同天线选择方法的比较示意图如图2所示。频谱效率为6bit/s/Hz,从6根发送天线选择4根进行发送时不同天线选择方法的比较示意图如图3所示。
可以看出,混合空间调制方法的误比特率低于传统的空间调制(SM)与后来扩展的空间调制(ESM)方法,且将传统天线选择(穷搜索算法)应用到混合空间调制(HSM-EDAS)后使系统性能得到较大提高。而本发明提出的基于QR分解的欧氏距离混合空间调制方法(RC-HSM-EDAS)保持了原穷搜索算法低误比特率的优点,且降低了算法时间复杂度。
采用本发明的天线选择方法(RC-HSM-EDAS)与基于欧式距离的混合空间调制方案(HSM-EDAS)的复杂度比较如表1所示。
表1
实验结果标明,本发明提出的基于QR分解的欧氏距离混合空间调制天线选择算法,避免了当发射天线总数较多、星座图较大时,传统的天线选择算法(比如穷搜索算法)需要对所有天线与符号组合都进行一次遍历的缺点,在保证系统性能不变的情况下,大大降低了算法复杂度,减小了算法运行时间。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
Claims (7)
1.一种基于QR分解的混合空间调制天线选择方法,其特征是,包括以下步骤:
(1)假设发送端总共有Ns根发送天线,从中选择Nt根发送数据,接收端有Nr根接收天线,则接收信号y=Hx+n,采用穷搜索算法实现最佳天线组合的选择方法;
(2)定义表示发送符号xi、xj中非零元素所在位置的序号的集合,即发送数据的天线序号的集合,Φ(l)表示对于给定天线组合l所有可能天线序号的集合;计算天线组合l对应的最小欧式距离,根据所述最小欧式距离以及发送数据的天线序号的集合重新确定最佳天线组合的选择方法;
根据最小欧式距离以及发送数据的天线序号的集合重新确定最佳天线组合的选择方法具体为:
其中,为天线组合l对应的最小欧式距离,xi,xj表示为xi=[0,...,s1,0,...,s1,0,...]T,xj=[0,...,s2,0,...,s2,0,...]T,其中s1,S是M维的正交幅度调制星座图所有符号的集合;
(3)定义一个上三角矩阵其中,使得D(l)中非零元素的最小值等价于天线组合l对应的最小欧式距离;
(4)对于上三角矩阵D(l),在i=j和i≠j时,分别对所述上三角矩阵D(l)进行简化处理;得到简化以后的最佳天线组合的选择方法。
2.如权利要求1所述的一种基于QR分解的混合空间调制天线选择方法,其特征是,所述步骤(1)中采用穷搜索算法实现最佳天线组合的选择方法具体为;
其中,xi=[0,...,s1,0,...,s1,0,...]T,表示发送端发送符号,Λ是所有可能发送符号的集合,表示天线组合l对应的Nt列服从独立同分布的标准高斯信道,Γ表示所有种可能的天线组合。
3.如权利要求1所述的一种基于QR分解的混合空间调制天线选择方法,其特征是,所述步骤(2)中,天线组合l对应的最小欧式距离为:
其中,表示发送符号xi、xj中非零元素所在位置的序号的集合,即发送数据的天线序号的集合,Φ(l)表示对于给定天线组合l所有可能天线序号的集合;表示天线组合l中对应第p或q根发送天线的信道向量。
4.如权利要求1所述的一种基于QR分解的混合空间调制天线选择方法,其特征是,所述步骤(3)中的上三角矩阵D(l)具体为:
其中,表示发送符号xi、xj中非零元素所在位置的序号的集合,即发送数据的天线序号的集合,Φ(l)表示对于给定天线组合l所有可能天线序号的集合;表示天线组合l中对应第p或q根发送天线的信道向量。
5.如权利要求4所述的一种基于QR分解的混合空间调制天线选择方法,其特征是,所述步骤(4)中,当i=j时,对所述上三角矩阵D(l)进行简化处理的具体方法为:
当i=j时,要求s1≠s2,并且对于QAM调制,任意两个星座点之间的最小距离可表示为此时,公式(4.a)简化为:
其中,表示发送符号xi中非零元素所在位置的序号的集合,即发送数据的天线序号的集合,Φ(l)表示对于给定天线组合l所有可能天线序号的集合;表示天线组合l中对应第p根发送天线的信道向量。
6.如权利要求4所述的一种基于QR分解的混合空间调制天线选择方法,其特征是,所述步骤(4)中,当i≠j时,对所述上三角矩阵D(l)进行简化处理的具体方法为:
将公式(4.b)表示为:
其中对做QR分解得到矩阵Q和R,其中,QHQ=I|Q|且R是一个(4×4)的上三角矩阵,则公式(4.b)初步简化为:
其中,siI和siQ分别表示si的实部和虚部,且M=N1N2,i=1,2;
由QAM星座点的对称性可知,因此公式(4.b)最终简化为:
其中,一个可分解M维正交幅度调制星座图可以分解为两个N1维和N2维的脉冲幅度调制星座图,M=N1N2,M1和M2分别代表N1维和N2维的脉冲幅度调制星座图x轴以上的部分。
7.如权利要求1所述的一种基于QR分解的混合空间调制天线选择方法,其特征是,所述步骤(4)中简化以后的最佳天线组合的选择方法为:
其中,
其中,表示发送符号xi中非零元素所在位置的序号的集合,即发送数据的天线序号的集合,Φ(l)表示对于给定天线组合l所有可能天线序号的集合;表示天线组合l中对应第p根发送天线的信道向量,siI和siQ分别表示si的实部和虚部,i=1,2,M1和M2分别代表N1维和N2维的脉冲幅度调制星座图x轴以上的部分。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410751818.7A CN105743618B (zh) | 2014-12-09 | 2014-12-09 | 一种基于qr分解的混合空间调制天线选择方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410751818.7A CN105743618B (zh) | 2014-12-09 | 2014-12-09 | 一种基于qr分解的混合空间调制天线选择方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105743618A CN105743618A (zh) | 2016-07-06 |
CN105743618B true CN105743618B (zh) | 2018-11-13 |
Family
ID=56239826
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410751818.7A Expired - Fee Related CN105743618B (zh) | 2014-12-09 | 2014-12-09 | 一种基于qr分解的混合空间调制天线选择方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105743618B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106911429A (zh) * | 2017-02-27 | 2017-06-30 | 电子科技大学 | 用于gsm通信系统的信号检测方法 |
CN110417448B (zh) * | 2018-04-27 | 2022-04-15 | 中兴通讯股份有限公司 | 选择天线的方法、设备及计算机可读存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102946267A (zh) * | 2012-08-31 | 2013-02-27 | 电子科技大学 | 一种基于qr分解的低复杂度mimo检测算法 |
CN103795503A (zh) * | 2014-02-28 | 2014-05-14 | 电子科技大学 | 一种qr分解检测方法 |
CN104022810A (zh) * | 2014-06-12 | 2014-09-03 | 电子科技大学 | 一种基于sm系统的发射端天线选择方法 |
-
2014
- 2014-12-09 CN CN201410751818.7A patent/CN105743618B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102946267A (zh) * | 2012-08-31 | 2013-02-27 | 电子科技大学 | 一种基于qr分解的低复杂度mimo检测算法 |
CN103795503A (zh) * | 2014-02-28 | 2014-05-14 | 电子科技大学 | 一种qr分解检测方法 |
CN104022810A (zh) * | 2014-06-12 | 2014-09-03 | 电子科技大学 | 一种基于sm系统的发射端天线选择方法 |
Non-Patent Citations (2)
Title |
---|
Antenna Selection in Spatial Modulation Systems;Rakshith Rajashekar;《IEEE COMMUNICATIONS LETTERS》;20130331;第17卷(第3期);第522页 * |
Reduced-Complexity Antenna Selection Schemes in Spatial Modulation;Zhiqiang Zhou,Ning Ge,Xiaokang Lin;《IEEE COMMUNICATIONS LETTERS》;20140131;第18卷(第1期);第14-16页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105743618A (zh) | 2016-07-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9634747B2 (en) | Apparatus and method of processing signal, and recording medium | |
CN104868944B (zh) | 一种自适应广义空间调制方法及其简化方法 | |
CN105162506B (zh) | 一种低复杂度的大规模天线选择方法 | |
CN105656529A (zh) | 基于广义空间调制系统的自适应调制方法 | |
CN104935370B (zh) | 一种用于mimo通信系统的空时联合调制的传输方法 | |
US11716131B2 (en) | Single input single output (SISO) physical layer key exchange | |
US11258649B2 (en) | Transmitter and receiver communication apparatus for non-coherent communication | |
CN106301496B (zh) | 基于天线选择和预编码的空间调制系统 | |
CN106788626B (zh) | 一种能够获得二阶发射分集的改进正交空间调制传输方法 | |
CN106209187A (zh) | 一种采用联合优化的自适应广义空间调制方法 | |
CN109167649B (zh) | 一种gsm-mbm系统低复杂度检测方法 | |
CN108988923A (zh) | 安全空间调制系统中基于信号泄露的天线选择方法 | |
CN105743618B (zh) | 一种基于qr分解的混合空间调制天线选择方法 | |
CN112929057A (zh) | 一种双重广义空间调制方法及系统 | |
CN111010220B (zh) | 基于能量效率的多用户多流的下行混合预编码方法及系统 | |
CN109831252B (zh) | 一种可见光通信中的多用户发送预编码设计方法 | |
Purwita et al. | Generalized time slot index modulation for LiFi | |
CN107493123B (zh) | 基于预编码辅助广义正交空间调制的低复杂度检测方法 | |
EP2269352B1 (en) | Techniques for multiple-subcarrier joint precoding | |
CN114124637B (zh) | 一种适用于降低ofdm系统峰均比的低复杂度方法 | |
CN106254039B (zh) | 一种广义空间调制系统的天线选择方法 | |
CN106788630B (zh) | 一种基于降低实时错误比特概率上界的功率分配方法 | |
CN107528807B (zh) | 一种基于相位判决的空间调制mqam信号检测方法 | |
EP4097645A1 (en) | Finite resolution decomposition of a matrix and matrix-vector multiplication | |
CN108809376B (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 | ||
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: 20181113 Termination date: 20191209 |