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

zhyerica

各位大牛，求助啊~~
我毕设的题目是交直流混联系统的潮流计算，自己编了交替迭代法的的PQ分解法程序，计算结果与书上算例结果误差最大有12%。
于是测试了纯交流部分的程序，也就是一般的PQ法，对照着现代电力系统分析P86上的例子，发现自己程序B1矩阵算到的结果跟书上的有差，致使第一次迭代时相角的修正量跟标准值就差了4%左右。但是查到我程序生成的B1矩阵跟自己手算的结果是一样的啊……
# h/ s& Q, {( V" Q$ d* r# h
附上相关部分代码 Ps 我得到的节点导纳矩阵和B2矩阵都是和书上的一样的，真心搞不懂了，求解答啊~~~
" r* g: h) B5 Z# ]2 g4 _& O: R0 t

function [Bp, Bpp] = makeB(Node,Branch)
' l, F) X3 B* u, d. |6 Z/ ]
%% define named indices into node, branch matrices
[PQ, PV, REF, ...
8 M* d. O* V9 t$ g4 n$ D    TRANS, CAPAC, LINE,...
/ T8 C$ i* Z/ g; L' n2 d    NODE_I, BUS_TYPE,PG,QG, PL, QL, V, THETA, ...
) U0 b) `2 U6 l. ~3 I- r    BRANCHNUM, NODE1, NODE2, R, X, BORK, BRANCH_TYPE, STATUS, P1, Q1, P2, Q2] = index;
8 l7 z- a, z+ N8 |
%% form index lists for slack, PV, and PQ buses
0 q8 ^4 g  s+ [; w! L7 G* Gref = find(Node(:, 2) == 3);   %% reference bus index
* V9 ]8 h8 a# i7 ?' h; ^9 Vpv  = find(Node(:, BUS_TYPE) == PV);    %% PV bus indices
, W7 {- N3 W- ~' Kpq  = find(Node(:, BUS_TYPE) == PQ);    %% PQ bus indices
0 ^3 @  T5 g& w# L. d+ g
" s+ F0 `' V" a6 ?%*************************************************************************%
" q8 L7 ~- @( f' u9 ?/ k%                              form Bp
; s4 N& G- _0 E5 |: A3 N%*************************************************************************%
3 u; Z0 b" b$ J( S
" R. y7 F  @- d9 F9 j" K5 E%找出非接地支路
0 H5 l2 K# C( E, Wtemp = find( Branch(:,NODE2)~=0 );
" Q5 g, O4 q) `- F9 J5 aNongroundingBranch = Branch(temp,:);
! s+ n- j% ^. F
%忽略线路充电电容
NongroundingBranch(:,BORK) = 0;
) \. _0 ?+ M' ]
%不考虑变压器非标准变比
5 {: F! V# c+ H' U) |7 ktemp = find( NongroundingBranch(:,BRANCH_TYPE)==TRANS );
) ]- f0 N( A/ \NongroundingBranch(temp,BORK) = 1;
  m' m2 a6 b) y" @. eY = makeY(Node,NongroundingBranch);
5 p  E$ J7 y3 t- R1 Z4 j9 [5 ]) ?1 wBp = imag( Y );
Bp = Bp([pq; pv], [pq; pv])  ;  % reduce Bp matrices
( G8 ~* ?4 }( g- t3 z+ N$ b% Y
8 Q! N# Y1 J# F- R" z9 w0 \! ~%*************************************************************************%
%                              form Bpp
' q; S) ^; A& b1 B# o& y' w  i%*************************************************************************%
Branch(:,4)=0;%忽略输电线路电阻
Y = makeY(Node,Branch);
Bpp = imag( Y );
Bpp = Bpp(pq, pq);    % reduce Bp matrices
) @: B9 Z' [' O
( c- y( c4 h+ ireturn;

zhyerica

自己顶~~~
/ \6 D& A2 {2 H% p4 h5 B拿那个5节点的例子，算出来的B1矩阵如下：
B1 =

7 u/ i# s) \! E   (1,1)      -6.5417
   (2,1)       3.9002
+ i- f6 B5 K) H3 [5 p   (3,1)       2.6415
   (1,2)       3.9002
* M0 }9 i( l, m$ k% t   (2,2)     -73.6789
' l4 ~$ k/ G- |7 u( h   (3,2)       3.1120
' M# V# T: z2 q. g/ B3 t' t   (4,2)      66.6667
   (1,3)       2.6415
   (2,3)       3.1120
   (3,3)     -39.0869
" I9 M7 x- P* D& U   (2,4)      66.6667
   (4,4)     -66.6667

wlm_28

网上有好像有这个程序

zhyerica

回复 3# wlm_28
8 v4 h8 r, K( ^" F+ t
; w  M3 y  F; B: w
    网上找到的PQ法B1、B2矩阵都是手动输入的，没看到程序实现的。主要是拿5节点的算出来跟书上也差挺多，B2是一样的，B1书上只给了因子表，但是从迭代第一步的结果上看应该就是B1的问题……