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

leegengda

PQ分解潮流程序中因子表的形成
程序代码如下：
节点导纳矩阵Y已经形成，BanNode是平衡节点，PVi是PV节点
) V, X& |) B) J% k
%形成因子表
( {* c3 l: q. w; E0 O. V/ ?function [ B B1 B2 ]=form_B(Y,BanNode,PVi)

B=imag(Y);
% ?& Z) m2 H4 FB1=B;
- k( y( U& ~& fB1(:,BanNode)=[];
B1(BanNode,:)=[];
& ?  r1 ]: m4 e! s
" L  k/ @2 I; O' L  HB2=B1;
B2(:,PVi)=[];
3 G; k+ n+ z" ]- G) xB2(PVi,:)=[];


6 [- t7 E. \; L! F9 G" z% }5 lend
$ p9 T7 [, l' I9 e* h
9 ~3 ~; B0 ~: S* Y7 q1 }请高手指教，这样形成的因子表有啥问题？