[发明专利]基于最大似然二维搜索的盲频偏估计方法

摘要

本发明公开了一种基于最大似然二维搜索的盲频偏估计方法，根据通信接收机接收到的信号信息，得到接收信号矩阵。初始化频偏搜索列表和码参数搜索列表。依据频偏搜索列表补偿接收信号矩阵得到频偏补偿信号矩阵。依据码参数搜索列表，对每个频偏补偿信号矩阵分别计算频偏‑码参数最大似然函数值。判决频偏‑码参数最大似然函数值中最大值所对应的频偏。由于本发明对接收信号采取了分组合并处理，使得频偏‑码参数最大似然函数的结果不依赖于数据起始位置，对于非完整性数据有很好的估计性能，因此本方法对同步误差有一定的容忍度。

主权项

一种基于最大似然二维搜索的盲频偏估计方法，其特征在于，包括步骤如下：(1)接收信号；通信接收端通过接收天线接收发射端发送的MIMO信号，得到接收信号矩阵Y；(2)初始化；2a)设置频偏搜索列表：依据频偏估计的精度要求设置频偏搜索步长，依据系统的频偏范围设置频偏搜索范围，由频偏搜索步长和频偏搜索范围得到频偏搜索列表，并将列表中第一个频偏作为当前频偏；2b)设置码参数搜索列表：依据发射天线数，将所有可能空时分组编码(STBC)所对应的码参数，即分组长度和每个分组中包含的信息符号数，进行列表，并将列表中第一组码参数作为当前码参数；(3)补偿接收信号矩阵：利用当前频偏对接收信号矩阵Y进行频偏补偿，得到频偏补偿信号矩阵(4)对频偏补偿信号矩阵进行预白化处理，得到预白化信号矩阵；(5)计算频偏‐码参数最大似然函数：5a)将预白化信号矩阵的实部和虚部并联，并进行分组合并处理，即根据当前码参数，将属于同分组的接收信号列向量并联，并计算其相关矩阵；5b)对步骤5a)中的相关矩阵做特征值分解，得到特征值矩阵；5c)根据当前码参数计算频偏‐码参数最大似然函数：其中Nr是接收天线数目，L为接收机接收信号的时隙数，ρk是步骤5a)中得到的特征值矩阵,l是当前码参数所对应的分组长度，n是当前码参数所对应的每个分组中包含的信息符号数；(6)记录当前频偏和当前码参数所对应的频偏‐码参数最大似然函数值；(7)判决是否已遍历码参数搜索列表，若满足，转入步骤(8)，否则，依据码参数搜索列表更新当前码参数，再转入步骤(5)；(8)判决是否已遍历频偏搜索列表，若满足，转入步骤(9)，否则，依据频偏搜索列表更新当前频偏，再转入步骤(4)；(9)取记录集中频偏‐码参数最大似然函数最大值所对应的频偏为判决频偏。