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

disp(U);

theta=zeros(5,1); for i=1:5

theta(i)=angle(U(i)); end

disp(theta); p=zeros(5,1); q=zeros(5,1); for i=1:5 for j=1:5

p(i)=p(i)+e(i)*(G(i,j)*e(j)-B(i,j)*f(j))+f(i)*(G(i,j)*f(j)+B(i,j)*e(j)); end end for i=1:5 for j=1:5

q(i)=q(i)+f(i)*(G(i,j)*e(j)-B(i,j)*f(j))-e(i)*(G(i,j)*f(j)+B(i,j)*e(j)); end end disp(p); 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 disp(s); disp(deta_s);