求助~~PQ法B1矩阵的形成

zhyerica

各位大牛，求助啊~~
我毕设的题目是交直流混联系统的潮流计算，自己编了交替迭代法的的PQ分解法程序，计算结果与书上算例结果误差最大有12%。
7 w0 T9 V8 |% f' w# W7 K于是测试了纯交流部分的程序，也就是一般的PQ法，对照着现代电力系统分析P86上的例子，发现自己程序B1矩阵算到的结果跟书上的有差，致使第一次迭代时相角的修正量跟标准值就差了4%左右。但是查到我程序生成的B1矩阵跟自己手算的结果是一样的啊……
9 r. n: ~6 Z5 p! B( ^+ l+ y
3 l" ^& ~  T: R. a6 ~/ `) e7 g, k0 F附上相关部分代码 Ps 我得到的节点导纳矩阵和B2矩阵都是和书上的一样的，真心搞不懂了，求解答啊~~~

" W" @  }4 G6 H8 z& L+ K9 a
2 O* V8 t6 f) Ofunction [Bp, Bpp] = makeB(Node,Branch)

%% define named indices into node, branch matrices
[PQ, PV, REF, ...
  ~% A9 K8 f3 ~) |2 w    TRANS, CAPAC, LINE,...
" o+ j! u2 V  l4 L" l    NODE_I, BUS_TYPE,PG,QG, PL, QL, V, THETA, ...
* K+ c1 A& X$ l8 A6 u0 Q8 ^    BRANCHNUM, NODE1, NODE2, R, X, BORK, BRANCH_TYPE, STATUS, P1, Q1, P2, Q2] = index;
: H3 a* X- j7 Y2 ~- ^
%% 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
6 p. t( X) Z% Q: N) epq  = find(Node(:, BUS_TYPE) == PQ);    %% PQ bus indices
% i7 j. Z6 y* J  @
; _7 n& a# h5 h6 f# {$ I" M%*************************************************************************%
- `6 U4 ]+ N% q/ ?$ E1 n" u% i%                              form Bp
  H# s; j+ @6 P6 O& i/ O%*************************************************************************%
/ m4 i$ h# d7 r" k3 o5 O/ N
9 ^* E  b/ w# {' D5 V%找出非接地支路
temp = find( Branch(:,NODE2)~=0 );
NongroundingBranch = Branch(temp,:);

%忽略线路充电电容
NongroundingBranch(:,BORK) = 0;
* C3 E7 Q, [" k2 Z
7 l5 c6 v# z& i' o1 n%不考虑变压器非标准变比
temp = find( NongroundingBranch(:,BRANCH_TYPE)==TRANS );
4 K5 ?, ]) x6 R) h3 Z4 s" DNongroundingBranch(temp,BORK) = 1;
Y = makeY(Node,NongroundingBranch);
Bp = imag( Y );
Bp = Bp([pq; pv], [pq; pv])  ;  % reduce Bp matrices

" r5 G: x: {& X, b: ?7 M%*************************************************************************%
* d5 O4 r& t/ Q; }/ {) T0 i%                              form Bpp
& x7 O9 y( B, w%*************************************************************************%
Branch(:,4)=0;%忽略输电线路电阻
Y = makeY(Node,Branch);
8 C0 k0 N0 z7 ~Bpp = imag( Y );
Bpp = Bpp(pq, pq);    % reduce Bp matrices
0 q& \" _5 y, R  \5 f4 A6 g
+ J6 h3 `' b" A" K7 A, B+ B6 Yreturn;

zhyerica

自己顶~~~
7 d# h  ?( a! `# \" _7 _  b4 W拿那个5节点的例子，算出来的B1矩阵如下：
. m* M/ T' x5 ?7 t3 Z) G, E7 KB1 =
6 w0 b1 u& h+ ]  W
   (1,1)      -6.5417
   (2,1)       3.9002
! L- I. [( V; ?& N8 m3 ]0 ?/ F% o   (3,1)       2.6415
   (1,2)       3.9002
   (2,2)     -73.6789
   (3,2)       3.1120
3 @; A/ W& \3 h$ J1 B. W   (4,2)      66.6667
) M, @  y3 l" C/ p5 ~; H, P   (1,3)       2.6415
4 |% D. p) b, P! ], U   (2,3)       3.1120
6 C! w5 c. S6 v. [: \. s6 S   (3,3)     -39.0869
   (2,4)      66.6667
   (4,4)     -66.6667

wlm_28

网上有好像有这个程序

zhyerica

回复 3# wlm_28

+ a- N& r6 o( X- {2 T$ ^/ @
    网上找到的PQ法B1、B2矩阵都是手动输入的，没看到程序实现的。主要是拿5节点的算出来跟书上也差挺多，B2是一样的，B1书上只给了因子表，但是从迭代第一步的结果上看应该就是B1的问题……