|
|
楼主 |
发表于 2013-9-15 11:11:54
|
显示全部楼层
一下是代码:
. ` Y2 o0 {- x6 f6 l5 C
' _2 ~2 ]; B* }3 s$ V( Ffunction mpc=case5_01' O3 p% {/ C, G1 P
% MATPOWER Case Format:Version 2
2 {0 D' K: ~9 N) Y: @- p4 Zmpc.version='2';
6 H9 [( T" x+ ]; Q/ Y* I" s1 G5 N8 T%%-----Power Flow Date-----%%& x, R2 l+ o: _+ ]' _# Y
%% system MVA base- x$ Y- P1 e7 K/ Z. o
mpc.baseMVA=100;( H/ e, y* f7 s# Y: M: D! F
%% bus date# @2 v" J% V$ q
% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin! W/ V1 }! J6 D( R4 U
mpc.bus=[7 _* x# B1 l C: L- H( Z
1 1 160 80 0 0 1 1 0 100 1 1.1 0.94;
* L8 s- H s- J9 N$ p 2 1 200 100 0 0 1 1 0 100 1 1.1 0.94;
6 h6 a. _8 ?3 |/ N 3 1 370 130 0 0 1 1 0 100 1 1.1 0.94;
# q- T m, C. @: `/ h) F& c4 `9 n 4 2 0 0 0 0 1 1.050 0 100 1 1.1 0.94;
" g- k% {" n' P5 x 5 3 0 0 0 0 1 1.050 0 100 1 1.1 0.94;1 a2 r0 {% E& y9 W8 Z3 t
];
& ?, Z! F7 @# r* i# \: T%% generator date! G! L# H: a" P1 ]7 f& Y- n' {- R
% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin
: A: v& ~' S) _5 Smpc.gen =[. _5 l5 b4 i8 I7 v8 W7 @9 y, p
4 500 0 99990 -9999 1.050 100 1 600 0;
/ r/ k; F( e% d/ b( g8 U- a0 J 5 0 0 99990 -9999 1.050 100 1 600 0;' i! ~$ m& l: \. M
];
4 r( ?7 ?7 s# X$ Q* b%% branch date
2 }5 ^% }$ s) m2 `0 P( K% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax ' v6 z$ F( P( H2 b) q
mpc.branch=[
! ]3 I! [+ D7 N5 G 2 1 0.04 0.25 0.5 0 0 0 0 0 1 -360 360;
* h' D6 k: H! _* L" Y 3 1 0.1 0.35 0 0 0 0 0 0 1 -360 360;
- ]+ [" K4 a% m9 F: N 3 2 0.08 0.3 0.5 0 0 0 0 0 1 -360 360;' q. E! Z4 w6 i7 `$ C9 Q4 I( k
3 5 0 0.03 0 0 0 0 1.05 0 1 -360 360;, {# j1 z. e8 [
2 4 0 0.015 0 0 0 0 1.05 0 1 -360 360;) f8 a* K( B4 J, k
];& `8 H+ }. B" e: b; T4 W# g
return; |
-
运行结果
|