第4章--微机继电保护的算法《电力系统微机保护》ppt课件.ppt

《第4章--微机继电保护的算法《电力系统微机保护》ppt课件.ppt》由会员分享，可在线阅读，更多相关《第4章--微机继电保护的算法《电力系统微机保护》ppt课件.ppt（86页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、微机继电保护技术微机继电保护技术第四章第四章 微机微机继电保保护的算法及其仿真的算法及其仿真2023/1/20Xian University of Science and Technology主要内容主要内容4.1 4.1 微机保护算法概述微机保护算法概述4.2 4.2 基于正弦信号的算法基于正弦信号的算法4.3 4.3 基于周期信号模型的傅立叶算法基于周期信号模型的傅立叶算法4.4 4.4 基于非周期信号模型的基于非周期信号模型的R-LR-L模型算法模型算法4.5 4.5 基于随机数学模型的递推最小二乘算法基于随机数学模型的递推最小二乘算法4.6 4.6 移相算法及滤序算法移相算法及滤序算法

2、4.7 4.7 相位比较器算法相位比较器算法4.8 4.8 突变量电流算法突变量电流算法4.9 4.9 微机保护算法选择微机保护算法选择2/862023/1/20Xian University of Science and Technology4.1 4.1 微机保护算法概述微机保护算法概述4.1.1 算法的含义算法的含义对采集的数据进行对采集的数据进行处理、分析、判断处理、分析、判断以便实现保护功能的方法以便实现保护功能的方法称为算法。称为算法。Xian University of Science and Technology4.1.2 微机保护算法类型微机保护算法类型1.1.直接计算信号特

3、征值的算法直接计算信号特征值的算法该算法直接由采样值经过某种运算，求出实际值，再与定值比该算法直接由采样值经过某种运算，求出实际值，再与定值比较，实现保护。较，实现保护。2.2.判断继电动作方程的算法判断继电动作方程的算法该类算法依据继电器的动作方程，将采样值直接代入动作方程，该类算法依据继电器的动作方程，将采样值直接代入动作方程，转换为运算式的判断。转换为运算式的判断。4.1.2微机保护算法类型微机保护算法类型3/862023/1/20Xian University of Science and Technology4.1 4.1 微机保护算法概述微机保护算法概述4.1.3 衡量算法的指标衡

4、量算法的指标1.1.算法的精度算法的精度指的是时计算出的结果信号与实际值的指的是时计算出的结果信号与实际值的逼近程度逼近程度。2.2.算法的速度算法的速度包括两个方面：包括两个方面：1 1、算法所要求的、算法所要求的采样点数采样点数；2 2、算法的、算法的运算工作量运算工作量。3.3.算法的精度与速度的关系算法的精度与速度的关系精度和速度总是精度和速度总是矛盾的矛盾的。若要计算精确，则往往要利用更多的。若要计算精确，则往往要利用更多的采样点和进行更多的计算工作量。采样点和进行更多的计算工作量。Xian University of Science and Technology4.1.3衡量算法的

5、指标衡量算法的指标4/862023/1/20Xian University of Science and Technology4.2 4.2 基于正弦信号的算法基于正弦信号的算法Xian University of Science and Technology4.2.1 4.2.1 半周内取最大绝对值算法半周内取最大绝对值算法1.1.算法原理算法原理要计算有效值，最直接的方法是在任意半周期内取各采样值的要计算有效值，最直接的方法是在任意半周期内取各采样值的绝对值，然后寻求最大值，将最大值除以绝对值，然后寻求最大值，将最大值除以 即得到有效值。即得到有效值。如图如图4-14-1所示，有：所示，有

6、：(4-2)其中，其中，为半周波内的采样值。为半周波内的采样值。4.2.1绝对值算法绝对值算法5/862023/1/20Xian University of Science and Technology4.2 4.2 基于正弦信号的算法基于正弦信号的算法Xian University of Science and Technology2.2.算法指标算法指标数据窗为数据窗为半个周波。半个周波。当采样时刻刚好落在最大值，则此算法本身当采样时刻刚好落在最大值，则此算法本身没有误差没有误差。但当采样时刻落在任意其它时刻，这种算法本身但当采样时刻落在任意其它时刻，这种算法本身会产生误差会产生误差。最大

7、绝对误差：最大绝对误差：(4-3)最大相对误差：最大相对误差：(4-4)4.2.1绝对值算法绝对值算法6/862023/1/20Xian University of Science and Technology4.2 4.2 基于正弦信号的算法基于正弦信号的算法Xian University of Science and Technology4.2.2 4.2.2 两点乘积算法两点乘积算法1.1.两点乘积算法原理两点乘积算法原理利用两个间隔一定的采样值计算利用两个间隔一定的采样值计算有效值及相位有效值及相位的方法。的方法。以电流为例，设以电流为例，设 和和 分别为两个电气角度相隔为分别为两个电

8、气角度相隔为 的采样的采样值（如图值（如图4-84-8所示），即：所示），即：(4-5)根据式（根据式（4-14-1）有：）有：(4-6)(4-7)4.2.2两点乘积算法两点乘积算法7/862023/1/20Xian University of Science and Technology4.2 4.2 基于正弦信号的算法基于正弦信号的算法Xian University of Science and Technology可以利用电压和电流的有效值测量阻抗，有两种方法：可以利用电压和电流的有效值测量阻抗，有两种方法：方法一方法一：利用公式可以得出测量阻抗的：利用公式可以得出测量阻抗的模值模值Z

9、Z 和和阻抗角阻抗角 。(4-9)同理电压也有上述关系：同理电压也有上述关系：(4-10)(4-11)(4-8)则有：则有：(4-13)(4-12)可得：可得：(4-14)(4-15)4.2.2两点乘积算法两点乘积算法8/862023/1/20Xian University of Science and Technology4.2 4.2 基于正弦信号的算法基于正弦信号的算法Xian University of Science and Technology方法二方法二：电流和电压都是正弦量，写成复数形式：电流和电压都是正弦量，写成复数形式参照式（参照式（4-64-6）、（）、（4-74-7）、

10、（）、（4-104-10）、（）、（4-114-11）有：）有：则有：则有：(4-16)式（式（4-164-16）的实部就是电阻分量）的实部就是电阻分量R R，虚部则为电抗分量，虚部则为电抗分量X X：(4-17)(4-18)4.2.2两点乘积算法两点乘积算法9/862023/1/20Xian University of Science and Technology4.2 4.2 基于正弦信号的算法基于正弦信号的算法Xian University of Science and Technology2 2 两点乘积算法特点两点乘积算法特点（1 1）算法本身所需的数据窗长度为工频的）算法本身所需的

11、数据窗长度为工频的1/41/4个周期个周期。（2 2）此算法基于正弦信号，因此要与带通滤波器）此算法基于正弦信号，因此要与带通滤波器配合使用配合使用。（3 3）算法本身与）算法本身与采样频率无关采样频率无关，因此对采样频率无特殊要求。，因此对采样频率无特殊要求。（4 4）算法本身）算法本身无误差无误差。（5 5）算法中要进行较多的乘除法，运算）算法中要进行较多的乘除法，运算工作量大工作量大。最后必须指出，两点乘积算法中的两点不必一定相隔最后必须指出，两点乘积算法中的两点不必一定相隔1/41/4个周个周期期，其数据窗长度可以取任意值，只是取任意值时的数学模型略微，其数据窗长度可以取任意值，只是取

12、任意值时的数学模型略微复杂一点而已。复杂一点而已。4.2.2两点乘积算法两点乘积算法10/864.2 4.2 基于正弦信号的算法基于正弦信号的算法4.2.3 4.2.3 三点乘积算法三点乘积算法1.1.三点乘积算法原理三点乘积算法原理利用连续利用连续三个采样点三个采样点的数据计算正弦信号特征量的方法的数据计算正弦信号特征量的方法对正弦电压和电流任意取三个连续时刻对正弦电压和电流任意取三个连续时刻 、对应的采样值分别为对应的采样值分别为 、和和 、。电压信号：电压信号：(4-21)(4-23)因为：因为：所以：所以：(4-22)同理有：同理有：4.2.3三点乘积算法三点乘积算法2023/1/20

13、Xian University of Science and TechnologyXian University of Science and Technology11/862023/1/20Xian University of Science and Technology4.2 4.2 基于正弦信号的算法基于正弦信号的算法Xian University of Science and Technology2.2.三点乘积算法特点三点乘积算法特点（1 1）算法）算法数据窗数据窗为为2 2T Ts s，因此，算法的延时与采样频率有关。，因此，算法的延时与采样频率有关。（2 2）算法本身）算法本身无

14、误差无误差。（3 3）作为基于正弦信号的算法，必须与滤波器）作为基于正弦信号的算法，必须与滤波器配合使用配合使用。（4 4）算法中要进行较多的乘除法，运算）算法中要进行较多的乘除法，运算工作量大工作量大。根据电压和电流得到根据电压和电流得到阻抗模值阻抗模值：(4-24)电阻分量：电阻分量：(4-26)(4-25)电抗分量：电抗分量：4.2.3三点乘积算法三点乘积算法12/862023/1/20Xian University of Science and Technology4.2 4.2 基于正弦信号的算法基于正弦信号的算法Xian University of Science and Tech

15、nology4.2.4 4.2.4 导数算法导数算法1.导数算法原理导数算法原理导数算法只需知道输入正弦量在某一时刻导数算法只需知道输入正弦量在某一时刻 的采样值及该时的采样值及该时刻对应的导数，即可算出刻对应的导数，即可算出有效值有效值和和相位相位。假设电压和电流为假设电压和电流为:(4-29)则电压和电流的导数为则电压和电流的导数为：(4-30)(4-27)(4-28)4.2.4导数算法导数算法13/862023/1/20Xian University of Science and Technology4.2 4.2 基于正弦信号的算法基于正弦信号的算法Xian University of

16、 Science and Technology在在 时刻有：时刻有：(4-34)(4-33)(4-31)(4-32)则可以利用正弦函数和余弦函数的平方和等于则可以利用正弦函数和余弦函数的平方和等于1 1得：得：(4-35)(4-36)4.2.4导数算法导数算法14/862023/1/20Xian University of Science and Technology4.2 4.2 基于正弦信号的算法基于正弦信号的算法Xian University of Science and Technology(4-37)于是被采样的电压和电流对应的阻抗为：于是被采样的电压和电流对应的阻抗为：(4-38)

17、以上分析表明，只要知道电压电流某点采样值和在该点的导数，以上分析表明，只要知道电压电流某点采样值和在该点的导数，就可以就可以求出电压电流的大小和相位，求出测量阻抗。就可以就可以求出电压电流的大小和相位，求出测量阻抗。为求导数，可取为求导数，可取 为两个相邻采样时刻为两个相邻采样时刻 和和 的中间的中间时刻（如图时刻（如图4-214-21所示）。然后用差分近似求导，则有：所示）。然后用差分近似求导，则有：(4-39)(4-40)4.2.4导数算法导数算法15/862023/1/20Xian University of Science and Technology4.2 4.2 基于正弦信号的算法

18、基于正弦信号的算法Xian University of Science and Technology电流、电压瞬时值用平均值替代，则有：电流、电压瞬时值用平均值替代，则有：(4-41)(4-45)(4-46)(4-42)另外求导数也可以用相邻三个采样点另外求导数也可以用相邻三个采样点 、，取取 ，则有：，则有：(4-43)(4-44)4.2.4导数算法导数算法16/862023/1/20Xian University of Science and Technology4.2 4.2 基于正弦信号的算法基于正弦信号的算法Xian University of Science and Technol

19、ogy2.2.导数算法特点导数算法特点（1 1）数据窗短，仅为两个或三个采样点，即时窗为一个或两）数据窗短，仅为两个或三个采样点，即时窗为一个或两个采样间隔，因此算法的时延与采样频率有关。个采样间隔，因此算法的时延与采样频率有关。（2 2）运算工作量与乘积算法相似。）运算工作量与乘积算法相似。（3 3）导数算法存在方法误差。）导数算法存在方法误差。（4 4）算法具有一定的抑制直流分量的能力，但对高次谐波具）算法具有一定的抑制直流分量的能力，但对高次谐波具有放大的作用，因此要求滤波器应有良好的滤除高频分量的能力有放大的作用，因此要求滤波器应有良好的滤除高频分量的能力。4.2.4导数算法导数算法1

20、7/862023/1/20Xian University of Science and Technology4.2 4.2 基于正弦信号的算法基于正弦信号的算法Xian University of Science and Technology4.2.5 半周绝对值积分算法半周绝对值积分算法1.1.半周绝对值积分算法原理半周绝对值积分算法原理算法的依据是一个正弦量在任意半个周期内绝对值的积分为一算法的依据是一个正弦量在任意半个周期内绝对值的积分为一常数常数S S，如图，如图4-284-28所示：所示：(4-48)积分面积积分面积S S为：为：(4-47)若已知积分面积若已知积分面积S S，则有：

21、，则有：4.2.5半周绝对值积分算法半周绝对值积分算法18/862023/1/20Xian University of Science and Technology4.2 4.2 基于正弦信号的算法基于正弦信号的算法Xian University of Science and Technology计算积分时可以用梯形法或矩形法。计算积分时可以用梯形法或矩形法。(4-49)如图如图4-294-29所示，设若干个小梯形面积之和为所示，设若干个小梯形面积之和为 ，则有：，则有：式中：式中：为为 时的采样值，时的采样值，N N为采样点数。为采样点数。梯形法求积分梯形法求积分：4.2.5半周绝对值积分算

22、法半周绝对值积分算法19/862023/1/20Xian University of Science and Technology4.2 4.2 基于正弦信号的算法基于正弦信号的算法Xian University of Science and Technology矩形法求积分：矩形法求积分：矩形法比梯形法公式较简洁，便于编程，但在相同的矩形法比梯形法公式较简洁，便于编程，但在相同的 下，下，精度较梯形法差。精度较梯形法差。(4-50)2.2.半周绝对值积分算法特点半周绝对值积分算法特点（1 1）数据窗长度为）数据窗长度为半个周期半个周期。（2 2）由于进行的是积分运算，故）由于进行的是积分运算

23、，故具有滤波功能具有滤波功能，对高频分量，对高频分量有拟制作用，但不能拟制直流分量有拟制作用，但不能拟制直流分量（3 3）本算法的精度与）本算法的精度与采样频率有关采样频率有关。（4 4）由于只有加法运算，计算）由于只有加法运算，计算工作量小工作量小。4.2.5半周绝对值积分算法半周绝对值积分算法20/862023/1/20Xian University of Science and Technology4.2 4.2 基于正弦信号的算法基于正弦信号的算法Xian University of Science and Technology4.2.6 4.2.6 频率算法频率算法指的是用采样值计算

24、交流信号频率的方法指的是用采样值计算交流信号频率的方法。(4-55)设采样频率为设采样频率为600Hz600Hz，则电压的采样值为：，则电压的采样值为：(4-54)则有：则有：(4-51)(4-52)(4-53)从而可求出：从而可求出：4.2.6频率算法频率算法21/862023/1/20Xian University of Science and Technology4.2 4.2 基于正弦信号的算法基于正弦信号的算法Xian University of Science and Technology(4-56)为使计算更为精确，采样信号最好用线电压。例如，对于一个为使计算更为精确，采样信号最

25、好用线电压。例如，对于一个频率为频率为45Hz45Hz的交流信号按的交流信号按 进行采样，取进行采样，取 ，得得 当：当：时，时，f=50Hz 当：当：时，时，f=050Hz之间之间 当：当：时，时，f=50100Hz之间之间 4.2.6频率算法频率算法22/862023/1/20Xian University of Science and Technology4.2 4.2 基于正弦信号的算法基于正弦信号的算法Xian University of Science and Technology4.2.7 4.2.7 算法中的平均值、差分值的误差分析算法中的平均值、差分值的误差分析从算法的原理中

26、我们知道，在实际应用的数值计算中，经常会从算法的原理中我们知道，在实际应用的数值计算中，经常会遇到通过采样值求瞬时值、微分值和积分值的。遇到通过采样值求瞬时值、微分值和积分值的。常用方法是：用常用方法是：用平均值代替瞬时值平均值代替瞬时值，用，用差分近似代替微分差分近似代替微分，用，用梯形法或矩形法则近似求积分梯形法或矩形法则近似求积分。设信号为设信号为:(4-57)4.2.7 算法误差分析算法误差分析经采样后，就可以得到两个采样值经采样后，就可以得到两个采样值 和和 的两个采样值的两个采样值（其中时刻位于（其中时刻位于 和和 采样时刻的中间，如图采样时刻的中间，如图4-364-36所示所示

27、。23/862023/1/20Xian University of Science and Technology4.2 4.2 基于正弦信号的算法基于正弦信号的算法Xian University of Science and Technology表达式：表达式：(4-58)(4-59)1 1 由平均值求瞬时值的误差分析由平均值求瞬时值的误差分析(4-60)显然，显然，平均值平均值与与瞬时值瞬时值之间也仅相差一个系数之间也仅相差一个系数 ，且这个系数只与角频率和采样时间有关。可以进行无误差地修正。且这个系数只与角频率和采样时间有关。可以进行无误差地修正。4.2.7 算法误差分析算法误差分析24/

28、862023/1/20Xian University of Science and Technology4.2 4.2 基于正弦信号的算法基于正弦信号的算法Xian University of Science and Technology通过式（通过式（4-604-60）可知，对于单一的纯正弦信号，可以由平均值）可知，对于单一的纯正弦信号，可以由平均值求出准确的瞬时值。具体计算公式为：求出准确的瞬时值。具体计算公式为：(4-61)式中，式中，为常数。为常数。4.2.7 算法误差分析算法误差分析25/862023/1/20Xian University of Science and Techno

29、logy4.2 4.2 基于正弦信号的算法基于正弦信号的算法Xian University of Science and Technology2 2 由差分值求微分值的误差分析由差分值求微分值的误差分析设信号表达式为：设信号表达式为：(4-62)于是：于是：显然，差分值与微分值之间也仅相差一个系数显然，差分值与微分值之间也仅相差一个系数 ，且这个系数只与角频率和采样时间有关。可以进行无误差地修正。且这个系数只与角频率和采样时间有关。可以进行无误差地修正。4.2.7 算法误差分析算法误差分析26/862023/1/20Xian University of Science and Technolo

30、gy4.2 4.2 基于正弦信号的算法基于正弦信号的算法Xian University of Science and Technology由此可以看出，由此可以看出，和和 二者之间仅相差一个系数二者之间仅相差一个系数 。这是一个与时刻。这是一个与时刻t t和初相角和初相角 无关的参数，仅与角频率无关的参数，仅与角频率 和采样和采样间隔间隔 有关。有关。因此，通过式因此，通过式（4-60）可知，对于单一的纯正弦信号，可以由可知，对于单一的纯正弦信号，可以由平均值求出准确的瞬时值。具体计算公式为平均值求出准确的瞬时值。具体计算公式为：(4-63)式中，式中，为常数。为常数。4.2.7 算法误差分析

31、算法误差分析27/862023/1/20Xian University of Science and Technology4.3 4.3 基于周期信号的算法基于周期信号的算法Xian University of Science and Technology4.3.1 4.3.1 周期函数的傅立叶级数与各次谐波关系周期函数的傅立叶级数与各次谐波关系傅氏算法假定被采样信号是一个周期性时间函数，除傅氏算法假定被采样信号是一个周期性时间函数，除基波基波外还外还含有不衰减的含有不衰减的直流分量直流分量和和各次谐波各次谐波。(4-65)设该周期信号为：设该周期信号为：(4-64)m m次谐波为：次谐波为：

32、其中其中 为直流分量，为直流分量，为为m m次谐波，次谐波，m m取自然数。取自然数。其中：其中：(4-67)(4-66)和和 分别为正弦分量和余弦分量系数。分别为正弦分量和余弦分量系数。4.3.1傅立叶级数与各次谐波关系傅立叶级数与各次谐波关系28/862023/1/20Xian University of Science and Technology4.3 4.3 基于周期信号的算法基于周期信号的算法Xian University of Science and Technology(4-69)(4-68)有了有了 和和 ，就可以计算出，就可以计算出m m次谐波的次谐波的幅值幅值和和相位相位

33、：m次谐波分量的复数形式为次谐波分量的复数形式为：(4-70)(4-71)将式将式（4-65）代入式代入式（4-62）中，有：中，有：4.3.1傅立叶级数与各次谐波关系傅立叶级数与各次谐波关系29/862023/1/20Xian University of Science and Technology4.3 4.3 基于周期信号的算法基于周期信号的算法Xian University of Science and Technology4.3.2 4.3.2 全波傅立叶算法全波傅立叶算法1.全波傅里叶算法的原理全波傅里叶算法的原理当已知周期函数当已知周期函数 时，其时，其m m次谐波的正、余弦系数

34、为：次谐波的正、余弦系数为：(4-74)(4-72)基波分量的正、余弦分量的系数为：基波分量的正、余弦分量的系数为：(4-75)(4-73)4.3.2全波全波傅立叶算法傅立叶算法30/862023/1/20Xian University of Science and Technology4.3 4.3 基于周期信号的算法基于周期信号的算法Xian University of Science and Technology矩形法计算正弦和余弦分量的系数矩形法计算正弦和余弦分量的系数.(4-77)(4-76)基波分量基波分量的正、余弦分量的系数为：的正、余弦分量的系数为：则，则，m m次次谐波分量谐

35、波分量的正、余弦分量的系数为：的正、余弦分量的系数为：其中，其中，为采样间隔，为采样间隔，为一个采样时间对应的相角为一个采样时间对应的相角。4.3.2全波全波傅立叶算法傅立叶算法31/862023/1/20Xian University of Science and Technology4.3 4.3 基于周期信号的算法基于周期信号的算法Xian University of Science and Technology式（式（4-764-76）和式（）和式（4-774-77）可写为另一种形式：）可写为另一种形式：(4-79)(4-78)则基波分量的则基波分量的复数形式复数形式为：为：(4-81

36、)(4-80)4.3.2全波全波傅立叶算法傅立叶算法32/862023/1/20Xian University of Science and Technology4.3 4.3 基于周期信号的算法基于周期信号的算法Xian University of Science and Technology计算出各基波电流、电压的正弦和余弦系数计算出各基波电流、电压的正弦和余弦系数 、后，后，可以计算可以计算阻抗阻抗。(4-84)(4-82)(4-83)4.3.2全波全波傅立叶算法傅立叶算法33/862023/1/20Xian University of Science and Technology4.3

37、 4.3 基于周期信号的算法基于周期信号的算法Xian University of Science and Technology2.2.全波傅立叶算法的滤波特性全波傅立叶算法的滤波特性傅氏方法也是一种傅氏方法也是一种滤波方法滤波方法。根据式（根据式（4-78）和（）和（4-79）基波分量在第）基波分量在第n点傅里叶正弦和余弦点傅里叶正弦和余弦项系数可表示为：项系数可表示为：(4-87)Z Z变换：变换：(4-88)4.3.2全波全波傅立叶算法傅立叶算法34/862023/1/20Xian University of Science and Technology4.3 4.3 基于周期信号的算法

38、基于周期信号的算法Xian University of Science and Technology因此，傅立叶正、余弦系数的因此，傅立叶正、余弦系数的单位冲激响应序列单位冲激响应序列为：为：(4-91)(4-90)正、余弦系数正、余弦系数频率特性频率特性如图如图4-414-41和图和图4-424-42所示。所示。图图4-414-41正弦项系数频率特性正弦项系数频率特性 图图4-424-42余弦项系数频率特性余弦项系数频率特性4.3.2全波全波傅立叶算法傅立叶算法35/862023/1/20Xian University of Science and Technology4.3 4.3 基于周

39、期信号的算法基于周期信号的算法Xian University of Science and Technology令电压输入信号为：令电压输入信号为：(4-91)(4-90)则第则第k k个采样值为：个采样值为：相对频率相对频率 的的幅频特性幅频特性定义为：定义为：不同初相角时傅立叶算法的幅频特性，如图不同初相角时傅立叶算法的幅频特性，如图4-434-43所示：所示：4.3.2全波全波傅立叶算法傅立叶算法36/862023/1/20Xian University of Science and Technology4.3 4.3 基于周期信号的算法基于周期信号的算法Xian University

40、of Science and Technology可见，此算法可以完全滤除直流分量，但不能滤除非整次谐波分量。可见，此算法可以完全滤除直流分量，但不能滤除非整次谐波分量。不同初相角时，对于基波、直流分量和各整次谐波分量的频率响应相同。不同初相角时，对于基波、直流分量和各整次谐波分量的频率响应相同。图图4-434-43全波傅立叶算法频率特性全波傅立叶算法频率特性4.3.2全波全波傅立叶算法傅立叶算法37/8638/862023/1/20Xian University of Science and Technology4.3 4.3 基于周期信号的算法基于周期信号的算法Xian Universit

41、y of Science and Technology4.3.3 4.3.3 半波傅立叶算法半波傅立叶算法1.1.半波傅立叶算法原理半波傅立叶算法原理仅用半周波的数据计算信号的仅用半周波的数据计算信号的幅值幅值和和相角相角。针对基波分量，计。针对基波分量，计算正弦和余弦系数的矩形方法如下：算正弦和余弦系数的矩形方法如下：(4-87)(4-88)4.3.3半波半波傅立叶算法傅立叶算法2023/1/20Xian University of Science and Technology4.3 4.3 基于周期信号的算法基于周期信号的算法Xian University of Science and T

42、echnology当当 N=12N=12时，上式为：时，上式为：(4-87)(4-88)同样，求出同样，求出 和和 后，就可以通过式（后，就可以通过式（4-824-82）式（）式（4-834-83）和）和式（式（4-844-84）求出信号的）求出信号的有效值有效值和和相角相角。2.2.全波傅里叶算法的滤波特性全波傅里叶算法的滤波特性根据式（根据式（4-95）和（）和（4-96）基波分量在第）基波分量在第n点傅里叶正弦和余弦点傅里叶正弦和余弦项系数可表示为：项系数可表示为：4.3.3半波半波傅立叶算法傅立叶算法39/862023/1/20Xian University of Science an

43、d Technology4.3 4.3 基于周期信号的算法基于周期信号的算法Xian University of Science and Technology(4-87)Z Z变换：变换：(4-88)(4-91)(4-90)因此，半波傅立叶正、余弦系数的因此，半波傅立叶正、余弦系数的单位冲激响应序列单位冲激响应序列为：为：4.3.3半波半波傅立叶算法傅立叶算法40/862023/1/20Xian University of Science and Technology4.3 4.3 基于周期信号的算法基于周期信号的算法Xian University of Science and Technol

44、ogy半波傅里叶算法正、余弦系数半波傅里叶算法正、余弦系数频率特性频率特性如图如图4-484-48和和4-494-49所示所示：图图4-48 4-48 正弦项系数频率特性正弦项系数频率特性 图图4-49 4-49 余弦项系数频率特性余弦项系数频率特性4.3.3半波半波傅立叶算法傅立叶算法41/862023/1/20Xian University of Science and Technology4.3 4.3 基于周期信号的算法基于周期信号的算法Xian University of Science and Technology图图4-504-50半波傅里叶算法频率特性半波傅里叶算法频率特性4.

45、3.3半波半波傅立叶算法傅立叶算法42/862023/1/20Xian University of Science and Technology4.3 4.3 基于周期信号的算法基于周期信号的算法Xian University of Science and Technology3.3.全波傅里叶算法与半波傅里叶算法的比较分析全波傅里叶算法与半波傅里叶算法的比较分析半周傅氏算法在故障后半周傅氏算法在故障后10ms即可进行计算，因而使保护的动作即可进行计算，因而使保护的动作时间减少了半个周期。时间减少了半个周期。从滤波效果来看，全波傅氏算法不仅能完全滤出各次谐波分量从滤波效果来看，全波傅氏算法不仅

46、能完全滤出各次谐波分量和稳定的直流分量，和稳定的直流分量，而且能较好地滤出高频分量而且能较好地滤出高频分量半波傅氏算法的滤半波傅氏算法的滤波效果不如全波算法，它不能滤去直流分量和偶次谐波。波效果不如全波算法，它不能滤去直流分量和偶次谐波。从精度来看，当故障发生半周后，半波算法即可计算出真值，从精度来看，当故障发生半周后，半波算法即可计算出真值，但精度差；全波傅氏算法在故障发生一周后才能计算出真值，精度但精度差；全波傅氏算法在故障发生一周后才能计算出真值，精度较半波好。较半波好。由图由图4-504-50可见，半波傅立叶算法可以完全滤除奇次谐波分量，可见，半波傅立叶算法可以完全滤除奇次谐波分量，但

47、不能滤除直流分量、偶次谐波分量和非整次谐波分量。不同初相但不能滤除直流分量、偶次谐波分量和非整次谐波分量。不同初相角时，对于基波和奇次谐波分量的频率响应相同。角时，对于基波和奇次谐波分量的频率响应相同。4.3.3半波半波傅立叶算法傅立叶算法43/864.4 4.4 基于非周期信号模型的基于非周期信号模型的R-LR-L模型算法模型算法4.4.1 4.4.1 线路的线路的R-LR-L模型模型 R-L R-L模型算法仅用于计算线路阻抗，用于距离保护。因而从模型算法仅用于计算线路阻抗，用于距离保护。因而从故障点到保护安装处的线路段可用一电阻和电感的串联来表示，故障点到保护安装处的线路段可用一电阻和电感

48、的串联来表示，即将输电线路等效为即将输电线路等效为R-LR-L简化模型。这样，在短路时，母线电压简化模型。这样，在短路时，母线电压u u和流过保护的电流和流过保护的电流i i与线路的电阻与线路的电阻R1R1和电感和电感L1L1之间有下列方程成之间有下列方程成立立 (4-104)式中：式中：、分别为故障点至保护安装处线路段的正序电阻和分别为故障点至保护安装处线路段的正序电阻和电感；电感；、分别为保护安装处的电压、电流（下面为了简分别为保护安装处的电压、电流（下面为了简便起见，省略掉时间符号便起见，省略掉时间符号t t）式（式（4-854-85）中的）中的u u、i i和和 都是可以测量和计算的，

49、未知数为都是可以测量和计算的，未知数为 和和 。解此微分方程有微分法和积分法两种途径。解此微分方程有微分法和积分法两种途径。4.4.14.4.144/862023/1/20Xian University of Science and Technology4.4 4.4 基于非周期信号模型的基于非周期信号模型的R-LR-L模型算法模型算法4.4.2 R-L4.4.2 R-L模型的微分算法模型的微分算法1.R-L1.R-L模型微分算法原理模型微分算法原理 如果在两个时刻如果在两个时刻 和和 分别测量分别测量 、和和 ，就可以得到，就可以得到两个独立的方程两个独立的方程 (4-105)(4-105)

50、(4-106)(4-106)式中，式中，D D表示表示 ，下标，下标“1 1”和和“2 2”分别表示测量时刻分别表示测量时刻 和和 。为了满足独立方程的需要，要求为了满足独立方程的需要，要求 。联立解以上两式，即可求得两个未知数联立解以上两式，即可求得两个未知数 和和 。(4-107)(4-107)(4-108)(4-108)4.4.4.4.2 245/862023/1/20Xian University of Science and Technology4.4 4.4 基于非周期信号模型的基于非周期信号模型的R-LR-L模型算法模型算法 在用计算机处理时，电流的导数可用差分来近似计算，最在用