CN101917363B - 一种多普勒频移的估算方法、装置 - Google Patents

Abstract

本发明实施例公开了一种多普勒频移的估算方法、装置，其中方法的实现包括：对同一长期演进LTE上行子帧的两个导频符号进行频域信道响应估计，得到所述两个导频符号的频域信道响应值；根据所述频域信道响应值，获取所述子帧内相差设定距离的两个解调参考信号DMRS的导频符号的频域信道响应的互相关值，获取所述子帧内相差所述设定距离的DMRS的导频符号的频域信道响应的自相关值；获取所述互相关值与自相关值的比值，根据所述比值获取多普勒频移值。采用对不同频率索引位置上的频域信道响应值进行自相关计算，可以消除噪声影响，提升了估算精度。对多普勒频移的估计不会局限于需要同样位置子载波上的导频符号，扩大了方法的应用范围。

Description

一种多普勒频移的估算方法、装置

技术领域

本发明涉及通信技术领域，特别涉及一种多普勒频移的估算方法、装置。

背景技术

在基于正交频分复用(Orthogonal Frequency DiVision Multiplex，OFDM)的无线移动系统中，可以利用信道的时变性引起的多普勒频移参数测量，在物理层进行物理层算法的优化，提高接收机的接收性能。在高层(物理层以上的层)进行资源调度的优化，从而最终达到提高系统容量。

受环境的影响，电磁波的信号采用散射的方式传播，信号经不同传播路径到达接收机。接收机接收到的信号是从多个路径收到的信号的矢量和。由于接收机的运动以及接收机周围地形的时变性，接收机收到的信号也呈现时变的特性。假定在移动的接收机的周围分布着密集的随机指向发射机阵列，电磁波的各个散射分量以相同的幅度到达接收机；另，由于各个散射分量的入射角度和相位彼此不同，现假定各个散射分量相互独立且均匀分布；于是所有散射分量的矢量和服从复高斯分布，幅度服从瑞利分布，这就是无线衰弱信道的克拉克(Clarke)模型。Clarke模型的多普勒功率谱及时域自相关函数如下：

$S\left(f\right)=\frac{1}{\mathrm{\π}{f}_d\sqrt{1-{\left(\frac{f-f_c}{f_d}\right)}^2}}$

r_h(Δt)＝J₀(2πf_dΔt)，其中f_d为最大多普勒频移，f_c为载波频率，J₀(2πf_dΔt)为第一类零阶贝塞尔(Bessel)函数。因此可以通过信道的时域自相关特性估计(或称估算)出无线信道的多普勒频移。

在基于OFDM的无线通信系统中，可以利用信道响应的频域自相关特性估算多普勒频移，多普勒频移参数估计的方案原理如下：

假定频域信道响应其中k为子载波索引，t为时间点，N为OFDM系统中快速傅氏变换(FastFourier Transform，FFT)的大小，h_l(t)为第l径的信道时域响应。

通过计算信道响应的频域相关，可以估计得到用户的最大多普勒频移值，如下：

$R_H\left(\mathrm{\Δt}\right)=E\left(H^{*}(k,t)H(k,t+\mathrm{\Δt})\right)=$

$=E(\frac{1}{\sqrt{N}}\underset{l=0}{\overset{N-1}{\mathrm{\Σ}}}{h}_l(t+\mathrm{\Δt})e^{\frac{-j2\mathrm{\πlk}}{N}}*\frac{1}{\sqrt{N}}\underset{l=0}{\overset{N-1}{\mathrm{\Σ}}}{h_l}^{*}\left(t\right)e^{\frac{j2\mathrm{\πlk}}{N}})$

$={\mathrm{Cr}}_h\left(\mathrm{\Δt}\right)$

按照Clarke模型，r_h(Δt)＝J₀(2πf_dΔt)，其中J₀(Δt)为第一类零阶的Bessel函数。

上式的推导中，假设不同径的时域信道响应不相关，且h_l(t)为广义平稳(Wide-sense stationary，WSS)平稳随机过程。

为了消除R_H(Δt)绝对幅度的影响，可以利用频域自相关和频域功率的比值，即相对值很方便的估计出最大多普勒频移。

$\frac{R_H\left(\mathrm{\Δt}\right)}{R_H\left(0\right)}=\frac{J_0\left(2\mathrm{\π}{f}_d\mathrm{\Δt}\right)}{J_0\left(0\right)}=J_0\left(2\mathrm{\π}{f}_d\mathrm{\Δt}\right)$

Bessel函数自变量在0到2.4之间为单调递减函数，可以利用Bessel函数的这个单调特性得到J_d。

即：

$f_d=\frac{1}{2\mathrm{\π\Δt}}{J_0}^{-1}\left(\frac{R_H\left(\mathrm{\Δt}\right)}{R_H\left(0\right)}\right)---\left(1\right)$

实际应用中，通过计算同一频域位置上的频域参考信号信道响应之间的相关性，可以很方便的估计得到最大多普勒频移值。为了减少求逆运算的复杂度，上式中可以采用频移值获取法进行Bessel函数的求逆计算。

采用上述方法，多普勒频移值估算精确度差且应用范围窄。

发明内容

本发明实施例要解决的技术问题是提供一种多普勒频移的估算方法、装置，提高估算精度，扩展应用范围。

为解决上述技术问题，本发明所提供的多普勒频移的估算方法实施例可以通过以下技术方案实现：

对同一长期演进LTE上行子帧的两个导频符号进行频域信道响应估计，得到所述两个导频符号的频域信道响应值；

根据所述频域信道响应值，获取所述子帧内相差设定距离的两个解调参考信号DMRS的导频符号的频域信道响应的互相关值，获取所述子帧内相差所述设定距离的DMRS的导频符号的频域信道响应的自相关值；

获取所述互相关值与自相关值的比值，根据所述比值获取多普勒频移值。

一种多普勒频移的估算方法，包括：

对同一全球微波互联接入上行时隙信号的两个符号内的导频符号分别进行频域信道响应估计，得到所述两个符号内的导频符号的频域信道响应值；

根据频域信道响应值，获取所述时隙信号内相差设定距离的两个导频符号的频域信道响应的互相关值，获取所述时隙信号内相差所述设定距离的导频符号的频域信道响应的自相关值；

获取所述互相关值与自相关值的比值，根据所述比值获取多普勒频移值。

一种多普勒频移的估算装置，包括：

响应值估计单元，用于对同一长期演进LTE上行子帧的两个导频符号进行频域信道响应估计，得到所述两个导频符号的频域信道响应值；

相关值计算单元，用于根据所述频域信道响应值，获取所述子帧内相差设定距离的两个解调参考信号DMRS的导频符号的频域信道响应的互相关值，获取所述子帧内相差所述设定距离的DMRS的导频符号的频域信道响应的自相关值；

比值获取单元，用于获取所述互相关值与自相关值的比值；

频移值获取单元，用于根据所述比值获取多普勒频移值。

一种多普勒频移的估算装置，包括：

估计响应值单元，用于对同一全球微波互联接入上行时隙信号的两个符号内的导频符号分别进行频域信道响应估计，得到所述两个符号内的导频符号的频域信道响应值；

获取相关值单元，用于根据频域信道响应值，获取所述时隙信号内相差设定距离的两个导频符号的频域信道响应的互相关值，获取所述时隙信号内相差所述设定距离的导频符号的频域信道响应的自相关值；

获取比值单元，用于获取所述互相关值与自相关值的比值；

获取频移值单元，用于根据所述比值获取多普勒频移值。

上述技术方案具有如下有益效果：采用对不同频率索引位置上的频域信道响应值进行自相关计算，从而估计出多普勒频移值，可以消除噪声影响，提升了估算精度。对多普勒频移的估计不会局限于需要同样位置子载波上的导频符号，扩大了方法的应用范围。

附图说明

为了更清楚地说明本发明实施例的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

图1为本发明实施例方法流程示意图；

图2为现有技术的LTE子帧结构示意图；

图3为本发明实施例另一方法流程示意图；

图4为本发明实施例Wimax块格式示意图；

图5为本发明实施例装置结构示意图；

图6为本发明实施例装置结构示意图；

图7为本发明实施例另一装置结构示意图；

图8为本发明实施例另一装置结构示意图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本发明实施例提供了一种多普勒频移的估算方法，如图1所示，该方法可以应用于长期演进(Long Term Evolution，LTE)网络系统，包括：

101：对同一长期演进LTE上行子帧的两个导频符号进行频域信道响应估计，得到上述两个导频符号的频域信道响应值。

102：根据上述频域信道响应值，获取上述子帧内相差设定距离的两个解调参考信号DMRS的导频符号的频域信道响应的互相关值，获取上述子帧内相差上述设定距离的DMRS的导频符号的频域信道响应的自相关值。

上述设定的距离可以是设定的RB的个数，例如1个RB或者2个RB；设定的距离，也可以包括多个不同的RB值。

因此，获取上述子帧内相差设定距离的两个解调参考信号DMRS的导频符号的频域信道响应的互相关值，可以是获取相差1个RB的两个DMRS的导频符号的频域信道响应的互相关值，也可以是获取相差2个RB的两个DMRS的导频符号的频域信道响应的互相关值；还可以是，获取相差1个RB的两个DMRS的导频符号的频域信道响应的互相关值，以及获取相差2个RB的两个DMRS的导频符号的频域信道响应的互相关值。

上述102的计算方法可以通过式(2)进行计算：

$R_H(\mathrm{\ΔK},\mathrm{\Δt})=E(H^{*}(k,t)H(k+\mathrm{\ΔK},t+\mathrm{\Δt})=$

$=E(\frac{1}{\sqrt{N}}\underset{l=0}{\overset{N-1}{\mathrm{\Σ}}}{h}_l(t+\mathrm{\Δt})e^{\frac{-j2\mathrm{\πlk}+\mathrm{\ΔK})}{N}}*\frac{1}{\sqrt{N}}\underset{l=0}{\overset{N-1}{\mathrm{\Σ}}}{h_l}^{*}\left(t\right)e^{\frac{j2\mathrm{\πlk}}{N}})$

$=E(\frac{1}{\sqrt{N}}\underset{l=0}{\overset{N-1}{\mathrm{\Σ}}}{h}_l(t+\mathrm{\Δt})e^{\frac{-j2\mathrm{\πl}\mathrm{\ΔK})}{N}}*\frac{1}{\sqrt{N}}\underset{l=0}{\overset{N-1}{\mathrm{\Σ}}}{h_l}^{*}\left(t\right)e^{\frac{j2\mathrm{\πlk}}{N}})$

$={\mathrm{Cr}}_h\left(\mathrm{\Δt}\right)*r_H\left(\mathrm{\ΔK}\right)---\left(2\right)$

式(2)中，ΔK为两个子载波在频域上相隔的子载波个数，Δt为两个子载波在时域上相隔的时间值。r_h(Δt)为时域自相关函数，r_H(ΔK)为频域自相关函数。C为常数，N为OFDM系统中快速傅氏变换(Fast Fourier Transform，FFT)的大小。DMRS为导频，由若干个子载波组成。自相关和互相关是统计值，由若干个子载波求平均得到。由于Δt＝0时互相关值与自相关值相等，令Δt＝0，根据上式可以计算自相关值。

103：计算上述互相关值与自相关值的比值，根据上述比值获取多普勒频移值。根据上述比值获取多普勒频移值的方法很多，例如查贝塞尔表等，对此本发明实施例不予限定。

上述102中设定距离可以为至少两个不同的值，则103计算上述互相关值与自相关值的比值包括：对应上述至少两个不同值计算各值对应的互相关值与自相关值的比值，计算上述各值对应的互相关值与自相关值的比值的平均值。对此，在本实施例的具体举例中有详细说明。

上述103的计算实现方式可以是通过计算式(3)然后查贝塞尔表得到多普勒频移值：

$f_d=\frac{1}{2\mathrm{\π\Δt}}{J_0}^{-1}\left(\frac{R_H(\mathrm{\ΔK},\mathrm{\Δt})}{R_H\left(\mathrm{\ΔK}\right)}\right)---\left(3\right)$

其中，J₀ ^-1J0-1是0阶贝塞尔函数，R_H(ΔK)为频域自相关函数。

如表1所示，为贝塞尔表的一个示例：

表1贝塞尔表示例

进一步地，可以对使用101～103的方法得到上述上行子帧属于同一LTE帧内的至少两个上行无线子帧的多普勒频移值求多普勒频移值的平均值。

可以分别计算出导频符号上子载波索引分隔不同大小的多普勒频移值，之后再进行联合判决。如计算出导频符号上相隔ΔK₁和ΔK₂时的多普勒频移值：

然后对这两个值进行联合判决得到第二多普勒频移值，可以更精确地估计出多普勒频移值。上述联合判决可以是：计算上述至少两个子帧的多普勒频移值的平均值。

采用本发明实施例方法，因为采用对不同频率索引位置上的频域信道响应值进行自相关计算，从而估计出多普勒频移值，可以消除噪声影响，提升了估算精度，

理由如下：由下式可见，噪声对自相关参数计算没有影响，假设不同频域上的干扰噪声功率不相关，即：

${\tilde{R}}_H(\mathrm{\ΔK},\mathrm{\Δt})=E(H^{*}(k,t)+N(k,t))(H(k+\mathrm{\ΔK},t+\mathrm{\Δt})+N(k+\mathrm{\ΔK},t+\mathrm{\Δt})=E(H^{*}(k,t)H(k+\mathrm{\ΔK},t+\mathrm{\Δt})$

$=R_H(\mathrm{\ΔK},\mathrm{\Δt})$

同理：

由此可见，对多普勒频移的估计不会局限于频域上同样位置子载波上的导频符号，扩大了方法的应用范围。

另外，本发明实施例使用相隔不同频率位置上的频域信道响应值进行互相关，并分别进行判决，最后对得到的每个多普勒频移值合并，进行联合判决，使得不需要太多的样点数即可获得理想的统计特性。

以下是上述方法实现中一个具体的举例：

以LTE系统中，由网络侧估计用户设备(User Equipment，UE)的多普勒频移为例，如图2所示，为LTE无线子帧的结构示意图，其中，

为每个上行时隙里包含的SC-FDMA符号数目；

为每个上行时隙里每个频域RB的大小，以每个资源块里包含的子载波数表示；为上行带宽配置，以

的倍数来表示；资源要素Resource element(k，l)中k为子载波索引，l为时域符号索引；

如图2所示：一个上行时隙(One uplink slot)T_slot；SC-FDMA符号(symbols)；

子载波(subcarriers)；

子载波(subcarriers)；资源要素(Resource element)资源块(Resource block)resource elements

图2中，黑色位置为导频子载波。如果LTE无线子帧为子帧内非跳频模式，UE同一子帧内不同时隙内导频符号位置相同；如果LTE无线子帧为子帧内跳频模式，UE同一子帧内不同时隙导频符号位置不同，但RB数相同，UE不同时隙内导频符号位置相隔一定的RB间隔。

对一个子帧内的导频符号分别进行频域信道响应估计，之后利用两个时隙的导频符号的频域信道响应值进行如下步骤的多普勒频移估计：

计算每个子帧两个时隙相差N个资源块(Resource Block，RB)位置的两个DMRS导频符号频域信道响应的互相关值，其中，N可配置，即可以任意配置相差RB的个数，例如1，2，3，4等。

以下为举例说明，具体说明相差1个RB和两个RB的情况：

计算每个子帧相差1个RB位置的第二个DMRS导频符号频域信道响应的自相关值。

计算每个子帧两个时隙相差2个RB位置的两个DMRS导频符号频域信道响应的互相关值。

计算每个子帧相差2个RB位置的第二个DMRS导频符号频域信道响应的自相关值。

连续N个子帧分别计算相差1个或2个RB时的互相关值的平均值和自相关值的平均值，然后并求得两个平均值的比值。

根据上述比值查Bessel表进行联合判决得到最终的多普勒频移值。联合判决方式可以是，对得到的多个多普勒频移值合并，然后求平均值。

本发明实施例还提供了另一种多普勒频移的估算方法，该方法可以用于全球微波互联接入(Worldwide Interoperability for Microwave Access，Wimax)等采用时隙信号传输的网络系统中，本实施例中的计算自相关值、互相关值、以及比值的公式可以参考图1对应实施例的公式，如图2所示，包括：

201：对同一Wimax上行时隙信号的两个符号内的导频符号分别进行频域信道响应估计，得到上述两个符号内的导频符号的频域信道响应值。

202：根据频域信道响应值，获取上述时隙信号内相差设定距离的两个导频符号的频域信道响应的互相关值，获取上述时隙信号内相差上述设定距离的导频符号的频域信道响应的自相关值。

上述202中，可以是获取至少两个上述时隙信号内相差设定距离的两个导频符号的频域信道响应的互相关值，获取至少两个上述时隙信号内相差上述设定距离的导频符号的频域信道响应的自相关值。

203：计算上述互相关值与自相关值的比值，根据上述比值获取多普勒频移值。上述根据比值获取多普勒频移值的方法可以是查贝塞尔表，对于根据比值获取频移值的具体方法本发明实施例不予限定。

若在202中计算了至少两个自相关值和互相关值，则203可以计算上述至少两个互相关值的和，计算上述至少两个自相关值的和，计算上述互相关值的和与上述自相关值的和的比值。

进一步地，在202中设定距离为至少两个不同值，则上述计算上述互相关值与自相关值的比值包括：对应上述不同值计算各值对应的互相关值与自相关值的比值，计算上述各值对应的互相关值与自相关值的比值的平均值；上述根据上述比值获取多普勒频移值包括：根据上述平均值获取多普勒频移值。

本实施例采用对不同频率索引位置上的频域信道响应值进行自相关计算，从而估计出多普勒频移值，所以可以消除噪声影响。对多普勒频移的估计不会局限于需要同样位置子载波上的导频符号，扩大了方法的应用范围。

另外，本发明实施例使用相隔不同频率位置上的频域信道响应值进行互相关，并分别进行判决，最后对得到的每个多普勒频移值合并，进行联合判决，使得不需要太多的样点数即可获得理想的统计特性。

以下是上述方法实现中一个具体的举例：

以Wimax上行系统中估计UE的多普勒频移为例，如图4所示为Wimax的块(tile)格式，一个tile有4个子载波，6个tile组成一个上行子信道，tile中有4个固定位置的导频子载波(pilor carrier)(图4中白色圆)，以及8个数据子载波(data carrier)(图4中黑色圆)。一个上行时隙(slot)子信道有48个数据子载波和24个导频子载波。

对一个时隙内的导频符号分别进行频域信道响应估计，之后利用两个导频符号的频域信道响应值进行如下步骤的多普勒频移估计。其中，时隙相差的子载波数量N可配置，即可以任意配置相差子载波的个数，例如4，8等，以下说明中具体说明相差4个子载波和8个子载波的情况。

计算每个时隙相差4个子载波位置的两个导频符号频域信道响应的互相关值。

计算每个时隙相差4个子载波位置的第二个导频符号频域信道响应的自相关值。

计算每个时隙相差8个子载波位置的两个导频符号频域信道响应的互相关值。

计算每个时隙相差8个子载波位置的第二个导频符号频域信道响应的自相关值。

连续N个时隙的信号分别计算相差4个或8个子载波时的互相关值的平均值和自相关值的平均值，然后并求得两个平均值的比值；

根据上述比值查Bessel表进行联合判决得到最终的多普勒频移值。判决方式可以是，对得到的多个多普勒频移值合并，然后求平均值。

本发明实施例还提供了一种多普勒频移的估算装置，如图5所示，包括：

响应值估计单元501，用于对同一长期演进LTE上行子帧的两个导频符号进行频域信道响应估计，得到上述两个导频符号的频域信道响应值；

相关值计算单元502，用于根据上述频域信道响应值，获取上述子帧内相差设定距离的两个解调参考信号DMRS的导频符号的频域信道响应的互相关值，获取上述子帧内相差上述设定距离的DMRS的导频符号的频域信道响应的自相关值；

比值获取单元503，用于获取上述互相关值与自相关值的比值；

频移值获取单元504，用于根据上述比值获取多普勒频移值。

可选地，上述相关值计算单元502，用于根据上述频域信道响应值，获取至少两个上述子帧内相差设定距离的两个解调参考信号DMRS的导频符号的频域信道响应的互相关值，获取至少两个上述子帧内相差上述设定距离的DMRS的导频符号的频域信道响应的自相关值；

上述比值获取单元503，用于计算上述至少两个互相关值的和，计算上述至少两个自相关值的和，计算上述互相关值的和与上述自相关值的和的比值。

可选地，上述设定距离为至少两个不同值；

则上述比值获取单元503，用于对应上述不同值计算各值对应的互相关值与自相关值的比值，计算上述各值对应的互相关值与自相关值的比值的平均值；

上述频移值获取单元504，用于根据上述平均值获取多普勒频移值。

进一步地，如图6所示，上述装置还包括：

平均值求取单元601，用于对上述频移值获取单元504获取的至少两个上行无线子帧的多普勒频移值求多普勒频移值的平均值。

本发明实施例还提供了另一种多普勒频移的估算装置，如图7所示，包括：

估计响应值单元701，用于对同一全球微波互联接入上行时隙信号的两个符号内的导频符号分别进行频域信道响应估计，得到上述两个符号内的导频符号的频域信道响应值；

获取相关值单元702，用于根据频域信道响应值，获取上述时隙信号内相差设定距离的两个导频符号的频域信道响应的互相关值，获取上述时隙信号内相差上述设定距离的导频符号的频域信道响应的自相关值；

获取比值单元703，用于获取上述互相关值与自相关值的比值；

获取频移值单元704，用于根据上述比值获取多普勒频移值。

可选地，上述获取相关值单元702，用于获取至少两个上述时隙信号内相差设定距离的两个导频符号的频域信道响应的互相关值，获取至少两个上述时隙信号内相差上述设定距离的导频符号的频域信道响应的自相关值；

上述获取比值单元703，用于计算上述至少两个互相关值的和，计算上述至少两个自相关值的和，计算上述互相关值的和与上述自相关值的和的比值。

可选地，上述设定距离为至少两个；

则上述获取比值单元703，用于对应上述不同值计算各值对应的互相关值与自相关值的比值，计算上述各值对应的互相关值与自相关值的比值的平均值；

上述获取频移值单元704，用于根据上述平均值获取多普勒频移值

进一步地，如图8所示，上述装置还包括：

求平均值单元801，用于对于上述上行时隙信号属于同一全球微波互联接入帧内的至少两个时隙信号的多普勒频移值求多普勒频移值的平均值。

本实施例采用对不同频率索引位置上的频域信道响应值进行自相关计算，从而估计出多普勒频移值，所以可以消除噪声影响。对多普勒频移的估计不会局限于需要同样位置子载波上的导频符号，扩大了方法的应用范围。

另外，本发明实施例使用相隔不同频率位置上的频域信道响应值进行互相关，并分别进行判决，最后对得到的每个多普勒频移值合并，进行联合判决，使得不需要太多的样点数即可获得理想的统计特性。

本领域普通技术人员可以理解实现上述实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件完成，上述的程序可以存储于一种计算机可读存储介质中，上述提到的存储介质可以是只读存储器，磁盘或光盘等。

以上对本发明实施例所提供的一种多普勒频移的估算方法、装置进行了详细介绍，本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上，本说明书内容不应理解为对本发明的限制。

Claims (16)

1.一种多普勒频移的估算方法，其特征在于，包括：

对同一长期演进LTE上行子帧的两个导频符号进行频域信道响应估计，得到所述两个导频符号的频域信道响应值；

根据所述频域信道响应值，获取所述子帧内相差设定距离的两个解调参考信号DMRS的导频符号的频域信道响应的互相关值，获取所述子帧内相差所述设定距离的DMRS的导频符号的频域信道响应的自相关值；

获取所述互相关值与自相关值的比值，根据所述比值获取多普勒频移值。

2.根据权利要求1所述方法，其特征在于，所述获取所述子帧内相差设定距离的两个解调参考信号DMRS的导频符号的频域信道响应的互相关值，获取所述子帧相差所述设定距离的DMRS的导频符号的频域信道响应的自相关值包括：

获取所述子帧内至少两个相差设定距离的两个解调参考信号DMRS的导频符号的频域信道响应的互相关值，获取所述子帧内至少两个相差所述设定距离的DMRS的导频符号的频域信道响应的自相关值；

所述获取所述互相关值与自相关值的比值包括：

计算所述至少两个互相关值的和，计算所述至少两个自相关值的和，计算所述互相关值的和与所述自相关值的和的比值。

3.根据权利要求1所述方法，其特征在于，所述设定距离为至少两个不同值；

所述获取所述互相关值与自相关值的比值包括：

对应所述至少两个不同值计算各值对应的互相关值与自相关值的比值，计算所述各值对应的互相关值与自相关值的比值的平均值；

所述根据所述比值获取多普勒频移值包括：

根据所述平均值获取多普勒频移值。

4.根据权利要求1至3任意一项所述方法，其特征在于，还包括：

对与所述上行子帧属于同一LTE帧内的至少两个上行无线子帧的多普勒频移值求多普勒频移值的平均值。

5.一种多普勒频移的估算方法，其特征在于，包括：

对同一全球微波互联接入上行时隙信号的两个符号内的导频符号分别进行频域信道响应估计，得到所述两个符号内的导频符号的频域信道响应值；

根据频域信道响应值，获取所述时隙信号内相差设定距离的两个导频符号的频域信道响应的互相关值，获取所述时隙信号内相差所述设定距离的导频符号的频域信道响应的自相关值；

获取所述互相关值与自相关值的比值，根据所述比值获取多普勒频移值。

6.根据权利要求5所述方法，其特征在于，

所述获取所述时隙信号内相差设定距离的两个导频符号的频域信道响应的互相关值，获取所述时隙信号内相差所述设定距离的导频符号的频域信道响应的自相关值包括：

获取所述时隙信号内至少两个相差设定距离的两个导频符号的频域信道响应的互相关值，获取至少两个所述时隙信号内相差所述设定距离的导频符号的频域信道响应的自相关值；

所述获取所述互相关值与自相关值的比值包括：

计算所述至少两个互相关值的和，计算所述至少两个自相关值的和，计算所述互相关值的和与所述自相关值的和的比值。

7.根据权利要求5所述方法，其特征在于，所述设定距离为至少两个不同值；

所述获取所述互相关值与自相关值的比值包括：

对应所述不同值计算各值对应的互相关值与自相关值的比值，计算所述各值对应的互相关值与自相关值的比值的平均值；

所述根据所述比值获取多普勒频移值包括：

根据所述平均值获取多普勒频移值。

8.根据权利要求5至7任意一项所述方法，其特征在于，还包括：

对于所述上行时隙信号属于同一全球微波互联接入帧内的至少两个时隙信号的多普勒频移值求多普勒频移值的平均值。

9.一种多普勒频移的估算装置，其特征在于，包括：

响应值估计单元，用于对同一长期演进LTE上行子帧的两个导频符号进行频域信道响应估计，得到所述两个导频符号的频域信道响应值；

相关值计算单元，用于根据所述频域信道响应值，获取所述子帧内相差设定距离的两个解调参考信号DMRS的导频符号的频域信道响应的互相关值，获取所述子帧内相差所述设定距离的DMRS的导频符号的频域信道响应的自相关值；

比值获取单元，用于获取所述互相关值与自相关值的比值；

频移值获取单元，用于根据所述比值获取多普勒频移值。

10.根据权利要求9所述装置，其特征在于，

所述相关值计算单元，用于根据所述频域信道响应值，获取至少两个所述子帧内相差设定距离的两个解调参考信号DMRS的导频符号的频域信道响应的互相关值，获取至少两个所述子帧内相差所述设定距离的DMRS的导频符号的频域信道响应的自相关值；

所述比值获取单元，用于计算所述至少两个互相关值的和，计算所述至少两个自相关值的和，计算所述互相关值的和与所述自相关值的和的比值。

11.根据权利要求9所述装置，其特征在于，所述设定距离为至少两个不同值；

所述比值获取单元，用于对应所述不同值计算各值对应的互相关值与自相关值的比值，计算所述各值对应的互相关值与自相关值的比值的平均值；

所述频移值获取单元，用于根据所述平均值获取多普勒频移值。

12.根据权利要求9至11任意一项所述装置，其特征在于，还包括：

平均值求取单元，用于对所述频移值获取单元获取的至少两个上行无线子帧的多普勒频移值求多普勒频移值的平均值。

13.一种多普勒频移的估算装置，其特征在于，包括：

估计响应值单元，用于对同一全球微波互联接入上行时隙信号的两个符号内的导频符号分别进行频域信道响应估计，得到所述两个符号内的导频符号的频域信道响应值；

获取相关值单元，用于根据频域信道响应值，获取所述时隙信号内相差设定距离的两个导频符号的频域信道响应的互相关值，获取所述时隙信号内相差所述设定距离的导频符号的频域信道响应的自相关值；

获取比值单元，用于获取所述互相关值与自相关值的比值；

获取频移值单元，用于根据所述比值获取多普勒频移值。

14.根据权利要求13所述装置，其特征在于，

所述获取相关值单元，用于获取至少两个所述时隙信号内相差设定距离的两个导频符号的频域信道响应的互相关值，获取至少两个所述时隙信号内相差所述设定距离的导频符号的频域信道响应的自相关值；

所述获取比值单元，用于计算所述至少两个互相关值的和，计算所述至少两个自相关值的和，计算所述互相关值的和与所述自相关值的和的比值。

15.根据权利要求13所述装置，其特征在于，所述设定距离为至少两个；

所述获取比值单元，用于对应所述不同值计算各值对应的互相关值与自相关值的比值，计算所述各值对应的互相关值与自相关值的比值的平均值；

所述获取频移值单元，用于根据所述平均值获取多普勒频移值。

16.根据权利要求13至15任意一项所述装置，其特征在于，还包括：

求平均值单元，用于对于所述上行时隙信号属于同一全球微波互联接入帧内的至少两个时隙信号的多普勒频移值求多普勒频移值的平均值。

Images