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

deta_U=zeros(3,1); pis=[-1.6; -2 ;-3.7; 5]; qis=[-0.8; -1; -1.3]; deta_p=zeros(4,1); deta_q=zeros(3,1); while (kp==0||kq==0) p=zeros(4,1); for i=1:4 for j=1:5

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

deta_p(i)=pis(i)-p(i); end for i=1:4

deta_p(i)=deta_p(i)/U(i); end

deta_theta=-inv(B1)*deta_p; for i=1:4

deta_theta(i)=deta_theta(i)/U(i); theta(i)=theta(i)+deta_theta(i); end

eps_theta=max(max(abs(deta_theta)),max(abs(deta_p))); if eps_theta<eps1 kp=1; else

kp=0; end

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

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

deta_q(i)=qis(i)-q(i); end for i=1:3

deta_q(i)=deta_q(i)/U(i); end

deta_U=-inv(B2)*deta_q; for i=1:3

U(i)=U(i)+deta_U(i); end

eps_U=max(max(abs(deta_U)),max(abs(deta_q))); if eps_U<eps1 kq=1;