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

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

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

[复制链接]

该用户从未签到

尚未签到

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

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

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

×
各位大牛,求助啊~~' b" u) v. f/ N1 h5 F
我毕设的题目是交直流混联系统的潮流计算,自己编了交替迭代法的的PQ分解法程序,计算结果与书上算例结果误差最大有12%。
6 M+ X; Y! w8 q, |& Q于是测试了纯交流部分的程序,也就是一般的PQ法,对照着现代电力系统分析P86上的例子,发现自己程序B1矩阵算到的结果跟书上的有差,致使第一次迭代时相角的修正量跟标准值就差了4%左右。但是查到我程序生成的B1矩阵跟自己手算的结果是一样的啊……4 j4 O" O, \% W6 U6 b/ B" d
4 O3 G/ ?+ _" t, {8 W6 G7 E: n: p
附上相关部分代码 Ps 我得到的节点导纳矩阵和B2矩阵都是和书上的一样的,真心搞不懂了,求解答啊~~~
- b) m+ D1 q% R, p, \; n) P- O+ {2 C; b' b+ a& z
) V( f7 f  g* D  P: J. n3 N
function [Bp, Bpp] = makeB(Node,Branch)7 ^3 n$ T9 |( B; k  y  j; X2 \
2 j5 E, l  e1 k/ l  d. a/ F; a
%% define named indices into node, branch matrices$ B3 ^! T3 g6 S
[PQ, PV, REF, ...; E* \9 r4 q0 t$ e3 m
    TRANS, CAPAC, LINE,...
1 a6 C& t$ _5 [    NODE_I, BUS_TYPE,PG,QG, PL, QL, V, THETA, ...
6 C3 S. W. j3 A- l7 O    BRANCHNUM, NODE1, NODE2, R, X, BORK, BRANCH_TYPE, STATUS, P1, Q1, P2, Q2] = index;
! B$ M, c" f* O$ e
& h) g# g! N, J# n0 j%% form index lists for slack, PV, and PQ buses
0 Q6 M3 H* v6 u! X9 Cref = find(Node(:, 2) == 3);   %% reference bus index/ e2 t% z; \- v, |" y5 M8 n  s
pv  = find(Node(:, BUS_TYPE) == PV);    %% PV bus indices7 q6 w% S8 h: K3 p7 D' I# B
pq  = find(Node(:, BUS_TYPE) == PQ);    %% PQ bus indices
( E0 l5 |& z' T( m0 P9 K
) u& s* C/ |- y# t( }3 r1 i%*************************************************************************%
9 Q  w5 X1 v5 K7 ~%                              form Bp
$ g" p- z. Z/ t: d1 _7 P0 ^' X' a%*************************************************************************%
9 \# \  _- C% c' [- d( Z  C1 o2 G/ O& m$ q
%找出非接地支路. L& W8 j) u3 O: A: M. M7 q4 }
temp = find( Branch(:,NODE2)~=0 );
& l9 A% ^9 _( H. bNongroundingBranch = Branch(temp,:);2 e0 v; @! q( ~9 o0 x. ?$ V

6 ^1 x, x1 `7 p$ W%忽略线路充电电容
! K5 H$ o5 V5 T+ ?+ g' I+ aNongroundingBranch(:,BORK) = 0;( S, u6 g8 `) z6 ]6 R4 q, w/ ]! L/ q0 G

" C" U# X& c- b( U% S%不考虑变压器非标准变比& _% r( b1 t1 \6 N, V2 |
temp = find( NongroundingBranch(:,BRANCH_TYPE)==TRANS );2 B0 p: l4 J+ o$ q( e9 v' H
NongroundingBranch(temp,BORK) = 1;
- p% {$ `/ s; `, P6 y: u  JY = makeY(Node,NongroundingBranch);
, g5 e% S& x+ j/ ^! ZBp = imag( Y );4 S2 M, _' W1 @* E6 t  x
Bp = Bp([pq; pv], [pq; pv])  ;  % reduce Bp matrices+ f0 F& a+ `: ?# o2 Q

- H- w# q6 W$ y3 q, |& P8 j%*************************************************************************%
6 L+ X  N4 n# ^1 J8 ]- F4 `" _$ ^%                              form Bpp6 G2 }- A) D6 N0 R2 R# [
%*************************************************************************%
0 n0 a) W% d- U/ A3 G6 B: |$ e* dBranch(:,4)=0;%忽略输电线路电阻5 p# X' ]! f: l+ @' N
Y = makeY(Node,Branch);
  [& Q0 ?" A' G3 l& T. sBpp = imag( Y );4 i: i0 l1 D2 o: n% @1 B
Bpp = Bpp(pq, pq);    % reduce Bp matrices/ z2 g( h0 ~& g* H! `8 ^0 l- m

& R* C9 u$ [, o* R, u, O% c$ Zreturn;
楼主热帖
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

该用户从未签到

尚未签到

 楼主| 发表于 2012-5-6 14:37:12 | 显示全部楼层
自己顶~~~
4 b$ [! u6 }/ @1 o( U拿那个5节点的例子,算出来的B1矩阵如下:
1 P. G. ]. X4 ^B1 =
- Z- }2 \4 V) w; \" A! d& |
1 W+ q3 V2 B3 ~' K2 |3 \6 g   (1,1)      -6.5417
5 P7 }- J5 H/ ~# p   (2,1)       3.9002
* a% t3 U5 i+ Q# m, a9 P, f/ f; G   (3,1)       2.6415' Q& q) n: M% i' s( T  K
   (1,2)       3.9002
: p5 M) e; Z, d1 j4 @! g) {; {   (2,2)     -73.6789  u; w; |6 s; v$ J; a% l
   (3,2)       3.1120
  _, @; J) V) P2 ~% T1 n! E   (4,2)      66.6667- x$ B0 L& ~- t5 k" K6 E9 R' u! c4 T  ^
   (1,3)       2.6415
  {1 w6 {* U- c; Y" J# Y   (2,3)       3.1120
9 o+ y6 C. C5 v5 W/ c   (3,3)     -39.0869! W6 c( l( j' z) R% x+ B
   (2,4)      66.6667. h0 }  ~+ N% ]' E
   (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
    7 o5 K- k1 D+ A, I1 ^3 n1 U4 r3 z1 h8 L- O$ t2 F) ^& A
    0 I7 g" p) f0 j- X. k7 {7 P: L8 B
        网上找到的PQ法B1、B2矩阵都是手动输入的,没看到程序实现的。主要是拿5节点的算出来跟书上也差挺多,B2是一样的,B1书上只给了因子表,但是从迭代第一步的结果上看应该就是B1的问题……
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
    您需要登录后才可以回帖 登录 | 立即加入

    本版积分规则

    招聘斑竹

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

    GMT+8, 2026-3-16 21:26

    Powered by Discuz! X3.5 Licensed

    © 2001-2025 Discuz! Team.

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