CN101915874B - 一种基于傅立叶变换的谐波检测方法 - Google Patents
一种基于傅立叶变换的谐波检测方法 Download PDFInfo
- Publication number
- CN101915874B CN101915874B CN 201010238721 CN201010238721A CN101915874B CN 101915874 B CN101915874 B CN 101915874B CN 201010238721 CN201010238721 CN 201010238721 CN 201010238721 A CN201010238721 A CN 201010238721A CN 101915874 B CN101915874 B CN 101915874B
- Authority
- CN
- China
- Prior art keywords
- harmonic
- signal
- method based
- sampling
- period
- 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
- 238000001514 detection method Methods 0.000 title abstract description 9
- 230000009466 transformation Effects 0.000 title abstract description 6
- 238000000034 method Methods 0.000 claims abstract description 60
- 230000001360 synchronised effect Effects 0.000 claims abstract description 20
- 238000001914 filtration Methods 0.000 claims abstract description 19
- 230000000737 periodic effect Effects 0.000 claims abstract description 5
- 238000005070 sampling Methods 0.000 claims description 28
- 239000000284 extract Substances 0.000 claims description 4
- 239000000203 mixture Substances 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 abstract description 5
- 238000012545 processing Methods 0.000 abstract description 3
- 238000010183 spectrum analysis Methods 0.000 abstract description 2
- 238000005259 measurement Methods 0.000 description 12
- 238000012360 testing method Methods 0.000 description 10
- 238000013528 artificial neural network Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000001228 spectrum Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 238000011960 computer-aided design Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000009977 dual effect Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000035800 maturation Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
Images
Landscapes
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
本发明公开了一种基于傅立叶变换的谐波检测方法，本方法包括以下几方面：(1)对被测信号进行过采样并进行模数转换，按所测信号估计的最大周期值，保证每组数据能采样到两个信号周期的数据；(2)用数字低通滤波器滤除基波以外的谐波成分；(3)用周期法求取基波周期；(4)一周期内均匀提取2*N(N-1为最高次谐波阶数)个数据点；(5)用同步采样FFT求取各谐波成分的参数。属于信号处理技术领域，也可用于其它信号的频谱分析。用该方法可得到非常高的谐波检测精度，实时性与准同步法的三到五个信号周期相比还能提高一倍左右，只需两个信号周期的时间，并且还能精确地测量出基波的频率。
Description
技术领域：
本发明涉及一种用于电力系统谐波检测和基波频率测量的方法，具体是一种基于同步采样FFT(快速傅立叶变换)的谐波参数检测方法，属于信号处理技术领域，也可用于其它信号的频谱分析。
背景技术：
电力系统谐波检测的主要方法有：基于模拟带通或带阻滤波器的谐波检测、基于瞬时无功功率的谐波检测、基于傅立叶变换的谐波检测、基于神经网络的谐波检测、基于小波变换的谐波检测。
基于模拟带通或带阻滤波器的谐波检测是早期的模拟谐波测量方法。基于瞬时无功功率理论的分析方法在解决谐波总量实时测量方面很有优势，但不能解决各次谐波成分含量的检测问题。
小波变换(Wavelet Transformation，WT)分析方法对波动谐波、快速变化的谐波检测有很大优势；但是WT并不能完全取代傅里叶变换，这是因为：一方面WT在稳态谐波检测方面并不具备理论优势，另一方面WT的理论和应用研究时间相对较短，WT应用在谐波测量方面尚处于初始阶段，还存在着许多不完善的地方，例如缺乏系统规范的最佳小波基的选取方法，缺乏构造频域行为良好(分频严格、能量集中)的小波函数以改善检测精度的规范方法。
神经网络(Neural Network，NN)在谐波检测中则具有计算量小，精度高，实时性好，抗干扰性好；但是NN的检测方法用于工程实际还有很多问题，例如：没有规范的NN构造方法，需要大量的训练样本，如何确定需要的样本数没有规范方法，NN的精度对样本有很大的依赖性。
另外基于小波变换和神经网络的谐波检测方法应用于工程实际中还存在着许多问题，都属于正在研究的新方法，研究和应用时间短，实现技术尚需完善，目前在工程应用中还未优先选用。
而现在应用最广泛的谐波检测方法是基于傅立叶变换的谐波检测方法。在稳态的谐波检测中，基于傅立叶变换的谐波检测方法具有精度较高，功能较多，使用方便，易于实现等优点，但它的实时性还有待改进。所以在新的谐波检测理论还没足够成熟时，研究改进基于傅立叶变换的谐波检测方法的实时性还有着很重要的意义。
造成FFT谐波检测误差的根源是由于采样不同步引起的频谱泄漏。解决频谱泄漏的办法主要有两类：一是同步误差一定的情况下，通过对采样数据的处理或测量结果的修正来减少测量误差，如准同步算法、加窗插值算法和准同步采样补偿法等；二是通过减少同步误差来减少测量误差，如双速率采样法、优化采样周期法等。
基于傅立叶变换的谐波检测方法中最经典的就是准同步采样法，准同步采样法通过迭代运算，在采样周期和信号周期不严格同步的情况下，仍能实现较高精度的谐波分析。不需要使用信号周期值，对采样起点无任何要求，只要求频率是相对稳定的，以牺牲时间换取精度，一般需要3至5个信号周期。当信号频率波动较大造成同步误差较大时，可通过增加迭代次数来提高测量精度。因此准同步采样法除了实时性相对较差外确有不少优点。
加窗插值算法可得到较高的测量精度，但它的实时性比准同步采样法还差，一般需要十个信号周期左右。
准同步采样补偿法、双速率采样法、采样周期优化法、非整周期采样法等这些方法要的数据仅为一个周期左右的数据，但是在测量前需要精确知道信号的基波周期。而在实际的测量中，基波周期可能会变化，测量前是不能精确知道。因此要使用这些方法就必须解决好基波周期的精确测量问题。
通过对现有方法的研究和利用现在DSP中AD转换的快速性，提出了一个实时性和精确度都较好的谐波检测方案。先用较高频率对信号进行过采样，再用低通滤波的方法精确测量出基波周期，然后再用直线似合的方法在一个整周期的时间内均匀求出2*N(N-1为最高次谐波阶数)个数据点，最后直接用FFT求出各谐波成分的参数，整个过程仅需两个信号周期。
发明内容：
所要解决的技术问题
针对准同步采样法实时性不理想这一不足，提出在保证较高的测量精度的同时提高实时性的解决方案。
技术方案(发明概述)
本发明的技术方案主要如下：
1、对被测信号进行过采样并进行模数转换，按所测信号估计的最大周期值，保证每组数据能采样到两个信号周期的数据。
2、设计IIR数字滤波器对数据进行滤波，将基波以外的谐波成分全部滤除。
在对数据进行滤波前，先用数字滤波器的采样频率对原始数据进行重采样，使用一个固定的频率进行重采样是为了使滤波器在设计时容易获得稳定系数和在滤波时能在较短的时间内稳定下来。
由于对数据进行滤波只是为了求基波周期，所以对滤波器的幅频特性和相频特性就没有太高的要求，只要滤波后基波的周期不变就可以了。最经典的数字滤波器为FIR滤波器和IIR滤波器，但在相同的性能指标下，FIR滤波器的阶数要比IIR滤波器高得多，为了减小计算量，这里选择IIR滤波器。
用MATLAB语言对IIR滤波器进行辅助设计，求出指定指标的滤波器系数。这里用椭圆滤波器模型设计IIR数字滤波器，因为它的阶数最小且容易获得稳定的滤波系统。
在设计滤波器时，为了使滤波器一直工作在较稳定的状态下，这里不是对每组数据单独进行滤波，而是把上一组数据滤波前后的M(M为滤波器长度)数据作为下一组的初始值，这样就相当于只是对一组很长的数据进行较稳定的滤波，避免了每组数据都需要等待一个较长的过渡期。
3、用周期法对滤波后的数据求取基波周期值。
为了减小上一组数据对下一组数据的影响，这里选择用第四个过零点减去第二个过零点来求取周期。第二个过零点前的时间为过渡时间，保证有一定时间的过渡期，使这种基波频率求取方法有更强的适应性。
4、对原始数据在一个基波周期的时间内均匀提取2*N(N-1为最高次谐波阶数)个点。
对起始点没有特别要求，只要保证其后有一个完整周期的数据就可以。当要提取的点不在原始数据上，就用直线拟合的方法对原始数据相邻的两点进行拟合。由于原始数据的采样频率非常高，所以用直线拟合的方法来求取仍能得到非常高的精度。用这个方法求取2*N个数据点，避免了软件同步采样时截断误差造成周期误差，减小了最终的频谱泄漏。
5、用同步采样法的傅立叶变换的谐波检测方法求取各谐波成分的值。
由于上面的2*N个点的数据是在一个整周期的时间里均匀求取的，它几乎达到理想的同步采样，因此用同步采样法来求取谐波参数引起的频谱泄漏就非常小。
有益效果
与准同步采样法相比，这种方法在满足较高精确度要求的同时，实时性与准同步的三到五个信号周期相比提高一倍左右。这种方法可以做到每两个信号周期更新一次数据，且能非常精确地测量出基波的周期。
附图说明
图1是本发明的谐波检测流程图；
具体实施方式
本发明提出了一种基于傅立叶变换谐波检测的新方法，以下结合实例详述，但不作为本发明的限定。本实例中，最高次谐波为63次，主要参数选择如下：
1、采样频率为：Fs＝50*128*40；
2、电网基波频率：在47Hz到53Hz之间波动；
3、每组数据的采样时间：0.044375s；
4、滤波器指标：Wp＝60Hz，Ws＝135Hz，Rp＝0.5db，Rs＝50db；
5、滤波器采样频率：50*128*2Hz；
6、FFT运算数据长度：2*N＝128。
本实例的处理流程框图如图1所示，信号经过A/D转换(模数转换)后送入DSP进行处理，滤波后求基波周期，在一个信号周期内均匀提取128个点，最后再进行FFT运算，得到各次谐波参数。
本实例中，用MATLAB进行滤波器辅助设计时求得的系数为：
a＝[1.0000 -3.9642 5.8940 -3.8955 0.9657]
b＝[0.0031 -0.0124 0.0186 -0.0124 0.0031]
滤波器的差分方程为：
y(n)＝b(1)*x(n)+b(2)*x(n-1)+...+b(5)*x(n-5)-a(2)*y(n-1)-...-a(5)*y(n-4)
数据经过滤波后得到的波形只剩下基波成分，这样就可以用周期法来求取基波的周期了，且得到非常高的精确度。
本实例中，信号基频取48Hz、50Hz和52Hz三个频率，被测信号的相关参数值和测量结果如下：
1、基波频率为52Hz检测结果
表1基波频率检测结果
| 实际基波频率 | 测量基波频率 | 误差(％) |
| 52.00000000000000 | 51.99998673122966 | 0.00002551686604 |
表2谐波幅值和初相位检测结果
2、基波频率为50Hz检测结果
表3基波频率检测结果
| 实际基波频率 | 测量基波频率 | 误差(％) |
| 50.00000000000000 | 50.00000142435692 | 0.00000284871385 |
表4谐波幅值和初相位检测结果
3、基波频率为48Hz检测结果
表5基波频率检测结果
| 实际基波频率 | 测量基波频率 | 误差(％) |
| 48.00000000000000 | 47.99999189650186 | 0.00001688228779 |
表6谐波幅值和初相位检测结果
结果表明，使用这种方法可以得到非常高的测量精确度，且只要约两个信号周期的时间，实时性也很好，还能非常精确地测量出基波的周期。
Claims (5)
1.一种基于傅立叶变换的谐波检测方法，其特征是：首先对被测信号进行过采样，再用数字低通滤波的方法求取信号周期，然后在一个周期的时间内均匀提取2*N个数据点，其中N-1为最高次谐波阶数，最后用同步采样FFT求出各谐波参数。
2.根据权利要求1所述的基于傅立叶变换的谐波检测方法，过采样是采用较高频率进行采样，采样频率＝F*2*N*K，其中F是信号的基波频率，N-1为最高次谐波阶数，K是采样倍率，K要求大于10，使相邻两点用直线拟合的方法求中间值时的误差很小，每组数据采样要保证大于两个信号基波周期。
3.根据权利要求1所述的基于傅立叶变换的谐波检测方法，数字低通滤波是指能将基波以外的谐波成分都滤除的滤波算法。
4.根据权利要求1所述的基于傅立叶变换的谐波检测方法，滤波求周期是指对只剩基波成分的波形用周期法求取基波周期。
5.根据权利要求1所述的基于傅立叶变换的谐波检测方法，均匀提取2*N个点是指在一个整周期的时间里用直线拟合插值的方法均匀求取2*N个数据点的值，N-1为最高次谐波阶数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201010238721 CN101915874B (zh) | 2010-07-20 | 2010-07-20 | 一种基于傅立叶变换的谐波检测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201010238721 CN101915874B (zh) | 2010-07-20 | 2010-07-20 | 一种基于傅立叶变换的谐波检测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101915874A CN101915874A (zh) | 2010-12-15 |
| CN101915874B true CN101915874B (zh) | 2013-01-09 |
Family
ID=43323435
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN 201010238721 Expired - Fee Related CN101915874B (zh) | 2010-07-20 | 2010-07-20 | 一种基于傅立叶变换的谐波检测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101915874B (zh) |
Families Citing this family (24)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103439567B (zh) * | 2011-10-21 | 2016-01-27 | 常州顺创电气科技有限公司 | 一种moa阻性电流在线监测系统的工作方法 |
| CN103267895B (zh) * | 2013-04-25 | 2015-04-22 | 国家电网公司 | 一种电网电流的谐波电流检测方法 |
| CN104251934B (zh) * | 2013-06-26 | 2018-08-14 | 华为技术有限公司 | 谐波分析方法和装置以及确定谐波间杂波的方法和装置 |
| CN103454495B (zh) * | 2013-09-13 | 2016-01-20 | 电子科技大学 | 自适应高精度快速频谱分析方法 |
| CN103576120B (zh) * | 2013-11-07 | 2016-09-28 | 国家电网公司 | 三次谐波分量准同步信息传输校验与自愈算法 |
| CN103760425A (zh) * | 2014-01-22 | 2014-04-30 | 湖南大学 | 一种基于时域准同步的介损角快速测量方法及其测量装置 |
| CN103869162B (zh) * | 2014-03-05 | 2017-02-08 | 湖南大学 | 一种基于时域准同步的动态信号相量测量方法 |
| CN103983849B (zh) * | 2014-05-07 | 2016-06-29 | 江苏天浩达科技有限公司 | 一种实时高精度的电力谐波分析方法 |
| CN104502676B (zh) * | 2015-01-09 | 2017-08-11 | 珠海派诺科技股份有限公司 | 一种无限脉冲响应的高精度瞬时信号检测方法 |
| CN104777356A (zh) * | 2015-03-10 | 2015-07-15 | 三峡大学 | 一种基于神经网络的实时高精度谐波检测方法 |
| CN105044459B (zh) * | 2015-07-21 | 2017-09-29 | 青岛艾诺智能仪器有限公司 | 一种谐波分析方法 |
| CN105067881A (zh) * | 2015-07-24 | 2015-11-18 | 淄博格瑞恩电气科技有限公司 | 一种基于谐波检测的用电安全监控方法 |
| CN105137183A (zh) * | 2015-09-15 | 2015-12-09 | 湖北工业大学 | 一种电力系统谐波分析方法及系统 |
| CN105388361B (zh) * | 2015-12-31 | 2018-01-23 | 武汉大学 | 双向插值同步化采样序列的fft电力谐波检测方法 |
| CN107085144B (zh) * | 2017-04-28 | 2019-08-20 | 珠海泰芯半导体有限公司 | 一种快速测量电力系统谐波的方法 |
| CN107271773B (zh) * | 2017-05-24 | 2020-09-04 | 中国矿业大学(北京) | 一种电力系统谐波的快速检测方法 |
| CN109030938B (zh) * | 2017-06-08 | 2021-05-11 | 许继集团有限公司 | 一种基于正弦滤波的抗谐波测频方法和装置 |
| CN110568254A (zh) * | 2019-09-27 | 2019-12-13 | 宁夏凯晨电气集团有限公司 | 一种准确检测故障电流中衰减直流分量参数的方法 |
| CN111007559B (zh) * | 2019-11-25 | 2023-09-15 | 中国辐射防护研究院 | 一种基于Hilbert-Huang变换的γ能谱平滑方法 |
| CN112291039B (zh) * | 2020-10-15 | 2022-03-22 | 天津大学 | 一种基于信道信号周期的数据编码传输方法 |
| CN113341226B (zh) * | 2021-06-21 | 2022-04-29 | 合肥美的暖通设备有限公司 | 谐波检测方法、装置、变频器及存储介质 |
| CN113533846B (zh) * | 2021-07-13 | 2024-08-06 | 国网江苏省电力有限公司营销服务中心 | 一种用电负荷数据高速同步采集计算方法和装置 |
| CN113608023B (zh) * | 2021-08-18 | 2024-08-20 | 国网湖南省电力有限公司 | 采样频率自适应的谐波检测方法及系统 |
| CN117169590B (zh) * | 2023-08-15 | 2024-03-29 | 嘉兴市科讯电子有限公司 | 一种基于软件变采样率的电力谐波分析的方法和装置 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1580791A (zh) * | 2004-05-20 | 2005-02-16 | 上海交通大学 | 基于自适应神经元在线整定的谐波检测方法 |
| CN101403774A (zh) * | 2008-11-07 | 2009-04-08 | 扬州中凌高科技发展有限公司 | 一种基于非同步采样的谐波分析方法 |
| JP2009264753A (ja) * | 2008-04-22 | 2009-11-12 | Yokogawa Electric Corp | 高調波成分測定装置 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0915273A (ja) * | 1995-06-30 | 1997-01-17 | Sony Tektronix Corp | 基本波成分周波数の測定方法及び測定装置 |
-
2010
- 2010-07-20 CN CN 201010238721 patent/CN101915874B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1580791A (zh) * | 2004-05-20 | 2005-02-16 | 上海交通大学 | 基于自适应神经元在线整定的谐波检测方法 |
| JP2009264753A (ja) * | 2008-04-22 | 2009-11-12 | Yokogawa Electric Corp | 高調波成分測定装置 |
| CN101403774A (zh) * | 2008-11-07 | 2009-04-08 | 扬州中凌高科技发展有限公司 | 一种基于非同步采样的谐波分析方法 |
Non-Patent Citations (2)
| Title |
|---|
| 任洪林等.基于滤波技术的电力网络基波频率测量方法.《电测与仪表》.2007,第44卷(第498期), * |
| 梅永 等.电力系统信号采集与谐波测量方法.《电测与仪表》.2008,第45卷(第513期), * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101915874A (zh) | 2010-12-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101915874B (zh) | 一种基于傅立叶变换的谐波检测方法 | |
| CN102288807B (zh) | 一种测量电网电压闪变的方法 | |
| CN103454497B (zh) | 基于改进加窗离散傅立叶变换的相位差测量方法 | |
| CN103869162B (zh) | 一种基于时域准同步的动态信号相量测量方法 | |
| CN102305891A (zh) | 一种电力系统低频振荡在线监测方法 | |
| CN107543962B (zh) | 主导间谐波频谱分布的计算方法 | |
| CN203287435U (zh) | 一种基于stm32f107vct6的微电网谐波与间谐波检测装置 | |
| CN105548739B (zh) | 一种避雷器运行状态信号处理方法 | |
| CN102435844A (zh) | 一种频率无关的正弦信号相量计算方法 | |
| CN102426416B (zh) | 控制对象频域幅值特性辩识方法和装置 | |
| CN103913631A (zh) | 用于数显表的电网频率测量方法和测量系统 | |
| CN103983849B (zh) | 一种实时高精度的电力谐波分析方法 | |
| CN102435849A (zh) | 使用数字滤波器提高频率测量精度的方法 | |
| CN109142863B (zh) | 一种电力系统测频方法及系统 | |
| CN103163372A (zh) | 采用hht的电力系统谐波分析方法 | |
| CN103529294A (zh) | 基于hht的光伏系统并网逆变器谐波检测系统及方法 | |
| CN108918964A (zh) | 一种稀疏性增强的谐波分析方法 | |
| CN109831184A (zh) | 一种带通滤波方法、系统、存储介质及终端设备 | |
| CN102253281A (zh) | 一种基于软件抗谐波干扰并自动换相测量电网频率的方法 | |
| CN102495250A (zh) | 一种基于Hilbert算法的准同步宽频无功电能表及其采样方法 | |
| CN112213560A (zh) | 一种基于z-adaline的高精度电网宽频信号测量方法 | |
| CN103969508A (zh) | 一种实时高精密的电力谐波分析方法及装置 | |
| CN101718816B (zh) | 基于四项系数Nuttall窗插值FFT的基波与谐波检测方法 | |
| CN104111373A (zh) | 一种用于智能变电站数字化电能的计量方法 | |
| CN104569675A (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 | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130109 |