具体实施方式

以下描述中，为了说明而不是为了限定，提出了诸如特定系统结构、技术之类的具体细节，以便透彻理解本申请实施例。然而，本领域技术人员应当清楚，在没有这些具体细节的其他实施例中也可以实现本申请。在其他情况中，省略对众所周知的系统、装置、电路以及方法的详细说明，以免不必要的细节妨碍本申请的描述。

下面通过具体实施例来说明本申请的技术方案。

实施例一

参照图1，示出了本申请实施例提供的一种频域相关系数的估计方法的步骤流程示意图，具体可以包括如下步骤：

S101、从为终端分配的子载波中获取长度最长的一段目标子载波。

需要说明的是，本方法可以应用于无线通信领域。具体地，本方法的执行主体可以是上述无线通信领域中具备信道估计等计算功能的通信设备，本申请实施例对通信设备的具体类型不作限定。

在本申请实施例中，可以根据系统为终端分配的子载波中长度最长的一段目标子载波来进行频域相关系数的估计。

在通信领域，频域上连续的12个子载波，或者时域上一个时隙，被称为一个RB。若终端分配到的RB是连续的，则可以直接通过计算PDP来获取频域相关系数。

因此，在本申请实施例中，在进行频域相关系数估计前，可以首先确定为终端分配的子载波是否是连续的。

在具体实现中，可以通过定义如下条件来确定为终端分配的子载波是否是连续的：

S(k+1)-S(k)＝1，k＝0,1,2,…,K-1……(1)

其中，S为子载波的绝对索引集合，K为该子载波的个数。

若针对某一终端分配的子载波满足上述式(1)设定的条件，可以确定系统为该终端分配的子载波是连续的，该终端分配的数量为K个的子载波即为目标子载波；否则，则认为该终端分配到的子载波是不连续的。

若终端分配到的子载波是不连续的，则可以首先获取该终端分配到的子载波中长度最长的一段目标子载波。

由于一个RB包括频域上的12个连续的子载波，那么，在根据式(1)确定终端分配到的子载波不连续时，该不连续的子载波必然是由多段连续的子载波组成的，只是这些多段子载波相互之间是不连续的。

示例性的，若为终端分配的子载波的个数为K个，该K个子载波可以包括多段，如k₁、k₂、……k_n段。虽然作为一个整体，上述K个子载波是不连续的，但从其中每段的角度来看，k₁、k₂、……k_n都可以是连续的。

在本申请实施例中，可以从多段子载波中获取长度最长的一段目标子载波。例如，目标子载波为第k_m段子载波，该段子载波的个数可以为K′。

S102、判断所述目标子载波的个数是否小于设定的子载波阈值。

通常，在终端子载波连续时，可以直接通过计算PDP来估计频域相关信息。但是，一般认为，只有在该连续的子载波的个数足够多时，最终计算得到的频域相关信息才准确。

因此，在本申请实施例中，对于不连续的子载波中的某一段连续的目标子载波，可以首先判断该段目标子载波的个数K′与设定的子载波阈值K_Thr之间的大小关系。

如果上述K′小于K_Thr，则可以执行S104。

S103、通过信道估计得到信道估计值。

在本申请实施例中，对于子载波不连续，且其中最长的一段连续的子载波的个数也不够多(小于设定的子载波阈值)时，无法直接通过计算PDP来估计频域相关信息。

针对此种情况，本申请实施例可以首先针对信道进行信道估计。

在本申请实施例中，可以采用多种信道估计算法进行信道估计，如MMSE算法、最小二乘准则(least-square，LS)信道估计算法等等，本申请实施例对此不作限定。

作为本申请实施例的一种示例，可以采用LS信道估计算法，计算信道估计值。

S104、采用所述信道估计值计算频域相关值。

然后，可以直接根据信道估计值计算频域相关值。

示例性的，可以采用如下公式，计算频域相关值r′_f：

其中，H(k)为信道估计值。

S105、基于所述频域相关值和预设的多个信道模型，估计频域相关系数。

在根据信道估计值计算出频域相关值后，本申请实施例可以将频域相关值与多个不同的信道模型的频域相关函数进行比较，从而确定最终的频域相关系数。

通常，某个信道模型下的频域相关函数是固定的。例如，在一个呈指数衰落的多径功率时延信道模型中，频域相关函数只与子载波之间的距离有关，频域相关函数的表达式可以如下式(3)表示：

其中，τ_max为最大多径时延，Δf为系统子载波间隔。

在具体实现中，可以分别确定计算得到的频域相关值与预设的多个不同的信道模型的频域相关函数的相似度，然后将相似度最大值对应的信道模型确定为目标信道模型，根据目标信道模型，估计最终的频域相关系数。

实施例二

参照图2，示出了本申请实施例提供的另一种频域相关系数的估计方法的步骤流程示意图，具体可以包括如下步骤：

S201、确定为终端分配的子载波是否连续。

需要说明的是，前述实施例一是针对为终端分配的子载波不连续，且从不连续的子载波中无法获取到一段个数足够多的连续的子载波这一情形来介绍如何进行频域相关系数的估计的。本实施例二将针对为终端分配的子载波连续，即，上述实施例一中的“长度最长的一段目标子载波”即为整段(K个)为终端分配的子载波。但该连续的子载波的个数较少(小于设定的子载波阈值)这一情形来介绍如何进行频域相关系数的估计的。

S202、判断连续的所述子载波的个数是否小于所述子载波阈值。

在本申请实施例中，可以采用前述实施例中的公式(1)来确定为终端分配的子载波是否连续。若该子载波连续，则可以继续判断连续的子载波的个数与设定的子载波阈值之间的大小关系。

如果连续的子载波个个数小于上述子载波阈值，则可以执行S203-S205，采用与前述实施例相同的方式，估算频域相关系数。

S203、通过信道估计得到信道估计值。

S204、采用所述信道估计值计算频域相关值。

S205、基于所述频域相关值和预设的多个信道模型，估计频域相关系数。

由于本实施例二的上述步骤S203-S205与前述实施例一中步骤S103-S105类似，可以相互参阅，本实施例对此不再赘述。

实施例三

参照图3，示出了本申请实施例提供的又一种频域相关系数的估计方法的步骤流程示意图，具体可以包括如下步骤：

S301、确定为终端分配的子载波是否连续。

S302、判断所述子载波的个数是否大于或等于所述子载波阈值。

S303、根据所述子载波，计算信道的功率时延分布。

需要说明的是，本实施例三是针对为终端分配的子载波连续，即，上述实施例一中的“长度最长的一段目标子载波”为整段(K个)为终端分配的子载波，且该连续的子载波的个数足够多这一情形来介绍如何进行频域相关系数的估计的。

对于此类情形，可以通过计算PDP的方式来进行频域相关系数的估计。

在本申请实施例中，可以首先计算PDP的长度。示例性的，可以采用如下公式(4)计算PDP的长度：

其中，N_PDP表示该PDP的长度，K为连续的子载波的个数。

然后，可以对信道估计值H(k)补零，得到：

并对上式进行快速傅里叶逆变换(inverse fast fourier transform，IFFT)，获得信道冲击响应(channel impulse response，CIR)：

h(n)＝IFFT(H_N),n＝0,1,2...,N_PDP-1……(6)

通过依次对CIR进行求模和平方运算，可以获得该信道的PDP：

pdp(n)＝|h(n)|²,n＝0,1,2...,N_PDP-1……(7)

需要说明的是，在计算得到信道的PDP后，还可以对该PDP进行抑噪处理，本申请实施例对此不再赘述。

S304、根据所述功率时延分布，估计所述频域相关系数。

针对终端分配到的子载波连续个数足够多这一情形，在计算得到PDP后，可以直接采用PDP进行频域相关系数的估计。

在具体实现中，可以首先根据在应用正交频分复用技术(orthogonal frequencydivision multiplexing，OFDM)进行调制时的快速傅里叶变换(fast fourier transform，FFT)点数，对PDP补零。

示例性的，可以采用如下公式(8)对PDP进行补零：

其中，N为OFDM调制时的FFT点数。

然后，通过FFT处理，可以获得频域相关性系数，即：

r_f(k)＝fft(pdp_N),k＝0,1,...,N-1……(9)

实施例四

参照图4，示出了本申请实施例提供的再一种频域相关系数的估计方法的步骤流程示意图，具体可以包括如下步骤：

S401、确定为终端分配的子载波是否连续。

S402、从组成不连续的所述子载波的多段子载波中获取长度最长的一段目标子载波。

S403、判断所述目标子载波的个数是否小于设定的子载波阈值。

需要说明的是，本实施例四是针对为终端分配的子载波不连续，但可从不连续的子载波中获取到一段个数足够多的连续的子载波这一情形来介绍如何进行频域相关系数的估计的。

由于本实施例四中上述步骤S401-S403与前述实施例一中步骤S101-S102类似，可以相互参阅，本实施例对此不再赘述。

若目标子载波的个数K′大于或等于子载波阈值K_Thr，则可以执行S404。

S404、根据所述目标子载波，计算信道的功率时延分布。

在本申请实施例中，对于终端子载波不连续，但能够从不连续的子载波中获取到一段个数足够多(个数为K′)的连续子载波的情形，可以通过计算PDP的方式，来进行频域相关系数的估计。

在本申请实施例中，可以首先计算PDP的长度。示例性的，可以采用如下公式(10)计算PDP的长度：

其中，N_PDP表示该PDP的长度，K′为获取到的连续的子载波的个数。

然后，可以对信道估计值H(k)补零，得到：

在完成上述步骤后，可以通过执行与前述实施例三中公式(5)、(6)、(7)类似的步骤，得到信道的PDP。

S405、根据所述功率时延分布，估计所述频域相关系数。

由于本实施例四中的上述步骤S405与前述实施例中步骤S304类似，可以相互参阅，本实施例对此不再赘述。

需要说明的是，上述实施例中各步骤的序号的大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本申请实施例的实施过程构成任何限定。

为了便于理解，下面结合具体的示例，对上述多种情形下的频域相关信息的估计过程作一介绍。

如图5所示，是本申请实施例提供的一种频域相关系数的估计方法的示意图。按照图5所示，在开始进行处理时，可以首先从为终端分配的子载波中获取长度最长的一段目标子载波。下面针对不同的情形，分别进行介绍。

(1).K个子载波连续。

整段(K个)为终端分配的子载波即为“长度最长的一段目标子载波”，对于K个子载波连续的情形，又可以根据子载波的个数K与设定的子载波阈值K_Thr之间的大小关系的不同，分为两种情形。即如下介绍的(1.1)和(1.2)。

(1.1).K≥K_Thr

对于K≥K_Thr，如图5所示，可以通过IFFT变换获得信道的CIR，在此基础上，可以计算信道的PDP，并通过抑噪处理后，进行FFT变换，获得频域相关系数。

(1.2).K＜K_Thr

对于K＜K_Thr，如图5所示，可以直接根据信道估计值计算频域相关值，并将其与多个信道模型进行比较，根据确定出的最相似的目标信道模型来获取频域相关系数。

(2).K个子载波不连续。

对于K个子载波不连续的情形，可以首先从不连续的子载波中获取到连续的且长度最长的一段子载波，该段子载波的个数为K′。根据上述个数K′与设定的子载波阈值K_Thr之间的大小关系的不同，分为两种情形。即如下介绍的(2.1)和(2.2)。

(2.1).K′≥K_Thr

对于K′≥K_Thr，如图5所示，可以通过IFFT变换获得信道的CIR，在此基础上，可以计算信道的PDP，并通过抑噪处理后，进行FFT变换，获得频域相关系数。

由图5可见，该情形的处理步骤与情形(1.1)的处理步骤相同。

(2.2).K′＜K_Thr

对于K′＜K_Thr，如图5所示，可以直接根据信道估计值计算频域相关值，并将其与多个信道模型进行比较，根据确定出的最相似的目标信道模型来获取频域相关系数。

由图5可见，该情形的处理步骤与情形(1.2)的处理步骤相同。

参照图6，示出了本申请实施例提供的一种频域相关系数的估计装置的示意图，具体可以包括目标子载波获取模块601、子载波个数判断模块602、信道估计模块603、频域相关值计算模块604、频域相关信息估计模块605，其中：

目标子载波获取模块601，用于从为终端分配的子载波中获取长度最长的一段目标子载波；

子载波个数判断模块602，用于判断所述目标子载波的个数是否小于设定的子载波阈值；

信道估计模块603，用于若所述目标子载波的个数小于所述子载波阈值，则通过信道估计得到信道估计值；

频域相关值计算模块604，用于采用所述信道估计值计算频域相关值；

频域相关信息估计模块605，用于基于所述频域相关值和预设的多个信道模型，估计频域相关系数。

在本申请实施例的一种可能的实现方式中，频域相关系数估计模块605具体用于：分别确定所述频域相关值与预设的多个信道模型的频域相关函数的相似度；将所述相似度最大值对应的信道模型确定为目标信道模型；根据所述目标信道模型，估计所述频域相关系数。

在本申请实施例的另一种可能的实现方式中，信道估计模块603具体用于：采用最小二乘准则信道估计算法，计算所述信道估计值。

在本申请实施例的另一种可能的实现方式中，频域相关函数可以表示为：

其中，τ_max为最大多径时延，Δf为系统子载波间隔。

在本申请实施例的另一种可能的实现方式中，子载波个数判断模块602具体用于若所述目标子载波的个数大于或等于所述子载波阈值，则根据所述目标子载波，计算信道的功率时延分布；根据所述功率时延分布，估计所述频域相关系数。

在本申请实施例的一种可能的实现方式中，频域相关信息估计模块605具体用于：根据所述功率时延分布，通过抑噪处理后，进行快速傅里叶变换，获得所述频域相关系数。

在本申请实施例的另一种可能的实现方式中，所述从为终端分配的子载波中获取长度最长的一段目标子载波，是基于如下公式获取：

S(k+1)-S(k)＝1，k＝0,1,2,…,K-1……(13)

其中，S为所述子载波的绝对索引集合，K为所述子载波的个数。

对于装置实施例而言，由于其与方法实施例基本相似，所以描述得比较简单，相关之处参见方法实施例部分的说明即可。

参照图7，示出了本申请一个实施例的一种通信设备的示意图。如图7所示，本实施例的通信设备700包括：处理器710、存储器720以及存储在所述存储器720中并可在所述处理器710上运行的计算机程序721。所述处理器710执行所述计算机程序721时实现上述频域相关信息的估计方法各个实施例中的步骤，例如图1所示的步骤S101至S106。或者，所述处理器710执行所述计算机程序721时实现上述各装置实施例中各模块/单元的功能，例如图6所示模块601至606的功能。

示例性的，所述计算机程序721可以被分割成一个或多个模块/单元，所述一个或者多个模块/单元被存储在所述存储器720中，并由所述处理器710执行，以完成本申请。所述一个或多个模块/单元可以是能够完成特定功能的一系列计算机程序指令段，该指令段可以用于描述所述计算机程序721在所述通信设备700中的执行过程。例如，所述计算机程序721可以被分割成目标子载波获取模块、子载波个数判断模块、信道估计模块、频域相关值计算模块、频域相关系数估计模块，各模块具体功能如下：

目标子载波获取模块，用于从为终端分配的子载波中获取长度最长的一段目标子载波；其中，任意一段子载波为连续的子载波；

子载波个数判断模块，用于判断所述目标子载波的个数是否小于设定的子载波阈值；

信道估计模块，用于若所述目标子载波的个数小于所述子载波阈值，则通过信道估计得到信道估计值；

频域相关值计算模块，用于采用所述信道估计值计算频域相关值；

频域相关系数估计模块，用于基于所述频域相关值和预设的多个信道模型，估计频域相关系数。

所述通信设备700可以是基站或终端，其中终端包括手机、平板或个人计算机。所述通信设备700可包括，但不仅限于，处理器710、存储器720。本领域技术人员可以理解，图7仅仅是通信设备700的一种示例，并不构成对通信设备700的限定，可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件，例如所述通信设备700还可以包括输入输出设备、网络接入设备、总线等。

所述处理器710可以是中央处理单元(Central Processing Unit，CPU)，还可以是其他通用处理器、数字信号处理器(Digital Signal Processor，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现成可编程门阵列(Field-Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。

所述存储器720可以是所述通信设备700的内部存储单元，例如通信设备700的硬盘或内存。所述存储器720也可以是所述通信设备700的外部存储设备，例如所述通信设备700上配备的插接式硬盘，智能存储卡(Smart Media Card，SMC)，安全数字(SecureDigital，SD)卡，闪存卡(Flash Card)等等。进一步地，所述存储器720还可以既包括所述通信设备700的内部存储单元也包括外部存储设备。所述存储器720用于存储所述计算机程序721以及所述通信设备700所需的其他程序和数据。所述存储器720还可以用于暂时地存储已经输出或者将要输出的数据。

本申请实施例还公开了一种通信设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如前述各个实施例所述的频域相关系数的估计方法。

本申请实施例还公开了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现如前述各个实施例所述的频域相关系数的估计方法。

本申请实施例还公开了一种计算机程序产品，当所述计算机程序产品在通信设备上运行时，使得所述通信设备执行前述各个实施例所述的频域相关系数的估计方法。

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