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

zhyerica

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

* m( [$ x1 S( V: t2 {# W附上相关部分代码 Ps 我得到的节点导纳矩阵和B2矩阵都是和书上的一样的，真心搞不懂了，求解答啊~~~
8 X0 @+ `% m- o. s3 i

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

%% define named indices into node, branch matrices
[PQ, PV, REF, ...
6 |5 h# A) j" \/ e2 S9 B6 G# A    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;
: w7 N5 ?8 `5 s' k: h
%% form index lists for slack, PV, and PQ buses
" P% h/ _) y# q" |ref = find(Node(:, 2) == 3);   %% reference bus index
; ~  ]# E) N3 J/ m; y" T1 B! L9 f- w6 mpv  = find(Node(:, BUS_TYPE) == PV);    %% PV bus indices
pq  = find(Node(:, BUS_TYPE) == PQ);    %% PQ bus indices

%*************************************************************************%
2 F7 o) s. ~" C8 I/ H%                              form Bp
%*************************************************************************%
( F5 u  u8 V* n/ o9 T! M2 Z
; j. j# G/ g7 K# J%找出非接地支路
4 a& Q3 x3 Y" n1 y; k8 ytemp = find( Branch(:,NODE2)~=0 );
NongroundingBranch = Branch(temp,:);
0 r; a' W/ W. n# G/ S3 N6 E
/ o6 T$ [( a/ a9 Q3 U& F7 O0 L%忽略线路充电电容
NongroundingBranch(:,BORK) = 0;
5 D0 `9 R; Y1 L- l2 y
%不考虑变压器非标准变比
, S6 S) C7 d3 x5 d$ `3 r, s/ qtemp = find( NongroundingBranch(:,BRANCH_TYPE)==TRANS );
NongroundingBranch(temp,BORK) = 1;
. N; O( @+ ?  h$ f) LY = makeY(Node,NongroundingBranch);
Bp = imag( Y );
3 |/ ~3 T9 o2 O4 n$ I+ v7 |Bp = Bp([pq; pv], [pq; pv])  ;  % reduce Bp matrices

5 e* `! L! C" C7 u  A/ R( F- }! F8 X%*************************************************************************%
%                              form Bpp
* M5 G- U* }7 m) K( o8 y; }$ F" R$ W%*************************************************************************%
- ^) S$ q6 ]& `7 s  GBranch(:,4)=0;%忽略输电线路电阻
Y = makeY(Node,Branch);
Bpp = imag( Y );
# R% k$ }( J7 I% b3 }Bpp = Bpp(pq, pq);    % reduce Bp matrices

return;

zhyerica

自己顶~~~
' Z2 @6 V' W5 o, t$ x. [拿那个5节点的例子，算出来的B1矩阵如下：
- g4 e6 s7 R0 P1 iB1 =
& x  E& n. s1 R
   (1,1)      -6.5417
6 H+ m# R/ A. X% f% _# Q   (2,1)       3.9002
   (3,1)       2.6415
   (1,2)       3.9002
8 b9 r& s4 f# g+ P   (2,2)     -73.6789
   (3,2)       3.1120
6 b2 Y1 E: f, c- v   (4,2)      66.6667
   (1,3)       2.6415
   (2,3)       3.1120
   (3,3)     -39.0869
   (2,4)      66.6667
   (4,4)     -66.6667

wlm_28

网上有好像有这个程序

zhyerica

回复 3# wlm_28
8 x3 z& Q9 K' c+ d. L

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