|
|
楼主 |
发表于 2013-9-15 11:11:54
|
显示全部楼层
一下是代码:2 r0 u2 @5 i% s
+ m2 s% H2 N2 |4 a/ l; q" {
function mpc=case5_01! u% b1 l9 N9 u$ m
% MATPOWER Case Format:Version 2: }5 S. s' x+ M/ l6 H. u$ ]
mpc.version='2';
$ o* t8 ~0 D2 j. k%%-----Power Flow Date-----%%5 Y; o7 x/ Q$ f E, m8 T) [
%% system MVA base" ^' Y- o7 r M7 {5 N( h- y( R
mpc.baseMVA=100;
% `/ v* _$ W6 y2 T/ G%% bus date
0 Q, R* F2 b3 k* i% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin
8 }: J; m( P) P1 w2 |0 v+ ^mpc.bus=[$ P, ?( F' i3 \ A2 e5 B/ H
1 1 160 80 0 0 1 1 0 100 1 1.1 0.94;. s5 }, N5 W% G k
2 1 200 100 0 0 1 1 0 100 1 1.1 0.94;
; W! S9 s; W! _) j/ \ 3 1 370 130 0 0 1 1 0 100 1 1.1 0.94;: l+ L7 q) |& y2 x
4 2 0 0 0 0 1 1.050 0 100 1 1.1 0.94;; O; p/ S, j8 o6 K) f3 Q; a
5 3 0 0 0 0 1 1.050 0 100 1 1.1 0.94;
* ?0 k, U& D6 s1 M1 N: L" ^# P5 g ];
9 r I3 ], S9 Y7 h9 a" P%% generator date( o' m5 G! ]$ `! A5 r- a0 L$ r: T
% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin
y+ p7 j7 |; ?mpc.gen =[3 {" W1 @5 D6 V
4 500 0 99990 -9999 1.050 100 1 600 0;
* `" T' T+ l. S: r: H 5 0 0 99990 -9999 1.050 100 1 600 0;
* _( Y' s$ e+ q; |9 `; ] ];4 |# G2 [4 O- T$ y9 G, k0 i7 P9 u7 q
%% branch date+ a. w" e+ O3 ?( x/ ?9 k
% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax
& F# }2 g: m o0 K4 @mpc.branch=[
4 j) X' X2 j, k+ W/ u4 u8 C& y 2 1 0.04 0.25 0.5 0 0 0 0 0 1 -360 360;
; @2 [2 F& Y! [# f. j* G' `( l 3 1 0.1 0.35 0 0 0 0 0 0 1 -360 360;* Y6 u1 d* v) u
3 2 0.08 0.3 0.5 0 0 0 0 0 1 -360 360;) d" O0 F6 A% ~6 {: G
3 5 0 0.03 0 0 0 0 1.05 0 1 -360 360;
) n$ F; _, ~3 U+ v M _$ a5 a8 g7 S7 [ 2 4 0 0.015 0 0 0 0 1.05 0 1 -360 360;
$ l* E& Y7 p5 c2 p v ];
% J# \% B8 Q7 P$ greturn; |
-
运行结果
|