CN101915874B - 一种基于傅立叶变换的谐波检测方法 - Google Patents

一种基于傅立叶变换的谐波检测方法 Download PDF

Info

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
sampling
data
period
Prior art date
Application number
CN 201010238721
Other languages
English (en)
Other versions
CN101915874A (zh
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 北海市深蓝科技发展有限责任公司
Priority to CN 201010238721 priority Critical patent/CN101915874B/zh
Publication of CN101915874A publication Critical patent/CN101915874A/zh
Application granted granted Critical
Publication of CN101915874B publication Critical patent/CN101915874B/zh

Links

Abstract

本发明公开了一种基于傅立叶变换的谐波检测方法,本方法包括以下几方面:(1)对被测信号进行过采样并进行模数转换,按所测信号估计的最大周期值,保证每组数据能采样到两个信号周期的数据;(2)用数字低通滤波器滤除基波以外的谐波成分;(3)用周期法求取基波周期;(4)一周期内均匀提取2*N(N-1为最高次谐波阶数)个数据点;(5)用同步采样FFT求取各谐波成分的参数。属于信号处理技术领域,也可用于其它信号的频谱分析。用该方法可得到非常高的谐波检测精度,实时性与准同步法的三到五个信号周期相比还能提高一倍左右,只需两个信号周期的时间,并且还能精确地测量出基波的频率。

Description

一种基于傅立叶变换的谐波检测方法
技术领域:
[0001] 本发明涉及一种用于电力系统谐波检测和基波频率测量的方法,具体是一种基于同步采样FFT (快速傅立叶变换)的谐波参数检测方法,属于信号处理技术领域,也可用于其它信号的频谱分析。
背景技术:
[0002] 电力系统谐波检测的主要方法有:基于模拟带通或带阻滤波器的谐波检测、基于瞬时无功功率的谐波检测、基于傅立叶变换的谐波检测、基于神经网络的谐波检测、基于小波变换的谐波检测。
[0003] 基于模拟带通或带阻滤波器的谐波检测是早期的模拟谐波测量方法。基于瞬时无 功功率理论的分析方法在解决谐波总量实时测量方面很有优势,但不能解决各次谐波成分含量的检测问题。
[0004] 小波变换(Wavelet Transformation, WT)分析方法对波动谐波、快速变化的谐波检测有很大优势;但是WT并不能完全取代傅里叶变换,这是因为:一方面WT在稳态谐波检测方面并不具备理论优势,另一方面WT的理论和应用研究时间相对较短,WT应用在谐波测量方面尚处于初始阶段,还存在着许多不完善的地方,例如缺乏系统规范的最佳小波基的选取方法,缺乏构造频域行为良好(分频严格、能量集中)的小波函数以改善检测精度的规范方法。
[0005] 神经网络(Neural Network, NN)在谐波检测中则具有计算量小,精度高,实时性好,抗干扰性好;但是NN的检测方法用于工程实际还有很多问题,例如:没有规范的NN构造方法,需要大量的训练样本,如何确定需要的样本数没有规范方法,NN的精度对样本有很大的依赖性。
[0006] 另外基于小波变换和神经网络的谐波检测方法应用于工程实际中还存在着许多问题,都属于正在研究的新方法,研究和应用时间短,实现技术尚需完善,目前在工程应用中还未优先选用。
[0007] 而现在应用最广泛的谐波检测方法是基于傅立叶变换的谐波检测方法。在稳态的谐波检测中,基于傅立叶变换的谐波检测方法具有精度较高,功能较多,使用方便,易于实现等优点,但它的实时性还有待改进。所以在新的谐波检测理论还没足够成熟时,研究改进基于傅立叶变换的谐波检测方法的实时性还有着很重要的意义。
[0008] 造成FFT谐波检测误差的根源是由于采样不同步引起的频谱泄漏。解决频谱泄漏的办法主要有两类:一是同步误差一定的情况下,通过对采样数据的处理或测量结果的修正来减少测量误差,如准同步算法、加窗插值算法和准同步采样补偿法等;二是通过减少同步误差来减少测量误差,如双速率采样法、优化采样周期法等。
[0009] 基于傅立叶变换的谐波检测方法中最经典的就是准同步采样法,准同步采样法通过迭代运算,在采样周期和信号周期不严格同步的情况下,仍能实现较高精度的谐波分析。不需要使用信号周期值,对采样起点无任何要求,只要求频率是相对稳定的,以牺牲时间换取精度,一般需要3至5个信号周期。当信号频率波动较大造成同步误差较大时,可通过增加迭代次数来提高测量精度。因此准同步采样法除了实时性相对较差外确有不少优点。
[0010] 加窗插值算法可得到较高的测量精度,但它的实时性比准同步采样法还差,一般需要十个信号周期左右。
[0011] 准同步采样补偿法、双速率采样法、采样周期优化法、非整周期采样法等这些方法要的数据仅为一个周期左右的数据,但是在测量前需要精确知道信号的基波周期。而在实际的测量中,基波周期可能会变化,测量前是不能精确知道。因此要使用这些方法就必须解决好基波周期的精确测量问题。
[0012] 通过对现有方法的研究和利用现在DSP中AD转换的快速性,提出了一个实时性和精确度都较好的谐波检测方案。先用较高频率对信号进行过采样,再用低通滤波的方法精确测量出基波周期,然后再用直线似合的方法在一个整周期的时间内均匀求出2*N(N-1为最高次谐波阶数)个数据点,最后直接用FFT求出各谐波成分的参数,整个过程仅需两个信 号周期。
发明内容:
[0013] 所要解决的技术问题
[0014] 针对准同步采样法实时性不理想这一不足,提出在保证较高的测量精度的同时提高实时性的解决方案。
[0015] 技术方案(发明概述)
[0016] 本发明的技术方案主要如下:
[0017] I、对被测信号进行过采样并进行模数转换,按所测信号估计的最大周期值,保证每组数据能采样到两个信号周期的数据。
[0018] 2、设计IIR数字滤波器对数据进行滤波,将基波以外的谐波成分全部滤除。
[0019] 在对数据进行滤波前,先用数字滤波器的采样频率对原始数据进行重采样,使用一个固定的频率进行重采样是为了使滤波器在设计时容易获得稳定系数和在滤波时能在较短的时间内稳定下来。
[0020] 由于对数据进行滤波只是为了求基波周期,所以对滤波器的幅频特性和相频特性就没有太高的要求,只要滤波后基波的周期不变就可以了。最经典的数字滤波器为FIR滤波器和IIR滤波器,但在相同的性能指标下,FIR滤波器的阶数要比IIR滤波器高得多,为了减小计算量,这里选择IIR滤波器。
[0021] 用MATLAB语言对IIR滤波器进行辅助设计,求出指定指标的滤波器系数。这里用椭圆滤波器模型设计IIR数字滤波器,因为它的阶数最小且容易获得稳定的滤波系统。
[0022] 在设计滤波器时,为了使滤波器一直工作在较稳定的状态下,这里不是对每组数据单独进行滤波,而是把上一组数据滤波前后的M(M为滤波器长度)数据作为下一组的初始值,这样就相当于只是对一组很长的数据进行较稳定的滤波,避免了每组数据都需要等待一个较长的过渡期。
[0023] 3、用周期法对滤波后的数据求取基波周期值。
[0024] 为了减小上一组数据对下一组数据的影响,这里选择用第四个过零点减去第二个过零点来求取周期。第二个过零点前的时间为过渡时间,保证有一定时间的过渡期,使这种基波频率求取方法有更强的适应性。
[0025] 4、对原始数据在一个基波周期的时间内均匀提取2*N(N_1为最高次谐波阶数)个点。
[0026] 对起始点没有特别要求,只要保证其后有一个完整周期的数据就可以。当要提取的点不在原始数据上,就用直线拟合的方法对原始数据相邻的两点进行拟合。由于原始数据的采样频率非常高,所以用直线拟合的方法来求取仍能得到非常高的精度。用这个方法求取2*N个数据点,避免了软件同步采样时截断误差造成周期误差,减小了最终的频谱泄漏。
[0027] 5、用同步采样法的傅立叶变换的谐波检测方法求取各谐波成分的值。
[0028] 由于上面的2*N个点的数据是在一个整周期的时间里均匀求取的,它几乎达到理想的同步采样,因此用同步采样法来求取谐波参数引起的频谱泄漏就非常小。
[0029] 有益效果
[0030] 与准同步采样法相比,这种方法在满足较高精确度要求的同时,实时性与准同步的三到五个信号周期相比提高一倍左右。这种方法可以做到每两个信号周期更新一次数据,且能非常精确地测量出基波的周期。
附图说明
[0031] 图I是本发明的谐波检测流程图;
具体实施方式
[0032] 本发明提出了一种基于傅立叶变换谐波检测的新方法,以下结合实例详述,但不作为本发明的限定。本实例中,最高次谐波为63次,主要参数选择如下:
[0033] I、采样频率为:Fs = 50*128*40 ;
[0034] 2、电网基波频率:在47Hz到53Hz之间波动;
[0035] 3、每组数据的采样时间:0. 044375s ;
[0036] 4、滤波器指标:ffp = 60Hz, Ws = 135Hz, Rp = O. 5db, Rs = 50db ;
[0037] 5、滤波器采样频率:50*128*2Hz ;
[0038] 6、FFT 运算数据长度:2*N = 128。
[0039] 本实例的处理流程框图如图I所示,信号经过A/D转换(模数转换)后送入DSP进行处理,滤波后求基波周期,在一个信号周期内均匀提取128个点,最后再进行FFT运算,得到各次谐波参数。
[0040] 本实例中,用MATLAB进行滤波器辅助设计时求得的系数为:
[0041] a= [1.0000 -3.9642 5.8940 -3.8955 0.9657]
[0042] b = [O. 0031 -O. 0124 O. 0186 -O. 0124 O. 0031]
[0043] 滤波器的差分方程为:
[0044] y (n) = b (I) *x (n)+b (2) *χ (η_1)+· · ·+b (5) *x (n_5)-a (2) *y (η_1)-· · ·-a (5) *y (η-4)
[0045] 数据经过滤波后得到的波形只剩下基波成分,这样就可以用周期法来求取基波的周期了,且得到非常高的精确度。[0046] 本实例中,信号基频取48Hz、50Hz和52Hz三个频率,被测信号的相关参数值和测量结果如下:
[0047] I、基波频率为52Hz检测结果
[0048] 表I基波频率检测结果
[0049]
Figure CN101915874BD00061
[0050] 表2谐波幅值和初相位检测结果 [0051]
Figure CN101915874BD00062
[0053] 表3基波频率检测结果
[0054]
Figure CN101915874BD00063
[0055] 表4谐波幅值和初相位检测结果
[0056]
Figure CN101915874BD00071
[0057] 3、基波频率为48Hz检测结果
[0058] 表5基波频率检测结果
[0059]
Figure CN101915874BD00072
[0060] 表6谐波幅值和初相位检测结果
[0061]
Figure CN101915874BD00073
[0062] 结果表明,使用这种方法可以得到非常高的測量精确度,且只要约两个信号周期的时间,实时性也很好,还能非常精确地测量出基波的周期。

Claims (5)

1. 一种基于傅立叶变换的谐波检测方法,其特征是:首先对被测信号进行过采样,再用数字低通滤波的方法求取信号周期,然后在一个周期的时间内均匀提取2*N个数据点,其中N-I为最高次谐波阶数,最后用同步采样FFT求出各谐波参数。
2.根据权利要求 I所述的基于傅立叶变换的谐波检测方法,过采样是采用较高频率进行采样,采样频率=F*2*N*K,其中F是信号的基波频率,N-I为最高次谐波阶数,K是采样倍率,K要求大于10,使相邻两点用直线拟合的方法求中间值时的误差很小,每组数据采样要保证大于两个信号基波周期。
3.根据权利要求I所述的基于傅立叶变换的谐波检测方法,数字低通滤波是指能将基波以外的谐波成分都滤除的滤波算法。
4.根据权利要求I所述的基于傅立叶变换的谐波检测方法,滤波求周期是指对只剩基波成分的波形用周期法求取基波周期。
5.根据权利要求I所述的基于傅立叶变换的谐波检测方法,均匀提取2*N个点是指在一个整周期的时间里用直线拟合插值的方法均匀求取2*N个数据点的值,N-I为最高次谐波阶数。
CN 201010238721 2010-07-20 2010-07-20 一种基于傅立叶变换的谐波检测方法 CN101915874B (zh)

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 CN101915874B (zh) 2010-07-20 2010-07-20 一种基于傅立叶变换的谐波检测方法

Country Status (1)

Country Link
CN (1) CN101915874B (zh)

Families Citing this family (17)

* Cited by examiner, † Cited by third party
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 中国矿业大学(北京) 一种电力系统谐波的快速检测方法
CN109030938A (zh) * 2017-06-08 2018-12-18 许继集团有限公司 一种基于正弦滤波的抗谐波测频方法和装置

Citations (3)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0915273A (ja) * 1995-06-30 1997-01-17 Sony Tektronix Corp 基本波成分周波数の測定方法及び測定装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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
Lobos et al. Real-time determination of power system frequency
CN104749498B (zh) 一种便携式局部放电检测及诊断装置
Routray et al. A novel Kalman filter for frequency estimation of distorted signals in power systems
CN103558436B (zh) 基于单相锁相环算法的检测电网电压幅值、频率和相角的方法
CN104698273B (zh) 谐波责任划分方法和系统
CN104502700B (zh) 电力信号的正弦参数测量方法和系统
CN101149425B (zh) 一种电子式电能表调试、校验自动化系统
de la O Serna et al. Instantaneous oscillating phasor estimates with Taylor $^ K $-Kalman filters
CN102650658B (zh) 一种时变非平稳信号时频分析方法
CN103630871B (zh) 电子式电能表动态特性测试装置
CN103323702B (zh) 复合电能质量扰动信号识别方法
CN102565542B (zh) 基于iec61850-9-2标准的容性设备介质损耗在线监测方法
CN102305891B (zh) 一种电力系统低频振荡在线监测方法
CN201993463U (zh) 数字化电能表检测装置
CN106841778B (zh) 基于pmu实现的次同步和超同步谐波参数的处理方法
CN103257280A (zh) 基于电场传感器的容性设备介损监测装置及监测方法
CN103401238B (zh) 一种基于总体测辨法的电力负荷建模方法
CN103353550A (zh) 一种测量电力系统信号频率及谐波参数的方法
CN103245832B (zh) 基于快速s变换的谐波时频特性参数估计方法及分析仪
CN107271768B (zh) 一种最小二乘拟合动态频率测量方法
CN101034117A (zh) 输电线路工频参数测量方法及装置
CN100557453C (zh) 基于相量测量技术的发电机一次调频实时定量分析方法
CN102565574A (zh) 船舶电能质量监测系统
CN104007346B (zh) 基于频率扫描的模拟量合并单元暂态延时测试方法
CN105046025B (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