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

yb(5,5)=(-33.33333i)/2; yb=yb+conj(yb'); k=0;

eps1=10^-4; unlineeps=1; G=real(yb); B=imag(yb); U=zeros(5,1);

e=[1; 1 ;1 ;1.05 ;1.05]; f=zeros(5,1);

pis=[-1.6; -2 ;-3.7; 5]; qis=[-0.8; -1; -1.3]; deta_p=zeros(4,1); deta_q=zeros(3,1); deta=zeros(8,1); deta_ef=zeros(8,1); deta_e0=zeros(5,1); deta_e1=zeros(5,1); deta_f0=zeros(5,1); deta_f1=zeros(5,1); p=zeros(4,1); q=zeros(3,1); for i=1:4 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

deta_p(i)=pis(i)-p(i); end for i=1:3 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

deta_q(i)=qis(i)-q(i); end

deta_UU=1.05*1.05-(e(4)*e(4)+f(4)*f(4)); jacobi=zeros(8,8); for i=1:4 for j=1:4

jacobi(2*i-1,2*j-1)=-G(i,j)*e(i)-B(i,j)*f(i); jacobi(2*i-1,2*j)=B(i,j)*e(i)-G(i,j)*f(i); end end for i=1:3 for j=1:4