|
|
楼主 |
发表于 2013-9-15 11:11:54
|
显示全部楼层
一下是代码:
0 z: F& [3 r( d9 `: a, q( Q5 I# j! h. J! u8 P1 b" C
function mpc=case5_01% u1 u, k3 y; |3 \ e! v
% MATPOWER Case Format:Version 2- _% `" Y( l/ z; K+ n% z1 t4 S
mpc.version='2';
1 G. l7 z+ ~0 E+ F" }%%-----Power Flow Date-----%%
- b5 P4 x: l4 ?& R2 e) K2 T# M6 p3 j%% system MVA base
8 G0 S( a' d# Z" l% tmpc.baseMVA=100;/ B. F' V0 W' d
%% bus date
8 ]# C/ Z1 F+ a5 _% p% X% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin
: l; U, Q5 I+ d7 C \, k+ m2 O# Gmpc.bus=[9 U- m9 X* Q! c' g% y) j3 t
1 1 160 80 0 0 1 1 0 100 1 1.1 0.94;0 m3 A8 O8 m3 B9 e
2 1 200 100 0 0 1 1 0 100 1 1.1 0.94;. g- ?7 y% O. U
3 1 370 130 0 0 1 1 0 100 1 1.1 0.94;8 R) F$ S1 m6 _; @
4 2 0 0 0 0 1 1.050 0 100 1 1.1 0.94;& {0 A) t5 W/ C. n
5 3 0 0 0 0 1 1.050 0 100 1 1.1 0.94;' o. A! ]. N5 m7 o. x7 e
];/ ~& E1 P% A* o
%% generator date# j1 w$ ]; k1 g% I) j- `
% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin ) a7 K' {& m! u% W& b/ A; C2 ~
mpc.gen =[8 v: R ^( V+ Q. h* f5 x+ t! @
4 500 0 99990 -9999 1.050 100 1 600 0;
6 N' h( H# [" B. F/ d7 U+ E0 g' } 5 0 0 99990 -9999 1.050 100 1 600 0;
6 c' \6 o* ~" p- j: M1 a- W) C ];: _0 _2 o% v" k
%% branch date
1 [) j! b+ `& k$ L% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax
! K# M4 O: u' n+ W! Pmpc.branch=[, r6 o' B! _, @" Z. Q3 q3 L
2 1 0.04 0.25 0.5 0 0 0 0 0 1 -360 360;
6 w/ j; |9 G8 H0 b4 i! [ 3 1 0.1 0.35 0 0 0 0 0 0 1 -360 360;0 k- D& a0 G. g; ^
3 2 0.08 0.3 0.5 0 0 0 0 0 1 -360 360;
; M1 _* _+ {' Y' G* E2 c 3 5 0 0.03 0 0 0 0 1.05 0 1 -360 360;
' c) G# v2 H" ~ 2 4 0 0.015 0 0 0 0 1.05 0 1 -360 360;
3 E, ~0 a% }0 s( `3 G+ q6 b ];
' K5 }/ B4 E+ |4 `! U5 s3 creturn; |
-
运行结果
|