|
|
楼主 |
发表于 2013-9-15 11:11:54
|
显示全部楼层
一下是代码:' m$ y8 {( L9 n' v; Y- P
& b$ g% ` H @) `" nfunction mpc=case5_013 Z; g, g4 H0 i
% MATPOWER Case Format:Version 2
; O8 P% D! D' r- ^mpc.version='2';4 k9 o; s! U1 K+ a$ H
%%-----Power Flow Date-----%%
4 ~, u, L G$ C2 G8 p d%% system MVA base! n$ {8 B& {$ H( Z8 ]# W, V, H
mpc.baseMVA=100;: B4 o& ~0 e: W8 }$ O( _; ?
%% bus date; A' c) s' f. E2 b- Y) s
% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin" W. e! c" j; F' ?( i) e3 K/ r, `
mpc.bus=[
* f0 @6 h3 ^& D$ V: t$ b/ D0 C 1 1 160 80 0 0 1 1 0 100 1 1.1 0.94;
$ a: K+ W) ~* F/ p7 h; I 2 1 200 100 0 0 1 1 0 100 1 1.1 0.94;5 x1 q8 F# p8 w
3 1 370 130 0 0 1 1 0 100 1 1.1 0.94;9 _0 O( ?; y1 H1 `$ V
4 2 0 0 0 0 1 1.050 0 100 1 1.1 0.94;' k' v P1 m* y e
5 3 0 0 0 0 1 1.050 0 100 1 1.1 0.94;
. Y! ?- Z1 y U* M! q ];
! m' ?# T% a2 B( ~& F%% generator date
+ W8 ]* a" L; p4 V+ d8 N% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin 5 S3 x, j. O' |6 G8 {7 q# S) D
mpc.gen =[$ ]- o7 O- W' |; v. b9 K& \, K
4 500 0 99990 -9999 1.050 100 1 600 0;
) Q; M; G4 }& C( o 5 0 0 99990 -9999 1.050 100 1 600 0;
# L2 W; k/ C4 t# D4 y9 y* x ];9 f, I" ~8 c, X
%% branch date
# J0 b) o6 B1 r( X. F- O" J7 R% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax
$ f. P* b. M+ V( o3 @& Impc.branch=[
% ~1 b5 K0 z) W3 k9 g 2 1 0.04 0.25 0.5 0 0 0 0 0 1 -360 360;3 z6 y- u: V8 C- x' r8 J
3 1 0.1 0.35 0 0 0 0 0 0 1 -360 360;
9 `' H% p2 j( t* `& @5 T x4 q7 t 3 2 0.08 0.3 0.5 0 0 0 0 0 1 -360 360;
! U% {( c( {1 f: }. h" x& v 3 5 0 0.03 0 0 0 0 1.05 0 1 -360 360;
% r" F( V8 {. a6 i( Q! K' _ N 2 4 0 0.015 0 0 0 0 1.05 0 1 -360 360;& U2 [' m) E5 l1 v- l5 R
];
& H) d. D. G6 Kreturn; |
-
运行结果
|