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

zhyerica

[讨论] 求助~~PQ法B1矩阵的形成

$ c/ u0 E: p/ H, T0 l( H& L; q求助, PQ分解法
各位大牛，求助啊~~
0 q: K0 W6 t1 m& G* ~我毕设的题目是交直流混联系统的潮流计算，自己编了交替迭代法的的PQ分解法程序，计算结果与书上算例结果误差最大有12%。
) j% S! k- B" q6 \于是测试了纯交流部分的程序，也就是一般的PQ法，对照着现代电力系统分析P86上的例子，发现自己程序B1矩阵算到的结果跟书上的有差，致使第一次迭代时相角的修正量跟标准值就差了4%左右。但是查到我程序生成的B1矩阵跟自己手算的结果是一样的啊……

8 T  w  U4 p) C9 P5 `附上相关部分代码 Ps 我得到的节点导纳矩阵和B2矩阵都是和书上的一样的，真心搞不懂了，求解答啊~~~

: F6 e* u# C- Y& _2 x& y/ T& N0 }
function [Bp, Bpp] = makeB(Node,Branch)

4 b6 Y! F- \8 j; T- S% d$ a4 q%% define named indices into node, branch matrices
[PQ, PV, REF, ...
    TRANS, CAPAC, LINE,...
    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
& c+ p+ ]- s! ~; ?) @8 ^1 p) {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
' y6 }: K3 c; |  l* s
# T) l  b6 Z/ I0 l' R%*************************************************************************%
8 {* F; m! Z/ u& L2 k. K. Y%                              form Bp
3 r$ _- p' Y" ]. Y& z3 v$ ?%*************************************************************************%

3 {0 Q) g* F. ]( s3 p5 N0 ~%找出非接地支路
temp = find( Branch(:,NODE2)~=0 );
  |: k- n" ^4 dNongroundingBranch = Branch(temp,:);

%忽略线路充电电容
NongroundingBranch(:,BORK) = 0;
  E* u. \# n1 z0 p
7 e) n: W+ ^: `1 ^( w2 ~/ H. m%不考虑变压器非标准变比
0 b& D8 c! s0 u9 |: utemp = find( NongroundingBranch(:,BRANCH_TYPE)==TRANS );
+ C! o7 B8 G( k. NNongroundingBranch(temp,BORK) = 1;
# Z, x  O. R; }, lY = makeY(Node,NongroundingBranch);
7 J; k  t3 A/ f5 l, b3 l* @" {Bp = imag( Y );
& M8 t8 G  v& Z! @! j7 {8 u4 KBp = Bp([pq; pv], [pq; pv])  ;  % reduce Bp matrices

; _, h: s+ p. w5 d; Z5 p%*************************************************************************%
%                              form Bpp
%*************************************************************************%
2 x/ c1 |2 Y" V: g. UBranch(:,4)=0;%忽略输电线路电阻
Y = makeY(Node,Branch);
  {" P- F& z- vBpp = imag( Y );
Bpp = Bpp(pq, pq);    % reduce Bp matrices

loverqi

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