1、微机继电保护原理微机继电保护原理 目 录第一章概述n第一节微机继电保护的发展历史及研究现状n第二节微机继电保护装置的特点和硬件组成第二章微机保护的数据采集系统n第一节 数据采集系统的作用与要求n第二节 模拟信号的预处理n第三节 采样定理n第四节 采用逐次逼近原理A/D芯片构成的数据采集系统n第五节 采用积分型A/D芯片构成的数据采集系统第三章微机保护装置的硬件电路原理n第一节 概述n第二节 管理微机系统的硬件电路原理n第三节 保护微机系统的硬件电路原理n第四节 微机保护装置的出口与信号电路 目录第四章 微机保护中数字滤波器的设计n第一节 概述n第二节 数字滤波器的实现和分析方法n第三节 简单滤
2、波器及其级联的分析方法n第四节 用零、极点配置法设计数字滤波器的方法第五章 微机保护的算法n第一节 微机保护的算法概述n第二节 基于正弦信号的算法n第三节 基于信号为周期函数模型的算法n第四节 微机保护中的阻抗算法n第五节 微机保护中的比相算法n第六节 微机保护中的滤序算法n第七节 微机保护中的复数求模值方法 目录第六章 输电线路的微机保护n第一节 输电线路的纵联保护n第二节 采用闭锁信号的高频保护n第三节 采允许锁信号的高频n第四节 高频方向保护中的方向元件n第五节 输电线路距离保护n第六节 故障类型判别与故障选相原理n第七节 微机保护中的工频变化量阻抗元件n第十节 微机保护装置中振荡与故障
3、的识别方法 目录 第七章 微机发电机变压器组保护 第二节 微机型发电机差动保护 第三节 发电机定子绕组单相接地保护 第四节 发电机的失磁保护 第五节 微机型变压器差动保护 第八章 微机保护的抗干扰设计 第一节 概述 第二节 干扰和干扰源 第三节 微机保护装置中的硬件抗干扰措施 第四节 微机保护装置中的软件抗干扰措施 第一章第一章 概述概述 第一节微机继电保护的发展历史及研究现状第一节微机继电保护的发展历史及研究现状英国剑桥大学的英国剑桥大学的P.G.Mclaran及其同事就提出用计算机构成电及其同事就提出用计算机构成电力系统继电保护的设想，并发表了力系统继电保护的设想，并发表了Sampling
4、 Techniques applied to derivation Letter 的文章。的文章。澳大利亚新南威尔士大学的澳大利亚新南威尔士大学的I.F.Morrison预测了输电线路计算机控制预测了输电线路计算机控制的前景。的前景。 美国西屋公司的美国西屋公司的G.D.Rockefeller发表年了发表年了利用数字计算机实现的利用数字计算机实现的故障保护故障保护的文章。的文章。1965年：年：1967年：年：1969年：年： 国际电子电气工程师学会教育委员会组织了一次世界性的计算国际电子电气工程师学会教育委员会组织了一次世界性的计算机继电保护研究班，对机继电保护研究班，对70年代以来的计算机
5、保护的研究成果进行年代以来的计算机保护的研究成果进行了总结和交流。到了总结和交流。到80年代中期计算机保护在电力系统中获得广泛年代中期计算机保护在电力系统中获得广泛应用。应用。 日本投入了一套以微处理机为硬件的控制与继电保护装置，全部代日本投入了一套以微处理机为硬件的控制与继电保护装置，全部代替了原有保护，大大减少了控制室的占地面积，并于替了原有保护，大大减少了控制室的占地面积，并于1980年发表年发表了试运行的结果。了试运行的结果。 1977年：年：1979年：年：1984年4月，华北电力大学研究的以MC6809CPU构成的MDP1型微机线路保护装置。这是我国研究成功的第一套微机线路保护装置
6、。并在河北某电厂投入试运行。我国微机保护的发展从硬件上看大体可分为三个阶段 以单CPU的8位微处理器构成的微机保护装置其主要特点为：1 采用8位微处理器MC6909CPU构成硬件系统。2 数据、地址、控制总线须引出插件外部。3 数据采集系统采用逐次逼近式A/D芯片构成。4 存储器的容量较小。5 保护的程序和定值都存在EPROM中，定值的修改十分不便。6 仅有软件时钟，当直流电源消失后，时钟停止运行，直流电源恢复后须重新校时。无接受GPS的接口。7 不具备数据远传功能。8 所有保护功能均由一个CPU处理，可靠性低。9 代表产品为WXB01WXB01型微机保护装置。第一阶段：第二个阶段：以多个8位
7、单片机组成的多微机系统。 其主要特点为：1 硬件为五个8位单片机组成的多微机系统。2 数据、地址、控制总线不引出插件。3 数据采集系统为VFC系统。4 保护装置的定值存在EEPROM中，定值修改十分方便。5 设有硬件时钟电路，装置直流电源消失后，依靠备用电池可使时钟继续运行。6 设计了与上位机通信的串行接口电路。7 具有液晶显示电路。调试方便。8 保护功能分散在各CPU中，可靠性高。9 代表产品为WXB11微机保护装置。其主要特点为：1硬件是以16位单片机构成的多微机系统。（80C196KB、M77）。2以日本三菱公司的M77芯片构成的微机系统，可以做到总线不引出芯片。3 保护装置的硬件设计除
8、了有硬件时钟外，装置还具备接受GPS全球定位系统的秒脉冲的接口。4 数据采集系统为VFC方式，最高转换频率为4MHZ。（VFC110）5 时钟芯片和EEPROM芯片均为串行方式。6具备较完善的通信网络，可应用于变电站综合自动化系统中。7 具有友好的人机接口界面。第三阶段以多个16位单片机组成的多微机系统。8 采用多层印刷电路板和表贴技术，增强了抗干扰能力。9 保护装置具备录波功能。其代表产品为CSL系列微机保护装置和LFP-900系列微机保护装置。二 微机保护的现状硬件方面：1采用32位DSP数字信号处理器或嵌入式系统构成的多微机系统。2 大屏幕液晶彩色显示，可显示汉字和图形。3 多层印刷电路
9、板，表贴技术、后插拔方式。软件方面：1 支持高级语言编程。人机界面采用WINGDONS系统2 自适应保护，智能体应用于保护中。3 人工神经网络在保护中的应用研究。4 模糊识别在继电保护中的应用研究。5 小波分析方法在继电保护中的应用研究。第二节微机继电保护装置的特点及硬件组成一微机保护装置的特点微机保护与常规保护相比具有以下特点：微机保护可以实现智能化。 微机保护可实现在线自检。 微机保护可 提供附加功能。 微机保护具有调试维护方便。 微机保护具有完善的网络通信功能 。微机保护便于采用一些新原理。微机保护的硬件可标准化。微机保护便于与数字传感器或光CT、光PT接口。二微机保护装置的硬件构成框图
10、DAS保护微机系统（1n个）开入电源 继电器开出开入管理微机系统GPS按键 LCD打印机逆变稳压电源 至通信网串口通信 二次电流电压DAS数据采集系统保护微机系统由单片机构成的微机系统。实现各种保护功能，根据保护对象的不同，保护功能可分布于16个插件中。管理微机系统由单片机构成的微机系统。主要功能是实现人机对话。例如，命令的输入，执行结果的输出显示或打印。开关量输入输入到微机保护装置中的开关量。在逻辑上仅有两种状态，1或0。一般是保护装置的投退压板，重合闸方式开关的接点位置等。开关量输出由微机保护的接口电路驱动继电器的部分。一般经光电隔离电路与继电器部分连接。按键一般为触摸按键，数量为516个
11、不等。GPS接受全球定位系统秒脉冲的接口。开入电源开关量输入部分的独立逆变稳压电源，24伏。电源装置的电源。为输入220V，输出为多个电压的逆变稳压电源。输出一般为5V、+12V（+15V）、-12V（-15V）、24V。第一节 数据采集系统的作用与要求第二章 微机保护的数据采集系统一 数据采集系统的作用将电流互感器二次、电压互感器二次的模拟信号转换为数字信号 。逐次逼近原理的A/D芯片构成的数据采集系统。采用VFC芯片构成的积分式数据采集系统。二 数据采集系统的设计要求 1模拟量。包括模拟量的数量，模拟量的性质，模拟量的动态范围。 2 数据采集系统的数据传送方式选择程序查询方式 中断方式 D
12、MA方式 转换完成？启动A/D转换读转换结果并存入指定地址单元调整地址指针指向下 一通道 通道号清零所有通道转换完成？程序查询方式的流程图 保护现场A/D中断请求状态复位读取A/D转换结果并存入循环存储区发第一通道转换命令地址更新通道号加1数据处理程序通道号清零发转换命令Y到最后通道？N采样中断程序 A/D转换完成中断程序 (2) 采样频率的选择与保护原理和采用的算法有关 （3）采样频率还与模拟量的数量，数据传送方式有关。 4 数据采集系统的精度对于A/D式数据采集系统，精度与A/D芯片的位数有关 。对VFC芯片构成的系统，精度与VFC的最高转换频率和保护中的算法有关。 3 采样频率的选择 （
13、1）采样频率的选择必须满足采样定理。 一 电压、电流变换器将CT、PT二次信号变换为适合A/D芯片量程的信号。 起隔离的作用。有利于防止干扰从二次回路侵入微机系统。在变换器的原、副边可加一屏蔽层。有利于抗干扰。二 模拟低通滤波器RR RC LR C CLR为了降低采样频率，满足采样定理，在采样之前先用一个模拟低 | )(|fH0/ ff通滤波器将频率高于采样频率一半的信号滤掉。 第二节 模拟信号的预处理 图23 A/D式数据采集系统的框图逐次逼近A/D芯片构成的典型数据采集系统， 变换器模拟低通采样/保持多路开关A/D变换器模拟低通采样/保持至微机系统 三 采样保持电路A1A2RAs输出A3模
14、拟输入逻辑输入逻辑参考Ch图26 采样保持电路原理图四 多个模拟通道的采样方式 1等间隔采样和不等间隔采样方式2跟踪采样和定位采样 3 多模拟通道的采样方式多通道同时采样方式 顺序采样方式 分组顺序采样方式 S/HS/HS/H多路转换开关A/D多路转换开关S/HA/DS/HS/HS/HS/H多路转换开关A/D S/HS/HS/HS/H多路转换开关A/D 第一组 第二组五 多路转换开关第四节采样定理 一 门函数的傅立叶变换)2(sin22sin)(2222cjeedtedtegFgFjjtjtjTTTg )(F 2 2 t 4 2 2 4 （a）门函数 (b)门函数的频谱1 )()2(sin2)
15、(nSSSTnncTtPF周期门函数阵列经过傅立叶变换 )(tPT )(TP t ST S 图211 门函数序列及其频谱二 冲激函数的傅立叶变换1)()()(dtetFtFtj)(t)(F图212 冲激函数及其频谱 t11将冲激函数变换到频域后它是一个连续频谱，对各频率成分的输出均为1。 对该冲激序列进行傅立叶变换可得：)()()(nSSTTntF )(tT )(Tt ST S 图213 冲激序列及其频谱三 矩形脉冲采样和冲激采样)(tf)(ts )(tfs)()()(tstftfs当用矩形脉冲对原始信号进行采样时 )()()()()(SnTTsnTtgtftPtftf当用冲激序列对原始信号进
16、行采样时 )()()()()(SnTsnTttfttftf 四 采样后信号的频谱分析)()()(tPtfFtfFTSS)()(TSPFF)()2(sin)()()2(sin)()2(sin*)(SnSSSnSSnSSSSnFncTnFncTnncFTF用矩形脉冲采样时，可对式（217）两边取傅立叶变换 可见，)(SF是)(F的反复，反复周期为S即将原始信号的频谱按s的间隔重复。其幅度被)2(sinSnc调制。当用冲激单位序列采样时，可对式（218）两边取傅立叶变换。)()()(ttfFtfFTSS)()(TSFF)(1)()(1)(*)(21nSSSnSnSSSnFTnFTnFF可见，经冲激序
17、列采样后的信号的频谱)(SF也是)(F反复周期为S即将原始信号的频谱按s的间隔重复。各频谱函数的幅度均相等。的反复。通过对以上两种采样后信号的频谱分析可知，连续信号经采样后，其频谱总是按采样频率S无限重复。当原始信号的最高频率为m时，如果采样频率mS2则采样后信号的频谱便不会发生混叠，否则将发生频谱混叠。 )(tf )(tT )(tfS t t t)(F )(T )(SF m m 第四节 采用逐次逼近原理的A/D芯片构成的数据采集系统 一 模数转换的一般原理nnBBBD2.222211RAUDUD/A转换器数码设定器控制器数字量输出 模拟输入比较器 图216逐次逼近A/D转换原理图四位A/D转
18、换的逼近过程示意图n第一次设定数码： 1 0 0 0n URUO URUO URU0 UR n第四次设定数码：1111 1101 1011 1001 0111 0101 0011 0001n 图217 四位A/D转换的逼近过程示意图数模转换的工作原理 USC-URabcRRR2R2R2R2R2RRF K3K2K1K0+ A1I 2R 2I3I4I)44332211443322112222(BBBBRUIBIBIBIBIRRUIR211221II 1234121III1348121IIIDRRURIUFRFSC 二模数转换芯片AD7874 S/HS/HS/H S/H控制逻辑内部时钟多路开关SAR
19、数据寄存器3V参考电压12位DAC INTRDCSCON 1inV CLK 2inV 3inV110 DD 4inVinREFoutREF 图219 AD7874内部逻辑电路图其主要特点如下 1 片内具有4路采样/保持放大器，可实现4通道同时采样；2 片内含有快速12位模数转换电路，每通道转换时间为8S3 输入电压范围为10V；4 5V供电电源；5 片上参考电压AD713AD 7874三利用AD7874芯片构成的数据采集系统 110 DD CON/CS/INT/-5VR1C1 C2-VS+VSR2 R3 +5V第五节采用积分型A/D芯片构成的数据采集系统一积分型AD转换原理 积分型AD又称为计
20、数式AD。设有一个计数器，输入计数器的脉冲信号的频率为f,计数的间隔为t，则在t时间内，加入计数器的脉冲数为N，则有：tfN 其一是用待转换的电压V控制时间间隔t，使t正比与输入电压V，而记数脉冲频率不变。在t时间间隔内，计得的脉冲数代表了输入电压的大小。 其二是用待转换的电压V控制计数脉冲的频率，使脉冲频率正比于输入电压，而记数的时间间隔不变。则记数结果代表了输入电压的大小。 积分型AD转换原理有电压时间型（V-t型），电压频率型（V-F型）两种方案。由上式可看出，在计数式模数转换系统中，可用两种方法实现。二 电荷平衡式V-F模数转换原理+VSt0t0TT t tA1A2单稳触发器VinCi
21、nt过零检测器RT CTRinKFout-VS+-IR Vint Vint设脉冲信号的周期为T，其中低脉冲的宽度为to,根据电荷平衡原理有如下关系： TRVtIiiR0因而：01tRIVTfiRiout 可见，输出脉冲信号的频率与输入电压的瞬时值成正比。三 两种电荷平衡式VFC芯片简介AD654OSCDRIVER +Vs CT CT -VsFout Dcom -Vin +Vin Fout 18+VsDcom 2-Vin3+ Vin 47CTCT-Vs561 AD654VFC芯片其主要参数如下：电源电压（Vs）：双端5VV18V；单端5V36V；500KHZ；模拟信号输入方式：负端电流输入方式，
22、 正或负端电压输入方式；模拟电压信号输入范围：单端供电方式时，0+Vs-4V；双端供电方式时，-Vs+Vs-4V；最大输入电流：1mA；输入阻抗：250M；输出方式及负载能力： 开路集电极输出，可驱动12个TTL负载；输出频率与输入电压关系：TioutRCVF10|满刻度输出频率：2 VFC110芯片Digtal CommonVFC110 Iin Vininput CommonAnalog Common VoutComparator in+Vs1413121110 9 8 FoutNC12 3+5Vref Out-Vs4 5 6 7EnableCos图225 VFC110管脚图One-Shot
23、Vref +15V5V112 1110 8 7 680RP5 FOUT50PF6 34 1325K214-15V -+-+ K 图226 满刻度为4MHZ的VFC110的接线图 四 由AD654VFC芯片构成的数据采集系统电路变换器浪涌吸收 器电压频率变 换 器光电隔离 器计数器iv, VFC计数器2i CT-5VRP1R1R2C1R3 RP2 R45V+Vs-Vs-+LBR6R5t Fout(KHZ)Vin(V)05-5-10500250VFC芯片输入电压与输出频率关系 五 由VFC110芯片构成的数据采集系统电路VFC110计数器光电隔离+15V5V5V地+15V-15V123456789
24、101112131412345678C21C51C11C31C41RW11RW21 六 VFC110芯片的电压与输出频率的关系f(MHZ)V10V4MHz 第三章 微机保护装置的硬件电路原理第一节 概述4U19英寸 变换器插件VFC插件CPU1插件CPU2插件CPU3插件CPU4插件人机接口插件开关量输入插件跳闸插件信号插件逻辑插件电源插件 第二节管理微机系统的硬件电路原理8031单片机按键电路液晶显示电路内部通信电路与上位机通信电路打印机接口电路硬件时钟电路EPROM SRAM EEPROM 8031单片机74LS37327256EPROM6264SRAM由8031单片机构成的8位单片机扩展
25、系统D0D7 A0-A15 P0口 ALE P2口 PSEN/ P3.6 P3.7 5V CS1 CS2/WE/ RD/ 二 按键电路 1 乒乓式按键电路并行口2 矩阵编码式按键电路 并 行口芯片0 1 2 012 5V 5V 三 并行打印机驱动电路并行口或非门或非门D1（3）D0（2）D2（4）D6（8）D7（31）打印机电源地（19-30）打印机5V（18）BUSY/RDY（11）选通线（1）打印地5V5V四 管理插件与保护插件的内部通信电路管理微机插件保护1微机插件保护2微机插件与门 与门5V5V第三节保护微机系统的硬件电路原理 8031单片机8253计数器电路定值写入电路内部通信电路开
26、关量输入电路开关量输出电路EPROM SRAM EEPROM8253可编程计数器芯片构成的计数器电路8253（3）8253（2）8253（1）Va IaIbVb Ic 3I0Vc(Vxl) 3U0A3/A2CS1/CS2/CS3/RD/WE/D0-D7 二定值存放电路7 4 LS1572817AEEPROM拨玛 开关 MEN 5V A0A6J45V 表32区号0123456789M0-M31111111011011100101110101001100001110110对应09区，四位拨码开关（M0、M1、M2、M3）的状态 24C16串行EEPROM构成的定值存放电路： 5VSCLSDA 三
27、开关量输出及开关量输入电路并行口1开关量输出电路 24V + 5V 7400 QDJ24V-CKJ1CKJ25VAXJ24V-74002 开关量输入电路并行口并行口 开入电源 24V外部接点开入电源 24V外部接点5V5V第四节 微机保护装置的出口和信号电路 一微机保护装置的跳闸出口电路J与门MC1413 1pin16pin 9pin V24V24 1 16 8 9电源地 2 3 415 14 13 5 6 7 12 11 10 QDJ2QDJ1QDJ3TXJCKJACKJBCKJCCKJRCKJQWXB-11WXB-11微机高压线路保护的跳闸出口电路微机高压线路保护的跳闸出口电路 V24 Q
28、DJ三取二闭锁GP、JL、LXGP、JL、LXGP、JL、LXGP、JL、LXGP、JL、LXGPJLLXGPCKJA U2 CKJB U2CKJC U2CKJR U2GJ U2GJ4 U2CKJA1 CKJB1 CKJC1CKJR1GJ1GJ41本地灯光信号中央信号+XMBDXJ1GI1 GJ41V24 V24V24V24CKJA2 CKJB2 CKJC2CKJR2 BDXJGJ12GJ22GJ32GJ U1FJ2 图317 信号插件电路图 第四章 微机保护中数字滤波器的设计第一节 概述有源滤波器无源滤波器R、L、C由R、L、C和运算放大器组成按频率特性划分低通滤波器高通滤波器带通滤波器带阻
29、滤波器模拟滤波器有源滤波器按实现方法划分非递归型滤波器按冲激响应划分有限冲激响应滤波器无限冲激响应滤波器数字滤波器递归型滤波器FFT型型数字滤波器的分类数字滤波器的优点 1 无物理器件，滤波器的性能稳定。 2 数字滤波器不受环境因素的影响。3 数字滤波器具有高度的灵活性。 4 数字滤波器不存在阻抗匹配的问题。5 数字滤波器可方便地做到分时复用。 S/HA/DT()滤波程序Y(n) X（t）X(n)数字滤波器的框图 第二节 数字滤波器的实现和分析方法一数字滤波器的冲激响应T（ .）)(t)(th设输入信号为冲激函数)(t，经过该系统后，其输出为)(th)(th称为滤波器的冲激响应。 当输入为冲激
30、信号时，滤波器的输出反映了滤波器的频率特性。所以冲激响应是滤波器特性的充分描述。一个任意时域函数)(tx，我们可以将其表示为无穷多个冲激函数之和。dtxtx)()()(当输入信号为)(tx时，该系统的输出为)(ty。则有： dtTxtxTty)()()()(对于一个时不变系统有： )()(11tthttT 所以有：dthxty)()()( 上式表明，只要我们知道了一个滤波器的冲激响应，就可以通过上式求出任意输入信号时滤波器的输出。该式是一个卷积分，可表示为：)(*)()(thtxty 上式也说明了在时域中滤波器的实现方法。即将输入信号与滤波器的冲激 响应进行卷积分即可实现。 在微机保护中，计算
31、机处理的是离散数字信号。因此数字滤波器应在离散域中实现。用)(ny)(nx)(nh分别表示与)(ty、)(tx 、)(th对应的离散数字信号：则：)(*)()(nhnxny二数字滤波器的实现方法在微机保护中实现数字滤波的方法。可归纳为以下三步：（1）确定数字滤波器的冲激响应)(th（2）按采样频率将)(th（3）用滤波系数与信号的采样值进行卷积，求出滤波器的输出 离散化，形成滤波系数。三数字滤波器的分析方法)().()(ZHZXZY)()()(ZXZYZHH(Z)称为系统的传递函数。 令传递函数中的Z为：STjeZ )(. | )(|)()(SSSTjTjTjeeHeHZH幅频特性为：| )(
32、|)(STjeHfF相频特性为：)()(STf 第三节简单滤波器及其级联的分析方法一差分滤波器差分滤波器的差分方程形式：)()()(knxnxny 式中：1k kZZXZXZY).()()(对上式作变换，得： 求得传递函数为 kZZXZYZH1)()()(令：STjeZ 幅频特性：|2sin|2)cos1 (2sin)cos1|sincos1|1| )(|22SSSSSSTjkTjTkTkTkTkTkjTkeeHSS相频特性：22)22()2(cos1sin)(111SSSSSSTkTktgtgTkctgtgTkTktgT| )(|fH 0m 0m 0m 0m 0ff 234 0 从而可求得滤
33、波器能滤除的谐波次数为： ).;2 , 1 , 0(.00kNmImIkNIm.2 , 02200sSSffmkTmfkfTk)( f0ff 2 1 二加法滤波器加法滤波器的差分方程形式：)()()(knxnxny1k式中： 对上式作变换，得：kZZXZXZY).()()(求得传递函数为：kZZXZYZH1)()()( 幅频特性：|2cos|2)cos1 (2sin)cos1|sincos1|1| )(|22SSSSSSTjkTjTkTkTkTkTkjTkeeHSS相频特性： SSSSSfkTTktgtgTkTktgT)2(cos1sin)(11滤波器能滤除的谐波次数为：)2.;2 , 1 ,
34、 0()21 (2).21 (.)21(00kNmImIkNIkNIm| )(|fH0ff0m0ff3 5 )( f 2 1 0 积分滤波器的差分方程为：三 积分滤波器)(.)2() 1()()(knxnxnxnxnyZ变换的形式为： kZZXZZXZZXZXZY).(.).().()()(21传递函数为： 1)1(2111.1)()()(ZZZZZZXZYZHkk滤波器的幅频特性： |2sin2) 1(sin| )(|SSTjTTkeHS 滤波器的相频特性为： SSfkTT)( 滤波器能滤除的谐波次数为： ).3 , 2 , 1 , 0(.) 1(0ImIIkNm| )(|fH0ff四加减滤
35、波器：)() 1(.)3()2() 1()()(knxnxnxnxnxnyk 该滤波器的特点是轮流进行加法与减法运算。其差分方程为：对应的变换形式为： kkZZXZZXZZXZZXZXZY).() 1(.).().().()()(321滤波器的传递函数为： 1)1(3211) 1(1) 1(.1)()()(ZZZZZZZXZYZHkkkkk为奇数和k为偶数时，滤波器的特性有所不同。 当k为奇数时 ：幅频特性： |2cos2) 1(sin| )(|SSTjTTkeHS相频特性： 22)(SSTkT.2 , 0) 1() 1(22) 1(0mNkfmfkfTksS当 时，滤波器的输出为零。 从而有
36、： .2 , 1 , 0(.10ImIIkNm其幅频特性和相频特性的形状与差分滤波器的形状类似。 当k为偶数时其幅频特性与相频特性的表达式： |2cos2) 1(cos| )(|SSTjTTkeHS幅频特性：SSfTkT)(相频特性：.2 , 1 , 0()21 () 1(2)21 ().21(10ImIkNIIkNm其幅频特性和相频特性的形状与加法滤波器的形状类似。 .25,23,2) 1() 1(22) 1(0mNkfmfkfTksS当 时，滤波器的输出为零 从而有： 五简单滤波器的级联)()()()()(121ZHzHZHZHZHmiim | )(| )(|1SSTjmiiTjeHeH)
37、()(1SmiiSTT设采样频率为600HZ。用差分滤波器、加法滤波器和积分滤波器进行级联。差分滤波器、加法滤波器、积分滤波器的传递函数分别为：)()()(321ZHZHZH、111)(ZZH221)(ZZH32131)(ZZZZH差分滤波器加法滤波器积分滤波器级联后滤波器的传递函数为： 642321211)1 ()1 ()1 ()(ZZZZZZZZZH滤波器对应的差分方程为： )6()4()2()()(nxnxnxnxny级联后滤波器能够完全滤除直流、3、6、9、12次谐波分量。 其幅频特性为： | )(|STjeH | )(|1STjeH | )(|2STjeH)(|3STjeH|2sin
38、|cos|4|2sin2sin|cos|2|2sin|2SSSSSSTTTTTT其相频特性为：SSSSSSSSTTTTTTTT322322)()()()(321如果我们把级联后的滤波器作为一个整体，分析它的幅频特性。可得 ：sSSSsSTjTjTjTjTjTTjTTjTeeeeHSSSS6sin6cos4sin4cos2sin2cos11)(642sssSSSTjTTTTTTeHS22)2sin6sin4(sin)6cos4cos2cos1 (| )(|滤波器的幅频特性 | )(|fH0ff第四节 用零、极点配置法设计数字滤波器的方法一 N阶线性时不变系统的差分方程和滤波器的零、极点N阶线性时
39、不变系统可用下面的常系数线性差分方程： ).()()()2() 1()() 1()()(102110inybinxaMnybnybnybNnxanxanxanyMiiNiiMN MN NNMMZaZaZaaZXZbZbZbZY 221102211)()1)(从而可得传递函数为： MiiiNiiiZbZaZH101)( 该传递函数的分子和分母可以表示为以下连乘积的形式： AZDZCAZHMiiNii1111)1 ()1 ()(MiiNiiDZCZ11)()(用sTjeZ代入得： MiiTjNiiTjTjDeCeAeHSSs11)()()()Im(Z)Re(Z Z平面上的单位圆iC称为滤波器的零点
40、； iD称为滤波器的极点； 二 用零、极点配置法设计非递归型数字滤波器在传递函数的公式（446）中，令：0ib，则： iNiiZaZH0)( 据此传递函数设计的数字滤波器称为非递归型数字滤波器。将该传递函数表示为连乘积的形式：)1 ()(11ZCAZHNiiiC为滤波器的零点，因该滤波器只有零点而没有极点，所以又称为全零点滤波器。 假定对m次谐波设置为零点，则滤波器的传递函数为： 2101010101011cos21)sincos1 ()sincos1)1 ()1 ()(00ZZTmZTmjZTmZTmjZTmZeZeZHsssssTjmTjmmss设采样频率为600HZ，则12,35NmsT
41、S 对直流分量、基波、二次谐波、直到11次谐波分别设置为零点，求出滤波器的传递函数如下。滤除直流、十二次谐波分量： 101)(ZZH滤除基波、十一次谐波分量： 21131)(ZZZH滤除二、十次谐波分量： 2121)(ZZZH滤除三、九次谐波分量： 231)(ZZH滤除四、八次谐波分量： 2141)(ZZZH滤除五、七次次谐波分量： 21531)(ZZZH滤除六次谐波分量： 161)(ZZH即有： )()(ZHZHmNm 设采样频率为600HZ。则：12,35NmsTS 用零、极点配置法设计两个数字滤波器。一个要求保留基波分量，而将直流和其他各整次谐波都滤除。另一各要求保留三次谐波，而将直流和
42、基波、二次、五次谐波滤除。1098764321218621212122113233231)31 ()1 ()1 ()31 ()1 ()1 ()1 ()1 ()(ZZZZZZZZZZZZZZZZZZZZZZZZH)10()9(3)8(2)7(3)6()4()3(3)2(2) 1(3)()(nxnxnxnxnxnxnxnxnxnxny保留基波分量的滤波器：滤波器的差分方程：22)sin32sin()2coscos31 (|2cos|2|2sin23sin|cos|2|2cos23cos|2sin|2| )(|sssssssssssTjTTTTTTTTTTTeHS滤波器的幅频特性：| )(|fH0f
43、f基波滤波器的幅频特性76543212121211221)31 ()1 ()31 ()1 ()(ZZZZZZZZZZZZZZZH保留三次谐波的数字滤波器： )7()6(2)5()4()3()2() 1(2)()(nxnnxnxnxnxnxnxny滤波器的差分方程： 2222)2sinsin3)2coscos31 ()sin32sin()2coscos31 (|2cos23cos|2sin|2| )(|ssssssssssTjTTTTTTTTTTTeHS滤波器的幅频特性： |)(|fH0ff滤波器对三次谐波的响应为：23)180sin90sin3()180cos90cos31 ()90sin31
44、80sin()180cos90cos31 (|45cos135cos|45sin2| )(|200200200200000900jeH三 用零、极点配置法设计递归型数字滤波器当式（446）的传递函数中，分母的系数0ib数字滤波器的当前最新输出不仅与历史上的输入有关，而且与历史上的输出有关。因此，这种滤波器成为递归型数字滤波器。递归型数字滤波器不仅 有零点、而且有极点。 时，构成的极点应设在单位圆内。 针对某次谐波设置一对极点为ssTjTjpepe、则极点对应的传递函数为：22111cos21)1 ()1 ()(ZpZTpZpeZpeZHSTjTjPSS 取23p，对基波设置为极点，采样频率为6
45、00HZ，6,35SSTmsT。则极点对应的传递函数为：2175. 05 . 11)(ZZZHp 在直流、二次、三次谐波处设置零点，对应的传递函数为： 54321221123321)1)(1)(1 ()(ZZZZZZZZZZHZ递归型滤波器的总传递函数为： 215432175. 05 . 1123321)(ZZZZZZZZH递归型数字滤波器所对应的差分方程为： )2(75. 0) 1(5 . 1)5()4(2)3(3)2(3) 1(2)()(nynynxnxnxnxnxnxny22)2sin75. 0sin5 . 1 ()2cos75. 0cos5 . 11 (|cos|2|2cos23cos
46、|2sin|2| )(|SSSSSSSSTjTTTTTTTTeHS该滤波器的幅频特性： 当866. 0p时，对基波分量，5 . 0p时，对基波分量，滤波器，p对极点频率的输出响应。 的大小直接影响到滤波器滤波器的输出响应为5.22。当的输出响应为1.52。可见第五章微机保护的算法第一节微机保护的算法概述算法是研究由若干个采样数据（已被量化）求取被测信号量值（包括有效值、相位、比值等）的方法。算法所研究的主要问题有两点。一是算法的计算精度。二是算法所用的数据窗 第二节基于正弦信号的算法一 半周内取最大绝对值算法二半周绝对值积分算法)(tu t 图52 半周绝对值积分算法示意图设积分的结果为S，则
47、：UtUtdtUdtatUSTmTmTm22|cossin| )sin(|202020 从而可求出有效值：22SU 在计算机中如何求积分呢？我们可以用求面积的方法求出S值。 )(tu S1 S2 S3 S4 S5 S6U0U1 U2U3U4U5U6STuuuuuuSSSSSSS)2|.2|2|(654321652110对于一个纯正弦信号，取绝对值后必有|60uu。 所以，上式可写为：STuuuS|)|.|(|510即有 SNKKTuS) |(120将S值代入式（54）中，得：KUTUUKNKNKSK*)(22)(12012022STK 式中： 三一阶导数算法(Mann_Morrison算法)
48、设：)cos()sin(cossin11111111tIitIitUutUumtmtmtmt2221212221212)(2)(IIiiUUuumttmtt则有： 从而可得： 21*)(2121ttuuU21*)(2121ttiiI 测量阻抗： 2 12122 1212ttttiiuuZ在计算机中如何求导数呢？通常是用差分来代替求导数。 )(ku ku1ku1k 1tkn 图54由采样值计算导数的示意图 211kktuuuSkktTuuu11t1时刻电压的导数为：1t时刻的电压1tu由平均值求得： 五 采样值积算法设：)sin(sin1111tIitUumtmt)(sin)sin()(sins
49、in122122TtItIiTtUtUummtmmt12ttT取11ttiu 、和 22ttiu 、两采样值的乘积：1 两采样值积算法)2cos(cos2).(sinsin.11111tIUttIUiummmmtt)22cos(cos2)(sin).(sin.11122TtIUTtTtIUiummmmtt取21ttiu 、和12ttiu 、两采样值乘积得：)2cos()cos(2)(sinsin.11121TtTIUTttIUiummmmtt)2cos()cos(2)sin()(sin.11112TtTIUtTtIUiummmmtt于是有： )2cos(cos2cos22.12211TtTIU
50、iuiummtttt)2cos(2coscos22.11221TtTIUiuiummttttTTiuiuiuiuIUttttttttmm212212211sincos)(cosTiuiuIUttttmmsinsin1221TTiiiiIttttm22122212sincos2如果都用电压或电流信号的采样值，则=0，此时可得 TTuuuuUttttm22122212sincos2 由于T是预先选定的值，所以sinT是常数。例如采样频率为600HZ，选定T=Ts，即t1和t2两个时刻仅相隔一个采样间隔，则sinT=0.5。当选择T=T/4，090T，则sinT=1，cosT=0，则有：222212
