电力系统稳态实验报告(12)

disp(s); disp(deta_s);

3.高斯赛德尔法

clear; clc;

yb=zeros(5,5);

yb(1,1)=(1.37874-6.26166i)/2;yb(1,2)=-0.62402+3.90015i;yb(1,3)=-0.75471+2.64150i; yb(2,2)=(1.45390-66.98082i)/2;yb(2,3)=-0.82987+3.11203i;yb(2,4)=63.49206i; yb(3,3)=(1.58459-35.73786i)/2;yb(3,5)=31.74603i; yb(4,4)=(-66.66667i)/2; yb(5,5)=(-33.33333i)/2; yb=yb+conj(yb'); G=real(yb); B=imag(yb);

s=[-1.6-0.8i;-2-i;-3.7-1.3i;5]; U0=[1;1;1;1.05;1.05]; U1=zeros(4,1); deta_U=zeros(4,1); eps1=10^-4; eps_U=1; theta=0; k=0;

while (eps_U>eps1) q=0; for j=1:5

q=q+U0(4)*conj(yb(4,j))*conj(U0(j)); end

q=imag(q);

s(4)=5+sqrt(-1)*q; p=0; for j=1:5

p=p+yb(4,j)*U0(j); end

p=p-yb(4,4)*U0(4);

U1(4)=(conj(s(4))/conj(U0(4))-p)/yb(4,4); theta=angle(U1(4));

U1(4)=1.05*(cos(theta)+sqrt(-1)*sin(theta)); for i=1:3 a=0; b=0; for j=1:i-1

a=a+yb(i,j)*U1(j); end

for j=i+1:4