[发明专利]基于快速S变换的谐波时频特性参数估计方法及分析仪

说明书

技术领域

本发明涉及一种电力系统谐波时频分析方法，尤其是涉及一种基于快速S变换的谐波时频特性参数估计方法

背景技术

谐波是电能质量问题的主要根源，其时频特征参数的准确分析是研究电网谐波问题的重要内容，是谐波潮流计算、谐波电能计量、绝缘设备监测、谐波补偿与抑制的最关键依据，对确保电网安全、经济运行具有重要意义。

现有谐波频域分析方法主要是FFT及其改进算法，但不能实现谐波时域、频域参数的全面分析，无法满足暂态或突变等特性的非平稳信号谐波信号快速、准确分析的要求。且，时域分析法主要是基于各种矢量变换（d-q变换、对称分量变换等）和瞬时无功功率理论的参数估计方法，其信息量少，计算量大。

谐波时频参数包括：各次谐波幅值、相位、起止时刻、持续时间等。目前现有的谐波时域、频域联合特征分析和估计方法有：

(1) 短时傅里叶变换 Dennis Gabor于1946年建立了短时傅里叶变换（Short Time Fourier Transform，STFT），STFT算法简单，在电网谐波谱分析、电压暂降扰动检测等方面取得了应用。但存在不可克服的缺陷：①同一信号不同窗函数的STFT结果相差很远；②信号分析的时、频分辨率在窗函数选定后无法改变，只能固定分辨率，无法兼顾高频信息和低频信息；③STFT的离散形式没有正交展开，很难实现高效算法。

（2）小波变换（Wavelet Transform，WT）是具有多分辨率特性的时频局部化分析方法。优势在于：①可以对信号的任意细节聚焦，适合处理突变信号②根据频带进行信息处理，结果较为稳定可靠。但仍存在诸多局限：①频域分辨率粗糙，可能导致各频带间严重的频率混叠；②引入尺度因子后，其结果仍不是真正意义上的时频谱；③只能定性分析而不能准确检测信号的幅值与谐波分量；④小波基函数无统一选择方法，分析结果受小波基函数的影响较大；⑤算法复杂，不适合嵌入式系统实现。

(3) S变换:1996年，Stockwell等人将连续WT和STFT结合，提出了S变换时频分析方法，将一维时间信号变成二维包含局部时频谱性质的时间-频率谱。S变换利用Gauss窗的高度与频率成正比、宽度随频率的增大而变小的特点，克服了STFT窗固定的缺陷，可获得某一时刻的频率信息和某一频率点的信号幅值信息，实现动态信号的时频特性分析。近年来，S变换备受重视，但处于发展中的S变换的不足也日渐显现。其中， S变换及其改进算法的时频矩阵信息量非常大，计算繁琐，难以嵌入式实现。 

发明内容

本发明所要解决的技术问题是，克服现有技术存在的上述缺陷，提供一种信息量

大，计算简单，易嵌入式的基于快速S变换的谐波时频特性参数估计方法及分析仪。

本发明解决其技术问题所采用的技术方案是，

本发明之基于快速S变换的谐波时频特性参数估计方法，包括以下步骤：

（1）低通滤波：将被测信号输入电力系统分析仪内，根据模数转换速率和谐波分析要求，确定被测信号进行谐波分析的最高次数和低通数字滤波器的电阻电容取值，将被测信号低通滤波处理；

（2）模数转换：将步骤（1）低通滤波好的被测信号经模数转换器高速模数转换为数字信号；

（3）快速S变换：将步骤（2）转换好的被测数字信号送入数字信号处理器内，完成构建离散多高斯窗，进行快速S变换，进行时频谐波幅值、相角、起止、突变时刻参数分析；

得到谐波时频特性参数结果。

进一步，步骤（3）中，所述快速S变换的方法包括以下步骤：

（1）FFT运算：将采样好的被测数字信号x(t)通过FFT运算求得到信号频谱;

信号的连续S变换S(τ，f )定义如下