|
|
楼主 |
发表于 2013-9-15 11:11:54
|
显示全部楼层
一下是代码:
$ D @% G8 ?" L$ i2 I- l
8 j: B; M% a7 u) m. Qfunction mpc=case5_01
3 E* o# D" c: _3 ~2 f! K, `) T+ G: p% MATPOWER Case Format:Version 2+ M5 X4 ?- ]- G+ Z7 a
mpc.version='2';: \1 R, v0 g. \3 K0 ^
%%-----Power Flow Date-----%%, p4 j( S6 J$ Y5 C/ Z/ ]9 k9 c4 ~
%% system MVA base
$ E* ~( O/ d& P6 G hmpc.baseMVA=100;
! f- ?7 F/ N6 o: \* M: g/ f( _ Z%% bus date
9 R! W3 r; n& z, M% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin
3 U5 |0 d! `" R3 Zmpc.bus=[& q2 M% L7 x4 L& S& O% c# u$ T) J
1 1 160 80 0 0 1 1 0 100 1 1.1 0.94;
! q5 M, |6 k) n* e2 }7 u+ U8 S 2 1 200 100 0 0 1 1 0 100 1 1.1 0.94;- m u% a# [ g" f" v
3 1 370 130 0 0 1 1 0 100 1 1.1 0.94;# {% w7 {; f; C4 k: g* f
4 2 0 0 0 0 1 1.050 0 100 1 1.1 0.94;
. x% A! V V- A* v* d 5 3 0 0 0 0 1 1.050 0 100 1 1.1 0.94;
6 q9 ]# j+ z, y3 G* [% | ];
6 p" R' P) V4 \8 a+ Q* k- g%% generator date
2 @% ^3 W: D8 B! C% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin
/ U6 n7 F7 m% i$ M8 t( q, kmpc.gen =[
) @* t3 o4 X/ }1 P9 x2 X 4 500 0 99990 -9999 1.050 100 1 600 0;
E3 M$ ~- ]/ t- [$ k* ]1 B 5 0 0 99990 -9999 1.050 100 1 600 0;7 |( M+ C/ l8 c9 T6 u
];( D+ P# y6 S7 A) V) C# l: G
%% branch date9 p/ h+ X0 Q* e4 P4 n: D" M) K7 n
% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax 6 I; {& ~) Y" s2 \8 T# c- S& k
mpc.branch=[
$ W% R/ o& c# j4 Z* m 2 1 0.04 0.25 0.5 0 0 0 0 0 1 -360 360;+ y+ W3 ^8 j3 n6 O N! g. V6 a
3 1 0.1 0.35 0 0 0 0 0 0 1 -360 360;
5 N* p" {+ Q, [1 v1 [ 3 2 0.08 0.3 0.5 0 0 0 0 0 1 -360 360;
* l+ \/ c3 f9 g 3 5 0 0.03 0 0 0 0 1.05 0 1 -360 360;
$ \- Y; S% y% x 2 4 0 0.015 0 0 0 0 1.05 0 1 -360 360;. k6 U: G5 s2 W$ E3 q) Y
];
- f1 p" b: H6 x! A. ~) dreturn; |
-
运行结果
|