CN111274534B - 一种改进全泄漏抑制的短时dft插值方法 - Google Patents
一种改进全泄漏抑制的短时dft插值方法 Download PDFInfo
- Publication number
- CN111274534B CN111274534B CN202010150216.1A CN202010150216A CN111274534B CN 111274534 B CN111274534 B CN 111274534B CN 202010150216 A CN202010150216 A CN 202010150216A CN 111274534 B CN111274534 B CN 111274534B
- Authority
- CN
- China
- Prior art keywords
- frequency
- spectrum
- dft
- short
- leakage
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 25
- 230000005764 inhibitory process Effects 0.000 title claims description 8
- 238000001228 spectrum Methods 0.000 claims abstract description 60
- 230000003595 spectral effect Effects 0.000 claims abstract description 26
- 230000001629 suppression Effects 0.000 claims abstract description 6
- 238000005070 sampling Methods 0.000 claims description 15
- 238000012935 Averaging Methods 0.000 claims description 3
- 230000000750 progressive effect Effects 0.000 claims description 3
- 230000006978 adaptation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Discrete Mathematics (AREA)
- Operations Research (AREA)
- Complex Calculations (AREA)
Abstract
本申请公开了一种改进全泄漏抑制的短时DFT插值方法,在本申请中首先基于最大旁瓣衰减窗函数频谱的旁瓣一致衰减特性,在高精度近似基础上推导得到了其线性比例及递推特性。根据这种特性能够在谐波分量未知前提下,将其泄漏干扰进行有效参数化近似,且仅与被测频率谱线位置的级数线性相关,使其能够适用于后续多谱线插值校正的解析过程。再次利用最大旁瓣衰减窗函数频谱的线性比例特性,采用多谱线方程组解析得到与单频率信号插值校正类似的显示表达式,进而突破了单频解析信号模型这一先决条件对插值算法的限制。最后,分别对实部和虚部采用相同解析方式,得到对应解析解,进而采用数学平均方式进一步提高短时CiR下频谱校正精度。
Description
技术领域
本申请涉及频谱校正技术领域,尤其涉及一种改进全泄漏抑制的短时DFT插值方法。
背景技术
目前常用的离散谱线幅值比值插值校正方法,必须在对应于待测频率附近的离散峰值谱线近似满足单频解析信号这一先决条件下,利用非线性拟合求取频率偏移量。当待测频率较低或对应离散频率接近奈奎斯特频率时,DFT(Discrete Fourier Transform,离散傅里叶变换)的频谱共轭对称特性使得来自待测信号负频率的“频谱短泄漏”分量量级急剧增大。负频率的频谱短泄漏干扰将导致基于单频解析信号模型的非线性拟合方法精度低。
为此,需要针对负频率的频谱短泄漏干扰进行有效抑制,其原理在利用正负频率共轭对称特性及任意H阶最小旁瓣衰减窗的高精度近似下的线性比例特性,通过谱线方程组解析得到频率校正的通用解析表达式。
因为离散采样信号经窗函数加权DFT后对应离散频率分布的不确定,以及窗函数频谱旁瓣衰减的非线性特性,使得无法完全消除其它频率分量所带来的频谱泄漏干扰,特别是短时CiR的情况下的频谱泄漏干扰。
目前,常用方法在于利用最大旁瓣衰减窗函数频谱的旁瓣一致衰减特性,将所有泄漏干扰以泰勒级数多项式近似,进而递推得到多谱线插值校正方法。但是,该方法的来自待测信号负频率的频谱短泄漏也引入多项式近似,明显对于低频率校正,此时占主导的负频率的频谱短泄漏会使得该方法存在偏差,影响含谐波的短时CiR下频谱校正精度。
发明内容
本申请提供了一种改进全泄漏抑制的短时DFT插值方法,以解决含谐波的短时CiR下频谱校正精度低的技术问题。
为了解决上述技术问题,本申请实施例公开了如下技术方案:
本申请实施例公开了一种改进全泄漏抑制的短时DFT插值方法,包括:对频谱进行离散采样,生成离散采样信号;
对所述离散采样信号进行加窗DFT,得到DFT频谱,建立所述DFT频谱对应的谐波正负频率的频谱长泄漏分量的泰勒级数多项式等效模型;
基于所述泰勒级数多项式等效模型,选择主瓣谱线,建立主瓣及附近局部峰值谱线的全泄漏频谱模型;
根据所述泰勒级数多项式等效模型和所述全泄漏频谱模型,计算得到频率估计的通用显示表达式;
对所述频率估计的通用显示表达式进行数学平均,得到全泄漏抑制的短时DFT插值的频率校正公式。
进一步的,对所述离散采样信号进行DFT,得到DFT频谱,建立所述DFT频谱对应的谐波正负频率的频谱长泄漏分量的泰勒级数多项式等效模型,包括:对所述离散采样信号进行加窗DFT,得到DFT频谱,其中,窗函数为H阶最大旁瓣衰减窗;
利用H阶最大旁瓣衰减窗函数频谱的旁瓣一致单调衰减特性,建立所述DFT频谱对应的谐波正负频率的频谱长泄漏分量的泰勒级数多项式等效模型。
进一步的,根据所述泰勒级数多项式等效模型和所述全泄漏频谱模型,计算得到频率估计的通用显示表达式,包括:
根据所述泰勒级数多项式等效模型中泰勒级数多项式的阶数,确定所述主瓣及附近局部峰值谱线数量;
分别建立所述峰值谱线的实部和虚部方程组,通过对所述方程组求解,得到频率估计的通用显示表达式,其中,所述频率估计的通用显示表达式包括实部解析结果和虚部解析结果。
进一步的,所述H阶最大旁瓣衰减窗的旁瓣渐进衰减速率为6(2H-1)dB/倍频程。
与现有技术相比,本申请的有益效果为:
在本申请中首先基于最大旁瓣衰减窗函数频谱的旁瓣一致衰减特性,在高精度近似基础上推导得到了其线性比例及递推特性。根据这种特性能够在谐波分量未知前提下,将其泄漏干扰进行有效参数化近似,且仅与被测频率谱线位置的级数线性相关,使其能够适用于后续多谱线插值校正的解析过程。进而,再次利用最大旁瓣衰减窗函数频谱的线性比例特性,采用多谱线方程组解析得到与单频率信号插值校正类似的显示表达式,进而突破了单频解析信号模型这一先决条件对插值算法的限制。最后,分别对实部和虚部采用相同解析方式,得到对应解析解,进而采用数学平均方式进一步提高含谐波的短时CiR下频谱校正精度。
应当理解的是,以上的一般描述和后文的细节描述仅是示例性和解释性的,并不能限制本申请。
附图说明
为了更清楚地说明本申请的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,对于本领域普通技术人员而言,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本申请实施例提供的一种改进全泄漏抑制的短时DFT插值方法流程示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
如图1所示,本申请实施例提供了一种改进全泄漏抑制的短时DFT插值方法,包括:对频谱进行离散采样,生成离散采样信号。
对所述离散采样信号进行加窗DFT,得到DFT频谱,建立所述DFT频谱对应的谐波正负频率的频谱长泄漏分量的泰勒级数多项式等效模型。具体包括:对所述离散采样信号进行加窗DFT,得到DFT频谱,其中,窗函数为H阶最大旁瓣衰减窗。利用H阶最大旁瓣衰减窗函数频谱的旁瓣一致单调衰减特性,建立所述DFT频谱对应的谐波正负频率的频谱长泄漏分量的泰勒级数多项式等效模型。
基于所述泰勒级数多项式等效模型,选择主瓣谱线,建立主瓣及附近局部峰值谱线的全泄漏频谱模型。
根据所述泰勒级数多项式等效模型和所述全泄漏频谱模型,计算得到频率估计的通用显示表达式。具体包括:根据所述泰勒级数多项式等效模型中泰勒级数多项式的阶数,确定所述主瓣及附近局部峰值谱线数量。分别建立所述峰值谱线的实部和虚部方程组,通过对所述方程组求解,得到频率估计的通用显示表达式,其中,所述频率估计的通用显示表达式包括实部解析结果和虚部解析结果。
对所述频率估计的通用显示表达式进行数学平均,得到全泄漏抑制的短时DFT插值的频率校正公式。
本实施例提供的改进全泄漏抑制的短时DFT插值方法采用的是旁瓣渐进衰减速率为6(2H-1)dB/倍频程的H阶最大旁瓣衰减窗,对应频谱的高精度近似为:
式(1)中,N为采样点数,v为待估频率偏移,ah为加权系数,且
式(1)表明为一个复数,其对应实部和虚部应具有相似的线性比例特性。
上述近似,使在保证其旁瓣一致渐进衰减特性及精度要求的前提下,最大旁瓣衰减窗频谱具有两个重要的线性比例递推特征:
将式(3)右侧以泰勒级数多项式展开,分别表示为:
其中,J∈N+,阶数J用以控制近似精度。
不失一般性,此处设信号s(n)仅含有f1和f2两个频率分量,且f1<f2。离散频谱局部最大值分别处于第l1和l2根谱线,则对应l1根DFT谱线值为:
利用频谱函数W正负对称特性及式(16),式(5)中第二大部分大括号内为谐波分量f2,其对应的正负频率W加权旁瓣衰减分量可分别近似为:
将式(6)带入对应l1根DFT谱线,即可得到全泄漏频谱模型为:
明显,该模型不仅考虑了共轭负频率的频谱泄漏分量;同时,在谐波参数未知前提下使得来自谐波分量的频谱泄漏干扰得以有效参数化,且仅与被测频率谱线位置l1的级数线性相关。
上述模型分为实部和虚部,对应为:
选择J+3根主瓣范围附近离散谱线,分别为:kp-1、kp、kp+1…kp+J+1,对应实部和虚部形成相同形式方程组CY=0。
以实部为例,对应系数矩阵分别表示为C=[C1 C2],YT=[Y1 Y2],其中:
Y2=[x0 x1…xJ-1] (12)
根据线性代数中的克莱姆法则,非零解的存在要求齐次线性方程组系数矩阵的行列式为零,即|C|=0。进而可得上述方程的解为:
仍以实部为例,其中Γ1和Γ2分别为:
Γ1=(2H+J-1)·[ΔJ+1SRe(l1)(H+kp+J)2-ΔJ+1SRe(kp-1)(H-kp)2+ΔJSRe(kp)(2H+J-1)(2H+J)]
Γ2=(2H+J-1)·ΔJ+2SRe(kp-1)
符号ΔJ表示J阶前向差分,具有以下性质:
同理可得虚部对应解析解
进而,在上述对应实部和虚部解析解基础上,再采用“数学平均”方式,即
根据离散采样性质,将上述方程的解换算为被测频率,即为频率校正公式,即:
其中,fS为采样频率,N为采样点数。
在本申请中首先基于最大旁瓣衰减窗函数频谱的旁瓣一致衰减特性,在高精度近似基础上推导得到了其线性比例及递推特性。
根据这种特性能够在谐波分量未知前提下将其泄漏干扰进行有效参数化近似(且仅与被测频率谱线位置l1的级数线性相关),使其能够适用于后续多谱线插值校正的解析过程;进而,再次利用最大旁瓣衰减窗函数频谱的线性比例特性,采用多谱线方程组解析得到与单频率信号插值校正类似的显示表达式(且仅与被测频率主瓣范围内的谱线、窗函数阶数H和近似阶数J相关),进而突破了单频解析信号模型这一先决条件对插值算法的限制;最后,分别对实部和虚部采用相同解析方式,得到对应解析解,进而采用“数学平均”方式进一步提高含谐波的短时CiR下频谱校正精度。
由于以上实施方式均是在其他方式之上引用结合进行说明,不同实施例之间均具有相同的部分,本说明书中各个实施例之间相同、相似的部分互相参见即可。在此不再详细阐述。
需要说明的是,在本说明书中,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的电路结构、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种电路结构、物品或者设备所固有的要素。在没有更多限制的情况下,有语句“包括一个……”限定的要素,并不排除在包括所述要素的电路结构、物品或者设备中还存在另外的相同要素。
本领域技术人员在考虑说明书及实践这里发明的公开后,将容易想到本申请的其他实施方案。本申请旨在涵盖本发明的任何变型、用途或者适应性变化,这些变型、用途或者适应性变化遵循本申请的一般性原理并包括本申请未公开的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的,本申请的真正范围和精神由权利要求的内容指出。
以上所述的本申请实施方式并不构成对本申请保护范围的限定。
Claims (2)
1.一种改进全泄漏抑制的短时DFT插值方法,其特征在于,包括:
对频谱进行离散采样,生成离散采样信号;
对所述离散采样信号进行加窗DFT,得到DFT频谱,其中,窗函数为H阶最大旁瓣衰减窗;
利用H阶最大旁瓣衰减窗函数频谱的旁瓣一致单调衰减特性,建立所述DFT频谱对应的谐波正负频率的频谱长泄漏分量的泰勒级数多项式等效模型;所述泰勒级数多项式等效模型为:
式中,v为待估频率偏移;
基于所述泰勒级数多项式等效模型,选择主瓣谱线,利用频谱函数W正负对称特性建立主瓣及附近局部峰值谱线的全泄漏频谱模型;所述全泄漏频谱模型为:
式中,l1表示为DFT谱线的数量;
根据所述泰勒级数多项式等效模型中泰勒级数多项式的阶数,确定所述主瓣及附近局部峰值谱线数量;
分别建立所述峰值谱线的实部和虚部方程组,通过对所述方程组求解,得到频率估计的通用显示表达式,其中,所述频率估计的通用显示表达式包括实部解析结果和虚部解析结果;
对所述频率估计的通用显示表达式进行数学平均,得到全泄漏抑制的短时DFT插值的频率校正公式。
2.根据权利要求1所述的短时DFT插值方法,其特征在于,所述H阶最大旁瓣衰减窗的旁瓣渐进衰减速率为6(2H-1)dB/倍频程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010150216.1A CN111274534B (zh) | 2020-03-06 | 2020-03-06 | 一种改进全泄漏抑制的短时dft插值方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010150216.1A CN111274534B (zh) | 2020-03-06 | 2020-03-06 | 一种改进全泄漏抑制的短时dft插值方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111274534A CN111274534A (zh) | 2020-06-12 |
CN111274534B true CN111274534B (zh) | 2023-11-21 |
Family
ID=71000523
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010150216.1A Active CN111274534B (zh) | 2020-03-06 | 2020-03-06 | 一种改进全泄漏抑制的短时dft插值方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111274534B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111884965A (zh) * | 2020-07-22 | 2020-11-03 | 云南电网有限责任公司电力科学研究院 | 基于全泄漏抑制的频谱校正方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101441233A (zh) * | 2008-12-02 | 2009-05-27 | 湖南海兴电器有限责任公司 | 基于Kaiser窗双谱线插值FFT的基波与谐波检测方法 |
CN103308766A (zh) * | 2013-05-15 | 2013-09-18 | 湖南大学 | 一种基于凯撒自卷积窗双谱线插值fft谐波分析方法及其装置 |
CN106018956A (zh) * | 2016-08-10 | 2016-10-12 | 北京妙微科技有限公司 | 一种加窗谱线插值的电力系统频率计算方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11372038B2 (en) * | 2016-11-11 | 2022-06-28 | University Of Tennessee Research Foundation | Magnetic field based wireless GMD/EMP-E3 impact monitoring device |
-
2020
- 2020-03-06 CN CN202010150216.1A patent/CN111274534B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101441233A (zh) * | 2008-12-02 | 2009-05-27 | 湖南海兴电器有限责任公司 | 基于Kaiser窗双谱线插值FFT的基波与谐波检测方法 |
CN103308766A (zh) * | 2013-05-15 | 2013-09-18 | 湖南大学 | 一种基于凯撒自卷积窗双谱线插值fft谐波分析方法及其装置 |
CN106018956A (zh) * | 2016-08-10 | 2016-10-12 | 北京妙微科技有限公司 | 一种加窗谱线插值的电力系统频率计算方法 |
Non-Patent Citations (3)
Title |
---|
Yong Shao."Time-Domain Interpolation Correction Method for Wide-Range Frequency Offset Power Grid".《2019 Chinese Control Conference (CCC)》.2019,全文. * |
李海峰 ; 刘俊勇 ; 徐斌 ; 付东丰 ; 王泽 ; .消除外部环境因素干扰的MOA在线监测参数修正方法.高电压技术.2018,(08),全文. * |
陈欣."电力系统谐波检测的研究及应用".《中国优秀硕士学位论文全文数据库信息科技辑》.2020,全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN111274534A (zh) | 2020-06-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6230106B1 (en) | Method of characterizing a device under test | |
CN109471095B (zh) | 一种基于快速迭代插值的fmcw雷达距离估计方法 | |
CN109100755B (zh) | 高精度gnss接收机射频前端群时延畸变的校正方法 | |
EP2378667A2 (en) | Method for compensating a frequency characteristic of an arbitrary waveform generator | |
US10444327B2 (en) | Method for determining parameters of a compression filter and associated multi-channel radar | |
KR101294771B1 (ko) | 크기 측정 데이터를 이용하는 필터 등화 | |
CN111274534B (zh) | 一种改进全泄漏抑制的短时dft插值方法 | |
KR100347699B1 (ko) | 전력스펙트럼밀도추정방법및장치 | |
Matania et al. | Novel approaches for the estimation of the spectrum background for stationary and quasi-stationary signals | |
Varlamov et al. | Research of the Hybrid Filter Bank Frequency Equalizer Transfer Function Coefficients Estimation Quality | |
Zhu et al. | Noise influence on estimation of signal parameter from the phase difference of discrete Fourier transforms | |
CN112213560A (zh) | 一种基于z-adaline的高精度电网宽频信号测量方法 | |
CN107315714B (zh) | 一种去卷积功率谱估计方法 | |
CN101308175A (zh) | 相位谱分析仪 | |
Ivanov et al. | Mitigation of Dispersion Distortions of Transionospheric Communication Channels when Total Electron Content Measurements are Corrupted with Stochastic Error | |
CN108801296B (zh) | 基于误差模型迭代补偿的传感器频响函数计算方法 | |
US20020123869A1 (en) | Simulation method and test arrangement for determining nonlinear signal distortion | |
Orović et al. | A class of highly concentrated time-frequency distributions based on the ambiguity domain representation and complex-lag moment | |
CN116299110A (zh) | 一种用于双指数波形式电磁脉冲测量的时域标定方法 | |
US20040189326A1 (en) | Vector-detecting apparatus and impedance measuring apparatus | |
CN114114231A (zh) | 一种基于fft-czt的雷达测距方法 | |
CN108008360B (zh) | 一种幅度加权的非线性调频波形设计方法 | |
Ferrero et al. | A fast frequency-domain interpolation method for the evaluation of the frequency and amplitude of spectral components | |
CN113702701B (zh) | 一种基于梳状波信号的幅相特性测量方法 | |
CN114265017B (zh) | 一种基于数字信号处理的相噪测量方法 |
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 |