哈工大电力系统分析短路及潮流计算实验上机程序

2*j,1);p=p+g*e-b*f;q=q+g*f+b*e;endW(2*k,x)=C(k,2)-Q(2*k-1,1)*p-Q(2*k)*q;W(2*k-1,x)=(C(k,3))^2-(Q(2*k-1,1)^2+Q(2*k)^2);End%%%================================================================%%%===============================================================%比较是否符合条件Max=0;fori=1:1:2*(A(1,1)-1)Max=max(abs(W(i,x)),Max);endifMax<A(5,1)break;end%%%=================================================================%%%=================================================================%求雅克比行列式fori=1:1:A(1,1)-1%===========================================forj=1:1:A(1,1)-1%===========================i~=jifi~=jk=C(i,1);g=real(Y(k,j));b=imag(Y(k,j));e=Q(2*k-1,1);f=Q(2*k,1);J(2*k-1,2*j-1)=-(g*e+b*f);J(2*k-1,2*j)=b*e-g*f;ifi~=3J(2*k,2*j)=-J(2*k-1,2*j-1);J(2*k,2*j-1)=J(2*k-1,2*j);elseJ(2*k,2*j-1)=0;J(2*k,2*j)=0;