微机保护算法仿真示例.pptxVIP

算法仿真示例;C4example1：数字滤波器设计;选择以下数字滤波器：

H1(z)=（1–Z-m/4)2

H2(z)=[（1–Z–7)/（1–Z–1)]2

m为每周期的采样点数，取m=96；

H1(z)是有限种激响应数字滤波器，为有梳状频谱特性的简单数字滤波器，主要用来初步消去整次谐波及直流分量的影响；

H2(z)为无限冲激响应数字滤波器，采用递归形式实现，是低通数字滤波器，能较好地消去高频分量的影响。(?！);分析：

H1(z)=（1–Z-m/4)2

结构：相减滤波单元的级联;

阶次:k=96/4=24

作用：使用相减滤波单元消除直流分量和

4，8，12，…l?4次谐波分量

类型：有限冲激响应数字滤波器

（FIR滤波器）;分析：

H2(z)=[（1–Z–7)/（1–Z–1)]2

结构：积分滤波单元的级联;

阶次:k=8

作用：使用积分滤波单元增加对高次谐波

的衰耗

类型：有限冲激响应数字滤波器

（FIR滤波器；论文误为IIR滤波器）

;N=96;f0=50;fs=N*f0;

a1=1;b1=[1zeros(1,23)-1];

%设置差分滤波单元的传递函数系数%k=N/m=24;(谐波阶次m=4)

%文献给出的积分滤波单元的传递函数

a2=1;b2=[ones(1,9)];

%设置积分滤波单元的传递函数系数k=8,

%另外注意：设置系数b时，b0=1;

f=0:1:fs/2;

h1=abs(freqz(b1,a1,f,fs));

H1=h1/max(h1);

%由传递函数系数确定传递函数的幅频特性;H11=H1.^2;

%级联差分滤波器

h2=abs(freqz(b2,a2,f,fs));

H2=h2/max(h2);

H22=H2.^2;

%级联积分滤波器

h=h1.*h2;H=h./max(h);

%计算级联滤波器传递函数的幅频特性

HH=H.^2;

;Figure

%级联差分、积分滤波器幅频特性

subplot(221);plot(f,H11);

xlabel(f/Hz);ylabel(H1);

subplot(223);plot(f,H22);

xlabel(f/Hz);ylabel(H2);

subplot(222);plot(f,HH);

xlabel(f/Hz);ylabel(H);

subplot(224);plot(f,20*log10(HH+10e-10));

axis([0,fs/2,-60,0]);xlabel(f/Hz));ylabel(H/dB);;分析：

H2(z)=[（1–Z–7)/（1–Z–1)]2

并未针对特定谐波滤波进行设置，最终结

果非常不理想：实际上二次谐波有最大输

出；且3、5、6、7等各次滤波也有很大输

出。

结论：

在N=96时，论文给出的滤波器不能作为工频

选频滤波器。;改进：

考虑积分滤波单元1消除l?3次谐波的影响:

m=3

阶次:k=N/m–1=31

考虑积分滤波单元2消除l?2次谐波的影响:

m=2

阶次:k=N/m–1=47

H2(z)=（1–Z–31)（1–Z–47)/（1–Z–1)2

;仿真m文件作以下修改：

N=96;f0=50;fs=N*f0;

a1=1;b1=[1zeros(1,23)-1];

%设置差分滤波单元的传递函数系数k=24

a2=1;b2=[ones(1,32)];

%设置积分滤波单元的传递函数系数k=31

a3=1;b3=[ones(1,48)];

%设置积分滤波单元的传递函数系数k=47

重新运行m文件，结果如下：

（参考c4e2.m）;滤波效果比较

原始输入信号

y0=sin(2*pi*f0*t)+0.5*sin(2*2*pi*f0*t)+0.5*

sin(3*2*pi*f0*t)+0.3*sin(4*2*pi*f0*t)+0.3*

sin(5*2*pi*f0*t);

滤波输出信号