发明内容
发明目的:为充分实现MIMO系统多天线的分集增益,利用UW可同时作CP和训练序列的优势,本发明提出了一种SC-MIMO系统中单流四天线开环分集收发方法,此装置可实现空间扩展分集,并且不影响UW序列的数据块结构。
技术方案:为实现上述目的,本发明在系统发射端,数据比特流首先经过补零操作,然后通过调制器,将调制后的数据符号分块,形成数据符号块,再以块为单位实施空间扩展,在每根天线上的数据符号块之间插入UW序列,同一根天线上插入的UW是相同的,不同的天线间插入的UW不同,UW的长度可根据系统的信道状况选择长或短两种。在接收端,按对应发射端所选择的空间扩展分集方案来确定接收端的具体操作,此操作可通过信令字段的指示进行选择但并不局限于此种指示方式。接收端通过相应的循环移位,时频转换,频域加减特定项等操作实现正确接收。
本发明的技术方案如下:
一种SC-MIMO系统中单流四天线开环分集收发方法,包括以下步骤:
1)在发射端,对待发送数据编码比特流后进行补零操作;
2)对补零后的数据比特流进行符号调制;
3)将调制后的数据符号流均匀分块,得到连续数据符号块;
4)对连续数据符号块进行空间扩展操作;
5)对经过空间扩展后的数据符号插入UW序列;
6)在接收端,对接收天线接收到的数据符号进行去首个UW操作,得到接收符号;
7)对接收符号进行均匀分块操作,得到接收数据符号块;
8)对接收数据符号块进行循环移位操作,得到时域信号;
9)对时域信号进行时频转换,得到频域信号;
10)在频域对频域信号进行加减特定项操作;
11)分别对接收信号的每一个频点进行频域均衡,得到频域的发送符号块估计;
12)将均衡后的频域符号估计进行频域时域转换,得到时域信号估计;
13)对时域信号估计进行处理,除去每块中的UW,获取空间流数据符号块估计;
14)对空间流数据符号块进行解块操作,恢复空间流数据符号估计;
15)对空间流数据符号估计进行解调,得到接收端的数据比特估计;
16)将数据比特估计进行去零操作,在接收端得到数据比特的估计值。
本发明的一个具体技术方案如下:
在发射端:
所述补零操作是在发射端待发送的数据编码比特流b(0)b(1)…b(Nbits-1)后补上Npad个零,Nbits为数据编码比特流的长度,Npad的数值为
Npad=Ndbpb×Nblock-Nbits
其中,Ndbpb表示每个数据符号块中包含的数据信息比特数,其值为Ndbpb=(Ndspb-G)×Ndbps,符号Ndspb表示每个数据符号块中包含的调制符号数,G表示插入的UW序列的长度,符号Ndbps表示每个星座调制符号包含的数据比特数;Nblock表示该编码数据比特流被划分成的块数,其值为计算式中的符号表示向上取整,
补零后的数据编码比特流变为:
b(0)b(1)…b(Nbits-1)b(Nbits)…b(Nbits+Npad-1),其中,
b(Nbits)…b(Nbits+Npad-1)为0;
所述符号调制是对补零后的数据编码比特b(0)b(1)…b(Nbits-1)b(Nbits)…b(Nbits+Npad-1)进行符号调制,调制后的数据符号流为s(0)s(1)…s(Nsym-1),Nsym表示调制后的符号流长度;
所述的将调制后的数据符号流均匀分块是将调制后的数据符号流均匀分成Nblock个块,每个块中包含K个数据调制符号,K=Ndspb-G,用符号si=[si(0),…,si(K-1)]表示数据符号流上的第i个数据符号块向量,i=0,1,…,Nblock;
所述空间扩展是以每四个连续数据符号块si,si+1,si+2,si+3,i=0,4,8…为一个组,对其进行空间扩展操作,假设经过空间扩展后第m根天线上待发送的第i个数据块用符号表示,则对于四根发射天线的各天线上的待发送数据如下:
天线一:
天线二:
天线三:
天线四:
其中,符号表示将K维的数据符号向量si末尾补零延长至Ndspb维后的数据符号向量,即上标*表示对矩阵或向量求共轭运算,矩阵Q是一个Ndspb×Ndspb维的置换矩阵,实现逆序循环移位操作,其具体定义如下:
其中Q2=0K×G,Q3=0G×K,Q4=IG×G;其中I和0分别表示单位阵和零矩阵;
(1)所述插入UW序列是对经过空间扩展后的数据符号插入UW序列,用um表示第m根天线上插入的UW序列,um=[um,0,...,um,G-1],长度为G,在保证保护间隔持续时长大于信道的最大时延扩展的情况下,G的值有两种选择,为G1和G2(且G1>G2),G1和G2是可由通信系统事先确定的两个正整数,若当前信道的最大时延扩展小于G2长度的保护间隔持续时长,则选择G=G2,当最大时延扩展较大,即最大时延扩展大于G2长度的保护间隔持续时长,选择G=G1,但仍然需要保证保护间隔持续时长大于信道最大时延扩展。在信令字段用1个比特来指示UW的长度;插入的具体方式为:首先在第m根天线上的待发送数据符号Nblock-1的最前端插入um,然后将每个待发送数据符号块的末尾G个零值用UW序列替换,完成插入UW后的待发送信号的符号序列用符号xm表示,则
其中 的定义为 是一个长度为Ndbps的数据符号块;
假设用符号r=[r(0)r(1)…r(Nsym+G(Nblock+1)-1)]表示接收天线上接收到的数据符号,在接收端,
所述丢弃首个UW序列是去掉数据符号r的前G个接收符号,所述前G个接收符号为[r(0)r(1)…r(G-1)],得到去掉了首个UW序列的符号y;
所述对接收符号进行均匀分块操作是将符号y均匀分成Nblock个块,每个块的长度为Ndspb,并用符号yi表示第i个接收数据符号块;
所述循环移位是对接收数据符号块yi+1,yi+2,i=0,4,8,…作循环移位操作,循环移位后的信号用符号y′i+1,y′i+2表示,则
y′i+1=Pyi+1
y′i+2=Pyi+2
其中矩阵P是一个N×N维的置换矩阵,实现循环移位操作,其具体定义如下:
其中P1=0(G+1)×(K-1),P2=I(G+1)×(G+1),P3=I(K-1)×(K-1),P4=0(K-1)×(G+1)
所述时频转换是对时域信号yi,y′i+1,y′i+2,yi+3,i=0,4,8,…作Ndspb点FFT,得到频域信号,即:
Yi=Wyi
Y′i+1=Wy′i+1
Y′i+2=Wy′i+2
Yi+3=Wyi+3
其中,W是一个Ndspb×Ndspb维的DFT阵,矩阵W的第(p,q)个元素为其中J为虚数单位,表示j2=-1;
所述频域加减特定项是在频域,将Y′i+1,Y′i+2,Yi+3分别减去和补充某些特定项,以方便后续操作,减去和补充特定项之后的频域信号用符号表示,则:
上述式中,符号Λm代表第m根发射天线到接收天线的信道循环矩阵Hm的频域对角化矩阵,即Λm=WHmWH,上标H表示对矩阵作共轭转置;
所述均衡为:定义则其中[Z1 Z2 Z3 Z4]T代表噪声项,上标T代表矩阵转置。分别对接收信号的每一个频点进行频域均衡,得到频域的发送符号块估计,即X1,X2,X3,X4的估计并用符号表示;
所述频域时域转换是将均衡后的频域符号块估计做Ndspb点FFT,得到时域信号估计即:
所述除去每块中的UW序列是将估计的信号最后的G个元素用零元素替换,即可得到发射端符号的估计即抽取出的前K个元素,得到即并由发射端的空间扩展操作对应可知,其中,i=0,4,8,…
所述解块操作是将恢复出的空间流数据符号块按照块的序号从小到大排列,恢复空间流数据符号即:
所述解调是对空间流数据符号估计进行解调,得到接收端的数据比特估计
所述去零是将接收端的数据比特估计去掉最后的Npad个比特值,最终接收端得到数据比特的估计值
本发明的另一个具体技术方案如下:
在发射端:
所述补零操作是在发射端待发送数据比特流后补上Npad个零,Npad的数值为:
Npad=Ndbpb×Nblock-Nbits
其中,Ndbpb表示每个数据符号块中包含的数据信息比特数,其值为Ndbpb=(Ndspb-G)×Ndbps,符号Ndspb表示每个数据符号块中包含的调制符号数,G表示插入的UW序列的长度,符号Ndbps表示每个星座调制符号包含的数据比特数;Nblock表示该编码数据比特流被划分成的块数,其值为Nblock计算式中的符号表示向上取整,补零后的数据比特流变为:b(0)b(1)…b(Nbits-1)b(Nbits)…b(Nbits+Npad-1),其中,b(Nbits)…b(Nbits+Npad-1)为0;
所述符号调制是对补零后的比特b(0)b(1)…b(Nbits-1)b(Nbits)…b(Nbits+Npad-1)进行符号调制,调制后的数据符号流为s(0)s(1)…s(Nsym-1),Nsym表示调制后的符号流长度;
所述的将调制后的数据符号流均匀分块是将调制后的数据符号流均匀分成Nblock个块,每个块中包含K个数据调制符号,K=Ndspb-G,用符号si=[si(0),...,si(K-1)]表示数据符号流上的第i个数据符号块向量,i=0,1,…,Nblock-1,
所述空间扩展操作是以每两个连续数据符号块si,si+1,i=0,2,4,…为一个组,对其进行空间扩展操作,四根发射天线的各天线上的待发送数据如下:
天线一:
天线二:
天线三:
天线四:
其中,符号表示将K维的数据符号向量si末尾补零延长至Ndspb维后的数据符号向量,即
所述的插入UW序列是用um表示第m根天线上插入的UW序列,um=[um,0,...,um,G-1],长度为G,G的值可为G1和G2两种,并通过信令字段的1个比特来指示,插入的具体方式为:首先在第m根天线上的待发送数据符号的最前端插入um,然后将每个待发送数据符号块的末尾G个零值用UW序列替换,完成插入UW后的待发送信号的符号序列用符号xm表示,则
令 的定义为 是一个长度为Ndbps的数据符号块;
假设用符号r=[r(0)r(1)…r(Nsym+G(Nblock+1)-1)]表示接收天线上接收到的数据符号,则
所述去首个UW是对接收数据符号r,去掉前G个接收符号,即去掉符号[r(0)r(1)…r(G-1)],得到符号y;
所述对接收符号进行均匀分块操作是将y均匀分成Nblock个块,每个块的长度为Ndspb,用符号yi表示第i个接收数据符号块;
所述循环移位操作是对接收数据符号块yi+1,i=0,2,4,…,作循环移位操作,循环移位后的信号用符号y′i+1表示,则
y′i+1=Pyi+1,
所述时频转换是对时域信号yi,y′i+1,i=0,2,4,…作Ndspb点FFT,得到频域信号,即:
Yi=Wyi
Y′i+1=Wy′i+1
W表示一个Ndspb×Ndspb维的DFT阵;
所述加减特定项操作是在频域,将Yi,Y′i+1分别减去和补充某些特定项,以方便后续操作,减去和补充特定项之后的频域信号用符号表示,则:
上述式中,符号Λm代表第m根发射天线到接收天线的信道循环矩阵Hm的频域对角化矩阵,即Λm=WHmWH,符号表示补零后的UW序列,即
所述频域均衡是:定义则其中代表噪声项。分别对接收信号的每一个频点进行频域均衡就得到频域的发送符号块估计,即得到X1,X2的估计,用符号表示;
所述频域时域转换是将均衡后的频域估计符号做Ndspb点FFT,得到时域信号估计即:
所述除去每个块中的UW是将估计的信号最后的G个元素用零元素替换,即可得到发射端符号的估计即抽取出的前K个元素,得到即并由发射端的空间扩展操作对应可知, 其中,i=0,2,4,…;
所述解块操作是将上述步骤中恢复出的空间流数据符号块按照块的序号从小到大排列,恢复空间流数据符号即:
所述解调是对空间流数据符号估计进行解调,得到接收端的数据比特估计
所述去零是将上述步骤中得到的数据比特估计去掉最后的Npab个比特值,最终接收端得到数据比特的估计值
有益效果:本发明提出的基于SC-MIMO系统的单流四发射天线开环分集收发装置可在不增加系统复杂度的情况下,充分利用MIMO系统多天线的分集增益和UW同时作为CP和训练序列的优势,克服传统的空间扩展操作无法与UW相结合的缺陷,弥补传统空间扩展无法实现空间分集增益的不足,本发明所提出的空间扩展方案在实现空间扩展功能,并与UW相结合的同时可提供空间分集增益,因此本发明可以明显改善系统性能。
具体实施方式
本发明涉及单载波多输入多输出(Single Carrier Multiple Input MultipleOutput,SC-MIMO)系统的单流四天线开环分集收发装置,假设用符号b(0)b(1)…b(Nbits)表示发射端待发送的数据比特流,长度为Nbits,两种单流四天线开环分集的收发装置的具体操作步骤如下:
方法一:
发送端的操作流程如图1所示,具体步骤如下:
(1)补零操作,在发射端待发送数据编码比特流后补上Npad个零,Npad的数值为
Npad=Ndbpb×Nblock-Nbits
其中,Ndbpb表示每个数据符号块中包含的数据信息比特数,其值为Ndbpb=(Ndspb-G)×Ndbps,符号Ndspb表示每个数据符号块中包含的调制符号数,G表示插入的UW的长度,符号Ndbps表示每个星座调制符号包含的数据比特数,例如对于QPSK调制,Ndbps为2;Nblock表示该编码数据比特流被划分成的块数,其值为符号表示向上取整。
补零后的数据比特流变为:
b(0)b(1)…b(Nbits-1)b(Nbits)…b(Nbits+Npad-1),其中,
b(Nbits)…b(Nbits+Npad-1)为0
(2)符号调制。对补零后的比特b(0)b(1)…b(Nbits-1)b(Nbits)…b(Nbits+Npad-1)进行符号调制,调制后的数据符号流为s(0)s(1)…s(Nsym-1),Nsym表示调制后的符号流长度。
(3)成块。将调制后的数据符号流均匀分成Nblock个块,每个块中包含K个数据调制符号,K=Ndspb-G.用符号si=[si(0),...,si(K-1)]表示数据符号流上的第i个数据符号块向量,i=0,1,…,Nblock-1.
(4)空间扩展。以每四个连续数据符号块si,si+1,si+2,si+3,i=0,4,8…为一个组,对其进行空间扩展操作,假设经过空间扩展后第m根天线上待发送的第i个数据块用符号表示,则对于四根发射天线的各天线上的待发送数据如下:
天线一:
天线二:
天线三:
天线四:
其中,符号表示将K维的数据符号向量si末尾补零延长至Ndspb维后的数据符号向量,即上标*表示对矩阵或向量求共轭。.
矩阵Q是一个Ndspb×Ndspb维的置换矩阵,实现逆序循环移位操作,其具体定义如下:
其中Q2=0K×G,Q3=0G×K,Q4=IG×G,其中
I和0分别表示单位阵和零矩阵。
(5)插入UW(即插入独特字Unique Word,UW).对经过空间扩展后的数据符号插入UW序列,用um表示第m根天线上插入的UW序列,um=[um,0,…,um,G-1],长度为G。在保证保护间隔持续时长大于信道的最大时延扩展的情况下,G的值有两种选择,为G1和G2(且G1>G2),G1和G2是可由通信系统事先确定的两个正整数,若当前信道的最大时延扩展小于G2长度的保护间隔持续时长,则选择G=G2,当最大时延扩展较大,即最大时延扩展大于G2长度的保护间隔持续时长,选择G=G1,但仍然需要保证保护间隔持续时长大于信道最大时延扩展。可在信令字段用1个比特来指示UW的长度。插入的具体方式为:首先在第m根天线上的待发送数据符号Nblock-1的最前端插入um,然后将每个待发送数据符号块的末尾G个零值用UW序列替换,完成插入UW后的待发送信号的符号序列用符号xm表示,则
其中 的定义为 是一个长度为
Ndbps的数据符号块.发射端的发送数据格式如图2所示。
下面以接收端配置一根接收天线为例说明接收端的具体操作步骤,操作流程如图4所示。若接收端配置多根接收天线,每一根天线上的操作步骤与一根天线配置上的操作相一致。假设用符号r=[r(0)r(1)…r(Nsym+G(Nblock+1)-1)表示接收天线上接收到的数据符号,接收端的具体操作步骤如下:
(1)丢弃首个UW。对接收数据符号r,首先去掉前G个接收符号,即去掉符号[r(0)r(1)…r(G-1)],称此操作为去首个UW,去掉首个UW的接收数据符号用y来表示.
(2)接收成块操作。在接收端,对接收符号y作与发射端相一致的成块操作,将y均匀分成Nblock个块,每个块的长度为Ndspb,用符号yi表示分成的第i个接收数据符号块。
(3)循环移位。对接收数据符号块yi+1,yi+2,i=0,4,8,…,作循环移位操作,循环移位后的信号用符号y′i+1,y′i+2表示,则
y′i+1=Pyi+1
y′i+2=Pyi+2
其中矩阵P是一个N×N维的置换矩阵,实现循环移位操作,其具体定义如下:
其中P1=0(G+1)×(K-1),P2=I(G+1)×(G+1),P3=I(K-1)×(K-1),P4=0(K-1)×(G+1)
(4)时频转换。对时域信号yi,y′i-1,y′i+2,yi+3,i=0,4,8,…作Ndspb点FFT,得到频域信号,即:
Yi=Wyi
Y′i+1=Wy′i+1
Y′i+2=Wy′i+2
Yi+3=Wyi+3
其中,W是一个Ndspb×Ndspd维的DFT阵,矩阵W的第(p,q)个元素为其中J为虚数单位,即j2=-1.
(5)频域加减特定项。在频域,将Y′i+1,Y′i+2,Yi+3分别减去和补充某些特定项,以方便后续操作,减去和补充特定项之后的频域信号用符号表示,则:
上述式中,符号Λm代表第m根发射天线到接收天线的信道循环矩阵Hm的频域对角化矩阵,即ΛmWHmWH,,上标H表示对矩阵作共轭转置,符号表示补零后的UW序列,即
(6)均衡.定义则其中代表噪声项,上标T代表矩阵转置。分别对接收信号的每一个频点进行频域均衡就可以得到频域的发送符号块估计,即得到X1,X2,X3,X4的估计,用符号 表示。
(7)频域时域转换。将均衡后的频域符号估计做Ndspb点FFT,得到时域信号估计即:
(8)除去每块中的UW.将估计的信号最后的G个元素用零元素替换,即可得到发射端符号的估计即抽取出的前K个元素,得到即并由发射端的空间扩展操作对应可知, 其中,i=0,4,8,…
(9)解块操作。将上述步骤中恢复出的空间流数据符号块按照块的序号从小到大排列,恢复空间流数据符号即:
(10)解调。对空间流数据符号估计进行解调,得到接收端的数据比特估计
(11)去零。将上述步骤中得到的数据比特估计去掉最后的Npad个比特值,最终接收端得到数据比特的估计值
方法二:
发送端的具体操作步骤如下:
(1)补零操作,在发射端待发送数据比特流后补上Npad个零,Npad的数值为
Npad=Ndbpb×Nblock-Nbits
其中,Ndbpb表示的意义与方案一中相同;Nblock同样表示数据比特流被划分成的块数,与方案一中的值不同,其值为补零后的数据比特流变为:b(0)b(1)…b(Nbits-1)b(Nbits)…b(Nbits+Npad-1),其中,b(Nbits)…b(Nbits+Npad-1)为0
(2)符号调制。补零后的比特b(0)b(1)…b(Nbits-1)b(Nbits)…b(Nbits+Npad-1)进行符号调制,调制后的数据符号流为s(0)s(1)…s(Nsym-1),Nsym表示调制后的符号流长度。
(3)成块。将调制后的数据符号流均匀分成Nblock个块,每个块中包含K个数据调制符号,K=Ndspb-G.用符号si=[si(0),…,si(K-1)]表示数据符号流上的第i个数据符号块向量,i=0,1,…,Nblock-1.
(4)空间扩展。以每两个连续数据符号块si,si+1,i=0,2,4,…为一个组,对其进行空间扩展操作.则对于四根发射天线的各天线上的待发送数据如下:
天线一:
天线二:
天线三:
天线四:
Q的定义与方案一中相同
(5)插入UW.对经过空间扩展后的数据符号插入UW序列,用um表示第m根天线上插入的UW序列,um=[um,0,...,um,G-1],长度为G,与方案一相同,G的值可为G1和G2两种,并通过信令字段的1个比特来指示。插入的具体方式为:首先在第m根天线上的待发送数据符号的最前端插入um,然后将每个待发送数据符号块的末尾G个零值用UW序列替换,完成插入UW后的待发送信号的符号序列用符号xm表示,则
其中 的定义为 是一个长度为
Ndbps的数据符号块.发送数据格式如图3所示。
下面以接收端配置一根接收天线为例说明接收端的具体操作步骤。若接收端配置多根接收天线,每一根天线上的操作步骤与一根天线配置上的操作相一致。假设用符号r=r(0)r(1)…r(Nsym+G(Nblock+1)-1)]表示接收天线上接收到的数据符号,接收端的具体操作步骤如下:
(1)丢弃首个UW.对接收数据符号r,去掉前G个接收符号,即去掉符号[r(0)r(1)…r(G-1),去掉首个UW的数据符号用y来表示.
(2)接收成块操作。在接收端,对接收符号y作与发射端相一致的成块操作,将y均匀分成Nblock个块,每个块的长度为Ndspb,用符号yi表示分成的第i个接收数据符号块。
(3)循环移位。对接收数据符号块yi+1,i=0,2,4,…,作循环移位操作,循环移位后的信号用符号y′i+1表示,则
y′i+1=Pyi+1
P的定义与方案一相同。
(4)时频转换。对时域信号yi,y′i+1,i=0,2,4,…作Ndspb点FFT,得到频域信号,即:
Yi=Wyi
Y′i+1=Wy′i+1
与方案一中相同,W表示一个Ndspb×Ndspb维的DFT阵.
(5)频域加减特定项。在频域,将Yi,Y′i+1分别减去和补充某些特定项,以方便后续操作,减去和补充特定项之后的频域信号用符号表示,则:
上述式中,符号Λm代表第m根发射天线到接收天线的信道循环矩阵Hm的频域对角化矩阵,即Λm=WHmWH,符号表示补零后的UW序列,即
(6)均衡.定义 则其中代表噪声项。分别对接收信号的每一个频点进行频域均衡就可以得到频域的发送符号块估计,即得到X1,X2的估计,用符号表示。
(7)频域时域转换。将均衡后的频域估计符号做Ndspb点FFT,得到时域信号估计即:
(8)除去每块中的UW.将估计的信号最后的G个元素用零元素替换,即可得到发射端符号的估计即抽取出的前K个元素,得到即由发射端的空间扩展操作对应可知, 其中,i=0,2,4,…
(9)解块操作。将上述步骤中恢复出的空间流数据符号块按照块的序号从小到大排列,恢复空间流数据符号即:
(10)解调。对空间流数据符号估计进行解调,得到接收端的数据比特估计
(11)去零。将上述步骤中得到的数据比特估计去掉最后的Npad个比特值,最终接收端得到数据比特的估计值