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

zhyerica

[讨论] 求助~~PQ法B1矩阵的形成
+ c" C; x% N+ S; H* I
3 X/ C* Y6 o8 H" |! X$ g/ J+ [0 M$ ^求助, PQ分解法
各位大牛，求助啊~~
# x! g. J' i4 ]( t2 f# e我毕设的题目是交直流混联系统的潮流计算，自己编了交替迭代法的的PQ分解法程序，计算结果与书上算例结果误差最大有12%。
于是测试了纯交流部分的程序，也就是一般的PQ法，对照着现代电力系统分析P86上的例子，发现自己程序B1矩阵算到的结果跟书上的有差，致使第一次迭代时相角的修正量跟标准值就差了4%左右。但是查到我程序生成的B1矩阵跟自己手算的结果是一样的啊……
# {5 [- h: _2 u8 r1 d
附上相关部分代码 Ps 我得到的节点导纳矩阵和B2矩阵都是和书上的一样的，真心搞不懂了，求解答啊~~~
& k1 `1 O0 T# s: @

: a6 _) L" Z6 l. Lfunction [Bp, Bpp] = makeB(Node,Branch)

, e/ h- {1 |' B( x' _) B%% define named indices into node, branch matrices
[PQ, PV, REF, ...
    TRANS, CAPAC, LINE,...
- U* L5 ^; n8 E  X8 R5 A% _4 Q8 n    NODE_I, BUS_TYPE,PG,QG, PL, QL, V, THETA, ...
    BRANCHNUM, NODE1, NODE2, R, X, BORK, BRANCH_TYPE, STATUS, P1, Q1, P2, Q2] = index;

%% form index lists for slack, PV, and PQ buses
ref = find(Node(:, 2) == 3);   %% reference bus index
pv  = find(Node(:, BUS_TYPE) == PV);    %% PV bus indices
pq  = find(Node(:, BUS_TYPE) == PQ);    %% PQ bus indices

" f+ g$ @+ {  d  \: B% `  a%*************************************************************************%
%                              form Bp
%*************************************************************************%

%找出非接地支路
temp = find( Branch(:,NODE2)~=0 );
NongroundingBranch = Branch(temp,:);

3 \! w. A9 C2 a) y" [# W%忽略线路充电电容
9 ~' o% _' y1 N8 ~NongroundingBranch(:,BORK) = 0;
, M* c8 l$ m- G5 B% U1 ~" a( `* A
' ~& K- O+ p" M" K" i2 I%不考虑变压器非标准变比
temp = find( NongroundingBranch(:,BRANCH_TYPE)==TRANS );
/ B. H" U) o5 D( NNongroundingBranch(temp,BORK) = 1;
Y = makeY(Node,NongroundingBranch);
Bp = imag( Y );
6 |, y' A0 O& q. B* Z' q! y* tBp = Bp([pq; pv], [pq; pv])  ;  % reduce Bp matrices

- I- v  B# S1 X%*************************************************************************%
%                              form Bpp
%*************************************************************************%
# ]& `8 a4 P( b5 ^' PBranch(:,4)=0;%忽略输电线路电阻
8 }; C* M7 h( ^: x1 @Y = makeY(Node,Branch);
Bpp = imag( Y );
2 a" `1 f, |1 N- NBpp = Bpp(pq, pq);    % reduce Bp matrices

loverqi

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