一种时间偏差估计方法及装置与流程

1.本技术涉及通信技术领域，尤其涉及一种时间偏差估计方法及装置。


背景技术：

2.在5g多用户-多输入多输出(multi-user multiple-input multiple-output，mu-mimo)通信场景下，由于基站与终端时钟振荡器随时间和温度变化产生频率漂移，基站与终端的定时同步存在随时间变化的偏差，也就是时间偏差，从而导致正交频分复用(orthogonal frequency division multiplexing，ofdm)系统频域信道估计不准。
3.现有技术中进行时偏校准的方法有两种，一是基于逆离散傅里叶变换(invert discrete fourier tranform，idft)将信道频域估计转换成时域信道估计，然后搜索峰值所在位置，计算其对应的时偏值。一种是基于现代谱估计的方法，如基于旋转不变技术的信号参数估计方法(estimation of signal parameters via rotational invariance techniques，esprit)及多信号分类法(multiple signal classification，music)等，主要通过子空间分解的方法区分信号子空间和噪声子空间，然后从频域信号空间中提取相位信息，从而获得时偏值。通过时偏值对原始信道估计结果进行校准。
4.然而，基于idft的时偏校准在多径场景中，信道各子径之间相互干扰，导致idft的峰值位置不能反映真实的时偏值，最终信道预测的结果精度低。基于现代谱估计的方法对多径的分辨能力较高，但运算的复杂度过高，实用性低。


技术实现要素：

5.本技术实施例提供一种时间偏差估计方法、基站、装置及存储介质，用以解决现有技术中运算复杂度和精度不能同时保证的缺陷，确保精度的同时降低运算的复杂度。
6.第一方面，本技术实施例提供一种时间偏差估计方法，包括：
7.根据滑动窗口确定基于第一信道估计周期的信道估计的原始信道频域响应矩阵对应的第一数据矩阵，第一数据矩阵为l
×
j的矩阵，l为滑动窗口的长度，j为滑动窗口在原始信道频域响应矩阵中滑动的次数，l和j均为正整数；
8.对第一数据矩阵进行滑动窗口离散傅里叶变换，确定第一子空间矩阵；
9.根据第一子空间矩阵进行时间偏差估计。
10.可选地，根据本技术一个实施例的时间偏差估计方法，对第一数据矩阵进行滑动窗口离散傅里叶变换，确定第一子空间矩阵，包括：
11.对第一数据矩阵中的每一行对应的行向量进行滑动窗口离散傅里叶变换，确定第二数据矩阵；
12.确定第二数据矩阵中每一列对应的列向量的模值；
13.根据模值超过第一预设门限值的列向量确定第一子空间矩阵。
14.可选地，根据本技术一个实施例的时间偏差估计方法，对第一数据矩阵中的每一行对应的行向量进行滑动窗口离散傅里叶变换，确定第二数据矩阵，包括：
15.根据公式一对第一数据矩阵的第一行对应的行向量进行点数为n的离散傅里叶变换，公式一的表达式为：
[0016][0017]
其中，xk为离散傅里叶变换结果，l为第一数据矩阵的行序号，k为频点序号，x(i)为采样信号，n为离散傅里叶变换的点数，n和j的取值相同；l、k和n均为正整数；
[0018]
在l大于1的情况下，根据公式二确定第一数据矩阵的第l行对应的行向量的离散傅里叶变换结果，公式二的表达式为：
[0019]
xk(l)＝e
j2πk/n
[xk(l-1)+x(l+n-1)-x(l-1)]
[0020]
根据公式一和公式二确定第二数据矩阵，第二数据矩阵的表达式为：
[0021][0022]
其中，v为第二数据矩阵。
[0023]
可选地，根据本技术一个实施例的时间偏差估计方法，根据第一子空间矩阵进行时间偏差估计，包括：
[0024]
根据第一子空间矩阵的除最后一行之外的前m-1个行向量和除第一行之外的后m-1个行向量确定第一矩阵；其中，m为第一子空间矩阵的行数；
[0025]
根据第一子空间矩阵的第m行对应的行向量和第一矩阵确定频谱矩阵；
[0026]
根据频谱矩阵确定基于第一信道估计周期的信道估计的时间偏差估计结果。
[0027]
可选地，根据本技术一个实施例的时间偏差估计方法，根据第一子空间矩阵的第m行对应的行向量和第一矩阵确定频谱矩阵的计算公式为：
[0028][0029]
其中，φ为频谱矩阵，ψ为第一矩阵，v为第一子空间矩阵的第m行对应的行向量进行共轭转置后的向量。
[0030]
可选地，根据本技术一个实施例的时间偏差估计方法，根据第一子空间矩阵进行时间偏差估计，包括：
[0031]
根据第二信道估计周期的时间偏差估计结果确定第一信道估计周期的谱峰搜索区间；其中，第二信道估计周期为第一信道估计周期的前一个信道估计周期；
[0032]
在第一信道估计周期的谱峰搜索区间的长度大于第二预设门限值的情况下，对空间谱函数进行黄金分割搜索确定第一信道估计周期的时间偏差估计结果。
[0033]
可选地，根据本技术一个实施例的时间偏差估计方法，根据频谱矩阵确定基于第一信道估计周期的信道估计的时间偏差估计结果之后，还包括：
[0034]
根据第一信道估计周期的时间偏差估计结果确定第三信道估计周期的谱峰搜索区间；其中，第三信道估计周期为第一信道估计周期的后一个信道估计周期，第一信道估计周期为信道估计开始后的第一个信道估计周期；
[0035]
在第三信道估计周期的谱峰搜索区间的长度大于第二预设门限值的情况下，对空间谱函数进行黄金分割搜索确定第三信道估计周期的时间偏差估计结果。
[0036]
可选地，根据本技术一个实施例的时间偏差估计方法，对空间谱函数进行黄金分割搜索确定第一信道估计周期的时间偏差估计结果，包括：
[0037]
以第二预设门限值为收敛精度，根据黄金分割搜索确定空间谱函数在第一信道估计周期的谱峰搜索区间内的极小值；
[0038]
根据极小值确定目标搜索区间，第一信道估计周期的时间偏差估计结果为目标搜索区间的中值。
[0039]
第二方面，本技术实施例还提供一种网络设备，包括：存储器，收发机，处理器：
[0040]
存储器，用于存储计算机程序；收发机，用于在处理器的控制下收发数据；处理器，用于读取存储器中的计算机程序并执行以下操作：
[0041]
根据滑动窗口确定基于第一信道估计周期的信道估计的原始信道频域响应矩阵对应的第一数据矩阵，第一数据矩阵为l
×
j的矩阵，l为滑动窗口的长度，j为滑动窗口在原始信道频域响应矩阵中滑动的次数，l和j均为正整数；
[0042]
对第一数据矩阵进行滑动窗口离散傅里叶变换，确定第一子空间矩阵；
[0043]
根据第一子空间矩阵进行时间偏差估计。
[0044]
可选地，根据本技术一个实施例的网络设备，对第一数据矩阵进行滑动窗口离散傅里叶变换，确定第一子空间矩阵，包括：
[0045]
对第一数据矩阵中的每一行对应的行向量进行滑动窗口离散傅里叶变换，确定第二数据矩阵；
[0046]
确定第二数据矩阵中每一列对应的列向量的模值；
[0047]
根据模值超过第一预设门限值的列向量确定第一子空间矩阵。
[0048]
可选地，根据本技术一个实施例的网络设备，对第一数据矩阵中的每一行对应的行向量进行滑动窗口离散傅里叶变换，确定第二数据矩阵，包括：
[0049]
根据公式一对第一数据矩阵的第一行对应的行向量进行点数为n的离散傅里叶变换，公式一的表达式为：
[0050][0051]
其中，xk为离散傅里叶变换结果，l为第一数据矩阵的行序号，k为频点序号，x(i)为采样信号，n为离散傅里叶变换的点数，n和j的取值相同；l、k和n均为正整数；
[0052]
在l大于1的情况下，根据公式二确定第一数据矩阵的第l行对应的行向量的离散傅里叶变换结果，公式二的表达式为：
[0053]
xk(l)＝e
j2πk/n
[xk(l-1)+x(l+n-1)-x(l-1)]
[0054]
根据公式一和公式二确定第二数据矩阵，第二数据矩阵的表达式为：
[0055]
[0056]
其中，v为第二数据矩阵。
[0057]
可选地，根据本技术一个实施例的网络设备，根据第一子空间矩阵进行时间偏差估计，包括：
[0058]
根据第一子空间矩阵的除最后一行之外的前m-1个行向量和除第一行之外的后m-1个行向量确定第一矩阵；其中，m为第一子空间矩阵的行数；
[0059]
根据第一子空间矩阵的第m行对应的行向量和第一矩阵确定频谱矩阵；
[0060]
根据频谱矩阵确定基于第一信道估计周期的信道估计的时间偏差估计结果。
[0061]
可选地，根据本技术一个实施例的网络设备，根据第一子空间矩阵的第m行对应的行向量和第一矩阵确定频谱矩阵的计算公式为：
[0062][0063]
其中，φ为频谱矩阵，ψ为第一矩阵，v为第一子空间矩阵的第m行对应的行向量进行共轭转置后的向量。
[0064]
可选地，根据本技术一个实施例的网络设备，根据第一子空间矩阵进行时间偏差估计，包括：
[0065]
根据第二信道估计周期的时间偏差估计结果确定第一信道估计周期的谱峰搜索区间；其中，第二信道估计周期为第一信道估计周期的前一个信道估计周期；
[0066]
在第一信道估计周期的谱峰搜索区间的长度大于第二预设门限值的情况下，对空间谱函数进行黄金分割搜索确定第一信道估计周期的时间偏差估计结果。
[0067]
可选地，根据本技术一个实施例的网络设备，根据频谱矩阵确定基于第一信道估计周期的信道估计的时间偏差估计结果之后，还包括：
[0068]
根据第一信道估计周期的时间偏差估计结果确定第三信道估计周期的谱峰搜索区间；其中，第三信道估计周期为第一信道估计周期的后一个信道估计周期，第一信道估计周期为信道估计开始后的第一个信道估计周期；
[0069]
在第三信道估计周期的谱峰搜索区间的长度大于第二预设门限值的情况下，对空间谱函数进行黄金分割搜索确定第三信道估计周期的时间偏差估计结果。
[0070]
可选地，根据本技术一个实施例的网络设备，对空间谱函数进行黄金分割搜索确定第一信道估计周期的时间偏差估计结果，包括：
[0071]
以第二预设门限值为收敛精度，根据黄金分割搜索确定空间谱函数在第一信道估计周期的谱峰搜索区间内的极小值；
[0072]
根据极小值确定目标搜索区间，第一信道估计周期的时间偏差估计结果为目标搜索区间的中值。
[0073]
第三方面，本技术实施例还提供一种时间偏差估计装置，包括：
[0074]
第一确定单元，用于根据滑动窗口确定基于第一信道估计周期的信道估计的原始信道频域响应矩阵对应的第一数据矩阵，第一数据矩阵为l
×
j的矩阵，l为滑动窗口的长度，j为滑动窗口在原始信道频域响应矩阵中滑动的次数，l和j均为正整数；
[0075]
第二确定单元，用于对第一数据矩阵进行滑动窗口离散傅里叶变换，确定第一子空间矩阵；
[0076]
估计单元，用于根据第一子空间矩阵进行时间偏差估计。
[0077]
可选地，根据本技术一个实施例的时间偏差估计装置，第二确定单元包括第一确定模块、第二确定模块和第三确定模块；
[0078]
第一确定模块用于对第一数据矩阵中的每一行对应的行向量进行滑动窗口离散傅里叶变换，确定第二数据矩阵；
[0079]
第二确定模块用于确定第二数据矩阵中每一列对应的列向量的模值；
[0080]
第三确定模块用于根据模值超过第一预设门限值的列向量确定第一子空间矩阵。
[0081]
可选地，根据本技术一个实施例的时间偏差估计装置，第一确定模块包括第一确定子模块、第二确定子模块和第三确定子模块；
[0082]
第一确定子模块用于根据下列公式一对第一数据矩阵的第一行对应的行向量进行点数为n的离散傅里叶变换，公式一的表达式为：
[0083][0084]
其中，xk为离散傅里叶变换结果，l为第一数据矩阵的行序号，k为频点序号，x(i)为采样信号，n为离散傅里叶变换的点数，n和j的取值相同；l、k和n均为正整数；
[0085]
第二确定子模块用于在l大于1的情况下，根据公式二确定第一数据矩阵的第l行对应的行向量的离散傅里叶变换结果，公式二的表达式为：
[0086]
xk(l)＝e
j2πk/n
[xk(l-1)+x(l+n-1)-x(l-1)]
[0087]
第三确定子模块用于根据公式一和公式二确定第二数据矩阵，第二数据矩阵的表达式为：
[0088][0089]
其中，v为第二数据矩阵。
[0090]
可选地，根据本技术一个实施例的时间偏差估计装置，估计单元包括第四确定模块、第五确定模块和第六确定模块；
[0091]
第四确定模块用于根据第一子空间矩阵的除最后一行之外的前m-1个行向量和除第一行之外的后m-1个行向量确定第一矩阵；其中，m为第一子空间矩阵的行数；
[0092]
第五确定模块用于根据第一子空间矩阵的第m行对应的行向量和第一矩阵确定频谱矩阵；
[0093]
第六确定模块用于根据频谱矩阵确定基于第一信道估计周期的信道估计的时间偏差估计结果。
[0094]
可选地，根据本技术一个实施例的时间偏差估计装置，根据第一子空间矩阵的第m行对应的行向量和第一矩阵确定频谱矩阵的计算公式为：
[0095][0096]
其中，φ为频谱矩阵，ψ为第一矩阵，v为第一子空间矩阵的第m行对应的行向量进
行共轭转置后的向量。
[0097]
可选地，根据本技术一个实施例的时间偏差估计装置，估计模块包括第七确定模块和第八确定模块；
[0098]
第七确定模块用于根据第二信道估计周期的时间偏差估计结果确定第一信道估计周期的谱峰搜索区间；其中，第二信道估计周期为第一信道估计周期的前一个信道估计周期；
[0099]
第八确定模块用于在第一信道估计周期的谱峰搜索区间的长度大于第二预设门限值的情况下，对空间谱函数进行黄金分割搜索确定第一信道估计周期的时间偏差估计结果。
[0100]
可选地，根据本技术一个实施例的时间偏差估计装置，估计单元还包括第九确定模块和第十确定模块；
[0101]
第九确定模块用于根据第一信道估计周期的时间偏差估计结果确定第三信道估计周期的谱峰搜索区间；其中，第三信道估计周期为第一信道估计周期的后一个信道估计周期，第一信道估计周期为信道估计开始后的第一个信道估计周期；
[0102]
第十确定模块用于在第三信道估计周期的谱峰搜索区间的长度大于第二预设门限值的情况下，对空间谱函数进行黄金分割搜索确定第三信道估计周期的时间偏差估计结果。
[0103]
可选地，根据本技术一个实施例的时间偏差估计装置，第八确定模块包括第四确定子模块和第五确定子模块；
[0104]
第四确定子模块用于以第二预设门限值为收敛精度，根据黄金分割搜索确定空间谱函数在第一信道估计周期的谱峰搜索区间内的极小值；
[0105]
第五确定子模块用于根据极小值确定目标搜索区间，第一信道估计周期的时间偏差估计结果为目标搜索区间的中值。
[0106]
第四方面，本技术实施例还提供一种计算机可读存储介质，计算机可读存储介质存储有计算机程序，计算机程序用于使处理器执行第一方面描述的时间偏差估计方法的步骤。
[0107]
第五方面，本技术实施例还提供一种时间偏差估计的装置，时间偏差估计的装置包括处理器和接口电路，接口电路用于接收计算机执行指令并传输至处理器，处理器运行计算机执行指令以实现第一方面描述的时间偏差估计方法中的步骤。
[0108]
可选地，该时间偏差估计的装置还包括存储器，存储器与处理器耦合，处理器用于实现第一方面描述的时间偏差估计方法中的步骤。
[0109]
可选地，该存储器用于存储程序指令和数据。该存储器与该处理器耦合，该处理器可以调用并执行该存储器中存储的程序指令，用于实现第一方面描述的时间偏差估计方法中的步骤。
[0110]
可选地，该时间偏差估计的装置还包括通信接口，该通信接口用于该时间偏差估计的装置与其他设备进行通信。
[0111]
可选地，该时间偏差估计的装置为芯片或芯片系统时，该通信接口包括该芯片或芯片系统上输入/输出接口、接口电路、输出电路、输入电路、管脚或相关电路等。该处理器体现为处理电路或逻辑电路。
[0112]
第六方面，本技术实施例还提供一种包含指令的计算机程序产品，指令被计算机执行时使得计算机执行第一方面描述的时间偏差估计方法的步骤。本技术实施例提供的时间偏差估计方法及装置，通过对信道频域响应矩阵进行滑动窗口傅里叶变换确定子空间矩阵，显著降低了多组离散傅里叶变换的运算复杂度，同时使用子空间近似方法代替现代谱估计方法中的特征值分解方法，保证运算精度的同时进一步降低了运算复杂度。
附图说明
[0113]
为了更清楚地说明本技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0114]
图1是现有技术提供的时偏随时间变化的曲线示意图；
[0115]
图2是本技术实施例提供的时间偏差估计方法的流程示意图之一；
[0116]
图3是本技术实施例提供的时间偏差估计方法的流程示意图之二；
[0117]
图4是本技术实施例提供的基站的结构示意图；
[0118]
图5是本技术实施例提供的时间偏差估计装置的结构示意图。
具体实施方式
[0119]
对于任意一个无线收发通信系统，无论是5g还是4g，接收端和发送端采用的晶振频率之间都会有物理上的偏差。接收端需要将这个偏差估计出来并对数据进行补偿。
[0120]
在5g多用户多输入多输出通信场景下，由于基站和终端之间的时钟晶振随着时间和温度的变化产生频率漂移，产生时间偏差从而导致正交频分复用系统频道信域估计不准确，信道预测能力下降。
[0121]
在用户高速移动时，基站与终端之间的无线信道迅速变化，基站的波束赋形权值和无线信道失配，这就需要利用多组历史信道估计结果来对未来时刻信道进行预测，从而获得更加准确的波束赋形权值，提升传输速率。
[0122]
但是由于时间偏差的存在，历史时刻的信道估计精准度不够，会增加信道预测的误差。
[0123]
图1是现有技术提供的时偏随时间变化的曲线示意图，如图1所示，由于时钟晶振随时间和温度变化产生频率漂移，基站和终端的定时同步存在随时间的偏差，即时偏。其中，横坐标为时间变化(snap index)，每个样点值表示对一次信道估计结果的快照，纵坐标为当前时刻信道估计结果中的时偏值(delay)。
[0124]
从图1中可以看出，时偏值随着时间的推移出现轻微的抖动，抖动范围在纳米级，与终端所选用的时钟周期相当。基站使用上行定时同步算法，在上行定时的跟踪阶段测量上行定时偏差，其本身具有一定的时偏校准能力。
[0125]
现有技术中上行定时同步算法主要是通过确定ofdm符号的起始点和结束点，以供调解使用，对时偏校准的精度要求不高，只需要在整数时钟周期量级。
[0126]
但是信道预测过程中对历史信道时偏校准的精度要求很高，需要在小数时钟周期量级，以符合信道预测的高精度要求。
[0127]
基于idft的时偏校准在多径场景中，信道各子径之间相互干扰，导致idft的峰值位置不能反映真实的时偏值，最终信道预测的结果精度低。基于现代谱估计的方法对多径的分辨能力较高，但运算的复杂度过高，实用性低。
[0128]
为解决上述技术问题，本技术实施例提供一种时间偏差估计方法及装置，在保证信道预测过程中时偏值估计精度的同时，降低运算的复杂度。
[0129]
为使本发明的目的、技术方案和优点更加清楚，下面将结合本发明中的附图，对本发明中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0130]
图2是本技术实施例提供的时间偏差估计方法的流程示意图之一，如图2所示，本技术实施例提供一种时间偏差估计方法，其执行主体可以为基站，该方法包括：
[0131]
步骤201、根据滑动窗口确定基于第一信道估计周期的信道估计的原始信道频域响应矩阵对应的第一数据矩阵，第一数据矩阵为l
×
j的矩阵，l为滑动窗口的长度，j为滑动窗口在原始信道频域响应矩阵中滑动的次数，l和j均为正整数。
[0132]
具体来说，获取到携带时偏的原始信道频域响应矩阵后，根据滑动窗口构造第一数据矩阵。
[0133]
选取某个天线阵列中某通道的原始信道频域响应矩阵为：
[0134]
h＝[x(0)，x(1)，x(2)，
…
，x(n
f-1)]
t
[0135]
其中x(k)表示第k个频点的信道系数，即第k个采样点的信道系数，nf表示总频点数，k＝0，1，2，
…
，n
f-1。
[0136]
滑动窗口为：xk＝[x(k)，x(k+1)，
…
，x(k+l-1)]
t
，
[0137]
其中，l为滑动窗口的长度，l的取值一般为不超过0.5(n
f-1)的最大整数。将xk称为观测向量，滑动窗口尽可能在全带宽内滑动，获取到观测向量x0，x1，
…
，x
j-1
，将多个观测向量排列成一个维度为l
×
j的矩阵s，即第一数据矩阵，其表达式为：
[0138][0139]
其中，j＝n-l+1，表示观测向量的个数。
[0140]
步骤202、对第一数据矩阵进行滑动窗口离散傅里叶变换，确定第一子空间矩阵。
[0141]
具体来说，对第一数据矩阵逐行进行滑动窗口离散傅里叶变换，来确定第一子空间矩阵。
[0142]
首先对矩阵s的第一行进行点数为n的离散傅里叶变换(discrete fourier tranform，dft)，得到表达式：
[0143][0144]
其中，xk为离散傅里叶变换结果，l为第一数据矩阵的行序号，l＝7，2，
…
，l，上标k为频点序号，x(i)为采样信号，n为离散傅里叶变换的点数，n和j的取值相同，l、k和n均为正
整数；
[0145]
迭代公式为：
[0146]
xk(l)＝e
j2πk/n
[xk(l-1)+x(l+n-1)-x(l-1)]
[0147]
根据迭代公式快速确定矩阵s的第l(1＞1)行的离散傅里叶变换结果，逐行存储在矩阵v中，v的表达式为：
[0148][0149]
对于矩阵v，求取每一列对应的列向量的模值或者模值的平方，选取模值或者模值的平方最大的m列向量对应的子矩阵作为近似子空间矩阵u’。
[0150]
m的取值可以通过设定一个能量百分比门限，选取模值的平方和超过该门限值的最小列数作为m的取值。
[0151]
m的取值同样可以通过设定一个门限值，选取模值或者模值的平方大于该门限值的向量数作为m的取值。
[0152]
对近似子空间矩阵u’进行正交化得到第一子空间矩阵u，正交化的方式可以为施密特正交化法，本技术实施例不对此作出限制。
[0153]
步骤203、根据第一子空间矩阵进行时间偏差估计。
[0154]
具体来说，得到正交化的第一子空间矩阵u后，根据该子空间矩阵进行时偏估计。
[0155]
根据上述方法将信号子空间和噪声子空间区分开，然后通过现代谱估计方法对时偏进行估计。本技术实施例涉及到的现代谱估计方法包括music算法和esprit算法，本技术中对两种算法进行了改进。
[0156]
在得到正交化第一子空间矩阵后，可以通过本技术提出的改进的music算法或者改进的esprit算法估计时偏，或者对改进后的两种算法结合进行时偏估计。
[0157]
本技术实施例提供的时间偏差估计方法，通过对信道频域响应矩阵进行滑动窗口傅里叶变换确定子空间矩阵，显著降低了多组离散傅里叶变换的运算复杂度，同时使用子空间近似方法代替现代谱估计方法中的特征值分解方法，进一步降低了时间偏差估计过程的运算复杂度。
[0158]
可选地，对第一数据矩阵进行滑动窗口离散傅里叶变换，确定第一子空间矩阵，包括：
[0159]
对第一数据矩阵中的每一行对应的行向量进行滑动窗口离散傅里叶变换，确定第二数据矩阵；
[0160]
确定第二数据矩阵中每一列对应的列向量的模值；
[0161]
根据模值超过第一预设门限值的列向量确定第一子空间矩阵。
[0162]
具体来说，对第一数据矩阵的第一行对应的行向量进行离散傅里叶变换，对第一行以外的行向量通过滑动窗口离散傅里叶变换的递推公式确定离散傅里叶变换结果，从而确定第二数据矩阵。通过第二数据矩阵每一列对应列向量的模值或者模值的平方，确定模值超出预设门限值的列向量组成近似空间矩阵，对近似空间矩阵进行正交化确定第一子空间矩阵。
[0163]
对离散傅里叶变换的基础公式做简单变换，可以得到：
[0164][0165]
其中，xk为离散傅里叶变换结果，l为第一数据矩阵的行序号，l＝1，2，
…
，l，上标k为频点序号，x(i)为采样信号，n为离散傅里叶变换的点数，n和j的大小相同。
[0166]
对于任意时间窗l-1(1＞1)时，上述公式改写为：
[0167][0168]
l时刻，公式改写为：
[0169][0170]
令p＝i+1，则：
[0171][0172]
整理后得到递推公式：
[0173]
xk(l)＝e
j2πk/n
[xk(l-1)+x(l+n-1)-x(l-1)]
[0174]
根据递推公式确定矩阵s的除第一行以外的每一行的行向量的离散傅里叶变换结果，得到第二数据矩阵，第二数据矩阵的表达式为：
[0175][0176]
对于矩阵v，求取每一列向量的模值或者模值的平方，选取模值或者模值的平方最大的m列向量对应的子矩阵作为近似子空间矩阵u’。
[0177]
m的取值可以通过设定一个能量百分比门限，选取模值的平方和超过门限值的最小列数作为m的取值。
[0178]
m的取值同样可以通过设定一个门限值，选取模值或者模值的平方大于门限值的列向量数作为m的取值。
[0179]
对近似子空间矩阵u’进行正交化得到第一子空间矩阵u，正交化的方式可以为施密特正交化法，本技术实施例不对此作出限制。
[0180]
本技术实施例提供的时间偏差估计方法，通过对信道频域响应矩阵进行滑动窗口傅里叶变换确定子空间矩阵，显著降低了多组离散傅里叶变换的运算复杂度，从而降低了时间偏差估计的运算复杂度。
[0181]
可选地，根据第一子空间矩阵进行时间偏差估计，包括：
[0182]
根据第一子空间矩阵的除最后一行之外的前m-1个行向量和除第一行之外的后m-1个行向量确定第一矩阵；其中，m为第一子空间矩阵的行数；
[0183]
根据第一子空间矩阵的第m行对应的行向量和第一矩阵确定频谱矩阵；
[0184]
根据频谱矩阵确定基于第一信道估计周期的信道估计的时间偏差估计结果。
[0185]
具体来说，根据第一子空间矩阵的前m-1行和后m-1行确定第一矩阵；根据第一子空间矩阵的第m行对应的行向量和第一矩阵确定频谱矩阵；根据频谱矩阵进行时偏估计。
[0186]
esprit算法是利用数据协方差矩阵信号子空间的旋转不变性来估计信号参数，要求阵列的几何结构存在不变性或者通过变换获得两个以上的相同子阵。只需要得到两个子阵间的旋转不变关系φ，即可对信号参数进行估计，本技术中涉及到的信号参数是指时偏值。
[0187]
将第一子空间矩阵u的最后一行对应的行向量，记为vh。将第一子空间矩阵u的前m-1行对应的矩阵记为u
↓
，将第一子空间矩阵u的后m-1行对应的矩阵记为u
↑
，其中m为第一子空间矩阵的行数。
[0188]
根据第一子空间矩阵的前m-1行和后m-1行确定第一矩阵，第一矩阵ψ的表达式为：
[0189][0190]
根据第一矩阵计算向量其表达式为：
[0191][0192]
根据向量和第一矩阵ψ确定频谱矩阵φ，其表达式为：
[0193][0194]
在现有的esprit算法中需要进行求逆运算来确定旋转不变关系φ，使得算法的运算复杂度过高。而本技术实施例中通过向量和矩阵的乘法规避了原有的求逆运算，降低了运算的复杂度。
[0195]
确定频谱矩阵φ后，对其进行特征值分解，得到复数特征值向量λ，计算m条路径的时延，其表达式为：
[0196][0197]
其中，angle(λ)表示括号内变量相位或者弧度，δf表示信道频域响应矩阵h的相邻子载波间隔。根据上述表达式确定m条路径的时延，选取其中最小的τ即首径时延作为第一信道估计周期的时间偏差估计结果。首径是指多径环境中最快到达的子径。
[0198]
本技术实施例提供的时间偏差估计方法，通过子空间近似方法代替现有esprit算法中所采用的特征值分解方法，降低了运算复杂度，同时对原有的矩阵求逆过程进行改进，通过基于正交矩阵的求逆规避将原本复杂的矩阵求逆问题转化为相对简单的矩阵与向量乘法，进一步降低了运算复杂度。
[0199]
可选地，根据第一子空间矩阵的第m行对应的行向量和第一矩阵确定频谱矩阵的计算公式为：
[0200][0201]
其中，φ为频谱矩阵，ψ为第一矩阵，v为第一子空间矩阵的第m行对应的行向量进行共轭转置后的向量。
[0202]
具体来说，根据上一实施例中的具体描述，可以得到频谱矩阵的计算公式为：
[0203][0204]
其中，φ为频谱矩阵，ψ为第一矩阵，v为第一子空间矩阵的第m行对应的行向量进行共轭转置后的向量。
[0205]
本技术实施例提供的时间偏差估计方法，对原本的esprit算法中的矩阵求逆过程进行改进，通过基于正交矩阵的求逆规避将原本复杂的矩阵求逆问题转化为相对简单的矩阵与向量乘法，进一步降低了运算复杂度。
[0206]
可选地，根据第一子空间矩阵进行时间偏差估计，包括：
[0207]
根据第二信道估计周期的时间偏差估计结果确定第一信道估计周期的谱峰搜索区间；其中，第二信道估计周期为第一信道估计周期的前一个信道估计周期；
[0208]
在第一信道估计周期的谱峰搜索区间的长度大于第二预设门限值的情况下，对空间谱函数进行黄金分割搜索确定第一信道估计周期的时间偏差估计结果。
[0209]
具体来说，在通过sdft得到正交化的第一子空间矩阵后，根据上一信道估计周期的时偏估计结果来确定当前信道估计周期的谱峰搜索区间，在谱峰搜索区间的长度大于第二预设门限值的情况下，对空间谱函数进行黄金分割搜索确定当前信道估计周期的时偏估计结果。
[0210]
music算法是基于子空间分解的算法，利用信号子空间和噪声子空间的正交性，构建空间谱函数，通过谱峰搜索来对信号参数进行估计，本技术中涉及到的信号参数是指时偏值。
[0211]
假设当前时刻是信道预测开始后的第n(n＞1)次信道估计，在上一次信道估计中时偏估计结果为τ
n-1
，则本次信道估计中的谱峰搜索区间为[τ
n-1-t，τ
n-1
+t|，其中t的取值可由实际情况决定。
[0212]
相对于现有技术中遍历搜索来进行信号参数的估计，本技术实施例中使用黄金分割算法来估计首径时延，确定时偏值。
[0213]
黄金分割算法改变原本搜索时中间节点的确定，中间节点不再是中间值或者插值得到，而是位于黄金分割点附近。通过黄金分割算法可以迅速确定空间谱函数的极值点。
[0214]
令a＝τ
n-1-t，b＝τ
n-1
+t，通过music算法确定的空间谱函数表达式为：
[0215][0216]
其中，向量ui为信号子空间矩阵u的第i列对应的列向量，向量e
τ
为频率导向矢量，为一个行向量，其表达式为：
[0217][0218]
其中，表示第k个频点对应的频率。
[0219]
在进行黄金分割搜索之前，令即黄金分割数，初始化标志位gold-flag为0。
[0220]
若满足迭代条件|a-b|＞δτ，其中δτ为预设门限值，其取值可根据需要的收敛精度确定，进入迭代搜索，每次迭代根据gold-flag标志位执行以下步骤：
[0221]
步骤301、若gold-flag＝0，执行：τ1＝a+(1-ρ)(b-a)，τ2＝a+ρ(b-a)，根据空间谱函数计算伪谱p1＝p(τ1)，p2＝p(τ2)；跳转至步骤304；
[0222]
步骤302、若gold-flag＝1，执行：p1＝p2，τ2＝a+ρ(b-a)，根据空间谱函数更新伪谱p2＝p(τ2)；跳转至步骤304；
[0223]
步骤303、若gold-flag＝2，执行：p2＝p1，τ1＝a+(1-ρ)(b-a)，根据空间谱函数更新伪谱p1＝p(τ1)；跳转至步骤304；
[0224]
步骤304、若p1＞p2，则a＝τ1，τ1＝τ2，置gold-flag＝1；若p1≤p2，则b＝τ2，τ2＝τ1，置gold-flag＝2；进入下一次迭代；
[0225]
步骤305、当|a-b|≤δτ时，迭代收敛，则输出迭代结果τ
est
＝(a+b)/2。
[0226]
最终的τ
est
即为首径时延，也就是最终的时间偏差估计结果。
[0227]
本技术实施例提供的时间偏差估计方法，通过前一信道估计周期的时偏估计结果确定当前信道估计周期的谱峰搜索区间，缩小了搜索范围，同时通过黄金分割搜索代替原有music算法中的遍历搜索，降低了搜索范围和搜索次数，在相同精度的情况下，显著降低了运算的复杂度。
[0228]
可选地，根据频谱矩阵确定基于第一信道估计周期的信道估计的时间偏差估计结果之后，还包括：
[0229]
根据第一信道估计周期的时间偏差估计结果确定第三信道估计周期的谱峰搜索区间；其中，第三信道估计周期为第一信道估计周期的后一个信道估计周期，第一信道估计周期为信道估计开始后的第一个信道估计周期；
[0230]
在第三信道估计周期的谱峰搜索区间的长度大于第二预设门限值的情况下，对空间谱函数进行黄金分割搜索确定第三信道估计周期的时间偏差估计结果。
[0231]
具体来说，以改进后的esprit算法输出的时偏值作为改进后music算法的初值，将两种算法进行结合来确定时偏值。
[0232]
对于信道预测开始后的第一次信道估计，将第一子空间矩阵u的最后一行对应的行向量，记为vh。将第一子空间矩阵u的前m-1行对应的矩阵记为u
↓
，将第一子空间矩阵u的后
m-1行对应的矩阵记为u
↑
，其中m为第一子空间矩阵的行数。
[0233]
根据第一子空间矩阵的前m-1行和后m-1行确定第一矩阵，第一矩阵ψ的表达式为：
[0234][0235]
根据第一矩阵计算向量其表达式为：
[0236][0237]
根据向量和第一矩阵ψ确定频谱矩阵φ，其表达式为：
[0238][0239]
确定频谱矩阵φ后，对其进行特征值分解，得到复数特征值向量λ，计算m条路径的时延，其表达式为：
[0240][0241]
其中，angle(λ)表示括号内变量相位或者弧度，δf表示信道频域响应矩阵h的相邻子载波间隔。根据上述表达式确定m条路径的时延，选取其中最小的τ即首径时延作为第一次信道估计的时间偏差估计结果，记为τ1。
[0242]
对于信道预测开始后的第2次信道估计，谱峰搜索区间为[τ
1-t，τ1+t]，其中t的取值可由实际情况决定。
[0243]
令a＝τ
1-t，b＝τ1+t，通过music算法确定的空间谱函数表达式为：
[0244][0245]
其中，向量ui为信号子空间矩阵u的第i列对应的列向量，向量e
τ
为频率导向矢量，为一个行向量，其表达式为：
[0246][0247]
其中，表示第k个频点对应的频率。
[0248]
在进行黄金分割搜索之前，令即黄金分割数，初始化标志位gold-flag为0。
[0249]
若满足迭代条件|a-b|＞δτ，其中δτ为预设门限值，其取值可根据需要的收敛精度确定。进入迭代搜索，每次迭代根据gold-flag标志位执行以下步骤：
[0250]
步骤401、若gold-flag＝0，执行：τa＝a+(1-ρ)(b-a)，τb＝a+ρ(b-a)，根据空间谱函数计算伪谱pa＝p(τa)，pb＝p(τb)；跳转至步骤404；
[0251]
步骤402、若gold-fag＝1，执行：pa＝pb，τb＝a+ρ(b-a)，根据空间谱函数更新伪谱pb＝p(τb)；跳转至步骤404；
[0252]
步骤403、若gold-flag＝2，执行：pb＝pa，τa＝a+(1-ρ)(b-a)，根据空间谱函数更新伪谱pa＝p(τa)；跳转至步骤404；
[0253]
步骤404、若pa＞pb，则a＝τa，τa＝τb，置gold-flag＝1；若pa≤pb，则b＝τb，τb＝τa，置gold-flag＝2；进入下一次迭代；
[0254]
步骤405、当|a-b|≤δτ时，迭代收敛，则输出迭代结果τ2＝(a+b)/2。
[0255]
最终的τ2即为首径时延，也就是第二次信道估计中的时间偏差估计结果。
[0256]
对于之后的每一次信道估计，都根据上一次信道估计的时偏值确定谱峰搜索区间，然后重复上述操作，即可得到信道预测开始后第n(n＞2)次信道估计中的时间偏差估计结果。
[0257]
本技术实施例提供的时间偏差估计方法，以改进后的esprt算法的时间偏差估计结果作为初值，与改进后的music算法进行结合，兼顾搜索精度和搜索次数，在精度相同的情况下，进一步降低了搜索次数。
[0258]
可选地，对空间谱函数进行黄金分割搜索确定第一信道估计周期的时间偏差估计结果，包括：
[0259]
以第二预设门限值为收敛精度，根据黄金分割搜索确定空间谱函数在第一信道估计周期的谱峰搜索区间内的极小值；
[0260]
根据极小值确定目标搜索区间，第一信道估计周期的时间偏差估计结果为目标搜索区间的中值。
[0261]
具体来说，以第二预设门限值作为收敛精度，根据黄金分割搜索确定空间谱函数在谱峰搜索区间内的极小值，根据该极小值可以缩小谱峰搜索区间到收敛精度，从而确定时偏值。
[0262]
黄金分割法属于一维搜索算法中的试探法，能使相邻两次搜索区间都具有相同的缩短率。通过前一信道估计周期的时偏值来确定当前信道估计周期的谱峰搜索区间，以缩小搜索范围。
[0263]
通过黄金分割法来确定空间谱函数在谱峰搜索区间内的极小值，对于空间谱函数p(τ)，谱峰搜索区间为[a，b]，第二预设门限值及收敛精度为δτ，其中δτ为预设门限值，其取值可根据需要的收敛精度确定。
[0264]
在进行黄金分割搜索之前，令即黄金分割数，初始化标志位gold-flag为0。
[0265]
若满足迭代条件|a-b|＞δτ，进入迭代搜索，每次迭代根据gold-flag标志位执行以下步骤：
[0266]
步骤501、若gold-flag＝0，执行：τ1＝a+(1-ρ)(b-a)，τ2＝a+ρ(b-a)，根据空间谱函数计算伪谱p1＝p(τ1)，p2＝p(τ2)；跳转至步骤504；
[0267]
步骤502、若gold-flag＝1，执行：p1＝p2，τ2＝a+ρ(b-a)，根据空间谱函数更新伪谱p2＝p(τ2)；跳转至步骤504；
[0268]
步骤503、若gold-flag＝2，执行：p2＝p1，τ1＝a+(1-ρ)(b-a)，根据空间谱函数更新伪谱p1＝p(τ1)；跳转至步骤504；
[0269]
步骤504、若p1＞p2，则a＝τ1，τ1＝τ2，置gold-flag＝1；若p1≤p2，则b＝τ2，τ2＝τ1，置gold-flag＝2；进入下一次迭代；
[0270]
步骤505、当|a-b|≤δτ时，迭代收敛，则输出迭代结果τ
est
＝(a+b)/2。
[0271]
最终的τ
est
即为首径时延，也就是最终的时间偏差估计结果。
[0272]
本技术实施例提供的时间偏差估计方法，通过黄金分割搜索代替原本music算法中的遍历搜索，显著减少了搜索次数，在相同精度的情况下，降低了时间偏差估计过程中的运算复杂度。
[0273]
下面结合一个具体的实施例对本技术提供的时间偏差估计方法进行说明，图3是本技术实施例提供的时间偏差估计方法的流程示意图之二，如图3所示，在获取到携带时偏的原始信道频域响应矩阵h后，通过本技术提供的子空间近似方法将信号子空间与噪声子空间区分开来。
[0274]
通过原始信道频域响应矩阵h构造第一数据矩阵，对第一数据矩阵进行sdft得到第二数据矩阵。取第二数据矩阵中模值超过预设门限值的列向量构成近似子空间矩阵，对近似子空间矩阵进行正交化得到子空间矩阵。
[0275]
通过改进后的esprit算法确定信道预测中第一次信道估计过程中的时偏值，基于该时偏值确定第二次信道估计过程中的谱峰搜索区间，通过改进后的music算法进行谱峰搜索，来确定第二次信道估计过程中的时偏值。
[0276]
之后每一次时偏估计都可以采用改进后的music算法进行，以前一次信道估计的时偏值确定当前信道估计中的谱峰搜索区间，基于黄金分割搜索对music算法确定的空间谱函数在对应的谱峰搜索区间确定极小值，从而不断缩小谱峰搜索区间至收敛精度，以缩小后的搜索区间的中值作为当前信道估计的时偏值。
[0277]
基于每一次信道估计过程中的时偏值对原始信道频域响应矩阵进行校准。
[0278]
图4是本技术实施例提供的网络设备的结构示意图，如图4所示，本技术实施例提供一种网络设备，包括存储器420，收发机400，处理器410：
[0279]
存储器420，用于存储计算机程序；收发机400，用于在处理器的控制下收发数据；处理器410，用于读取存储器中的计算机程序并执行以下操作：
[0280]
根据滑动窗口确定基于第一信道估计周期的信道估计的原始信道频域响应矩阵对应的第一数据矩阵，第一数据矩阵为l
×
j的矩阵，l为滑动窗口的长度，j为滑动窗口在原始信道频域响应矩阵中滑动的次数，l和j均为正整数；
[0281]
对第一数据矩阵进行滑动窗口离散傅里叶变换，确定第一子空间矩阵；
[0282]
根据第一子空间矩阵进行时间偏差估计。
[0283]
具体来说，收发机400，用于在处理器410的控制下接收和发送数据。
[0284]
其中，在图4中，总线架构可以包括任意数量的互联的总线和桥，具体由处理器410代表的一个或多个处理器和存储器420代表的存储器的各种电路链接在一起。总线架构还可以将诸如外围设备、稳压器和功率管理电路等之类的各种其他电路链接在一起，这些都是本领域所公知的，因此，本文不再对其进行进一步描述。总线接口提供接口。收发机400可以是多个元件，即包括发送机和接收机，提供用于在传输介质上与各种其他装置通信的单元，这些传输介质包括无线信道、有线信道、光缆等传输介质。处理器410负责管理总线架构和通常的处理，存储器420可以存储处理器410在执行操作时所使用的数据。
[0285]
处理器410可以是中央处埋器(cpu)、专用集成电路(application specific integrated circuit，asic)、现场可编程门阵列(field-programmable gate array，fpga)
或复杂可编程逻辑器件(complex programmable logic device，cpld)，处理器也可以采用多核架构。
[0286]
可选地，对第一数据矩阵进行滑动窗口离散傅里叶变换，确定第一子空间矩阵，包括：
[0287]
对第一数据矩阵中的每一行对应的行向量进行滑动窗口离散傅里叶变换，确定第二数据矩阵；
[0288]
确定第二数据矩阵中每一列对应的列向量的模值；
[0289]
根据模值超过第一预设门限值的列向量确定第一子空间矩阵。
[0290]
可选地，对第一数据矩阵中的每一行对应的行向量进行滑动窗口离散傅里叶变换，确定第二数据矩阵，包括：
[0291]
根据公式一对第一数据矩阵的第一行对应的行向量进行点数为n的离散傅里叶变换，公式一的表达式为：
[0292][0293]
其中，xk为离散傅里叶变换结果，l为第一数据矩阵的行序号，k为频点序号，x(i)为采样信号，n为离散傅里叶变换的点数，n和j的取值相同；l、k和n均为正整数；
[0294]
在l大于1的情况下，根据公式二确定第一数据矩阵的第l行对应的行向量的离散傅里叶变换结果，公式二的表达式为：
[0295]
xk(l)＝e
j2πk/n
[xk(l-1)+x(l+n-1)-x(l-1)]
[0296]
根据公式一和公式二确定第二数据矩阵，第二数据矩阵的表达式为：
[0297][0298]
其中，v为第二数据矩阵。
[0299]
可选地，根据第一子空间矩阵进行时间偏差估计，包括：
[0300]
根据第一子空间矩阵的除最后一行之外的前m-1个行向量和除第一行之外的后m-1个行向量确定第一矩阵；其中，m为第一子空间矩阵的行数；
[0301]
根据第一子空间矩阵的第m行对应的行向量和第一矩阵确定频谱矩阵；
[0302]
根据频谱矩阵确定基于第一信道估计周期的信道估计的时间偏差估计结果。
[0303]
可选地，根据第一子空间矩阵的第m行对应的行向量和第一矩阵确定频谱矩阵的计算公式为：
[0304][0305]
其中，φ为频谱矩阵，ψ为第一矩阵，v为第一子空间矩阵的第m行对应的行向量进行共轭转置后的向量。
[0306]
可选地，根据第一子空间矩阵进行时间偏差估计，包括：
[0307]
根据第二信道估计周期的时间偏差估计结果确定第一信道估计周期的谱峰搜索区间；其中，第二信道估计周期为第一信道估计周期的前一个信道估计周期；
[0308]
在第一信道估计周期的谱峰搜索区间的长度大于第二预设门限值的情况下，对空间谱函数进行黄金分割搜索确定第一信道估计周期的时间偏差估计结果。
[0309]
可选地，根据频谱矩阵确定基于第一信道估计周期的信道估计的时间偏差估计结果之后，还包括：
[0310]
根据第一信道估计周期的时间偏差估计结果确定第三信道估计周期的谱峰搜索区间；其中，第三信道估计周期为第一信道估计周期的后一个信道估计周期，第一信道估计周期为信道估计开始后的第一个信道估计周期；
[0311]
在第三信道估计周期的谱峰搜索区间的长度大于第二预设门限值的情况下，对空间谱函数进行黄金分割搜索确定第三信道估计周期的时间偏差估计结果。
[0312]
可选地，对空间谱函数进行黄金分割搜索确定第一信道估计周期的时间偏差估计结果，包括：
[0313]
以第二预设门限值为收敛精度，根据黄金分割搜索确定空间谱函数在第一信道估计周期的谱峰搜索区间内的极小值；
[0314]
根据极小值确定目标搜索区间，第一信道估计周期的时间偏差估计结果为目标搜索区间的中值。
[0315]
在此需要说明的是，本发明实施例提供的上述网络设备，能够实现上述方法实施例所实现的所有方法步骤，且能够达到相同的技术效果，在此不再对本实施例中与方法实施例相同的部分及有益效果进行具体赘述。
[0316]
图5是本技术实施例提供的时间偏差估计装置的结构示意图，如图5所示，本技术实施例提供一种时间偏差估计装置，包括：
[0317]
第一确定单元501，用于根据滑动窗口确定基于第一信道估计周期的信道估计的原始信道频域响应矩阵对应的第一数据矩阵，第一数据矩阵为l
×
j的矩阵，l为滑动窗口的长度，j为滑动窗口在原始信道频域响应矩阵中滑动的次数，l和j均为正整数；
[0318]
第二确定单元502，用于对第一数据矩阵进行滑动窗口离散傅里叶变换，确定第一子空间矩阵；
[0319]
估计单元503，用于根据第一子空间矩阵进行时间偏差估计。
[0320]
可选地，第二确定单元包括第一确定模块、第二确定模块和第三确定模块；
[0321]
第一确定模块用于对第一数据矩阵中的每一行对应的行向量进行滑动窗口离散傅里叶变换，确定第二数据矩阵；
[0322]
第二确定模块用于确定第二数据矩阵中每一列对应的列向量的模值；
[0323]
第三确定模块用于根据模值超过第一预设门限值的列向量确定第一子空间矩阵。
[0324]
可选地，第一确定模块包括第一确定子模块、第二确定子模块和第三确定子模块；
[0325]
第一确定子模块用于根据公式一对第一数据矩阵的第一行对应的行向量进行点数为n的离散傅里叶变换，公式一的表达式为：
[0326]
[0327]
其中，xk为离散傅里叶变换结果，l为第一数据矩阵的行序号，k为频点序号，x(i)为采样信号，n为离散傅里叶变换的点数，n和j的取值相同；l、k和n均为正整数；
[0328]
第二确定子模块用于在l大于1的情况下，根据公式二确定第一数据矩阵的第l行对应的行向量的离散傅里叶变换结果，公式二的表达式为：
[0329]
xk(l)＝e
j2πk/n
[xk(l-1)+x(l+n-1)-x(l-1)]
[0330]
第三确定子模块用于根据公式一和公式二确定第二数据矩阵，第二数据矩阵的表达式为：
[0331][0332]
其中，v为第二数据矩阵。
[0333]
可选地，估计单元包括第四确定模块、第五确定模块和第六确定模块；
[0334]
第四确定模块用于根据第一子空间矩阵的除最后一行之外的前m-1个行向量和除第一行之外的后m-1个行向量确定第一矩阵；其中，m为第一子空间矩阵的行数；
[0335]
第五确定模块用于根据第一子空间矩阵的第m行对应的行向量和第一矩阵确定频谱矩阵；
[0336]
第六确定模块用于根据频谱矩阵确定基于第一信道估计周期的信道估计的时间偏差估计结果。
[0337]
可选地，根据第一子空间矩阵的第m行对应的行向量和第一矩阵确定频谱矩阵的计算公式为：
[0338][0339]
其中，φ为频谱矩阵，ψ为第一矩阵，v为第一子空间矩阵的第m行对应的行向量进行共轭转置后的向量。
[0340]
可选地，估计模块包括第七确定模块和第八确定模块；
[0341]
第七确定模块用于根据第二信道估计周期的时间偏差估计结果确定第一信道估计周期的谱峰搜索区间；其中，第二信道估计周期为第一信道估计周期的前一个信道估计周期；
[0342]
第八确定模块用于在第一信道估计周期的谱峰搜索区间的长度大于第二预设门限值的情况下，对空间谱函数进行黄金分割搜索确定第一信道估计周期的时间偏差估计结果。
[0343]
可选地，估计单元还包括第九确定模块和第十确定模块；
[0344]
第九确定模块用于根据第一信道估计周期的时间偏差估计结果确定第三信道估计周期的谱峰搜索区间；其中，第三信道估计周期为第一信道估计周期的后一个信道估计周期，第一信道估计周期为信道估计开始后的第一个信道估计周期；
[0345]
第十确定模块用于在第三信道估计周期的谱峰搜索区间的长度大于第二预设门限值的情况下，对空间谱函数进行黄金分割搜索确定第三信道估计周期的时间偏差估计结
果。
[0346]
可选地，第八确定模块包括第四确定子模块和第五确定子模块；
[0347]
第四确定子模块用于以第二预设门限值为收敛精度，根据黄金分割搜索确定空间谱函数在第一信道估计周期的谱峰搜索区间内的极小值；
[0348]
第五确定子模块用于根据极小值确定目标搜索区间，第一信道估计周期的时间偏差估计结果为目标搜索区间的中值。
[0349]
在此需要说明的是，本发明实施例提供的上述时间偏差估计装置，能够实现上述以基站为执行主体的时间偏差估计方法实施例所实现的所有方法步骤，且能够达到相同的技术效果，在此不再对本实施例中与方法实施例相同的部分及有益效果进行具体赘述。
[0350]
需要说明的是，本技术实施例中对单元/模块的划分是示意性的，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式。另外，在本技术各个实施例中的各功能单元/模块可以集成在一个处理单元/模块中，也可以是各个单元/模块单独物理存在，也可以两个或两个以上单元/模块集成在一个单元/模块中。上述集成的单元/模块既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
[0351]
集成的单元/模块如果以软件功能单元/模块的形式实现并作为独立的产品销售或使用时，可以存储在一个处理器可读取存储介质中。基于这样的理解，本技术的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)或处理器(processor)执行本技术各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：u盘、移动硬盘、只读存储器(read-only memory，rom)、随机存取存储器(random access memory，ram)、磁碟或者光盘等各种可以存储程序代码的介质。
[0352]
另一方面，本技术实施例还提供一种计算机可读存储介质，计算机可读存储介质存储有计算机程序，计算机程序用于使处理器执行上述各实施例提供的时间偏差估计方法，包括：
[0353]
根据滑动窗口确定基于第一信道估计周期的信道估计的原始信道频域响应矩阵对应的第一数据矩阵，第一数据矩阵为l
×
j的矩阵，l为滑动窗口的长度，j为滑动窗口在原始信道频域响应矩阵中滑动的次数，l和j均为正整数；对第一数据矩阵进行滑动窗口离散傅里叶变换，确定第一子空间矩阵；根据第一子空间矩阵进行时间偏差估计。
[0354]
计算机可读存储介质可以是处理器能够存取的任何可用介质或数据存储设备，包括但不限于磁性存储器(例如软盘、硬盘、磁带、磁光盘(mo)等)、光学存储器(例如cd、dvd、bd、hvd等)、以及半导体存储器(例如rom、eprom、eeprom、非易失性存储器(nand flash)、固态硬盘(ssd))等。
[0355]
需要说明的是：计算机可读存储介质可以是处理器能够存取的任何可用介质或数据存储设备，包括但不限于磁性存储器(例如软盘、硬盘、磁带、磁光盘(mo)等)、光学存储器(例如cd、dvd、bd、hvd等)、以及半导体存储器(例如rom、eprom、eeprom、非易失性存储器(nand flash)、固态硬盘(ssd))等。
[0356]
另一方面，本技术实施例还提供一种包含指令的计算机程序产品，该计算机程序产品在计算机上运行时，使得计算机执行上述各实施例提供的时间偏差估计方法。
[0357]
另一方面，本技术实施例还提供一种芯片，包括处理器。该处理器用于读取并运行存储器中存储的计算机程序，以执行本技术实施例提供的时间偏差估计方法中由终端设备执行的相应操作和/或流程。可选地，该芯片还包括存储器，该存储器与该处理器通过电路或电线与存储器连接，处理器用于读取并执行该存储器中的计算机程序。进一步地，该芯片还包括通信接口，处理器与该通信接口连接。通信接口用于接收需要处理的数据和/或信息，处理器从该通信接口获取该数据和/或信息，并对该数据和/信息进行处理。该通信接口可以是输入输出接口。
[0358]
另外需要说明的是：本技术实施例中术语“和/或”，描述关联对象的关联关系，表示可以存在三种关系，例如，a和/或b，可以表示：单独存在a，同时存在a和b，单独存在b这三种情况。字符“/”一般表示前后关联对象是一种“或”的关系。
[0359]
本技术的说明书和权利要求书中的术语“第一”、“第二”、“目标”等是用于区别类似的对象，而不用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便本技术的实施例能够以除了在这里图示或描述的那些以外的顺序实施，且“第一”、“第二”、“目标”等所区分的对象通常为一类，并不限定对象的个数，例如第一对象可以是一个，也可以是多个。
[0360]
本技术实施例中术语“多个”是指两个或两个以上，其它量词与之类似。
[0361]
本技术实施例提供的技术方案可以适用于多种系统，尤其是5g系统。例如适用的系统可以是全球移动通讯(global system of mobile communication，gsm)系统、码分多址(code division multiple access，cdma)系统、宽带码分多址(wideband code division multiple access，wcdma)通用分组无线业务(general packet radio service，gprs)系统、长期演进(long term evolution，lte)系统、lte频分双工(frequency division duplex，fdd)系统、lte时分双工(time division duplex，tdd)系统、高级长期演进(long term evolution advanced，lte-a)系统、通用移动系统(universal mobile telecommunication system，umts)、全球互联微波接入(worldwide interoperability for microwave access，wimax)系统、5g新空口(new radio，nr)系统等。这多种系统中均包括终端设备和网络设备。系统中还可以包括核心网部分，例如演进的分组系统(evloved packet system，eps)、5g系统(5gs)等。
[0362]
本技术实施例涉及的终端设备，可以是指向用户提供语音和/或数据连通性的设备，具有无线连接功能的手持式设备、或连接到无线调制解调器的其他处理设备等。在不同的系统中，终端设备的名称可能也不相同，例如在5g系统中，终端设备可以称为用户设备(user equipment，ue)。无线终端设备可以经无线接入网(radio access network，ran)与一个或多个核心网(core network，cn)进行通信，无线终端设备可以是移动终端设备，如移动电话(或称为“蜂窝”电话)和具有移动终端设备的计算机，例如，可以是便携式、袖珍式、手持式、计算机内置的或者车载的移动装置，它们与无线接入网交换语言和/或数据。例如，个人通信业务(personal communication service，pcs)电话、无绳电话、会话发起协议(session initiated protocol，sip)话机、无线本地环路(wireless local loop，wll)站、个人数字助理(personal digital assistant，pda)等设备。无线终端设备也可以称为系统、订户单元(subscriber unit)、订户站(subscriber station)，移动站(mobile station)、移动台(mobile)、远程站(remote station)、接入点(access point)、远程终端
设备(remote terminal)、接入终端设备(access terminal)、用户终端设备(user terminal)、用户代理(useragent)、用户装置(user device)，本技术实施例中并不限定。
[0363]
本技术实施例涉及的网络设备，可以是基站，该基站可以包括多个为终端提供服务的小区。根据具体应用场合不同，基站又可以称为接入点，或者可以是接入网中在空中接口上通过一个或多个扇区与无线终端设备通信的设备，或者其它名称。网络设备可用于将收到的空中帧与网际协议(internet protocol，ip)分组进行相互更换，作为无线终端设备与接入网的其余部分之间的路由器，其中接入网的其余部分可包括网际协议(ip)通信网络。网络设备还可协调对空中接口的属性管理。例如，本技术实施例涉及的网络设备可以是全球移动通信系统(global system for mobile communications，gsm)或码分多址接入(code division multiple access，cdma)中的网络设备(basetransceiver station，bts)，也可以是带宽码分多址接入(wide-bandcode division multiple access，wcdma)中的网络设备(nodeb)，还可以是长期演进(long term evolution，lte)系统中的演进型网络设备(evolutional node b，enb或e-nodeb)、5g网络架构(nextgeneration system)中的5g基站(gnb)，也可以是家庭演进基站(homeevolved node b，henb)、中继节点(relaynode)、家庭基站(femto)、微微基站(pico)等，本技术实施例中并不限定。在一些网络结构中，网络设备可以包括集中单元(centralized unit，cu)节点和分布单元(distributed unit，du)节点，集中单元和分布单元也可以地理上分开布置。
[0364]
网络设备与终端设备之间可以各自使用一根或多根天线进行多输入多输出(multi input multi output，mimo)传输，mimo传输可以是单用户mimo(single user mimo，su-mimo)或多用户mimo(multiple user mimo，mu-mimo)。根据根天线组合的形态和数量，mimo传输可以是2d-mimo、3d-mimo、fd-mimo或massive-mimo，也可以是分集传输或预编码传输或波束赋形传输等。
[0365]
本领域内的技术人员应明白，本技术的实施例可提供为方法、系统、或计算机程序产品。因此，本技术可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本技术可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器和光学存储器等)上实施的计算机程序产品的形式。
[0366]
本技术是参照根据本技术实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机可执行指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机可执行指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0367]
这些处理器可执行指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的处理器可读存储器中，使得存储在该处理器可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0368]
这些处理器可执行指令也可装载到计算机或其他可编程数据处理设备上，使得在
计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0369]
显然，本领域的技术人员可以对本技术进行各种改动和变型而不脱离本技术的精神和范围。这样，倘若本技术的这些修改和变型属于本技术权利要求及其等同技术的范围之内，则本技术也意图包含这些改动和变型在内。