%短路計算
clc;clear?all;
SNG=50;%發電機容量
E2=10.5;%發電機電動勢
Xd2=0.2;%發電機次暫態電抗
X1=0.4;%單位長度正序電抗
X0=3*X1;%單位長度零序電抗
L=80;%線路長度
SB=100;%選擇基準容量為100MVA
UN=110;%額定電壓為110kv
Uav=1.05*UN;
UB=Uav;?%基準電壓取額定電壓1.05倍
IB=SB/sqrt（3）/UB;%計算基準電流
SNT=50;%變壓器容量
Uk=10.5;%變壓器短路電壓百分數
XG1=Xd2*SB/SNG;%發電機正序電抗標幺值
XG2=XG1;%發電負序電抗
XT1=Uk/100*SB/SNT;%變壓器正序電抗標幺值
XT2=XT1;%變壓器負序電抗
XT0=XT1;%變壓器零序電抗


????disp（‘---------------------------‘）
????Zd=Z2*（Z0+3*Zg）/（Z2+Z0+3*Zg）;
????M=sqrt（3）*sqrt（1-Z2*（Z0+3*Zg）/（Z2+Z0+3*Zg）^2）;
????If1=j*Uf0/（Z1+Zd）;
????If2=-If1*Z0/（Z2+Z0）;
????If0=-If1*Z2/（Z2+Z0）;???
end
if?Type_shortage==3
????disp（‘短路類型：兩相短路‘）
????disp（‘---------------------------‘）
????Zd=Z2+Zf;
????M=sqrt（3）;
????If1=j*Uf0/（Z1+Zd）;
????If2=-If1;
????If0=0;
end
if?Type_shortage==4
????disp（‘短路類型：三相短路‘）
????disp（‘---------------------------‘）
????Zd=0;
????M=1;
????If1=j*Uf0/（Z1+Zd）;
????If2=0;
????If0=0;
end

If=M*If1;??%根據M值及正序電流計算流入大地的短路電流?
disp（‘正序電流‘）
disp（If1）
disp（‘負序電流‘）
disp（If2）
disp（‘零序電流‘）
disp（If0）
disp（‘短路電流‘）
disp（If）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-10-10?15:39??源碼\
?????文件????????1277??2018-10-13?16:02??源碼\shortage123.m
?????文件??????790016??2017-03-28?22:46??源碼\論文.doc
?????目錄???????????0??2018-10-10?15:39??演示\
?????文件????26502656??2017-03-28?22:05??演示\10716論文.avi
?????文件??????????64??2018-10-10?16:34??說明.txt
?????目錄???????????0??2018-10-10?15:38??需求\
?????文件???????11604??2017-03-13?17:08??需求\功能說明.docx