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

zhyerica

[讨论] 求助~~PQ法B1矩阵的形成
$ z& M! h/ h( t5 I+ @$ O
0 f6 A- B+ ?3 L7 |2 X2 x* p8 W求助, PQ分解法
各位大牛，求助啊~~
我毕设的题目是交直流混联系统的潮流计算，自己编了交替迭代法的的PQ分解法程序，计算结果与书上算例结果误差最大有12%。
于是测试了纯交流部分的程序，也就是一般的PQ法，对照着现代电力系统分析P86上的例子，发现自己程序B1矩阵算到的结果跟书上的有差，致使第一次迭代时相角的修正量跟标准值就差了4%左右。但是查到我程序生成的B1矩阵跟自己手算的结果是一样的啊……
/ i9 }6 P1 w+ {) `. [/ @
; J" C, q! g; v3 B" y, }* f' f附上相关部分代码 Ps 我得到的节点导纳矩阵和B2矩阵都是和书上的一样的，真心搞不懂了，求解答啊~~~

  i( {9 N7 N5 ]  p9 ^0 c
function [Bp, Bpp] = makeB(Node,Branch)
! e' C$ Y: S/ ~
$ `  z' N/ N% c, @4 I  e# I%% define named indices into node, branch matrices
) m3 `+ n4 H$ e( F5 @, K; }2 T[PQ, PV, REF, ...
- _$ U7 x3 W+ D! t7 Y" z7 E" K: Q8 E& j    TRANS, CAPAC, LINE,...
  E$ x6 l. o2 G. o1 g" ?    NODE_I, BUS_TYPE,PG,QG, PL, QL, V, THETA, ...
    BRANCHNUM, NODE1, NODE2, R, X, BORK, BRANCH_TYPE, STATUS, P1, Q1, P2, Q2] = index;
: p) [0 i, s! m# L6 ]( Y5 D  k
6 b# I& b3 O$ n6 {%% 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
) N5 k. z5 g: M# L8 I  [- a: _. c( Gpq  = find(Node(:, BUS_TYPE) == PQ);    %% PQ bus indices
3 t' }5 n0 r# O) l
( i3 {  f3 T5 I5 e%*************************************************************************%
! h0 A9 F# S( f6 A%                              form Bp
8 a( e& x' m' A% v2 d9 a%*************************************************************************%
* S+ a2 p1 ?* h+ m+ b# _+ K5 R# y
%找出非接地支路
temp = find( Branch(:,NODE2)~=0 );
* S/ Q- K+ H/ H" B% k- WNongroundingBranch = Branch(temp,:);
6 L; C3 [2 h/ ]. z. k" z% r' X+ N
%忽略线路充电电容
NongroundingBranch(:,BORK) = 0;
/ w; \, V& V% j0 E  X/ f2 R
%不考虑变压器非标准变比
temp = find( NongroundingBranch(:,BRANCH_TYPE)==TRANS );
NongroundingBranch(temp,BORK) = 1;
Y = makeY(Node,NongroundingBranch);
  ]2 p# l& Z6 f; W& bBp = imag( Y );
Bp = Bp([pq; pv], [pq; pv])  ;  % reduce Bp matrices
! I; g( i$ a  W2 z2 C3 O- @" A  m/ y
%*************************************************************************%
# Z4 m" Z& u% L- ?' `" u8 A. _%                              form Bpp
%*************************************************************************%
Branch(:,4)=0;%忽略输电线路电阻
Y = makeY(Node,Branch);
6 u' L1 _2 p- T3 \* fBpp = imag( Y );
8 u; F' \, i: p7 ]) l% [4 c5 o6 gBpp = Bpp(pq, pq);    % reduce Bp matrices

loverqi

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