CN112685689A - 一种基于IpDFT的频率阶跃单频信号频率估计方法 - Google Patents
一种基于IpDFT的频率阶跃单频信号频率估计方法 Download PDFInfo
- Publication number
- CN112685689A CN112685689A CN202011406494.5A CN202011406494A CN112685689A CN 112685689 A CN112685689 A CN 112685689A CN 202011406494 A CN202011406494 A CN 202011406494A CN 112685689 A CN112685689 A CN 112685689A
- Authority
- CN
- China
- Prior art keywords
- frequency
- signal
- ipdft
- dft
- estimation method
- 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 65
- 239000000126 substance Substances 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 abstract description 7
- 238000005070 sampling Methods 0.000 description 7
- 238000004088 simulation Methods 0.000 description 5
- 238000012360 testing method Methods 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Landscapes
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
本发明公开了一种基于IpDFT的频率阶跃单频信号频率估计方法,用于频率阶跃单频信号的频率估计。利用本发明公开的方法可以通过方程得到DFT单元与阶跃变化频率的关系,使用六个不同的DFT单元来消除符号转换的影响,从而实现对例如FSK的频率阶跃变化的单频信号进行高精度的频率估计。
Description
技术领域
本发明涉及一种基于IpDFT的信号频率估计的方法,具体涉及一种用于频率阶跃单频信号的频率估计方法。
背景技术
近几十年来,频率的快速准确估计被广泛应用于电力系统、雷达/声呐和通信中。在过去几十年中,出现了大量的参数方法(时域方法)或非参数方法(频域方法)。用快速傅里叶变换(FFT)简化的插值离散傅里叶变换(IpDFT)是目前最常用的非参数方法之一。对于单频信号,IpDFT方法完全解决了栅栏效应引起的估计误差。在IpDFT方法中,默认的信号模型是时不变信号,这意味着信号的幅度、相位和频率在采样周期内保持不变。然而,时变信号在实际场景中也具有研究意义。例如,QAM信号、PSK信号和FSK信号可以看作是载波信号,其振幅、相位和频率随不同采样周期传输的信息而变化。在电力系统中,相位估计的标准对时变信号给出了详细的测试要求,包括幅值、相位和频率的阶跃变化。目前,基于IpDFT的时变信号频率估计的研究仍处于初级阶段。
发明内容
针对以上问题,本发明提出了一种基于IpDFT的用于频率阶跃变化的单频信号的频率估计方法。用线性方程给出了DFT单元与阶跃变化频率的关系,使用六个不同的DFT单元来消除符号转换的影响。
具有阶跃变化参数的单频信号可以描述为:
步骤一,用长度为N的矩形窗口给x(n)加窗,得到一个N点采样时间序列x(n),其中n=0,1,L,N。
步骤二,重写ω0=2πl0/N=2π(k0+δ0)/N,其中,l0为获得的信号周期数,δ0和k0分别l0的分数和整数部分。
步骤三,让符号转换发生在第L离散点上
步骤四,对x(n)在位置k处进行DFT变换,得到
步骤五,基于DFT的最大幅度指标可以进行粗频率估计。当已知任意六个不同的DFT单元,我们可以根据步骤四得出以下线性方程:
Xk=Wkη
步骤六,根据步骤五的方程,可估计η
虽然我们的方法可以使用任何六个不同的DFT-bins,但是在实际应用中推荐使用六个最大值的DFT-bins以获得最佳的估计结果。
有益效果:本发明可以通过方程得到DFT单元与阶跃变化频率的关系,使用六个不同的DFT单元来消除符号转换的影响,从而实现对例如FSK的频率阶跃变化的单频信号进行高精度的频率估计。
附图说明
图1是本发明用于对有阶跃频率的单频信号频率估计提供的方法流程图;
图2是频率不变,振幅和相位变化时改变不同参数本发明方法与现有方法比较图;
图3是信号在幅度阶跃变化时本发明方法与同类型方法比较图;
图4是信号在相位阶跃变化时本发明方法与同类型方法比较图;
图5是信号在频率阶跃变化时本发明方法与同类型方法比较图。
具体实施方式
下面结合附图对本发明的技术方案做进一步的详细说明:
如图1所示为本发明用于对具有阶跃变化参数的单频信号频率估计提供的方法流程图,频率估计的方法,包括以下步骤:
步骤一,用长度为N的矩形窗口给x(n)加窗,得到一个N点采样时间序列 x(n),其中n=0,1,L,N。
步骤二,我们可以重写ω0=2πl0/N=2π(k0+δ0)/N,其中,l0为获得的信号周期数,δ0∈[-0.5,0.5]和分别为l0的分数和整数部分。求 x(n)的N点DFT,其N点DFT可以表示为:
步骤三,让符号转换发生在样本L上
步骤四,在位置k处进行DFT变换,得到
Xk=Wkη
其中,η=[p1+p2,q1+q2,u1+u2,v1+v2,λ1+λ2,λ1λ2]T,且
步骤六,根据步骤五的方程,可估计η
为了进一步说明该方法,做了以下仿真和实验,测试它与同类型方法相比的优越性和针对具体实例的鲁棒性。
对于所有仿真,设置N=128,符号转换发生在第64个样本上(L=64),对于每个参数,执行3000此运行来评估统计属性。采用均方误差来评价所提出的估计量和其他方法的性能。均方误差的计算公式为:
首先,考虑频率保持不变,振幅和相位变化的情况。假设 且保持频率不变。图2中(a)描述了当信噪比为40dB,L=64,l0从0.5变换到3.5时,本方法与现有方法均方误差的比较。可以看出本方法和UFE 方法可以得到最小的均方误差,即-90dB。图2中(b)描述了当l0=1.06,L=64, 信噪比从0到40dB时,本方法与现有方法均方误差的比较。可以看出本方法和UFE 方法的均方误差随着信噪比的增加而线性减小,这意味着只有这两种方法才能在幅值和相位阶跃变化时正确估计频率。图2中(c)描述了当l0=1.06,信噪比为40dB, L从1到64时,本方法与现有方法均方误差的比较。本方法和UFE方法能提供最高的估计准确率。
表1、表2、表3是在频率阶跃变化但振幅、相位不变(U1=1,U2=1)情况下改变不同参数时,本方法与现有方法的均方误差比较。如表格所示,只有本发明提出的方法才能正确估计跳变点前后的频率。本方法能达到最小均方误差,并提供最佳的估计精度。
图3、图4是信噪比为40dB时信号在幅度和相位阶跃变化的仿真结果。其中,幅度变化步长为0.1,相位变化步长为π/18。阶跃变化发生在t=0.2s并在t=0.3s释放。图3结果表明,即使在阶跃变化发生和释放时,本方法法仍能跟踪参考频率。与现有方法相比,无论信号类型如何,本方法都能准确地估计出任何L值下的频率。图5 结果显示了本方法在信噪比为40dB时频率阶跃信号的性能。频率尖峰出现在t=1秒,并在t=2s时释放。当跳变点L≤88时,本算法可以精确估计ω1,当跳变点L≥40时,本算法可以精确估计ω2,当跳变点L∈[41,88]时,本方法可以同时精确估计ω1和ω2。对于IpDFT方法,跳频将影响整个128点采样周期。但是跳频对本方法的影响较小。结合如图3、图4和图5所示,需要突出的一个重要点是,当信号跳变的参数时,本发明的方法不需要收敛到期望值所需的持续时间。
具体实验部分,将在实际的5G射频一致性测试场景下测试所提出方法的鲁棒性。试验台平台建立在Aeroflex PXI模块上。在设置中,符号速率为160K的8-FSK信号的基带I/Q波形由运行在PC上的Aeroflex IQ Creator软件创建。调制器输出信号随后加载到矢量信号发生器(VSG)的非易失性存储器中。在低通滤波之后,I/Q信号波形被馈送到正交调制器以在下行链路N1工作频带(2110mhz)生成载波调制FSK 信号。矢量信号分析仪(VSA)用于以环回方式检测发射的RF信号,并将其下变频为载波频率为2110mhz的基带。基带中的采样频率为fs=17.6MHz,在这种情况下,基带中八个载波信号的归一化频率(ωi,i=1,L,8)为0.0277,0.0562,0.0848,0.1133, 0.1419,0.1705,0.1900和0.2276。我们设置N=128,L=64。表7为本方法与其他方法均方误差结果比较。
仿真和实验结果表明,本发明提出的方法能够实现高精度的频率估计。
表1信噪比为40dB,L为64时li(i=1,2)和均方误差的关系
表2 l1=1.5,l2=0.51且L为64时信噪比和均方误差的关系
表3 l1=1.5,l2=0.51且信噪比为40dB时L和均方误差的关系
表4信噪比为40dB且L为64时li(i=1,2)和均方误差的关系
表5 l1=1.5,l2=0.5且L为64时信噪比和均方误差的关系
表6 l1=1.5,l2=0.5且信噪比为40dB时L和均方误差的关系
表7不同8FSK信号符号组合情况下频率估计的偏差和均方误差
Claims (3)
2.根据权利要求1所述的一种基于IpDFT的频率阶跃单频信号频率估计方法,其特征在于,ω0=2πl0/N=2π(k0+δ0)/N,其中,l0为获得的信号周期数,δ0和k0分别l0的分数和整数部分。
3.根据权利要求1所述的一种基于IpDFT的频率阶跃单频信号频率估计方法,其特征在于,所述六个不同的DFT单元取前六个最大值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011406494.5A CN112685689B (zh) | 2020-12-03 | 2020-12-03 | 一种基于IpDFT的频率阶跃单频信号频率估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011406494.5A CN112685689B (zh) | 2020-12-03 | 2020-12-03 | 一种基于IpDFT的频率阶跃单频信号频率估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112685689A true CN112685689A (zh) | 2021-04-20 |
CN112685689B CN112685689B (zh) | 2024-03-22 |
Family
ID=75446054
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011406494.5A Active CN112685689B (zh) | 2020-12-03 | 2020-12-03 | 一种基于IpDFT的频率阶跃单频信号频率估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112685689B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1663211A (zh) * | 2003-03-27 | 2005-08-31 | 都科摩欧洲通信技术研究所有限公司 | 估计多个信道的设备和方法 |
CN1921466A (zh) * | 2006-09-08 | 2007-02-28 | 清华大学 | 利用旋转技术解决ofdm截取位置跳变的信道估计方法 |
US20200200806A1 (en) * | 2017-06-09 | 2020-06-25 | Ecole Polytechnique Federale De Lausanne (Epfl) | Method for estimating synchrophasors during static and dynamic conditions |
-
2020
- 2020-12-03 CN CN202011406494.5A patent/CN112685689B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1663211A (zh) * | 2003-03-27 | 2005-08-31 | 都科摩欧洲通信技术研究所有限公司 | 估计多个信道的设备和方法 |
CN1921466A (zh) * | 2006-09-08 | 2007-02-28 | 清华大学 | 利用旋转技术解决ofdm截取位置跳变的信道估计方法 |
US20200200806A1 (en) * | 2017-06-09 | 2020-06-25 | Ecole Polytechnique Federale De Lausanne (Epfl) | Method for estimating synchrophasors during static and dynamic conditions |
Non-Patent Citations (1)
Title |
---|
郭履翔: "基于DFT的实正弦信号频率估计算法的研究", 中国优秀硕士学位论文全文数据库(电子期刊), no. 6, 15 June 2020 (2020-06-15), pages 136 - 12 * |
Also Published As
Publication number | Publication date |
---|---|
CN112685689B (zh) | 2024-03-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104749432B (zh) | 基于聚焦s变换的多分量非平稳信号瞬时频率估计方法 | |
US10298338B2 (en) | Method for evaluating quality of radio frequency signals for stellite navigation system | |
CN110113278B (zh) | 一种基于全数字接收机的调制方式识别方法 | |
CN107085140B (zh) | 基于改进的SmartDFT算法的非平衡系统频率估计方法 | |
CN107204840B (zh) | 基于dft和迭代校正的正弦信号频率估计方法 | |
CN103323667B (zh) | 贝塞尔函数与虚拟阵列相结合的sfm信号的参数估计方法 | |
CN109342813A (zh) | 一种基于dft和二分法的正弦信号频率估计方法 | |
US20080151772A1 (en) | Display methods and apparatus for transmission characteristics | |
Fu et al. | Parameter Measurement of $ M $-Ary PSK Signals With Finite Rate of Innovation | |
EP3490206A1 (en) | Modulation method, apparatus applicable to ovxdm system and ovxdm system | |
CN108572352B (zh) | 一种基于欠采样的相位编码信号的参数估计方法 | |
CN112685689A (zh) | 一种基于IpDFT的频率阶跃单频信号频率估计方法 | |
US9294324B2 (en) | Signal generator and signal generation method with cyclic prefix generation | |
CN113810326A (zh) | 时偏估计的方法、装置、电子设备及存储介质 | |
CN111539323B (zh) | 一种循环前缀线性调频信号的频率估计方法与装置 | |
US20180278344A1 (en) | Error measurement method using a time-variant stopband test signal | |
CN102185811B (zh) | 一种载波频率估计方法 | |
CN112883787B (zh) | 一种基于频谱匹配的短样本低频正弦信号参数估计方法 | |
CN113406386B (zh) | 一种基于数字下变频的信号频率精确估计方法 | |
US10379162B1 (en) | System for performing modulation analysis without using a modulated signal | |
CN112671528B (zh) | 一种基于dft的符号跳变单频复指数信号频率估计方法 | |
CN108169715B (zh) | 一种同相和正交通道相位不平衡度确定方法及系统 | |
CN101977173A (zh) | Ofdm信号的带宽盲估计方法 | |
CN114338317A (zh) | 一种色噪声中的载波频率获取方法 | |
CN111522000A (zh) | 一种基于OFDM-chirp波形的目标检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |