CN101718816A - 基于四项系数Nuttall窗插值FFT的基波与谐波检测方法 - Google Patents
基于四项系数Nuttall窗插值FFT的基波与谐波检测方法 Download PDFInfo
- Publication number
- CN101718816A CN101718816A CN200910154682A CN200910154682A CN101718816A CN 101718816 A CN101718816 A CN 101718816A CN 200910154682 A CN200910154682 A CN 200910154682A CN 200910154682 A CN200910154682 A CN 200910154682A CN 101718816 A CN101718816 A CN 101718816A
- Authority
- CN
- China
- Prior art keywords
- harmonic
- interpolation
- coefficient
- 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.)
- Granted
Links
- 238000001514 detection method Methods 0.000 title abstract description 9
- 238000001228 spectrum Methods 0.000 claims abstract description 10
- 238000000034 method Methods 0.000 claims description 24
- 230000003595 spectral effect Effects 0.000 claims description 22
- 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
- 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
Images
Landscapes
- Complex Calculations (AREA)
- 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.3125,a1=0.46875,a2=0.1875,a3=0.03125。
上述步骤3)中,计算基波和各次谐波插值系数的显式表达式为:
式中,δm为第m次谐波的插值系数,为第m次谐波相邻最大和次大谱线的幅值之比。
上述步骤4)中,计算基波和各次谐波的频率、幅值和相位的插值运算公式为:
fm=(km+δm)fs/N (3)
本发明首次提出基于显式表达式直接计算插值系数的四项系数Nuttall窗FFT基波与谐波检测方法,减小了计算量,提高了运算稳定性,提高了加窗插值FFT的分析精度。本发明的优点是:1、本发明的计算量小,运算稳定性好,由于采用了显式的插值系数计算,无需求解多次方程,在噪声和干扰下也有很好的适用性,易于单片机或嵌入式系统实现;2、本发明的检测精度高,由于四项系数Nuttall窗的最大旁瓣为-61dB,每个旁瓣的衰减速度为42dB,因此基波和谐波的检测精度高。
附图说明
图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.3125,a1=0.46875,a2=0.1875,a3=0.03125。
上述步骤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.4390、δ3=0.6586、δ4=0.8781、δ5=0.0976、δ6=0.3171、δ7=0.5366、δ8=0.7562、δ9=0.9757。
(4)最后通过插值运算公式(3)、(4)和(5)得到基波和各次谐波的频率、幅值和相位,结果如表1所示。
表1实施例检测结果
综上所述,本发明基于四项系数Nuttall窗插值FFT的基波与谐波检测方法可以直接计算基波和各次谐波的插值系数,计算量较小,运算稳定性好,基波和谐波的检测精度高。
Claims (4)
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.根据权利要求1所述的方法,其特征在于:步骤2)中,所述的四项系数Nuttall窗函数的构造方法为:
其中a0=0.3125,a1=0.46875,a2=0.1875,a3=0.03125。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009101546820A CN101718816B (zh) | 2009-11-23 | 2009-11-23 | 基于四项系数Nuttall窗插值FFT的基波与谐波检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009101546820A CN101718816B (zh) | 2009-11-23 | 2009-11-23 | 基于四项系数Nuttall窗插值FFT的基波与谐波检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101718816A true CN101718816A (zh) | 2010-06-02 |
CN101718816B CN101718816B (zh) | 2011-05-18 |
Family
ID=42433420
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009101546820A Expired - Fee Related CN101718816B (zh) | 2009-11-23 | 2009-11-23 | 基于四项系数Nuttall窗插值FFT的基波与谐波检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101718816B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102033164A (zh) * | 2010-11-16 | 2011-04-27 | 哈尔滨工业大学 | 一种计算电信号的基波分量采样信号序列的方法和系统 |
CN102253282A (zh) * | 2011-06-10 | 2011-11-23 | 中国矿业大学 | 纳托尔窗函数连续频谱内插电力谐波参数获取方法 |
CN102305886A (zh) * | 2011-05-31 | 2012-01-04 | 浙江大学 | 电网电压谐波畸变及不平衡时基波电压同步信号检测方法 |
CN103454495A (zh) * | 2013-09-13 | 2013-12-18 | 电子科技大学 | 自适应高精度快速频谱分析方法 |
CN104391178A (zh) * | 2014-12-05 | 2015-03-04 | 国家电网公司 | 一种基于Nuttall窗的时移相位差稳态谐波信号校正方法 |
CN105353215A (zh) * | 2015-12-23 | 2016-02-24 | 合肥工业大学 | 基于Nuttall窗四谱线插值FFT的谐波检测方法 |
-
2009
- 2009-11-23 CN CN2009101546820A patent/CN101718816B/zh not_active Expired - Fee Related
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102033164A (zh) * | 2010-11-16 | 2011-04-27 | 哈尔滨工业大学 | 一种计算电信号的基波分量采样信号序列的方法和系统 |
CN102033164B (zh) * | 2010-11-16 | 2012-09-26 | 哈尔滨工业大学 | 一种计算电信号的基波分量采样信号序列的方法和系统 |
CN102305886A (zh) * | 2011-05-31 | 2012-01-04 | 浙江大学 | 电网电压谐波畸变及不平衡时基波电压同步信号检测方法 |
CN102253282A (zh) * | 2011-06-10 | 2011-11-23 | 中国矿业大学 | 纳托尔窗函数连续频谱内插电力谐波参数获取方法 |
CN102253282B (zh) * | 2011-06-10 | 2014-07-16 | 中国矿业大学 | 纳托尔窗函数连续频谱内插电力谐波参数获取方法 |
CN103454495A (zh) * | 2013-09-13 | 2013-12-18 | 电子科技大学 | 自适应高精度快速频谱分析方法 |
CN103454495B (zh) * | 2013-09-13 | 2016-01-20 | 电子科技大学 | 自适应高精度快速频谱分析方法 |
CN104391178A (zh) * | 2014-12-05 | 2015-03-04 | 国家电网公司 | 一种基于Nuttall窗的时移相位差稳态谐波信号校正方法 |
CN104391178B (zh) * | 2014-12-05 | 2017-12-15 | 国家电网公司 | 一种基于Nuttall窗的时移相位差稳态谐波信号校正方法 |
CN105353215A (zh) * | 2015-12-23 | 2016-02-24 | 合肥工业大学 | 基于Nuttall窗四谱线插值FFT的谐波检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101718816B (zh) | 2011-05-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101701984B (zh) | 基于三项系数Nuttall窗插值FFT的基波与谐波检测方法 | |
CN101701982B (zh) | 基于加窗插值fft的电力系统谐波检测方法 | |
CN101603985B (zh) | 高准确度正弦信号测量方法 | |
CN106483374B (zh) | 一种基于Nuttall双窗全相位FFT的谐波间谐波检测方法 | |
CN103454497B (zh) | 基于改进加窗离散傅立叶变换的相位差测量方法 | |
CN101718816B (zh) | 基于四项系数Nuttall窗插值FFT的基波与谐波检测方法 | |
CN103245832B (zh) | 基于快速s变换的谐波时频特性参数估计方法及分析仪 | |
CN102288807B (zh) | 一种测量电网电压闪变的方法 | |
CN103575984A (zh) | 基于凯塞窗双谱线插值fft的谐波分析方法 | |
CN101261292A (zh) | 基于5项Rife-Vincent(I)窗双谱线插值FFT的基波与谐波检测方法 | |
CN103197141A (zh) | 一种测量电力系统信号频率及谐波参数的方法 | |
CN104714075B (zh) | 一种电网电压闪变包络参数提取方法 | |
CN104062528A (zh) | 基于汉宁乘积窗的信号谐波分析方法及系统 | |
CN109541312A (zh) | 一种新能源汇集地区次同步谐波检测方法 | |
CN110133564A (zh) | 一种电流传感器相位特性的同步测试方法及装置 | |
CN103399204A (zh) | 一种基于Rife-Vincent(II)窗插值FFT的谐波与间谐波检测方法 | |
CN105137181A (zh) | 基于Nuttall-Kaiser组合窗双谱线插值的谐波分析方法 | |
CN102590598B (zh) | 基于多阈值比较的超声波信号过零点预测方法 | |
CN102868402A (zh) | 一种测试模数转换器主要性能指标的测试方法 | |
CN103795411A (zh) | 基于五项最大旁瓣衰减窗三谱线插值测试sfdr的方法 | |
CN103941090A (zh) | 基于谱线能量插值的谐波测量方法 | |
CN108776263A (zh) | 基于高阶汉宁自卷积窗及改进插值算法的谐波检测方法 | |
CN109581045B (zh) | 一种满足iec标准框架的间谐波功率计量方法 | |
CN103197143A (zh) | 基于汉宁窗fft算法与遍历滤波的谐波、间谐波检测方法 | |
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 |