基于牛顿迭代法的高精度快速开方算法.pdfVIP

第28卷第3期 电力自动化设备 V01．28No．3 ElectricPowerAutomation Mar．2008 2008年3月 Equipment o 基于牛顿迭代法的高精度 快速开方算法 张小呜1，一，李永新1 (1．南京理工大学机械工程学院，江苏南京210094； 2．江苏工业学院计算机科学与x-4r--．Y_系，江苏常州213164) 摘要：针对全波傅氏算法中的开平方运算较耗时影响微机保护瞬动性问题，提出一种选取接近定 点数开方真值的牛顿迭代初值的方法。该方法利用数字信号处理器(DSP)的移位指令、256个单元 的查表技术和DSP的硬件乘法器，通过1次查表和1．2次乘法运算，就能快速确定迭代误差小于 2’9的迭代初值。在TIDSP集成开发平台上。运行牛顿迭代开平方汇编程序。运算结果表明：该算法 达到2娟以上迭代精度．且占用内存小。非常适合带有硬件乘法器的嵌入式微处理器实现。 关键词：开平方；迭代；微机保护；DSP 744 中图分类号：TM 文献标识码：A 基于交流采样的微机保护测量算法中最常用的 束条件满足： 是并联补偿全波傅氏算法[1】。该算法计算谐波幅值 占=lXn+l-X。I≤2。16 (1) 时，要进行谐波实部平方和虚部平方的开平方运算。 2牛顿迭代初值选取新算法 因嵌入式微处理器无专门的开方指令。目前常见的开 平方算法主要有牛顿迭代(Newton．Raphson)算法[21 假设具有图l数据格式的被开方数A最大为 和逐位循环算法C3-6]。在微处理器应用场合，通常用1FFFFFF．FFFF H．用下列二进制通式表示： 牛顿迭代法实现开方运算。牛顿迭代法每次产生1 A=口。a，l…alao．口一la一2…口一m+l口一。、(2) 个开方数．直到前后2次迭代的误差小于希望的精 将式(2)展开为以2为底的权展开和形式[12】： 度为止。若取固定的迭代初值且被开方数较大时， 迭代非常慢。一般10次以上才能达到精度。若迭代 a一2×2—2+…+n—m+1×2一刖。1+口一m×2一肌(3) 初值越接近开方结果，则迭代收敛的速度越快[7．¨]。 其中，O≤乃≤24，1≤m≤16，口。=‘0’或‘1’，口。= 这里研究了定点数直接牛顿迭代开方初值的快 ‘0’或‘l’。 速选定算法，减少牛顿迭代次数，提高定点数牛顿迭 将式(3)中最高含‘1的二进制位数n的权力作为 代开方结果的速度．使定点数开方牛顿迭代快速算 公因子提到公式最左边．即在a。=‘1’的情况下，有 法能充分发挥定点CPU快速处理定点数的优势。 A=2”×(1+口。1×2-1+…+o】×2—4“+口o×2一“+口一1× 1 被开方定点数数据格式定义 、／万一=2n／2×(1+口。一lx2—1+…+口lx2-n+I+ 微机保护系统测量的电流、功率等变化范围较 aox2咄+Q一1x2-I-n+口一2x2-2-n+…+ 大，为了达到较高的牛顿迭代开方精度1．5X10～，设 口一。+l×2一”+1一”+口一。×2一”””)172 (4) 置被开方定点数数据格式为8Byte，其中6Byte(即 从公式(4)可看出，、／丁的值由2个乘积项组 3个字)为整数部分，2 成，若牛顿迭代初值仅取第1项2n