重提问 PQ分解法B1矩阵该如何求啊

zhyerica

[讨论] 求助~~PQ法B1矩阵的形成
3 o7 }# d9 L# T/ i
* D( j0 {5 H3 O, ^( F- y求助, PQ分解法
各位大牛，求助啊~~
我毕设的题目是交直流混联系统的潮流计算，自己编了交替迭代法的的PQ分解法程序，计算结果与书上算例结果误差最大有12%。
5 F( U1 y3 C5 {3 K5 v. k% x于是测试了纯交流部分的程序，也就是一般的PQ法，对照着现代电力系统分析P86上的例子，发现自己程序B1矩阵算到的结果跟书上的有差，致使第一次迭代时相角的修正量跟标准值就差了4%左右。但是查到我程序生成的B1矩阵跟自己手算的结果是一样的啊……
/ @$ O$ ~* f$ I9 z! ^- A4 U
附上相关部分代码 Ps 我得到的节点导纳矩阵和B2矩阵都是和书上的一样的，真心搞不懂了，求解答啊~~~


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

%% 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;
  I0 W) U0 Q7 y+ }
1 M" ?( R. X5 U6 Z* |0 l- t%% form index lists for slack, PV, and PQ buses
# q$ z* C& K, [, d' H" v# Zref = find(Node(:, 2) == 3);   %% reference bus index
) t  @* S: r6 v* Ipv  = find(Node(:, BUS_TYPE) == PV);    %% PV bus indices
  K- ^, u$ b# Lpq  = find(Node(:, BUS_TYPE) == PQ);    %% PQ bus indices
  ?5 e% e3 n0 `% V- f& s
% @3 @+ n' x; y$ ~%*************************************************************************%
%                              form Bp
( R% f8 ~" G+ ~' y%*************************************************************************%
. L' J1 p3 H: m5 w; m5 w
%找出非接地支路
temp = find( Branch(:,NODE2)~=0 );
8 ^6 V+ z5 u8 L2 W. Q, VNongroundingBranch = Branch(temp,:);
% Q- N; Z/ U5 a" Q  _
) t0 [) w9 y2 P%忽略线路充电电容
NongroundingBranch(:,BORK) = 0;

2 u" m! u* z% p2 i" A4 M%不考虑变压器非标准变比
temp = find( NongroundingBranch(:,BRANCH_TYPE)==TRANS );
6 S5 j6 Y8 f/ w( N# z- t2 C2 {NongroundingBranch(temp,BORK) = 1;
) [9 l7 v7 S0 q. Y' T9 Z4 SY = makeY(Node,NongroundingBranch);
$ G9 U4 ~+ X2 t( ABp = imag( Y );
, \& Q* u9 P: [' gBp = Bp([pq; pv], [pq; pv])  ;  % reduce Bp matrices
6 y5 {' ~) y: c) _4 o
%*************************************************************************%
%                              form Bpp
) B( z" A$ A+ u% a4 \% S%*************************************************************************%
Branch(:,4)=0;%忽略输电线路电阻
Y = makeY(Node,Branch);
& E4 {2 j3 D+ Q; C) P! ^/ A9 t. PBpp = imag( Y );
3 b' Z/ n+ c# S( ?/ jBpp = Bpp(pq, pq);    % reduce Bp matrices

loverqi

这些基础的我都不会，那书我们也学的，。。。要学习