|
|
楼主 |
发表于 2013-9-15 11:11:54
|
显示全部楼层
一下是代码:
4 ~" H5 p Z6 h% T7 ~( C1 c# e& l: }
# b/ D( P$ d. t/ [: ?0 @3 ]) pfunction mpc=case5_014 f1 p+ G$ K/ g3 H, v% }
% MATPOWER Case Format:Version 2
9 C! ~) K! R9 h$ I6 u w+ ^! d+ @mpc.version='2';0 Z: A( s) X! R# @
%%-----Power Flow Date-----%%/ t9 G- k# O; z# Z, [
%% system MVA base# s* A3 F5 ^1 A' g
mpc.baseMVA=100;
7 A5 e V- D, ? N' q8 h, L%% bus date2 Z9 G! d* Q# O6 a1 M& ^4 b
% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin% |) F& t9 ~1 ?7 m" ]& ^
mpc.bus=[# ~/ u" P4 ]" z* h6 O
1 1 160 80 0 0 1 1 0 100 1 1.1 0.94;1 p% a7 d4 [! I; K- r( F9 D( @
2 1 200 100 0 0 1 1 0 100 1 1.1 0.94;
7 U" x5 Q/ C0 W+ Z8 w 3 1 370 130 0 0 1 1 0 100 1 1.1 0.94;
* g q) \: y- x' m+ U2 T 4 2 0 0 0 0 1 1.050 0 100 1 1.1 0.94;" Y! q) u1 o5 }/ |% h
5 3 0 0 0 0 1 1.050 0 100 1 1.1 0.94;
@: V1 O+ D' F5 k( } ];9 y; O6 T, ^7 d) N0 m
%% generator date
9 c/ n! \/ P$ w/ o P% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin 1 @1 A) v1 \6 _ u) M3 W
mpc.gen =[& C+ s& q& O7 x6 I4 p
4 500 0 99990 -9999 1.050 100 1 600 0;
1 p7 W) X4 E" Z 5 0 0 99990 -9999 1.050 100 1 600 0;4 X7 D4 E6 n" C- I, m0 E7 {' B8 X2 S
];
& z @: e+ M( J: I" S! s%% branch date
a0 R* U9 l6 r# N J. f% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax
4 C6 y8 j2 K6 |+ N0 q `* Gmpc.branch=[
' c% m! U. J9 A2 Y 2 1 0.04 0.25 0.5 0 0 0 0 0 1 -360 360;+ G" c1 w/ ?! l2 B D. x! T
3 1 0.1 0.35 0 0 0 0 0 0 1 -360 360;
, }( B8 ~% ~0 a& W9 ^4 \0 Q 3 2 0.08 0.3 0.5 0 0 0 0 0 1 -360 360;
- ~1 J. v8 G) C 3 5 0 0.03 0 0 0 0 1.05 0 1 -360 360;
1 T9 H& C c. Y5 n$ z 2 4 0 0.015 0 0 0 0 1.05 0 1 -360 360;
1 e6 L9 S3 F4 x5 C ];
9 H9 M% x& Z' Yreturn; |
-
运行结果
|