|
|
|
马上加入，结交更多好友，共享更多资料，让你轻松玩转电力研学社区！
您需要 登录 才可以下载或查看，没有账号？立即加入
×
各位大牛，求助啊~~
+ m' J& y: O: ` A1 f8 d" [我毕设的题目是交直流混联系统的潮流计算，自己编了交替迭代法的的PQ分解法程序，计算结果与书上算例结果误差最大有12%。
( c5 Y: N3 r8 Y/ i I1 K4 [5 L! ?于是测试了纯交流部分的程序，也就是一般的PQ法，对照着现代电力系统分析P86上的例子，发现自己程序B1矩阵算到的结果跟书上的有差，致使第一次迭代时相角的修正量跟标准值就差了4%左右。但是查到我程序生成的B1矩阵跟自己手算的结果是一样的啊……. {9 x) e; J: \0 e2 F; o
% v: `- l4 w- `1 T; g- G/ Q* e$ N附上相关部分代码 Ps 我得到的节点导纳矩阵和B2矩阵都是和书上的一样的，真心搞不懂了，求解答啊~~~! H& m+ a" Z8 j
3 ~5 J$ e1 K5 F1 x' S) z" r) O% c$ N+ K
function [Bp, Bpp] = makeB(Node,Branch)
Z4 l- I7 \% d; R6 W. a
& W& J. q& h1 v5 f7 N$ g%% define named indices into node, branch matrices) @" K& Y: g2 [0 v9 ]- p, _+ C
[PQ, PV, REF, ...
9 Q7 A; C) d/ k* h9 { TRANS, CAPAC, LINE,...
/ o( M" I" ~) p g/ z% G% T NODE_I, BUS_TYPE,PG,QG, PL, QL, V, THETA, ...
' t- F+ R3 ]! x! O BRANCHNUM, NODE1, NODE2, R, X, BORK, BRANCH_TYPE, STATUS, P1, Q1, P2, Q2] = index;
9 ~' z) `! `/ e2 u( t9 M' `
) k! v+ j9 \, ]3 F%% form index lists for slack, PV, and PQ buses. J. u1 ?1 e9 q- Z1 L
ref = find(Node(:, 2) == 3); %% reference bus index3 _* y9 m4 m* B2 [
pv = find(Node(:, BUS_TYPE) == PV); %% PV bus indices3 i5 Y# z! M. }) g* t7 `7 m
pq = find(Node(:, BUS_TYPE) == PQ); %% PQ bus indices
. K$ o% M; V- w& Q9 E6 A9 y; T0 b; N; H3 z
%*************************************************************************%5 |5 f/ M! O& b9 S. F
% form Bp a2 C/ K# W0 g' y
%*************************************************************************%
1 E* N, p# ^3 H J0 g, V' t6 P$ N% y' ^& d9 R1 ~- w+ t& W; M
%找出非接地支路0 A. S- R7 _9 N; O" m3 @- L
temp = find( Branch(:,NODE2)~=0 );
/ u+ d! o8 G7 F' |NongroundingBranch = Branch(temp,:);/ p- {4 w8 d2 N: {; ^
7 _, ?) B$ J3 t$ ]) f! q% e6 T%忽略线路充电电容
$ v1 g8 w; Z/ Y) x5 l. v9 @NongroundingBranch(:,BORK) = 0;5 E" G5 y0 X, v* x% L. ]+ G
9 w" b8 `+ k8 Y) v% l4 X
%不考虑变压器非标准变比+ j# u$ C6 n( N5 H7 q
temp = find( NongroundingBranch(:,BRANCH_TYPE)==TRANS );( a7 Z F' q7 I4 @+ ?# ~: b
NongroundingBranch(temp,BORK) = 1;
3 _' p" P7 l2 R( P Y& VY = makeY(Node,NongroundingBranch);5 [# l2 @6 o6 L8 C7 I. Y
Bp = imag( Y );4 L7 `* |" S/ V
Bp = Bp([pq; pv], [pq; pv]) ; % reduce Bp matrices
- A6 _3 x: l2 C: M7 J1 K S9 S7 l- P+ A" V# u8 C* E
%*************************************************************************%
, |' v! Z n4 q+ d: V/ B2 t% form Bpp
. u: u- I- q0 P" ~8 i% ~( I%*************************************************************************%% [" | w5 P, }1 i) _5 ~
Branch(:,4)=0;%忽略输电线路电阻1 E S6 U* b4 i* ~5 R
Y = makeY(Node,Branch);
0 a5 e- f. n" `% O# u3 }( hBpp = imag( Y );
6 Z/ Q- S' }* A+ `/ fBpp = Bpp(pq, pq); % reduce Bp matrices- b+ v1 z4 |- S; y
* P) I5 b* w, b
return; |
|
赣公网安备36040302000210号 )|网站地图
狗仔卡
发表于 2012-5-6 11:04:01
提升卡
置顶卡
沉默卡
喧嚣卡
变色卡
显身卡