CN106053936A - 一种获取电学信号瞬时频率的方法及系统 - Google Patents
一种获取电学信号瞬时频率的方法及系统 Download PDFInfo
- Publication number
- CN106053936A CN106053936A CN201610434916.7A CN201610434916A CN106053936A CN 106053936 A CN106053936 A CN 106053936A CN 201610434916 A CN201610434916 A CN 201610434916A CN 106053936 A CN106053936 A CN 106053936A
- Authority
- CN
- China
- Prior art keywords
- frequency
- value
- electrical signal
- primary iteration
- instantaneous
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/02—Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
本发明涉及一种获取电学信号瞬时频率的方法及系统,灵活运用短时傅里叶变换方法,确定初始迭代频率,并将所述初始迭代频率应用于用于获取瞬时频率的迭代算法之中;本发明能够提供更加精确的初始迭代频率,加速迭代过程的收敛速度,并且降低窗口长度和采样频率对于迭代结果的影响。本方法对信号的阶跃变化敏感,但是对于下列三种干扰因素不敏感:电压或电流信号被其他信号调制、包含衰减的直流偏执和噪声,可见,本方法非常适用于电力系统这种不稳定状态情况下的瞬时频率的跟踪。
Description
技术领域
本发明涉及电力领域,尤其涉及一种获取电学信号瞬时频率的方法及系统。
背景技术
电力系统中,电压和电流信号的频率值是最重要的状态参量和操作指标,它也是电力系统发生孤岛现象的标示。频率值是电力系统的安全、稳定运行的一个至关重要的参数,能够反映整个系统的电能质量和安全情况,过高或过低的频率都会对电力系统带来危害,比如频率降低0.1Hz意味着800MW的电力缺失。在过去的几十年里,在电力工程领域的频率测量方面,研究者做了大量的研究工作,这些研究工作在时域领域收获颇丰,但是在变换域,能够用于频率跟踪的主要有传统傅里叶变换、最小均方误差法和短时傅里叶变换,而这些方法在使用过程中都存在一定缺陷。
最小均方误差法中,利用一个矩阵的伪逆运算可以确定信号频率。利用泰勒展开式构造的简化表达式可以使得最小均方误差法能够应用于在线模式和非在线模式。最小均方误差法的测量精度由每个采样间隔内的迭代次数和硬件的计算能力决定,这种方法相较于传统傅里叶变换能够解决非稳定状态下的信号测量问题,但是这种方法在非稳定状态中性能较差,因为泰勒展开后的算法需要更长的数据,而较长的数据窗口包含更多的不良数据,故最小均方误差法迭代速度较慢,并且受窗口长度和采样频率的影响很大。
短时傅里叶变换(SDFT:Short-time DFT)是一种基于窗口的运算,要求在每个周期内包含整数个采样点。如果一个窗口对应一个周期,包含了非整数个采样点,SDFT将产生错误,而这种情况是电力传输系统中普遍存在的。在相对稳定的状态下借助最小二乘法,SDFT能够连续的估计出瞬时频率、抑制谐波和平滑噪声,而对于非平稳状态,简单的SDFT算法无能为力。
可见,不论是最小均方误差法还是短时傅里叶变换都不宜直接应用于电力系统瞬时频率的获取,有必要综合最小均方误差法和短时傅里叶变换的优势,研究出一种适用于非稳态频率测量的准确高效的方法。
发明内容
为此,本发明提供了一种获取电学信号瞬时频率的方法及系统,通过灵活运用短时傅里叶变换方法,确定初始迭代频率,并将所述初始迭代频率应用于用于获取瞬时频率的迭代算法之中,所述初始迭代频率能够加速迭代过程的收敛速度,并且降低窗口长度和采样频率对于迭代结果的影响。
一种获取电学信号瞬时频率的方法,包括以下步骤:
S1、预先设定好采样频率,并根据所述采样频率对电学信号进行采样;
S2、获取不同窗口长度下的频率标定值,所述频率标定值fi=f/Li,Li为不同的数据长度,f为采样信号频率,fi为不同数据长度下的频率标定值;
S3、按照短时傅立叶变换算法(SDRT)获取采样时刻的不同窗口长度下的频率计算值;
S4、计算S2和S3中相同数据窗口长度下频率标定值与频率计算值的差异值,确定最小差异值对应的频率计算值;
S5、将最小差异值对应的频率计算值作为采样时刻的初始迭代频率;
S6、根据所述初始迭代频率计算采样时刻的瞬时频率。
优选的,S6包括以下步骤:
S61.根据所述初始迭代频率获取所述采样时刻的第一次频率估计值
S62、用所述第一次频率估计值代替所述初始迭代频率,重复步骤S61获取所述采样时刻的第二次频率估计值
S63、按照S62中所述方法依次获取下一次的频率估计值,当获取频率估计值的次数达到预设次数或相邻频率估计值之间的差值不大于预设阈值时,停止获取下一次的频率估计值;
S64、将最后一次获取的频率估计值作为所述采样时刻的瞬时频率值。
优选的,S61中包括如下步骤:
S611、获取过渡参数X1和X2,其中A1和是所述电学信号的初始振幅和初始相位;
S612、根据所述过渡参数获取所述第一次频率估计值其中, 若或者有否则
优选的,在进行采样之前,让所述电学信号通过低通滤波器,所述低通滤波器用于将频率高于f1Hh的信号滤除,f1是系统的标称频率,Hh是输入信号包含的最高谐波次数。
优选的,所述电学信号包括电压信号或电流信号。
一种获取电学信号瞬时频率的系统,包括采样频率设定模块、采样模块、初始迭代频率获取模块和瞬时频率计算模块,
所述采样频率设定模块用于电学信号的初始频率设定采样频率f,
所述采样模块根据采样频率f对电学信号进行采样,并将采样结果传输至初始迭代频率获取模块,
所述初始迭代频率获取模块用于获取初始迭代频率,并将所述初始迭代频率传输至瞬时频率计算模块,由所述瞬时频率计算模块计算出电学信号的瞬时频率。
优选的,所述初始迭代频率获取模块包括频率标定值获取单元,频率计算值获取单元和初始迭代频率获取单元,
所述频率标定值获取单元用于根据采样信号频率和窗口长度获取不同窗口长度下的频率标定值,并将所述频率标定值传输至初始迭代频率获取单元;所述频率计算值获取单元用于根据短时傅里叶变换算法获取不同窗口长度下的频率计算值,并将所述频率计算值传输至初始迭代频率获取单元;所述初始迭代频率获取单元用于根据所述频率标定值和频率计算值获取初始迭代频率。
优选的,所述采样频率设定模块设定的采样频率满足公式,f≥2Hhf1,Hh是所述电学信号中包含的最高谐波次数,f1是系统的标称频率,。
优选的,还包括低通滤波器,所述低通滤波器用于将频率高于Hhf1的信号在进入采样模块之前滤除。
本发明的有益效果是:
本发明提供了一种获取电学信号瞬时频率的方法及系统,灵活运用短时傅里叶变换方法,确定初始迭代频率,并将所述初始迭代频率应用于用于获取瞬时频率的迭代算法之中;本发明能够提供更加精确的初始迭代频率,加速迭代过程的收敛速度,并且降低窗口长度和采样频率对于迭代结果的影响。本方法对信号的阶跃变化敏感,但是对于下列三种干扰因素不敏感:电压或电流信号被其他信号调制、包含衰减的直流偏执和噪声,可见,本方法非常适用于电力系统这种不稳定状态的瞬时频率的跟踪。
附图说明
图1是一种获取电学信号瞬时频率的方法的流程图;
图2是一种获取电学信号瞬时频率的系统结构图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作进一步地详细描述。
在一个实施例中,一种获取电学信号瞬时频率的方法如图1所示,包括以下步骤:
S1、预先设定好采样频率,并根据所述采样频率对电学信号进行采样;
S2、获取不同窗口长度下的频率标定值,所述频率标定值fi=f/Li,Li为不同的数据长度,f为采样信号频率,fi为不同数据长度下的频率标定值;
S3、按照短时傅立叶变换算法(SDRT)获取采样时刻的不同窗口长度下的频率计算值;
S4、计算S2和S3中相同数据窗口长度下频率标定值与频率计算值的差异值,确定最小差异值对应的频率计算值;
S5、将最小差异值对应的频率计算值作为采样时刻的初始迭代频率;
S6、根据所述初始迭代频率计算采样时刻的瞬时频率。
具体地,S6包括以下步骤:
S61.根据所述初始迭代频率获取所述采样时刻的第一次频率估计值
S62、用所述第一次频率估计值代替所述初始迭代频率,重复步骤S61获取所述采样时刻的第二次频率估计值
S63、按照S62中所述方法依次获取下一次的频率估计值,当获取频率估计值的次数达到预设次数或相邻频率估计值之间的差值不大于预设阈值时,停止获取下一次的频率估计值;
S64、将最后一次获取的频率估计值作为所述采样时刻的瞬时频率值。
具体地,S61中包括如下步骤:
S611、获取过渡参数X1和X2,其中A1和是所述电学信号的初始振幅和初始相位;
S612、根据所述过渡参数获取所述第一次频率估计值其中, 若或者有否则
具体地,在进行采样之前,让所述电学信号通过低通滤波器,所述低通滤波器用于将频率高于f1Hh的信号滤除,f1是系统的标称频率,Hh是输入信号包含的最高谐波次数。
具体地,所述电学信号包括电压信号或电流信号。
本实施例中,采样频率为60*32Hz,本发明的瞬时频率获取结果如表一所示:
表一初始迭代频率获取结果
输入频率为测试频率,带有符号“*”为本方法确定的初始迭代频率,比如,第一行输入频率为55Hz,窗口长度为L=35,采用短时傅立叶变换方法的频率估计值为55.1432Hz,这个值是第一行中所有数值与表达式60×32/L差值的最小值。接下来可以选择这个频率估计值为初始迭代频率,显而易见,本方法获得的初始迭代频率与输入频率已经非常接近,使用本方法的初始迭代频率能够获得较快的迭代速度,提升瞬时频率的获取效果。
第二个实施例,如图2所示,一种获取电学信号瞬时频率的系统,包括采样频率设定模块、采样模块、初始迭代频率获取模块和瞬时频率计算模块,
所述采样频率设定模块用于电学信号的初始频率设定采样频率f,
所述采样模块根据采样频率f对电学信号进行采样,并将采样结果传输至初始迭代频率获取模块,
所述初始迭代频率获取模块用于获取初始迭代频率,并将所述初始迭代频率传输至瞬时频率计算模块,由所述瞬时频率计算模块计算出电学信号的瞬时频率。
具体地,所述初始迭代频率获取模块包括频率标定值获取单元,频率计算值获取单元和初始迭代频率获取单元,
所述频率标定值获取单元用于根据采样信号频率和窗口长度获取不同窗口长度下的频率标定值,并将所述频率标定值传输至初始迭代频率获取单元;所述频率计算值获取单元用于根据短时傅里叶变换算法获取不同窗口长度下的频率计算值,并将所述频率计算值传输至初始迭代频率获取单元;所述初始迭代频率获取单元用于根据所述频率标定值和频率计算值获取初始迭代频率。
具体地,所述采样频率设定模块设定的采样频率满足公式,f≥2Hhf1,Hh是所述电学信号中包含的最高谐波次数,f1是系统的标称频率,。
具体地,还包括低通滤波器,所述低通滤波器用于将频率高于Hhf1的信号在进入采样模块之前滤除。
以上所揭露的仅为本发明较佳实施例而已,当然不能以此来限定本发明之权利范围,因此依本发明权利要求所作的等同变化,仍属本发明所涵盖的范围。
Claims (9)
1.一种获取电学信号瞬时频率的方法,其特征在于,包括以下步骤:
S1、预先设定好采样频率,并根据所述采样频率对电学信号进行采样;
S2、获取不同窗口长度下的频率标定值,所述频率标定值fi=f/Li,Li为不同的数据长度,f为采样信号频率,fi为不同数据长度下的频率标定值;
S3、按照短时傅立叶变换算法(SDRT)获取采样时刻的不同窗口长度下的频率计算值;
S4、计算S2和S3中相同数据窗口长度下频率标定值与频率计算值的差异值,确定最小差异值对应的频率计算值;
S5、将最小差异值对应的频率计算值作为采样时刻的初始迭代频率;
S6、根据所述初始迭代频率计算采样时刻的瞬时频率。
2.根据权利要求1所述的一种获取电学信号瞬时频率的方法,其特征在于,S6包括以下步骤:
S61.根据所述初始迭代频率获取所述采样时刻的第一次频率估计值
S62、用所述第一次频率估计值代替所述初始迭代频率,重复步骤S61获取所述采样时刻的第二次频率估计值
S63、按照S62中所述方法依次获取下一次的频率估计值,当获取频率估计值的次数达到预设次数或相邻频率估计值之间的差值不大于预设阈值时,停止获取下一次的频率估计值;
S64、将最后一次获取的频率估计值作为所述采样时刻的瞬时频率值。
3.根据权利要求2所述的一种获取电学信号瞬时频率的方法,其特征在于,S61中包括如下步骤:
S611、获取过渡参数X1和X2,其中A1和是所述电学信号的初始振幅和初始相位;
S612、根据所述过渡参数获取所述第一次频率估计值其中, 若或者有否则
4.根据权利要求1~3中任意一项所述的一种获取电学信号瞬时频率的方法,其特征在于,在进行采样之前,让所述电学信号通过低通滤波器,所述低通滤波器用于将频率高于f1Hh的信号滤除,f1是系统的标称频率,Hh是输入信号包含的最高谐波次数。
5.根据权利要求4所述的一种获取电学信号瞬时频率的方法,其特征在于,所述电学信号包括电压信号或电流信号。
6.一种获取电学信号瞬时频率的系统,其特征在于,包括采样频率设定模块、采样模块、初始迭代频率获取模块和瞬时频率计算模块,
所述采样频率设定模块用于电学信号的初始频率设定采样频率f,
所述采样模块根据采样频率f对电学信号进行采样,并将采样结果传输至初始迭代频率获取模块,
所述初始迭代频率获取模块用于获取初始迭代频率,并将所述初始迭代频率传输至瞬时频率计算模块,由所述瞬时频率计算模块计算出电学信号的瞬时频率。
7.根据权利要求6所述的一种获取电学信号瞬时频率的系统,其特征在于,所述初始迭代频率获取模块包括频率标定值获取单元,频率计算值获取单元和初始迭代频率获取单元,
所述频率标定值获取单元用于根据采样信号频率和窗口长度获取不同窗口长度下的频率标定值,并将所述频率标定值传输至初始迭代频率获取单元;所述频率计算值获取单元用于根据短时傅里叶变换算法获取不同窗口长度下的频率计算值,并将所述频率计算值传输至初始迭代频率获取单元;所述初始迭代频率获取单元用于根据所述频率标定值和频率计算值获取初始迭代频率。
8.根据权利要求6或7所述的一种获取电学信号瞬时频率的系统,其特征在于,所述采样频率设定模块设定的采样频率满足公式,f≥2Hhf1,Hh是所述电学信号中包含的最高谐波次数,f1是系统的标称频率,。
9.根据权利要求8所述的一种获取电学信号瞬时频率的系统,其特征在于,还包括低通滤波器,所述低通滤波器用于将频率高于Hhf1的信号在进入采样模块之前滤除。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610434916.7A CN106053936B (zh) | 2016-06-17 | 2016-06-17 | 一种获取电学信号瞬时频率的方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610434916.7A CN106053936B (zh) | 2016-06-17 | 2016-06-17 | 一种获取电学信号瞬时频率的方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106053936A true CN106053936A (zh) | 2016-10-26 |
CN106053936B CN106053936B (zh) | 2019-05-14 |
Family
ID=57168490
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610434916.7A Active CN106053936B (zh) | 2016-06-17 | 2016-06-17 | 一种获取电学信号瞬时频率的方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106053936B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107092579A (zh) * | 2017-03-20 | 2017-08-25 | 东南大学 | 一种基于ffb改进的sdft频率估计方法 |
CN111896802A (zh) * | 2020-06-30 | 2020-11-06 | 成都四威功率电子科技有限公司 | 一种频率自适应的采样方法 |
CN112505413A (zh) * | 2020-11-25 | 2021-03-16 | 华中科技大学 | 一种时频分析方法和系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4310800A (en) * | 1979-11-07 | 1982-01-12 | General Motors Corporation | Digital frequency measuring circuitry |
CN101603984A (zh) * | 2009-07-01 | 2009-12-16 | 湖南大学 | 电气信号频率的数字化实时检测方法 |
CN101858938A (zh) * | 2009-12-18 | 2010-10-13 | 西安电子科技大学 | 基于自适应滤波原理的瞬时频率测量方法 |
CN104330622A (zh) * | 2014-10-29 | 2015-02-04 | 广东电网有限责任公司电力科学研究院 | 电力系统中正弦波信号的频率测量方法及系统 |
CN104849548A (zh) * | 2015-06-01 | 2015-08-19 | 海南大学 | 一种电力系统瞬时频率监测方法及系统 |
-
2016
- 2016-06-17 CN CN201610434916.7A patent/CN106053936B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4310800A (en) * | 1979-11-07 | 1982-01-12 | General Motors Corporation | Digital frequency measuring circuitry |
CN101603984A (zh) * | 2009-07-01 | 2009-12-16 | 湖南大学 | 电气信号频率的数字化实时检测方法 |
CN101858938A (zh) * | 2009-12-18 | 2010-10-13 | 西安电子科技大学 | 基于自适应滤波原理的瞬时频率测量方法 |
CN104330622A (zh) * | 2014-10-29 | 2015-02-04 | 广东电网有限责任公司电力科学研究院 | 电力系统中正弦波信号的频率测量方法及系统 |
CN104849548A (zh) * | 2015-06-01 | 2015-08-19 | 海南大学 | 一种电力系统瞬时频率监测方法及系统 |
Non-Patent Citations (1)
Title |
---|
郭少琨 等: "瞬时频率估计方法对比研究", 《国外电子测量技术》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107092579A (zh) * | 2017-03-20 | 2017-08-25 | 东南大学 | 一种基于ffb改进的sdft频率估计方法 |
CN111896802A (zh) * | 2020-06-30 | 2020-11-06 | 成都四威功率电子科技有限公司 | 一种频率自适应的采样方法 |
CN112505413A (zh) * | 2020-11-25 | 2021-03-16 | 华中科技大学 | 一种时频分析方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN106053936B (zh) | 2019-05-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110579682B (zh) | 一种故障录波数据暂态同源对比方法及装置 | |
Su et al. | Power harmonic and interharmonic detection method in renewable power based on Nuttall double‐window all‐phase FFT algorithm | |
CN102435844B (zh) | 一种频率无关的正弦信号相量计算方法 | |
CN102338827B (zh) | 用于电网波形畸变分析和电力谐波参数自动监测方法 | |
CN102901874A (zh) | 基于相位与时标测量的电子式互感器绝对延时检测方法 | |
CN106053936A (zh) | 一种获取电学信号瞬时频率的方法及系统 | |
CN106199183A (zh) | 一种实现次同步振荡在线辨识告警的pmu和方法 | |
CN106980044B (zh) | 一种适应风电接入的电力系统谐波电流估计方法 | |
CN103941088A (zh) | 一种基于三相信号的电力系统频率快速测量方法 | |
CN110244151B (zh) | 一种同步相量测量方法、系统及相关组件 | |
CN104635109A (zh) | 一种从故障录波数据中直接计算提取电网暂态非周期分量幅值和时间常数的高精度方法 | |
CN102253282A (zh) | 纳托尔窗函数连续频谱内插电力谐波参数获取方法 | |
CN105224811A (zh) | 一种基于反馈迭代式频率跟踪的pmu动态数据处理方法 | |
CN109557371A (zh) | 一种用于配电网相量测量的同步授时守时方法 | |
CN106092156A (zh) | 交流伺服串行通信编码器位置反馈脉冲分频输出系统和方法 | |
CN105510706A (zh) | 一种高精度欠采样测频方法 | |
CN110417031B (zh) | 一种自动发电控制系统频率偏差系数分段整定方法 | |
CN103995180A (zh) | 一种考虑不等式约束的电力系统频率估计方法 | |
CN104849548B (zh) | 一种电力系统瞬时频率监测方法及系统 | |
Zhuang et al. | Four harmonic analysis and energy metering algorithms based on a new cosine window function | |
CN106324403B (zh) | 一种基于状态识别的继电保护适应性滤波算法 | |
CN116243042A (zh) | 一种配电网的电压暂降检测方法 | |
CN109061299A (zh) | 一种基于径向基函数神经网络的电力系统谐波检测方法 | |
CN113076628B (zh) | 适用于现代大电网频率安全指标的解析方法与系统 | |
CN108362940A (zh) | 一种半周波畸变信号的动态谐波频率提取方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |