三机九节点潮流计算(共4页).docx

精选优质文档-倾情为你奉上clc;clear;close all;n1=9;n2=9;B1=load(b1.txt,ascii);B2=load(b2.txt,ascii);m=0;r=0;for i=1:n1 if B1(i,8)=1 m=m+1; endendr=n1-m-1;e=B1(1:n1,2);f=B1(1:n1,3);ps=B1(1:n1,4)-B1(1:n1,6);qs=B1(1:n1,5)-B1(1:n1,7);Y=zeros(n1);for j=1:n2 p=B2(j,2);q=B2(j,3);J=zeros(2*n1); Y(p,q)=-1./(B2(j,4)+1i*B2(j,5)*B2(j,7); Y(q,p)=Y(p,q); Y(p,p)=Y(p,p)+1/(B2(j,4)+1i*B2(j,5)+1i*B2(j,6); Y(q,q)=Y(q,q)+1/(B2(j,