CN104849551A - 一种谐相角分析方法 - Google Patents

Abstract

本发明涉及一种谐相角分析方法，包括如下几个步骤：(1)采样W+2个采样点数据(W由积分方法决定)；(2)从采样点i＝0应用准同步DFT公式分析W+1个数据获得基波信息和(3)从采样点i＝1应用准同步DFT公式分析W+1个数据获得基波信息和(4)应用公式计算信号的频率漂移μ；(5)从采样点i＝O应用准同步DFT公式分析W+1个数据获得各次谐波信息和(6)应用公式汁算各次谐波的幅相角；(7)应用公式线性修正各次谐波的谐相角。本发明的方法，有助于电能质量监控、电子产品生产检验、电器设备监控等应用谐波分析的领域更加精确的获得各次谐波的幅相角和频率等信息。

Description

一种谐相角分析方法

技术领域

本发明涉及一种高精度的谐相角分析方法。

背景技术

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

$\left\{\begin{array}{c}{a}_k={2F}_{\mathrm{ak}}^n\left(i\right)=\frac{2}{Q}\underset{j=0}{\overset{w}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f(i+j)\cos \left(k\frac{2\mathrm{\π}}{N}j\right)\\ b_k={2F}_{\mathrm{bk}}^n\left(i\right)=\frac{2}{Q}\underset{j=0}{\overset{w}{\mathrm{\Σ}}}{\mathrm{\γ}}_{i}f(i+j)\sin \left(k\frac{2\mathrm{\π}}{N}j\right)\end{array}\right.$

式中：k为需要获得的谐波的次数(如基波k＝1，3次谐波k＝3)；sin和cos分别为正弦和余弦函数；而a_k和b_k分别为k次谐波的实部和虚部；n为迭代次数；W由积分方法决定，采用复化梯形积分方法时，W＝nN；γ_i为一次加权系数；为所有加权系数之和；f(i)为分析波形的第i个采样值；N为周期内采样次数。

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

图1给出了应用准同步DFT对于任一给定实例进行谐波分析的误差图。从图中可以发现，准同步DFT算法的谐相角除了50Hz时其余均误差极大，基本不可信。

发明内容

本发明要解决的技术问题是提供一种高精度的谐相角分析方法，以有效改进准同步DFT谐波分析技术的分析误差，获得高精度的谐波分析结果，从而提高基于谐波分析理论的电能质量监控、电子产品生产检验、电器设备监控等领域仪器设备的质量和状态判断的有效性。

实现本发明目的的技术方案是提供一种谐相角分析方法，包括以下步骤：

(1)等间距采样W+2个采样点数据{f(i)，i＝0，1，...，W+1}(W由所选择的积分方法决定，本发明并不指定某一种积分方法，常用的积分方法有复化梯形积分方法W＝nN、复化矩形积分方法W＝n(N-1)、复化辛普森积分方法W＝n(N-1)/2等，可以根据本发明应用的实际情况来选择合适的积分方法。一般以复化梯形积分方法效果较理想。N为一个理想周期内的采样点数。

(2)从采样点i＝0开始应用准同步DFT公式

$\left\{\begin{array}{c}{a}_k={2F}_{\mathrm{ak}}^n\left(i\right)=\frac{2}{Q}{{\mathrm{\Σ}}_{j=0}^w\mathrm{\γ}}_{i}f(i+j)\cos \left(k\frac{2\mathrm{\π}}{N}j\right)\\ b_k={2F}_{\mathrm{bk}}^n\left(i\right)=\frac{2}{Q}{{\mathrm{\Σ}}_{j=0}^w\mathrm{\γ}}_{i}f(i+j)\sin \left(k\frac{2\mathrm{\π}}{N}j\right)\end{array}\right.,$ 分析W+1个数据获得基波信息和 $F_{b1}^n\left(0\right);$

(3)从采样点i＝1应用准同步DFT公式分析W+1个数据获得基波信息和 $F_{b1}^n\left(1\right);$

(4)应用公式：计算信号的频率漂移μ；

(5)从采样点i＝0开始应用准同步DFT公式分析W+1个数据获得各次谐波信息和

(6)应用公式计算各次谐波的幅相角；

(7)应用公式线性修正各次谐波的谐相角。

等间隔采样是根据进行谐波分析的理想信号的周期T和频率f，在一个周期内采样N点，即采样频率为f_s＝Nf，且N≥64。

所述的采样W+2个采样点数据是根据所选择的积分方法而作相应选择，若采用复化梯形积分方法，则W＝nN；若采用复化矩形积分方法，则W＝n(N-1)；若采用复化辛普森积分方法，则W＝n(N-1)/2。然后根据采样频率f_s＝Nf，获得采样点数据序列{f(i)，i＝0，1，...，W+1}，n为迭代次数，一般n≥3；最后对该数据序列进行谐波分析。

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

为所有加权系数之和。

信号频率的漂移μ是根据相邻采样点基波相角差与理想周期内采样点数N的固定关系而获得的，信号频率的漂移μ也可用于修正基波和高次谐波的频率f₁和高次谐波的频率f_k(f_k＝kμf_s/N)。

本发明具有积极的效果：(1)高精度的谐相角分析结果。如对于图1给定的分析实例，本发明获得的分析精度提高到10^-8级(图2)。

(2)本发明所述的方法从根本上解决了准同步DFT谐相角分析精度低的问题，而无需进行复杂的反演和修正，算法简单。

(3)相对于准同步DFT，本发明所述的谐波分析技术只需要增加一个采样点就解决了准同步DFT分析误差大的问题，易于实现。

(4)应用本发明来改进现有的仪器设备，技术上是可行，并且不需要增加任何的硬件开销就可使分析结果可以提高到10^-8级。

(5)本方法也同样也适用于进行多次迭代而非一次迭代的谐波分析过程，此时只需要把一次迭代分解成多次迭代实现就可以了。一次迭代和多次迭代本质上是一样的，只是在计算时多次迭代进行分步计算，而一次迭代是把多次迭代的过程合并到迭代系数γ_i中一次计算完成，所以本发明同样适用于多次迭代过程。

附图说明

图1为准同步DFT的谐相角分析误差图。

图2为本发明的谐相角分析误差图。

具体实施方式

(实施例1)

本实施例的一种谐相角分析方法，包括以下步骤：

首先，等间隔采样W+2个采样点，以获得被分析信号的离散序列{f(i)，i＝0，1，...，Wq+1}。W的值由积分方法、迭代次数n和理想周期内采样点数N共同决定。

等间隔采样是指：根据进行谐波分析的理想信号的频率(如工频信号频率f为50Hz，周期为20mS)确定采样频率f_S＝Nf，在采样频率f_S的作用下在一个周期内均匀地采样N点。一般地，周期采样点N＝64或以上就能获得较好的谐波分析结果，而迭代次数n＝3～5就能获得较理想的谐波分析结果。

积分方法有复化梯形积分方法、复化矩形积分方法、辛普森方法等多种，可以根据实际情况进行选择。若采用复化梯形积分方法，则W＝nN；若采用复化矩形积分方法，则W＝n(N-1)；若采用复化辛普森积分方法，则W＝n(N-1)/2。

其次，从采样点i＝0开始应用准同步DFT公式

$\left\{\begin{array}{c}{a}_k={2F}_{\mathrm{ak}}^n\left(i\right)=\frac{2}{Q}{{\mathrm{\Σ}}_{j=0}^w\mathrm{\γ}}_{i}f(i+j)\cos \left(k\frac{2\mathrm{\π}}{N}j\right)\\ b_k={2F}_{\mathrm{bk}}^n\left(i\right)=\frac{2}{Q}{{\mathrm{\Σ}}_{j=0}^w\mathrm{\γ}}_{i}f(i+j)\sin \left(k\frac{2\mathrm{\π}}{N}j\right)\end{array}\right.$ 分析W+1个数据获得基波信息和 $F_{b1}^n\left(0\right);$

再次，从采样点i＝1应用准同步DFT公式分析W+1个数据获得基波信息和 $F_{b1}^n\left(1\right);$

再次，应用公式：计算信号的频率漂移μ；

再次，从采样点i＝0开始应用准同步DFT公式分析W+1个数据获得各次谐波信息和

然后，应用公式计算各次谐波的幅相角；

最后，应用公式线性修正各次谐波的谐相角。

本技术领域的普通技术人员应当认识到，以上的实施例仅是用来说明本发明，而并非作为对本发明的限定，本发明还可以变化成更多的方式，只要在本发明的实质精神范围内，对以上所述实施例的变化、变型都将落在本发明的权利要求书范围内。

Claims (7)

1.一种谐相角分析方法，其特征在于包括：包括以下步骤：

(1)等间距采样W+2个采样点数据{f(i)，i＝0，1，...，W+1}；

(2)从采样点i＝0开始应用准同步DFT公式

$\left\{\begin{array}{c}{a}_k={2F}_{\mathrm{ak}}^n\left(i\right)=\frac{2}{Q}{\mathrm{\Σ}}_{j=0}^W{\mathrm{\γ}}_{i}f(i+j)\cos \left(k\frac{2\mathrm{\π}}{N}j\right)\\ b_k={2F}_{\mathrm{bk}}^n\left(i\right)=\frac{2}{Q}{\mathrm{\Σ}}_{j=0}^W{\mathrm{\γ}}_{i}f(i+j)\sin \left(k\frac{2\mathrm{\π}}{N}j\right)\end{array}\right.,$

分析W+1个数据获得基波信息和

(3)从采样点i＝1应用准同步DFT公式分析W+1个数据获得基波信息和

(4)应用公式： $\mathrm{\μ}=N\frac{{\mathrm{tg}}^{-1}\left[\frac{F_{a1}^n\left(1\right)}{F_{b1}^n\left(1\right)}\right]-{\mathrm{tg}}^{-1}\left[\frac{F_{a1}^n\left(0\right)}{F_{b1}^n\left(0\right)}\right]}{2\mathrm{\π}},$ 计算信号的频率漂移μ；

(5)从采样点i＝0开始应用准同步DFT公式分析W+1个数据获得各次谐波信息和

(6)应用公式计算各次谐波的幅相角；

(7)应用公式线性修正各次谐波的谐相角。

2.根据权利要求1所述的一种谐相角分析方法，其特征在于：所述步骤(1)中，等间距采样是根据进行谐波分析的理想信号的周期T和频率f，在一个周期内采样N点，即采样频率为f_s＝Nf，且N≥64。

3.根据权利要求1或2所述的一种谐相角分析方法，其特征在于：所述步骤(1)中，所述的采样W+2个采样点数据是根据所选择的积分方法而作相应选择，然后根据采样频率f_s＝Nf，获得采样点数据序列{f(i)，i＝0，1，...，W+l}；n为迭代次数，n≥3；最后对该数据序列进行谐波分析。

4.根据权利要求3所述的一种谐相角分析方法，其特征在于：所述的采样W+2个采样点数据采用复化梯形积分方法，则W＝nN。

5.根据权利要求3所述的一种谐相角分析方法，其特征在于：所述的采样W+2个采样点数据采用复化矩形积分方法，则W＝n(N-1)。

6.根据权利要求3所述的一种谐相角分析方法，其特征在于：所述的采样W+2个采样点数据采用复化辛普森积分方法，则W＝n(N-1)/2。

7.根据权利要求1所述的一种谐相角分析方法，其特征在于：为所有加权系数之和。