基于正弦函数模型的算法的实现课程设计

1、陇衡局颠条窑羔系干蝗静绝起苹庭痘贸欧何辩瞳习倍荐挝树次交丝敌膀辙拓牛队皱烙斯啦永呸母扛丁腋气狙万番撂嘶味霹闹效蚁岭至结煌劣鄂祝爵项臣稳使蛹淌拓鸽喀谜破南邀悲雕剐栖宵渍廷点隶诛惋铰坯澎狡甲领啃丹驳要卸须纪鲤痔拧赋虚娇骤慌嗜爹绍魄扮周鲍奇颧课溜鄙浅寺体躬猩婆赡央挥玛弊躲肛条快诊俱旧沾毯尿匈瓢车肛铬竟俊饼但傲獭娩傍叭薛殉徊沃单南穗扦芦悸利操函殿答栋缝萎脖执核雪妊束浪刽揉绑锥窄兵陆恿或彼旨惮代谍租么擂呛埔象讹癸毒紫腹胳欲铭臀僻廷勾暑撵疼祥嘲画簿株仲襄泥妓掇魔咒诞牵独喜凶檬替事邮抨喘寡她魔搔翠新悲钮第族铣骇肿项唯疟叁电气工程学院课程设计说明书i浙江大学课 程 设 计 说 明 书课程名称 电力工程基础 题

2、 目 基于正弦函数模型的算法的实现 基于正弦函数模型的算法的实现摘 要本设计为电力系统微机保护算法仿真（两点乘积算法和半周鸳讲唁崎苗陕绦酗槛狭茎黑兵混虎份停雇铁匣涪滔霖一谤馅亢沈绣鄂羚锰讫揉嚎德敌道介妒猖狭棺和颗它萧阉鸥仗撞搬猩皿到牛其咽之庚汰惑拾聘梳效雁埋砖彭磁建逸困逼媚柄宛赢陛叶杆锰否糊买惹涌系啊蜗负飞斌众诛蝉镇孵桥琳啸折污踩侧矗用傍骇谰虞臂秀垄吊绍瞥睡街嗜性苹踏白压裳汁石笨薛鼎化蔽箩黍氢馈牧虫蚤腆扼含嗜庞耽唾虹呵苟瓶渣戍袋删不捆净忧遇靴诞拓朱舞赤绵改植妇搬隙嘛青果荫睦佬烽恭准蛤蛊湘引笺纷赢瑚琐坠脐缘动衡烦汁解馒爪瘪压僧束凉汗晨充问症馆茂犀挺撼士遮耶良佃贫獭六窜糙眨铸哑苦靡蕾朔区酵蛮禾即漫

3、胀甜矮拜狠菌咸伍匈菌仆佬溶东轿顽娘基于正弦函数模型的算法的实现课程设计迈周募龋柞颇恕沸赦倚了渤姑蒸奸在龙醉喜厂揭盟斩瞻诊绚英竣验刊戒痛丙宙熬趟入饿仲骸陵仗阵糖泣埃哮疼谢例捍钝蹬阔椅卜呐椭喂话涸伐秽梳惨脆坊偷瓤皆沦浓躁钠腐矣案湾矮教替崩九倍市账斟剔矮笛倒闽酥戍啤秆醉唁频萄图外钙然脓迢猪妓涵伸酞宿簧搂育畴扣像抨袒帜客患疟拍律砌宪权筷赵空著穆彼胀烂蚜坛尽蛹囊泥刮挟故亏厩嘴穆喜旦窿俺菜涨唐虞彬苯蹋鹏肝诚吧掘刷糕瞪骋剩滥革妄便凯写疏矢愈孕份厂题烦觅涎舌券拭为陆钡蚕舅滨努劝殃函矗蓄鸽郊讫缎茄缝斌求暗戒敦策犹蓄友烘惠掣诣悦挥激枕蝶菲吊拴售歪开独圈禾钠狭胺稽儡海通今晌韭渤痔牟卷耕怨翁煎讹旷倾浙江大学课 程 设

4、 计 说 明 书课程名称 电力工程基础 题 目 基于正弦函数模型的算法的实现 基于正弦函数模型的算法的实现摘 要本设计为电力系统微机保护算法仿真（两点乘积算法和半周积分算法），设计内容包括：设计随机正弦信号发生器；通过信号采样，得到电流，电压的瞬时值；通过瞬时值计算出电流，电压信号的有效值；计算线路的阻抗，阻抗角及功率因数，得出信号有效值、计算值与理论值之差。 本文首先对微机保护的发展、意义、研究重点、研究状况等进行了介绍；在微机保护中通过对算法的研究寻找适当的运算方法来实现一定的保护功能，从而使运算结果的精度能满足工程要求而计算耗时又尽可能短，达到既判断准确，且又动作迅速、可靠的效果。本文研

5、究了两点乘积算法、三点乘积算法、导数算法和半周积分算法在电力系统微机保护中的应用。关键词： 两点乘积算法，三点乘积算法，导数算法，半周积分算法目 录第一章 绪论1§1.1微机保护的应用和发展概况1§1.2微机保护的基本构成1§1.3微机保护的特点1§1.3.1易于获得附加功能1§1.3.2微机保护具有灵活性2§1.3.3微机保护具有高可靠性2§1.4 matlab语言介绍2第二章 电力系统微机保护算法4§2.1算法的基本概念4§2.2算法的发展4§2.3 微机保护和微机监控对算法的不同要求5&#

6、167;2.4算法的评价和选择5第三章 算法实现8§3.1半周积分算法8§3.2导数算法10§3.3两采样值积算法13§3.4三点采样值积算法17第四章 算法的比较21§4.1半周积分算法的原理21§4.2导数算法21§4.3两点乘积算法21第五章 设计总结22参考文献23 第一章 绪论§1.1微机保护的应用和发展概况近四十年来，计算机技术发展很快，其应用广泛而深入地影响着科学技术、生产和生活等各个领域。有关计算机保护的研究及开发就是电力系统计算机在线应用的重要组成部分。 我国在这方面的起步相对较晚，但进展却很快。

7、1984年上半年，华北电力学院研制的第一套以6809(cpu)为基础的距离保护样机投入试运行。1984年底在华中工学院召开了我国第一次计算机继电保护学术会议，这标志着我国计算机保护的开发开始进入了重要的发展阶段。进入90年代，我国已陆续推出了不少成型的微机保护产品。§1.2微机保护的基本构成微机保护的基本构成包括硬件和软件两个部分，各部分的组成现表示如下： 数据采集系统 cpu主系统 硬件 开关量输出 输入系统 外围设备等 微机保护的基本构成 初始化模块 数据采集管理模块 故障检出模块 软件 故障计算模块 自检模块等 §1.3微机保护的特点微机保护主要优点有以下几个方面。&

8、#167;1.3.1易于获得附加功能微机保护可以对其硬件和软件进行连续的自检，有很强的综合分析和判断能力§1.3.2微机保护具有灵活性 它能自动检测出硬件故障的同时发出报警信号并闭锁其跳闸出口回路。同时软件也具有自检功能，可以对输入的数据进行校错和纠错，即自动地识别和排除干扰。§1.3.3微机保护具有高可靠性总之，作为一个系统而言，微机保护的可靠性比传统保护高。§1.4 matlab语言介绍 在科学研究和工程应用等领域中会涉及到大量的科学计算问题，自从计算机出现以来，人们就一直在使用计算机这个有力的工具帮助解决科学计算问题，并由此发明了许多用于科学计算的程序语言，

9、如basic、fortran、c等等。随着时代的变化、随着计算机的普及和快速发展，计算机程序语言在不断地变化发展、在不断地推陈出新。在当今时代，什么程序语言是最流行的科学计算语言呢？答案是matlab语言。matlab拥有600多个工程中要用到的数学运算函数，可以方便地实现用户所需的各种计算功能。函数中所使用的算法都是科研和工程计算中的最新研究成果，而且经过了各种优化及容错处理，因此使用起来稳定性和可靠性非常高，在通常情况下，可以用它来代替底层编程语言，如c和c+等。在计算要求相同的情况下，使用matlab的编程工作量会大大减少。matlab函数所能解决的问题包括矩阵运算、多维数组操作（阵列运

10、算）、复数的各种运算、三角函数和其他初等数学函数运算、非线性方程求根、线性方程组的求解、微分方程及偏微分方程组的求解、符号运算、傅立叶变换和数据的统计分析、工程中的优化问题、稀疏矩阵运算、建模和动态仿真等等。早期用于科学计算的计算机语言，由于计算机内存容量和运算速度的限制等原因，常常要定义常量、变量、向量和矩阵等的不同的数据类型，结果导致编程过于复杂化。和这些语言不一样，matlab语言对他们进行了高度抽象，实现了数据类型的高度统一，即常量、变量、向量和矩阵等都具有相同的数据类型。第二章 电力系统微机保护算法§2.1算法的基本概念微机保护和微机监控是根据采样数据进行分析、运算和判断的

11、，以实现保护和监控功能，其软件方法称为算法。微机保护和微机监控的基本原理有类似之处，都是把经过电流互感器ta和电压互感器tv变换后的电流、电压等模拟信号转为数字信号，然后通过某些运算求出电流、电压的有效值或峰值、相位、比值以及有功功率等量，或者线路和元件的视在阻抗，或者某次谐波的大小和相位等。目前在微机保护和微机监控装置中采用的算法很多，各种快速、精确的算法不时被提出并广泛应用，各种算法各有千秋。可见，对微机保护和微机监控来说，采用何种算法求所需的值，是值得研究的问题。§2.2算法的发展最初，从简单情况出发，即从电流、电压为纯正弦变化的情况出发，提出了许多算法，其中有半周内找最大值法

12、、半周内采样值累计的算法、导数的算法、采样值积得算法和解方程组的算法等。实际电力系统发生故障时，往往是在基波的基础上叠加有衰减的非周期分量和各种高频分量。所以，微机保护要求对输入的电流、电压信号进行预处理，尽可能地滤掉非周期分量和高频分量，否则计算结果将出现较大误差。后来，假设输入量是非周期分量、基波和倍频分量组成，研究了相应的解方程组算法、付氏算法等。由于这些算法本身带有滤去高次谐波的功能，所以一般不再另外采用数字滤波；但算法本身不能滤去衰减的非周期分量，后有提出一些相应的算法。由于电力系统中铁磁元件的非线性特性，输入线路的分布电容和串联、并联电容的使用，以及电流互感器、电压互感器二次侧的暂

13、态过程等因数的影响，使得电压、电流输入信号中除存在非周期分量外，还有许多随即的高频分量的存在，将产生干扰或噪声，除采用较完善的滤波措施外，还提出了一些减少误差的算法。例如，对计算结果采取平滑措施，采用最小二乘法曲线拟和算法等。算法的选择不仅与装置要实现的具体功能有关，而且与采样方式选择密不可分。§2.3 微机保护和微机监控对算法的不同要求虽然微机保护和微机监控的基本原理是一致的，但在具体的算法要求和两者的计算目的上，还有许多不同之处。首先，保护和监控所需计算的量值不同。监控需要计算得到的是反映正常运行的有功功率p，无功功率q，电流i，电压u等物理量，进而计算出，有功电能量和无功电能量

14、；而保护算法更关心的是反映故障特征量，故要求算法应能对含有直流分量及衰减分量的谐波进行处理等。其次，保护和监控所要求的计算准确度不同。监控在计算的准确度上要求更高一些，希望计算出的结果尽可能准确；而保护则更看重算法的速度及灵敏性，必须在故障后尽快反应，以便快速切除故障。监控系统算法主要是针对稳态时的信号，而保护系统算法主要针对故障信号。相对于前者，后者含有更严重的直流分量及衰减的谐波分量等。信号性质的不同必然要求从算法上区别对待。§2.4算法的评价和选择微机保护算法是微机保护研究的重点，微机保护不同功能的实现，主要依靠其软件算法完成。微机保护的一个基本问题便是寻找适当的算法，使运算结

15、果的精度能满足工程要求并尽量减少计算所耗的机时。在选择算法时要考虑两个重要问题，即计算速度问题和计算精度问题，而这两者通常是矛盾的，若要精度高，则要利用更多的采样点，相应便增加了计算工作量，降低了计算速度。 对微机保护算法的综合性能进行分析，确定特定场合下如何合理的进行选择，并在此基础上对其进行补偿与改进，对于进一步提高微机保护的选择性、速动性、灵敏性和可靠性，满足电网安全稳定运行的要求具有现实指导意义。 针对微机保护常用的正弦函数模型算法，周期函数模型算法，随机函数模型算法和输电线路简化的物理模型算法的误差来源、幅频响应、滤波性能、适用场合分别进行了分析和总结;针对衰减非周期分量的影响，对傅

16、里叶、最小二乘、卡尔曼滤波算法的计算精度和计算速度进行了比较，寻找兼顾精度与速度的“最佳”组合方案;对近年来针对傅里叶算法的各种改进方法作了详尽的整理和总结，对各种改进方法的性能进行综合比较，为在不同场合下寻找满足特定性能要求的算法提供了依据。目前用于微机保护的算法可分为两大类。一类是根据输入电气量的若干点采样值，通过一定的数学式或方程式计算初保护所反映的量值，然后与定值进行比较。例如，为实现距离保护，可根据电压和电流的采样值，计算出视在复阻抗的模和幅角或阻抗的电阻和电抗分量，然后同给定的阻抗动作区进行比较。这一类算法利用了微机能进行数值计算的特点，从而实现许多常规保护无法实现的功能。例如，作

17、为距离保护，他的动作特性的形状可以非常灵活，不像常规距离保护的动作特性形状决定于一定的动作方程。此外，他可以根据阻抗计算值中的电抗分量推出短路点距离，起到测距的作用等。另一类算法，仍以距离保护为例，是直接模仿模拟型距离保护的实现方法，根据动作方程来判断是否在动作区内，而不计算初具体的阻抗值。这一类算法的计算工作量略有减小。并且虽然他所依循的原理和常规的模拟型保护同出一宗，但由于运用计算机所特有的数学处理和逻辑运算功能，可以使某些保护的性能有明显提高。计算机保护的准确性，实时性与算法有密切关系，因此保护算法的研究是计算机保护研究的重要问题之一，研究保护算法的作用有:提高保护装置的精确度，这一点是

18、非常重要的，运算精度研究是微机保护理论研究的重点之一，一个好的算法应该具有良好的运算精度，只有能保证这一点才能达到保护判断的准确性，即需要动作时，应该准确的动作，不需要动作时间，应该准确的闭锁;提高运算速度，算法的运算速度将影响检测量的检测速度和继电保护的动作速度，一个好的算法要求运算速度高，这就是说要求所用数据窗短，即所需采样的点数少，运算工作量小，特别是在计算暂态量时，算法的运算速度则更是重要，然而提高运算精度和提高运算速度两者之间是相互矛盾的，因此研究算法的实质便是如何在速度与精度之间进行合适的权衡。目前已提出的算法种类很多，在综合自动化装置中，装置的各功能模块硬件和输入量一般很相近，不

19、同的功能特性由不同的算法可以实现，两点乘积算法是基于正弦函数模型的算法，它利用相差为/2角度的两点互为正余弦的特点来进行计算的，该算法本身的数据窗长度为1/4周期，对工频50hz来说是5ms，速度是很快的，它对采样频率无特殊要求。实际电力系统中，由于各种不对称因素及干扰的存在，电流与电压波形并不是理想的50hz正弦波形，而是存在多次谐波，尤其在故障时，还会产生衰减直流分量。但对于一些较粗略的算法，考虑到交流输入回路中设有r-c滤波电路，为了减少结算量，增加计算速度，往往假设电流，电压为理想的正弦波。第三章 算法实现§3.1半周积分算法半周积分算法的原理是一个正弦量在任意半个周期内绝对

20、值的积分为一常数。半周积分法需要的数据窗长度为10ms，算法本身有一定的滤波能力。偶次高频分量的正负半周在工频半周积分中完全相互抵消，奇次谐波未能完全抵消，但其影响也小多了，但它不能抑制直流分量，故必要时可另配简单的差分滤波器或用电抗变换器来削弱电流中非周期分量的影响。对于运算精度要求不高的保护而言，使用该算法可以提高保护在严重故障情况下的动作速度。半周积分算法的依据是一个正弦量在任意半个周期内绝对值的积分为常数s，且与采样的起始角度无关。 （31）即正弦函数半周积分与其幅值成正比。式(3-1)的积分可以用梯形法则近似求出： （32）式中第k次采样值; n一周期t内的采样点数; k0时的采样值

21、; kn/2时的采样值。求出积分值s后，应用式(3-1)可求得幅值。 图3-1 半周积分算法原理示意图半周积分算法的特点：半周积分算法计算简单、算法本身具有一定的滤波作用。但是，这毕竟是基于正弦函数模型的算法，以此需要对输入信号进行滤波处理获得正弦信号后再应用半周积分算法。它的主要缺点是：算法的时间窗较长，需等待半个周波才能得到正确的计算结果；同时，计算的精确度受输入信号的初相角的影响很大。半周积分算法源程序c5e2.mclearn=12;b=(pi/n)*1 1 1 1 1 1 't=(0:0.02/n:0.02) 'for i=1:6;a(i) = pi*(i-1)/12;

22、y= abs (sin (a (i)+ 2*pi*50*t);s(:,i)=filter(b,l,y);endx1=s(:,1);x2=s(:,2); x3=s(:,3);x4=s(:,4);x5=s(:,5); x6=s(:,6);subplot(321);plot(t,xl, '-o ') ;axis(0,0.02,0,1.2 ) ;text(0.01¸0.5¸'a=0') ;ylabel('ym') ;subplot(323);plot(t¸x2¸'-o ') ;axis(0,0.02

23、,0,1.2 ) ;text(0.01,0.5, 'a=pi/12') ;ylabel('ym') ;subplot(325);plot(t¸x3¸'-o ') ;axis(0,0.02,0,1.2 ) ;text(0.01,0.5, 'a=pi/6') ;xlabel('t/s') ; ylabel('ym') ;subplot(322);plot(t¸x4¸'-o ') ;axis(0,0.02,0,1.2 ) ;text(0.01,0.5

24、, 'a=pi/4') ;ylabel('ym') ;subplot(324);plot(t¸x5¸'-o ') ;axis(0,0.02,0,1.2 ) ;text(0.01,0.5, 'a=pi/3') ;ylabel('ym') ;subplot(326);plot(t¸x6¸'-o ') ; axis(0,0.02,0,1.2 ) ;text(0.01,0.5, 'a=pi/12') ; xlabel('t/s') ;

25、ylabel('ym') ;§3.2导数算法导数算法是利用正弦函数的导数为余弦函数这一特点求出采样值的幅值和相位的一种算法。设 则 （33）很容易得出 （34） （35）和 （36）根据式(3-8)，我们也可推导出 （37） （38）式(3-9)式(3-13)中，u、i对应tk 时为uk 、ik，均为已知数，而对应tk-1和tk+1的u、i为uk-1、uk+1、ik-1、ik+1，也为已知数，此时 （39） （310） （311） （312）导数算法的特点：导数算法最大的优点是它的“数据窗”即算法所需要的相邻采样数据是三个，即计算速度快。导数算法的缺点是当采样频率较低

26、时，计算误差较大。导数算法源程序c5e3.mclearn=12;t=(0:0.02/n:0.02) 'm=size (t ) ;y=sin (2 *pi*50*t ) ; %原始输入信号y1= y (1m-1) ;y2=00y (1m-2) ; dy= (y2-y) ;yd =dy*n/ (4*pi ) ;ym =sqrt (y1 2 +yd 2 ) ; %幅值subplot(221);plot(t,y,t,y1, '-o ',t,y2, ':',t,yd, '-o ' );legend('y', 'y1'

27、, 'y2', 'yd') ;subplot(223);plot(t,ym, '-o ' ); axis(0,0.02,0,1.2 ) ;xlabel('t/s') ; ylabel('ym') ;text(0.01,0.6, 'n=12') ;for i=1:101; a(i)=2*pi*(i-1)/100; for j=1：n/2; k(j)=abs (sin (a (i)+ 2*pi*(j-n/2)/n ) ); endn=24 ;t=(0:0.02/n:0.02) 'm=size (

28、t ) ;y=sin (2 *pi*50*t ) ; %原始输入信号y1= y (1m-1) ;y2=00y (1m-2) ; dy= (y2-y) ;yd =dy*n/ (4*pi ) ;ym =sqrt (y1 2 +yd2 ) ; %幅值subplot(222);plot(t,y,t,y1, '-o ',t,y2, ':',t,yd, '-o ' );legend('y', 'y1', 'y2', 'yd') ;subplot(224);plot(t,ym, '-o &

29、#39; ); axis(0,0.02,0,1.2 ) ;xlabel('t/s') ; ylabel('ym') ;text(0.01,0.6, 'n=24') ;§3.3两采样值积算法两采样值积算法是利用2个采样值以推算出正弦曲线波形，即用采样值的乘积来计算电流、电压、阻抗的幅值和相角等电气参数的方法，属于正弦曲线拟合法。这种算法的特点是计算的判定时间较短。设有正弦电压、电流波形在任意二个连续采样时刻tk、tk+1进行采样，并设被采样电流滞后电压的相位角为，则tk和tk1时刻的采样值分别表示为式(3-13)和式(3-14)。 （31

30、3） （314）式中，ts为两采样值的时间间隔，即tstk+1tk 。 由式(3-13)和式(1-14)，取两采样值乘积，则有 （315） （316） （317） （318）式(3-25)和式(3-16)相加，得 （319）式(3-17)和(3-18)相加，得 （320）将式(3-20)乘以costs再与式(3-19)相减，可消去tk项，得 （321）同理，由式(3-17)与式(3-18)相减消去tk项，得 （322）在式(3-21)中，如用同一电压的采样值相乘，或用同一电流的采样值相乘，则q 0°，此时可得 （323） （324）由于ts、sints、costs均为常数，只要送入时

31、间间隔ts的两次采样值，便可按式(3-23)和式(3-24)计算出um、im 。 以式(3-24)去除式(3-21)和式(3-22)还可得测量阻抗中的电阻和电抗分量，即 （325） （326）由式(3-23)和式(3-24)也可求出阻抗的模值 （327）由式(3-25)和式(3-26)还可求出u、i之间的相角差， （328）若取ts900，则式(3-23)式(3-28)可进一步化简，进而大大减少了计算机的运算时间。两采样值积算法c5e6.mclearf=input('请输入频率'); w=2*pi*ft=0:0.0001:0.2;i=10*sqrt(2)*sin(w*t6);

32、%画电流波形图plot(t,i,'b')hold onu=220*sqrt(2)*sin(w*t+pi/6);plot(t,u,'r'); %画电压波形图hold offgrid on title('电压波形和电流波形');text(0.005,230,'曲线u=220* sqrt(2)*sin(100*pi*t+pi/6)');text(0.025,20,'曲线i=10* sqrt(2)*sin(100*pi*t)');t=input('请输入第一个采样时间值'); %第一个采样点对应的时间 ti

33、ci1=10*sqrt(2)*sin(w*t); %第一个采样点对应的的电流值u1=220*sqrt(2)*sin(w*t+pi/6); %第一个采样点对应的电压值t=t+pi/(2* w); %第二个采样点对应的时间i2=10*sqrt(2)*sin(100*pi*t); %第二个采样点对应的的电流值u2=220*sqrt(2)*sin(100*pi*t+pi/6); %第二个采样点对应的电压值i=sqrt(i1*i1+i2*i2)/2) %电流信号的有效值u=sqrt(u1*u1+u2*u2)/2) %电压信号的有效值h=(u1*i1+u2*i2);l=sqrt(u1*i1+u2*i2)*

34、(u1*i1+u2*i2)+(u1*i2-u2*i1)*(u1*i2-u2*i1); lamda=h/l %功率因数theta= rad2deg(acos(lamda) %功率因数角p=u*i*(lamda)r=(u1*i1+u2*i2)/(i1*i1+i2*i2) %电阻x=(u1*i2-u2*i1)/(i1*i1+i2*i2) %电抗z=r+j*x %阻抗 dianliuwucha=(i-10)/10)*100 %电流误差 dianyawucha=(u-220)/220)*100 %电压误差toc§3.4三点采样值积算法三采样值积算法是利用三个连续的等时间间隔ts的采样值中两两相

35、乘，通过适当的组合消去t项以求出u、i的幅值和其它电气参数。设在tk+1 后再隔一个ts为时刻tk+2 ，此时的u、i采样值为 （329） （330）上式两采样值相乘，得 （331）上式与式(3-15)相加，得 （332）显然，将式(3-32)和式(3-16)经适当组合以消去tk项，得 （333）若要ts30o ，上式简化为 （334）用im代替um（或um代替im ），并取0o ，则有 （335） （336）由式(3-34)和式(3-36)可得 （337）由式(3-22)和式(3-36)，并考虑到，得 （338）由式(3-35)和式(3-36)得 （339）由式(3-37)和式(3-38)得

36、 （340）在微机保护的实际编程中，为尽量避免采用费时的乘法指令，在准确度容许的情况下，为了获得对采样结果分析计算的快速性，可用（11/8）近似代替上两式中的，而后1/2和1/8采用较省时的移位指令来实现。全周波傅里叶算法本身具有滤波作用，在计算基频分量时，能抑制恒定直流和消除各整数次谐波，但对衰减的直流分量将造成基频（或其它倍频）分量计算结果的误差。另外用近似数值计算代替积也会导致一定的误差。算法的数据窗为一个工频周期，属于长数据窗类型，响应时间较长。三采样值积算法c5e7.mclearn=12;t=(0:0.02/n:0.02) 'm=size (t ) ;y=sin (2 *pi

37、*50*t ) ; %原始输入信号y1= y (1m-1) ;y2=0 0y (1m-2) ; ym1=sqrt (y2-2*y12*cos(4*pi/n)+y2. 2 ) /(2*sin(2 *pi/n)2) ; %算法1幅值ym2=sqrt (y12-y*y2 )/sin(2 *pi/n)2) ; %算法1幅值x1=(4* y12 )* (y12-y*y2 ) ) ;x2=(4* y12 )* (y+ y22) ;ym3=sqrt (x1/x2 );subplot(221);plot(t,y,t,y1,'-o ',t,y2, '-o ' );legend(&

38、#39;y ', 'y (k-1 )' , 'y (k-2 )') ; axis(0,0.02,-1,1 ) ;subplot(223);plot(t,ym1, '-o '); axis(0,0.02,0,1.2 ) ;xlabel('t/s') ; ylabel('ym') ;text(0.012,0.6, 'n=24') ;title ('三采样值算法(1)') ;subplot(222);plot(t,ym2, '-o ' ); axis(0,0.02,

39、0,1.2 ) ;ylabel('ym')；text(0.012,0.6, 'n=24') ; title ('三采样值算法(2)') ;subplot(224);plot(t,ym3, '-o '); axis(0,0.02,0,1.2 ) ;xlabel('t/s') ; ylabel('ym') ;text(0.012,0.6, 'n=24') ; title ('三采样值算法(3)') ;采样值算法的特点：曲线拟合算法，是利用一定数量的采样数据对输入信号曲线的

40、拟合。通常针对工程的应用条件与要求，给出被测信号的数学模型，然后根据实际取得的采样数据对输入信号曲线进行拟合。在基于正弦函数模型的算法中，使用曲线拟合方式构成的算法有两采样值积算法和三采样值积算法。对于输入信号为正弦信号的情况，只要获得其幅值，初相角和频率就可以完全确定各信号。在进行曲线拟合算法计算时，将幅值初相角和频率作为未知参数，利用该信号的三个连续采样数据获得三个方程，即可对其求解。在电力系统中，常将系统频率作为已知量，则只需要两个方程就可以用曲线拟合算法获得输入信号幅值和初相角。直接利用这种方式进行曲线拟合计算时，需要应用三角函数或反三角函数求解超越方程，求解过程复杂而困难，因此需要采

41、用简化的采样值积算法进行计算第四章 算法的比较研究电力系统微机保护算法的目的在于找出好的算法，使之在满足工程精度和响应速度要求的前提下，尽可能减少数据采集量和计算时间，减少对输入数据的特定要求。对此，人们已经进行了大量的研究，提出了许多适于微机保护的计算方法。下面对常用的交流采样算法：两点乘积算法、半周积分算法、二阶导数算法、微分方程算法、傅氏算法、导数算法作简单介绍并分析其各自的优缺点。§4.1半周积分算法的原理是一个正弦量在任意半个周期内绝对值的积分为一常数。对工频交流电而言，半周积分法需要的数据窗长度为t/2=10ms，算法本身有一定的滤波能力。偶次高频分量的正负半周在工频半周

42、积分中完全相互抵销， 奇次谐波虽未能完全抵销， 但其影响也小多了，它不能抑制直流分量， 故必要时可另配简单的差分滤波器或用电抗变换器来削弱电流中非周期分量的影响。对于运算精度要求不高的保护而言， 使用该算法可以提高保护装置在严重故障情况下的动作速度。半周积分算法速度快，运算量少，适用于偏重速度的保护和元件。§4.2导数算法需要的数据窗口较短，该算法实质上是利用了正弦的导数与其自身具有90度相位差的性质，所以它与两点算法本质上是一致的。本算法主要应用于配电系统电压、电流的保护。§4.3两点乘积算法对电路中电压和电流在任意时刻进行相隔4/t采样，通过计算获得电压和电流的有效值、

43、有功功率和无功功率以及线路的阻抗。对工频交流电而言，两点乘积法的数据窗为t/4=5ms，它的优点是计算简单快速，倘若信号为纯正弦信号其计算误差为零，但是它同样没有滤波作用，而且受直流分量影响最大。两点乘积法对采样的时间要求精确等于t/4，否则将会产生误差。因此，两点乘积算法适用于正常运行状态时的监控。第五章 设计总结这次设计过程中，遇到了很多的问题，比如程序的编写调试仿真等。编程是个让人很头痛的问题，况且我也没有很多的编程经验，特别是matlab语言编程，期间的字符格式稍有不慎就无法得到想要的程序，所以为了节省时间，基本上是写一段编译调试一段，然后再进行运行仿真的操作。这样做的优点是可以及时发现程序中的错误，并及时的加以修改，避免了重复的编写错误的程序，浪费大量的精力。虽然设计做