基于最小均方误差均衡的混合载波阶数选择方法与流程

本发明属于无线通信领域。

背景技术：

随着科技的进步，移动通信的信道环境变得越来越复杂，在时域与频域同时衰落的时频双选信道下，传统的单载波体制系统和多载波体制系统都难以获得很好的性能，而基于加权分数傅里叶变换的混合载波体制系统却能够在这种恶劣的信道下获得相对较好的效果。对于混合载波体制系统，其误码性能的好坏主要取决于变换阶数的选择，因而亟需一种变换阶数的选择方法来选择最优的变换阶数，进一步提升系统性能。

技术实现要素：

本发明是为了解决混合载波系统由于变换阶数难以选择，而导致系统误码性能高的问题，本发明提供了一种基于最小均方误差均衡的混合载波阶数选择方法。

基于最小均方误差均衡的混合载波阶数选择方法，该方法包括如下步骤：

步骤一、采用接收天线接收当前数据块信号y(i)，并对当前数据块信号y(i)处理，获得当前数据块信号的平均功率

其中，n为数据块信号的长度，且n为正整数，y(i)^*表示y(i)的共轭，

步骤二、根据噪声功率谱密度n0和系统带宽b，计算噪声功率pz＝n0b，其中，n0＝-174dbm/hz；

步骤三、根据步骤一得到的当前数据块信号的平均功率py和步骤二中得到的噪声功率pz，计算有用信号平均功率ps，其中，ps＝py-pz；

步骤四、根据当前数据块信号y(i)中的导频序列，估计出信道状态信息htl；

步骤五、采用mmse均衡方法对信道状态信息htl进行处理，获得均衡矩阵g；

步骤六、判断当前混合载波阶数α的值是否大于4，若结果为是，执行步骤十一；若结果为否，执行步骤七，混合载波阶数α的初始值为0；

步骤七、根据加权分数傅里叶变换矩阵f^α、信道状态信息htl和均衡矩阵g计算信号残余干扰分布矩阵

步骤八、根据加权分数傅里叶变换矩阵f^α和均衡矩阵g计算噪声干扰分布矩阵

步骤九、根据有用信号平均功率ps、噪声功率pz、信号残余干扰分布矩阵和噪声干扰分布矩阵计算当前数据块信号y(i)中每比特所受到的干扰功率u^α(n)；n＝0,1,…n-1；

步骤十、通过干扰功率u^α(n)计算当前数据块信号y(i)中每比特所受干扰功率的方差v(α)＝var(u^α(n))，其中，α＝α+0.1，执行步骤六；

步骤十一：在获得的所有方差v(α)中挑选出最小值，该方差v(α)的最小值所对应的混合载波阶数α为最优阶数，且该最优阶数α作为下一个数据块信号的调制阶数，从而完成最优阶数α的选择。

步骤四中所述的信道状态信息htl的增益矩阵表达式为：

其中，hn,l表示时变信道第l个抽头在第n个采样时刻的信道增益系数，n∈[0,n-1]，l∈[0,l-1]。

步骤五中，所述的采用mmse均衡方法对信道状态信息htl进行处理，获得均衡矩阵g的具体过程，采用如下公式二实现：

其中，i表示单位矩阵，表示htl的共轭转置。

步骤七中，所述的根据加权分数傅里叶变换矩阵f^α、信道状态信息htl和均衡矩阵g计算信号残余干扰分布矩阵的具体过程，采用如下公式三实现:

其中，i表示单位矩阵。

步骤八中，所述的根据加权分数傅里叶变换矩阵f^α和均衡矩阵g计算噪声干扰分布矩阵的具体过程，采用如下公式四实现:

步骤九中，所述的根据有用信号平均功率ps、噪声功率pz、信号残余干扰分布矩阵和噪声干扰分布矩阵计算当前数据块信号y(i)中每比特所受到的干扰功率u^α(n)的具体过程，采用如下公式五实现：

原理说明：采用信道估计方法获得信道的信道状态信息；根据接收到的信号计算有用信号和噪声的平均功率；根据获得的信道状态信息、有用信号功率和噪声功率计算每比特信号所受干扰的功率；改变变换阶数，同时利用最小均方误差均衡公式和信道状态信息计算每一变换阶数下干扰功率分布的方差；比较所有计算出的方差，其中最小值所对应的变换阶数，即为所要选择的最优变换阶数。

本发明带来的有益效果是，本发明通过联合考虑信道状态信息和变换阶数对干扰功率分布的影响，为基于最小均方误差均衡的混合载波系统提供了一种变换阶数选择方法，通过该方法，可以为混合载波系统选择出最优的变换阶数，进一步优化系统性能。

具体实施方式

具体实施方式一：本实施方式所述的基于最小均方误差均衡的混合载波阶数选择方法，该方法包括如下步骤：

步骤一、采用接收天线接收当前数据块信号y(i)，并对当前数据块信号y(i)处理，获得当前数据块信号的平均功率

其中，n为数据块信号的长度，且n为正整数，y(i)^*表示y(i)的共轭，i为变量；

步骤二、根据噪声功率谱密度n0和系统带宽b，计算噪声功率pz＝n0b，其中，n0＝-174dbm/hz；

步骤三、根据步骤一得到的当前数据块信号的平均功率py和步骤二中得到的噪声功率pz，计算有用信号平均功率ps，其中，ps＝py-pz；

步骤四、根据当前数据块信号y(i)中的导频序列，估计出信道状态信息htl；

步骤五、采用mmse均衡方法对信道状态信息htl进行处理，获得均衡矩阵g；

步骤六、判断当前混合载波阶数α的值是否大于4，若结果为是，执行步骤十一；若结果为否，执行步骤七，混合载波阶数α的初始值为0；

步骤七、根据加权分数傅里叶变换矩阵f^α、信道状态信息htl和均衡矩阵g计算信号残余干扰分布矩阵

步骤八、根据加权分数傅里叶变换矩阵f^α和均衡矩阵g计算噪声干扰分布矩阵

步骤九、根据有用信号平均功率ps、噪声功率pz、信号残余干扰分布矩阵和噪声干扰分布矩阵计算当前数据块信号y(i)中每比特所受到的干扰功率u^α(n)；n＝0,1,…n-1；

步骤十、通过干扰功率u^α(n)计算当前数据块信号y(i)中每比特所受干扰功率的方差v(α)＝var(u^α(n))，其中，α＝α+0.1，执行步骤六；

步骤十一：在获得的所有方差v(α)中挑选出最小值，该方差v(α)的最小值所对应的混合载波阶数α为最优阶数，且该最优阶数α作为下一个数据块信号的调制阶数，从而完成最优阶数α的选择。

本实施方式中，本发明主要在接收端对所接收的当前数据块信号进行最优阶数提取，且提取完成后可通过回传链路传给发送端，该最优阶数作为发送端发送给接收端的下一数据块的调制阶数，也即最优阶数作为接收端所接收的下一数据块的调制阶数。

在接收端接收信号之前，发送端对所要发送的信号进行处理：

步骤一一，将待发送信号进行调制，并分成长度为n的数据块；n为正整数；

步骤一二，长度为n的数据块进行-α阶加权分数傅里叶变换获得调制后的信号；

其中，α阶的加权分数傅里叶变换的矩阵表示形式如下：

f^α＝w0(α)i+w1(α)f+w2(α)pi+w3(α)pf(公式六)，

式中：{w0(α)，w1(α)，w2(α)，w3(α)}为加权系数，具体表达式为：

i为单位阵，f为离散傅里叶变换矩阵，f归一化表达式为：

其中：q＝e^-j2π/n，p为置换矩阵，p的具体形式为：

步骤一三，对调制后的信号添加长度为l的循环前缀，再插入导频序列，并送入信道进行传输，发至接收端。

具体实施方式二：本实施方式与具体实施方式一所述的基于最小均方误差均衡的混合载波阶数选择方法的区别在于，步骤四中所述的信道状态信息htl的增益矩阵表达式为：

其中，hn,l表示时变信道第l个抽头在第n个采样时刻的信道增益系数，n∈[0,n-1]，l∈[0,l-1]。

具体实施方式三：本实施方式与具体实施方式一所述的基于最小均方误差均衡的混合载波阶数选择方法的区别在于，步骤五中，所述的采用mmse均衡方法对信道状态信息htl进行处理，获得均衡矩阵g的具体过程，采用如下公式二实现：

其中，i表示单位矩阵，表示htl的共轭转置。

本实施方式中，所述的mmse(minimummeansquareerror，最小均方误差)均衡方法，通过这种均衡方法，可以有效抑制均衡过程中的噪声放大问题，使均衡后的信号与原始信号之间的方差达到最小。

具体实施方式四：本实施方式与具体实施方式一所述的基于最小均方误差均衡的混合载波阶数选择方法的区别在于，步骤七中，所述的根据加权分数傅里叶变换矩阵f^α、信道状态信息htl和均衡矩阵g计算信号残余干扰分布矩阵的具体过程，采用如下公式三实现:

其中，i表示单位矩阵。

本实施方式中，经过均衡和变换后所恢复出来的信号与原始信号之间必然会存在一些失真，而这些失真可以看作是由信号间的残余干扰造成的，因而用代表恢复信号的矩阵f^αghtlf^-α减去代表原始信号的矩阵i即可得到信号残余干扰分布矩阵。

具体实施方式五：本实施方式与具体实施方式一所述的基于最小均方误差均衡的混合载波阶数选择方法的区别在于，步骤八中，所述的根据加权分数傅里叶变换矩阵f^α和均衡矩阵g计算噪声干扰分布矩阵的具体过程，采用如下公式四实现:

本实施方式中，信号在通过信道时，会受到加性高斯白噪声的影响。这种噪声的功率本来是均匀分布的，但在经过均衡和变换后，其功率不再呈现均匀分布，而是变为如公式四所示的分布形式。

具体实施方式六：本实施方式与具体实施方式一所述的基于最小均方误差均衡的混合载波阶数选择方法的区别在于，步骤九中，所述的根据有用信号平均功率ps、噪声功率pz、信号残余干扰分布矩阵和噪声干扰分布矩阵计算当前数据块信号y(i)中每比特所受到的干扰功率u^α(n)的具体过程，采用如下公式五实现：

本实施方式中，信号所受到的干扰由两部分组成，一部分是信号的残余干扰，另一部分是噪声干扰。对于信号残余干扰分布矩阵和噪声干扰分布矩阵其第n行的2范数分别表示了信号残余干扰和噪声干扰在信号第n比特上的分布情况。因而，只要用信号功率乘以信号残余干扰功率矩阵第n行的2范数再加上噪声功率与噪声干扰分布矩阵第n行2范数的乘积，即可获得信号第n比特上所受的干扰总功率。