全相位FFT算法在谐波测量中的应用概要.docx

全相位FFT算法在谐波测量中的应用概要.docx

《全相位FFT算法在谐波测量中的应用概要.docx》由会员分享，可在线阅读，更多相关《全相位FFT算法在谐波测量中的应用概要.docx（12页珍藏版）》请在冰豆网上搜索。

全相位FFT算法在谐波测量中的应用概要

电测与仪表

ElectricalMeasurement&Instrumentation

总第49卷第554期2012年

第2期

Vol.49No.554

Feb.2012

引言

在电力系统中,由于各种非线性电力设备负荷（如各种晶闸管整流装置、变频装置等因素的影响,电网频率通常会发生偏移,从而产生各次谐波分量,使得电网污染日益严重,这已成为影响电能质量的主要公害之一,对电力系统的安全和经济运行造成了极大的影响[1]

。

所以,实时测量电网中的谐波含量以确切掌握电网谐波的实际状况,对于防止谐波危害、维护电网的安全运行是十分必要的[2]

。

目前常用的最直接的电网谐波参数的分析方法

是应用FFT算法。

为解决不同步采样而造成的频谱泄漏,越来越多的人更关注各种算法来解决这种非同步

采样的问题,文献[3-7]详细的介绍了加窗插值法,并

逐步改进,但此方法不足之处在于对谐波的相角测量精度不高。

频谱校正理论的不断发展与完善,使谐波测量在精度上达到了一定高度,文献[8]中提到的频谱校正方法比较多,许多频谱校正理论方法已经能准确地自动识别出各谐波成分,并自动校正其频率、幅值和相位[9-10]

。

本采用的全相位FFT算法（以下简称apFFT,它具有“相位不变性”,能够自动搜索峰值,无需任何校正措施,具有很优良的抑制频谱泄漏特性,大大的提高了电力系统谐波测量的精度[11]

。

1

全相位分析方法的原理

全相位FFT算法在谐波测量中的应用

付贤东1

康喜明2

卢永杰1

尹星

1

（1.三峡大学电气与新能源学院,湖北宜昌443002;2.青海省电力公司海西供电公司修试所,青海海西816000

摘要:

传统FFT算法的结果可以通过一些算法实现频谱校正,如FFT插值法、比值法等,但在对谐波参数测量时存在一定的误差,精度有限,影响谐波分析结果的准确性。

本文在现有的离散频谱校正方法基础上,提出一种全相位FFT算法。

该算法实现自动搜索各个谐波峰值,能够有效提高对电网谐波频率、相位和幅值的测量精度。

通过matlab仿真工具,将全相位FFT算法与比值法测量结果进行比较,结果表明该算法简单实用,精度高,特别是对间谐波的检测,能有效防止频谱泄露,抗干扰能力强。

关键词:

全相位FFT;数据预处理;相位不变性;间谐波中图分类号:

TM935

文献标识码:

A

文章编号:

1001-1390（201202-0019-04

ApplicationofAllPhaseFFTAlgorithminHarmonicMeasurement

FUXian-dong1

KANGXi-ming2

LUYong-jie1

YINXing

1

（1.CollegeofElectricalEngineering&NewEnergy,ChinaThreeGorgesUniversity,Yichang443002,Hubei,China.

2.InstituteofOverhauling&Experiment,HaiXiElectricPowerCompany,Haixi816000,Qinghai,ChinaAbstract:

TheresultofthetraditionalFFTalgorithmcanbecorrectedthroughsomemethods,suchastheFFTwindowinterpolationalgorithm,theratiomethodandsoon.However,thetestaccuracyislimitedandthereisonelargermeasurementerrorthataffecttheaccuracyofharmonicanalysis.Basedonthemethodofdiscretespectrumcorrection,thispaperstudiesanewcorrectionmethod-AllphaseFFTalgorithm.Thisalgorithmcanautomaticallysearchforthepeakofeachharmonic,anditcaneffectivelyimprovethepoweroftheharmonicfrequency,phaseandamplitudeaccuracy.Inthispaper,withthematlabtool,comparingthetwomethods,theresultsshowthattheAllphaseFFTalgorithmissimple,practicalandhighprecision,especiallyforthedetectionofinter-harmonics.Itcaneffectivelypreventspectralleakageandhasastronganti-interferenceability.

Keywords:

all-phaseFFT,pre-processingofthedata,phaseinvariant,inter-harmonics19--

电测与仪表

ElectricalMeasurement&Instrumentation

总第49卷第554期2012年第2期

Vol.49No.554

Feb.2012

全相位分析方法的根本思想是为解决信号数据截断所产生误差而提出的,该方法主要考虑了数据段中心样本点的所有可能截断组合,因此它能有效的减小栅栏效应和泄漏现象。

假设并设原信号的离散序列为x（n,输入数据为{…,x（-N+1,…,x（0,…,x（N-1,…},若把所有的分段都考虑,则可形成如下分段x0~xN-1,有:

x0:

x（0,…,x（N-2,x（N-1x1:

x（-1,x（0,…,x（N-3,x（N-2……

xN-1:

x（-N+1,…,x（-1,x（0

若将各分段数据以x0为中心对齐,这样对信号进行截断,就充分考虑到了所有情况,再对所有情况做FFT变换,最后将所有结果利用窗函数加权求和,具体过程如图1所示。

在如上所示的一组数据中信号的离散序列为:

x（n=ej（2πk0n/N+φ0

（1

则DFT为:

XN（k=12π

X（ejω·W（ejωω=

2πkN

=1

2π

2πδ（ω-ω0X（ejωejφ0

]·sin（ωN/2

sin（ω/2e

-j1

2

n-1ω

ω=

2πkN

=

sin（π（k-k0

sin（π（k-k0/N

e

-j（N-1（k-k0π/Ne

jφ0

（2

式中k=0,1,…,（N-1。

由式（2可知（k-k0为频率偏离值,所以DFT相位与频率偏离值有关,将上述信号

经全相位处理后输出,根据文献[9]中推导过程知apFFT输出为:

YN（k=1

Nsin2

（π（k-k0

sin2（π（k-k0/N

e

jφ0

（3

从公式知apFFT相位为初相位φ0,与频率偏离值无关,它具有“相位不变性”,所以apFFT的相位-频率

特性是一条水平线,即其群延时为0,在任一频率点测相位都可以。

2apFFT的数据处理过程

apFFT的算法是先对数据进行全相位处理,然后进行传统的FFT运算。

经全相位处理后的信号,可由原来的带有跳变的信号变为平滑连续的信号,使得周期延拓后不会出现信号的跳变。

本文中数据处理如下:

（1构成一个N点的汉宁窗;

（2汉宁窗对自己求卷积,得到（2N-1点的卷积窗;（3求（2N-1点的卷积窗的和;

（4将卷积窗的每一项除以卷积窗的和,得到（2N-1点的归一化卷积窗;

（5将数据的1:

（2N-1项和归一化卷积窗相乘,得到加窗的（2N-1项;延时N,将数据的（N+1:

（3N-1项和归一化卷积窗相乘,同样得到加窗的（2N-1项;（6将第1项和（N+1项,第2项和（N+2项...第（N-1

项和第（2N-1项相加,将延时N后的数据按同样的方式相加,得到经过全相预处理的N点序列。

实现过程框图如图2所示。

3实验仿真

为了验证算法的准确性,假定电网波形可用表达为:

x（t=∑n

Ancos2πnf0t+φn∑∑（4

其中,基波、谐波、间谐波的幅值按实际电网信号（这里选择公网谐波电压的限制范围内,在电网标称

电压为0.38kV的情况下

的特点设置,各种频率信号的相位任意设定,具体参数如表1所示。

为了能更清楚的说明此方法的优越性,本文将文献[8-9]所讲的比值法与本方法在相同的条件下通过

图1N阶apFFT频谱分析过程基本框图Fig.1N-orderapFFTspectrumanalysisofthe

basicblockdiagram

x（-N+1x（-1x（0x（1x（N-1图2apFFT法系统结构图

Fig.2

SystemstructurediagramoftheapFFTmethod

20--

电测与仪表

ElectricalMeasurement&Instrumentation

总第49卷第554期2012年

第2期

Vol.49No.554

Feb.2012

matlab仿真比较,设采样频率为3000Hz,数据长度为1024,两种方法均采用汉宁加窗,在无噪和加噪（信噪

比为40dB情况下,得到各次谐波、间谐波幅值、相位、频率与给定值的比较。

具体的数据如表2所示。

4仿真结果分析

表2中的所有实验结果可以看到:

（1在无噪时,apFFT法测的频率、相位和幅值基本上与给定值一致,最大的相角偏差仅为0.0002°,比值法测的最大相角为1.58614°;在相角测量精度方面,apFFT法甚至比值法

高出6个数量级（如在间谐波频率为25Hz时,

apFFT相位测量误差仅为0.00001°,而比值法的测量误差为1.58614°;（2在加噪时,apFFT法测量的频率、相位和幅值精度比传统比值法基本高出1~2个数量级;（3从分析偶次谐波（k=2,4

时,比值法估计精度有所下降（如偶次谐波相位偏差达到了0.5°左右,而apFFT则不然,其频率、相位和幅值基本上无偏差。

表1

各谐波参数

Tab.1

Theparametersoftheharmonic

表2谐波实测结果分析表

Tab.2

Measuredresultsofharmonicanalysistable

分析测试结果,我们不难得出如下结论:

（1从频率、幅值和相位三个参数来看,无论是在无噪还是加噪情况下,用apFFT法测得的电网谐波的精度都高于比值法,其中在测试相角精度上apFFT优势明显;

（2对间谐波参数的测量,加噪的情况下apFFT法对频率、相角和幅值的测量精度依然很高;（3apFFT法减少了频率间的干涉现现象,对间谐波的测量精度高;

（4apFFT法抗噪能力强,具有实用优势。

5

结束语

本文提出的apFFT算法具有较高的谐波参数估计精度,尤其是对间谐波参数的准确测量,这对谐波补偿装置的设计具有十分重要的意义。

另外两个电网要并网,需要电压信号同频同相同幅,测出同一时刻的两路的样点精确相位是首要的,apFFT法无疑在其中可发挥很大的优势。

不足之处是当信号中包含噪声时,如何提高谐波参数测量精度还值得作进一步的研究。

/Hz/V/（°251.86105038051000.622016517.5302000.936025013.340

315

1.55

80

/Hz

2550100165200250315apFFT24.9999950.00000

99.99999164.99999200.00000

249.99999314.99999-0.000010-0.00001-0.000010-0.00001-0.0000125.0023049.9974899.99999164.99999200.00000249.99999314.999990.00230-0.00252-0.00001-0.000010.00000-0.00001-0.00001apFFT（

25.0003949.99999100.00156164.99999199.99922249.99984314.99984

0.00039-0.000010.00156-0.00001-0.00078-0.00016-0.00016

25.0023849.9974899.97053164.99904199.98791249.99898314.99958

0.00238

-0.00252

-0.02947

-0.00096

-0.01219

-0.00102

-0.00042

/（°1052030604080apFFT9.999985.0000019.9999930.00000

59.9999840.00000

79.99999-0.000020-0.000010-0.000020-0.000018.413865.1427321.5419030.0036160.5328140.0025780.10298-1.586141.142730.54190.002610.532810.002570.10298apFFT（

9.937545.0003319.8924730.0001859.9568240.0098080.03958

-0.062460.00033-0.107530.00018-0.043180.009800.03958

8.044465.1429821.5939230.0073560.6358840.0058980.109970.044460.142981.59392

0.0007350.635880.005890.10997

/V1.86380

0.6217.50.9313.31.55apFFT1.86001379.999990.6200017.499990.9300013.299991.549990.00001-0.000010.00000-0.000010.00000-0.00001-0.000011.86121380.341830.6230317.512010.9317913.309011.551020.00121

0.341830.003030.012010.001790.009010.00102apFFT（

1.85688380.000770.6216117.500370.9297713.299971.54993-0.003120.000770.001610.000370.00977-0.00003-0.000071.86187

380.342530.6226417.512110.9318613.309221.551450.001870.342530.002640.012110.001860.009220.00145

21--

电测与仪表

ElectricalMeasurement&Instrumentation

总第49卷第554期2012年第2期

Vol.49No.554

Feb.2012

参考文献

[1]吴国乔,王兆华.离散频谱的全相位校正法[J].数据采集与处理,2005,20（3:

286-290.

WUGuo-qiao,WANGZhao-hua.AllPhaseCorrectionMethodforDiscreteSpectrum[J].JournalofDataAcquisition&Processing,2005,20（3:

286-290.

[2]肖雁鸿,毛筱.电力系统谐波测量方法综述[J].电网技术,2002,26（6:

61-64.

XIAOYan-hong,MAOXiao.ASurveyonMeasuringMethodforHarmonicsinPowerSystem[J].PowerSystemTechnolog,2002,26（6:

61-64.

[3]张建秋,沈毅.非整周期采样算法的误差分析[J].哈尔滨工业大学学报,1995,27（5:

99-103.

ZHANGJian-qiu,SHENYi.TheAnalysisoftheNon-Integer-PeriodSamplingAlgorithm[J].JournalofHarbinInstituteofTechnology,1995,27（5:

99-103.

[4]潘文,钱俞寿.基于加窗差值FFT的电力谐波测量理论（II双插值FFT理论[J]电工技术学报,1994,2（5:

53-56.

PanWen,QianYu-shou.PowerHarmonicsMeasurementBasedonWindowsandInterpolatedFFT（IIDualInterpolatedFFTAlgorithms[J].TransactionsofChinaElectrotechnicalSociety,1994,2（5:

53-56.[5]任祖华.基于窗函数的FFT谐波参数估计算法[J].电测与仪表,2010,47（5:

8-11.

RENZu-hua.TheHarmonicParameterEstimationAlgorithmBasedontheFFTWindowFunctions[J].ElectricalMeasurement&Instrumentation,2010,47（5:

8-11.

[6]李朝生,周进有,李先允.应用改进FFT算法实现风电场谐波的检测[J].电测与仪表,2010,47（11:

12-16.

LIChao-sheng,ZHOUJin-you,LIXian-yun.UsinganImprovedFFTAlgorithmtoDetectHarmonicinaWindPowerPlant[J].ElectricalMeasurement&Instrumentation,2010,47（11:

12-16.

[7]庞浩,李东霞.应用FFT进行电力系统谐波分析的改进算法[J].中国电机工程学报,2003,23（6:

50-54.

PANGHao,LIDong-xia.AnImprovedAlgorithmforHarmonicAnalysisofPowerSystemUsingFFTTechnique[J].ProceedingsoftheCSEE,2003,23（6:

50-54.[8]丁康,张晓飞.频谱校正理论的发展[J].振动工程学报,2002,13（1:

14-22.

DingKang,ZHANGXiao-fei.AdvancesinSpectrumCorrectionTheory

[J].JournalofVibrationEngineering.2002,13（1:

14-22.

[9]黄翔东.全相位数字信号处理[D].天津大学博士生论文,2006.12.1.HUANGXiang-dong.All-PhaseDigitalSignalProcessing[D].TianjinUniversitydoctoralthesis,2006.12.1.

[10]谢明,丁康.频谱分析的校正方法[J].振动工程学报,1994;7

（2:

172-179.

XIEMing,DINGKang.RectifyingTechniqueofSpectrumAnalysis[J].JournalofVibrationEngineering,1994;7（2:

172-179.

[11]王兆华,候正信,苏飞.全相位FFT频谱分析[J].通信学报,2003,24

（11A:

6-19.

WANGZhao-hua,HOUZheng-xin,SUFei.AllPhaseSpectrumAnalysis[J].JournalofChinaInstituteofCommunications,2003,24（11A:

6-19.作者简介:

收稿日期:

2011-07-02

（田春雨编发

付贤东（1984-,男,湖北武汉人,硕士研究生,研究方向为电力电子技术。

Email:

fxdaky@

康喜明（1984-,男,内蒙古呼和浩特和林格尔县人,青海省海西供电公司职员,助理工程师。

卢永杰（1984-,男,陕西西安人,硕士研究生,研究方向为电力电子及高压直流输电技术。

尹星（1987-,女,江苏泰州人,硕士在读,主要从事电力系统谐波的分析计算的研究。

22--