一种基于DFT的符号跳变单频复指数信号频率估计方法技术

本发明专利技术公开了一种基于DFT的符号跳变单频复指数信号频率估计方法，用于M

【技术实现步骤摘要】
一种基于DFT的符号跳变单频复指数信号频率估计方法


[0001]本专利技术涉及一种基于DFT的信号频率估计的方法，具体涉及一种用于 M
‑
QAM符号跳变单频复指数信号的频率估计算法。

技术介绍

[0002]随着高速无线通信容量需求的快速增长，频谱效率高的M元正交振幅调制 (M
‑
QAM)与相干接收机和数字信号处理相结合引起了人们广泛关注。由于频率偏移的影响严重降低了系统性能，因此，需要一个频率偏移估计器来保证系统的可靠性。在实际应用中，快速傅里叶变换(FFT)由于低复杂度，常被用作DFT的快速算法。当频率估计采用一定长度的FFT时，估计的精度会受到频率泄露和栅栏效应的影响。为获得更高精度，频域插值算法被广泛应用，现已有几种插值DFT(IpDFT)方法，但它们只能应用于不需要对其进行数据调制的纯正弦信号。然而，在许多应用中，未知数据被调制在信号上，如M
‑
QAM信号。在这种情况下，符号跳变是随机发生的。当在DFT周期发生符号跳变时，DFT系数可能会严重失真，导致传统的IpDFT方法产生较大估计误差。

技术实现思路

[0003]针对以上问题，本专利技术提出了一种基于DFT的用于符号跳变的单频复指数信号的频率估计方法。用方程给出了FFT单元与频率偏移的关系，使用三个不同的DFT单元来消除符号跳变的影响。
[0004]采样的M
‑
QAM信号s(n)序列可以描述为：
[0005]x(n)＝s(n)+q(n)
[0006]＝U
m exp(jω0n)+q(n)
[0007]其中为调制在信号上的第m(m∈[0,M
‑
1))个未知QAM 符号，A
m
和分别对应符号的幅度和相位。ω0＝2πΔfT
s
为信号频率，q(n)为方差为的加性高斯白噪声。设Δf是发射与本地振荡器之间的频率偏移，T
s
是采样频率。本专利技术所用的频率估计方法步骤如下：
[0008]步骤一，用长度为N的矩形窗口给x(n)加窗，得到一个N点采样时间序列 x(n)，其中n＝0,1,L,N，L为离散点个数；
[0009]步骤二，重写ω0＝2πl0/N＝2π(k0+δ0)/N，其中，l0为获得的正弦信号周期数，δ0和k0分别为归一化频率的分数和整数部分步骤三，让符号跳变发生在第L离散点上
[0010]步骤四，对x(n)在位置k处进行DFT变换，得到
[0011][0012]其中
[0013]步骤五，根据X(k)的最大值所在位置估计其中X(k)为S(k)在有噪情况下的观测值：令令以及代入步骤四得到
[0014][0015]步骤六，当已知任意三个不同的DFT单元S(k1)、S(k2)和S(k3)时，可用步骤五的式子得到三个线性方程。
[0016]步骤七，根据以下公式估计和
[0017][0018]其中，上标代表是估计的值。
[0019]步骤八，通过步骤七的式子计算和因为步骤五中令以及所以反代换得到ω0、和的估计值。
[0020]如前所述，任何三种不同的DFT单元都可以在所提出的方法中使用。在实际应用中，可以优先采用最大振幅X(k0)及其相邻的DFT单元X(k0+1)和 X(k0‑
1)，因为频谱泄露引起的噪声和干扰对这些DFT单元影响很小。
[0021]有益效果：本专利技术可以通过方程得到DFT单元与频率偏移的关系，使用三个不同的DFT单元来消除符号跳变的影响，从而实现对例如M
‑
QAM的符号跳变的单频复指数信号进行高精度的频率估计。
附图说明
[0022]图1是本专利技术用于对有符号跳变的单频复指数信号频率估计提供的方法流程图；
[0023]图2是信噪比为40dB，L为10，改变l0时本专利技术方法与同类型方法比较图；
[0024]图3是l0为1.1，L为10，改变信噪比时本专利技术方法与同类型方法比较图；
[0025]图4是信噪比为40dB，l0为1.1，改变L时本专利技术方法与同类型方法比较图。
具体实施方式
[0026]下面结合附图对本专利技术的技术方案做进一步的详细说明：
[0027]如图1所示为本专利技术用于对有符号跳变的的单频复指数信号频率估计提供的方法流程图，如图1所示，用于蓝牙信号频率偏移测试的方法，包括以下步骤：
[0028]步骤一，用长度为N的矩形窗口给x(n)加窗，得到一个N点采样时间序列 x(n)，其中n＝0,1,L,N。
[0029]步骤二，我们可以重写ω0＝2πl0/N＝2π(k0+δ0)/N，其中，l0为获得的正弦信号周期数，δ0和k0分别为归一化频率的分数和整数部分。求x(n)的N点 DFT，其N点DFT可以表示为：
[0030][0031]其中，S(k)是序列s(n)的DFT，Q(k)是噪声q(n)的DFT，
[0032]步骤三，让符号跳变发生在样本L上
[0033][0034]其中，和是任意两个不同的未知符号。
[0035]步骤四，在位置k处进行DFT变换，得到
[0036][0037]其中，
[0038]步骤五，举例说明，例：令令以及则可以把步骤四的式子重写为
[0039][0040]步骤六，当已知任意三个不同的DFT单元S(k1)、S(k2)和S(k3)时，可用步骤五的式
子得到如下线性方程：
[0041][0042]步骤七，根据以下公式估计和
[0043][0044]其中，上标代表是估计的值。
[0045]步骤八，通过步骤七的式子计算和因为步骤五中令以及所以ω0、和可用下式估计为：
[0046][0047]为了进一步说明该方法，做了以下仿真和实验，测试它与同类型方法相比的优越性和针对具体实例的鲁棒性。
[0048]对于所有仿真设置N＝128，L＝10，使用4
‑
QAM符号即 {1+j,
‑
1+j,1
‑
j,
‑1‑
j}，采用最大振幅X(k0)及其相邻DFT单元X(k0+1)和 X(k0‑
1)的三条谱线进行插值。图2
‑
图4为改变信噪比、l0、L三个不同参数时本专利技术所提供的方法与已有同类型方法的比较，其中(a)表示{1+j,
‑
1+j}， (b)表示{1+j,
‑1‑
j}，(c)表示{1+j,1
‑
j}。图2显示了信噪比为40dB，L 为10，改变l0时本专利技术方法与同类型方法性能比较。图3显示了l0为1.1，L为10，改变信噪比时本专利技术方法与同类型方法性能比较。图4显示了信噪比为40dB， l0为1.1，改变L时本专利技术方法与同类型方法性能比较。
[0049]在现实世界的IEEE802.1X射频一致性测试场景下，测试该方法的鲁棒性。试验台平台建立在Aerofle本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于DFT的符号跳变单频复指数信号频率估计方法，其特征在于，包括如下步骤：步骤1：采集N点离散时间信号s(n)，其中n＝0,1,L,N，L为离散点个数；步骤2：让符号跳变发生在第L离散点上其中，和是任意两个不同的未知符号，ω0为信号频率，A
m
和分别对应符号的幅度和相位，m＝1,2；步骤3：对s(n)在位置k处进行DFT变换，得到其中，步骤4：令以及代入公式(2)得到步骤5：以任意三个不同的已知的DFT单元S(k1)、S(k2)和S(k3)，分别根据公式(3)成立三个线性方程；步骤6：根...

【专利技术属性】
技术研发人员：王开，蒋威，雍婷，周海峰，刘珊，钱思雨，
申请(专利权)人：东南大学，
类型：发明
国别省市：