谐波测量标准推荐算法的研究.docxVIP

谐波测量标准推荐算法的研究 0 非平稳信号的谐波分析算法 随着非线性负荷的增加和能源电子能源在能源系统中的广泛应用，系统中的波形和间波形变得越来越复杂。它们的存在会对电网的环境产生较大的污染,影响通信设备,造成元器件的发热甚至损坏,还会导致波形畸变,造成继电保护误动作,甚至可能引发闪变与谐振,因此谐波与间谐波的准确检测对提高电能质量具有重要意义。 谐波分析方法一直是学者关注的热点[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17],并取得了不少理论成果。平稳信号谐波分析主要采用快速傅里叶变换(fast Fourier transform,FFT)法、加窗法、插值法、频谱校正等;对于非平稳信号的谐波分析有小波变换、希尔伯特黄变换(Hilbert-Huang transform,HHT)、神经网络、Prony法、局域均值分解(local mean decomposition,LMD)算法等。 为了更好地指导工程实践,IEEE和IEC分别制订了相关国际标准,即IEEE Std.519—1992和国际标准IEC 61000-4-7。在这些标准中,傅里叶分析作为建议的标准信号处理工具,取得电压和电流波形的谐波含量。为提高精度,FFT算法要求时域采样点数N足够大,N越大FFT的维数越高,计算量就越大。但在执行时遇到一个困难:如果每个周期采样点数取2的幂(如256),10周期内的采样点数就不是2的幂(2 560),用DSP无法编FFT程序;如果10周期采样点数取2的幂(如2 048),每个周期采样点数就不是2的幂(204.8)——电压、电流、有功功率、无功功率的计算不方便。 研究电力系统的谐波不难发现:N维FFT求出的N个频域系数大部分是不必要的,有用的仅是其中的几十个系数,因为一般情况下电力系统谐波测量只要求50次以下谐波。频域点数(H)远小于时域点数(N)的FFT已经引起了研究者的关注,文献论证了Hlog2N时直接离散傅里叶变换(discrete Fourier transform,DFT)比标准FFT计算量少。修剪算法是在标准FFT的计算树上修剪掉无用的分枝,从而减少计算量,不过对于数千点的FFT,这种修剪规则比较复杂。本文提出一种降阶FFT算法可以解决上述问题。它规则简单,计算结果与标准FFT完全相同,而计算量和内存占用量少。 1 减少fft算法 1.1 谐波次数的确定 可以证明,N维DFT可以表示成矩阵形式 为了提高精度,一般要求N足够大。随着DSP的使用N已经取到1 024、2 048,甚至4 096。相比之下,电力系统中要分析的谐波次数比N小得多,一般只是几十次。设需要分析的谐波最高次数为H-1,那么式(1)中只需要计算X的前H个分量。定义2个正整数M和K,它们满足 显然要分析0~H次谐波只需计算 1.2 降阶fft算法的计算步骤 依次取向量y的第(iK+1)个元素构成向量y0,取向量y的第(iK+2)个元素构成向量y1,取y的第(iK+M)个元素构成向量yM-1。yk∈RM×1(k=0,1,2,???,K-1)。这样式(6)就可以写成降阶形式 式中每个F0yk(k=0,1,2,???,M-1)都可以用标准M维FFT来计算,而Dk是对角矩阵, 所以降阶FFT算法的计算步骤是: 2)计算ρk=D0ηk,k=0,1,2,???,K-1,共(K-1)?(M-1)次复数乘法运算,K(M-1)次复数加法运算。 3)计算,共(K-1)次复数加法运算,1次乘法运算。 降阶FFT的总计算量是:复数乘法运算0.5KM log2M+(K-1)(M-1)次,复数加法运算KM log2M+(K-1)M次。在时域采样点数都是KM点的前提下,标准FFT的总计算量是:复数乘法0.5KM log2(KM)次,复数加法KM log2(KM)次。一般K4时,降阶FFT的计算量会比标准FFT少,且前者占用的内存单元也比后者少得多。 如果K=M=32,采样1 024个点,计算0~31次谐波,则标准FFT的总计算量是:复数乘法5 120次,复数加法10 240次;降阶FFT的总计算量是:复数乘法3 521次,复数加法6 149次。同时降阶FFT在计算过程中只需暂存32个中间结果,而标准FFT在计算过程中需要暂存1 024个中间结果,仅此一项降阶FFT就节省了近1 000个RAM。 1.3 权系数的变换 非同步采样时,为了克服频谱泄漏,通常采用加窗FFT算法。加窗FFT同样可以降阶计算,以Hanning窗为例,典型的Hanning窗是 一般取α0=1,α1=-1,w(n)称为权系数。令 则加Hanning窗后的计算公式变成 依次取向量z的第(iK+1)个元素构成向量z0,取向量z的第(iK+2)个元素构成向量z1,取z的第(iK+M)