|
楼主 |
发表于 2013-9-15 11:11:54
|
显示全部楼层
一下是代码:
+ s x% o( E" i$ Y% l/ |4 k+ `$ F
! h& N) s# B& s/ V& bfunction mpc=case5_01
4 k0 U; H5 p/ e$ _3 _, N3 V$ b% MATPOWER Case Format:Version 21 }3 a5 I1 w: x
mpc.version='2';6 e# o3 t" P( {4 d o
%%-----Power Flow Date-----%%5 v# i+ P0 x; H2 G- X
%% system MVA base8 |! S5 D6 S0 N2 E# b$ W0 x
mpc.baseMVA=100;, d, D, ?- D- o$ a- F
%% bus date
/ \) k! e9 u- q% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin5 p7 n0 e5 x9 A' C; M9 s
mpc.bus=[
$ G6 |1 s( `" @+ G$ G 1 1 160 80 0 0 1 1 0 100 1 1.1 0.94;7 Z/ _. Q" k5 p- d2 \' q/ y
2 1 200 100 0 0 1 1 0 100 1 1.1 0.94;& c0 F! V0 Z/ D# G
3 1 370 130 0 0 1 1 0 100 1 1.1 0.94;
9 _; v- u; i3 T; C8 k 4 2 0 0 0 0 1 1.050 0 100 1 1.1 0.94;" e6 m: J) {7 l, o
5 3 0 0 0 0 1 1.050 0 100 1 1.1 0.94; S2 d- M( W! V9 ~1 J
];
' Y, p; w8 h& l& G* K. ?$ ~' U%% generator date
/ `1 ^! Z9 k& k5 ^% G" }" x% ]% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin # R8 H& U- D" l6 I( v5 t
mpc.gen =[! e W( V: e0 b& q
4 500 0 99990 -9999 1.050 100 1 600 0;
' r( V' V. Y/ _9 ~ 5 0 0 99990 -9999 1.050 100 1 600 0;
9 ~- O1 G" l% X ];
2 a/ Q6 v+ ~& ~8 }* r/ u) L%% branch date" O+ y: q& H) `8 c
% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax 6 p5 O, B& J/ F4 H4 c8 S3 t
mpc.branch=[+ U& I7 }6 h* h
2 1 0.04 0.25 0.5 0 0 0 0 0 1 -360 360;
1 a, ]8 t! [2 l3 B' z 3 1 0.1 0.35 0 0 0 0 0 0 1 -360 360;: P! Y$ t' c& d* O
3 2 0.08 0.3 0.5 0 0 0 0 0 1 -360 360;" {7 ]6 e+ _, \# G1 Q. T/ h1 V
3 5 0 0.03 0 0 0 0 1.05 0 1 -360 360;
9 z: U; A, Z) n7 \* g 2 4 0 0.015 0 0 0 0 1.05 0 1 -360 360;2 [" }/ B) G: g& u% S* U
];
6 K6 {7 H/ r, J; P! q$ X3 Areturn; |
-
运行结果
|