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

zhyerica

各位大牛，求助啊~~
我毕设的题目是交直流混联系统的潮流计算，自己编了交替迭代法的的PQ分解法程序，计算结果与书上算例结果误差最大有12%。
$ P% S  x( y1 s: G; i4 R于是测试了纯交流部分的程序，也就是一般的PQ法，对照着现代电力系统分析P86上的例子，发现自己程序B1矩阵算到的结果跟书上的有差，致使第一次迭代时相角的修正量跟标准值就差了4%左右。但是查到我程序生成的B1矩阵跟自己手算的结果是一样的啊……

9 q7 N% F- E1 n% _& u) P& z1 ~附上相关部分代码 Ps 我得到的节点导纳矩阵和B2矩阵都是和书上的一样的，真心搞不懂了，求解答啊~~~

5 c! I3 u' Y/ l+ F) S
. i. J, T; F  t4 H, Vfunction [Bp, Bpp] = makeB(Node,Branch)
; X5 K8 L, U$ G: _* f* `* ]5 U
%% define named indices into node, branch matrices
[PQ, PV, REF, ...
4 e# x3 I  |+ A' {9 B1 R' T/ t    TRANS, CAPAC, LINE,...
6 N0 u8 [( L) X9 t5 \    NODE_I, BUS_TYPE,PG,QG, PL, QL, V, THETA, ...
    BRANCHNUM, NODE1, NODE2, R, X, BORK, BRANCH_TYPE, STATUS, P1, Q1, P2, Q2] = index;

%% form index lists for slack, PV, and PQ buses
3 j; Q2 N" r, ~2 Kref = find(Node(:, 2) == 3);   %% reference bus index
pv  = find(Node(:, BUS_TYPE) == PV);    %% PV bus indices
3 N$ P& D+ d+ w1 vpq  = find(Node(:, BUS_TYPE) == PQ);    %% PQ bus indices

%*************************************************************************%
%                              form Bp
% M+ i$ _. j/ L: J6 p5 ~%*************************************************************************%
7 b2 |; h# |9 b8 `$ x* K8 `% J# g2 M
%找出非接地支路
! S* b4 ?- o( U) vtemp = find( Branch(:,NODE2)~=0 );
NongroundingBranch = Branch(temp,:);

%忽略线路充电电容
8 s6 |! W5 q& @& d$ j  l: MNongroundingBranch(:,BORK) = 0;
- H3 E" A+ ^  M
4 J' _* K. o6 w3 |4 ?/ @! U6 i) H& z%不考虑变压器非标准变比
temp = find( NongroundingBranch(:,BRANCH_TYPE)==TRANS );
2 X: y7 ]( N% |; V. x2 _NongroundingBranch(temp,BORK) = 1;
/ j( g) ~) u6 n- v$ eY = makeY(Node,NongroundingBranch);
" U4 r/ {& \: f( s" HBp = imag( Y );
, ?  S7 x1 w( BBp = Bp([pq; pv], [pq; pv])  ;  % reduce Bp matrices
& E& e3 H* ]! C# v, D
& T1 g( \: o+ ?* L%*************************************************************************%
%                              form Bpp
6 w7 ~! E; \; Q+ `%*************************************************************************%
( m- ?3 }6 T7 d/ M5 sBranch(:,4)=0;%忽略输电线路电阻
+ d' f! q# @9 YY = makeY(Node,Branch);
) J) X1 D# j, E( YBpp = imag( Y );
9 r- N7 L9 p& Q; \* X( t% l8 p3 YBpp = Bpp(pq, pq);    % reduce Bp matrices
' K+ J$ P6 [& o% s( D( p
3 o* ]9 ^* ?5 R) U+ g$ \7 l# Breturn;

zhyerica

自己顶~~~
( K+ `& S. [: ]4 n9 F8 S6 U拿那个5节点的例子，算出来的B1矩阵如下：
7 B0 _, m# H9 \B1 =
6 }$ C2 C* L& T- X
   (1,1)      -6.5417
   (2,1)       3.9002
; W" Z* o5 ?# g; l6 D7 u  ]% P   (3,1)       2.6415
   (1,2)       3.9002
& W. p$ P- @& R+ c& z! J; ?  `2 a   (2,2)     -73.6789
   (3,2)       3.1120
   (4,2)      66.6667
   (1,3)       2.6415
   (2,3)       3.1120
& B" B; ?. t, Z- j7 Q   (3,3)     -39.0869
) H: M' @8 @3 [/ C3 I   (2,4)      66.6667
- [  v& Q( L! y   (4,4)     -66.6667

wlm_28

网上有好像有这个程序

zhyerica

回复 3# wlm_28
7 M2 Z5 t& m6 }0 Y1 `. q' r5 `

% [5 d" L$ @6 M0 ~% E# n0 C2 r+ N    网上找到的PQ法B1、B2矩阵都是手动输入的，没看到程序实现的。主要是拿5节点的算出来跟书上也差挺多，B2是一样的，B1书上只给了因子表，但是从迭代第一步的结果上看应该就是B1的问题……