电力系统微机保护中开平方运算的一种新的快速算法.pdfVIP

2005 第十届全国保护和控制学术研讨会 力系统微机保护中开平方运算的一种新的快速算法 曹海欧,袁宇波,李 澄 (江苏省电力科学研究院，南京 2 10036) 摘要：在分析了查表法以及迭代法的基础上，提出了一种新的开平方运算方法－综合法。该运 算方法综合了查表法以及迭代法，运算过程中只需要一次牛顿迭代计算，并且据精度要求可以 定量的确定表格的大小。既具有运算速度快的特点，又具有精度高、占用内存小的特点。解决 了长期以来开平方算法存在耗时长、精度低、存贮数据值范围难以确定的问题。 关键词：开平方；牛顿迭代法；查表法 0 概述 在微机保护装置中，基于交流采样的电力系统微机测量和保护技术应用十分广泛,测量保护中常常 要求计算基波电压、电流的有效值，功率和功率因数以及各谐波分量。比如利用全波傅氏算法计算电压 量、电流量的有效值时，首先根据全波计算公式计算出电压电流量的实部、虚部，再通过对实部和虚部 的平方和进行开平方运算才能得到电流电压量的有效值，在各种微处理 中都没有直接的开平方运算指 令，在微机保护软件设计时，开平方运算的实现主要有以下几种方法： a.采用数值计算法：主要是采用迭代方法，但是选用迭代法存在收敛和发散的问题，另外，即使算 法是收敛的，收敛的速度也会不一样，有的结果只要几次迭代就可以达到所需的精度，有时候迭代几十 次也未必达到需要的精度。 b.调用微处理 系统的函数库：现在的单片机和DSP 的C 编译 都提供了丰富的运行支持库函数， 在C 程序和汇编程序中都可以直接调用它们，从而实现非线性运算。这样做虽然方便，但是库函数经过 编译之后，目标文件十分庞大，而且运算量很大，不能满足实时性高的系统要求。 c.查表法：根据自变量的范围和精度的要求制作一张表格，在计算的时候通过输入值确定表的地址， 再有地址得到相应的结果。查表法的优点在于运算快速，能够满足实时性的要求。但是随着精度要求的 提高，或者输入范围的扩大，所需制作的表格也就越大，极大消耗了系统有限的资源。 本论文中把查表法和迭代法结合起来，充分利用查表法的快速性和迭代法的占用资源小的特点，称 之为综合法。总的思路为：首先对自变量进行移位缩短被开平方数的取值范围，然后在被开方数的取值 范围之内，制作一张表格，根据精度要求，表格内存储了相隔一定间距的点的开方值，把表格事先存放 在微处理 系统中的FLASH 中。计算时，先通过查表法找出计算点所在的区间位置，然后再通过一次牛 顿迭代计算出开平方结果。 1 牛顿迭代法[1] 所谓迭代法就是递归的应用某个公式来决定未知量，并使之逐步逼近真实值的一种方法。对于方程 f(x)=0 构造一个迭代式： xk 1j (xk ) （其中k=0 、1、2 ···） (1) 选择初值x 反复代人该迭代式，得到一个迭代序列 x ，若k 时，x 的极限存在且等于x* ，则称 0 k k 546 曹海欧等 电力系统微机保护中开平方运算的一种新的快速算法 x j (x ) j x 为迭代函数。 k 1 k 是收敛的，否则是发散的， ( ) 牛顿迭代法是一种 单迭代法，在混和法中，我们采用牛顿迭代法。开平方函数f (x )x2 c 0 的根的 牛顿迭代公式为: 1 c X 1 [X ] (2) k k 2 X k 对公式(2)配方之后得： X 1 [X c ]2 (3) k 1 c k 2X