傅里叶算法的采样电流计算

傅里叶算法的采样电流计算

*******广西大学*******摘要：微机继电保护是用数学运算的方法实现故障的测量、分析和判断的。通过全波傅立叶算法可用于求出各次谐波分量的幅值和相角，并具有一定的滤波作用。本文探讨了傅氏算法在电力系统中的应用。介绍了全波傅立叶算法的基本原理。通过仿真验证了该算法的实用性。关键词：微机继电保护；电力系统；算法引言在微机保护装置中，首先要对反映被保护设备的电气量模拟量进行采集，然后对这些采集的数据进行数字滤波，再对这些经过数字滤波的数字信号进行数学运算、逻辑运算，并进行分析判断，最终输出跳闸命令、信号命令或计算结果，以实现各种继电保护功能。这种对数据进行处理、分析、判断以实现保护功能的方法称为算法。目前广泛采用全波傅氏算法和最小二乘法作为电力系统微机保护提取基波分量的算法。傅立叶算法可用于求出各谐波分量的幅值和相角，所以它在微机保护中作为计算信号幅值的算法被广泛采用。实际上，傅立叶算法也是一种滤波方法。分析可知，全周傅氏算法可有效滤除恒定直流分量和各正次谐波分量。傅里叶算法原理一个周期函数满足狄里赫利条件，就可以将这个周期函数分解为一个级数，最为常用的级数是傅里叶级数，傅氏算法的基本思路来自傅里叶级数，

即一个周期性函数可以分解为直流分量、基波分量及各次谐波的无穷级数，如L”I=z[bcos（nwt）+asm（nwt）]TOC\o"1-5"\h\z（t） n 1n 1 （1.1）n=0式中"1表示基波角频率；an和bn分别是各次谐波的正弦和余弦的幅值，其中比较特殊的有：b表示直流分量，a,b表示基波分量正、余弦项的幅0 1 1值。根据傅氏级数的原理，可以求出a、b分别为2T・•,a=fTisin（nwt）dtnT0（t） 1 （1.2）\o"CurrentDocument"『 2T・ /b=fTicos(nwt)dt（1.3）nT0(t) 1（1.3）于是n次谐波电流分量可表示为i(t)=bcos(nwt)+asln(nwt)n n 1n 1 (1.4)据此可求出n次谐波电流分量的有效值和相角为I=官nYba=arctann（1.5）n a（1.5）n其中a、b可用梯形积分法近似求出为

1 NT・二1 NT・二2knn_.a= [2zism ]k=1八1「・bn=N[l0+(1.6), 2kmicos +1]k-i(1.7)式中N——基波信号1周期采样点数，k——第k次采样值l0,lN k=0和k=N时的采样值求出基波分量（n=1）的实部和虚部％4，即可求出信号的幅值。当米样频率为600Hz时，取wT=3Oo（N=12），基波正、余弦的系数如下表所示，于是可得到式（1.6）和（1.7）的采样计算公式为i(1.8)（1.9）a- [(i+1-1-1)+3(i+1-1-1)+2(i-1)](1.8)（1.9）1 121 5 7 11 % 2 4 8 10 3 9,1 b= [i+1—i—i+1+1)+,3(i—i—i+1)—2i]1 120 2 4 8 10 12 1 5 7 10 6湾由i,i,i，…，ik—0』,2,...,12□-工1中012 12 时刻的采样值。基波正弦和余弦的系数（N=12时）k012345678gto1112sln{k—)N0121楚212aL逅2-1.迈L■70cost止叩-)N1g21I0£-7互212L-7aL2V321实例")=100sin(wt+：6)=100sin(314t+?利用前面叙述的傅里叶算法进行计算，采样周期为12点（N=12）,则间隔时间为At=002S。12所以，i=100sin(314Atk+E(t) 6=100sin(314x%2k)k=0,1,2,...,1211、 ，:+T)(—=0.1666671)6计算得到采样值表:k0123456h5086.58926599.99998686.64232950.0919240.132721-49.862009k789101112h-86.509486-99.999775-860721742-50.229703-0.29198649.723891把采样值分别代入式(1.8)和式(1.9)得a=86.701959七=42.642633由匕,«的值代入式(1.5)得'I=68.321387<a=26.189348)=0.145496HL1i(t)=寸2x68.321387sin(3141+0.145496n)=96.621032sin(3141+0.145496n)i(t)结语本文通过对全波傅立叶算法原理的介绍，并通过仿真验证了12点全波傅立叶算法在电力系统数据采集中的应用。仿真结果表明：全波傅立叶算法在电力系统中具有一定的实用性可求出各谐波分量的幅值和相角，而且，可以有效地滤除恒定直流分量和各