CN101247375A - 一种进行载波频率偏移估计的方法及装置 - Google Patents
一种进行载波频率偏移估计的方法及装置 Download PDFInfo
- Publication number
- CN101247375A CN101247375A CNA2007100792496A CN200710079249A CN101247375A CN 101247375 A CN101247375 A CN 101247375A CN A2007100792496 A CNA2007100792496 A CN A2007100792496A CN 200710079249 A CN200710079249 A CN 200710079249A CN 101247375 A CN101247375 A CN 101247375A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msup
- carrier frequency
- frequency offset
- mover
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000013508 migration Methods 0.000 title 1
- 230000005012 migration Effects 0.000 title 1
- 230000004044 response Effects 0.000 claims abstract description 21
- 238000007476 Maximum Likelihood Methods 0.000 claims abstract description 14
- 239000000969 carrier Substances 0.000 claims description 17
- 230000008569 process Effects 0.000 claims description 16
- 125000004122 cyclic group Chemical group 0.000 claims description 13
- 230000008030 elimination Effects 0.000 claims description 8
- 238000003379 elimination reaction Methods 0.000 claims description 8
- 238000012804 iterative process Methods 0.000 claims description 6
- 238000006467 substitution reaction Methods 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 5
- 101100272279 Beauveria bassiana Beas gene Proteins 0.000 claims 2
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 claims 1
- 230000005540 biological transmission Effects 0.000 abstract description 12
- 230000001788 irregular Effects 0.000 abstract description 6
- 230000035939 shock Effects 0.000 abstract 1
- 238000006243 chemical reaction Methods 0.000 description 19
- 238000010586 diagram Methods 0.000 description 11
- 238000005070 sampling Methods 0.000 description 11
- 238000012545 processing Methods 0.000 description 10
- 238000004364 calculation method Methods 0.000 description 8
- 238000012549 training Methods 0.000 description 7
- 230000010363 phase shift Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000005562 fading Methods 0.000 description 4
- 239000000654 additive Substances 0.000 description 3
- 230000000996 additive effect Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000021615 conjugation Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 241000695274 Processa Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
Images
Landscapes
- Mobile Radio Communication Systems (AREA)
Abstract
本发明公开了一种进行载波频率偏移估计的方法及装置。基站对接收信号进行二维联合最大似然估计,得到信道冲击响应和载波频率偏移的向量泛数;根据最小二乘原则以及所述向量泛数对载波频率偏移进行估计,适用于任意子载波分配方案下的载波频率偏移进行估计,对用户间的子载波分配方式无要求,不再像现有技术中必须为用户分配连续子载波,对子载波的规则分配和非规则分配均适用。由于本发明提供的实施例支持子载波的无规则分配方式,因此,***可以根据信道情况为用户分配较好的子载波即传输信道,以获得更多的频率分集增益,进一步降低误码率,增大传输速率。
Description
技术领域
本发明涉及移动通信领域,特别是指一种进行载波频率偏移(CarrierFrequency Offset,CFO)估计的方法及装置。
背景技术
正交频分复用(Orthogonal Frequency Division Multiplexing,OFDM)传输技术不仅具有较高的频谱效率,而且能够有效对抗频率选择性衰落,适合于宽带传输,被广泛应用于各种无线宽带通信***中。OFDM是一种特殊的多载波调制方式,其主要思想就是在频域内将总的信道划分成多个子信道,不同的子信道使用不同的子载波进行调制,各子载波之间彼此重叠且保持相互正交,并行进行传输。这样能够有效对抗由于多径时延造成的频率选择性衰落。
由于OFDM信号的各子载波是相互正交的,因此,OFDM***能够有效利用频谱资源,大大提高数据传输的性能和速率。但是,实际应用中,由于载波频率偏移及多普勒效应,破坏了OFDM信号的正交性,造成子载波之间的相互干扰,导致整个***的性能严重下降,因此,对载波频率偏移进行估计、进而根据载波频率偏移估计值进行补偿是非常必要的。
正交频分复用接入(Orthogonal Frequency Division Multiplexing Access,OFDMA)是以OFDM为基础的一种多址接入技术,OFDMA***与OFDM***具有相同的优势,如高频谱利用率、支持高速传输和良好的抗多径衰落能力。OFDMA***与OFDM***相同,同样对载波频率偏移非常敏感,频率载波频率偏移及多普勒效应使子载波之间的正交性遭到破坏,导致OFDMA***的性能变差,
OFDMA***中,由于是多用户接入,需要对子载波进行划分,将其划分成若干个子信道,然后将这些子信道分配给各用户使用。划分子信道的方式主要有三种。第一种方式是基于子带进行划分,将整个带宽划分成几个连续的子带,每个子带上的子载波与一个子信道相对应,这种方式实现简单,并且相邻子信道干扰较少,但是在利用频率分集方面效果较差。第二种方式是交织的子载波分配,即每个子信道的子载波等间隔分布在频域上,这种方式虽然能够获得比较好的频率分集增益,但是却会导致子信道之间的干扰较为敏感。第三种方式是对第二种方式的改进,不同子信道包含的子载波以伪随机的方式交织,每个子信道的子载波随机地分布在频域上,通过不同的基站使用不同的序列码来减少基站之间的干扰。
目前,OFDMA***中,多基于子载波连续分配的情况进行载波频率偏移估计。将整个子载波划分成若干个子带,并且每个子带的子载波是连续的,不同用户在接入时占用不同的子带。用户在连续两个OFDMA符号的相同子带位置上发送相同的训练符号。3个用户的情况下,如图1所示,每个用户的训练符号都在频域上占用了连续的子载波。用户终端收到时域信号后,通过傅立叶变换,可恢复出训练符号的频域值。由于受到载波频率偏移的影响,该训练符号的频域值与实际发送的值产生偏差,且连续两个OFDMA符号中的序列符号的偏差值是不一样的,这样,可根据序列符号在连续两个OFDMA符号中所占用的子载波和相位关系来对每个用户的载波频率偏移进行估计。由于该处理方式必须为每个用户分配连续的子载波,这就使得子载波的分配极不灵活,无法有效利用频率分集来提高OFDMA***的性能,并且由于每个用户需要占用连续两个OFDMA符号发送训练符号,使得***开销比较大。而对于子载波不连续分配的情况,目前还没有相应的载波频率偏移估计方案。
另外,在OFDMA***的上行链路中,基站需要同时处理多个用户的上行信号,会产生用户间干扰,使得性能进一步恶化,因此,根据载波频率偏移估计值进行载波频率偏移补偿就是非常重要的。一种处理方式是,基站对收到的信号直接乘以每个用户的载波频率偏移估计值,进行载波频率偏移补偿,但是,这样会存在多址间干扰,为了减小多址间干扰,需要在每个用户的子带之间***保护间隔,以减少干扰,而且这种处理方式只适用于为每个用户分配连续子载波的情况,使得子载波的分配极不灵活,无法有效利用频率分集来提高OFDMA***的性能。由于不同用户所处位置不同、移动速度不同等因素,使得不同用户与基站之间的载波频率偏移和多普勒效应各不相同,这对载波频率偏移补偿提出了更高的要求。另一种处理方式采用基站反馈法进行载波频率偏移补偿,即基站根据收到的信号对各个用户的上行链路进行载波频率偏移估计,然后向各用户终端返回载波频率偏移估计值,用户终端根据收到的载波频率偏移估计值,对发射信号的载波频率进行调整,以实现载波频率偏移补偿,减少由于载波频率偏移而造成的干扰。虽然,目前的方案能够对载波频率偏移进行补偿,但没有考虑用户间干扰,使得载波频率偏移补偿效果不佳。
发明内容
有鉴于此,本发明实施例提供一种进行载波频率偏移估计的方法及装置,灵活分配子载波,以有效利用频率分集来提高***性能。进一步地,本发明实施例提供了载波频率偏移补偿方案,有效消除用户间干扰。
本发明实施例提供的进行载波频率偏移估计的方法包括:基站对接收信号进行二维联合最大似然估计,得到信道冲击响应和载波频率偏移的向量泛数;根据最小二乘原则和所述向量泛数对载波频率偏移进行估计。
本发明实施例提供的进行载波频率偏移估计的装置包括载波频率偏移估计单元,用于对接收信号进行二维联合最大似然估计,得到信道冲击响应和载波频率偏移的向量泛数,根据最小二乘原则及所述向量泛数对载波频率偏移进行估计,得到载波频率偏移估计值。
本发明提供的实施例中,基站对接收信号进行二维联合最大似然估计,得到信道冲击响应和载波频率偏移的向量泛数,根据最小二乘原则及向量泛数对载波频率偏移进行估计,适用于任意子载波分配方案下的载波频率偏移进行估计,对用户间的子载波分配方式无要求,不再像现有技术中必须为用户分配连续子载波,对子载波的规则分配和非规则分配均适用。由于本发明提供的实施例支持子载波的无规则分配方式,因此,***可以根据信道情况为用户分配较好的子载波即传输信道,以获得更多的频率分集增益,进一步降低误码率,增大传输速率。
附图说明
图1为现有技术中训练符号频域分布示意图;
图2为本发明实施例中用户占用子载波频域分布示意图;
图3为本发明实施例中任意子载波分配方式下用户占用子载波情况示意图;
图4为本发明实施例中载波频率偏移估计流程图;
图5为本发明实施例中基站接收机结构示意图;
图6为本发明实施例中CFO估计单元结构示意图;
图7为本发明实施例中载波频率偏移补偿流程图;
图8为本发明实施例中CFO补偿单元结构示意图;
图9为本发明实施例中干扰消除单元结构示意图。
具体实施方式
本发明提供的实施例中,基站对接收信号进行二维联合最大似然估计,得到信道冲击响应和载波频率偏移的向量泛数,根据最小二乘原则及向量泛数对载波频率偏移进行估计。
下面结合附图对本发明提供的实施例的具体实现方案进行详细描述。
fc为载波中心频率,ΔF为子载波间隔,K个用户占用2M个数据子载波,2M个数据子载波被划分成Q个无交集的子信道,可记为Mq,其中,q=0,1,…,Q-1。每个子信道是为用户分配子载波的最小单位。划分得到的集合满足 且当q≠q′时,Mq∩Mq′为空集,即2M个数据子载波均被划分至子信道中,并且各子信道中包含的子载波不得有重复。每个用户所占用的子载波的数量和位置可根据各用户所处的环境进行最佳分配。如图2所示,可任意为用户分配子载波,分配给用户的子载波可以不是连续的,3个用户的情况下,只要保证各用户占用的子载波没有重合即可,一个用户占用的子载波无需连续,大大提高了子载波分配的灵活性,进而获得更多的频率分集增益。图3为5个接入用户的子载波分配示意,其中,实线为数据子载波,不同的箭头形状表示相应子载波被不同用户占用,虚线表示虚拟子载波,该子载波不进行传输。
第k个用户占用p(k)个子载波,占用的子载波序号 每个用户的信号经子载波映射后扩展成N点的频域信号,X(k)(n)为第k个用户发送的、且经四相移频键控(Quaternary Phase Shift Keying,QPSK)或正交幅度调制(Quadrature Amplitude Modulation,QAM)调制的频域信号,其中n∈{1,2,…,N}。X(k)(n)只在第k个用户所占用的子载波上有值,其他位置均为0,具体可参见图2。第k个用户的载波频率分配表为 其中本发明提供的实施例中,A(k)中取值为1的位置可连续,也可不连续。X(k)(n)经N点快速傅立叶逆变换(Inverse Fast Fourier Transform,IFFT)变换后变为x(k)(n)。无线信道是多径衰落信道,可采用时域冲击响应来描述,这样,第k个用户与基站之间的信道时域冲击响应h(k)=[h(k)(0),h(k)(1),…,h(k)(L(k))]T,这是对无线信道的表述,其中L(k)为冲击响应长度,对应于载波频率分配表A(k),信道时域冲击响应存在取值的位置在子载波的任意位置上,可以连续,也可以不连续。
理想情况下,即没有载波频率偏差、没有时间同步误差、没有基站的接收机噪声等情况下,第k个用户的信号x(k)(n)经过信道后,基站得到的时域接收信号可表示为用户信号与信道时域冲击响应的卷积,可表示为
当存在载波频率偏移时,ε(k)为第k个用户的归一化载波频率偏移,即载波频率偏移与子载波间隔ΔF之比,在该归一化载波频率偏移下,每个时域采样间隔所对应的相位偏移量ω(k)=2πε(k)/N,第n个采样点所对应的相位偏移量为nω(k),因此,基站收到的第k个用户的第n个采样点的信号为
由于多用户同时接入,因此,基站收到的第n个采样点的信号是所有用户的第n个采样点的信号及噪声的叠加,即基站收到的第n个采样点的信号为K个用户发送的信号与噪声之和,可表示为:
其中,z(n)为高斯白噪声。
将公式(1)表示为矩阵形式,得到
其中,r=[r(0),r(1),…,r(N-1)]T,用于表示基站收到的时域信号在N个采样点位置上的值; 用于表示第k个用户时域信号不同采样点位置所对应的由载波频率偏移而产生的相位偏移量,这是因为相位偏移量与采用点位置和载波频率偏移均有关;Ak=[Ak]p,q=x(k)(|p-q|N),为公式(1)中第k个用户发送的信号x(k)(n-l)的矩阵表示形式,p为Ak的行数,q为Ak的列数; 为第k个用户所对应的信道时域冲击响应,需要在h(k)的基础上填充Ncp-L(k)个0;z=[z(0),z(1),…,z(N-1)]T,用于表示基站接收信号所对应的加性高斯白噪声。
公式(2)可进一步表示为:
r=Q(ω)ξ+z (3)
其中,Q(ω)=[Γ(ω1)A1,Γ(ω2)A2,…,Γ(ωK)AK],用于表示已经预先包含了载波频率偏移但未经信道发送的信号,ω即为载波频率偏移;ξ=[ξ1,ξ2,…,ξK]T,用于表示与所有用户相对应的信道时域冲击响应。通过Ak能够表示出收发双方约定的导频信号。由于分配给用户终端的子载波可以连续也可以不连续,因此,导频信号所占用的子载波相应地连续或不连续,且分配给用户终端的子载波不连续时,分配给用户终端的各子载波之间的间隔可没有特定规律。
对于基站而言,载波频率偏移信息ω和信道信息h(k)均是未知的,因此,利用二维联合最大似然估计可得:
在实际***中,通常情况下认为一个用户的ω在一帧数据长度内的变化是非常小的,因此,可假定为固定值,又由于公式(4)中 均是未知的,如果直接计算,计算量将会非常大,因此,首先假设是已知的,则只有是未知的,此时,要使公式(4)的等号右边得到最小值,则根据最小二乘原则得到:
另外,如果用户终端在发射信号之前为了减少码间干扰而增加了循环前缀,则基站收到信号后,在进行载波频率偏移估计之前,会先去除循环前缀。
根据公式(6)就可得到每个用户的载波频率偏移估计。但是,由于公式(6)需要计算的是包含K个载波频率偏移的表达式求最值的问题,计算复杂度非常大,因此本发明实施例通过迭代算法得到K个载波频率偏移,具体处理流程如图4所示:
步骤402:k=1时,将作为变量,除以外的其余均取初始值,如均取0,将的取值代入公式(6),得到的取值。得到取值的具体处理过程可为:根据载波频率偏移可能取值的范围,依次以步长为单位,搜索到能够使最大的
步骤403:令k=k+1,将作为变量,除以外,的取值均为各自的更新值,的取值均取初始值,将除以外的其余的取值代入公式(6),得到的取值。得到取值的具体处理过程可为:根据载波频率偏移可能取值的范围,确定初始值,然后依次以步长为单位,搜索到能够使最大的
步骤404:判断k是否等于K,如果是,则继续执行步骤405;否则,返回执行步骤403。
步骤405:判断是否收敛,即各的取值是否逐渐趋向于一个取值,如果是,则完成载波频率偏移的估计;否则,继续执行步骤406。判断是否收敛可通过本次迭代过程得到的与上次迭代过程得到的之间的差值是否小于设定值来实现,如果本次迭代过程得到的与上次迭代过程得到的之间的差值小于设定值,则表明收敛。如果本次迭代为第一次迭代,则上次迭代过程得到的可取初始值。
步骤406:令k=1,将作为变量,除以外的其余均取上次迭代过程中得到的各的值,即的取值为上次迭代过程得到的值,将除外的其余的取值代入公式(6),得到的取值。得到取值的具体处理过程可为:根据载波频率偏移可能取值的范围,确定初始值,然后依次以步长为单位,搜索到能够使最大的
步骤407:令k=k+1,将作为变量,除以外,的取值均为各自的更新值,的取值为上次迭代过程得到的值,将除以外的其余的取值代入公式(6),得到的取值。得到取值的具体处理过程可为:根据载波频率偏移可能取值的范围,确定初始值,然后依次以步长为单位,搜索到能够使最大的
步骤408:判断k是否等于K,如果是,则返回执行步骤405;否则,返回执行步骤407。
如图4所示的迭代过程,虽然需要多步迭代才能够得到收敛的,但每次的运算量都比较小,采用迭代法进行载波频率偏移估计通常比直接进行最大似然进行载波频率偏移估计要大大节省运算时间,通过仿真得到二者所需运算时间的比值大概为1∶15,而且,一般只需3到5次的循环过程就能够得到收敛的,因此,节省了运算时间的同时还大大降低了运算量。
图5为本发明实施例中基站接收机结构示意图,如图5所示,基站接收机包括接收天线、下变频单元、模数(AD)变换单元、串并变换单元、傅立叶变换单元、CFO估计单元、CFO补偿单元、并串变换单元和星座点解映射单元。其中,接收天线用于接收用户终端发射的射频信号,并向下变频单元提供该射频信号;下变频单元用于对收到的射频信号进行下变频得到基带信号,并向AD变换单元提供该基带信号;AD变换单元用于对收到的基带信号进行A/D变换,得到离散数字信号,并向串并变换单元提供该离散数字信号;串并变换单元用于将串行数据变换为并行数据、即采样点信号,并向傅立叶变换单元和CFO估计单元分别提供采样点信号;傅立叶变换单元用于对收到的采样点信号进行傅立叶变换得到频域信号,并向CFO补偿单元提供该频域信号,傅立叶变换可通过离散傅立叶变换(Discrete Fourier Transform,DFT)运算或快速傅立叶变换(Fast Fourier Transform,FFT)运算实现;CFO估计单元用于根据收到的接收信号和包含在接收信号中、收发双方约定的导频信号对载波频率偏移进行估计,然后向CFO补偿单元提供载波频率偏移估计值;CFO补偿单元用于根据收到的载波频率偏移估计值对收到的频域信号进行载波频率偏移补偿,还原出用户终端发射的信号,然后向并串变换单元提供该还原信号;并串变化单元用于对收到的信号进行并串变换,得到串行数据,然后向星座点解映射单元提供该串行数据;星座点解映射单元用于对收到的串行数据进行星座点解映射。AD变换单元与串并变换单元之间可进一步包括循环前缀去除单元,用于对收到的离散数字信号去除循环前缀,并向串并变换单元提供去除循环前缀的离散数字信号。CFO估计单元的具体处理为:对接收信号进行二维联合最大似然估计,得到信道冲击响应和载波频率偏移的向量泛数,根据最小二乘原则及所述向量泛数对载波频率偏移进行估计,得到载波频率偏移估计值。
图6为本发明实施例中CFO估计单元结构示意图,如图6所示,CFO估计单元包括生成器和估计器,其中,生成器用于生成载波频率偏移的可能取值范围及步长,并向估计器提供该取值范围及步长;估计器用于对接收信号进行二维联合最大似然估计,得到信道冲击响应和载波频率偏移的向量泛数,根据最小二乘原则,向量泛数、及载波频率偏移的可能取值范围和步长,对载波频率偏移进行估计,得到载波频率偏移估计值,然后向CFO补偿单元提供载波频率偏移估计值。理论上,载波频率偏移的取值范围可从负无穷到正无穷,但是,生成器可根据实际***的实现方式及对载波频率偏移估计值的精度要求,确定载波频率偏移的可能取值范围及步长,以降低计算复杂度。
根据以上描述可见,本发明提供的实施例可以适用于任意子载波分配方案下的载波频率偏移进行估计,对用户间的子载波分配方式没有任何要求,不再像现有技术中必须为用户分配连续子载波,对子载波的规则分配和非规则分配均适用。由于本发明提供的实施例支持子载波的无规则分配方式,因此,***可以根据信道情况为用户分配较好的子载波即传输信道,以获得更多的频率分集增益,进一步降低误码率,增大传输速率。
得到载波频率偏移估计值后,可进一步根据该载波频率偏移估计值进行载波频率偏移补偿。
如果不考虑噪声的影响,第k个用户的信号经过信道后,基站收到的信号可表示为:
用户终端在发射信号之前增加循环前缀时,基站去除循环前缀后,接收信号与用户信号之间的关系仍然如公式(7)所示。
这样,公式(1)可表示为: 其中,r(n)为K个用户发送的信号与噪声之和,经过傅立叶变换后得到:
其中,r=[r(0),r(1),…,r(N-1)]T, 为相位偏移量,即第k个用户的信号不同采样点位置所对应的由载波频率偏移而产生的相位偏移量,z=[z(0),z(1),…,z(N-1)]T,Y(k)=DFT(y(k)),C(k)=DFT(c(k)),Z=DFT(z)。
公式(8)可进一步表示为:
其中,C′(k)=DFT(c′(k)), 是c(k)相位偏移量的补偿量,通过对c(k)进行共轭得到。
利用载波频率分配表为 和R,第k个用户的信号Y′(k)可表示为:
Y′(k)=A(k)R (11)
将公式(11)、公式(9)均代入公式(10)可得: 最后得到:
公式(12)中,为第k个用户的信号,为用户间干扰,是加性噪声。通过公式(12)可见,该公式并不能消除多用户干扰(Multi-User Interference,MUI)。本发明提供的实施例中,利用多用户检测和干扰消除理论,多次循环以进一步消除用户间干扰,Y(k),j表示第j次迭代时第k个用户的频域信号,具体处理流程如图7所示:
步骤701:j=0时, 由于C′(i)与C(i)共轭,因此,通过A(i)R与补偿向量C′(i)的卷积可消除C(i)中载波频率偏移的影响,实现载波频率偏移补偿,其中,i=1,2,…,K,从而对接收的K个用户中每个用户的信号进行载波频率偏移补偿。
单纯地进行了载波频率偏移补偿之后,用户间干扰仍然存在,因此,可进一步消除用户间干扰,具体处理过程见以下描述。
步骤702:令j=j+1,根据得到的Y(i),j-1消除用户间干扰,即对于第k个用户而言,根据公式(9)可得 Y(i),j-1为j-1次迭代过程中接收的第k个用户的已经进行过载波频率偏移补偿的信号,Y′(k),j为本次、即第j次迭代过程中接收的第k个用户的未经载波频率偏移补偿的信号。
步骤703:根据得到的Y′(k),j重新估计第k个用户的接收信号,根据公式(12)的第一部分可得
步骤704:判断Y(k),j是否收敛,如果是,则继续执行步骤705;否则,返回执行步骤702,进行再一次的迭代。判断Y(k),j是否收敛可采用期望最小化的方式得到,即E{||Y(k),j-Y(k),j-1||2}<Δ,其中,Δ为收敛精度,如Δ为1e-5、1e-6、1e-7等期望的值,如果满足条件时、即E{||Y(k),j-Y(k),j-1||2}<Δ成立时,则认为Y(k),j收敛,不再继续迭代过程;否则,继续进行迭代,即返回执行步骤702,直至得到满足条件的Y(k),j。
步骤705:对第k个用户的信号进行输出。
以上过程是针对K个用户中的每个用户进行的。
通常,经过4至5次的迭代过程就能够得到收敛的Y(k),j。在上述迭代过程中引入了圆周卷积,一般认为圆周卷积的复杂性高于傅立叶变换处理,从表面上看增加了迭代过程的复杂性。但是在步骤702和步骤703中,Y(k),j和A(k)Y′(k),j中最多只有p(k)个非零元素,并且当归一化载波频率偏移ε(k)较小时,如|ε(k)|<0.5,C(k)与C′(k)中的大多数元素均非常小,可以用零来代替。如果C′(k)={C′(k)(0),C′(k)(1),…,C′(k)(N-1)},则 由此可见,序列C′(k)中,C′(k)(N/2)附近的元素的值都接近于0,因此,W为序列C′(k)中非零元素的数量,W为奇数时,可采用序列 来替代C′(k),采用序列 替代C(k);W为偶数时,可采用序列 来替代C′(k),采用序列 替代C(k)。这样,大大降低了圆周卷积的复杂性。C(k)(m)中非零的取值只有W个,因此,步骤702中只需要做W次乘法,又由于第k个用户只在p(k)个子载波上传输数据,其他子载波上没有该用户的数据,因此,圆周卷积的乘法数大约是W·p(k),W很小时,运算量较小。
图8为本发明实施例中CFO补偿单元结构示意图,如图8所示,CFO单元对接收信号r进行傅立叶变换得到R,对所述r中的相位偏移量的补偿量进行傅立叶变换得到补偿向量C′(i),根据补偿向量C′(i)和载波频率偏移估计值对收到的频域信号进行载波频率偏移补偿后,还可进一步用于根据经过频率偏移补偿的信号进行用户间干扰的消除,这样,CFO补偿单元包括补偿单元和干扰消除单元,其中,补偿单元用于对接收信号r进行傅立叶变换得到R,对所述r中的相位偏移量的补偿量进行傅立叶变换得到补偿向量C′(i),根据所述补偿向量C′(i)和收到的载波频率偏移估计值对每个用户的信号进行载波频率偏移补偿,并向干扰消除单元提供经过载波频率偏移补偿的信号;干扰消除单元用于将收到的经过载波频率偏移补偿的信号经过迭代过程得到满足设定条件的已经进行过载波频率偏移补偿的信号。补偿单元收到的信号为r。
当干扰消除单元采用图7所示的迭代过程消除用户间干扰时,如图9所示,干扰消除单元包括消除单元、估计单元和判决单元,其中,消除单元用于根据收到的用户信号消除用户间干扰,并向估计单元提供消除用户间干扰的信号;估计单元用于根据收到的信号重新估计用户信号,并向判决单元提供该用户信号;判决单元用于确定收到的用户信号满足设定条件时,对用户信号进行输出,确定收到的用户信号不满足设定条件时,向消除单元提供该用户信号。消除单元第一次收到的用户信号来自于补偿单元,后续收到的用户信号则来自于判决单元。
因此,本发明实施例中提供的载波频率偏移补偿方案,不仅能够实现载波频率偏移的补偿,还能够消除用户间干扰。
以上提供的方案可适用于任何多载波***,如OFDM***、OFDMA***等。
本发明提供的实施例中,基站对接收信号进行二维联合最大似然估计,得到信道冲击响应和载波频率偏移的向量泛数,根据最小二乘原则及向量泛数对载波频率偏移进行估计,适用于任意子载波分配方案下的载波频率偏移估计,对用户间的子载波分配方式无要求,不再像现有技术中必须为用户分配连续子载波,对子载波的规则分配和非规则分配均适用。由于本发明提供的实施例支持子载波的无规则分配方式,因此,***可以根据信道情况为用户分配较好的子载波即传输信道,以获得更多的频率分集增益,进一步降低误码率,增大传输速率。
另外,本发明提供的实施例中,使用迭代的方式对各用户的载波频率偏移进行估计,大大降低了运算复杂度,缩短了运算时间。根据发明实施例提供的方案,用户的训练序列仅需要一个符号,节省了***开销。
此外,本发明提供的实施例中,基站根据载波频率偏移进行载波频率偏移补偿时,充分考虑了用户间干扰,不仅能够实现载波频率偏移的补偿,还能够消除用户间干扰,获得更准确的用户信号估计值。在消除用户间干扰过程中引入了参数W,降低了圆周卷积的复杂度,弥补了采用圆周卷积相对复杂的缺点,易于硬件的实现。本发明实施例提供的载波频率偏移补偿方案不仅适用于用户低速移动的情况,同样也适用于用户高速移动的情况。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (18)
1. 一种进行载波频率偏移估计的方法,其特征在于,该方法包含:
基站对接收信号进行二维联合最大似然估计,得到信道冲击响应和载波频率偏移的向量泛数;
根据最小二乘原则及所述向量泛数对载波频率偏移进行估计。
2. 根据权利要求1所述的方法,其特征在于,所述接收信号为r=Q(ω)ξ+z,其中,Q(ω)为已经预先包含了载波频率偏移但未经信道发送的信号,ω为载波频率偏移,ξ为与所有用户相对应的信道冲击响应,z为基站接收信号所对应的噪声。
c1、k等于用户数量K时,继续执行步骤d1,k小于用户数量K时,返回执行步骤b1;
g1、k等于用户数量K时,则返回执行步骤d1,k小于用户数量K时,返回执行步骤f1。
7. 根据权利要求1至6任一所述的方法,其特征在于,
所述基站接收信号,之前进一步包括:用户终端在发射信号之前加循环前缀;
所述基站对载波频率偏移进行估计,之前进一步包括:基站去除循环前缀。
8. 根据权利要求1至6任一所述的方法,其特征在于,所述信道冲击响应存在取值的位置在子载波的任意位置上。
9. 根据权利要求1至6任一所述的方法,其特征在于,所述对载波频率偏移进行估计,之后进一步包括,基站根据得到的载波频率偏移估计值进行载波频率偏移补偿:
对接收信号r进行傅立叶变换得到R;
对所述r中的相位偏移量的补偿量进行傅立叶变换得到补偿向量C′(i);
根据所述载波频率偏移估计值和所述补偿向量C′(i)进行载波频率偏移补偿。
10. 根据权利要求9所述的方法,其特征在于,所述根据载波频率偏移估计值和补偿向量C′(i)进行频率偏移补偿为 其中,A(i)为分配给用户终端的载波频率分配表,i=1,2,…,K。
11. 根据权利要求9所述的方法,其特征在于,所述根据载波频率偏移估计值和补偿向量C′(i)进行载波频率偏移补偿包括以下迭代过程:
a2、j=0时, 其中,j为迭代次数;
b2、j=j+1时,对于第k个用户,第j次迭代过程中接收的第k个用户的未经载波频率偏移补偿的信号 其中,Y(i),j-1为第j-1次迭代过程中接收的第k个用户的已经进行过载波频率偏移补偿的信号;
c2、根据所述Y′(k),j重新估计接收的第k个用户的信号
d2、确定Y(k),j满足设定条件时,对第k个用户的信号进行输出,确定Y(k),j不满足设定条件时,返回执行步骤b2。
12. 根据权利要求11所述的方法,其特征在于,所述Y(k),j满足的设定条件,是:E{||Y(k),j-Y(k),j-1||2}<Δ。
13. 根据权利要求11所述的方法,其特征在于,所述C(i)或C′(k)中包含W个非零元素。
14. 一种进行载波频率偏移估计的装置,其特征在于,该装置包括载波频率偏移估计单元,用于对接收信号进行二维联合最大似然估计,得到信道冲击响应和载波频率偏移的向量泛数,根据最小二乘原则及所述向量泛数对载波频率偏移进行估计,得到载波频率偏移估计值。
15. 根据权利要求14所述的装置,其特征在于,所述载波频率偏移估计单元包括生成器和估计器,其中,
所述生成器用于生成载波频率偏移的取值范围及步长,并向估计器提供该取值范围及步长;
所述估计器用于对接收信号进行二维联合最大似然估计,得到信道冲击响应和载波频率偏移的向量泛数,根据最小二乘原则、所述向量泛数、及所述载波频率偏移的取值范围和步长,对载波频率偏移进行估计,得到载波频率偏移估计值。
16. 根据权利要求14或15所述的装置,其特征在于,所述装置进一步包括:载波频率偏移补偿单元,用于对接收信号r进行傅立叶变换得到R,对所述r中的相位偏移量的补偿量进行傅立叶变换得到补偿向量C′(i),根据所述补偿向量C′(i)和收到的载波频率偏移估计值对收到的信号进行载波频率偏移补偿。
17. 根据权利要求16所述的装置,其特征在于,所述载波频率偏移补偿单元包括:补偿单元和干扰消除单元,其中,
所述补偿单元用于对接收信号r进行傅立叶变换得到R,对所述r中的相位偏移量的补偿量进行傅立叶变换得到补偿向量C′(i),根据所述补偿向量C′(i)和收到的载波频率偏移估计值对每个用户的信号进行载波频率偏移补偿,并向所述干扰消除单元提供经过载波频率偏移补偿的信号;
所述干扰消除单元用于将收到的经过载波频率偏移补偿的信号经过迭代过程得到满足设定条件的已经进行过载波频率偏移补偿的信号。
18. 根据权利要求17所述的装置,其特征在于,所述干扰消除单元包括:消除单元、估计单元和判决单元,其中,
所述消除单元用于根据收到的信号消除用户间干扰,并向估计单元提供消除用户间干扰的信号;
所述估计单元用于根据收到的信号重新估计用户终端的信号,并向判决单元提供该信号;
所述判决单元用于确定收到的信号满足设定条件时,对用户终端的信号进行输出,确定收到的用户信号不满足设定条件时,向消除单元提供该信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007100792496A CN101247375B (zh) | 2007-02-13 | 2007-02-13 | 一种进行载波频率偏移估计的方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007100792496A CN101247375B (zh) | 2007-02-13 | 2007-02-13 | 一种进行载波频率偏移估计的方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101247375A true CN101247375A (zh) | 2008-08-20 |
CN101247375B CN101247375B (zh) | 2012-08-08 |
Family
ID=39947582
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2007100792496A Active CN101247375B (zh) | 2007-02-13 | 2007-02-13 | 一种进行载波频率偏移估计的方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101247375B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010057382A1 (zh) * | 2008-11-21 | 2010-05-27 | 北京天碁科技有限公司 | 一种用于无线通信***的频偏估计方法及装置 |
CN102857465A (zh) * | 2012-07-25 | 2013-01-02 | 北京邮电大学 | 一种ofdm***同步误差的补偿方法及装置 |
CN103780518A (zh) * | 2012-10-25 | 2014-05-07 | 意法-爱立信有限公司 | 频率偏移估计方法及装置 |
CN102224718B (zh) * | 2008-11-21 | 2014-12-17 | 爱立信电话股份有限公司 | 对于时间和/或空间相关信号的盲频率偏移估计 |
CN106559375A (zh) * | 2016-12-05 | 2017-04-05 | 北京邮电大学 | Ofdm***中多普勒频偏估计方法及装置和消除方法及装置 |
CN107819716A (zh) * | 2017-10-31 | 2018-03-20 | 锐捷网络股份有限公司 | 一种基于频域的频偏补偿方法及设备 |
CN110247873A (zh) * | 2019-06-17 | 2019-09-17 | 中山大学 | 干扰消除交织ofdma上行链路cfo估计方法 |
CN110417699A (zh) * | 2019-05-30 | 2019-11-05 | 北京邮电大学 | 一种基于机器学习的ofdm***符号定时同步的方法 |
CN112272959A (zh) * | 2018-04-17 | 2021-01-26 | 株式会社Ntt都科摩 | 用户终端 |
CN112910821A (zh) * | 2021-02-01 | 2021-06-04 | 成都爱瑞无线科技有限公司 | 一种增加载波频率偏移估计可靠性和估计范围的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100518157C (zh) * | 2004-08-17 | 2009-07-22 | 中兴通讯股份有限公司 | 一种移动通信***频偏补偿的方法和装置 |
-
2007
- 2007-02-13 CN CN2007100792496A patent/CN101247375B/zh active Active
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010057382A1 (zh) * | 2008-11-21 | 2010-05-27 | 北京天碁科技有限公司 | 一种用于无线通信***的频偏估计方法及装置 |
CN101404518B (zh) * | 2008-11-21 | 2012-07-04 | 北京天碁科技有限公司 | 一种用于无线通信***的频偏估计方法及装置 |
CN102224718B (zh) * | 2008-11-21 | 2014-12-17 | 爱立信电话股份有限公司 | 对于时间和/或空间相关信号的盲频率偏移估计 |
CN102857465A (zh) * | 2012-07-25 | 2013-01-02 | 北京邮电大学 | 一种ofdm***同步误差的补偿方法及装置 |
CN102857465B (zh) * | 2012-07-25 | 2015-06-17 | 无锡北邮感知技术产业研究院有限公司 | 一种ofdm***同步误差的补偿方法及装置 |
CN103780518A (zh) * | 2012-10-25 | 2014-05-07 | 意法-爱立信有限公司 | 频率偏移估计方法及装置 |
CN106559375A (zh) * | 2016-12-05 | 2017-04-05 | 北京邮电大学 | Ofdm***中多普勒频偏估计方法及装置和消除方法及装置 |
CN106559375B (zh) * | 2016-12-05 | 2019-08-02 | 北京邮电大学 | Ofdm***中多普勒频偏估计方法及装置和消除方法及装置 |
CN107819716A (zh) * | 2017-10-31 | 2018-03-20 | 锐捷网络股份有限公司 | 一种基于频域的频偏补偿方法及设备 |
CN107819716B (zh) * | 2017-10-31 | 2021-01-05 | 锐捷网络股份有限公司 | 一种基于频域的频偏补偿方法及设备 |
CN112272959A (zh) * | 2018-04-17 | 2021-01-26 | 株式会社Ntt都科摩 | 用户终端 |
CN112272959B (zh) * | 2018-04-17 | 2024-03-01 | 株式会社Ntt都科摩 | 用户终端 |
CN110417699A (zh) * | 2019-05-30 | 2019-11-05 | 北京邮电大学 | 一种基于机器学习的ofdm***符号定时同步的方法 |
CN110247873A (zh) * | 2019-06-17 | 2019-09-17 | 中山大学 | 干扰消除交织ofdma上行链路cfo估计方法 |
CN112910821A (zh) * | 2021-02-01 | 2021-06-04 | 成都爱瑞无线科技有限公司 | 一种增加载波频率偏移估计可靠性和估计范围的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101247375B (zh) | 2012-08-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101247375B (zh) | 一种进行载波频率偏移估计的方法及装置 | |
CN109391579B (zh) | 用于共同相位误差及载波间干扰估测及补偿的方法 | |
US6826240B1 (en) | Method and device for multi-user channel estimation | |
US6362781B1 (en) | Method and device for adaptive antenna combining weights | |
JP5497800B2 (ja) | 直交周波数分割多重無線通信システムにおけるキャリア周波数オフセットを補償するための方法及び装置 | |
CN1838656B (zh) | 上行链路载波频率同步和天线加权向量估计的装置和方法 | |
EP1473862B1 (en) | Apparatus and method for transmitting training symbol groups in an OFDM communications system using multiple antennas | |
CN102113286B (zh) | 用于多载波***中的ici消除的迭代信道估计方法和设备 | |
JP4198428B2 (ja) | 無線伝送装置 | |
US8588204B2 (en) | Efficient channel estimation method using superimposed training for equalization in uplink OFDMA systems | |
CN101577692A (zh) | 一种正交频分复用***的信道估计方法和装置 | |
CN102113253B (zh) | 用于多载波***中的ici消除的无线电信道模型 | |
JP7204950B2 (ja) | チャネル及び位相雑音の同時推定のための巡回パイロットシーケンス | |
US8411773B2 (en) | Simplified equalization scheme for distributed resource allocation in multi-carrier systems | |
KR20170016653A (ko) | 비직교 파형을 갖는 다중 반송파 시스템에서 간섭 성분을 구성하고 채널을 추정하는 방법 및 장치 | |
CN100484117C (zh) | 交织ofdma上行链路***的载波频偏校正方法 | |
CN101247383A (zh) | 一种进行载波频率偏移补偿的方法及装置 | |
US20060251037A1 (en) | Apparatus and method for performance improvement of channel estimation in broadband wireless access communication system | |
Fan et al. | Channel Estimation and Interference Cancellation for OFDM Systems Based on Total Least Squares Solution. | |
KR20170089738A (ko) | 무선 통신 시스템에서 위상 에러를 추정 및 보상하는 방법 및 장치 | |
CN101997805A (zh) | 导频符号的处理方法和装置 | |
CN101335551A (zh) | 基于dft-s-gmc***多天线分集方案的sinr估计方法 | |
Thabet et al. | Synchronization error reduction using guard-band allocation for wireless communication systems | |
Stemick et al. | OFDM-FDMA scheme for the uplink of a mobile communication system | |
CN101958874A (zh) | 基于角分复用的双正交频分多址蜂窝*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |