1.一种用于NB‑IoT网络中已经连接到服务小区的用户设备UE的联合时频同步方法，以获得最优时频偏移估计，所述方法包括以下阶段：‑阶段1：相对于包含频域中的期望资源元素的预先计算的参考NB‑IoT子帧Bm，解调接收的包含所述频域中的资源元素的NB‑IoT子帧Rm，Bm的预先计算的资源元素是从所述服务小区的信息中导出的，并且将解调的NB‑IoT子帧Wm的资源元素存储在矩阵W中，其中，m是等于或大于0的自然数，并且呈现第m个处理的子帧；

‑阶段2：计算由 定义的度量Vm

(f，t)，其中，Vm(f，t)表示第m个处理的子帧的度量，Wk,l是第m个处理的子帧Wm的矩阵W的元素，k表示所述矩阵W的行数，l表示所述矩阵W的列数；并计算关于预定义的离散频率偏移值f∈F和离散时间偏移值t∈T的区域的函数 和pl(f)＝exp(2πjτlf)；其中F是频率偏移区域，T是时间偏移区域，NC是子帧中的子载波数，NS是子帧中的符号数，Δfs是子载波间距，τl是子帧内第l个OFDM符号的开始的时间偏移；

‑阶段3：在多个解调的NB‑IoT子帧Wm上以非相干方式递归组合所述度量Um(f，t)，产生新的度量Um(f，t)；

‑阶段4：基于当前处理的解调的NB‑IoT子帧Wm的子帧能量，递归计算多个解调的NB‑IoT子帧Wm上的参考值Sm；

‑阶段5：将非相干组合的新度量Um(f，t)与所述参考值Sm进行比较，

‑如果Um(f，t)/Sm的归一化度量值不超过阈值KThreshold(m)，则递归重复所述方法的阶段

1至5，直到解调了预定义的子帧最大数量Mmax，其中，Mmax是大于0的自然数，否则‑阶段6：如果所述Um(f，t)/Sm的归一化度量值超过所述阈值KTreshold(m)，则基于Um(f，t)的非相干组合度量值来计算关于所述离散频率偏移值f∈F和所述离散时间偏移值t∈T的arg‑max值，获得所述时频偏移估计的最优值，假设组合的所述子帧Wm的数量没有超过Mmax。

2.根据权利要求1所述的用于NB‑IoT网络中已经连接到服务小区的UE的联合时频同步方法，其中，对于所有Ns＝14个OFDM符号的所述离散频率偏移值f∈F的所有复指数被预先计算并存储在查找表中。

3.根据权利要求1所述的用于NB‑IoT网络中已经连接到服务小区的UE的联合时频同步方法，其中，对于所有Nc＝12个OFDM符号的所述离散时间偏移值t∈T的所有复指数被预先计算并存储在查找表中。

4.根据权利要求1所述的用于NB‑IoT网络中已经连接到服务小区的UE的联合时频同步方法，其中，所述度量Vm(f，t)的计算以两阶段方法进行，其中，在第一阶段，首先基于所述离散频率偏移值f∈F计算中间值ck(f)，在第二阶段，基于所述第一阶段的中间值ck(f)以及基于离散值t∈T和f∈F计算所述度量Vm(f，t)。

5.根据权利要求1所述的用于NB‑IoT网络中已经连接到服务小区的UE的联合时频同步方法，其中，在锚载波的子帧#5上传输的窄带主同步信号在阶段1中被解调，而所述预先计算的参考NB‑IoT子帧Bm是恒定的，与LTE帧号无关。

6.根据权利要求1所述的用于NB‑IoT网络中已经连接到服务小区的UE的联合时频同步方法，其中，在对于每个偶数LTE帧号的锚载波的子帧9上传输的窄带辅同步信号在阶段1中被解调，而所述参考NB‑IoT子帧Bm的预先计算取决于小区ID。

7.根据权利要求1所述的用于NB‑IoT网络中已经连接到服务小区的UE的联合时频同步方法，其中，在窄带物理广播信道NPBCH中传递的主信息比特MIB有效载荷信息在阶段1中被解调，而所述参考NB‑IoT子帧Bm的预先计算在同步过程开始时执行，并且取决于所述NPBCH的资源元素的重新编码。

8.根据权利要求1所述的用于NB‑IoT网络中已经连接到服务小区的UE的联合时频同步方法，其中，窄带唤醒信号NWUS在阶段1中被解调，而所述参考NB‑IoT子帧Bm的预先计算取决于NWUS所关联的小区ID和寻呼时机的时间。