PQ分解潮流程序中因子表的形成

leegengda

PQ分解潮流程序中因子表的形成
程序代码如下：
节点导纳矩阵Y已经形成，BanNode是平衡节点，PVi是PV节点

%形成因子表
function [ B B1 B2 ]=form_B(Y,BanNode,PVi)

B=imag(Y);
) f" [' P9 I  q% c! R$ g7 [B1=B;
B1(:,BanNode)=[];
/ L  _3 A' L+ z9 f" ?. qB1(BanNode,:)=[];
/ p* \0 W2 O' G- i) C; R& Z* I( b
2 ^) u+ x) g( f1 r7 pB2=B1;
9 o& m' G( T$ DB2(:,PVi)=[];
" _2 n# w9 Z. i! g! K" i7 gB2(PVi,:)=[];

- f& T/ X9 e: V* }# v
end
6 [# r! f6 b: `7 q5 i- r; |
请高手指教，这样形成的因子表有啥问题？