继电保护仿真作业

继电保护技术与应用

xlabel('t/s');ylabel('三相电压波形u/V');legend('ua','ub','uc');

subplot(212);

plot(t,ia,'k',t,ib,'y:',t,ic,'r-'); xlabel('t/s');ylabel('三相电流波形i/A');legend('ia','ib','ic');

%滤波

for e=5:3001

t1(1,e-4)=t(1,e);

Ia(1,e-4)=(ia(1,e)+sqrt(3)*ia(1,e-1)+2*ia(1,e-2)+sqrt(3)*ia(1,e-3)+ia(1,e-4))/6;

Ib(1,e-4)=(ib(1,e)+sqrt(3)*ib(1,e-1)+2*ib(1,e-2)+sqrt(3)*ib(1,e-3)+ib(1,e-4))/6;

Ic(1,e-4)=(ic(1,e)+sqrt(3)*ic(1,e-1)+2*ic(1,e-2)+sqrt(3)*ic(1,e-3)+ic(1,e-4))/6;

Ua(1,e-4)=(ua(1,e)+sqrt(3)*ua(1,e-1)+2*ua(1,e-2)+sqrt(3)*ua(1,e-3)+ua(1,e-4))/6;

Ub(1,e-4)=(ub(1,e)+sqrt(3)*ub(1,e-1)+2*ub(1,e-2)+sqrt(3)*ub(1,e-3)+ub(1,e-4))/6;

Uc(1,e-4)=(uc(1,e)+sqrt(3)*uc(1,e-1)+2*uc(1,e-2)+sqrt(3)*uc(1,e-3)+uc(1,e-4))/6;

end

figure(2),subplot(1,2,1);

plot(t,ua);xlabel('t/s');ylabel('ua/V');legend('滤波前ua');

subplot(1,2,2);

plot(t1,Ua);xlabel('t/s');ylabel('ua/V');legend('滤波后Ua');

figure(3),subplot(1,2,1); plot(t,ia);legend('滤波前

ia');xlabel('t/s');ylabel('ia/A');

subplot(1,2,2);

plot(t1,Ia);legend('滤波后Ia'); xlabel('t/s');ylabel('Ia/A'); %两点乘积算法 for n=1440:1679

u1=Ua(n-1)/2+Ua(n-2)/2; u2=Ua(n+29)/2+Ua(n+28)/2;

i1=Ia(n-1)/2+Ia(n-2)/2; i2=Ia(n+29)/2+Ia(n+28)/2;

r(n-1439,1)=(u1*i1+u2*i2)/(i2^2+i1^2);

x(n-1439,1)=(u1*i2-u2*i1)/(i2^2+i1^2);

end

figure(8),plot(x,'k');legend('电抗x');

4 参考文献

[1] 张保会，尹项根. 电力系统继电保护[M].中国电力出版社.2005.5 130-131 [2] 杨奇逊，黄少锋.微型机继电保护基础[M].中国电力出版社.2007.8

[3]陈皓.微机保护原理及算法仿真[M].中国电力出版社.2007.1 159-166.

[4]关敬欢.电力系统继电保护现状与发展探讨[J].现代商贸工业，2009年，第18期：303-304.

[5]刘平，刘鹏龙.基于Matlab的电力系统故障分析与仿真[A]. 2007系统仿真技术及其应用学术研讨会[C]. 中国自动化学会,2007年8月，427-429.

[6]胡文丽.基于录波信息的故障特征量分析方法的研究[D].华北电力学,2004年. [7]霍兰茹.继电保护通用分析和仿真校验系统的研究[D].华北电力大学,2009年. [8] 朱广伟.微机继电保护在企业供电系统中的应用及发展趋势[A].第4届中国金属协会青年学术年会论文集[C].