信息工程学院(课程设计)
end
disp('迭代次数为');disp(Times) for no=1:n
Vn(no)=B2(no,3);
str1 = sprintf('节点%d的电压是%d',no); disp(str1);disp(Vn(no)); end
2.潮流计算输入参数说明： n=4 n1=4
isb=4 pr=0.00001
B1=[1 2 0.1+0.4i 0.3056i 1 0;1 3 0+0.3i 0 1.1 0;1 4 0.12+0.5i 0.0382i 1 0;2 4 0.08+0.4i 0.02826i 1 0]
%各支路参数:首末节点号，支路阻抗，对地容抗，变压器变比，是否含变压器 B2=[0 -0.3-0.18i 1 0 0 2;0 -0.55-0.13i 1 0 0 2;0 0.5+0i 1 1.1 0 3;0 0 1 1.05 0 1]
%无功补偿容量，各节点注入功率，初始电压，节点电压给定值和相角，节点类型，\表示平衡节点，
X=[1 0;2 0;3 0;4 0] %\表示PQ节点，\表示PV节点
%X=input('节点号和对地参数:X=');
3.不对称短路的m语言源程序：
%求发电机参数的标幺值,计算公式：X=Xd1*(SB/SGN) clear
Sn=120;SB=120;Xdc1=0.9;Xdc2=0.45; XG1b=Xdc1*(SB/Sn);
disp('一.发电机1的电抗值 XG1b='),disp(XG1b) XG2b=Xdc2*(SB/Sn);
disp('发电机2的电抗值 XG2b='),disp(XG2b)
%求变压器T的各序等值电路电抗的参数，计算公式：XT=Uk%/100*(SB/STN) clear
ST1=60;ST2=60;SB=120;Uk1=10.5;Uk2=10.5; XT1=(Uk1/100)*(SB/ST1);
disp('二.变压器T的各序电抗 XT1='),disp(XT1) XT2=(Uk2/100)*(SB/ST2); disp('XT2='),disp(XT2)
XT0=XT1; %由于变压器是静止电器，所以各序参数相等 disp('XT0='),disp(XT0)
%输电线1的各序等值电路的电抗标幺值计算，计算公式：
46
信息工程学院(课程设计)
Xl1=Xl2=x0*(SB/Uav^2),Xl0=3*Xl1 clear
SB=120;x0=0.4;L=105;Uav=115; Xl1b=x0*L*(SB/Uav^2);
disp('三.输电线的各序电抗 Xl1b='),disp(Xl1b) Xl2b=x0*L*(SB/Uav^2); disp('Xl2b='),disp(Xl2b) Xl0b=3*x0*L*(SB/Uav^2); disp('Xl0b='),disp(Xl0b)
%负荷1的各序等值电路的电抗标幺值。计算公式：X1L=X1*(SB/Sn) clear
SB=120;Sn=60;X11=1.2;X12=0.35; X1L1b=X11*(SB/Sn);
disp('四.负荷1各序参数 X1L1b='),disp(X1L1b) X1L2b=X12*(SB/Sn);
disp('XlL2b='),disp(X1L2b)
%负荷2的各序等值电路的电抗标幺值。计算公式：X2L=X1*(SB/Sn) SB=120;Sn=40;X21=1.2;X22=0.35; X2L1b=X21*(SB/Sn);
disp('负荷2各序参数 X2L1b='),disp(X2L1b) X2L2b=X22*(SB/Sn);
disp('X2L2b='),disp(X2L2b)
%单相接地短路时的短路电流计算
Uf0=0.95;Zkk1=0.83;Zkk2=0.44;Zkk0=0.78; SB=120;UB=115;
If1=SB/(UB*sqrt(3))*Uf0/(Zkk1+Zkk2+Zkk0); disp('单相短路电流的正序分量If1='),disp(If1) If=3*If1;
disp('单相短路电流If='),disp(If)
%两相短路时的短路电流计算
Uf0=0.95;Zkk1=0.83;Zkk2=0.44;Zkk0=0.78; SB=120;UB=115;
If1=SB/(UB*sqrt(3))*Uf0/(Zkk1+Zkk2);
disp('两相短路电流的正序分量If1='),disp(If1) If=sqrt(3)*If1;
disp('两相短路电流If='),disp(If)
%两相接地短路时的短路电流计算
Uf0=0.95;Zkk1=0.83;Zkk2=0.44;Zkk0=0.78; SB=120;UB=115;
IB=SB/(UB*sqrt(3));
47
信息工程学院(课程设计)
Zzeta=Zkk2*Zkk0/(Zkk2+Zkk0);
m(1,1)=sqrt(3)*sqrt(1-(Zkk2*Zkk0/((Zkk2+Zkk0)*(Zkk2+Zkk0)))); If1=SB/(UB*sqrt(3))*Uf0/(Zkk1+Zzeta);
disp('两相接地短路电流的正序分量If1='),disp(If1) If=m(1,1)*If1;
disp('两相接地短路电流If='),disp(If)
48