|
|
楼主 |
发表于 2013-9-15 11:11:54
|
显示全部楼层
一下是代码: X) m% i0 W M. ^0 h! \
1 I2 H) B: V, P- [& h, S. ]/ p
function mpc=case5_01
, S: }3 N) p) X8 V$ ^- O0 ?2 A% MATPOWER Case Format:Version 2
/ r' X! t5 C8 e! j- cmpc.version='2';
s+ J+ Z* ~/ g' E" i: E- G%%-----Power Flow Date-----%%$ y3 Q [. P+ k0 h: \2 G
%% system MVA base- I4 u- ~4 u7 l& p# F( u
mpc.baseMVA=100;# ?6 S2 Y' U0 ]# u7 b% Z
%% bus date5 A4 [" c8 b$ Y! c
% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin) @" ~0 p! B, y7 x0 w1 Y1 d
mpc.bus=[6 a4 G) E: ~$ |1 b
1 1 160 80 0 0 1 1 0 100 1 1.1 0.94;
) s! r3 U y4 z! _$ j# A+ \' { 2 1 200 100 0 0 1 1 0 100 1 1.1 0.94;/ U) d) o& p; L1 |3 v1 @
3 1 370 130 0 0 1 1 0 100 1 1.1 0.94;
& }6 `0 G2 E+ ?7 x" Y8 R: y 4 2 0 0 0 0 1 1.050 0 100 1 1.1 0.94;8 G) t; M* T) a* ^6 w
5 3 0 0 0 0 1 1.050 0 100 1 1.1 0.94;) K5 P9 n9 n* ^, G! i1 |. N6 ]
];
1 U* F: J% o6 f8 C0 Z%% generator date
7 }2 `& W" I8 r% i+ E m- l% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin ; N* a0 I' g$ I+ S+ `
mpc.gen =[
# ~! }1 f Y; I' D 4 500 0 99990 -9999 1.050 100 1 600 0;8 s1 _& Z) ^* ^0 f7 r) I `' `! Q
5 0 0 99990 -9999 1.050 100 1 600 0;' x* U8 K N2 W+ }/ f0 y! x) A
];* z5 N( G# G% L% g
%% branch date& l, u9 S* \- T- L0 B! a! q
% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax
8 m0 ]4 n& ~: Z6 `9 R, [5 ?. Gmpc.branch=[* s+ |, z( c/ P$ C5 ^0 {3 Y4 s
2 1 0.04 0.25 0.5 0 0 0 0 0 1 -360 360;
* v5 ^. B4 W: ] J, {' B 3 1 0.1 0.35 0 0 0 0 0 0 1 -360 360;$ \: s! O, n# p/ @7 ]0 i
3 2 0.08 0.3 0.5 0 0 0 0 0 1 -360 360;
6 T+ I: u( g$ }7 X 3 5 0 0.03 0 0 0 0 1.05 0 1 -360 360;
2 V8 p- S& j1 s+ k 2 4 0 0.015 0 0 0 0 1.05 0 1 -360 360;
2 {( s7 J( A$ L& k2 a0 m; ~ ];5 u2 H3 `7 o% y; ]
return; |
-
运行结果
|