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

zhyerica

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

求助, PQ分解法
各位大牛，求助啊~~
  b& {$ {0 F5 `' K" }5 b  M! ^" p我毕设的题目是交直流混联系统的潮流计算，自己编了交替迭代法的的PQ分解法程序，计算结果与书上算例结果误差最大有12%。
于是测试了纯交流部分的程序，也就是一般的PQ法，对照着现代电力系统分析P86上的例子，发现自己程序B1矩阵算到的结果跟书上的有差，致使第一次迭代时相角的修正量跟标准值就差了4%左右。但是查到我程序生成的B1矩阵跟自己手算的结果是一样的啊……
: C2 k1 W, c) I! ?% a
附上相关部分代码 Ps 我得到的节点导纳矩阵和B2矩阵都是和书上的一样的，真心搞不懂了，求解答啊~~~


function [Bp, Bpp] = makeB(Node,Branch)
$ C- R4 @2 O5 ^' a4 g; H9 e
%% 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;
% L2 ]7 s) T, I( c/ }% r' I
) C/ T% g3 ^+ ?%% form index lists for slack, PV, and PQ buses
ref = find(Node(:, 2) == 3);   %% reference bus index
& J8 N$ ^4 ?: Q  N' A$ Z8 fpv  = find(Node(:, BUS_TYPE) == PV);    %% PV bus indices
pq  = find(Node(:, BUS_TYPE) == PQ);    %% PQ bus indices

%*************************************************************************%
1 Y2 T' Q4 u! S# R3 z% X, p5 ]+ C%                              form Bp
%*************************************************************************%

  ^. b* \! D4 H: c! Q, d+ r%找出非接地支路
temp = find( Branch(:,NODE2)~=0 );
NongroundingBranch = Branch(temp,:);
4 y/ h4 B, X: @! |/ S6 L5 R
- v9 t- c' Q# ^  W%忽略线路充电电容
NongroundingBranch(:,BORK) = 0;

% t- |0 a( w/ p%不考虑变压器非标准变比
" a- B* j# T3 |* c' Ctemp = find( NongroundingBranch(:,BRANCH_TYPE)==TRANS );
  N  b, N! p" _# cNongroundingBranch(temp,BORK) = 1;
: y- I% ~" H" y9 dY = makeY(Node,NongroundingBranch);
9 t% ?; _/ [9 O0 M/ J, DBp = imag( Y );
' m+ m+ U  f# R( K- N5 Y& K% hBp = Bp([pq; pv], [pq; pv])  ;  % reduce Bp matrices

%*************************************************************************%
$ b" m" E; J  t%                              form Bpp
% Z% q, O' i# k' h+ L%*************************************************************************%
Branch(:,4)=0;%忽略输电线路电阻
$ h: f5 @7 G* y& A' _& }Y = makeY(Node,Branch);
* r. O  U9 S/ m; V' ^Bpp = imag( Y );
) ?% K) r$ P# n" d' JBpp = Bpp(pq, pq);    % reduce Bp matrices

loverqi

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