专利名称：基于Kaiser窗双谱线插值FFT的基波与谐波检测方法
技术领域：
本发明涉及一种用于信号中基波与谐波参数检测的方法。具体是一种基于 Kaiser窗双谱线插值FFT (快速傅立叶变换)的基波与谐波参数检测方法，属于 信号处理技术领域，也可用于其它信号的谱分析。
背景技术：
以电力信号检测为例，随着电力电子技术的发展，非线性设备在电力系统 中得到广泛应用，非线性设备带来的谐波问题对电力系统安全、稳定、经济运 行构成的威胁日趋严重。基波与谐波的精确测量可为电网电能计量、谐波潮流 计算、设备入网检测、电力系统谐波补偿与抑制等提供科学依据。
基波与谐波分量检测方法一般有(1)基于频域分析的FFT方法，其特点 是电网频率波动时，检测精度较低；(2)滤波器检测谐波方法，其特点是难以 获得理想频率特性；(3)基于瞬时无功功率理论的方法，其特点是计算量大， 处理复杂；(4)基于神经网络理论和小波变换的方法，其特点是计算量大，难 以在嵌入式系统中实现。
基于FFT的基波、谐波分析算法，易于在DSP (数字信号处理器)上实现， 是当今应用最广泛的一种谐波分析方法。在电力系统中，电网电压、电流畸变 导致基波频率变化，由于非同步采样造成FFT算法存在频谱泄露和栅栏效应等 问题，使得基波与谐波检测的精度受到影响。减少FFT算法的频谱泄露和栅栏 效应影响、提高电测量中基波和各谐波检测精度是电测量信号分析和电能质量 管理中的难题。
已有专利文件"测量工厂谐波的方法和测量仪"(200310105446.2)，"计量 工厂谐波的方法和计量仪"(200310105446.2)，"电力系统谐波定量计量方法和计量仪"(98110414.2)等。其发明的目的是精确定量地测量电力系统中谐波电 压、电流的动态特性。已有技术的不足之处是已有技术提供了一些基波与谐 波测量的电气设备设计方法，但由于谐波具有多样性、随机性和多态性等特征， 基被频率波动造成的频谱泄露和栅栏效应依然存在，因此实时、高精度的基波 与谐波测量分析难以实现。
经文献检索发现，对信号采用加窗插值的法，可有效实现实时、高精度的 基波与谐波测量。目前，国内外学者提出了一系列加窗插值FFT算法，但均属 于采用基本窗函数和广义余弦窗函数对信号加权，对于动态信号分析效果受到 窗函数固定旁边性能的制约。而Kaiser窗可定义一组可调的窗函数，其主瓣能 量和旁瓣能量的比例近乎最大，且可自由选择主瓣宽度和旁瓣高度之间的比重。 针对此特点，何益宏等学者提出以实时无功理论为基础，采用Kaiser窗对信号 取样(何益宏，卓放等.Kaiser窗在谐波电流检测中的应用，电网技术，2003， 27 (1): 9-12);吴峻等学者提出基于Kaiser窗的频率测量法(吴峻，吴崇昊. 基于Kaiser窗的改进傅氏测频算法，江苏电机工程，2008， 27 (4): 17-19)， 以上文献均未涉及基于Kaiser窗插值FFT处理的相关处理方法。
发明内容
为克服已有技术、方法的不足，本发明的目的在于提供基于Kaiser窗双谱 线插值FFT方法，该方法能快速、准确检测信号中基波和各谐波分量。 设包含多项整数次谐波的时域被测信号为
式中，//表示谐波的项数；Q代表谐波的次数；当/2=1时，^=1，力、A、 ^分别 为基波的频率(Hz)、幅值(V)和初相位(。)；当好1时，Q、 A、 ^分别为第/2项谐波的次数、幅值(V)和初相位(。)。信号;c(f)经过采样率为厶的数据采集系统后，得 到离散序列
根据傅立叶变换的定义，信号< )的连续傅立叶变换为
= I;义("K"一 (3)
在嵌入式系统上或计算机上实现时，数据总是有限长的，即相当于信号JC(")
被一个窗函数vK")("二o,i，…w-i)截短为w点长序列。加窗后，信号频谱由原
来的谱线变为以/o为中心向两边扩展的连续谱，信号频谱旁瓣之间的相互干扰，
能量泄露到整个频带，即为频谱泄露。
对％(/)以离散频率间隔△/ = / 进行离散化，可得到y(A:A/)
(yt二0，l,…w-l，代表离散频谱线的序号)。由于电网频率并非恒定，存在一 定波动，加上被测信号中可能存在间谐波(即非整数次谐波)，因而难以严格做到 同步采样。非同步采样造成峰值频率很难正好位于离散谱线频点上，即
(&为非整数)，导致的对信号峰值点的观测偏差，即为栅栏效应。
基于以上分析，本发明提出的基于Kaiser窗的基波与谐波检测方法，其特 征在于对被测信号进行模数转换和低通滤波处理，采用可自由选择主瓣与旁 瓣衰减之间的比重、设计实现灵活、旁瓣电平低、衰减速度快的Kaiser窗对信 号进行加窗FFT运算，运用双谱线插值算法检测被测信号中的基波和各谐波量。 它包括下列步骤
(1)模数转换与低通滤波
高速模数转换器将被测信号(模拟量)40转换为数字量，采样频率为A， 该数字量作为基波与谐波检测处理的输入量。根据谐波分析精度要求，确定需检测的最高次谐波次数及频率/皿，设计截止频率为稍大于/皿的低通数字滤 波器，输入量经过数字低通滤波器后，得到去除高频信号分量的输入量JC(")。
<formula>formula see original document page 6</formula>(4)
式中，//为所含最高谐波的次数。
(2)构建离散Kaiser窗(根据设计精度要求选择参数P ) Kaiser窗是一种最优化窗，其离散表达式为
<formula>formula see original document page 6</formula>(5)
式中，Io(灼是第1类变形零阶贝塞尔函数，"是窗函数的形状参数，为任意正数 值；W为FFT运算的数据长度，也是窗函数的长度("-l,2，…W-l)。 Kaiser窗的连续频谱函数，其频域表示为
<formula>formula see original document page 6</formula>
(6)
图2分别给出了 "=
时Kaiser窗的幅频特性图，其 中，^=0， Kaiser窗己退化为矩形窗。当y -ll时，旁瓣衰减速率达到18分贝/ 倍频程，旁瓣电平达到了-82.3分贝。随着yS值增大，Kaiser窗函数频谱的旁瓣 峰值越小，渐近衰减速率越大，同时，主瓣相应也加宽。因此，具体y9值可根 据实际精度要求选定，体现了应用Kaiser窗设计实现的灵活性。 (3)基于Kaiser窗的FFT运算
根据检测精度的要求，确定FFT运算的数据长度W。为便于嵌入式系统实现，本发明考虑只含有单次谐波存在的情况，将去除高频分量的输入信号量I(")
进行加Kaiser窗FFT运算，可得到加窗后信号的离散傅立叶变换的表达式为
<formula>formula see original document page 7</formula>
式中，r/产&A/为第A项谐波的峰值频率，ff(0是Kaiser窗的连续频谱函数。 (4)双谱线插值算法 不失一般性，设需要测量的为第《^/Z)项谐波，为简单起见，忽略其余各次 谐波对第/项谐波的泄漏影响，此时，式(7)变为
<formula>formula see original document page 7</formula>(8)
由于栅栏效应，第/项谐波的峰值频率A:,A/很难正好位于抽样频点上，即 & 一般不是整数。设在峰值频点附近抽样得到的幅值最大和次最大谱线分别为 禾tl &， fe^ A ^,.2 (=&+1)，这两条谱线的幅值分别为Y尸P^'沐i4/)l和 力4^V'(A:。A/)l，设Xy2-yi)/(y2+yi)、"-^ 一 0.5，，可知or取值范围为[-0.5， 0.5]，可得以《为自变量，以义为因变量的函数关系/1=/ 0)，有
<formula>formula see original document page 7</formula>
其反函数为a吐"(A)。由A可求出参数"，则被测信号的频率计算公式为
<formula>formula see original document page 7</formula>
被测信号的幅值计算公式为<formula>formula see original document page 8</formula>被测信号的初相位计算公式
A = arg[I二 (W1] +冗/ 2 - arg,2;r(A: - <) / 7V)〗(i2)
由此可计算出被测信号基波与谐波的真实频率和各项参数。 本发明所述的模数转换是指通过对被测信号(可以是电压或电流信号等)
进行高速模数转换，将被测信号转换为数字量。
本发明所述的低通滤波是指能将被测信号中高频信号进行滤除的滤波算
法。根据模数转换速率和谐波分析要求，确定被测信号进行谐波分析的最高次
数和低通数字滤波器的各参数，对被测信号进行低通滤波处理，可以消除高次
谐波对谐波分析结果的影响。
本发明所述的信号数据长度iV是进行一次FFT变换的数据长度，7V根据谐 波检测精度和计算机或嵌入式系统设备的运行速度综合考虑确定。
本发明所述的Kaiser窗是一类最优化窗函数，可自由选择主瓣宽度和旁瓣 高度之间的比重。在非同步采样情况下，不加窗的FFT运算将存在严重的频谱 泄露和栅栏效应，检测精度低，加Kaiser窗双谱线插值FFT算法的频谱泄露和 栅栏效应能显著降低，能提高基波与谐波检测精度。
本发明所述的双谱线插值FFT算法是指在非同步采样情况下，离散谱线与 真实频率点存在偏差，通过寻找真实频率点附近的两根峰值谱线，采用多项式 拟合的方法求取真实谱线处的频率值、幅值和初相位。
本发明与己有技术相比有以下优点1、本发明可以对被测信号的基波与谐 波分量进行快速、准确检测；2、本发明可以克服被测信号频率波动的影响，实现频率准确测量；3、本发明便于嵌入式系统实现，可以连续、长期对被测信号
进行检测。
图1是本发明的基波与谐波测量原理框图 图2是"=
时Kaiser窗的离散时域特性图，可见，^值越 大，Kaiser窗越窄。图3分别给出了 "=
。 Kaiser窗的频谱幅度函数为
<formula>formula see original document page 10</formula> (14)
令w-2兀A:/W，可得
<formula>formula see original document page 10</formula> (15)
进行多项式拟合逼近，得到基于Kaiser窗的计算公式为a = =4.26201515义+ 0.50674586;i5 + 0.22516483A5
十0.13687018义7 (16)
频率计算公式为
/0 = &4/" = (" + &+0単/ (17)
幅值计算公式为
j =(力"2) (3.81160858+ 0，79465173 2 + 0.08673112a4
十0.006796術6) (18)
相位计算公式为
& = arg[Z"M/)] + W2 - (AT —1)(" - 0.5);r/W (19)
本实施例的基波与谐波参数检测误差结果如下
表1 基波~21次谐波幅值相对误差(％)
谐波次数1234567
误差9.2 E-28,6 E-29.2 E-28.7 E-29.1 E-28.8 E-29.1 E-2
谐波次数891011121314
误差8.8 E-29.1E-29.0 E-29.1 E-28.9 E-29.1 E-29.0 E-2
谐波次数15161718192021
误差9.0 E-29.0 E-29.0 E-29.1 E-29.0 E-29.1 E-29.0 E-2
表2基波~21次谐波初相位相对误差(％)
谐波次数1234567
误差6.0 E-59.3 E-54.9 E-63.9 E-69.0 E—72.3E-62.0E-7
谐波次数891011121314
误差4.2 E-62.3 E-62.6 E-653 E-51.8 E-61.2 E-51.8 E-5
谐波次数15161718192021
误差5.0 E-62.3 E-61.3E-62.犯一42.8 E—4—4.4E—56.7 E—5
本实施例还采用了其他经典窗进行了实验，结果表明，Kaiser窗具有良好的频谱泄漏抑制作用，采用基于Kaiser窗双谱线插值FFT的基波与谐波检测方法 能够在较短数据长度的情况下对谐波进行准确分析，有利于高精度谐波分析算 法的嵌入式系统实现。
权利要求
1. 一种用于信号处理的基于Kaiser窗双谱线插值FFT的基波与谐波参数检测方法，其特征是对被测信号进行模数转换和低通滤波处理，采用可自由选择主瓣与旁瓣衰减之间的比重、设计实现灵活、旁瓣电平低、衰减速度快的Kaiser窗对信号进行加窗FFT处理，运用双谱线插值算法检测被测信号中的基波和各谐波量。
2. 根据权利要求1所述的模数转换是指通过对被测信号(可以是电压或电流信 号等)进行高速模数转换，将被测信号转换为数字量。
3. 根据权利要求1所述的低通滤波是指能将被测信号中高频信号进行滤除的滤 波算法。根据模数转换速率和谐波分析要求，确定被测信号进行谐波分析的最 高次数和低通数字滤波器的各参数，对被测信号进行低通滤波处理，可以消除 高次谐波对谐波分析结果的影响。
4. 根据权利要求1所述的Kaiser窗是可自由选择主瓣与旁瓣衰减之间的比重的 一类最优化窗函数，在非同步采样情况下，不加窗的FFT运算将存在严重的频 谱泄露和栅栏效应，检测精度低，加Kaiser窗双谱线插值FFT算法的频谱泄露 和栅栏效应能显著降低，能有效提高基波与谐波检测精度，且设计实现灵活。
5. 根据权利要求1述的双谱线插值FFT算法是指在非同步采样情况下，离散谱 线与真实频率谱线存在偏差，通过寻找真实频率点附近的两根峰值谱线，采用 多项式拟合的方法，求出实用的插值计算式，由此求取真实谱线处的频率值、 幅值和初相位。
全文摘要
本发明公开了一种基于Kaiser窗双谱线插值FFT(快速傅立叶变换)的基波与谐波检测方法对被测信号进行模数转换和低通滤波处理，采用可自由选择主瓣与旁瓣衰减之间的比重、设计实现灵活、旁瓣电平低、衰减速度快的Kaiser窗对信号进行加窗FFT运算，运用双谱线插值算法检测被测信号中的基波和各谐波量。包括下列方面(1)模数转换与数字低通滤波；(2)基于Kaiser窗的FFT运算；(3)双谱线插值和多项式拟合求取基波与谐波参数。
文档编号G01R23/16GK101441233SQ20081014378
公开日2009年5月27日 申请日期2008年12月2日 优先权日2008年12月2日
发明者周良璋, 张红瑛, 博 曾, 李聪聪, 杨步明, 毛群辉, 和 温, 滕召胜, 一 王, 王璟珣, 贺静丹, 高云鹏, 魏双双 申请人:湖南海兴电器有限责任公司