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

zhyerica

[讨论] 求助~~PQ法B1矩阵的形成
% A2 H/ Q3 c* U. F) n
求助, PQ分解法
各位大牛，求助啊~~
$ q3 `. {) q8 _7 {我毕设的题目是交直流混联系统的潮流计算，自己编了交替迭代法的的PQ分解法程序，计算结果与书上算例结果误差最大有12%。
于是测试了纯交流部分的程序，也就是一般的PQ法，对照着现代电力系统分析P86上的例子，发现自己程序B1矩阵算到的结果跟书上的有差，致使第一次迭代时相角的修正量跟标准值就差了4%左右。但是查到我程序生成的B1矩阵跟自己手算的结果是一样的啊……

6 d' _6 ^+ {, ~$ h8 k附上相关部分代码 Ps 我得到的节点导纳矩阵和B2矩阵都是和书上的一样的，真心搞不懂了，求解答啊~~~
& m- E. Z0 u$ H4 }3 J

- e" n5 [& s# T$ ~function [Bp, Bpp] = makeB(Node,Branch)
; w$ A3 n0 C+ a
%% define named indices into node, branch matrices
  |+ j& U3 {3 e* V  ^: O; C1 L[PQ, PV, REF, ...
    TRANS, CAPAC, LINE,...
    NODE_I, BUS_TYPE,PG,QG, PL, QL, V, THETA, ...
5 K  R6 R  N7 y. I    BRANCHNUM, NODE1, NODE2, R, X, BORK, BRANCH_TYPE, STATUS, P1, Q1, P2, Q2] = index;
+ M8 ~' ?" x& {4 {9 a! E  A" |1 y
* t+ o# m/ a1 w/ |% }; L- P%% form index lists for slack, PV, and PQ buses
" N8 ?# ^3 S2 t6 `2 Zref = find(Node(:, 2) == 3);   %% reference bus index
9 `9 Q) k# N: r( Spv  = find(Node(:, BUS_TYPE) == PV);    %% PV bus indices
pq  = find(Node(:, BUS_TYPE) == PQ);    %% PQ bus indices

%*************************************************************************%
" _! j. z2 o/ }7 M7 ?%                              form Bp
, z; k7 V8 r! I; c%*************************************************************************%

%找出非接地支路
% w/ h; t; P: }6 Ktemp = find( Branch(:,NODE2)~=0 );
NongroundingBranch = Branch(temp,:);
) j- I" X# \4 ?8 X( Z% j7 g9 D
% m$ D/ {" _2 C' s0 H%忽略线路充电电容
NongroundingBranch(:,BORK) = 0;
7 ^5 I3 c" b* k
%不考虑变压器非标准变比
; @( G# p; ]. O# g6 [& v6 R; Otemp = find( NongroundingBranch(:,BRANCH_TYPE)==TRANS );
NongroundingBranch(temp,BORK) = 1;
Y = makeY(Node,NongroundingBranch);
Bp = imag( Y );
* X) W8 ~$ ^, b* b. v- v+ ^Bp = Bp([pq; pv], [pq; pv])  ;  % reduce Bp matrices

%*************************************************************************%
/ f( u3 b; ]4 m%                              form Bpp
%*************************************************************************%
Branch(:,4)=0;%忽略输电线路电阻
1 E6 B" O* m7 p' X3 pY = makeY(Node,Branch);
Bpp = imag( Y );
. p" b" i, n9 ?$ h0 B0 H  ?Bpp = Bpp(pq, pq);    % reduce Bp matrices

loverqi

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