基于复频谱共轭插值的电网信号参数估计方法及系统技术方案

本发明专利技术公开了一种基于复频谱共轭插值的电网信号参数估计方法及系统，此估计方法中的电网信号参数包括频率、幅值和相位，其中频率估计方法具体为：S1、获取N点采样序列x(n)，并对N点采样序列x(n)进行傅立叶变换得到N根DFT谱线V[0]，V[1]，...，V[N

【技术实现步骤摘要】
基于复频谱共轭插值的电网信号参数估计方法及系统


[0001]本专利技术主要涉及电网
，具体涉及一种基于复频谱共轭插值的电网信号参数估计方法及系统。

技术介绍

[0002]随着大量具有强随机的分布式电源和负荷接入电网，强随机性使电网信号受到强噪声干扰和信号参数快速波动，因此需要在很短观测时间内完成测量才有精确测量的可能性。这给现有基波参数估计方法带来严峻挑战，关系着保护控制系统响应的快速性和准确性，同时会影响电能计量的准确性。因此在强噪声和短观测窗条件下实现电网基波参数准确测量对确保电力系统安全稳定运行有重要意义。
[0003]国内外研究工作者提出了多种电网参数估计算法，从数学原理方面可分为两大类：参数方法和非参数方法。参数方法包括最大似然估计、Prony算法、MUSIC算法和矩阵束算法等，需要准确的信号频率成份个数做前置输入才能实现参数的精确估计。非参数方法主要指离散傅里叶变换(Discrete Fourier transform，DFT)类算法，因计算量小和抗噪性强被广泛应用于电网参数估计。然而，DFT类算法的测量精确度受到频率泄漏和栅栏效应的制约。为解决上述问题，加窗插值DFT(Windowed Interpolated DFT，WIpDFT)算法被提出。
[0004]经典的WIpDFT算法使用频谱中最大的两根或三根谱线的幅值进行插值计算实现频率、幅值和相位参数测量。但是，在现有算法推导过程中，负频率的影响均被忽略不计。但是，当采样信号的周波个数较少时，基波参数的测量精确度会因为这个近似操作而降低。
[0005]近年来，考虑负频谱影响的WIpDFT方法已经被提出，进一步提升了频率、幅值和相位参数的估计性能。例如：例如使用复频谱进行加窗三谱线插值，即使信号的测量时间在2.5个周期内，也能实现精确的频率估计。例如提出了基于尺度因子的IpDFT算法来减小负频谱的干扰。例如使用信号频谱的实部和虚部来进行插值计算，消除了负频谱对频率估计的影响。然而，在强噪声和短观测窗长同时存在时，上述算法的测量性能并不理想。

技术实现思路

[0006]本专利技术要解决的技术问题就在于：针对现有技术存在的技术问题，本专利技术提供一种精准快速的基于复频谱共轭插值的电网信号参数估计方法及系统。
[0007]为解决上述技术问题，本专利技术提出的技术方案为：
[0008]一种基于复频谱共轭插值的电网信号参数估计方法，电网信号参数包括频率、幅值和相位，其中频率估计方法具体为：
[0009]S1、获取N点采样序列x(n)，并对N点采样序列x(n)进行傅立叶变换得到N根DFT谱线V[0]，V[1]，...，V[N
‑
1]；
[0010]S2、获取N根DFT谱线中幅度最大的谱线，取其索引值为k
m
，然后以预定准则确定索引值k0、k1和k2；
[0011]S3、基于索引值k0、k1、k2和过程系数，得到电网频率。
[0012]优选地，在步骤S2中，幅度最大谱线的位置预先结合采样频率f
s
和采样点数N计算得到，即k
m
＝50/(f
s
/N)。
[0013]优选地，步骤S1的具体过程为：
[0014]N点采样序列对应的DFT谱线为：
[0015][0016]式中，ω
k
＝(2π/N)k，k∈[0,N
‑
1]，A0和φ0分别为电压或电流信号的幅值和相位。
[0017]优选地，步骤S2中的预定准则为：
[0018][0019]优选地，步骤S3中，基于如下基波频率估计公式得到电网频率；
[0020][0021]式中，f
s
为采样频率，为角频率，为过程参数。
[0022]优选地，步骤S3中，得到过程参数的公式如下：
[0023][0024]优选地，其中幅值和相位的估计方法中，选择模值最大的DFT谱线V[k1]进行幅值和相位估计。
[0025]优选地，得到幅值估计值的具体公式为：
[0026][0027]式中P0的估计值通过频率估计值得到。
[0028]优选地，得到相位估计值的具体公式为：
[0029][0030]式中P0的估计值通过频率估计值得到。
[0031]本专利技术还公开了一种基于复频谱共轭插值的电网信号参数估计系统，包括存储器和处理器，所述存储器上存储有计算机程序，所述计算机程序在被处理器运行时执行如上所述方法的步骤。
[0032]与现有技术相比，本专利技术的优点在于：
[0033]本专利技术基于矩形窗和频谱共轭特性设计了新的三谱线插值DFT电网信号参数估计方法，该方法推导过程考虑了基波的负频谱成份，利用频谱的复数共轭特性进行插值处理，
消除频谱泄漏和栅栏效应，实现参数的快速精确测量。
附图说明
[0034]图1为本专利技术的短采样窗口长度条件下的频率估计性能图。
[0035]图2为本专利技术的短采样窗口长度条件下的幅值估计性能图。
[0036]图3为本专利技术的短采样窗口长度条件下的相位估计性能图。
[0037]图4为本专利技术的不同信噪比条件下的频率估计性能图。
[0038]图5为本专利技术的不同信噪比条件下的幅值估计性能图。
[0039]图6为本专利技术的不同信噪比条件下的相位估计性能图。
[0040]图7为本专利技术的不同采样点条件下的频率估计性能图。
[0041]图8为本专利技术的不同采样点条件下的幅值估计性能图。
[0042]图9为本专利技术的不同采样点条件下的相位估计性能图。
[0043]图10为本专利技术的实验平台环境下的频率估计性能图。
[0044]图11为本专利技术的实验平台环境下的幅值估计性能图。
[0045]图12为本专利技术的估计方法在实施例的流程图。
具体实施方式
[0046]以下结合说明书附图和具体实施例对本专利技术作进一步描述。
[0047]如图12所示，本专利技术实施例的基于复频谱共轭插值的电网信号参数估计方法，其中电网信号参数包括频率、幅值和相位，频率估计方法具体为：
[0048]S1、获取N点采样序列x(n)，并对N点采样序列x(n)进行傅立叶变换得到N根DFT谱线V[0]，V[1]，...，V[N
‑
1]；
[0049]S2、获取N根DFT谱线中幅度最大的谱线，取其索引值为k
m
，然后以预定准则确定索引值k0、k1和k2；
[0050]S3、基于索引值k0、k1、k2和过程系数，得到电网频率。
[0051]具体地，电力系统中电压或电流信号通常建模为被高斯白噪声污染的实值正弦信号：
[0052]x(t)＝A0cos(2πf0t+φ0)+ε(t)
ꢀꢀꢀ
(1)
[0053]式(1)中，A0，f0和φ0分别为电压或电流信号的幅值、频率和相位；t为单位为秒的连本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于复频谱共轭插值的电网信号参数估计方法，其特征在于，电网信号参数包括频率、幅值和相位，其中频率估计方法具体为：S1、获取N点采样序列x(n)，并对N点采样序列x(n)进行傅立叶变换得到N根DFT谱线V[0]，V[1]，...，V[N
‑
1]；S2、获取N根DFT谱线中幅度最大的谱线，取其索引值为k
m
，然后以预定准则确定索引值k0、k1和k2；S3、基于索引值k0、k1、k2和过程系数，得到电网频率。2.根据权利要求1所述的基于复频谱共轭插值的电网信号参数估计方法，其特征在于，在步骤S2中，幅度最大谱线的位置预先结合采样频率f
s
和采样点数N计算得到，即k
m
＝50/(f
s
/N)。3.根据权利要求1或2所述的基于复频谱共轭插值的电网信号参数估计方法，其特征在于，步骤S1的具体过程为：N点采样序列对应的DFT谱线为：式中，ω
k
＝(2π/N)k，k∈[0,N
‑
1]，A0和φ0分别为电压或电流信号的幅值和相位。4.根据权利要求3所述的基于复频谱...

【专利技术属性】
技术研发人员：熊德智，邓汉钧，李庆先，杨帅，欧阳洁，谢伟峰，刘茜，
申请(专利权)人：国网湖南省电力有限公司供电服务中心计量中心国家电网有限公司，
类型：发明
国别省市：