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

zhyerica

各位大牛，求助啊~~
5 k* n5 g4 R- Y: @3 S2 }  M1 c我毕设的题目是交直流混联系统的潮流计算，自己编了交替迭代法的的PQ分解法程序，计算结果与书上算例结果误差最大有12%。
于是测试了纯交流部分的程序，也就是一般的PQ法，对照着现代电力系统分析P86上的例子，发现自己程序B1矩阵算到的结果跟书上的有差，致使第一次迭代时相角的修正量跟标准值就差了4%左右。但是查到我程序生成的B1矩阵跟自己手算的结果是一样的啊……
: u( Z% E8 _1 e4 s; d7 j
  Y2 u4 l! Q. a3 s2 B* c附上相关部分代码 Ps 我得到的节点导纳矩阵和B2矩阵都是和书上的一样的，真心搞不懂了，求解答啊~~~
: W+ G, a* G) [

function [Bp, Bpp] = makeB(Node,Branch)
* j: _/ b! s$ M2 Q  J6 n4 Y/ p
%% define named indices into node, branch matrices
[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;

; U: ~. n6 V9 M1 [%% 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
pq  = find(Node(:, BUS_TYPE) == PQ);    %% PQ bus indices
0 Q7 I* D$ \" }5 Q" }, i3 Z( F
( J; {# ?9 L: T7 i4 ^%*************************************************************************%
%                              form Bp
. \" c* P% H$ V0 ^! Q2 ]%*************************************************************************%
' S, k; Q! ^6 e" |# H, g
1 W( l7 t, p0 r# C8 B: x6 `%找出非接地支路
temp = find( Branch(:,NODE2)~=0 );
+ d& X* N8 R  d! v( I; }; n: MNongroundingBranch = Branch(temp,:);

' _) G0 U$ x1 {2 U%忽略线路充电电容
NongroundingBranch(:,BORK) = 0;

, @7 C& H! ]. ?2 I3 q' d%不考虑变压器非标准变比
temp = find( NongroundingBranch(:,BRANCH_TYPE)==TRANS );
NongroundingBranch(temp,BORK) = 1;
* B7 v- w! h7 L$ Z- QY = makeY(Node,NongroundingBranch);
Bp = imag( Y );
Bp = Bp([pq; pv], [pq; pv])  ;  % reduce Bp matrices
/ ]2 Y) h1 n3 o$ I% c
%*************************************************************************%
- y; @4 z% R& b  r; C' y%                              form Bpp
* a$ Y( ]# l4 E5 ]% s8 S6 \%*************************************************************************%
Branch(:,4)=0;%忽略输电线路电阻
% e/ c) v) H; k1 g) I0 S; \4 mY = makeY(Node,Branch);
Bpp = imag( Y );
/ B0 O  D. P; V0 kBpp = Bpp(pq, pq);    % reduce Bp matrices
7 {% F4 g: A0 s3 v. n9 N
return;

zhyerica

自己顶~~~
拿那个5节点的例子，算出来的B1矩阵如下：
: ^# [7 b. Z4 s/ h9 aB1 =

   (1,1)      -6.5417
   (2,1)       3.9002
- K1 V' y6 @' J   (3,1)       2.6415
& O1 i7 t/ W8 j; w1 X! z8 p   (1,2)       3.9002
   (2,2)     -73.6789
   (3,2)       3.1120
   (4,2)      66.6667
   (1,3)       2.6415
0 d; ]" e# x& S   (2,3)       3.1120
4 C  x( f4 R5 ~: H5 H   (3,3)     -39.0869
& a/ E$ Q9 R3 c$ Z' A7 U, Q   (2,4)      66.6667
   (4,4)     -66.6667

wlm_28

网上有好像有这个程序

zhyerica

回复 3# wlm_28
; _# B. @9 @1 t; O* P
1 [8 z4 ~0 U4 Q% i, B, x
    网上找到的PQ法B1、B2矩阵都是手动输入的，没看到程序实现的。主要是拿5节点的算出来跟书上也差挺多，B2是一样的，B1书上只给了因子表，但是从迭代第一步的结果上看应该就是B1的问题……