极坐标下的雅可比矩阵形成

329698914

for m=1:Npq
        for n=1:Nbus
            pt(n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n))); %由节点电压求得的PQ节点注入有功功率
7 [4 i  i6 Q! b) O1 }            qt(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n))); %由节点电压求得的PQ节点注入无功功率
        end
& P$ Q; s* j+ ~. x2 y5 N, K! v        Unbalance(2*m-1)=p(m)-sum(pt); %计算PQ节点有功功率不平衡量
        Unbalance(2*m  )=q(m)-sum(qt); %计算PQ节点无功功率不平衡量
end %[Unbalance]是节点不平衡量矩阵
. \- E3 X* l$ `: q) w5 Y6 q
0 K& ?2 o7 L( c9 Z; j8 j    for m=1:Npq
- I$ ~9 _8 Q% B" P        for n=1:Nbus
/ `! Y, j' q# X5 s            h0(n)= u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));
            n0(n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));
% l* z9 E" r3 a5 N1 K            j0(n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));
9 Z- E$ W' @6 `7 W; n- ]            l0(n)=-u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));
# w: V7 W, H6 K0 @        end
9 R5 D/ |6 `* A/ j- f; `) Z        H(m,m)=sum(h0)-u(m)^2*(G(m,m)*sin(delt(m)-delt(m))-B(m,m)*cos(delt(m)-delt(m)));
        N(m,m)=sum(n0)+u(m)^2*(G(m,m)*cos(delt(m)-delt(m))+B(m,m)*sin(delt(m)-delt(m)))-2*u(m)^2*G(m,m);
        J(m,m)=sum(j0)+u(m)^2*(G(m,m)*cos(delt(m)-delt(m))+B(m,m)*sin(delt(m)-delt(m)));
        L(m,m)=sum(l0)+u(m)^2*(G(m,m)*sin(delt(m)-delt(m))-B(m,m)*cos(delt(m)-delt(m)))+2*u(m)^2*B(m,m);
% A/ T# d! c# m% [( n" D        Jacobi(2*m-1,2*m-1)=H(m,m);
        Jacobi(2*m-1,2*m  )=N(m,m);
        Jacobi(2*m  ,2*m-1)=J(m,m);
) P" t4 V: n' U! _        Jacobi(2*m  ,2*m  )=L(m,m);
    end %计算m=n情况下的Jacobi矩阵中的子矩阵元素
: O/ R! ]: i& q) I: P  ^) N+ ~; O   
* q7 l5 ~9 V4 Z2 |( u& E4 e9 n7 T  {    for m=1:Npq
- d# ]& R0 B, d. I: ~% o) @, r/ ?& R        for n=1:Npq
            if m==n
            else
                H(m,n)=-u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));
5 U) _) X. o5 I5 k" Q- C5 i9 {4 C. T                J(m,n)= u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));
                N(m,n)=-J(m,n);
& Z# T. F! t7 B' ?9 |                L(m,n)= H(m,n);
2 Z% d# Y9 |7 J# t3 F                Jacobi(2*m-1,2*n-1)=H(m,n);
* r9 D- D9 f# W4 q& q                Jacobi(2*m-1,2*n  )=N(m,n);
3 B: \: w, y3 q" L, j/ x9 M# h7 L                Jacobi(2*m  ,2*n-1)=J(m,n);
2 v/ l' S. V  M4 M( H                Jacobi(2*m  ,2*n  )=L(m,n);
: G( W. w9 O0 }/ P- _( p+ M) N, ?1 P4 ]) u            end
        end
    end %计算m≠n情况下的Jacobi矩阵中的子矩阵元素

! r. A9 _$ I/ H问题：１、在算对角元素时候，为什么不能用前面pt(n)和qt(n)代呢
　　　２、对角元素H(m,m),N(m,m),J(m,m),L(m,m)中还有u(m)^2*(G(m,m)*sin(delt(m)-delt(m))-B(m,m)*cos(delt(m)-delt(m)));等，书上不是就u(m)^2*G(m,m)＋Ｑ（i）?
( a4 e' s- h' B4 D觉得很疑问？谢谢大家回答我！