ts1=('以下是短路点'); ft=NUm2str(D(r,1)); ts2=('各时刻的结果'); dn=strcat(ts1,ft,ts2); disp(dn); for k=1:n Vb(k)=V0(k,1)-Idb*Z(k,D(r,1)); end for i1=1:nl p=B(i1,1);q=B(i1,2); if p~=0&B(i1,8)==0 if B(i1,6)==0 k=B(i1,5); V(i1)=Vb(p)-Vb(q)./k else k=1./B(i1,5); V(i1)=k*Vb(p)-Vb(q); end else V(i1)=1-Vb(q); end Ib(i1)=V(i1)./B(i1,3); end disp('0s时短路点的电流为'); disp(Idb); disp('0s时各节点的电压标么值为'); disp(Vb);
disp('0s时各支路的电流标么值为'); disp(Ib) for a=1:N Is=zeros(1,nl);Isy=zeros(1,nl);Idzs=0; s(a)=input('请输入短路后几秒（若要终止本程序请输入-1）s='); if s(a)==-1 disp('本程序已被终止，谢谢使用！'); break end for k=1:nl if B(k,1)==0&B(k,8)==1 Sn(k)=B(k,7); cf=Idb./Ib(k);Z1(k)=cf*Z(D(r,1),D(r,1));X(k)=abs(imag(Z1(k))); Xjs1(k)=X(k)*Sn(k)./Sb; if s>4|Xjs1(k)>|B(k,9)==0 Is(k)=1./Xjsl(k); if B(k,9)==0 Isy(k)=Is(k)*Sb./(sq3*vb(D(r,1))); tsl=('无穷大电源'); ft=num2str(k); ts2=('对短路点的计算电抗Xjsl，对应的电流标么值，有名值分别为'); dn=strcat,ts2);disp(dn); disp(Xjsl(k));disp(Is(k));disp(Isy(k)); else Isy(k)=Is(k)*Sn(k)./(sq3*vb(D(r,1))); tsl=('等值电源'); ft=num2str(k); ts2=('对短路点的计算电抗Xjsl，对应的电流标么值，有名值分别为'); dn=strcat(tsl,ft,ts2);
disp(dn); disp(Xjsl(k));disp(Is(k));disp(Isy(k)); end else Xjs=Xjs1(k); tsl=('等值电源'); ft=num2str(k); ts2=('对短路点的计算电抗Xjsl为'); dn=strcat(tsl,ft,ts2);disp(dn); disp(Xjs) Is(k)=input('查与计算电抗Xjs对应的电流标么值I='); Isy(k)=Is(k)*Sn(k)./(sq3*vb(D(r,1))); disp('对应的电流标么值、有名值（单位：KA）分别为');disp(Is(k));disp(Isy(k)); end Idzs=Idzs+Isy(k); end end Id(a)=Idzs; sj=num2str(s(a));ts=('秒对应的短路点短路总电流为（单位：kA）');jg=strcat(sj,ts); disp(jg) disp(Id(a)); end plot(s,Id);xlabel('时间（秒）');ylabel('短路电流大小的标么值');title('短路电流曲线'); end otherwise break end
思考题： 1.计算短路电流的目的是什么 答：选择电气设备,计算各种保护的定值,制定运行方案等 2. 在计算短路电流时，是否用到了节点阻抗阵中的所有元素据此，可以如何简化短路计算程序 答：不需要，可以利用高斯消元法消去部分中间节点和联络节点，化简成等效网络再进行计算。
实验四 潮流计算 1. 用牛顿拉夫逊法求解如下非线性方程组，注意保留迭代过程中的中间值 22??3??x1?x2?5x1?1?0(0) x??2?3? 2????x1?x2?3x2?3?0clc clear all syms x1 x2 f1=x1.^2+x2.^2-5*x1+1;f2=x1.^2+x2.^2-3*x2-3; J(1,1)=diff(f1,x1);J(1,2)=diff(f1,x2); J(2,1)=diff(f2,x1);J(2,2)=diff(f2,x2); % Xn=input('X0=') Xn=[3;3]; a=;b=1;t=0; %迭代次数 while (b) fn1=subs(f1,{x1 x2},{Xn(1,1) Xn(2,1)}); fn2=subs(f2,{x1 x2},{Xn(1,1) Xn(2,1)}); fn(1,1)=fn1;fn(2,1)=fn2; Js(1,1)=subs(J(1,1),{x1 x2},{Xn(1,1) Xn(2,1)});Js(1,2)=subs(J(1,2),{x1 x2},{Xn(1,1) Xn(2,1)}); Js(2,1)=subs(J(2,1),{x1 x2},{Xn(1,1) Xn(2,1)});Js(2,2)=subs(J(2,2),{x1 x2},{Xn(1,1) Xn(2,1)}); Ja=-1*Js; Xns=double(Ja\\fn)
昆明理工大学电力工程学院学生实验报告
