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

leegengda

PQ分解潮流程序中因子表的形成
5 j5 `+ \# ^! b3 o' ?) A程序代码如下：
节点导纳矩阵Y已经形成，BanNode是平衡节点，PVi是PV节点
/ j( ^1 x6 P4 K  F1 `5 Z! F" a
- J. \. d; ^1 U! E%形成因子表
* r6 _9 ^6 u; D, F( kfunction [ B B1 B2 ]=form_B(Y,BanNode,PVi)
7 B9 }: y/ v- c1 y5 G2 ?9 v5 t0 W
/ R- f: K' v. W/ d% TB=imag(Y);
B1=B;
B1(:,BanNode)=[];
. a7 o$ \& @+ uB1(BanNode,:)=[];
( e/ m! r# F& C: o! e# G
B2=B1;
B2(:,PVi)=[];
* L- e( j$ `# s+ U" \7 VB2(PVi,:)=[];

0 i  M  {" b9 m3 ^0 X" p
end
" A1 j8 m* @# ^
1 v2 L# n0 }; \3 |! I5 w请高手指教，这样形成的因子表有啥问题？