一种基于宽带ofdm系统频变多普勒频移的估计方法
【技术领域】
[0001] 本发明涉及通信技术领域，更具体地，涉及一种基于宽带OFDM系统频变多普勒频 移的估计方法。
【背景技术】
[0002] 随着移动通信技术和互联网技术的迅猛发展，产生了很多种宽带无线通信技术， 而这些宽带无线通信技术的降临使得用户对宽带业务的高需求成为可能。OFDM技术作为宽 带无线通信技术中的一种核心技术，不仅能够有效地对抗多径干扰，而且具有频谱利用率 高和数据传输能力强的优点，OFDM技术在实际应用中应用极为广泛。
[0003] OFDM系统是一种子载波间保持相互正交的多载波通信系统，系统对子载波间的正 交性要求很高，当系统中存在频率偏差时，子载波间正交性将会遭到破坏，导致子载波间干 扰ICI的产生，从而降低系统的性能。在无线通信系统中，当收发两端存在相对运动时，系 统会产生多普勒效应，接收端的接收信号频率与发送信号之间会产生频率偏移现象。因此， 多普勒频移会对系统性能造成严重的影响。
[0004] 在宽带OFDM系统中，系统的多普勒频移并不像窄带OFDM系统那样可以认为在频 域上是不变的，相反地，宽带OFDM系统中的多普勒频移随着频率的变化而变化。因此，在宽 带OFDM系统中研究多普勒频移的相关问题已经变得十分迫切了。
【发明内容】
[0005] 本发明的目的在于克服现有窄带OFDM多普勒频移估计方法存在的不足，提出了 一种基于宽带OFDM系统频变多普勒频移的估计方法，本发明不仅能够有效地估计得到宽 带OFDM系统所有子载波上的多普勒频移，而且能够获得较好的估计性能。
[0006] 为了实现上述目的，本发明的技术方案如下：
[0007] -种基于宽带OFDM系统频变多普勒频移的估计方法，具体步骤如下：
[0008] a)设计包含两块导频簇的导频分布；
[0009] b)估计多普勒频移向量W (0, 0);
[0010] C)通过IFFT计算多普勒频移的复指数函数卬）：
[0011] d)估计出导频簇位置的多普勒频移；
[0012] e)利用得到的多普勒频移，通过线性拟合得到所有子载波上的多普勒频移。
[0013] 进一步地，所述步骤a)中，根据最小均方误差MMSE准则，导频序列可以选取为：
[0014]
[0015]
[0016] 其中，kp_T，…，kp+Q+T是导频簇的导频序号，固定值E为多普勒频移估计的功率， r (k) e [0, 1]为一个均匀分布的随机变量，k表示导频序号，T、Q均为正整数；导频序列表 示为 S (pP+kp-T)，…，S (pP+kp+Q+T)。
[0017] 此外，导频序号可以表不为kpl_T,…，kpl+Q+T和k p2_T,…，kp2+Q+T的连续导频组 成两块导频簇。其中，Q为一个正整数。也就是说，导频符号的数量为2*(2T+1+Q)。由 于两块导频簇的分布情况相同，且估计方法也相同，为了研究方便，在接下来的估计算法 的描述中，只描述其中一块导频簇的估计算法，另外，统一假设这块导频簇的导频序号为 kp-T，…，kp+Q+T，其对应的多普勒频移则用fd表示。导频序列可以表示为s(pP+k p-T)，… ,s (pP+kp+Q+T) 〇
[0018] 进一步地，所述步骤b)中，多普勒频移向量W(0,0)的表达式如下：
[0019]
[0020] 其中，ft表示第t个子载波上的多普勒频移，是某个不确定的子载波上的多普勒频 移，t e (1，2, ···，《，P表示子载波个数，w (i)表示第i个子载波上的多普勒频移的复指数 函数的 FFT 变换，i = 1，2, ..·，Ρ。
[0021] 假设多普勒频移的大小小于半个子载波的间隔，即满足：
[0022]
[0023] 多普勒频移向量W(0, 0)的主要部分或者主瓣在频域上会集中在两端的位置，即 集中在子载波序号〇,…，T与P-T，…，P-I的两段区间。其中，T为一正整数。那么，多普勒 频移向量W(0,0)可以近似为：
[0024]
[0025] 所述步骤c)中，在接收端，把导频序号为kp，…，kp+Q的子载波通过FFT解调输出 信号组成一个向量，可以表TK为：
[0026]
[0027] 上式中，V(P)为噪声向量，w(p) = [w(P-T),…，w(P-l), w(0),…，wW]T是由多 普勒频移向量W(0,0)中的非零元素组成。E为一个（Q+1)X(2T+1)的矩阵，矩阵的表达式 为：
[0028]
[0029] 其中，当满足Q+1多2T+1，且矩阵E满足列满秩时，多普勒频移参数向量W(p)可以 被估计为：
[0030]
[0031] 其中，1*??)为w(p)的估计，E+矩阵为E矩阵的伪逆。
[0032] 然后通过IFFT运算可以得到复指数函数f y27rU，表达式如下所示：
[0033]
[0034] 其中，fd表示导频簇位置对应的多普勒频移，是一个具体的子载波上的多普勒频 移。
[0035] 所述步骤d)中，€_在频域两端的误差会比频域的中间部分更大。所以，应该选 择贫?)的中间部分的样本估计多普勒频移。估计公式为：
[0036]
[0037] 其中，λ为一正整数，且λ〈Ρ/2,符号'Ζ'表示取相位角操作。
[0038] 所述步骤e)中，假设采用这种方法估计得到的各组多普勒频移数据为：
[0039] t = {tl, t2, ···, tn}, fd = {fdl, fd2, ···, fdn}
[0040] 其中，t表示子载波序号，fd表示多普勒频移的估计值，n表示数据量的大小。
[0041] 接着，采用最小二乘法对得到的数据做线性拟合，拟合后的表达式可以表示如 下：
[0042] fd = at+b
[0043] 其中，a与b分别为：
[0044]
[0045]
[0046] 这样，各个子载波上的多普勒频移都估计出来了。
[0047] 此外，还可以通过加窗的方法来改善W(0, 0)的估计性能，得到多普勒频移估计的 改进算法，具体过程如下：
[0048] 接收信号经过加窗处理后可得：
[0049]
[0050] 其中，f(p，t)为一窗函数，可选择为Kaiser窗或Hamming窗。可得：
[0051]
[0052]
[0053] 印）已经变为,印）的FFT输出向量则变为：
[0054]
[0055]
[0056]
[0057] 与现有的技术相比，本发明能够实现宽带OFDM系统频变多普勒频移的估计，并获 得良好的估计性能。
【附图说明】
[0058] 图1是本发明实现的流程图。
[0059] 图2是导频符号在OFDM符号中的分布图。
[0060]图3是第一块导频簇的复指数函数^)估计示意图。
[0061]图4是第二块导频簇的复指数函数cX/)估计示意图。
[0062] 图5是拟合后的多普勒频移估计相对误差曲线。
【具体实施方式】
[0063] 下面结合附图对本发明做进一步的描述，但本发明的实施方式并不限于此。
[0064] 本发明的思想为：利用两块导频簇估计得到两块导频簇位置的多普勒频移，再根 据估计得到的多普勒频移采用线性拟合的方法估计得到所有子载波上的多普勒频移，整体 的流程图如图1所示。
[0065] 本发明的具体实现过程如下：
[0066] 首先，在发送端设计如图2所示的包含两块导频簇的导频符号分布。根据最小均 方误差MMSE准则，可以设计导频序列如下式：
[0067]
[0068]
[0069] 其中，固定值E为多普勒频移估计的功率，r(k) e [0, 1]为一个均匀分布的随机 变量，T为一正整数，导频簇的长度都应不小于2T+1，kp表示导频簇位置的子