|
楼主 |
发表于 2013-9-15 11:11:54
|
显示全部楼层
一下是代码:/ N2 V2 J) G4 y( V
: c( z! _! Y$ a( f2 j# P5 d
function mpc=case5_01) C% _+ \7 X4 s) N
% MATPOWER Case Format:Version 29 g1 M- M# r9 l! x
mpc.version='2';
% v8 S! I% c! S( U; L+ ~. L%%-----Power Flow Date-----%%
$ Q) s. v2 n8 S1 Y%% system MVA base
2 ~& H$ R6 U6 S" S9 ]' Tmpc.baseMVA=100;( |) |$ h# L7 A2 w( G/ u, J
%% bus date
7 b$ F* ^5 W2 M6 P7 P* L7 A% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin
0 t) J* r- j/ ?8 ~( A8 ~- P8 Hmpc.bus=[. ^6 d# e, A" m9 d
1 1 160 80 0 0 1 1 0 100 1 1.1 0.94;
2 Y; |1 {6 X3 d, X 2 1 200 100 0 0 1 1 0 100 1 1.1 0.94;
0 Q' Y, N9 [* t/ R7 M0 l 3 1 370 130 0 0 1 1 0 100 1 1.1 0.94;
9 t% m. v* e; q* \, R# L- k 4 2 0 0 0 0 1 1.050 0 100 1 1.1 0.94;2 t9 G: m: S9 g
5 3 0 0 0 0 1 1.050 0 100 1 1.1 0.94;/ x2 _1 v, [$ P
];
+ v5 r+ u; D6 C5 x$ E: n; t+ H%% generator date
6 V. B5 W! d6 k1 Q% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin
9 V: ~' V( P& z4 G$ r3 c6 n2 }. gmpc.gen =[; g8 m- l/ ~$ y1 ^# K1 Z6 t
4 500 0 99990 -9999 1.050 100 1 600 0;+ V+ T# S7 o7 \8 J! |. @) H
5 0 0 99990 -9999 1.050 100 1 600 0;
+ n# P" b. q/ o( V3 d5 J4 y! W3 E+ R ];, R' I8 d2 c) t! `# f
%% branch date
3 C; H/ c2 s9 F" t% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax - E I4 {/ M& ^, S2 M
mpc.branch=[
" x5 B- e) ~8 g 2 1 0.04 0.25 0.5 0 0 0 0 0 1 -360 360;
6 o6 w7 y: E1 t3 R$ L/ \$ u. B 3 1 0.1 0.35 0 0 0 0 0 0 1 -360 360;5 @2 c2 t" G: l; J* [% o
3 2 0.08 0.3 0.5 0 0 0 0 0 1 -360 360;9 t' Q: Y: m3 \" ~* [9 D
3 5 0 0.03 0 0 0 0 1.05 0 1 -360 360;
3 L; e- H/ P% m 2 4 0 0.015 0 0 0 0 1.05 0 1 -360 360;$ ^, q2 q! }' b }/ f& ]0 _
];
$ H$ u% a) B0 Z. q dreturn; |
-
运行结果
|