CN106771591A - 一种复杂电力谐波的参数估计方法 - Google Patents
一种复杂电力谐波的参数估计方法 Download PDFInfo
- Publication number
- CN106771591A CN106771591A CN201710023967.5A CN201710023967A CN106771591A CN 106771591 A CN106771591 A CN 106771591A CN 201710023967 A CN201710023967 A CN 201710023967A CN 106771591 A CN106771591 A CN 106771591A
- Authority
- CN
- China
- Prior art keywords
- spectral line
- frequency
- signal
- harmonic
- frequency spectrum
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
Abstract
本发明公开了一种复杂电力谐波的参数估计方法,首先将复杂谐波进行离散化,其次对离散信号进行加窗截断,然后获取频域内多根谱线的幅值,根据谱线的理论幅值和实际幅值,将谐波参数估计问题转化为含约束条件的非线性规划问题,最终求解非线性规划问题,得到谐波参数。本发明精度高,计算量小,能较好的克服不同谐波频谱泄露干扰而带来的问题。
Description
技术领域
本发明属于电能质量分析与控制领域,特别涉及了一种复杂电力谐波的参数估计方法。
背景技术
实际电网的复杂谐波包含丰富的间谐波成分,且间谐波有一些不同于整数次谐波的特征,如果简单的使用FFT方法将会出现很多问题,这也是复杂平稳谐波分析最大的挑战。相对于整数次谐波,间谐波有如下特点:
1、间谐波频率不是基波频率的整数倍,甚至很难确定其大概周期,所以要实现对间谐波分量的同步采样基本是不可能的。
2、间谐波和整数次谐波/间谐波的频谱可能靠的很近,这就要求复杂谐波的分析方法必须具有很高的分辨率。
3、间谐波的强度很弱,使其更容易受频谱泄露的影响,尤其是当间谐波和整数次谐波靠的很近的时候,这种影响就更为明显。
在IEC规定下,谱线间的频率间隔为Δf=fs/N=5Hz,Hanning窗的主瓣宽度为20Hz,对于任意的相邻的整数次谐波而言,这已经足以将其分辨开来。再加之Hanning的加权作用,谐波分量的泄露的影响已经非常小,故高精度的估计整数次谐波分量的参数是较容易实现的。不过对于间谐波而言,考虑到间谐波上述的三个特点,FFT的局限性就非常严重,具体体现在以下三方面:
1、对于整数次谐波,选取合适的采样频率和采样点数,通过加窗减小频谱泄露的影响,再利用插值法克服栅栏效应,即可实现整数次谐波参数的高精度估计,但对于间谐波而言,由于间谐波的频率不是基波的整数倍,故即使在同步采样的时候(对整数次谐波而言),也很难实现对间谐波的高精度分析。
2、由于间谐波的幅值仅为整数次谐波的百分之几或更小,于是临近的比较强的谐波分量的频谱泄露有可能导致出现虚假间谐波分量或者间谐波谱峰消失,当间谐波和谐波间隔比较小的时候,这种现象会更突出。
3、FFT分析时,为了将信号的两个比较接近频率分量分辨开来,必须提高频率分辨率,这可以通过增加采样时间或者采用主瓣比较窄的窗函数来实现。
增加采样时间确实可以达到高的分辨率,但对于间谐波而言,单纯的增加采样时间是不合适的,这是因为实际的间谐波频谱是随时间变化的,且具有一定随机性,比如由电弧产生的间谐波分量,所以增加时长的后果可能导致前后分析的间谐波频谱发生变化,使结果失去意义。
选取具有窄的主瓣的窗函数是也可以提高频率分辨率,但是提高的程度是比较受限的,不可能达到任意的分辨率,另一方面,提高分辨率和降低频谱泄露效应也是矛盾的:如果从降低频谱泄露的角度看,应该选择旁瓣电平低衰减快的窗函数;但旁瓣衰减快的窗函数同时其主瓣宽度也会增大,这时必然就降低了频率分辨率,同时主瓣间互相干扰的程度也增大,这会对参数估计造成更严重的误差。
经过检索,在国内尚无从非线性规划的角度进行谐波分析的思路,也无将谐波分析在频域内转化成一个有约束的非线性规划问题的研究。
发明内容
为了解决上述背景技术提出的技术问题,本发明旨在提供一种复杂电力谐波的参数估计方法,克服传统FFT法分析复杂谐波所存在的弊端,将非线性规划问题引入谐波分析,开拓了复杂谐波分析的思路。
为了实现上述技术目的,本发明的技术方案为:
一种复杂电力谐波的参数估计方法,包括以下步骤:
(1)稳态谐波信号x(t)包含间谐波,对信号x(t)进行离散化,得到数字序列x(n);
(2)选择窗函数,对数字序列x(n)进行加窗截断,得到加窗信号的频谱;
(3)根据步骤(2)得到的加窗信号的频谱,计算主瓣内谱峰处某谱线的理论幅值和实际幅值yi:
yi=A W(δ-i)+ε
上式中,A为谐波信号的幅值,W(*)为窗函数的连续频谱函数,δ为由非同步采样引起的频率偏移量,ε为由于其他分量泄露引起的未知干扰量,下标i表示主瓣内谱线的索引,i=0表示主瓣内最高谱峰处谱线,最高谱峰处谱线左侧由近及远的各条谱线依次用下标i=-1,i=-2,…,i=-I表示,最高谱峰处谱线左侧由近及远的各条谱线依次用下标i=1,i=2,…,i=I表示,I为最高谱峰处谱线左右两侧分别具有的谱线条数;
(4)根据谱线的理论幅值和实际幅值yi,将谐波参数估计问题转化为含约束条件的非线性规划问题:
s.t 0≤δ<1,A>0
上式中,min表示最小值,H为目标函数,ki为各条谱线的加权因子;
(5)求解步骤(4)的非线性规划问题,得到谐波参数A和δ。
进一步地,在步骤(1)中,设稳态谐波信号则x(t)离散后的的数字序列其中,A为信号幅值,f0为信号频率,为信号初相,t表示连续时刻,fs为离散化的采样频率,n表示离散序列。
进一步地,在步骤(2)中,所述加窗信号的频谱X(k):
X(k)=AejθW(k-k0')
其中,W(*)为窗函数的连续频谱函数,N为采样点数,k是频谱所对应的数字位置,k0'=f0/Δf,Δf为频率分辨率,设频谱最高谱峰点为k0,则k0=k0'+δ,δ为由非同步采样引起的频率偏移量。
进一步地,在步骤(2)中,所述窗函数采用Hanning窗。
采用上述技术方案带来的有益效果:
本发明通过快速傅里叶变换得到频域内多根谱线的幅值,在频域内将谐波参数估计问题转化成一个有约束的非线性最优化问题,当谱间干扰较小时,该方法和比值法参数估计一致,当谐波分量之间干扰较强时,该方法能综合利用多根谱线之间的信息,精度高,计算量小,能较好的克服不同谐波频谱泄露干扰而带来的问题。本发明具有较强的通用性,可以统一目前的多谱线加权插值方法,能同时求出谐波分量的幅值的频率,且不受所加窗函数和谐波分量之间间隔的影响。
附图说明
图1为本发明的基本流程图。
具体实施方式
以下将结合附图,对本发明的技术方案进行详细说明。
如图1所示,一种复杂电力谐波的参数估计方法,具体过程如下。
设一稳态谐波信号为的形式,以采样频率fs=1/Ts对x(t)离散化,得数字序列:
其中,A为信号幅值,f0为信号频率,为信号初相,fs为离散化的采样频率,Ts为采样周期。
选择窗函数w(n),对信号x(n)加窗截断,只考虑在正频率f0附近的频谱,则加窗信号x(n)w(n)的频谱为:
X(k)=AejθW(k-k0')
其中,W(*)为窗函数的连续频谱函数,N为采样点数,k是频谱所对应的数字位置,k0'=f0/Δf,Δf为频率分辨率,设频谱最高谱峰点为k0,则k0=k0'+δ,δ为由非同步采样引起的频率偏移量。
为叙述方便,这里以主瓣内谱峰处相邻的三根谱线为例进行分析,三根谱线的对应的幅值分别为y1=|X(k1)|、y0=|X(k0)|以及y-1=|X(k-1)|,设δ为由非同步采样的引起的频率偏移量,如果所加的窗函数为Hanning窗,则无其他分量干扰时:
当无其他分量干扰时,上述三个等式严格相等,但由于间谐波的存在,三个等式是无法严格满足相等的,即三根谱线实际值应该是:
其中,ε1、εmax、ε2是由于其他分量泄露引起的未知干扰量。
从最优化的角度看,谐波参数估计问题转化成如下含约束非线性规划问题:
s.t 0≤δ<1
A>0
其中,k1、k-1、k0是各根谱线处的加权因子。
可见,这种思路的核心是通过Hanning窗频谱的形式实现三根谱线的拟合。
当无其余分量泄露的干扰的时候,上面三项平方每一项均为0,故H恒等于0,此时优化问题变成两组独立的方程组:
通过这两个方程组可以解出幅值A和频率偏移量δ。这和比值法是一致的。
当k1:k-1:k0=1:1:1或k1:k-1:k0=1:2:1,且无其余分量泄露的干扰时,该方法转化成一般的三谱线插值方法。可见,该方法可以统一目前的多谱线插值。
当存在由于其他谐波分量的泄露而导致的干扰的时,求解上述优化问题即可得到需要估计分量的的幅值和频率偏差。
实施例仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明保护范围之内。
Claims (4)
1.一种复杂电力谐波的参数估计方法,其特征在于,包括以下步骤:
(1)稳态谐波信号x(t)包含间谐波,对信号x(t)进行离散化,得到数字序列x(n);
(2)选择窗函数,对数字序列x(n)进行加窗截断,得到加窗信号的频谱;
(3)根据步骤(2)得到的加窗信号的频谱,计算主瓣内谱峰处某谱线的理论幅值和实际幅值yi:
yi=AW(δ-i)+ε
上式中,A为谐波信号的幅值,W(*)为窗函数的连续频谱函数,δ为由非同步采样引起的频率偏移量,ε为由于其他分量泄露引起的未知干扰量,下标i表示主瓣内谱线的索引,i=0表示主瓣内最高谱峰处谱线,最高谱峰处谱线左侧由近及远的各条谱线依次用下标i=-1,i=-2,…,i=-I表示,最高谱峰处谱线左侧由近及远的各条谱线依次用下标i=1,i=2,…,i=I表示,I为最高谱峰处谱线左右两侧分别具有的谱线条数;
(4)根据谱线的理论幅值和实际幅值yi,将谐波参数估计问题转化为含约束条件的非线性规划问题:
s.t 0≤δ<1,A>0
上式中,min表示最小值,H为目标函数,ki为各条谱线的加权因子;
(5)求解步骤(4)的非线性规划问题,得到谐波参数A和δ。
2.根据权利要求1所述一种复杂电力谐波的参数估计方法,其特征在于:在步骤(1)中,设稳态谐波信号则x(t)离散后的的数字序列其中,A为信号幅值,f0为信号频率,为信号初相,t表示连续时刻,fs为离散化的采样频率,n表示离散序列。
3.根据权利要求1所述一种复杂电力谐波的参数估计方法,其特征在于:在步骤(2)中,所述加窗信号的频谱X(k):
X(k)=AejθW(k-k0')
其中,W(*)为窗函数的连续频谱函数,N为采样点数,k是频谱所对应的数字位置,k0'=f0/Δf,Δf为频率分辨率,设频谱最高谱峰点为k0,则k0=k0'+δ,δ为由非同步采样引起的频率偏移量。
4.根据权利要求3所述一种复杂电力谐波的参数估计方法,其特征在于:在步骤(2)中,窗函数采用Hanning窗。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710023967.5A CN106771591B (zh) | 2017-01-13 | 2017-01-13 | 一种复杂电力谐波的参数估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710023967.5A CN106771591B (zh) | 2017-01-13 | 2017-01-13 | 一种复杂电力谐波的参数估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106771591A true CN106771591A (zh) | 2017-05-31 |
CN106771591B CN106771591B (zh) | 2019-05-10 |
Family
ID=58948209
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710023967.5A Expired - Fee Related CN106771591B (zh) | 2017-01-13 | 2017-01-13 | 一种复杂电力谐波的参数估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106771591B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108107393A (zh) * | 2017-12-06 | 2018-06-01 | 中国矿业大学 | 一种谐波分析中谱峰可信度判断方法 |
CN108107269A (zh) * | 2017-12-07 | 2018-06-01 | 中国矿业大学 | 一种谐波分析中幅值求解方法 |
CN109375023A (zh) * | 2018-12-03 | 2019-02-22 | 晶晨半导体(上海)股份有限公司 | Emi测试方法及系统 |
CN111308198A (zh) * | 2020-03-10 | 2020-06-19 | 国网江苏省电力有限公司扬州供电分公司 | 一种基于Hanning窗的加窗插值DFT的谐波测量装置及测量方法 |
CN112180163A (zh) * | 2020-09-30 | 2021-01-05 | 南方电网科学研究院有限责任公司 | 一种基于fpga的电力谐波检测系统及方法 |
CN112485522A (zh) * | 2020-12-09 | 2021-03-12 | 国网四川省电力公司电力科学研究院 | 基于电能数据感知的平顶窗函数同步相量测量方法及装置 |
CN112578345A (zh) * | 2020-11-20 | 2021-03-30 | 福瑞泰克智能系统有限公司 | 一种雷达遮挡检测方法、装置、设备及存储介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101894190A (zh) * | 2010-07-14 | 2010-11-24 | 北京交通大学 | 牵引变电站谐波发射水平估计方法 |
US20100309952A1 (en) * | 2008-10-16 | 2010-12-09 | Advantest Corporation | Distortion identification apparatus, test system, recording medium and distortion identification method |
CN104820949A (zh) * | 2015-04-23 | 2015-08-05 | 深圳供电局有限公司 | 一种配电网运行方式综合评价的方法及系统 |
CN104897960A (zh) * | 2015-06-15 | 2015-09-09 | 中南民族大学 | 基于加窗四谱线插值fft的谐波快速分析方法及系统 |
CN105353215A (zh) * | 2015-12-23 | 2016-02-24 | 合肥工业大学 | 基于Nuttall窗四谱线插值FFT的谐波检测方法 |
CN105372493A (zh) * | 2014-08-31 | 2016-03-02 | 盛吉高科(北京)科技有限公司 | 基于三条dft复数谱线的信号幅值和相位测量方法 |
-
2017
- 2017-01-13 CN CN201710023967.5A patent/CN106771591B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100309952A1 (en) * | 2008-10-16 | 2010-12-09 | Advantest Corporation | Distortion identification apparatus, test system, recording medium and distortion identification method |
CN101894190A (zh) * | 2010-07-14 | 2010-11-24 | 北京交通大学 | 牵引变电站谐波发射水平估计方法 |
CN105372493A (zh) * | 2014-08-31 | 2016-03-02 | 盛吉高科(北京)科技有限公司 | 基于三条dft复数谱线的信号幅值和相位测量方法 |
CN104820949A (zh) * | 2015-04-23 | 2015-08-05 | 深圳供电局有限公司 | 一种配电网运行方式综合评价的方法及系统 |
CN104897960A (zh) * | 2015-06-15 | 2015-09-09 | 中南民族大学 | 基于加窗四谱线插值fft的谐波快速分析方法及系统 |
CN105353215A (zh) * | 2015-12-23 | 2016-02-24 | 合肥工业大学 | 基于Nuttall窗四谱线插值FFT的谐波检测方法 |
Non-Patent Citations (2)
Title |
---|
董海波等: "畸变电源电压下的谐波检测方法与仿真研究", 《电气传动》 * |
金广厚等: "基于最有数学规划方法的谐波污染管理", 《电力自动化设备》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108107393A (zh) * | 2017-12-06 | 2018-06-01 | 中国矿业大学 | 一种谐波分析中谱峰可信度判断方法 |
CN108107393B (zh) * | 2017-12-06 | 2019-11-05 | 中国矿业大学 | 一种谐波分析中谱峰可信度判断方法 |
CN108107269A (zh) * | 2017-12-07 | 2018-06-01 | 中国矿业大学 | 一种谐波分析中幅值求解方法 |
CN108107269B (zh) * | 2017-12-07 | 2019-11-05 | 中国矿业大学 | 一种谐波分析中幅值求解方法 |
CN109375023A (zh) * | 2018-12-03 | 2019-02-22 | 晶晨半导体(上海)股份有限公司 | Emi测试方法及系统 |
CN111308198A (zh) * | 2020-03-10 | 2020-06-19 | 国网江苏省电力有限公司扬州供电分公司 | 一种基于Hanning窗的加窗插值DFT的谐波测量装置及测量方法 |
CN111308198B (zh) * | 2020-03-10 | 2021-09-24 | 国网江苏省电力有限公司扬州供电分公司 | 一种基于Hanning窗的加窗插值DFT的谐波测量方法 |
CN112180163A (zh) * | 2020-09-30 | 2021-01-05 | 南方电网科学研究院有限责任公司 | 一种基于fpga的电力谐波检测系统及方法 |
CN112578345A (zh) * | 2020-11-20 | 2021-03-30 | 福瑞泰克智能系统有限公司 | 一种雷达遮挡检测方法、装置、设备及存储介质 |
CN112485522A (zh) * | 2020-12-09 | 2021-03-12 | 国网四川省电力公司电力科学研究院 | 基于电能数据感知的平顶窗函数同步相量测量方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN106771591B (zh) | 2019-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106771591A (zh) | 一种复杂电力谐波的参数估计方法 | |
CN101403774B (zh) | 一种基于非同步采样的谐波分析方法 | |
CN101441233A (zh) | 基于Kaiser窗双谱线插值FFT的基波与谐波检测方法 | |
CN108918964B (zh) | 一种稀疏性增强的谐波分析方法 | |
CN107085140B (zh) | 基于改进的SmartDFT算法的非平衡系统频率估计方法 | |
CN103575984A (zh) | 基于凯塞窗双谱线插值fft的谐波分析方法 | |
CN105137180B (zh) | 基于六项余弦窗四谱线插值的高精度谐波分析方法 | |
CN103308804A (zh) | 基于快速k-s变换电能质量扰动信号时频参数提取方法 | |
CN105353215A (zh) | 基于Nuttall窗四谱线插值FFT的谐波检测方法 | |
CN106018956B (zh) | 一种加窗谱线插值的电力系统频率计算方法 | |
CN110095650A (zh) | 基于五项Rife-Vincent(I)窗的四谱线插值FFT的复杂谐波检测分析方法 | |
CN109946512A (zh) | 一种改进频域插值的动态功率分析方法 | |
CN109855852A (zh) | 一种基于短时校正傅里叶变换的信号包络提取方法 | |
CN103941090A (zh) | 基于谱线能量插值的谐波测量方法 | |
Carpinelli et al. | A new advanced method for an accurate assessment of harmonic and supraharmonic distortion in power system waveforms | |
CN103245830B (zh) | 一种结合ar谱估计与非线性优化的间谐波检测方法 | |
CN103543331B (zh) | 一种计算电信号谐波和间谐波的方法 | |
CN108390698A (zh) | 一种基于插值fft算法的电力线载波参数测量方法 | |
CN108107269B (zh) | 一种谐波分析中幅值求解方法 | |
Bracale et al. | On some spectrum estimation methods for analysis of nonstationary signals in power systems. Part I. Theoretical aspects | |
CN107315714B (zh) | 一种去卷积功率谱估计方法 | |
Zhang et al. | Frequency shifting and filtering algorithm for power system harmonic estimation | |
CN109709378A (zh) | 瞬变电信号的频率和幅值自适应算法 | |
CN113295923A (zh) | 一种基于改进s变换的VFTO信号频谱分析方法 | |
CN110320400B (zh) | 准同步采样和改进能量算子的电压闪变包络参数提取方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB02 | Change of applicant information |
Address after: 221116 No. 1 University Road, copper mountain, Jiangsu, Xuzhou Applicant after: China University of Mining & Technology Address before: 221116 Xuzhou University Road, Jiangsu, No. 1 Applicant before: China University of Mining & Technology |
|
CB02 | Change of applicant information | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190510 Termination date: 20210113 |
|
CF01 | Termination of patent right due to non-payment of annual fee |