微机继电保护傅立叶算法研究

《微机继电保护傅立叶算法研究》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、BatchDoc-Word文档批量处理工具微机继电保护傅立叶算法研究张迪，王维庆新疆大学电气工程学院,乌鲁木齐（830008）E-mail:zhangdi2434@yahoo.cn摘要：全波傅里叶算法可以计算信号谐波分量，半波傅里叶算法把运算的数据窗减少了一半，运算量比全波算法减少，但二者都会受到衰减直流分量的影响。本文在传统傅里叶算法的的基础上，每次求和时把衰减直流分量减去，从而去除了衰减直流分量的影响。不论全波还是半波，相角都不能直接算出，需要对计算出的值修正后才能得到正确结果。此外联合Matlab

2、和VisualBasic可以方便的对算法仿真。关键词：全波傅里叶算法,半波傅里叶算法,衰减直流分量,1．引言全波傅里叶算法以傅里叶级数为基础，可以准确计算出故障函数中的各次谐波，但是传统算法有两个明显缺点。（1）衰减直流分量的影响。（2）相角如果直接套用公式，结果错误。针对这两大缺点，本文提出两种不同的解决方法。第一，滤去衰减直流分量。挑出故障函数的三个周期的最大值和最小值（共六个），用三个最大值计算出上包罗曲线的三个参数A,B,C,从而得到上包罗曲线函数，同理也可以利用最小值算出下包罗曲线，再把A,B

3、,C同a,b,c相加求和除以二，就得到了衰减直流分量的三个参数，就相当于得到衰减直流分量函数。那么在用傅里叶算法求和时从y中把衰减直流分量减掉，就可以达到将其滤掉的目的。通过Matlab的仿真结果来看，滤波效果还是很明显的。第二，相角的问题，直接套公式的结果绝大部分是错的，所以把期望值和计算值在-2π到2π做了比较，通过曲线发现其中有一定规律性，按照他的差值规律，将计算值同期望值作拟合，最后得到很好的拟合结果，通过仿真结果可以看出。传统半波傅立叶算法也有以上的两个缺点，虽然计算时间短于全波，但是以不能计

4、算偶次谐波的幅值为代价。用Matlab编程时为了很直观的看出计算结果和仿真图形，参考利用VB和Matlab混合编程的文章之后，本文选取了一种通用的方法，即把Matlab作为VB的一个ActiveX控件加入到VB中，实现在VB中调用Matlab，使用VB界面录入数据，后台计算交给Matlab，主要实现对全波和半波的图形仿真比较，体现各自算法的特点。2．全波傅立叶算法编程修正2.1全波傅立叶算法原理传统傅立叶算法的原理是基于对周期信号的傅立叶分解.并把对连续信号的积分改成对离散量的求和,公式如下[7]-4-

5、BatchDoc-Word文档批量处理工具nU(t)=∑Uksin(kωt+φk)k=1(1)-4-BatchDoc-Word文档批量处理工具资助来源：教育部博士点基金:大型风力发电机组电能品质监测与智能控制的研究(20060755001)-4-BatchDoc-Word文档批量处理工具-4-BatchDoc-Word文档批量处理工具2ak=TT∫0u(t)cos(kωt)dt(2)-4-BatchDoc-Word文档批量处理工具2bk=TT∫0u(t)sin(kωt)dt(3)-4-Batc

6、hDoc-Word文档批量处理工具U=a2+b2(4)-4-BatchDoc-Word文档批量处理工具kkφ=arctank−bk(5)-4-BatchDoc-Word文档批量处理工具kak-4-BatchDoc-Word文档批量处理工具其中U(t)是信号,k是谐波次数,T是信号的周期,公式离散化后得:Uk是谐波幅值,φk是谐波初相角。把-4-BatchDoc-Word文档批量处理工具-4-BatchDoc-Word文档批量处理工具Nak=2∑ucos(nk2π)(6)-4-Batch

7、Doc-Word文档批量处理工具nNn=1N-4-BatchDoc-Word文档批量处理工具Nbk=2∑usin(nk2π)(7)-4-BatchDoc-Word文档批量处理工具nNn=1N-4-BatchDoc-Word文档批量处理工具2πfN为一个周波的采样次数,连续函数中相角ωt离散化后为n2πN,因为ω离散化后为-4-BatchDoc-Word文档批量处理工具。N2.2传统全波算法首先,根据以上的原理利用Matlab编写了仿真程序。电源函数为：y=20+20*sin(100*pi*t

8、)+8*sin(200*pi*t+pi/3)+10*sin(300*pi*t)+2*sin(400*pi*t+pi/6)+5*sin(500*pi*t+pi/4)+sin(600*pi*t)；在Matlab中运行传统傅里叶算法,且无衰减直流分量的程序结果如图1。图1无衰减直流分量的频率图Fig1thefrequencymapofDCcomponentwithoutdecaying从图上可以看出傅立叶算法很准确的算出各谐波的幅值且不受直流分量