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

deta(i)=deta_p((i+1)/2)-deta_pp((i+1)/2); end

for i=2:2:6

deta(i)=deta_q(i/2)-deta_qq(i/2); end

deta(8)=deta_UU+deta_UUU; deta_ef=-inv(jacobi)*deta; for i=1:4

deta_e1(i)=deta_ef(2*i-1); deta_f1(i)=deta_ef(2*i); end

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

deta_ppp(i)=deta_ppp(i)+deta_e1(i)*(G(i,j)*deta_e1(j)-B(i,j)*deta_f1(j))+deta_f1(i)*(G(i,j)*deta_f1(j)+B(i,j)*deta_e1(j));

end end

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

deta_qqq(i)=deta_qqq(i)+deta_f1(i)*(G(i,j)*deta_e1(j)-B(i,j)*deta_f1(j))-deta_e1(i)*(G(i,j)*deta_f1(j)+B(i,j)*deta_e1(j));

end end

unlineeps=max(max(abs(deta_e1-deta_e0)),max(abs(deta_f1-deta_f0))); unlineeps=max(unlineeps,max(abs(deta_qq-deta_qqq))); unlineeps=max(unlineeps,max(abs(deta_pp-deta_ppp))); for i=1:4

deta_e0(i)=deta_e1(i); deta_f0(i)=deta_f1(i); end k=k+1; end for i=1:4

e(i)=e(i)+deta_e0(i); f(i)=f(i)+deta_f0(i); end for i=1:5

U(i)=e(i)+sqrt(-1)*f(i); end