|
|
楼主 |
发表于 2013-9-15 11:11:54
|
显示全部楼层
一下是代码:
( |( y7 E( P S$ p) v# l
$ B1 s4 [8 \; e) Kfunction mpc=case5_01
' j4 b! q3 O: \9 k/ ?/ W4 f% MATPOWER Case Format:Version 2' H& O& e7 t2 d/ P' R
mpc.version='2';
6 ]0 \) O; P1 M: r%%-----Power Flow Date-----%%
( k$ F: r' K& D( v2 J( x6 V%% system MVA base! K3 ~ [7 i2 W3 s, z& ?
mpc.baseMVA=100;! j/ ?: ~ D5 w
%% bus date
. E+ o" Y+ q) R# V6 F+ K% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin/ C0 [. K! ]: y/ d& L( n% J
mpc.bus=[# j9 d, p" Z9 a p4 D6 z
1 1 160 80 0 0 1 1 0 100 1 1.1 0.94;
$ Z% ^ ?* c" l, A' T- L5 l7 A5 q) @ 2 1 200 100 0 0 1 1 0 100 1 1.1 0.94;
$ t7 L9 F @1 m: e" k 3 1 370 130 0 0 1 1 0 100 1 1.1 0.94;4 O9 t5 r8 D, i9 A0 w" @, L% w
4 2 0 0 0 0 1 1.050 0 100 1 1.1 0.94;% w- G) c7 p, g& V: W/ p# n
5 3 0 0 0 0 1 1.050 0 100 1 1.1 0.94;
s2 z& w# K2 p" R$ j( M ];$ }$ }+ _5 f! j1 V* ]! M) L
%% generator date3 f& O( \# s/ v5 A- }
% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin
- t7 D% B1 Z: L! |+ s0 o" C& ?7 nmpc.gen =[" H) Q6 |& A7 D2 b4 K
4 500 0 99990 -9999 1.050 100 1 600 0;
/ U- O" [7 {5 W5 f: [! H* a 5 0 0 99990 -9999 1.050 100 1 600 0;
/ J. }/ [0 U% ?: \7 k6 M8 ?5 l( s: Q- ? ];
! w0 B: F. o' }+ [8 F9 g% F%% branch date3 Y- @ L+ M( A; e6 x/ z
% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax
/ |% Z) I) T' o0 W/ umpc.branch=[
& h8 _" f$ C% F, u 2 1 0.04 0.25 0.5 0 0 0 0 0 1 -360 360;
3 G1 Q* n+ Z' a 3 1 0.1 0.35 0 0 0 0 0 0 1 -360 360;3 z Y9 i, i( U; m: j" V# z
3 2 0.08 0.3 0.5 0 0 0 0 0 1 -360 360;
y7 {& G4 i! E" @" h 3 5 0 0.03 0 0 0 0 1.05 0 1 -360 360; i: {: D' _0 @ a
2 4 0 0.015 0 0 0 0 1.05 0 1 -360 360;
. G. p2 s/ |0 V: b l ];8 k" B7 ~! o2 t4 V1 K. `
return; |
-
运行结果
|