本发明涉及一种适用于正交频分复用(ofdm)系统的基于导频叠加的信道估计方法和检测技术，可不使用专门的导频子载波，提高数据传输速率，属于移动通信系统中的信号处理领域。
背景技术：
ofdm具有传输效率高、实现简单、易与多输入多输出(mimo)结合的优点，是目前第四代(4g)和第五代(5g)蜂窝移动通信系统的主流空中接口技术。为了在接收端进行相干解调，需要估计信道。基于导频(pilot)辅助的信道估计方法需要将接收端已知的导频符号在频域上周期性地插入到数据当中，与数据一起传输。在这种方法中，由于导频需要占用专门的子载波，因而降低了数据的传输速率。为了提高系统的频谱效率，可以将导频叠加到数据上形成所谓的导频叠加(superimposedtraining)方式。为了将叠加导频产生的干扰分散到各个数据符号上，可以先对数据进行预编码(precoding)，然后将导频插入位置处的编码后信号删除，形成dnst(datanullingsuperimposedtraining)方案。通过引入干扰控制因子，pdst(partial-datasuperimposedtraining)方案不删除编码后的信号，得到了更好的性能。但是，虽然增加了控制因子，但导频与数据之间的相互干扰仍然存在，对系统性能仍有较大影响。
技术实现要素：
本发明的目的是为了解决上述问题，提供一种导频叠加信道估计方法，针对ofdm系统，将需发送的信号进行预编码，然后通过移位选择合适的导频叠加位置，最小化导频与数据之间的相互干扰，从而有效降低接收端的误码率。
为了实现上述目的，本发明采用如下的技术方案：一种导频叠加信道估计方法，包括以下步骤：
步骤1)、发送端将需要发送的数据符号组成的向量d乘以一个n×n酉矩阵w进行预编码，得到n×1维列向量s；
步骤2)、计算s中位置为{i,i df,i 2df,k,i (np-1)df}处元素的功率和λi。由于i取值为{0,1,2,k,np-1}，所以λi共有np个值；
步骤3)、发送端选择最小的
步骤4)、发送端将s在{i*,i* df,i* 2df,k,i* (np-1)df}处的元素乘以α，并将np个导频叠加到该位置处，得到ofdm的频域发送信号x；
步骤5)、发送端将x经过反傅里叶变换(idft)变换到时域，然后加循环前缀，调制到载波发送出去。接收端通过相反过程得到经过多径信道的接收信号y；
步骤6)、接收端利用y在{i*,i* df,i* 2df,k,i* (np-1)df}处的信号估计信道，然后抵消导频后对数据进行检测。
其中，n表示ofdm的子载波个数；np表示导频的个数；df＝n/np表示导频间隔；0＜α＜1表示干扰控制因子。
有益效果：
1、本发明的方法能提高信道估计精度，从而降低误码率；
2、本发明的方法将导频叠加在数据上，不需要额外的导频子载波，可以提高系统的频谱效率。
3、本发明的方法适用范围广，除了ofdm，还适用于下一代通信系统的多载波技术，包括滤波器组多载波(fbmc)、广义频分复用(gdmc)和通用滤波多载波(ufmc)，等等。
附图说明
图1是本发明的计算机仿真结果图。
具体实施方式
下面结合附图对本发明的具体实施方式进行描述，以便本领域的技术人员更好的理解本发明。需要特别注意的是，在以下的描述中，由于对一些已知的技术和功能的详细描述可能会淡化本发明的内容，这些描述在这里将被忽略。
考虑一个ofdm系统，有n个子载波，每个子载波带宽为bsub，总带宽为b＝nbsub。需要传输的数据符号用向量可表示为d＝[d(0),d(1),...,d(n-1)]t，其中d(k)表示调制在第k个子载波上的数据符号，[·]t表示向量转置。d(k)可以采用正交相移键控(qpsk)或其它的调制方式。为了将叠加导频产生的干扰分散到各个数据符号上，d乘以一个n×n酉矩阵w进行预编码，得到s＝wd。w可取沃尔什哈达码(walshhadamard)矩阵，或其它酉矩阵，但需满足i＝whw。这里，i表示n×n的单位阵，[·]h表示共轭转置。假设np个相同功率和间隔的导频叠加到数据上用于接收端的信道估计。根据采样定理，导频的间隔df＝n/np需小于信道的最大时延扩展τmax，即：
其中，t＝1/b表示采样周期。为了计算位置为{i,i df,i 2df,k,i (np-1)df}处信号的功率和λi，我们定义一个n×n的对角矩阵j(i)，其第k个对角元素可表示为：
其中，i＝0,1,...,df-1；n＝0,1,...,np-1。于是，我们得到导频初始位置位移为i时，λi的计算公式为：
λi＝[j(i)s]h[j(i)s][公式三]
本发明选择最小的
为了控制导频与数据之间的相互干扰，我们对导频叠加处的s乘以实数因子α，其取值范围为0＜α＜1。最终得到的频域发送信号为：
x＝(i-αj(i*))wd j(i*)p(i*)[公式五]
p是n×1的向量，其第k个元素表示为：
其中，p(n)表示第n个导频符号。最后，发送端将x经过反傅里叶变换(idft)变换到时域，然后加长度为ncp的循环前缀，调制到载波发送出去。接收端通过下变频等相反过程得到经过多径信道的接收信号。
假设信道在连续的m个ofdm符号内保持不变，则第m个接收到的频域信号为：
ym＝hxm vm[公式七]
其中，vm表示方差为
其中，h(k)表示第k个子载波上的信道频域响应(cfr)。假设
相应的时域信道冲激响应(cir)为：
其中，f1表示进过裁剪的dft矩阵，即f1的第p行第q列的元素满足
其中，q＝0,1,...,ncp-1；
变换到频域可得到每个子载波上的cfr，即：
其中，f2表示进过裁剪的dft矩阵，即f2的第p行第q列的元素满足
其中，p＝0,1,...,n-1；q＝0,1,...,ncp-1。
最后，第m个ofdm符号上传输的数据由下式解调：
其中，n×n的对角矩阵gm表示信道的均衡矩阵，即其第k个对角元素为：
为了提高数据检测的精度，可以进行迭代检测。此时，第l次迭代的结果可由第l－1次得到：
其中，
图1给出了本发明的仿真结果，仿真参数如下：m取10，n取64，qpsk调制，np＝8，α＝0.71。w选择64×64的walshhadamard阵。多径瑞利衰落信道(rayleigh)有8根径(path)，时延分别为{0,1,k,7}，第l根径的平均功率为
根据以上描述，可以得到一种适用于ofdm的导频叠加信道估计方法为：
步骤1)、发送端将需要发送的数据符号组成的向量d乘以一个n×n酉矩阵w进行预编码，得到n×1维列向量s；；
步骤2)、计算s中位置为{i,i df,i 2df,k,i (np-1)df}处元素的功率和λi。由于i取值为{0,1,2,k,np-1}，所以λi共有np个值；
步骤3)、发送端选择最小的
步骤4)、发送端将s在{i*,i* df,i* 2df,k,i* (np-1)df}处的元素乘以α，并将np个导频叠加到该位置处，得到ofdm的频域发送信号x；
步骤5)、发送端将x经过反傅里叶变换(idft)变换到时域，然后加循环前缀，调制到载波发送出去。接收端通过相反过程得到经过多径信道的接收信号y；
步骤6)、接收端利用y在{i*,i* df,i* 2df,k,i* (np-1)df}处的信号估计信道，然后抵消导频后对数据进行检测。
以上所述即使本发明的实施方法，应当指出：对于本技术领域的普通技术人员来说，再不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
1.一种导频叠加信道估计方法，其特征在于，该方法包括以下步骤：
1)发送端将需要发送的数据符号组成的向量d乘以一个n×n酉矩阵w进行预编码，得到n×1维列向量s；
2)计算s中位置为{i,i df,i 2df,k,i (np-1)df}处元素的功率和λi；由于i取值为{0,1,2,k,np-1}，所以λi共有np个值；
3)发送端选择最小的
4)发送端将s在{i*,i* df,i* 2df,k,i* (np-1)df}处的元素乘以α，并将np个导频叠加到该位置处，得到ofdm的频域发送信号x；
5)发送端将x经过反傅里叶变换变换到时域，然后加循环前缀，调制到载波发送出去，接收端通过相反过程得到经过多径信道的接收信号y；
6)接收端利用y在{i*,i* df,i* 2df,k,i* (np-1)df}处的信号估计信道，然后抵消导频后对数据进行检测；
其中，n表示ofdm的子载波个数；np表示导频的个数；df＝n/np表示导频间隔；0＜α＜1表示干扰控制因子。
技术总结