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

zhyerica

[讨论] 求助~~PQ法B1矩阵的形成
* d/ _  y) J8 {4 r! F7 y
求助, PQ分解法
, N' a5 [' X& e各位大牛，求助啊~~
6 l/ a0 y( I' k2 F0 c  Y" |我毕设的题目是交直流混联系统的潮流计算，自己编了交替迭代法的的PQ分解法程序，计算结果与书上算例结果误差最大有12%。
1 {* I; J) A& s& h8 A于是测试了纯交流部分的程序，也就是一般的PQ法，对照着现代电力系统分析P86上的例子，发现自己程序B1矩阵算到的结果跟书上的有差，致使第一次迭代时相角的修正量跟标准值就差了4%左右。但是查到我程序生成的B1矩阵跟自己手算的结果是一样的啊……
( S2 ?; p) \- B: W& E6 N
# T) L# e4 w/ u% ^, A6 o$ s  N3 ~附上相关部分代码 Ps 我得到的节点导纳矩阵和B2矩阵都是和书上的一样的，真心搞不懂了，求解答啊~~~

/ v' g. I" @! v
function [Bp, Bpp] = makeB(Node,Branch)

%% define named indices into node, branch matrices
- [( ^  j, z' _3 d7 J) y[PQ, PV, REF, ...
    TRANS, CAPAC, LINE,...
" u' e! }  E0 e, |% }+ F/ y    NODE_I, BUS_TYPE,PG,QG, PL, QL, V, THETA, ...
" a" q0 n2 j) q6 X    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
/ h; R, j; B/ l* h; I% R5 Tpq  = find(Node(:, BUS_TYPE) == PQ);    %% PQ bus indices
- K1 |( ^3 n/ `
%*************************************************************************%
/ G* ]. T1 K; Y1 k( D+ v%                              form Bp
%*************************************************************************%

%找出非接地支路
$ C* N+ \6 u, R4 _7 R7 Ctemp = find( Branch(:,NODE2)~=0 );
NongroundingBranch = Branch(temp,:);
7 {1 {0 _% C7 F' E, Z# J
0 ^+ A5 `$ }4 r2 D%忽略线路充电电容
& w# R4 E# j  z, S9 Q0 N0 ]! @NongroundingBranch(:,BORK) = 0;

8 K' }; `+ d& K. H; E: P( t% i%不考虑变压器非标准变比
3 r: W  n2 j' y4 W5 s- E1 ~temp = find( NongroundingBranch(:,BRANCH_TYPE)==TRANS );
" }4 ]) M9 o# zNongroundingBranch(temp,BORK) = 1;
Y = makeY(Node,NongroundingBranch);
- l- k+ e! q' |3 O8 EBp = imag( Y );
* s( R6 ^7 |+ I$ Y6 s. B3 UBp = Bp([pq; pv], [pq; pv])  ;  % reduce Bp matrices
- d1 j+ H4 C- t& `
( l' x9 c* U8 b9 F/ u; S%*************************************************************************%
; G, p" X. ^- P$ ~1 X%                              form Bpp
%*************************************************************************%
Branch(:,4)=0;%忽略输电线路电阻
% T$ q0 x3 z" M- _0 IY = makeY(Node,Branch);
Bpp = imag( Y );
% g. k  L% g0 a9 ]- f1 I: h' BBpp = Bpp(pq, pq);    % reduce Bp matrices

loverqi

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