数字式继电保护基础.pptVIP

4 、三采样值积算法 三采样值积算法是利用三个连续的等时间间隔 T S 的采样值 中两两相乘，通过适当的组合消去 ωt 项以求出 u 、 i 的幅值和其 它电气参数。 设在 t k +1 后再隔一个 T S 为时刻 t k +2 ，此时的 u 、 i 采样值为 ) ( sin S k m 3 T 2 t U u ? ? ? (8-34) ) sin( ? ? ? ? ? ? s k m 3 T 2 t I i (8-35) 上式两采样值相乘，得 )] cos( [cos ? ? ? ? ? ? ? ? s k m m 3 3 T 4 t 2 I U 2 1 i u (8-36) 上式与式 (8-20) 相加，得 )] cos( cos cos [ ? ? ? ? ? ? ? ? ? ? s k s m m 3 3 1 1 T 2 t 2 T 2 2 2 I U 2 1 i u i u 显然，将式 (8-37) 和式 (8-21) 经适当组合以消去 ωt k 项，得 s 2 s 2 2 3 3 1 1 m m T 2 T 2 i u 2 i u i u I U ? ? ? sin cos cos ? ? ? 若要 ωT s ＝ 30 o ，上式简化为 ) ( cos 2 2 3 3 1 1 m m i u i u i u 2 I U ? ? ? ? 用 I m 代替 U m （或 U m 代替 I m ），并取 ＝ 0 o ，则有 ) ( 2 2 2 3 2 1 m u u u 2 U ? ? ? ) ( 2 2 2 3 2 1 m i i i 2 I ? ? ? (8-40) (8-41) 由式 (8-39) 和式 (8-41) 可得 2 2 2 3 2 1 2 2 3 3 1 1 m m i i i i u i u i u I U R ? ? ? ? ? ? ? cos (8-42) 由式 (8-27) 和式 (8-41) ，并考虑到，得 2 2 2 3 2 1 1 2 2 1 m m i i i i u i u I U X ? ? ? ? ? ? sin (8-43) 由式 (8-40) 和式 (8-41) 得 2 2 2 3 2 1 2 2 2 3 2 1 m m i i i u u u I U z ? ? ? ? ? ? (8-44) 由式 (8-42) 和式 (8-43) 得 (8-45) 2 2 3 3 1 1 2 2 1 1 i u i u i u i u i u arctg ? ? ? ? ? 三采样值积算法的数据窗是 2 T s 。 从精确角度看，如果输入 信号波形是纯正弦的，这种算法没有误差，因为算法的基础是 考虑了采样值在正弦信号中的实际值。 三、傅里叶算法（傅氏算法） 1. 全周波傅里叶算法 全周波傅里叶算法是采用由 cos nω 1 t 和 sin nω 1 t ( n =0 ， 1 ， 2 … ) 正弦函数组作为样品函数，将这一正弦样品函数与待分 析的时变函数进行相应的积分变换，以求出与样品函数频率 相同的分量的实部和虚部的系数。进而可以求出待分析的时 变函数中该频率的谐波分量的模值和相位。 根据傅里叶级数，我们将待分析的周期函数电流信号 i(t) 表示为 ? ? t n I t n I I t i 1 1 n ns 1 1 n nc 0 ? ? sin cos ? ? ? ? ? ? ? ? ? 式中 n —— n 次谐波（ n =1 ， 2 ， … ）； I 0 —— 恒定电流分量； I nc 、 I ns —— 分别表示 n 次谐波的余弦分量电流和正弦 分量电流的幅值。 (8-46) 当我们希望得到 n 次谐波分量时，可用 cos nω 1 t 和 sin nω 1 t 分别乘式 (8-46) 两边，然后在 t 0 到 t 0 ＋ T 积分，得到 dt t n t i T 2 I T t t 1 nc 0 0 ? ? ? ? cos ) ( dt t n t i T 2 I T t t 1 ns 0 0 ? ? ? ? sin ) ( (8-47) (8-48) 每工频周期 T 采样 N 次，对式 (8-47) 和式 (8-48) 用梯形法 数值积分来代替，则得 N n 2 k i N 2 I N 1 k k nc ? cos ? ? ? N n 2 k i N 2 I N 1 k k ns ? sin ? ? ? (8-49) (8-50) 式中 k 、 i k —— 第 k 采样及第 k 个采样值 电流 n 次谐波幅值（最大值）和相位（余弦函数的初相）分别为 2 nc 2 ns nm I I I ? ? (8-51) (8-52) nc ns n I I arctg ?