OpenFile = uigetfile('*.*');
All_D = textread(OpenFile);
BusNum = All_D(1,1); BranchNum = All_D(1,2);
L = 0;
for j = 3:BranchNum + 3
if All_D(j,1) == 0
break;
end
L = L + 1; M=L+2;
end
LineI = All_D( 3:M,3 ); LineJ = All_D( 3:M,4 ); LineX = All_D( 3:M,5); LineX0 = All_D(3:M,6);
M=M+2;L=0;
for j = M:BranchNum + M
if All_D(j,1) == 0
break;
end
L = L + 1;
end
N=M+L-1;
Linei = All_D( M:N,1); LineT = All_D(M:N,2); LineF = All_D(M:N,3);
if size( LineI ) >= 1
TempY1 = 1 ./ LineX.* 1i ; TempB1 = -1 ./ LineT.* 1i ; TempB2 = -1 ./ LineF.* 1i ;
k=1;
for j=1:BranchNum
if LineX0(k)>0|| LineX0(k)<0
TempY0 (k)= 1 / LineX0(k)* 1i;
end
k=k+1;
end
Y1 = sparse( [ LineI;LineJ ],[ LineJ;LineI ],[ TempY1,TempY1],22,22);
Y0 = sparse( [ LineI;LineJ ],[ LineJ;LineI ],[ TempY0,TempY0],22,22);
Y0 = Y0 + diag( sum( -Y0 ) );Y=Y0;
Y1 = Y1 + diag( sum( -Y1 ) ); Y2=Y1;
B1=sparse( Linei,Linei ,TempB1,22,22);
B2=sparse( Linei,Linei ,TempB2,22,22);
Y1 =Y1+B1;Y2=Y2+B2;
end
k=1;
for j=1:BusNum
if Y0(k,k)==0
Y0(k,k)=1i;
end
k=k+1;
end
k=1;
for j=1:L
Y0( Linei(k), : ) = 0;
Y0(: , Linei(k) ) = 0;
Y0( Linei(k), Linei(k) ) = 1* 1i;
k=k+1;
end
Z0=inv(Y0); Z1=inv(Y1); Z2=inv(Y2);
a=zeros(22,1);
a(e)=1;b1=Z1*a;b2=Z2*a;b0=Z0*a;
n=eye(22,22);
n=sum(n);
n=n';
T=[1,1,1;-0.5-0.866i,-0.5+0.866i,1;-0.5+0.866i,-0.5-0.866i,1];
p=-0.5-0.866i;q=-0.5+0.866i;