专利名称：搜索最佳载波频率偏移校正值的方法
技术领域：
本发明属于移动通信系统领域，更具体地，涉及在移动通信系统的上行链路的接 收机处对载波频率偏移(CFO，Carrier Frequency Offset)进行补偿时，针对每一个单用户 搜索其最佳CFO校正值的方法。
背景技术：
^t JE M ^ ^ it & \ (OFDMA, Orthogonal Frequency Division MultipleAccess)系统的上行链路中，信号在传输时由于多普勒效应等因素会造成载 波频率的偏移，接收端处不同用户的不同CFO会导致多用户干扰(MUI，Multi-User Interference)。为了抑制MUI，需要对上行链路的CFO做出估计，进而在基站端使用某种方 式，利用估计出的上行链路CFO来恢复子载波之间的正交性。由于每一个单用户的发射信 号受到不同的CFO的影响，而又必须通过CFO校正来抑制MUI，所以该CFO校正过程是一个 重要的，且具有挑战性的任务。针对上述问题，存在一种在OFDMA系统上行链路的接收机处 对CFO进行补偿的方法。该方法包括计算针对每一个单用户的最佳CFO校正值，所述最佳 CFO校正值是使针对所述每一个单用户的平均信号干扰比(SIR，Signal-to-Interference Ratio)最大的CFO校正值；利用所获得的每一个单用户的最佳CFO校正值，对每一个单用 户的CFO进行补偿。上述方法的一个重要步骤是计算针对每一个单用户的最佳CFO校正值。所述最 佳CFO校正值是使所述每一个单用户的平均SIR最大的CFO校正值。所述最佳CFO校正值 位于所有用户CFO的最小值min{￡l，ε2，…εΜ}(其中ε m表示某一个单用户m的归一化 CFO值，M表示所有用户的数目)到最大值max { ε ” ε 2，…εΜ}之间。针对某一个单用户
m，传统搜索其最佳CFO校正值的方法是按以下过程进行的，具体流程如图2所示从所
有用户CFO的最小值Hiinh1, ε2,…ε Μ}开始，以某一固定步长L递增，直到增加到所有 用户 CFO 的最大值 max { ε 工，ε 2，... ε J 为止，以 min { ε ” ε 2，··· ε M}, min{ ε ε 2，... eM}+L、min{ ε ” ε 2,…eM}+2L、->max{ ε ε 2，…ε J为该用户平均信号干扰比函数 SIRw的自变量，依次计算上述自变量对应SIRw的值，比较所得函数值的大小，则其中的最 大值所对应的自变量为要搜索的最佳CFO校正值^w)= argmaxSIR(w)。但是，当
所有用户CFO的最大值max { ε ε 2，…ε Μ}与最小值Hiinh1, ε 2，…ε J差值较大时， 传统搜索方法的复杂度较高。
发明内容
本发明的目的是解决搜索最佳CFO校正值时传统方法复杂度较高的问题，提出一 种在移动通信系统的上行链路的接收机处搜索最佳载波频率偏移校正值的方法。为了实现上述目的，本发明的技术方案是一种在移动通信系统的上行链路的接 收机处搜索最佳载波频率偏移校正值的方法，其特征在于，包括如下步骤
3
步骤1 针对某一个单用户m，确定搜索其最佳CFO校正值的区间[a，b]，其中 a、b分别表示所有用户CFO的最小值、最大值，即a = min { ε ε 2，…ε Μ}，b = max { ε ε2, ."εΜ}，其中，￡m，m=l，2，.··，M，表示某一个单用户m的归一化CFO值，M表示所有 用户的数目；步骤2:令、=￡1+(1-3)03-￡1)，、=￡1+3 03-￡1)，计算、、、的值，其中，β = 0. 618 ；步骤3 判断、与、的差值是否达到精度要求。设要求的精度为1，若It^t2I <1， 搜索停止，此时的值即为要搜索的针对某一个单用户m的最佳CFO校正值^w)否则，
2bo ’
若Itft2I彡1，转步骤4;步骤4 以、、t2为自变量，分别计算对应的平均信号干扰比函数SIRw，记为SIRw (、)、SIRw (t2)，判断 SIRw (、) ^ SIR(m) (t2)是否成立，若成立，则置 a = t1; b = b, ^ = t2，t2 = a+β (b-a)，然后转步骤 3 ；否则，若 SIRw U1) > SIR(m) (t2)，则置 a = a, b = t2, t2
= ^a1 = a+(l-β ) (b-a)，然后转步骤3，当I t「t21 < 1时，搜索停止，此时^^的值即为 要搜索的针对某一个单用户m的最佳CFO校正值g㈨；步骤5 重复步骤1至步骤4，直到搜索出所有用户的最佳CFO校正值41)、 €)、…、I)为止。优选地，所述上行链路为正交频分多址接入系统的上行链路。本发明的有益效果本发明通过步骤4，利用迭代的方式逐步缩小搜索范围，搜索 移动通信系统的上行链路的接收机处的最佳载波频率偏移校正值。该方法与传统搜索方法 相比，大大降低了算法复杂度。
图1是对CFO进行补偿的OFDMA系统的上行链路工作原理图。图2是用传统方法搜索最佳CFO校正值的流程图。图3是用本发明方法搜索最佳CFO校正值的流程图。附图标记说明子载波映射1，N点离散傅里叶逆变换(IDFT) 2，添加保护间隔模块 3，D/A转换4，上变频5，下变频6，A/D转换7，搜索最佳CFO校正值8，频率偏移补偿9，去除 保护间隔10，N点离散傅里叶变换(DFT) 11，子载波逆映射12，其中N表示每个OFDM符号中 存在的子载波个数。
具体实施例方式下面将结合附图，给出本发明的具体实施例。需要说明的是实施例中的参数并不 影响本发明的一般性。为了便于对具体实施例理解，先对CFO进行补偿的OFDMA系统的上行链路的发射 机部分和接收机部分的工作原理进行说明，具体工作原理如图1所示。发射机部分包括子载波映射1、IDFT 2、添加保护间隔3、D/A转换4、上变频5，具体工作过程为考虑具有M个用户的OFDMA系统，其中每一个单用户通过独立的多径信道 来与基站进行通信。假定在每一个OFDM符号中存在N个子载波，将这N个子载波平均分配 给M个用户，这样，每一个单用户具有P = N/M个子载波。在第k个子载波处针对第m个
单用户的信息符号表示为;^w)，ke Γω，其中Γm是分配给用户m的子载波集合。于是
M
IJr^ = ·^,!,···,"-1丨，并且对于i乒j，_ 针对交织的子载波分配方案，将
m=\丄；_| Ily =U。
rm 定义为rm = {m-l+rM|r = 0，…，卩-丄}。在OFDMA系统中，保护间隔的长度等于Ng个采样，并且假定其长于最大信道延迟 扩展。在发射机处经过离散傅里叶逆变换和保护间隔插入之后，第m个单用户的时域序列
jcf)由下式给出
j2nnkJm) _ γ γΠ N -Nr ^ η ^ N_l. Χη -乙 kermAk e , g接收机部分包括下变频6、A/D转换7、搜索最佳CFO校正值8、频率偏移补偿9、去 除保护间隔10、DFT11、子载波逆映射12，具体工作过程为在通过衰落信道、下变频和A/D 转换之后，第m个单用户的信号由以下等式给出其中，“*”表示线性卷积，^w)是第m个单用户的信道冲激响应。假定仅对于
n = 0,…，L-I是非零的，其中L是最大信道延迟扩展。通过考虑CFO和加性噪声，接收到 的基带信号由以下等式表示
M 2πεΜη^ =N,-Ng 彡 η 彡 Ν-1·
m=\其中，εω，πι=1，…，M表示归一化的第m个单用户的CF0，并且ζη是加性白高斯 噪声。对接收到的基带信号rn，首先搜索每一个单用户的最佳载波频率偏移CFO校正值，该 CFO校正值是使所述每一个单用户的平均SIR最大的CFO校正值，它位于所述所有用户CFO 的最小值min{ ε ε 2，…ε J到最大值max { ε ε 2，…ε J之间。为了校正CF0，在进行离
散傅里叶变换之前，针对某一个单用户m，将所接收到的信号rn乘以时域序列
e Ν ,
2πε{,η)η
得到。0 ^ n^ N-I0对。进行离散傅里叶变换之后，得到频域信号Y，子
yn _ rne‘^n
载波逆映射后在Y中选择针对用户m的子载波，作为单用户m的输出信号。在对CFO进行补偿的OFDMA系统的上行链路中，根据上述发射机部分和接收机部 分的工作过程，可以推导出在第k个子载波处针对第m个单用户，接收到的所需信号功率与 干扰信号功率比值的表达式，即在第k个子载波处针对第m个单用户平均信号干扰比函数
sir4w)(O的表达式为
5'm . 2π{ε- ) sin v m J SIR^CO = ^--T^
y y sin2^ei-t + gj-k)
'=I q^T^q^k sin ---
N其中，t 代表函数的自变量，t e [min{ ε 丨，ε 2，- εΜ}, max{ ε ε 2,…ε Μ}]
2
^eTi.,]，//，表示对;进行离散傅里叶变换后所得向量的第屮个
元素，Ε[.]表示求平均。特别地，如果子载波分配方案采用交织分配，则VA，ferw, SIR(Zw)⑴=SIRf⑴=SIR⑷⑴。其中，SIRw⑴表示某一个单用户m的平均信号干扰 比函数。本实施例的仿真参数为用户数M = 4，子载波数N = 64的对CFO进行补偿的OFDMA 系统的上行链路，子载波分配方案为交织分配。所有用户的CFO值[ε ” ε2, ε 4，ε J =
。采用传统的搜索方法时，步长L = 0. 001，采用本发明方法时，要 求的精度1 =0. 001。在计算用户平均信号干扰比的值时，取一个属于用户m的子载波k，根据上述第k 个子载波处针对第m个单用户平均信号干扰比函数SIR^)的表达式，依次计算其对应
sir4w)(/)的值，在本实施例中，由于子载波分配方案采用的是交织分配，所以的 值等于该用户平均信号干扰比的值，记为SlRw (t)。下面结合具体仿真参数，对传统搜索方法步骤与本发明方法步骤进行说明。在本实施例中，使用传统搜索方法搜索某一个单用户m最佳CFO校正值的步骤为 从所有用户CFO的最小值min{ ε ε 2，ε 3，ε 4} = 0. 08开始，以固定步长L = 0. 001递增， 直到增加到所有用户CFO的最大值max { ε ε 2, ε 3, ε J = 0. 16为止，分别令t = 0. 08，t =0. 08+0. 001，…，t = 0. 16，依次计算其对应SIR(m) (t)的值，比较所得函数值的大小，则
其中的最大值所对应的自变量为要搜索的最佳CFO校正值<w)。分别令m= l，m = 2，…， m = 4，重复上述步骤，直到求出所有用户的最佳CFO校正值^1)，·^2)，…,‘4)为止。本发明方法流程图如图3所示。在本实施例中，使用本发明方法搜索所有用户最 佳CFO校正值的步骤为步骤1 针对某一个单用户m，确定搜索其最佳CFO校正值的区间 [a, b], a =
min{ ε ” ε 2，ε 3，ε J = 0. 08, b = max{ ε 工，ε 2，ε 3，ε 4} = 0. 16 ；步骤2 令、=ει+(1-β) (b-a)，t2 = a+β (b-a)，计算、、t2 的值。其中，β = 0. 618 ；步骤3 令要求的精度1 =0.001。若|trt2| < 1，搜索停止，此时^^的值即为
要搜索的针对某一个单用户m的最佳CFO校正值否则，若It^t2I彡1，转步骤4;步骤4 以、、、为自变量，计算SIRw U1^SIRw (t2)的值，判断SIRw (、)彡SIRw
6(t2)是否成立，若成立，则置a = t1 b = b，、= t2，t2 = a+β (b-a)，然后转步骤3 ；否则， ^ SIR(m) (t^ > SIRw(t2)，则置 a = a, b = t2，t2 = t1 、= a+(1_ β ) (b_a)，然后转步骤
3。当|trt2| <1时，搜索停止，此时的值即为要搜索的针对某一个单用户m的最佳 CFO校正值^w)；步骤5 分别重复步骤1到步骤4，直到求出所有用户的最佳CFO校正值 ^),·.·，44)为止。下面分析传统搜索法以及本发明的方法复杂度，以验证本发明的性能。针对某一个单用户m，计算其平均信号干扰比函数SIRw (t)时，实数乘法次数为
b — a
N，实数除法次数为N+1。对于传统搜索方法，每一个单用户要计算次平均信号干扰比
b — Cib — Ci
函数，则总的实数乘法次数为# χ——X M，总的实数除法次数为(iV +1) X —— χ M。
LL
对于本发明的方法，令k为本发明方法的迭代次数，针对每一个单用户，第一次迭代时，需
计算、和t2以及两次平均信号干扰比函数，以后每次迭代需计算、或者t2，以及一次平
均信号干扰比函数。这样，总的实数乘法次数为(k+l)(N+l)M，总的实数除法次数为(k+1)
(N+1) M0根据以上分析，在本实施例中，代入具体的仿真参数，使用传统搜索
法，总的实数乘法次数为64x °·^"0·08 χ 4 = 20480，总的实数除法次数为
0.001
(64 +1)χ aiQQQ0108 x4 = 20800 ；使用本发明的方法，所需迭代次数为7次，则总的
实数乘法次数为(7+1) X (64+1) X4 = 2080，总的实数除法次数为(7+1) X (64+1) X4 = 2080。可见，与传统搜索方法相比，本发明的方法能够大大降低算法复杂度。以上实例仅为本发明的优选例子而已，本发明的使用并不局限于该实例，凡在本 发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范 围之内。
权利要求
一种在移动通信系统的上行链路的接收机处搜索最佳载波频率偏移校正值的方法，其特征在于，包括如下步骤步骤1针对某一个单用户m，确定搜索其最佳CFO校正值的区间[a，b]，其中a、b分别表示所有用户CFO的最小值、最大值，即a＝min{ε1，ε2，…εM}，b＝max{ε1，ε2，…εM}，其中，εm，m＝1，2，…，M，表示某一个单用户m的归一化CFO值，M表示所有用户的数目；步骤2令t1＝a+(1 β)(b a)，t2＝a+β(b a)，计算t1、t2的值，其中，β＝0.618；步骤3判断t1与t2的差值是否达到精度要求。设要求的精度为l，若|t1 t2|＜l，搜索停止，此时的值即为要搜索的针对某一个单用户m的最佳CFO校正值否则，若|t1 t2|≥l，转步骤4；步骤4以t1、t2为自变量，分别计算对应的平均信号干扰比函数SIR(m)，记为SIR(m)t1)、SIR(m)(t2)，判断SIR(m)(t1)≤SIR(m)(t2)是否成立，若成立，则置a＝t1，b＝b，t1＝t2，t2＝a+β(b a)，然后转步骤3；否则，若SIR(m)(t1)＞SIR(m)(t2)，则置a＝a，b＝t2，t2＝t1，t1＝a+(1 β)(b a)，然后转步骤3，当|t1 t2|＜l时，搜索停止，此时的值即为要搜索的针对某一个单用户m的最佳CFO校正值步骤5重复步骤1至步骤4，直到搜索出所有用户的最佳CFO校正值为止。FSA00000286777500011.tif,FSA00000286777500012.tif,FSA00000286777500013.tif,FSA00000286777500014.tif,FSA00000286777500015.tif,FSA00000286777500016.tif,FSA00000286777500017.tif
2.根据权利要求1所述的方法，其特征在于，所述的上行链路为正交频分多址接入系 统的上行链路。
全文摘要
本发明公开了一种搜索最佳载波频率偏移校正值的方法，涉及一种在移动通信系统的上行链路的接收机处对载波频率偏移进行补偿时，针对每一个单用户搜索其最佳载波频率偏移校正值的方法。针对在移动通信系统的上行链路的接收机处搜索最佳载波频率偏移校正值时，传统搜索方法复杂度较高的问题，本发明通过利用迭代的方式逐步缩小搜索范围，搜索最佳载波频率偏移校正值。该方法与传统搜索方法相比，大大降低了算法复杂度。
文档编号H04L27/26GK101958872SQ201010294508
公开日2011年1月26日 申请日期2010年9月28日 优先权日2010年9月28日
发明者李少谦, 杨佳, 王刚, 肖悦, 雷霞 申请人:电子科技大学