|
|
楼主 |
发表于 2013-9-15 11:11:54
|
显示全部楼层
一下是代码:
- h+ L5 t9 w4 E
( l, g" M, s- tfunction mpc=case5_01
' y: q9 @& A, T- r" a9 z# K% MATPOWER Case Format:Version 2; S$ Q4 |, U+ |7 b- [7 Y( b' n( E
mpc.version='2';
% r2 S3 ^; y7 b8 k6 @%%-----Power Flow Date-----%%
7 D1 h4 A& B+ A" }" n%% system MVA base
; W3 v* W0 z% c1 Cmpc.baseMVA=100;
. \6 k5 f L u* @* }%% bus date
8 s0 _" s" R8 l8 u% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin/ d9 c0 ^2 Y+ D, l, a3 @! Y
mpc.bus=[
7 t/ N1 l$ j, J9 [: q! C/ O2 w 1 1 160 80 0 0 1 1 0 100 1 1.1 0.94;: B, |5 ~% r3 J
2 1 200 100 0 0 1 1 0 100 1 1.1 0.94;
5 x3 l! _2 Q( F 3 1 370 130 0 0 1 1 0 100 1 1.1 0.94;
8 N" f0 m1 x. l5 y/ t, R s 4 2 0 0 0 0 1 1.050 0 100 1 1.1 0.94;
( h( s+ o0 g/ l- d, P 5 3 0 0 0 0 1 1.050 0 100 1 1.1 0.94;2 x6 I& S) P1 V+ i7 A) U( g. e! @) l; f- o
];$ M! f* d4 \; u3 X' G! A2 b
%% generator date
3 X- E7 [4 L9 I0 p. \1 u% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin ( [9 |& s% H: C0 r) ]3 y+ H1 j( y
mpc.gen =[& r( p4 t9 N6 E9 W) k
4 500 0 99990 -9999 1.050 100 1 600 0;( R# Y- B. P* ~9 i: k9 V
5 0 0 99990 -9999 1.050 100 1 600 0; }/ e2 y2 {1 L. k
];
$ L% @- v( c2 |%% branch date
% w* u+ B y" F( d( V' {% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax
z: F0 u' X; H8 b- ^: F |mpc.branch=[% c- m0 h( U) e3 L5 ?
2 1 0.04 0.25 0.5 0 0 0 0 0 1 -360 360;) i! D) o% j- a$ @$ j$ L2 Y5 u" |4 G
3 1 0.1 0.35 0 0 0 0 0 0 1 -360 360;
% C& L- f) w9 F. Z: r1 p6 k 3 2 0.08 0.3 0.5 0 0 0 0 0 1 -360 360;
3 \0 K/ L1 B6 m! ]- J 3 5 0 0.03 0 0 0 0 1.05 0 1 -360 360;
5 n$ v2 E3 _: g# ?+ J1 A 2 4 0 0.015 0 0 0 0 1.05 0 1 -360 360;* S! D( l4 p* x3 B
];2 u/ N' A# A% U4 e5 W# {# K3 m
return; |
-
运行结果
|