CN105911341B - 一种谐波无功功率的测量方法 - Google Patents
一种谐波无功功率的测量方法 Download PDFInfo
- Publication number
- CN105911341B CN105911341B CN201610247907.7A CN201610247907A CN105911341B CN 105911341 B CN105911341 B CN 105911341B CN 201610247907 A CN201610247907 A CN 201610247907A CN 105911341 B CN105911341 B CN 105911341B
- Authority
- CN
- China
- Prior art keywords
- harmonic wave
- discrete
- reactive power
- hanning
- convolution
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R21/00—Arrangements for measuring electric power or power factor
- G01R21/001—Measuring real or reactive component; Measuring apparent energy
- G01R21/003—Measuring reactive component
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
本发明公开了一种谐波无功功率的测量方法,包括(1)构建长度为N的p阶离散Hanning自卷积窗序列;(2)利用得到的p阶离散Hanning自卷积窗序列处理含有谐波分量的离散信号x(n),得到加窗后的离散频谱;(3)利用三谱线插值确定谱线的准确位置:以及(4)通过上述得到的谱线准确位置,即可得到谐波的频率、幅值和相角,进而可直接计算无功功率,完成无功功率的测量。本发明的方法可以快速精确地测量无功功率,而且测量效率高,计算量和步骤简化,精度高。
Description
技术领域
本发明涉及电力自动化领域的电能测量技术,具体涉及一种谐波无功功率的测量方法。
背景技术
无功功率控制是电网节能的重要措施,在提高电能质量、降低网损、提高电网运行的稳定性和安全性等方面具有极为重要的意义。合理的无功功率控制及补偿的前提是进行准确的无功测量,同时无功电量也是衡量供电系统电能质量的重要指标之一,因此无功功率测量具有重要意义。
随着非线性负载的广泛使用,电网中的谐波越来越丰富。关于谐波无功功率测量的方法有很多,例如可以使用Hilbert数字滤波器来进行测量,或通过基于Hibert变换的FFT插值重构的方法进行无功功率测量。但是,对于使用Hilbert数字滤波器测量在非同步采样时,会存在较大误差,测量精度不高,而采用基于Hibert变换的FFT插值重构的方法时,其步骤多,计算量比较大,效率以及精度都不高。
发明内容
针对现有技术的以上缺陷或改进需求,本发明提供了一种谐波无功功率的测量方法,其基于Hanning自卷积窗和三谱线插值方式,可以快速精确地测量无功功率,而且测量效率高,计算量和步骤简化,精度高。
为实现上述目的,按照本发明,提供一种谐波无功功率的测量方法,包括如下具体步骤:
(1)构建长度为N的p阶离散Hanning自卷积窗序列,其中所述自卷积窗序列通过长度为M的Hanning窗做离散自卷积得到,具体包括将p个长度为M的Hanning窗序列做离散自卷积,得到长度为pM-p+1的序列,然后在尾部补零,从而得到长度为 N=pM的Hanning自卷积窗序列;
(2)利用上述得到的p阶离散Hanning自卷积窗序列处理含有谐波分量的离散信号x(n),得到加窗后的离散频谱,具体为
式中,H为谐波的个数,f1、fh、Ah、分别表示第1次谐波的频率、第h次谐波的频率、幅值和相角,Δf=fs/N,WHp为离散频谱,kh为峰值频点,fs为采样频率。
(3)利用三谱线插值确定谱线的准确位置,具体包括如下过程:
首先,通过中间参数其中,
kha为kh附近抽样得到的幅值最大的谱线,谱线kha-1和kha+1分别为与kha相邻的两根谱线,
进一步得到:
计算δ=v-1(α),对上式进行拟合,进而求解得到δ。
因为α是可以通过求解得到的,通过拟合可以得到δ和α的关系(即δ=v-1(α)),就可以得到δ。
(4)通过上述得到的谱线准确位置,即可得到谐波的频率、幅值和相角,进而可直接计算无功功率,完成无功功率的测量。
作为本发明的进一步优选,所述谐波的信号频率为fh=(kha+δ)Δf。
作为本发明的进一步优选,所述谐波的幅值通过对三谱线进行加权求得,即幅值为:
Ah=N-1(y1+2y2+y3)g(δ),
其中,N-1为长度N的倒数,
作为本发明的进一步优选,所述谐波的相位为
作为本发明的进一步优选,所述含有谐波分量的离散信号x(n)为通过含有谐波分量的待处理信号x(t)以采样频率fs均匀采样得到的离散时间信号。
作为本发明的进一步优选,所述p阶Hanning自卷积窗通过p-1个Hanning窗进行自卷积运算得到。
作为本发明的进一步优选,由p-1个长度为M的Hanning窗自卷积得到长度为N 的p阶Hanning自卷积窗的离散频谱WHp为:
其中 k为离散点0,…,N-1。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,具有以下有益效果:
(1)本发明中,无功功率计算所需要的谐波的频率、幅值和相位通过采用基于Hanning自卷积窗,利用Hanning自卷积窗的优越性能可以抑制信号的频谱泄漏,使得计算过程简洁优化;
(2)本发明中,谐波的频率、幅值和相位的计算同时利用幅值最大的谱线及其相邻2根谱线进行三谱线插值,以确定谱线的准确位置,进而得到谐波的频率、幅值和相角,计算过程简单、精度高;
(3)本发明中利用所求的谐波的频率、幅值和相位参数直接求解无功功率,简化了算法的步骤,易于在嵌入式系统中实现。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
按照本发明一个实施例的基于Hanning自卷积窗和三谱线插值的谐波无功功率测量算法,其步骤如下:
(1)构建长度为N的p阶离散Hanning自卷积窗序列,自卷积窗序列通过长度为 M的Hanning窗做离散自卷积得到。
具体地,同将p个长度为M的Hanning窗序列做离散自卷积,得到长度为pM-p+1 的序列,在尾部补零得到长度为N=pM的序列。
其中,长度为M的Hanning窗的离散表达式为
定义p阶Hanning自卷积窗为p-1个Hanning窗进行自卷积运算,即
由p-1个长度为M的Hanning窗自卷积得到长度为N的p阶Hanning自卷积窗,其离散频谱为
其中k为离散点0,…,N-1。
(2)设含有谐波分量的信号x(t)以采样频率fs均匀采样得到离散时间信号为
式中H表示谐波的个数。h=1时,f1、A1、分别为基波的频率、幅值和初相位;当h≠1时,fh、Ah、分别表示为第h次谐波的频率、幅值和相角。
用p阶Hanning自卷积窗处理信号x(n),得到加窗后的离散频谱为
式中Δf=fs/N。
(3)由于栅栏效应,kh一般不是整数。设在峰值频点kh附近抽样得到的幅值最大的谱线为kha,其左边谱线为kha-1,右边谱线为kha+1。记δ=kh-kha。
在非同步采样的情况下,最大谱线以及其左右谱线都具有较大值,为了充分利用这些谱线信息,采用三谱线插值。令
y1=|X(kha-1)|
y2=|X(kha)|
y3=|X(kha+1)|
设中间参数得到
记α=v(δ),其反函数记为δ=v-1(α)。计算δ=v-1(α),对上式进行拟合,进而求解得到δ。
因为α是可以通过求解得到的,通过拟合可以得到δ和α的关系(即δ=v-1(α)),就可以得到δ。
本实施例中优选采用最小二乘法多项式拟合,可以优选通过在Matlab中用函数polyfit实现。
(4)根据步骤(3)的结果,对谐波的频率、幅值以及相位角分别进行计算,具体为:
信号频率为
fh=(kha+δ)Δf
幅值通过对三谱线进行加权求得
记
得到幅值计算公式
Ah=N-1(y1+2y2+y3)g(δ)
相位修正公式
最后,根据上述谐波的频率、幅值以及相位角三个参数,直接即可计算得到无功功率。例如可以根据Budeanu计算无功功率。
本方法中,采用Hanning自卷积窗,具有良好的抑制频谱泄漏的性能,并且能够根据需要选取自卷积的阶数p。另外,采取了三谱线的插值算法,即利用了幅值最大的谱线以及其相邻的两根谱线,充分利用了频谱的信息,能够在频率波动的情况下准确计算谐波参数。另外,根据Budeanu关于无功功率的定义直接计算无功功率,简化了算法流程,增加了执行效率。
而且,通过窗函数自卷积过后的补零运算,使最终FFT计算长度为2的整数次幂,便于使用基于2的FFT算法,便于在DSP等嵌入式设备中实现。
在一个具体实施例中,设置采样频率为3200Hz,采样的点数N=1024。(2)根据阶数不同,选取不同长度的离散Hanning窗。当p=1时,M=1024;p=2时,M=512; p=4时,M=256。(3)离散Hanning窗做p阶自卷积,并且通过补零后得到长度为N 的自卷积窗。(4)加窗,即利用自卷积窗截断采样得到的谐波信号。(5)对加窗后的谐波信号进行FFT变换,得到信号的频谱。(6)对频谱进行三谱线插值,得到谐波的频率、幅值和相角。(6)根据Budeanu关于谐波的计算公式计算谐波。
p=1,M=1024时
δ=0.01281132α5-0.07313921α3+0.66660636α
g(δ)=0.12388464δ4+0.52603680δ2+1.33333957
p=2,M=512时
δ=0.01371699α5-0.08706146α3+1.04093202α
g(δ)=0.00029875δ4+0.00221128δ2+0.00908164
p=4,M=256时
δ=0.01857964α5-0.12973263α3+1.80839104α
g(δ)=(0.00210426δ4+0.02860162δ2+0.20613957)*10-5
通过实验可知,当频率波动时,得到的结果精度很高,完全满足测量的要求。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种谐波无功功率的测量方法,包括如下具体步骤:
(1)构建长度为N的p阶离散Hanning自卷积窗序列,其中所述自卷积窗序列通过长度为M的Hanning窗做离散自卷积得到,具体包括:将p个长度为M的Hanning窗序列做离散自卷积,得到长度为pM-p+1的序列,然后在尾部补零,从而得到长度为N=pM的Hanning自卷积窗序列;
(2)利用上述得到的p阶离散Hanning自卷积窗序列处理含有谐波分量的离散信号x(n),得到加窗后的离散频谱,具体为
式中,H为谐波的个数,f1、fh、Ah、分别表示第1次谐波的频率、第h次谐波的频率、幅值和相角,Δf=fs/N,WHp为离散频谱,kh为峰值频点,fs为采样频率;
(3)利用三谱线插值确定谱线的准确位置,具体包括如下过程:
首先,确定中间参数其中,
y1=|X(kha-1)|
y2=|X(kha)|
y3=|X(kha+1)|,
kha为kh附近抽样得到的幅值最大的谱线,谱线kha-1和kha+1分别为与kha相邻的两根谱线,
进一步得到:
对上式进行拟合,得到δ和α的关系,即δ=v-1(α),进而计算得到δ;
(4)通过上述得到的谱线准确位置,即可得到谐波的频率、幅值和相角,进而可直接计算无功功率,完成无功功率的测量。
2.根据权利要求1所述的一种谐波无功功率的测量方法,其中,所述谐波的信号频率为fh=(kha+δ)Δf。
3.根据权利要求1所述的一种谐波无功功率的测量方法,其中,所述谐波的幅值通过对三谱线进行加权求得,即幅值为:
Ah=N-1(y1+2y2+y3)g(δ),
其中,N-1为长度N的倒数,
4.根据权利要求1所述的一种谐波无功功率的测量方法,其中,所述谐波的相位为
5.根据权利要求1-4中任一项所述的一种谐波无功功率的测量方法,其中,所述含有谐波分量的离散信号x(n)为通过含有谐波分量的待处理信号x(t)以采样频率fs均匀采样得到的离散时间信号。
6.根据权利要求1-4中任一项所述的一种谐波无功功率的测量方法,其中,所述p阶Hanning自卷积窗通过p-1个Hanning窗进行自卷积运算得到。
7.根据权利要求1-4中任一项所述的一种谐波无功功率的测量方法,其中,由p-1个长度为M的Hanning窗自卷积得到长度为N的p阶Hanning自卷积窗的离散频谱WHp为:
其中k为离散点0,…,N-1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610247907.7A CN105911341B (zh) | 2016-04-20 | 2016-04-20 | 一种谐波无功功率的测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610247907.7A CN105911341B (zh) | 2016-04-20 | 2016-04-20 | 一种谐波无功功率的测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105911341A CN105911341A (zh) | 2016-08-31 |
CN105911341B true CN105911341B (zh) | 2018-08-31 |
Family
ID=56747778
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610247907.7A Expired - Fee Related CN105911341B (zh) | 2016-04-20 | 2016-04-20 | 一种谐波无功功率的测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105911341B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107525987B (zh) * | 2017-08-01 | 2020-06-19 | 合肥工业大学 | 基于自适应序列全相位dft的同期并网参数测量方法 |
CN108107393B (zh) * | 2017-12-06 | 2019-11-05 | 中国矿业大学 | 一种谐波分析中谱峰可信度判断方法 |
CN108107269B (zh) * | 2017-12-07 | 2019-11-05 | 中国矿业大学 | 一种谐波分析中幅值求解方法 |
CN109490627A (zh) * | 2018-12-25 | 2019-03-19 | 新疆大学 | 基于混合卷积窗的新型电力系统谐波分析方法 |
CN109709378B (zh) * | 2019-01-21 | 2021-02-19 | 国网安徽省电力有限公司六安供电公司 | 瞬变电信号的频率和幅值自适应算法 |
CN114966373A (zh) * | 2022-05-24 | 2022-08-30 | 江苏科技大学 | 一种集成电路模数转换芯片参数测试方法与测试系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5262715A (en) * | 1991-07-22 | 1993-11-16 | Landis & Gyr Betriebs Ag | Arrangement for the measurement of reactive power or reactive energy |
CN103308766A (zh) * | 2013-05-15 | 2013-09-18 | 湖南大学 | 一种基于凯撒自卷积窗双谱线插值fft谐波分析方法及其装置 |
CN103543325A (zh) * | 2012-07-11 | 2014-01-29 | 索尼公司 | 功耗管理设备、功耗管理系统、功耗管理方法以及程序 |
CN103339518B (zh) * | 2011-02-04 | 2015-12-16 | 富士通株式会社 | 电源排插和功率测量方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5957397B2 (ja) * | 2013-03-05 | 2016-07-27 | 株式会社東芝 | 状態推定装置、状態推定方法及び状態推定プログラム |
-
2016
- 2016-04-20 CN CN201610247907.7A patent/CN105911341B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5262715A (en) * | 1991-07-22 | 1993-11-16 | Landis & Gyr Betriebs Ag | Arrangement for the measurement of reactive power or reactive energy |
CN103339518B (zh) * | 2011-02-04 | 2015-12-16 | 富士通株式会社 | 电源排插和功率测量方法 |
CN103543325A (zh) * | 2012-07-11 | 2014-01-29 | 索尼公司 | 功耗管理设备、功耗管理系统、功耗管理方法以及程序 |
CN103308766A (zh) * | 2013-05-15 | 2013-09-18 | 湖南大学 | 一种基于凯撒自卷积窗双谱线插值fft谐波分析方法及其装置 |
Non-Patent Citations (2)
Title |
---|
一种用于分析电网谐波的多谱线插值算法;吴静 等;《中国电机工程学报》;20060430;第26卷(第8期);第55-60页 * |
基于卷积窗的电力系统谐波误差估计与数值模拟;张介秋 等;《中国电机工程学报》;20041231;第24卷(第12期);第34-37页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105911341A (zh) | 2016-08-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105911341B (zh) | 一种谐波无功功率的测量方法 | |
CN104502700B (zh) | 电力信号的正弦参数测量方法和系统 | |
CN102435844B (zh) | 一种频率无关的正弦信号相量计算方法 | |
CN104502706B (zh) | 电力信号的谐波幅值测量方法和系统 | |
CN104535836B (zh) | 电力信号的基波频率测量方法和系统 | |
CN104897961B (zh) | 基于互乘法窗函数的三谱线插值fft谐波分析方法及系统 | |
CN104897960B (zh) | 基于加窗四谱线插值fft的谐波快速分析方法及系统 | |
CN104502698B (zh) | 电力信号的频率测量方法和系统 | |
CN107478896A (zh) | 一种基于级联广义积分器的频率自适应谐波电流检测方法 | |
CN103353550A (zh) | 一种测量电力系统信号频率及谐波参数的方法 | |
CN104459320B (zh) | 电力信号的谐波相位测量方法和系统 | |
CN109831184B (zh) | 一种带通滤波方法、系统、存储介质及终端设备 | |
CN105811408B (zh) | 面向mmc换流站的谐波能量传递函数建模方法 | |
CN104502675B (zh) | 电力信号的基波幅值测量方法和系统 | |
CN104459321B (zh) | 电力信号的基波相位测量方法和系统 | |
CN106980044B (zh) | 一种适应风电接入的电力系统谐波电流估计方法 | |
CN102636693A (zh) | 一种结合fft与非线性最小二乘的谐波分析算法 | |
CN103236691A (zh) | 基于复仿射数学理论的三相不平衡潮流计算方法 | |
CN103018555A (zh) | 一种高精度的电力参数软件同步采样方法 | |
CN105486921A (zh) | 凯撒三阶互卷积窗三谱线插值的谐波与间谐波检测方法 | |
CN105334381B (zh) | 一种交流有功功率的测量方法和装置 | |
CN102495285B (zh) | 对称窗函数功率重心估计电力谐波参数的方法 | |
CN106018960A (zh) | 一种基于压缩传感的同步相量测量方法 | |
CN110068729B (zh) | 一种信号相量计算方法 | |
CN104111373A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180831 Termination date: 20200420 |
|
CF01 | Termination of patent right due to non-payment of annual fee |