微机继电保护算法
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
微机继电保护算法三个分析题
1设计一个减法数字滤波器，并编程验证其滤波效果MATLAB程序1
clc;
clear;
a1=1;b1=[1 0 0 0 0 0 0 0 -1];
f=0:1:450;
h1=abs(freqz(b1,a1,f,1200));
H1=h1/max(h1);
plot(f,H1);
xlabel('f/Hz');ylabel('H1');
title('滤波器的幅频特性')
MATLAB程序2
clc;
clear;
a1=1;b1=[1 0 0 0 0 0 0 0 -1];
f=0:1:450;
h1=abs(freqz(b1,a1,f,1200));
H1=h1/max(h1);
plot(f,H1);
xlabel('f/Hz');ylabel('H1');
N=24;
t1=(0:0.02/N:0.04);
m=size(t1);
Va=100*sin(2*pi*50*t1);
Va1=25+100*sin(2*pi*50*t1)+40*sin(6*pi*50*t1)+20*sin(6*pi*100*t1); Y=zeros(1,8);
for jj=9:m(2)
Y(jj)=(Va1(jj)-Va1(jj-8))/1.414;
end
subplot(231);
plot(t1,Va,'-ro');
xlabel('t/s');ylabel('v/Va');
title('基波电压Va')
subplot(232);
plot(t1,Va1,'-bs');
xlabel('t/s');ylabel('v/Va1');
title('叠加直流分量和3、6次谐波后的电压Va') subplot(233);
plot(t1,Y,'-g*');
xlabel('t/s');ylabel('v/V');
title('经滤波后的输出电压Y')
subplot(234);
plot(t1,Va,'-ro',t1,Va1,'-bs',t1,Y,'-g*');
xlabel('t/s');ylabel('v/Y');
grid on
title('滤波效果的仿真波形')
2.推导三采样乘积算法，并给出当采样频率f=1200HZ时，参数计算公式及验证MATLAB程序
clc;
clear;
N=24;
t1=(0:0.02/N:0.04);
m=size(t1);
Va=100*sin(2*pi*50*t1);
Ia=50*sin(2*pi*50*t1-pi/6);
for jj=4:m(2)
R(jj)=((Va(jj-1)*Ia(jj-1)+Va(jj-3)*Ia(jj-3)-sqrt(3)*Va(jj-2)*Ia(jj-2))/(Ia(jj-1)^2+Ia(jj-3 )^2-sqrt(3)*Ia(jj-2)^2));
X(jj)=(sqrt(((Va(jj-1)*Ia(jj-3)+Va(jj-3)*Ia(jj-1))*(Va(jj-1)*Ia(jj-3)+Va(jj-3)*Ia(jj-1))-( sqrt(3)*(Va(jj-1)*Ia(jj-2))+(Ia(jj-2)*Ia(jj-1))*(Va(jj-1)*Ia(jj-2))+(Ia(jj-2)*Ia(jj-1)))-(sq rt(3)*(Va(jj-3)*Ia(jj-2)+Va(jj-2)*Ia(jj-3))*(Va(jj-3)*Ia(jj-2)+Va(jj-2)*Ia(jj-3)))))/(Ia(jj -1)^2+Ia(jj-3)^2-sqrt(3)*Ia(jj-2)^2));
Z(jj)=sqrt(((Va(jj-1)*Va(jj-1)+Va(jj-3)*Va(jj-3)-sqrt(3)*Va(jj-2)*Va(jj-2))/(Ia(jj-1)*Ia( jj-1)+Ia(jj-3)*Ia(jj-3)-sqrt(3)*Ia(jj-2)*Ia(jj-2))));
O(jj)=180/pi*atan(sqrt(((Va(jj-1)*Ia(jj-3)+Va(jj-3)*Ia(jj-1))*(Va(jj-1)*Ia(jj-3)+Va(jj-3 )*Ia(jj-1))-(sqrt(3)*(Va(jj-1)*Ia(jj-2))+(Ia(jj-2)*Ia(jj-1))*(Va(jj-1)*Ia(jj-2))+(Ia(jj-2)*I a(jj-1)))-(sqrt(3)*(Va(jj-3)*Ia(jj-2)+Va(jj-2)*Ia(jj-3))*(Va(jj-3)*Ia(jj-2)+Va(jj-2)*Ia(jj-3)))))/(Va(jj-1)*Ia(jj-1)+Va(jj-3)*Ia(jj-3)-sqrt(3)*Va(jj-2)*Ia(jj-2)));
end
subplot(231);
plot(t1,Va,'-ro',t1,Ia,'--bo');
title （'测量到的电压和电流量'）
subplot(232);
plot(t1,R,'-bo');
ylabel('R');
title （'计算得到的电阻值'）
subplot(233);
plot(t1,X,'-bo');
ylabel('X');
title （'计算得到的电抗值'）
subplot(234);
plot(t1,Z,'-bo');
ylabel('Z');
title （'计算得到的阻抗值'）
subplot(235);
plot(t1,O,'-bo');