重提问 PQ分解法B1矩阵该如何求啊

zhyerica

[讨论] 求助~~PQ法B1矩阵的形成
6 y5 H9 H. J7 c8 \4 k
( M, Y$ I5 d' ]" @* `求助, PQ分解法
各位大牛，求助啊~~
" D* c7 ?- }0 Q1 r) Y  f我毕设的题目是交直流混联系统的潮流计算，自己编了交替迭代法的的PQ分解法程序，计算结果与书上算例结果误差最大有12%。
' R$ X) i! |8 m) O; l1 X( s  ~于是测试了纯交流部分的程序，也就是一般的PQ法，对照着现代电力系统分析P86上的例子，发现自己程序B1矩阵算到的结果跟书上的有差，致使第一次迭代时相角的修正量跟标准值就差了4%左右。但是查到我程序生成的B1矩阵跟自己手算的结果是一样的啊……

附上相关部分代码 Ps 我得到的节点导纳矩阵和B2矩阵都是和书上的一样的，真心搞不懂了，求解答啊~~~
( }" j# L9 M+ a5 J: ]7 [( v1 \

function [Bp, Bpp] = makeB(Node,Branch)
1 E, z, p: I2 C# g' Z
%% define named indices into node, branch matrices
1 m  D% z  Y. b[PQ, PV, REF, ...
    TRANS, CAPAC, LINE,...
& O; k4 H& f& i8 q" `    NODE_I, BUS_TYPE,PG,QG, PL, QL, V, THETA, ...
    BRANCHNUM, NODE1, NODE2, R, X, BORK, BRANCH_TYPE, STATUS, P1, Q1, P2, Q2] = index;

5 ]/ G/ j% p4 R& O( Y6 f%% form index lists for slack, PV, and PQ buses
ref = find(Node(:, 2) == 3);   %% reference bus index
9 f! F* ?4 E7 }; R# z3 Y3 `) fpv  = find(Node(:, BUS_TYPE) == PV);    %% PV bus indices
pq  = find(Node(:, BUS_TYPE) == PQ);    %% PQ bus indices
8 z1 i; }! o* @
# m# B. D( e* ^7 {%*************************************************************************%
  C8 E! H+ V, L5 _2 @%                              form Bp
%*************************************************************************%

%找出非接地支路
5 r& c- Z& R) S- T0 o! Itemp = find( Branch(:,NODE2)~=0 );
NongroundingBranch = Branch(temp,:);

: y: A+ C0 l6 z/ N0 x%忽略线路充电电容
$ K: \% }4 o/ wNongroundingBranch(:,BORK) = 0;
6 w; o' x1 L+ b5 W: u2 U7 e: T5 \& o
- |3 G1 x- K8 ?( S+ ^/ O) y- i%不考虑变压器非标准变比
" ~' ~8 Q) W- l+ g7 z; Xtemp = find( NongroundingBranch(:,BRANCH_TYPE)==TRANS );
NongroundingBranch(temp,BORK) = 1;
Y = makeY(Node,NongroundingBranch);
" b" f, n$ }2 |& TBp = imag( Y );
Bp = Bp([pq; pv], [pq; pv])  ;  % reduce Bp matrices
& e( A# n# A7 n( F7 J1 v6 T, o
%*************************************************************************%
%                              form Bpp
8 v: k1 O) u% H5 l5 Y$ ~%*************************************************************************%
Branch(:,4)=0;%忽略输电线路电阻
0 Y' h+ V- I. m3 d: ZY = makeY(Node,Branch);
. D' o/ f  r0 `" D3 F- QBpp = imag( Y );
Bpp = Bpp(pq, pq);    % reduce Bp matrices

loverqi

这些基础的我都不会，那书我们也学的，。。。要学习