数字化继电器中递推DFT算法的几点实际问题.docVIP

数字化继电器中 递推DFT算法的几点实际问题 作者：Hatem A. DarwishMagdy Fikri 翻译：测控06-3班 张维 摘要：在本文中作者将会呈现数字化继电中可逆DFT算法的现实问题，并着重于DFT滤波器的不稳定性导致的放大错误。在估计幅度和相位放大的判断方面也很重要。本文针对新的放大错误的解决方案进行深入测试，，利用DSP为基础实现结果校正。本文将对继电保护工程师尤其是关注与数字化继电保护工程师有帮助。 关键词：可逆DFT，向量估计，DSP，微机继电保护 1.引言 现代数字化继电器是目前被公认的用于电力系统保护高效化和自动化的核心装置。继电器被设计有包括控制、测量、模拟、干预和事件记录的功能。因此，继电器中的测量数值算法应当尤其被关注，特别是在以微机和DSP实现的强电控制的继电气设计中。 DFT滤波器在各种参考书和工业报告中都被广泛用于基波和谐波测量。DFT滤波器主要有两种形式：非递推算法和递推算法。与非递推算法相比，运算量的减少使递推算法得到普遍的使用。有部分技术文章已经强调了一下在输入信号中递归算法对于动态干扰和系统干扰有抑制作用。虽然普遍认为两种算法都有相同的表现，但其实只是理论上的。事实是，已经证明，DFT滤波器被用于浮点微处理时会偏离理想结果，实际操作中输入信号的动态范围和频域会产生一系列相位和幅度的放大错误。 98年上半年开始这种现象在很多的电力系统继电保护实验室运用浮点DSP板研究项目中出现。实验结果是，即使是理想的的50HZ正弦输入信号也会产生一个极低频率的放大误差，并随着频率的正负偏离增大，误差也随之增加。此后这个问题虽然被解决，但却缺乏一个合理的证明[3]。在[3]中，放大错误的原因被归结为数学运算错的积累，与递推滤波器中反馈[5]-[7]相似。但是，递推DFT的等值框图和带反馈递推滤波器的内在区别减弱了这一证明的有效性。本文着重于强调在继电安全方面对这一错误分析和补救原则。 本文中已指明递推DFT的几种可能的放大误差，实际说明这些放大误差的计算和应用的建议，用等效的Simulink模型证明实验结果，并提出和评价六个可供选择的解决方案。 2.递推DFT的数学公式 为了提取???入信号的基波成分（Xj），可用以下几个公式[1]-[3]达成一个完整周期的递推DFT滤波器。 其中，N是周期数，Sj和Cj是DFT表达式的正弦和余弦部分。输入信号的有效值Fj和相角XX由以下公式给出： 相量Fj在复平面上是固定的，所以用滤波器的估计值简化了了电力系统相量之间的相关性。式子Cj和Sj分别是相量的虚部和实部，而相应的滤波器的正弦系数代表了相应的信号。所以，F’j可以由下式计算： 在这里i是复数算子，应注意到，只要Xj和Xj-N相等则F’j与F’j-1也相等。而如果输入的有偏离信号50Hz的偏离量（Δf）而取样始终保持不变，则根据文献[1]相量F’(50+Δf)与原始F’(50)关系为： [5]式可证明F’(50+Δf)的变化范围与Δf频域有限周期有关。另外，也可以观察到ΔF的量和相因子ΔΦ的变化。关系式为： 目前为止所得出的结果大部分应用于以ΔΦ直接测量ΔF。但是，没有人注意到其产生的有限周期量ΔF，因为它相对较小。这个量会在浮点微处理器的实时操作中被放大到使数字继电器降低可靠性。 3.递归DFT失稳的实验记录 递归DFT使用模拟[8][9]进行试验和使用，试验设置参数被应用于理想输入信号模式的产生，这个DFT测试与报告[10][11]中的测试方法类似。试验设置包括DSP板DS1003，该板带有I/O输入输出板DS2201[12][13]。之所以选择DS1003是因为它带有TMS320C40的DSP芯片，TMS320C40的机器周期为6M每秒，支持浮点运算，内置40位的累加器，32位寄存器和16位储存地址，其卓越的运算容量可概述为[14]-[17]。而一般通用的I/O板DS2201有20个分布式的模拟输入通道，5个12位A/D转换器，8个模拟输出端口，16个通用数字I/O口，5个PWM通道等等。A/D和D/A转换器使用10伏的参考电压，所有输入模拟信号都在DSP程序中被数字化并缩小比例至原始信号的十分之一，也就是说一个峰值为10伏的输入信号将内部表示为浮点数1.0。应注意到，浮点运算导致舍入误差的可能性会比整数运算的大，这是因为在整数运算的比例因数会小心地选择来避免舍入运算，而在浮点运算中却是不同的。浮点运算中这这种现象的细节和整数与浮点数之间的关系的说明将在4部分和5部分中说明。 (a)刚输入信号 (b)输入信号10分钟 (c)输入信号1小时 图1. Δf=0.5Hz量的放大误差实验数据 输入的正弦信号峰值为3.3伏，频率接近49.5Hz。而使用这种信号后，可用图1中的递推DFT记录到一个有限