|
|
楼主 |
发表于 2013-9-15 11:11:54
|
显示全部楼层
一下是代码: Z5 K4 j4 T+ N% n
0 Z/ h. ^! s+ E, g' B3 r* I
function mpc=case5_01" g9 K' I J& M8 P
% MATPOWER Case Format:Version 25 F, K3 V4 ]2 B# e5 }: ~5 h' [
mpc.version='2';
+ h' s! L7 g4 |! V p%%-----Power Flow Date-----%%# t) p2 w, R/ U
%% system MVA base; e$ b; d% i8 f1 G
mpc.baseMVA=100;
6 L' n X C/ ~0 Q, e$ S; S%% bus date8 b. |8 f$ J2 S: Y2 |$ i" X
% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin
' V J2 z* ?6 m: J, R% _! e. Bmpc.bus=[% l4 X/ I0 @+ q. `
1 1 160 80 0 0 1 1 0 100 1 1.1 0.94;" I7 o6 Z4 T# @4 y2 @6 c/ Q
2 1 200 100 0 0 1 1 0 100 1 1.1 0.94;
, q9 F& i2 P- E: N; ~+ Z 3 1 370 130 0 0 1 1 0 100 1 1.1 0.94; D1 M2 x# z' Y5 _
4 2 0 0 0 0 1 1.050 0 100 1 1.1 0.94;5 R7 y- |$ p- x! I" D
5 3 0 0 0 0 1 1.050 0 100 1 1.1 0.94;. V6 R* j% I. E% n- m- ~& g; A
];% g5 J0 F b& i/ r* a% _
%% generator date
; ?/ j A( h! @9 m; {4 K# a' g8 Z* s% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin 5 [* K, l3 Y4 k; u
mpc.gen =[
' b1 I4 M+ V- D: s1 F( k 4 500 0 99990 -9999 1.050 100 1 600 0;" H( P- U/ S- ^9 B, }4 i
5 0 0 99990 -9999 1.050 100 1 600 0;6 t# E/ ^6 ]- s7 ~4 u# O0 Z
];& n6 X( ]- [# S2 _5 r1 p \
%% branch date9 I( w7 D- D9 b k
% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax ; w; O- e2 a( U4 J
mpc.branch=[3 }1 r& E7 C3 ~ n/ Z2 R# n
2 1 0.04 0.25 0.5 0 0 0 0 0 1 -360 360;
s1 H1 \; Q$ l$ O 3 1 0.1 0.35 0 0 0 0 0 0 1 -360 360;
) C" M8 z) G, i! n d 3 2 0.08 0.3 0.5 0 0 0 0 0 1 -360 360;
- x, t7 v g) P 3 5 0 0.03 0 0 0 0 1.05 0 1 -360 360;
+ K& q6 ~! g$ c1 G- R( m 2 4 0 0.015 0 0 0 0 1.05 0 1 -360 360;3 G6 B/ @0 @8 }0 a2 t8 }/ `
];
; I1 r1 E8 p/ S, U3 preturn; |
-
运行结果
|