|
|
楼主 |
发表于 2013-9-15 11:11:54
|
显示全部楼层
一下是代码:
1 S' |6 F9 j8 d9 Q- j0 ^8 m' R) n1 T
function mpc=case5_01
! [* c: p+ y3 F; l4 R% MATPOWER Case Format:Version 2
2 B r9 S( @2 a; g* k1 M$ fmpc.version='2';: `( L4 u4 v' N6 r" D1 p( ] d
%%-----Power Flow Date-----%%7 r: T O$ j$ ?! g5 a5 i
%% system MVA base* H% q9 F% c8 z+ m$ L8 G0 v% k2 x
mpc.baseMVA=100;
. z3 A6 E' x# h, ~( |%% bus date' Z1 i& u- T* a, R2 Q+ l( T/ e
% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin" }4 V5 B7 J: _: |0 \
mpc.bus=[
' Q, q* ]2 q1 I* v) }3 Q, `6 j 1 1 160 80 0 0 1 1 0 100 1 1.1 0.94;
4 u" T, K/ O3 K }( K5 E 2 1 200 100 0 0 1 1 0 100 1 1.1 0.94;. X! Y! z* V2 J' x1 E, \
3 1 370 130 0 0 1 1 0 100 1 1.1 0.94;' i$ |: r) O9 L( N M0 Y( j
4 2 0 0 0 0 1 1.050 0 100 1 1.1 0.94;
7 o# R3 T9 u9 g# |8 T/ K; j' C( G 5 3 0 0 0 0 1 1.050 0 100 1 1.1 0.94;' ^6 q+ w; J& B$ d0 C& {, p
];
" k9 ?; T7 {3 _6 i! D! J1 F%% generator date7 O7 Y/ m, j5 }! M: v) `
% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin / {' w2 q7 c: g* O# b2 U! n
mpc.gen =[
6 P0 x; X- ^. h$ c 4 500 0 99990 -9999 1.050 100 1 600 0;
! Q* ^2 y5 e. n8 r3 Q 5 0 0 99990 -9999 1.050 100 1 600 0;
! X) z' b* a2 e( w- Y* w8 y; ] ];
- m9 x6 q) T* W% r. c% y%% branch date
( r) i |' g3 [) t. l$ ~% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax ; {# X8 W4 ?7 h
mpc.branch=[
3 w( \# _2 z: @$ K 2 1 0.04 0.25 0.5 0 0 0 0 0 1 -360 360;7 D; y9 } z: i) ]
3 1 0.1 0.35 0 0 0 0 0 0 1 -360 360;6 {: N4 _1 D$ I9 K I6 N
3 2 0.08 0.3 0.5 0 0 0 0 0 1 -360 360;
& u$ D1 a% R& W. u2 A 3 5 0 0.03 0 0 0 0 1.05 0 1 -360 360;' v$ m ?8 w' l* S: f5 c
2 4 0 0.015 0 0 0 0 1.05 0 1 -360 360;& e+ s- G7 ]+ w
];3 F% y3 d# m6 I. V% N _
return; |
-
运行结果
|