技术领域
本发明涉及一种载波频偏估测方法,特别是涉及一种利用内插法求得两 段具有几乎重复特性的信号,再根据这两段信号所求的延迟相关值进行载波 频偏估测的方法。
背景技术
传统在时间轴上补偿载波频偏的
算法利用延迟相关的方法来估计此频偏 的大小,这种方法有个最基本的要求,就是传送的前导信号中必须至少有两 段是相同的信号,可能是两个或多个相同的OFDM symbol或者是一个OFDM symbol中由两段或更多段相同信号组成。有了重复的特性,就可以利用延迟 相关的方法来估算载波频偏,其基本想法可用下列数学式子来说明:
其中,z:延迟相关函数,表示累积N次相同信号
位置的取样点乘积。
rn:收到的信号在第n点的取样值。
rn+D:收到的信号在第n+D点的取样值。
sn:传送的信号在第n点的取样值。
sn+D:传送的信号在第n+D点的取样值。
fΔ:载波
频率偏移。
Ts:取样周期。
D:表示延迟的长度(在没有噪声与频偏下,隔了D点的取样值会是一样 的)。
*:表示取共轭复数。
由此相关函数可以发现该函数的
相位里面就包含我们所要的载波频偏, 只要取出此相关函数的相位,并做适当的normalized即可得到频偏大小.运 算如下:
在802.16e的规格里面,前导信号(preamble)在连续时间域上有3段重 复的特性,因此延迟相关的技术就可以被采用。但不同的是,这三段重复的 信号在经取样后取得的值并不具有如连续时间域上具有完美的重复特性。换 句话说重复信号的相对应位置点上的信号值会有一个
相位差,这导致了延迟 相关的结果会多出一个非理想的相位(undesired phase term),该相位会因 为不同的sector的前导信号而有所不同。请参考图1a与图1b。图1a与图 1b表示一个从振幅响应来看具有三段重复特性的信号,以不同的取样速率所 得到的离散时间域的数据值的示意图。在图1a中,若我们将这
模拟信号等间 距的取9点信号下来,则可以观察到取出来的9点信号中每段重复信号的相 对应位置的取样值是一样的。但在图1b中,若改取10点信号下来,则可以 发现在每个重复信号里面的相对应位置的信号值几乎不太一样。在802.16e 标准中的前导信号,由于我们采用的FFT size是1024,512以及128点,这 三个值都不是3的整数倍,会造成如图1b中,三段重复信号里面的取样点并 非完全相同。使得在进行延迟相关的运算时会有非理想的相位出现。就数学 上的式子来观察也可以得到如前述的结果:对于完美重复特性的信号进行延 迟相关的操作得出来的表示式如下:
而采用非完美重复特性的信号所做出来的延迟相关结果如下所示:
从表示式(1)与(2)就可以发现,不完美重复特性的信号会有Θn相位 出现,因此必须先将这个值消除后才能将求得更正确的
频率偏移值。
发明内容
本发明的目的是将原本前导信号的非完美重复
片段通过一个简单的内插
滤波器,重建出重复片段信号,再应用延迟相关(delay correlation)的运算 就可以得到相当准确的频偏估计。
本发明提供一种载波频偏估测方法,用以估测接收信号的载波频偏值, 包括接收前导信号,该前导信号包括第一段重复信号与第二段重复信号;对 该第一段重复信号进行内插取样程序以求得第一信号;根据该第一信号与该 第二段重复信号取得延迟相关值,并根据该延迟相关值估测该载波频偏值。
在另一
实施例中,该载波频偏估测方法更包括取得该延迟相关值的
角度 值;将该角度值乘上预定值;将该角度值与该预定值的积减去该接收信号的 相位旋转值以求得差值,并根据该差值求得该载波频偏值。
本发明更提供一种载波频偏估测系统,用以估测接收信号的载波频偏值, 该接收信号包括第一段重复信号与第二段重复信号,该系统包括共轭信号产 生单元、第一乘法器以及载波频偏估测器。该共轭信号产生单元,包括内插 滤波器与运算单元,接收该第一段重复信号,并通过该内插滤波器与该运算 单元以产生对应于该第二段重复信号的第一共轭信号。该第一乘法器,计算 该第二段重复信号与该第一共轭信号的第一乘积值。该载波频偏估测器,根 据该乘积值并累积一段时间后估测该载波频偏值。
在一实施例中,该载波频偏估测器更包括相位计算单元、第二乘法器以 及减法器。该相位计算单元根据该乘积值取得相位值。该第二乘法器,求得 该相位值与预定值的第二乘积值。该减法器,接收该第二乘积值并减去该接 收信号的相位旋转值以求得差值,并根据该差值求得该载波频偏值。
附图说明
图1a与图1b表示一个从振幅响应来看具有三段重复特性的信号,以不 同的取样速率所得到的离散时间域的数据值的示意图。
图2为可消除载波频偏估测时所产生的非理想相位的一种载波频偏估测 系统的功能方
块图。
图3为可消除载波频偏估测时所产生的非理想相位的另一种载波频偏估 测系统的功能方块图。
图4为根据本发明的另一种载波频偏估测系统的功能方块图。
图5为
输入信号r(n)经过图4中内插滤波器48的一实施例处理后所取 得的数据值的示意图。
图6为根据本发明的另一种载波频偏估测系统的功能方块图。
图7为输入信号r(n)经过图5中内插滤波器68的一实施例处理后所取 得的数据值的示意图。
附图符号说明
21、41、61-延迟单元
22、42、62-运算单元
23、43、63-算数平均单元
26、46、66-相位计算单元
28、47、67-相位消除单元
48、68-内插滤波器
24、27、29、44a、44b、44c、39a、39b、39c、64a、64b、64c-乘法 器
31-第一延迟单元
32-第二延迟单元
33a-第一运算单元
33b-第二运算单元
35a-第一算数平均单元
35b-第二算数平均单元
37a-第一相位计算单元
37b-第二相位计算单元
D1、D2、D3、D4、D5、D6、D7、D8、D9、D10、D41、D42、D51、D52、D61、D62、D71、 D72、D11、D12、D21、D22、D31、D32-数据值。
具体实施方式
图2为可消除载波频偏估测时所产生的非理想相位的一种载波频偏估测 系统的功能方块图。延迟单元21接收输入信号r(n),并延迟D个取样点, 其中D表示接收信号中重复出现的信号所占的长度。换言之,延迟单元21为 一个存储装置(移位寄存器一shift register或内存memory),输入信号会从 r(n)一直累积到r(n+D-1)的值,当第r(n+D)的值进入该延迟单元后,第r(n) 的值就会输出来。运算单元22接收延迟单元21输出的信号,并求得该信号 的共轭复数。乘法器24将信号r(n+D)与运算单元22输出的信号r(n)相乘 以得到信号x(n)。举例来说,若rn表示输入信号r(n)在第n点的取样值,则 rn+D则表示延迟D个取样点后的输入取样值,亦即为输入信号r(n)在第n+D 点的取样值。则信号x(n)在第n点的数据值xn即可表示如下:
xn=rnr*n+D
算数平均单元23接收多个信号x(n)的数据值,并根据所求得算数平均 数,产生延迟相关值(delay correlation value)z。载波频偏估测器25则 根据延迟相关值z估测载波频偏值在本实施例中,载波频偏估测器25更 包括相位计算单元26、乘法器27与29,以及相位消除单元28。相位计算单 元26接收并求得延迟相关值z的相位角∠z。在本发明的一实施例802.16e 系统,由于802.16e前导信号刚好有个特性,即三个sector的前导信号其延 迟相关出来的相位里面的非理想相位分别落在0度、及+/-120度附近,所以 在CFO estimate里面只需要乘上3就可以大致将此非理想相位消除。因此乘 法器27将延迟相关值z的相位角∠z乘3,并在相位消除单元28拿掉±2π项。 乘法器29则将相位消除单元28的
输出信号乘上补偿值以求得载波频 偏值其中D为输入信号中重复的长度,Ts为输入信号的取样频率。
延迟相关值z可由下列方式求得:
在802.16e标准中,位于不同的sector的前导信号其Ф会不同,分别是 0度+/-120度,所以延迟相关值z的相位角∠z=-2πfΔDTs+Ф可写成下列三个 式子:
∠z1=-2πfΔDTs
∠z2=-2πfΔDTs+120
∠z3=-2πfΔDTs-120
将上述三个式子分别乘上3后再相加,便可消除+/-360度的项。其中运 算单元22也可放在延迟单元的前面或后面,但补偿值须调整为
图3为可消除载波频偏估测时所产生的非理想相位的另一种载波频偏估 测系统的功能方块图。在图3中的载波频偏估测系统,将前导信号中的三段 重复信号都列入考虑,除计算第一段与第二段的第一延迟相关值,同时也计 算第一段与第三段的第二延迟相关值,其表示式如下:
z12表示第一延迟相关值,z13表示第二延迟相关值。在表示式中,Ф与Θ分 别表示第一延迟相关值与第二延迟相关值中因前导信号的三段重复信号的不 完美重复特性(|sn|≠|sn+D|≠|sn+2D|),所造成的相位偏移。且根据实验发现这两 个角度(Ф与Θ)刚好是互成反相,换句话说,Ф=-Θ。因此只要将这两种 延迟相关的相位取出后相加就可以得到载波频偏值。其数学表示式如下:
接着取出z的相位,并根据z的相位求得到载波频偏值如下所示:
第一延迟单元31接收输入信号r(n),并延迟D个取样点,其中,D表示 接收信号中重复出现的信号所占的长度。换言之,第一延迟单元31为用以延 迟三段重复的信号中的第一段重复信号。第一运算单元33a接收第一延迟单 元31输出的信号,并求得该信号的共轭复数。乘法器39a将信号r(n+D)与 第一运算单元33a输出的信号r(n)相乘以得到信号x12(n)。第一算数平均单 元35a接收多个信号x12(n)的数据值,并根据所求得算数平均数,产生延迟 相关值z12。第一相位计算单元37a接收并求得延迟相关值z的相位角∠z12。 第二延迟单元32接收输入信号r(n),并延迟2D个取样点,其中D表示接收 信号中重复出现的信号所占的长度。换言之,延迟单元32为用以延迟三段重 复的信号中的第一、二段重复信号。第二运算单元33b接收第二延迟单元32 输出的信号,并求得该信号的共轭复数。乘法器39b将信号r(n+2D)与第二 运算单元33b输出的信号r(n)相乘以得到信号x13(n)。第二算数平均单元35b 接收多个信号x13(n)的数据值,并根据所求得算数平均数,产生延迟相关值 z13。第二相位计算单元37b接收并求得延迟相关值z的相位角∠z13。加法器 38将相位角∠z12与∠z13相加以求得相位角∠z。在本发明的另一实施例中, 可先通过乘法器将第一延迟相关值与第二延迟相关值相乘,在通过相位计算 单元,如第一相位计算单元37a或第二相位计算单元37b,求得载波偏移值 的相位角∠z。乘法器39c则将载波偏移值的相位角∠z乘上补偿值 以求得载波频偏值其中,第一运算单元33a也可放在第一延迟单元31的 前面或后面,第二运算单元33b也可放在第一延迟单元32的前面或后面,但 补偿值须调整为
图4为根据本发明的另一种载波频偏估测系统的功能方块图。延迟单元 41接收输入信号r(n),并延迟D个取样点,其中,D表示接收信号中重复出 现的信号所占的长度。换言之,延迟单元41为用以延迟三段重复的信号中的 第一段重复信号。内插滤波器48是针对第一段重复信号做动作,并以内插法 的方式求得以三倍取样频率所求得的过度取样信号。
为更清楚说明内插滤波器48的运作方式,请参考图5。图5为输入信号 r(n)经过图4中内插滤波器48的一实施例处理后所取得的数据值的示意图。 因为要取得与第二段重复信号【D5,D6,D7】有更完美重复特征的第一信号, 我们在数据值D1至D4之间,以内插法的方式求得每个三等分点上的数据值, 即D11、D12、D21、D22等等。接着内插滤波器48会根据第二段重复信号【D5,D6, D7,】的取样点位置,自该等数据值中选择对应的数据值,以得到第一信号。 在本实施例中,第一信号由数据值D12、D22、D32所组成。
另外,滤波器的位置除了可以设计在图4所展示的位置,也可以将这个 装置放在延迟单元41之前,或放在另一路,如图6所示68的地方,(有多种 组合,延迟单元与运算单元可在同一或不同路,滤波器又可在延迟单元或运 算单元的前或后)图7是针对第二段重复信号做内插取样,说明内插滤波器 68的运作方式,请参考图7。图7为输入信号r(n)经过图6中内插滤波器68 的一实施例处理后所取得的数据值的示意图。因为要取得与第一段重复信号 【D1,D2,D3,D4】有更完美重复特征的第一信号,我们在数据值D4至D8之间, 以内插法的方式求得每个三等分点上的数据值,即D41、D42、D51、D52、D61、D62 等等。接着内插滤波器68会根据第一段重复信号【D1,D2,D3,D4】的取样点 位置,自该等数据值中选择对应的数据值,以得到第一信号。在本实施例中, 第一信号由数据值D41、D51、D61以及D71所组成。
运算单元42接收内插滤波器48输出的第一信号,并求得第一共轭信号, 该第一共轭信号由该第一信号中数据值的共轭复数值所组成。乘法器44a将 信号r(n+D)与运算单元42输出的信号相乘以得到信号x(n)。算数平均单 元43接收多个信号x(n)的数据值,并根据所求得算数平均数,并根据该算 数平均数产生延迟相关值z。载波频偏估测器45则根据延迟相关值z估测载 波频偏值在本实施例中,载波频偏估测器45更包括相位计算单元46、 乘法器44b与44c,以及相位消除单元47。相位计算单元46接收并求得延迟 相关值z的相位角∠z。由于802.16e前导信号刚好有个特性,即三个sector 的前导信号其延迟相关出来的相位里面的非理想相位分别落在0度、及 +/-120度附近,所以在CFO estimate里面只需要乘上3就可以大致将此非 理想相位消除。因此,乘法器44b将延迟相关值z的相位角∠z乘3,并在相 位消除单元47消去±2π。乘法器44c则将相位消除单元47的输出信号乘上 补偿值以求得载波频偏值其中,D为输入信号中重复的长度,Ts 为输入信号的取样频率。
虽然本发明已以具体实施例揭露如上,然其仅为了易于说明本发明的技 术内容,而并非将本发明狭义地限定于该实施例,任何所属技术领域中具有 通常知识者,在不脱离本发明的精神和范围内,当可作些许的更动与润饰, 因此本发明的保护范围当视本发明的
申请专利范围所界定者为准。