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

zhyerica

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

附上相关部分代码 Ps 我得到的节点导纳矩阵和B2矩阵都是和书上的一样的，真心搞不懂了，求解答啊~~~
, u  Q8 @7 k3 j: h  ^) [
+ b$ ~8 ~+ T) M! ^9 d
& n5 T% U9 Y* _! x1 @; w6 Mfunction [Bp, Bpp] = makeB(Node,Branch)
/ |' y$ k/ A$ c2 G
' R/ H# V; Y" f- B. F/ D%% define named indices into node, branch matrices
  X1 p( U% s7 B( _; L. l[PQ, PV, REF, ...
6 i3 @2 }8 @- h1 G    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;

& B: j# O0 X0 Y7 S* @; S%% form index lists for slack, PV, and PQ buses
5 i% k3 c- s! L: Y& b- xref = 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

%*************************************************************************%
% g* \5 N! Y, U+ H%                              form Bp
%*************************************************************************%

%找出非接地支路
temp = find( Branch(:,NODE2)~=0 );
1 W8 u' E% b! D8 h% d0 @& n5 S  WNongroundingBranch = Branch(temp,:);

; P' r1 Z( O, T3 ]  f%忽略线路充电电容
NongroundingBranch(:,BORK) = 0;
. q% [) D6 F# V% T( q0 W
$ X( v, @% s1 t" h  z%不考虑变压器非标准变比
temp = find( NongroundingBranch(:,BRANCH_TYPE)==TRANS );
NongroundingBranch(temp,BORK) = 1;
6 A/ l' T2 N' \/ A1 R$ r, T1 BY = makeY(Node,NongroundingBranch);
( j* A% y! [# b/ s# `/ H" FBp = imag( Y );
Bp = Bp([pq; pv], [pq; pv])  ;  % reduce Bp matrices

# S/ G$ b) @: D' v+ R& r8 B%*************************************************************************%
9 s: T. q$ n1 j3 _%                              form Bpp
/ X( p$ |- q  D2 V%*************************************************************************%
, G: Q& j1 w7 b& ^( i% E+ F8 ]Branch(:,4)=0;%忽略输电线路电阻
+ w+ \+ ~+ Z0 {8 b+ vY = makeY(Node,Branch);
Bpp = imag( Y );
Bpp = Bpp(pq, pq);    % reduce Bp matrices

return;

zhyerica

自己顶~~~
拿那个5节点的例子，算出来的B1矩阵如下：
B1 =

7 D' i) _' l& u3 B. n) h- O8 j   (1,1)      -6.5417
   (2,1)       3.9002
   (3,1)       2.6415
& k# C+ X: y6 v  l   (1,2)       3.9002
   (2,2)     -73.6789
   (3,2)       3.1120
( W" J/ |( z# O   (4,2)      66.6667
! e  n5 B9 ?. N9 C2 q   (1,3)       2.6415
   (2,3)       3.1120
1 S" Y* f$ p  `# L) @* W$ R6 H, O   (3,3)     -39.0869
/ Q8 F- w6 Z0 h- q8 `3 ~9 M   (2,4)      66.6667
* R; o7 V* H7 v. M8 ^   (4,4)     -66.6667

wlm_28

网上有好像有这个程序

zhyerica

回复 3# wlm_28
6 t7 L6 ?1 l+ g5 ?: Z& ?
( E( T8 m7 }+ W! \3 d% t
    网上找到的PQ法B1、B2矩阵都是手动输入的，没看到程序实现的。主要是拿5节点的算出来跟书上也差挺多，B2是一样的，B1书上只给了因子表，但是从迭代第一步的结果上看应该就是B1的问题……