CN101701984B - 基于三项系数Nuttall窗插值FFT的基波与谐波检测方法 - Google Patents
基于三项系数Nuttall窗插值FFT的基波与谐波检测方法 Download PDFInfo
- Publication number
- CN101701984B CN101701984B CN2009101546816A CN200910154681A CN101701984B CN 101701984 B CN101701984 B CN 101701984B CN 2009101546816 A CN2009101546816 A CN 2009101546816A CN 200910154681 A CN200910154681 A CN 200910154681A CN 101701984 B CN101701984 B CN 101701984B
- Authority
- CN
- China
- Prior art keywords
- harmonic
- coefficient
- interpolation
- fft
- amplitude
- 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 26
- 230000003595 spectral effect Effects 0.000 claims abstract description 27
- 238000001228 spectrum Methods 0.000 claims description 5
- ODKSFYDXXFIFQN-UHFFFAOYSA-M argininate Chemical compound [O-]C(=O)C(N)CCCNC(N)=N ODKSFYDXXFIFQN-UHFFFAOYSA-M 0.000 claims description 3
- 238000005070 sampling Methods 0.000 abstract description 3
- 238000009472 formulation Methods 0.000 abstract 2
- 239000000203 mixture Substances 0.000 abstract 2
- 238000001514 detection method Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Landscapes
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
本发明公开了一种基于三项系数Nuttall窗插值FFT的基波与谐波检测方法。通过模数转换器得到被测信号的采样数据;对采样数据进行加三项系数Nuttall窗FFT运算;对基波和各次谐波在FFT谱线中搜索幅值谱线的最大和相邻次大谱线,根据相邻谱峰的幅值之比按显式表达式直接计算基波和各次谐波的插值系数;最后通过插值运算得到基波和各次谐波的频率、幅值和相位。本发明首次提出基于显式表达式直接计算插值系数的三项系数Nuttall窗FFT基波与谐波检测方法。本方案的优点是提供一种计算量小、分析精度高的加窗插值FFT基波与谐波检测方法。
Description
技术领域
本发明涉及一种信号中基波与谐波检测领域,尤其涉及一种基于三项系数Nuttall窗插值FFT的基波与谐波检测方法。
背景技术
以电力信号检测为例,电力***谐波影响电力设备的正常运行,对基波和谐波进行准确检测是电网和电力设备安全稳定运行的基本要求。加窗插值FFT(快速傅里叶变换)方法是电力***谐波检测的常用方法。在非同步采样时,傅里叶变换存在频谱泄漏和栅栏效应。加窗插值FFT通过加窗运算抑制频谱泄漏,通过插值运算消除栅栏效应的影响。加窗插值FFT方法常用的余弦窗函数有Hanning窗、Blackman窗、Blackman-Harris窗、Rife-Vincent窗、最优余弦窗等。
已有专利文件“基于Nuttall窗双峰插值FFT的基波与谐波检测方法”(200710035653.3),其发明的目的在于提供基于Nuttall窗的双峰插值FFT方法,该方法应用四项三阶Nuttall窗,采用双峰插值算法,不具有显式直接计算的插值系数计算公式,计算量大。
已有技术的不足之处是,基于项数少的窗函数的加窗插值FFT方法计算量小,具有显式的插值系数计算公式,但分析精度低;基于项数多的窗函数的加窗插值FFT方法分析精度高,但需要求解多次方程,无法通过显式表达式直接计算插值系数,计算量大。例如,Hanning窗为两项系数余弦窗,其插值公式为显式的计算公式,简单明了,可以直接计算,计算量小,计算稳定性好,但分析精度较低;Blackman-Harris窗和四项三阶Nuttall窗函数的项数都多于Hanning窗,分析精度较高,但在插值系数计算时需要求解多次方程或多项式拟合逼近,无法直接计算插值系数,计算量大,在噪声和干扰下可能导致方程的有效解不存在。
发明内容
本发明的目的针对上述技术的不足,提供一种基于三项系数Nuttall窗插值FFT的基波与谐波检测方法,它具有显式的插值系数计算公式,计算量小,计算稳定性好,基波和谐波的检测精度高。
基于三项系数Nuttall窗插值FFT的基波与谐波检测方法包括如下步骤:
1)通过模数转换器以采样频率fs将被测电压电流信号从模拟信号转化为数字信号,得到N点长度的采样数据;
2)构造N点长度的三项系数Nuttall窗函数,对N点长度的采样数据进行加三项系数Nuttall窗FFT运算,得到FFT谱线X(k),k=0,1,...,N,其中FFT运算的数据长度为N;
3)对基波和各次谐波在FFT谱线中搜索幅值谱线的最大和相邻次大谱线,根据相邻谱峰的幅值之比按显式表达式直接计算基波和各次谐波的插值系数;
4)通过插值运算得到基波和各次谐波的频率、幅值和相位。
上述步骤2)中,三项系数Nuttall窗函数的构造方法为:
其中a0=0.375,a1=0.5,a2=0.125。
上述步骤3)中,计算基波和各次谐波插值系数的显式表达式为:
式中,δm为第m次谐波的插值系数, 为第m次谐波相邻最大和次大谱线的幅值之比。
上述步骤4)中,计算基波和各次谐波的频率、幅值和相位的插值运算公式为:
fm=(km+δm)fs/N (3)
本发明首次提出基于显式表达式直接计算插值系数的三项系数Nuttall窗FFT基波与谐波检测方法,减小了计算量,提高了运算稳定性,提高了加窗插值FFT的分析精度。本发明的优点是:1、本发明的计算量小,运算稳定性好,由于采用了显式的插值系数计算,无需求解多次方程,在噪声和干扰下也有很好的适用性,易于单片机或嵌入式***实现;2、本发明的检测精度高,由于三项系数Nuttall窗的最大旁瓣为-47dB,每个旁瓣的衰减速度为30dB,因此基波和谐波的检测精度高。
附图说明
图1为基于三项系数Nuttall窗插值FFT的基波与谐波检测方法的框图。
具体实施方式
基于三项系数Nuttall窗插值FFT的基波与谐波检测方法包括如下步骤:
1)通过模数转换器以采样频率fs将被测电压电流信号从模拟信号转化为数字信号,得到N点长度的采样数据;
2)构造N点长度的三项系数Nuttall窗函数,对N点长度的采样数据进行加三项系数Nuttall窗FFT运算,得到FFT谱线X(k),k=0,1,...,N,其中FFT运算的数据长度为N;
3)对基波和各次谐波在FFT谱线中搜索幅值谱线的最大和相邻次大谱线,根据相邻谱峰的幅值之比按显式表达式直接计算基波和各次谐波的插值系数;
4)通过插值运算得到基波和各次谐波的频率、幅值和相位。
上述步骤2)中,三项系数Nuttall窗函数的构造方法为:
其中a0=0.375,a1=0.5,a2=0.125。
上述步骤3)中,计算基波和各次谐波插值系数的显式表达式为:
式中,δm为第m次谐波的插值系数, 为第m次谐波相邻最大和次大谱线的幅值之比。
上述步骤4)中,计算基波和各次谐波的频率、幅值和相位的插值运算公式为:
fm=(km+δm)fs/N (3)
实施例1
以某电器工作时的电流谐波检测为例,设该电器的电流表达式为其设定值如表1所示。应用本发明基于三项系数Nuttall窗插值FFT的基波与谐波检测方法测量其基波和2~9次谐波(但本发明不局限于2~9次谐波),本实施例中通过Matlab仿真软件来说明其实施过程。
(1)通过模数转换器得到该电器电流的采样数据,其中模数转换器采用美信MAX125CEAX集成电路芯片,采样频率fs=10kHz,数据长度N=2048。Matlab软件在原始信号中增加80dB的高斯白噪声以表示其测量噪声。
(2)构造2048点三项系数Nuttall窗,对采样数据加窗,并进行2048点FFT变换,得到1024个FFT谱线,记为X(0)、X(1)、...、X(1023)。
(3)对基波和各次谐波在谱线FFT中搜索幅值谱线的最大和相邻次大谱线,得到:k1=10、k2=20、k3=30、k4=40、k5=51、k6=61、k7=71、k8=81、k9=91;根据式(2)计算基波和各次谐波的插值系数,得到:δ1=0.2195、δ2=0.4398、δ3=0.6585、δ4=0.8777、δ5=0.0977、δ6=0.3167、δ7=0.5363、δ8=0.7565、δ9=0.9757。
(4)最后通过插值运算公式(3)、(4)和(5)得到基波和各次谐波的频率、幅值和相位,结果如表1所示。
表1 实施例检测结果
综上所述,本发明基于三项系数Nuttall窗插值FFT的基波与谐波检测方法可以直接计算基波和各次谐波的插值系数,计算量小,运算稳定性好,基波和谐波的检测精度高。
Claims (1)
1.一种基于三项系数Nuttall窗插值FFT的基波与谐波检测方法,其特征在于包括如下步骤:
1)通过模数转换器以采样频率fs将被测电压电流信号从模拟信号转化为数字信号,得到N点长度的采样数据;
2)构造N点长度的三项系数Nuttall窗函数,对N点长度的采样数据进行加三项系数Nuttall窗FFT运算,得到FFT谱线X(k),k=0,1,...,N,其中FFT运算的数据长度为N;
3)对基波和各次谐波在FFT谱线中搜索幅值谱线的最大和相邻次大谱线,根据相邻谱峰的幅值之比按显式表达式直接计算基波和各次谐波的插值系数;
4)通过插值运算得到基波和各次谐波的频率、幅值和相位;
步骤2)中,所述的三项系数Nuttall窗函数的构造方法为:
其中a0=0.375,a1=0.5,a2=0.125;
步骤3)中,所述的计算基波和各次谐波插值系数的显式表达式为:
步骤4)中,所述的基波和各次谐波的频率、幅值和相位的插值运算公式为:
fm=(km+δm)fs/N (3)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009101546816A CN101701984B (zh) | 2009-11-23 | 2009-11-23 | 基于三项系数Nuttall窗插值FFT的基波与谐波检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009101546816A CN101701984B (zh) | 2009-11-23 | 2009-11-23 | 基于三项系数Nuttall窗插值FFT的基波与谐波检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101701984A CN101701984A (zh) | 2010-05-05 |
CN101701984B true CN101701984B (zh) | 2011-05-18 |
Family
ID=42156904
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009101546816A Expired - Fee Related CN101701984B (zh) | 2009-11-23 | 2009-11-23 | 基于三项系数Nuttall窗插值FFT的基波与谐波检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101701984B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102243272A (zh) * | 2011-04-01 | 2011-11-16 | 重庆大学 | 一种高精度采样数据同步的谐波分析方法 |
CN102331526B (zh) * | 2011-06-10 | 2013-07-17 | 中国矿业大学 | 汉宁窗函数连续频谱内插电力谐波参数获取方法 |
CN102253282B (zh) * | 2011-06-10 | 2014-07-16 | 中国矿业大学 | 纳托尔窗函数连续频谱内插电力谐波参数获取方法 |
CN104062528A (zh) * | 2014-07-04 | 2014-09-24 | 武汉大学 | 基于汉宁乘积窗的信号谐波分析方法及*** |
CN104391178B (zh) * | 2014-12-05 | 2017-12-15 | 国家电网公司 | 一种基于Nuttall窗的时移相位差稳态谐波信号校正方法 |
CN104749433A (zh) * | 2015-01-04 | 2015-07-01 | 云南电网公司电力科学研究院 | 一种准确计算分析谐波及电能质量的方法 |
CN104833937B (zh) * | 2015-05-21 | 2017-08-11 | 湖南大学 | 一种基于mir‑rsd高精度余弦窗插值fft算法的谐波测量通道校准方法 |
CN105137183A (zh) * | 2015-09-15 | 2015-12-09 | 湖北工业大学 | 一种电力***谐波分析方法及*** |
CN105137181A (zh) * | 2015-09-15 | 2015-12-09 | 福州大学 | 基于Nuttall-Kaiser组合窗双谱线插值的谐波分析方法 |
CN105353215A (zh) * | 2015-12-23 | 2016-02-24 | 合肥工业大学 | 基于Nuttall窗四谱线插值FFT的谐波检测方法 |
CN105717359A (zh) * | 2016-02-19 | 2016-06-29 | 云南电网有限责任公司电力科学研究院 | 一种谐波分析算法 |
CN107271002A (zh) * | 2017-06-19 | 2017-10-20 | 重庆邮电大学 | 一种快速高精度的频谱校正插值算法 |
CN109782063A (zh) * | 2018-10-23 | 2019-05-21 | 国网安徽省电力有限公司芜湖供电公司 | 一种基于纳托尔自卷积窗三谱线插值fft的动态间谐波分析方法 |
CN110007148B (zh) * | 2019-03-28 | 2021-03-16 | 东南大学 | 一种基于离散频谱相位和幅值综合内插的单频信号频率估计方法 |
CN114659791B (zh) * | 2022-02-28 | 2023-07-04 | 广东机电职业技术学院 | 汽轮机故障检测方法、***、装置和存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6014620A (en) * | 1995-06-21 | 2000-01-11 | Telefonaktiebolaget Lm Ericsson | Power spectral density estimation method and apparatus using LPC analysis |
CN101113995A (zh) * | 2007-08-29 | 2008-01-30 | 湖南大学 | 基于Nuttall窗双峰插值FFT的基波与谐波检测方法 |
CN101261618A (zh) * | 2008-04-14 | 2008-09-10 | 湖南大学 | 三角自卷积窗加权算法 |
-
2009
- 2009-11-23 CN CN2009101546816A patent/CN101701984B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6014620A (en) * | 1995-06-21 | 2000-01-11 | Telefonaktiebolaget Lm Ericsson | Power spectral density estimation method and apparatus using LPC analysis |
CN101113995A (zh) * | 2007-08-29 | 2008-01-30 | 湖南大学 | 基于Nuttall窗双峰插值FFT的基波与谐波检测方法 |
CN101261618A (zh) * | 2008-04-14 | 2008-09-10 | 湖南大学 | 三角自卷积窗加权算法 |
Non-Patent Citations (2)
Title |
---|
张介秋.基于卷积窗的电力***谐波误差估计与数值模拟.《中国电机工程学报》.2004,第24卷(第12期),34-37. * |
李明.时频原子方法在间谐波分析中的应用.《电网技术》.2009,第33卷(第17 期),81-85. * |
Also Published As
Publication number | Publication date |
---|---|
CN101701984A (zh) | 2010-05-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101701984B (zh) | 基于三项系数Nuttall窗插值FFT的基波与谐波检测方法 | |
CN101701982B (zh) | 基于加窗插值fft的电力***谐波检测方法 | |
CN103454497B (zh) | 基于改进加窗离散傅立叶变换的相位差测量方法 | |
CN101603985B (zh) | 高准确度正弦信号测量方法 | |
CN106483374B (zh) | 一种基于Nuttall双窗全相位FFT的谐波间谐波检测方法 | |
CN103245832B (zh) | 基于快速s变换的谐波时频特性参数估计方法及分析仪 | |
CN102288807B (zh) | 一种测量电网电压闪变的方法 | |
CN101718816B (zh) | 基于四项系数Nuttall窗插值FFT的基波与谐波检测方法 | |
CN103308766A (zh) | 一种基于凯撒自卷积窗双谱线插值fft谐波分析方法及其装置 | |
CN103575984A (zh) | 基于凯塞窗双谱线插值fft的谐波分析方法 | |
CN101261292A (zh) | 基于5项Rife-Vincent(I)窗双谱线插值FFT的基波与谐波检测方法 | |
CN109782063A (zh) | 一种基于纳托尔自卷积窗三谱线插值fft的动态间谐波分析方法 | |
CN103197141A (zh) | 一种测量电力***信号频率及谐波参数的方法 | |
CN104714075B (zh) | 一种电网电压闪变包络参数提取方法 | |
CN109946512B (zh) | 一种改进频域插值的动态功率分析方法 | |
CN103399204A (zh) | 一种基于Rife-Vincent(II)窗插值FFT的谐波与间谐波检测方法 | |
CN110133564A (zh) | 一种电流传感器相位特性的同步测试方法及装置 | |
CN105137181A (zh) | 基于Nuttall-Kaiser组合窗双谱线插值的谐波分析方法 | |
CN102590598B (zh) | 基于多阈值比较的超声波信号过零点预测方法 | |
CN103795411A (zh) | 基于五项最大旁瓣衰减窗三谱线插值测试sfdr的方法 | |
CN102998528B (zh) | 动态正弦畸变信号中频率分量的提取方法及其实现装置 | |
CN108776263A (zh) | 基于高阶汉宁自卷积窗及改进插值算法的谐波检测方法 | |
CN103197143A (zh) | 基于汉宁窗fft算法与遍历滤波的谐波、间谐波检测方法 | |
CN202102051U (zh) | 一种电网间谐波测量仪 | |
CN103543331A (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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110518 Termination date: 20111123 |