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

else

kq=1; end

if max(eps_theta,eps_U)>eps1 k=k+1; end end for i=1:5

U(i)=U(i)*(cos(theta(i))+sqrt(-1)*sin(theta(i))); end disp(U); disp(theta); p=zeros(5,1); q=zeros(5,1); for i=1:5 for j=1:5

p(i)=p(i)+abs(U(i))*abs(U(j))*(G(i,j)*cos(theta(i)-theta(j))+B(i,j)*sin(theta(i)-theta(j))); end end disp(p); for i=1:5 for j=1:5

q(i)=q(i)+abs(U(i))*abs(U(j))*(G(i,j)*sin(theta(i)-theta(j))-B(i,j)*cos(theta(i)-theta(j))); end end disp(q);

y0=[-4i;-8i;-4i;0;0]; s=zeros(5,5); for i=1:5 for j=1:5

s(i,j)=U(i)*(conj(U(i))*conj(y0(i))+(conj(U(i))-conj(U(j)))*conj(-yb(i,j))); end end for i=1:5 s(i,i)=0; end

deta_s=zeros(5,5); for i=1:5 for j=1:5

deta_s(i,j)=s(i,j)+s(j,i); end end