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

zhyerica

[讨论] 求助~~PQ法B1矩阵的形成
9 @! |8 d! `* k
求助, PQ分解法
, ]" t% N- h( |7 ]9 v! J各位大牛，求助啊~~
我毕设的题目是交直流混联系统的潮流计算，自己编了交替迭代法的的PQ分解法程序，计算结果与书上算例结果误差最大有12%。
于是测试了纯交流部分的程序，也就是一般的PQ法，对照着现代电力系统分析P86上的例子，发现自己程序B1矩阵算到的结果跟书上的有差，致使第一次迭代时相角的修正量跟标准值就差了4%左右。但是查到我程序生成的B1矩阵跟自己手算的结果是一样的啊……
* n$ L( }, [; T0 J- m* B
8 `( o* P7 _) O/ Y附上相关部分代码 Ps 我得到的节点导纳矩阵和B2矩阵都是和书上的一样的，真心搞不懂了，求解答啊~~~
0 \! X( A/ W3 W$ x! Q% n

2 p% W/ Z, ^$ C) L7 Xfunction [Bp, Bpp] = makeB(Node,Branch)

%% define named indices into node, branch matrices
9 }7 {  ~( U7 z5 W$ l[PQ, PV, REF, ...
  c8 w. k( k  J5 G    TRANS, CAPAC, LINE,...
% i. P  T' V& a1 Q: @: s$ F: x( R    NODE_I, BUS_TYPE,PG,QG, PL, QL, V, THETA, ...
    BRANCHNUM, NODE1, NODE2, R, X, BORK, BRANCH_TYPE, STATUS, P1, Q1, P2, Q2] = index;
" E; {! l8 E9 ~3 ^" w0 ]* w! _
/ k' b$ r/ C# b9 q5 C/ s8 |%% 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
. @) ~+ K% Q4 a4 Lpq  = find(Node(:, BUS_TYPE) == PQ);    %% PQ bus indices

%*************************************************************************%
%                              form Bp
%*************************************************************************%

%找出非接地支路
temp = find( Branch(:,NODE2)~=0 );
NongroundingBranch = Branch(temp,:);
9 s$ q9 j/ `# B
4 o; t2 c% A1 V* P& k%忽略线路充电电容
  T; I0 Q- }1 Q) g3 UNongroundingBranch(:,BORK) = 0;
2 W6 b6 |: r, |6 n
%不考虑变压器非标准变比
temp = find( NongroundingBranch(:,BRANCH_TYPE)==TRANS );
, X+ l. D" d: X8 P4 A* ?8 `+ ^) c# bNongroundingBranch(temp,BORK) = 1;
Y = makeY(Node,NongroundingBranch);
; P- z& X2 `2 x" O5 J' C0 GBp = imag( Y );
( ]( W- g  g3 M% \- pBp = Bp([pq; pv], [pq; pv])  ;  % reduce Bp matrices
" m8 I5 G/ ^; m1 h
%*************************************************************************%
%                              form Bpp
%*************************************************************************%
' X% S( R5 N6 B( TBranch(:,4)=0;%忽略输电线路电阻
3 J  M$ V7 `" j. z' `Y = makeY(Node,Branch);
Bpp = imag( Y );
Bpp = Bpp(pq, pq);    % reduce Bp matrices

loverqi

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