CN101113995A - 基于Nuttall窗双峰插值FFT的基波与谐波检测方法 - Google Patents

Abstract

本发明公开了一种基于Nuttall窗双峰插值FFT(快速傅立叶变换)的基波与谐波检测方法：对被测信号进行模数转换和低通滤波处理，采用旁瓣电平低、衰减速度快的Nuttall窗对信号进行加窗FFT运算，运用双峰插值算法检测被测信号中的基波和各谐波量。包括下列方面：(1)模数转换与数字低通滤波；(2)基于Nuttall窗的FFT运算；(3)双峰插值和多项式拟合求取基波与谐波参数。

Description

基于Nuttall窗双峰插值FFT的基波与谐波检测方法

所属技术领域

本发明涉及一种用于信号中基波与谐波参数检测的方法。具体是一种基于Nuttall窗双峰插值FFT(快速傅立叶变换)的基波与谐波参数检测方法，属于信号处理技术领域，也可用于其它信号的谐波分析。

背景技术

以电力信号检测为例，随着电力电子技术的发展，非线性设备在电力系统中得到广泛应用，非线性设备带来的谐波问题对电力系统安全、稳定、经济运行构成的威胁日趋严重。基波与谐波的精确测量可为电网电能计量、谐波潮流计算、设备入网检测、电力系统谐波补偿与抑制等提供科学依据。

基波与谐波分量检测方法一般有：(1)基于频域分析的FFT方法，其特点是电网频率波动时，检测精度较低；(2)滤波器检测谐波方法，其特点是难以获得理想频率特性；(3)基于瞬时无功功率理论的方法，其特点是计算量大，处理复杂；(4)基于神经网络理论和小波变换的方法，其特点是计算量大，难以在嵌入式系统中实现。

基于FFT的基波、谐波分析算法，易于在DSP(数字信号处理器)上实现，是当今应用最广泛的一种谐波分析方法。在电力系统中，电网电压、电流畸变导致基波频率变化，由于非同步采样造成FFT算法存在频谱泄露和栅栏效应等问题，使得基波与谐波检测的精度受到影响。减少FFT算法的频谱泄露和栅栏效应影响、提高电测量中基波和各谐波检测精度是电测量信号分析和电能质量管理中的难题。

已有专利文件“测量工厂谐波的方法和测量仪”(200310105446.2)，“计量工厂谐波的方法和计量仪”(200310105446.2)，“电力系统谐波定量计量方法和计量仪”(98110414.2)等。其发明的目的是精确定量地测量电力系统中谐波电压、电流的动态特性。已有技术的不足之处是：已有技术提供了一些基波与谐波测量的电气设备设计方法，但由于谐波具有多样性、随机性和多态性等特征，基波频率波动造成的频谱泄露和栅栏效应依然存在，因此实时、高精度的基波与谐波测量分析难以实现。

发明内容

为克服已有技术的不足，本发明的目的在于提供基于Nuttall窗的双峰插值FFT方法，该方法能快速、准确检测信号中基波和各谐波分量。

设包含多项整数次谐波的时域被测信号为

其中，M为所含谐波的项数；r_m代表谐波的次数；当m＝0时，r₀＝1，f₀、A₀和₀分别为基波的频率、幅值和初相角；当m≠0时，r_m、A_m和_m分别为第m项谐波的次数、幅值和初相角。信号x(t)经过采样率为f_s的数据采集系统后，得到离散序列

根据傅立叶变换的定义，信号x(n)的连续傅立叶变换为

$X\left(f\right)=\underset{n=-\∞}{\overset{\∞}{\mathrm{\Σ}}}x\left(n\right)e^{-j2\mathrm{\πfn}}$

在嵌入式系统上或计算机上实现时，数据总是有限长的，即相当于信号x(n)被一个窗函数w(n)(n＝0，1，…N-1)截短为N点长序列。加窗后，信号频谱由原来的谱线变为以f₀为中心向两边扩展的连续谱，信号频谱旁瓣之间的相互干扰，能量泄露到整个频带，即为频谱泄露。

对X(f)以离散频率间隔Δf＝f_s/N进行离散化，可得到X(kΔf)(k＝0，1，…N-1，代表离散频谱线的序号)。由于电网频率并非恒定，存在一定波动，加上被测信号中可能存在间谐波(即非整数次谐波)，因而难以严格做到同步采样。非同步采样造成峰值频率很难正好位于离散谱线频点上，即f₀＝k₀Δf(k₀为非整数)，导致的对信号峰值点的观测偏差，即为栅栏效应。

基于以上分析，本发明提出的基于Nuttall窗的基波与谐波检测方法，其特征在于：对被测信号进行模数转换和低通滤波处理，采用旁瓣电平低、衰减速度快的Nuttall窗对信号进行加窗FFT运算，运用双峰插值算法检测被测信号中的基波和各谐波量。它包括下列步骤：

(1)模数转换与低通滤波

高速模数转换器将被测信号(模拟量)x(t)转换为数字量，采样频率为f_s，该数字量作为基波与谐波检测处理的输入量。根据谐波分析精度要求，确定需检测的最高次谐波次数及频率f_MAX，设计截止频率为稍大于f_MAX的低通数字滤波器，输入量经过数字低通滤波器后，得到去除高频信号分量的输入量x(n)。

式中，M为所含最高谐波的次数。

(2)构建离散四项三阶Nuttall窗

Nuttall窗是一种余弦组合窗，其离散表达式为

$w\left(n\right)=\underset{g=0}{\overset{G-1}{\mathrm{\Σ}}}{(-1)}^g{b}_g\cos (2\mathrm{\πn}\·g/N)$

式中，G为Nuttall窗函数的项数；N为FFT运算的数据长度；也是窗函数的长度(n＝1，2，…N-1)；b_g应满足以下约束条件

$\underset{g=0}{\overset{G-1}{\mathrm{\Σ}}}{b}_g=1$

$\underset{g=0}{\overset{G-1}{\mathrm{\Σ}}}{(-1)}^g{b}_g=1$

四项三阶Nuttall窗的系数分别是b₀＝0.338946、b₁＝0.481973、b₂＝0.161054、b₃＝0.018027，其离散时域特性如图2所示。图3给出了四项三阶Nuttall窗的幅频特性，四项三阶Nuttall窗具有优良的旁瓣特性，旁瓣衰减速率达到30分贝/倍频程，旁瓣电平达到了-83分贝，采用四项三阶Nuttall窗对输入信号进行加窗处理。

(3)基于Nuttall窗的FFT运算

根据检测精度的要求，确定FFT运算的数据长度N。为便于嵌入式系统实现，本发明考虑只含有单次谐波存在的情况，将去除高频分量的输入信号量x(n)进行加窗FFT运算，可得到加窗后信号的离散傅立叶变换的表达式为

其中，W(·)是w(n)的频谱函数。

(4)双峰插值算法

非同步采样情况下，离散峰值频率f₀＝k₀Δf不处于离散谱线频点上，设需要计算的各频率峰值点k₀附近幅值最大和次最大的谱线分别k₁和k₂。显然，k₁≤k₀≤k₂(k₂＝k₁+1)，这两条谱线的幅值分别为y₁＝|X(k₁Δf)|和y₂＝|X(k₂Δf)|，设 $\mathrm{\β}=\frac{y_2-y_1}{y_2+y_1}$ 、α＝k₀-k₁-0.5，可知α取值范围为[-0.5，0.5]，可得以α为自变量，以β为因变量的函数关系β＝h(α)

$\mathrm{\β}=\frac{\left|W\left(\frac{2\mathrm{\π}(-\mathrm{\α}+0.5)}{N}\right)\right|-\left|W\left(\frac{2\mathrm{\π}(-\mathrm{\α}-0.5)}{N}\right)\right|}{\left|W\left(\frac{2\mathrm{\π}(-\mathrm{\α}+0.5)}{N}\right)\right|+\left|W\left(\frac{2\mathrm{\π}(-\mathrm{\α}-0.5)}{N}\right)\right|}$

其反函数为α＝h^-1(β)。由β可求出参数α，则被测信号的频率计算公式为

f₀＝k₀Δf＝(α+k₁+0.5)Δf

被测信号的幅值计算公式为

$A=\frac{A_1\left|W\left(\frac{2\mathrm{\π}(k_1-k_0)}{N}\right)\right|+A_2\left|W\left(\frac{2\mathrm{\π}(k_2-k_0)}{N}\right)\right|}{\left|W\left(\frac{2\mathrm{\π}(k_1-k_0)}{N}\right)\right|+\left|W\left(\frac{2\mathrm{\π}(k_2-k_0)}{N}\right)\right|}$

被测信号的初相位计算公式

＝arg[X(kΔf)]+π/2-arg[W2π(kΔf-k₀)T_s)]

由此可计算出被测信号基波与谐波的真实频率和各项参数。

本发明所述的模数转换是指通过对被测信号(可以是电压或电流信号等)进行高速模数转换，将被测信号转换为数字量。

本发明所述的低通滤波是指能将被测信号中高频信号进行滤除的滤波算法。根据模数转换速率和谐波分析要求，确定被测信号进行谐波分析的最高次数和低通数字滤波器的各参数，对被测信号进行低通滤波处理，可以消除高次谐波对谐波分析结果的影响。

本发明所述的信号数据长度N是进行一次FFT变换的数据长度，N根据谐波检测精度和计算机或嵌入式系统设备的运行速度综合考虑确定。

本发明所述的Nuttall窗是一类余弦组合窗函数，在非同步采样情况下，不加窗的FFT运算将存在严重的频谱泄露和栅栏效应，检测精度低，加Nuttall窗FFT算法的频谱泄露和栅栏效应能显著降低，能提高基波与谐波检测精度。

本发明所述的双峰插值FFT算法是指在非同步采样情况下，离散谱线与真实频率点存在偏差，通过寻找真实频率点附近的两根峰值谱线，采用多项式拟合的方法求取真实谱线处的频率值、幅值和初相位。

本发明与已有技术相比有以下优点：1、本发明可以对被测信号的基波与谐波分量进行快速、准确检测；2、本发明可以克服被测信号频率波动的影响，实现频率准确测量；3、本发明便于嵌入式系统实现，可以连续、长期对被测信号进行检测。

附图说明

图1是本发明的基波与谐波测量原理框图

图2是四项三阶Nuttall窗离散时域特性图

图3是四项三阶Nuttall窗幅频特性图

图4是FFT双峰插值算法示意图

图5是本发明的基波与谐波测量流程框图

具体实施方式

本发明提出了一种基于Nuttall窗双峰插值FFT的基波与谐波检测方法。以下结合附图作详述，但不作为本发明的限定。

本实施例为可用于实现电力系统谐波电能计量的三相多功能谐波电能表的谐波参数检测部分。

设定本实施例中谐波检测最高次数为21次，模数转换器采用TDK公司生产的71M6513，CPU采用ADI公司生产的BF533处理器，主要参数选择如下：

(1)采样速率：f_s＝2520；

(2)电网基波频率：在50Hz左右波动；

(3)FFT运算数据长度：N＝512。

本实施例的处理流程框图如图1所示，信号经过A/D转换(模数转换)后送入DSP(BF533)进行处理，BF533完成构建离散四项三阶Nuttall窗，进行加窗插值FFT运算，得到各次谐波参数。

本实施例中离散四项三阶Nuttall窗的表达式如下

$w\left(n\right)=\underset{g=0}{\overset{G-1}{\mathrm{\Σ}}}{(-1)}^g{b}_g\cos (2\mathrm{\πn}\·g/N)$

式中，G为Nuttall窗函数的项数；N为FFT运算的数据长度，也是窗函数的长度(n＝1，2，…N-1)；bg应满足以下约束条件

$\underset{g=0}{\overset{G-1}{\mathrm{\Σ}}}{b}_g=1$

$\underset{g=0}{\overset{G-1}{\mathrm{\Σ}}}{(-1)}^g{b}_g=1$

四项三阶Nuttall窗的系数分别是b₀＝0.338946、b₁＝0.481973、b₂＝0.161054、b₃＝0.018027，其离散时域特性如图2所示。图3给出了四项三阶Nuttall窗的幅频特性，四项三阶Nuttall窗具有优良的旁瓣特性，旁瓣衰减速率达到30dB/倍频程，旁瓣电平达到了-83dB。

本实施例中双峰插值算法的示意图如图4所示，k₁和k₂分别表示真实频率点附近的两根谱线，以基频f_s＝50Hz为例，则k₁和k₂将位于40～60Hz频率范围内，找出该范围内的两根峰值谱线，则真实频率点位于两根谱线之间。两根峰值谱线的幅值分别为y₁＝|X(k₁Δf)|和y₂＝|X(k₂Δf)|，设 $\mathrm{\β}=\frac{y_2-y_1}{y_2+y_1}$ 、α＝k₀-k₁-0.5，可知α取值范围为[-0.5，0.5]。

Nuttall窗的频谱幅度函数为

$W\left(w\right)=\underset{m=0}{\overset{M-1}{\mathrm{\Σ}}}{(-1)}^m\frac{b_m}{2}[W_R(w-\frac{2\mathrm{\π}}{N}m)+W_R(w+\frac{2\mathrm{\π}}{N}m)]$

式中， $W_R\left(w\right)=\frac{\sin \frac{\mathrm{Nw}}{2}}{\sin \frac{w}{2}}{e}^{-j\frac{N-1}{2}w}.$

令 $w=\frac{2\mathrm{\π}}{N}k,$ 可得

$W\left(\frac{2\mathrm{\π}}{N}k\right)=\sin \mathrm{\πk}\·e^{-\mathrm{j\πk}}\left[\underset{m=0}{\overset{M-1}{\mathrm{\Σ}}}{(-1)}^m\frac{b_m}{2}\frac{\sin \frac{2\mathrm{\πk}}{N}}{\sin \frac{\mathrm{\π}(k-m)}{N}\sin \frac{\mathrm{\π}(k+m)}{N}}\right]$

将k＝-α±0.5代入上式，因|-α±0.5|≤1，且N一般较大，可得

$\left|W\left(\frac{2\mathrm{\π}(-\mathrm{\α}\±0.5)}{N}\right)\right|\≈\left|\sin \mathrm{\π}(-\mathrm{\α}\±0.5)\right[\underset{m=0}{\overset{M-1}{\mathrm{\Σ}}}{(-1)}^m\frac{b_m}{\mathrm{\π}}\frac{N(-\mathrm{\α}\±0.5)}{{(-\mathrm{\α}\±0.5)}^2-m^2}\left]\right|$

进行多项式拟合逼近，得到基于四项三阶Nuttall窗的计算公式为

α＝H(β)＝2.95494514β+0.17671943β³+0.09230694β⁵

频率计算公式为

f₀＝k₀Δf＝(α+k₁+0.5)Δf

幅值计算公式为

$A=\frac{(y_1+y_2)}{N}(3.20976143+0.9187393{\mathrm{\α}}^2+0.14734229{\mathrm{\α}}^4)$

相位计算公式为

₀＝arg[X(k_iΔf)]+π/2-π(α-0.5(-1)ⁱ)

式中，i＝1，2。

本实施例的基波与谐波参数检测误差结果如下：

表1基波～21次谐波幅值相对误差(％)

谐波次数	1	2	3	4	5	6	7
								误差	-0.0002	0.0003	0.0002	0.0009	0.0002	-0.0001	-0.0001
谐波次数	8	9	10	11	12	13	14
								误差	-0.0009	0.0001	0.0001	-0.0001	-0.0006	0.0001	0.0004
谐波次数	15	16	17	18	19	20	21
								误差	-0.0002	0.0005	-0.0000	0.0005	-0.0002	-0.0008	-0.0001

表2基波～21次谐波初相位相对误差(％)

谐波次数	1	2	3	4	5	6	7
								误差	0.0209	0.0409	-0.0013	0.0039	-0.0004	-0.0041	0.0026
谐波次数	8	9	10	11	12	13	14
								误差	-0.0100	-0.0019	-0.0347	0.0291	-0.0105	0.0097	-0.0254
谐波次数	15	16	17	18	19	20	21
								误差	0.0010	0.0178	0.0001	-0.0132	0.0307	-0.0017	0.0010

Claims (5)

1.一种用于信号处理的基于Nuttall窗双峰插值FFT的基波与谐波参数检测方法，其特征是：对被测信号进行模数转换和低通滤波处理，采用旁瓣电平低、衰减速度快的Nuttall窗对信号进行加窗FFT处理，运用双峰插值算法检测被测信号中的基波和各谐波量。

2.根据权利要求1所述的模数转换是指通过对被测信号(可以是电压或电流信号等)进行高速模数转换，将被测信号转换为数字量。

3.根据权利要求1所述的低通滤波是指能将被测信号中高频信号进行滤除的滤波算法。根据模数转换速率和谐波分析要求，确定被测信号进行谐波分析的最高次数和低通数字滤波器的各参数，对被测信号进行低通滤波处理，可以消除高次谐波对谐波分析结果的影响。

4.根据权利要求1所述的Nuttall窗是一类余弦组合窗函数，在非同步采样情况下，不加窗的FFT运算将存在严重的频谱泄露和栅栏效应，检测精度低，加Nuttall窗FFT算法的频谱泄露和栅栏效应能显著降低，能有效提高基波与谐波检测精度。

5.根据权利要求1述的双峰插值FFT算法是指在非同步采样情况下，离散谱线与真实频率谱线存在偏差，通过寻找真实频率点附近的两根峰值谱线，采用多项式拟合的方法求取真实谱线处的频率值、幅值和初相位。

Images