半波傅氏算法及其改进算法的实现

目录
第一章半波傅氏算法 (1)
第二章半波傅氏算法的误差分析 (3)
第三章滤除衰减非周期分量的新算法 (4)
第四章仿真计算 (7)
第五章结论 (9)
参考文献 (10)
半周傅氏算法及其改进算法的实现
摘要
提出一种利用半波傅氏算法消除衰减非周期分量对基波分量影响的快速算法，新算法的数据窗是半个周期的采样值加两个采样点，而其滤波效果远远优于半波傅氏算法。

该算法理论上可以完全消除任意衰减时间常数τ的非周期分量对基波分量的影响。

通过大量的仿真试验表明，新算法滤除衰减非周期分量能力强，计算简单，速度快，具有实际应用价值。

大多数微机保护算法的计算可视为对交流信号中参数的估算过程，对算法性能的评价也取决于其是否能在较短数据窗中，从信号的若干采样值中获得基波分量或某次谐波分量的精确估计值。

目前广泛采用全波傅氏算法和最小二乘算法作为电力系统微机保护提取基波分量的算法。

但由于半波傅氏算法只用半个周期的采样数据，响应快，但滤波能力相对较弱，故只能用于保护切除出口或近处故障。

为使保护快速动作，选择数据窗较短的快速算法就成为关键。

从衰减非周期分量对半波傅氏算法的影响分析入手，提出新的计算方法，可完全滤除衰减非周期分量及奇次谐波分量，以提高其滤波能力。

关键词:微机保护衰减非周期分量半波傅氏算法快速算法
第一章 半波傅氏算法
为了分析衰减非周期分量对半波傅氏算法的影响，设电力系统故障电流有如下形式：
(1-1)
式中 I m (n),φn 分别为n 次谐波的幅值和初相角。

因半波傅氏算法不能滤除偶次谐波，所以设式(1)中n 为奇数，则所得的n 次谐波分量的实部模值an 和虚部模值bn 的时域表达式分别为：
(1-2)
(1-3)
式中 T 为基波分量的周期；ω为基波分量的角频率，ω＝2π/T 。

在计算机上实现时，是对离散的采样值进行计算。

用离散采样值表示的半波傅氏算法为：
(1-4)
(1-5)
式中 k 表示从故障开始时的采样点序号；N 为每个周期的采样点数。

n 次谐波的幅值I m (n)和初相角φn 为：
(1-6)
(1-7) 假设暂不考虑输入信号(如式(1-1)的形式)中的衰减非周期分量，根据式(1-4)、式(1-5)利用半波傅氏算法得到的理论值为：
(1-8)
(1-9)
第二章半波傅氏算法的误差分析
如果输入信号中包含衰减非周期分量，将使半波傅氏算法的计算结果产生误差，具体分析如下：
(2-1)
(2-2)
令
(2-3)
(2-4)
由式(2-1)、式(2-2)可知，当输入信号中包含有衰减非周期分量时，I0≠0，α≠0,则w a≠0,w b≠0。

从而看出，n次谐波的实部和虚部与理论值相比，存在误差w a 和w b。

因此，消除w a和w b是将半波傅氏算法应用于快速保护的关键之一。

第三章滤除衰减非周期分量的新算法
为了全部使用故障后的采样值，取k≥N/2，同时，为了使新算法的推导更趋于精确，下面以时域形式介绍新算法的推导过程。

a.取第一个数据窗，使t∈［0，T/2］，利用半波傅氏算法有：
(3-1) 令
(3-2) 则式(3-1)可以简化为：
(3-3)
，取第2个数据窗，使t∈
b.取延时ΔT为一个采样周期时间T
s
［ΔT,(T/2)+ΔT］,有：
(3-4) 令
(3-5) 在理论上，移动的数据窗大小(即ΔT)可任意确定，但为了提高算法的计算
较合适。

一旦确定了每个周期的采样速度以达到快速计算的目的，ΔT选取为T
s
点数N，ΔT也就随之确定。

同时，若谐波次数n和延时ΔT确定,k
a ,k
b
就成为
两个常数。

则式(3-4)可化简为：
(3-6)
c.延时2ΔT,取第3个数据窗，使t∈［2ΔT,(T/2)+2ΔT］，有：
(3-7) 由式(3-3)、式(3-6)、式(3-7)可以看出，3个方程组中只有5个未知数，
而为了校正衰减非周期分量对半波傅氏算法的影响，只要计算出w
a 和w
b
的值，
即可对半波傅氏算法由于衰减非周期分量引起的误差进行校正，式中的未知数A，B和e-αΔT只需作为中间变量，没有必要求出。

其计算过程如下：利用式(3-3)、式(3-6)、式(3-7)，先消除A，B两个中间变量。

令：
Q＝a
n ′-k
a
a
n
+k
b
b
n
(3-8)
R＝b
n ′-k
a
b
n
-k
b
a
n
(3-9)
X＝a
n ″-2k
a
a
n
′+a
n
(3-10)
Y＝b
n ″-2k
a
b
n
′+b
n
(3-11)
这里的Q，R，X，Y值可根据采样值实时计算出。

所以由式(3-8)～式(3-11)得：
w a /w
b
＝X/Y (3-12)
w b Q-w
a
R＝k
b
(w2a+w2b) (3-13)
由式(3-12)和式(3-13)得：
(3-14)
式(3-14)是由于衰减非周期分量对半波傅氏算法产生的影响数据。

则由式(2-1)和式(2-2)可得，消除衰减非周期分量对半波傅氏算法影响的校正量a
c和
n
c应为：
b
n
分析新算法的整个计算过程可知，半个周期后第3个采样间隔的计算量较大，但其计算时间仅约80μs，完全能够满足实时控制的要求。

第四章 仿真计算
通过设置下列输入信号：
i(t)＝50e -t/τ+50sin(ω1t+φ1)+ 15sin(3ω1t)+10sin(5ω1t)
对新算法进行仿真计算，并与半波傅氏算法和全波傅氏算法进行了比较，其结果见表1。

这里取τ＝30ms,ω1＝100π,φ1＝30°，其对n 次谐波分量的计算程序流程图如图1。

图4-1 程序流程图
开始 输入信号数据
去第1数据窗，利用半波傅氏公式计算n a 和n b 延时S T ，去第2数据窗，再计算'n a 和'n b 延时2S T ，去第3数据窗，再计算"n a 和"n b 计算中间变量Q 、R 、X 、Y 数值 计算a w 和b w
输出计算结果 n n nc w a a -= wb b b n nc -=
开始
表4-1 仿真计算结果
从仿真计算的输入信号可以看出，本算例输入信号中含有衰减非周期分量的初值为100%基波幅值，之所以设置这样大的衰减非周期分量初值(在实际中属于比较严重情况)，就是为了人为增大衰减非周期分量对滤波算法的影响，来检验新算法的有效性。

从表1可见，通过与全波傅氏算法和半波傅氏算法的比较，本文提出的新算法具有很高的计算精度。

算法
幅值
相角 计算值
误差值 /(%)
计算值
误差值 /(%)
全波傅氏算法 57.41243 14.82486 28.96597 -3.44676 半波傅氏算法 99.43623 98.87245 27.50386 -8.32044 新算法
49.99940
0.00000
29.99998 0.00004
第五章结论
本文在分析衰减非周期分量对半波傅氏算法产生的影响的基础上，介绍了一种新算法，不仅保留了原来半波傅氏算法的功能，又增添了对衰减非周期分量的过滤作用。

新算法所采用的数据窗仅为半个周期的采样值加两个采样点，计算简单，速度快，精度高；同时其滤除衰减非周期分量的能力又不受衰减非周期分量时间常数大小的限制。

特别适合于需要快速动作的继电保护
参考文献
[1]．杨奇逊，黄少锋．微型机继电保护基础（第3版）．北京：中国电力出版社，2007．
[2]．许正亚．变压器及中低压网络数字式保护．北京：中国水利水电出版社，2004．
[3]．陈皓．微机保护原理及算法仿真．北京：中国电力出版社，2007．
[4]．杨新民，杨隽琳．电力系统微机保护培训教材（第2版）．北京：中国电力出版社，2008．
[5]．罗士萍．微机保护实现原理及装置．北京：中国电力出版社，2001
[6]．丁书文，张承学，龚庆武，等．半波傅氏算法的改进——一种新的微机保护交流采样快速算法．电力系统自动化，1999，23(5)：18-20．。