第四章 微机继电保护算法课件

1、1第四章 微机保护算法一、概述一、概述被测信号采样值算法的核心问题：算法的核心问题：根据采样值计算出表征被保护根据采样值计算出表征被保护对象运行特点的物理量，如：电流、电压等信号对象运行特点的物理量，如：电流、电压等信号的有效值、相位、视在阻抗、序分量、基波分量、的有效值、相位、视在阻抗、序分量、基波分量、谐波分量等。谐波分量等。被测信号量值A/D转换算法1.1.保护算法的定义保护算法的定义 对采样数据进行分析、运算与判断，以实对采样数据进行分析、运算与判断，以实现各种继电保护功能的方法称为算法。现各种继电保护功能的方法称为算法。2RjXkZKZsetsetKZZ动作与边界条件：2.2.保护算

2、法的分类保护算法的分类由采样值求出保护所反映的量值由采样值求出保护所反映的量值以距离保护为例说明：以距离保护为例说明：电流I电压U测量阻抗Z=U/I与整定值比较保护算法保护算法采样值采样值优点：优点：容易实现常规保护无法实现的功能。如距离保护的特性可以很灵活。3270arg90270arg90kkKsetsetKsetkksetZZUI ZZZUI Z 或 -（a）测量阻抗在圆上 （b）测量阻抗在圆内 （c）测量阻抗在圆外 动作与边界条件：RjXoZsetRjXoZsetRjXoZsetZKZset+ZKZKZsetZKZK ZsetZset+ZKZKZset+ZKZK Zset（a）（b）（

3、c）根据继电器动作方程进行判断根据继电器动作方程进行判断 电流、电压相量电流I电压U阻抗继电器动作方程阻抗继电器动作特性采样值采样值保护算法保护算法特点：特点：不计算出具体的阻抗值。43.3.衡量算法的指标衡量算法的指标算法的速度算法的速度算法所要求的采样点数（数据窗）算法所要求的采样点数（数据窗）算法的运算量算法的运算量算法的精度算法的精度精度与速度之间的关系：精度与速度之间的关系：精度精度 数据窗长度数据窗长度增加，计算量增加，计算量算法的滤波性能算法的滤波性能 有些具有数字滤波的功能，有些则需要配以数有些具有数字滤波的功能，有些则需要配以数字滤波一起工作。字滤波一起工作。研究算法的实质研

4、究算法的实质：如何在：如何在速度速度和和精度精度两方面进行权衡两方面进行权衡 5二、假定输入为正弦函数的算法 )sin(2)(0tIti)sin(2)(0tUtu设i1、i2和u1、u2分别为两个相邻采样时刻tK和tK+1的采样值（tK+1=tK+T），则有： ）（1 )sin(2)(01KKtItii)sin(2)sin(2)(00112TtItItiiKKK）（2 )sin()cos(2)cos()sin(200TtITtIKK基于如下假设：基于如下假设：输入信号为纯正弦量，因此采用该类算法要获得比较理想的结果，必须与数字滤波器配合使用。 1.两点乘积算法设输入信号为：）（3 )sin(2

5、)(01KKtUtuu)sin(2)sin(2)(00112TtUtUtuuKKK）（4 )sin()cos(2)cos()sin(200TtUTtUKK求电流有效值求电流有效值I I由（1）和（2）式可得：）（5 )cos(2)sin()cos(012KtITTii6由（1）和（5）式可得：）（6 )(sin)cos(2222122212TTi iiiI求电压有效值求电压有效值U U方法与求电流有效值相同，可求得：）（7 )(sin)cos(2222122212TTuuuuU求阻抗（求阻抗（R R、X X）根据电流I和电压U求阻抗R、X的公式为：2sinsinIUIIUX2coscosIUI

6、IUR先求 和 ，将式（1）（4）两两相乘可得：sinUIcosUI1 1002sin()sin()KKu iUITT）（8 )22cos(cos0KTUI)22cos(cos)sin()sin(201010122KKKTUITTUIiu）（9 )222cos(cos0TTUIK)sin()sin(2)sin()sin(20000121KKKKTTTUITTUIiu）（10 )22cos()cos(0TTTUIK)sin()sin(2)sin()sin(20001012TTTUITTUIiuKKKK）（ 11 )22cos()cos(0TTTUIK7由式（10）和式（11）可求得：)222co

7、s()22cos(cos2002211TTTUIiuiuKK）（ 13 )22cos()cos(2cos20TTTUIK)22cos(2)cos()cos(01221TTTTUIiuiuK）（14 )22cos(2)cos(cos20TTTUIK由式（13）、（14）可求得：）（ 15 )(sin2)cos()(cos212212211TTiuiuiuiuUI由式（8）（11）可进一步求得：)cos()cos(1221TTUIiuiu)sin(sin2TUI）（12 )sin(2sin1221TiuiuUI即：由式（6）、（12）和（15）可求得：)cos(2)(sin2)sin(2sinsi

8、n212221212212TiiiiTTiuiuIUIIUX）（ 16 )cos(2)sin()(2122211221Ti iiiTiuiu8)cos(2)(sin2)(sin2)cos()(coscos21222122122122112TiiiiTTTiuiuiuiuIUIIUR）（17 )cos(2)cos()(21222112212211Ti iiiTiuiuiuiu特点特点数据窗仅为很短的一个采样间隔（两个采样点）；算式较复杂。当 时，公式可简化为：2T22221iiI22221uuU22211221iiiuiuX22212211iiiuiuR92.导数算法设输入信号为：)sin(2)

9、(0tIti)sin(2)(0tUtu设t1时刻电流、电压信号的瞬时值为： ）（ 1 )sin(2011tIi）（2 )sin(2011tUu求电流有效值I、电压有效值U对式（1）、（2）求导，可得： ）（3 )cos(2011tIi）（4 )cos(2011tUu由式（1）（4）可得：212122iiI212122uuU求阻抗（R、X）将电压电流表示成复数形式： 1112iIji1112uUju（5）（6）10由式（5）、（6）可得：1111111 1112212211111uu iiujuuiuiUZRjXjiiiIjiii 求电流、电压信号的导数基本思想：用差分近似求导。下面以电流信号为

10、例进行说明：如下图所示，电流信号在t1时刻的采样值i1和导数值i1可以用与t1时刻相邻的两个连续采样时刻tK和tK+1的采样值iK和iK+1近似计算，即：211KKiiiTiiiKK11特点： 数据窗仅为很短的一个采样间隔（两个采样点）； 求导数将放大高频分量； 差分近似求导数，要求有较高的采样频率。113.半周积分算法 基本思想：一个正弦信号在任意半周内，其绝对值积分（求面积）为常数S。 由上式可得：00)sin(2)sin(2dttIdttISItI22)cos(20SI42积分值S与积分起点的初相角无关 求面积S面积S可以采用梯形法近似求得：SNKKTiS120特点 数据窗长度为10ms

11、； 具有一定的滤除高频分量的能力； 不能抑制直流分量； 适用于要求不高的电流、电压保护中，可以采用差分滤波器滤除信号中的非周期分量。124.平均值、差分值的误差分析在实际应用中，常采用平均值平均值代替瞬时值瞬时值，用差分值差分值近似代替微分微分，用梯形法则近似求积分。当输入信号为纯正弦纯正弦信号时，用平均值平均值可以求出准确的瞬准确的瞬时值时值，用差分差分也可以求出准确的微分值准确的微分值。设信号为：)sin()(tXtxm设x(t)的两个采样值为x(n)和x(n+1),有：)2/(sin)(smTtXnx)2/(sin) 1(smTtXnx由平均求瞬时值由平均求瞬时值)2/(sin)2/(s

12、in212) 1()(smsmTtXTtXnxnx)2cos()()2cos()sin(ssmTtxTtX结论：结论：平均值x(n)+x(n+1)/2与瞬时值x(t)之间仅差一个系数 ，该系数与时刻t和初相角无关，仅与角频率和采样间隔Ts有关。)2cos(sT)1()(2) 1()()2/cos(1)(nxnxKnxnxTtxps对于单一纯正弦信号，由平对于单一纯正弦信号，由平均值求瞬时值的公式为：均值求瞬时值的公式为：13由差分值求微分值由差分值求微分值)sin()(tXtxm)cos()(tXdttdxm)2/(sin)2/(sin1)() 1(smsmssTtXTtXTTnxnxdttd

13、xTTTtXTTtXTsssmssms)()2sin(2)2sin()cos(2)2sin()cos(2结论：结论：差分值x(n+1)-x(n)/Ts与微分值dx(t)/dt之间仅差一个系数 该系数与时刻t和初相角无关，仅与角频率和采样间隔Ts有关。)2sin(2ssTT对于单一纯正弦信号，由差对于单一纯正弦信号，由差分值求微分值的公式为：分值求微分值的公式为：)() 1()() 1()2/sin(2)(nxnxKnxnxTdttdxcs当Ts足够小时，sin(Ts/2)越接近于Ts/2，Kc也越接近于1/Ts。14三、突变量电流算法 1.基本原理线路发生故障时，短路如图（a）所示。对于系统结

14、构不发生变化的线性系统，利用叠加原理可以得到（b）和（c）两个分解图。由叠加原理可得：)()()(tititikLm故障后的测量电流故障后的测量电流负荷电流负荷电流故障电流分量故障电流分量)()()(tititiLmk可求得故障电流分量为：对于正弦信号，在时间上间隔整周的两个瞬时值，其大小相等，即：)()(TtitiLLT:工频信号的周期工频信号的周期因此：)()()(TtititiLmk在非故障阶段测量电流等于负荷电流，即：)()(TtiTtimL)()()(Ttititimmk则故障分量电流为：短路前后的电流波形示意图短路前后的电流波形示意图15故障分量电流的采样值计算公式为：(1) Nk

15、kkiii当系统正常运行时，（1）式的输出为0；当系统刚发生故障的一周期内，（1）式输出的是纯故障分量；当负荷电流发生变化时，（1）式也有输出。因此（1）式反映的是电流的变化电流的变化，称为电流“突变量突变量”。当系统频率发生变化时，ik和ik-N对应电流波形的相位将有一个差值，当k在电流过零附近时，由于电流变化较快，不大的引起的不平衡电流较大，因此常采用下式求突变量电流。(2) 2NkNkNkkkiiiii说明：（2）式对应的突变量的存在时间不是20ms，而是40ms。如果由于频率偏移，造成ik和ik-N之间有一个相角差，则ik-N和ik-2N的相角差也应基本相同，（2）式右侧中的两项可以部

16、分抵消。因此采用（2）式可以补偿由于频率偏离产生的不平衡电流。162.频率变化的影响设一个工频周期的采样点数为N，分析电网实际频率偏离50Hz时对突变量计算公式（2）： 的影响。NkNkNkkkiiiii2以A相电流为例，设：)sin()(tItima）（3 )23cos()2cos(2sin2 )23cos(2sin)2cos(2sin2 )2(sin)(sin)(sin)sin( )2()()()()(TtTtTITtTTtTITtITtITtItITtiTtiTtititimmmmmmaaaaa（3）取最大值的条件是： 0)2cos(0)23cos(。，二者的结果是一致的或TtTt）（可

17、得：由12220)2cos(kTtTt）（）（因此：,.3 ,2, 1 ,0 12223kTkTt）（4 )sin()2sin(2) 12(2cos2sin2 )23cos(2sin2)(TTITkTITtTItimmma17Nkkkiii按公式（1）： 计算突变量受频率变化的影响为： ）（5 2sin2 )2cos(2sin2 )(sin)sin( )()()(TITtTITtItITtititimmmmaaa按公式（1）和（2）计算突变量的最大相对误差如下表所示。 f(Hz)最大相对误差484949.55050.55152（1）式的误差（%）25.0712.566.2806.2812.56

18、25.07（2）式的误差（%）6.231.580.3900.391.586.23结论结论：采用公式（2）计算突变量时，系统频率变化的影响要小得多。18四、傅里叶算法 000)cos()sin()(nnntnbtnatxN次谐波正弦项系数N次谐波余弦项系数基波角频率根据傅里叶级数和三角函数的正交性，可求出系数：）（ 1 )sin()(200TndttntxTa）（2 )cos()(200TndttntxTb因此，x(t)中的n次谐波分量可以表示为：1.基本原理 基本思想：假定被采样的模拟信号是一个周期性时间函数，可以通过傅里叶级数展开，表示为： )cos()sin()(00tnbtnatxnnn

19、同时，x(t)中的n次谐波分量又可以表示为： )sin(2)(0nnntnXtx)cos(sin2)sin(cos200tnXtnXnnnn比较xn(t)的两个表达式可得：nnnXacos2nnnXbsin2因此可以求n次谐波的幅值和相位：222nnnbaXnnnabtg1192.a2.an n和和b bn n的特点分析的特点分析nnTnXdttntxTacos2 )sin()(200 sin2 )cos()(200nnTnXdttntxTb从上式可以得出：采用傅氏算法求出的n次谐波分量xn(t)的正弦项系数an和bn是xn(t)的初始相角n的函数。也就是说， an和bn的值与积分开始时刻xn

20、(t)的相角有关。由于x(t)是周期函数，因此，可以得到计算an和bn的更一般的表达式为：TndttnttxTta0011)sin()(2）（TndttnttxTtb0011)cos()(2）（上式中若t1=0，即假定取从故障开始起的一个周期来积分，当t10时，x(t+ t1)将相对于时间坐标的零点向左平移，相当于积分从故障后t1开始。结论：结论：a an n超前超前b bn n9090改变t1不会改变n次谐波分量的有效值，但初始相角会改变。203.求阻抗（R、X） 设一个正弦信号为： 因此，正弦信号x(t)可表示为向量形式： )sin(2)(0tXtx可表示为向量形式：sincos.jXXX

21、eXj x(t)的正弦项系数和余弦项系数为：cos21Xa sin21Xb 2211.bjaXeXj将正弦电流、电压信号表示为向量形式：2211IIjbjaIeI2211.UUjbjaUeU212111111111.)()(IIIIUUIIUUbajbajbajbajbaIUjXRZ212111111111)()(IIIUIUIUIUbabaabjbbaa21211111)(IIIUIUbabbaaR21211111)(IIIUIUbabaabX电阻、电抗计电阻、电抗计算公式为：算公式为：214.离散傅里叶算法每工频周期采样N点，利用梯形法则可以求得： 122sin()NnKKaxnKNN12

22、2cos()NnKKbxnKNN特点：数据窗为一个工频周期，即20ms；运算量大，N次乘法和加法；抑制恒定直流分量和整数次谐波分量，对由于非周期分量引起的低频分量抑制能力较差。半波傅里叶算法的正弦项系数和余弦项系数的计算式为：2142sin()NnKKaxnKNN2142cos()NnKKbxnKNN特点：数据窗较短，为10ms；计算量较小，N/2次乘法和加法；不能滤除恒定直流分量和偶次谐波分量。225.5.傅里叶算法举例傅里叶算法举例电压：电流：全波傅里叶算法)3sin(*2100t)sin(*25t输入信号中无谐波分量输入信号中无谐波分量23)33sin(*2 . 0*2100)32sin

23、(*1 . 0*2100)3sin(*2100ttt)3sin(*2 . 0*25)2sin(*1 . 0*25)sin(*25ttt电压：电流：输入信号中有谐波输入信号中有谐波分量分量24电压：电流：半波傅里叶算法)3sin(*2100t)sin(*25t输入信号中无谐波输入信号中无谐波分量时分量时25)33sin(*2 . 0*2100)32sin(*1 . 0*2100)3sin(*2100ttt)3sin(*2 . 0*25)2sin(*1 . 0*25)sin(*25ttt电压：电流：输入信号中有谐波输入信号中有谐波分量（分量（2次和次和3次）时次）时26)33sin(*2 . 0*

24、2100)3sin(*2100tt)3sin(*2 . 0*25)sin(*25tt电压：电流：输入信号中有谐波分量（输入信号中有谐波分量（3次）时次）时276.6.衰减直流分量对傅氏算法的影响及补偿方法衰减直流分量对傅氏算法的影响及补偿方法设输入信号为：)()cos()sin()(000tiCetnbtnaCetxtnnnt=1/=1/， 为衰减时间常数为衰减时间常数设每周波采样N点，即采样间隔为Ts=T/N，则第m次采样值为：mmmmTmiCriCexs式中：sTer因此，对上式两端在一周期内积分，用矩形积分矩形积分近似，有：mNmkskmNmkskmNmkskkmNmkskTiTrCTi

25、CrTx1111)(考虑到交流分量在一个周期内的积分为0，即 ，因此有：01mNmkskTimNmkskmNmkskTrCTx11）（）（1 /1111NkkNmmNmkkmNmkkmNmkkrrxrxC数据窗右移一个采样间隔）（）（2 /112NkkNmmNmkkrrxCr由（2）/（1）可得：）（3 /112mNmkkmNmkkxxr因此根据采样值，利用（1）和（3）式，可以计算出衰减直流分量的初始值C和时间常数。28考虑n次谐波分量的正弦项系数an在第m个采样点处的计算式为:mmmiCrx)2sin()(2)2sin(2)(11NkNmkNmkNkNmknNnkiCrNNnkxNmaNk

26、NmkNkNmkNnkiNNnkCrN11)2sin(2)2sin(2NkNmkNkkNmNnkiNNnkrCrN11)2sin(2)2sin(2 /11）（NkkNmmNmkkrrxC)()2sin(21maNnkrCrNnNkkNm故障信号中故障信号中n n次谐波的正弦项系数次谐波的正弦项系数将 带入上式，可得：)()2sin()/(2)(111maNnkrrxNmanNkkNkkmNmkkn)/ )2sin(211NkkNkkArNnkrNK)()(1maxKmanmNmkkAn因此可求得：(4) )()(1mNmkkAnnxKmama同理可求得：(5) )()(1mNmkkBnnxKm

27、bmb)/ )2cos(211NkkNkkBrNnkrNK因此，当衰减时间常数已知时，修正系数KA和KB可以离线求得。297.递归式离散傅里叶算法根据离散傅里叶算法的计算式，考虑n次谐波分量的正弦项系数 在第 个采样点处的计算式为：nal）（1 )2sin(2)(1lNlKKnNnKxNla在第 个采样点的计算式为：1l）（2 )2sin(2) 1(12lNlKKnNnKxNla比较式（1）和（2）可得：NNlnxNlnxNlalaNllnn2)1(sin2) 1(sin2)() 1(11）（ 3 2) 1(sin)(2)(11NlnxxNlaNlln同理可以求得n次谐波分量的余弦项系数 的递

28、推表达式为： nbNNlnxNlnxNlblbNllnn2)1(cos2) 1(cos2)() 1(11）（4 2) 1(cos)(2)(11NlnxxNlbNlln特点：计算量小，每次只需要一次乘法和一次加减法运算；需要考虑累积误差对算法精度的影响。30）（6 )2sin(2)(1NkNmknNnkxNma考虑n次谐波分量的正弦项系数an和余弦项系数bn在第m个采样点处的计算式为:）（7 )2cos(2)(1NkNmknNnkxNmban和bn在第m+1个采样点处的计算式为:）（8 )2sin(2) 1(11NkNmknNnkxNma）（9 )2cos(2) 1(11NkNmknNnkxNm

29、b下面分析下面分析a an n(m)(m)、b bn n(m)(m)、a an n(m+1)(m+1)、b bn n(m+1)(m+1)之间的关系：之间的关系： )2sin(2)2sin(2)(21NkNmkNmnNnkxNNnxNma（6）式可变形为： 22) 1(sin22sin221NkNmkNmNnNknxNNnxN )2sin(2) 1(cos)2cos(2) 1(sin22sin221NkNmkNmNnNknNnNknxNNnxNNkNmkNkNmkNmNnNknxNNnNknxNNnxN221)2sin(2) 1(cos2)2cos(2) 1(sin22sin2式式式式11211

30、1222cos()sin (1)sin2sin22222cos()sinsin2(1)cos()Nk m NmmkNk mNmnknxn kxnxnNNNnxnkxna mnNNNN （）（）（）31)2cos()2cos(2) 1(cos)2sin(2112nxnxNknxNnNmmNkNmk1111)2sin(2)2sin() 1(2cos)2sin(2mnmNkNmkxNnNNnmbxNnkxNnN因此，可得：）（ 10 )2sin(2)2sin() 1()2cos() 1()2sin(2)(11mnnNmnxNnNNnmbNnmaNnxNma同理，根据（7）式可求得：）（ 11 )2c

31、os(2)2sin() 1()2cos() 1()2cos(2)(11mnnNmnxNnNNnmaNnmbNnxNmb根据（10）、（11）式可以求得：）（ 12 )2sin()()2cos()() 1(NnmbNnmamannn）（ 13 )-(2)2cos()()2sin()() 1(11NmmnnnxxNNnmbNnmamb32两种递归算法仿真分析）（ 1 2) 1(sin)(2)() 1(11NlnxxNlalaNllnn）（2 2) 1(cos)(2)() 1(11NlnxxNlblbNllnn）（3 )2sin()()2cos()() 1(NnmbNnmamannn）（4 )-(2

32、)2cos()()2sin()() 1(11NmmnnnxxNNnmbNnmamb电压：电流：)3sin(*2100t)sin(*25t输入信号：输入信号：递归算法二递归算法二递归算法一递归算法一递归算法二：递归算法一：33计算阻抗计算阻抗电压：电流：)3sin(*2100t)sin(*25t输入信号：输入信号：递归算法二递归算法二递归算法一递归算法一34衰减直流分量的影响衰减直流分量的影响输入信号：输入信号：)3sin(*210010002. 0tet算法二算法二算法一算法一时间常数已知的情况，r按下式计算：sTer35衰减直流分量的影响衰减直流分量的影响输入信号：输入信号：)3sin(*2

33、10010002. 0tet算法二算法二算法一算法一时间常数未知的情况，r按下式计算：mNmkkmNmkkxxr112/36五、解微分方程算法）（1 11dtdiLiRu1.基本原理 对于一般的输电线路，在短路情况下，线路分布电容产生的影响主要表现为高频分量，采用低通滤波器将高频分量滤除，就可以忽略线路分布电容的影响，因此，输电线路等效为R-L模型。在短路时，下列方程成立，即： 上式中：R1、L1分别为故障点至保护安装处线路段的正序电阻和电感； u、i为保护安装处的电压和电流。对于相间短路相间短路时，应采用u和i ，如AB相间短路时，取为uab和ia-ib对于单相接地短路单相接地短路时，取相电

34、压相电压和相电流加零序补偿电流相电流加零序补偿电流，以A相为例，（1）式可改写为：(2) )3()3(0101dtiKidLiKiRuxaraa（2）式中，Kr、Kx分别为电阻和电感分量的零序补偿系数，可用下式求出：1103rrrKr1103LLLKx其中：r0、r1、L0、L1分别为输电线路每公里的零序和正序电阻和电感。37dtdi11111DLiRu21212DLiRuD表示dtdi211221121DiDiDuDuR211212211DiDiiuiuL求得：2.短数据窗算法采用插值法可求得电流、电压信号在t1和t2时刻的值为： 221211KKKKiiiiii、221211KKKKuuu

35、uuu、采用差分近似求导求得： sKKsKKTiiDTiiD1211、)()()()(111111111KKKKKKKKKKKKKKKKiiiiiiiiiiuuiiuuRsKKKKKKKKKKKKKKKKTiiiiiiiiiiuuiiuuL)()()()(111111111以（1）式为例，两个不同采样时刻t1和t2分别测量u、i和 ，得到两个独立的方程： 3.长数据窗算法采用插值法可求得电流、电压信号在t1和t2时刻的值为： 采用差分近似求导求得： 2221211KKKKiiiiii、2221211KKKKuuuuuu、sKKsKKTiiDTiiD12211、384.算法的稳定性分析21122

36、1121DiDiDuDuR211212211DiDiiuiuL实质就是分析R1和L1的计算公式会不会出现 的情况。00当在出口附近短路时，分子将趋近于0，因此，如果分母出现两个非常接近的数相减，就会出现 的情况，从而导致算式的不稳定，出现很大的误差。00为便于分析便于分析，假设电流和电流的导数都是正弦正弦的，即：11sinmIi )sin(11DmDD)sin(12mIi)sin(12DmDD上式中：1为t1时刻电流的相角，D为电流的导数超前电流的角度，为t2滞后t1的角度。)sin(sin)sin()sin(11112112DDmmDIDiDi)sin(cos)cos(sinsin)sin(

37、)sincoscos(sin111111DDDmmDI)cos(sinsin)sin(sincos1111DDmmDI(1) sinsinDmmDI同理可求得：(3) )sin(sin2112DmmDUDuDu(2) sinsin1221mmUIiuiu电压超前电流的角度电压超前电流的角度对分母的分析对分母的分析从（1）式可以看出：分母的值与与t1时刻电流的相角1无关；在相间短路时，电流的导数总是超前于电流90，即D=90，带入（1）式可得：sin2112mmDIDiDi394.算法特点仅用于计算线路阻抗，应用于距离保护中； 不受电网频率变化的影响；不需要滤除非周期分量；具有分布电容的长线路，

38、将对算法产生误差；差分近似求导带来的误差。222121121iiiuiuR2221122111iiiuiuLX上式与两点乘积算法一样。因此，为了提高分母的数值，以便提高算法的稳定性，常采用长数据窗算法 )sin(sin2112DmmDUDuDusinsin1221mmUIiuiu因此，越接近90，分母的值越大。当=90时，D1=i2，D2=-i1，有：对电感计算公式的分析对电感计算公式的分析电感L的计算公式中的分子为：当金属性短路时， 90，因此上式同分母一样，其值与1无关。对电阻计算公式的分析对电阻计算公式的分析电阻R的计算公式中的分子为：当金属性短路时， 很小，可能出现两个相近的数相减。因此，电阻分量的计算相对误差一般要比电抗分量的误差大。 )sin(D40五、最小二乘方算