一种基于MIR‑RSD高精度余弦窗插值FFT算法的谐波测量通道校准方法技术

本发明专利技术公开了一种基于MIR‑RSD高精度余弦窗插值FFT算法的谐波测量通道校准方法，包括数据采集模块、对数据进行加窗插值FFT运算的谐波参数检测模块、用于谐波测量通道参数调整的校准模块，其中谐波参数检测模块采用基于MIR‑RSD高精度余弦窗的插值FFT算法，将时域信号加MIR‑RSD窗并进行离散傅里叶变换得到离散频谱Xw(k)，找到离散频谱中频率fN附近的幅度最大谱线k1及其左边谱线k2、右边谱线k3，相应的频谱幅度分别为│X(k1)│、│X(k2)│和│X(k3)│，计算加权比值系数β，从而求出非整数部分频率成分δ，再根据δ值进而求出幅值、频率和相角作为标准量与被校准的谐波测量通道进行校准。

【技术实现步骤摘要】
一种基于MIR-RSD高精度余弦窗插值FFT算法的谐波测量通道校准方法
本专利技术涉及信号相量测量领域，具体是一种基于MIR-RSD高精度余弦窗插值FFT算法的谐波测量通道校准方法。
技术介绍
众所周知，电力系统谐波测量与分析已经成为电力系统领域发展的一个重要研究方向。近几十年来，随着各种电力电子装置(主要是一些非线性设备)的广泛应用，使得电网的谐波污染日益严重，电能质量下降。因此，实时测量和分析电网谐波分量，对于防止谐波危害、提高电能质量十分必要。可是在谐波测量中产生的误差直接影响到测量精度，不能达到预期防止谐波的效果，降低了工作效率。于是，对谐波分析结果进行校准必不可少。而谐波校准装置应运而生，可满足对谐波分析结果的校准。它采用DSP芯片作为处理器，处理信息量大、速度快，能够满足精度要求和实时性要求。它作为标准仪器，用来检测谐波分析结果是否准确。本专利技术采用了三谱线加窗插值FFT算法对信号进行处理，该算法通过分析加窗信号傅里叶变换的频域表达式，利用谐波频点附近的三根频谱的幅值确定谐波谱线的准确位置，进而得到谐波的幅值、频率及相位，推导的三谱线插值修正算法能够进一步提高谐波分析的准确性。
技术实现思路
本专利技术的目的是为了针对现有技术存在的不足，提供了一种操作简单、运算速度快、精度更高的谐波校正方法。本专利技术提出的解决方案为：将时域信号加窗并进行离散傅里叶变换得到离散频谱Xw(k)，找到离散频谱中频率fN附近的幅度最大谱线k1及其左边谱线k2、右边谱线k3，相应的频谱幅度分别为│X(k1)│、│X(k2)│和│X(k3)│，计算加权比值系数β，从而求出非整数部分频率成分δ，再根据δ值进而求出幅值、频率和相角作为标准量与被校准的谐波测量通道进行校准。本专利技术提出的基于MIR-RSD高精度余弦窗插值FFT算法的谐波测量通道校准方法包括以下步骤：步骤一：信号源发出信号，同时送至数据采集模块和被校准的谐波测量通道；步骤二：数据采集模块获取信号并对其模数转换，获得信号模数转换后的数据；步骤三：谐波参数检测模块采用基于MIR-RSD高精度余弦窗三谱线加权插值FFT算法对信号模数转换后的数据进行谐波分析，获得谐波测量结果L1；步骤四：读取被校准的谐波测量通道的谐波测量结果L2，并以谐波参数检测模块的谐波测量结果L1为真值，计算被校准的谐波测量通道的谐波测量结果L2的绝对误差、标准差；步骤五：根据被校准的谐波测量通道的谐波测量结果L2的绝对误差、标准差，确定被校准的谐波测量通道的系统误差分量，计算被校准的谐波测量通道在各次谐波测量时的系统误差校准系数，从而完成被校准的谐波测量通道的校准工作。所述的方法，基于MIR-RSD高精度余弦窗插值FFT算法实现谐波准确分析，其工作步骤包括：步骤一：选择长度为N的窗函数w(n)，n∈[0，N-1]，对信号模数转换后的数据x(n)进行加窗处理，获得加窗后的数据xw(n)＝x(n)w(n)，并对其进行离散傅里叶变换，得到离散频谱Xw(k)，k∈[0，N-1]；步骤二：考虑到非同步采样和干扰的影响，信号的第i次谐波频率fi在离散频谱所对应的位置k0将偏离离散频谱Xw(k)的谱线位置，即包含整数部分频率成分k1及非整数部分频率成分δ＝k0-k1，找到离散频谱中频率fN附近的幅度最大谱线k1及其左边谱线k2、右边谱线k3，相应的频谱幅度分别为│X(k1)│、│X(k2)│和│X(k3)│；步骤三：计算加权比值系数β步骤四：设窗函数w(n)的频谱函数为W(k)，忽略频谱中负频率部分量的影响后的加窗信号的离散频谱函数为则β可改写为构造函数β＝g(δ)，通过多项式拟合，δ值可由反函数δ＝g-1(β)求得；步骤五：根据δ值可得第i次谐波频率频率值fi根据δ值及窗函数的频谱函数可得第i次谐波幅值Ai根据δ值及窗函数的频谱函数可得第i次谐波相角值i所述的方法，窗函数w(n)的时域形式为这里H≥2为窗函数系数，且H的取值由ah的表达式决定，其中所述的方法，MIR-RSD高精度余弦窗被称为MaximumImageinterferenceRejectionwindowswithRapidSidelobeDecayrate，具有最大抗干扰能力和最大旁瓣衰减率，且计算速度快、精度高。所述的方法，校准系数的求取过程，对谐波参数检测模块和被校准的谐波测量通道同时进行M次参数测量，每次测量在基波频率f0下以Δ＝0.1为步长来实施，得到一系列成对的测量结果，即L21、L11，L22、L12，…，L2M、L1M，对M个点采用最小二乘法拟合得到方程L1＝γL2+τ其中γ和τ为校准系数，则L2可以被校准为那么第i次谐波的校准系数的求取，只需改变频率值，即在fi＝if0的条件下，重复上述步骤即可。本专利技术基于MIR-RSD高精度余弦窗插值FFT算法的谐波测量通道校准方法，具有操作简单、计算精准等特点，采用MIR-RSD高精度余弦窗完成插值FFT算法，计算快速、准确，能保证实现准确校准。附图说明图1是本专利技术处理流程的原理框图；图2是本专利技术中基于MIR-RSD高精度余弦窗插值FFT算法实现流程图。具体实施方式本专利技术提出了一种基于MIR-RSD高精度余弦窗插值FFT算法的谐波测量通道的校准方法。以下结合附图作详细说明：本实施例的处理流程的原理框图如图1所示，信号经过包含对其进行降压、滤波、模数转换的数据采集模块后转化为数字量，并送入配备DSP芯片的谐波参数检测模块，在这里将完成对转换后的数据加窗插值FFT运算，得到精确计算后的真值L1，读取被校准的谐波测量通道的谐波测量结果L2，并以谐波参数检测模块的谐波测量结果L1为真值，对系统进行M次测量，得到关于L1、L2的M组数据，根据L1、L2的M组数据，采用最小二乘法拟合计算被校准的谐波测量通道在各次谐波测量时的校准系数，从而完成被校准的谐波测量通道的校准工作。图1中的数据采集模块主要是由调理电路组成，调理电路包括电压互感器、电流互感器、归一化电路和低通滤波器，其中电压互感器、电流互感器将三相标准源产生的信号转换成小电压、小电流信号进入归一化电路和低通滤波器，采用归一化电路的目的是使同一个被测量的信号流经不同互感器以后，仍然能够输出同一电平的信号，以便于系统进行处理，该电路可由一个数字电位计和运算放大器组成的增益可调的放大电路组成，而低通滤波器则是为了滤除高频信号，以确定所测信号的最高频率，因为在加窗插值FFT算法中要采用最高频率来确定采样点的数量。信号经过前端采集模块后进入谐波参数检测模块，由DSP芯片和相应辅助电路组成，DSP芯片型号为TMS320F2812，该模块通过MIR-RSD高精度余弦窗能精确计算出各次谐波参数，得到测量结果L1，与被校准的谐波测量通道的谐波测量结果L2做最小二乘法拟合得到校准系数，从而完成被校准的谐波测量通道的校准工作。下面以求取基波频率f的校准系数为例，具体阐述校准系数的求取过程：设基波频率f0＝50Hz，采用三相标准源HBS1030在基波频率下对频率进行9次测量，其中步长Δ＝0.1，9次频率取值分别为49.6、49.7、49.8、49.9、50.0、50.1、50.2、50.3、50.4，信号同时送至数据采集模块和被校准的谐波测量通道后，获得9组数据如表一本文档来自技高网...

【技术保护点】
一种基于MIR‑RSD高精度余弦窗插值FFT算法的谐波测量通道校准方法，包括获取信号并对其模数转换的数据采集模块、对信号模数转换后的数据进行加窗插值FFT运算的谐波参数检测模块、用于谐波测量通道参数调整的校准模块，其工作步骤包括：步骤一：信号源发出信号，同时送至数据采集模块和被校准的谐波测量通道；步骤二：数据采集模块获取信号并对其模数转换，获得信号模数转换后的数据；步骤三：谐波参数检测模块采用基于MIR‑RSD高精度余弦窗三谱线加权插值FFT算法对信号模数转换后的数据进行谐波分析，获得谐波测量结果L1；步骤四：读取被校准的谐波测量通道的谐波测量结果L2，并以谐波参数检测模块的谐波测量结果L1为真值，对系统进行M次测量，得到关于L1、L2的M组数据；步骤五：根据L1、L2的M组数据，采用最小二乘法拟合计算被校准的谐波测量通道在各次谐波测量时的校准系数，从而完成被校准的谐波测量通道的校准工作。

【技术特征摘要】
1.一种基于MIR-RSD高精度余弦窗插值FFT算法的谐波测量通道校准方法，包括获取信号并对其模数转换的数据采集模块、对信号模数转换后的数据进行加窗插值FFT运算的谐波参数检测模块、用于谐波测量通道参数调整的校准模块，其工作步骤包括：步骤一：信号源发出信号，同时送至数据采集模块和被校准的谐波测量通道；步骤二：数据采集模块获取信号并对其模数转换，获得信号模数转换后的数据；步骤三：谐波参数检测模块采用基于MIR-RSD高精度余弦窗三谱线加权插值FFT算法对信号模数转换后的数据进行谐波分析，其工作步骤包括：步骤1：选择长度为N的窗函数w(n)，n∈[0，N-1]，对信号模数转换后的数据x(n)进行加窗处理，获得加窗后的数据xw(n)＝x(n)w(n)，并对其进行离散傅里叶变换，得到离散频谱Xw(k)，k∈[0，N-1]，窗函数w(n)的时域形式为这里H≥2为窗函数系数，且H的取值由ah的表达式决定，其中步骤2：考虑到非同步采样和干扰的影响，信号的第i次谐波频率fi在离散频谱所对应的位置k0将偏离离散频谱Xw(k)的谱线位置，即包含整数部分频率成分k1及非整数部分频率成分δ＝k0-k1，找到离散频谱中频率fN附近的幅度最大谱线k1及其左边谱线k2、右边谱线k3，相应的频谱幅度分别为│X(k1)│、│X(k2)│和│X(k3)│；步骤3：计算加权比值系数β步骤4：...

【专利技术属性】
技术研发人员：温和，陈洲，黎福海，滕召胜，郭斯羽，胡亮，张军号，张政，
申请(专利权)人：湖南大学，
类型：发明
国别省市：湖南;43