CN106053936B - 一种获取电学信号瞬时频率的方法及*** - Google Patents
一种获取电学信号瞬时频率的方法及*** Download PDFInfo
- Publication number
- CN106053936B CN106053936B CN201610434916.7A CN201610434916A CN106053936B CN 106053936 B CN106053936 B CN 106053936B CN 201610434916 A CN201610434916 A CN 201610434916A CN 106053936 B CN106053936 B CN 106053936B
- Authority
- CN
- China
- Prior art keywords
- frequency
- sampling
- electrical signal
- value
- obtaining
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000005070 sampling Methods 0.000 claims description 71
- 238000004364 calculation method Methods 0.000 claims description 28
- 230000007704 transition Effects 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 3
- 238000009499 grossing Methods 0.000 abstract 1
- 238000005259 measurement Methods 0.000 description 4
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
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.根据权利要求1所述的一种获取电学信号瞬时频率的方法,其特征在于,还包括一种获取电学信号瞬时频率的***,包括采样频率设定模块、采样模块、初始迭代频率获取模块和瞬时频率计算模块,
所述采样频率设定模块用于电学信号的初始频率设定采样频率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 CN106053936A (zh) | 2016-10-26 |
CN106053936B true CN106053936B (zh) | 2019-05-14 |
Family
ID=57168490
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610434916.7A Expired - Fee Related CN106053936B (zh) | 2016-06-17 | 2016-06-17 | 一种获取电学信号瞬时频率的方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106053936B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107092579A (zh) * | 2017-03-20 | 2017-08-25 | 东南大学 | 一种基于ffb改进的sdft频率估计方法 |
CN111896802B (zh) * | 2020-06-30 | 2023-03-31 | 成都四威功率电子科技有限公司 | 一种频率自适应的采样方法 |
CN112505413B (zh) * | 2020-11-25 | 2021-09-14 | 华中科技大学 | 一种时频分析方法和*** |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4310800A (en) * | 1979-11-07 | 1982-01-12 | General Motors Corporation | Digital frequency measuring circuitry |
CN101603984B (zh) * | 2009-07-01 | 2011-04-06 | 湖南大学 | 电气信号频率的数字化实时检测方法 |
CN101858938A (zh) * | 2009-12-18 | 2010-10-13 | 西安电子科技大学 | 基于自适应滤波原理的瞬时频率测量方法 |
CN104330622B (zh) * | 2014-10-29 | 2017-06-09 | 广东电网有限责任公司电力科学研究院 | 电力***中正弦波信号的频率测量方法及*** |
CN104849548B (zh) * | 2015-06-01 | 2018-05-25 | 海南大学 | 一种电力***瞬时频率监测方法及*** |
-
2016
- 2016-06-17 CN CN201610434916.7A patent/CN106053936B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN106053936A (zh) | 2016-10-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chen et al. | Comparative study of harmonic and interharmonic estimation methods for stationary and time-varying signals | |
Su et al. | Power harmonic and interharmonic detection method in renewable power based on Nuttall double‐window all‐phase FFT algorithm | |
Zhang et al. | Analysis of white noise on power frequency estimation by DFT-based frequency shifting and filtering algorithm | |
WO2016138764A1 (zh) | 一种提升pmu同步相量测量精度的方法 | |
CN108414833B (zh) | 一种信号分量频率的精确估计方法 | |
Zhang et al. | Improved smoothing frequency shifting and filtering algorithm for harmonic analysis with systematic error compensation | |
CN106018956B (zh) | 一种加窗谱线插值的电力***频率计算方法 | |
CN106053936B (zh) | 一种获取电学信号瞬时频率的方法及*** | |
Chen et al. | An efficient Prony-based solution procedure for tracking of power system voltage variations | |
Chang et al. | Measurement techniques for stationary and time-varying harmonics | |
CN111984920B (zh) | 次/超同步谐波参数识别方法、装置、设备和介质 | |
CN106501602B (zh) | 一种基于滑窗频谱分离的基波参数测量方法 | |
CN103941088A (zh) | 一种基于三相信号的电力***频率快速测量方法 | |
CN110967658A (zh) | 一种基于数字微差法的模拟量输入合并单元校验仪溯源的方法 | |
WO2018090783A1 (zh) | 电机定子电阻的在线辨识方法、装置和电机控制*** | |
Ma et al. | Harmonic and interharmonic analysis of mixed dense frequency signals | |
CN103543331B (zh) | 一种计算电信号谐波和间谐波的方法 | |
CN109557365B (zh) | 电网电压频率和相角检测方法及***、单相电网检测*** | |
Zhang et al. | Frequency shifting and filtering algorithm for power system harmonic estimation | |
Arabaci et al. | A novel motor speed calculation method using square wave speed sensor signals via fast Fourier transform | |
CN105044459A (zh) | 一种谐波分析方法 | |
CN104849548B (zh) | 一种电力***瞬时频率监测方法及*** | |
CN104330644A (zh) | 检测电力***中正弦波信号幅值的方法和*** | |
CN105137198A (zh) | 一种基于Nuttall窗-五点变换FFT的介损测量新方法 | |
KR101997633B1 (ko) | Teo 및 desa를 이용한 자동 동기화 파라미터 측정 장치 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190514 |
|
CF01 | Termination of patent right due to non-payment of annual fee |