设为首页收藏本站|繁體中文 快速切换版块

 找回密码
 立即加入
搜索
查看: 1357|回复: 3

[讨论] 求助~~PQ法B1矩阵的形成

[复制链接]

该用户从未签到

尚未签到

发表于 2012-5-6 11:04:01 | 显示全部楼层 |阅读模式

马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!

您需要 登录 才可以下载或查看,没有账号?立即加入

×
各位大牛,求助啊~~2 A/ |# `( Z( w
我毕设的题目是交直流混联系统的潮流计算,自己编了交替迭代法的的PQ分解法程序,计算结果与书上算例结果误差最大有12%。: G; l7 t2 x% Q' p
于是测试了纯交流部分的程序,也就是一般的PQ法,对照着现代电力系统分析P86上的例子,发现自己程序B1矩阵算到的结果跟书上的有差,致使第一次迭代时相角的修正量跟标准值就差了4%左右。但是查到我程序生成的B1矩阵跟自己手算的结果是一样的啊……$ t" J2 X0 j7 \* v- W2 }: |

: 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, ...1 R6 ], ~% p5 }. n6 X" N
    TRANS, CAPAC, LINE,...; u2 Y% F5 u" l" I0 o
    NODE_I, BUS_TYPE,PG,QG, PL, QL, V, THETA, ...- Y3 h' F  {5 p
    BRANCHNUM, NODE1, NODE2, R, X, BORK, BRANCH_TYPE, STATUS, P1, Q1, P2, Q2] = index;7 V" g- m  B: y5 |0 [# ~
% K6 W) X, P& C4 @) r
%% form index lists for slack, PV, and PQ buses. c. F$ h. F& l6 Z2 H8 [/ w
ref = find(Node(:, 2) == 3);   %% reference bus index6 m' b2 t3 S4 c1 c: t
pv  = find(Node(:, BUS_TYPE) == PV);    %% PV bus indices
, b8 r( {) a; J5 opq  = find(Node(:, BUS_TYPE) == PQ);    %% PQ bus indices$ \0 H8 V- w) ^% b, t0 ]0 r

/ y3 P7 H7 x% ]# X%*************************************************************************%
8 u5 B8 f5 ]* a2 V( E. B* c%                              form Bp% |  ?# L7 S  P; y1 K4 E/ x2 @
%*************************************************************************%
: M$ x- F, A7 B* ^" u& t- t* {( I6 O0 k" m0 x
%找出非接地支路
; \; 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 [( B! R) t+ E) s3 W+ h+ i3 f1 d- R
%忽略线路充电电容
  R7 o8 ]6 y8 ]0 LNongroundingBranch(:,BORK) = 0;* S, N( n+ ~" ?: E* \
+ l6 z0 @4 m* ?2 ~, i- `
%不考虑变压器非标准变比
* h  G8 V: l) M' T+ y; t# Vtemp = find( NongroundingBranch(:,BRANCH_TYPE)==TRANS );7 R1 J4 K0 p/ q1 a& `4 V, [6 E
NongroundingBranch(temp,BORK) = 1;
5 |% s& x: \9 v9 P5 a! FY = makeY(Node,NongroundingBranch);$ t- G$ A. Z% ]  q+ `1 R
Bp = imag( Y );" P% T  ^0 T1 r5 t+ [2 c& b2 c3 S
Bp = Bp([pq; pv], [pq; pv])  ;  % reduce Bp matrices6 l0 e2 C0 h1 v+ w0 w
" q; q. w! b  s* K% s
%*************************************************************************%
% E& N4 R8 J, D. {0 F%                              form Bpp
6 z; G* z2 i% A) f: m%*************************************************************************%5 e2 E/ V* N, M4 V8 C* X6 {9 C
Branch(:,4)=0;%忽略输电线路电阻2 Z3 z2 [+ \, E% e1 q" z2 l. n
Y = makeY(Node,Branch);! e, l' p# ~% a' q
Bpp = imag( Y );- D7 Y, ]1 I8 U' M$ t. O8 ^$ X
Bpp = Bpp(pq, pq);    % reduce Bp matrices
7 ?- J* X. q+ [: A4 c
+ a$ ]8 J& |% ^& n/ O% s% Kreturn;
楼主热帖
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

该用户从未签到

尚未签到

 楼主| 发表于 2012-5-6 14:37:12 | 显示全部楼层
自己顶~~~
  Q" i+ F) Z9 N+ S拿那个5节点的例子,算出来的B1矩阵如下:
; z. q4 |4 b4 l4 x( [% b# l) PB1 =7 U- O: j/ N5 B3 D

/ 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' O: [' ?6 B! E1 b2 q0 z* D- Q
   (2,2)     -73.6789& @3 U, t. x6 i" ^! s; s
   (3,2)       3.1120
; Y* u$ D, w# z3 t4 k   (4,2)      66.6667
: D8 Y4 f1 h: }( p# k   (1,3)       2.64155 m1 M4 f% K/ C; x3 @- v
   (2,3)       3.1120
$ [6 N/ n" B+ U" [! ~   (3,3)     -39.08696 N  [( @$ y+ q, n3 ?) ]: X7 s( U3 B
   (2,4)      66.6667
% \  H0 E" e0 U* I. C5 R1 F   (4,4)     -66.6667
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    无聊
    2017-9-26 15:28
  • 签到天数: 21 天

    连续签到: 1 天

    [LV.4]偶尔看看III

    累计签到:21 天
    连续签到:1 天
    发表于 2012-5-6 20:37:46 | 显示全部楼层
    网上有好像有这个程序
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

    该用户从未签到

    尚未签到

     楼主| 发表于 2012-5-7 10:04:07 | 显示全部楼层
    回复 3# wlm_28
      v0 i6 a3 {9 F# w/ W0 H
    3 D3 S: Y- J! S3 C" `/ q- k0 h2 R2 `7 n- t' V; \
        网上找到的PQ法B1、B2矩阵都是手动输入的,没看到程序实现的。主要是拿5节点的算出来跟书上也差挺多,B2是一样的,B1书上只给了因子表,但是从迭代第一步的结果上看应该就是B1的问题……
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
    您需要登录后才可以回帖 登录 | 立即加入

    本版积分规则

    招聘斑竹

    小黑屋|手机版|APP下载(beta)|Archiver|电力研学网 ( 赣ICP备12000811号-1|赣公网安备36040302000210号 )|网站地图

    GMT+8, 2026-3-16 20:15

    Powered by Discuz! X3.5 Licensed

    © 2001-2025 Discuz! Team.

    快速回复 返回顶部 返回列表