短路计算matlab程序文档格式.docx
- 文档编号:15161918
- 上传时间:2022-10-28
- 格式:DOCX
- 页数:15
- 大小:17.76KB
短路计算matlab程序文档格式.docx
《短路计算matlab程序文档格式.docx》由会员分享,可在线阅读,更多相关《短路计算matlab程序文档格式.docx(15页珍藏版)》请在冰豆网上搜索。
Positive;
Negative;
Zero;
%runpowerflow:
opt=mpoption('
ENFORCE_Q_LIMS'
1);
results=runpf(case9gs,opt);
Vp=[results.bus(:
8).*(cos(pi/180*results.bus(:
9))+j*sin(pi/180*results.bus(:
9)))];
clc;
%Menu:
Bucle=1;
while(Bucle==1)
disp('
________________________________________________________________________________________'
)
%%%%%%Programtocalculatethesymetricalandunsymetricalshortcircuitfaults%%%%%%'
);
1)Symetricalfault(threephaseshortcircuitfault)'
2)Unsymetricalfault(Asinglephasetogroundfault)'
0)Forexittheprogram'
obtion=input('
Choosetheoptionyouwanttobecalculatedbytypingitsnumber:
'
'
s'
%Calculationofthefault
switch(obtion)
%%
case'
1'
disp('
YouHavechoose:
SYMETRICALFAULT'
'
Pleasefollowthemenu'
%Askforthefault'
sbusnumberandtheRfvalue:
*******************************************************************************'
*DATA*'
bus=input('
Please,numberthebuswherethethreephaseshortcircuitfaultoccur.Busnumber:
%bus=bus+1;
%Sincenodezeroisalsoanalyzed,thebusnumberXwillhaveaX+1index
Rf=input('
PleasewritethevalueofthearcresistanceRfperphase.Rf='
%Calculation:
If=zeros(8,1);
If(bus)=Vp(bus)/(Zpos(bus,bus)+Rf);
Vs=Vp-Zpos*If;
%ShowingtheResults:
*RESULTS*'
FaultCurrent=sprintf('
Theshortcircuitcurrentsis:
%gA.withangle%g?
Iref*abs(If(bus)),180/pi*angle(If(bus)));
disp(FaultCurrent);
figure();
compass(If(bus),'
blue'
holdon
compass(If(bus)*a,'
red'
compass(If(bus)*a^2,'
green'
text(-abs(If(bus)),-abs(If(bus)),'
ShortCircuitcurrent'
disp('
TheVoltagesatallthebusesinthesystemare:
_____________________________________________'
BusVAngle[Degree]'
VoltageTokke=sprintf('
Tokke%gV%g?
Uref*abs(Vs
(1)),180/pi*angle(Vs
(1)));
disp(VoltageTokke)
VoltageVinje=sprintf('
Vinje%gV%g?
Uref*abs(Vs
(2)),180/pi*angle(Vs
(2)));
disp(VoltageVinje)
VoltageSonga=sprintf('
Songa%gV%g?
Uref*abs(Vs(3)),180/pi*angle(Vs(3)));
disp(VoltageSonga)
VoltageVemork=sprintf('
Vemork%gV%g?
Uref*abs(Vs(4)),180/pi*angle(Vs(4)));
disp(VoltageVemork)
VoltageRjukan=sprintf('
Rjunkan%gV%g?
Uref*abs(Vs(5)),180/pi*angle(Vs(5)));
disp(VoltageRjukan)
VoltageFlesaker=sprintf('
Flesaker%gV%g?
Uref*abs(Vs(6)),180/pi*angle(Vs(6)));
disp(VoltageFlesaker)
VoltageTveiten=sprintf('
Tveiten%gV%g?
Uref*abs(Vs(7)),180/pi*angle(Vs(7)));
disp(VoltageTveiten)
VoltageRod=sprintf('
Rod%gV%g?
Uref*abs(Vs(8)),180/pi*angle(Vs(8)));
disp(VoltageRod)
figure()
subplot(2,4,1)
compass(Vs
(1),'
b'
compass(Vs
(1)*a,'
r'
compass(Vs
(1)*a^2,'
g'
text(-abs(Vs
(1)),-abs(Vs
(1)),'
Tokke`sVoltagevectors'
subplot(2,4,2)
compass(Vs
(2));
compass(Vs
(2)*a,'
compass(Vs
(2)*a^2,'
text(-abs(Vs
(2)),-abs(Vs
(2)),'
Vinje`sVoltagevectors'
subplot(2,4,3)
compass(Vs(3));
compass(Vs(3)*a,'
compass(Vs(3)*a^2,'
text(-abs(Vs(3)),-abs(Vs(3)),'
Songa`sVoltagevectors'
subplot(2,4,4)
compass(Vs(4));
compass(Vs(4)*a,'
compass(Vs(4)*a^2,'
text(-abs(Vs(4)),-abs(Vs(4)),'
Vemork`sVoltagevectors'
subplot(2,4,5)
compass(Vs(5));
compass(Vs(5)*a,'
compass(Vs(5)*a^2,'
text(-abs(Vs(5)),-abs(Vs(5)),'
Rjukan`sVoltagevectors'
subplot(2,4,6)
compass(Vs(6));
compass(Vs(6)*a,'
compass(Vs(6)*a^2,'
text(-abs(Vs(6)),-abs(Vs(6)),'
Flesaker`sVoltagevectors'
subplot(2,4,7)
compass(Vs(7));
compass(Vs(7)*a,'
compass(Vs(7)*a^2,'
text(-abs(Vs(7)),-abs(Vs(7)),'
Tveiten`sVoltagevectors'
subplot(2,4,8)
compass(Vs(8));
compass(Vs(8)*a,'
compass(Vs(8)*a^2,'
text(-abs(Vs(8)),-abs(Vs(8)),'
Rod`sVoltagevectors'
*CUR
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 短路 计算 matlab 程序