|
楼主 |
发表于 2013-9-15 11:11:54
|
显示全部楼层
一下是代码:4 e4 E! I0 Q3 `. K. X
: ^& P+ Y; z$ I1 `" s. n" l
function mpc=case5_01: [$ y5 k, ^+ ~/ u6 K
% MATPOWER Case Format:Version 2
* F* N7 d3 g1 Q8 U+ C8 E7 }mpc.version='2';/ s6 ]; Z" ^7 Z- u6 E4 u. O
%%-----Power Flow Date-----%%
; I0 S# p" R9 U* k, \, i: b%% system MVA base3 n; G; Z" Q+ [) ]+ ~
mpc.baseMVA=100;+ N7 Q# z: x) g4 I
%% bus date
$ V3 Z4 r4 g! A' C9 s% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin7 Y+ q! j$ x& [6 I5 B% k7 F
mpc.bus=[. J1 T! w2 t& Z8 o+ f
1 1 160 80 0 0 1 1 0 100 1 1.1 0.94;/ R; ]1 J8 D% M
2 1 200 100 0 0 1 1 0 100 1 1.1 0.94;$ c$ C3 ]- w, \& v8 a- p
3 1 370 130 0 0 1 1 0 100 1 1.1 0.94;6 h7 M/ I8 d6 l3 p
4 2 0 0 0 0 1 1.050 0 100 1 1.1 0.94;0 \, U' O# d4 J% w& e, w+ }
5 3 0 0 0 0 1 1.050 0 100 1 1.1 0.94;
4 P& ? A0 [+ W9 T1 c; \ ];
9 t X5 n% w2 m0 t' O. S) z2 p%% generator date6 T; L6 Q% H# z# O
% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin 9 t: }3 J3 k% X: B2 h) K
mpc.gen =[, v9 K! Q+ m" f1 Y, Q% q4 m' ]
4 500 0 99990 -9999 1.050 100 1 600 0;
! _( p* T, e) @: M 5 0 0 99990 -9999 1.050 100 1 600 0;" O8 i( W& k# f
];
6 _: o7 m' i" h2 Q, v1 b+ V%% branch date
9 T2 z5 o9 F3 k. n! I/ c% n% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax
* G4 Y; n- ^; {5 ?' r" W+ g1 }mpc.branch=[
, N6 e( b! a; m8 b6 w2 i& ^+ T; q( ?) S, b, \ 2 1 0.04 0.25 0.5 0 0 0 0 0 1 -360 360;
4 w7 p+ X1 [# t! }9 r* e4 A 3 1 0.1 0.35 0 0 0 0 0 0 1 -360 360;
" Q+ H7 U9 l" i, M! h 3 2 0.08 0.3 0.5 0 0 0 0 0 1 -360 360;
& t$ b$ H; W# B5 Q' Z1 m+ b* I 3 5 0 0.03 0 0 0 0 1.05 0 1 -360 360;
8 t; T% L; y6 f# c3 P6 ]) {' ^ 2 4 0 0.015 0 0 0 0 1.05 0 1 -360 360;" t. I" G1 y Y9 k2 y/ P& y
];
# h% n+ z7 {! s2 c0 Ureturn; |
-
运行结果
|