|
|
楼主 |
发表于 2013-9-15 11:11:54
|
显示全部楼层
一下是代码:
4 i5 O, H- A. [8 x/ ^6 l* D# G
0 ]# U3 W/ f# K3 F" [0 M9 r0 Z$ Gfunction mpc=case5_017 Q2 a& _+ C0 O9 g
% MATPOWER Case Format:Version 2 U; K+ ?! {0 e
mpc.version='2';
2 Y6 Z/ z# Z" b4 i$ b( U+ {4 Z%%-----Power Flow Date-----%%
8 T/ [3 d. m5 n% E8 f%% system MVA base. c$ e0 @' z& J) J& d
mpc.baseMVA=100;& |2 |6 a8 `' \) R9 B
%% bus date1 o' M# E6 Q1 U6 {8 e( \( A
% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin
8 W. ]' E4 ^; m& dmpc.bus=[
3 K) j I0 [; d$ ]! j 1 1 160 80 0 0 1 1 0 100 1 1.1 0.94;, }9 p0 ?+ R) k) D/ S
2 1 200 100 0 0 1 1 0 100 1 1.1 0.94;
" Y% }8 [8 ~$ A* w6 o 3 1 370 130 0 0 1 1 0 100 1 1.1 0.94;& H! Y8 j K- A5 u
4 2 0 0 0 0 1 1.050 0 100 1 1.1 0.94;4 L! O! ]4 t" w% O$ \' J$ c
5 3 0 0 0 0 1 1.050 0 100 1 1.1 0.94;
; @* Y5 m. N/ b# F5 O+ n ];
2 c3 p7 S8 D1 X, M%% generator date4 c# q1 I/ K6 B# y0 ]# K2 v0 H
% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin
7 [# t. O, q! |% ~mpc.gen =[. y, d6 j( ? l
4 500 0 99990 -9999 1.050 100 1 600 0;/ j. J1 p' U; s& e( K9 H& K
5 0 0 99990 -9999 1.050 100 1 600 0;6 l3 y7 G7 ~7 m s5 j
];1 Q7 \0 S9 a: j2 Y: W; j
%% branch date: q0 X9 H3 T; V7 I
% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax / R9 X% s. A$ v8 `0 P/ e
mpc.branch=[( i! c+ Q4 N. o. s9 q7 L& F- Z
2 1 0.04 0.25 0.5 0 0 0 0 0 1 -360 360;4 C$ h; W x+ Q3 i) F
3 1 0.1 0.35 0 0 0 0 0 0 1 -360 360;
9 l3 U( I/ H/ x h( | 3 2 0.08 0.3 0.5 0 0 0 0 0 1 -360 360;
2 P# F4 x0 c9 r3 R: }: m 3 5 0 0.03 0 0 0 0 1.05 0 1 -360 360;8 T5 V; B ?2 E3 h
2 4 0 0.015 0 0 0 0 1.05 0 1 -360 360;2 E3 K/ C- G5 j( ~$ b- ~
];
% |0 w% P! H1 m- m# z! Ireturn; |
-
运行结果
|