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

zhyerica

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

: n& b, n2 W2 t* b% y  T- O附上相关部分代码 Ps 我得到的节点导纳矩阵和B2矩阵都是和书上的一样的，真心搞不懂了，求解答啊~~~
9 Y7 v- N9 R9 P( E
, ?: f* s& Y* y: P  ?  N
- l6 P7 `# U# v3 Z0 y# Rfunction [Bp, Bpp] = makeB(Node,Branch)
+ T% h  q+ o/ W" ?& B% r
% q5 T- y: a! m' |%% define named indices into node, branch matrices
, G/ ~/ Q, I' {[PQ, PV, REF, ...
    TRANS, CAPAC, LINE,...
    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
ref = find(Node(:, 2) == 3);   %% reference bus index
pv  = find(Node(:, BUS_TYPE) == PV);    %% PV bus indices
, b8 r( {) a; J5 opq  = find(Node(:, BUS_TYPE) == PQ);    %% PQ bus indices

/ y3 P7 H7 x% ]# X%*************************************************************************%
8 u5 B8 f5 ]* a2 V( E. B* c%                              form Bp
%*************************************************************************%
: M$ x- F, A7 B* ^" u
%找出非接地支路
; \; p' H& p/ A* M$ A8 h* ]temp = find( Branch(:,NODE2)~=0 );
8 |( {0 o; T3 [! o3 }/ YNongroundingBranch = Branch(temp,:);
6 b) C9 a: n1 _& S* l8 [
%忽略线路充电电容
  R7 o8 ]6 y8 ]0 LNongroundingBranch(:,BORK) = 0;

%不考虑变压器非标准变比
* h  G8 V: l) M' T+ y; t# Vtemp = find( NongroundingBranch(:,BRANCH_TYPE)==TRANS );
NongroundingBranch(temp,BORK) = 1;
5 |% s& x: \9 v9 P5 a! FY = makeY(Node,NongroundingBranch);
Bp = imag( Y );
Bp = Bp([pq; pv], [pq; pv])  ;  % reduce Bp matrices

%*************************************************************************%
% E& N4 R8 J, D. {0 F%                              form Bpp
6 z; G* z2 i% A) f: m%*************************************************************************%
Branch(:,4)=0;%忽略输电线路电阻
Y = makeY(Node,Branch);
Bpp = imag( Y );
Bpp = Bpp(pq, pq);    % reduce Bp matrices
7 ?- J* X. q+ [: A4 c
+ a$ ]8 J& |% ^& n/ O% s% Kreturn;

zhyerica

自己顶~~~
  Q" i+ F) Z9 N+ S拿那个5节点的例子，算出来的B1矩阵如下：
; z. q4 |4 b4 l4 x( [% b# l) PB1 =

/ N. [* P* z2 ^7 v2 q8 U   (1,1)      -6.5417
4 ?; q; ~' B( p   (2,1)       3.9002
0 v) K% J, L# R5 t0 t   (3,1)       2.6415
% i0 G3 e5 S9 I9 C   (1,2)       3.9002
   (2,2)     -73.6789
   (3,2)       3.1120
; Y* u$ D, w# z3 t4 k   (4,2)      66.6667
: D8 Y4 f1 h: }( p# k   (1,3)       2.6415
   (2,3)       3.1120
$ [6 N/ n" B+ U" [! ~   (3,3)     -39.0869
   (2,4)      66.6667
% \  H0 E" e0 U* I. C5 R1 F   (4,4)     -66.6667

wlm_28

网上有好像有这个程序

zhyerica

回复 3# wlm_28
  v0 i6 a3 {9 F# w/ W0 H
3 D3 S: Y- J! S3 C
    网上找到的PQ法B1、B2矩阵都是手动输入的，没看到程序实现的。主要是拿5节点的算出来跟书上也差挺多，B2是一样的，B1书上只给了因子表，但是从迭代第一步的结果上看应该就是B1的问题……