소개글
Matlab 을 이용한 전력 고장계산
목차
1.내용
2.예시
3.수행 및 결과화면
4,코드
본문내용
Matlab 을 이용한 전력 고장계산
내용
모선(Bus)사고 와 선로(Line)사고 시에 사고 종류를 입력받아
1선 지락 (single line ground fault), 2선 지락 (double line ground fault), 단락 사고 (Line to line fault), 3상평형사고 (three phase fault)
각 사고 시에 고장 전류, 전압 및 Line flow를 구해 고장을 해석하기 위한 프로그램
예시 계통
발전기 데이터
<중 략>
ia=zeros(q,q);
ib=zeros(q,q);
ic=zeros(q,q);
for k=1:q
l=data1(k,1);
m=data1(k,2);
ia(l,m)=[1,1,1]*[i0f(l,m);i1f(l,m);i2f(l,m)];
ib(l,m)=[1,a^2,a]*[i0f(l,m);i1f(l,m);i2f(l,m)];
ic(l,m)=[1,a,a^2]*[i0f(l,m);i1f(l,m);i2f(l,m)];
end
iF=iF(2)+iF(3);
fprintf(`\n\n고장전류> magnitude=%.4f degree=%.2f\n`,abs(iF),angle(iF)*180/pi)
fprintf(`\n고장전압\n`)
fprintf(`N0. Bus voltage A voltage B voltage C angle A angle B angle C \n`)
for k=1:ysize+1
fprintf(` %d %f %f %f %.2f %.2f %.2f\n`,k ,abs(vabcf(k,1)),abs(vabcf(k,2)),abs(vabcf(k,3)),angle(vabcf(k,1))*180/pi,angle(vabcf(k,2))*180/pi,angle(vabcf(k,3))*180/pi)
end
fprintf(`\nLine current flow\n`)
fprintf(`\n from bus to bus phase current A phase current B phase current C\n`)
for k=1:ysize+1
for l=1:ysize+1
if abs(ib(k,l))|abs(ia(k,l))~=0
fprintf(` %d %d %f %f %f \n`,k,l,abs(ia(k,l)),abs(ib(k,l)),abs(ic(k,l)));
참고 자료
없음