CN104853367A - 一种终端移动速度估计方法和装置 - Google Patents

Abstract

本发明提供一种终端移动速度估计方法和系统，包括：选取采样间隔内符合条件的两列频域解调参考信号位置估计得到的频域信道，计算噪声估计值；引入修正因子，结合所述噪声估计值计算所述频域信道在相应采样间隔内的自相关函数值；引入限制因子，判断所述自相关函数值的模是否小于所述限制因子，如果是，则将所述自相关函数值作为有效数据，与所有之前得到的相同采样间隔内的自相关函数值作迭代平均，求得终端移动速度；如果不是，则沿用上个周期得到的自相关函数值求得的终端移动速度。本发明适用范围广泛，可以比较稳定准确地估计出终端移动速度。

Description

一种终端移动速度估计方法和装置

技术领域

本发明涉及无线通信技术领域，尤其涉及一种终端移动速度估计方法和装置。

背景技术

长期演进（Long Term Evolution，LTE）无线通信技术以正交频分复用（Orthogonal Frequency Division Multiplexing，OFDM）技术为基础。根据OFDM标准，数据分别通过许多相互正交的子载波传输，而且有一个保护间隔，来有效克服多径干扰问题。在通信系统中，终端的快速移动会产生较大的多普勒频移，尤其在多径环境下会造成信号幅度的快速衰落和信号幅度的迅速变化，进而破坏接收信号子载波之间的正交特性，大幅度影响系统性能。

LTE系统采用多输入多输出（Multiple-Input Multiple-Output，MIMO）技术，结合高阶调制、频选调度、功率控制等技术，以达到更高的系统吞吐量和频谱效率。这些技术都需要基于事先掌握的终端无线信道信息，其中终端移动速度信息十分重要。基站端可以根据终端当前的移动速度，来对信道估计方式和信号检测相关算法进行自适应调整，完成自适应天线模式切换和自适应编码调制等；终端移动速度信息还可用于终端服务基站的切换优化，降低掉话概率；另外，还可以帮助基站优化负载均衡和协同传输等。因此，为了改善系统性能，估计终端的移动速度十分必要。

当前已存在很多速度估计方案：其中根据下行链路信道信息（主要包括秩指示RI和预编码矩阵指示PMI）变化的概率来估计终端移动速度的方法，所估速度非常粗略，并且只能在具有RI/PMI反馈的情况下才能应用。根据发射信号的采样计算复信道的自相关函数的估计值，结合发射信号到达角的角度差的估计方法极易受到噪声干扰，尤其是当终端速度比较低时，在信噪比条件不好的情况下，该方法的性能会大打折扣，而且本方法还必须得到发射信号到达角的角度差信息。

因此，根据以上描述可以看出，现有技术中的速度估计方法受噪声影响较大，估算速度粗略，易导致结果不够精确的技术问题。

发明内容

（一）要解决的技术问题

本发明提供一种终端移动速度估计方法和装置，以解决现有技术中受噪声影响大，估算结果不够精确的技术问题。

（二）技术方案

为解决上述技术问题，本发明提供一种终端移动速度估计方法，包括：

选取采样间隔内符合条件的两列频域解调参考信号位置估计得到的频域信道，计算噪声估计值；

引入修正因子，结合所述噪声估计值计算所述频域信道在相应采样间隔内的自相关函数值；

引入限制因子，判断所述自相关函数值的模是否小于所述限制因子，如果是，则将所述自相关函数值作为有效数据，与所有之前得到的相同采样间隔内的自相关函数值作迭代平均，求得终端移动速度；如果不是，则沿用上个周期得到的自相关函数值求得的终端移动速度。

进一步地，所述引入修正因子，结合所述噪声估计值计算所述频域信道在相应采样间隔内的自相关函数值包括：

利用如下公式计算采样间隔t内的自相关函数值

${\widetilde{\mathrm{\ρ}}}_t=\frac{\frac{1}{K}\underset{k=1}{\overset{K}{\mathrm{\Σ}}}({\tilde{H}}_2\left(k\right)\·{\tilde{H}}_1^{*}\left(k\right))}{\frac{1}{K}\underset{k=1}{\overset{K}{\mathrm{\Σ}}}({\tilde{H}}_1\left(k\right)\·{\tilde{H}}_1^{*}\left(k\right))-\mathrm{mod}\mathrm{ify}\_\mathrm{factor}\·P_{\mathrm{estimated}\_\mathrm{moise}}},$

其中 ${\tilde{H}}_1=\{{\tilde{H}}_1\left(k\right),k=\mathrm{1,2},...,K\},$ ${\tilde{H}}_2=\{{\tilde{H}}_2\left(k\right),k=\mathrm{1,2},...,K\}$ 分别为采样间隔t内两列频域解调参考信号位置估计得到的频域信道序列，K为分配给该终端的子载波数，modify_factor为修正因子，P_{estimated_noise}为噪声估计值。

进一步地，所述将所述自相关函数值作为有效数据，与所有之前得到的相同采样间隔内的自相关函数值作迭代平均，求得终端移动速度包括：

将第n个子帧计算得到的与前n-1个子帧做平滑，得到

${\stackrel{\‾}{\mathrm{\ρ}}}_t^n=\frac{{\widetilde{\mathrm{\ρ}}}_t^n+(n-1)\·{\stackrel{\‾}{\mathrm{\ρ}}}_t^{n-1}}{n};$

根据信道相关值与终端移动速度之间的函数映射关系，通过查表法得到所对应的终端移动速度v_t。

进一步地，

所述修正因子取值为0.12；

和/或，所述限制因子取值为1.8。

进一步地，所述方法还包括：

分别计算不同采样间隔t₁和t₂所对应的终端移动速度v_t1和v_t2，求平均得到最终估计速度：

另一方面，本发明还提供一种终端移动速度估计系统，包括：

选取单元，用于选取采样间隔内符合条件的两列频域解调参考信号位置估计得到的频域信道，计算噪声估计值；

修正单元，用于引入修正因子，结合所述噪声估计值计算所述频域信道在相应采样间隔内的自相关函数值；

限制单元，用于引入限制因子，判断所述自相关函数值的模是否小于所述限制因子，如果是，则将所述自相关函数值作为有效数据与所有之前得到的相同采样间隔内的自相关函数值作迭代平均并输出；如果不是，则输出上个周期得到的自相关函数值；

速度估计单元，用于利用所述限制单元输出的自相关函数值求得终端移动速度。

进一步地，所述修正单元还用于：

利用如下公式计算采样间隔t内的自相关函数值

${\widetilde{\mathrm{\ρ}}}_t=\frac{\frac{1}{K}\underset{k=1}{\overset{K}{\mathrm{\Σ}}}({\tilde{H}}_2\left(k\right)\·{\tilde{H}}_1^{*}\left(k\right))}{\frac{1}{K}\underset{k=1}{\overset{K}{\mathrm{\Σ}}}({\tilde{H}}_1\left(k\right)\·{\tilde{H}}_1^{*}\left(k\right))-\mathrm{mod}\mathrm{ify}\_\mathrm{factor}\·P_{\mathrm{estimated}\_\mathrm{moise}}},$

其中 ${\tilde{H}}_1=\{{\tilde{H}}_1\left(k\right),k=\mathrm{1,2},...,K\},$ ${\tilde{H}}_2=\{{\tilde{H}}_2\left(k\right),k=\mathrm{1,2},...,K\}$ 分别为采样间隔t内两列频域解调参考信号位置估计得到的频域信道序列，K为分配给该终端的子载波数，modify_factor为修正因子，P_{estimated_noise}为噪声估计值。

进一步地，

所述限制单元还用于：利用如下公式将第n个子帧计算得到的与前n-1个子帧做平滑，得到

${\stackrel{\‾}{\mathrm{\ρ}}}_t^n=\frac{{\widetilde{\mathrm{\ρ}}}_t^n+(n-1)\·{\stackrel{\‾}{\mathrm{\ρ}}}_t^{n-1}}{n};$

所述速度估计单元还用于：根据信道相关值与终端移动速度之间的函数映射关系，通过查表法得到所对应的终端移动速度vt。

进一步地，

所述修正因子取值为0.12；

和/或，所述限制因子取值为1.8。

进一步地，所述系统还包括：

平均单元，用于分别计算不同采样间隔t₁和t₂所对应的终端移动速度v_t1和v_t2，求平均得到最终估计速度：

（三）有益效果

可见，在本发明提供的一种移动速度估计方法和系统中，利用了频域解调参考信号的频域信道信息，在计算自相关函数时引入了修正因子和限制因子两个因子，能够有效地减少噪声的影响，并且使得自相关函数的计算能够适用于不同的信道条件，适用范围广泛，可以比较稳定准确地估计出终端移动速度。本发明实施例利用大间隔信道响应估计的速度对小间隔信道响应对应的速度进行修正的方式，能够进一步提高终端移动速度估计的精度。

附图说明

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

图1是本发明实施例终端移动速度估计方法的基本流程示意图；

图2是本发明实施例1终端移动速度估计方法的流程示意图；

图3是本发明实施例1中0.5ms采样间隔下频域信道抽样示意图；

图4是本发明实施例1中0.5ms和1ms采样间隔下不同速度对应的自相关函数示意图；

图5是本发明实施例1中1ms采样间隔下频域信道抽样示意图；

图6是本发明实施例采样间隔为0.5ms，4RB，SCME信道中不同SNR下的速度估计结果；

图7是本发明实施例采样间隔为0.5ms，4RB，ITU-R信道中不同SNR下的速度估计结果；

图8是本发明实施例采样间隔为0.5ms，10RB，SCME信道中不同SNR下的速度估计结果；

图9是本发明实施例采样间隔为0.5ms，10RB，ITU-R信道中不同SNR下的速度估计结果；

图10是本发明实施例采样间隔为1ms，4RB，SCME信道中不同SNR下的速度估计结果；

图11是本发明实施例采样间隔为1ms，4RB，ITU-R信道中不同SNR下的速度估计结果；

图12是本发明实施例采样间隔为1ms，10RB，SCME信道中不同SNR下的速度估计结果；

图13是本发明实施例采样间隔为1ms，10RB，ITU-R信道中不同SNR下的速度估计结果；

图14是本发明实施例终端移动速度估计系统的基本结构示意图；

图15是本发明实施例终端移动速度估计系统的一个优选结构示意图。

具体实施方式

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

本发明实施例首先提供一种终端移动速度估计方法，参见图1，包括：

步骤101：选取采样间隔内符合条件的两列频域解调参考信号位置估计得到的频域信道，计算噪声估计值。

步骤102：引入修正因子，结合所述噪声估计值计算所述频域信道在相应采样间隔内的自相关函数值。

步骤103：引入限制因子，判断所述自相关函数值的模是否小于所述限制因子，如果是，则将所述自相关函数值作为有效数据，求得终端移动速度；如果不是，则沿用上个周期得到的自相关函数值求得的终端移动速度。

可见，在本发明实施例提供的一种移动速度估计方法中，利用了频域解调参考信号的频域信道信息，在计算自相关函数时引入了修正因子和限制因子两个因子，能够有效地减少噪声的影响，并且使得自相关函数的计算能够适用于不同的信道条件，适用范围广泛，可以比较稳定准确地估计出终端移动速度。

优选地，引入修正因子，结合所述噪声估计值计算所述频域信道在相应采样间隔内的自相关函数值可以包括：

利用如下公式计算采样间隔t内的自相关函数值

${\widetilde{\mathrm{\ρ}}}_t=\frac{\frac{1}{K}\underset{k=1}{\overset{K}{\mathrm{\Σ}}}({\tilde{H}}_2\left(k\right)\·{\tilde{H}}_1^{*}\left(k\right))}{\frac{1}{K}\underset{k=1}{\overset{K}{\mathrm{\Σ}}}({\tilde{H}}_1\left(k\right)\·{\tilde{H}}_1^{*}\left(k\right))-\mathrm{mod}\mathrm{ify}\_\mathrm{factor}\·P_{\mathrm{estimated}\_\mathrm{moise}}},$

其中 ${\tilde{H}}_1=\{{\tilde{H}}_1\left(k\right),k=\mathrm{1,2},...,K\},$ ${\tilde{H}}_2=\{{\tilde{H}}_2\left(k\right),k=\mathrm{1,2},...,K\}$ 分别为采样间隔t内两列频域解调参考信号位置估计得到的频域信道序列，K为分配给该终端的子载波数，modify_factor为修正因子，P_{estimated_noise}为噪声估计值。

优选地，将所述自相关函数值作为有效数据，与所有之前得到的相同采样间隔内的自相关函数值作迭代平均，求得终端移动速度可以包括如下步骤：

将第n个子帧计算得到的与前n-1个子帧做平滑，得到

${\stackrel{\‾}{\mathrm{\ρ}}}_t^n=\frac{{\widetilde{\mathrm{\ρ}}}_t^n+(n-1)\·{\stackrel{\‾}{\mathrm{\ρ}}}_t^{n-1}}{n};$

根据信道相关值与终端移动速度之间的函数映射关系，通过查表法得到所对应的终端移动速度v_t。

优选地，修正因子的取值可以为0.12，限制因子的取值可以为1.8。在实际应用中，可以对修正因子和限制因子作出调整。在测试的过程中，当所估速度普遍高于实际速度时，应首先适当提高修正因子；如果所估速度还较高，则可以提高限制因子。反之，当所估速度普遍低于实际速度时，则应首先适当减小修正因子；如果所估速度还较低，则可以减小限制因子。

优选地，方法还可以包括：分别计算不同采样间隔t₁和t₂所对应的终端移动速度v_t1和v_t2，求平均得到最终估计速度：例如，可以首先计算0.5ms采样间隔下的信道自相关系数，另外还在1ms采样间隔的信道信息可用的情况下计算1ms采样间隔下的信道自相关系数，得到所估速度之后对0.5ms采样间隔下估计得到的速度进行修正，使得所估终端速度更加准确，这种利用大间隔信道响应估计的速度对小间隔信道响应对应的速度进行修正的方式，能够进一步提高终端移动速度估计的精度。

实施例1：

下面以在上行OFDM系统中对终端移动速度进行速度估计为例，来详细说明本发明一个实施例的具体实现过程，参见图2：

步骤201：选取0.5ms采样间隔内符合条件的两列频域解调参考信号位置估计得到的频域信道，计算噪声估计值，进入步骤202和步骤206。

本步骤中，首先判断当前子帧内目标用户是否有资源分配，若有，则选取基站端第一条接收天线上接收到的信号，抽取该子帧中DMRS位置估计得到的频域信道。如图3所示，本发明实施例1所利用的是上行OFDM系统中物理上行共享信道（Physical Uplink Shared Channe，PUSCH）上的DMRS导频信道信息，对所估DMRS信道进行连续抽样，则抽样间隔为0.5ms。由于现有通信系统在资源调度的时候在时域上是以子帧为单位的，因此每一个调度周期都可以得到间隔0.5ms的两列所估DMRS信道序列，分别表示为和K为分配给该终端的子载波数。

步骤202：引入修正因子，结合噪声估计值计算所述频域信道在0.5ms采样间隔内的自相关函数值。

本步骤中，结合噪声估计值P_{estimated_noise}，由下式计算采样间隔为0.5ms的自相关函数值：

${\widetilde{\mathrm{\ρ}}}_{0.5}=\frac{\frac{1}{K}\underset{k=1}{\overset{K}{\mathrm{\Σ}}}({\tilde{H}}_2\left(k\right)\·{\tilde{H}}_1^{*}\left(k\right))}{\frac{1}{K}\underset{k=1}{\overset{K}{\mathrm{\Σ}}}({\tilde{H}}_1\left(k\right)\·{\tilde{H}}_1^{*}\left(k\right))-\mathrm{mod}\mathrm{ify}\_\mathrm{factor}\·P_{\mathrm{estimated}\_\mathrm{moise}}},$

其中modify_factor为修正因子，建议取值为0.12。

步骤203：引入限制因子，判断自相关函数值的模是否小于限制因子，如果是，进入步骤204；如果不是，进入步骤205。

本步骤中，为了消除某次或者某几次格外离谱的数据对于所估速度的影响，对求得的做一个判断，当满足下式时：

我们认为此数据是有效数据，进入步骤204；否则舍弃，进入步骤205。在此建议限制因子limit_factor取值为1.8。

步骤204：将自相关函数值作为有效数据，与所有之前得到的相同采样间隔内的自相关函数值作迭代平均，求得终端移动速度，进入步骤213。

本步骤中，将当前第n个子帧计算得到的与前(n-1)个子帧做平滑，得到如下式所示：

${\stackrel{\‾}{\mathrm{\ρ}}}_{0.5}^n=\frac{{\widetilde{\mathrm{\ρ}}}_{0.5}^n+(n-1)\·{\stackrel{\‾}{\mathrm{\ρ}}}_{0.5}^{n-1}}{n};$

然后根据信道相关值与UE速度之间的函数映射关系，通过查表，由对应得到UE速度v_0.5，该表可以根据图4中的曲线建立，图中实线对应0.5ms采样间隔所对应的UE移动速度值，虚线对应1ms采样间隔所对应的UE移动速度值。

步骤205：沿用上个周期得到的自相关函数值，求得终端移动速度。

本步骤中，当自相关函数值的模不小于限制因子时，则认为当前周期所估移动速度与上个周期得到的移动速度相同，可以沿用上个周期得到的自相关函数值查表求得的UE速度v_0.5，进入步骤213。

步骤206：判断当前子帧与上一个子帧的DMRS是否有重复资源部分，如果有，进入步骤207，如果没有，进入步骤208。

步骤207：选取1ms采样间隔内符合条件的两列频域解调参考信号位置估计得到的频域信道，计算噪声估计值。

本步骤中，由于当前子帧与上一个子帧的DMRS具有重复资源部分，所以可以分别将这些RB上两个子帧内的第一个时隙的DMRS信号提取出来，记为R_m(k)，m＝1,3，则R₁(k)与R₃(k)间隔为1ms，抽样示意图如图5所示，进入步骤209。

步骤208：沿用上一个周期估计得到的UE速度v₁，进入步骤213。

步骤209-步骤212中对1ms采样时间间隔内终端移动速度的估计步骤，与步骤202-步骤205中针对0.5ms采样时间间隔内终端移动速度的估计步骤相同，步骤211、212得到UE速度v₁，进入步骤213。

步骤213：对v_0.5和v₁做平均，得到最终估计速度v_final，即

$v_{\mathrm{final}}=\frac{v_{0.5}+v_1}{2}$

至此，则完成了本发明实施例1对于终端移动速度估计方法的全过程。

在此对本发明实施例在TD-LTE系统上行链路级平台中仿真实现，仿真参数如表1所示。

表1仿真参数设置

对于信道模型，SCME是无直射径的信道建模，而ITU-R是有直射径的信道建模。该仿真中，修正因子取值0.12，限制因子取值1.8。仿真结果参见图6-图13。由上述附图可以看出，SNR高的时候仿真结果更加接近理论值，因此精度更高；1ms采样时间间隔下的仿真数据优于0.5ms间隔，大间隔信道响应估计的速度对小间隔信道响应的修正具有提高速度估计准确性的效果；SCME信道模型仿真结果相较ITU-R的结果更好。

本发明实施例还提供一种终端移动速度估计系统，参见图14，包括：

选取单元1401，用于选取采样间隔内符合条件的两列频域解调参考信号位置估计得到的频域信道，计算噪声估计值；

修正单元1402，用于引入修正因子，结合所述噪声估计值计算所述频域信道在相应采样间隔内的自相关函数值；

限制单元1403，用于引入限制因子，判断所述自相关函数值的模是否小于所述限制因子，如果是，则将所述自相关函数值作为有效数据与所有之前得到的相同采样间隔内的自相关函数值作迭代平均并输出；如果不是，则输出上个周期得到的自相关函数值；

速度估计单元1404，用于利用所述限制单元输出的自相关函数值求得终端移动速度。

优选地，修正单元1402还可以用于：

利用如下公式计算采样间隔t内的自相关函数值

${\widetilde{\mathrm{\ρ}}}_t=\frac{\frac{1}{K}\underset{k=1}{\overset{K}{\mathrm{\Σ}}}({\tilde{H}}_2\left(k\right)\·{\tilde{H}}_1^{*}\left(k\right))}{\frac{1}{K}\underset{k=1}{\overset{K}{\mathrm{\Σ}}}({\tilde{H}}_1\left(k\right)\·{\tilde{H}}_1^{*}\left(k\right))-\mathrm{mod}\mathrm{ify}\_\mathrm{factor}\·P_{\mathrm{estimated}\_\mathrm{moise}}},$

其中 ${\tilde{H}}_1=\{{\tilde{H}}_1\left(k\right),k=\mathrm{1,2},...,K\},$ ${\tilde{H}}_2=\{{\tilde{H}}_2\left(k\right),k=\mathrm{1,2},...,K\}$ 分别为采样间隔t内两列频域解调参考信号位置估计得到的频域信道序列，K为分配给该终端的子载波数，modify_factor为修正因子，P_{estimated_noise}为噪声估计值。

优选地，限制单元1403还可以用于：利用如下公式将第n个子帧计算得到的与前n-1个子帧做平滑，得到

${\stackrel{\‾}{\mathrm{\ρ}}}_t^n=\frac{{\widetilde{\mathrm{\ρ}}}_t^n+(n-1)\·{\stackrel{\‾}{\mathrm{\ρ}}}_t^{n-1}}{n};$

速度估计单元1404还可以用于：根据信道相关值与终端移动速度之间的函数映射关系，通过查表法得到所对应的终端移动速度v_t。

优选地，修正因子取值可以为0.12；和/或，限制因子取值可以为1.8。在实际应用中，可以对修正因子和限制因子作出调整。在测试的过程中，当所估速度普遍高于实际速度时，应首先适当提高修正因子；如果所估速度还较高，则可以提高限制因子。反之，当所估速度普遍低于实际速度时，则应首先适当减小修正因子；如果所估速度还较低，则可以减小限制因子。

优选地，系统还可以包括：平均单元1501，如图15，用于分别计算不同采样间隔t₁和t₂所对应的终端移动速度v_t1和v_t2，求平均得到最终估计速度： $v_{\mathrm{final}}=\frac{v_{t1}+v_{t2}}{2}.$

可见，本发明实施例具有如下有益效果：

在本发明实施例提供的一种移动速度估计方法和系统中，利用了频域解调参考信号的频域信道信息，在计算自相关函数时引入了修正因子和限制因子两个因子，能够有效地减少噪声的影响，并且使得自相关函数的计算能够适用于不同的信道条件，适用范围广泛，可以比较稳定准确地估计出终端移动速度。本发明实施例利用大间隔信道响应估计的速度对小间隔信道响应对应的速度进行修正的方式，能够进一步提高终端移动速度估计的精度。

最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (10)

1.一种终端移动速度估计方法，其特征在于，包括：

选取采样间隔内符合条件的两列频域解调参考信号位置估计得到的频域信道，计算噪声估计值；

引入修正因子，结合所述噪声估计值计算所述频域信道在相应采样间隔内的自相关函数值；

引入限制因子，判断所述自相关函数值的模是否小于所述限制因子，如果是，则将所述自相关函数值作为有效数据，与所有之前得到的相同采样间隔内的自相关函数值作迭代平均，求得终端移动速度；如果不是，则沿用上个周期得到的自相关函数值求得的终端移动速度。

2.根据权利要求1所述的终端移动速度估计方法，其特征在于，所述引入修正因子，结合所述噪声估计值计算所述频域信道在相应采样间隔内的自相关函数值包括：

利用如下公式计算采样间隔t内的自相关函数值

${\widetilde{\mathrm{\ρ}}}_t=\frac{\frac{1}{K}\underset{k=1}{\overset{K}{\mathrm{\Σ}}}({\tilde{H}}_2\left(k\right)\·{\tilde{H}}_1^{*}\left(k\right))}{\frac{1}{K}\underset{k=1}{\overset{K}{\mathrm{\Σ}}}({\tilde{H}}_1\left(k\right)\·{\tilde{H}}_1^{*}\left(k\right))-\mathrm{mod}\mathrm{ify}\_\mathrm{factor}\·P_{\mathrm{estimated}\_\mathrm{moise}}},$

其中 ${\tilde{H}}_1=\{{\tilde{H}}_1\left(k\right),k=\mathrm{1,2},...,K\},$ ${\tilde{H}}_2=\{{\tilde{H}}_2\left(k\right),k=\mathrm{1,2},...,K\}$ 分别为采样间隔t内两列频域解调参考信号位置估计得到的频域信道序列，K为分配给该终端的子载波数，modify_factor为修正因子，P_{estimated_noise}为噪声估计值。

3.根据权利要求2所述的终端移动速度估计方法，其特征在于，所述将所述自相关函数值作为有效数据，与所有之前得到的相同采样间隔内的自相关函数值作迭代平均，求得终端移动速度包括：

将第n个子帧计算得到的与前n-1个子帧做平滑，得到

${\stackrel{\‾}{\mathrm{\ρ}}}_t^n=\frac{{\widetilde{\mathrm{\ρ}}}_t^n+(n-1)\·{\stackrel{\‾}{\mathrm{\ρ}}}_t^{n-1}}{n};$

根据信道相关值与终端移动速度之间的函数映射关系，通过查表法得到所对应的终端移动速度v_t。

4.根据权利要求1所述的终端移动速度估计方法，其特征在于：

所述修正因子取值为0.12；

和/或，所述限制因子取值为1.8。

5.根据权利要求1至4中任一项所述的终端移动速度估计方法，其特征在于，所述方法还包括：

分别计算不同采样间隔t₁和t₂所对应的终端移动速度v_t1和v_t2，求平均得到最终估计速度：

6.一种终端移动速度估计系统，其特征在于，包括：

选取单元，用于选取采样间隔内符合条件的两列频域解调参考信号位置估计得到的频域信道，计算噪声估计值；

修正单元，用于引入修正因子，结合所述噪声估计值计算所述频域信道在相应采样间隔内的自相关函数值；

限制单元，用于引入限制因子，判断所述自相关函数值的模是否小于所述限制因子，如果是，则将所述自相关函数值作为有效数据与所有之前得到的相同采样间隔内的自相关函数值作迭代平均并输出；

如果不是，则输出上个周期得到的自相关函数值；

速度估计单元，用于利用所述限制单元输出的自相关函数值求得终端移动速度。

7.根据权利要求6所述的终端移动速度估计系统，其特征在于，所述修正单元还用于：

利用如下公式计算采样间隔t内的自相关函数值

${\widetilde{\mathrm{\ρ}}}_t=\frac{\frac{1}{K}\underset{k=1}{\overset{K}{\mathrm{\Σ}}}({\tilde{H}}_2\left(k\right)\·{\tilde{H}}_1^{*}\left(k\right))}{\frac{1}{K}\underset{k=1}{\overset{K}{\mathrm{\Σ}}}({\tilde{H}}_1\left(k\right)\·{\tilde{H}}_1^{*}\left(k\right))-\mathrm{mod}\mathrm{ify}\_\mathrm{factor}\·P_{\mathrm{estimated}\_\mathrm{moise}}},$

其中 ${\tilde{H}}_1=\{{\tilde{H}}_1\left(k\right),k=\mathrm{1,2},...,K\},$ ${\tilde{H}}_2=\{{\tilde{H}}_2\left(k\right),k=\mathrm{1,2},...,K\}$ 分别为采样间隔t内两列频域解调参考信号位置估计得到的频域信道序列，K为分配给该终端的子载波数，modify_factor为修正因子，P_{estimated_noise}为噪声估计值。

8.根据权利要求7所述的终端移动速度估计系统，其特征在于：

所述限制单元还用于：利用如下公式将第n个子帧计算得到的与前n-1个子帧做平滑，得到

${\stackrel{\‾}{\mathrm{\ρ}}}_t^n=\frac{{\widetilde{\mathrm{\ρ}}}_t^n+(n-1)\·{\stackrel{\‾}{\mathrm{\ρ}}}_t^{n-1}}{n};$

所述速度估计单元还用于：根据信道相关值与终端移动速度之间的函数映射关系，通过查表法得到所对应的终端移动速度v_t。

9.根据权利要求6所述的终端移动速度估计系统，其特征在于：

所述修正因子取值为0.12；

和/或，所述限制因子取值为1.8。

10.根据权利要求6至9中任一项所述的终端移动速度估计系统，其特征在于，所述系统还包括：

平均单元，用于分别计算不同采样间隔t₁和t₂所对应的终端移动速度v_t1和v_t2，求平均得到最终估计速度：