NF=1;
n=4;
nl=6;
Sb=100;
D=[4,0];
vb=[10.5 115 6.3 37];
V0=[1;1;1;1];
B=[0 1 0.15i 0 1 0 70 1 1;
0 2 6.67i 0 1 0 0 0 2;
0 3 0.71i 0 1 0 15 1 1;
1 2 0.26i 0 1 0 0 0 2;
2 3 1.11i 0 1 1 0 0 2;
3 4 0.53i 0 1 1 0 0 2];
N=2;
m=0;Z=zeros(n);sq3=sqrt(3);s=zeros(N,1);
for k1=1:nl
p=B(k1,1);q=B(k1,2);
if B(k1,6)==0
k=1./B(k1,5);
else k=B(k1,5);
end
if p==0
if q>m %追加接地树支
Z(q,q)=B(k1,3);m=m+1;
else %追加接地连支
for i1=1:m
Z(i1,m+1)=-Z(i1,q); Z(m+1,i1)=- Z(q,i1);
end
Z(m+1,m+1)=Z(q,q)+B(k1,3);
for i1=1:m
for j=1:m
Z(i1,j)=Z(i1,j)-Z(i1,m+1)*Z(m+1,j)./Z(m+1,m+1);
end
Z(i1,m+1)=0;
end
for i1=1:m+1
Z(m+1,i1)=0;
end
end
else if q>m %追加不接地树支
for i1=1:m
Z(i1,q)=Z(i1,p)*k;Z(q,i1)=Z(p,i1)*k;
end
Z(q,q)=k^2*Z(p,p)+k^2*B(k1,3);