[发明专利]一种基于改进不完全S变换的基波和谐波参数估计方法

说明书

本发明公开了一种基于改进不完全S变换的基波和谐波参数估计方法，该方法包括数据采集及预处理模块(1)、改进不完全S变换计算模块(2)、基波和谐波参数估计模块(3)等3个功能模块；数据采集和预处理模块(1)对采集的数据进行规格化处理；改进不完全S变换计算模块(2)对规格化数据进行改进不完全S变换运算；基波和谐波参数估计模块(3)首先利用改进不完全S变换的行幅值向量中心段均值求得基波和谐波的幅值，然后根据幅值采用辛格函数的反函数值以及频谱对应基波和谐波次大谱线和谱峰频率坐标差值的极性，得到基波和谐波的频率和相位校正值，进而求得基波和谐波的频率及相位；本发明为电力系统基波和谐波参数的检测提供了一种实时精确的手段。

技术领域

本发明涉及电力系统电能质量信号分析和参数估计领域，特别涉及一种基波和谐波参数估计方法。

背景技术

电力系统谐波是一种危害很大的电能质量问题，为有针对性治理谐波，需要快速准确估计其参数。

已有的谐波参数估计方法常采用FFT(快速傅里叶变换)进行分析，但FFT存在栅栏效应和频谱泄漏，只有在严格同步采样的基础上才能得到精确的参数值。为克服在非同步采样情况下FFT方法参数估计的不足，在FFT方法基础上产生了一些频谱分析的校正方法，其原理都是根据FFT分析的加窗类型，以及FFT的相关性质，确定信号的频谱分布函数，结合实际信号的频谱分布，反推出信号的精确参数。但相关方法对于谱峰检测并没有给出明确方法，并且仅适用于稳态信号分析。

S变换在电能质量信号检测和参数估计领域已有众多应用。不完全S变换的求解仅针对FFT谱峰进行运算，大大减少了运算时间，且在运算中考虑了谱峰两侧主要谱线的作用，在幅值计算中消除了FFT频谱泄漏造成的误差，可以得到信号的精确幅值。同时，在不完全S变换的计算过程中含有对谱峰精确检测方法，形成了完整应用技术。

发明内容

本发明提供一种基于改进不完全S变换的基波和谐波参数估计方法，通过仅针对基波和谐波谱峰频率坐标进行S变换的后续计算，以及利用不完全S变换在幅值估计对频谱泄漏的修正特性，以解决基波和谐波参数估计的准确性和实时性问题。

一种基于改进不完全S变换的基波和谐波参数估计方法，包括数据采集及预处理模块(1)、改进不完全S变换计算模块(2)、基波和谐波参数估计模块(3)，所述3个功能模块依次连接实现电力系统基波和谐波参数的检测，其数据采集及预处理模块(1)根据基波和谐波参数检测指令提取缓冲区的数据段，并对提取的数据进行规格化处理，改进不完全S变换计算模块(2)对预处理后的数据做改进不完全S变换，得到对应基波和谐波的不完全S变换的复数矩阵，基波和谐波参数估计模块(3)首先根据改进不完全S变换计算模块(2)的计算结果求其模矩阵，接着通过计算模矩阵行幅值向量中心段均值求得基波和谐波的幅值，然后采用误差控制的逐步逼近法求得辛格函数对应基波和谐波幅值的反函数值绝对值，并结合频谱序列中对应基波和谐波的次大谱线频率坐标与谱峰频率坐标差值的极性，确定基波和谐波的频率校正值，进而得到基波和谐波的频率值，最后依据频率校正值求得相位校正值，进而得到基波和谐波的相位值。

所述的一种基于改进不完全S变换的基波和谐波参数估计方法，其数据采集及预处理模块(1)所进行的规格化处理为将采集的数据换算为实际值再除以额定值。

所述的一种基于改进不完全S变换的基波和谐波参数估计方法，其改进不完全S变换计算模块(2)通过功率谱包络的动态测度检测基波和谐波频率坐标值，动态测度标准化阈值根据信号噪声强度设定为0.15~0.45之间，所谓标准化为动态测度值除以数据长度。

所述的一种基于改进不完全S变换的基波和谐波参数估计方法，其改进不完全S变换计算模块(2)针对检测到的基波和谐波频率分别采用不同的窗宽系数进行不完全S变换的高斯窗函数的加窗运算。

所述的一种基于改进不完全S变换的基波和谐波参数估计方法，其基波和谐波参数估计模块(3)的计算步骤如下：