资源描述：
《电力系统短路电流计算-matlab程序.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、%电力系统极坐标下的牛顿-拉夫逊法潮流计算disp('电力系统极坐标下的牛顿-拉夫逊法潮流计算:');clearn=input('请输入结点数：n=');n1=input('请输入PV结点数：n1=');n2=input('请输入PQ结点数：n2=');isb=input('请输入平衡结点：isb=');pr=input('请输入精确度：pr=');K=input('请输入变比矩阵看:K=');C=input('请输入支路阻抗矩阵：C=');y=input('请输入支路导纳矩阵：y=');U=input('请输入结点电压矩阵：U=');S=input('请输入各
2、结点的功率：S=');Z=zeros(1,n);N=zeros(n1+n2,n2);L=zeros(n2,n2);QT1=zeros(1,n1+n2);form=1:nforR=1:nC(m,m)=C(m,m)+y(m,R);ifK(m,R)~=0C(m,m)=C(m,m)+1/((K(m,R)*C(m,R))/(K(m,R)-1));C(R,R)=C(R,R)+1/((K(m,R)^2*C(m,R))/(1-K(m,R)));C(m,R)=C(m,R)*K(m,R);C(R,m)=C(m,R);endendendform=1:nforR=1:nifm~=RZ(
3、m)=Z(m)+1/C(m,R);endendendform=1:nforR=1:nifm==RY(m,m)=C(m,m)+Z(m);elseY(m,R)=-1/C(m,R);endendenddisp('结点导纳矩阵:');disp(Y);disp('迭代中的雅克比矩阵:');G=real(Y);B=imag(Y);O=angle(U);U1=abs(U);k=0;PR=1;P=real(S);Q=imag(S);whilePR>prform=1:n2UD(m)=U1(m);endform=1:n1+n2forR=1:nPT(R)=U1(m)*U1(R)*(G
4、(m,R)*cos(O(m)-O(R))+B(m,R)*sin(O(m)-O(R)));endPT1(m)=sum(PT);PP(m)=P(m)-PT1(m);PP1(k+1,m)=PP(m);endform=1:n2forR=1:nQT(R)=U1(m)*U1(R)*(G(m,R)*sin(O(m)-O(R))-B(m,R)*cos(O(m)-O(R)));endQT1(m)=sum(QT);QQ(m)=Q(m)-QT1(m);QQ1(k+1,m)=QQ(m);endPR1=max(abs(PP));PR2=max(abs(QQ));PR=max(PR1,PR
5、2);form=1:n1+n2forR=1:n1+n2ifm==RH(m,m)=U1(m)^2*B(m,m)+QT1(m);elseH(m,R)=-U1(m)*U1(R)*(G(m,R)*sin(O(m)-O(R))-B(m,R)*cos(O(m)-O(R)));endendendform=1:n1+n2forR=1:n2ifm==RN(m,m)=-U1(m)^2*G(m,m)-PT1(m);elseN(m,R)=-U1(m)*U1(R)*(G(m,R)*cos(O(m)-O(R))+B(m,R)*sin(O(m)-O(R)));endendendform=1:
6、n2forR=1:n1+n2ifm==RJ(m,m)=U1(m)^2*G(m,m)-PT1(m);elseJ(m,R)=U1(m)*U1(R)*(G(m,R)*cos(O(m)-O(R))+B(m,R)*sin(O(m)-O(R)));endendendform=1:n2forR=1:n2ifm==RL(m,m)=U1(m)^2*B(m,m)-QT1(m);elseL(m,R)=-U1(m)*U1(R)*(G(m,R)*sin(O(m)-O(R))-B(m,R)*cos(O(m)-O(R)));endendendJJ=[HN;JL];disp(JJ);PQ=[P
7、P';QQ'];DA=-inv(JJ)*PQ;DA1=DA';form=1:n1+n2OO(m)=DA1(m);endform=n:n1+n2+n2UU1(m-n1-n2)=DA1(m);endUD2=diag(UD);UU=UU1*UD2;form=1:n1+n2O(m)=O(m)+OO(m);endform=1:n2U1(m)=U1(m)+UU(m);endform=1:n1+n2o(k+1,m)=180/pi*O(m);endform=1:n2u(k+1,m)=U1(m);endk=k+1;endform=1:nb(m)=U1(m)*cos(O(m));c
8、(m)=U1(m)*si
