clear;
y (1,2)=1/(0。06+0。18i); y (1,3)=1/(0。06+0。18i); y (1,4)=1/(0。04+0。12i);
y(1,5)=1/(0。02+0。06i);
y(2,3)=1/(0。01+0。
03i);y(2,5)=1/(0。08+0。24i);
y(3,4)=1/(0。08+0。24i);
y(4,5)=0;
for i=1:5
for j=i:5
y(j,i)=y(i,j);
end
end
Y=0;
for i=1:5
for j=1:5
if i~=j
Y(i,j)=-y(i,j);
end
end
end for i=1:5
Y(i,i)=sum(y(i,:));
end
Y
G=real(Y);
B=imag(Y);
S(1)=0。
2+0。2i;
S(2)=-0。45-0。15i;
S(3)=-0。4-0。05i;
S(4)=-0。6-0。1i;
S(5)=0;
P=real(S);
Q=imag(S);
U=ones(1,5);
U(5)=1。
06;
e=zeros(1,5);
ox=ones(8,1);
fx=ones(8,1);
count=0
while max(fx)>1e-5
for i=1:4
for j=1:4
H(i,j)=0;N(i,j)=0;M(i,j)=0;L(i,j)=0;oP(i)=0;oQ(i)=0;
end
end
for i=1:4
for j=1:5
oP(i)=oP(i)-U(i)*U(j)*(G(i,j)*cos(e(i)-e(j))+B(i,j)*sin(e(i)-e(j)));
oQ(i)=oQ(i)-U(i)*U(j)*(G(i,j)*sin(e(i)-e(j))-B(i,j)*cos(e(i)-e(j)));
end
oP(i)=oP(i)+P(i); oQ(i)=oQ(i)+Q(i);
end
fx=[oP,oQ]';
%当i~=j时候求H,N,M,L 如下:
for i=1:4
for j=1:4
if i~=j H(i,j)=-U(i)*U(j)*(G(i,j)*sin(e(i)-e(j))-B(i,j)*cos(e(i)-e(j)));
N(i,j)=-U(i)*U(j)*(G(i,j)*cos(e(i)-e(j))+B(i,j)*sin(e(i)-e(j)));
L(i,j)=H(i,j); M(i,j)=-N(i,j);
end
end
end
H,N,M,L
%当i=j 时H,N,M,L如下:
for i=1:4
for j=1:5
if i~=j
H(i,i)=H(i,i)+U(i)*U(j)*(G(i,j)*sin(e(i)-e(j))-B(i, j)*cos (e(i)-e(j)));
N(i,i)=N(i,i)-U(i)*U(j)*(G(i, j)*cos(e(i)-e(j))+B(i,j)*sin(e(i)-e(j))); M(i,i)=M(i,i)-U(i)*U(j)*(G(i,j)*cos(e(i)-e(j))+B(i,j)*sin(e(i)-e(j))); L(i,i)=L(i,i)-U(i)*U(j)*(G(i,j)*sin(e(i)-e(j))-B(i,j)*cos(e(i)-e(j)));
end
end
N(i,i)=N(i,i)-2*(U(i))^2*G(i,i);
L(i,i)=L(i,i)+2*(U(i))^2*B(i,i);
end
J=[H,N;M,L]
ox=-((inv(J))*fx);
for i=1:4
oe(i)=ox(i);
oU(i)=ox(i+4)*U(i);
end
for i=1:4
e(i)=e(i)+oe(i);
U(i)=U(i)+oU(i);
end
count=count+1;
end
ox,U,e,count
for j=1:5
P(i)=U(i)*U(j)*(G(i,j)*cos(e(i)-e(j))+B(i,j)*sin(e(i)-e(j)))+P(i);
Q(i)=U(i)*U(j)*(G(i,j)*sin(e(i)-e(j))-B(i,j)*cos(e(i)-e(j)))+Q(i);
end
S(5)=P(5)+Q(5)*sqrt(-1);
S
I=Y*U'。
全部