本发明涉及无线通信技术领域,更具体地,涉及一种数据传输方法。
背景技术:
滤波器组多载波系统(fbmc)作为未来无线通信的一个关键技术,是通信系统研究中的一大热点。相较于正交频分复用技术(ofdm),fbmc使用原型滤波器代替矩形窗,不同于使用矩形窗时信号的突变,经过原型滤波器的信号存在滚降过程,是平缓变化的,因此使用fbmc技术时信号会具有更好的频谱特性。然而作为多载波调制的一种,fbmc系统中的发送信号也是由各个载波分量叠加产生,因而难以避免峰均功率比(papr)过高的问题。现有的降低papr的方法中,dft-s-fbmc技术实现简单,不存在处理时延问题,但是因为在操作过程中对于每个符号段是独立处理,符号段相加时会产生新的峰值信号,所以papr无法显著降低,fb-s-fbmc技术虽然能够显著降低papr,但是实现复杂且处理时延问题比较严重。
技术实现要素:
针对现有技术的缺陷,提出了采用相邻符号段之间合作发送数据的方法来降低papr,提出了一种主要的合作方法以及该方法的两种特殊变形。
为实现上述目的,本发明提供一种数据传输方法,包括:将待发送的数据符号划分成数据符号段并编序号,并对待发送的数据符号段按照其序号进行处理,包括加窗、fft,并采用滤波器组技术将处理过的数据符号映射到多个子载波上,所述对待发送的数据符号段按照其序号处理使得数据符号段映射到多个子载波上后,对于偶数序号数据符号段,偶数序号子载波传输实数部分,奇数序号子载波传输虚数部分,对于奇数序号数据符号段,偶数序号子载波传输虚数部分,奇数序号子载波传输实数部分,实现相邻数据符号段之间合作发送数据。
可选地,假设子载波总数为m,滤波器组的符号周期为t,t时刻内发送ms个数据符号,将待发送的数据符号序列记为{a0,a1,a2,...},所述待发送的数据符号序列要经过长度为ms的滑动窗口,每隔ms/2点滑动截取数据符号,得到在第n个数据符号段xn上的符号为::
其中,每个序号的数据符号段的长度为ms,第n个数据符号段xn的第k个元素记为xn[k]。
可选地,所述对待发送的数据符号段按照其序号处理包括:
用w表示长度为ms的窗矢量,对第n个数据符号段xn加窗后得到rn:
rn=diag(w)·xn
其中,diag函数的作用是将矢量转换成对角矩阵形式,符号·表示矩阵的乘法运算,第n个数据符号段xn的第k个元素xn[k]加窗后对应rn的第k个元素rn(k),通过ms点离散傅里叶变换将rn转换到频域形式rn:
rn=fft(rn)
fft函数表示傅里叶变换,将rn第m个元素的取值记为rn[m]。对于第n个数据符号段,将经过进一步调制后得到的en中第m个元素记为en[m],当n为偶数可以表示为:
当n为奇数时,en[m]可以表示为:
其中,
将ms个元素en[m]映射到m个子载波上,第m个子载波上的信息为in[m];
in[m]通过ifft转换到时域后与长度为km的滤波器相乘,第i个时刻的发送信号sn[i]可以表示为:
其中,g[i]表示滤波器的冲激响应,长度为km,
将所有经过上述步骤得到的sn[i]进行依次位移
可选地,可供选择的窗包括但不限于w1,w2,它们的第k个元素w1[k]和w2[k]可以分别表示为:
可选地,将ms个元素en[m]映射到m个子载波上,第m个子载波上的信息in[m]可以表示为:
或者:
或者:
可选地,所述对待发送的数据符号段按照其序号处理包括:
将xn拆分为an,bn,cn,dn四个子段,满足an和bn的长度为
当n为偶数时的构造方式为:
当n为奇数时构造方式为:
定义:
其中,
通过ms点离散傅里叶变换将en转换到频域形式en:
en=fft(rn)
fft函数表示傅里叶变换,将第m个元素的取值记为en[m];
将ms个元素en[m]映射到m个子载波上,第i个时刻的第m个子载波上的信息为in[m];
in[m]通过ifft转换到时域后与长度为km的滤波器相乘,第i个时刻的发送信号sn[i]可以表示为:
其中,g[i]表示滤波器的冲激响应,长度为km,
将所有经过上述步骤得到的sn[i]进行依次位移
可选地,所述对待发送的数据符号段按照其序号处理包括:
将实虚交错的符号序列xn拆分为an,bn,cn,dn四个子段,满足an和bn的长度为
对于偶数数据符号段,
对于奇数数据符号段:
可选地,将ms个元素en[m]映射到m个子载波上,第m个子载波上的信息in[m]可以表示为:
或者:
或者:
总体而言,通过本发明所构思的以上技术方案与现有技术相比,具有以下有益效果:
由于滤波器的滚降特性,其下降过程具有缓慢、平滑的特点,符号段之间相互重叠,若对符号段进行独立处理,相加传输时会产生新的峰值信号,所设计的方法采用相邻符号段之间合作发送数据的方法降低papr,不会产生新的峰值信号,同时处理复杂度和延迟较低。本发明提出了一种主要方法以及该方法的两种变形。
附图说明
图1为本发明提供的数据传输方法所用发射机的具体实现框架示意图;
图2为本发明提供的对数据符号进行处理所用的两种窗矢量函数曲线示意图,图2a为窗矢量w1对应的函数曲线,图2b为窗矢量w2对应的函数曲线;
图3为本发明提供的数据传输方法与dft-s-fbmc对应的发送信号papr性能曲线。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
dft-s-fbmc是近些年提出的一种有效降低发送信号峰均功率比的技术,核心思想是将符号序列先经过一个dft变换转换到频域后再输入到fbmc通信系统的发送端。本发明在该技术的基础上进行了改进,添加了对符号序列的预处理,从而进一步降低发送信号的峰均功率比。
本发明提出了一种低峰均功率比(papr)的数据发送接收技术,是基于dft-s-fbmc技术的一种变形技术,将待发送数据分段、加窗后再进行dft,映射到载波上后,对于偶数序号数据符号段,偶数序号载波传输实数部分,奇数序号载波传输虚数部分,奇数序号反之。
本发明的核心思想是在载波上发送频域符号,且数据段之前载波之间发送的数据是实虚交错的,并且通过数据段之间的滑动重叠以及窗函数的作用,可以实现相邻数据符号段之间合作发送数据。
在一个示例中,发射机的具体实现框图如图1所示:信号经过滑动窗口分段截取后加窗,再通过fft转换到频域,进行实虚交错提取后映射到载波上,再进行基于滤波器组的调制系统发送端的一系列调制步骤(包括ifft、滤波器组滤波、错位相加等)。
为了实现上述目的,本发明提出了一种数据传输方法,假设子载波总数为m,fbmc的复数符号周期为t,t时刻内发送ms个数据符号,将待发送的数据符号序列记为{a0,a1,a2,...}。
在一个示例中,本发明提供的数据传输方法发送过程包括了以下步骤:
(1)子载波总数为m,滤波器组的符号周期为t,t时刻内发送ms个数据符号,将待发送的数据符号序列记为{a0,a1,a2,...},待发送的数据符号序列要经过长度为ms的滑动窗口,每隔ms/2点滑动截取数据符号,得到在数据符号段上的符号为:将待发送的数据符号序列分为一个个长度为ms的数据符号段{x0,x1,x2,...},其中第n个数据符号段xn的第k个元素记为xn[k]:
其中,每个序号的数据符号段的长度为ms,第n个数据符号段xn的第k个元素记为xn[k]。
对待发送的数据符号按照符号段编序号,每个数据段包含个ms符号,且相邻数据段之间重叠ms/2个符号,具体如下所示:
(2)用w表示长度为ms的窗矢量,对第n个数据符号段xn加窗后得到rn:
rn=diag(w)·xn
其中,diag函数的作用是将矢量转换成对角矩阵形式,符号·表示矩阵的乘法运算,第n个数据符号段xn的第k个元素xn[k]加窗后对应rn的第k个元素rn(k),
将可供选择的两种窗矢量w记为w1,w2,它们的第k个元素w1[k]和w2[k]可以分别表示为:
具体地,上述两种窗矢量函数曲线如图2所示,其中图2a为窗矢量w1对应的函数曲线,图2b为窗矢量w2对应的函数曲线。
(3)通过ms点离散傅里叶变换将rn转换到频域形式rn:
rn=fft(rn)
fft函数表示傅里叶变换,将第m个元素的取值记为rn[m]。
(4)对于第n个数据符号段,将经过进一步调制后得到的en中第m个元素记为en[m],当n为偶数可以表示为:
当n为奇数时,en[m]可以表示为:
其中,
(5)将ms个元素en[m]映射到m个子载波上,第m个子载波上的信息in[m]可以表示为:
或者:
或者:
(6)in[m]通过ifft转换到时域后与长度为km的滤波器相乘,第i个时刻的发送信号sn[i]可以表示为:
其中,g[i]表示滤波器的冲激响应,长度为km,
(7)将所有经过上述步骤得到的sn[i]进行依次位移
其中,位移
在另一个示例中,本发明提供的数据传输方法还可以有如下两种变形形式:
变形1:
将待传输复数符号序列依次分为长度为ms的数据符号段{x0,x1,x2,...},将第n个数据符号段xn经过拆分,再构造后,得到新的长度为ms的序列rn,将rn经过ms点dft转换到频域后映射到载波上,然后进行采样、ifft、滤波、位移相加,得到最终的发送信号,具体可包括如下步骤:
(1)得到数据符号段{x0,x1,x2,...},每段长度为
其中,这里的数据符号分分断方式如下所示:
(2)将xn拆分为an,bn,cn,dn四个子段,满足an和bn的长度为
(3)分为n为偶数、奇数两种情况对数据符号段进行再构造,得到长度为ms的en,将其中的第k个元素记为en[k],当n为偶数时的构造方式为:
当n为奇数时构造方式为:
定义:
其中,
(4)通过ms点离散傅里叶变换将en转换到频域形式en:
en=fft(en)
fft函数表示傅里叶变换,将第m个元素的取值记为en[m]。
(5)与之前的操作一致。
(6)与之前的操作一致。
(7)与之前的操作一致。
需要说明的是,该变形方案相比于之前的方案,不同之处在于在fft之前就对序列进行再构造,使得fft的结果直接满足实虚交错的条件不需要再进行实部和虚部的提取。
变形2
下面,讨论一种特例情况,即发送的符号序列是实虚交错的,那么在拆分和再构造的时候,需要做出一些调整:当n为偶数时,拆分后的子段需要满足cn,dn为实数an,bn的第一个元素为虚数;当n为奇数时,则需要满足cn,dn为虚数,an,bn的第一个元素为实数,在进行序列再构造的时候,也做出了调整,具体包括如下步骤:
(1)数据符号的分段方法与变形1中步骤(1)的操作一致,得到实虚交错的xn。
(2)将实虚交错的符号序列xn拆分为an,bn,cn,dn四个子段,将每个子段内的元素表示为a[k],b[k],c[k],d[k],其中an和bn的长度为
n为偶数
n为奇数
(3)分为n为偶数、奇数两种情况对数据符号段进行再构造,得到长度为ms的en,将其中的第k个元素记为en[k]:
n为偶数时构造方式为:
当n为奇数时构造方式为:
定义:
其中,
(4)与变形1中的操作一致。
(5)与变形1中的操作一致。
(6)与变形1中的操作一致。
(7)与变形1中的操作一致。
以下为本发明提供的一个更为具体的实施例:
实施例1
本实施例的特定条件为,子载波总数m为1024个,t时刻内发送的数据符号数ms为256个,待发送的符号序列ak长度为10368,其中{a1,a2,...,a128},{a10239,a10240,...,a10368}为0起信号过渡作用,{a129,a130,...,a10238}为复数符号,所使用滤波器的重叠因子k=4。在该具体条件下,实施步骤可以表示为:
(1)将待发送的数据符号序列分为80个长度为256的数据符号段{x0,x1,...,x80},其中第n个数据符号段的第k个元素记为xn[k]:
xn[k]=a128(n-1)+k,0≤k≤255,(n≠0)
(2)用w表示长度为256的窗矢量,对第n个数据符号段xn加窗后得到:
rn=diag(w)·xn
其中,diag函数的作用是将矢量转换成对角矩阵形式,将可供选择的两种窗矢量w记为w1,w2,它们的第k个元素w1[k]和w2[k]可以分别表示为:
这两种窗函数的函数曲线可参见图1所示。
(3)通过256点离散傅里叶变换将rn转换到频域形式rn:
rn=fft(rn)
fft函数表示傅里叶变换,将第m个元素的取值记为rn[m]。
(4)对于第n个数据符号段,将经过进一步调制后得到的en中第m个元素记为en[m],当n为偶数可以表示为:
当n为奇数时,en[m]可以表示为:
其中,
(6)将256个元素en[m]映射到1024个子载波上,第m个子载波上的信息in[m]表示为:
(7)in[m]通过ifft转换到时域后与长度为4096的滤波器相乘,第i个时刻的发送信号sn[i]可以表示为:
其中,g[i]表示滤波器的冲激响应,长度为4096,
(8)将所有经过上述步骤得到的sn[i]进行依次位移512累加,如图2中所示,得到最终的发送信号s[i]。
在如下的仿真实例中,g[i]取phydas滤波器(参考文献:filterbankdesignformulticarriertransmissionandspectrumsensing)或sine滤波器(g[i]=sin(π(i-1)/1024),i=1,2,...,1024)的冲激响应函数,记为g1[i]和g2[i],窗取w1或w2,仿真四种情况下发送信号的papr性能,并与采用传统dft-s-fbmc调制的情况进行对比,ccdf曲线对比见图3。
图3中的横坐标表示的是发送信号峰均功率比的db值,纵坐标表示仿真出来的数据中比该值大的概率。黑色实线表示传统的dft-s-fbmc,方形标注曲线表示使用窗w1和phydas滤波器,圆形曲线表示使用窗w1和sine滤波器,三角形标注曲线表示使用窗w2和phydas滤波器,星形标注曲线表示使用窗w2和sine滤波器。可以看到,无论采用哪种组合,papr的性能都相比于传统的dft-s-fbmc有了明显改善,其中采用窗w2,sin滤波器的情况下性能最优,相比于传统方案提升了2.4db左右。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。