[发明专利]一种准同步DFT幅值线性修正系数的计算方法

说明书

本发明涉及一种谐波分析技术，具体为一种准同步DFT幅值线性修正系数的计算方法；包括使用准同步DFT线性修正公式对准同步DFT的谐波分析幅值进行线性修正时，线性修正系数M和β_i的计算步骤为：步骤一：Δu在[‑30,30]范围内间隔0.1取值Δμ_i，并按照公式和A₁(Δu)＝产生L＝601组A_1k数据，并记录为的形式；步骤二，确定正整数M＝10；步骤三、按照Δμ_i和M生成数组Δ；步骤四、按照公式B＝(Δ^TΔ)^‑1Δ^TA求解B；步骤五、求均方差步骤六、逐一增加M值到20，重复上述步骤三至五并记录对应的G和B；步骤七、使G最小的M和β_i即为准同步DFT的幅值线性修正系数。

技术领域

本发明涉及一种谐波分析技术，具体为一种准同步DFT幅值线性修正系数的计算方法。

背景技术

谐波分析技术在电能质量监控、电子产品生产检验、电器设备监控等众多领域应用广泛，是进行电网监控、质量检验、设备监控的重要技术手段。目前谐波分析应用最广泛的技术是离散傅里叶变换(DFT)和快速傅里叶变换(FFT)。准同步采样技术和DFT技术相结合的谐波分析技术能够提高谐波分析的精度，其算式为：

式中：k为需要获得的谐波的次数(如基波k＝1，3次谐波k＝3 )；sin和cos分别为正弦和余弦函数；而ak和bk分别为k次谐波的实部和虚部；A_k为幅值；n为迭代次数；W＝nN；γ_i为一次加权系数；为所有加权系数之和；f(i)为分析波形的第i个采样值；N为周期内采样次数。

一次迭代系数γ_i由积分方法、理想周期采样点N和迭代次数n 决定，具体推导过程参见文献【戴先中.准同步采样应用中的若干问题[J].电测与仪表,1988,(2):2-7】。

在工程应用中，谐波分析总是进行有限点的采样和难以做到严格意义的同步采样。这样，在应用准同步DFT进行谐波分析时，就会存在由于截断效应导致的长范围泄漏和由于栅栏效应导致的短范围泄漏，使得分析结果精度不高，甚至不可信。

图1给出了应用准同步DFT对于任一给定实例进行谐波分析的误差图。从图中可以发现，准同步DFT算法的谐波幅值误差极大，基本不可信。

为了解决信号频偏时准同步DFT算法的谐波幅值误差较大的问题，可以采用线性修正的方法对谐波幅值进行修正，具体方法是：使用准同步DFT线性修正公式

对准同步DFT的谐波分析幅值进行线性修正。

发明内容

本发明提供一种准同步DFT幅值线性修正系数的计算方法，来确定准同步DFT的幅值线性修正系数M和β_i。理想的修正系数能够有效减小准同步DFT谐波分析技术的分析误差，获得高精度的谐波分析结果，从而提高基于谐波分析理论的电能质量监控、电子产品生产检验、电器设备监控等领域仪器设备的质量和状态判断的有效性。

本发明的技术方案为：一种准同步DFT幅值线性修正系数的计算方法，包括使用准同步DFT线性修正公式对准同步DFT的谐波分析幅值进行线性修正时，线性修正系数M和β_i的计算步骤为：

步骤一：Δμ在[-30,30]范围内间隔0.1取值Δμ_i，并按照公式

和产生L＝601组A_1k数据，并记录为的形式；

步骤二、确定正整数M＝10；