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

zhyerica

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

附上相关部分代码 Ps 我得到的节点导纳矩阵和B2矩阵都是和书上的一样的，真心搞不懂了，求解答啊~~~
- b) m+ D1 q% R, p

function [Bp, Bpp] = makeB(Node,Branch)

%% define named indices into node, branch matrices
[PQ, PV, REF, ...
    TRANS, CAPAC, LINE,...
1 a6 C& t$ _5 [    NODE_I, BUS_TYPE,PG,QG, PL, QL, V, THETA, ...
6 C3 S. W. j3 A- l7 O    BRANCHNUM, NODE1, NODE2, R, X, BORK, BRANCH_TYPE, STATUS, P1, Q1, P2, Q2] = index;
! B$ M, c" f* O$ e
& h) g# g! N, J# n0 j%% form index lists for slack, PV, and PQ buses
0 Q6 M3 H* v6 u! X9 Cref = find(Node(:, 2) == 3);   %% reference bus index
pv  = find(Node(:, BUS_TYPE) == PV);    %% PV bus indices
pq  = find(Node(:, BUS_TYPE) == PQ);    %% PQ bus indices
( E0 l5 |& z' T( m0 P9 K
) u& s* C/ |- y# t( }3 r1 i%*************************************************************************%
9 Q  w5 X1 v5 K7 ~%                              form Bp
$ g" p- z. Z/ t: d1 _7 P0 ^' X' a%*************************************************************************%
9 \# \  _- C% c' [
%找出非接地支路
temp = find( Branch(:,NODE2)~=0 );
& l9 A% ^9 _( H. bNongroundingBranch = Branch(temp,:);

6 ^1 x, x1 `7 p$ W%忽略线路充电电容
! K5 H$ o5 V5 T+ ?+ g' I+ aNongroundingBranch(:,BORK) = 0;

" C" U# X& c- b( U% S%不考虑变压器非标准变比
temp = find( NongroundingBranch(:,BRANCH_TYPE)==TRANS );
NongroundingBranch(temp,BORK) = 1;
- p% {$ `/ s; `, P6 y: u  JY = makeY(Node,NongroundingBranch);
, g5 e% S& x+ j/ ^! ZBp = imag( Y );
Bp = Bp([pq; pv], [pq; pv])  ;  % reduce Bp matrices

- H- w# q6 W$ y3 q, |& P8 j%*************************************************************************%
6 L+ X  N4 n# ^1 J8 ]- F4 `" _$ ^%                              form Bpp
%*************************************************************************%
0 n0 a) W% d- U/ A3 G6 B: |$ e* dBranch(:,4)=0;%忽略输电线路电阻
Y = makeY(Node,Branch);
  [& Q0 ?" A' G3 l& T. sBpp = imag( Y );
Bpp = Bpp(pq, pq);    % reduce Bp matrices

& R* C9 u$ [, o* R, u, O% c$ Zreturn;

zhyerica

自己顶~~~
4 b$ [! u6 }/ @1 o( U拿那个5节点的例子，算出来的B1矩阵如下：
1 P. G. ]. X4 ^B1 =
- Z- }2 \4 V) w; \" A! d& |
1 W+ q3 V2 B3 ~' K2 |3 \6 g   (1,1)      -6.5417
5 P7 }- J5 H/ ~# p   (2,1)       3.9002
* a% t3 U5 i+ Q# m, a9 P, f/ f; G   (3,1)       2.6415
   (1,2)       3.9002
: p5 M) e; Z, d1 j4 @! g) {; {   (2,2)     -73.6789
   (3,2)       3.1120
  _, @; J) V) P2 ~% T1 n! E   (4,2)      66.6667
   (1,3)       2.6415
  {1 w6 {* U- c; Y" J# Y   (2,3)       3.1120
9 o+ y6 C. C5 v5 W/ c   (3,3)     -39.0869
   (2,4)      66.6667
   (4,4)     -66.6667

wlm_28

网上有好像有这个程序

zhyerica

回复 3# wlm_28
7 o5 K- k1 D+ A, I1 ^3 n1 U

    网上找到的PQ法B1、B2矩阵都是手动输入的，没看到程序实现的。主要是拿5节点的算出来跟书上也差挺多，B2是一样的，B1书上只给了因子表，但是从迭代第一步的结果上看应该就是B1的问题……