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

jacobi(2*i,2*j-1)=B(i,j)*e(i)-G(i,j)*f(i); jacobi(2*i,2*j)=G(i,j)*e(i)+B(i,j)*f(i); end end

for i=1:2:7 jacobi(i,i)=0; jacobi(i,i+1)=0; for j=1:5

jacobi(i,i)=jacobi(i,i)-(G((i+1)/2,j)*e(j)-B((i+1)/2,j)*f(j)); jacobi(i,i+1)=jacobi(i,i+1)-(G((i+1)/2,j)*f(j)+B((i+1)/2,j)*e(j)); end

jacobi(i,i)=jacobi(i,i)-G((i+1)/2,(i+1)/2)*e((i+1)/2)-B((i+1)/2,(i+1)/2)*f((i+1)/2); jacobi(i,i+1)=jacobi(i,i+1)+B((i+1)/2,(i+1)/2)*e((i+1)/2)-G((i+1)/2,(i+1)/2)*f((i+1)/2); end

for i=2:2:6

jacobi(i,i-1)=0; jacobi(i,i)=0; for j=1:5

jacobi(i,i-1)=jacobi(i,i-1)-(G(i/2,j)*f(j)-B(i/2,j)*e(j)); jacobi(i,i)=jacobi(i,i)-(G(i/2,j)*e(j)-B(i/2,j)*f(j)); end

jacobi(i,i-1)=jacobi(i,i-1)+B(i/2,i/2)*e(i/2)-G(i/2,i/2)*f(i/2); jacobi(i,i)=jacobi(i,i)+G(i/2,i/2)*e(i/2)+B(i/2,i/2)*f(i/2); end

jacobi(8,7)=-2*e(4); jacobi(8,8)=-2*f(4); while (unlineeps>eps1) deta_pp=zeros(4,1); for i=1:4 for j=1:5

deta_pp(i)=deta_pp(i)+deta_e0(i)*(G(i,j)*deta_e0(j)-B(i,j)*deta_f0(j))+deta_f0(i)*(G(i,j)*deta_f0(j)+B(i,j)*deta_e0(j));

end end

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

deta_qq(i)=deta_qq(i)+deta_f0(i)*(G(i,j)*deta_e0(j)-B(i,j)*deta_f0(j))-deta_e0(i)*(G(i,j)*deta_f0(j)+B(i,j)*deta_e0(j));

end end

deta_UUU=-(deta_e0(4)*deta_e0(4)+deta_f0(4)*deta_f0(4));