WO2020220804A1 - 一种可变采样率的傅里叶分析方法 - Google Patents

一种可变采样率的傅里叶分析方法 Download PDF

Info

Publication number
WO2020220804A1
WO2020220804A1 PCT/CN2020/076318 CN2020076318W WO2020220804A1 WO 2020220804 A1 WO2020220804 A1 WO 2020220804A1 CN 2020076318 W CN2020076318 W CN 2020076318W WO 2020220804 A1 WO2020220804 A1 WO 2020220804A1
Authority
WO
WIPO (PCT)
Prior art keywords
frequency
sampling
fourier analysis
sampling rate
signal
Prior art date
Application number
PCT/CN2020/076318
Other languages
English (en)
French (fr)
Inventor
曹继伟
李立毅
王凯思源
张成明
Original Assignee
哈尔滨工业大学
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 哈尔滨工业大学 filed Critical 哈尔滨工业大学
Publication of WO2020220804A1 publication Critical patent/WO2020220804A1/zh
Priority to US17/152,844 priority Critical patent/US20210141854A1/en

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/14Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/14Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • G06F17/141Discrete Fourier transforms
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/16Spectrum analysis; Fourier analysis
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/34Testing dynamo-electric machines
    • G01R31/343Testing dynamo-electric machines in operation

Definitions

  • Figure 2 is the flow chart of searching for the optimal sampling rate

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Engineering & Computer Science (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)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Measurement Of Resistance Or Impedance (AREA)

Abstract

一种可变采样率的傅里叶分析方法,包括以下步骤:S100通过对比初始设置的采样频率对所要分析的信号进行初步采样,进而采用傅里叶分析分析其基波频率及基波幅值;S200通过分析得到的基波幅值进行初步判定,确定满足整周期截断的采样频率,重新对信号进行采样;S300采样后的信号对其进行傅里叶分析,利用三谱线法确定最佳采样频率,即,既满足整周期截断,又满足无频谱泄露的最佳采样频率;S400重新对信号进行采样,得到各次谐波的频率及幅值组成。可变采样率的傅里叶分析方法可以采用较小的采样率实现高精度、高频信号的采集及分析,减小了采集系统的成本及加快了分析速度。

Description

一种可变采样率的傅里叶分析方法 技术领域
本发明属于高频信号测试领域,特别是涉及一种可变采样率的傅里叶分析方法。
背景技术
傅里叶分析最常见的即为离散傅里叶变换(DFT)和快速傅里叶变换(FFT)。DFT是重要的谐波分析工具,可以针对复杂信号的采样序列进行数学变换,分离基波信号与各谐波信号。通常为了保证频率分辨率,需要增加采样序列长度N,当N值较大时,DFT需要进行N2次复数乘法运算,所需的时间过长,若要保证良好的实时性则会对大幅增加对硬件的要求。
FFT算法采用蝶形运算方式,能在较短的时间内实现谐波检测,但频率分辨率低,且要求同步采样和整数周期截断。若能保证精准的同步采样,FFT针对谐波的测量可达到的精度很高;但对于间谐波的测量则取决于频率分辨率。通常,商业功率分析仪遵从IEC标准而设置的窗长度与采样率只能满足(1~10)Hz的频率分辨率,对于频率非整数的间谐波测量往往不够。非同步采样条件下,傅里叶方法固有的频谱泄漏和栅栏效应会大幅增加幅值和频率测量的误差,而在实际工程应用中实现严格的同步采样是难以实现的,因此如何减小频谱泄漏和栅栏效应是国内外学者研究的重点。
傅里叶方法的分析结果受频谱泄漏和栅栏效应影响较大,而且二者是相辅相成的,只有当被测频率分量恰好与频率轴单元重合时,才能得到精确的分析结果。通常,现有的谐波分析仪器和分析方法在低频领域会使用同步采样和栅栏效应的双重局限来实现该前提,但实则是为采样条件增加了强约束,缩小了采样率的有效作用空间,因此在高频领域的分析会受到采样率的限制。
为了使被测频率尽量与频率轴单元点重合,通常的谐波分析算法会要求满足同步采样和整周期截断两个条件,其中严格的同步采样是采样频率是所有频率分量的整数倍,否则一旦采样频率对某一次谐波形成非同步采样,则会对傅里叶分析的所有谱线结果产生影响。而实际工程中由于信号中包含谐波与间谐波,频率分量种类繁多,且各分量频率未知,很难实现严格的同步采样;同时整周期截断需要前期仪器的滤波来确定基波周期来计算截断时间窗的长度,在采取PWM或存在噪声的条件下,波形的过零点存在大幅振动,会导致基波周期的测量失准,而若加入滤波器的使用又会对主分量幅值产生一定的削减。
不发生频谱泄漏的充分必要条件即为被测频率与频率单元重合,同步采样与整周期截断都是为实现该前提而衍生出来的条件。对于工作频率达到几百甚至数千赫兹的高速电 机,如果为实现同步采样求取各谐波次数的最小公倍数,则可能会超出目前硬件采样率的极限。
因此本发明专利针对这种情况,提出一种新型的可变采样率的傅里叶分解方法,从而提高傅里叶分解方法的精度,有效减小对于高频信号尤其是谐波信号的频谱泄露和栅栏效应。
发明内容
本发明为解决高速电机的高频信号的分析采集与分析问题,提出一种可变采样率的傅里叶分析方法,从而提高傅里叶分解方法的精度,有效减小对于高频信号尤其是谐波信号的频谱泄露和栅栏效应。
本发明通过以下技术方案实现:一种可变采样率的傅里叶分析方法,所述傅里叶分析方法包括以下步骤:
S100通过对比初始设置的采样频率对所要分析的信号进行初步采样,进而采用傅里叶分析分析其基波频率及基波幅值;
S200通过分析得到的基波幅值进行初步判定,确定满足整周期截断的采样频率,重新对信号进行采样;
S300采样后的信号对其进行傅里叶分析,利用三谱线法确定最佳采样频率,即,既满足整周期截断,又满足无频谱泄露的最佳采样频率;
S400重新对信号进行采样,得到各次谐波的频率及幅值组成。
进一步的,在步骤S100中,具体的,通过公式f=np/60直接计算运行频率f,以基频估计值
Figure PCTCN2020076318-appb-000001
的2n倍为初始采样频率进行采样,
进一步的,所述步骤S100前还包括以下步骤:
步骤S000获得被测高速电机反馈的转速信号。
进一步的,所述最佳采样频率具有多个数值,且呈现周期性变化。
进一步的,所述最佳采样频率的周期性变化,随频率升高而周期逐渐增加。
本发明的有益效果在于:本发明设计了一种可变采样率的傅里叶分析方法,能够根据转速信号快速估计初始采样率,随后根据计算算法能够快速确定最优化采样率,能够有效减小高频信号的采集误差及分析误差,能够将栅栏效应和频谱泄露减小至0,有效提高傅里叶分析算法的分析精度。同时,采用该种算法可以采用较小的采样率实现高精度、高频信号的采集及分析,减小了采集系统的成本及加快了分析速度。
附图说明
图1为可变采样率FFT分解的计算流程图;
图2为最优化采样率搜索流程图;
图3为采样频率变化与基波幅值与频率的关系示意图;
图4为采样频率的变化与基频谱线之间的关系;
图5为采样频率变化时所呈现的周期性规律;
图6为主旁瓣分布的一种情况的示意图;
图7为主旁瓣分布的另一种情况的示意图;
图8为基于二分法的搜索原理示意图。
具体实施方式
下面将结合本发明实施例中的附图对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
参照图1所示,本发明提供了一种可变采样率的傅里叶分析方法,傅里叶分析方法包括以下步骤:
S100通过对比初始设置的采样频率对所要分析的信号进行初步采样,进而采用傅里叶分析分析其基波频率及基波幅值;
S200通过分析得到的基波幅值进行初步判定,确定满足整周期截断的采样频率,重新对信号进行采样;
S300采样后的信号对其进行傅里叶分析,利用三谱线法确定最佳采样频率,即,既满足整周期截断,又满足无频谱泄露的最佳采样频率;
S400重新对信号进行采样,得到各次谐波的频率及幅值组成。
在本部分优选实施例中,在步骤S100中,具体的,通过公式f=np/60直接计算运行频率f,以基频估计值
Figure PCTCN2020076318-appb-000002
的2n倍为初始采样频率进行采样,
在本部分优选实施例中,所述步骤S100前还包括以下步骤:
步骤S000获得被测高速电机反馈的转速信号。
在本部分优选实施例中,在步骤S300中,采样频率根据所需要判别的谐波频率及谐波幅值对象进行最优化调整,从而确定最佳采样频率。
在本部分优选实施例中,最佳采样频率具有多个数值,且呈现周期性变化。
在本部分优选实施例中,最佳采样频率的周期性变化,随频率升高而周期逐渐增加。
具体的,参照图1所示,根据被测高速电机反馈的转速信号,由公式f=np/60直接计算运行频率f,为了能够同时保证整周期截断和序列长度的要求,可以以基频估计值
Figure PCTCN2020076318-appb-000003
的2n倍为初始采样频率进行采样,同时最终的采样序列长度也要控制为N=2p(p>n),即信号周期数为2p-n,记作N p;而对于m次谐波,在序列中的周期数为m×N p。其中需要注意的是,基波在一个周期中被采样了2n个点,则m次谐波一个周期只被采样了2n/m个点,当谐波次数较大时,可能会无法满足采样定理,因此在设定采样频率时,需要根据电机特点对n进行适当调节,保证满足下式:
Figure PCTCN2020076318-appb-000004
式中,m max为所需分析的最高谐波次数。
Figure PCTCN2020076318-appb-000005
则实际基频应该对应第2 p-n条谱线,则此时m次谐波的频率则对应第m×2 p-n条谱线,即满足了与频率单元重合的条件。但由于转速反馈的估计值存在误差,一定会出现频谱泄漏结果,为了最大限度地减少频谱泄漏,本方法通过改变采样率来获得不同的信号样本序列,利用这些样本序列进行傅里叶分析,并寻找幅值结果的最大值。当幅值达到最大值时,频谱泄漏几乎全部消除,而经过不断修正的采样率,也逐渐逼近了最优采样率f sop=f 1×2 p-n
结合图1来说明,图1中三谱线分析是用来取代FFT的低计算量方法,由前文的分析可知,仅需针对基频实现优化,即可最大程度完成所有谐波分量的频谱泄漏的消除。因此,在每次更改采样率之后,仅需关注基频分量的幅值与频率情况,而无需进行FFT完成全频域的计算。三谱线分析会基于之前获得的基频信息来估计新分析结果中的谱线位置。
Figure PCTCN2020076318-appb-000006
式中[]取整数,可以四舍五入。求得对应谱线的幅值A *(k),其计算方法如下所示:
Figure PCTCN2020076318-appb-000007
Figure PCTCN2020076318-appb-000008
Figure PCTCN2020076318-appb-000009
然而
Figure PCTCN2020076318-appb-000010
与实际基频f 1的差异可能会导致取整函数的结果出现差错,使得对应谱线位置的次序相差1,即恰好计算的是主瓣两侧的谱线。为了避免出现这种差错,通常会计算第k条谱线左右相邻的谱线幅值,并找到三者中的最大值作为基频幅值的估计值。后续步骤中若再次改变采样率时,计算的谱线位置未必出错,因此并不是每次都需要计算三条谱线的幅值。为此,设置一个三谱线分析的开启阈值:
Figure PCTCN2020076318-appb-000011
式中A t为开启阈值,当计算得到的谱线幅值小于A t,即其旁瓣与主瓣几乎相等,则说明该谱线一定不是主瓣,需要另行计算寻找主瓣。而主瓣对应的频率单元即为新得到的基频分析结果:
Figure PCTCN2020076318-appb-000012
根据新得到的基频值f 1 *,按照之前的设定定义新的采样率2 n×f 1 *
参照图3所示,图3以8000Hz的采样频率采样一组基波频率为500Hz,幅值为10V;含有三次谐波(频率1500Hz,幅值为3.3V)、五次谐波(频率为2500Hz,1.7V);取n=4,p=9,对其进行分析的结果。图3可见各频率分量的幅值都会随着采样率的变化展现出一种类抛物线的波动规律,且谐波的波动周期明显较基波小很多;且当基波取得峰值时,各谐波也一定可以达到峰值,反之则不成立。若继续扩大采样率的变化范围,可以发现幅值的波动周期会随着采样率的增加而逐渐增大。图4说明利用8kHz采样基频为500Hz的信号,序列长度为512对应的是32个周期,此时基频对应的为频率轴上的第32根谱线;而其相邻的峰值由于也同样消除了频谱泄漏的影响,必定也使得基频与频率单元重合,而由于采样率的改变,频率单元也会发生相应的变化,故左右两侧的峰值对应的基频恰好为第33条谱线和第31条谱线。同理,基于以上论述,我们也可以得知,与8kHz相邻的左侧波谷值所对应的基频,应当恰好处第32与第33条谱线中央,而与其相邻的左侧波谷值所对应的基频,应当恰好处第31与第32条谱线中央。图5说明了不同的采样频率所处的位置。
参照图2、图6和图7所示,图6中左侧谱线幅值大于右侧谱线幅值,说明基频的真 实值在主瓣与左侧旁瓣之间;同理,图7中情况对应的则是基频的真实值在主瓣与右侧旁瓣之间,则通过这种方式,可以判断出2 n×f 1 *与f sop的方位关系。除此之外,主瓣与旁瓣的幅值还能用来进一步修正得到的基频测量值,对于图6中的情况,修正判断为:
Figure PCTCN2020076318-appb-000013
Figure PCTCN2020076318-appb-000014
对于图7中的情况,则变更频率为:
Figure PCTCN2020076318-appb-000015
Figure PCTCN2020076318-appb-000016
式中f 1 **为经过修正的基频值,A *(k),A *(k-1),A *(k+1)分别为主瓣和左右侧旁瓣的谱线幅值,λ l,λ r分别为向左侧和向右侧修正采样率时的修正系数。
参照图2和图8所示,经过三谱线法的修正后,f 1 **和基频真实值间的偏差已经很小了,但若继续使用三谱线法优化采样率,最后极有可能会不收敛,因为三谱线法的优化较为粗略,幅度较大,会导致最终优化值在实际值周围震荡,但无法收敛。因此,考虑使用二分法来快速锁定最佳采样率的精确位置。
二分法是一种适合在大数据量区间中搜索的快速方法,借助二分区域的原理可以以指数级速度缩小区域,如示意图8所示,其计算次数n 0取决于搜索区间L=[F a,F b]和计算精度e:
Figure PCTCN2020076318-appb-000017
式中的取整函数需要正向取整。由于经过两次修正后的采样率与最优采样率极为接近,因此可以大幅压缩二分法的搜索区间,从而减少搜索次数。如前,经过修正的2 n×f 1 **必定分布在f s c与f s d之间,因此可以通过使f s c与f s d二者相减,计算搜索区间。但该搜索区间过大,而2 n×f 1 **与最佳采样率f sop的具体差值未知,不宜盲目缩小区间,而利用三谱线修正基频的原理可以进一步缩小搜索区间。再次使用三谱线法分析以2 n×f 1 **采样得到 的新序列,同样按照上述步骤,根据旁瓣幅值判断采样率修正的方向,然后利用主瓣、旁瓣的幅值比例求取基频的修正系数λ:
Figure PCTCN2020076318-appb-000018
式中A **(k),A **(k-1),A **(k+1)为再次三谱线分析的主瓣、旁瓣幅值结果。为确保最优采样率f sop在区间中,将基频的修正量扩大至两倍。
则对应的另一区间端点为:
Figure PCTCN2020076318-appb-000019
搜索区间确定好后,根据二分法原理取两端点采样率和的二分之一为新采样率,并利用三谱线方法计算基于新采样率的新信号序列的主瓣幅值A op,并计算A op与两端点采样率对应的主瓣幅值的差值为:
Figure PCTCN2020076318-appb-000020
Δ a=A op-A a                 (16)
Δ b=A op-A b                 (17)
根据Δ a和Δ b的大小,可以确定更新区间的端点,并重新赋值:
Figure PCTCN2020076318-appb-000021
设置计算精度为10 -4级,当新区间两端点对应采样率的主瓣幅值结果间的差值Δ ab=|A a-A b|小于该精度时,即可认为二分法搜索完成。

Claims (5)

  1. 一种可变采样率的傅里叶分析方法,其特征在于,所述傅里叶分析方法包括以下步骤:
    S100通过对比初始设置的采样频率对所要分析的转速信号进行初步采样,进而采用傅里叶分析分析其基波频率及基波幅值;
    S200通过分析得到的基波幅值进行初步判定,确定满足整周期截断的采样频率,重新对信号进行采样;
    S300采样后的信号对其进行傅里叶分析,利用三谱线法确定最佳采样频率,即,既满足整周期截断,又满足无频谱泄露的最佳采样频率;
    S400重新对信号进行采样,得到各次谐波的频率及幅值组成。
  2. 根据权利要求1所述的一种可变采样率的傅里叶分析方法,其特征在于,在步骤S100中,具体的,通过公式f=np/60直接计算运行频率f,以基频估计值
    Figure PCTCN2020076318-appb-100001
    的2n倍为初始采样频率进行采样,
  3. 根据权利要求2所述的一种可变采样率的傅里叶分析方法,其特征在于,所述步骤S100前还包括以下步骤:
    步骤S000获得被测高速电机反馈的转速信号。
  4. 根据权利要求2所述一种可变采样率的傅里叶分析方法,其特征在于,所述最佳采样频率具有多个数值,且呈现周期性变化。
  5. 根据权利要求3所述一种可变采样率的傅里叶分析方法,其特征在于,所述最佳采样频率的周期性变化,随频率升高而周期逐渐增加。
PCT/CN2020/076318 2019-04-28 2020-02-23 一种可变采样率的傅里叶分析方法 WO2020220804A1 (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US17/152,844 US20210141854A1 (en) 2019-04-28 2021-01-20 Fourier Analysis Method with Variable Sampling Frequency

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201910349613.9A CN109975611A (zh) 2019-04-28 2019-04-28 一种可变采样率的傅里叶分析方法
CN201910349613.9 2019-04-28

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US17/152,844 Continuation US20210141854A1 (en) 2019-04-28 2021-01-20 Fourier Analysis Method with Variable Sampling Frequency

Publications (1)

Publication Number Publication Date
WO2020220804A1 true WO2020220804A1 (zh) 2020-11-05

Family

ID=67086847

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2020/076318 WO2020220804A1 (zh) 2019-04-28 2020-02-23 一种可变采样率的傅里叶分析方法

Country Status (3)

Country Link
US (1) US20210141854A1 (zh)
CN (1) CN109975611A (zh)
WO (1) WO2020220804A1 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113447713A (zh) * 2021-06-25 2021-09-28 南京丰道电力科技有限公司 一种基于傅式快速高精度的电力系统频率测量方法及装置
CN113466552A (zh) * 2021-07-14 2021-10-01 南京海兴电网技术有限公司 一种定间隔采样下频率跟踪方法

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109975611A (zh) * 2019-04-28 2019-07-05 哈尔滨工业大学 一种可变采样率的傅里叶分析方法
CN111736014B (zh) * 2020-07-20 2020-11-24 中国空气动力研究与发展中心低速空气动力研究所 一种频率信号采集和数据处理方法及装置
CN112817250A (zh) * 2020-12-28 2021-05-18 深圳市翌日科技有限公司 传感器数据采集方法及电路
CN113608023B (zh) * 2021-08-18 2024-08-20 国网湖南省电力有限公司 采样频率自适应的谐波检测方法及系统
CN115265691B (zh) * 2022-09-26 2022-12-09 中国空气动力研究与发展中心设备设计与测试技术研究所 一种科氏流量计振动频率跟踪方法及系统
CN116359605B (zh) * 2023-04-21 2023-09-12 中国计量科学研究院 一种基于二次加权的谐波信号分析方法
CN117169590B (zh) * 2023-08-15 2024-03-29 嘉兴市科讯电子有限公司 一种基于软件变采样率的电力谐波分析的方法和装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2232494B (en) * 1989-05-24 1993-01-27 Marconi Instruments Ltd Frequency spectrum determination using interpolation
CN105353215A (zh) * 2015-12-23 2016-02-24 合肥工业大学 基于Nuttall窗四谱线插值FFT的谐波检测方法
CN105866543A (zh) * 2016-03-31 2016-08-17 济南大学 一种消除基波、谐波对间谐波检测干扰的间谐波检测方法
CN106053937A (zh) * 2016-06-21 2016-10-26 电子科技大学 一种基于fft+ft的基波测频方法
CN109975611A (zh) * 2019-04-28 2019-07-05 哈尔滨工业大学 一种可变采样率的傅里叶分析方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI376507B (en) * 2008-06-06 2012-11-11 Univ Ishou Method of optimizing spectrum of power analyzer

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2232494B (en) * 1989-05-24 1993-01-27 Marconi Instruments Ltd Frequency spectrum determination using interpolation
CN105353215A (zh) * 2015-12-23 2016-02-24 合肥工业大学 基于Nuttall窗四谱线插值FFT的谐波检测方法
CN105866543A (zh) * 2016-03-31 2016-08-17 济南大学 一种消除基波、谐波对间谐波检测干扰的间谐波检测方法
CN106053937A (zh) * 2016-06-21 2016-10-26 电子科技大学 一种基于fft+ft的基波测频方法
CN109975611A (zh) * 2019-04-28 2019-07-05 哈尔滨工业大学 一种可变采样率的傅里叶分析方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
WANG, KAISIYUAN: "Research on the Harmonics Analyzing Algorithm of High-Frequency Electrical Power Signal", CHINA MASTER'S THESES FULL-TEXT DATABASE, 15 January 2019 (2019-01-15), ISSN: 1674-0246 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113447713A (zh) * 2021-06-25 2021-09-28 南京丰道电力科技有限公司 一种基于傅式快速高精度的电力系统频率测量方法及装置
CN113447713B (zh) * 2021-06-25 2023-03-07 南京丰道电力科技有限公司 一种基于傅式快速高精度的电力系统频率测量方法及装置
CN113466552A (zh) * 2021-07-14 2021-10-01 南京海兴电网技术有限公司 一种定间隔采样下频率跟踪方法
CN113466552B (zh) * 2021-07-14 2024-02-02 南京海兴电网技术有限公司 一种定间隔采样下频率跟踪方法

Also Published As

Publication number Publication date
CN109975611A (zh) 2019-07-05
US20210141854A1 (en) 2021-05-13

Similar Documents

Publication Publication Date Title
WO2020220804A1 (zh) 一种可变采样率的傅里叶分析方法
CN111222088B (zh) 一种改进的平顶自卷积窗加权电力谐波幅值估计方法
CN102539915A (zh) 时延傅立叶变换测频法精确计算电力谐波参数方法
CN101261292A (zh) 基于5项Rife-Vincent(I)窗双谱线插值FFT的基波与谐波检测方法
CN108427031B (zh) 基于多项式拟合及非干扰区域划分的间谐波检测方法
CN102818930B (zh) 一种高精度快速计算电力谐波参数的方法
CN103353550A (zh) 一种测量电力系统信号频率及谐波参数的方法
CN105486921A (zh) 凯撒三阶互卷积窗三谱线插值的谐波与间谐波检测方法
CN106226590A (zh) 一种电力系统同步相量测量方法
CN104502707A (zh) 一种基于三次样条插值的电力系统同步相量测量方法
Borkowski et al. Comparison of sine-wave frequency estimation methods in respect of speed and accuracy for a few observed cycles distorted by noise and harmonics
CN110954746A (zh) 一种基于四项Nuttall余弦窗的六插值FFT算法
CN114460527A (zh) 关联度延拓Hibert移相电子式互感器校验仪溯源方法和系统
CN109541304A (zh) 基于六项最小旁瓣窗插值的电网高次弱幅值谐波检测方法
CN105372492B (zh) 基于三条dft复数谱线的信号频率测量方法
CN106053936B (zh) 一种获取电学信号瞬时频率的方法及系统
CN117169590A (zh) 一种基于软件变采样率的电力谐波分析的方法和装置
CN117234176A (zh) 一种基于扫频试验的控制回路频域辨识方法
Serov et al. Comparative Analysis of Techniques for Reducing the Signal Parameters Measurement Error for the Quadrature Demodulation Application
Serov et al. Application of simulation modeling to determine the parameters of electrical signals by using a quadrature demodulator
CN105372493B (zh) 基于三条dft复数谱线的信号幅值和相位测量方法
CN114487589A (zh) 电网宽频信号自适应测量方法、装置及系统
Yue et al. Modified algorithm of sinusoid signal frequency estimation based on Quinn and Aboutanios iterative algorithms
CN113129912A (zh) 一种单音信号的检测方法
Ardeleanu et al. Fundamental frequency estimation based on mean values

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 20799226

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 20799226

Country of ref document: EP

Kind code of ref document: A1

32PN Ep: public notification in the ep bulletin as address of the adressee cannot be established

Free format text: NOTING OF LOSS OF RIGHTS PURSUANT TO RULE 112(1) EPC (EPO FORM 1205A DATED 10/05/2022)

122 Ep: pct application non-entry in european phase

Ref document number: 20799226

Country of ref document: EP

Kind code of ref document: A1