CN102986292A - 用于估计载波频率偏移的方法与装置 - Google Patents

用于估计载波频率偏移的方法与装置 Download PDF

Info

Publication number
CN102986292A
CN102986292A CN201180014043XA CN201180014043A CN102986292A CN 102986292 A CN102986292 A CN 102986292A CN 201180014043X A CN201180014043X A CN 201180014043XA CN 201180014043 A CN201180014043 A CN 201180014043A CN 102986292 A CN102986292 A CN 102986292A
Authority
CN
China
Prior art keywords
signal
time
carrier frequency
preorder
rotation
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.)
Pending
Application number
CN201180014043XA
Other languages
English (en)
Inventor
冯灏泓
孙素梅
何晋强
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Agency for Science Technology and Research Singapore
Original Assignee
Agency for Science Technology and Research Singapore
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Agency for Science Technology and Research Singapore filed Critical Agency for Science Technology and Research Singapore
Publication of CN102986292A publication Critical patent/CN102986292A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2657Carrier synchronisation
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2614Peak power aspects
    • H04L27/262Reduction thereof by selection of pilot symbols
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2614Peak power aspects
    • H04L27/2621Reduction thereof using phase offsets between subcarriers
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2689Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation
    • H04L27/2692Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation with preamble design, i.e. with negotiation of the synchronisation sequence with transmitter or sequence linked to the algorithm used at the receiver
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W56/00Synchronisation arrangements
    • H04W56/0035Synchronisation arrangements detecting errors in frequency or phase
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/14Relay systems
    • H04B7/15Active relay systems
    • H04B7/155Ground-based stations
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W24/00Supervisory, monitoring or testing arrangements
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W72/00Local resource management
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W88/00Devices specially adapted for wireless communication networks, e.g. terminals, base stations or access point devices
    • H04W88/02Terminal devices
    • H04W88/04Terminal devices adapted for relaying to or from another terminal or user

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Mobile Radio Communication Systems (AREA)
  • Synchronisation In Digital Transmission Systems (AREA)

Abstract

本发明公开了一种用于估计通信系统(100)中的载波频率偏移的方法(500)。在所述实施例中,该通信系统(100)包括第一通信装置(120)、第二通信装置(122)、以及中继(110)。该方法(500)包括:在该第一通信装置(120)处,以多个不同的第一旋转角度旋转(530)第一产生的信号区块以形成对应的第一前序信号集合(410);在该第二通信装置(122)处,以多个不同的第二旋转角度旋转(530)第二产生的信号区块以形成对应的第二前序信号集合(420);从该第一通信装置(120)和该第二通信装置(122)中的每个通信装置,将该各相应第一前序信号集合(410)和第二前序信号集合(420)作为时域信号发射(540)至该中继(110);在该第一通信装置(120)处,接收来自该中继(110)的经重传的时域信号,该经重传的信号为该第一前序信号集合(410)与该第二前序信号集合(420)的组合;以及基于该接收到的经重传的时域信号来估计(550)该信道频率偏移。还公开了一种用于估计通信系统中的载波频率偏移的装置以及集成电路。

Description

用于估计载波频率偏移的方法与装置
技术领域
本发明涉及一种用于估计载波频率偏移的方法与装置,尤其但是并非排他地,本发明涉及双向中继。 
背景技术
用于将协作分集引进无线传输的中继技术具有通信量和范围扩展的潜能。具体而言,双向中继是能够在单向中继中所需时间的仅一半时间内支持两个源节点之间的数据分组交换的高效网络传输方案。在双向中继中,两个源节点使用相同时槽和带宽发送独立数据流。然而,由于共享频谱资源,所以由任一源节点所传达的信号不仅交叉传送至另一源节点,而且还作为干扰经由中继节点传送回其本身。因此,当设计双向中继系统中用于频率同步化的前序信号时,考虑干扰。 
一种广泛应用于点对点传输中频率同步化的前序信号具有IEEE802.11a/g/n标准中采用的周期性结构。然而,当将此前序信号应用于双向中继设置时,该前序信号产生来自源节点的传输,一旦该传输在到达中继节点之前彼此混合,则可能无法区分该传输。当发生此状况时,当中继节点将此接收到的复合信号重新引导回两个源节点时,可能不再能够可靠地估计该两个源节点之间的载波频率偏移(CFO)。因此,目前用于点对点传输或者单向中继传输的标准前序信号可能并不能够直接应用于双向中继传输。 
已知为中继网络而设计的训练序列并且其可以用于CFO估计,但是这些训练序列通常用于达到信道估计的目的。例如,用于单向中继的基于恒定振幅零自相关(CAZAC)序列的周期性前序信号可以用于估计目的地与多个中继站之间的CFO。然而,在这种情况下,由于存在由CFO产生的载波间干扰,所以以高SNR观察到CFO估计中的误差极限(error floor)。同样地,此策略是为单向中继而设计的,并且当应用于双向中继时,性能降低。 
本发明的一个目的在于提供一种解决现有技术中的至少一个问题的用于估计CFO的方法与装置,和/或向公众提供有用的选择。 
发明内容
在本发明的具体表达中,提供一种用于在通信装置处估计载波频率偏移的方法,所述方法包括: 
产生信号区块; 
以多个不同旋转角度旋转所述区块以形成对应的第一前序信号集合; 
将所述第一前序信号集合作为时域信号发射至中继; 
接收来自所述中继的经重传的时域信号,所述经重传的时域信号为所述第一前序信号集合与来自另一通信装置的第二前序信号集合的组合;以及 
基于所述接收到的经重传的时域信号来估计所述信道频率偏移。 
所述方法可以进一步包括:将优化调制应用至所述第一前序信号集合以形成所述时域信号。优选地,所述信号区块的各个信号对应于用于发射所述第一前序信号集合的多个次载波中的各相应的次载波。 
有利地,所述方法可以进一步包括: 
排列多个组合,每个组合将所述旋转信号区块的布置与多个群集相关联; 
用所述多个次载波中的各相应次载波来调制所述多个组合中的每个组合;以及 
从所述经调制的多个组合中选择出具有最佳组合的被选信号,所述最佳组合使所述对应信号的峰值平均功率比最小化;其中所述被选信号为所述时域信号。 
优选地,所述方法可以进一步包括:确定由所述多个次载波所使用的多个调制符号以形成所述第一前序信号集合中的每一个。所述方法可以进一步包括:通过应用所述多个不同旋转角度中的一个旋转角度,来移位各个所述信号的频率。可选地,在所述通信装置处旋转所述产生的区块的步骤可以包括比例缩放所述信号中的每个信号的振幅。 
在变型例中,所述多个不同旋转角度中的一个旋转角度可以由所述多个次载波的数目以及所述时域信号的预定长度获得。可选地,所述多个不同旋转角度中的一个旋转角度可以由先前信号区块的角度获得。 
有利地,所述多个不同旋转角度中的一个旋转角度可以使用所述载波频率偏移估计的克拉美-罗界来确定。可选地,所述克拉美-罗界可以为近似值。优选地,所述多个不同旋转角度中的一个旋转角度可具有介于0与π之间并且包括0及π的值。 
在第二变型例中,所述多个不同旋转角度中的一个旋转角度可以根据其他通信装置的数目来获得。 
可选地,所产生的区块可以为IEEE 802.11前序信号。优选地,所述第一前序信号集合中的一个前序信号可以由旋转先前旋转信号区块获得。可选地,所述时域信号可以使用正交频分复用来发射。 
优选地,所述时域信号为非周期性的。优选地,所述经重传的时域信号可包括由所述中继从另一装置重传的训练信号。优选地,所述载波频率偏移介于所述时域信号与所述接收到的经重传的时域信号之间。 
有利地,估计所述载波频率偏移的步骤可以包括:对所述接收到的经重传的时域信号进行线性滤波。在又一变型例中,估计所述载波频率偏移的步骤可以进一步包括:在所述线性滤波信号上执行相关。 
在本发明的第二具体表达中,提供一种用于在第一通信装置处估计载波频率偏移的中继方法,所述方法包括: 
在中继处,接收从所述第一通信装置发射的第一时域信号,以及从第二通信装置发射的第二时域信号; 
其中所述第一时域信号包括第一前序信号集合,通过以对应的多个不同的第一旋转角度旋转第一产生的信号区块来形成所述第一前序信号集合;并且 
其中所述第二时域信号包括第二前序信号集合,通过以对应的多个不同的第二旋转角度旋转第二产生信号区块来形成所述第二前序信号集合;以及 
从所述中继将经重传的时域信号重传至所述第一通信装置,以允许所述第一通信装置基于所述经重传的时域信号来估计所述载波频率偏移,所述经重传的时域信号为所述第一前序信号集合与所述第二前序信号集合的组合。 
在本发明的第三具体表达中,提供一种用于估计通信系统中的载波频率偏移的方法,所述通信系统包括:第一通信装置、第二通信装置以及中继,所述方法包括: 
在所述第一通信装置处,以多个不同的第一旋转角度旋转第一产生的信号区块,以形成相应第一前序信号集合; 
在所述第二通信装置处,以多个不同的第二旋转角度旋转第二产生的信号区块,以形成相应第二前序信号集合; 
从所述第一通信装置和所述第二通信装置中的每个通信装置将所述各相应的第一前序信号集合和第二前序信号集合作为时域信号发射至所述中继; 
在所述第一通信装置处,接收来自所述中继的经重传的时域信号,所述经重传的信号为所述第一前序信号集合与所述第二前序信号集合的组合;以及 
基于所述接收到的经重传的时域信号来估计所述信道频率偏移。 
在第三具体表达中,有利地,所述多个不同的第一旋转角度的起始角度与所述多个不同的第二旋转角度的起始角度可以相差π。优选地,所述多个不同的第一旋转角度的所述起始角度为0。 
本发明还涉及一种用于执行上文所论述的任何方法或者在优选实施例中所阐述的那些方法的设备或者通信装置。具体而言,在本发明的第四具体表达中,提供一种通信装置,包括: 
处理器,其被配置以产生信号区块并且以多个不同的旋转角度旋转所述区块,以形成对应的第一前序信号集合;以及 
发射机,其被配置以将所述第一前序信号集合作为时域信号发射至中继; 
接收机,其被配置以接收来自所述中继的经重传的时域信号,所述经重传的时域信号为所述第一前序信号集合与来自另一通信装置的第二前序信号集合的组合;并且 
其中所述处理器进一步被配置以基于所述接收到的经重传的时域信号来估计所述信道频率偏移。 
在本发明的第五具体表达中,提供一种用于通信装置的集成电路,包括: 
处理单元,其被配置以产生信号区块并且以多个不同旋转角度旋转所述区块,以形成对应的第一前序信号集合; 
接口,其被配置以将所述第一前序信号集合作为时域信号发射至中继,并且进一步被配置以接收来自所述中继的经重传的时域信号,所述经重传的时域信号为所述第一前序信号集合与来自另一通信装置的第二前序信号集合的组合;并且 
其中所述处理单元进一步被配置以基于所述接收到的经重传的时域信号来估计所述信道频率偏移。 
能够从所描述的实施例理解到的是,所述方法与装置可以: 
-产生并且发射具有低峰值平均功率比(PAPR)的前序信号; 
-降低估计CFO的复杂度; 
-当估计CFO时,不同频率估计技术的使用具有稳健性; 
-产生精确的CFO估计; 
-具有接近最佳的CFO估计性能; 
-使用滤波器去除已知频率分量,这使CFO估计性能接近克拉美-罗界; 
-在不经受性能损失或者估计精度损失的情况下,使用具有低复杂度的估计器电路; 
-对中继站数目或者各站中的天线数目无任何限制; 
-将廉价的无失真电路用于CFO估计;并且 
-具有带有足以进行PAPR优化的自由度的前序信号结构。 
附图说明
现将参照附图仅以举例的方式来阐述本发明的优选实施例,附图中: 
图1为根据优选实施例的具有两个源节点以及中继节点的通信系统的示意图; 
图2为图1的源节点的发射部分的示意图; 
图3为图1的源节点的接收部分的示意图; 
图4为从图2的发射部分产生并且发射出的经旋转的前序信号的两个集合的示意图; 
图5为估计图1的通信系统中的CFO的方法的流程图; 
图6为在图2的发射部分处产生并且发射前序信号的方法的流程图; 
图7为随着CFO改变,不同数目的CFO估计区块的最小CFO估计误差值的曲线图; 
图8为图5的前序信号中的一个前序信号的时域波形的曲线图; 
图9为图5的前序信号中的一个前序信号的频域波形的曲线图; 
图10为随着SNR改变并且在应用方案1的信道条件的情况下,图5的CFO估计中的均方误差(MSE)的曲线图; 
图11为随着SNR改变并且在应用方案2的信道条件的情况下,图5的CFO估计中的MSE的曲线图;以及 
图12为随着SNR改变并且在应用方案3的信道条件的情况下,图5的CFO估计中的MSE的曲线图。 
具体实施方式
在本说明书中可能使用以下符号。大写粗体字母和小写字母分别表示矩阵和向量。除非另有说明,矩阵和向量中的所有索引均从零开始。符号 
Figure DEST_PATH_GDA00002455764000061
表示回旋。 表示具有零均值并且具有协方差矩阵 
Figure DEST_PATH_GDA00002455764000063
的多元高斯分布。 
系统模型 
图1示出根据优选实施例的通信系统100。通信系统100包括中继节点110和两个源节点,即源1 120和源2 122。中继节点110和源节点各自均能够进行双向中继通信。换言之,源1 120能够向中继节点110发射信号并且从中继节点110接收信号。类似地,源2 122也能够向中继节点110发射信号并且自中继节点110接收信号。 
变量SA表示节点A,而hAB表示自节点A至节点B的信道。rA表示在节点A处接收到的信号。A和B可以取值0、1和2,A和B取值0、1和2分别表示节点与中继110、源1 120以及源2 122相关联。中继节点110和源节点S1与S2使用表示为f0、f1和f2的载波频率来发射。源节点S1和S2可以使其载波频率f1和f2与中继节点110的载波频率,即f0对准。f0用作通用频率并且可以使用测距技术来实现源节点S1和S2处的的频率对准。 
在能够进行双向中继通信的系统中,可以在两个时槽上为S1和S2执行测距。在第一时槽(即,时间1)中,对于第n个离散时间取样,S1和S2经由信道h10,n和h20,n同时将其分别表示为x1,n和x2,n的数据包发射至中继节点110。当在中继节点110处接收到所发射的数据包x1,n和x2,n时,数据包重叠以形成接收到的信号r0,n。在第二时槽(即,时间2)中,中继节点110对之前在时间1期间接收到的信号r0,n进行比例缩放。随后,中继节点110将该经比例缩放的信号重传回两个源S1和S2。经重传的信号经由两个独立返回信道h01,n和h02,n而传输,以分别返回到S1和S2处。在S1和S2处,经重传的信号分别被接收为r1,n和r2,n。中继节点110可以被称为“响应机”,因为其通过将信号重传回源来“响应”自源接收到的信号。 
信道h10,n和h01,n可以不同,同样地,h20,n和h02,n也可以不同。所有信道可以使所发射的信号与加性高斯白噪声(AWGN)一起经受频率选择性衰落。当 执行测距时,在中继节点110处可以存在相对于频率f0的小残余频率差f1-f0和f2-f0。 
对于第n个离散时间取样而言,在时间2中S1接收到的离散信号可以以数学方式表示为: 
r 1 , n = α r 11 , n + α e j 2 π ( f 2 - f 1 ) n r 21 , n + u 1 , n - - - ( 1 )
其中 
r 11 , n ≡ h 11 , n ⊗ x 1 , n ,
r 21 , n ≡ h 21 , n ⊗ x 2 , n .
r11,n表示源自S1的重传回S1的分量信号,其包括源自S1的数据包x1,n。r21,n表示源自S2的分量信号,其包括经由中继向前发射至S1的数据包x2,n。应注意的是,也能够为S2在时间2中接收到的中继信号r2,n写出类似于方程式1的表达式。 
r1,n表示对于第n个离散时间取样而言在S1处接收到的信号。在r1,n中三个分量混合到一起。第一分量αr11,n表示源自S1的自中继节点110重传回S1的信号分量。该第一分量包括在时间1中由S1发射至中继节点110的消息x1,n,并且该第一分量可以被视为经由包括自S1至中继节点110并且自中继节点110回到S1的信道的复合信道h11,n而传输。α表示缩放比例因子,并且应注意的是,第一分量不具有载波频率偏移(CFO)。 
第二分量 
Figure DEST_PATH_GDA00002455764000074
表示源自S2的现经由中继节点110发射至S1的信号分量。该第二分量包括在时间1中S2已发送至中继节点110的消息x2,n,并且该第二分量可以被视为经由包括自S2至中继节点110并且自中继节点110至S1的信道的复合信道h21,n而传输。可以看出,第二分量经历缩放比例因子α,并且显著地经受载波频率偏移f2-f1。此CFO的量与在不存在中继节点110并且从S2至S1进行直接传输情况下的CFO的量相同。 
第三分量u1,n表示有色高斯噪声,第三分量u1,n的相关性为非时变的。 
源节点S1和S2
现参照图2,图2图示出图1的源节点S1和S2的发射部分200。发射部分200包括:处理器220,其被配置为从起始前序信号230产生包括前序信号 430、440的时域信号;以及天线210,其被配置为将前序信号430、440发射至中继节点110。可以在处理器220中产生前序信号430、440并且可以使用下文将描述的方法510来发射前序信号430、440。 
起始前序信号230为预定的,并且其可以被储存在发射部分200内的存储器中且随后被提供至处理器220。可选地,也可以使用算法在处理器220内产生起始前序信号230。起始前序信号230也可以取IEEE 802.11a/g/n标准中所定义的前序信号的值而预定。与物理层相关的IEEE 802.11a/g/n规范的内容均以引用的方式并入本文,即分别为:IEEE Std 802.11a-1999,Part11:Wireless LAN Medium Access Control(MAC)and Physical Layer(PHY)Specifications-High-speed Physical Layer in the 5GHz Band,IEEE,1999,IEEEStd 802.11g-2003,Part 11:Wireless LAN Medium Access Control(MAC)and Physical Layer(PHY)Specifications-Amendment 4:Further Higher Data Rate Extension in the 2.4GHz Band,IEEE,2003以及IEEE Std 802.11n-2009,Part 11:Wireless LAN Medium Access Control(MAC)and Physical Layer(PHY)Specifications Amendment 5:Enhancements for Higher Throughput。 
图3示出了图1的源节点S1和S2的接收部分300。接收部分300包括:天线310,其被配置为接收来自中继节点110的信号 
Figure DEST_PATH_GDA00002455764000081
以及处理器320,其被配置为根据接收到的信号 
Figure DEST_PATH_GDA00002455764000082
来估计载波频率偏移φ。在源节点S1处φ=f2-f1,而在源节点S2处φ=f1-f2。处理器320进一步包括逐块线性滤波器330和频率估计器340。逐块线性滤波器330从接收到的信号 
Figure DEST_PATH_GDA00002455764000083
中去除已知的频率分量。在步骤560中更详细地描述此滤波器330。因此,保留包括载波频率偏移φ的信号分量,并且可以使用频率估计器340来估计该载波频率偏移φ。频率估计器340可以采取基本相关器电路的形式。 
本领域技术人员应理解的是,尽管在此说明书中使用两个单独的天线210、310来描述源节点S1和S2的发射部分200和接收部分300,但是可以在源节点S1和S2中使用既能够发射又能够接收的单个天线来实施天线210、310。类似地,尽管描述了两个处理器220、320,但是应理解的是,可以使用单个处理器来产生经旋转的前序信号430、440以及估计载波频率偏移φ。 
使用前序信号的CFO估计的概述 
现参照图4和图5,可以使用前序信号来估计存在于方程式1的第二分量中的CFO f2-f1。图4示出在时域中,在时间1期间,分别从源节点S1和S2的发射部分产生并且发射的经旋转的前序信号430和440的两个集合410和420。图5示出估计图1的通信系统100中的CFO的方法500。 
集合410和420各自包括NBLK个前序信号,其分别为430和440。将各个前序信号430、440均旋转一角度。各个前序信号430、440均具有长度为L的取样,其中TS为取样间隔。逐个循序地发射各个集合410、420的前序信号430、440。应注意的是,各个连续的前序信号旋转的角度随时间而改变,并且对于S1和S2而言,各个连续的前序信号旋转的角度与紧接其先前的前序信号分别相差角度θ1和θ2。 
估计CFO的方法500可以分为两个部分。第一部分发生于时间1,并且涉及在源节点S1和S2中的每个源节点处产生将发射至中继节点110的经旋转的前序信号430和440。随后,第二部分发生于时间2中,并且涉及在S1和S2处接收回包含经旋转的前序信号430和440的信号,随后在接收到的信号上执行CFO估计。 
在步骤510中,在时间1期间,分别从源节点S1和S2产生并且发射经旋转的前序信号430和440。步骤510可以进一步包括下文将结合图6更详细描述的步骤520至步骤540。应注意的是,可选地,在执行方法500之前可以离线执行步骤520至步骤536。在源节点S1和S2处自起始前序信号230产生分别表示为x1,n和x2,n的前序信号430、440。x1,n和x2,n各自包含具有L个取样的NBLK个区块并且在时域中。将区块旋转角度kθ1和kθ2分别应用至x1,n和x2,n中的每一个,使得, 
x 1 , n + kL = e jk θ 1 x 1 , n , k = 0,1,2 , . . . , N BLK - 1 , - - - ( 2 )
x 2 , n + kL = e jk θ 2 x 2 , n , k = 0,1,2 , . . . , N BLK - 1 ,
k为区块索引,而θ1和θ2分别为在S1和S2中使用的区块旋转角度。 
如在图4中所示,经旋转的前序信号430、440的各个集合410、420分别包括NBLK个前序信号区块,其中各个前序信号430、440均为L个取样长。两个源的第一前序信号区块分别标示为 和 
Figure DEST_PATH_GDA00002455764000094
在时间1中,分别从源节点S1和S2发射前序信号430、440的两个集合410、420。 
在步骤550中,在时间2中,在源节点S1和S2处估计CFO。中继节点110将从源节点S1和S2发射的前序信号430、440重传回源节点S1和S2。将重传执行为来自中继节点110的广播,并且以与在时间1中当从S1和S2产生并且发射信号时相同的区块旋转格式在S1和S2处接收该重传。 
应注意的是,如在方程式1中所示,在S1和S2中所接收回的信号为从源节点S1和S2发射的信号的组合。在源节点S1与S2处接收到的信号 
Figure DEST_PATH_GDA00002455764000101
可以表示为: 
Figure DEST_PATH_GDA00002455764000102
Figure DEST_PATH_GDA00002455764000103
其中 
Figure DEST_PATH_GDA00002455764000104
Figure DEST_PATH_GDA00002455764000105
ρ 1 = e j φ 1 , ρ 2 = e j φ 2 ,
NCFO≤Nmax, 
Figure DEST_PATH_GDA00002455764000108
和 
Figure DEST_PATH_GDA00002455764000109
分别表示从源节点S1和S2发射的包括前序信号430和440的集合410和420的信号。 
Figure DEST_PATH_GDA000024557640001010
和 分别表示应用于 和 
Figure DEST_PATH_GDA000024557640001013
的旋转。 
Figure DEST_PATH_GDA000024557640001014
表示存在于 
Figure DEST_PATH_GDA000024557640001015
中的高斯噪声,并且其具有均值零和方差 
Figure DEST_PATH_GDA000024557640001016
Nmax为可用于CFO估计的区块的最大数目,该最大数目视实际实施例而预定并且比前序信号区块总数目NBLK小。NBLK为循环前缀去除和两个源节点S1和S2之间的时序未对准的结果。NCFO为用于CFO估计的实际区块数目,其可以不超过Nmax。另外,应注意的是,ρ2和ρ2表示在时间2中在相应的源S1和S2处感知到的信号 
Figure DEST_PATH_GDA000024557640001017
和 
Figure DEST_PATH_GDA000024557640001018
的相应CFO。 
Figure DEST_PATH_GDA000024557640001019
中的CFO在S1处已知,同样地, 
Figure DEST_PATH_GDA000024557640001020
中的CFO在S2处已知。 
因此,在源节点S1处, 
Figure DEST_PATH_GDA000024557640001021
Figure DEST_PATH_GDA000024557640001024
φ11, 
φ2=2π(f2-f1)L+θ2, 
和φ2未知,而φ1和 
Figure DEST_PATH_GDA00002455764000112
已知。NCFO已知。应注意的是,尽管消息x1,n已知,但是因为信道h11,n对于S1未知,所以 
Figure DEST_PATH_GDA00002455764000113
未知。 
在源节点S2处, 
Figure DEST_PATH_GDA00002455764000114
Figure DEST_PATH_GDA00002455764000115
Figure DEST_PATH_GDA00002455764000116
Figure DEST_PATH_GDA00002455764000117
φ12, 
φ2=2π(f1-f2)L+θ1
Figure DEST_PATH_GDA00002455764000118
和φ1未知,而φ2和 
Figure DEST_PATH_GDA00002455764000119
已知。NCFO已知。类似地,因为信道h12,n未知,所以 
Figure DEST_PATH_GDA000024557640001110
未知。在两个源节点S1和S2处,φ1和φ2表示所感知的区块旋转角度,而θ1和θ2分别为源节点S1和S2的实体区块旋转角度。实体区块旋转角度θ1和θ2表示发射所经历的实体CFO。随着前序信号410和420人工地将旋转θ1和θ2引入前序信号区块430和440,所感知的区块旋转φ1和φ2反映区块旋转θ1和θ2的组合效应和实体CFO。φ2为包括CFO的变量。 
随后,可以在接收到的信号 
Figure DEST_PATH_GDA000024557640001111
上执行线性滤波560,其后是由经线性滤波的信号进行CFO的频率估计570。在本说明书中,下文将更详细地描述这些步骤。 
前序信号产生和发射(步骤520至步骤540) 
随后,通过使用图6来描述前序信号430、440的产生。图6示出在图1的源节点S1和S2的发射部分200处产生和发射前序信号430、440的方法510。 
在520中,在S1和S2中提供起始前序信号230,并且确定区块旋转角度θ1和θ2以在S1和S2进行相应使用。在S1和S2中提供的起始前序信号230可以相同,或者它们可以不同。克拉美-罗界(CRB)可以用作确定θ1和θ2的最佳值的准则。可以使用CRB显示θ1=0和θ2=π的值为最佳值。因此,可以在两个源处将旋转角度θ1和θ2设定为相距π弧度。 
CFO的CRB反映所估计的CFO中的估计误差的最小统计可达值。在基于方程式3和方程式4估计CFO中的CRB可以通过以下方程式给出: 
Figure DEST_PATH_GDA000024557640001112
其中 
Figure DEST_PATH_GDA00002455764000121
并且 
Figure DEST_PATH_GDA00002455764000122
其中★为克洛涅克积算子。 
对于给定 
Figure DEST_PATH_GDA00002455764000123
而言,可以通过将方程式5代入方程式7来获得CFO估计f2-f1的CRB。举例而言,从S1的角度看,CFO估计的CRB将为: 
Figure DEST_PATH_GDA00002455764000124
其中 
Figure DEST_PATH_GDA00002455764000125
根据方程式5,φ2=2π(f2-f1)L+θ2。由于θ2为确定的并且事先已知,因此CFO f2-f1可以根据f2-f1=(mod(φ22,2π)-π)/(2πL)从φ2的估计获得,其中mod(A,2π)为A在2π上的模运算。注意,CFO必须落入|f1-f2|≤0.5/L,以无任何相位模糊性地估计该CFO。如在方程式7中所给出的CRB为变量 
Figure DEST_PATH_GDA00002455764000126
φ1和φ2的联合函数。因此,找出CRB的最小值作为φ1和φ2的孤立函数可能涉及非常复杂的解法。 
通过使用CRB近似值可以使该解法易于处理。在该近似值中,将变量φ1和φ2与所有其他变量解耦。此近似值利用以下知识: 
方程式3中的噪声协方差 
Figure DEST_PATH_GDA00002455764000127
为带式矩阵。可以使分块对角近似值为: 
Figure DEST_PATH_GDA00002455764000128
如果信道h01平坦衰落,则方程式11将提供精确等值。 
因此,方程式7的CRB可以简化为以下方程式的近似CRB: 
Figure DEST_PATH_GDA00002455764000129
Figure DEST_PATH_GDA000024557640001210
λ(φ21)为CFO估计误差的均方误差(MSE),并且为由以下方程式所给出的非负函数: 
其中, g 1 = [ 1 , ρ 1 , ρ 1 2 , · · · , ρ 1 N CFO - 1 ] T , g 2 = [ 1 , ρ 2 , ρ 2 2 , · · · , ρ 2 N CFO - 1 ] T , ρ 1 = e j φ 1 并且  ρ 2 = e jφ 2
现参照图7,图7示出当CFOφ21改变时,不同数目NCFO个CFO估计区块的最小CFO估计误差,即λ(φ21)的值的图表。对于任何NCFO,当差值为-π或π时,函数λ(φ21)偏小。这可以表明,稳健选择是设计前序信号以使得此差值接近-π或π。 
在方程式12中可以看出,通过将方程式11应用于方程式7,现将区块旋转角度φ1和φ2与 
Figure DEST_PATH_GDA00002455764000136
和 
Figure DEST_PATH_GDA00002455764000137
分离。这允许独立于 和 而最小化近似CRB,并且可以获得使非负函数λ(φ21)最小化的特定值φ21。 
如图7中所示,对于任何NCFO,当φ21=0时,λ(φ21)在其最大值处,并且方程式12中的有效CRB成为: 
Figure DEST_PATH_GDA000024557640001310
= ∞
因为在时间1中分别从两个源S1和S2传达的前序信号430和440在时间2中返回该源时不能被区分开,所以λ(φ21)在其最大值处。 
根据方程式5和方程式6,S1和S2的所感知的旋转角度之间的差异,即φ21,由以下方程式组组成: 
Figure DEST_PATH_GDA000024557640001312
由于实体CFO,即±(f1-f2)可以由于测距而减小,所以θ12的小值产生φ12的小值。因此,为了避免获得φ21=0,区块旋转角度θ1和θ2不应为θ21。换言之,S1和S2不应均发送周期性前序信号。在源节点S1和S2处发送周期性前序信号为最差状况回落条件。当不涉及近似值时,方程式15将适用于方程式7。这是因为当φ21=0时,矩阵 
Figure DEST_PATH_GDA000024557640001313
成为奇异的并且方程式7的CRB(φ2)成为无限的。 
还可以看出,方程式13中所给出的下界对应于当在时间1中S1不发送信号时的情况。这表明还原至单向中继的情况,其中两个源S1和S2轮流发射消息,各个轮回持续2个时槽,即在时间1和时间2中。在单向中继的情况中,中继节点110在时间2中中继其在时间1中接收的消息。 
另外,如果信道为平坦衰落,则下界将为点对点传输中有效CFOφ2的CRB。这意味着,S1和S2可以能够执行无干扰通信。由于 
Figure DEST_PATH_GDA00002455764000141
因此,当φ21=π并且NCFO为奇数时,可以获得下界。 
使用方程式16,当θ21=π时,确定φ21=π。作为测距的结果,可以忽略实际CFO f2-f1。因此,为了使常规前序信号的修改最小化,S1和S2的实体区块旋转角度可以为: 
θ1=0且θ2=π(18) 
在这种情况下,分配θ1和θ2,使得θ1和θ2在介于0与π之间的范围上均等地间隔,即θ1和θ2间隔开π弧度,并且可以达到CRB下界的估计性能。应注意的是,对于中继节点的数目和各个中继节点中的天线的数目不强加限制。 
应注意的是,尽管上文使用S1作为示例来描述CFO的CRB和近似CRB,但是应理解的是,对于S2,可以类似地获得CFO的CRB和近似CRB。 
现返回图6,在步骤530中,在时域中,源节点S1和S2的起始前序信号230分别旋转了区块旋转角度θ1和θ2,并且受制于频谱规则以形成经旋转的前序信号430和440。所应用的频谱规则包括将频域信号的次载波的功率限定在频谱屏蔽以下。应注意的是,当将频率移位应用至相应的频域信号时,可以了解时域信号区块的旋转。 
使用一具体示例,在两个源节点S1与S2处,所提供的起始前序信号230具有IEEE 802.11a/g/n设计。将称为短前序信号的用于CFO估计的周期性训练序列用作起始前序信号230,并且每个该周期性训练序列均被构建为10个16取样的区块。换言之,NBLK=10并且L=16。通信系统100使用正交频分复用(OFDM),其中各个符号包括64个取样。 
如方程式18中取区块旋转角度θ1和θ2。在这种情况下,可以说,由单位振幅因子来比例缩放起始前序信号230,并且区块旋转后的前序信号430类似于起始前序信号230。在源节点S2处,使用区块旋转角度θ2=π并且将起始前序信号230旋转θ2=π。 
作为一示例,随后使PAPR最小化,使得前序信号满足IEEE802.11a/g频谱屏蔽。应注意的是,可以存在大量可以满足频谱屏蔽的备选信号,并且 具有最低PAPR的信号如方程式19和方程式20中所给出的。通过以下方程式组分别给出在S1和S2处的优化的前序信号: 
Figure DEST_PATH_GDA00002455764000151
及 
Figure DEST_PATH_GDA00002455764000152
其中, 是确保 
Figure DEST_PATH_GDA00002455764000154
与 
Figure DEST_PATH_GDA00002455764000155
均消耗相同的平均功率的正规化因子。 
具体而言,在本示例中,可以如以下所述来优化前序信号。IEEE802.11a/g频谱屏蔽阻挡对应于64离散傅立叶变换(DFT)的总数为64的频率频段外的频率频段{0,27,28,29,...,35,36,37}。应注意的是,用于IEEE 802.11a/g的频谱屏蔽比用于IEEE 802.11n的频谱屏蔽更严格。 
对于S2而言,如下文将在方程式26中所示,在使用离散傅立叶变换(DFT)将 
Figure DEST_PATH_GDA00002455764000156
转换为时域之后,作为选择区块旋转的结果,仅占用16个频率位置{2,6,10,14,18,22,26,30,34,38,42,46,50,54,58,62}。在位置{30,34}处存在与频谱屏蔽的重叠并且消除这些重叠位置。因此,能够保留在{2,6,10,14,18,22,26,38,42,46,50,54,58,62}处覆盖16个频率频段中的14个的子集。 
在步骤532中,在源节点S1和S2中的每个处,确定每个次载波的调制集。这涉及确定可以用于调制的群集。在本示例中,与在IEEE 802.11a/g/n中由短训练序列所使用的调制集相同,两个群集以 
Figure DEST_PATH_GDA00002455764000157
的形式可用。因此,调制集的大小为Nmod=2。 
在步骤534中,在源节点S1和S2中的每个处,由频域中的Nmod个可能群集中的每一个群集来调制Nsub个次载波中的每个次载波。S1中的负载频域信号反映于方程式19中,并且当负载12个次载波时,Nsub=12。因此,在S1中将存在 
Figure DEST_PATH_GDA00002455764000161
个可能群集。 
在S2中,所负载的频域信号反映于方程式20中,并且当负载14个次载波时,Nsub=14。因此,可以用于排列的自由度导致总数为214-1个的可能群集(而非214个,因为如果在每一取样中可以通过正负号转换来由一个取样设计产生另一个取样设计,则64取样设计的两个集合实质上相等)。 
当在S1和S2中调制Nsub个次载波中的每个次载波时,将尚未消除的负载频域信号(即 
Figure DEST_PATH_GDA00002455764000162
或者 
Figure DEST_PATH_GDA00002455764000163
其中k=0,1,2,...,63)转换为时域。这可以通过在 
Figure DEST_PATH_GDA00002455764000164
或者 上执行64点DFT来完成。在将频域信号转换为时域时还使用上取样和内插法,以便俘获发生在取样之间的时域信号值。对于每一个完成的调制,计算峰值平均功率比(PAPR)。 
在步骤536中,在各个源节点S1和S2处,将 
Figure DEST_PATH_GDA00002455764000166
个组合中具有最低PAPR的一个组合选择为最佳调制。通过这样做,使传输的PAPR最小化,因此所选择的组合产生优化调制。在源节点S1处,获得的具有优化调制的时域信号为 
Figure DEST_PATH_GDA00002455764000167
形成各个L取样区块。在本示例中,时域信号具有L=16取样的周期,并且可以使用以下分配对 
Figure DEST_PATH_GDA00002455764000168
分组以形成L=16取样区块: 
x ~ 1 , n = ( e jπ ) n - m 16 x ~ 1 , m , n = 0,1,2 , . . . , 63 , m = mod ( n , 16 ) - - - ( 21 )
这导致 
DFT ( x ~ 1,0 x ~ 1,1 x ~ 1,2 · · · x ~ 1,63 - 1 ) = DFT ( ( - 1 ) 0 x ~ 1 [ BLK ] ( - 1 ) 1 x ~ 1 [ BLK ] ( - 1 ) 2 x ~ 1 [ BLK ] ( - 1 ) 3 x ~ 1 [ BLK ] ) , - - - ( 22 )
其中 
x ~ 1 [ BLK ] = x ~ 1,0 x ~ 1,1 x ~ 1,2 · · · x ~ 1 , L - 1 , L = 16 , - - - ( 23 )
则在S1处的前序信号430的集合410可以表达为: 
Figure DEST_PATH_GDA00002455764000172
此处,NBLK为4并且P1为表示功率的变量,在S1处将以该功率来发射前序信号430。接着,对于 
Figure DEST_PATH_GDA00002455764000173
其中N为DFT大小并且L为各区块的长度,通过以角度θk=π/2k旋转紧接先前的区块来产生第一前序信号区块之后的NBLK个前序信号430区块中的每一个。这等同于将各个紧接先前的区块乘以 
Figure DEST_PATH_GDA00002455764000174
因此,以方程式24为例,由于N=64,L=16,则 
类似地,在源节点S2处,在调制之后,从 
Figure DEST_PATH_GDA00002455764000176
的64点DFT获得优化时域信号 
Figure DEST_PATH_GDA00002455764000177
再一次,用于将取样分组为大小L=16的区块的分配为: 
x ~ 2 , n = ( e jπ ) n - m 16 x ~ 2 , m , n = 0,1,2 , . . . , 63 , m = mod ( n , 16 ) - - - ( 25 )
这导致 
DFT ( x ~ 2 , 0 x ~ 2 , 1 x ~ 2,2 · · · x ~ 2 , 63 - 1 ) = DFT ( ( - 1 ) 0 x ~ 2 [ BLK ] ( - 1 ) 1 x ~ 2 [ BLK ] ( - 1 ) 2 x ~ 2 [ BLK ] ( - 1 ) 3 x ~ 2 [ BLK ] ) , - - - ( 26 )
其中 
x ~ 2 [ BLK ] = x ~ 2 , 0 x ~ 2 , 1 x ~ 2 , 2 · · · x ~ 2 , L - 1 , L = 16 , - - - ( 27 )
则在S2处的前序信号440的集合420可以表达为: 
x ~ 2 = 64 52 P 2 ( - 1 ) 0 x ~ 2 [ BLK ] ( - 1 ) 1 x ~ 2 [ BLK ] ( - 1 ) 2 x ~ 2 [ BLK ] · · · ( - 1 ) N BLK - 1 x ~ 2 [ BLK ] , - - - ( 28 )
P2为表示功率的变量,在S2处将以该功率来发射前序信号430。接着,对于 
Figure DEST_PATH_GDA00002455764000183
其中N为DFT大小并且L为各区块的长度,通过以角度θk=π/2k旋转紧接先前的区块来产生第一前序信号区块之后的NBLK个前序信号440区块中的每一个。这等同于将各紧接先前的区块乘以 
Figure DEST_PATH_GDA00002455764000184
注意,如果来自方程式24的 
Figure DEST_PATH_GDA00002455764000185
的第一区块取为 
Figure DEST_PATH_GDA00002455764000186
则破坏频谱屏蔽。因此,为 
Figure DEST_PATH_GDA00002455764000187
使用不同于 
Figure DEST_PATH_GDA00002455764000188
的个别设计。 
应注意的是,尽管已描述步骤520至步骤536,使得S1和S2并行地执行其处理,但是可选地,S1和S2可以以非并行的方式执行步骤520至步骤536。另外,S1和S2可以逐个操作。 
在本示例中,在S1中具有根据方程式19的设计的负载频域信号具有2.24dB的PAPR。应注意的是,方程式19的设计类似于IEEE 802.11a/g/n中用于常规前序信号的设计。在S2中具有方程式20的设计的负载频域信号具有2.20dB的PAPR。因此,方程式20的前序信号设计可以具有较低PAPR的优点。这也可以说明,存在充分自由度来支持PAPR优化。 
另外,尽管已使用最低PAPR作为优化准则来描述步骤534和步骤536,但是设想,可以采用利用可用自由度的一些其他优化准则,例如,时域信号的自相关和/或交叉相关的最小化。 
参照图8和图9,图8图示 
Figure DEST_PATH_GDA00002455764000189
的时域波形,而图9图示 
Figure DEST_PATH_GDA000024557640001810
的频域波形。两个波形均为在区块旋转角度为θ2=π的情况下的源节点S2的波形。在图8 中, 
Figure DEST_PATH_GDA00002455764000191
的时域波形与在类似条件下但是在区块旋转内产生的对应时域波形形成对照。在图9中,可以看出,由于 
Figure DEST_PATH_GDA00002455764000192
的波形落入频谱屏蔽内,所以 
Figure DEST_PATH_GDA00002455764000193
的波形与频谱屏蔽要求达到完全一致。 
返回图6,在步骤540中,在S1处的前序信号430的集合410由S1发射。同样地,在S2处的前序信号440的集合420由S2发射。同时执行来自S1和S2的发射。 
线性滤波(步骤560) 
来自S1和S2的分别包括前序信号430的集合410和前序信号的440集合420的时域信号的发射发生在时间1中。在时间2中,在S1与S2处均自中继节点110接收信号 
Figure DEST_PATH_GDA00002455764000194
在步骤560中,在接收到的信号 上执行线性滤波。双向中继通信系统与点对点传输系统之间的差异在于:与后者系统中的一个频率音调不同,在前者系统中处理两个频率音调。如从方程式3可以看出,接收到的信号 
Figure DEST_PATH_GDA00002455764000196
包括两个分别以 
Figure DEST_PATH_GDA00002455764000197
和 
Figure DEST_PATH_GDA00002455764000198
表示的频率音调。 
在各个源节点S1和S2处,一个频率为事先已知,并且此已知频率可以使用定制滤波器来去除。因此,在步骤570中执行CFO估计之前,此滤波器可以执行自干扰减轻。以源节点S1为例,如方程式29中所定义的简单逐块滤波器 
Figure DEST_PATH_GDA00002455764000199
可以用以从由方程式3、方程式4和方程式5所定义的接收到的信号 
Figure DEST_PATH_GDA000024557640001910
中去除已知频率分量φ1。 
Figure DEST_PATH_GDA000024557640001911
由于 
Figure DEST_PATH_GDA000024557640001912
所以,可以将滤波器输出给出为 
Figure DEST_PATH_GDA000024557640001913
Figure DEST_PATH_GDA00002455764000201
应注意的是,滤波可以具有再成形并且进一步着色被感知的噪声频谱的潜能。这可以导致CFO估计性能的损失。然而,如下文将示出,通过检查来自滤波信号的包括CFO估计的分量(即,包括φ2的分量)的克拉美-罗界(CRB)是否大于来自滤波前的信号的分量的CRB,估计CFO的方法500可以不经受估计性能的损失。使用方程式7可以将φ2的CRB计算为: 
Figure DEST_PATH_GDA00002455764000202
其中 
Figure DEST_PATH_GDA00002455764000203
使用方程式8和方程式30可以证明 
Figure DEST_PATH_GDA00002455764000204
因此,使用方程式7,可以证明 由于φ2的CRB,即 
Figure DEST_PATH_GDA00002455764000206
不大于来自滤波前信号的CRB,即CRB(φ2),所以可以说,逐块线性滤波器 
Figure DEST_PATH_GDA00002455764000207
去除已知频率,即φ1而不包括CRB并且因此在CFO估计中不产生损失。 
估计CFO频率(步骤570) 
在步骤570中,根据线性滤波的信号来估计CFO的频率。在将滤波应用于接收到的信号后,在经滤波的信号中可以仅留下一个频率音调,并且因此可以使用本领域技术人员所已知的用于点对点传输中的CFO估计的任何技术来执行估计CFO,例如,使用基于最大似然(ML)的估计器。 
在存在有色噪声的情况下,现在可以将CFO估计的问题描述为估计单个音调。应注意的是,滤波器 
Figure DEST_PATH_GDA00002455764000208
可以具有无自干扰的优点。因此,可以使用基本相关器。 
φ2,EST表示CFO估计,并且 
Figure DEST_PATH_GDA00002455764000211
代表方程式31的滤波输出。使用基本相关器的可以具有非常简单的优点。 
模拟结果 
使用由产生和发射前序信号的方法510所获得之前序信号430和440来实施模拟。这些模拟使用线性滤波560以去除已知频率。在模拟中,所有信道分接点经历独立的瑞雷衰落,其中,其量值由指数幂延迟分布 
Figure DEST_PATH_GDA00002455764000212
调整,其中n为分接点索引,并且τrms为均方根延迟扩展。在模拟中使用描述延迟的扩展上升度的信道参数的三个集合: 
情况1:Lh=8,其中τrms=1; 
情况2:Lh=16,其中τrms=5;以及 
情况3:Lh=16,其中τrms=∞。 
情况3的信道参数对应于均匀幂延迟分布模型并且作为最差延迟扩展情况以进行比较。 
为简便起见,假设两个源节点S1和S2以相等的功率传达,即P=P1=P2,并且在所有接收机处的AWGN的方差相等, 
Figure DEST_PATH_GDA00002455764000213
即 
Figure DEST_PATH_GDA00002455764000214
信噪比定义为 
Figure DEST_PATH_GDA00002455764000215
设定在中继处所应用的缩放比例因子α以使在时间2中的中继处的总发射功率P0保持与在时间1中的总发射功率相等即P0=P1+P2=2P。后者允许与最大发射功率要求相一致,该最大发射功率要求限制对其他共信道用户所产生的干扰。 
图10为图示随着SNR改变,并且在f1=0.001,f2=-0.002以及NCFO=5的情况下,并且在应用情况1的情况下的CFO估计中的MSE的图表。图11为图示随着SNR改变,并且在f1=0.001,f2=-0.002以及NCFO=5的情况下,并且在应用情况2的情况下的CFO估计中的MSE的图表。图12为图示随着SNR改变,并且在f1=0.001,f2=-0.002以及NCFO=5的情况下,并且在应用情况3的情况下的CFO估计中的MSE的图表。在这三个图表中,所使用的估计器为基本相关器(即,诸如方程式34的相关器)或者基于最大似然(ML)的估计器。在所有情况下,使用估计CFO的方法500。还呈现图示自S1和S2角度看的克拉美--罗界性能的曲线,正如图示当使用周期性前序信号时的克拉美-罗界性能的曲线。 
如从曲线可以看出的,在所有三个图表中使用周期性前序信号的CRB性能执行得最差,并且通过以前序信号430、440来替换周期性前序信号而使MSE减少了大于28倍。还评估由方程式12所导出的近似CRB,但是该近似CRB未显示在图10、图11和图12中。对于最不分散信道(即,具有情况1的条件的信道)而言,在近似CRB情况下的性能与在用方程式7所计算的精确CRB情况下的性能之间存在小于1%的差异,而对于最分散信道(即,具有情况3的条件的信道)而言,在近似CRB情况下的性能与在用方程式7所计算的精确CRB情况下的性能之间存在稍微大于1%的差异。因此,可以说,在前序信号430、440设计中使用近似CRB是适当的。另外,可以观察到在S1和S2处的MSE性能几乎相同。这可能是因为在S1和S2处以相等的功率发射,并且还因为所有信道在统计上均相等,这相应地保持了系统对称。 
在所使用的估计器之间进行比较,图10、图11和图12图示出,在高SNR值处基本相关器能够产生约比CRB的MSE高1.25倍的较低的MSE。另外,可以看出,τrms值越小,基本相关器的MSE越接近CRB性能的MSE。具体而言,相对于其中τrms=∞的图12而言,基本相关器的MSE比CRB性能的MSE高1.25倍。然而,在其中τrms=1的图10中,基本相关器的MSE比CRB性能的MSE高1.1倍。值得注意的是,尽管基本相关器结构简单并且其不具有信道信息,但是其仍然可以产生接近CRB性能的性能。 
另外,在图10、图11和图12中,在低SNR水平处,通过在所有三种情况中基本相关器均产生较低的MSE来将基本相关器图示为比基于ML的方案优越。对于图10、图11和图12而言,高SNR水平处的曲线之间的狭窄性能间隙表明,使用方法510所产生并且发射的前序信号对在确定CFO估计中所使用的估计器的选择可以具有稳健性。 
在本说明书中,交替使用术语“前序信号”和“前序信号区块”来代表前序信号430和/或440。交替使用术语“源”和“源节点”来代表通信系统100的源节点120、122。 
尽管上文已参照双向中继情况(其中,仅存在两个源节点)描述了方法500,但是设想,可以将方法500应用于涉及K个(大于一个)源的多向中继情况。在这种情况下,用起始前序信号230的预定集合来配备各个源。分配应用至K个起始前序信号230集合的区块旋转角度θ1、θ2、……、θK,使得这些区块旋转角度在介于0与π的范围上均等间隔。 
尽管已详细描述本发明的示例性实施例,但是本领域技术人员读者将清楚,在不脱离本发明的范围的情况下,许多变型是可能的。 
另外,本领域技术人员将理解,可以将源和/或中继节点实现为诸如移动电话的移动装置和/或诸如基地站的固定装置。同样地,设想,可以将源和/或中继节点实现为集成电路或者芯片上系统方案。 

Claims (47)

1.一种用于在通信装置处估计载波频率偏移的方法,所述方法包括:
产生信号区块;
以多个不同旋转角度旋转所述区块以形成对应的第一前序信号集合;
将所述第一前序信号集合作为时域信号发射至中继;
接收来自所述中继的经重传的时域信号,所述经重传的时域信号为所述第一前序信号集合与来自另一通信装置的第二前序信号集合的组合;以及
基于所述接收到的经重传的时域信号来估计所述信道频率偏移。
2.根据权利要求1所述的用于估计载波频率偏移的方法,进一步包括:
将优化调制应用至所述第一前序信号集合以形成所述时域信号。
3.根据权利要求2所述的用于估计载波频率偏移的方法,其中所述信号区块的各个信号对应于用于发射所述第一前序信号集合的多个次载波中的各相应的次载波。
4.根据权利要求3所述的用于估计载波频率偏移的方法,其中应用所述优化调制的步骤包括:
排列多个组合,每个组合将所述旋转信号区块的布置与多个群集相关联;
用所述多个次载波中的各相应次载波来调制所述多个组合中的每个组合;以及
从所述经调制的多个组合中选择出具有最佳组合的被选信号,所述最佳组合使所述对应信号的峰值平均功率比最小化;其中所述被选信号为所述时域信号。
5.根据前述任一项权利要求所述的用于估计载波频率偏移的方法,其中旋转所述区块的步骤包括:
确定由所述多个次载波所使用的多个调制符号以形成所述第一前序信号集合中的每一个。
6.根据权利要求4或5所述的用于估计载波频率偏移的方法,其中在所述通信装置处旋转所产生的信号区块的步骤进一步包括:
通过应用所述多个不同旋转角度中的一个旋转角度,来移位所述信号中的每个的频率。
7.根据权利要求6所述的用于估计载波频率偏移的方法,其中在所述通信装置处旋转所产生的信号区块的步骤包括:比例缩放所述信号中的每个的振幅。
8.根据前述任一项权利要求所述的用于估计载波频率偏移的方法,其中所述多个不同旋转角度中的一个旋转角度由所述多个次载波的数目以及所述时域信号的预定长度获得。
9.根据前述任一项权利要求所述的用于估计载波频率偏移的方法,其中所述多个不同旋转角度中的一个旋转角度由先前信号区块的角度获得。
10.根据前述任一项权利要求所述的用于估计载波频率偏移的方法,其中所述多个不同旋转角度中的一个旋转角度是使用所述载波频率偏移估计的克拉美-罗界而确定的。
11.根据权利要求10所述的用于估计载波频率偏移的方法,其中所述克拉美-罗界为近似值。
12.根据权利要求10或11所述的用于估计载波频率偏移的方法,其中所述多个不同旋转角度中的一个旋转角度具有介于0与π之间并且包括0及π的值。
13.根据前述任一项权利要求所述的用于估计载波频率偏移的方法,其中所述多个不同旋转角度中的一个旋转角度是根据其他通信装置的数目而获得的。
14.根据前述任一项权利要求所述的用于估计载波频率偏移的方法,其中所产生的区块为IEEE 802.11前序信号。
15.根据前述任一项权利要求所述的用于估计载波频率偏移的方法,其中所述第一前序信号集合中的一个前序信号是由旋转先前旋转信号区块而获得。
16.根据前述任一项权利要求所述的用于估计载波频率偏移的方法,其中所述时域信号使用正交频分复用来发射。
17.根据前述任一项权利要求所述的用于估计载波频率偏移的方法,其中所述时域信号为非周期性的。
18.根据前述任一项权利要求所述的用于估计载波频率偏移的方法,其中所述经重传的时域信号包括由所述中继从另一装置重传的训练信号。
19.根据前述任一项权利要求所述的用于估计载波频率偏移的方法,其中所述载波频率偏移介于所述时域信号与所述接收到的经重传的时域信号之间。
20.根据前述任一项权利要求所述的用于估计载波频率偏移的方法,其中估计所述载波频率偏移的步骤包括:对所述接收到的经重传的时域信号进行线性滤波。
21.根据权利要求20所述的用于估计载波频率偏移的方法,其中估计所述载波频率偏移的步骤进一步包括:在所述线性滤波信号上执行相关。
22.一种用于在第一通信装置处估计载波频率偏移的中继方法,所述方法包括:
在中继处,接收从所述第一通信装置发射的第一时域信号,以及从第二通信装置发射的第二时域信号;
其中所述第一时域信号包括第一前序信号集合,通过以对应的多个不同的第一旋转角度旋转第一产生的信号区块来形成所述第一前序信号集合;并且
其中所述第二时域信号包括第二前序信号集合,通过以对应的多个不同的第二旋转角度旋转第二产生的信号区块来形成所述第二前序信号集合;以及
从所述中继将经重传的时域信号重传至所述第一通信装置,以允许所述第一通信装置基于所述经重传的时域信号来估计所述载波频率偏移,所述经重传的时域信号为所述第一前序信号集合与所述第二前序信号集合的组合。
23.一种用于估计通信系统中的载波频率偏移的方法,所述通信系统包括:第一通信装置、第二通信装置以及中继,所述方法包括:
在所述第一通信装置处,以多个不同的第一旋转角度旋转第一产生的信号区块,以形成对应的第一前序信号集合;
在所述第二通信装置处,以多个不同的第二旋转角度旋转第二产生的信号区块,以形成对应的第二前序信号集合;
从所述第一通信装置和所述第二通信装置中的每个通信装置将所述各相应的第一前序信号集合和第二前序信号集合作为时域信号发射至所述中继;
在所述第一通信装置处,接收来自所述中继的经重传的时域信号,所述经重传的信号为所述第一前序信号集合与所述第二前序信号集合的组合;以及
基于所述接收到的经重传的时域信号来估计所述信道频率偏移。
24.根据权利要求23所述的用于估计载波频率偏移的方法,其中所述多个不同的第一旋转角度的起始角度与所述多个不同的第二旋转角度的起始角度相差π。
25.根据权利要求24所述的用于估计载波频率偏移的方法,其中所述多个不同的第一旋转角度的所述起始角度为0。
26.一种通信装置,包括:
处理器,其被配置以产生信号区块并且以多个不同的旋转角度旋转所述区块,以形成对应的第一前序信号集合;以及
发射机,其被配置以将所述第一前序信号集合作为时域信号发射至中继;
接收机,其被配置以接收来自所述中继的经重传的时域信号,所述经重传的时域信号为所述第一前序信号集合与来自另一通信装置的第二前序信号集合的组合;并且
其中所述处理器进一步被配置以基于所述接收到的经重传的时域信号来估计所述信道频率偏移。
27.根据权利要求26所述的通信装置,其中所述发射机进一步被配置以:
将优化调制应用至所述第一前序信号集合以形成所述时域信号。
28.根据权利要求27所述的通信装置,其中所述信号区块的各信号对应于用于发射所述第一前序信号集合的多个次载波中的各相应的次载波。
29.根据权利要求28所述的通信装置,其中所述发射机进一步被配置以:
排列多个组合,每个组合将所述旋转信号区块的布置与多个群集相关联;
用所述多个次载波中的各相应次载波来调制所述多个组合中的每个组合;以及
从经调制的多个组合中选择出具有最佳组合的被选信号,所述最佳组合使所述对应信号的所述峰值平均功率比最小化;其中所述被选信号为所述时域信号。
30.根据权利要求26至权利要求29中任一项所述的通信装置,其中所述处理器进一步被配置以:
确定由所述多个次载波所使用的多个调制符号以形成各个所述第一前序信号集合。
31.根据权利要求29或30所述的通信装置,其中所述处理器进一步被配置以:
通过应用所述多个不同旋转角度中的一个旋转角度,来移位所述信号中的每个信号的频率。
32.根据权利要求31所述的通信装置,其中所述处理器进一步被配置以:
比例缩放所述信号中的每个信号的振幅。
33.根据权利要求26至权利要求32中任一项所述的通信装置,其中所述多个不同旋转角度中的一个旋转角度由所述多个次载波的数目以及所述时域信号的预定长度获得。
34.根据权利要求26至权利要求33中任一项所述的通信装置,其中所述多个不同旋转角度中的一个旋转角度由先前信号区块的角度获得。
35.根据权利要求26至权利要求34中任一项所述的通信装置,其中所述多个不同旋转角度中的一个旋转角度使用所述载波频率偏移估计的克拉美-罗界来确定。
36.根据权利要求35所述的通信装置,其中所述克拉美-罗界为近似值。
37.根据权利要求36或36所述的通信装置,其中所述多个不同旋转角度中的一个旋转角度具有:介于0与π之间并且包括0及π的值。
38.根据权利要求26至权利要求37中任一项所述的通信装置,其中所述多个不同旋转角度中的一个旋转角度是根据其他通信装置的数目而获得的。
39.根据权利要求26至权利要求38中任一项所述的通信装置,其中所产生的区块为IEEE 802.11前序信号。
40.根据权利要求26至权利要求39中任一项所述的通信装置,其中所述第一前序信号集合中的前序信号是由旋转先前旋转信号区块而获得的。
41.根据权利要求26至权利要求40中任一项所述的通信装置,其中所述时域信号使用正交频分复用来发射。
42.根据权利要求26至权利要求41中任一项所述的通信装置,其中所述时域信号为非周期性的。
43.根据权利要求26至权利要求42中任一项所述的通信装置,其中所述经重传的时域信号包括由所述中继从另一装置重传的训练信号。
44.根据权利要求26至权利要求43中任一项所述的通信装置,其中所述载波频率偏移介于所述时域信号与所述接收到的经重传的时域信号之间。
45.根据权利要求26至权利要求44中任一项所述的通信装置,其中所述处理器进一步被配置以对所述接收到的经重传的时域信号进行线性滤波。
46.根据权利要求45所述的通信装置,其中所述处理器进一步被配置以在所述线性滤波信号上执行相关。
47.一种用于通信装置的集成电路,包括:
处理单元,其被配置以产生信号区块并且以多个不同旋转角度旋转所述区块,以形成对应的第一前序信号集合;
接口,其被配置以将所述第一前序信号集合作为时域信号发射至中继,并且进一步被配置以接收来自所述中继的经重传的时域信号,所述经重传的时域信号为所述第一前序信号集合与来自另一通信装置的第二前序信号集合的组合;并且
其中所述处理单元进一步被配置以基于所述接收到的经重传的时域信号来估计所述信道频率偏移。
CN201180014043XA 2010-01-26 2011-01-26 用于估计载波频率偏移的方法与装置 Pending CN102986292A (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
SG201000591 2010-01-26
SG201000591-6 2010-01-26
PCT/SG2011/000036 WO2011093798A1 (en) 2010-01-26 2011-01-26 Method and device for estimating carrier frequency offset

Publications (1)

Publication Number Publication Date
CN102986292A true CN102986292A (zh) 2013-03-20

Family

ID=44319599

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201180014043XA Pending CN102986292A (zh) 2010-01-26 2011-01-26 用于估计载波频率偏移的方法与装置

Country Status (5)

Country Link
US (1) US20120300644A1 (zh)
CN (1) CN102986292A (zh)
SG (1) SG182720A1 (zh)
TW (1) TW201203897A (zh)
WO (1) WO2011093798A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103986535A (zh) * 2014-05-29 2014-08-13 国家电网公司 一种测试接收机可接受频率偏移的装置

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013013697A1 (en) * 2011-07-22 2013-01-31 Alcatel Lucent Machine type communications in a radio network
CN104066147A (zh) * 2013-03-19 2014-09-24 中兴通讯股份有限公司 基于下行探测参考信号搜索网络节点的方法、装置及设备
CN103701733B (zh) * 2013-09-28 2017-03-01 河北工业大学 一种td‑lte中继系统频偏估计的方法
US10560302B2 (en) 2017-08-28 2020-02-11 Indian Institute of Technology Kharagpur Method and system for joint training sequences design for correlated channel and frequency offsets estimation

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008069555A1 (en) * 2006-12-05 2008-06-12 Electronics And Telecommunications Research Institute Apparatus and method for channel estimation and synchronization in ofdm/ofdma relay system
US20090041095A1 (en) * 2007-08-10 2009-02-12 Samsung Electronics Co., Ltd. Method and apparatus for fine frequency synchronization in wireless broadband (wibro) system without using global positioning system (gps) receiver
CN101483622A (zh) * 2008-01-11 2009-07-15 株式会社Ntt都科摩 在双向中继网络中信道估计的方法、设备和系统
CN101515917A (zh) * 2009-03-25 2009-08-26 东南大学 基于双向中继的多用户无线通信系统及方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8958493B2 (en) * 2004-03-31 2015-02-17 Infineon Technologies Ag Operation for backward-compatible transmission
EP1937006A1 (en) * 2006-12-22 2008-06-25 Siemens Networks GmbH & Co. KG Multi-antenna relay station with two-way channel
US8208522B2 (en) * 2008-03-07 2012-06-26 Nokia Corporation System and methods for receiving OFDM symbols having timing and frequency offsets
CN101800616B (zh) * 2009-02-10 2012-11-21 富士通株式会社 数据中继装置、通信装置和方法
US8488539B2 (en) * 2009-07-16 2013-07-16 Ralink Technology Corp. Method of generating preamble sequence

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008069555A1 (en) * 2006-12-05 2008-06-12 Electronics And Telecommunications Research Institute Apparatus and method for channel estimation and synchronization in ofdm/ofdma relay system
US20090041095A1 (en) * 2007-08-10 2009-02-12 Samsung Electronics Co., Ltd. Method and apparatus for fine frequency synchronization in wireless broadband (wibro) system without using global positioning system (gps) receiver
CN101483622A (zh) * 2008-01-11 2009-07-15 株式会社Ntt都科摩 在双向中继网络中信道估计的方法、设备和系统
CN101515917A (zh) * 2009-03-25 2009-08-26 东南大学 基于双向中继的多用户无线通信系统及方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
GONGPU WANG ET AL.: "Joint Frequency Offset and Channel Estimation Methods for Two-Way Relay Networks", 《IEEE GLOBECOM 2009》, 31 December 2009 (2009-12-31) *
L. B. THIAGARAJAN ET AL.: "CARRIER FREQUENCY OFFSET AND CHANNEL ESTIMATION IN SPACE-TIME NON-REGENERATIVE TWO-WAY RELAY NETWORK", 《2009 IEEE》, 31 December 2009 (2009-12-31) *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103986535A (zh) * 2014-05-29 2014-08-13 国家电网公司 一种测试接收机可接受频率偏移的装置

Also Published As

Publication number Publication date
TW201203897A (en) 2012-01-16
WO2011093798A1 (en) 2011-08-04
US20120300644A1 (en) 2012-11-29
SG182720A1 (en) 2012-08-30

Similar Documents

Publication Publication Date Title
CN104158785B (zh) 用于确定无线移动通信系统中用户设备的位置的方法
JP4057603B2 (ja) 送信方法、受信方法、送信装置、受信装置
RU2365055C2 (ru) Точное вхождение в синхронизм
CN102176707B (zh) 多载波系统中使用信道估计的码元定时校准
KR100707287B1 (ko) 무선 통신을 위한 편파 상태 기술들
CN100556017C (zh) 在多入多出正交频分复用系统中为同步而发射前导的方法
US7672221B2 (en) Radio receiver and radio signal receiving method
CN101689905B (zh) 在无线通信系统中用于确认和Category 0比特的发射分集
CN101194481B (zh) Ofdm系统中的导频发送方法和设备
CN101027864A (zh) 多天线无线通信系统的符号定时检测方法
CN101958865A (zh) 解调参考信号的生成方法及装置
CN102986292A (zh) 用于估计载波频率偏移的方法与装置
CN101075999B (zh) 室内ofdm系统toa训练符号构建及toa估计方法与装置
CN101223753A (zh) 用于精细ofdm符号同步的方法和同步器以及用于ofdm符号接收的方法/接收器
CN101447962B (zh) 同步信息的发送与接收方法及装置、同步系统
CN108737317B (zh) 广义混合载波频选信道传输方法
CN100562002C (zh) 一种正交频分复用调制系统中保护间隔的填充方法及其通信系统
CN108347399B (zh) 一种前导信号的生成和接收方法以及接收装置
CN104244398B (zh) 基于ofdm调制的微功率无线通信系统同步信号的生成方法
CN101291311A (zh) 多输入多输出正交频分复用系统的同步实现方法及装置
CN102404071B (zh) 同步训练符号块编码调制方法和装置
CN103546416A (zh) Ofdm符号同步方法
Zoltowski et al. Complementary codes based channel estimation for MIMO-OFDM systems
CN101789816B (zh) 一种ofdma多天线系统的符号检测方法及装置
CN102761505A (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
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20130320