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

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节点注入有功功率
            qt(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n))); %由节点电压求得的PQ节点注入无功功率
        end
        Unbalance(2*m-1)=p(m)-sum(pt); %计算PQ节点有功功率不平衡量
" U5 Y% `% I" W8 k        Unbalance(2*m  )=q(m)-sum(qt); %计算PQ节点无功功率不平衡量
end %[Unbalance]是节点不平衡量矩阵
! I+ l! B5 T, B' @8 j. f
9 y5 S2 d, \) v( s, P9 |    for m=1:Npq
        for n=1:Nbus
! X" B2 Y9 F# l3 m, e% F' ^5 H* q  W            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)));
            j0(n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));
! B& k  Y+ l* V' K+ N            l0(n)=-u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));
4 y" o7 l6 g. r        end
5 d7 ^% ~7 T# U8 r4 L5 L        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);
# v& K8 r6 H8 r        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);
        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);
2 H( W/ n0 v2 u( E        Jacobi(2*m  ,2*m  )=L(m,m);
8 H% Q! f& x% [! V7 ^, r    end %计算m=n情况下的Jacobi矩阵中的子矩阵元素
( l) l6 W* R. s2 D+ A: y- G   
+ x9 M$ U6 j* f& N5 d    for m=1:Npq
        for n=1:Npq
2 w8 {5 E. h3 S3 B2 y            if m==n
            else
9 {- O4 F# u. b. f* s                H(m,n)=-u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));
                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);
/ y" Q0 m6 h% w% N                L(m,n)= H(m,n);
                Jacobi(2*m-1,2*n-1)=H(m,n);
                Jacobi(2*m-1,2*n  )=N(m,n);
                Jacobi(2*m  ,2*n-1)=J(m,n);
                Jacobi(2*m  ,2*n  )=L(m,n);
& s7 S6 G+ B2 N  _) @            end
        end
9 H3 G0 ]; @1 c. N+ N2 ~    end %计算m≠n情况下的Jacobi矩阵中的子矩阵元素
; `3 S+ {) X' R' n, {' C# @
问题：１、在算对角元素时候，为什么不能用前面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）?
觉得很疑问？谢谢大家回答我！